Sub-picosecond photo-induced displacive phase transition in two-dimensional MoTe$_2$

TL;DR

This study demonstrates a sub-picosecond, laser-induced phase transition in monolayer MoTe₂ from a semiconducting to a topological phase, driven by electronic excitations that soften lattice vibrations, enabling ultrafast control of topological properties.

Contribution

First-principles calculations reveal a purely electronic, ultrafast phase transition mechanism in 2D MoTe₂ triggered by specific photon energies, advancing control over topological phases.

Findings

Transition occurs within sub-picosecond timescale.

Triggered by photons with energy over 1.96 eV.

Carrier density needed is 3.4×10^{14} cm^{-2}.

Abstract

Photo-induced phase transitions (PIPTs) provide an ultrafast, energy-efficient way for precisely manipulating the topological properties of transition-metal ditellurides, and can be used to stabilize a topological phase in an otherwise semiconducting material. Using first-principles calculations, we demonstrate that the PIPT in monolayer MoTe$_2$ from the semiconducting 2H phase to the topological 1T$'$ phase can be triggered purely by electronic excitations that soften multiple lattice vibrational modes. These softenings, driven by a Peierls-like mechanism within the conduction bands, lead to structural symmetry breaking within sub-picosecond timescales, which is shorter than the timescale of a thermally driven phase transition. The transition is predicted to be triggered by photons with energies over $1.96$\,eV, with an associated excited carrier density of $3.4\times10^{14}$\,cm$^{-2}$, which enables a controllable phase transformation by varying the laser wavelength. Our results provide insight into the underlying physics of the phase transition in 2D transition-metal ditellurides, and show an ultrafast phase transition mechanism for manipulation of the topological properties of 2D systems.

Tables

Table 1. Table 1: Excited carrier concentration n 𝑛 n and the maximum amount of excitation N 𝑁 N by a laser pulse with a fluence of 100 mJ/cm 2 , as well as their corresponding excitation percentage in the total number of valence electrons N tot subscript 𝑁 tot N_{\textrm{tot}} , at different laser excitation energies.

Energy (eV)	1.58	1.96	2.34	2.42	2.63
$n$ ($\times {10}^{14}$ cm-2)	0.9	3.4	10.1	11.8	13.2
$n/N_{\text{tot}}$ (%)	0.8	2.9	8.5	9.9	11.1
$N$ ($\times {10}^{14}$ cm-2)	7.2	12.2	16.0	17.2	17.6
$N/N_{\text{tot}}$ (%)	6.0	10.2	13.5	14.5	14.8

Equations

$${\Gamma }_{\text{2H}}={\text{E}}^{\prime \prime }{\text{+A}}_1^{\prime }{\text{+E}}^{\prime }{\text{+A}}_2^{\prime \prime }.$$

$$V(u)=a+\frac{b}{2}{u}^2+\frac{c}{4}{u}^4,$$

Full text

A photo-induced phase transition (PIPT), resulting from cooperative electron-lattice interactions through transiently changing the electronic states of the solid by photoexcitations Nasu (2004), is completely different from thermally or pressure-induced phase transitions. A PIPT not only gives access to different phases in the solid to explore anomalous properties absent in the ground state phase Bisoyi and Li (2016), but also enables precisely controllable phase transitions towards target structures with desired physical properties at high speeds Wu and Liou (2011). Thus, the initial discovery of a PIPT in organic charge transfer crystals has triggered great interest in a variety of fields Koshihara et al. (1990). Compared with the vast number of materials exhibiting thermodynamic phase transitions, there are still very few solids undergoing a PIPT, and many investigations currently focus on bulk crystals possessing 1D correlated electron chains in which PIPTs may be easily triggered due to inherent instabilities caused by electron-electron and electron-lattice interactions in these quasi-1D systems Nasu (2004); Hanamura (1997).

