We have calculated the optical absorption for InGaNAs and GaNSb using the band anticrossing (BAC) model and a self-consistent Green’s function (SCGF) method. In the BAC model, we include the interaction of isolated and pair N levels with the host matrix conduction and valence bands. In the SCGF approach, we include a full distribution of N states, with non-parabolic conduction and light-hole bands, and parabolic heavy-hole and spin-split-off bands. The comparison with experiments shows that the first model accounts for many features of the absorption spectrum in InGaNAs; including the full distribution of N states improves this agreement. Our calculated absorption spectra for GaNSb alloys predict the band edges correctly but show more features than are seen experimentally. This suggests the presence of more disorder in GaNSb alloys in comparison with InGaNAs.
Keywords:Dilute nitride semiconductors; Self-consistent Green’s function; Optical absorption; Band anticrossing model
The substitution of a small fraction x of nitrogen atoms, for group V elements in conventional III-V semiconductors such as GaAs and GaSb, strongly perturbs the conduction band (CB) of the host semiconductor. The band structure of dilute nitride alloys has been widely investigated . We have recently developed  a SCGF approach to calculate the density of states (DOS) near the conduction band edge (CBE) in these alloys.
One way to test the accuracy of this model is to look at optical absorption spectra for dilute nitride samples, where we expect to see features related to the N states present in the samples. The absorption spectrum arises from transitions between valence and conduction band states. It provides knowledge of the energy gap in semiconductors, and also gives significant information about the band structure of materials. Experimental measurements of absorption spectra can be used to benchmark band structure calculations. In this paper, we investigate two different materials: InyGa1-yNxAs1-x, for which the band structure has been widely studied and many of the features are well established and GaN xSb1-x, for which much less information has been reported in the literature.
We consider two different models for the band structure of dilute nitride alloys, firstly a five-level band anticrossing (BAC) model, including the host semiconductor CB and valence bands, isolated N and pair N-N states and, secondly, the linear combination of isolated nitrogen states (LCINS) model [3,4], which allows for interaction between N states on nearby sites as well as inhomogeneous broadening and produces a distribution of N state energies. In the LCINS model, the band structure of the alloys is calculated using a SCGF approach .
For InyGa1-yNxAs1-x alloys, we find that the BAC model reproduces the main features in the absorption spectrum, in agreement with previous work [5,6]. However this model shows some additional features which are related to the N and N-N state energies, reflecting that in the BAC model, we have ignored the actual distribution of localised states. Including the LCINS distribution of N states in InyGa1-yNxAs1-x using the SCGF approach  removes the additional features found in the BAC calculations and gives absorption spectra that are in very good agreement with experimental data.
We then apply our methods to GaN xSb1-x, where much less information is known theoretically and experimentally. The overall width of the optical spectrum can be well fitted by our models for the absorption spectrum. Both the BAC and LCINS models account for the absorption edge of GaNSb alloys, supporting the presence of a band anti-crossing interaction in these alloys. However, the five-level BAC model gives more features than are seen experimentally in the absorption spectrum. Including a distribution of localised state energies, obtained by modifying those calculated for GaNAs, makes the calculated absorption spectra smoother and gives better agreement with experimental data but still shows some discrepancies around the localised state peak energies. These results suggest the presence of more disorder in GaNSb samples than in InGaNAs. This disorder may be due to sample inhomogeneities or due to an intrinsically broader distribution of N states in GaNSb than in InGaNAs.
The remainder of this paper is organised as follows. In the ‘Methods’ Section we first provide an overview of optical absorption calculation, followed by a description of the band structure models used for dilute nitride alloys. The theoretical results for InGaNAs and GaNSb are presented and compared with experiment in the ‘Results and discussion’ Section. Finally, we summarise our conclusions in the last section.
The absorption spectrum α(E) describes the rate of absorption of photons with energy per unit distance and can be described using a ‘single-electron’ approximation. In this approach, the absorption spectrum for allowed transitions between valence band v and conduction band c states is given by [7-10]
where is the photon energy, e and m0 are the electron charge and mass, c is the speed of light and nr is the refractive index. f is the Fermi-Dirac distribution function. Here, we assume a filled valence and empty conduction band, so that fv-fc=1. The matrix element in Equation 1 can be written for transitions between valence p- and conduction s-like zone centre states as 
where is the interaction energy, and the momentum interband matrix element, p, can be estimated from experiment as
where is the conduction band effective mass, Eg is the band gap (between conduction and valence band) and Δso is the spin-orbit-splitting energy. Therefore, Equation 1 can be written as
where E=Eck-Evk is the transition energy between conduction (Eck) and valence (Evk) states with wavevector k, Jcv(E) is the joint density of states, and we ignore for now the energy dependence of Mb.
The hole-electron interaction can be included assuming Elliot’s theory, which applies to parabolic and nondegenerate bands. According to Elliot’s model , the absorption spectrum is modified because of the hole-electron interaction through a multiplicative function Fex given by 
where κ is the static dielectric constant, h is the Planck constant, is the reduced mass, and , where mh and ml are the heavy-hole and light-hole effective masses.
Including the hole-electron interaction, the total absorption spectrum αtot(E) can then be obtained as
where αlh, αhh and αso are the absorption spectra for transitions from the light-hole (LH), heavy-hole (HH) and spin-orbit split-off (SO) bands to the conduction band, respectively.
The effect of the incorporation of N in (In)GaNAs alloys can be described in different ways. We investigate here how the model chosen influences the calculated alloy absorption spectrum. We first present a simple model including isolated and pair N states using the BAC model. This model includes the nonparabolicity of the conduction and light-hole and split-off bands and the interaction between the split-off and conduction bands. In the second model, we then include the full LCINS distribution of localised states using the SCGF model.
Optical absorption of dilute nitride alloys in five-level BAC model
Here we first consider a simpler model, including isolated and pair N states and their interaction with the host semiconductor conduction, valence and spin-orbit split-off bands. The conventional BAC model treats the host III-V conduction band dispersion as a parabolic band . Test calculations that we have undertaken show that the inclusion of band non-parabolicity can strongly modify the calculated absorption spectra due to the change in the joint density of states caused by the non-parabolicity. In order to treat the host matrix conduction band nonparabolicity and the effect of N on the alloy conduction band dispersion, we construct a 5×5 Hamiltonian. This includes the Kane nonparabolicity of the host matrix conduction band, due to interactions with the light-hole and split-off bands, and treats the effects of N using a three-level BAC model , including isolated and pair N states. This Hamiltonian is given as
where Δso indicates the spin-orbit-splitting energy, and Ev0 is the energy of the valence band maximum. The energy of the isolated N levels (EN), N-N pair states (ENN) and the conduction band edge (Ec0) are assumed to vary with composition, x, and temperature, T, as 
The interaction parameters are assumed to vary as and , where the concentrations of single N and N-N pair states xN and xNN, respectively, are determined from the total N concentration x as xNN=6x2 and xN=x-2xNN. The chosen values of the above parameters for (In)GaNAs and GaNSb are given in Tables 1 and 2, respectively [16-18].
Calculating the eigenvalues of the matrix in Equation 9 gives the dispersion for five bands, namely, the light-hole and split-off valence bands and three conduction bands denoted by El, Em and Eu for lower, middle and upper bands, respectively. However, the five-band model of Equation 9 overestimates the LH nonparabolicity and omits the heavy-hole band dispersion. Figure 1 displays the band dispersion for In 0.04Ga0.96N0.01As0.99 where we have included a parabolic heavy-hole band, and a nonparabolic light-hole band calculated using the six-band Luttinger-Kohn (LK) valence band Hamiltonian . In this model, the band dispersion of the light-hole band is given by
Figure 1. The band dispersion for In 0.04Ga0.96N0.01As0.99 in five-level BAC model. Arrows show the possible optical transition from spin-orbit split-off band (solid arrows), light hole (dashed arrows) and heavy hole (dash-dotted arrows) to conduction sub-bands.
where and are the parabolic light-hole and spin-orbit bands, respectively, and
The fractional Γ character, defined as the contribution of the host matrix CB states to a given level, can be calculated by finding the eigenvectors of Equation 9 and is given by
Here, for each conduction sub-band, we use the appropriate energy E, as shown in Figure 1 by El, Em and Eu.
The joint density of states for transitions from the valence band vi to the conduction band cf, times the host matrix Γ character of the conduction state, determines the absorption strength. It is given by
where is the energy separation between the CB cf and the valence band vi states. The absorption spectrum for transitions between bands vi and cf is given, using Equation 4, by
Figure 1 represents the band dispersion and all possible optical transitions for In 0.04Ga0.96N0.01As0.99. The total absorption spectrum is given by the summation of nine individual absorption components, shown in Figure 1, for transitions from VBs vi (including LH, HH and SO bands) to CBs cf calculated using the 5-level BAC model (i.e.El, Em and Eu) as
where αSO-l, αSO-m and αSO-u are the absorption spectra from split-off band to lower, middle and upper sub-bands, respectively, and with a similar notation used for transitions from the HH and LH bands.
Optical absorption of dilute nitride alloys in the LCINS model
In order to include the full distribution of N states, we need to use the Green’s function method in the framework of the LCINS model. In our recent work , we showed that the conduction band Green’s function could be written as
where the (complex) energy shift of each localised state j is given by
We can solve Equations 19 and 20 self-consistently by an iterative method  to calculate ΔEj(E). The density of states per unit volume for the CB, projected onto the host matrix conduction states, is given by
The joint density of states between the CB host matrix components and valence band vi can then be obtained using
where Evi is the energy of the valence band vi which, as in the previous section, can be the LH, HH, or split-off (SO) band. We take into account the nonparabolicity of the LH band given by Equation 13, but assume parabolic heavy-hole and split-off bands, as we find that the parabolic split-off band dispersion in the relevant energy range is very close to that obtained when nonparabolicity effects are also included.
Having the joint DOS the optical absorption spectrum can be calculated using
The Green’s function given by Equation 19 ignores the nonparabolicity of the host semiconductor CB. In order to consider the Kane non-parabolicity, the Green’s function given by Equation 19 is modified to include the conduction-valence band interaction
Results and discussion
Here, we present the absorption spectra calculated using the five-level BAC and LCINS model and compare them with experiments. Perlin et al. [5,6,22-24] measured the optical absorption spectra for In0.04Ga0.96N0.01As0.99 and In0.08Ga0.92N0.015As0.985 and compared them with GaAs absorption data. Turcotte et al. [16,25] recently measured the optical absorption spectrum of GaNxAs1-x and InyGa1-yNxAs1-x for several values of x and y. Here, we calculate the absorption spectra for In0.04Ga0.96N0.01As0.99 and compare them with Skierbiszewski measurements at different temperatures.
Five-level model for InyGa1-yNxAs1-x
The interaction between the InGaAs valence and conduction bands and isolated and pair N states in InyGa1-yNxAs1-x can be described using Equation 9. The band structure parameters for In yGa1-yAs are taken to vary with In composition, y, and temperature, T, as shown in Tables 1 and 3. Also, the energy and the interaction of isolated and pair N states are taken to vary with In composition as given in Table 1. Figure 1 shows the calculated band structure of In 0.04Ga0.96N0.01As0.99 where the three conduction sub-bands (Eu(k), Em(k) and El(k)) are determined as the eigenvalues of Equation 9. Also, we consider the lowest eigenvalue of Equation 9 as the split-off band energy (ESO). The non-parabolic light-hole (ELH) is given by Equation 13, and the heavy-hole band (EHH) has been taken to be parabolic.
Table 3. Electrical and optical parameters in GaAs
The fractional Γ character of the conduction sub-bands is also required in order to calculate the joint density of states between the Γ-like conduction band components and the valence bands. Figure 2 shows the calculated Γ character of the CB for In 0.04Ga0.96N0.01As0.99, obtained using Equation 15. It is observed that in the lower sub-band, fΓ has its maximum value at the CBE and decreases toward zero at the top of the lowest band. It increases again from zero to a maximum value of around 0.4 and goes back to zero in the middle band. Then in the upper band, it increases gradually from its minimum at the bottom of the upper band, approaching an approximately constant value around E=2.1 eV.
Figure 2. The fractionalΓ character for InyGa1-yNxAs1-x withy=4% andx=1%. The fractional Γ character calculated using the five-level BAC model for conduction sub-bands, at T=10 K.
Figure 3 shows the calculated contributions of the different transitions to the total absorption spectrum. The solid, dashed and dotted lines in this figure represent the contributions due to transitions from the LH, HH and SO bands, respectively, to the conduction sub-bands. The red, blue and green lines indicate transition to the upper, middle and lower conduction sub-bands, respectively. The summation of these nine transitions is shown by the black dash-dotted line in this figure. Multiplying this by Fex gives the total absorption spectrum shown by the brown circles in this figure.
Figure 3. The absorption spectrum for InyGa1-yNxAs1-x, withy=4% andx=1% calculated using the five-level model atT=10 K. The contribution of the transitions between three valence bands and three conduction sub-bands are shown. The solid arrows designate the transitions from HH and LH bands to lower (l), middle (m) and upper (u) conduction sub-bands. The dashed arrows indicate the transitions from spin-orbit split-off band to conduction sub-bands.
The usual BAC model predicts a gap in the DOS [2,27] of (In)GaNAs alloys. However, it is clear from Figure 3 that the joint DOS for different transitions overlaps and fills this gap. Therefore, no gap is found in the absorption spectrum in (In)GaNAs alloys when using the BAC model.
Figure 4 compares the calculated absorption spectrum using the five-level BAC model with that measured and calculated by Skierbiszewski . The black dots in this figure are the experimental results for the absorption spectrum of In 0.04Ga0.96N0.01As0.99 at T=10 K. The solid black line shows the calculated absorption coefficient using the two-level BAC model  with constant VNc=2.7 eV and EN=1.65 eV. This line shows some discrepancies with the experimental data, especially around the transition to the upper sub-band of the BAC model. The arrows in this figure indicate the different transitions from the split-off, heavy- and light-hole band edges to the lower and upper sub-band edges, in the BAC model. The dashed red line in this figure displays the absorption spectrum calculated using the five-level BAC model, which shows much better agreement with the experimental measurements. This model still shows some steps corresponding to the transitions from the HH and LH band to the lower, middle and upper conduction sub-bands (see Figure 3), whereas the experiment shows a much smoother absorption spectrum and has only one pronounced step around E=1.85 eV. This is due to the fact that in the five-band model, we have considered isolated and pair N states. Considering the full distribution of N states makes the calculated absorption spectrum smoother, in better agreement with the experimental data.
Figure 4. The absorption spectrum for InyGa1-yNxAs1-x withy=4% andx=1%. The red dashed line and blue squares represent the absorption spectrum calculated using the five-level BAC model and the SCGF including the LCINS distribution of localised states (obtained for GaN 0.012As0.988), respectively. The black dots and line show the measured and calculated spectrum by Skierbiszewski  at temperature 10 K.
LCINS approach for InyGa1-yNxAs1-x
In order to calculate the absorption spectrum using the LCINS model, we first calculate the Green’s function for the CB, given by Equation 24. The inset in Figure 5 shows histograms of the distribution of localised states for GaN xAs1-x with x=0.84% and x=1.2%, calculated using the LCINS approach [3,4]. This figure shows that the LCINS distributions for x=0.84% and x=1.2% are very similar. This implies that the LCINS distribution for GaN xAs1-x with x=1.0% can be approximated by the one for x=1.2%. However, the calculated CBE at x=1.2% (indicated by E-) is about 30 meV lower than the CBE for x=0.84%. In the five-band BAC model, we assumed that the energy gap in In yGa1-yAs is given by Eg,GaAs-1.33y+0.27y2. Also, the BAC model parameters in Table 1 suggest that including 4% In in In 0.04Ga0.96N0.01As0.99 reduces the interaction VN(E) by about 5%, which we ignore in the LCINS model. Because we do not have the LCINS distribution for In 0.04Ga0.96N0.01As0.99, we approximate it here by the LCINS distribution calculated for GaN 0.012As0.988.
Figure 5. The density of state for InyGa1-yNxAs1-x, withy=4% andx=1%. Density of states calculated using the SCGF approach and LCINS distribution of N states at room temperature. Inset shows calculated distribution of N cluster state energies at low temperature, weighted by their interactions with the conduction band edge state for GaN xAs1-x with x=0.84% and x=1.2%.
We then solve Equations 20 and 24 self-consistently . Figure 5 shows the CB DOS calculated by Equation 21 for In 0.04Ga0.96N0.01As0.99 alloy. We observe that use of the LCINS distribution of states inhibits the gap predicted by the BAC model in the DOS of (In)GaNAs alloys.
The blue squares in Figure 4 show the calculated absorption spectrum at T=10 K including the full LCINS distribution of N states, which is compared with the absorption spectrum measured by Skierbiszewski . Clearly, the sharp steps that we saw in the five-level BAC model disappear due to the inclusion of the distribution of localised states. This gives a better overall agreement with the experimental data. The remaining discrepancies between the calculated and experimental data may be partly due to the fact that we have approximated the N distribution by the one that was previously calculated for GaN 0.012As0.988.
The room temperature absorption coefficient calculated from the SCGF method including the full LCINS distribution of N states is shown in Figure 6. The solid black line in this figure displays the absorption spectrum measured by Skierbiszewski  at T=300 K. The red circles in this figure show the absorption spectrum calculated in the LCINS model, where Fex is calculated using Equation 5. The blue diamonds here indicate the result when we consider Fex=1. The dashed blue and red lines in this figure show the optical absorption calculated using the five-level BAC model, assuming Fex=1 and given by Equation 5, respectively. Here, we again observe that the results calculated using the LCINS model have lower values in comparison with those calculated by the five-level BAC model. This is because of the differences in the band non-parabolicity that we have considered for the valence bands, in the LCINS and five-level LCINS models. Figures 4 and 6 suggest that Fex might have a stronger temperature dependence than what we have considered in our calculation. The temperature dependence is considered only in the static dielectric constant as shown in Table 3 in calculating exciton Rydberg energy in Equation 6.
Figure 6. Room temperature absorption spectrum for In 0.04Ga0.96N0.01As0.99 calculated using LCINS and five-level models. The absorption spectrum calculated using the LCINS (red circles and blue diamonds) and five-level BAC (dashed lines) models. The red and blue curves display the results with and without including the electron-hole interaction. The solid black line shows the experimental data.
The absorption spectrum for GaN xSb1-x
We can apply our method to calculate the band structure and absorption spectrum of other dilute nitride alloys. Here, we extend our calculations to investigate the absorption spectrum of GaNSb. The room temperature band gap of GaSb is about 725 meV, around half that of GaAs. Lindsay et al. [17,28] have reported that N-related defect levels lie close to the CBE in GaNSb and therefore strongly perturb the lowest conduction states in this alloy. The band gap and optical properties in GaN xSb1-x have been shown to be strongly affected and highly sensitive to the distribution of the nitrogen atoms. Lindsay et al.  found that there is a wide distribution of N levels lying close to and below the CBE. The higher-lying N states push the CBE down in energy, as in GaAs, but the large number of lower energy N states are calculated to mix in strongly with the conduction band edge states, severely disrupting the band edge dispersion in GaNSb.
Here, we first investigate the band structure and optical absorption spectra of GaN xSb1-x in the five-level BAC model and compare the results with the absorption spectra measured by Veal et al.  and Jefferson et al. . We then apply the SCGF method to GaN xSb1-x in the Section ‘LCINS model for GaN xSb1-x’. As the LCINS distributions have not yet been calculated for these alloys, we modify those calculated for GaNAs alloys and use them in our calculations.
Five-level BAC model for GaN xSb1-x
When a single Sb atom is replaced by N in GaSb, the N atom introduces a localised state with energy EN. However, a GaNSb alloy can also contain clusters of N atoms, such as N-N nearest neighbour pairs as well as larger clusters that introduce states in the band gap of GaSb. Table 2 contains the band parameters that we use for GaN xSb1-x, including the isolated N state energies, EN and N pair state energies, ENN relative to the valence band maximum energy and the BAC interaction parameters βN and βNN. As shown in this table, isolated N states are calculated to be less than 0.1 eV above the conduction band minimum, while the N pair states have energies that lie in the GaSb band gap. The calculated energy gap of GaN xSb1-x depends strongly on the assumed N distribution, reflecting that N cluster states introduce a series of defect levels close to the CBE in this alloy. In addition, the interaction parameters (βN and βNN) in GaN xSb1-x are calculated to be about 20% larger than for GaNAs alloys.
Figure 7 shows the conduction and valence band dispersion, calculated using the five-level BAC model given by the Hamiltonian of Equation 9. The solid lines in this figure show the conduction sub-bands. Here, we include the isolated and pair N states and their interaction with the GaSb conduction and valence bands, as explained in the previous section. Since EN is very close to the GaSb CBE and ENN just below it, we observe that the lower sub-band (El) is almost flat and located within the GaSb band gap. The band edge minimum for this band is 0.39 eV and its maximum energy is 0.45 eV. This implies that substitution of only x=1.2% N by Sb in GaSb reduced the energy gap rapidly from 725 to 390 meV. This value for the band gap of GaN xSb1-x with x=1.2% is very close to that which was previously measured  and calculated using k.p and ab initio pseudopotential  calculations.
Figure 7. The band dispersion of GaN 0.012Sb0.988 calculated by the five-level BAC model. The solid lines display the conduction bands including upper (Eu), middle (Em) and lower (El) sub-bands. The dashed lines show the spin-orbit split-off (ESO), light-hole (ELH) and heavy-hole (EHH) bands.
The middle sub-band (Em) lies between 0.55 and 0.78 eV, and the upper sub-band (Eu) minimum is close to 1.0 eV. The blue dashed line in Figure 7 shows the non-parabolic spin-orbit split-off band (ESO), calculated by the lowest eigenvalue of Equation 9. Also, the non-parabolicity of the light-hole band (ELH) has been taken into account using Equation 13, while we assumed that the heavy-hole (EHH) band has a parabolic dispersion.
Given the band dispersion, we can calculate the optical absorption as described earlier for InGaNAs. The red dashed line (with circles) in Figure 8 shows the absorption spectrum of GaN 0.012Sb0.988 calculated using the five-level BAC model. Green and brown solid lines in this figure show the absorption spectra measured by Mudd et al.  for GaN xSb1-x with x=1.18% and x=1.22%. Our calculated absorption edge is in good agreement with these experiments. However, there are two sharp steps in the calculated spectra corresponding to transitions from the light- and heavy-hole bands to the middle and upper sub-bands. We observe that the experimental absorption spectrum, α(E), increases from zero at the CBE, to about 3×103 cm -1 at energy E=0.55 eV, that is the calculated band edge of the middle sub-band. After this point the slope of the spectra decreases up to E=1 eV, the minimum of the upper conduction sub-band. Then, due to the transition from valence bands to the upper conduction sub-band, the magnitude of the calculated absorption spectra increases rapidly. The spin-orbit splitting energy, Δso, is 0.76 eV. Therefore, transitions from the split-off band commence at 1.15 eV, where we see a small increase in the calculated absorption spectrum due to transitions from the spin-orbit split off band to the lowest conduction sub-band El. The calculated optical absorption using a 5×5 k.p Hamiltonian, accounts well for the absorption edge. Wang et al.  have also measured the absorption edge of GaN xSb1-x with x=0.3%, x=0.7% and x=1.4%, with the measured band edge energies in very good agreement with those calculated by the five-level model of this paper.
Figure 8. The absorption spectrum of GaN 0.012Sb0.988. The absorption spectrum for Gan xSb1-x with x=1.2%, calculated using the SCGF model and including the distribution of localised states (blue diamonds), in comparison with experimental data measured for GaN xSb1-x with x=1.18% and x=1.22% (solid lines). The results calculated by the five-level BAC model are displayed by the red line with circles.
We note however that there are two sharp features in the calculated results that are not observed in experiments. This could be due to the fact that we included only isolated and pair N states in this model and ignored the distribution of N states and their inhomogeneous broadening. The results for the calculated absorption spectrum using a five-level BAC model suggest that we need to include the full distribution of N states in optical absorption calculations.
LCINS model for GaN xSb1-x
It has been shown that the calculated electronic structure of GaN xSb1-x strongly depends on the assumed distribution of N atoms . Therefore, in order to calculate an accurate band dispersion for this alloy, we need to have the distribution of localised states. Unfortunately, such a distribution has not been calculated for GaN xSb1-x. However, we expect that the distribution of N states in GaN xSb1-x should have a similar general form to the LCINS distribution that Lindsay et al.  have calculated for GaN xAs1-x and for GaN xP1-x alloys . Therefore, here, we consider the LCINS distribution of N states in GaN xAs1-x and, with some small modifications, use that for GaN xSb1-x alloys.
In the previous section, we have seen that the calculated energy of an isolated N state EN is at about 0.82 eV. So, we first need to shift the LCINS distribution of GaN 0.012As0.988 to locate the highest peak at this energy. The dashed red line in the inset of Figure 9 displays the LCINS distribution of N states, weighted by , calculated for GaN xAs1-x with x=1.2%, and shifted down in energy by 888 meV. This distribution can be approximated by three Gaussian distributions, each corresponding to different N environments. It is observed in this figure that if we align the main peak at E=0.82 eV, the lowest peak corresponding to pair N-N states is located at 0.55 eV, which is higher than the values that we considered for ENN in the BAC model. Therefore, we shift the Gaussian distributions corresponding to pairs and larger clusters of N states down by a further 70 meV. Moreover, the BAC model parameters in Table 2 suggest that in GaNSb, the interaction parameters, βN and βNN, are 20% stronger than in GaNAs. Therefore, we multiply the N LCINS values by 1.44 to account for this difference. The blue solid line in the inset of Figure 9 presents the distribution of N states that we consider for GaN 0.012Sb0.988 in our calculation.
Figure 9. The DOS of GaN 0.012Sb0.988, calculated using the SCGF method and including the LCINS distribution given in the inset. Inset displays the distribution of N states assumed for GaN 0.012Sb0.988 (blue line), in comparison with the LCINS distribution of GaN 0.012As0.988, shifted down in energy by 888 meV (dashed red line). The zero of energy is taken to be at the top of the GaSb valence band.
Having the distribution of N states, we are able to calculate the Green’s function for GaN 0.012Sb0.988 using Equations 19 and 20, self-consistently. Also, the density of CB states can be calculated using Equation 21. Figure 9 shows the DOS of GaN 0.012Sb0.988 calculated by the SCGF method and including the distribution of localised states shown by the solid blue line in the inset of Figure 9. The gaps corresponding to isolated and pair N states are clearly observed in this plot. Also, at energies around 0.65 eV, the DOS has a small gap that is related to the higher cluster of N states.
We can also calculate the absorption spectrum using the SCGF model. The blue line with diamonds in Figure 8 shows the calculated absorption coefficient using this method. As expected, this method shows a better agreement with experiments than the result of the five-level BAC model (shown by red circles in this plot).
For the considered N distribution, this calculation suggests more gaps in the DOS of GaN 0.012Sb0.988 compared to GaN 0.012As0.988. However, experimental data indicate that there are fewer features in the GaNSb absorption spectra than in the GaNAs ones. This could be due to inhomogeneities in the samples investigated experimentally, either due to fluctuations in the N composition in the experimental samples or because of intrinsic differences between the short-range N ordering in GaNSb and in InGaNAs samples.
Recent work by Mudd et al.  has shown that the composition dependence of the energy gap in GaN xSb1-x is well described using a three-level model including interactions between the host matrix band edge and the N isolated states and N-N pair states. The energy gap calculated using the LCINS model is also determined primarily by these interactions. The energy gap calculated here using the SCGF and LCINS method is consistent with experiment for the N composition x=1.2% which we consider and should closely follow the theoretical energy gap results presented in  as a function of N composition x.
In this paper, we presented an analysis of the optical absorption spectra of dilute nitride alloys, calculated using the band structure model presented in our earlier work . We have considered two different models to calculate the absorption spectra in InGaNAs and GaNSb alloys and compared our results with experimental measurements. We note however that there are some discrepancies between experimental data in similar samples that make quantitative comparison difficult.
Two models have been considered to calculate the absorption spectrum in these materials: a five-level BAC model and a LCINS-based model. The five-level BAC model included isolated and pair N states and their interactions with the host semiconductor valence and conduction bands. The results of this model for InGaNAs alloys give an overall good agreement with experiments, and predict accurate absorption edge for these alloys. However, the results of the five-level BAC model include several additional features not seen experimentally, supporting the need to consider a full distribution of N state energies in the electronic structure calculations.
We therefore extended our calculations to include the LCINS distribution using the SCGF approach presented in . The calculated absorption spectra using this approach for InGaNAs provide very good agreement with experiments, supporting the validity of the LCINS approach to describe dilute nitride conduction band structure.
Our calculated absorption spectra for GaNSb alloys fit well with experiments at the absorption edge , and predict the correct band gap in these alloys. However, the absorption spectrum calculated in the BAC model contains features associated with individual transitions to lower and upper sub-bands in the model that are not seen in the measured absorption spectra. Taking the distribution of localised states into account reduces the impact of these features and gives results more similar to experimental absorption. But we still see some dips in our calculated spectra that are not seen in any experiment. We conclude that the distribution of N states in the GaNSb alloys studied are different from that for InGaNAs samples. We conclude that further work is required to address and resolve why more structure is found in the calculated absorption spectra compared to what is observed in the experimentally measured spectra.
The authors declare that they have no competing interests.
SF proposed the SCGF approach to study the band structure of dilute nitride alloys. EOR suggested to apply this method to calculate the absorption spectrum. All calculations have been carried out by MS. All authors helped in drafting the manuscript. All authors read and approved the final manuscript.
This work was supported by the Science Foundation Ireland (06/IN.1/I90; 10/IN.1/I2994; 07/IN.1/I1810). The authors thank Tim Veal for providing measured value of absorption spectrum data for GaNSb samples.
J Phys: Condens Matter 2004, 16:3257-3276. Publisher Full Text
Semiconductor Sci Technol 2002, 17(8):803. Publisher Full Text
Appl Phys Lett 2000, 76:1279. Publisher Full Text
Phys Rev 1964, 133(2A):553-563. Publisher Full Text
Phys Rev B 1995, 52(3):1463. Publisher Full Text
Phys Rev B 1985, 32(2):1090. Publisher Full Text
Phys Rev 1957, 108(6):1384. Publisher Full Text
J Appl Phys 1982, 53(10):123-181. Publisher Full Text
Phys Rev Lett 1999, 82(6):1221-1224. Publisher Full Text
Phys E: Low-dimensional Syst Nanostructures 2006, 32:249-253. Publisher Full Text
Turcotte S, Beaudry J-N, Masut RA, Desjardins P, Bentoumi G, Leonelli R: Experimental investigation of the variation of the absorption coefficient with nitrogen content in GaAsN and GaInAsN grown on GaAs (001).
J Appl Phys 2008, 104(8):083511-083511. Publisher Full Text
Semiconductor Sci Technol 2009, 24:033001. Publisher Full Text
Solid State Commun 1999, 112:443-447. Publisher Full Text
J Appl Phys 2001, 89(11):5815. Publisher Full Text
Appl Phys Lett 2005, 87(13):132101-132101. Publisher Full Text
Appl Phys Lett 1998, 73:3703. Publisher Full Text
Appl Phys Lett 2000, 76(17):2409. Publisher Full Text
Turcotte S, Larouche S, Beaudry J-N, Martinu L, Masut RA, Desjardins P, Leonelli R: Evidence of valence band perturbations in GaAsN/GaAs (001): combined variable-angle spectroscopic ellipsometry and modulated photoreflectance investigation.
Phys Status Solidi B 2011, 248:1176-1179. Publisher Full Text
Appl Phys Lett 2006, 89(11):111921-111921. Publisher Full Text
Appl Phys Lett 2006, 88(15):152109-152109. Publisher Full Text
Appl Phys Lett 2013, 103:042110. Publisher Full Text
J Appl Phys 2009, 105(1):014904-014904. Publisher Full Text