Corrigendum: Band structure of Ge 1 − x Sn x alloy: a full-zone 30-band k  ·  p model (2019 New J. Phys. 21 073037)

In this corrigendum, we present corrections to the paper ‘ Band structure of Ge - x 1 Sn x alloy: a full-zone 30-band · k p model ’ ( 2019 New J. Phys. 21 073037 ) by Zhigang Song et al . Two typos are found in the P 1 and Q 2 expressions in table 3, which have been corrected as following. Abstract A full-zone 30-band k · p model is developed as an ef ﬁ cient and reliable tool to compute electronic band structure in Ge 1 − x Sn x alloy. The model was ﬁ rst used to reproduce the electronic band structures in Ge and α -Sn obtained with empirical tight binding and ab initio methods. Input parameters for the 30-band k · p model are carefully calibrated against prior empirical predications and experimental data. Important material properties such as effective mass for electrons and holes, Luttinger parameters, and density of states are obtained for Ge 1 − x Sn x alloy with the composition range 0  <  x  <  0.3. The 30-band k  ·  p model that requires far less computing resources is a necessary capability for optimization of sophisticated devices made from Ge 1 − x Sn x alloy with a large parameter space to explore.


Introduction
The field of Si photonics has seen impressive growth since the early vision in 1990s [1]. The huge infrastructure of the global Si electronics industry is expected to benefit the fabrication of highly sophisticated Si photonic devices at costs that are lower than those currently required for compound semiconductors. Following the discovery of high speed SiGe electronics [2], the landscape of Si-based optical platform quickly expanded to include Ge and SiGe alloys. Among all the necessary components that make up the complete set of the Si photonics, efficient light sources such as LEDs and lasers are the most challenging because Si, Ge, and their alloys are all indirect bandgap materials. In order to improve the emission efficiency, effort has been directed towards developing direct bandgap material by incorporating yet another group-IV element, Sn, into the mix [3]. Early theoretical studies suggested that the group-IV alloy SiGeSn may possess a direct bandgap since Sn has negative bandgap at Γ-point [4]. But the progress in developing actual light sources has been painfully slow due to the grand technical challenges in growing high quality SiGeSn materials with sufficient Sn to turn the materials into direct bandgap. Si, with its band gap at Γ valley sitting 2.28 eV above its indirect band gap at the X-valley [5], naturally requires high Sn composition to pull its Γ-valley bandgap down to form direct-bandgap SiSn alloy. The difference for Ge, on the other hand, is smaller-only 0.14 eV between its indirect L-valley and direct Γ valley [6]. As a result, relatively small amount of Sn composition, estimated to be around 6%-8% [7], is sufficient for GeSn alloy to turn direct. Mid infrared (IR) electrical-injection lasers based on direct-bandgap GeSn as active region with latticed matched SiGeSn barrier have been predicted [8,9]. To date optically pumped mid-IR GeSn lasers have been demonstrated [10][11][12]. There is reason to be optimistic that the grand prize this material system has to offer, i.e. electrically pumped GeSn laser diodes, is finally within reach and the community will soon be facing the task of device performance optimization. This task is complicated by the fact that many material parameters remain unknown for a wide range of alloy compositions that need to be extracted from electronic band structures across the entire Brillion zone (BZ).
The electronic band structure of bulk Ge 1−x Sn x alloy can be calculated using empirical pseudopotential method [13][14][15], empirical tight binding (ETB) [16] and ab initio [17,18] approaches. These methods are known to be accurate in predicting electronic band structures that agree well with experiments [18], but the demand on computing time and resource makes them impractical to be employed in optimization of sophisticated photonic Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence.
Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. devices that involve heterostructures and/or nanostructure such as superlattices, quantum wells and quantum dots, especially when a large parameter space needs to be explored. In comparison, the 8-band k·p model is far more efficient in calculating the band structure near the Γ-point of the BZ, which works well for direct bandgap III-V semiconductors since essentially all electronic and optical processes occur close to the proximity of the Γpoint. However, for GeSn alloys where the Γ and L valleys are not separated far apart in energy, electrons are subject to inter-valley scattering such that both valleys are populated. It is therefore necessary to obtain full BZ band structure with a model that is computationally efficient. The 30-band k · p model developed as an expansion of the simpler 8-band approach is capable of calculating the full BZ band efficiently and has been used successfully in reproducing full-zone band structures for Si, Ge, SiGe alloy [19], α-Sn [20] and other III-V structures [21].
In this article, we present a 30-band k·p model that can be used to map out full-zone electronic energy band structures for Ge 1−x Sn x alloys where Sn composition varies from 0 to 0.3. This composition range covers most GeSn alloys currently under development for device applications. The model is anchored at reproducing the previously known band structures of Ge and α-Sn that are obtained using either ETB or ab initio method. This step establishes a set of reliable input parameters for the 30-band k · p model at the two extreme compositions: x=0, 1. The parameters needed for any Ge 1−x Sn x alloy in 0<x<0.3 are then obtained in the spirit of virtual crystal approximation (VCA) with either linear or quadratic interpolations between the two extreme points at the composition x that is adjusted with available experimental results. The electronic band structures calculated with the 30-band k·p model allow us to extract many electronic and optical properties that are important to electronic and photonic devices of all kinds made from Ge 1−x Sn x . In this work we present the results of effective masses along different crystalline directions, Luttinger parameters, and density of states (DOS) at Γ and L valleys for relaxed Ge 1−x Sn x alloys spanning across the alloy composition range 0.0<x<0.3 at room temperature. In addition, this model can be used to pinpoint the alloy composition x for Ge 1−x Sn x to transition from indirect to direct bandgap and the prediction x=7.25% is in good agreement with experimental measurement [22]. It should be pointed out that for the purpose of demonstration we choose to study relaxed GeSn at room temperature, but this model can easily be extended to include the effect of strain as well as at other temperatures. This study illustrates the efficiency and effectiveness of the 30-band k· p model in calculating the electronic band structures of the Ge 1−x Sn x alloy, allowing for extraction of practically all electronic and optical properties associated with its energy band for device simulation and optimization.

