## Acceptor and donor impurity levels in hexagonal-diamond silicon

Marc Túnica, Alberto Zobelli, and Michele Amato ∗

Université Paris-Saclay, CNRS, Laboratoire de Physique des Solides, 91405 Orsay, France

(*michele.amato@universite-paris-saclay.fr)

(Dated: October 22, 2024)

## Abstract

Recent advances in the characterization of hexagonal-diamond silicon (2H-Si) have shown that this material possesses remarkably different structural, electronic and optical properties as compared to the common cubic-diamond (3C) polytype. Interestingly, despite the wide range of physical properties analyzed, to date no study has investigated impurity energy levels in 2H-Si. Here, we present results of ab initio DFT simulations to describe the effect of p- and n-type substitutional doping on the structural and electronic properties of hexagonal-diamond Si (2H-Si). We first provide a detailed analysis of how a given impurity can assume a different local symmetry depending on the host crystal phase. Then, by studying neutral and charged dopants, we carefully estimate donors and acceptors transition energy levels in 2H-Si and compare them with the cubic-diamond (3C) case. In the case of acceptors, the formation energy is always lower in 2H-Si and is associated with a shallower charge transition level with respect to 3C-Si. On the other hand, donors prefer the cubic phase and have transition energies smaller with respect to 2H-Si. Finally, by employing a simple model based on the 2H/3C band offset diagram, we prove the physical validity of our findings and we show how holes can be used to stabilize the 2H-Si phase. Overall, the described doping properties represent a robust starting point for further theoretical and experimental investigations.

## I. INTRODUCTION

Recent developments and advances in semiconductor growth methods demonstrated that the hexagonal-diamond polytype of Si, Ge and SiGe alloys (2H in the Rasmdell notation [1]) can be fabricated in the form of nanowires [2-6]. This crystal structure presents a C v 3 symmetry (P63/mmc space group) with three equal length bonds in the c-plane (0001) and a larger one along the [0001] direction of the hexagonal crystal cell. Compared to the 3C cubic-diamond phase (Fd3m group symmetry) where four equal length bonds induce a Td symmetry, in the 2H phase the atomic ordering along the ⟨ 111 ⟩ cubic direction changes from an ABC to an ABAB stacking.

The main reason for studying these novel group IV polytypes lies in the possibility to observe unexplored electronic and optical properties with respect to those of their 3C counterparts. To give an example, Ge, which in its cubic form is an indirect semiconductor, has a direct band gap in its hexagonal polytype [7, 8]. Yet, hexagonal SiGe alloys [9] and quantum wells [10] have recently attracted a great attention because of their direct band gap and the experimentally observed light emission that can be tailored via the Ge concentration and quantum confinement effects [3, 10, 11].

In the case of Si, while the cubic phase has an indirect band gap of 1.12 eV with a conduction band minimum close to the X point, the 2H structure presents an indirect band gap of about 0.95 eV and a conduction band minimum at the L point [12, 13]. The direct band gap is also reduced from 3.40 eV to 1.63 eV [12]. Indeed, 2H-Si presents a larger optical absorption spectra in the visible range than 3C-Si [14], making it a good candidate for solar panel materials. Furthermore, 3C/2H Si junctions are also of interest for photovoltaics due to their intrinsic tendency to promote charge carrier separation [8, 14, 15].

It is well known that most technologies based on carrier transport in semiconductors (e.g. electronics, embedded systems, light emitting devices and solar cells) become functional only if the material can be appropriately doped. Impurities play a crucial role as they allow to modulate the carrier concentration around the Fermi level. Furthermore, in different type of semiconductor compounds, doping can be employed to stabilize one polytype with respect to the others [16-19]. For instance, Cai et al. [20] recently demonstrated that the direct band gap 2H-Ge phase can be stabilized by incorporating electrons into the 3C-Ge crystal. Finally, most of the experimental grown nanowires can contain significant concentration of extrinsic dopants due to the substrate atomic contamination. This is, for instance, the case of 2H-Ge and 2H-SiGe alloys [3] fabricated with the crystal transfer methods which usually present a high concentration of As atoms ( ∼ 10 18 cm -3 ).

Although the behaviour of impurities in 3C-Si has been widely discussed, to date no study has investigated impurity energy levels in 2H-Si. The experimental characterization of dopants in this polytype still remains challenging because of the difficulties in its fabrications methods. Ab-initio theoretical approaches can be hence an efficient way to investigate this question. A few studies focused on the analysis of the stability of single-point and neutral substitutional defects in 2H-Si [21-23] but they were essentially limited to the structural stability of defects. Here, we present a comprehensive study carried out via ab initio DFT simulations on the effect of p- and n-type substitutional defects on the structural and electronic properties of 2H-Si. We first provide a detailed analysis on how a given impurity can assume a different local symmetry depending on the host crystal phase. Then, by studying neutral and charged dopants, we carefully estimate donors and acceptors thermodynamic transition levels in 2H-Si and compare them with the 3C-Si case. Finally, by employing a simple model based on the 2H/3C band offset diagram, we show how the carrier concentration can be used to stabilize one phase with respect to the other.

## II. METHODOLOGY

All simulations were performed by spin-polarized ab initio DFT as implemented in the SIESTA code [24], using the generalized gradient approximation (GGA) with the Perdew, Burke, Ernzerhof (PBE) functional [25]. Since the diameter of the 2H nanowires grown experimentally is usually larger than 100 nm, we considered bulk supercells as an accurate and reliable approximation for such systems. Indeed at such large diameters, quantum confinement effect on the band gap can be neglected [26-28]. Norm-conserving pseudopotentials of Trouillier-Martins type were chosen to replace the core electrons. Valence electrons were instead represented using a doubleζ (DZ) polarized basis set for Si, and doubleζ plus two polarization orbitals (DZP) basis sets for impurities. We considered p-type (B, Al, Ga, In) and n-type (P, Sb, As) as common substitutional dopants in 2H-Si. Lattice parameters and band gap of pristine bulk 3C- and 2H-Si were calculated by sampling the Brillouin zone using a 8 × × 8 8 Monkhorst-Pack k-points grid [29]. Table I presents the converged results of the lattice parameter and the band gap calculations for both the phases as well as their experimental values. As expected, when compared to experiments, it is clear that both quantities are affected by the error associated to the GGA exchange-correlation functional. We used a denser grid of 25 × 25 × 25 k-points to calculate the dielectric constant of both phases. This quantity was needed for the charge correction scheme described below [30]. The values of the in-plane and out-of-plane components, reported in Table I, are in perfect agreement with other DFT calculations [31, 32] and only slightly overestimate the experimental values [33]. For the doped systems, a supercell of 512 (a 4 × × 4 4 replica of the 8-atoms conventional supercell with orthogonal cell vectors) and of 500 atoms (a 5 × × 5 5 replica of the primitive cell) were used for the 3C and 2H phases, respectively. This choice has guaranteed the convergence of the formation energy value for both neutral and charged impurities as it eliminated the spurious interaction between defect periodical images (see Supplemental Materials (SM) [34]). Such doped supercells correspond to a dopant concentration between 10 19 cm -3 and 10 20 cm -3 , a value that is one order of magnitude larger with respect to the impurity concentration usually measured in 2H-samples [3]. A2 × × 2 2 k-grid and a real space grid cut-off of 500 Ry were considered to ensure the proper convergence of the results. Self-consistency field cycles were converged with a tolerance of 10 -5 eV. Atomic positions were optimized until the forces were less than 10 -2 eV/Å and the stress lower than 10 -1 GPa. Relaxation was performed adopting a conjugate gradient algorithm that calculates DFTforces through the Hellman-Feynman theorem. Charged structures were relaxed starting from

