Research Article  Open Access
Jacky Even, Laurent Pedesseau, Eric Tea, Samy Almosni, Alain Rolland, Cédric Robert, JeanMarc Jancu, Charles Cornet, Claudine Katan, JeanFrançois Guillemoles, Olivier Durand, "Density Functional Theory Simulations of Semiconductors for Photovoltaic Applications: Hybrid OrganicInorganic Perovskites and III/V Heterostructures", International Journal of Photoenergy, vol. 2014, Article ID 649408, 11 pages, 2014. https://doi.org/10.1155/2014/649408
Density Functional Theory Simulations of Semiconductors for Photovoltaic Applications: Hybrid OrganicInorganic Perovskites and III/V Heterostructures
Abstract
Potentialities of density functional theory (DFT) based methodologies are explored for photovoltaic materials through the modeling of the structural and optoelectronic properties of semiconductor hybrid organicinorganic perovskites and GaAs/GaP heterostructures. They show how the properties of these bulk materials, as well as atomistic relaxations, interfaces, and electronic bandlineups in small heterostructures, can be thoroughly investigated. Some limitations of available standard DFT codes are discussed. Recent improvements able to treat manybody effects or based on densityfunctional perturbation theory are also reviewed in the context of issues relevant to photovoltaic technologies.
1. Introduction
Photovoltaic (PV) solar electricity is one of the key technologies of the 21st century to reduce the world’s reliance on fossil fuels for energy generation. Reduced costs and higher conversion efficiencies are of crucial importance to make PVbased technologies economically more competitive. The quest for quality and improved performances for future solar cells has attracted a tremendous research effort over the last decade towards the semiconductor heterostructures, nanostructured materials, and thin films. Various approaches are used ranging from highcost/highperformance IIIV technology, multiple junctions, and concentrator systems to the lowcost thin films [1]. Obviously, the design of novel and/or efficient PV devices requires a deep understanding of underlying material’s properties including chemical composition, mechanical, electrical, and optical properties, which could be achieved from stateoftheart ab initio approaches. Such knowledge is also desirable to reach lowcost PV cells composed of earthabundant elements based materials capable of full recycling.
Among the large panel of available theoretical approaches, the densityfunctional theory (DFT) has become overwhelmingly popular. This success greatly relies on the fact that no input adjustable parameters are needed and that efficient numerical codes exist and a strong versatility especially in the description of the ground state properties of semiconductors and metals. Increase in computing power has afforded further capabilities in system’s size that DFT methods can handle. However, the current limit does not yet reach the window 10000–10 million atoms involved in the active device region of PV cells. Indeed, modern semiconductor optoelectronic devices have a feature size of few nanometers, including nanostructures like quantum wells [2] or quantum dots [3]. Such systems consist of involved two (2D) and/or three (3D) dimensional geometries composed of multiple materials and/or alloys. Nevertheless, DFT methods are useful to get insight into physical phenomena of each component part of the device, that is, separate materials or small heterostructures. For example, a quantitative design of PV requires a reliable prediction of the electronic bandgaps, bandlineups, and effective masses.
It is well known that DFT based on the local density approximation (LDA) [4] or generalized gradient approximation (GGA) [5] does not reproduce accurately the excited states of compounds. Conversely, hybrid methods that include a fraction of HartreeFock exchange may circumvent the bandgap problem, but their performances strongly depend on the material of interest. Interestingly, the HSE06 hybrid functional proposed by Heyd et al. [6] is a good alternative to compute bandgaps, band offsets, or alloys properties [7, 8]; even so, it fails to reproduce the directindirect crossover in GaAsP alloys [8]. Alternatively, reliable results can be obtained from manybody perturbation theory (MBPT), especially within the GW approach (GW, where G stands for Green’s function and W for the screened potential) that may be used in a perturbative scheme [9–11] or selfconsistently [12, 13]. For charge separation, key quantities are the ionization potential and electron affinity, respectively, of the donor and acceptor materials, as they control the relative alignment of electron and hole levels. The drop of the interface’s potential in heterostructures can be estimated by DFT within a superlattice (SL) approximation, which consequently offers a nice estimate of the bandlineup. For a complex stack, a GW treatment is beyond reach, but the DFT potential drop at the interface can be efficiently corrected by the GW eigenvalues obtained for the bulk valence band states [11, 14, 15]. A good estimation of alloys electronic properties is also a difficult task in DFT simulation. Indeed, semiconductors, even conventional ones, undergo a significant bowing of the bandgap; that is, the bandgap energy follows , where is a bowing parameter. Supercell approaches for statistically random alloys are more adequate than virtual crystal approximations. For special alloy compositions, equivalent results can be obtained from DFT with special quasirandom structures (SQS) (small supercells) that reproduce mixing enthalpies and atomic correlations of very large supercells [8]. Noteworthily, SQS models include chemical mixing, strain, and atomic relaxation effects.
In addition, semiempirical methods for mechanical or electronic properties studies, as the valence force field (VFF) and the tight binding approximation [16] or elasticity and the k.p method [17], require accurate electronic parameters as input. They can be derived from DFT calculations or experiment. For example, quantitative estimates of electromechanical tensors of bulk materials are achieved from densityfunctional perturbation theory (DFPT) [18–20]. An efficient use of the “” theorem, where only byproducts of a firstorder perturbation calculation are required, provides second and thirdorder derivatives of the total energy, provided that atomicdisplacement variables have been eliminated. For secondorder derivatives, various physical responses of insulating crystals can be obtained, including elastic constants, linear piezoelectric tensors, and linear dielectric susceptibility, as well as tensor properties related to internal atomic displacements like Born charges and phonons [20, 21]. For thirdorder components, related to physical properties such as nonlinear electrical susceptibilities, nonlinear elasticity, or photoelastic and electrostrictive effects, finite difference techniques and symmetry analysis must be associated with the DFPT method [22]. Finally, excitonic and transport properties can also be studied by DFT methods, but they are beyond the scope of the present paper.
This paper aims at illustrating, with two examples, some of the DFT’s potentialities for PV technologies. First, we will focus on hybrid organic/inorganic perovskites (HOP) that open new routes for solar cells (Section 3.1). In particular, a direct optical transition and isotropic activity are predicted for a model 3D system. Secondly, the investigation of structural, elastic, and electronic properties will be illustrated on IIIV semiconductor heterostructures for GaP/Si pseudosubstrates (Section 3.2). Computational details are given in Section 2.
2. Computational Details
Total energy DFT calculations were carried out using the ABINIT code [20] within LDA or GGA. A planewave basis set with an energy cutoff of 340 and 680 eV for compounds investigated, respectively, in Sections 3.1 and 3.2, was used to expand the electronic wavefunctions. The reciprocal space integration was performed over MonkhorstPack grids [23]. The energy was computed from the linear response method and convergence is accurately reached with tolerance on the residual potential which stems from differences between the input and output potentials. Pseudopotentials were constructed either from the Fritz Haber Institute (FHI) format [24] or using the HartwigsenGoedeckerHutter (HGH) scheme [25]. The dielectric properties were studied from the densityfunctional perturbation theory (DFPT) based on the linear response theory implemented in SIESTA code [26]. All manybody GW calculations were performed at the G_{0}W_{0} level [9]. The quasiparticles energies were converged with respect to the energy cutoffs, the number of kpoints in the Brillouin Zone (BZ), and the number of bands used to compute the dielectric function and the selfenergy. The plasmonpole model was used to describe the dynamic dependence of the screening function [20].
3. Results
3.1. Optoelectronic Properties of Hybrid Organic/Inorganic Perovskites for High Efficiency and LowCost Solar Cells
Hybrid organic/inorganic materials have attracted increasing interest over the past decade due to their potential applications. In particular, hybrid nanomaterials are expected to offer a way towards enhanced performances of optoelectronic devices [27]. The versatility of the organic part affords the possibility of fine tuning material’s optoelectronic properties [28]. For example, it has been shown that the optical spectra of HOP can be easily tailored by varying the organic cation, which improves optical efficiencies and the tuning of the emission wavelength. Among them, selfassembled layered (2D) HOP have recently shown enhanced nonlinear optical properties in microcavities [29]. Moreover, following pioneering recent works [30–32], 3D HOP based on relatively small organic cations have also been shown to drastically improve the photonic conversion in dye sensitized solar cells (DSSC) [33–38] and expected to afford efficiencies up to 20% in the near future [39]. It is even predicted to open “a new era and a new avenue of research and development for lowcost solar cells … likely to push the absolute power conversion efficiency toward that of CIGS (20%) and then toward and beyond that of crystalline silicon (25%)” [40]. Our theoretical work on this rapidly evolving topic shows that these compounds can be considered now as a new class of semiconductors [41]. Indeed, compared to alternative strategies based on inorganic semiconductor quantum dots (QDs) or extremely thin absorbers coated upon the internal surface of a mesoporous TiO_{2} electrode, 3D HOP offer different benefits. Among those, one can note the ease and low temperature of synthesis, the tailoring of electronic properties by chemical substitution, hole transport, their high stability in dry air, and the formation of pnlike heterojunctions on TiO_{2}.
DFT calculations on prototypes 3D or 2D HOP combined with symmetry analysis of the band edge Bloch states show that the ordering of the band edge states is found reversed compared to tetrahedrally bonded semiconductor structures [28, 41]. Moreover, from the computation of Kane’s energy parameters, we have explained the underlying mechanism of the transverse electric (TE) optical activity of 2D HOP resulting from a subtle interplay between the electronic structure and exciton binding energies [28]. Interestingly, it has also been shown that the spinorbit coupling (SOC) induces a large splitting of the conduction bands, in comparison with those of the valence bands of cubic semiconductors. The importance of SOC was also put forward for 3D HOP [41]. To illustrate this point and the understanding gained from DFT methodologies for hybrid compounds, we consider one of the 3D HOP proposed for PV devices. It belongs to the CH_{3}NH_{3}PbX_{3} family (where X is a halogen atom) and exhibits a disordered average cubic phase (space group Pm3m) at room temperature. Structural disorder is associated with both the rotation of ammonium cations and tilt of the lead halide octaedra. CH_{3}NH_{3}PbI_{3} presents an ordered structure at low temperature, which is orthorhombic (space group Pnma), with a cell doubling when compared to the room temperature phase [41, 45] (Figure 1). The electronic degeneracy lifting is associated with the low temperature symmetry breaking and strain into the unit cell: Ǻ, Ǻ, and Ǻ [45]. The low temperature cell strain is larger along the axis. The theoretical study of 3D HOP predicts an isotropic absorption at high temperature for these materials used in the active zone of solar cells [41]. The triply degenerate CB without SOC of 3D HOP states is indeed associated with the vectorial representation of a cubic simple group. For the low temperature orthorhombic Pnma phase of the 3D CH_{3}NH_{3}PbI_{3 }HOP, the number of CB states is doubled (Figure 2(a)). The fundamental electronic transitions for the CH_{3}NH_{3}PbI_{3 }compound with and without SOC are found at the DFTLDA level equal to 0.4 (Figure 2(b)) and 1.4 eV (Figure 2(a)), respectively. These values are a little bit smaller than the values computed at the DFTGGA level [41]. Even though the fundamental transitions calculated without SOC for the low temperature Pnma structure (1.4 eV) compare nicely with the values obtained experimentally (1.5 eV, [30]), the agreement is fortuitous and stems from large error cancellations. SOC effect indeed strongly reduces the bandgap (Figure 2), without modifying the main character of the optical transitions [41]. At the same time, it is known that a DFT ground state computation systematically underestimates the excited states. This deficiency can be solved with the inclusion of manybody effects, namely, the GW selfenergy correction for the bandgap and the subsequent resolution of Bethe Salpeter’s equation to account for the excitonic features. However, such calculations are far beyond available computational resources for large systems. In the corresponding double group including spinors, the CB of the high temperature cubic phase is split by SOC into twofold degenerate states and fourfold degenerate states [41]. A similar effect is predicted for the lowtemperature phase (Figure 2(b)). The SOC splitting (on the order of 1.0 eV) is much larger than the one usually encountered in the valence band (VB) of cubic conventional semiconductors. The CB minimum is associated with the twofold degenerate and odd spinorbit splitoff (SO) states, leading for symmetry reasons to a strong and isotropic optical transition with the even VB states [41].
(a)
(b)
(a)
(b)
These findings differ from those derived for analogue 2D HOP with larger organic cations. In fact, the 2D densities of states obtained from 4FPEPI ((C_{5}H_{11}NH_{3})_{2}PbI_{4}) [28] and that obtained from (C_{18}H_{37}NH_{3})_{2}PbI_{4}, whose structure is shown in Figure 3 [46], show quite different features. The bandgap remains located at the Γpoint are direct and associated with only three active Bloch states: a nondegenerate level for the valenceband maximum (VBM) and two nearly doubly degenerate levels for the conductionband minimum (CBM12) ( Figure 4, [28]). Moreover, due to the 2D character of these layered compounds, no energy dispersion occurs along the ΓX direction (characterizing the stacking axis in real space), illustrating the PbIPb bond breaking and the dielectric mismatch between inorganic and organic sheets. The fundamental transition of hybrid organic/inorganic layered (2D) perovskites is indeed expected to display a TE character, and it could be even further enhanced by orientational disorder introduced by the organic layer [28]. Contrarily, in the CH_{3}NH_{3}PbI_{3} HOP, an almost isotropic optical activity is expected from symmetry that should also be enhanced at room temperature with the disorder observed for the cubic phase [41]. The computation of the relative imaginary dielectric permittivity is computationally involved with such large systems. It is then limited in this work to a calculation without SOC, for the 3D HOP and for a 2D HOP compound with a shorter alkyl chain (C_{5}H_{11}NH_{3})_{2}PbI_{4} [41, 47]. Relative imaginary dielectric permittivity computed at the GGAlevel presented in Figure 5 for CH_{3}NH_{3}PbI_{3} (a) and (C_{5}H_{11}NH_{3})_{2}PbI_{4} (b) confirms our earlier theoretical predictions based on symmetry analyses [28, 41]. The results are given for light polarizations along the three crystallographic axes (1–3) or for a polycrystalline sample in the case of CH_{3}NH_{3}PbI_{3} (a) and TE or TM polarizations for (C_{5}H_{11}NH_{3})_{2}PbI_{4} (b).
(a)
(b)
(a)
(b)
3.2. IIIV Semiconductors Heterostructures for GaP/Si Pseudosubstrates
High efficiency PV technology is based on IIIV materials growth on expensive Ge or GaAs substrates [1, 48, 49]. With the combination of bulk GaInP and GaAs, 30% conversion efficiencies are attained. Higher conversion efficiencies than 40% can be reached using triple GaInP/Ga(In)As/Ge junctions in solar cells. Efficient tandem cells based on Si substrates have been recently proposed [50–52]. They are based on alloyed heterostructures GaAsP(N)/GaP(N)/Si doublepnjunctions. The key point concerns a coherent growth of IIIV semiconductors on a Si substrate. An initial deposition of thin GaP layer (with a lattice mismatch equal to 0.36%) and a subsequent overgrowth of the dilutednitride GaAsPN compound latticematched to Si are a possible solution to obtain a 1.7 eV/1.1 eV tandem cell. Elaborated growth strategies have however to be used in order to overcome the problems related to extended structural defects formation like antiphase domains and MicroTwins [50–52].
GaAsP(N)/GaP(N) material combinations are good candidates to reach optimal absorption at 1.7 eV [53, 54]. Accurate band structure calculations of strained Ga(NAsP) heterostructures have been performed within the framework of the extendedbasis sp3d5 tightbinding model [16]. DFT calculations are helpful to yield comparisons and band parameters for semiempirical approaches. They are however limited for the simulation of ultrashort superlattices or binary alloys. Figure 6 shows the atomic positions of GaAs(4)/GaP(4) SLs grown along the direction. The GaAs layer is strained onto a GaP substrate where the DFT lattice parameter is of Ǻ, to be compared to the experimental value of 5.44 Ǻ at K [55]. After relaxation, the lattice constant of GaP is recovered almost everywhere along the direction in the GaP layer, whereas the GaAs layer is strongly strained along the direction with an average lattice parameter Ǻ. For GaAs bulk, Ǻ in our DFT simulations, in very good agreement with the experimental value of 5.65 Ǻ, at low temperature [55]. It should be noted that fairly corresponds to the value estimated from the linear elasticity theory using experimental elastic constants ( GPa and GPa).
DFT offers the possibility to study the structural properties beyond classical elasticity, especially close to the interfaces between materials [11, 15]. Figure 7 represents the variation of the atomic interlayer spacing close to the PGaAs interface of the relaxed stack shown in Figure 6. The same stack is simulated using the semiempirical VFF approach [16]. VFF is optimized to describe semiconductor heterostructures at room temperature. Thus, the VFF atomic interlayer spacing values are then shifted to slightly higher values when compared to the DFT results. Both methods show nevertheless that bond stretching undergoes a small fluctuation close to the interface. This fluctuation only extends over roughly one monolayer in the two layers apart from the interface. The layer spacing fluctuation can be interpreted by using an atomic charge analysis. Indeed, the Hirshfeld population analysis [20, 42] is based on the separation of the electron density in real space close to atoms. Conceptually simple, it should be however used with some caution [43]. It was used to describe qualitatively the variation of atomic polarizabilities and atomic charges in IIIV semiconductor nanostructures with heteropolar bonds [44]. Figure 8 shows the variation of the atomic charges in shortperiod GaAs(4)/GaP(4) SLs. Atomic charges calculated for P (−0.15 e) and As (−0.18 e) atoms correspond to the ones obtained in GaP and GaAs bulk, respectively. This also applies for Ga atoms ( e or +0.18 e), except for those located exactly at the interface (+0.16 e). The stretching of the PGa and GaAs bonds connected to the Ga atom at the interface can thus be related to their reduced polarizabilities. This effect is well reproduced by the VFF computations. It should be stressed that advanced empirical force field models are now developed beyond the VFF method in order to get a fine agreement with DFT results for elastic constants and LOTO phonon splitting [56].
The drop of the interface’s potential between materials in heterostructures is a very important parameter for PV structure’s design. Modern theory of bandlineup in semiconductors relies on the model solid approach of van de Walle and Martin [57–59]. The calculation of the valence bandlineup between two semiconductors A (GaAs) and B (GaP) is based on the construction of a superlattice with A/B pseudomorphic interfaces as presented in Figures 6 and 9. The drop at the A/B interface and the planar averaged ab initio potential is then calculated [57]. To complete the computation of valence band discontinuity , one has also to perform the band calculations for bulk and shift the averaged ab initio potentials and valence band maxima, and . The valence band discontinuity is then calculated by transitivity . The conduction band energy levels are represented schematically in Figure 9, but they are not used in the valence bandlineup computation. The valence bandlineup must also be corrected to account for strain and SOC effects [59] or for the valence band selfenergy corrections at the G_{0}W_{0} level [11, 15]. The oneshot G_{0}W_{0} selfenergy correction is known to imperfectly correct the bandgap calculated at the LDA level. The quality of the G_{0}W_{0} approximation is indeed closely tied to the quality of the LDA eigenfunctions at the starting point. The quasiparticle selfconsistent GW (QSGW) method [12] results in stateoftheart predictions of excitedstate properties. However, the additional corrections associated with QSGW selfconsistency mostly affect conduction band levels in conventional semiconductors [12].
We apply in the following the model solid approach for the GaAs/GaP interface and would like to point out that a similar LDA + G_{0}W_{0} study was already presented for the InAs/GaP interface taking into account strain and SOC [15]. The LDA + G_{0}W_{0} study of unstrained bulk GaAs and GaP (Figure 10) yields as expected a direct bandgap for GaAs and an indirect bandgap for GaP. A detailed inspection of the G_{0}W_{0} corrections (Table 1) demonstrates that they are required to recover the right order for the X and L bandgaps in GaP. In addition, the effects are more pronounced for the two compounds at Γ than at X. This point is very important for GaAsP(N)/GaP(N) SL [53] and (In)GaAs/GaP quantum dots [60], where the observed optical transitions result from a fine competition between the Γ and X conduction band states. The DFT simulations carried out only at the LDA level for superlattices are then not expected to yield accurate predictions especially for the conduction subbands. As DFT + LDA fails to reproduce correctly the electron effective masses [12], the quantum confinement is not indeed reproduced correctly for small nanostructures.

