Temperature-driven topological quantum phase transitions in a phase-change material Ge2Sb2Te5.

The Ge2Sb2Te5 is a phase-change material widely used in optical memory devices and is a leading candidate for next generation non-volatile random access memory devices which are key elements of various electronics and portable systems. Despite the compound is under intense investigation its electronic structure is currently not fully understood. The present work sheds new light on the electronic structure of the Ge2Sb2Te5 crystalline phases. We demonstrate by predicting from first-principles calculations that stable crystal structures of Ge2Sb2Te5 possess different topological quantum phases: a topological insulator phase is realized in low-temperature structure and Weyl semimetal phase is a characteristic of the high-temperature structure. Since the structural phase transitions are caused by the temperature the switching between different topologically non-trivial phases can be driven by variation of the temperature. The obtained results reveal the rich physics of the Ge2Sb2Te5 compound and open previously unexplored possibility for spintronics applications of this material, substantially expanding its application potential.

The Ge2Sb2Te5 compound is a phase-change material (PCM) which has been long used in optical memory devices such as re-writable DVD-RAM and is also a leading candidate for next generation non-volatile electronic memory known as phase-change random-access memory (PC-RAM)12. Recently, a new concept of nanostructured PCMs has been developed based on [GeTe][Sb2Te3] short-period superlattices, referred to as interfacial phase-change materials345 among which Ge2Sb2Te5 (n = 2, m = 1) is regarded as a prototype conventional PCM. Initially it was assumed that in this type of materials the switching between memory states is due to the amorphous-crystalline phase-transition of the separate relatively thick superlattice sublayers456. However, latter it was demonstrated that the superlattice kept functioning while the GeTe sublayer thickness was narrowed down to 2–3 GeTe bilayers and thus it was concluded that the temperature-induced phase-change occurred within the crystalline state, rather than between amorphous and crystalline phases, as was verified with transmission electron microscopy3.

At present, it is established that Ge2Sb2Te5, i.e. [GeTe][Sb2Te3] system can adopt four different hexagonal layered structures in which the primary bonds in different layers are aligned according to the ordering of Ge, Sb and Te layers7. They are Kooi structure, experimentally the most stable phase of Ge2Sb2Te58, which is formed by nonuple layer (NL) building blocks, where GeTe bilayers are incorporated into the Sb2Te3 quintuple layer (QL); Petrov structure9, in which the [GeTe]2 blocks are sandwiched between Sb2Te3 QLs with Ge-Te-Te-Ge sequence within the block; Inverted-Petrov structure with Te-Ge-Ge-Te sequence; and so-called Ferro structure, in which the atomic layer sequence in the [GeTe]2 block has Ge-Te-Ge-Te order, i.e. like in ferroelectric bulk GeTe. The relative stability of these four structures depends on temperature. Earlier it was shown by means of ab-initio calculations of enthalpy as a function of temperature10 the Kooi phase has the lowest enthalpy at 0 K, in agreement with earlier density functional theory (DFT) calculation7 and experiment8. However, upon raising the temperature the enthalpy of the Kooi structure increases and above ≈125 K the Ferro phase becomes the most stable one10. The metastable Petrov and Inverted Petrov phases are involved into atomic mechanism of the phase transition between stable Kooi and Ferro phases10. Despite the results of ref. 10 correctly reproduce the sequence and stability of the Ge2Sb2Te5 crystal structures they underestimate the transition temperature which according to experiments is higher than the room temperature.

Despite the advances in comprehension of the crystal structure and atomic mechanisms of crystal phase transformations the electronic structure of Ge2Sb2Te5 phases is currently not fully understood. For the stable low-temperature NL structured Kooi phase the narrow-gap band insulator (BI) phase has been obtained theoretically7111213. It is in contrast to related NL-structured Ge2Bi2Te5 and Sn2Sb(Bi)2Te5 compounds, which are shown to be topological insulators (TIs)1214. The [GeTe][Sb2Te3] compound with opposite n, m indices (n = 1, m = 2) has also been predicted to be TI15. The calculations for Ferro and Petrov structures also predicted the BI phase, while another metastable structure, the Inverted Petrov structure is shown to possess the Dirac semimetal quantum phase7.

