Photoinduced multistage phase transitions in Ta2NiSe5

‡ 1 ‡ 1,2 ‡ 3 1 1 1 Q. M. Liu , D.Wu ∗ , Z. A. Li , L. Y. Shi, Z. X. Wang, S. J. Zhang, 1 1 3 3 1 1, 4 T. Lin, T. C. Hu, H. F. Tian, J. Q. Li, T. Dong, and N. L. Wang∗ 1 International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China 2 Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, China 3 Beijing National Laboratory for , Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 4 Collaborative Innovation Center of Quantum Matter, Beijing, China

Abstract: Utrafast control of material physical properties represents a rapid developing field in con- densed matter physics. Yet, accessing to the long-lived photoinduced electronic states is still in its early stage, especially with respect to an insulator to metal phase transition. Here, by combing transport mea- surement with ultrashort photoexcitation and coherent phonon , we report on photoinduced multistage phase transitions in Ta2NiSe5. Upon excitation by weak pulse intensity, the system is triggered to a short-lived state accompanied by a structural change. Further increasing the excitation intensity beyond a threshold, a photoinduced steady new state is achieved where the resistivity drops by more than four orders at temperature 50 K. This new state is thermally stable up to at least 350 K and exhibits the lattice structure different from any of the thermally accessible equilibrium states. Transmission electron microscopy reveals an in-chain Ta atom displacement in the photoinduced new structure phase. We also found that nano-sheet samples with the thickness less than the optical penetration depth are required for attaining a complete transition.

∗ Corresponding author: [email protected]; [email protected]

Ultrashort laser pulse not are only a powerful tool to excite and probe non-equilibrium electronic processes in transient states, but also emerges as an useful method to induce phase transitions that may or may not be thermally accessible. The latter has drawn increasing attention because it enables ultrafast optical manipulation and control over material properties [1–13]. Among various photoinduced (PI) phase transitions and related phenomena, the ultrafast switching from an insulating to a metallic state is particularly attractive, for it has high potential for device applications. While some phase transitions are transient, in the sense that they recover on very short time scales, e.g. picoseconds (ps), others are permanent or meta-stable and require something other than time to reset them. Up to now, the PI insulator to steady metallic state is achieved only in few systems exhibiting competing orders or broken symmetry phases, e.g. manganites perovskites [3, 4, 9, 13] and a charge density wave (CDW) crystal of 1T-TaS2[8], with their PI states only surviving at temperatures far below the room temperature. PI nonthermal insulator-metal phase transition was also reported in a strongly correlated material VO2 [14, 15], however, a recent study questioned the claim and found no evidence for a hidden transient Mott-Hubbard non-thermal phase in VO2 [16]. Here we report on distinct PI multistage phase transitions with insulator-to-metal transition (IMT) characteristics in Ta2NiSe5 -the compound is currently known as an important candidate of excitonic insulator (EI) [17–19]. We show that in Ta2NiSe5, the new steady phase induced by the strong laser pulses has a lattice structure different from any of the initial equilibrium states and is thermally stable up to at least 350 K. Ta2NiSe5 is a two-dimensional (2D) van der Waals compound with typical in-plane quasi-one-dimensional (1D) substructure composed of parallel Ta and Ni atomic chains along a-axis of the lattice (Fig. 1a). Upon cooling, the system shows a second-order phase transition at Tc = 326 K accompanied by a slight shear distortion of those atomic 1D chains, that reduces the lattice symmetry from orthorhombic to monoclinic [19]. Our start point is the observation of the sharp resistance drop of ultrathin Ta2NiSe5 crystals induced by ultrafast laser pulses. The samples, with thicknesses of ∼ 30 nm (less than the penetration depth ∼ 55 nm at λ = 800 nm), are exfoliated from the bulk crystals (See Methods. Images of typical devices are shown in Supplementary Figure 1). To induce the resistance switching, laser pulses from a Ti:sapphire amplifier system with 800 nm wavelength, 35 fs duration and 2 1 kHz repetition were firstly used as the “writing” pulses with a fluence of ∼ 3.5 mJ/cm . As shown in Fig. 1b, after writing pulses, the resistance drops approximately four orders of magnitude at 50 K thereafter retains in a low resistance (LR) state. This PI-LR state is thermally stable up to at least 350 K by the measurement. Unlike the pristine resistivity showing a transition at Tc of 326 K, no anomaly is indicated in the temperature dependent resistivity down to 1.8 K for PI-LR state (Supplementary Figure 2). The activation energy of the PI-LR state is tiny and estimated as 4.3 meV, twenty times smaller than 2 the one of pristine state. To induce the PI-LR state, a threshold fluence ∼ 2.5 mJ/cm is required. We remark that, due to the very high resistance at low temperature, we could not measure the value of the resistance for the pristine nano-thickness sample below 50 K. It is estimated that the resistance can drop by 8 orders at 4 K by such ultrafast laser pulse excitations. The PI-LR state is repeated from device to device, no matter the sample growth batches. To gain insight into the microscopic of the PI-LR state, we investigated the time-resolved photo-excitations by pump-probe spectroscopy. In the experiments, the energies of pump and probe 2 2 pulse are kept low at ∼ 50 µJ/cm and ∼ 3 µJ/ cm respectively, to ensure minimal disturbance of those states (see Methods). Fig. 2a depicts the photoinduced relative change of transient reflectivity ΔR/R(t). The rising time is about 100 fs. It can be clearly recognized that the decay process of photoinduced reflectivity change ΔR/R(t) is superimposed by several coherent oscillations [20]. The reflectivity change could be well reproduced by two-exponential decay processes[21, 22], a fast one in the time scale of 0.5 -0.6 ps and a slow one with ∼ 10 ps. Presence of rapid and slow decay dynamics after

