Abstract
Motivated by the research and analysis of new materials for photovoltaics and by the possibility of tailoring their optical properties for improved solar energy conversion, we have focused our attention on the (GaAs)_{1−x}Ge_{2x} series of alloys. We have investigated the structural properties of some (GaAs)_{1−x}Ge_{2x} compounds within the localdensity approximation to densityfunctional theory, and their optical properties within the Quasiparticle Selfconsistent GW approximation. The QSGW results confirm the experimental evidence of asymmetric bandgap bowing. It is explained in terms of violations of the octet rule, as well as in terms of the order–disorder phase transition.
Keywords:
Photovoltaics; III–V IVdoped alloys; Bandgap bowing; Order–disorder phase transition; DFT; Quasiparticle Selfconsistent GWIntroduction
The design of semiconductors with controlled bandgaps E_{G}, unit cell parameters, and low defect concentration is the ultimate aim in several important areas of industrial applications—electronics, photonics, lightemitting devices, and photovoltaics (PV). Efficient collection of solar energy requires materials that absorb light from different portions of the solar spectrum, followed by efficient conversion into electrons and holes at p–n junctions. A natural approach to the design of new semiconductors is to alloy two materials with similar lattice parameters but different bandgaps. For example, Ge (E_{G} = 0.67 eV [1] at 300 K) and GaAs (E_{G} = 1.43 eV [1] at 300 K) have very similar lattice parameters, 5.649 and 5.66 Å, respectively [2,3]. There is thus the appealing possibility that (GaAs) _{1−x}Ge_{2x} alloys with intermediate bandgaps can be realized, in particular one characterized by a direct gap, 1 <E_{G} < 1.4 eV (i.e., the average between the bandgaps of Ge and GaAs), which corresponds to the maximum efficiency solar cell for a single bandgap material [4]. Indeed, several theoretical [512] and experimental [1318] papers have been published on studies of metastable alloys between III–V and IV group semiconductors, formally (III–V) _{1−x}(IV)_{2x} compounds.
A group of such mixed single crystal metastable semiconductors covering a wide composition range was synthesized by vapor phase deposition techniques. Noreika et al. [17] deposited (GaAs)_{1−x}Si_{x} on GaAs(111) by the reactive rf sputtering technique, and they reported an optical bandgap at room temperature of about 1.28 eV for (GaAs)_{0.45}Si_{0.55}. Baker et al. [18] measured the Raman spectra of (GaSb)_{1−x}Ge_{x} alloys, and found both GaSb and Gelike optical modes. The Gelike mode frequency depends on the alloy’s composition within about 40 cm^{−1}, whereas the GaSblike mode does not.
(III–V)–IV alloys such as (GaAs)_{1−x}Ge_{2x}, characteristically display a large negative, Vshaped bowing of E_{G} as a function of the alloy composition x. A minimum value of ~0.5 eV was detected by Barnett et al. [13] at a Ge concentration of about 35%, corresponding to the critical value (x_{c}) for phase transition between an ordered zincblende (ZB) and a disordered diamond structure. In the ordered GaAsrich phase, Ga and As preferentially form donor–acceptor pairs, whereas in the Gerich phase, they are randomly distributed in the alloy forming a mixture of n type (As in Ge) and p type (Ga in Ge) semiconductors. This phase transition has been put forward to explain [19] the large bowing. Several models have been developed for the ZB → diamond phase transition [8,1926]. The stochastic model by Kim and Stern [22] well reproduces this phase transition along the <100> direction at x_{c} = 0.3. However, it poorly describes the growth along the <111> direction, with accumulation of Ge on alternate {111} planes. In general, kinetic models seem to be more appropriate descriptors of the ZB—diamond phase transition than thermodynamic ones: the latter do not take into account the nonuniqueness of the critical composition x_{c} as a function of kinetic growth; they require as input the critical concentration at which the transition takes place, but no restrictions on the formation of Ga–Ga and As–As bonds are imposed. Other models based on the percolation method [26] lead to ZB → diamond transition at x_{c} = 0.57; percolation theory also does not account for different growth conditions. Rodriguez et al. [8] reported that the growth direction and avoidance of “bad bond” formation (i.e., Ga–Ga and As–As bonds) (longrange order, LRO) effects are the main factors responsible for atomic ordering in (GaAs)_{1−x}Ge_{2x} alloys. According to the same model, E_{G} is influenced only by nearest neighbor (NN) atomic interactions (shortrange order, SRO) effects. In an extension of the stochastic model of growth along the <100> direction, Holloway and Davis [23,24] formulated a model for alloys grown in the <100> and <111> directions. SRO effects are common for these two directions. In contrast, the impact of LRO is quite different: a tendency to convert to <111> As growth is predicted [24] as a consequence of the instability of the growth in the <111> Ga direction. In a previous paper [25], the same authors note that the transition from ZB to diamond does not affect the energy gap of (GaAs)_{1−x}Ge_{2x}: this model predicts a critical concentration for the order–disorder transition with x_{c} = 0.75, without any dependence on the method of growth. SRO and LRO effects on the electronic properties of many other IVdoped III–V alloys have also been compared by combining the special quasirandom structures (SQS) and the simulatedannealing (SA) methods for cells of various sizes in conjunction with an empirical pseudopotential approach [27]. In particular, the direct bandgaps of ideal random Al_{1−x}Ga_{x} As, Ga_{1−x} In_{x}P, and Al_{1−x} In_{x} As alloys were studied. SRO effects are reported to increase the optical bowing of the direct bandgap.
Surface faceting has also been detected in these systems, reported to take place with a subsequent phase separation between the GaAsrich ZB and Gerich diamond region during the growth on (001)oriented GaAs substrates [15]. A direct consequence attributable to the faceting is the bandgap narrowing of such (GaAs)_{1−x}Ge_{2x} (0 <x < 0.22) alloy layers grown by lowpressure metal–organic vapor phase epitaxy. A similar phenomenon has been reported only once previously, for InAs_{y}Sb_{1−y} grown by molecular beam epitaxy grown at low temperature (T_{g}) [28]. It has been also demonstrated that growth temperature [16] strongly affects the nature of the alloy. (GaAs)_{1−x}Ge_{2x} layers have been epitaxially grown on GaAs (100) substrates at different temperatures. Transmission electron microscopy analysis revealed that at T_{g} = 550°C, Ge separated from GaAs into domains of ~100 Å. Singlephase alloys are detected differently at T = 430°C.
In spite of considerable recent research in novel complex materials for photovoltaics, the relationship between chemical and optical properties of III–V–IV alloys and similar materials is still unknown, and is a matter of current debate. In the present paper, we investigate the chemical nature of the bowing in (GaAs)_{1−x}Ge_{2x} alloys. In particular, we theoretically investigated the structural and optical properties of four different intermediate structured compounds that range between “pure” GaAs and “pure” Ge (x_{Ge} = 0.25, 0.50 (two samples), 0.75).
Computational Details
We performed calculations by using densityfunctional theory (DFT), within both the localdensity approximation (LDA) [29,30] and the generalized gradient approximation (GGA) of Perdew and Wang [3133]. We used Blöchl’s allelectron projectoraugmented wave (PAW) method [34,35], with PAW potentials with d electrons in the semicore for both Ga and Ge. Cutoff energies of 287 and 581 eV were set as the expansion and augmentation charge of the plane wave basis. The force convergence criterion for these models was 0.01 eV/Å. The initial (GaAs)_{1−x}Ge_{2x} models consisting of eight atoms were optimized with a 10^{3} Γ centered kpoints sampling scheme.
All the total energy calculations were also performed with the generalized fullpotential LMTO method of Ref. [36]. Calculated structural properties and heats of reaction predicted by the two methods were almost identical, indicating that the results are well converged.
The thermodynamic stability of these alloys was calculated as the ΔE products–reactants of the equation:
It is expected that LDA and GGA predict reasonable heats of reaction of the type in Eq. 1, since reactions involve rearrangement of atoms on a fixed (zincblende) lattice, and there is a large cancelation of errors. Optical properties are much less well described. The LDA is well known to underestimate semiconductor bandgaps, and moreover, the dispersion in the conduction band is poor. In Ge, the LDA gap is negative and Γ_{1c} is lower than L_{1c} in contradiction to experiment. Also the ΓX dispersion is often strongly at variance with experiment: in GaAsX_{1c−}Γ_{1c} is about twice the experimental value of 0.48 eV.
When considering (GaAs)_{1−x}Ge_{2x} alloys, any of the three points (Γ, X, L) may turn out to be minimumgap points, so all must be accurately described. Thus, the LDA is not a suitable vehicle for predicting optical properties of these structures.
It is widely recognized that the GW approximation of Hedin [37] is a much better predictor of semiconductor optical properties. The GW approximation is a perturbation theory around some noninteracting Hamiltonian H_{0}; thus the quality of the GW result depends on the quality of H_{0}. It is also important to mention that for reliable results, care must be taken to use an allelectron method [38]. We adopt here a particularly reliable allelectron method, where not only the eigenfunctions are expanded in an augmented wave scheme, but the screened coulomb interaction W and the selfenergy Σ = iGW are represented in a mixed plane wave and molecular orbital basis [39,40]. All core states are treated at the Hartree–Fock level.
GW calculations in the literature usually take H_{0} from the LDA; thus, we may call this the G^{LDA}W^{LDA} approximation. There are many limitations to G^{LDA}W^{LDA}; see e.g., Ref. [41]. In particular, the G^{LDA}W^{LDA} gap for GaAs is 1.33 eV. The Quasiparticle Selfconsistent GW (QSGW) approximation, recently developed by one of us [42], overcomes most of these limitations. Semiconductor energy band structures are well described with uniform reliability. Discrepancies with experimental semiconductor bandgaps are small and highly systematic (e.g., for most semiconductors [41]), and the origin of the error can be explained in terms of ladder diagrams missing in the random phase approximation (RPA) to the polarizability Π(r, r′,ω) [43]. While standard QSGW would be sufficient for this work, we can do a little better by exploiting our knowledge of the small errors originating from the missing vertex in Π. In principle ladder diagrams can be included explicitly via the Bethe–Salpeter equation, but it is very challenging to do. It has never been done in the QSGW context except in a very approximate manner [43]. On the other hand, in sp semiconductors, the consequences of this vertex are well understood. The RPA results in a systematic tendency for the dielectric constant, ε_{∞}, to be underestimated. The error is very systematic: to a very good approximation ε_{∞} is too small by a universal factor of 0.8, for a wide range of semiconductors and insulators [44]. This fact, and the fact that quasiparticle excitations are predominantly controlled by the static limit of W, provides a simple and approximate remedy to correct this error: we scale Σ (more precisely ) by 0.8. While such a postprocessing procedure is admittedly ad hoc, the basis for it is well understood and the scaling results in a very accurate ab initio scheme for determining energy band structures (to within ~0.1 eV when the effect of zeropoint motion on bandgaps is taken into account) and effective masses for essentially any semiconductor. Here, we adopt this scaling procedure to refine our results to this precision. In any case corrections are small, and our conclusions do not depend in any way on this scaling. Results for GaAs and Ge are shown in Table 1.
Table 1. Left, LDA calculated bandgaps (LMTO [36], Spin–Orbit effects included) for Γ, XL points for GaAs and Ge Right, QSGW bandgaps for the same points in Ge and GaAs (eV, 0 K), compared with measured values at 0 K. The selfenergy was scaled by a factor 0.8, as described in the text. Raw (unscaled) QSGW levels are slightly larger than experiment
Results
We performed preliminary calculations at the DFT level of GaAs and Ge; Table 2 lists the main structural optimized parameters of the two most stable polymorphs of GaAs, zincblende (ZB, group 216, F43m, Z = 4) and wurtzite (WZ, group 186, P 63 mc, Z = 2) and of Ge in its cubic form (group 227, Fd3m, Z = 8).
Table 2. The energy difference (ΔE, per unit, eV) between ZB and WZ polymorphs of GaAs, lattice constant, a, and bulk moduli B(GPa) of GaAs (ZB) and Ge (diamond)
As seen from Table 2, the LDA generates structural properties closer to experiment than GGA in this context. Thus, we use LDA to study structural properties. Both Ga–As and Ge–Ge bond lengths are 2.43 Å in their most stable polymorph.
ZB–GaAs is constituted by interpenetrating fcc sublattices of cations (Ga) and anions (As). The diamond lattice of Ge may be thought of as the ZB structure with Ge occupying both cation and anion sites. Here, we consider 8atom (GaAs)_{1−x}Ge_{2x} compounds that vary the Ge composition, including pure GaAs (x = 0) to x= 0.25 (2 Ge atoms), x = 0.50 (4 Ge atoms), x = 0.75 (6 Ge atoms) (see Fig. 1), and finally pure Ge (x = 1).
Figure 1. Four (GaAs)_{1−x}Ge_{2x} models investigated. [Ga, small gray; As, large white; Ge, large black]
At first, we performed an analysis of the Ge dimer in bulk GaAs, at site positions (1/4, 1/4,1/2) and (0,1/2,3/4). We denote this as “alloy model I”, the dimer in an 8atom GaAs cell with lattice vectors (100), (010), and (001) before relaxation. It can be considered a highly concentrated molecular substitutional Ge_{2} defect in GaAs, for which we predict stability owing to the donor–acceptor selfpassivation mechanism.^{1} The first layer of I consists only of As; the second and the third layers (along [001]), are Ge–Ga, and Ge–As, respectively. The fourth is pure Ga. Then, the overall sequence is a repeated “sandwichlike” structure, ···/As/Ge–Ga/Ge–As/Ga/···. The bond lengths were calculated to be 2.38 (Ga–Ge), 2.42 (Ge–Ge), 2.44 (Ga–As), and 2.47 Å (Ge–As)—relatively small variations around the calculated values in bulk Ge and GaAs (2.43 Å). This is perhaps not surprising as the electronic structure can roughly be described in terms of nearly covalent twocenter bonds [electronegativity χ = 1.81, 2.01, and 2.18, for Ga, Ge, and As, respectively (http://www.webelements.com webcite)]. In the alloy I the number of “bad bonds” [7,8], i.e., the number of III–IV and IV–V nearest neighbors, is 12, or 37.5% of the total. According to the Bader analysis [5759], in the pure host, the difference in electronegativity is responsible for charge transfer from cation to anion. In the alloy formation process, the introduction of Ge reduces the ionic character of the GaAs bond, while increasing the ionic character of the Ge–Ge bond. When a Ge dimer is inserted in GaAs, 0.32 electrons are transferred away from Ge_{Ga} site, while Ge_{As} gains 0.21 electrons The charge deficit on Ga, is reduced from 0.6 electrons in bulk GaAs to 0.47e, while the charge excess on As is reduced from 0.6 to 0.5e. ΔE for reaction Eq. (1) was 0.55 eV, and the optimized lattice parameter was a = 5.621 Å. We have also considered Ge donors (Ge_{Ga}) and acceptors (Ge_{As}) in the pure 8atom GaAs host cell, separately. The formation energy has been computed according to the Zhang–Northrup formalism [60]. In particular, we calculate ΔE to be 1.03 eV for Ge_{Ga} and 0.84 for Ge_{As}. The sum of the single contributions (1.87 eV) is larger than the heat of formation of the dimer, structure I (0.55 eV). Two reasons explain this difference in energy. First, in the I model alloy, at least one correct bond III–V is formed while in the separate Ge_{Ga} (IV–V) and Ge_{As} (IV–III) cases only bad bonds are formed. The isolated Ge_{Ga} is a donor; the isolated Ge_{As} is an acceptor. Neither is stable in their neutral charged state. We have tested it in a previous analysis (see Footnote) where we calculated +1 and −1 as the most stable charged state for Ge_{Ga} and Ge_{As}, for almost the range of the electronic chemical potential, μ_{e}. These two charged states are indeed formally isoelectronic with the host GaAs. That the stabilization energy 1.87–0.55 = 1.33 eV is only slightly smaller than the host GaAs bandgap establishes that the pair is stabilized by a selfpassivating donor–acceptor mechanism.
We considered two alternative structures for the x_{Ge} = 0.50 case. In the IIa structure, Ge atoms are substituted for host atoms at (1/2, 0, 1/2), (1/2, 1/2, 0), (3/4, 3/4, 1/4), and (3/4, 1/4, 3/4); then the lattice was relaxed. It results in a stacking ···/Ga–Ge/Ge–As/··· along <001>. The three cubic directions are no longer symmetryequivalent: the optimized lattice parameters were found to be a = 5.590 Å, b = c = 5.643 Å. The four intralayer bond lengths were calculated to be Ga–Ge (2.39 Å), Ge–As (2.48 Å), Ge–Ge (2.42 Å), and Ga–As (2.44). Because of the increased amount of Ge, structure IIa was less polarized than I, as confirmed by the slightly more uniform bond lengths. In IIa alloy, the number of “bad bonds” is 16 (i.e., 50%) and ΔE rises to 0.72 eV. In the IIb structure, Ge atoms are substituted for host atoms at (1/4, 1/4, 1/4), (1/4, 3/4, 3/4), (3/4, 3/4, 1/4), and (3/4, 1/4, 3/4). This structure consists of a stack of pure atomic layers, ···/Ga/Ge/As/Ge···, and thus it contains only nearest neighbors of the (Ga–Ge) and (Ge–As) type: thus all bonds are “bad bonds” in this IIb compound. Bond lengths were calculated to be 2.40 Å and 2.49 Å, respectively, and optimized lattice parameters were a = c = 5.682, b = 5.560 Å. In this structure, ΔE = 1.40 eV, almost double that of IIa with identical composition. It supports the picture [7,8] that III–IV and IV–V bonds are less stable than their III–V, IV–IV counterparts.
This result confirmed findings of an analysis of the substitutional defect Ge in GaAs (see Footnote). In that case, we checked the stability of Ge_{2} dimers (donor–acceptor pair formation) versus isolated Ge couples (ntype Ge + ptype Ge) in GaAs matrix in GaAs supercells. We calculate the energy reaction Ge_{2}:GaAs → Ge_{Ga}:GaAs + Ge_{As}:GaAs to be positive, with ΔE = 0.39 eV, and interpret this as the gain of one III–IV (Ga–Ge) and one IV–V (Ge–As) bond and the loss of one IV–IV (Ge–Ge) bond. (Note that the IIb structure corresponds to the high concentration limit of isolated couples.) According to phase transition theory, the symmetry lowering for the two intermediate systems is the fingerprint of an ordered–disordered phase transition [61]. The calculated deviation from the ideal cubic case (c/a = 1) is 0.94 and 2.15% for IIa and IIb models, respectively, confirming energetic instability for the IIb alloy.
The last model, III, [Ge] = 0.75, consists of pure Ge except that Ga is substituted at (0, 0, 0) and As at (1/4,1/4,1/4). The calculated bond lengths were 2.39, 2.43, 2.45, and 2.48 Å for Ga–Ge, Ge–Ge, Ga–As, and Ge–As, respectively. Cubic symmetry is restored: the optimized lattice parameter (a = 5.624 Å) is nearly identical to structure I. Similarly, ΔE is almost the same as I (~0.54 eV). Indeed, I and III are formally the same model with the same concentration (25%) of Ge in GaAs (I) and GaAs in Ge (III) and the same number of bad bonds, 12. By analogy to model I, we have also calculated the formation energy of a single substituted Ge in the cell. We have also made a preliminary calculation of the stability of isolated Ga acceptors (Ga_{Ge}) and As donors (As_{Ge}) versus that of the substitutional molecular GaAs_{Ge2} in Ge pure host (a supercell of 64 atoms see Footnote); for such concentrations (x_{GaAs} = 0.0312 = 1/32 and x_{Ga(As)} = 0.0156 = 1/64), the molecular substitutional GaAs_{Ge2} is only stabilized by 0.057 eV with respect to the separate couple acceptor–donor. This small stabilization for GaAs_{Ge2} compared to isolated Ga_{Ge} and As_{Ge} confirms the expected similar probability of finding a mixture of n type and p type semiconductors in the “disordered” Gerich phase. For reference states needed to balance a reaction, we used the most stable polymorph, of the elemental compounds i.e., orthorhombic Ga and rhombohedral As [62]. Garich and Asrich conditions have been considered, corresponding to μ_{Ga(As)} = μ_{Ga(As)}^{bulk}, respectively. In the case of the 8atom cells, the formation energy for Ga_{Ge} and As_{Ge} are 0.26 eV and 0.58 eV, respectively. Thus, the model III alloy stabilizes the isolated III and V substitutionals by 0.30 eV (ΔE(Ga_{Ge}) + ΔE(As_{Ge})−ΔE_{III}). This energy is 0.4 eV less than the host bandgap, indicating that the stabilization energy is a little more complicated than a simple selfpassivating donor–acceptor mechanism, as we found for the Ge molecule in GaAs.
Collecting ΔE for the different systems containing equal numbers of Ge cations and anions, we find an almost exactly linear relationship between ΔE and the number of bad bonds, as Fig. 2 shows. This striking result confirms that the electronic structure of these compounds is largely described in terms of independent twocenter bonds. For stoichiometric compounds, it suggests an elementary model Hamiltonian for the energetics of any alloy with equal numbers of Ge cations and anions: ΔE = 0.54, 0.72, 1.40 eV for N = 12, 16, 32, respectively, where N is the number of bad bonds.
Figure 2. Heat of formation (ΔE) of the alloy models versus the number of “bad bonds”
Even if small variations are expected in the lattice parameter of the alloys, the Vegard’s law:
where a_{(GaAs)1 −xGe2x}, a_{Ge}, a_{GaAs} are the lattice parameters of the final alloy and its components, respectively, represents a useful tool for predicting a trend in terms of lattice parameter variation for our alloy models. We have thus tested the predicted versus calculated values for the lattice parameter. In particular, for models I and III, we have used the calculated lattice parameters, a (LDA). On the other hand, since for the reducedsymmetry models IIa and IIb is a ≠ b = c and a = c ≠ b, respectively, we have approximated the lattice parameter as the cubic root of the volume of each of the x_{Ge} = 0.5 cell models. Table 3 reports the experimental, theoretically predicted, and calculated (LDA) lattice parameters, based on Vegard’s Law. These values showed the almost perfect matching between GaAs and Ge lattice parameter and at the same time the marked deviation of model IIb from the trend, thus confirming the main contribution of “bad bonds” to the final instability of the alloy.
Table 3. Lattice parameters: a_{exp} obtained by Eq. (2) using experimental lattice parameters; a_{theor} calculated from Eq. (2), but with optimized lattice parameters at the PAW/LDA level for Ge and GaAs; a_{calc} the PAW/LDA optimized lattice parameters for models I, IIa, IIb, and III. (Italic is for values extrapolated as^{3}√V.)
We have performed QSGW calculations on the optimized structures (I, IIa, IIb, and III) and also for the pure GaAs and Ge 8atom cells. Table 4 shows the QSGW bandgaps for the Γ and R points. In Fig. 3, we report the electronic structure for all the models considered (pure GaAs, Ge, and the intermediate alloy models). In the simple cubic supercell, R corresponds to L of the original ZB lattice; Γ has both X and Γ points folded in. It is evident that there is a pronounced bowing at both Γ and L, as also shown in Fig. 4, where we report the QSGW bandgap as function of Ge concentration.
Table 4. QSGW bandgaps for Γ and R in ordered (GaAs)_{1−x}Ge_{2x} alloys
Figure 3. Electronic structure for the considered systems, GaAs (first, left up), Ge (last, right bottom), and the four intermediate alloys I, IIa, IIb, III
Figure 4. QSGW calculated bowing of the bandgap at Γ and R versus different concentration of Ge atoms
Figure 5 summarizes the relationship between DFT (a, lattice constant) and QSGW (E_{G}, bandgap). From left to right, we report the bandgap for Ge, III, IIb, IIa, I, and pure GaAs. Once more we note the marked discontinuity for IIb from the general trend.
We finally remark on the lack of definitive analysis of the state of the alloy. In particular, the validity of two viewpoints, a probabilistic growth model based on a layerbylayer deposition that rejects highenergy bond formation (As–As, Ga–Ga) [22] and a thermodynamic equilibrium based on an effective Hamiltonian (but which is able to describe electronic states [63]) needs to be assessed. In thermodynamic models [19], the Boltzmann weight can always be realized regardless of the measurement time; in probabilistic models [21,22] further equilibration after the atomic deposition is not possible. It is apparent from the results of Figs 4 and 5 that our theoretical alloy models, and optical properties in particular, depend sensitively on the arrangement of atoms in the alloy.
Conclusions
We have focused on the class of (III–V)_{1−x}IV_{2x} alloys, as candidate new materials with applications relevant to photovoltaics. Previous experiments reported an asymmetric (nonparabolic) bowing of the bandgap as a function of the concentration of the III–V and IV constituents in the alloy.
We have built and optimized 8atom (GaAs)_{1−x}Ge_{2x} ordered compounds, with x ranging from 0 to 1, as elementary models of alloys. For these systems, we have thus employed DFT to determine structural properties and reaction energies, and QSGW to study optical properties. For the more diluted and more concentrated Ge models, I and III, we have predicted stabilizing clustering effects accompanied by a lowering of the products–reactants excess energy. These two systems are symmetric and additionally characterized by an almost identical lattice parameter. In other words, the calculated excess energy of the two intermediate models (IIa and IIb, x_{Ge} = 0.5), clearly showed that the octet rule violation has lead to the final instability of the alloys. In particular, the larger the number of III–IV and IV–V bonds, the larger the instability of the model. We detected a linear relationship between formation energy and number of bad bonds in the alloys. The relevance of this result stems by the fact that for stoichiometric compounds an elementary model Hamiltonian for the energetics of any alloy with equal numbers of Ge cations and anions as function of the number of bad bonds can be developed.
Our QSGW calculations confirm the bowing of the alloy both at the Γ and L points. We also detected direct relationships between optical and mechanical properties: a diminished cohesion for the intermediate alloys (IIb is even almost metallic), with a sensitive reduction in the bandgap was clearly coupled with an increase in lattice parameter and with a reduced symmetry of these two structures. The reduction in symmetry for the intermediate alloys is also considered the fingerprint of an ordered–disordered phase transition for our alloy models.
Footnotes
^{1} We preliminarily performed calculations on the stability of substitutional Ge donor (Ge_{Ga}), acceptor (Ge_{As}), and Ge pairs in GaAs. We have both compared the stability of Ge_{2} molecule and 2Ge isolated in a 64atom supercell GaAs host. Similarly, we calculated the stability of As_{Ge}, Ga_{Ge}, and GaAs_{Ge2} still in a 64atom supercell Ge host. For sake of consistency, these calculations were performed at the same level of theory of present calculations (PAW/LDA), with same cutoff, reduced kpoint sampling (4^{3}Γ centered), and force convergence threshold which is reduced up to 0.05 eV/Å.
Acknowledgments
This research was supported by a Grant from KAKENHI (#21245004) and the Global COE Program [Chemical Innovation] from the Ministry of Education, Culture, Sports, Science, and Technology of Japan.
Open Access
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
References

Kittel C: Introduction to Solid State Physics. Wiley, New York; 2005.

Levinstein M: Handbook Series on Semiconductor Parameters vol 1, 2. World Scientific, London; 1999.

J. Phys.: Condens. Matter. 2002, 14:9579.
COI number [1:CAS:528:DC%2BD38XovVaksrc%3D]; Bibcode number [2002JPCM...14.9579W]
Publisher Full Text 
Streetman BG, Banerjee S: Solid State electronic Devices. Prentice Hall, New Jersey; 2000.

Phys. Rev. B. 2002, 66:075131. Publisher Full Text

Newman KE, Dow JD, Bunker B, Abels LL, Raccah PM, Ugur S, Xue DZ, Kobayashi A:
Phys. Rev. B. 1989, 39:657.
COI number [1:CAS:528:DyaL1MXhtVCisr8%3D]; Bibcode number [1989PhRvB..39..657N]
Publisher Full Text 
Phys. Rev. B. 1993, 47:1889.
COI number [1:CAS:528:DyaK3sXhtlKltbk%3D]; Bibcode number [1993PhRvB..47.1889O]
Publisher Full Text 
Rodriguez AG, NavarroContreras H, Vidal MA:
Phys. Rev. B. 2001, 63:115328.
Bibcode number [2001PhRvB..63k5328R]
Publisher Full Text 
Surf. Sci.. 1992, 267:87.
COI number [1:CAS:528:DyaK38XksFKhs78%3D]; Bibcode number [1992SurSc.267...87I]
Publisher Full Text 
Phys. Rev. B. 1993, 47:16336.
Bibcode number [1993PhRvB..4716336O]
Publisher Full Text 
Superlattices and Microstruct. 1985, 1:275.
COI number [1:CAS:528:DyaL28XkvVGns7g%3D]; Bibcode number [1985SuMi....1..275N]
Publisher Full Text 
Bowen MA, Redfield AC, Froelich DV, Newman KE, Allen RE, Dow JD:
J. Vac. Sci. Technol. B. 1983, 1(3):747.
COI number [1:CAS:528:DyaL3sXkvFChu7c%3D]
Publisher Full Text 
Barnett SA, Ray MA, Lastras A, Kramer B, Greene JE, Raccah PM, Abels LL:
Electron. Lett.. 1982, 18:891.
COI number [1:CAS:528:DyaL38XmtVGru7w%3D]; Bibcode number [1982ElL....18..891B]
Publisher Full Text 
Alferov ZhI, Zhingarev MZ, Konnikov SG, Mokan II, Ulin VP, Umanskii VE, Yavich BS:

Norman AG, Olson JM, Geisz JF, Moutinho HR, Mason A, AlJassim MM, Vernon M:
Appl. Phys. Lett.. 1999, 74:1382.
COI number [1:CAS:528:DyaK1MXhsFWrtbY%3D]; Bibcode number [1999ApPhL..74.1382N]
Publisher Full Text 
Banerjee I, Chung DW, Kroemer H:
Appl. Phys. Lett.. 1985, 46:494.
COI number [1:CAS:528:DyaL2MXhsVWisLc%3D]; Bibcode number [1985ApPhL..46..494B]
Publisher Full Text 
J. Appl. Phys.. 1974, 45:3690.
COI number [1:CAS:528:DyaE2cXls1Sgs7k%3D]; Bibcode number [1974JAP....45.3690N]
Publisher Full Text 
Baker SH, Bayliss SC, Gurman SJ, Elgun N, Bates JS, Davis EA:
J. Phys.: Condens. Matter. 1993, 5:519.
COI number [1:CAS:528:DyaK3sXhtlaltro%3D]; Bibcode number [1993JPCM....5..519B]
Publisher Full Text 
Phys. Rev. B. 1983, 27:7495.
COI number [1:CAS:528:DyaL3sXktlWjtb8%3D]; Bibcode number [1983PhRvB..27.7495N]
Publisher Full Text 
Rodriguez AG, NavarroContreras H, Vidal MA:
Appl. Phys. Lett.. 2000, 77:2497.
COI number [1:CAS:528:DC%2BD3cXntlSqu70%3D]; Bibcode number [2000ApPhL..77.2497R]
Publisher Full Text 
Stern EA, Ellis F, Kim K, Romano L, Shah SI, Greene JE:
Phys. Rev. Lett.. 1985, 54:905.
COI number [1:CAS:528:DyaL2MXhsFyns74%3D]; Bibcode number [1985PhRvL..54..905S]
PubMed Abstract  Publisher Full Text 
Phys. Rev. B. 1985, 32:1019.
COI number [1:CAS:528:DyaL2MXkvFKltb0%3D]; Bibcode number [1985PhRvB..32.1019K]
Publisher Full Text 
Phys. Rev. B. 1987, 35:2767.
COI number [1:CAS:528:DyaL2sXhs1Klsrk%3D]; Bibcode number [1987PhRvB..35.2767D]
Publisher Full Text 
Phys. Rev. B. 1987, 35:3823.
COI number [1:CAS:528:DyaL2sXhvVKktLw%3D]; Bibcode number [1987PhRvB..35.3823H]
Publisher Full Text 
Phys. Rev. Lett.. 1984, 53:830.
COI number [1:CAS:528:DyaL2cXlt1Srsr4%3D]; Bibcode number [1984PhRvL..53..830H]
Publisher Full Text 
Fiz. Tekh. Poluprovodn.. 1982, 16:890.
Sov. Phys. Semicond. 16, 570 (1982)

Phys. Rev. B. 1995, 51:10462.
COI number [1:CAS:528:DyaK2MXlsVSgtrc%3D]; Bibcode number [1995PhRvB..5110462M]
Publisher Full Text 
Seong Y, Norman AG, Ferguson IT, Booker GR:
J. Appl. Phys.. 1993, 73:8227.
COI number [1:CAS:528:DyaK3sXlt1ymu7w%3D]; Bibcode number [1993JAP....73.8227S]
Publisher Full Text 
Phys. Rev. B. 1981, 23:5048.
COI number [1:CAS:528:DyaL3MXktFejurk%3D]; Bibcode number [1981PhRvB..23.5048P]
Publisher Full Text 
Phys. Rev. Lett.. 1980, 45:566.
COI number [1:CAS:528:DyaL3cXltlyhtrY%3D]; Bibcode number [1980PhRvL..45..566C]
Publisher Full Text 
Burke K, Perdew JP: Y. Wang. Electronic Density Functional Theory: Recent Progress and New Directions, Plenum Press, New York; 1998.

Perdew JP: Electronic Structure of Solids 91. Akademie Verlag, Berlin; 1991.

Perdew JP, Chevary JA, Vosko SH, Jackson KA, Pederson MR, Singh DJ, Fiolhais C:
Phys. Rev. B. 1992, 46:6671.
COI number [1:CAS:528:DyaK38XlvFyks7c%3D]; Bibcode number [1992PhRvB..46.6671P]
Publisher Full Text 
Phys. Rev. B. 1994, 50:17953.
Bibcode number [1994PhRvB..5017953B]
Publisher Full Text 
Phys. Rev. B. 1999, 59:1758.
COI number [1:CAS:528:DyaK1MXkt12nug%3D%3D]; Bibcode number [1999PhRvB..59.1758K]
Publisher Full Text 
Methfessel M, van Schilfgaarde M, Casali RA: Electronic Structure and Physical Properties of Solids: The Uses of the LMTO Method, Lecture Notes in Physics 535. SpringerVerlag Berlin; 2000:114.

Phys. Rev.. 1965, 139:A796.
Bibcode number [1965PhRv..139..796H]
Publisher Full Text 
GomezAbal R, Li X, Scheffler M, AmbroschDraxl C:
Phys. Rev. Lett.. 2008, 101:106404.
Bibcode number [2008PhRvL.101j6404G]
PubMed Abstract  Publisher Full Text 
Solid State Commun.. 2002, 121:461.
COI number [1:CAS:528:DC%2BD38XhsFagsbY%3D]; Bibcode number [2002SSCom.121..461K]
Publisher Full Text 
Kotani T, van Schilfgaarde M, Faleev SV:
Phys. Rev. B. 2007, 76:165106.
Bibcode number [2007PhRvB..76p5106K]
Publisher Full Text 
van Schilfgaarde M, Kotani T, Faleev SV:
Phys. Rev. B. 2006, 74:245125.
Bibcode number [2006PhRvB..74x5125V]
Publisher Full Text 
van Schilfgaarde M, Kotani T, Faleev SV:
Phys. Rev. Lett.. 2006, 96:226402.
Bibcode number [2006PhRvL..96v6402V]
PubMed Abstract  Publisher Full Text 
Shishkin M, Marsman M, Kresse G:
Phys. Rev. Lett.. 2007, 99:246403.
COI number [1:STN:280:DC%2BD1c%2FnslSqsw%3D%3D]; Bibcode number [2007PhRvL..99x6403S]
PubMed Abstract  Publisher Full Text 
Chantis AN, van Schilfgaarde M, Kotani T:
Phys. Rev. Lett.. 2006, 96:086405.
Bibcode number [2006PhRvL..96h6405C]
PubMed Abstract  Publisher Full Text 
Viña L, Logothetidis S, Cardona M:
Phys. Rev. B. 1984, 30:1979.
Bibcode number [1984PhRvB..30.1979V]
Publisher Full Text 
Lautenschlager P, Garriga M, Logothetidis S, Cardona M:
Phys. Rev. B. 1987, 35:9174.
COI number [1:CAS:528:DyaL2sXlt1Wgtbk%3D]; Bibcode number [1987PhRvB..35.9174L]
Publisher Full Text 
Yeh CY, Lu ZW, Froyen S, Zunger A:
Phys. Rev. B. 1992, 46:10086.
COI number [1:CAS:528:DyaK38Xms1Wqtrw%3D]; Bibcode number [1992PhRvB..4610086Y]
Publisher Full Text 
J. Phys.: Condens. Matter. 2003, 15:L197.
COI number [1:CAS:528:DC%2BD3sXjsVans7g%3D]; Bibcode number [2003JPCM...15L.197W]
Publisher Full Text 
Phys. Rev. B. 1997, 55:4027.
COI number [1:CAS:528:DyaK2sXhtlGhsL4%3D]; Bibcode number [1997PhRvB..55.4027K]
Publisher Full Text 
Hellwege KH, Madelung O: Landolt–Börnstein New Series Group III. Springer, Berlin; 1982. PubMed Abstract  Publisher Full Text

Singh J: Physics of Semiconductors and their Heterostructures. McGraw & Hill, New York; 1993.

Arabi H, Pourghazi A, Ahmadian F, Nourbakhsh Z:
Phys B: Condens Matter. 2006, 373:16.
COI number [1:CAS:528:DC%2BD28XhtVGitLw%3D]; Bibcode number [2006PhyB..373...16A]
Publisher Full Text 
CRC Handbook of Chemistry & Physics. 1997. PubMed Abstract

Phys. Rev. B. 1994, 49:4710.
COI number [1:CAS:528:DyaK2cXit1OnsLc%3D]; Bibcode number [1994PhRvB..49.4710M]
Publisher Full Text 
BautistaHernandez A, PerezArrieta L, Pal U:
J. F. Rivas Silva, Rev. Mex. Fis.. 2003, 49:9.
COI number [1:CAS:528:DC%2BD3sXhvFWruro%3D]

Henkelman G, Arnaldsson A, Jónsson H:
Comput. Mater. Sci.. 2006, 36:354. Publisher Full Text

Sanville E, Kenny SD, Smith R, Henkelman G:
J. Comp. Chem.. 2007, 28:899.
COI number [1:CAS:528:DC%2BD2sXjtFCit78%3D]
Publisher Full Text 
Tang W, Sanville E, Henkelman G:
J. Phys.: Condens. Matter. 2009, 21:084204.
Bibcode number [2009JPCM...21h4204T]
Publisher Full Text 
Phys. Rev. Lett.. 1991, 67:2339.
COI number [1:CAS:528:DyaK3MXms1Ohs7s%3D]; Bibcode number [1991PhRvL..67.2339Z]
PubMed Abstract  Publisher Full Text 
Zh. Eksp. Teor. Fiz.. 1937, 7:19.
COI number [1:CAS:528:DyaA1cXislGk]

Phys. Rev. B. 1996, 54:16676.
COI number [1:CAS:528:DyaK2sXnt1Wr]; Bibcode number [1996PhRvB..5416676M]
Publisher Full Text 
BarYam Y, Kandel D, Domany E:
Phys. Rev. B. 1990, 41:12869.
Bibcode number [1990PhRvB..4112869B]
Publisher Full Text