Theoretical model
The 30-band k·p model Hamiltonian matrix can be written as [19]: The 15 energy-level parameters of various Γ states from O h and the SOC strength introduced by Cardona and Pollak [23] can be obtained from theoretical methods such as ETB, MBJLDA and HSE06. The other ten parameters, P 1 , P 2 , P 3 , P 4 , Q 1 , Q 2 , R 1 , R 2 , T 1 , T 2 , used in the present 30-band k · p model are the matrix elements of the linear momentum operator p, such as l l , representing the coupling strength between different bands as shown in figure 1, which can be optimized by nonlinear methods.
Based on above 30-band k · p model, effective mass can be calculated by the finite difference methods, where the second and mixed derivatives are evaluated using five-point stencil. In addition, to accelerate the speed of computation and save computation resources in the DOS calculation, sampling BZ methods [24] is used. The kpoints and corresponding weights are obtained through the Vienna ab initio simulation package (VASP ) [25]. Besides, in order to improve the accuracy of the DOS at Γ-valley, we increase more k-point around Γ point.
All details of the ab initio calculation, effective mass calculation and parameters optimization are presented in the appendix section.

Results and discussion
The starting point of our calculation is the electron structure of single crystal Ge and α-Sn. The previous application of the 30-band k · p model to obtain the Ge band structure [19] was based on parameters at liquid helium temperature. As we know the band structure of Ge has a strong dependence on temperature, for instance, the Ge band gaps at L-valley and Γ-valley are 0.66 eV and 0.80 eV, respectively, at room temperature which are significantly lower than those at 0.747 eV and 0.89 eV at that of liquid helium [5,6]. Since essentially all GeSn-based devices are developed for room temperature operation, we recalculated the Ge band structure at room temperature using the ETB(sp 3 d 5 s * ) method and then extracted parameters from it to feed into the 30band k·p model in order to reproduce the band structure again. For convenience, we assume that the the valence band maximum is at potential zero and that all other values are referenced to it in all calculations. Results from the two different methods are shown in figure 2 where the two band structures are nicely matched globally across the full BZ. Table 1 lists eigenvalues, energy gaps, effective masses and Luttinger parameters at highsymmetry points in BZ that are obtained from experiments, ETB and the 30-band k · p methods at 300 K for Ge. The good overall agreement between them clearly validates the effectiveness of the 30-band k·p model.
For α-Sn it has been reported that ab initio calculation could result in underestimating the band gap because of its negative band gap [20], we therefore chose to employ two amended methods involving hybrid functional to improve the ab initio calculation. Two such methods are used to obtain the electronic band structure of α-Sn. One is based on modified Becke-Johnson local density approximation (MBJLDA) [26] and result is shown (red curve) in figure 3(a). The other is proposed by Heyd, Scuseria, and Ernzerhof (HSE06) [27] and its result is shown (red curve) in figure 3(b). Table 2 also lists the band gap results as table 1. From the two band structures, we have extracted two sets of parameters to feed into the 30-band k·p model. In both cases, they were able to reproduce the electronic band structure shown (blue curve) in figure 3 that matches the two ab initio hybrid functional results reasonably well, except near the X-valley where a substantial discrepancy can be observed. While this discrepancy remains to be resolved in future investigation, it does not affect the material properties that we aim to extract because Γ-valley and L-valley in conduction band are the major concerns for GeSn alloys as they sit significantly lower than the X-valley in energy. It is not difficult to see from figure 3 that although the two different hybrid functional methods have produced largely matching electronic band structure over the entire BZ, there are some inconsistencies. In the case of HSE06 method, for instance, the effective mass calculated from the Γ 7 conduction band is positive at Γ-point as shown in figure 3(b) (green circle), but its experimental value as measured turns out to be negative [5,6] which agrees with what was obtained from the MBJLDA method. The Similar discrepancy has also been found in HgTe and ScPtBi topological material with the HSE06 approach [28]. The general agreement between the two methods and the possibility of local inconsistency with the HSE06 method lead us with confidence to follow the MBJLDA method for the subsequent calculation of the electronic band structures of Ge 1−x Sn x alloys and extraction of their properties.
The VCA [19] is used to extend the 30-band k · p model to Ge 1−x Sn x alloy. Although the crystal structure of Ge 1−x Sn x alloy does not belong to group O h due to the broken centrosymmetry, the group O h model is used in order to simplify the optimization of parameters to a reasonable level [19]. We have used quadratic interpolation between Ge and α-Sn to come up with the coupling strengths (table 3) as well as the necessary energy levels of different bands and SOC splitting at Γ point (table 4) as input parameters for the 30-band k·p model for each alloy composition of Ge 1−x Sn x . Bowing effect in both conduction band minimum and SOC in valence band is considered. For the coupling matrix elements, we are mainly concerned with those that couple with the conduction and valence bands. These input parameters are optimized in order for the k · p model to produce bandgap results that best match those in [17] and to yield effective masses close to what were obtained in [14].
Having done the optimization, the band structure of Ge 1−x Sn x alloy is calculated by diagonalizing the alloy Hamiltonian and depicted in figure 4 with three different Sn composition. Besides, we are more interesting in the band gap variation with Sn composition. Figure 5 shows the dependence of bandgap at L-valley, E g (L) (solid blue), and at Γ-valley, E g (Γ) (solid red) on the Sn composition x, calculated using the 30-band model. The   figure 5) [29][30][31][32] in terms of bandgap values and rate of bandgap decrease with Sn composition. The crossover of L and Γ-valley occurs at the Sn composition x=7.25% according to our 30-band k·p model, which is in good agreement with the empirical predication of x=7.18% [17] and the experimental result of x≈7.8% [22]. Some of the discrepancies can be accounted for by the different bandgaps used in the calculations for Ge as shown in table 1. The electronic band structure obtained from 30-band model allows for extraction of important material parameters such as effective mass of electrons at L and Γ-valleys, Luttinger parameters and DOS effective mass, that directly influence electronic and optical processes in Ge 1−x Sn x alloy. For such an immature material system Table 2. Effective masses and energy gaps obtained from literature, MBJLDA, HSE06 and k·p methods for α-Sn. All energies are in eV and the effective mass unit is m 0 that is the mass of free electron. Corresponding literature data are from [5].  Table 3. Matrix elements of the linear momentum p used in the present 30-band k · p model. The parameters are listed in a.u. as previous SiGe alloy [19] and eV nm for calculation.