As an emerging platform for PIPTs, 2D and layered transition-metal ditellurides exhibit polymorphisms with distinct physical properties Wang et al. (2017a); Fei et al. (2018); Zhang et al. (2019a). Recent experimental evidence for a PIPT in MoTe2 and WTe2 shows that structural phase transitions often involve abrupt changes in the electronic structure such as the emergence of novel topological states Cho et al. (2015); Sie et al. (2019); Yang et al. (2017). For instance, in few-layer MoTe2, an irreversible transition from a semiconducting hexagonal 2H phase to a topological distorted octahedral 1T*′* phase occurs under laser irradiation, which can be used to fabricate an ohmic heterophase homojunction with accurate control of micron-patterning in a desired area Cho et al. (2015). The phase transition is irreversible because the 1T*′* phase formed in the electronically excited state is metastable, and to return to the ground-state 2H phase, the system has to overcome a potential barrier of around 0.71 eV Kolobov et al. (2016); Krishnamoorthy et al. (2018). The laser-induced 1T*′* MoTe2 is stable for more than one week in ambient conditions without any protection Tan et al. (2018). Moreover, the structural distortion from the 2H to the 1T*′* phase results in an intrinsic band inversion between Te $p$ and Mo $d$ bands Qian et al. (2014); Choe et al. (2016), enabling ultrafast manipulation of the topological character of MoTe2. The topological phase also displays gate-tunable superconductivity, providing a new potential platform to realize Majorana bound modes Sajadi et al. (2018); Fatemi et al. (2018).

Despite all this progress, the microscopic nature of the photo-induced phase transition remains unclear. Triggered by external stimulation by light, the system undergoes changes in temperature, strain, electronic excitation, chemical state, and lattice vibrational modes. The transition may depend on one of these factors or on a combination of them Cho et al. (2015); Kretschmer et al. (2017). One proposal is that Te vacancies created by irradiation trigger the local phase transition Cho et al. (2015); Wang et al. (2017b); Yoshimura et al. (2018). Other theoretical and experimental studies have argued that accumulated heat is a main driving force for the phase transition Wang et al. (2017b); Tan et al. (2018). Moreover, a strain-induced phase transition has been observed in monolayer MoTe2 Duerloo et al. (2014); Song et al. (2016), thus laser-induced thermal strain may also contribute to the observed phase transition. Yet another explanation is that it is the electronic excitation that plays a critical role in the phase transition Kolobov et al. (2016); Krishnamoorthy et al. (2018). Overall, the transition mechanism is still debated.

In this work, we demonstrate that the phase transition of monolayer MoTe2 can be triggered by photoexcitation of carriers alone. Using first-principles calculations, we study the excited-state properties of 2H MoTe2 at photoexcitation energies $E$ of 1.58 eV (785 nm), 1.96 eV (633 nm), 2.34 eV (532 nm), 2.42 eV (514 nm), and 2.63 eV (473 nm) by exciting all available electrons consistent with these energies. We find that a soft mode displacive phase transition occurs for $E>1.96$ eV. Photoexcitation modifies the electronic potential, resulting in a complete softening of the lattice vibrations that leads to spontaneous structural distortions with new equilibrium positions for the atoms. This structural transition takes place within one picosecond, shorter than the timescale for the photoexcited electrons to transfer their energy to the lattice. We also rule out a thermally-driven or a strain-induced phase transition under laser irradiation by comparing the thermodynamic stability of 2H and 1T*′* MoTe2 monolayers. Our findings not only reveal the potential origin of the PIPT in 2D transition-metal ditellurides, but also provide new insight into the observed dynamic transitions in their 3D counterparts, another question that remains unclear Sie et al. (2019); Zhang et al. (2019b); Kim et al. (2017).

.1 Lattice mode softening

We start by examining atomic displacements in the excited state to see whether (and how) a phase transition may be triggered. With $D_{3h}$ point group, the irreducible representations of the vibrational modes in the 2H phase at the $\Gamma$ point are expressed as Yamamoto et al. (2014)

[TABLE]

The calculated phonon dispersion is shown in Fig. 1, and the calculated lattice mode frequencies are consistent with previous reports Kan et al. (2015); Guo et al. (2015).