the neutral defect optimized structure. A compensating jellium background was added to avoid divergences due to the defect-defect long-range electrostatic interaction between charged periodic cells.

For all the considered impurities, we calculated their formation energy, Eform , in a given crystal phase (2H or 3C). This quantity is defined as the energy required to replace a host atom with a neutral or charged defect. In the Zhang and Northrup formalism [35], it reads as follows:

$$E _ { f o r m } = E _ { t o t a l } ^ { D } - E _ { t o t a l } ^ { p r i s t } + \mu _ { 0 } - \mu _ { D } + q ( \mathfrak { e } _ { V B M } + \mu _ { e } ) + E _ { c o r r },$$

where E D total is the total energy of the doped system, E prist total is the total energy of the pristine system, µ 0 is the chemical potential of Si in a given phase while µ D is that of the defect. µ 0 is here defined as the total ground energy per atom of a pristine Si cell whereas µ D is calculated as the total energy of the isolated impurity atom. However, it is worth to note that there is no univocal criteria to calculate the chemical potential of an impurity as this quantity strongly depends on the synthesis conditions. As a consequence, the calculated value of the formation energy can assume a wide range of values depending on the way µ D is defined. Nevertheless, since we are interested in the relative stability of a dopant among the two phases, the chemical potential µ D erases if we refer to the difference between formation energies, ∆ E f orm 2 H -3 C = E H 2 f orm -E C 3 f orm rather than in Eform value in itself (this also holds because the difference of µ 0 among the two phases can be neglected). Yet, referring back to Eq. (1), the term q ( ε VBM + µ e ) is the charge multiplied by the electron chemical potential, µ e , using as energy reference the valence band maximum, ε VBM . µ e varies from 0 to the electronic band gap. Finally, the term Ecorr corresponds to the Freysoldt-Neugebauer-Van de Walle's finite-size correction [30] that cancels artificial interactions between charged defects and that has been proven to give very accurate results for different type of defects [36]. As other proposed schemes [36-39], this approach implies the solution of the Poisson equation for a model system and the alignment of the electrostatic potentials to derive the formation energy correction. In particular, solving the Poisson equation for the periodic model potential needs the calculation of the dielectric tensor profile of the material [39], whose static components are presented in Table I.

By calculating formation energies, we can also estimate the thermodynamic transition level between two charge states, ε ( q 1/ q 2), which corresponds to the specific value of the electron chemical potential for which the two states have the same value of Eform . This quantity can be calculated

as:

$$\varepsilon$$

$$& \varepsilon ( q _ { 1 } / q _ { 2 } ) = \frac { E _ { t o t a l } ^ { D } ( q _ { 1 } ) + E _ { c o r r } ( q _ { 1 } ) - E _ { t o t a l } ^ { D } ( q _ { 2 } ) - E _ { c o r r } ( q _ { 2 } ) } { q _ { 2 } - q _ { 1 } } \\ & - \varepsilon _ { V B M }.$$

For each dopant the formation energy plots are reported for different values of q (-1, 0 and +1) and they are linear functions of the chemical potential as shown in Eq. 1. The intersection of the formation energy lines for two different charge states q 1 and q 2 is the thermodynamic transition level calculated through Eq. 2

Transition levels ε ( q 1 / q 2 ) correspond to thermal ionization energies that can be measured by deep level transient spectroscopy (DLTS) [40]. We employed Eq. (2) to calculate both acceptors, ε A ( q 1 / q 2 ) and donor transition levels ε D q ( 1 / q 2 ) . It is worth to note that transition levels are not affected by the dopant chemical potential which assumes the same value for the charged and neutral impurity and hence cancels out in Eq. (2).

In the following analysis, Eq. (2), which depends on the valence band maximum of each phase taken separately, allows to compare the doping with a given impurity of the 2H and the 3C phase of Si. On the other hand, to give an accurate quantitative estimation of transition levels, in particular in the case of the phase coexistence [4, 5, 14], we align ε ( q 1/ q 2) to a common reference level, i.e. the vacuum level Φ (see SM for the details of the calculations [34]). Indeed, many works have demonstrated that, even if within a semi-local approach, energy levels of dopants generally remain preserved if they are referred to a common well defined potential [41-43]. This choice improves the quantitative comparison among the two phases. On the basis of this consideration, Eq. (2) can be modified as:

$$& \varepsilon ( q _ { 1 } / q _ { 2 } ) = \frac { E _ { t o t a l } ^ { D } ( q _ { 1 } ) + E _ { c o r r } ( q _ { 1 } ) - E _ { t o t a l } ^ { D } ( q _ { 2 } ) - E _ { c o r r } ( q _ { 2 } ) } { q _ { 2 } - q _ { 1 } } \\ & - ( \varepsilon _ { V B M } - \Phi ) = \overline { \varepsilon } ( q _ { 1 } / q _ { 2 } ) - \overline { \varepsilon } _ { V B M }.$$

where the charge transition level ε ( q 1 / q 2 ) and the ε VBM are referred to the vacuum potential.

## III. STRUCTURAL PROPERTIES

Westart our discussion by a detailed structural analysis of the relaxed doped structures. Table II and III report the first-neighbors bond lengths of group III and group V impurities in the cubic and hexagonal crystals. In the 3C-Si, both neutral donors and acceptors maintain the tetrahedral

symmetry ( Td ) forming four first-neighbor bonds of equal length. Variations in bond lengths compared with the pristine system are uniform and depend on the mismatch between the covalent radius of the impurity and that of the host element, as already reported by Amato et al. [21, 22]. The bond contractions and expansions discussed here agree with previous results. As for smaller defects, such as B and P, the first neighbors impurity bonds are compressed in respect to those of the pristine cell (see first block of Table II and III) while for all the other reported defects (Al, Ga, In, Sb and As in Table II and III) they tend to expand since their covalent radius is larger than those of Si. More in general, we can say that impurities in the 3C-Si phase induce an expansion (contraction) of their first neighbors bonds and, consequently, there is a noticeable contraction (expansion) between the first and second neighbors atomic layers that gradually fades within the crystal. The charged impurities also preserve the cubic Td symmetry. As for B -, the extra electron does not induce any appreciable change in the bond distances (see first block of Table II). Al -, Ga -, and As -first neighbors distances are only slightly reduced compared to the neutral case (see Table II). These tiny bond length variations are proportional to the impurity size, being more pronounced for larger dopants. Therefore, adding an extra electron to large radius acceptors (Al, Ga, In) in the 3C-Si host lattice induces a larger contraction compared to the neutral case, presumably because the lattice has fewer degrees of freedom to adjust to the extra electron configuration, resulting in a more significant contraction. The analysis of donors shows that, when positive charges are added, bonds tend to contract as well, but without any clear general trend (see Table III).