FIG. 1. Resistivity switching of Ta2NiSe5 by the 35-fs laser pulse at 800 nm. a, The layered crystal structure of Ta2NiSe5. A quasi one-dimensional structure is formed by a single Ni chain and two Ta chains along a-axis. b, The four-probe resistance of the pristine (green) and PI-LR state (pink) after a 2 writing excitation ∼3.5 mJ/cm . Inset is the schematic of the sample and experimental configuration. The length of dotted red line indicates the time delay between the pump and probe pulses

excitation has been observed in many systems. In general, the number of photoexcited hot electrons at zero time delay is related to the amplitude of ΔR/R. Those excited high-energy hot electrons release and transfer their energy to lattice through the emission of optical phonons. The sub-picosecond decay would be mainly attributed to hot electron relaxation via the energy exchange with strongly coupled phonons, the ∼ 10 ps decay process would be related to the energy exchange with other phonons or dephasing process. The fast decay process is in accordance well with the electronic pre-thermalization process revealed by tr-ARPES experiments [23, 24]. The coherent oscillations arise from the coherent phonon excitations. We extract the coherent phonon spectra by fast Fourier transformation from transient reflectivity spectra after subtracting the electronic decay background, as shown in Fig. 2b. For pristine state at 50K, five phonon modes of 1.0-, 2.0-, 3.0-, 3.7-, 4.0-THz are detected. These modes are identified as Ag symmetries in earlier studies [25], arising from displacive excitation of coherent phonons [26]. Remarkable differences can be found between the spectra of pristine state and PI-LR state at low temperatures. As shown in Fig.2b, the 2.0-and 3.7-THz modes are absent for PI-LR state at 50K. Besides, the peak frequencies of the remaining modes of PI-LR state are a little higher comparing with the ones of pristine state respectively. These changes in phonon

FIG. 2. Pump-probe response spectra. a, Transient photoinduced reflectivity at various states. The signal is made up of the electronic response and the coherent oscillations. The gray lines represent fits to the measured data. b, the corresponding FFTs from transient reflectivity spectra after substracting the incoherent components. c, the temperature evolution of the coherent phonon spectra of pristine state (upper panel) and the PI-LR state (lower panel). In the experiment, the light polarizations of pump and probe are set parallel and perpendicular to a-axis respectively. number and frequency imply a PI lattice structural phase transition. More information of coherent modes evolution are recognized from the temperature dependent spectra (Fig. 2c). For the spectra of pristine sample, the 2.0 THz mode is present solely in the low temperature monoclinic phase, and the modes 3.7-and 4.0 THz gradually merge into one 3.7-THz across Tc accompanying the structural transition from monoclinic to orthorhombic phase. While for the PI-LR state, the temperature-dependent spectra show only three recognizable modes throughout the measured temperature range, irrespective of thermal cycles. Although the PI-LR phase does not show a phase transition in the measured temperature range, the coherent phonon modes in the PI-LR phase are very similar to those in the high temperature pristine phase in Ta2NiSe5, therefore, it is expected that the PI-LR phase is in the similar lattice symmetry but with different atomic construction comparing to the pristine high temperature orthorhombic phase. The structural transformation was directly clarified by transmission electron microscopy (TEM). Fig. 3a depicts the perspective view of Ta2NiSe5 atomic model. An appropriate [110]-type zone-axis (Fig. 3b) of Ta2NiSe5 was selected for TEM structural characterization, by which atomically resolved images can be readily obtained (see Methods). Fig. 3c shows the electron diffraction pattern and the corresponding morphologies of a [110]-oriented pristine nano-flakes at room temperature. The stripe-like domain in TEM morphology and the splitting or elongation of (3-30)-type reflections in diffraction pattern are arising from the twinning formation of monoclinic phase [27]. Upon heating across Tc, the lattice structure transforms into orthorhombic phase along with the disappearance of both the stripe-like domain and the splitting/elongation of (3-30) reflections (see Supplementary Figure 3). The high-angle annual dark-field (HAADF) atomic imaging experiments were performed. Fig. 3d shows an atomically resolved HAADF image of the [110]-oriented pristine state at room temperature, in which the HAADF intensities scale with the atomic number of constituent elements: the brightest dots correspond to Ta (Z = 73) columns and the least bright dots correspond to Ni (Z = 28) columns. A periodic structural unit is framed and false-colored in Fig. 3d, and its enlargement is shown in Fig. 3e together with the corresponding positions of Ta columns in the [110] projection. Note that it is hard to differentiate the HAADF atomic images between the low temperature monoclinic and high temperature orthorhombic phases of pristine Ta2NiSe5 in the limitation of TEM experiments, consistent with the fact that the two pristine phases have quite tiny differences in crystallographic parameters as revealed by X-ray experiments [28]. Similar TEM structural characterizations were carried out for writing-pulse excited Ta2NiSe5 nano-flakes. Fig. 3f shows the TEM morphology and the corresponding electron diffraction of 2 [110]-oriented Ta2NiSe5 after 800 nm pulses (at ∼ 3.5 mJ/cm ). Clearly, the twinning related stripe contrast in TEM morphology and the splitting of reflections in diffraction pattern are both absent, in stark contrast to the results (Fig. 3c) for the pristine state (See Supplementary Figure 4 for more [010]- oriented). Fig. 3g shows a representative [110]-oriented atomically resolved HAADF image and the Fig. 3h is the enlargement of the periodic structural unit frame correspondingly together with the Ta-columns extracted in [110]-projection. The clean and uniform atomic periods in Fig. 3g suggest that there no observable PI vacancies or defects in the induced hidden state. Despite the greater similarity in the chain-like structure along a-axis between Fig. 3d and 3g, careful inspection of their respective atomic arrangements unveils notable Ta-lattice shear distortion comparing Fig. 3e with Fig. 3h. The colored lines and arrows in both Fig. 3e and Fig. 3h mark the Ta lattice displacement, we can recognize clearly a small and a large shear motion respectively between different Ta-chains. The small one can be considered as a shear motion of the Ta against the Ni chains, with a counter Ta-Ta displacement measured to be about 0.5 ±(0.05) Å along < 110 >-type direction for the PI-LR phase. While the large shear motion recorded a relative displacement of about 1.6 ±(0.05) Å along < 110 >-type direction, which occurs between the two nearest-neighbouring Ta-chains, where the weakest Ta-Se bonds exist. Considering the large shear motion of Ta lattice, geometrically, it corresponds closely to a a/2 a-axis sliding between the nearest-neighbouring Ta-chains (see Supplementary Figure 5 for details). It would be interesting to compare the above coherent phonon excitations and structural characterization with the equilibrium structural transition revealed by Raman scattering studies [29–32]. Raman studies have revealed significant role played by the 2-THz mode in the structural phase transition. The mode has a B2g symmetry in the high temperature orthorhombic phase but evolves into an Ag symmetry in the low-temperature monoclinic phase. It corresponds to the in-chain Ta lattice oscillations or shear motion of TaSe6 octahedra along the a-axis[29]. An instability of this mode as the origin of phase transition was suggested by density function theory based calculations[29] and indeed observed in a recent Raman measurement[30]. In the present study, the optical pump-probe experiment detects only the Ag coherent phonon modes. In the pristine state below Tc, five Ag modes were detected at 50 K and two of them (2.0-and 4.0-THz) disappeared due to evolving into B2g modes above Tc, the result is in accordance with the Raman measurement. On the other hand, the 2-THz mode is also absent in the PI-LR phase. In fact, the pump fluence dependent pump-probe measurement in the monoclinic phase presented below indicates a clear weakening and disappearance of the 2-THz mode upon increasing fluence. It is expected that an offset displacement gradually develops in the direction of the in-chain phonon mode, and the oscillatory motion would be on top of that. Increasing the pump field would lead to a much larger offset displacement, similar to what has been observed in WTe2 and MoTe2 [33, 34]. Above a threshold pump field, the lattice would not recover to its original state and a new equilibrium