In this work, we focus on electronic structure of stable low-temperature Kooi and high-temperature Ferro structures. We start with ordered Kooi structure which was also considered in the previous works712. The equilibrium crystal structure of the Kooi phase, obtained using VASP code (see details in the Method section) is shown in Fig. 1(a). The calculated bulk band structure demonstrates the insulating state with a gap of 25.2 meV in the middle of the Γ-A direction (Fig. 1(b), red lines) which is trivial band insulator in terms of the topological invariant that is in agreement with earlier results71112.

Figure 1

(a) Unit cell of the low-temperature Ge2Sb2Te5 Kooi phase. (b) Bulk band structure of ordered Kooi phase with equilibrium crystal structure parameters (red solid lines – VASP calculation, blue dashed lines – ABINIT calculation). Electronic structure of the ordered Kooi structure with λ/λ0 = 1.4: (c) bulk band structure near the A-gap with indication of the weights of Te and Sb orbitals, (d) surface electronic structure. (e) Bulk band structure of Kooi-Ge2Sb2Te5 with Ge/Sb mixing, and (f) its magnified view near the A-gap with indication of the weights of Te and Ge/Sb p orbitals (VCA-ABINIT calculation). (f) Surface electronic spectrum for Kooi-Ge2Sb2Te5 with Ge/Sb mixing. Red and blue circles demote the positive and negative sign of the in-plane spin, respectively.

The obtained small band-gap value may be an indication that the system is close to the topological quantum phase transition (TQPT) and can be converted into the topological phase by increasing spin-orbit interaction strength. We artificially increased the spin-orbit interaction strength λ in the ordered equilibrium structure and found that it leads to shift of the gap towards the A point along with its narrowing. Upon further increasing the spin-orbit interaction strength the system has gone through the critical point of the TQPT (at λ/λ0 ≈ 1.2), the gap becomes inverted achieving at λ/λ0 = 1.4 a width of 76 meV at the A point (Fig. 1(c)). This result is in line with the fact that the TI phase was predicted for Ge2Bi2Te512 in which the atomic spin-orbit coupling strength is larger owing to the larger atomic mass of Bi as compared to Sb. The bulk band inverted topology should manifest itself in formation of the spin-polarized Dirac state at the surface. As can be seen in Fig. 1(d) the gapless surface state with typical for related NL-structured Ge2Bi2Te5 and Sn2Sb(Bi)2Te5 TIs dog-leg dispersion1214 arises.

However, it is known from experiments that stable Ge2Sb2Te5 contains mixed Ge/Sb atomic layers16. Earlier, the influence of Ge/Sb mixing on electronic structure of the Ge2Sb2Te5 Kooi phase was considered within ordered 2 × 2 supercell approach11. It was demonstrated that invariant depends on relative concentration of Ge in the inner and outer Ge/Sb layers of NL (which was varied from 0 to 100% in steps of 25% in that model, containing 4 atoms in an atomic plane). It was shown that the material is trivial BI when the Ge atoms completely occupy inner layers (as in Fig. 1(a)) while it is TI at 50/50 mixing in the Ge/Sb layers.

