1. Introduction
Brain-inspired computing [
1,
2,
3], light- or electrically-controlled reconfigurable optical metamaterials over a wide spectral range from the visible to THz [
4,
5,
6], holographic and nonvolatile displays of high resolution [
7,
8], integrated photonic circuits [
9,
10] and many other applications need phase-change materials (PCMs) of the next generation with better energy efficiency, wider temperature and spectral range for reliable operation compared to current flagship PCMs as Ge
2Sb
2Te
5, doped Sb
2Te, etc. Gallium tellurides seem to be promising candidates to achieve the necessary requirements because of higher melting and crystallization temperatures, combined with a low switching power and fast switching rate [
11,
12,
13]. Simultaneously, Ga
2Te
3 and non-stoichiometric compositions appear to be atypical PCMs, characterized by regular tetrahedral structure, the absence of Peierls distortion in the crystalline phase and metavalent bonding [
14,
15,
16]. The sp
3 gallium hybridization in cubic and amorphous Ga
2Te
3 [
16,
17] is also different from conventional p-bonding in flagship PCMs, raising a question of the phase-change mechanism. In addition, gallium tellurides exhibit a number of unexpected and highly unusual phenomena as nanotectonic compression [
15] or viscosity anomaly just above melting [
18,
19]. The nanotectonic compression involves simultaneous co-crystallization of the stable ambient and metastable high-pressure (HP) forms on usual heating of glassy g-GaTe
3 above the glass transition temperature
. The appearance of metallic HP-polymorphs seems to be beneficial for the PCM performance because of higher optical and electric contrast, accompanied by lower power consumption and a possibility of multilevel writing [
16]. The viscosity anomaly for Ga-Te melts in the vicinity of sesquitelluride composition [
18,
19] appears to be more significant than the observed “double kink” in liquid germanium tellurides and other PCM alloys [
20,
21,
22,
23], allowing to distinguish between two contrasting models: (1) fragile-to-strong transition [
20,
21,
22], and incipient liquid-liquid immiscibility [
23]. Recent high-energy X-ray and neutron diffraction measurements [
24] have shown an enhanced small-angle scattering below the scattering vector
≲ 0.4 Å
-1, related to dense metallic liquid droplets in a semiconducting Ga
2Te
3 melt and correlating with a non-monotonic viscosity behavior, which can be quantitatively described by the Taylor model [
25] for two-phase emulsions. The incipient transient immiscibility within a semiconductor – metal transition in liquid telluride PCMs is an interesting topic for further studies.
The main goal of the current report is to unravel the atomic structure of amorphous Ga
2Te
5 PLD film using high-energy X-ray diffraction, supported by first-principles simulations, and compare it with the crystal structure of tetragonal gallium pentatelluride [
26,
27]. In contrast to cubic gallium sesquitelluride Ga
2Te
3, which is a congruently melting compound in the Ga-Te binary system [
28,
29,
30], tetragonal Ga
2Te
5 is stable over limited temperature range and exhibits peritectic decomposition before melting. Consequently, the relationship between amorphous material, obtained by nearly instantaneous freezing of highly excited fragments, particles, liquid globules, etc., existing in the laser-induced plasma (plume) [
31], and metastable crystal is expected to be complex, leaving room for various intermediate configurations and states. A deep insight into the atomic structure and respective electronic, optical, thermal and other properties is a key for rational design of the next-gen PCMs and new functional materials for photovoltaics, thermoelectric, DNA sensing and energy storage applications [
32,
33,
34,
35,
36,
37].
3. Thermal Properties and Crystallization on Heating
The obtained Ga
2Te
5 PLD thin films appear to be amorphous and vitreous according to XRD and DSC measurements. Typical DSC traces of Ga
2Te
5 PLD and bulk glassy g-GaTe
3 are shown in
Figure 1(e). The endothermic glass transition temperature
increases with gallium content from 448 K (g-GaTe
3, 25 at.% Ga) to 491 K (g-Ga
2Te
5, 28.57 at.%) accompanied by exothermic crystallization. Bulk g-GaTe
3 shows two intense thermal features, peaked at 492 and 602 K, and a very weak intermediate effect at 547 K. The 492 K peak corresponds to primary Te crystallization (both usual trigonal
and high-pressure monoclinic
forms), while the 602 K feature is giving rise to cubic (
) and high-pressure rhombohedral (
) Ga
2Te
3, suggesting the nanotectonic contraction in a viscous supercooled liquid [
15].
In contrast, g-Ga
2Te
5 PLD film shows a single narrow and intense exothermic effect peaked at 545 K with a crystallization onset at
≈ 535 K. The
in situ XRD measurements with a typical DSC heating rate of 10 K min
-1 were used to reveal the nature of crystallizing phase(s). Surprisingly, the first weak Bragg peaks at the scattering angles
= 30.5° and 63.2° (the Co
incident radiation) have appeared just in the vicinity of
at ≈483 K,
Figure 1(a). They correspond to the (002) and (042) reflections of tetragonal Ga
2Te
5,
[
26], which was reported to be stable only between 673 and 768 K,
Figure 1(d). These unexpected results suggest the existence of certain Ga
2Te
5 motifs in the vitreous PLD film evolving into nano-crystallites on heating in a viscous supercooled liquid. Further crystallization advances in the vicinity of
when the remaining Bragg peaks of tetragonal Ga
2Te
5 become visible and grow over the 513 ≲
≲ 653 K temperature range,
Figure 1(b). In addition to the majority gallium pentatelluride
polymorph, trigonal tellurium,
[
57], was also detected as a minority phase,
Figure 1(c,f). Cooling the crystallized sample, the observed phases remain intact, specifically tetragonal Ga
2Te
5, metastable at room temperature.
Even more surprisingly, gallium pentatelluride appears to be perfectly stable after 15 months at room temperature,
Figure S1 (the Supplementary Information), in contrast to bulk Ga
2Te
5, transforming into cubic Ga
2Te
3 and trigonal tellurium within several weeks [
28]. In other words, a controlled crystallization of the amorphous Ga
2Te
5 PLD film yields a high quality stable tetragonal crystal promising for photovoltaics, thermoelectrics, energy storage and memory applications [
32,
33,
34,
35,
36,
37]. On the contrary, slow cooling or fast quenching of molten Ga
2Te
5 gives a polycrystalline mixture of cubic gallium sesquitelluride and trigonal Te,
Figure S1, fully consistent with the Ga-Te phase diagram,
Figure 1(d).
4. Electric and Optical Properties
The measured electrical conductivity of bulk crystalline Ga
2Te
5 is shown in
Figure 2(a). In contrast to previously reported results [
58], the electrical conductivity follows the Arrhenius relationship over the entire temperature range
where
is the pre-exponential factor,
the conductivity activation energy,
and
have the usual meaning. Nevertheless, the derived activation energy
= 0.227±0.003 eV is identical to that for intrinsic conductivity in tetragonal single-crystal Ga
2Te
5, measured in the direction perpendicular to the
-axis [
58]. The conductivity pre-factor,
= 17±3 S cm
-1, is located at the lower limit of electronic transport regime over the extended states [
59], that is,
= 0.45 eV is approximately the electrical bandgap.
The electrical conductivity for glassy g-Ga
2Te
5 was obtained by interpolation of the available data for amorphous and glassy
alloys, 0 ≤
≤ 0.4 [
15,
60,
61],
Figure 2(b). The room-temperature conductivity appears to be a decreasing exponential function of the gallium content
,
where
is the conductivity of amorphous Te, and
< 0 is a constant. In other words, the electronic conductivity of disordered Ga-Te materials essentially depends on tellurium concentration. The interpolated g-Ga
2Te
5 conductivity is lower by two orders of magnitude than
of bulk crystalline pentatelluride. The estimated conductivity activation energy,
= 0.41 eV, and the pre-exponential factor,
≈ 220 S cm
-1, suggest the electrical bandgap of glassy polymorph is about 0.8 eV.
Figure 3 shows the optical properties of Ga
2Te
5 PLD films. The absorption measurements reveal the fundamental optical absorption edge below the incident wavelength λ ≲ 1.2 μm accompanied by distinct interference fringes indicating a homogeneous nature and uniform thickness of the PLD film. The observed fringes allow both the refractive index
and the film thickness to be calculated using the Swanepoel method [
62]. In addition, the well-defined interference over the used spectral range and rather thick PLD film (2.7 μm) enable also the refractive index dispersion
to be estimated.
Two approaches are usually applied to represent the refractive index dispersion [
63]. The Cauchy approximation of the derived data
is given in
Figure 3(c)
where
,
, and
are constants which are characteristic of any given material. Since the Cauchy equation is inappropriate in a region of anomalous dispersion [
63], the Sellmeier approach is often used, considering the existence in an optical material of dipole oscillators with a resonance frequency
:
where
and
are two characteristic constants, and
is the speed of light. Usually, the Sellmeier equation is written with a series of terms to account for different resonance frequencies over an extended domain, that is,
,
, etc.
The Sellmeier coefficients
and
allow the normal dispersion of optical glasses to be calculated over a wide spectral range. In our case, we were limited to the original Sellmeier equation (5) with the following coefficients:
= 3.5017 and
= 0.3992,
Figure 3(d). The higher order terms were inaccessible caused by insufficient spectral range and experimental uncertainty. Nevertheless, the two approaches describe reasonably well the derived
values.
The optical absorption results were also used to calculate the optical bandgap
applying the Tauc relation [
64]:
where
is the absorption coefficient,
the photon energy, and
≈ 10
5 cm
-1 eV
-1 is a constant.
As expected, the derived bandgap
= 0.98±0.02 eV,
Figure 3(b), for glassy gallium pentatelluride was found to be smaller than that for g-Ga
2Te
3, 1.20 eV [
16], supporting a predominant role of the tellurium content on electronic and optical properties of Ga-Te alloys. Simultaneously, the optical
appears to be comparable with the electrical counterpart,
= 0.82 eV.
Thermally annealed and crystallized Ga
2Te
5 PLD film exhibits more complicated optical absorption,
Figure 4. The absorbance below
≲ 1 μm shows a distinct blue-shift, while the low-energy absorbance becomes more intense and mostly loses interference fringes indicating less homogeneous material in both the chemical composition and thickness. Taking into account the existence of (nano)crystallites in the annealed PLD film, additional scattering corrections were applied simultaneously with the usual reflection corrections. The refractive index
of g-Ga
2Te
5 was used for these calculations.
The Mie theory of light scattering for turbidity
measurements and the wavelength exponent
were employed for the scattering corrections [
65,
66]
where
is the scattering path length, and
and
are the intensities of the incident and transmitted beam, respectively. The turbidity depends on several parameters
where
is the scattering coefficient,
the particle number density,
the average volume of the particle, and
is the wavelength exponent. Combining Eqs. (8) and (9), one obtains
where
is the optical density, and
is a constant. The values of
and
were obtained plotting
vs.
(
Figure S2), which allows both the turbidity
and
to be determined and the scattering corrections to be calculated. The theoretical Heller wavelength exponent
[
65] yields the average particle size
, which appears to be
≈ 110 nm for c-Ga
2Te
5,
Figure S2. The derived
value is consistent with the size of crystallites, obtained from the XRD linewidth
, yielding
> 50 nm. The final absorbance corrected for reflection and scattering is shown in
Figure 4(a).
The derived optical absorption coefficient α, presented in
Figure 4(b), exhibits two optical processes above and below the incident photon energy
≈ 1.3 eV. Basically, the overall shape of α is reminiscent of that in silicon or carbon [
67,
68], which is related to direct optical transition at high
and indirect optical absorption at low photon energies [
68,
69,
70,
71,
72]. Assuming direct optical transition in crystallized Ga
2Te
5 above 1.3 eV
the direct optical bandgap was found to be
= 1.36±0.03 eV,
Figure 4(c). The constant
in Eq (11) is given by
, where
is the electron charge,
is a reduced electron and hole effective mass [
69].
The optical absorption plotted as
vs. photon energy
, Eq. (7), yields the indirect optical bandgap
= 0.40±0.03 eV,
Figure 4(d). The derived value is consistent with the electrical bandgap of c-Ga
2Te
5,
= 0.45 eV,
Figure 2(a).
The experimental data for tetragonal Ga
2Te
5, obtained using the conductivity and Hall effect measurements, are strongly anisotropic and changing over a wide range between 0.46 and 1.79 eV [
58]. The calculated
values are also variable, 0.86 ≤
≤ 1.7 [
72,
73], depending on the applied simulation method. Nevertheless, the results of electrical and optical measurements show a reasonable contrast between amorphous (SET) and crystalline (RESET) states for Ga
2Te
5.
5. High-Energy X-Ray Diffraction
The high-energy X-ray diffraction data in
-space are shown in
Figure 5. In contrast to fast cooled Ga
2Te
5 in a thin-walled silica tube, mostly consisting of cubic Ga
2Te
3 and trigonal tellurium with some vitreous fraction (
Figure S1), the splat quenching of tiny Ga
2Te
5 droplets yields essentially glassy material with non-negligible nanocrystals of cubic gallium sesquisulfide,
Figure 5(a). The spontaneous Ga
2Te
3 crystallization is consistent with the Ga-Te phase diagram,
Figure 1(d), related to peritectic decomposition of Ga
2Te
5 above 768 K. The obtained Ga
2Te
5 PLD films are fully vitreous with a distinct glass transition temperature at 491 K,
Figure 1(e). The X-ray structure factor
of g-Ga
2Te
5 PLD appears to be intermediate between bulk glassy GaTe
4 and Ga
2Te
3 PLD film,
Figure 5(b), suggesting structural similarities and revealing a systematic evolution of vitreous
materials with increasing gallium content
independently on preparation techniques.
In particular, we observe an emerging and growing first sharp diffraction peak (FSDP), also shifting to lower
with increasing
from
= 0.94±0.01 Å
-1 (GaTe
4,
= 0.2) to 0.86±0.01 Å
-1 (Ga
2Te
3,
= 0.4),
Figure 5c). The isolated FSDPs were obtained using the subtraction procedure [
74,
75]. The FSDP systematics (position
and amplitude
) reveals monotonic nearly linear trends as a function of
,
Figure S3, for both bulk glasses (0.17 ≤
≤ 0.25) and vitreous PLD films (0.2857 ≤
≤ 0.40), suggesting structural similarities on the short- and intermediate-range scale.
Distinct high-
oscillations, clearly visible over the extended
-range for the Ga
2Te
5 interference function
-
,
Figure 5(d), enable high real-space resolution for atomic pair-distribution
and total correlation
functions after usual Fourier transform
where
is the experimental number density,
the Lorch window function, and
= 30 Å
-1.
The derived
for g-Ga
2Te
5 PLD film is shown in
Figure 6. The asymmetric feature between 2.4 and 3.2 Å corresponds to Ga-Te and Te-Te nearest neighbors (NN). A Gaussian fitting (
Table 1) yields tetrahedral gallium coordination at 2.64 Å, consistent with the Ga-Te coordination numbers and NN distances in crystalline and glassy gallium tellurides [
14,
15,
16,
26,
27,
28,
77,
78]. On the contrary, the Te-Te atomic pairs in glassy Ga
2Te
5 are markedly shorter (2.80 Å) than those in tetragonal gallium pentatelluride (3.027 Å) [
26]. Nevertheless, this Te-Te NN distance is typical for amorphous and trigonal tellurium [
57,
79] and Te-rich binary and ternary glasses [
15,
77,
78,
80]. The partial
correlation functions for tetragonal gallium pentatelluride are compared in
Figure 6(b) with experimental
for g-Ga
2Te
5. We note both similarities and differences for the two materials.
The crystalline counterpart is stable over a narrow temperature range from 673 to 768 K [
28,
29,
30],
Figure 1(d). In contrast to layered Al
2Te
5 and In
2Te
5 [
27], gallium pentatelluride has a 3D structure consisting of infinite chains, parallel to the
axis, formed by edge-sharing ES-GaTe(II)
4 tetrahedra,
Figure 6(c). Every four ES-GaTe(II)
4 entities, belonging to the neighboring chains, are linked together by Te(I) species, which are located in the center of squares, formed by Te(II).
These square-planar Te
5 units (crosses) are presumably mostly missing in glassy Ga
2Te
5 since the Te-Te NN coordination number is significantly lower in the glass,
= 1.01±0.08 (
Table 1), compared to that in the crystal,
= 1.6 = ⅕×4 + ⅘×1; the average Te-Te coordination for the Te(I) and Te(II) species,
Figure 6(c). The ES-GaTe
4 tetrahedra seem existing in the glass network, since a weak shoulder at 3.39±0.02 Å is observed for the asymmetric second neighbor peak, centered at ≈4.3 Å. The short Ga-Ga second neighbor correlations, characteristic of ES-units in tetragonal Ga
2Te
5, are located at 3.424 Å,
Figure 6(b). However, a weak average Ga-Ga weighting factor,
= 0.02964 vs.
= 0.28405 or
= 0.68630, and, consequently, a small amplitude of this feature as well as the truncation ripples, related to a finite
-range of the Fourier transform, enable only a rough estimation of
= 1.2±0.4, compared to
= 2. A deep insight into the atomic structure of vitreous Ga
2Te
5 PLD films yields first-principles molecular dynamics.
6. First-Principles Molecular Dynamics
Simulated FPMD X-ray structure factor
and pair-distribution function
for glassy Ga
2Te
5 in comparison with experimental results are shown in
Figure 7(a,b). We note that the GGA approximation with hybrid PBE0 functional describes well the experimental data as it was reported earlier [
15,
16,
44,
45,
46]. The positions and amplitudes of the diffraction features in both
- and
-space are reproduced.
The calculated partial structure factors
are displayed in
Figure 7(c). As expected, the main contribution to the FSDP comes from the Ga-Ga partial
. The simulated
,
Figure 7(d), reveal complicated short- and intermediate-range order.
The asymmetric Ga-Те NN correlations are peaked at 2.62 Å and suggest at least two contributions with slightly different bond lengths. The Ga-Te coordination number is consistent with the experiment,
= 3.97,
Table 1, assuming a tetrahedral gallium local environment. In addition to homopolar Te-Te bonds at 2.80 Å, a weak Ga-Ga NN feature at 2.42 Å was also found. The amplitude of this peak is too small to be observed experimentally,
Figure 7(b). The Ga-Ga second neighbors between 3 and 4.5 Å have a bimodal distribution. The shoulder at ≈3.35 Å indicates the ES-units, while the main contribution at 3.92 Å is related to corner-sharing CS-entities. Consequently, the fraction of ES-GaTe
4,
= 0.45, is significantly lower than that in tetragonal polymorph,
= 1, that is, only ES-GaTe
4 are present. The experimental value,
= 0.6±0.2, is reasonably consistent with the FPMD result. Basically, the experimental and FPMD structural parameters were found to be similar or identical,
Table 1.
The Ga and Te local coordination distributions are presented in
Figure 8. The tetrahedral gallium coordination
contains negligible number of Ga-Ga homopolar pairs (
= 1). In contrast, tellurium has multiple coordination environments
, 1 ≤
≤ 4, but only two-fold Te
2F (50.5%) and three-fold Te
3F (47.6%) coordinated species appear to be abundant. The tellurium forms reveal a large variability in Te-Te bonds, 0 ≤
≤ 4, from pure heteropolar Te-Ga coordination (
= 0) to fully homopolar environment (
=
). We should, however, note a small fraction of Te
4F species (1.87%) and a negligible number of Te
5 square-planar entities (0.23%), the only form of tellurium subnetwork in tetragonal gallium pentatelluride,
Figure 6(c). This result is coherent with the reduced Te-Te coordination number
,
Table 1.
The geometry of GaTe
4 units yields either Te-Ga-Te bond angles or the orientational order parameter
[
81,
82].
Figure 9(a,c) shows the calculated
bond angle distribution for g-Ga
2Te
5 in comparison with tetragonal Ga
2Te
5 and cubic Ga
2Te
3, respectively. We note a broad and slightly asymmetric
function, centered at 103.3±0.3°, for the PLD film (
Figure S4). The Te-Ga-Te angles in the two crystalline references, characterizing both distorted ES-GaTe
4 tetrahedra in gallium pentatelluride and regular CS-units in cubic sesquitelluride, are located within the glassy envelope but not reproducing it by simple broadening. Nevertheless, the tetrahedral geometry in tetragonal Ga
2Te
5 seems to be closer to that in the glass.
The connectivity of GaTe
4 tetrahedra is given by the Ga-Te-Ga triplets or the respective
distributions,
Figure 9(b,d). A difference to the crystalline connectivity in this case is even more significant but a remote resemblance to connected ES-entities in the tetragonal crystal also exists.
The order parameter
[
81,
82] is often used to evaluate the polyhedral topology and distinguish between tetrahedral and non-tetrahedral local geometry of four-fold coordinated GaTe
4 entities
where
is the Te-Ga-Te angle of the central gallium atom with its nearest Te neighbors
and
. The average value of
changes between 0 for an ideal gas and
= 1 for a regular tetrahedral network. The
probability distribution function is shown in
Figure 10(a). Asymmetric
is peaked at
= 0.93 and decreases sharply both ways to high and low
. Usually, the tetrahedral limit is set at
≥ 0.8 [
44,
83]. The
integration within these limitations shows that 97% of GaTe
4 units belong to tetrahedral geometry. The remaining entities (3%) presumably are defect octahedral species GaTe
4 with two missing Te neighbors around the central Ga atom, characterized by 0.4 <
< 0.8. The regular defect octahedron (
and
) has
= ⅝.
Two-fold and three-fold coordinated tellurium can explain the asymmetric shape of the Ga-Te NN peak. The calculated Ga-Te
2F and Ga-Te
3F distances are presented in
Figure 10(b). The two distributions are broad and asymmetric but have slightly different maxima. The Ga-Te
3F bonds are longer (a broad maximum at 2.71±0.03 Å) compared to Ga-Te
2F, peaked at 2.63±0.02 Å. Similar difference,
= 0.12±0.01 Å, was reported in monoclinic Ga
2S
3 [
84] with ordered distribution of gallium vacancies. The ratio of Ga-Te
3F to Ga-Te
2F bond populations,
= 1.75, was also found to be similar to the expected stoichiometric ratio
= 2 for tetrahedral Ga species, having the formal oxidation state Ga
3+ and Ga
2Te
3 alloy composition.
The connectivity analysis shows that all Ga and Te species are connected. The analysis of Te-Te connectivity reveals a different size of
fragments,
Figure 11(a). Tellurium monomers (
= 1), that is, Te atoms with only heteropolar Te-Ga bonds, and dimers (
= 2) represent a relative majority, 55%, of all
fragments. The remaining fragments can be divided into two groups: (i) 3 ≤
≤ 6 (see the inset in
Figure 11(a)), and (ii) oligomeric chains,
= 15 for the used 210-atom simulation box. Group (i) represents remnants of square-planar Te
5 units in tetragonal Ga
2Te
5, the inset in
Figure 11(b), confirmed also by bond angle distribution
. Group (ii) is similar to chains in trigonal tellurium,
[
57], supported by a contribution at about
≈ 103°,
Figure 11(b). The two groups appear as a consequence of limited thermal stability of tetragonal gallium pentatelluride and peritectic reaction Ga
2Te
5 ⇌ Ga
2Te
3 +2Te above 768 K.
The intermediate-range order in glassy and amorphous materials is often described by ring statistics, that is, by population of
rings in case of gallium tellurides. The ring population
[
54] in glassy Ga
2Te
5 (this work) and Ga
2Te
3 [
16] PLD films in comparison with crystalline references: tetragonal Ga
2Te
5 [
26], cubic [
17] and rhombohedral [
85] Ga
2Te
3 is shown in
Figure 12. The
population was found to be different for the two Ga
2Te
5 forms. The dominant 8-membered rings in tetragonal polymorph are hardly populated in the PLD film. On the contrary, the most populated
= 5 entities in g-Ga
2Te
5 are absent in c-Ga
2Te
5. The peritectic nature of tetragonal crystal seems to be related to this difference.
Gallium sesquitelluride Ga
2Te
3 is a congruently melting compound,
Figure 1(d). As it was reported earlier [
16], the ring statistics in g-Ga
2Te
3 represents a disordered mixture of
for the ambient and high-pressure polymorphs,
Figure 12(b,d), related to nanotectonic contraction in a viscous supercooled melt.
Microscopic voids and cavities in amorphous Ga-Te alloys, obtained using the Dirichlet–Voronoi tessellation [
56], are displayed in
Figure 13(b). The fraction of voids
, normalized to the volume of the FPMD simulation box, was found to be nearly invariant, 27 ≤
≤ 29%, over the gallium content
between 0.20 (bulk g-GaTe
4) and 0.40 (PLD g-Ga
2Te
3). This is coherent with a small change in the number density, ≈2% over the same composition range. Typical cavity radius varies between 0.2 and 4 Å, slightly increasing with
,
Figure 13(a).
The total electronic density of states (eDOS) is shown in
Figure 14 and appears to be typical for glassy and crystalline chalcogenides [
15,
16,
59,
86,
87,
88]. The valence band (VB) consists of three sub-bands between the Fermi level
and −16 eV. The upper part, roughly centered at −3 eV, mostly consists of Te 5p and Ga 4p states, and also includes non-negligible d-electron contributions, as it revealed by the eDOS projections (pDOS) on Ga and Te atomic pseudo-wave functions. The middle-energy sub-band, centered at −8 eV, essentially contains Te 5p and Ga 4s electron states, while the lower part, peaked at −13 eV, has an s-character, populated by Te 5s electrons together with Ga s-, p- and d-electron contributions. The derived eDOS and pDOS are similar to those in g-Ga
2Te
3 PLD film [
16] and suggest sp
3 gallium hybridization having also d-electron contribution.
The inverse participation ratio IPR [
89,
90] allows localized (large IPR → 1) and extended electron states (small IPR ≈
, where
is the number of atoms in the simulation box) to be distinguished
where
is a single-particle Kohn-Sham eigenfunction. The calculated IPRs, derived using the projections of
onto an atomic basis set and the atomic orbital coefficients, are shown in
Figure 14(a), plotted together with the eDOS. As it was reported earlier [
16,
46], a higher electron localization appears at the band tails (the top of the valence and the bottom of the conduction bands), consistent with the theories of disordered semiconductors [
59]. The remaining electron states in the vicinity of the bandgap are delocalized (extended). Deeper states of the lower-energy sub-bands, participating in the covalent bonding, are localized even more strongly.
The derived GGA/PBE0 bandgap
= 0.80 eV appears to be smaller than the experimental optical bandgap
= 0.98 eV,
Figure 3(b), but nearly identical to the interpolated electrical counterpart
= 0.82 eV,
Figure 2(a). The main contribution of the Te 5p electron states to the upper part of the valence band and at the bottom of the conduction band is also consistent with a dominant role played by Te on the electronic conductivity of Ga-Te alloys,
Figure 2(b).