The three modes E*′′, E′* and A${}_{2}^{\prime\prime}$ correspond to the three primary lattice distortions along the pathway from the 2H to the 1T*′* phase, so could play a role in triggering the phase transition (for the three primary lattice distortions, see “Crystal Structures” in Supplementary Information). We calculate the potential energy as a function of displacement amplitude of these vibrations under different photoexcitation energies and corresponding carrier densities to investigate potential structural instabilities. The results are summarized in Fig. 2(a)-(c), with the energy of the excited 1T*′* phase set to zero for comparing the relative stability between the 2H and 1T*′* phases. In the ground state, the energy of the 2H phase is $80$ meV lower than that of the 1T*′* phase. However, as the system is optically excited, the 2H phase becomes energetically less favorable than the 1T*′* phase due to population inversion Kolobov et al. (2016); Krishnamoorthy et al. (2018).

For the E*′′* mode (corresponding to the opposite in-plane motion of top and bottom Te layers), the Te atoms vibrate harmonically around their equilibrium positions before photoexcitation. With increasing excitation energy (and corresponding excited carrier density), the potential energy surface gradually flattens, and a double-well potential forms at $E>1.96$ eV. Thus, a distorted octahedral structure can form as atoms move to new equilibrium positions [Fig. 2(d)]. A similar situation occurs for the A${}_{2}^{\prime\prime}$ mode (an out-of-plane displacement of Mo atoms and an opposite-direction displacement of Te atoms). For $E>1.96$ eV, the Mo atoms can bounce up or down to the two minimum energy sites of the new double-well potential, while the Te atoms move in opposite directions, leading to an out-of-plane distorted variant of 2H MoTe2 [Fig. 2(f)]. The evolution of the E*′* mode under irradiation is slightly different. This mode is composed of in-plane displacements with the single Mo layer and the two Te layers moving opposite to each other. Increasing $E$ lowers the energy barrier for positive displacements along the eigenvector of the E*′* mode. The barrier becomes a deep well at 2.63 eV, leading to the atomic distortion shown in Fig. 2(e), with distorted in-plane atomic positions yielding a twisted octahedral coordination around Mo atoms. It should be noted that at 2.63 eV this mode dominates, but at lower energies the potential along the E*′* mode still exhibits a barrier whereas both E*′′* and A${}_{2}^{\prime\prime}$ modes show double-well potential energy surfaces, and therefore the latter two dominate in that regime. All three distortions of the corresponding mode eigenvectors are along the phase transition pathway from 2H to 1T*′* MoTe2. The full transition from 2H to 1T*′* MoTe2 probably occurs along a higher-dimensional path in configuration space in which photoexcitation lowers the energy barrier from 0.77 eV to 0.08 eV Krishnamoorthy et al. (2018), but our results suggest that the softening of these key modes can initiate the PIPT with a purely electronic mechanism.

.2 Critical excitation energy $E_{C}$

The transition from a single-well potential to a double-well potential shown in Fig. 2(a)-(c) indicates a displacive phase transition, which involves a phonon frequency that falls to zero at a critical excitation energy $E_{C}$ and corresponding excited carrier density. These zero-frequency soft modes drive the crystal instability along the displacements of the corresponding mode eigenvectors, leading to spontaneous symmetry breaking Dove (1993, 1997). We examine the soft-mode behaviour under optical excitation in the framework of Landau theory for the E*′′* and A${}_{2}^{\prime\prime}$ modes that soften at lower energies. For a soft mode displacive phase transition, the order parameter can be chosen as the amplitude of the distortion of the soft mode eigenvector $u$ Dove (1993). The lattice potential energy, equivalent to the free energy in the Landau theory Kittel (1976); Zeks and Blinc (1974), can be expressed by a fourth order even polynomial expansion of the order parameter $u$

[TABLE]