We took into account the Ge/Sb mixing in the calculations within virtual crystal approximation (VCA) using ABINIT code. First we checked the band structure for the ordered phase with ABINIT code and found it in a good agreement with the VASP result (Fig. 1(b), dashed blue lines). The calculated spectrum has a gap of 29.5 meV in the middle of the Γ-A direction. Next, according to the experiment16, we constructed Ge0.56/Sb0.44 and Ge0.44/Sb0.56 virtual atoms for the outer and inner Ge/Sb layers of NL, respectively. We find that at normal, zero-pressure conditions the minimum gap of 61 meV appears at the A point (Fig. 1(e)) and, being inverted (Fig. 1(f)) results in non-trivial topological invariant . This signifies that the Kooi structure possesses the TI quantum phase. Thereby we can conclude that Ge/Sb mixing in the Ge2Sb2Te5 Kooi structure effectively increases the spin-orbit interaction. As can be seen in Fig. 1(g) presenting the surface band structure, calculated for 5 NL slab of the Kooi phase with Ge/Sb mixing, the topological surface state with the Dirac point at the Fermi level arises in the spectrum.

The distinctive feature of the Ferro structure as opposed to other Ge2Sb2Te5 phases is the lack of inversion symmetry owing to Ge-Te-Ge-Te layer sequence in the [GeTe]2 block (Fig. 2 (a)). It results in the spin-orbit splitting of the bulk energy spectrum (Fig. 2(b)) which resembles the Rashba-like band splitting in the bulk GeTe17. The presented spectrum has a tiny gap of 11 meV in the A–H direction at k ≈ 0.1 Å−1 away from the A point. Near the gap the spin texture differs from Rashba spin-helical picture and demonstrates almost collinear spin alignment perpendicular to the A–H (k) direction with sizable S component of the same sign both below and above the gap (Fig. 2(c)). Away from the H–A–L plane at k = π/c ± 0.015 Å−1 (where k = π/c is the H–A–L plane) the gap closes forming a pair of the Weyl nodes (Fig. 2(d)) that is distinctive feature of the topological Weyl semimetal (TWS) phase.

Figure 2

(a) Unit cell of the high-temperature Ferro-Ge2Sb2Te5 structure. (b) Bulk band spectrum; red/blue circles indicate positive/negative sign of in-plane, S spin components. (c) Magnified view of the spectrum in the vicinity of the A point in -H–A–H direction with indication of S and S spin components (middle panels) and spin-resolved constant energy contours taken below and above the gap (left and right panels for -H–A and A–H directions, respectively); here arrows show in-plane spin direction and red/blue color correspond to positive/negative sign of S. (d) The dependence of the gap width on k|| at k = π/c − 0.015 and k = π/c + 0.015 Å−1; the zero-gap points correspond to the Weyl nodes.

In principle, TWS can be realized by breaking either time-reversal or inversion symmetry of topological Dirac semimetals18. In this regard, the existence of the TWS in the Ferro phase can be understood from the comparison with the intermediate Inverted Petrov structure. The Inverted Petrov structure possessing the Dirac semimetal phase differs from the Ferro structure only in the atomic layer sequence within [GeTe]2 block so that the latter one is inversion asymmetric.

The Weyl fermions in the bulk are predicted to provide realization of the chiral anomaly, giving rise to a negative magnetoresistance under parallel electric and magnetic fields, the semi-quantized anomalous Hall effect, unusual optical conductivity, non-local transport and local non-conservation of the current19202122232425. At the surface the bulk band topology should manifest itself via formation of unusual surface states which form disjoint Fermi arcs which connect the projections of the pairs of Weyl nodes onto the surface Brillouin zone2627. The Fermi arc surface states are predicted to show unconventional quantum oscillations in magneto-transport, as well as unusual quantum interference effects in tunneling spectroscopy27282930.

