LIGHT-EMITTING OR LIGHT-ABSORBING COMPONENT
20220135878 · 2022-05-05
Inventors
- Silvana BOTTI (EINDHOVEN, NL)
- Friedhelm BECHSTEDT (EINDHOVEN, NL)
- Jozef Everardus Maria HAVERKORT (EINDHOVEN, NL)
- Erik Petrus Antonius Maria Bakkers (Eindhoven, NL)
- Elham FADALY (EINDHOVEN, NL)
- Alain DIJKSTRA (EINDHOVEN, NL)
Cpc classification
H01L33/34
ELECTRICITY
H01L31/035227
ELECTRICITY
B82Y20/00
PERFORMING OPERATIONS; TRANSPORTING
C01P2002/76
CHEMISTRY; METALLURGY
H01L33/16
ELECTRICITY
H01L31/036
ELECTRICITY
H01L31/028
ELECTRICITY
H01L33/24
ELECTRICITY
International classification
Abstract
The invention relates to a light-emitting component comprising a light-emitting section consisting of a Hex-Si.sub.1−xGe.sub.x compound material, said Hex-Si.sub.1−xGe.sub.x compound material having a direct band gap for emitting light.
The invention also pertains to a light-absorbing component comprising a light-absorbing section consisting of a Hex-S.sub.1−xGe.sub.x compound material, said Hex-Si.sub.1−xGe.sub.x compound material having a direct band gap for absorbing light.
Claims
1-9. (canceled)
10. A light-emitting component comprising a light-emitting section comprising a Hex-Si.sub.1−xGe.sub.x compound material, the Hex-Si.sub.1−xGe.sub.x compound material having a direct band gap for emitting light.
11. The light-emitting component according to claim 10, wherein the Hex-Si.sub.1−xGe.sub.x compound material is structured to emit light with a B-coefficient for radiative emission of 0.7×10.sup.−10 cm.sup.3/s <B.sub.rad<11×10.sup.−10 cm.sup.3/s at 300K.
12. The light-emitting component according to claim 10, wherein the Hex-Si.sub.1−xGe.sub.x compound material is structured to emit light with a B-coefficient for radiative emission of 0.7×10.sup.−10 cm.sup.3/s <B.sub.rad<8.3×10.sup.−10 cm.sup.3/s at 300 K.
13. The light-emitting component according to claim 10, wherein the Hex-Si.sub.1−xGe.sub.x compound material is structured to emit light between 1.8 μm for x=0.65 and 3.5 μm for x=1.0.
14. The light-emitting component according to claim 10, wherein the Hex-Si.sub.1−xGe.sub.x compound material comprises strained quantum well structures of a different composition of the hexagonal Si.sub.1−xGe.sub.x compound material structured to emit light between 1.5 μm and 7.0 μm.
15. The light-emitting component according to claim 10, wherein the Hex-Si.sub.1−xGe.sub.x compound material exhibits direct band gap emission with a sub-nanosecond recombination lifetime.
16. The light-emitting component according to claim 10, wherein the Hex-Si.sub.1−xGe.sub.x compound material exhibits a linear dependence of the photoluminescence intensity versus the excitation power.
17. The light-emitting component according to claim 10, wherein x of the Hex-Si.sub.1−xGe.sub.x compound material is defined with 0.2<x<1.0, or with 0.6<x<1.0, or with 0.2<x<0.99, or with 0.2<x<0.9, or with 0.6<x<0.9, or with 0.6<x<0.99.
18. The light-emitting component according to claim 10, comprising a monolithic structure including a Cub-Si substrate provided with the Hex-Si.sub.1−xGe.sub.x compound material as the light-emitting section.
19. A light-absorbing component comprising a light-absorbing section including a Hex-Si.sub.1−xGe.sub.x compound material having a direct band gap for absorbing light.
20. The light-absorbing component according to claim 19, wherein x of the Hex-Si.sub.1−xGe.sub.x compound material is defined with 0.2<x<1.0, or with 0.6<x<1.0, or with 0.2<x<0.99, or with 0.2<x<0.9, or with 0.6<x<0.9, or with 0.6<x<0.99.
21. The light-absorbing component according to claim 19, comprising a monolithic structure including a Cub-Si substrate provided with the Hex-Si.sub.1−xGe.sub.x compound material as the light-absorbing section.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0023] Notwithstanding any other forms which may fall within the scope of the methods and devices as set forth in the Summary, specific embodiments of the methods and devices will now be described, by way of example, and with reference to the accompanying drawings in which:
[0024]
[0025]
[0026]
[0027]
[0028]
[0029]
[0030]
[0031]
[0032]
[0033]
[0034]
[0035]
[0036]
[0037]
[0038]
[0039]
[0040]
[0041]
[0042]
[0043]
[0044]
[0045]
[0046]
[0047]
[0048]
DETAILED DESCRIPTION OF THE DRAWINGS
[0049] Fluidly merging the fields of integrated electronics and photonics using a single materials family and CMOS compatible process flow is the major goal of distributed information technologies. Silicon (Si) crystallized in the usual cubic (diamond) structure has dominated the electronics industry for more than half a century. However, cubic silicon (Cub-Si), germanium (Cub-Ge) and Si.sub.1−xGe.sub.x-alloys are all indirect band gap semiconductors that cannot emit light efficiently.
[0050]
[0051] The band structure of Cub-Si, presented in
[0052] As shown in
[0053] To investigate how the direct band gap energy can be tuned by alloying Ge with Si, the band structures of Hex-Si.sub.1−xGe.sub.x (for 0<x<1) were calculated using ab initio density functional theory (DFT) and a cluster expansion method for isostructural hexagonal binary alloys, using so-called ab-initio calculations detailed in the paragraph below. Selected results, presented in
[0054] As to the ab-initio calculations it should be noted that according to the invention all calculations were performed using density functional theory (DFT) as implemented in the Vienna Ab Initio Simulation Package (VASP) A with the projector augmented wave method. A plane-wave cutoff of 500 eV was used and Ge 3d electrons were included as valence electrons. Brillouin zone integrations were carried out using 12×12×6 ┌-centered k-point grids for lonsdaleite Ge and 12×6×6 ┌-centered k-point grids for Si-Ge, ensuring a convergence of total energies to 1 meV/atom. For structural calculations, the PBEsol exchange-correlation potential was used, together with a convergence threshold of 1 meV/A on Hellmann-Feynman forces. The modified Becke-Johnson exchange potential in combination with local density approximation (MBJLDA) was preferred for electronic structures and optical properties, as it guarantees band gaps in excellent agreement with experiments available in the prior art. Spin-orbit coupling were included in all calculations.
[0055] Alloys are studied using a cluster expansion method for isostructural lonsdaleite binary alloys. For the cluster expansion, the macroscopic alloy is divided into clusters of 8 atoms obtained from the primitive wurtzite (WZ) unit cell. In this way, it is possible to study 46 different structures ranging from pure Ge to pure Si. This method becomes more accurate with increasing size of the clusters, and it was verified that the thermodynamic averages are not significantly modified by performing calculations with 16 atom clusters. The radiative lifetime (τ.sub.rail) at temperature (T) is calculated using the formula:
[0056] where A.sub.cvk denotes the radiative recombination rate for vertical optical transitions between a conduction state |ck>and a valence state |vk>, with one-particle energies ε.sub.ck and ε.sub.vk, and Fermi occupation functions ƒ.sub.ck and ƒ.sub.vk and w.sub.k the k-point weight. In order to reproduce experimental conditions, n=10.sup.19 cm.sup.−3 charge carriers due to n-doping in the conduction band were included, and modified accordingly the chemical potential of electrons. The radiative recombination rate is given by:
[0057] where n.sub.eƒƒ is the refractive index of the effective medium (here set approximately to the experimental value for cubic Ge for which n.sub.eƒƒ=5). The squares of the momentum matrix elements can be either averaged over all directions corresponding to the emission of unpolarized light, as in Eq. (A2), or only the in-plane component is considered, for light polarized perpendicularly to the wire axis. Denser k-point grids were necessary to calculate lifetimes (72x72x36 for lonsdaleite Ge and 24x12x12 for Si-Ge).
[0058]
[0059] Remarkably, the radiative lifetimes of Hex-Si.sub.1−xGe.sub.x alloys are significantly lower than that of pure Hex-Ge, for which the lowest energy transition is dipole forbidden at the ┌-point. This observation can be traced to the reduced symmetry in the random Hex-Si.sub.1−xGe.sub.x alloys, which leads to mixing of Ge s-states into the lowest conduction band wave function.
[0060] It has been found, that the calculated lifetimes of the Hex-Si.sub.1−xGe.sub.x alloys are approaching those of III-V semiconductors, such as GaAs. As to a comparison with group III-V semiconductors, the measured lifetime of
[0061] Hex-Si.sub.0.2Ge.sub.0.8 at low temperature is very comparable to the recombination lifetimes reported in prior art literature for III-V compound semiconductors, which are generally of the order of 1 ns. Prior art literature reported a temperature independent lifetime of 1 ns in core/shell GaAs/AlGaAs nanowires, very similar to the yet unpassivated Hex-Si.sub.0.2Ge.sub.0.8 nanowire shells according to the invention. The comparison of the quenching ratio of the integrated photoluminescence intensity when increasing the temperature from 4 K to 300 K, compares quite favorable for Hex-SiGe where this ratio varies between a factors of 15-100 as shown in
[0062] Examples according to the invention relate to Ge-rich alloys of Hex-Si.sub.1−xGe.sub.x as they combine a direct band gap, strong optical transitions and wavelength tunability. Here, it is demonstrated experimentally that Ge-rich alloys of Hex-Si.sub.1−xGe.sub.x are indeed direct gap semiconductors, observe strong emission, and a temperature independent nanosecond radiative lifetime. The results are shown to be in remarkable quantitative agreement with theoretical predictions.
[0063] In
TABLE-US-00001 TABLE S1 Extracted Lattice Parameters of the Hex-Si.sub.1−xGe.sub.x samples: Wurtzite lattice parameters of all measured Hex-Si.sub.1−xGe.sub.x samples with corresponding error-values extracted from XRD measurements. Ge-Content Lattice parameter (a) (Å) Lattice parameter (c) (Å) 1.0 3.9855 ± 0.0003 6.5772 ± 0.0003 0.92 3.9789 ± 0.0001 6.5542 ± 0.0001 0.86 3.9649 ± 0.0005 6.5431 ± 0.0004 0.75 3.9505 ± 0.0008 6.5257 ± 0.0001 0.63 3.9206 ± 0.0000 6.4790 ± 0.0005
[0064] In an example of the method according the invention Ge-rich Si.sub.1−xGe.sub.x alloys are grown around a thin (˜35 nm diameter) WZ gallium arsenide (GaAs) core that is lattice matched to Ge as shown in
[0065]
[0066] For all measured samples at least three individual hexagonal reflections have been measured. For the pure Hex-Ge sample, the azimuth was varied to enhance the fidelity of the extracted lattice parameters. In addition, a cubic GaAs substrate reflection has always been used as an anchor in reciprocal space to correct for any possible alignment offsets. From the measured symmetric reflections, as shown in the full series in
[0067] To accurately determine the peak positions, all RSMs were corrected according to the peak positions of the cubic GaAs-substrate reflections to eliminate any angular alignment offsets. Then a 2D-Gauss fit was performed on the data-sets in q-space before gridding, to reduce the influence of possible artefacts coming from the gridding-routine. For plotting the dataset, the irregularly spaced q-coordinates, as measured and transformed from the angular-space, have been gridded into a regularly spaced q-coordinate system.
[0068] The combined results from the XRD measurements can be found in table S1 where the measured lattice parameters are given for each measured Ge-concentration. For all samples the influence of the WZ-GaAs core material on the Si.sub.1−xGe.sub.x lattice parameter can be neglected because of the fact that a relatively thin GaAs core (around 35 nm) is surrounded by a thick (several 100 nm) Si.sub.1−xGe.sub.x shell. Hence, the crystalline properties of the Hex-Si.sub.1−xGe.sub.x shell dominate the whole structure. Furthermore, Hex-Ge and WZ-GaAs are nearly lattice matched (see lattice parameter of WZ-GaAs which implies that basically no strain in the shell is expected for the samples with high Ge-concentrations (>60%) as also confirmed by FEM-simulation. This is an important aspect since it confirms the high fidelity of the found lattice parameters, especially for the lattice parameter of pure Hex-Ge. The errors given in table S1 consider the accuracy of defining the peak position with a 2D-fit as described, as well as the scattering of the individual lattice parameter values extracted from the evaluation of multiple peaks. The instrumental resolution can be neglected for the error estimation, since the contribution to the given errors will be much smaller than the total error-values.
[0069] The XRD measurements have been carried out at the Deutsches-Elektronen-Synchrotron (DESY) in Hamburg, at the high-resolution-diffraction beamline P08. For the diffraction experiments a high precision 6-circle diffractometer has been used, the photon energy was set to 15 keV with a corresponding wavelength of 0.8266 Å. The energy was carefully chosen to ensure a high photon flux while still being able to access higher-indexed reflections, needed for the precise measurements of the lattice parameters. The x-ray beam has been shaped by a slit system and the resulting spot-size on the sample was 200 μm (horizontal)×100 μm (vertical), a size sufficient to illuminate a few thousands of wires at once. For measuring the scattered signal coming from the wires, a Dectris—“Mythen” 1D X-ray detector has been used; this detector offers a high dynamic range and, due to the small pixel size (50 μm), an increased angular resolution in 2θ, compared to most 2D detectors. For the conversion of the measured angular coordinates to reciprocal space coordinates and all further data processing, such as 2D-peak-fitting and post-processing for plotting, the freely available software library “Xrayutilities” in combination with Python 3.6 has been used.
[0070]
[0071] These results are in good agreement with the TEM measurements performed for the same samples, as shown
[0072] The structural quality of the crystals was investigated by transmission electron microscopy (TEM). Two different sample preparation methods were used. In the standard axial analysis, nanowires were mechanically transferred to a holey carbon TEM grid. Concerning the cross-section TEM studies, nanowires were prepared using a Focused Ion Beam (FIB). In both cases, high resolution TEM and Scanning TEM analyses were conducted using a JEM ARM200F probe-corrected TEM operated at 200 kV. For the chemical analysis, Electron Dispersive Xray (EDX) spectroscopy measurements were carried out using the same microscope equipped with a 100 mm.sup.2 EDX silicon drift detector. TEM lamellae were prepared in a FEI Nova Nanolab 600i Dual beam system. For this, the nanowires were initially transferred with the aid of a Kleindiek nano-manipulator from the growth substrate to a piece of Si and then arranged to lie parallel to each other. These nanowires were covered with electron- and ion-beam induced metal deposition to protect them during the procedure. The lamella was cut out by milling with 30 kV Ga ions and thinned down with subsequent steps of 30, 16, and 5 kV ion milling in order to minimize the Ga-induced damage in the regions imaged with TEM.
[0073] By exploring the optical properties of the Hex-Si.sub.1−xGe.sub.x nanowires probed using power and temperature dependent photoluminescence spectroscopy, shown in
[0074]
[0075] As to the fitting using the Lasher-Stern-Würfel (LSVV) model, the observed photoluminescence spectra of Hex-Ge and Hex-SiGe all consist out of a single peak. The observation of a single photoluminescence peak is attributed to a band-to-band (BtB) recombination. The absence of excitonic effects at low temperatures is due to an As-doping level of 9.10.sup.18 cm.sup.−3 as deduced by Atom Probe Tomography shown in
[0076] To accurately establish whether the observed photoluminescence is due to BtB recombination, the experimental spectra were fitted with the Lasher-Stern-Würfel (LSVV) model. This model, that predicts the shape of a photoluminescence peak, is derived from the Planck-Einstein radiation law and is given by:
[0077] In this equation Δμ is the splitting of the quasi-fermi levels and a(E) is the absorptivity. In modelling the absorptivity, parabolic bands have been assumed. Corrections for an Urbach tail and an excitation dependent Burstein-Moss shift have been made in analogy in the prior art. Both the temperature dependent and the excitation power dependent photoluminescence measurements were fitted as shown in
[0078] The high quality fits by the LSW model unambiguously show that the observed photoluminescence is exclusively due to BtB recombination. It is of paramount importance for the analysis that measured recombination lifetimes are due to BtB recombination and not due to e.g. an impurity or defect related optical transition. According to an advantageous effect according to the invention, it is noted that the deduced carrier temperature exceeds 700 K at the highest excitation densities.
[0079]
[0080] As to the temperature dependence of the fundamental band gap, although the temperature dependence of the fundamental band gap is most often described by the Varshni equation, the Vina equation provides a more accurate description for elevated temperatures
[0081] in which a is a constant, b represents the strength of the electron-phonon interaction, and θ is the Debye temperature of the material. For the band gap of Hex-Ge the Vina equation is fitted in
[0082] The shrinkage of the Si.sub.0.20Ge.sub.0.80 band gap, which is displayed in
[0083]
[0084]
[0085]
[0086] The radiative lifetime as well as the radiative emission efficiency of Hex-Si.sub.0.20Ge.sub.0.80 are deduced. It is important to note that the measured decay lifetime is determined by the fastest recombination process, which can be either radiative or non-radiative in nature. It is therefore crucial to choose experimental conditions in which the measured recombination lifetime is exclusively governed by pure radiative recombination. This can be achieved at low temperature, since non-radiative processes are commonly thermally activated and therefore negligible.
[0087] As to the temperature dependence of the integrated photoluminescence intensity a detailed Arrhenius analysis of the temperature dependence of the integrated PL is provided as presented in
[0088] Following the LSW analysis as described above, it was concluded that the photoluminescence spectrum can be explained by Band-to-Band (BtB) recombination with only a minor influence of the acceptor related transition. As a consequence, the rate equation model was limited to a three level system incorporating the conduction band, the valence band and a “killer defect” which is characterized by an activated non-radiative recombination lifetime. The one-center model in the classification of Reshchikov was then used, which is explained in more detail by a configuration coordinate diagram. In this one-center model, the internal quantum efficiency (η.sub.int) for radiative emission varies with temperature according to the ratio of the radiative recombination rate divided by the total recombination rate by η.sub.int=τ.sub.r.sup.−1/(τ.sub.r.sup.−1+τ.sub.nr.sup.−1(T)). The low excitation data collected at 68 W/cm.sup.2, which are presented in
[0089] The excellent quality of the Arrhenius fit provides evidence that the non-radiative recombination into the yet unknown killer defect can indeed be explained by an activated non-radiative recombination rate.
[0090] The temperature dependence of the photoluminescence intensity can thus be expressed as
[0091] in which the photoluminescence quenching rate into the non-radiative center is given by
In most semiconductors, different non-radiative recombination centers exist which feature e.g. activation energy E.sub.A, E.sub.B and quenching rates R.sub.A, R.sub.B resulting in
[0092] It is instructive to perform this analysis to three different generations of Hex-SiGe samples which are specified in table S2 and whose Arrhenius plots are shown in
when both non-radiative channels are fully activated (room temperature). The first quenching mechanism seems to have disappeared in sample B which was grown at a higher temperature. In sample B, we only observe photoluminescence quenching above a temperature of 100 K, which is again tentatively attributed to be at least partially due to surface recombination. The activation energy E.sub.B=4±5 meV is tentatively explained by the de-trapping from localized states due to alloy fluctuations in the Hex-SiGe nanowire shell. Once the carriers are de-trapped, they will quickly diffuse to the nanowire surface where they recombine non-radiatively. In sample A, both quenching mechanisms have disappeared as
at an excitation density of 36 kW/cm.sup.2, thus showing that sample A remains in the radiative limit up to 220 K. The quality of sample A is probably higher due its thick hex-SiGe shell which reduces the amount of surface recombination as well as by its length which reduces the influence of re-evaporating arsenic (As) and gallium (Ga) from unwanted growth on the substrate. To be completely sure, we have regrown sample A resulting in an identical temperature dependence as the first grown sample.
[0093]
TABLE-US-00002 TABLE S2 Growth parameters of Hex-SiGe samples with increasing quality: Here listed are the input temperature of the MOVPE reactor for the SiGe-shell growth and the resulting geometry of the nanowires that are presented in FIGS. 4B and 4B. nanowire Growth nanowire Core nanowire Temp Ge-content Shell Radius Diameter Length Sample # (° C.) (%) (nm) (nm) (μm) Sample A 700 79 650 175 8 Sample B 700 80 400 35 2.5 Sample C 650 75 150 35 2.5
[0094]
[0095]
[0096]
[0097] As to the temperature dependence of the radiative lifetime, as shown in
[0098] Non-radiative recombination features an activated behavior at low temperature which is governed by τ.sub.nr.sup.−1(T)=τ.sub.n,r0.sup.−1e.sup.−E.sup.
in which E.sub.t is the trapping level, E.sub.i is the intrinsic Fermi level and τ.sub.p,0 is the SRH lifetime for minority holes. At higher temperature, the SRH lifetime is expected to decrease with T.sup.−1/2 due to the fact that both τ.sub.n0 and τ.sub.p0 are inversely proportional to the thermal velocity. It is concluded that it is clearly not possible to interpret the observed temperature independent recombination lifetimes as being due to non-radiative recombination.
[0099] Auger recombination is discussed, which might be expected due to the high n-doping by unintentional arsenic (As) incorporation during growth. The Auger rate includes two different processes, the nnp-Auger process in which the excess energy is transferred to an electron and the npp-Auger process in which the excess energy is transferred to a hole. With the present invention, high n-doping exists due to As incorporation during growth, resulting in a doping concentration n.sub.0. It is expected that the nnp-Auger process will be most important in the n-doped Hex-SiGe samples. The Auger coefficients are however temperature dependent, which results in a T-dependent recombination lifetime, which is not consistent with the observations. Most importantly, as shown in the inset of
[0100] As to the possibility of radiative recombination, the radiative lifetime for an intrinsic semiconductors increases with T.sup.3/2 showing sub nanosecond radiative lifetimes at low temperature which increase to more than a microsecond at room temperature. For a degenerately doped semiconductor, the radiative lifetime is expected to be temperature independent since the B-coefficient for radiative recombination is proporuonai to
in which L is the spontaneous radiative recombination rate. It can be easily seen that for a degenerate semiconductor p∝T.sup.3/2, L∝T.sup.3/2 while n becomes temperature independent. Both the B-coefficient for radiative recombination rate and the radiative lifetime are thus expected to be independent of temperature.
[0101] Photoluminescence lifetime measurements for all three samples are presented in
[0102] In order to again obtain the correct statistics, photoluminescence lifetime measurements were performed on more than 60 different nanowires taken from sample A at 4 K and at 300 K. The data are displayed in
[0103] Moreover, measurements were performed at high excitation density to saturate non-radiative processes and to maintain the radiative limit up to increased temperature. Typical results from time-resolved luminescence measurements on a single wire from the Si.sub.0.20Ge.sub.0.80 sample are presented in
[0104]
[0105] To be sure that the data for an individual wire are representative, more than 60 individual wires swiped from high crystal quality sample A were analyzed, which are presented in
[0106] As to the excitation power dependence of the integrated photoluminescence intensity it is noted, that at low excitation density, Δn<n.sub.0, the nonradiative, radiative and Auger processes all yield a linear dependence of the PL-intensity versus excitation power with a slope of unity. However, this simplified analysis assumes that the non-radiative recombination centers are not being saturated.
[0107] Since no deviation from a linear behavior was observed, the data suggest that, even if non-radiative recombination centers would be present, it is possible to saturate them with excitation power. This suggests that no non-radiative recombination centers in the bulk of the material, implying that the radiative limit is reached. It is noted that this argument applies both for Δn<n.sub.0 and Δn>n.sub.0.
[0108] At high excitation density, Δn>n.sub.0, a prior art analysis is used, according to which analysis the total carrier generation rate G should be equal to the total carrier recombination rate by
G=An+Bn.sup.2+Cn.sup.3 (A7)
[0109] in which An is the Shockley-Read-Hall nonradiative recombination rate, Bn.sup.2 is the radiative recombination rate and Cn.sup.3 is the Auger nonradiative rate. At high excitation density (which is above 500 W/cm.sup.2 for Hex-Ge as shown by bandfilling in
[0110] In the plot of the integrated photoluminescence intensity versus excitation density, Eq. A6 yields a slope of two for non-radiative recombination (provided that the non-radiative recombination centers are not being saturated, see above), a slope of unity for radiative recombination and a slope of 2/3 for Auger recombination. It is noted that a decrease of the PL-intensity at the highest excitation power was not observed, providing a first indication that Auger recombination losses are not yet dominant in this material.
[0111] For the Hex-Si.sub.0.20Ge.sub.0.80 sample, a clear boundary between the Δn<n.sub.0 and the Δn>n.sub.0 regime could not be established due to the added complication of alloy broadening. Most probably, the Si.sub.0.20Ge.sub.0.80 alloy will be composed out of compositional pockets in which either Δn<n.sub.0 or Δn>n.sub.0 applies. The observation of a slope of exactly unity, as shown in the inset of
[0112] As to the radiative efficiency and B-coefficient of Hex-SiGe, in order to compare the radiative emission strength of Hex-SiGe with other well-known direct band gap semiconductors like e.g. GaAs or InP, the radiative emission rate is compared at room temperature which is most relevant for device applications. By making the comparison at 300 K, excitonic effects as well as effects due to carrier localization in the Hex-SiGe alloy are not relevant anymore. The key parameter to compare the radiative efficiency of a semiconductor is the B-coefficient which is a recombination rate, corrected for the doping density.
[0113] The radiative rate per volume of a semiconductor R.sub.rad can be expressed in terms of the B-coefficient, n- and p-type doping concentration n.sub.0 and p.sub.0 and the number of excited electron-hole pairs Δn=Δp. For a highly n-doped semiconductor, which yields n.sub.0>>Δn, R.sub.rad can be expressed as:
R.sub.rad=B.sub.rad(n.sub.0Δn)(p.sub.0+Δ.sub.P)≈B.sub.radn.sub.0Δ.sub.p (A9)
[0114] The experimentally observed radiative lifetime τ.sub.rad is determined by the recombination rate per volume R.sub.rad and the number of excited electron hole pairs Δn=Δp such that τ.sub.rad=Δp/R.sub.rad. Combining this result with equation A9 gives a definition for the B-coefficient of:
[0115] In which τ.sub.rad is the radiative lifetime at 300 K and n.sub.0 is the activated donor density. To determine the B.sub.rad coefficient the determined values are evaluated for τ.sub.rad and the doping density n.sub.0.
[0116] The measured photoluminescence-lifetimes show a spread over different wires as shown in
[0117] It is argued whether the measured photoluminescence decay time at 300 K, is equal to the radiative lifetime. The main supportive argument is provided by
[0118] A second point of concern might be whether the degenerate limit Δn<n.sub.0 is still valid. The main evidence for this point is that, for most wires, an excitation power is measured independent photoluminescence decay time in the same excitation range as in the inset of
[0119] Since the measurements were performed in the radiative limit and the carriers accumulate in the direct minimum at low temperature, the observation of a direct band gap emission with a sub-nanosecond recombination lifetime was concluded. The data reveals conclusive evidence for Hex-Si.sub.1−xGe.sub.x (0.65<x<1) as being a new class of direct band gap semiconductors with a large optical matrix element according to an aspect of the invention. Subsequently the radiative transition rate of Hex-SiGe was compared with other direct band gap semiconductors. The radiative transition rate R.sub.rad is quantified by R.sub.rad=B.sub.rad.Math.n.Math.p in which n and p are the electron and hole densities and B.sub.rad is the coefficient for radiative recombination, which is directly related to the transition dipole moments. The coefficient B.sub.rad can be deduced from a measurement of the pure radiative lifetime, τ.sub.rad by
in which n.sub.0 is the activated donor density.
[0120] The donor density n.sub.0 has been estimated using two techniques, the first of which is atom probe tomography, shown in
[0121] Now combining the upper bound for the donor density of 9.Math.10.sup.18 cm.sup.−3 with the upper bound of 1.6 ns for the radiative lifetime, a lower bound is obtained for the B-coefficient of 0.7.Math.10.sup.−10 cm.sup.3/s, which is roughly 2× smaller than the B-coefficient of InP. Using the lower limits for n.sub.0 and τ.sub.rad an upper limit of 11.Math.10.sup.−10 cm.sup.3/s is found for the B-coefficient, which is 9× larger as for InP. A comparison of B-coefficients of different III-V materials, Cub-Si and Hex-Si.sub.0.2Ge.sub.0.8 is made in table S3 below. Extracting the B-coefficient and thus the transition matrix elements is of great importance for possible device applications of Hex-SiGe for e.g. lasers, modulators, detectors and LEDs which all critically depend on the strength of the light-matter interaction.
[0122] As to the radiative efficiency and B-coefficient of Hex-SiGe, 0.7.Math.10.sup.−10 cm.sup.3/s<B.sub.rad<11.Math.10.sup.−10 cm.sup.3/s at 300 K is obtained, which is comparable in magnitude to GaAs and InP and almost 5 orders of magnitude larger than for Cub-Si as shown in Table S3 below. Hex-Si.sub.1−xGe.sub.x is thus a fully silicon compatible semiconductor with a radiative emission strength comparable to a direct band gap III-V semiconductor.
TABLE-US-00003 TABLE S3 Radiative coefficients of Hex-Si.sub.0.2Ge.sub.0.8, GaAs, InP and Cubic-Si: Listed are the B-coefficients for Hex-Si.sub.0.2Ge.sub.0.8 as calculated and described above S8 as to the radiative efficiency and B-coefficient of Hex-SiGe, and values for GaAs, InP and Cub-Si as known in the prior art. Hex-Si.sub.0.2Ge.sub.0.8 GaAs InP Cub-Si 0.7 .Math. 10.sup.−10 cm.sup.3/s − 11 .Math. 10.sup.−10 cm.sup.3/s 3.5 .Math. 10.sup.−10 cm.sup.3/s 1.2 .Math. 10.sup.−10 cm.sup.3/s 4.73 .Math. 10.sup.−15 cm.sup.3/s
[0123]
[0124] Hereto, time-correlated single photon counting measurements have been performed on single Si.sub.0.2Ge.sub.0.8 wires. The wires have been mechanically transferred onto a silicon wafer with a chromium (15 nm), Gold (300 nm) and SiO.sub.x (12 nm) top layer to act as a back mirror. This approach enhances the measured intensity and avoids potential optical signals emitted by the wafer. The samples with transferred Si.sub.0.2Ge.sub.0.8 wires were mounted in an Oxford Instruments HiRes2 helium flow cryostat and were excited with a 1030 nm, NKT ONEFIVE Origami femto-second pulsed laser with a 40 MHz repetition rate. The photoluminescence signal was measured in a backscattering geometry using a 36X gold coated cassegrain objective which focused the excitation laser to a spot of ˜3 μm. The laser was filtered out of the PL signal using a 1350 nm long pass filter. Using an achromatic lens the PL signal was then focused onto a SM2000 single mode fiber and fed to a Single Quantum superconducting-nanowire-single-photon-detector which was optimized for a >35% quantum efficiency at 1800 nm and a >15% quantum efficiency at 2000 nm. The 1350 nm long pass filter in combination with the SM2000 fiber defined a spectral interval of 1350 nm to ˜2300 nm over which PL was integrated. The time correlations between a laser pulse and a detection event were measured and counted using a PicoQuant PicoHarp 300 module. The overall instrumental response function (IRF) shows a FWHM of 96 ps with a decay time of τ.sub.IRF=21 ps which is the minimum observable decay time of the system. All measurements presented in
[0125] According to an aspect of the invention the direct nature of the band gap of Hex-Si.sub.1−xGe.sub.x has been established and the size of the direct band gap can be tuned via compositional engineering.
[0126]
[0127]
[0128] Direct band gap Hex-Si.sub.1−xGe.sub.x opens a pathway towards tight monolithic integration of Hex-Si.sub.1−xGe.sub.x light sources with passive cubic Si-photonics circuitry on the same chip. This will reduce stray capacitances thereby increasing performance and reducing energy consumption which is important for green information and communication technologies. According to an aspect to the invention on a conventional Cub-Si substrate of a light-emitting or light-absorbing component a Hex-Si.sub.1−xGe.sub.x compound structure is provided as the light-emitting or light-absorbing section. Possible integration routes are strain-induced transformation of Hex-Si.sub.1−xGe.sub.x, for instance by a dielectric (i.e. SiO.sub.x or SiN.sub.x) strain envelope, or alternatively by template-assisted selective area growth of the hexagonal phase.
[0129] According to an another aspect of the invention, the usable wavelength region of the Hex-Si.sub.1−xGe.sub.x compound material for emitting light can be extended towards a range between 1.4 μm and 7.0 μm using strain. This includes the application of strain on the Hex-Si.sub.1−xGe.sub.x compound material by growing strained quantum well structures using a different composition of the hexagonal Si.sub.1−xGe.sub.x compound material in each layer of the quantum well structure.
[0130]
[0131] According to a method step according to the invention, the Group 11 element catalyst seeds, here Au catalyst seeds were deposited in nano disks arrays arrangement on a Group III-V compound semiconductor substrate, here a GaAs (111).sub.B substrate via an electron beam lithography technique. The growth was performed at a reactor flow of 8.2 standard litres per minute (slm) and a reactor pressure of 50 mbar. For the GaAs nanowires, the growth template was annealed at a set thermocouple temperature of 635° C. under an AsH.sub.3 flow set to a molar fraction of χ.sub.AsH3=6.1×10.sup.−3 mols. Then, the growth was performed at a set temperature of 650° C. with trimethylgallium (TMGa) and Arsine (AsH.sub.3) as material precursors set to molar fractions of χ.sub.TMGa=1.9×10.sup.−5 mol, χ.sub.AsH3=4.55×10.sup.−5 mol, respectively, resulting in a V/III ratio of 2.4.
[0132] After the growth of the GaAs core nanowires, in a further method step the GaAs core nanowires are chemically treated with a cyanide based solution to remove the Au catalyst particles to avoid gold contamination in the SiGe shells, (see,
[0133]
[0134] As to the crystal quality of the WZ GaAs nanowire Cores grown with an example of a method according the invention,
[0135]
[0136]
[0137]
[0138] Finally,
[0139] For the APT measurements, individual nanowires (nanowires) were isolated from a forest of nanowires as described previously with a Kleindiek nano-manipulator inside a FEI Nova Nanolab 600i Dual beam. APT analyses were carried out in a LEAP 4000X-HR from Cameca. The system is equipped with a laser generating picosecond pulses at a wavelength of 355 nm. The experimental data were collected at laser or voltage pulse rates between 65-125 kHz with laser pulse energies between 5-10 pJ or pulse fractions between 25-27.5%. No significant differences between laser and voltage pulses are seen aside from a slightly higher compression of the core in laser pulsed mode and a lower quality of the mass spectra in voltage pulsed mode. During the analysis the sample is kept at a base temperature of 20 K in a vacuum of ˜2.Math.10.sup.−11 mbar.
[0140]
[0141]
[0142]
[0143] The calculated band structure of hexagonal Ge is shown in
[0144] The emission blue shifts with increasing Si concentration. The experimentally obtained emission energies are included in
[0145] An aspect of the invention an optical device is characterized by hexagonal Si.sub.1−xGe.sub.x compounds with 0.6<x<0.9. This material compound is important for integration of photonic functionalities in the Silicon industry. This material can be used to fabricate light-emitting diodes (LEDs), lasers, and detectors, and can be integrated in Si technology or can be used as active device in passive optical circuitry. Application of such devices is in logic chips, telecommunication, chemical sensing, IR imaging etc.
[0146] The efficient light emission has been experimentally verified by performing time-resolved photoluminescence lifetime measurements, yielding photoluminescence lifetimes of the order of 0.5-0.8 ns for an detailed aspect of an optical device according to the invention being characterized by Si.sub.0.78Ge.sub.0.22.
[0147] Further aspects of the invention pertain to optoelectronic applications and/or optoelectronic products comprising an light emitting body or component comprised of hexagonal Si.sub.1−xGe.sub.x compounds, which more in particular exhibit a direct band gap at Brillouin zone center for 0.6<x<1.0.
[0148] Additionally according to a further aspect of the invention, the optical emitting body or component comprised of hexagonal Si.sub.1−xGe.sub.x compounds may exhibit large matrix elements for hexagonal Si.sub.1−xGe.sub.x compounds with 0.2<x<1.0.
[0149] In particular optical emitting component comprised of hexagonal Si.sub.1−xGe.sub.x compounds with the particular range 0.6<x<0.8 exhibit large oscillator strengths and a direct fundamental gap 0.6-0.8 eV, which corresponds to emission wave length 1.5-2.0 μm.
[0150] The method or process according to the invention allows fabricating or manufacturing hexagonal Si.sub.1−xGe.sub.x crystals throughout the whole range of stoichiometry of the alloy. The hexagonal crystal phase of the Si.sub.1−xGe.sub.x crystal is achieved via the “Crystal Transfer” technique where the hexagonal crystal structure of the Si.sub.1−xGe.sub.x is adopted from a hexagonal material template. According to an aspect of the method according to the invention, this is achieved by utilizing non-tapered and single-crystalline wurtzite (hexagonal) GaAs nanowires (nanowires) that act as a template and then overgrow them with an epitaxial layer of Si.sub.1−xGe.sub.x using Metal Organic Vapor Phase Epitaxy (MOVPE) technique.
[0151] An experimental observation of efficient light emission for an example of an optoelectronic compound according to the invention comprised of Si.sub.0.78Ge.sub.0.22 as evidenced by the 0.7 ns photoluminescence (PL) lifetime is shown in the graph of
[0152] In an aspect of the method according to the invention the GaAs nanowires are grown by the Vapor-Liquid-Solid (VLS) mechanism utilizing gold (Au) catalyst seeds deposited in nano disks arrays arrangement on a GaAs (111).sub.B substrate via electron beam lithography. The growth template in
[0153] Eventually, the GaAs nanowire core is overgrown with the Si.sub.1−xGe.sub.x shell by introducing the suitable gas precursors of the shell growth (Germane (G.sub.eH.sub.4) and Disilane (Si.sub.2H.sub.6)) in the reactor, see
[0154] In an example of the method according to the invention the epitaxy process was performed at a reactor flow of 8.2 standard litres per minute (slm). The growth template was annealed at a surface temperature of 526° C. monitored by argus under A.sub.sH.sub.3 flow set to a molar fraction of χ.sub.AsH3=6.1×10.sup.−3. Then, the growth was performed at a surface temperature of 546° C. with TMGa, and A.sub.sH.sub.3 as material precursors set to molar fractions of χ.sub.TMGa=1.9×10.sup.−5, χ.sub.AsH3=4.55×10.sup.−5, respectively, achieving a V/III ratio of 2.4. The Si.sub.1−xGe.sub.x crystals were grown at a surface temperature of 560° C. at a molar fraction of χ.sub.SiGe=1.55×10.sup.−4.
[0155] In
[0156]
[0157] For a better understanding of the vapor-solid shell growth mechanism as well as for future optical characterization experiments it is important to assess the homogeneity of the Ge distribution along the entire shell. Therefor the Ge content is determined in the SiGe shell of the samples grown. To quantify the Ge content, the samples were prepared by focused ion beam (FIB) for cross-sectional TEM studies. The cross-sectional bright-field (BF) TEM image of , ⋄, Δ and ∇ in
[0158]
[0159] To further reveal the SiGe shell growth mechanism, it is useful to plot the resulting incorporated Ge content in the shell, as determined by EDX spectroscopy, as a function of the Ge content in the precursor gas flow in the MOVPE chamber, for all the heterostructures investigated. This is shown in
[0160] A characteristic feature of the cubic Si/SiGe system is the appearance of different growth modes as the Ge content; consequently, the lattice strain increases in the SiGe shell and it is observed for the hexagonal phase as well. For most of the samples (for a Ge content of less than 80 atom % grown at 600° C. and of less than 60 atom % grown at 700° C., indicated by solid symbols in
[0161] The growth kinetics have been studied in more detail and especially the growth rate as a function of time and temperature. The linear dependence of the SiGe thickness with respect to the growth time indicates that the growth rate is roughly constant (
[0162] In addition, the study of the growth rate was performed with a series of samples grown for the whole range of SiGe stoichiometries and for the two different temperatures of 600 and 700° C. In
[0163] An aspect of the method according to the invention pertains to the growth of single-crystalline defect-free SiGe with or on a hexagonal diamond crystal structure. According to the invention the growth temperature leads to significant differences in the resulting SiGe layer morphology: SiGe layers grown at lower temperatures exhibit uniform layer-by-layer F-M growth, while at higher temperatures and for the same Ge content, island-based S-K growth dominates. The maximum stoichiometry achieved for defect-free smooth layer thin film growth was at 77 atom % Ge, significantly above the predicted direct band gap transition at 65 atom % Ge. Moreover, it is observed that layer growth continues to this high value significantly more at the growth temperature of 600° C. Therefor this allows for the utilization of hexagonal SiGe in optoelectronic applications.
[0164] Examples of the GaP/Si/SiGe core/multishell nanowires according the invention were developed in a low-pressure (50 mbar) Aixtron Close Coupled Showerhead (CCS) MOVPE reactor. In this case though, the Si shell thickness was kept to a minimum of 10-12 nm. This could potentially act as sacrificial buffer layer to trap any P or Ga species from the GaP core from diffusing into the SiGe shell.
[0165] Directly after the growth of the Si shell and in the same growth run, the SiGe shell was grown. This was done by switching off the Si.sub.2H.sub.6 precursor gas and gradually lowering the temperature of the MOVPE growth chamber to either 600 or 700° C. Once the aimed temperature was reached and stabilized, Si.sub.2H.sub.6 and GeH.sub.4 were introduced to the chamber and their flows were adjusted to suit the desired stoichiometry. Si.sub.2H.sub.6 molar flow was modified from 2.87×10.sup.−7 at the low end to 1.00×10.sup.−4 at the high end, whereas GeH.sub.4 was modified from 3.66×10.sup.−6 at the low end to 3.33×10.sup.−4 at the high end in order to achieve the desired ratio. Hydrogen (H.sub.2) was used as a carrier gas for the precursors and the total flow into the rector was 8.2 L/min. At the end of the SiGe shell growth, the precursor flows were terminated, the heating elements were switched off and the chamber was allowed to cool down to room temperature.
[0166] For the Transmission Electron Microscopy, TEM, studies, two different sample preparation methods were used. In the standard axial analysis, nanowires were mechanically transferred to a holey carbon TEM grid. Concerning the cross-section TEM studies, nanowires were prepared using FIB. In both cases, HRTEM and STEM analyses were conducted using a JEM ARM200F aberration-corrected TEM operated at 200kV. For the chemical analysis, EDX measurements were carried out using the same microscope equipped with a 100 mm.sup.2 EDX silicon drift detector.
[0167] TEM Lamellae. TEM lamellae were prepared in a FEI Nova Nanolab 600i Dualbeam system. For this, the nanowires were initially transferred with the aid of a nanomanipulator from the growth substrate to a piece of Si and then arranged to lie parallel to each other. These nanowires were covered with electron- and ion-beam induced metal deposition to protect them during the procedure. The lamella was cut out by milling with 30 kV Ga ions and thinned down with subsequent steps of 30, 16, and 5 kV ion milling in order to minimize the Ga-induced damage in the regions imaged with TEM.
[0168] In another example according to the invention, high-quality defect-free lonsdaleite Si and Ge can be grown on hexagonal nanowire substrates. These hexagonal phases of group-IV semiconductors exhibit improved electronic and optical properties for optoelectronic applications. While lonsdaleite Si is a well-characterized indirect semiconductor, experimental data and reliable calculations on lonsdaleite Ge are scarce and not consistent regarding the nature of its gap. Using ab initio density-functional theory described above, accurate structural, electronic, and optical properties for hexagonal Ge can be determined. Given the well-known sensitivity of electronic-structure calculations for Ge to the underlying approximations, the performance of several exchange-correlation functionals, including meta-GGA and hybrid functionals can be tested. It is validated for cubic Ge, obtaining atomic geometries and band structures in agreement with available experimental data. The same approach is applied to predict electronic and optical properties of lonsdaleite Ge. According to a further aspect of the invention, an optoelectionic device comprised of lonsdaleite Ge as a direct semiconductor with only weakly dipole-active lowest optical transitions, small band gap, huge crystal-field splitting, and strongly anisotropic effective masses. The unexpectedly small direct gap and the oscillator strengths of the lowest optical transitions are explained in terms of symmetry and back-folding of energy bands of the diamond structure.
[0169] The integration of a material featuring efficient interaction with light into Si technology is of high technological interest. In fact, the copper interconnect between transistors on a chip has become a bigger challenge than reducing transistor size. A possible solution to this critical bottleneck are optical interconnects. Si in the diamond structure is an indirect band-gap material and cannot be used for this purpose. Several attempts to obtain light emission from Si in the telecommunication band had only limited success.
[0170] Si and Ge, despite their chemical similarity, are very different from an optical perspective. Ge, one of the most important and widely used semiconductors, crystallizes like Si in the cubic diamond structure (space group Fd
[0171] Besides the thermodynamically stable diamond structure, other metastable allotropes of Ge have been explored for optical applications. Ge in the hexagonal lonsdaleite structure (space group P63/mmc), occasionally also called wurtzite Ge, is attracting increasing attention as a promising material for optoelectronics. In the lonsdaleite phase, the Ge atoms feature the same tetrahedral nearest-neighbor coordination as in the cubic diamond structure, but, instead of an ABC stacking of adjacent Ge bilayers along the threefold symmetry axis, the lonsdaleite phase is characterized by an AB stacking. Therefore, reference is made to cubic Ge in the diamond structure also as 3C-Ge and to hexagonal Ge in the lonsdaleite structure as 2H-Ge (see
[0172] In the hexagonal Brillouin zone (BZ) of the lonsdaleite structure, the L point of the diamond-structure BZ that lies on the cubic [111] axis is mapped onto the G point. Therefore, the lowest conduction-band minimum (CBM) at the L point of cubic Ge is folded onto the G point rendering 2H-Ge a direct-gap semiconductor. Comparing hexagonal Si and Ge, it seems that breaking the k-selection rule is easier in Ge, since theoriginal and the backfolded conduction band are energetically very close. The exact ordering of the lowest conduction bands at G is extremely important, as the electron radiative lifetime of the material strongly depends on the symmetry of these states.
[0173] Together with optical emission or absorption measurements for photon energies comparable with the size of the band gap, accurate electronic-structure calculations can provide detailed answers concerning conduction-band ordering and the strength of optical transitions. Despite the fact that Ge is an elemental material, the experience acquired with calculations of diamond-structure Ge proves that Ge is a difficult system for accurate band-structure studies. On the one hand, it is essential to account for spin-orbit coupling (SOC) and to treat the shallow Ge 3d shell as valence electrons. On the other hand, the approximations used to describe exchange and correlation (XC) contributions to the electron-electron interaction significantly influence the k-space position of the lowest conduction-band minimum and the size of the direct and indirect gaps.
[0174] Applying density-functional theory (DFT), a Kohn-Sham (KS) band structure obtained within the local-density approximation (LDA) or any flavor of the generalized gradient approximation (GGA) is not sufficient, as the system is erroneously predicted to be metallic in this case. Moreover, the band gap is extremely sensitive to the value of the lattice constant. Therefore, an accurate description of the atomic geometry is indispensable. More sophisticated approaches, beyond semi-local XC functionals, are needed to obtain reliable quasiparticle states.
[0175] In fact, 3C-Ge band structures have also been computed approximating the XC self-energy within Hedin's GW approximation, and it is expected this state-of-the-art approach for excited states to work also for lonsdaleite Ge. However, the computational cost coming along with Green's function calculations is very high and would likely make their application to more complex systems (such as alloys, doped or defective crystals, surfaces, or interfaces) unfeasible. It is remarked that also the empirical-pseudopotential method (EPM), widely used for 3C-Ge, can be helpful at a reduced computational cost. However, empirical approaches that very accurately reproduce experimental data for 3C-Ge would need additional assumptions to be reliably applied to 2H-Ge.
[0176] A careful analysis of the electronic and optical properties of 2H-Ge is presented, with a particular focus on the choice of accurate and computationally efficient XC functionals for ground-state and excited-state calculations. The functionals are first tested against experimental data for 3CGe. They are then used for a careful analysis of the electronic and optical properties of 2H-Ge. In view of potential optoelectronic applications, special interest is given to conduction-band ordering, direct and indirect band gaps, band splittings, effective masses, optical transition strengths, and radiative lifetimes.
[0177] All calculations were performed with the Vienna Abinitio Simulation Package (VASP) with the projector-augmented wave (PAVV) method and a plane-wave cutoff of 500 eV. The shallow Ge 3d electrons were explicitly included as valence electrons. BZ integrations were carried out using 12×12×12 (3C-Ge) or 12×12×6 (2H-Ge) r-centered k points (unless otherwise stated), ensuring a convergence of total energies to 1 meV/atom. Atomic geometries and elastic properties were calculated with (semi-)local XC functionals, using the LDA as well as the GGA parametrizations PBE, PBEsol (a modified version of the PBE functional optimized for solids), and AM05.
[0178] The ground-state atomic structures, the isothermal bulk modulus B.sub.0, and its pressure derivative B′.sub.0 were determined by a series of fixed-volume relaxations and a subsequent fit of the resulting energy-over-volume curve to the Vinet equation of state (EOS). The internal cell parameters were relaxed until the Hellmann-Feynman forces drop below 1 meV/Å. It was found that the inclusion of SOC has essentially no impact on the lattice parameters and only a minor effect on the elastic constants. This observation is in line with general conclusions for other simple solids and zincblende-type semiconductors.
[0179] It is well known that KS band structures calculated in the LDA or GGA significantly underestimate all band gaps and interband transition energies. Quasiparticle calculations in the state-of-the-art GW approximation, on the other hand, are challenging and computer-time consuming for Ge, due to the necessity to include SOC, to account for the 3d electrons, and to calculate the full dynamical screening.
[0180] What is more, the GW quasiparticle band structures need to be computed self-consistently to overcome the problem of the negative fundamental gap in the LDA/GGA starting electronic structure of both 3C-Ge and 2H-Ge (see Sections III and IV). One reason for the negative gaps is the overestimation of the p-d repulsion. This is a direct consequence of the underestimated binding energy of the Ge 3d electrons within LDA or GGA, which pushes the p-like valence-band maximum (VBM) towards higher energies. An improved description of the localized d states can be achieved within the DFT+U method with a Hubbard parameter U for the 3d electrons which comes at the price of introducing a tunable parameter. The DFT+U method was tested in the Dudarev approach using a small but reasonable value U=1.3 eV, which reproduces the lattice constant of 3C-Ge and is in rough agreement with the picture of an atomic Coulomb integral of about U.sup.atom=15 eV screened by the bulk Ge dielectric constant.
[0181] The HSE06 hybrid functional was used with a fraction a=0.25 of short-range Fock exchange and an inverse screening length μ=0.2 Å.sup.−1 to calculate reliable band structures for cubic and hexagonal Ge. It has been shown that the HSE06 functional yields reasonable indirect and direct gaps for Ge and many other sp semiconductors. The most important contribution to the gap opening within the GW approach is due to the screened-exchange part of the electronic self-energy. The Coulomb hole, the second contribution to the GW self-energy, mainly influences the absolute position of the one-particle energies. In the HSE06 functional, the fraction a of Fock exchange simulates the important non-locality feature of the self-energy and the screening of the electron-electron interaction by an average dielectric constant of 1/a.
[0182] As a computationally cheap alternative to hybrid functionals, we also consider the meta-GGA functional MBJLDA of Tran and Blaha that is based on the modified Becke-Johnson (MBJ) exchange functional. The MBJLDA functional does not only give reasonable band gaps for 3C-Ge but also for other semiconductors. The strongly reduced computational cost allows for the application of the MBJLDA functional to more complex systems. In particular in the context of potential optoelectronic applications of 2H-Ge, also strained, disordered, or defective systems with larger supercells become computationally accessible. Moreover, both the hybrid and the meta-GGA functional allow for an easy inclusion of SOC.
[0183] Having in mind optoelectronic applications (e.g. lasing), the global optical emission properties of 2 H-Ge near the fundamental absorption edge can be characterized by the optical transition matrix elements of the near-edge transitions and the radiative lifetime of the material. Here, the optical transition matrix elements are calculated in the longitudinal gauge. They are given as matrix elements <ck|p|vk> of the momentum operator p between conduction band c and valence band v at a given k point.
[0184] The optical matrix elements at the G point can be linked to characteristic quantities from k:p perturbation theory introducing the average of the squared momentum matrix element over spin-orbit degenerate states i, j=1,2 in the conduction and valence bands at the zone center,
Then, the Kane energy reads
and the (dimensionless) optical oscillator strength
where ⊥/∥ stands for light polarized perpendicular/parallel to the c axis of the lonsdaleite structure. For 3C-Ge, these two directions are equivalent due to the isotropy of the material.
[0185] The radiative lifetime τ at temperature T, as a global measure for the light-emission properties of a material, is given by the thermally averaged recombination rate:
where A.sub.cvk denotes the radiative recombination rate for vertical optical transitions between a conduction state |ck> and a valence state |vk> with the one-particle energies ε.sub.ck and ε.sub.vk and k-point weight w.sub.k. The radiative recombination rate reads
with n.sub.eff the refractive index of the effective medium consisting of the Ge sample and its environment (set to 1 in the following). The squares of the momentum matrix elements are averaged over all directions corresponding to the emission of unpolarized light. Eq. (5) is given in the independent-(quasi)particle approximation, i.e., neglecting excitonic effects, which, however, can be easily taken into account. In Eq. (4), it is assumed that the thermalization of electrons and holes after their injection is more efficient than the radiative (or nonradiative) recombination. Whereas the convergence of the radiative lifetimes with the number of bands is very fast, we need 72×72×72 (3C-Ge) or 72×72×36 (2HGe) k points to sample the BZ with sufficiently high accuracy.
[0186] As to the diamond-structure of germanium, the lattice constant, elastic properties, and cohesive energy of 3C-Ge have been calculated with various XC correlation functionals (see Table S4). Comparing with experimental values, the expected tendencies are visible: The LDA tends to overbind, whereas the inclusion of gradient corrections, in particular within PBE, leads to an underestimation of the strength of the chemical bonds. The functionals PBEsol and AM05 yield the best agreement with experiment. However, they still slightly overestimate the experimental lattice constant. Further improvement can be obtained by DFT+U calculations (see Table S4) at the price of introducing the adjustable parameter U. Also the isothermal bulk modulus B.sub.0, its pressure derivative B′.sub.0, and the cohesive energy are consistent with experiment.
[0187] Subsequently, the band structure of 3C-Ge including SOC was calculated using the PBEsol, HSE06, and MBJLDA functionals at the PBEsol lattice constant (see
TABLE-US-00004 TABLE S4 Lattice constant a.sub.0, isothermal bulk modulus B.sub.0, its pressure derivative B′.sub.0, and cohesive energy E.sub.coh of 3C—Ge. Experimental (Exp.) values are given for comparison. Method a.sub.0 (Å) B.sub.0 (GPa) B′.sub.0 E.sub.coh (eV/at.) LDA 5.626 71.8 4.92 4.63 PBE 5.760 58.7 5.01 3.73 PBEsol 5.673 67.3 4.89 4.15 AM05 5.677 65.7 4.82 3.92 PBEsol + U 5.652 68.7 4.86 4.15 AM05 + U 5.656 67.2 4.87 3.91 Exp. 5.658.sup.a 75.0.sup.b 3.85.sup.c 5.652.sup.d 64.7.sup.e 5.0(1).sup.e 77(4).sup.f 4.3(1.0).sup.f .sup.aX-ray diffraction at T = 298.15 K [65]. .sup.bObtained from ultrasonic measurements of the elastic moduli C.sub.11 and C.sub.12 at ambient pressure at T = 298.15 K [66] using the relation B.sub.0 = (C.sub.11 + 2C.sub.12)/3. .sup.cFrom Ref. [67]. .sup.dX-ray diffraction at T = 10 K [68]. .sup.eFrom fitting an EOS to room-temperature experimental data for various pressures [69]. .sup.fFrom fitting the Vinet EOS to room-temperature experimental data for various pressures [70].
[0188] The states at the high-symmetry points of the BZ are labeled according to the double-group notation. The small difference between the PBEsol and the experimental lattice constant corresponds to an isotropic tensile strain of <0.4%. The volume deformation potential for the direct gap (the most volume-sensitive band-to-band transition) amounts to −9.0 eV (MBJLDA) implying that differences in the direct gap due to the discrepancy in the lattice constant are smaller than 0.1 eV. For the sake of comparability, all calculations of electronic and optical properties presented in the following are based on the PBEsol lattice constant.
[0189] GGA functionals like PBEsol yield a negative KS band gap for 3C-Ge in contradiction to experimental findings (cf. Table S5) which is why they are unsuitable for the description of the electronic structure of this material. The band ordering, band energies, and spin-orbit splittings Dso obtained with the more sophisticated HSE06 functional agree well with experimental results. Comparing HSE06 and MBJLDA band structures close to the fundamental gap, we find similar indirect (Γ.sub.8v.sup.+.fwdarw.L.sub.6c.sup.+) and direct (Γ.sub.8v.sup.+.fwdarw.Γ.sub.7c.sup.−) band gaps. Also the spin-orbit splittings of the p states are much the same. Further away from the band-gap region, the discrepancy between HSE06 and MBJLDA band energies increases. This is, however, not a crucial problem here, since we are mostly interested in optoelectronic properties that are governed by the electronic structure in the vicinity of the band gap. In particular, the ordering of the Γ.sub.7c.sup.− and L.sub.6c.sup.+ conduction-band minima is correct, only their energy distance is slightly underestimated compared to experiment (independent of temperature). Note that GW corrections on top of HSE06 or MBJLDA band structures are known to overestimate the gaps.
[0190] In Table S5, the electron and hole effective masses of relevant band extrema are compiled. Besides the band masses at the G point, also the masses of the conduction-band minimum at the L point parallel and perpendicular to the L-G line are given. The masses have been derived from the corresponding HSE06 and MBJLDA band structures. The HSE06 masses are in excellent agreement with experimental values. The MBJLDA functional slightly overestimates the experimental band masses which is in line with previous observations and the generally lower band widths in the MBJLDA calculation compared to the HSE06 calculation (cf.
[0191] Considering the findings for 3C-Ge, we rely on the PBEsol functional for the structural properties of 2H-Ge. The HSE06 and MBJLDA functionals will be used to study the electronic and optical properties of lonsdaleite Ge. This strategy is corroborated by the fact that both allotropes of Ge feature similar chemical bonding properties, i.e. they are both insulators with tetrahedral coordination. Therefore, the performance of the functionals should be largely transferable.
TABLE-US-00005 TABLE S5 Band energies and spin-orbit induced band splittings Δ.sub.so of 3C—Ge at high-symmetry points of the BZ calculated with the PBEsol, HSE06, and MBJLDA functionals at the PBEsol lattice constant. Experimental low- and room-temperature values are provided for comparison. All values in eV. Experiment State PBEsol HSE06 MBJLDA low temperature room temperature Γ.sub.7v.sup.+ −0.292 −0.314 −0.270 −0.297.sup.a Γ.sub.8v.sup.+ 0.000 0.000 0.000 Γ.sub.7c.sup.− −0.174 0.678 0.705 0.887(1).sup.a, 0.898(1).sup.b 0.805(1).sup.b Γ.sub.6c.sup.− 2.296 2.969 2.564 Γ.sub.8c.sup.− 2.508 3.189 2.760 Γ.sub.6c.sup.− + 0.200.sup.a Δ.sub.so(Γ.sub.8v.sup.+ − Γ.sub.7v.sup.+) 0.292 0.341 0.270 0.297.sup.a Δ.sub.so(Γ.sub.8c.sup.−− Γ.sub.6c.sup.−) 0.212 0.220 0.196 0.200.sup.a L.sub.6v.sup.− −1.585 −1.719 −1.530 L.sub.4v.sup.− + L.sub.5v.sup.− −1.401 −1.519 −1.359 L.sub.6v.sup.− + 0.228.sup.a L.sub.6c.sup.+ 0.005 0.675 0.653 0.744(1).sup.b 0.6643.sup.c Δ.sub.so(L.sub.4v+5v.sup.− − L.sub.6v.sup.−) 0.184 0.200 0.171 0.228.sup.a X.sub.5v −3.155 −3.441 −3.047 X.sub.5c 0.585 1.177 1.142 .sup.aSchottky-barrier electroreflectance at 10 K [72]. .sup.bMagnetoabsorption at 1.5 K and 293 K [73]. .sup.cOptical absorption-edge fine structure at 291 K [74].
TABLE-US-00006 TABLE S6 Effective electron and hole masses of 3C—Ge in units of the free electron mass m. The VBM at Γ splits into a heavy hole (m.sub.h.sup.hh), light hole (m.sub.h.sup.lh), and spin-orbit split-off hole (m.sub.h.sup.so). The heavy-holeand light-hole masses are averaged over the Γ-X and Γ-L directions. The masses of the conduction-band minimum at L are given both parallel (m.sub.e.sup.∥) and perpendicular (m.sub.e.sup.⊥) to the L-Γ direction. Mass HSE06 MBJLDA Exp. m.sub.h.sup.so(Γ.sub.7v.sup.+) 0.097 0.122 0.095(7).sup.a m.sub.h.sup.lh(Γ.sub.8v.sup.+) 0.043 0.059 0.0438(30) B∥[100].sup.b 0.0426(20) B∥[111].sup.b 0.0430(30) B∥[110].sup.b m.sub.h.sup.hh(Γ.sub.8v.sup.+) 0.203 0.233 0.284(1) B∥[100].sup.b 0.376(1) B∥[111].sup.b 0.352(4) B∥[110].sup.b m.sub.e(Γ.sub.7c.sup.−) 0.034 0.047 0.0380(5).sup.a m.sub.e.sup.∥(L.sub.6c.sup.+) 1.573 1.728 1.588(5).sup.c, 1.59.sup.d m.sub.e.sup.⊥(L.sub.6c.sup.+) 0.090 0.096 0.08152(8).sup.c, 0.0823.sup.d .sup.aPiezomagnetoreflectance at 30 K [75]. .sup.bCyclotron resonance at 4 K with magnetic field B oriented in various directions [76]. .sup.cCyclotron resonance at 1.4 K [77]. .sup.dMagnetophonon resonance at 120 K [78].
[0192] As to the atomic geometry and bonding lonsdaleite Germanium, the positions of the four atoms in the unit cell of the lonsdaleite structure are defined by the hexagonal lattice constants a and c, as well as the dimensionless internal cell parameter u. In Table S7, the results of the calculations of the structural properties of 2H-Ge with various XC functionals are compiled. As already discussed for 3C-Ge, a consistent over- and underestimation of bond lengths depending on the choice of the functional is observed. These tendencies are also reflected in the values for the cohesive energy E.sub.coh and the bulk modulus B.sub.0.
[0193] On average, the bond lengths of 2H-Ge are slightly reduced in comparison to the 3C-Ge bond length d=√{square root over (3)}/4a.sub.0. This can be understood in terms of a detailed analysis of the atomic geometry. The Ge-Ge bonds parallel to the c axis (bond length d.sub.∥=μc and those not parallel to the c axis (bond length d.sub.⊥=√{square root over (.sup.a.sup.
[0194] It can be deduced from Table S7 that c/a>(c/a).sub.ideal and u<u.sub.deal. The relation d.sub.⊥<d<d.sub.∥ holds, resulting in tetrahedral that are slightly elongated along the c axis. For instance, we find d.sub.⊥=2.451 Å and d.sub.∥=2.468 Å for the PBEsol functional. The average bond length d.sub.av=2.455 Å of 2H-Ge is only slightly smaller than the 3C-Ge bond length d=2.457 Å. These findings are in line with the empirical rule of Lawaetz for III-V semiconductor compounds, which states that materials with c/a>(c/a).sub.ideal favor a zincblende ground-state structure, or for an elemental materials as Ge, the diamond structure. The computed values of u nearly follow the relation u=1/4+1/3(c/a).sup.−2 indicating that the deformation of bonding tetrahedra in 2H-Ge can be explained to a good share by deviations of the bond angles from the ideal value.
[0195] In summary, a relatively strong hexagonal crystal deformation in 2H-Ge is observed, which is characterized by large (c/a>(c/a).sub.ideal and u-u.sub.ideal, despite the presence of covalent bonds. The calculated lattice parameters of 2H-Ge agree very well with the available experimental data. However, experimental structural parameters are scarce in the literature and, in some cases, have been obtained from nanostructured and potentially strained samples.
TABLE-US-00007 TABLE S7 Structural and elastic properties of 2H—Ge. Structural parameters a, c, and u as wells as the isothermal bulk modulus B.sub.0, its pressure derivative B′.sub.0, and cohesive energy E.sub.coh. Available experimental data are given for comparison. B.sub.0 E.sub.coh Method a (Å) c (Å) c/a u (GPa) B′.sub.0 (eV/at.) LDA 3.962 6.539 1.6504 0.3742 71.9 5.00 4.61 PBE 4.058 6.692 1.6492 0.3744 59.1 4.74 3.71 PBEsol 3.996 6.590 1.6492 0.3744 67.6 4.81 4.14 AM05 3.999 6.594 1.6490 0.3745 66.0 4.95 3.91 PBEsol + 3.980 6.568 1.6503 0.3743 69.0 4.80 4.13 U Am05 + 3.982 6.568 1.6503 0.3743 67.4 4.81 3.89 U Exp 3.96.sup.a 6.57.sup.a 1.659.sup.a 3.9878 6.5776 .sub. 1.6494.sup.b (20).sup.b (3).sup.b .sup.aPotentially strained samples from micro-indentation [79]. .sup.bRoom-temperature x-ray diffraction of unstrained crystalline nanowires [80].
TABLE-US-00008 TABLE S8 Band energies of 2H—Ge computed with different XC functionals. The Γ.sub.9v.sup.+ is used as energy zero. The crystal-field and spin-orbit splitting parameters Δ.sub.cf, Δ.sub.so.sup.∥, and Δ.sub.so.sup.⊥ for the VBM have been calculated from the band energies. Band energies from an empirical-pseudopotential model (EPM) are given for comparison. All energies in eV. State HSE06 MBJLDA EPM (Ref. [33]) Γ.sub.7v−.sup.+ −0.484 −0.433 −0.490 Γ.sub.7v+.sup.+ −0.134 −0.120 −0.129 Γ.sub.9v.sup.+ 0.000 0.000 0.000 Γ.sub.8c.sup.− 0.286 0.298 0.310 Γ.sub.7c.sup.− 0.614 0.632 0.766 U.sub.5c 0.615 0.620 Δ.sub.cf 0.288 0.270 Δ.sub.so.sup.∥ 0.329 0.282 Δ.sub.so.sup.⊥ 0.320 0.274
[0196] The band structure of 2H-Ge including SOC has been calculated with the HSE06 and MBJLDA functionals for the PBEsol atomic structure (see
[0197] 2H-Ge is found to be a direct-gap semiconductor with a band gap of 0.286 eV (HSE06) or 0.298 eV (MBJLDA). Note that the precise gap values are highly sensitive to lattice strain due to the huge deformation potentials of the gap-forming states. The results agree well with previously calculated gaps of 0.31 eV (empirical-pseudopotential method) or 0.32 eV (HSE06 calculation), being slightly higher than a GW gap of 0.23 eV obtained by the prior art. As for 3CGe, the HSE06 and MBJLDA energies of the near-gap states match excellently. Deviations occur for states further away from the band-gap region (see
[0198] The conduction-band minima of 3C-Ge are the four L.sub.6c.sup.+ states. The L.sub.6c.sup.+state in [111] direction is backfolded to the Γ.sub.8c.sup.− state in 2H-Ge and becomes the conduction-band minimum of 2H-Ge. The other L points of 3C-Ge are mapped onto a point between M and L on the U line of the hexagonal BZ. In the ideal lonsdaleite structure, the backfolded L points lie at ⅔
[0199] The Γ.sub.8c.sup.− state in 2H-Ge is downshifted by 0.4 eV compared to the cubic L.sub.6c.sup.+ state which cannot be understood by simple folding arguments. It is, however, in agreement with the behavior of Si going from diamond to lonsdaleite structure, whereas it is in clear contrast with the small band-gap opening in biatomic semiconductors when the structure changes from zincblende to wurtzite. The cubic Γ.sub.7c.sup.− state coincides with the Γ.sub.7c.sup.− state of the second conduction band in the lonsdaleite structure.
[0200] Due to the presence of inversion symmetry in the lonsdaleite structure, the valence bands in 2H-Ge do not show a spin-orbit-induced splitting of the k dispersion along the ┌-M line as it occurs in wurtzite semiconductors. Following k p theory, we can write the energy splittings at ┌ as
[0201] These formulas allow to extract the crystal-field splitting Δ.sub.cf and the spin-orbit splitting parameters parallel and perpendicular to the c axis, Δ.sub.so.sup.∥, and Δ.sub.so.sup.⊥. The band ordering Γ.sub.9c.sup.+>Γ.sub.7v.sup.+>Γ.sub.7v.sup.− at the top of the valence bands which we find for 2H-Ge is in line with the ordering observed in wurtzite semiconductors (except AlN and ZnO for which the band ordering is Γ.sub.7v.sup.+>Γ.sub.9 v>σ.sub.7v−). Note that the subscript indices 7v± represent bands of the same symmetry; the symbols ±merely serve to distinguish between the upper and the lower state. They are not to be confused with superscript parity indices.
[0202] The crystal-field splitting has been extracted from a calculation without SOC and used Eq. (6) to compute the spin-orbit splitting parameters from the band splittings of the calculation including SOC. The resulting values are compiled in Table S8. In particular, the direction-averaged spin-orbit splitting Δ.sub.so=(Δ.sub.so.sup.∥+2Δ.sub.so.sup.⊥)/3 compares well to the spin-orbit splitting of the VBM in 3C-Ge (cf. Table S5). The crystal-field splitting in 2H-Ge is much larger than for III-V compounds that crystallize in zincblende or wurtzite structure under ambient conditions. The large crystal-field splitting for 2H-Ge is, however, in accordance with the significant deformation of the bonding tetrahedra, as indicated by the increase of c/a (see Table S7) with respect to its ideal value. The large Δ.sub.cf shifts the Γ.sub.9v.sup.+ level toward higher energies and, hence, explains the observed small direct gap. It is emphasized that the quasicubic approximation Δ.sub.so.sup.∥=Δso.sup.⊥ that was used in the prior art is not valid for 2H-Ge and leads to splitting parameters at variance with the values.
[0203] The effective masses of the band edges at ┌ are compiled in Table S9. The small electron mass with almost vanishing anisotropy for the Γ.sub.7c.sup.−conduction band of 2H-Ge is of the order of magnitude of the Γ.sub.7c.sup.−mass of 3C-Ge in Table S6. The mass tensor at the conduction-band minimum Γ.sub.8c.sup.−, on the other hand, is highly anisotropic with a large mass along the hexagonal c axis and a small mass in the plane perpendicular to it. These values qualitatively agree with the longitudinal and transverse masses m.sub.e.sup.∥(L.sub.6c.sup.+) and m.sub.e.sup.⊥(L.sub.6c.sup.+) at the L.sub.6c.sup.+minimum of 3C-Ge. The strong direction dependence of the Γ.sub.8c.sup.− conduction-band dispersion is consistent with the identification of the band symmetry. The hole masses also exhibit strong asymmetries, especially for the Γ.sub.9v.sup.+ and Γ.sub.7v.sup.+ bands.
TABLE-US-00009 TABLE S9 Effective electron and hole masses of 2H—Ge in units of the free electron mass m. The masses are given for several directions in the BZ. The VBM at G splits into a heavy hole (m.sub.h.sup.hh), light hole (m.sub.h.sup.lh), and split-off hole (m.sub.h.sup.so). Mass Direction HSE06 MBJLDA m.sub.h.sup.so(Γ.sub.7v−.sup.+ Γ .fwdarw. M 0.252 0.325 Γ .fwdarw. A 0.044 0.053 m.sub.h.sup.lh(Γ.sub.7v+.sup.+) Γ .fwdarw. M 0.079 0.101 Γ .fwdarw. A 0.085 0.120 m.sub.h.sup.hh(Γ.sub.9v.sup.+) Γ .fwdarw. M 0.055 0.074 Γ .fwdarw. A 0.463 0.526 m.sub.e(Γ.sub.8c.sup.−) Γ .fwdarw. M 0.076 0.089 Γ .fwdarw. A 0.997 1.088 m.sub.e(Γ.sub.7c.sup.−) Γ .fwdarw. M 0.038 0.052 Γ .fwdarw. A 0.033 0.042
[0204] The oscillator strengths of optical transitions between the three uppermost valence and two lowest conduction bands of 2H-Ge are given in Table S10. Transitions that are dipole forbidden due to group-theoretical arguments are indicated by horizontal lines. These symmetry considerations corroborate the identified band ordering at the ┌ point (Γ.sub.7c.sup.−>Γ.sub.8c.sup.−). In
TABLE-US-00010 TABLE S10 Optical transitions between valence and conduction bands of 2H—Ge characterized by transition energy, optical transition matrix element, Kane energy, and oscillator strength. Transitions that are dipole forbidden by symmetry are indicated by horizontal lines. The values have been calculated with the HSE06 and MBJLDA functionals for light polarized perpendicular and parallel to the c axis. Transition Optical transition energy matrix element Kane energy Oscillator strength Transition Method ε.sub.ck − ε.sub.vk (eV) p.sup.⊥ (ℏ/a.sub.B) p.sup.∥ (ℏ/a.sub.B) E.sub.p.sup.⊥ (eV) E.sub.p.sup.∥ (eV) f.sup.⊥ f.sup.∥ Γ.sub.9v.sup.+ .fwdarw. Γ.sub.8c.sup.− HSE06 0.286 6.48 .Math. 10.sup.−3 — 2.29 .Math. 10.sup.−3 — 8.00 .Math. 10.sup.−3 — MBJLDA 0.298 5.92 .Math. 10.sup.−3 — 1.91 .Math. 10.sup.−3 — 6.39 .Math. 10.sup.−3 — Γ.sub.7v+.sup.+ .fwdarw. Γ.sub.8c.sup.− HSE06 0.419 — — — — — — MBJLDA 0.418 — — — — — — Γ.sub.7v−.sup.+ .fwdarw. Γ.sub.8c.sup.− HSE06 0.770 — — — — — — MBJLDA 0.730 — — — — — — Γ.sub.9v.sup.+ .fwdarw. Γ.sub.7c.sup.− HSE06 0.614 0.447 — 10.9 — 17.7 — MBJLDA 0.632 0.394 — 8.44 — 13.4 — Γ.sub.7v+.sup.+ .fwdarw. Γ.sub.7c.sup.− HSE06 0.748 0.384 0.388 8.01 8.21 10.7 11.0 MBJLDA 0.752 0.343 0.330 6.42 5.92 8.54 7.87 Γ.sub.7v−.sup.+ .fwdarw. Γ.sub.7c.sup.− HSE06 1.098 0.214 0.661 2.49 23.8 2.27 21.7 MBJLDA 1.065 0.178 0.603 1.72 19.8 1.62 18.6
[0205] Lonsdaleite Ge, being a direct semiconductor with a very weak lowest optical transition, exhibits significant variations in luminescence and absorption in comparison to cubic Ge. The effects can be expected to be stronger than for Si and SiGe alloys. For a clear illustration of the global light-emission properties, the radiative lifetime t as a function of temperature is shown in
[0206] The large gap difference between 3C-Ge and 2H-Ge and its consequence for the temperature-dependent band populations explain the huge difference of τ by several orders of magnitude for low temperatures. (Note that for the thermalization of electrons and holes in Ge nanocrystals similar curves have been published in the prior art.) Manipulation of the atomic structure of 2H-Ge by straining or alloying, for instance, may lead to an inversion of the ┌.sub.8c.sup.− and ┌.sub.7c.sup.−conduction states which is likely to drastically improve the light-emission properties of lonsdaleite Ge, thus providing a vast playground for engineering its optoelectronic performance.
[0207] The lonsdaleite (2H) phase of Ge, which can be grown using hexagonal III-V nanowire templates, is considered a good candidate for Si on-chip optical interconnects and Sicompatible quantum light sources, thanks to its predicted direct band gap. Since experimental data and reliable calculations on 2H-Ge are scarce and often inconsistent, we first established the computational approach for efficient predictive ab initio calculations in this work. The performance of several XC functionals of DFT was benchmarked, including meta-GGA and hybrid functionals, to calculate the experimentally and theoretically well known structural and electronic properties of diamond-structure (3C) Ge. In a second step, these functionals were used to predict the structural, electronic, and optical properties of lonsdaleite Ge.
[0208] The atomic structure of 2H-Ge was computed with the PBEsol functional which is shown to yield excellent lattice parameters for the well studied cubic phase of Ge. The electronic structures of cubic and lonsdaleite Ge were calculated with the HSE06 hybrid functional and the MBJLDA meta-GGA, finding consistent results with both approaches, and an excellent agreement with the available experimental data. The ┌.sub.8c.sup.− conduction-band minimum of lonsdaleite Ge results from the backfolding of the L point of diamond-structure Ge onto the ┌ point of the hexagonal BZ, while the ┌.sub.7c.sup.− conduction-band state, that is derived from the lowest conduction band at ┌ of cubic Ge, is pushed towards higher energies. The energetic ordering of the three highest valence bands is ┌.sub.9c.sup.+>┌.sub.7v.sup.+>┌.sub.7v.sup.−. While the spin-orbit splittings of the hexagonal and cubic phase are similar, a huge crystal-field splitting is observed in 2H-Ge. The crystalfield splitting is responsible for the small ┌.sub.9c.sup.+.fwdarw.┌.sub.8c.sup.− band gap of only about 0.3 eV. The second conduction-band minimum ┌.sub.7c.sup.− is higher in energy by about 0.3 eV. The calculated electron and hole effective masses of cubic Ge are in good agreement with values in literature. Consequently, reliable effective masses for electrons and holes in 2H-Ge is predicted.
[0209] The dipole-allowed and dipole-forbidden optical transitions between the uppermost valence bands and lowest conduction bands near the G point and their polarization dependence is consistent with the symmetry identification of the bands. It appears that lonsdaleite Ge is a semiconductor with a direct fundamental gap in the infrared which exhibits a non-vanishing but small optical oscillator strength only for ordinary light polarization. The optical transitions to the second lowest conduction band instead are dipole allowed with large oscillator strengths. It is noticed that the distance between the first and second conduction band, as well as the size of the band gap, appear to be sensitive to the structural parameters. Consequently, a careful investigation of the luminescence properties, including their time dependence, and the absorption edge, also considering effects of strain, are suggested to further clarify the optical and optoelectronic properties of the promising new material lonsdaleite Ge.