where the harmonic prefactor $b$ is equal to the square of the soft mode phonon frequency $\omega^{2}$ Dove (1993); Monserrat (2018). Before the phase transition, there is a single-well potential, corresponding to $b>0$ ($c$ is a small positive value with a much weaker $E$ dependence). When $b$ becomes negative, $V(u)$ becomes a double-well function with minima at $u=\pm\sqrt{-b/c}$ and a local maximum value at $u=0$. Therefore, the phonon softening can be described quantitatively by calculating the phonon frequencies of the two symmetric E*′′* and A${}_{2}^{\prime\prime}$ modes as a function of photoexcitation energy. As shown in Fig. 3(a), at 1.58 eV both modes soften slightly, and drop significantly at 1.96 eV. The phonon frequency becomes imaginary for $E>1.96$ eV, confirming again that a soft mode phase transition occurs at $E>1.96$ eV. It should be noted that at $E_{C}=1.96$ eV, the potential energy surface of the E*′* mode is still a single well. The energy surface in the subspace spanned by both soft modes is calculated to further evaluate the spontaneous lattice distortions. Fig. 3(b) shows the energy surface mapped on the plane spanned by the two soft modes E*′′* and A${}_{2}^{\prime\prime}$ at the photoexcitation energy of 2.34 eV. These results suggest that the A${}_{2}^{\prime\prime}$ mode dominates the distortion, thus the PIPT is probably triggered along that distortion.

In the Landau theory, the phase transition occurs at $b=0$. For temperature-driven ferroelectric/magnetic phase transitions, $b\propto(T-T_{0})$, giving the well-known Curie-Weiss law Kittel (1976). In the detailed analysis of the potential well of the A${}_{2}^{\prime\prime}$ mode, we find $b\propto(E_{C}-E)^{\gamma}$, where $E_{C}=1.96$ eV and $\gamma=0.24$. Since $b=\omega^{2}$, we fit the mode frequencies as a function of $E$ by $\omega\propto(E_{C}-E)^{\gamma/2}$, as shown in Fig. 3(a). The obtained $E_{C}=1.96$ eV is in good agreement with the change of electron distribution [Fig. 4(b)], from which one can see that the electron localization function also changes completely by $E\geq 1.96$ eV. A similar $E_{C}$ is also predicted for the E*′′* mode in Fig. 3(a). Thus our results suggest that the lattice distortions occur spontaneously after photoexcitation. Our prediction of photo-induced lattice mode softenings, marked by clear changes in both Raman and infrared spectra (E*′′* is Raman active and A${}_{2}^{\prime\prime}$ is infrared active), can be verified experimentally by ultrafast coherent phonon spectroscopy Zhang et al. (2019b) or second harmonic generation Song et al. (2018).

For a thermally induced phase transition, $b=\tilde{b}(T-T_{0})$, where $\tilde{b}$ is a positive constant Dove (1993), the potential well changes gradually with temperature towards the phase transition. In contrast, in a PIPT, the potential well change is determined by the number of excited carriers, which is related to the excitation energy $E$ and the density of states (DOS). For instance, when $E$ is only slightly larger than the band gap, few electrons are excited to the conduction bands, leading to a negligible change of the potential well. Thus in the photo-induced phase transition, $b$ has a much weaker $E$ dependence than the $T$ dependence in a thermally induced phase transition, leading to a much smaller $\gamma$.

.3 Photo-induced changes in the electronic state

We next investigate the electronic state to rationalize the phonon softening upon photoexcitation discussed above. When a solid undergoes a structural phase transition, the ionic positions change because the symmetry of the electronic potential changes Wall et al. (2012). Light-induced electronic excitations in the 2H lattice may be responsible for such a change in the electronic potential. Fig. 4(a) shows the 3D electron localization function of 2H MoTe2 at different photoexcitation energies. As shown in Fig. 4(b), the electron distribution around Mo atoms changes from isolated triangles (dashed red line) at 0 eV to connected hexagons (dashed red line) at 1.58 eV. As the excitation energy increases further to 1.96 eV, the electrons in the Mo layer gather between top and bottom Te atoms, and the in-plane triangles around Mo atoms appear once again, but with a reversed orientation compared to the ground state, indicating an electronic potential change. As the excitation energy increases to 2.34 eV, the inverted-triangle region of the electron density shrinks.

