1. Introduction
Gallium-aluminium nitride alloys may be promising materials for optical applications, particularly for light emission in the ultraviolet range. Indeed, the forbidden band of these wide gap semiconductor alloys seems to be suitable to many applications in this spectral range. Moreover, the residual deformation in these layers, due to the lattice mismatch with the substrate (generally sapphire), may be changed by the aluminium content of the alloy.
Basic studies on these ternary solid solutions were recently developed. The literature concerning Raman spectra from Ga1−XAlXN crystals is poor: we note only the paper published by Hayashi et al. Reference Hayashi, Itoh, Sawaki and Akasaki[1],who found a one-mode behaviour for the polar phonons in a very small compositional range (x≤0.15). In this communication, we present a report on the dynamical properties of Ga1−XAlXN crystals (0<x≤1) investigated by Raman spectroscopy with a special attention to the behaviour of long wavelength polar phonons.
2. Samples and experiments
Eight Ga1−XAlXN layers (with 0.16≤x≤0.84), together with GaN layers, have been used for the present study. The 2 μm thick samples were grown by MOVPE on thin AlN buffer layers previously deposited on a (0001) sapphire substrate; the aluminium content x of the layers was determined from a measurement of the lattice constant Reference Clur, Briot, Rouvière, Andenet, Le Vaillant, Gil, Aulombard, Demangeot, Frandon and Renucci[2]. Several AlN layers (x = 1) were also available.Unfortunately, they did not yield any useful results.
All the spectra have been recorded at room temperature in a backscattering geometry, along or perpendicular to the c-axis of the wurtzite crystal. Micro-Raman measurements with a lateral resolution of about 1 μm, allowing incidence on the layer edge, were performed using a Dilor set-up. The light source was the 488 nm line of a Ar+ laser.
3. Experimental results
Examples of Raman spectra, recorded in
and configurations, are given in Figure 1a and Figure 1b respectively (x′ and y′ represent two orthogonal axes perpendicular to the c axis). For low aluminium content, sharp phonon lines are clearly evidenced and most of them obey the selection rules for backscattering geometry in wurtzite crystals, as evidenced in Figure 2, for x = 0.45 solid solution. So the vibrational modes can be unambiguously followed in this range. However, it should be pointed out that the E1(LO) phonon is found in forbidden configuration, as it was also observed for pure GaN layers. An additional weak structure peaks in the intermediate range between TO and LO phonons, around 640 cm−1, on the high frequency side of the line corresponding to the E2 phonon.
For x>0.5,the intensity of Raman spectra decreases and the peaks become broader and broader. Moreover, we observe for x>0.7 a complete relaxation of selection rules, probably due to the polycrystalline structure of the layers, rendering the assignment of the experimental features difficult. In spite of the lack of unambiguous signature in the whole range of composition, a plot of the measured phonon frequencies in the whole composition range is reported in Figure 3. For the extrapolation up to x=1, we used the experimental phonon frequencies in AlN crystals from McNeil et al. Reference Mcneil, Grimsditch and French[3], which can be compared to the values recently calculated by Karch et al. Reference Karch, Portisch, Bechstedt, Pavone and Strauch[4].
Our frequency variation is in good agreement with previous results obtained in a narrower compositional range by Hayashi et al. Reference Hayashi, Itoh, Sawaki and Akasaki[1]. The continuous evolution observed in Figure 3 seems to be very satisfactory. Note the crossing of the E1(TO) and E2 modes, which is clearly found at x = 0.35. In spite of its uncertain signature, the wide feature evidenced near 660 cm−1 for high aluminium content was assigned to an E2 mode, because it is very close to the corresponding mode in AlN Reference Mcneil, Grimsditch and French[3]. Our data seem to support one-mode behaviour for the LO and TO vibration modes. On the contrary, our results concerning the E2 phonon are quite ambiguous: weak and wide features (which did not clearly obey the usual selection rules) in the spectra might suggest a possible two mode behaviour.
The best fit of our data gives these analytical formulae for the dependence of the A1(LO) and A1(TO) peaks position on Al content, valid in the case where x>0.05 :
Phonon frequencies values are given in cm−1.
Concerning the measured phonon frequencies, the experimental uncertainty is 1 or 2 cm−1 up to x = 0.55. For x>0.55, this value increases up to 10 cm−1.
In the Raman spectra of the solid solutions, a weak dip was observed (see Figure 1) for the whole range of composition, for parallel polarizations of incident and scattered light. Its origin seems to be related to the alloy and not to any crystalline defect, because its frequency is increasing linearly with the aluminum content (see Figure 3). This feature suggests an interference between a discrete mode and a continuum of excitations.
The discrete excitation coupled to the hypothetical continuum might be the “silent” B1 mode in the ternary solid solutions, according to a linear frequency variation between the phonon frequencies of the pure binary compounds, calculated in ref. Reference Karch, Portisch, Bechstedt, Pavone and Strauch[4] and Reference AZUHATA, MATSUNAGA, SHIMADA, YOSHIDA, SOTA, SUZUKI and NAKAMURA[5]. Activation of the silent mode in Raman spectra might proceed from a q-dependent process likely to occur in the alloys because of the breakdown of translational symmetry.
4. Presentation of the model
The purpose of the following model Reference GROENEN, CARLES, LANDA, GUERRET, FONTAINE and GENDRY[6] is to calculate the frequency and the intensity of the LO q = 0 polar phonons in the wurtzite Ga1-XAlXN solid solutions, as a function of the aluminium content. It should be noted that the non-polar E2 modes cannot be calculated in the frame of this model.
In the calculation of coupled vibration modes in ternary solid solutions,the first step is the estimation of impurity mode frequencies. The later is usually estimated from the average frequency for q=0 optical phonons and from its frequency shift, incoming from the extensive or compressive strain suffered by the impurity cluster in the host matrix Reference LANDA, CARLES, RENUCCI and Engström[7]. In the present case, the phonon dispersion, as calculated by Karch et al. Reference Karch, Portisch, Bechstedt, Pavone and Strauch[4] for AlN and by Azuhata et al. Reference AZUHATA, MATSUNAGA, SHIMADA, YOSHIDA, SOTA, SUZUKI and NAKAMURA[5] for GaN, was taken into account in the estimation of impurity modes, to obtain better agreement with experimental results. Due to the concavity of the LO branches at the centre of Brillouin zone, both impurity modes were shifted towards lower frequencies: we have obtained 600 cm−1 and 655 cm−1, respectively, for Ga in AlN and for Al in GaN (instead of the values 631 cm−1 and 684 cm−1, when only q=0 optical phonons are taken into account).
The second step is to define the “mechanical» frequencies ωTi(x), for each uncoupled oscillator (i=1 or 2, for AlN or GaN respectively) in the solid solution, as a function of x: a linear variation of these frequencies will be assumed, between the frequency of the impurity mode and that of the q=0 TO phonon in the pure crystal.
The third step is to define an average complex dielectric constant, for the coupled oscillators i = 1, 2 in the ternary solid solution, as:
where
and are related to the LO-TO splitting of the oscillators in their corresponding pure crystals. Hence the invariable effective charge of each oscillator in the ternary solution was assumed to be invariable. In the expression of ε(ω,x), we introduce empirical damping constants γi of oscillators.
It is known,from the Hon and Faust theory Reference HON and FAUST[8] for two oscillators coupled by the macroscopic electric field, that Im(−ε−1) gives the expected contribution of the LO phonons to the Raman spectra. So the last step is to calculate this function for 0<x<1. The calculated frequencies of A1(LO) phonons are compared with the experimental values in Figure 4a. Peak values of the calculated function Im(−ε−1), corresponding to the expected intensity of LO phonon lines,are plotted in Figure 4b.
5. Comparison of the calculated results with the experimental data
Good agreement is observed between the frequency measured for the LO phonon (A1 or E1) and the calculated LO branch corresponding to the AlN-like oscillator.The other LO branch corresponds to the weak structure, which is hardly evidenced in the 620-640 cm−1 range for solid solutions with low aluminium content. These experimental observations are completely explained by the very low intensity of the LO phonon for the GaN-like oscillator, compared to the AlN-like oscillator, as calculated by our model and illustrated in Figure 3b. So the so-called one-mode behaviour is found for the LO phonons in GaXAl1−XN. It was previously interpreted in the case of zincblende III-V ternary solid solutions as a two mode-behaviour, with a quasi complete intensity transfer from one oscillator to the other Reference LANDA, CARLES, RENUCCI and Engström[7].
6. Summary
The frequencies of long wavelength optical phonons of Ga1−XAlXN solid solutions have been measured by Raman spectroscopy in the whole compositional range. These data give evidence for a one-mode behaviour of the A1 and E1 modes. The frequencies of the LO phonons were calculated from a model based on a dielectric formalism and were compared to the experimental results. Finally, a Fano interference effect observed in the Raman spectra is discussed.