Compared to the 3C-Si case, the 2H-Si host lattice presents a lower symmetry due to the anisotropy along the [0001] hexagonal direction ( C v 3 symmetry). In the pristine structure, for a given Si atom, three equal length first-neighbor bonds in the c-(0001) plane can be defined. The fourth bond is longer and perpendicular to the plane, breaking the tetrahedral symmetry. This arrangement of the atoms increases the degrees of freedom for the structural relaxation around the defect. As shown in the first block of Table II, the local symmetry of the B 0 impurity is C v 3 : there exist three equivalent bonds in the plane and a longer one along the [0001] crystal direction. This is a consequence of three valence electrons of the acceptor trying to keep the three-fold coordination. The difference between the average of the three in-plane bonds and the out-of-plane one is larger than that of the 2H pristine case. As for B -, we observe the same symmetry but slightly more pronounced: an equivalence of the three bonds in the plane and a slight expansion of the

perpendicular one with respect to the neutral case. Looking at the second block of Table II, we can recognize a similar behavior in the case of Al: both Al 0 and Al -keep a clear C v 3 symmetry with three equal length bonds in the plane and a longer one along the c hexagonal direction. As for Ga 0 and In 0 (third and fourth blocks of Table II), though it is not very clear, it seems that the local symmetry is closer to the Td cubic one. Two bonds in the plane are larger than the third one, and even larger than the perpendicular interatomic distance, with this effect being more pronounced for In due to its larger atomic size. When a charge is added, the tendency toward Td symmetry becomes more evident, especially for In (fourth block of Table II). Possible factors that could explain this deviation from the expected bond length pattern include local distortions, impurity size mismatch as well as electronic wave function symmetry breaking. On the other hand, donors in 2H-Si induce similar contraction and expansion of the local structure, but with a different impurity symmetry. As reported in the first block of the Table III, P 0 has a C v 3 symmetry with a difference between the in-plane and the out-of-the-plane bond lengths larger than those of the pure 2H-Si. When a charge is added, the out-of-plane bond length is slightly reduced, while all the in-plane bonds maintain the same value. There is almost no change in the symmetry between neutral and charged cases, keeping a pronounced C v 3 symmetry. The same behavior can be recognized also for As and Sb (second and third block of Table III). It is worth to note that the difference between the in-plane and the out-of-plane bond length average is strictly related to the impurity ionic radius being lower for Sb and larger for P.

## IV. FORMATION ENERGY

The local structural analysis cannot provide alone a full picture of the impurity stability in different phases. Formation energy calculations of both neutral and charged dopants are required to analyze this question in depth. Figures 1 and 2 report the formation energy, Eform , calculated with Eq. (1), as a function of the chemical potential relative to the maximum of the valence band for p- and n-type dopants, respectively. Numerical values associated to the graphical plots of Figures 1 and 2 are provided in Table IV and Table V.

Looking at Figure 1, one can see, as expected, that negatively charged acceptors are more stable than neutral impurities because of their tendency to attract an electron from the host crystal to create a hole in the valence band. Interestingly, for both the neutral and charged cases, their

formation energy is always lower in the 2H phase with respect to the 3C one (see Figure 1). This is even clearer if one looks at the first column of Table IV where the difference in formation energy, ∆ E H 2 -3 C for a given impurity in the 3C and in the 2H phase is reported for both the charged and the neutral cases. The more ∆ E H 2 -3 C is negative, the more the impurity will prefer the 2H phase with respect to the 3C-one.

This behavior can be ascribed to the clear preference for acceptors to be in the 2H phase as a consequence of its stronger tendency to induce a local C v 3 symmetry and a three-fold coordination around the impurity (see Section III). The same explanation still holds when charged impurities are considered as this phase preference does not change (see Figure 1).

Similarly to the case of acceptors, Figure 2 shows that positively charged donors are more stable than neutral impurities because the additional donor electron, occupying an antibonding state, is only weakly bonded to the impurity and tends to be delocalized in the crystal. Furthermore, while neutral donors do not present any clear phase preference as ∆ E H 2 -3 C is close to zero, charged donors reveal a marked tendency to be more stable in the 3C-phase (see Figure 2 and the first two columns of Table V). So for charged donors the 3C-host lattice is the ideal crystal environment to form four equal length covalent bonds with a marked Td symmetry.

Our results clearly show that, close to the 2H/3C homo-junction [4, 5, 21], the ∆ E f orm 2 H -3 C of charged donors is positive while for charged acceptors it is negative. To align the Fermi levels of the two phases a space charge will be generated and thus the dopants will be necessarly charged at the interface. It is hence expected that a gradient in the dopant concentration at the interface between the two phases can be build up with donors accumulating on the 3C-Si side and acceptors in 2H-Si.

## V. THERMODYNAMIC TRANSITION LEVELS

The calculation of the formation energies of neutral and charged impurities allows to employ Eq. (2) to calculate thermodynamic transition levels, ε ( q 1 / q 2 ) , which are essential for the experimental characterization of the defect. These charge transition energies correspond to the crossing point where the charged and the neutral defect formation energies coincide, as observed in Figures 1 and 2. It is clear from the figures that, for acceptor dopants, the difference in the defects formation energy between the two phases changes in the order of a few meV with the state of charge of the impurities ( ∆ ε 2 H -3 C A also reported in the sixth column of Table IV). More negative

values of ∆ ε 2 H -3 C A correspond to a transition level that is shallower in the 2H phase. For B impurities the transition energy is positive, those of Ga and Al are practically zero (less than 2 meV), while that of In is negative. This means that by increasing the size of the defect, the difference between the transition energies between 2H and 3C will be more negative and hence the defect level will be shallower in 2H-Si. In the case of smaller radius impurities, ∆ ε 2 H -3 C A is positive and the defect level is shallower in 3C-Si. In summary, if one looks at the second and sixth column of Table IV, it is possible to conclude that charged acceptors are always more stable in the 2H phase and present charge transition levels shallower with respect to the 3C-Si. The magnitude of these energy shift depends from the ionic radius of the dopant.

A similar analysis can be performed for donors. The difference of the transition energies between the two phases is about hundreths of meV (see sixth column of Table V and Fig. 2) and are always negative. This means that the donor dopant levels will be always shallower in 3C-Si with respect to 2H-Si. As in the case of acceptors, even here we can recognize a relation with the ionic radius of the dopant. Indeed, ∆ ε 2 H -3 C D decreases as a function of the radius mismatch. These findings agree with the hydrogenic model of donors that considers four covalent bonds and an unpaired one whose energy relative to the bottom of the conduction band is directly proportional to the effective mass and inversely proportional to the dielectric function. As reported in Table I, the in-plane component of the dielectric function is lower for 2H- than for 3C-Si. This means that the transition level of donors in 2H-Si will be deeper with respect to those of 3C-Si. The discussed behavior, with the appropriate differences, was already observed in polytypes of GaAs, GaP, and InP [44].

Nevertheless a note of caution is needed when looking at the transition levels shown in Figure 1 and Figure 2. Indeed, some of these levels are resonant with the crystal host band states. This is directly related to the use of total energies obtained from semi-local functionals which cannot provide accurate ionization energies for shallow dopants [45, 46] but can only qualitatively reproduce the experimental chemical trend among different type of dopants, as clearly shown in the last column of Table IV and Table V. More advanced methods, such as quasi-particle G0W0 calculations, are known to provide accurate results but they are associated with a very high computational cost. Although in our case the absolute values of the transition energies for each phase taken separately may not be correct, the calculated semi-quantitative variation between the two phases is still valid and insightful. To further improve the comparison of the transition energy levels between the two phases, we calculated them with respect to a common reference point, i.e. the vacuum level, Φ ,