FIG. 3. TEM structural characterization of both pristine and laser-treated Ta2NiSe5. a, Perspective view of Ta2NiSe5 atomic model. b, Projection view along [110] zone-axis. c, For pristine Ta2NiSe5, [110]-oriented electron diffraction pattern and morphology, inset shows the zoom-in part of the splitting of (3-30) reflections due to the twinning formation. d, atomically-resolved high angle annual dark-field image taken along [110] zone-axis corresponding to (c), e, Enlarged part outline in (d), with a schematic of Ta atom columns corresponding to the brightest contrast. Descriptions for f, g, h are the same as (c,d,e) but for hidden state Ta2NiSe5. Red line and arrows in (h) mark the small relative displacement of Ta atoms with respect to the one in (e). The blue dashed arrows denote the locate where large Ta-atoms relative displacement occurred. Scale bars in (c) and (f) correspond to 500 nm, in (d) and (g) correspond to 0.5 nm. False colors scales with the image intensities for the purpose of visual clarity. state is reached eventually. The in-chain Ta atoms displacement found by TEM agrees well with the lattice instability of B2g mode proposed by the Raman work. Nevertheless, the large displacement of Ta atoms revealed by TEM measurement yields evidence that the PI new state is different from any of already known phases in pristine Ta2NiSe5. As the Ta lattice sliding has an intimate correlation with the excitonic insulator state in Ta2NiSe5 [30, 31], it can be expected that the pristine exciton condensation will be destroyed as soon as the Ta lattice transformation occurred. Note that previous report has revealed a high-pressure induced a/2 inter-layer-sliding transition in Ta2NiSe5 with the lattice transformed from monoclinic to orthorhombic symmetry along with a IMTs electronic transition [35]. As well, in pristine Ta2NiSe5, the electronic transition at 326 K is accompanied concomitantly by a lattice transition with the Ta sublattice slightly sheared along a-axis. It is therefore not surprising that the Ta lattice transformation will come along with an electronic phase transition in our observations in Ta2NiSe5 system. Recent structure investigation combined with laser-pumping and first-principle calculation study on layered van der Waals CDW material 1T-TaS2 suggested that the interlayer electronic ordering can play an important role in electronic phase transitions[36, 37]. The metal-insulator transition is driven not by the 2D order itself, but by the vertical ordering or stacking of the 2D CDWs. Among 13 different stackings, there exist two exceptionally stable stacking configurations, one being an insulator and the other being a metal. They have a small energy difference, and their competition is responsible for the metal-insulator transition in 1T-TaS2 [37]. The stacking scenario was used also to explain photoinduced insulator-metal phase transition in 1T-TaS2 compound. However, the observed in-chain displacement of Ta-Ta atoms in Ta2NiSe5 in the photoinduced LR state indicated that the distortion appears more complicated than a simple stacking fault proposed in Ref.[37]. To generate the PI-LR state in Ta2NiSe5, a pulse fluence beyond a threshold is required. Fig. 4a shows a shot-by-shot dependence of resistance measured at 50 K. It records a steep drop of resistivity following the excitations beyond a threshold fluence. From the figure, the threshold is estimated near ∼ 2.5 mJ/cm2, Each subsequent shot will drive the resistance drop into a new steady step value until another shot incident. With the shot number increasing, the resistivity drops quasi-exponentially to a local minimum value. Higher density excitations will drive the resistance to smaller values till reaching a minimum constant. This step-by-step type behavior is closely coincident with the phase percolation scenario, as it was found in transition metal oxide La2/3Ca1/3MnO3 [9]. That means the resistance change does not simply arise from the number of absorbed photons. The writing pulses can induce a range of resistivity values, indicating the multi-phase coexistence of pristine and PI-LR phases, dependent on the shot number and intensity. Thus, we can imagine an exponential growth of the excited volume fraction of PI-LR phase with increasing shot number or intensity. We also found that the same PI-LR phase can be driven by laser pulses with longer pulse duration of 100 fs but at higher fluence, suggesting that the pulse field plays a dominant role in the PI effect. We remark that, to achieve the PI-LR state, the delicate substrate environment has to be taken into account. In Ta2NiSe5, the PI-LR state is always detected when the sample’s thickness is thinner than the light penetration depth. However, for thicker bulk samples the 2.0-and 3.7-THz modes are usually recovered after the writing pulses ceased (see Supplementary Figure 6), with a high resistance state being kept as well. The observations are understandable. For the exfoliated thin flake sample transferred on the substrate with thickness smaller than the optical penetration depth, once the PI-LR state is triggered, the underlying lattice would be completely transformed into a new phase that is weakly affected from the substrate. Whereas for a thick bulk sample, the underneath unexcited or incompletely transited lattice may inevitably exert a restoring strain to the upper excited layers, resulting in a possible recovery of the 2.0-and 3.7-THz modes when the pump fluence is below the damage threshold. Switching of electronic properties more easily in thinner superlattice samples by femtosecond laser pulse has also been reported in other phase change materials, for example, in Ge2Sb2Te5 in terms of strain effect [38]. For relative weak pulse energies, the temporal photoresponse spectra behave in a different rule. Fig. 4b presents the coherent phonon spectra in a pristine nano-thickness sample with an increase of pump fluence up to 2.2 mJ/cm2, The original oscillations in the time domain is presented in Supplementary Figure 7. Initially at weak pump fluence, the primitive five coherent modes are recorded. With the pump fluence increasing, the 2.0-THz mode tends degraded gradually and disappear completely above ∼ 550 µJ/ cm2, while the 3.7-and 4.0 THz modes get merged into a broad one, suggesting a PI structural phase transition. However, this transition is transient, it doesn’t yield a permanent reduction of DC resistivity and the coherent phonon spectra will recover to the initial state when the pump fluence is tuned back to small values (See Supplementary Figure 8). With the excitation intensity increasing, the three low frequency modes show a clear red shift and an asymmetric broadening to the lower-frequency side, the phenomena can be ascribed to coherent anharmonic oscillations under high density photoexcitations [39, 40]. The red shift of the peak frequency can be considered as an electronic softening effect relating with the transient photoexcited carriers. While for the asymmetric broadening modes, the effect is typically referred as the nonthermal modification to the lattice potential-energy surface at high intensity excitations [41, 42]. This broadening modification is a transient effect that only lasts for a short time before thermal equilibrium is attained between electrons and lattice background. We performed the time-domain trace analysis for the modification effect. Fig. 4c shows the Fourier transformed spectra obtained after cutting the initial time delay Td at two different positions for a pump fluence 1.5 mJ/ cm2. As can be seen, if we perform Fourier transformation of the pump-probe waveform after the first 2 ps (Td = 2 ps), the red shift and broadening of the phonon modes approach to disappear. The inset shows

