Optical absorption of dilute nitride alloys using self-consistent Green’s function method

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.


Background
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 [1]. We have recently developed [2] 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: In y Ga 1-y N x As 1-x , for which the band structure has been widely studied and many of the features are well established and GaN x Sb 1−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 [2].
For In y Ga 1-y N x As 1-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 In y Ga 1-y N x As 1-x using the SCGF approach [2] 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 x Sb 1−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 http://www.nanoscalereslett.com/content/9/1/51 of GaNSb alloys, supporting the presence of a band anticrossing 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.

Methods
The absorption spectrum α(E) describes the rate of absorption of photons with energy E = ω 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][8][9][10] where ω is the photon energy, e and m 0 are the electron charge and mass, c is the speed of light and n r is the refractive index. f is the Fermi-Dirac distribution function. Here, we assume a filled valence and empty conduction band, so that f v − f c = 1. The matrix element in Equation 1 can be written for transitions between valence p-and conduction s-like zone centre states as [9] where E p = (2m 0 / 2 )p 2 is the interaction energy, and the momentum interband matrix element, p, can be estimated from experiment as where m * e is the conduction band effective mass, E g 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 = E ck − E vk is the transition energy between conduction (E ck ) and valence (E vk ) states with wavevector k, J cv (E) is the joint density of states, and we ignore for now the energy dependence of M b . The hole-electron interaction can be included assuming Elliot's theory, which applies to parabolic and nondegenerate bands. According to Elliot's model [11], the absorption spectrum is modified because of the hole-electron interaction through a multiplicative function F ex given by [8] where and R y is the exciton Rydberg energy, given by [12,13] where κ is the static dielectric constant, h is the Planck 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 http://www.nanoscalereslett.com/content/9/1/51 semiconductor conduction, valence and spin-orbit splitoff bands. The conventional BAC model treats the host III-V conduction band dispersion as a parabolic band [14]. 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 threelevel BAC model [15], including isolated and pair N states. This Hamiltonian is given as where so indicates the spin-orbit-splitting energy, and E v0 is the energy of the valence band maximum. The energy of the isolated N levels (E N ), N-N pair states (E NN ) and the conduction band edge (E c0 ) are assumed to vary with composition, x, and temperature, T, as [15] and The interaction parameters are assumed to vary as  Tables 1 and 2, respectively [16][17][18].
Calculating the eigenvalues of the matrix in Equation  9 gives the dispersion for five bands, namely, the lighthole and split-off valence bands and three conduction bands denoted by E l , E m and E u 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.04 Ga 0.96 N 0.01 As 0.99 where we have included a parabolic heavy-hole band, E HH = E v0 − 2 k 2 /(2m h ) and a nonparabolic light-hole band calculated using the six-band Luttinger-Kohn (LK) valence band Hamiltonian [21]. In this model, the band dispersion of the light-hole band is given by where E l = E v0 − 2 k 2 /(2m l ) and E so = E v0 − so − 2 k 2 /(2m so ) 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 cal- N-N interaction [17] β NN 3.39 (eV) http://www.nanoscalereslett.com/content/9/1/51 culated 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 E l , E m and E u .
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 E f ,i = ω = E vi − E cf 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.04 Ga 0.96 N 0.01 As 0.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. E l , E m and E u ) 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 [2], 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 [2] to calculate E j (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 E vi 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 heavyhole 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][23][24]  The interaction between the InGaAs valence and conduction bands and isolated and pair N states in In y Ga 1-y N x As 1-x can be described using Equation 9. The band structure parameters for In y Ga 1−y As 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.04 Ga 0.96 N 0.01 As 0.99 where the three conduction sub-bands (E u (k), E m (k) and E l (k)) are determined as the eigenvalues of Equation 9. Also, we consider the lowest eigenvalue of Equation 9 as the split-off band energy (E SO ). The non-parabolic light-hole (E LH ) is given by Equation 13, and the heavy-hole band (E HH ) has been taken to be parabolic. 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.04 Ga 0.96 N 0.01 As 0.99 , obtained using Equation 15. It is observed that in the lower subband, 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.   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 [5]. The black dots in this figure are the experimental results for the absorption spectrum of In 0.04 Ga 0.96 N 0.01 As 0.99 at T = 10 K. The solid black line shows the calculated absorption coefficient using the two-level BAC model [6] with constant V Nc = 2.7 eV and E N = 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 subbands (see Figure 3), whereas the experiment shows a  The absorption spectrum for In y Ga 1-y N x As 1-x with y = 4% and x = 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.012 As 0.988 ), respectively. The black dots and line show the measured and calculated spectrum by Skierbiszewski [5] at temperature 10 K. http://www.nanoscalereslett.com/content/9/1/51 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.