calculated from a Si slab model [47, 48] as described by Eq. (3) (see SM for more details about the procedure [34]). Indeed, even when the theoretical description is improved with beyond-DFT methods, band edges in semiconductors undergo significant shifts but charge transition levels of defects remain practically unaffected with respect to their semi-local energy positions [41-43]. Thus the calculated ∆ ε 2 H -3 C A and ∆ ε 2 H -3 C D presented in Table IV and Table V and graphically shown in Fig. 3 are a valid reference for future experiments aiming to measure the variation in the ionization energy among the two phases.

In agreement with other studies [8, 13, 15], we found that the 2H/3C Si junction presents a type-II alignment, as shown by the difference between hexagonal and cubic valence band maximum (VBM), ∆ VBM H 2 -3 C ≈ 0 154 eV and between hexagonal and cubic conduction band mini-. mum (CBM), ∆ CBM H 2 -3 C ≈ 0 036 eV. More accurate calculations [13] with the Heyd-Scuseria-. Ernzerhof hybrid functional (HSE06) [49] and the modified Becke-Johnson potential combined with local-density approximation correlation (mBJLDA) [50] functional report a valence band offset of ∆ VBM H 2 -3 C mBJLDA = 0 26 eV and . ∆ VBM H 2 -3 C HSE = 0 31 eV. Calculations . for the hexagonal/cubic homojunction supercell using the Local Density Approximation (LDA) with the GW correction [14] estimate a value of ∆ VBM H 2 -3 C GW = 0 25 eV, slightly larger than the one we ob-. tained. On the other hand, in the same work [14] a value of ∆ CBM H 2 -3 C GW = 0 20 eV was re-. ported, while Keller et al. [13] estimate conduction band offsets to be ∆ CBM H 2 -3 C MBJLDA = 0 07 eV . and ∆ CBM H 2 -3 C HSE = 0 04 eV, which is closer to our GGA results. . It is possible to conclude that, despite the use of different computational approaches and alignment procedures, the type-II character of the band offset is observed in various studies [8, 13-15]. The 3C-VBM lies well below the 2H-VBM, with differences on the order of hundreds of meV, similar to the magnitude of the calculated | ∆ E f orm q 2 H -3 C , | . When comparing the two phases, acceptor transition levels of both materials (referred to the vacuum level) are energetically lower in 3C-Si compared to those in 2H-Si. For donors, the band offset will lead to a closer alignment of the energy levels, although their exact relative positions may slightly differ slightly depending on the computational approach used.

## VI. BAND-OFFSET MODEL OF DOPING

In order to further confirm our DFT findings and to depict a full physical picture of the dopant stability in 2H and 3C-Si, we build up a simple model following the arguments employed by Dalpian et al. [17, 18] to explain the stability of III-V polytypes and recently adopted by Cai et

al. [20] to study 2H-Ge. Ideal doping can be represented by simply adding electrons and holes to the pristine supercell of 2H- and 3C-Si taking into account total energy correction due to finite size effects [30]. Figure 4 reports the difference in the total energy between the two phases, ∆ E total 2 H -3 C per atom as a function of the carrier concentration. The results show that introducing electrons tends to stabilize the 3C phase while introducing holes is a way to make the 2H phase more stable. These results are a consequence of the type II band offset between the two phases: since 2H-Si has a higher valence band maximum than 3C (see Fig. 3), it will cost less energy to create holes in the 2H-phase than in the 3C-one. On the other hand, since 3C-Si has a lower conduction band minimum than 2H-Si (see Fig. 3), it will cost less energy to add electrons in the 3C-Si than in 2HSi. Moreover, as the 2H/3C conduction band offset is very tiny, the slope of the electron doping curve (red line of Fig. 4) is found to be small while that of the hole incorporation is larger (green line of Fig. 4). Interestingly, this fully confirms our DFT results for impurity doped supercells which show that i) acceptors are more stable in the 2H phase and have shallower transition energy levels if compared to 3C-Si, ii) donors prefer instead to occupy cubic crystal site with shallower transition energies if compared to 2H-Si.

## VII. CONCLUSION

In conclusion, this work presents DFT calculations of the structural and electronic properties of p- and n-type impurities in 2H-Si, comparing them with the 3C-Si case. We find that acceptors tend to be more stable in the hexagonal phase because of its stronger tendency to induce a local C v 3 symmetry and three-fold coordination around the impurity. This behaviour is valid for all the considered acceptors (B, Al, Ga and In). In particular, charged acceptors are always more stable in the 2H phase and present a shallower charge transition level with respect to the 3C-Si crystal host.

On the other hand, we find that neutral donors (P, As, Sb) have no phase preference while positively charged donors have lower formation energies in the cubic phase than in the hexagonal one. This can be ascribed to the fact that for charged donors the 3C-host lattice is the ideal crystal environment to form four equal length covalent bonds with a marked Td symmetry. Furthermore, donors levels are shallowers in 3C-Si than in 2H-Si.

The physical validity of the doped supercell calculations is reinforced by a simple model based on the 2H/3C band offset diagram through which we show how holes can be used to stabilize the

2H-Si phase while electrons prefer the 3C-Si crystal phase.

All these results assume a particular importance for more accurate calculations beyond the DFT framework and for future experiments. On the one hand, these results could be further assessed by hybrid-DFT or GW based calculations that could confirm the difference in transition level positions among 2H and 3C-Si by providing numerical values to be directly compared with experiments. On the other hand, being a theoretical estimation of transition energy levels in 2H-Si, these findings represent a valid theoretical reference for future measurements of the variation in the ionization energy of a given dopant among the two phases. Moreover, in the case of 2H/3C homojunctions [4, 5, 21], because of the different stability of charged dopants, a gradient in dopant concentration at the interface between the two phases may be experimentally detected with donors accumulating on the 3C-Si side and acceptors in 2H-Si.

## ACKNOWLEDGMENTS

All the authors acknowledge the ANR AMPHORE project (ANR-21-CE09-0007) and the ANR TULIP (ANR-24-CE09-5076). Part of the high-performance computing resources for this project were granted by the Institut du développement et des ressources en informatique scientifique (IDRIS) under the allocations AD010914974 and AD010915077 via GENCI (Grand Equipement National de Calcul Intensif). The authors also acknowledge the use of the CERES high performance computer cluster at the Laboratoire de Physique des Solides of the Université Paris-Saclay (Orsay, France).

- [1] L. S. Ramsdell, Studies on silicon carbide, Am. Mineral. 32 (1947).
- [2] L. H. Galvão Tizei and M. Amato, Electronic structure and optical properties of semiconductor nanowires polytypes, Eur. Phys. J. B 93 , 1 (2020).
- [3] E. M. Fadaly, A. Dijkstra, J. R. Suckert, D. Ziss, M. A. Van Tilburg, C. Mao, Y. Ren, V. T. van Lange, K. Korzun, S. Kölling, et al. , Direct-bandgap emission from hexagonal Ge and SiGe alloys, Nature 580 , 205 (2020).
- [4] L. Vincent, G. Patriarche, G. Hallais, C. Renard, C. Gardès, D. Troadec, and D. Bouchier, Novel heterostructured Ge nanowires based on polytype transformation, Nano Lett. 14 , 4828 (2014).
- [5] J. Tang, J.-L. Maurice, F. Fossard, I. Florea, W. Chen, E. Johnson, M. Foldyna, L. Yu, and P. R.