FIG. 4. Pulse fluence threshold and Fluence dependent PI state properties. a, Shot-to-shot resistivity change of Ta2NiSe5 ultrathin sample at 50 K. Each shot includes a sequence of 5 pulses. b, Fluence dependent coherent phonon spectra for moderate fluences at 50 K. The energy of probe pulse kept constant in a small value of ∼ 3 µJ/ cm2. The dashed are guidelines. c, the red-shift effect of the 3.0-THz mode as a function of the time delay Td for incident 1.5 mJ/ cm2. The inset shows an exponential approximation recovery to the initial frequency, where the error bar is the standard deviation evaluated when extracting the peak frequency from the coherent phonon spectra.

the change of 3 THz mode of Fourier transformed spectra after cutting different time delays Td. A more detailed analysis of the transient phonon spectra as a function of cutting time delay Td at various excitation intensity is presented in Supplementary Figure 9. We now discuss implications of the results. First, our experiments clearly demonstrate multi-staged PI phase transitions. At small excitation fluence, the system is triggered to a transient or short-lived state accompanied by a structural change. Upon increasing the pump fluence beyond a threshold, the sample is driven to a stable phase showing high conductivity. This phase has a crystal structure different from the original one with observable sliding of Ta-atom, which is also robust against thermal excitations up to at least 350 K. Obviously, this phase must correspond to one of the pronounced minima in the free energy landscape. The energy barrier is high enough to prevent the system from thermally recovering to the pristine phase. In principle, the PI atomic vacancies might provide such a stabilized energy barrier, like a case of irreversible IMTs from insulating 2H-MoTe2 to metallic 1T’-MoTe2 induced by high-intensity photo-irradiations[43]. However, this possibility can be excluded here in Ta2NiSe5 system, because our atomically resolved TEM experiments indicate a homogeneous atom arrangement with negligible vacancy or disorder in the PI state of Ta2NiSe5. Therefore, the PI multistage phase transitions in Ta2NiSe5 must be a result of self-organization through different trajectories from photoexcited states. Second, our measurements provide useful information about exciton formation in Ta2NiSe5 system. The binding energy of exciton appears to be very sensitive to the lattice distortion. The phase of Ta2NiSe5 below Tc = 326 K has been widely considered as a prototype excitonic insulator [17–19, 44–46]. As is known, for a small gap semiconductor, an excitonic instability occurs only if the exciton binding energy is larger than the electronic band energy gap. Then, the system would form a new phase through spontaneous formation and Bose condensation of excitons in which the exciton binding energy would become smaller than the band energy gap. Our measurement reveals the presence of a non-thermal stable phase with subtle difference in structure relative to the pristine Ta2NiSe5. The absence of a phase transition suggests very small exciton binding energy in this phase. Therefore, the subtle structural change could induce a large change in the exciton binding energy. Nevertheless, the driving mechanism for the PI phase transitions may not be directly linked to the exciton interaction. Although it is not easy to solve the issue conclusively at present stage, a comparative study to the pure semiconductor Ta2NiS5, an isostructural sister compound to Ta2NiSe5 without exhibiting exciton insulator ground state [17], would be particularly helpful. To this end, we performed similar measurement on Ta2NiS5. Indeed, we observed similar PI-LR state in Ta2NiS5 (See Supplementary Figure 10). The results indicate that the exciton interaction is not a prerequisite for the formation of PI-LR state. Finally, our work reveals a significantly different photoinduced effect from reported work for the compound. The issue of photoinduced phase transition in Ta2NiSe5 has been heavily debated with controversial results. Several reports denied a possible of photoinduced phase transition [23, 47, 48], however some recent reports support that transition [24, 25, 49]. For example, in earlier tr-ARPES and ultrafast pump-probe studies on Ta2NiSe5 [23, 47], Mor et al. reported that the photon absorption 2 saturated at a critical fluence of Fc =0.2 mJ/cm upon photoexcitation, and above the saturation fluence the band energy gaps widened, reflecting an increase of the exciton condensate density in the EI phase in Ta2NiSe5 [23]. In accord with the tr-ARPES results, their pump-probe measurement also indicated a saturation phenomenon above a critical excitation fluence [47]. Based on those results, they argued that Ta2NiSe5 would exhibit a blocking mechanism when pumped in the near-infrared regime, preventing a nonthermal structural phase transition. Apparently, our experimental observation of the PI-LR phase in Ta2NiSe5 is in sharp contrast to those results. The strong NIR laser pulses can definitely induce a structural phase transition. It is possible that the peak electric field of laser pulse used by Mor et al. was too small to drive a transition (See Supplementary Table 1). Another possibility is that the samples they used in their tr-ARPES study are thick bulk samples. As we mentioned above, the sharp switching is achieved currently only in thin flake samples. We note that more recent tr-ARPES measurements using an 800 nm amplified laser system also indicated a photo-induced transient insulator-to-metal transitions in Ta2NiSe5, though a stable state was not observed [24, 49]. Those results are consistent with our observation of transient structural change at moderate excitations which can recovers to the initial state once the pump fluence is tuned back to a small value. Our present work brings new perspective to the control and manipulation of physical properties in Ta2NiSe5. The realization of ultrafast photoinduced irreversible phase transition with a dramatic resistivity change will open up new prospects for designing novel functional photonic and electronic devices.