Figures 4(c) and (d) show a comparison of the DOS calculated for the static lattice and for the lattice distorted along the E*′′* and A${}_{2}^{\prime\prime}$ phonon modes, with the displacement amplitudes corresponding to the potential minimum at 2.34 eV. As expected from the discussion in Fig. 2 and 3, the A${}_{2}^{\prime\prime}$ mode distortion leads to a larger change in the DOS: there is a significant redistribution of density of states in the conduction band towards lower energies. A similar picture emerges for the E*′′* mode, but in that case the changes are small. In the ground state the conduction band is empty and therefore a redistribution of the density of states towards lower energies does not affect the total energy of the system. However, upon photoexcitation, the conduction bands become occupied, and it then becomes energetically favourable to undergo a lattice distortion to red shift the conduction bands and lower the overall energy of the system. This is analogous to the Peierls mechanism, and provides a microscopic picture for the driving force behind the mode softening induced by photoexcitation. We note that the DOS around the valence band maximum also changes, and when those states are occupied in the ground state, the distortion leads to an energy increase as expected. However, under photoexcitation those states are unoccupied (holes), and the DOS redistribution has a minor effect on the total energy.

.4 Excited carrier concentrations

As discussed above, for a given photoexcitation energy we consider that the number of incoming photons is as large as necessary to excite all available electrons that can be excited at the given laser energy. Therefore, the excited carrier concentration becomes a function of the excitation energy. Table 1 shows the excited carrier concentrations $n$ under different excitation energies. The excited carrier density at 1.96 eV is $3.4\times 10^{14}$ cm*-2*, corresponding to an excited carrier density of 2.9%. We also estimate the maximum amount of electron-hole pairs generated by a femtosecond laser with a fluence of 100 mJ/cm2 assuming that no saturation occurs [for the calculations of the maximum amount of electron-hole pairs, see “Within 20 fs of Laser Excitation” in Supplementary Information]. At 1.96 eV, excitation of up to 10.2% is possible for a laser pulse with a fluence of 100 mJ/cm2 (see Table 1), and thus the excited carrier concentrations in our calculations are experimentally accessible. As a PIPT can take place for excitations over 2.9%, a laser pulse with a minimum fluence of 28 mJ/cm2 is needed to induce a phase transition at 1.96 eV. This is consistent with the experimental fact that with insufficient laser fluence, the phase transformation cannot take place even for long irradiation times Cho et al. (2015); Tan et al. (2018).

The excited carrier population is related to the DOS Bernardi et al. (2014). Optical transitions are dominated by peaks in the DOS, and the dominant peaks in MoTe2 are marked by arrows in Fig. 5(c), where the optical transitions correspond to the excitation energies of 2.52 eV, 2.74 eV, and 2.96 eV. In 2D materials, reduced electronic screening typically leads to strong excitonic effect. The Coulomb interaction of excited electrons and holes redistributes oscillator strength, giving rise to tightly bound excitons with eigenvalues lower than the band gap. By including electron-hole interactions, the optical transition energies between the DOS peaks are reduced to 2.22 eV, 2.45 eV, and 2.68 eV [shown in parentheses in Fig. 5(c)]. With laser energies larger than 2.22 eV, the concentrations of excited carriers become higher, corresponding to strong optical excitation.

.5 Timescale of the phase transition

The timescale after laser excitation is divided into three regimes: (1) excited electrons thermalization (in 20 fs), (2) soft mode phase transition (117-292 fs), and (3) lattice heating (within hundreds of ps).

In the first stage, the temperature of the electronic system increases by photoexcitation and the excited electrons thermalize within 20 fs at a high electronic temperature, while the lattice remains cold. In photoexcited MoTe2, carrier-carrier scattering is responsible for redistributing carrier energy, resulting in a thermalized distribution function of carriers Shah (2013). The carrier relaxation time can be calculated from the quasiparticle self-energy Kotani and van Schilfgaarde (2010); Peng et al. (2018a). The excited electrons thermalize in 20 fs [Fig. 5(b)], whereas the electron-hole recombination occurs on the order of picoseconds [see “Within 20 fs of Laser Excitation” in Supplementary Information]]. Thus all the excited electrons are in thermal equilibrium before recombining or interacting with phonons.