i Cabarrocas, Natural occurrence of the diamond hexagonal structure in silicon nanowires grown by a plasma-assisted vapour-liquid-solid method, Nanoscale 9 , 8113 (2017).

- [6] H. I. T. Hauge, M. A. Verheijen, S. Conesa-Boj, T. Etzelstorfer, M. Watzinger, D. Kriegner, I. Zardo, C. Fasolato, F. Capitani, P. Postorino, et al. , Hexagonal silicon realized, Nano Lett. 15 , 5855 (2015).
- [7] C. Rödl, J. Furthmüller, J. R. Suckert, V. Armuzza, F. Bechstedt, and S. Botti, Accurate electronic and optical properties of hexagonal germanium for optoelectronic applications, Phys. Rev. Mater. 3 , 034602 (2019).
- [8] T. Kaewmaraya, L. Vincent, and M. Amato, Accurate estimation of band offsets in group iv polytype junctions: A first-principles study, J. Phys. Chem. C 121 , 5820 (2017).
- [9] H. I. T. Hauge, S. Conesa-Boj, M. A. Verheijen, S. Koelling, and E. P. Bakkers, Single-crystalline hexagonal silicon-germanium, Nano Lett. 17 , 85 (2017).
- [10] W. H. Peeters, V. T. van Lange, A. Belabbes, M. C. van Hemert, M. M. Jansen, R. Farina, M. A. van Tilburg, M. A. Verheijen, S. Botti, F. Bechstedt, et al. , Direct bandgap quantum wells in hexagonal silicon germanium, Nat. Comm. 15 , 5252 (2024).
- [11] P. Borlido, J. R. Suckert, J. Furthmüller, F. Bechstedt, S. Botti, and C. Rödl, From pseudo-direct hexagonal germanium to direct silicon-germanium alloys, Phys. Rev. Mater. 5 , 114604 (2021).
- [12] C. Rödl, T. Sander, F. Bechstedt, J. Vidal, P. Olsson, S. Laribi, and J.-F. Guillemoles, Wurtzite silicon as a potential absorber in photovoltaics: Tailoring the optical absorption by applying strain, Phys. Rev. B 92 , 045207 (2015).
- [13] M. Keller, A. Belabbes, J. Furthmüller, F. Bechstedt, and S. Botti, First-principles insight in structureproperty relationships of hexagonal Si and Ge polytypes, Phys. Rev. Mater. 7 (2023).
- [14] M. Amato, T. Kaewmaraya, A. Zobelli, M. Palummo, and R. Rurali, Crystal phase effects in si nanowire polytypes and their homojunctions, Nano Lett. 16 , 5694 (2016).
- [15] A. Belabbes, S. Botti, and F. Bechstedt, Band lineup at hexagonal si x ge 1- x/si y ge 1- y alloy interfaces, Phys. Rev. B 106 , 085303 (2022).
- [16] V. Heine, C. Cheng, and R. J. Needs, The preference of silicon carbide for growth in the metastable cubic form, J. Am. Cer. Soc. 74 , 2630 (1991).
- [17] G. M. Dalpian, Y. Yan, and S. H. Wei, Impurity-induced phase stabilization of semiconductors, Appl. Phys. Lett. 89 (2006).
- [18] G. M. Dalpian and S.-H. Wei, Hole-mediated stabilization of cubic gan, Phys. Rev. Lett. 93 , 216401 (2004).

- [19] A. Liu, Z. Xi, M. Li, J. Yang, L. Qi, and R. Goldman, Influence of excess silicon on polytype selection during metal-mediated epitaxy of gan nanowires, Appl. Phys. Lett. 125 (2024).
- [20] X. Cai, H. X. Deng, and S. H. Wei, Carrier-stabilized hexagonal ge, Phys. Rev. B 103 , 245202 (2021).
- [21] M. Amato, S. Ossicini, E. Canadell, and R. Rurali, Preferential positioning, stability, and segregation of dopants in hexagonal si nanowires, Nano Lett. 19 , 866 (2019).
- [22] M. Amato, T. Kaewmaraya, and A. Zobelli, Extrinsic Doping in Group IV Hexagonal-Diamond-Type Crystals, J. Phys. Chem. C 124 , 17290 (2020).
- [23] L. Sun, M. R. G. Marques, M. A. L. Marques, and S. Botti, Point defects in hexagonal silicon, Phys. Rev. Mater. 5 , 64605 (2021).
- [24] J. M. Soler, E. Artacho, J. D. Gale, A. García, J. Junquera, P. Ordejón, and D. Sánchez-Portal, The siesta method for ab initio order-n materials simulation, J. Condens. Matter Phys. 14 , 2745 (2002).
- [25] J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77 , 3865 (1996).
- [26] M. Amato, M. Palummo, and S. Ossicini, Band structure analysis in sige nanowires, Mater. Sci. Eng.: B 177 , 705 (2012).
- [27] M. Amato, M. Palummo, and S. Ossicini, Segregation, quantum confinement effect and band offset for [110] sige nws, Phys. Status Solidi B 247 , 2096 (2010).
- [28] N. S. Mohammad, Understanding quantum confinement in nanowires: basics, applications and possible laws, J. Phys.: Cond. Matter 26 , 423202 (2014).
- [29] H. J. Monkhorst and J. D. Pack, Special points for brillouin-zone integrations, Phys. Rev. B 13 , 5188 (1976).
- [30] C. Freysoldt, J. Neugebauer, and C. G. Van De Walle, Fully Ab initio finite-size corrections for charged-defect supercell calculations, Phys. Rev. Lett. 102 , 016402 (2009).
- [31] S. Baroni and R. Resta, Ab initio calculation of the macroscopic dielectric constant in silicon, Phys. Rev. B 33 , 7017 (1986).
- [32] F. Zheng, J. Tao, and A. M. Rappe, Phys. Rev. B 95 (2017).
- [33] H. H. Li, Refractive index of silicon and germanium and its wavelength and temperature derivatives, J. Phys. Chem. Ref. Data 9 , 561 (1980).
- [34] See Supplemental Material at URL-will-be-inserted-by-publisher for the details of the vacuum levels calculations as well as the convergence of the dopant formation energy with the supercell size. Supplemental Material also contains Refs. [47,48].

