-
Polarization spontaneous and piezo: fundamentals and their implementation in ab initio calculations
Authors:
Pawel Strak,
Pawel Kempisty,
Konrad Sakowski,
Jacek Piechota,
Izabella Grzegory,
Eva Monroy,
Agata Kaminska,
Stanislaw Krukowski
Abstract:
Fundamental properties of spontaneous and piezo polarization are reformulated and critically reviewed. It was demonstrated that Landau definition of polarization as a dipole density could be used to the infinite systems. The difference between the bulk polarization and surface polarity are distinguished thus creating clear identification of both components. The local model of spontaneous polarizat…
▽ More
Fundamental properties of spontaneous and piezo polarization are reformulated and critically reviewed. It was demonstrated that Landau definition of polarization as a dipole density could be used to the infinite systems. The difference between the bulk polarization and surface polarity are distinguished thus creating clear identification of both components. The local model of spontaneous polarization was created and used to calculate spontaneous polarization as the electric dipole density. It was shown that the proposed local model correctly predicts c-axis spontaneous polarization values of the nitride wurtzite semiconductors. It was also shown that the proposed model predicts zero polarization in the plane perpendicular to the c-axis, in accordance with symmetry requirements. In addition, the model results are in accordance with polarization equal to zero for zinc blende lattice. These data confirm the basic correctness of the proposed model. The spontaneous polarization values obtained for all wurtzite III nitrides (BN, AlN, GaN and InN) are in basic agreement with the earlier calculations using Berry phase and slab models of Bernardini et al. {Bernardini et al. Phys Rev B 56 (2001) R10024 & 63 (2001) 193201} but not with Dreyer et al. {Dreyer et al. Phys. Rev X 6 (2016) 021038}. Wurtzite nitride superlattices ab initio calculations were performed to derive polarization-induced fields in the coherently strained lattices showing good agreement with the polarization values. The strained superlattice data were used to determine the piezoelectric parameters of wurtzite nitrides obtaining the values that were in basic agreement with the earlier data. Zinc blende superlattices were also modeled using ab initio calculations showing results that are in agreement with the absence of polarization of all nitrides in zinc blende symmetry.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Augmentation of the Electron Counting Rule with Ising Model
Authors:
Karol Kawka,
Pawel Kempisty,
Konrad Sakowski,
Stanislaw Krukowski,
Michal Bockowski,
David Bowler,
Akira Kusaba
Abstract:
On semiconductor growth surfaces, surface reconstructions appear. Estimation of the reconstructed structures is essential for understanding and controlling growth phenomena. In this study, the stability of a mixture of two different surface reconstructions is investigated. Since the number of candidate structures is enormous, the structures sampled by Bayesian optimization are analyzed. As a resul…
▽ More
On semiconductor growth surfaces, surface reconstructions appear. Estimation of the reconstructed structures is essential for understanding and controlling growth phenomena. In this study, the stability of a mixture of two different surface reconstructions is investigated. Since the number of candidate structures is enormous, the structures sampled by Bayesian optimization are analyzed. As a result, the local electron counting (EC) rule alone was found to be insufficient to explain such stability. Then, augmenting the EC rule, a data-driven Ising model is proposed. The model allows the evaluation of the whole enormous number of candidate structures. The approach is expected to be useful for theoretical studies of such mixtures on various semiconductor surfaces.
△ Less
Submitted 8 February, 2024;
originally announced February 2024.
-
Coulomb contribution to Shockley-Read-Hall (SRH) recombination
Authors:
Konrad Sakowski,
Pawel Strak,
Pawel Kempisty,
Jacek Piechota,
Izabella Grzegory,
Piotr Perlin,
Eva Monroy,
Agata Kaminska,
Stanislaw Krukowski
Abstract:
Defect-mediated nonradiative recombination, known as Shockley-Read-Hall (SRH) recombination is reformulated. The introduced model considers Coulomb attraction between charged deep defect and the approaching free carrier, showing that this effect may cause considerable increase of the carrier velocity approaching the recombination center. The effect considerably increases the carrier capture rates.…
▽ More
Defect-mediated nonradiative recombination, known as Shockley-Read-Hall (SRH) recombination is reformulated. The introduced model considers Coulomb attraction between charged deep defect and the approaching free carrier, showing that this effect may cause considerable increase of the carrier velocity approaching the recombination center. The effect considerably increases the carrier capture rates. It is demonstrated that in the typical semiconductor device or semiconductor medium, the SRH recombination cannot be neglected at low temperatures. The SRH is more effective in the case of low doped semiconductors. Effective screening by mobile carrier density could reduce the effect, leading to SRH rate increase.
△ Less
Submitted 26 January, 2024; v1 submitted 18 October, 2023;
originally announced October 2023.
-
Macrosteps dynamics and the growth of crystals and epitaxial layers
Authors:
Stanislaw Krukowski,
Konrad Sakowski,
Paweł Strak,
Paweł Kempisty,
Jacek Piechota,
Izabella Grzegory
Abstract:
Step pattern stability of the vicinal surfaces during growth was analyzed using various surface kinetic models. It was shown that standard analysis of the vicinal surfaces provides no indication on the possible step coalescence and therefore could not be used to elucidate macrostep creation during growth. A scenario of the instability, leading go macrostep creation was based on the dynamics of the…
▽ More
Step pattern stability of the vicinal surfaces during growth was analyzed using various surface kinetic models. It was shown that standard analysis of the vicinal surfaces provides no indication on the possible step coalescence and therefore could not be used to elucidate macrostep creation during growth. A scenario of the instability, leading go macrostep creation was based on the dynamics of the step train. The critical is step motion at the rear of the train which leads to double and multiple step creation. The condition is that the step density ratio in and out of the train lower than 2 prevents double step formation irrespective of the kinetics. For higher step density ratio low density of the step promotes single step stability. Fast step kinetics from lower terrace stabilizes the single steps slow (high barrier) is promoting step coalescence. The incorporation kinetics from upper terrace role is close to neutral. The creation of double step creates slow the step in front to accelerate and catch the previous double step while those behind catch up the double step creating multistep structure. The multistep are not mobile as the alimentation leads to emission of single step which moves forward. The final structure consist of macrosteps and superterraces with the number of single steps moving forward. Thus the single step motion is essential crystal growth mode despite the presence of the macrosteps. The macrostep are prone to creation of the overhangs which results from surface dynamics coupling to impingement from the mother phase. The angular preferential access of the bulk material to the macrostep edge, leads to diffusive instability. Therefore it is expected that harmful influence of the macrosteps by creation of inclusions and dislocation is stronger during growth from the liquid phase.
△ Less
Submitted 11 January, 2023; v1 submitted 22 March, 2022;
originally announced March 2022.
-
Polarization doping ab initio verification of the concept charge conservation and nonlocality
Authors:
Ashfaq Ahmad,
Pawel Strak,
Pawel Kempisty,
Konrad Sakowski,
Jacek Piechota,
Yoshihiro Kangawa,
Izabella Grzegory,
Michal Leszczynski,
Zbigniew R. Zytkiewicz,
Grzegorz Muziol,
Eva Monroy,
Stanislaw Krukowski,
Agata Kaminska
Abstract:
In this work, we study the emergence of polarization doping in AlxGa1-xN layers with graded composition from a theoretical viewpoint. We demonstrate that the charge conservation law applies for fixed and mobile charges separately, leading to nonlocal compensation phenomena involving bulk fixed and mobile charge and polarization sheet charge at the heterointerfaces. The magnitude of the effect allo…
▽ More
In this work, we study the emergence of polarization doping in AlxGa1-xN layers with graded composition from a theoretical viewpoint. We demonstrate that the charge conservation law applies for fixed and mobile charges separately, leading to nonlocal compensation phenomena involving bulk fixed and mobile charge and polarization sheet charge at the heterointerfaces. The magnitude of the effect allows obtaining technically viable mobile charge density for optoelectronic devices without impurity doping (donors or acceptors). Therefore, it provides an additional tool for the device designer, with the potential to attain high conductivities: high carrier concentrations can be obtained even in materials with high dopant ionization energies, and the mobility is not limited by scattering at ionized impurities.
△ Less
Submitted 22 March, 2022; v1 submitted 2 June, 2021;
originally announced June 2021.
-
Comment to "Topological insulators and superconductors" by Xiao-Liang Qi and Shou-Cheng Zhang Reviews of Modern Physics 83, 1057 (2011)
Authors:
Stanisław Krukowski,
Paweł Kempisty,
Paweł Strak,
Konrad Sakowski
Abstract:
The paper reports on the solution in the semi-space corresponding to topological surface state. It was shown that the exponential surface state results from selection of the energy outside the allowed energy range. In the result, the kinetic energy of the state is negative proving that the proposed state is not physically correct solution.
The paper reports on the solution in the semi-space corresponding to topological surface state. It was shown that the exponential surface state results from selection of the energy outside the allowed energy range. In the result, the kinetic energy of the state is negative proving that the proposed state is not physically correct solution.
△ Less
Submitted 20 February, 2019;
originally announced February 2019.
-
Comment to "Quantum Spin Hall Effect and Topological Phase Transition in HGTe Quantum Wells" by B. Andrei Bernevig, Taylor L. Hughes and Shou-Cheng Zhang Science 314, 1757 (2006)
Authors:
Stanisław Krukowski,
Paweł Kempisty,
Paweł Strak,
Konrad Sakowski
Abstract:
The manuscript contains Comment to "Quantum Spin Hall Effect and Topological Phase Transitio in HGTe Quantum Wells" by B. Andrei Bernevig, Taylor L. Hughes and Shou-Cheng Zhang Science 314, 1757 (2006). The text proves that the quantum state, identified by these authors as the surface state located inside HgTe well is incorrect physically. The kinetic energy of the state is negative. The state can…
▽ More
The manuscript contains Comment to "Quantum Spin Hall Effect and Topological Phase Transitio in HGTe Quantum Wells" by B. Andrei Bernevig, Taylor L. Hughes and Shou-Cheng Zhang Science 314, 1757 (2006). The text proves that the quantum state, identified by these authors as the surface state located inside HgTe well is incorrect physically. The kinetic energy of the state is negative. The state cannot be identified as proper eigenfunction (vector) in the Hilbert space. Therefore the state could be used as a primer in the topological state theory, developed in the following publications.
△ Less
Submitted 13 February, 2019;
originally announced February 2019.
-
Adsorption of ammonia at GaN(0001) surface in the mixed ammonia/hydrogen ambient - a summary of ab initio data
Authors:
Pawel Kempisty,
Stanislaw Krukowski
Abstract:
Adsorption of ammonia at NH3/NH2/H covered GaN(0001) surface was analyzed using results of ab initio calculations. The whole configuration space of partially NH3/NH2/H covered GaN(0001) surface was divided into zones differently pinned Fermi level: at Ga broken bond state for dominantly bare surface (region I), at VBM for NH2 and H covered (region II), and at CBM for NH3 covered surface (region II…
▽ More
Adsorption of ammonia at NH3/NH2/H covered GaN(0001) surface was analyzed using results of ab initio calculations. The whole configuration space of partially NH3/NH2/H covered GaN(0001) surface was divided into zones differently pinned Fermi level: at Ga broken bond state for dominantly bare surface (region I), at VBM for NH2 and H covered (region II), and at CBM for NH3 covered surface (region III). The extensive ab intio calculations show validity of electron counting rule (ECR) for all mixed coverage, for bordering these three regions. The adsorption was analyzed using newly identified dependence of the adsorption energy on the charge transfer at the surface. For region I and II ammonia adsorb dissociatively, disintegrating into H adatom and HN2 radical for large fraction of vacant sites while for high coverage the ammonia adsorption is molecular. The dissociative adsorption energy strongly depends on the Fermi level at the surface (pinned) and in the bulk (unpinned) while the molecular adsorption energy is determined by bonding to surface only, in accordance to the recently published theory. The molecular adsorption is determined by the energy of covalent bonding to the surface. Ammonia adsorption in region III (Fermi level pinned at CBM) leads to unstable configuration both molecular and dissociative which is explained by the fact that Ga-broken bond sites are doubly occupied by electrons. The adsorbing ammonia brings 8 electrons to the surface, necessitating transfer of the electrons from Ga-broken bond state to Fermi level, energetically costly process. Adsorption of ammonia at H-covered site leads to creation of NH2 radical at the surface and escape of H2 molecule. The process energy is close to 0.12 eV, thus not large, but the inverse process is not possible due to escape of the hydrogen molecule.
△ Less
Submitted 24 May, 2014;
originally announced May 2014.
-
General aspects of the vapor growth of semiconductor crystals - a study based on DFT simulations of the NH3/NH2 covered GaN(0001) surface in hydrogen ambient
Authors:
Pawel Kempisty,
Pawel Strak,
Konrad Sakowski,
Stanislaw Krukowski
Abstract:
Vapor growth of semiconductors is analyzed using recently obtained dependence of the adsorption energy on the electron charge transfer between the surface adsorbed species and the bulk [Krukowski et al. J. Appl. Phys. 114 (2013) 063507, Kempisty et al. ArXiv 1307.5778 (2013)]. Ab initio calculations were performed to study the physical properties of GaN(0001) surface in ammonia-rich conditions, i.…
▽ More
Vapor growth of semiconductors is analyzed using recently obtained dependence of the adsorption energy on the electron charge transfer between the surface adsorbed species and the bulk [Krukowski et al. J. Appl. Phys. 114 (2013) 063507, Kempisty et al. ArXiv 1307.5778 (2013)]. Ab initio calculations were performed to study the physical properties of GaN(0001) surface in ammonia-rich conditions, i.e. covered by mixture of NH3 molecules and NH2 radicals. The Fermi level is pinned at valence band maximum (VBM) and conduction band minimum (CBM) for full coverage by NH3 molecules and NH2 radicals, respectively. For the crossover content of ammonia of about 25% monolayer (ML), the Fermi level is unpinned. It was shown that hydrogen adsorption energy depends on the doping in the bulk for the unpinned Fermi level, i.e. for this coverage. Surface structure thermodynamic and mechanical stability criteria are defined and compared. Mechanical stability of the coverage of such surfaces was checked by determination of the desorption energy of hydrogen molecules. Thermodynamic stability analysis indicates that initally equilibrium hydrogen vapor partial pressure steeply increases with NH3 content to attain the crossover NH3/NH2 coverage, i.e. the unpinned Fermi level condition. For such condition the entire range of experimentally accessible pressures belongs showing that vapor growth of semiconductor crystals occurs predominantly for unpinned Fermi level at the surface, i.e. for flat bands. Accordingly, adsorption energy of most species depends on the doping in the bulk that is basis of the possible molecular scenario explaining dependence of the growth and the doping of semiconductor crystals on the doping in the bulk
△ Less
Submitted 19 November, 2013;
originally announced November 2013.
-
Energy of adsorption at semiconducting surfaces with Fermi level differently pinned -- ab initio study
Authors:
Paweł Kempisty,
Paweł Strak,
Konrad Sakowski,
Stanisław Krukowski
Abstract:
It is shown that adsorption energy at semiconductor surfaces critically depends on the charge transfer to or from the adsorbed species. For the processes without the charge transfer, such as molecular adsorption of closed shell systems, the adsorption energy is determined by the bonding only. In the case involving charge transfer, such as open shell systems like metal atoms or the dissociating mol…
▽ More
It is shown that adsorption energy at semiconductor surfaces critically depends on the charge transfer to or from the adsorbed species. For the processes without the charge transfer, such as molecular adsorption of closed shell systems, the adsorption energy is determined by the bonding only. In the case involving charge transfer, such as open shell systems like metal atoms or the dissociating molecules, the energy attains different value for the Fermi level differently pinned. The DFT simulation of adsorption of ammonia at hydrogen covered GaN(0001) confirms these predictions: the molecular adsorption is independent on the coverage while the dissociative process energy varies by several electronvolts.
△ Less
Submitted 22 August, 2013; v1 submitted 22 July, 2013;
originally announced July 2013.
-
Fermi level influence on the adsorption at semiconductor surfaces - ab initio simulations
Authors:
Stanislaw Krukowski,
Pawel Kempisty,
Pawel Strak
Abstract:
Chemical adsorption of the species at semiconductor surfaces is analyzed showing the existence of the two contributions to adsorption energy: bond creation and charge transfer. It is shown that the energy of quantum surface states is affected by the electric field at the surface, nevertheless the potential contribution of electron and nuclei cancels out. The charge transfer contribution is Fermi l…
▽ More
Chemical adsorption of the species at semiconductor surfaces is analyzed showing the existence of the two contributions to adsorption energy: bond creation and charge transfer. It is shown that the energy of quantum surface states is affected by the electric field at the surface, nevertheless the potential contribution of electron and nuclei cancels out. The charge transfer contribution is Fermi level independent for pinned surfaces. Thus for Fermi level pinned at the surface, the adsorption energy is independent on the Fermi energy i.e. the doping in the bulk. The DFT simulations of adsorption of hydrogen at clean GaN(0001) and silicon at SiC(0001) surfaces confirmed independence of adsorption energy on the doping in the bulk. For the Fermi level nonpinned surfaces the charge contribution depends on the position of Fermi level in the bulk. Thus adsorption energy is sensitive to change of the Fermi energy in the bulk, i.e. the doping. The DFT simulations of adsorption of atomic hydrogen at 0.75 ML hydrogen covered GaN(0001) surface confirmed that the adsorption energy may be changed by about 2 eV by the doping change from n- to p-type.
△ Less
Submitted 12 June, 2013;
originally announced June 2013.
-
Foundations of ab initio simulations of electric charges and fields at semiconductor surfaces within slab models
Authors:
Stanisław Krukowski,
Paweł Kempisty,
Paweł Strąk
Abstract:
Semiconductor surfaces were divided into charge categories, i.e. surface acceptor, donor and neutral ones that are suitable for simulations of their properties within a slab model. The potential profiles, close to the charged surface states, accounting for explicit dependence of the point defects energy, were obtained. A termination charge slab model was formulated and analyzed proving that two co…
▽ More
Semiconductor surfaces were divided into charge categories, i.e. surface acceptor, donor and neutral ones that are suitable for simulations of their properties within a slab model. The potential profiles, close to the charged surface states, accounting for explicit dependence of the point defects energy, were obtained. A termination charge slab model was formulated and analyzed proving that two control parameters of slab simulations exist: the slope and curvature of electric potential profiles which can be translated into a surface and volumetric charge density. The procedures of slab model parameter control are described and presented using examples of the DFT simulations of GaN and SiC surfaces showing the potential profiles, linear or curved, depending on the band charge within the slab. It was also demonstrated that the field at the surface may affect some surface properties in a considerable degree proving that verification of this dependence is obligatory for a precise simulation of the properties of semiconductor surfaces.
△ Less
Submitted 9 February, 2013;
originally announced February 2013.
-
High efficiency UV emitters - theoretical investigation of GaN/AlN heterostructures
Authors:
Pawel Strak,
Pawel Kempisty,
Stanislaw Krukowski
Abstract:
Density functional theory simulations were used to obtain physical properties of GaN/AlN system. Combination of these two compounds into multiquantum well (MQW) structure will induce strong electrostatic effect leading to emergence of high magnitude dipole layers at the AlN/GaN interfaces, which were first postulated by {Tersoff Phys. Rev. B 30(8) pp.4874 (1984)} and already identified in GaN/InN…
▽ More
Density functional theory simulations were used to obtain physical properties of GaN/AlN system. Combination of these two compounds into multiquantum well (MQW) structure will induce strong electrostatic effect leading to emergence of high magnitude dipole layers at the AlN/GaN interfaces, which were first postulated by {Tersoff Phys. Rev. B 30(8) pp.4874 (1984)} and already identified in GaN/InN by {Romanowski et al. J. Phys. Chem C 114. 14410 (2010)}. When combining GaN and AlN in to a heterostructure a spatial projection of wavefunctions indicate that valence band offset between states becomes of order of 0.85 V. Systematic analysis of influence of number of Ga atomic layers on the properties of wells have shown that for thickness up to 4 Ga layers, GaN behave as carriers locating potential minimum rather, while for larger thickness it is a standard quantum well. In all cases wells has strongly localized quantum states close to valance band maxima (VBM) and conduction band minimum (CBM). The calculated oscillator strength values rapidly decreases for the well thickness in excess of 8 Ga layers (~21 Å) which indicates that wells for UV emitters should be much thinner than these based on InGaN/GaN systems. The Quantum Confined Stark Effect (QCSE) related changes of the transition energy in function of the geometric arrangement were also obtained.
△ Less
Submitted 7 November, 2012; v1 submitted 29 August, 2012;
originally announced August 2012.
-
On the nature of Surface States Stark Effect at clean GaN(0001) surface
Authors:
Paweł Kempisty,
Stanisław Krukowski
Abstract:
Recently developed model allows for simulations of electric field influence on the surface states. The results of slab simulations show considerable change of the energy of quantum states in the electric field, i.e. Stark Effect associated with the surface (SSSE - Surface States Stark Effect). Detailed studies of the GaN slabs demonstrate spatial variation of the conduction and valence band energy…
▽ More
Recently developed model allows for simulations of electric field influence on the surface states. The results of slab simulations show considerable change of the energy of quantum states in the electric field, i.e. Stark Effect associated with the surface (SSSE - Surface States Stark Effect). Detailed studies of the GaN slabs demonstrate spatial variation of the conduction and valence band energy revealing real nature of SSSE phenomenon. It is shown that long range variation of the electric potential is in accordance with the change of the energy of the conduction and valence bands. However, at short distances from GaN(0001) surface, the valence band follows the potential change while the conduction states energy is increased due to quantum overlap repulsion by surface states. It is also shown that at clean GaN(0001) surface Fermi level is pinned at about 0.34 eV below the long range projection of the conduction band bottom and varies with the field by about 0.31 eV due to electron filling of the surface states.
△ Less
Submitted 30 May, 2012;
originally announced May 2012.
-
Polarization and polarization induced electric field in nitrides - critical evaluation based on DFT studies
Authors:
Pawel Strak,
Pawel Kempisty,
Konrad Sakowski,
Stanislaw Krukowski
Abstract:
Density Functional Theory (DFT) calculations were used to evaluate polarity of group III nitrides, such as aluminum nitride (AlN), gallium nitride (GaN) and indium nitride (InN) providing physically sound quantitative measure of polarity of these materials. Two different approaches to polarization of nitride semiconductors were assessed and the conclusions have been used to develop models. It was…
▽ More
Density Functional Theory (DFT) calculations were used to evaluate polarity of group III nitrides, such as aluminum nitride (AlN), gallium nitride (GaN) and indium nitride (InN) providing physically sound quantitative measure of polarity of these materials. Two different approaches to polarization of nitride semiconductors were assessed and the conclusions have been used to develop models. It was shown that Berry phase formulation of the electron related polarization component provides a number of various solutions, different for various selection of the simulated volume. The electronic part gives saw-like pattern for polarization. A total number of these solutions, related to well known scaling of the geometric phase, is equal to the number of valence electrons in the system. Summation with similar pattern for ionic part provides several polarization values. Standard dipole density formulation depends on the selection of the simulation volume in periodic continuous way. Using a condition of continuous embedding into the infinite medium, and simultaneously, the zero surface charge representation at crystal boundary provides to physically sound solution. This solution is corresponding to maximal and minimal polarization values and also corresponds to different physical termination of the crystal surfaces, either bare or covered by complementary atoms. This change leads to polarization and electric field reversal. The polarization and related fields in finite size systems were obtained.
△ Less
Submitted 16 November, 2011; v1 submitted 20 October, 2011;
originally announced October 2011.
-
Polarization analysis and Density Functional Theory (DFT) simulations of the electric field in InN/GaN multiple quantum wells (MQWs)
Authors:
Z. Romanowski,
P. Kempisty,
K. Sakowski,
S. Krukowski
Abstract:
Results of the first ab initio simulations of InN/GaN multiquantum well (MQW) system are presented. The DFT results confirm the presence of the polarization charge at InN/GaN interfaces, i.e. at polar InN/GaN heterostructures. These results show the potential jumps which is related to the presence of dipole layer at these interfaces. An electrostatic polarization analysis shows that the energy m…
▽ More
Results of the first ab initio simulations of InN/GaN multiquantum well (MQW) system are presented. The DFT results confirm the presence of the polarization charge at InN/GaN interfaces, i.e. at polar InN/GaN heterostructures. These results show the potential jumps which is related to the presence of dipole layer at these interfaces. An electrostatic polarization analysis shows that the energy minimum condition can be used to obtain the field in InN/GaN system, employing standard polarization parameters. DFT results are in good agreement with polarization data confirming the existence of electric field leading to separation of electron and holes in QWs and emergence of Quantum Confined Stark Effect (QCSE).
△ Less
Submitted 16 November, 2009;
originally announced November 2009.