The Ferro phase can adopt six different surface terminations depending on Sb2Te3 QL and GeTe bilayers sequence near the cleavage plane. These terminations should differ in bending of the surface potential owing to polarity in the GeTe bilayers (Ge+0.4, Te−0.4). The geometries of the terminations are shown schematically in insets in Fig. 3. The calculated surface electronic spectra (Fig. 3, odd columns) demonstrate that all surface terminations hold the trivial Rashba-split surface states resulting from splitting off from bulk bands due to the band-bending effect which is negative for QL-GeTe-GeTe-, GeTe-GeTe-QL-, and GeTe-QL-GeTe- terminations (first column in Fig. 3) and positive for terminations shown in the third column. Besides these states, as can be seen in the Fermi surface maps (Fig. 3, even columns), each surface termination holds the Fermi arcs connecting the Weyl’s pairs. In most cases they connect the Weyl nodes within each - pair while in case of TeGe-QL-TeGe- termination the arcs connect points of neighboring pairs via the hole-like Rashba surface state which crosses the Fermi level twice in the - direction, however does it once in - and touches the Weyl nodes on the conduction band side. Similar effect of the Weyl nodes reconnecting has been observed recently in the TWS phase arising in BiTeI under pressure31.

Figure 3

Surface electronic structure for different Ferro-Ge2Sb2Te5 surface terminations (odd columns) which structure is shown in insets and corresponding Fermi surfaces (even columns).

Small red points on Fermi surfaces mark positions of the Weyl nodes. k and k axis directed along - and -, respectively.

In summary, on the basis of ab-initio calculations we provide an important ingredient to the physics of prospective phase-change material Ge2Sb2Te5 demonstrating that temperature-induced structural phase transformation is accompanied by the quantum topological phase transition from TI phase in the low-temperature Kooi crystal structure to TWS phase in the high-temperature Ferro structure. We also demonstrate that the Ge/Sb mixing in the low-temperature structure is crucial for formation of the TI phase. Together with earlier predicted Dirac semimetal phase for intermediate metastable Inverted Petrov structure the TI and Weyl semimetal phases form a rich topological family realized in the same material and switching between the topological phases is ensured by the temperature. The Dirac surface states of TI, prospective for spintronic applications as well as exotic bulk and surface electronic states of TWS providing an ideal platform for many novel physical phenomena, such as negative magnetoresistance, anomalous quantum Hall effect and chiral magnetic effect can be realized in Ge2Sb2Te5, substantially expanding the application potential of this material. Since Ge2Sb2Te5 compound naturally has a polycrystal structure with randomly oriented crystallites the utilization of the predicted surface states of the topological quantum phases requires precise epitaxial growth of the Ge2Sb2Te5 films.

Methods

Electronic structure calculations were carried out within density functional theory (DFT). For bulk band structure calculations we used the Vienna Ab Initio Simulation Package (VASP)3233. The interaction between the ion cores and valence electrons was described by the projector augmented-wave (PAW) method3435. Relativistic effects, including spin-orbit interaction (SOI), were taken into account. For this calculations, the PBE exchange-correlation functional36 was used and DFT-D3 van der Walls correction37 was applied for accurate structure optimization. To treat the disordered Kooi phase we employed a virtual crystal approximation (VCA) as implemented in the ABINIT code38, where the averaged potential of a virtual atom occupying a site in the Ge/Sb sublattice is defined as a mixture VV = xVG + (1 − x)VS of Ge (VG) and Sb (VS) pseudopotentials. In ABINIT calculations we used GGA-PBE Hartwigsen-Goedecker-Hutter (HGH) relativistic norm-conserving pseudopotentials which include the SOI39. For surface electronic structure calculations for the Weyl phase, first the results of VASP calculations were used in the WANNIER90 code40 to construct tight-binding model. The chosen basis consists of six spinor p-type orbitals for each atom: , , , , , . The low-lying s orbitals are not taken into consideration. Surface tight-binding model is derived from the bulk one with inclusion of band-bending effects obtained from direct surface calculations within DFT. The surface spectrum has been calculated within surface Green function approach4142. The invariant is calculated from the parity of occupied electronic states at the time-reversal invariant points of the bulk Brillouin zone43.

Additional Information

How to cite this article: Eremeev, S. V. et al. Temperature-driven topological quantum phase transitions in a phase-change material Ge2Sb2Te5. Sci. Rep. 6, 38799; doi: 10.1038/srep38799 (2016).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.