- [35] S. B. Zhang and J. E. Northrup, Chemical potential dependence of defect formation energies in GaAs: Application to Ga self-diffusion, Phys. Rev. Lett. 67 , 2339 (1991).
- [36] C. Freysoldt, B. Grabowski, T. Hickel, J. Neugebauer, G. Kresse, A. Janotti, and C. G. Van de Walle, First-principles calculations for point defects in solids, Rev. Mod. Phys. 86 , 253 (2014).
- [37] I. Dabo, B. Kozinsky, N. E. Singh-Miller, and N. Marzari, Electrostatics in periodic boundary conditions and real-space corrections, Phys. Rev. B 77 , 115139 (2008).
- [38] Y. Kumagai and F. Oba, Electrostatics-based finite-size corrections for first-principles point defect calculations, Phys. Rev. B 89 , 195205 (2014).
- [39] M. H. Naik and M. Jain, Coffee: corrections for formation energy and eigenvalues for charged defect simulations, Comp. Phys. Comm. 226 , 114 (2018).
- [40] R. K. Willardson, E. R. Weber, and M. Stavola, Identification of Defects in Semiconductors (Academic Press, 1998).
- [41] A. Alkauskas and A. Pasquarello, Band-edge problem in the theoretical determination of defect energy levels: The o vacancy in zno as a benchmark case, Phys. Rev. B 84 , 125206 (2011).
- [42] H.-P. Komsa, P. Broqvist, and A. Pasquarello, Alignment of defect levels and band edges through hybrid functionals: Effect of screening in the exchange term, Phys. Rev. B 81 , 205118 (2010).
- [43] A. Alkauskas, P. Broqvist, and A. Pasquarello, Defect energy levels in density functional calculations: Alignment and band gap problem, Physical Rev. Lett. 101 , 046405 (2008).
- [44] G. Giorgi, M. Amato, S. Ossicini, X. Cartoixà, E. Canadell, and R. Rurali, Doping of iii-v arsenide and phosphide wurtzite semiconductors, J. Phys. Chem. C (2020).
- [45] L.-W. Wang, Density functional calculations of shallow acceptor levels in si, J. Appl. Phys. 105 (2009).
- [46] G. Zhang, A. Canning, N. Grønbech-Jensen, S. Derenzo, and L.-W. Wang, Shallow impurity level calculations in semiconductors using ab initio methods, Phys. Rev. Lett. 110 , 166404 (2013).
- [47] J. Junquera, M. Zimmer, P. Ordejón, and P. Ghosez, First-principles calculation of the band offset at BaO/BaTiO3 and SrO/SrTiO3 interfaces, Phys. Rev. B 67 , 155327 (2003).
- [48] M. Bertocchi, M. Amato, I. Marri, and S. Ossicini, Phys. Status Solidi C 14 , 1700193 (2017).
- [49] J. Heyd, G. E. Scuseria, and M. Ernzerhof, Hybrid functionals based on a screened Coulomb potential, The Journal of Chemical Physics 118 , 8207 (2003).
- [50] F. Tran and P. Blaha, Accurate Band Gaps of Semiconductors and Insulators with a Semilocal Exchange-Correlation Potential, Phys. Rev. Lett. 102 , 226401 (2009).
- [51] B. J. Isherwood and C. A. Wallace, Measurement of the Lattice Parameter of Silicon, using a Double-

diffraction Effect, Nature 212 , 173 (1966).

- [52] K. Godwod, R. Kowalczyk, and Z. Szmid, Application of a precise double X-ray spectrometer for accurate lattice parameter determination, Phys. Status Solidi (a) 21 , 227 (1974).
- [53] J. W. Precker and M. A. da Silva, Experimental estimation of the band gap in silicon and germanium from the temperature-voltage curve of diode thermometers, Am. J. Phys. 70 , 1150 (2002).
- [54] J. Low, M. Kreider, D. Pulsifer, A. Jones, and T. Gilani, Band Gap Energy in Silicon, Am. J. Undergrad. Res. 7 (2008).
- [55] H. S. Ahn, S. W. Kim, G. S. Lee, K. H. Kim, J. H. Lee, D. H. Ha, Y. T. Chun, and S. Ryu, Optical property of hexagonal (2H) silicon crystal, Semicond. Sci. Technol. 36 , 095023 (2021).
- [56] Y. Zhang, Z. Iqbal, S. Vijayalakshmi, and H. Grebel, Stable hexagonal-wurtzite silicon phase by laser ablation, Appl. Phys. Lett. 75 , 2758 (1999).
- [57] J. M. Besson, E. H. Mokhtari, J. Gonzalez, and G. Weill, Electrical properties of semimetallic silicon III and semiconductive silicon IV at ambient pressure, Phys. Rev. Lett. 59 , 473 (1987).
- [58] F. Fabbri, E. Rotunno, L. Lazzarini, N. Fukata, and G. Salviati, Visible and Infra-red Light Emission in Boron-Doped Wurtzite Silicon Nanowires, Sci. Rep. 4 , 3603 (2014).
- [59] E. Burstein, G. Picus, B. Henvis, and R. Wallis, Absorption spectra of impurities in silicon-I: GroupIII acceptors, J. Phys. Chem. Solids 1 , 65 (1956).
- [60] M. Steger, A. Yang, D. Karaiskaj, M. L. W. Thewalt, E. E. Haller, J. W. Ager, M. Cardona, H. Riemann, N. V. Abrosimov, A. V. Gusev, A. D. Bulanov, A. K. Kaliteevskii, O. N. Godisov, P. Becker, and H.J. Pohl, Shallow impurity absorption spectroscopy in isotopically enriched silicon, Phys. Rev. B. 79 , 205210 (2009).
- [61] D. Karaiskaj, J. A. H. Stotz, T. Meyer, M. L. W. Thewalt, and M. Cardona, Impurity Absorption Spectroscopy in 28 Si: The Importance of Inhomogeneous Isotope Broadening, Phys. Rev. Lett. 90 , 186402 (2003).
- [62] C. Xiao, D. Yang, X. Yu, L. Xiang, and D. Que, Determination of the Boron and Phosphorus Ionization Energies in Compensated Silicon by Temperature-Dependent Luminescence, Silicon 9 , 147 (2017).
- [63] N. Q. Vinh, B. Redlich, A. F. G. van der Meer, C. R. Pidgeon, P. T. Greenland, S. A. Lynch, G. Aeppli, and B. N. Murdin, Time-Resolved Dynamics of Shallow Acceptor Transitions in Silicon, Phys. Rev. X. 3 , 11019 (2013).
- [64] A. Fargi, N. Hizem, A. Kalboussi, and A. Souifi, Electrical Analysis of Indium Deep Levels Effects on Kink Phenomena of Silicon NMOSFETs, World J. Nano Sci. Eng. 4 , 9 (2014).

- [65] G. J. Parker, S. D. Brotherton, I. Gale, and A. Gill, Measurement of concentration and photoionization cross section of indium in silicon, J. Appl. Phys. 54 , 3926 (1983).
- [66] G. F. Cerofolini, G. U. Pignatel, E. Mazzega, and G. Ottaviani, Supershallow levels in indium-doped silicon, J. Appl. Phys. 58 , 2204 (1985).
- [67] C. E. Jones and G. E. Johnson, Deep level transient spectroscopy studies of trapping parameters for centers in indium-doped silicon, J. Appl. Phys. 52 , 5159 (1981).
- [68] G. Picus, E. Burstein, and B. Henvis, Absorption spectra of impurities in silicon-II: Group-V donors, J. Phys. Chem. Solids 1 , 75 (1956).
- [69] H. Jiang and Y.-C. Shen, Ionization potentials of semiconductors from first-principles, J. Chem. Phys. 139 (2013).
- [70] I. Marri, M. Amato, M. Bertocchi, A. Ferretti, D. Varsano, and S. Ossicini, Phys. Chem. Chem. Phys. 22 , 25593 (2020).