Matrix elements
with much to be explored, linear extrapolation of these parameters is often used, but their accuracies are obviously unreliable, and therefore not recommended [33]. To the best of our knowledge, these important parameters have not been reported previously. Figure 6 shows 111 [14]. As the Sn composition exceeds the crossover point, the bandgap of Ge 1−x Sn x alloy becomes direct, the availability of Luttinger parameters enables the use of simpler 8-band k · p model as a convenient method to   [17]. The scatter points are experimental data [29][30][31][32]. The zoom inset shows the indirect-direct crossover at the Sn composition of 7.25% for our 30-band k·p model compared with that of 7.18% from [17].
investigate various phenomena that take place around the Γ point such as the vertical optical transitions associated with optical absorption and emission processes. Luttinger parameters versus the Sn composition are shown in figure 7 where a strong nonlinear dependence is revealed.
At last, DOS in conduction band and valence bands, which are key features determining a range of electronic and optical properties that impact device performance, are calculated for Ge 1−x Sn x alloy with different compositions. Figures 8(a) and (b) show the DOS in L and Γ-valley in the conduction band as a function of the electron energy from their respective band edge, respectively. It is not difficult to see that the L-valley DOS is much higher than that in Γ-valley as expected. Figure 8(c) shows the combined conduction band DOS which is the summation of L and Γ-valley DOS. It can be seen that as the Ge 1−x Sn x alloy makes the transition from indirect to direct bandgap (x>7.25%), the total conduction band DOS around its band edge reduces, positively