Methods Sample growth: The Ta2NiSe5 single crystals were synthesized by traditional chemical vapor transportation (CVT) technology. The quartz tube loaded with the elements (99.99% Ta power, 99.99% Ni power and 99.999% Se pellets) in a stoichiometric ratio together with additional iodine (2 % of total ◦ mass) as the transport agent was vacuumed and flamed-sealed. Then the quartz tube was heated at 750 C ◦ -880 C for 10 days in a two-zone furnace, followed by a natural cooling process. The obtained single crystals (normally in size of 5×2×0.1 mm) were checked by the resistance, electron diffraction energy spectrum and infrared spectroscopy. Device fabrication and temporal spectra of optical excitations: The thin Ta2NiSe5 samples used in the experiments were mechanically exfoliated from the single crystals using the sticky tape and then transferred onto the cleaned sapphire substrates. The thickness of the films (10-50 nm in thickness) was measured by the KLA-Tencor’s P-6 stylus profiler. The electrodes were fabricated by ultra-violet lithography using the bilayer photoresists (PMGI and AZ1500) and followed by the electron beam evaporation deposited with 10 nm Ti and 60 nm Au successively. The typical image of the devices taken from the Olympus BX51 Microscope. The temporal resistance under photoexcitations were in-situ recorded by the standard four-electrode method using a Keithley 2400 Source-Meter. A sample holder of the Liquid-He flow optical cryostat (Oxford Company) was used. A Ti:Sapphire Amplifier (800 nm, 35 fs, 1 kHz-repetition rate) was supplied for both the writing and pump-probe pulses. The shot-number of the writing pulse was controlled by a shutter. The temperature dependence of resistance R(T) were also measured before and after the optical experiments by PPMS (Quantum Design Company). For optical pump-probe experiment, the two-electrodes geometry are normally used. The diameter of the writing pulse spot was set normally in the range of 0.2 -1 mm, fully covering the channel area between the electrodes. For pure optical pump-probe measurement, the diameter of the pump spot was 200 300 µm, more than two times larger than the one of probe pulse (90 µm). The photoinduced phase transition can be achieved by both the front and backside incident. TEM experiment: A mechanical exfoliation method with scotch tape was used to obtain the Ta2NiSe5 nanoflakes, which were placed on 1000 mesh copper grids for transmission electron microscopy measurements. For the laser-treatment, the nanoflakes were excited at temperature 50 K by a set of 30 2 pulses (800 nm, 35 fs, ∼3.5 mJ/cm ). In situ heating TEM observations were performed on a JEOL JEM-2100F microscope operated at 200 kV. The use of in situ heating TEM holder (model 652, Gatan Inc.) allows specimen temperature to be varied from room temperature to 750 K. Atomically resolved high-angle annual dark-field (HAADF) images were obtained using a probe-corrected JEOL ARM-200C TEM. TEM bright-field imaging and selected area diffraction patterns were recorded using both ARM-200C and JEM-2100F microscopes.