In the second stage, the symmetry breaking displacements in 2H MoTe2 take place within 292 fs. As the laser pulse creates electronic excitations, the new potential energy surfaces in Fig. 2(a)-(c) are formed immediately. The photo-induced lattice distortions, as a result of phonon softening, occur within the timescale of lattice vibrations. This can be estimated by the lattice vibration periods from 117 fs (A${}_{2}^{\prime\prime}$ mode) to 292 fs (E*′′* mode). Such an estimate agrees well with previous molecular dynamics simulations Kolobov et al. (2016) and pump-probe spectroscopy Zhang et al. (2019b). The structural distortion can take place before the electrons transfer their energy to the lattice by carrier-phonon interaction (1-100 ps Shah (2013)).

In the third stage, the hot carriers are cooled by transferring energy to the lattice, converting electronic photoexcitation energy to heat in hundreds of picoseconds through interaction with various phonons.

.6 Competing mechanisms

One proposed mechanism for the PIPT in MoTe2 is that after the excited electrons transfer their energy to the lattice through electron-phonon coupling, atomic reorganization toward the 1T*′* phase may occur as a result of local heat Wang et al. (2017b) or thermal strain Duerloo et al. (2014); Song et al. (2016). We next demonstrate that these thermal effects alone cannot trigger the phase transition.

As the thermalized electronic system transfers its energy to the lattice, local heat is generated. The estimated local temperature during the laser-irradiation is 670 K Cho et al. (2015). As the temperature increases, the phonons in the 1T*′* MoTe2 phase are more highly excited than the phonons in the 2H MoTe2 phase, because the 1T*′* phase has more low-energy phonons than the 2H phase [Fig. 6(a)]. Thus the entropy of the 1T*′* phase becomes higher than that of the 2H phase as temperature increases (the entropy increases with the occupancy). It is thereby possible for the 2H structure to transform into 1T*′* MoTe2.

Thermodynamic stability at finite temperature can be described by the difference in Helmholtz free energy van Setten et al. (2007); Peng et al. (2017, 2018b) [for the calculations of thermodynamic stability, see “Thermal Effects (about hundreds of ps)” in Supplementary Information]. At 0 K, the free energy difference between the 2H phase and the 1T*′* phase is just 72 meV per rectangular unit cell, indicating that the stability of the 1T*′* phase is quite close to that of the 2H phase. As shown by the blue curve in Fig. 6(b), at temperatures higher than 930 K, 1T*′* MoTe2 becomes thermodynamically more stable than 2H MoTe2. This is consistent with the fact that the pure 1T*′* phase appears in a temperature range between 930 K (for Te deficiency) and 1170 K (for Te excess) Keum et al. (2015). However, the temperature to facilitate the phase transition, 930 K, is significantly higher than the laser-induced temperature of 670 K Cho et al. (2015). Even after taking thermal strain into consideration [red curve in Fig. 6(b)], the phase transition temperature of 810 K is still higher than the laser-induced temperature [for the role of thermal strain, see “Thermal Effects (about hundreds of ps)” in Supplementary Information]. Therefore the laser-induced thermal effect on the phase transition has probably been overestimated in earlier studies Cho et al. (2015); Wang et al. (2017b).

In fact, previous experiments have suggested that the heating effect alone cannot induce the structural transformation Tan et al. (2018). For pure 2H MoTe2, no significant change occurs after annealing; whereas laser irradiation leads to a mixed 2H + 1T*′* phase. In the laser-irradiated mixed phase, the 2H phase disappears completely after annealing, leaving a pure 1T*′* crystal Tan et al. (2018). Therefore, electronic excitation may drive spontaneous lattice distortions, while the heating effect accelerates the phase transition via thermal displacements, in which the intermediate states decay to the 1T*′* phase [for the role of thermal displacement, see “Thermal Effects (about hundreds of ps)” in Supplementary Information].