TABLE I. Calculated GGA lattice parameter, bandgap energy, and static dielectric constants components for 3C- and 2H-Si as well as the corresponding experimental values.

|       | a (Å)               | c (Å)               | E GAP (eV)          | ε xx            | ε zz            |
|-------|---------------------|---------------------|---------------------|-----------------|-----------------|
| 3C-Si | 5.45 (5.43 a )      | /                   | 0.58 (1.12-1.23 b ) | 13.07 (11.7 c ) | 13.07 (11.7 c ) |
| 2H-Si | 3.84 (3.82-3.84 d ) | 6.34 (6.26-6.34 d ) | 0.46 (0.92 e )      | 12.89           | 13.90           |

a From Ref. [51] and [52].

b From Ref. [53] and [54].

c From Ref. [33].

d From Ref. [6], [55], [56], and [57]. Values summarized in Ref. [13].

e From Ref. [58].

TABLE II. Impurity first neighbor's bond distances for acceptors in the 3C-Si and 2H-Si host lattices after relaxation. For the 2H phase d 1, d 2, and d 3 correspond to distances calculated in the c-(0001) plane while d 4 lies along the [0001] hexagonal direction. In parenthesis the contraction (negative) or expansion (positive) with respect to the pristine Si distances are given.

|    | d 1 (Å)   | d 2 (Å)   | d 3 (Å)   | d 4 (Å)   |
|----|-----------|-----------|-----------|-----------|
| 3C | 2 . 079   | 2 . 079   | 2 . 079   | 2 . 079   |
|    | (-11.9%)  | (-11.9%)  | (-11.9%)  | (-11.9%)  |
| 2H | 2 . 076   | 2 . 076   | 2 . 076   | 2 . 092   |
|    | (-11.9%)  | (-11.9%)  | (-11.9%)  | (-11.8%)  |
| 3C | 2 . 079   | 2 . 079   | 2 . 079   | 2 . 079   |
|    | (-11.9%)  | (-11.9%)  | (-11.9%)  | (-11.9%)  |
| 2H | 2 . 073   | 2 . 073   | 2 . 073   | 2 . 095   |
|    | (-12%)    | (-12%)    | (-12%)    | (-11.6%)  |
| 3C | 2 . 427   | 2 . 427   | 2 . 427   | 2 . 427   |
|    | (2.9%)    | (2.9%)    | (2.9%)    | (2.9%)    |
| 2H | 2 . 420   | 2 . 427   | 2 . 427   | 2 . 430   |
|    | (2.8%)    | (3.0%)    | (3.0%)    | (2.5%)    |
| 3C | 2 . 426   | 2 . 426   | 2 . 426   | 2 . 426   |
|    | (2.8%)    | (2.8%)    | (2.8%)    | (2.8%)    |
|    | 2 . 421   | 2 . 421   | 2 . 421   | 2 . 433   |
| 2H | (2.8%)    | (2.8%)    | (2.8%)    | (2.6%)    |
|    | 2 . 470   | 2 . 470   | 2 . 470   | 2 . 470   |
| 3C | (4.7%)    | (4.7%)    | (4.7%)    | (4.7%)    |
|    | 2 . 461   | 2 . 473   | 2 . 473   | 2 . 463   |
| 2H | (4.5%)    | (5.0%)    | (5.0%)    | (3.9%)    |
| 3C | 2 . 464   | 2 . 464   | 2 . 464   | 2 . 464   |
|    | (4.4%)    | (4.4%)    | (4.4%)    | (4.4%)    |
|    | 2 . 462   | 2 . 463   | 2 . 463   | 2 . 466   |
| 2H | (4.6%)    | (4.6%)    | (4.6%)    | (4.0%)    |
| 3C | 2 . 581   | 2 . 581   | 2 . 581   | 2 . 581   |
|    | (9.4%)    | (9.4%)    | (9.4%)    | (9.4%)    |
|    | 2 . 572   | 2 . 586   | 2 . 586   | 2 . 569   |
| 2H | (9.2%)    | (9.8%)    | (9.8%)    | (8.4%)    |
| 3C | 2 . 571   | 2 . 571   | 2 . 571   | 2 . 571   |
|    | (9.0%)    | (9.0%)    | (9.0%)    | (9.0%)    |
|    | 2 . 571   | 2 . 573   | 2 . 573   | 2 . 573   |
| 2H | (9.2%)    | (9.2%)    | (9.2%)    | (8.5%)    |

TABLE III. Impurity first neighbor's bond distances for donors in the 3C-Si and 2H-Si host lattices after relaxation. For the 2H phase d 1, d 2, and d 3 correspond to distances calculated in the c-(0001) plane while d 4 lies along the [0001] hexagonal direction. In parenthesis the contraction (negative) or expansion (positive) with respect to the pristine Si distances are given.

|    | d 1 (Å)   | d 2 (Å)   | d 3 (Å)   | d 4 (Å)   |
|----|-----------|-----------|-----------|-----------|
| 3C | 2 . 36    | 2 . 36    | 2 . 36    | 2 . 36    |
|    | (-0.1 %)  | (-0.1 %)  | (-0.1 %)  | (-0.1 %)  |
| 2H | 2 . 35    | 2 . 35    | 2 . 35    | 2 . 47    |
|    | (-0.2 %)  | (-0.2 %)  | (-0.2 %)  | (-0.1 %)  |
| 3C | 2 . 36    | 2 . 36    | 2 . 36    | 2 . 36    |
|    | (-0.1 %)  | (-0.1 %)  | (-0.1 %)  | (-0.1 %)  |
| 2H | 2 . 35    | 2 . 35    | 2 . 35    | 2 . 37    |
|    | (-0.2 %)  | (-0.2 %)  | (-0.2 %)  | (0.0 %)   |
| 3C | 2 . 47    | 2 . 47    | 2 . 47    | 2 . 47    |
|    | (4.7 %)   | (4.7 %)   | (4.7 %)   | (4.7 %)   |
|    | 2 . 44    | 2 . 44    | 2 . 44    | 2 . 45    |
| 2H | (5.0 %)   | (5.0 %)   | (5.0 %)   | (3.9 %)   |
| 3C | 2 . 46    | 2 . 46    | 2 . 46    | 2 . 46    |
|    | (4.4 %)   | (4.4 %)   | (4.4 %)   | (4.4 %)   |
|    | 2 . 43    | 2 . 43    | 2 . 43    | 2 . 45    |
| 2H | (4.6 %)   | (4.6 %)   | (4.6 %)   | (4.0 %)   |
| 3C | 2 . 57    | 2 . 57    | 2 . 57    | 2 . 57    |
|    | (9.1 %)   | (9.1 %)   | (9.1 %)   | (9.1 %)   |
|    | 2 . 57    | 2 . 57    | 2 . 57    | 2 . 58    |
| 2H | (9.1 %)   | (9.1 %)   | (9.1 %)   | (8.9 %)   |
| 3C | 2 . 57    | 2 . 57    | 2 . 57    | 2 . 57    |
|    | (9.0 %)   | (9.0 %)   | (9.0 %)   | (9.0 %)   |
|    | 2 . 57    | 2 . 57    | 2 . 57    | 2 . 58    |
| 2H | (9.1 %)   | (9.1 %)   | (9.1 %)   | (8.9 %)   |

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00451v2-figure-4902ee12-a1ac-4541-8310-d9dbcb6f7ba5.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00451v2-figure-a4472957-6df4-4578-a946-59b3a3050b6e.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00451v2-figure-f9488941-7293-4cb3-bce3-5658bae3efc8.jpg)