Acknowledgements We would like to thank Profs. H. X. Yang and Z. Xu for their help in experiment and useful discussions. This work was supported by National Natural Science Foundation of China (No. 11888101, 11974019), the National Key Research and Development Program of China (No. 2017YFA0302904, 2017YFA0300303).

Author contributions D.W., Q.M.L., L.Y.S., Z.X.W.. performed crystal growth, device fabrication, and electrical property measurements. Q.M.L., D.W., L.Y.S., Z.X.W., S.J.Z., T.L., T.C.H., T.D., N.L.W., performed optical experiments and analyzed the data. Z.A.L., H.F.T., J.Q.L., performed TEM experiments and analyzed the data. D.W., Q.M.L., N.L.W. wrote the manuscript, with contributions from all the authors. N.L.W. and D.W. conceived the project. All authors have approved the final version of the manuscript. These authors contributed equally: Q.M.Liu, D.Wu, Z.A.Li.

References: [1] Claudio Giannetti, Massimo Capone, Daniele Fausti, Michele Fabrizio, Fulvio Parmigiani, and Dragan Mihailovic, “Ultrafast optical spectroscopy of strongly correlated materials and high-temperature superconductors: a non-equilibrium approach,” Advances in Physics 65, 58–238 (2016). [2] K. Nasu (Ed.), “Photoinduced phase transitions,” World Scientific,Singapore (2004). [3] N. Takubo, Y. Ogimoto, M. Nakamura, H. Tamaru, M. Izumi, and K. Miyano, “Persistent and reversible all-optical phase control in a manganite thin film,” Physical Review Letters 95, 017404 (2005). [4] Manfred Fiebig, Kenjiro Miyano, Yoshinori Tokura, and Yasuhide Tomioka, “Visualization of the local insulator-metal transition in Pr0.7Ca0.3MnO3,” Science 280, 1925–1928 (1998). [5] D. Fausti, R. I. Tobey, N. Dean, S. Kaiser, A. Dienst, M. C. Hoffmann, S. Pyon, T. Takayama, H. Takagi, and A. Cavalleri, “Light-induced in a stripe-ordered cuprate,” Science 331, 189–191 (2011). [6] Hirohiko Ichikawa, Shunsuke Nozawa, Tokushi Sato, Ayana Tomita, Kouhei Ichiyanagi, Matthieu Chollet, Laurent Guerin, Nicky Dean, Andrea Cavalleri, Shin-ichi Adachi, Taka-hisa Arima, Hiroshi Sawa, Yasushi Ogimoto, Masao Nakamura, Ryo Tamaki, Kenjiro Miyano, and Shin-ya Koshihara, “Transient photoinduced ’hidden’ phase in a manganite,” Nature materials 10, 101–105 (2011). [7] Ryusuke Matsunaga, Naoto Tsuji, Hiroyuki Fujita, Arata Sugioka, Kazumasa Makise, Yoshinori Uzawa, Hirotaka Terai, Zhen Wang, Hideo Aoki, and Ryo Shimano, “Light-induced collective pseudospin precession resonating with higgs mode in a superconductor,” Science 345, 1145–1149 (2014). [8] L. Stojchevska, I. Vaskivskyi, T. Mertelj, P. Kusar, D. Svetin, S. Brazovskii, and D. Mihailovic, “Ultrafast switching to a stable hidden quantum state in an electronic crystal,” Science 344, 177–180 (2014). [9] Jingdi Zhang, Xuelian Tan, Mengkun Liu, S. W. Teitelbaum, K. W. Post, Feng Jin, K. A. Nelson, D. N. Basov, Wenbin Wu, and R. D. Averitt, “Cooperative photoinduced metastable phase control in strained manganite films,” Nature Materials 15, 956–960 (2016). [10] D N Basov, R D Averitt, and D Hsieh, “Towards properties on demand in quantum materials,” Nature Materials 16, 1077–1088 (2017). [11] T. F. Nova, A. S. Disa, M. Fechner, and A. Cavalleri, “Metastable ferroelectricity in optically strained SrTiO3,” Science 364, 1075–1079 (2019), arXiv:1812.10560. [12] Xian Li, Tian Qiu, Jiahao Zhang, Edoardo Baldini, Jian Lu, Andrew M. Rappe, and Keith A. Nelson, “Terahertz field–induced ferroelectricity in quantum paraelectric SrTiO3,” Science 364, 1079–1082 (2019). [13] A. S. McLeod, Jingdi Zhang, M. Q. Gu, F. Jin, G. Zhang, K. W. Post, X. G. Zhao, A. J. Millis, W. B. Wu, J. M. Rondinelli, R. D. Averitt, and D. N. Basov, “Multi-messenger nanoprobes of hidden magnetism in a strained manganite,” Nature Materials 19, 397–404 (2020), arXiv:1910.10361. [14] Bradley J. Siwick Vance R. Morrison, Robert. P. Chatelain, Kunal L. Tiwari, Ali Hendaoui, Andrew Bruh´acs, Mohamed Chaker, “A photoinduced metal-like phase of monoclinic VO2 revealed by ultrafast electron diffraction,” Science 346, 445 (2014). [15] Cooke, and Bradley J. Siwick, “How optical excitation controls the structure and properties of vanadium dioxide,” Proceedings of the National Academy of Sciences U. S. A 116, 450–455 (2019), arXiv:1805.01430. [16] Martin R. Otto, Laurent P. Ren´Emond, Mohamed Chaker, David G. e de Cotret, David A. Valverde-Chavez, Kunal L. Tiwari, Nicolas Luciana Vidas, Daniel Schick, El´ıas Mart´ınez, Daniel Perez-Salinas, Alberto Ramos-´Alvarez, Simon Cichy, Sergi Batlle-Porro, Allan S. Johnson, Kent A. Hallman, Richard F. Haglund, and Simon Wall, “Does VO2 Host a Transient Monoclinic Metallic Phase?” Physical Review X 10, 031047 (2020). [17] Steven A. Sunshine and James A. Ibers, “Structure and Physical Properties of the New Layered Ternary Chalcogenides Ta2NiS5 and Ta2NiSe5 ,” Inorganic Chemistry 24, 3611–3614 (1985). [18] Y. Wakisaka, T. Sudayama, K. Takubo, T. Mizokawa, M. Arita, H. Namatame, M. Taniguchi, N. Katayama, M. Nohara, and H. Takagi, “Excitonic Insulator State in Ta2NiSe5 Probed by Photoemission Spectroscopy,” Phys. Rev. Lett. 103, 026402 (2009). [19] T. Kaneko, T. Toriyama, T. Konishi, and Y. Ohta, “Orthorhombic-to-monoclinic phase transition of Ta2NiSe5 induced by the Bose-Einstein condensation of excitons,” Phys. Rev. B 87, 035121 (2013). [20] H. J. Zeiger, J. Vidal, T. K. Cheng, E. P. Ippen, G. Dresselhaus, and M. S. Dresselhaus, “Theory for displacive excitation of coherent phonons,” Phys. Rev. B 45, 768–778 (1992). [21] Daniel Werdehausen, Tomohiro Takayama, Marc H¨oppner, Gelon Albrecht, Andreas W. Rost, Yangfan Lu, Dirk Manske, Hidenori Takagi, and Stefan Kaiser, “Coherent order parameter oscillations in the ground state of the excitonic insulator Ta2NiSe5,” Science Advances 4, 1–8 (2018). [22] O. V. Misochko, Muneaki Hase, K. Ishioka, and M. Kitajima, “Observation of an amplitude collapse and revival of chirped coherent phonons in bismuth,” Phys. Rev. Lett. 92, 197401 (2004). [23] Selene Mor, Marc Herzog, Denis Goleˇz, Philipp Werner, Martin Eckstein, Naoyuki Katayama, Minoru Nohara, Hide Takagi, Takashi Mizokawa, Claude Monney, and Julia St¨ahler, “Ultrafast electronic band gap control in an excitonic insulator,” Phys. Rev. Lett. 119, 086401 (2017). [24] Kozo Okazaki, Yu Ogawa, Takeshi Suzuki, Takashi Yamamoto, Takashi Someya, Shoya Michimae, Mari Watanabe, Yangfan Lu, Minoru Nohara, Hidenori Takagi, Naoyuki Katayama, Hiroshi Sawa, Masami Fujisawa, Teruto Kanai, Nobuhisa Ishii, Jiro Itatani, Takashi Mizokawa, and Shik Shin, “Photo-induced semimetallic states realised in electron hole coupled insulators,” Nature Communications 9, 1–6 (2018). [25] Takeshi Suzuki, Yasushi Shinohara, Yangfan Lu, Mari Watanabe, Jiadi Xu, Kenichi L. Ishikawa, Hide Takagi, Minoru Nohara, Naoyuki Katayama, Hiroshi Sawa, Masami Fujisawa, Teruto Kanai, Jiro Itatani, Takashi Mizokawa, Shik Shin, and Kozo Okazaki, “Detecting electron-phonon couplings during photo-induced phase transition,” (2020), arXiv:2002.10037. [26] H. J. Zeiger, J. Vidal, T. K. Cheng, E. P. Ippen, G. Dresselhaus, and M. S. Dresselhaus, “Theory for displacive excitation of coherent phonons,” Physical Review B 45, 768–778 (1992). [27] F. J. Di Salvo, C. H. Chen, R. M. Fleming, J. V. Waszczak, R. G. Dunn, S. A. Sunshine, and James A. Ibers, “Physical and structural properties of the new layered compounds Ta2NiS5 and Ta2NiSe5,” Journal of The Less-Common Metals 116, 51–61 (1986). [28] Jian Yan, Ruichun Xiao, Xuan Luo, Hongyan Lv, Ranran Zhang, Yan Sun, Peng Tong, Wenjian Lu, Wenhai Song, Xuebin Zhu, and Yuping Sun, “Strong Electron−phonon Coupling in the Excitonic Insulator Ta2NiSe5,” Inorganic Chemistry 58, 9036–9042 (2019). [29] Alaska Subedi, “Orthorhombic-to-monoclinic transition in Ta2NiSe5 due to a zone-center optical phonon instability,” Phys. Rev. Materials 4, 083601 (2020). [30] Min-Jae Kim, Armin Schulz, Tomohiro Takayama, Masahiko Isobe, Hidenori Takagi, and Stefan Kaiser, “Phononic soft mode behavior and a strong electronic background across the structural phase transition in the excitonic insulator Ta2NiSe5,” Phys. Rev. Research 2, 042039 (2020). [31] Pavel A. Volkov, Mai Ye, Himanshu Lohani, Irena Feldman, Amit Kanigel, Kristjan Haule, and Girsh Blumberg, “Critical charge fluctuations and quantum coherent state in excitonic insulator Ta2NiSe5,” arXiv (2020), 2007.07344. [32] Kwangrae Kim, Hoon Kim, Jonghwan Kim, Changil Kwon, Jun Sung Kim, and B. J. Kim, “Direct observation of excitonic instability in Ta2NiSe5,” (2020), arXiv:2007.08212. [33] Edbert J. Sie, Clara M. Nyby, C. D. Pemmaraju, Su Ji Park, Xiaozhe Shen, Jie Yang, Matthias C. Hoffmann, B. K. Ofori-Okai, Renkai Li, Alexander H. Reid, Stephen Weathersby, Ehren Mannebach, Nathan Finney, Daniel Rhodes, Daniel Chenet, Abhinandan Antony, Luis Balicas, James Hone, Thomas P. Devereaux, Tony F. Heinz, Xijie Wang, and Aaron M. Lindenberg, “An ultrafast symmetry switch in a weyl semimetal,” Nature 565, 61 (2019). [34] M. Y. Zhang, Z. X. Wang, Y. N. Li, L. Y. Shi, D. Wu, T. Lin, S. J. Zhang, Y. Q. Liu, Q. M. Liu, J. Wang, T. Dong, and N. L. Wang, “Light-Induced Subpicosecond Lattice Symmetry Switch in MoTe2,” Phys. Rev. X 9, 021036 (2019). [35] Akitoshi Nakano, Kento Sugawara, Shinya Tamura, Naoyuki Katayama, Kazuyuki Matsubayashi, Taku Okada, Yoshiya Uwatoko, Kouji Munakata, Akiko Nakao, Hajime Sagayama, Reiji Kumai, Kunihisa Sugimoto, Naoyuki Maejima, Akihiko Machida, Tetsu Watanuki, and Hiroshi Sawa, “Pressure-induced coherent sliding-layer transition in the excitonic insulator Ta2NiSe5,” IUCrJ 5, 158–165 (2018). [36] Quirin Stahl, Maximilian Kusch, Florian Heinsch, Gaston Garbarino, Norman Kretzschmar, Kerstin Hanff, Kai Rossnagel, Jochen Geck, and Tobias Ritschel, “Collapse of layer dimerization in the photo-induced hidden state of 1T-TaS2,” Nature Communications 11, 1–7 (2020) [37] Sung-Hoon Lee, Jung Suk Goh, and Doohee Cho, “Origin of the insulating phase and first-order metal-insulator transition in 1t−tas2,” Phys. Rev. Lett. 122, 106404 (2019). [38] Janne Kalikka, Xilin Zhou, Eric Dilcher, Simon Wall, Ju Li, and Robert E. Simpson, “Strain-engineered diffusive atomic switching in two-dimensional crystals,” Nature Communications 7 (2016), 10.1038/ncomms11983. [39] P. Tangney and S. Fahy, “Calculations of theA1 phonon frequency in photoexcited tellurium,” Physical Review Letters 82, 4340 (1999). [40] Muneaki Hase, Masahiro Kitajima, Shin-ichi Nakashima, and Kohji Mizoguchi, “Dynamics of coherent anharmonic phonons in bismuth using high density photoexcitation,” Phys. Rev. Lett. 88, 067401 (2002). [41] Roman Yusupov, Tomaz Mertelj, Viktor V Kabanov, Serguei Brazovskii, Primoz Kusar, Jiun-Haw Chu, Ian R Fisher, and Dragan Mihailovic, “Coherent dynamics of macroscopic electronic order through a symmetry breaking transition,” Nature Physics 6, 681(2010). [42] German Sciaini, Maher Harb, Sergei G Kruglik, Thomas Payer, Christoph T Hebeisen, Frank-J. Meyer zu Heringdorf, Mariko Yamaguchi, Michael Horn-von Hoegen, Ralph Ernstorfer, and R J Dwayne Miller, “Electronic acceleration of atomic motions and disordering in bismuth,” Nature 458, 56 (2009). [43] Suyeon Cho, Sera Kim, Jung Ho Kim, Jiong Zhao, Jinbong Seok, Dong Hoon Keum, Jaeyoon Baik, Duk-Hyun Choe, K. J. Chang, Kazu Suenaga, Sung Wng Kim, Young Hee Lee, and Heejun Yang, “Phase patterning for ohmic homojunction contact in MoTe2,” Science 349, 625–628 (2015). [44] Koudai Sugimoto, Satoshi Nishimoto, Tatsuya Kaneko, and Yukinori Ohta, “Strong Coupling Nature of the Excitonic Insulator State in Ta2NiSe5 ,” Phys. Rev. Lett. 120, 247602 (2018). [45] Giacomo Mazza, Malte R¨osner, Lukas Windg¨atter, Simone Latini, Hannes H¨ubener, Andrew J. Millis, Angel Rubio, and Antoine Georges, “Nature of Symmetry Breaking at the Excitonic Insulator Transition: Ta2NiSe5,” Phys. Rev. Lett. 124, 197601 (2020). [46] Matthew D. Watson, Igor Markovievre, Michael Merz, Amir A. Haghighirad, and Philip D. C. Edgar Abarca Morales, Patrick Le FKing, “Band hybridization at the semimetal-semiconductor transition of Ta2NiSe5 enabled by mirror-symmetry breaking,” Phys. Rev. Research 2, 013236 (2020). [47] Selene Mor, Marc Herzog, Johannes Noack, Naoyuki Katayama, Minoru Nohara, Hide Takagi, Annette Trunschke, Takashi Mizokawa, Claude Monney, and Julia Stahler, “Inhibition of the photoinduced structural phase transition in the excitonic insulator Ta2NiSe5,” Physical Review B 97, 115154 (2018). [48] Daniel Werdehausen, Tomohiro Takayama, Gelon Albrecht, Yangfan Lu, Hidenori Takagi, and Stefan Kaiser, “Photo-excited dynamics in the excitonic insulator Ta2NiSe5,” Journal of Physics Condensed Matter 30 (2018), 10.1088/1361-648X/aacd76. [49] Paolo Andrich, Hope M. Bretscher, Prachi Telang, Anupam Singh, Luminita Harnaga, A. K. Sood, and Akshay Rao, “Ultrafast melting and recovery of collective order in the excitonic insulator Ta2NiSe5,” arXiv , 3–10 (2020), arXiv:2007.03368.