LDA + G_{0}W_{0} studies have been recently applied directly to very shortperiod SL instead of using the corrections for bulk materials [14, 61]. Major error in that case however lies in the assumption that the selfenergy in the central part of the superlattice can already be identified with the selfenergy of the bulk. The GW selfenergy contains indeed the huge image effect due to the dielectric discontinuities at the interfaces of shortperiod SL. This can be partially corrected using a semiempirical approach based on the dielectric constant discontinuity [61]. In this work, we present LDA + G_{0}W_{0} studies for bulk GaAs and GaP, but taking into account strain effects. The computed G_{0}W_{0 }valence band shifts are equal to −0.08 eV and −0.25 eV for unstrained GaAs and GaP, respectively. The interplay of strain and manybody effects can be investigated for bulk materials. LDA + G_{0}W_{0} calculations were performed for a biaxially strained GaAs on GaP (001) (Figure 10). In order to understand the effect of the strain state, the strain tensor was decomposed into a hydrostatic component and a pure biaxial component . The energy shifts and definitions of deformation potentials are given in Figure 11. In addition, the interpretation of the electronic band diagram shown in Figure 10 must take into account the band folding effect associated with the doubling of the primitive cell from the rhombohedral system to the tetragonal one. Electronic energy bandgaps deduced for strained GaAs from LDAG_{0}W_{0} at the and Γ points are almost equal. Calculations based on experimental deformation potentials predict an indirect bandgap at point (Table 2, [53, 62]). The G_{0}W_{0} valence band shift for GaAs is equal to −0.19 eV.