Another proposed mechanism for the PIPT in MoTe2 is that laser irradiation can create Te vacancies, which can trigger a local phase transition Cho et al. (2015); Wang et al. (2017b); Yoshimura et al. (2018). Consistent with this picture, a recent study has demonstrated that the phase transition involves a set of successive dynamics such as Te vacancy diffusion and ordering to generate the 1T*′* phase Si et al. (2019), which is similar to the Ge-Sb-Te phase-change memory materials Zhang et al. (2012). A vacancy-induced phase transition is related to a change in chemical composition, which is analogous to the composition-dependent phase transition in Mo1-xWxTe2 Lv et al. (2017). We emphasize that our work demonstrates that a phase transition can occur even in the absence of Te vacancies, and instead be driven by photoexcited carriers alone. Indeed, a recent experimental study has reported an electric-field induced structural transition in 2H MoTe2 Zhang et al. (2019a), supporting our conclusion that changes in electron density alone can induce a phase transition without the need for Te vacancies. Similar photo-induced transitions have recently been predicted in the perovskites BaTiO3 and PbTiO3 Paillard et al. (2019).

.7 Summary

In summary, we demonstrate that a sub-picosecond phase transition in monolayer MoTe2 from the 2H phase to a distorted 1T*′* phase can be triggered by purely electronic excitation, even when the lattice is still cold. Such electronic excitation changes the electron density, resulting in a soft mode displacive phase transition. The phase transition from the 2H phase to the 1T*′* phase can be induced by a critical excitation energy of 1.96 eV, at which 2.9% of the carriers are excited. We also show that this phase transition mechanism can be understood in the frame of Landau theory, and is quite different from a temperature-driven phase transition. The microscopic picture of the phase transition is analogous to a Peierls distortion, with a red shift of the conduction bands that are occupied upon photoexcitation that lowers the overall energy of the system. Our findings shed new light on the microscopic origin of PIPT in 2D transition-metal ditellurides, and are expected to motivate both fundamental and applied studies of ultrafast phase transitions in these new class of materials for topological switching and neuromorphic computing.

.8 Methods

All calculations are performed based on density functional theory (DFT) using the Vienna ab-initio simulation package (VASP) Kresse and Furthmüller (1996). The exchange-correlation energy is calculated within the generalized gradient approximation (GGA) in the Perdew-Burke-Ernzerhof (PBE) parameterization Perdew et al. (1996). In the excited state simulations, we mimick optical excitation by promoting valence electrons from high-lying valence band states to low-lying conduction band states. This method, typically called $\Delta$ self-consistent field ($\Delta$SCF) Jones and Gunnarsson (1989); Görling (1996); Hellman et al. (2004), introduces noninteracting electron-hole pairs by changing the occupation numbers of the Kohn-Sham orbitals and is less computationally demanding than more advanced approaches such as constrained DFT Mauri and Car (1995) and excited-state force calculations within $GW$+BSE Ismail-Beigi and Louie (2003) but has been successfully used in studying PIPTs Li et al. (2011); Kolobov et al. (2016); Chen et al. (2017). The occupied conduction band is determined by energy conservation at photoexcitation energies $E$ of 1.58 eV (785 nm), 1.96 eV (633 nm), 2.34 eV (532 nm), 2.42 eV (514 nm), and 2.63 eV (473 nm) (see “Computational Details” in Supplementary Information). We choose these laser excitation energies because the corresponding wavelengths are experimentally feasible. We keep these occupancies fixed throughout the calculations. This scenario corresponds to the time domain of interest, namely, after the laser-excited carriers populate the states near the band edges but before these carriers relax to the band edges via electron-phonon interactions. Electron-hole recombination occurs on the picoseconds timescale, which is much larger than a single phonon period and therefore the number of excited electrons can be reasonably assumed to remain constant as in the calculation. It should be noted that, for simplicity, we use the PBE band structures to determine the number of the excited electrons under optical excitation. The $GW_{0}$ band structure correction is an almost rigid upshift of the PBE conduction bands of 0.69 eV, and after a simple scissor shift, the required laser excitation energies are still experimentally feasible. In subsequent work, it would be interesting to see how bound excitations (excitons) modify the picture.

Data Availability

The raw data used in this study are available upon reasonable request from the corresponding authors.