LCINS approach for In y Ga 1-y N x As 1-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 x As 1−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 x As 1−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 y Ga 1−y As is given by E g,GaAs − 1.33y + 0.27y 2 . Also, the BAC model parameters in Table 1  We then solve Equations 20 and 24 self-consistently [2]. Figure 5 shows the CB DOS calculated by Equation 21 for In 0.04 Ga 0.96 N 0.01 As 0.99 alloy. We observe that use of the  Figure 5 The density of state for In y Ga 1-y N x As 1-x , with y = 4% and x = 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 x As 1−x with x = 0.84% and x = 1.2%.

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 [5]. 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.012 As 0.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 [5] at T = 300 K. The red circles in this figure show the absorption spectrum calculated in the LCINS model, where F ex is calculated using Equation 5. The blue diamonds here indicate the result when we consider F ex = 1. The dashed blue and red lines in this figure show the optical absorption calculated using the five-level BAC model, assuming F ex = 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 F ex 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.

The absorption spectrum for GaN x Sb 1−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 x Sb 1−x have been shown to be strongly affected and highly sensitive to the distribution of the nitrogen atoms. Lindsay et al. [28] 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 x Sb 1−x in the five-level BAC model and compare the results with the absorption spectra measured by Veal et al. [20] and Jefferson et al. [29]. We then apply the SCGF method to GaN x Sb 1−x in the Section 'LCINS model for GaN x Sb 1−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 x Sb 1−x
When a single Sb atom is replaced by N in GaSb, the N atom introduces a localised state with energy E N . 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 x Sb 1−x , including the isolated N state energies, E N and N pair state energies, E NN 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 x Sb 1−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 x Sb 1−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 E N is very close to the GaSb CBE and E NN just below it, we observe that the lower subband (E l ) 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 x Sb 1−x with x = 1.2% is very close to that which was previously measured [20] and calculated using k.p [28] and ab initio pseudopotential [30] calculations.
The middle sub-band (E m ) lies between 0.55 and 0.78 eV, and the upper sub-band (E u ) minimum is close to 1.0 eV. The blue dashed line in Figure 7 shows the non-parabolic spin-orbit split-off band (E SO ), calculated by the lowest eigenvalue of Equation 9. Also, the non-parabolicity of the light-hole band (E LH ) has been taken into account using Equation 13, while we assumed that the heavy-hole (E HH ) 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.012 Sb 0.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. [31] for GaN x Sb 1−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 × 10 3 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 E l . The calculated optical absorption using a 5 × 5 k.p Hamiltonian, accounts well for the absorption edge. Wang et al. [32] have also measured the absorption edge of GaN x Sb 1−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. 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 x Sb 1−x
It has been shown that the calculated electronic structure of GaN x Sb 1−x strongly depends on the assumed distribution of N atoms [28]. 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 x Sb 1−x . However, we expect that the distribution of N states in GaN x Sb 1−x should have a similar general form to the LCINS distribution that Lindsay et al. [3] have calculated for GaN x As 1−x and for GaN x P 1−x alloys [17]. Therefore, here, we consider the LCINS distribution of N states in GaN x As 1−x and, with some small modifications, use that for GaN x Sb 1−x alloys.
In the previous section, we have seen that the calculated energy of an isolated N state E N is at about 0.82 eV. So, we first need to shift the LCINS distribution of GaN 0.012 As 0.988 to locate the highest peak at this energy. The dashed red line in the inset of Figure 9  √ N c , calculated for GaN x As 1−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 E NN 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.012 Sb 0.988 in our calculation.
Having the distribution of N states, we are able to calculate the Green's function for GaN 0.012 Sb 0.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.012 Sb 0.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.012 Sb 0.988 compared to GaN 0.012 As 0.988 [2]. 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. [31] has shown that the composition dependence of the energy gap in GaN x Sb 1−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 [31] as a function of N composition x.

Conclusions
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 [2]. 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 [2]. 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 [31], 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. http://www.nanoscalereslett.com/content/9/1/51