Conclusion
In conclusion, we present a 30-band k · p model that can be used efficiently to calculate the electronic band structures for relaxed Ge 1−x Sn x alloy of various compositions. This model was first validated to show consistency in calculating the electronic band structures of the α-Sn and Ge by ab initio and ETB methods, respectively. Two sets of input parameters were optimized for band structure calculation at various Sn compositions: (1) coupling matrix elements and (2) eigen energy levels and SO splitting at Γ point for various Sn  compositions, 0<x<0.3. Based on this model, the indirect-to-direct crossover is determined to take place at the Sn composition of 7.25% which is good agreement with previous predication and experimental result. Bowing parameters that describe the relationship between bandgaps at Γ and L-valley and Sn composition are obtained. The calculated band structure allows for extraction of electron effective mass at Γ valley, longitudinal and transverse electron effective mass at L-valley, as well as HH and LH effective masses along different crystalline directions. From these HH and LH effective masses, Luttinger parameters that characterize the valence band around Γ-valley in Ge 1−x Sn x alloy can be established. Finally, DOS in both conduction and valence band as well as their associated effective-mass DOS can be calculated as well. The 30-band k· p model fulfills the need for an efficient and effective tool that can be used in calculating the electronic band structure of Ge 1−x Sn x alloy across its entire BZ, establishing material parameters, and optimizing device performance. the lattice constant. The results were nearly the same.
At last, we gave a more clear description about the process of optimization 30-band · k p model parameters as shown in [19], which could be treated as an optimization process that includes two steps. First, 30-band k · p parameters of single crystal Ge and Sn are optimized. The target of this step was that the band structure calculated by 30-band k · p model with all optimized parameters should agree well with the results from ETB and ab initio. Obviously, it was a multiple variables optimization problem and can be solved by the hill climbing optimization method. The figures 2 and 3 confirmed the validity of our 30-band k·p optimized parameters. Secondly, 30-band k · p parameters of GeSn alloy are optimized in the frame of VCA where the the Hamiltonian of alloy can be written : = are from the first step. Although here the VCA was implemented between Ge at room temperature and α-Sn at 0 K, similar method has been employed in [16,34]. Now the optimization target was that the band gap variation profile at Γand L-valley should match with the description from empirical analytical formula [17] according to the experimental results at room temperature and the effective mass variation profile of electron, HH and LH at Γ-valley should be consistent with the trend in [14]. Since 30-band k·p model was very easy to solve numerically, Ge 1−x Sn x alloys with 100 different Sn composition between 0 and 0.3 are used to optimize 30-band k · p parameters. Noting that the k · p model is not related with microscopic configurations. The curves calculated by 30-band model with optimize parameters were consistent with the dashed lines in figure 5. Again, it confirmed the validity of our methods and parameters. The hill climbing optimization method played a key role in the whole process.