FIG. 1. Formation energy diagram as a function of the chemical potential for acceptors a) In, b) Ga, c) Al, and d) B. For each dopant the plots are reported for two values of q (-1 and 0) and they are linear functions of the chemical potential as shown in Eq. 1. The chemical potential is set to 0 at the valence band maximum. The intersection of the formation energy lines for two charge states q = 0 and q = -1 is the thermodynamic transition level calculated through Eq. 2. The continuous line corresponds to the 3C phase and the dashed line to the 2H. The low energy dark pink region corresponds to the valence band.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00451v2-figure-e792b9ae-a22d-4c4b-aa39-0fc9ad2fc0e5.jpg)

TABLE IV. Acceptors: columns two and three report the difference in the formation energy between the 2H and 3C phases for neutral ( ∆ E 2H-3C 0 form , ) and charged impurities ( ∆ E 2H-3C -1 form , ). The third and fourth columns correspond to the transition energy levels for the 2H and 3C phases calculated Eq. (2) between the charge state q = 0 and q = -1. The fifth column reports the difference in acceptor transition energies, ∆ ε 2H-3C A , between the 2H and 3C phases. Negative (positive) values correspond to shallower (deeper) states in the hexagonal phase. The last column lists the experimental transition levels for the cubic phase.

|    |   ∆ E 2H-3C 0 , form (eV) |   ∆ E 2H-3C - 1 , form (eV) |   ε 2H A ( 0 / - 1 ) (meV) |   ε 3C A ( 0 / - 1 ) (meV) |   ∆ ε 2H-3C A (meV) | ε 3C, exp A ( 0 / - 1 ) (meV)   |
|----|---------------------------|-----------------------------|----------------------------|----------------------------|---------------------|---------------------------------|
| B  |                     -0.24 |                       -0.23 |                         -7 |                        -15 |                   8 | 30, 42-45 a                     |
| Al |                     -0.25 |                       -0.25 |                         16 |                         15 |                   0 | 55-57 b                         |
| Ga |                     -0.25 |                       -0.25 |                         33 |                         34 |                   0 | 58-65 c                         |
| In |                     -0.25 |                       -0.26 |                         52 |                         63 |                 -11 | 150-180 d                       |

a From Ref. [59-62].

b From Ref. [59, 63].

c From Ref. [59, 63].

d From Ref. [59, 64-67].

TABLE V. Donors: columns two and three report the difference in the formation energy between the 2H and 3C phases for neutral ( ∆ E 2H-3C 0 form , ) and charged impurities ( ∆ E 2H-3C + 1 form , ). The third and fourth columns correspond to the transition energy levels for the 2H and 3C phases calculated with Eq. (2) between the charge state q = +1 and q =0. The fifth column reports the difference in donors transition energies, ∆ ε 2H-3C D , between the 2H and 3C phases. Positive (negative) values correspond to shallower (deeper) states in the hexagonal phase. The last column lists the experimental transition levels for the cubic phase.

|    |   ∆ E 2H-3C 0 , form (eV) |   ∆ E 2H-3C + 1 , form (eV) |   ε 2H D (+ 1 / 0 ) (meV) |   ε 3C D (+ 1 / 0 ) (meV) |   ∆ ε 2H-3C D (meV) | ε 3C, exp D (+ 1 / 0 ) (meV)   | ε 3C, exp D (+ 1 / 0 ) (meV)   |
|----|---------------------------|-----------------------------|---------------------------|---------------------------|---------------------|--------------------------------|--------------------------------|
| P  |                      0.01 |                        0.27 |                       392 |                       647 |                -145 | 1076 -                         | 1078                           |
| As |                      0.02 |                        0.27 |                       394 |                       642 |                -138 | 1071                           | b                              |
| Sb |                      0.01 |                        0.26 |                       413 |                       658 |                -135 | 1081                           | c                              |

a From Ref. [60-63, 68].

b From Ref. [68].

c From Ref. [63, 68].

FIG. 2. Formation energy diagram as a function of the chemical potential for donors a) Sb, b) As, and c) P. For each dopant the plots are reported for two values of q (0 and +1) and they are linear functions of the chemical potential as shown in Eq. 1. The chemical potential is set to 0 at the valence band maximum. The intersection of the formation energy lines for two charge states q = +1 and q = 0 is the thermodynamic transition level calculated through Eq. 2. The continuous line corresponds to the 3C phase and the dashed line to the 2H. The first high energy pink region corresponds to the hexagonal conduction band while the high energy pink region (darker) to the 3C conduction band.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00451v2-figure-ca093e06-ceb0-4a1c-a433-20130ae07807.jpg)

FIG. 3. Transition energy levels of acceptors, ε A , (red lines) and donors, ε D , (blue lines) in 3C-Si (left) and 2H-Si (right) as calculated in Eq. (3). It is worth to note that ε A and ε D are evaluated referring to the vacuum potential Φ taken as zero of the energy. The DFT-GGA calculated value of ionization potential and electronic affinity for 3C-Si is in good agreement with previous theoretical works using the same approach [48, 69, 70].

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00451v2-figure-66b18390-123f-4e0c-8fed-3bc9ba1c43b8.jpg)

FIG. 4. Difference in the total energy, ∆ E total 2 H -3 C , per atom between 2H-Si and 3C-Si as a function of the carrier concentration per atom for electrons (red square symbols line) and hole (green diamond symbols line).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00451v2-figure-61f61f25-f7d3-4e4f-8201-082e22fd2dd1.jpg)
Marc Túnica Alberto Zobelli Michele Amato
Acceptor and donor impurity levels in hexagonal-diamond silicon
Recent advances in the characterization of hexagonal-diamond silicon (2H-Si)
have shown that this material possesses remarkably different structural,
electronic, and optical properties as compared to the common cubic-diamond (3C)
polytype. Interestingly, despite the wide range of physical properties
analyzed, to date no study has investigated impurity energy levels in 2H-Si.
Here, we present results of ab initio DFT simulations to describe the effect of
p- and n-type substitutional doping on the structural and electronic properties
of hexagonal-diamond Si (2H-Si). We first provide a detailed analysis of how a
given impurity can assume a different local symmetry depending on the host
crystal phase. Then, by studying neutral and charged dopants, we carefully
estimate donors and acceptors transition energy levels in 2H-Si and compare
them with the cubic-diamond (3C) case. In the case of acceptors, the formation
energy is always lower in 2H-Si and is associated with a shallower charge
transition level with respect to 3C-Si. On the other hand, donors prefer the
cubic phase and have transition energies smaller with respect to 2H-Si.
Finally, by employing a simple model based on the 2H/3C band offset diagram, we
prove the physical validity of our findings and we show how holes can be used
to stabilize the 2H-Si phase. Overall, the described doping properties
represent a robust starting point for further theoretical and experimental
investigations.