The planar averaged and fully averaged ab initio potentials are calculated independently in each crystal bulk, for example, and , and for the GaAs(8)/GaP(8) SL (Figure 9, [20, 64]). The shifted planar averaged potentials for bulk GaP and strained bulk GaAs are shown in Figure 12, respectively, in green and red. The potential curve of GaP has been downshifted by −0.66 eV to fit into the respective part far from interface of the potential curve of the GaAs(8)/GaP(8) SL. The potential curve of strained GaAs has been upshifted by 0.62 eV. The effective lineup of the LDA averaged potentials for GaP and strained GaAs can then be deduced: . Taking into account the energy shifts for bulk and , the valence bandlineup in a GaAs/GaP heterostructure is finally computed at the LDA level:
Adding the valence band G_{0}W_{0} corrections, one finally ends up with a valence bandlineup value of . The GaAs/GaP valence bandlineup calculated from experiment [55] including strain effects and without SOC is equal to 0.55 eV. Another study using a similar method at the LDA level with a somewhat smaller plane wave cutoff (220 eV) yields a comparative value of 0.61 eV [65]. As three independent DFT calculations are needed in the model solid approach (Figure 9), such an agreement is satisfying.
Finally, it is possible to calculate at the LDA level the valence band spinorbit splitting of GaAs and GaP using pseudopotentials constructed from the HGH scheme [25]. Theoretical results are in good agreement with experiment [55] for GaAs 0.35 eV (exp. 0.34 eV) and GaP 0.04 eV (exp. 0.08 eV). Stateoftheart DFT calculations have been recently developed in order to fully take into account SOC in the singleparticle Green function G and the screened interaction W of the GW approximation [66]. It is possible in such a way to correctly predict the reverse band ordering of IIVI semiconductor alloys containing Hg atoms. These compounds are important for the design of topological insulators [67, 68].
4. Conclusions
We have demonstrated that DFT simulations yield relevant information on PV materials. Quantitative design of PV devices is beyond the possibility of nowadays DFT codes but greatly benefits from reliable physical predictions of structural and optoelectronic properties of bulk materials. Studies of small heterostructures may give also additional information on atomistic relaxations, interfaces, and electronic bandlineups. DFT methods are efficiently used in connection with atomistic simulation tools, like the VFF and tight binding codes. Much effort has been done recently by the DFT community to take into account precisely manybody effects for bulk and small nanostructures or to develop a densityfunctional perturbation theory. These improvements are also promising for the future of PV materials.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was performed using HPC resources from GENCI CINES and IDRIS 2013 (c2013096724). It is also supported by Agence Nationale pour la Recherche (PEROCAI Project ANR201004 and MENHIRS Project ANR (Grant no. ANR2011PRGE00701)).
References
 A. Feltrin and A. Freundlich, “Material considerations for terawatt level deployment of photovoltaics,” Renewable Energy, vol. 33, no. 2, pp. 180–185, 2008. View at: Publisher Site  Google Scholar
 M. Guézo, S. Loualiche, J. Even et al., “Ultrashort, nonlinear, optical time response of Fedoped InGaAs/InP multiple quantum wells in 1.55μm range,” Applied Physics Letters, vol. 82, no. 11, pp. 1670–1672, 2003. View at: Publisher Site  Google Scholar
 C. Cornet, C. Labbé, H. Folliot et al., “Timeresolved pump probe of 1.55μm InAs/InP quantum dots under high resonant excitation,” Applied Physics Letters, vol. 88, no. 17, Article ID 171502, 2006. View at: Publisher Site  Google Scholar
 J. P. Perdew and Y. Wang, “Accurate and simple analytic representation of the electrongas correlation energy,” Physical Review B: Condensed Matter and Materials Physics, vol. 45, Article ID 13244, 1992. View at: Publisher Site  Google Scholar
 J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1996. View at: Google Scholar
 J. Heyd, G. E. Scuseria, and M. Ernzerhof, “Hybrid functionals based on a screened Coulomb potential,” Journal of Chemical Physics, vol. 118, no. 8207, Article ID 219906, 2006. View at: Google Scholar
 A. Wadehra, J. W. Nicklas, and J. W. Wilkins, “Band offsets of semiconductor heterostructures: a hybrid density functional study,” Applied Physics Letters, vol. 97, Article ID 092119, 2010. View at: Google Scholar
 J. W. Nicklas and J. W. Wilkins, “Accurate ab initio predictions of IIIV directindirect band gap crossovers,” Applied Physics Letters, vol. 97, Article ID 091902, 13 pages, 2010. View at: Google Scholar
 M. P. Surh, S. G. Louie, and M. L. Cohen, “Quasiparticle energies for cubic BN, BP, and BAs,” Physical Review B: Condensed Matter and Materials Physics, vol. 43, no. 11, pp. 9126–9132, 1991. View at: Publisher Site  Google Scholar
 N. Chimot, J. Even, H. Folliot, and S. Loualiche, “Structural and electronic properties of BAs and B_{x}Ga_{1x}As, BxIn_{1x}As alloys,” Physica B: Condensed Matter, vol. 364, pp. 263–272, 2005. View at: Publisher Site  Google Scholar
 P. Prodhomme, F. FontaineVive, A. V. D. Geest, P. Blaise, and J. Even, “Ab initio calculation of effective work functions for a TiN/HfO_{2}/SiO_{2}/Si transistor stack,” Applied Physics Letters, vol. 99, no. 2, Article ID 022101, 2011. View at: Publisher Site  Google Scholar
 M. van Schilfgaarde, T. Kotani, and S. Faleev, “Quasiparticle selfconsistent GW theory,” Physical Review Letters, vol. 96, no. 22, Article ID 226402, 2006. View at: Publisher Site  Google Scholar
 J. Vidal, F. Trani, F. Bruneval, M. A. L. Marques, and S. Botti, “Effects of electronic and lattice polarization on the band structure of delafossite transparent conductive oxides,” Physical Review Letters, vol. 104, no. 13, Article ID 136401, 2010. View at: Publisher Site  Google Scholar
 R. Shaltaf, G.M. Rignanese, X. Gonze, F. Giustino, and A. Pasquarello, “Band offsets at the Si/SiO_{2} interface from manybody perturbation theory,” Physical Review Letters, vol. 100, no. 18, Article ID 186401, 2008. View at: Publisher Site  Google Scholar
 L. Pedesseau, J. Even, A. Bondi et al., “Theoretical study of highly strained InAs material from firstprinciples modelling: application to an ideal QD,” Journal of Physics D: Applied Physics, vol. 41, no. 16, Article ID 165505, 2008. View at: Publisher Site  Google Scholar
 C. Robert, M. Perrin, C. Cornet, J. Even, and J. M. Jancu, “Atomistic calculations of Ga(NAsP)/GaP(N) quantum wells on silicon substrate: band structure and optical gain,” Applied Physics Letters, vol. 100, no. 11, Article ID 111901, 2012. View at: Publisher Site  Google Scholar
 C. Cornet, A. Schliwa, J. Even et al., “Electronic and optical properties of InAsInP quantum dots on InP(100) and InP(311) B substrates: theory and experiment,” Physical Review B: Condensed Matter and Materials Physics, vol. 74, no. 3, Article ID 035312, 2006. View at: Publisher Site  Google Scholar
 S. Baroni, P. Gianozzi, and A. Testa, “Elastic constants of crystals from linearresponse theory,” Physical Review Letters, vol. 59, p. 2662, 1987. View at: Publisher Site  Google Scholar
 X. Gonze and J.P. Vigneron, “Densityfunctional approach to nonlinearresponse coefficients of solids,” Physical Review B: Condensed Matter and Materials Physics, vol. 39, no. 18, pp. 13120–13128, 1989. View at: Publisher Site  Google Scholar
 X. Gonze, B. Amadon, P.M. Anglade et al., “ABINIT: firstprinciples approach to material and nanosystem properties,” Computer Physics Communications, vol. 180, no. 12, pp. 2582–2615, 2009. View at: Publisher Site  Google Scholar
 J. Even, F. Doré, C. Cornet, L. Pedesseau, A. Schliwa, and D. Bimberg, “Semianalytical evaluation of linear and nonlinear piezoelectric potentials for quantum nanostructures with axial symmetry,” Applied Physics Letters, vol. 91, no. 12, Article ID 122112, 2007. View at: Google Scholar
 L. Pedesseau, C. Katan, and J. Even, “On the entanglement of electrostriction and nonlinear piezoelectricity in noncentrosymmetric materials,” Applied Physics Letters, vol. 100, no. 3, Article ID 031903, 2012. View at: Publisher Site  Google Scholar
 H. J. Monkhorst and J. D. Pack, “Special points for Brillouinzone integrations,” Physical Review B: Condensed Matter and Materials Physics, vol. 13, no. 12, pp. 5188–5192, 1976. View at: Publisher Site  Google Scholar
 M. Fuchs and M. Scheffler, “Ab initio pseudopotentials for electronic structure calculations of polyatomic systems using densityfunctional theory,” Computer Physics Communications, vol. 119, no. 1, pp. 67–98, 1999. View at: Publisher Site  Google Scholar
 C. Hartwigsen, S. Goedecker, and J. Hutter, “Relativistic separable dualspace Gaussian pseudopotentials from H to Rn,” Physical Review B: Condensed Matter and Materials Physics, vol. 58, no. 7, pp. 3641–3662, 1998. View at: Google Scholar
 J. M. Soler, E. Artacho, J. D. Gale et al., “The SIESTA method for ab initio orderN materials simulation,” Journal of Physics: Condensed Matter, vol. 14, no. 11, p. 2745, 2002. View at: Publisher Site  Google Scholar
 D. B. Mitzi, S. Wang, C. A. Feild, C. A. Chess, and A. M. Guloy, “Conducting layered organicinorganic halides containing (110)oriented perovskite sheets,” Science, vol. 267, no. 5203, pp. 1473–1476, 1995. View at: Google Scholar
 J. Even, L. Pedesseau, M.A. Dupertuis, J.M. Jancu, and C. Katan, “Electronic model for selfassembled hybrid organic/perovskite semiconductors: reverse band edge electronic states ordering and spinorbit coupling,” Physical Review B: Condensed Matter and Materials Physics, vol. 86, no. 20, Article ID 205301, 2012. View at: Publisher Site  Google Scholar
 Y. Wei, J. S. Lauret, L. Galmiche, P. Audebert, and E. Deleporte, “Strong excitonphoton coupling in microcavities containing new fluorophenethylamine based perovskite compounds,” Optics Express, vol. 20, no. 9, pp. 10399–10405, 2012. View at: Publisher Site  Google Scholar
 A. Kojima, K. Teshima, Y. Shirai, and T. Miyasaka, “Organometal halide perovskites as visiblelight sensitizers for photovoltaic cells,” Journal of the American Chemical Society, vol. 131, no. 17, pp. 6050–6051, 2009. View at: Publisher Site  Google Scholar
 M. M. Lee, J. Teuscher, T. Miyasaka, T. N. Murakami, and H. J. Snaith, “Efficient hybrid solar cells based on mesosuperstructured organometal halide perovskites,” Science, vol. 338, no. 6107, pp. 643–647, 2012. View at: Publisher Site  Google Scholar
 H. S. Kim, C. R. Lee, J. H. Im et al., “Lead iodide perovskite sensitized allsolidstate submicron thin film mesoscopic solar cell with efficiency exceeding 9%,” Scientific Reports, vol. 2, article 591, 2012. View at: Google Scholar
 J. H. Heo, S. H. Im, J. H. Noh et al., “Efficient inorganicorganic hybrid heterojunction solar cells containing perovskite compound and polymeric hole conductors,” Nature Photonics, vol. 7, pp. 487–492, 2013. View at: Google Scholar
 J. Burschka, N. Pellet, S. Moon et al., “Sequential deposition as a route to highperformance perovskitesensitized solar cells,” Nature, vol. 499, no. 7458, pp. 316–319, 2013. View at: Publisher Site  Google Scholar
 M. Liu, M. B. Johnston, and H. J. Snaith, “Efficient planar heterojunction perovskite solar cells by vapour deposition,” Nature, vol. 501, pp. 395–398, 2013. View at: Publisher Site  Google Scholar
 J. Even, L. Pedesseau, J. M. Jancu, and C. Katan, “DFT and k · p modelling of the phase transitions of lead and tin halide perovskites for photovoltaic cells,” Physica Status Solidi—Rapid Research Letters, vol. 8, pp. 31–35, 2014. View at: Publisher Site  Google Scholar
 S. D. Stranks, G. E. Eperon, G. Grancini et al., “Electronhole diffusion lengths exceeding 1 micrometer in an organometal trihalide perovskite absorber,” Science, vol. 342, no. 6156, pp. 341–344, 2013. View at: Publisher Site  Google Scholar
 G. Xing, N. Mathews, S. Sun et al., “Longrange balanced electron and holetransport lengths in organicinolongrange balanced electron and holetransport lengths in organicinorganic CH_{3}NH_{3}PbI_{3},” Science, vol. 342, no. 6156, pp. 344–347, 2013. View at: Publisher Site  Google Scholar
 N. G. Park, “Organometal perovskite light absorbers toward a 20% efficiency lowcost solidstate mesoscopic solar cell,” Journal of Physical Chemistry Letters, vol. 4, pp. 2423–2429, 2013. View at: Publisher Site  Google Scholar
 H. J. Snaith, “Perovskites: the emergence of a new era for lowcost, highefficiency solar cells,” Journal of Physical Chemistry Letters, vol. 4, no. 21, pp. 3623–3630, 2013. View at: Publisher Site  Google Scholar
 J. Even, L. Pedesseau, J. M. Jancu, and C. Katan, “Importance of spinorbit coupling in hybrid organic/inorganic perovskites for photovoltaic applications,” Journal of Physical Chemistry Letters, vol. 4, no. 17, pp. 2999–3005, 2013. View at: Publisher Site  Google Scholar
 F. L. Hirshfeld, “Bondedatom fragments for describing molecular charge densities,” Theoretica Chimica Acta, vol. 44, no. 2, pp. 129–138, 1977. View at: Publisher Site  Google Scholar
 S. Saha, R. K. Roy, and P. W. Ayers, “Are the Hirshfeld and Mulliken population analysis schemes consistent with chemical intuition?” International Journal of Quantum Chemistry, vol. 109, no. 9, pp. 1790–1806, 2009. View at: Publisher Site  Google Scholar
 A. Krishtal, P. Senet, and C. van Alsenoy, “Origin of the sizedependence of the polarizability per atom in heterogeneous clusters: the case of AlP clusters,” Journal of Chemical Physics, vol. 133, no. 15, Article ID 154310, 2010. View at: Publisher Site  Google Scholar
 T. Baikie, Y. Fang, J. M. Kadro et al., “Synthesis and crystal chemistry of the hybrid perovskite (CH_{3}NH_{3})PbI_{3} for solidstate sensitised solar cell applications,” Journal of Materials Chemistry A, vol. 1, no. 18, pp. 5628–5641, 2013. View at: Publisher Site  Google Scholar
 D. G. Billing and A. Lemmerer, “Synthesis, characterization and phase transitions of the inorganicorganic layered perovskitetype hybrids [(C_{n}H_{2n+1}NH_{3})_{2}PbI_{4}] ($n=12$, 14, 16 and 18),” New Journal of Chemistry, vol. 32, pp. 1736–1746, 2008. View at: Publisher Site  Google Scholar
 D. G. Billing and A. Lemmerer, “Synthesis, characterization and phase transitions in the inorganicorganic layered perovskitetype hybrids [(C_{n}H_{2n+1}NH_{3})_{2}PbI_{4}], ($n=4$, 5 and 6,” Acta Crystallographica B, vol. 63, pp. 735–747, 2007. View at: Publisher Site  Google Scholar
 K. A. Bertness, S. R. Kurtz, D. J. Friedman, A. E. Kibbler, C. Kramer, and J. M. Olson, “29.5%efficient GaInP/GaAs tandem solar cells,” Applied Physics Letters, vol. 65, no. 8, pp. 989–991, 1994. View at: Publisher Site  Google Scholar
 J. F. Geisz and D. J. Friedman, “IIINV semiconductors for solar photovoltaic applications,” Semiconductor Science and Technology, vol. 17, no. 8, pp. 769–777, 2002. View at: Publisher Site  Google Scholar
 W. Guo, A. Bondi, C. Cornet et al., “Thermodynamic evolution of antiphase boundaries in GaP/Si epilayers evidenced by advanced Xray scattering,” Applied Surface Science, vol. 258, no. 7, pp. 2808–2815, 2012. View at: Publisher Site  Google Scholar
 T. Nguyen Thanh, C. Robert, W. Guo et al., “Structural and optical analyses of GaP/Si and (GaAsPN/GaPN)/GaP/Si nanolayers for integrated photonics on silicon,” Journal of Applied Physics, vol. 112, Article ID 053521, 2012. View at: Google Scholar
 T. Quinci, J. Kuyyalil, T. N. Thanh et al., “Defects limitation in epitaxial GaP on bistepped Si surface using UHVCVDMBE growth cluster,” Journal of Crystal Growth, vol. 380, pp. 157–162, 2013. View at: Publisher Site  Google Scholar
 C. Robert, A. Bondi, T. N. Thanh et al., “Room temperature operation of GaAsP(N)/GaP(N) quantum well based lightemitting diodes: effect of the incorporation of nitrogen,” Applied Physics Letters, vol. 98, no. 25, Article ID 251110, 2011. View at: Publisher Site  Google Scholar
 S. Almosni, C. Robert, T. Nguyen Thanh et al., “Evaluation of InGaPN and GaAsPN materials latticematched to Si for multijunction solar cells,” Journal of Applied Physics, vol. 113, Article ID 123509, 2013. View at: Google Scholar
 S. L. Chuang, Physics of Optoelectronic Devices, Wiley Interscience, New York, NY, USA, 1995.
 P. Han and G. Bester, “Interatomic potentials for the vibrational properties of IIIV semiconductor nanostructures,” Physical Review B: Condensed Matter and Materials Physics, vol. 83, Article ID 174304, 2011. View at: Publisher Site  Google Scholar
 C. G. van de Walle and R. M. Martin, “Theoretical study of band offsets at semiconductor interfaces,” Physical Review B: Condensed Matter and Materials Physics, vol. 35, no. 15, pp. 8154–8165, 1987. View at: Publisher Site  Google Scholar
 C. G. van de Walle, “Band lineups and deformation potentials in the modelsolid theory,” Physical Review, vol. 39, p. 1871, 1989. View at: Publisher Site  Google Scholar
 C. G. van de Walle and J. Neugebauer, “Universal alignment of hydrogen levels in semiconductors, insulators and solutions,” Nature, vol. 423, no. 6940, pp. 626–628, 2003. View at: Publisher Site  Google Scholar
 C. Robert, C. Cornet, P. Turban et al., “Electronic, optical, and structural properties of (In,Ga)As/GaP quantum dots,” Physical Review B: Condensed Matter and Materials Physics, vol. 86, no. 20, Article ID 205316, 2012. View at: Publisher Site  Google Scholar
 C. Mitra, B. Lange, C. Freysoldt, and J. Neugebauer, “Quasiparticle band offsets of semiconductor heterojunctions from a generalized marker method,” Physical Review B: Condensed Matter and Materials Physics, vol. 84, no. 19, Article ID 193304, 2011. View at: Publisher Site  Google Scholar
 R. G. Dandrea and A. Zunger, “Firstprinciples study of intervalley mixing: ultrathin GaAs/GaP superlattices,” Physical Review B: Condensed Matter and Materials Physics, vol. 43, no. 11, pp. 8962–8989, 1991. View at: Publisher Site  Google Scholar
 O. Madelung, Ed., Semiconductors: Intrinsic Properties of Group IV Elements and IIIVIIVI and IVII Compounds, vol. 22 of LandoltBörnstein New Series, Group III, part A, Springer, Berlin, Germany, 1987.
 A. Baldereschi, S. Baroni, and R. Resta, “Band offsets in latticematched heterojunctions: a model and firstprinciples calculations for GaAs/AlAs,” Physical Review Letters, vol. 61, p. 734, 1988. View at: Publisher Site  Google Scholar
 M. di Ventra, M. Peressi, and A. Baldereschi, “Chemical and structural contributions to the valenceband offset at GaP/GaAs heterojunctions,” Physical Review B: Condensed Matter and Materials Physics, vol. 54, no. 8, p. 5691, 1996. View at: Google Scholar
 R. Sakuma, C. Friedrich, T. Miyake, S. Blügel, and F. Aryasetiawan, “GW calculations including spinorbit coupling: application to Hg chalcogenides,” Physical Review B: Condensed Matter and Materials Physics, vol. 84, Article ID 085144, 2011. View at: Publisher Site  Google Scholar
 O. A. Pankratov, S. V. Pakhomov, and B. A. Volkov, “Supersymmetry in heterojunctions: bandinverting contact on the basis of Pb_{1x}Sn_{x}Te and Hg_{1x}Cd_{x}Te,” Solid State Communications, vol. 61, no. 2, pp. 93–96, 1987. View at: Google Scholar
 C. L. Kane and J. E. Moore, “Topological insulators,” Physics World, vol. 24, no. 2, pp. 32–36, 2011. View at: Google Scholar
Copyright
Copyright © 2014 Jacky Even et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.