Persistent order due to transiently enhanced nesting in an electronically excited charge density wave.

Non-equilibrium conditions may lead to novel properties of materials with broken symmetry ground states not accessible in equilibrium as vividly demonstrated by non-linearly driven mid-infrared active phonon excitation. Potential energy surfaces of electronically excited states also allow to direct nuclear motion, but relaxation of the excess energy typically excites fluctuations leading to a reduced or even vanishing order parameter as characterized by an electronic energy gap. Here, using femtosecond time- and angle-resolved photoemission spectroscopy, we demonstrate a tendency towards transient stabilization of a charge density wave after near-infrared excitation, counteracting the suppression of order in the non-equilibrium state. Analysis of the dynamic electronic structure reveals a remaining energy gap in a highly excited transient state. Our observation can be explained by a competition between fluctuations in the electronically excited state, which tend to reduce order, and transiently enhanced Fermi surface nesting stabilizing the order.

Ordered states in condensed matter emerge in thermodynamic equilibrium below a critical temperature Tc if competing thermal fluctuations are weak compared with the energy gain driving the order. This delicate balance depends critically on the dimensionality of a system, where low-dimensional structures introduce instabilities towards broken symmetry ground states, but on the same footing enhance fluctuations. Charge density waves (CDWs) form due to cooperative effects of electrons and the crystal lattice and represent widely studied model systems of such broken symmetry ground states. Under equilibrium conditions CDW formation can be described phenomenologically by effective mean field models1, but on a microscopic level a description is far from trivial, because the electronic and lattice contributions can hardly be disentangled due to their intrinsic coupling23. Optical excitation of CDW systems into non-equilibrium states faster than the characteristic time scales of electrons or phonons breaks this link and has opened opportunities to investigate the interactions underlying CDW formation4567891011. In addition, tailored excitations of specific resonances in the mid-infrared have been reported to lead to novel, non-equilibrium material properties that are not accessible in thermal equilibrium1213141516.

The formation of a CDW is driven by an instability of the electronic system to a spatially periodic perturbation. Particularly in quasi one-dimensional systems parallel parts of the Fermi surface (FS) are nested by an ordering vector CDW, and lead to a divergence of the electronic Lindhard susceptibility1. Electron-phonon (e-ph) coupling imprints this ordering tendency on the lattice and freezes a soft phonon mode into a periodic lattice distortion, leading to a complex many-body problem. While this is a widely considered explanation for CDW formation117, momentum-dependent e-ph coupling may modify this picture31819. This coupled charge and lattice periodicity 1/qCDW creates an energy gap 2Δ in the electronic structure due to Bragg scattering at the modulated charge density. Thereby the periodic lattice distortion and the CDW are stabilized by the strength of e-ph coupling and the diverging susceptibility, that is, the FS topology. Optical excitation creates fluctuations of the electron density, which on average reduce the electronic charge modulation and 2Δ. Simultaneously, coherent lattice dynamics may drive the system towards a high-symmetry state and decrease the lattice distortion, further reducing the CDW modulation and 2Δ. Such ultrafast suppression in 2Δ has been widely observed in previous experiments employing femtosecond time- and angle-resolved photoemission spectroscopy (trARPES) after femtosecond laser excitation5620 in agreement with the theory21, which emphasizes the impact of optically induced fluctuations on the interaction leading to CDW formation. However, an enhanced driving force for CDW formation maintaining order in the non-equilibrium state was so far not observed.

Here we demonstrate a hindered melting of the CDW state after femtosecond laser excitation in the CDW model system RTe3 (R=Dy, Ho) through a residual gap 2Δ in a transient optically excited state. This residual 2Δ is fluence independent and is a consequence of improved FS nesting conditions in the excited state, for which we provide evidence through a complete analysis of 2Δ as a function of two independent momentum directions k, k, binding energy E and pump-probe delay t. This improved nesting increases the interactions underlying CDW formation which compete with incoherent fluctuations induced by the excitation. This surprising observation suggests a pathway to control material properties under non-equilibrium conditions relevant for a wide range of material classes with broken symmetry ground states.

Results

Static FS of RTe3

We begin by analyzing the electronic structure of HoTe3 in thermal equilibrium using high-resolution laser-ARPES at hν=7 eV, which is representative for the series of rare-earth tritellurides. Its quasi-2D electronic structure is determined by square nets of Te-planes, which are stacked between buckled RTe layers that allow to tune the CDW temperature and gap size through chemical pressure22. Figure 1b (left) shows the continuous, metallic FS of HoTe3 for T>Tc, which agrees well with tight binding (TB) calculations (solid lines) of the Te planes taking into account the overlap of Te-5p and 5p orbitals (refs 17, 18), for details see Supplementary Note 1 and Supplementary Fig. 1. The curvature of the diamond-shaped FS in the (k, k)-plane is determined by the ratio of orbital overlap along and perpendicular to the chains of p/p-orbitals in the TB model described by t⊥/t|| (ref. 17), see Fig. 1a. For T<Tc (Fig. 1b, right) CDW formation results in shadow bands (green dashed lines), which are translated along k by the nesting vector qCDW≈0.7c* (ref. 17) (Fig. 1a), and a CDW energy gap 2Δ opens in the vicinity of (k, k)=(0.15, 0.30) Å− where the FS is nested due to overlapping main and shadow bands (see Supplementary Fig. 1). With decreasing k the nesting gradually weakens and residual metallic pockets appear on the FS, where spectral weight is transferred into the shadow bands. Note that the imperfect nesting results in a k-dependent shift of the gap centre along the FS, while the full gap value 2Δ remains constant (see Supplementary Fig. 1). In the heavier members of RTe3, starting from TbTe3, also a second perpendicular CDW transition occurs along a* at a lower T (refs 18, 22). Here we concentrate on the first CDW transition with the larger 2Δ, which is well separated on the FS.

Figure 1

Time-dependent Fermi surface map.

(a) Experimental scheme and TB model. The Te-5p orbitals in the Te-layers of RTe3 (bottom) are coupled by t|| and t⊥ describing interactions parallel and perpendicular to the Te chains and lead to a diamond-shaped Fermi surface in reciprocal space, which we excess by our position-sensitive time-of-flight photoelectron spectrometer (pTOF) (top). Nesting leads to shadow bands translated along qCDW (dashed lines). (b) Static FS of HoTe3 for T=300 K>Tc (left) and T=180 K

Transient FS and band structure

Using femtosecond trARPES, see Fig. 1a, we investigate the electronic structure after fs laser excitation in DyTe3, which is very similar to HoTe3 (refs 17, 18). Figure 1c–e shows the gapped FS region (red box in Fig. 1b) probed by hνprobe=6.0 eV after optical excitation with hνpump=1.5 eV pump pulses at an absorbed fluence of F=0.27 mJ cm−2 for selected time delays. A movie of the transient FS is also available (Supplementary Movie 1). We observe filling of the gapped region starting from the metallic pocket. At 200 fs a nearly ungapped FS is found, which is accompanied by a shift of spectral weight back to the TB main band. This provides direct evidence for an optically driven gap-closing transition. At a first glance the transient state at 200 fs is very similar to the situation at T>Tc and consistent with the TB model623.

Our electron time-of-flight spectrometer24 allows extracting spectra along arbitrary in-plane momenta and we first analyze the trARPES intensity I(E, k1, t) with k1 perpendicular to the FS contour, see Fig. 1c and Supplementary Movie 2. Figure 2 shows an occupied band with CDW gap before excitation. At t=0 fs the unoccupied band at the top of the CDW gap is populated, which enables a direct determination of the full transient 2Δ. Both bands and the gap energy agree well with the interacting TB model (see Supplementary Note 1). Simultaneously, the occupied band is depleted, which implies hole excitation. With increasing pump-probe delay t the lower and upper bands move towards each other as the CDW order recedes and 2Δ is reduced. At t=100 fs the lower band reaches EF and at t=200 fs, E(k1) is well described by the TB main band in agreement with the FS analysis above.

Figure 2

Time-dependent CDW band structure.

(a–d) trARPES data along k1 for various pump-probe delays. Main and shadow TB bands are shown as solid and dashed lines, and symbols show lower (blue) and upper (red) band dispersion of the TB model including CDW gap, where the symbol sizes represent the spectral weight. (e) Spectra integrated around the region marked in b as a function of pump-probe delay. Red and blue symbols are peak positions of upper and lower CDW bands, respectively, obtained from a fitting procedure. (f) Peak shift of upper and lower band. Shown is the absolute value of the shift. Error bars represent 95% confidence intervals of the fits.

In Fig. 2e we plot trARPES spectra at kF and clearly observe periodic variations of 2Δ. The peak positions of the lower and upper CDW bands are determined by Lorentzian line fits of I(E, t), for details see Supplementary Note 2 and Supplementary Fig. 2. Before excitation, the lower CDW band occurs at E−EF=−0.12(1) eV (blue markers). Beginning at t=0 fs, we monitor the upper CDW band at E−EF=0.34(3) eV (red markers) and determine 2Δ=0.46(4) eV from the energy difference of the upper and lower CDW peaks. With increasing t both lower and upper CDW bands shift towards EF. At t=200 fs the lower (upper) band reaches its maximum (minimum) and 2Δ is smallest before the bands shift away from EF again at larger delays. Remarkably, two distinct peaks are observed in the spectra at all times (Supplementary Fig. 2), providing direct evidence for 2Δ remaining finite in the transient metallic state.

Figure 3 shows 2Δ(t) (see Supplementary Note 3 and Supplementary Fig. 3) for various fluences F. For F≤0.12 mJ cm−2 a relatively small decrease is followed by an oscillation with a period of ∼0.5 ps. This coherent response in a weakly perturbative regime lasts for several picoseconds and corresponds to the excitation of the CDW amplitude mode162325262728. For higher fluence in a strongly perturbative regime F≥0.18 mJ cm−2, corresponding to ∼220 meV per unit cell (see Supplementary Note 4), which is comparable to the electronic energy gain upon CDW formation δE∼250 meV per unit cell estimated from the TB model (see Supplementary Note 1 and Supplementary Fig. 1), we find a stronger initial decrease of 2Δ, followed by a recovery of 2Δ overlaid by anharmonic and strongly damped oscillations. These oscillations are a fingerprint of the coherent atomic rearrangements during the collapse and recovery of the CDW order in the transient potential energy surface. The minimal value of 2Δ saturates at 170 meV independent of fluence for F>0.12 mJ cm−2. This previously unobserved residual 2Δ in the transient metallic state demonstrates the persistence of electronic order in a highly perturbed system, in contrast to equilibrium conditions at T>Tc (ref. 18). This incomplete suppression of the CDW state is an unexpected result because at the highest fluence an excess energy of almost twice the energy gained by gapping the electronic structure is injected and still the gap persists (Fig. 3).

Figure 3

Time-dependent CDW gap size.

Time-dependent CDW gap size 2Δ determined from the separation of upper and lower CDW band as a function of pump-probe delay for various fluences.

Asymmetric gap suppression

Our experimental approach further allows for the analysis of the symmetry of 2Δ. A closer look at the transient peak positions in Fig. 2e reveals a larger pump-induced shift of the upper than of the lower CDW peak, indicating an asymmetric decrease of 2Δ. From the relative shift of both bands in Fig. 2f we find a maximal shift of 110 meV for the lower band (blue), while the upper band (red) changes by 60% more up to 190 meV. The momentum dependence of I(E, k2, t) along k2 (see Fig. 1c) parallel to the FS contour is shown in Figs 4a,b for t=0 and 200 fs and exhibits lower and upper CDW bands dispersing along k2. A full movie is available as Supplementary Movie 3. This relates to a shift of the centre of 2Δ along the FS because 2Δ(k) is only centred around EF for perfect nesting17 near (k, k)=(0.17, 0.29) Å−1. To study the transient state, it is helpful to look at the pump-induced changes where the gap is minimal (t=200 fs). The pump-induced peak shifts (see Supplementary Note 5 and Supplementary Fig. 4) |E(k2)up,down(200 fs)−E(k2)up,down(0 fs)| are depicted in Fig. 4c and exhibit a pronounced dependence on the FS position: With decreasing k, the transient shift of the upper band decreases, while the shift of the lower band increases, leading to a smaller asymmetry of the gap reduction. Finally, at (k, k)∼(0.17, 0.29) Å−1 both peak shifts become equal and the CDW gap decreases symmetrically around its centre. These observations manifest in a transient reduction of the dispersion of the gap centre with k2, as seen in Fig. 4d.

Figure 4

Gap collapse along the Fermi surface.

(a,b) Electron dispersion along k2 for t=0 fs and t=200 fs, respectively. Markers denote energy positions of upper (red) and lower (blue) CDW bands. (c) Peak shift between t=0 fs and =200fs of upper and lower CDW band as function of position on the FS determined by fits, with error bars as 95% confidence intervals. Lines and shaded areas are linear fits and 95% confidence bands, respectively. (d) Centre of the CDW gap as a function of position on the FS at t=0 fs (open symbols) and at t=200 fs (filled symbols). Solid lines are linear fits. The inset shows the relative change of slope Δm/m as function of fluence. (e,f) Tight binding bands for t⊥=0.35 eV (e) and t⊥=0.25 eV (f). The grey areas marks the regions of good nesting where main (solid) and shadow (dashed) band lines overlap.

Discussion

Such a transient change of the gap centre indicates an ultrafast modification of the nesting condition. The nesting vector and the curvature of the bands (Fig. 1) determine the band crossing of main and shadow bands (Fig. 2) and thus 2Δ(k). Changes of qCDW alone can be ruled out to be driving an asymmetric closing of 2Δ since this would lead to a constant shift along the FS (see Supplementary Note 6 and Supplementary Fig. 5) and since time-resolved diffraction29 does not resolve such changes. In contrast, a transient modification of the electronic band dispersion, which can be parametrized with modified TB parameters, can explain the observed modification of the gap dispersion (Supplementary Fig. 5). The experimentally observed smaller gap dispersion can be reproduced by a reduction of t⊥ leading to a smaller curvature of the TB FS, illustrated in Figs 4e,f. This results in a slower deviation from perfect nesting of main and shadow bands along the FS, and hence an improved nesting, as shown by the increase of the grey area in Fig. 4f. In the limit of a vanishing curvature, the whole FS would be nested (we note that to maintain nesting for large changes of t⊥, also a slight change of qCDW is required). To capture the experimental change in gap dispersion a reduction of t⊥ as large as 25% in the TB model is necessary (see Supplementary Note 7 and Supplementary Fig. 6). This scenario also explains the symmetric gap reduction at the point of perfect nesting ((k, k)∼(0.17, 0.29) Å−1) where 2Δ is symmetric to EF. We conclude to observe a trend towards enhanced nesting, which was not observed as a function of temperature1718 and must be unique for non-equilibrium conditions.

The gap dynamics is linked to ion motion, as evidenced by the coherent oscillations of 2Δ. We speculate that such directed ion motion in the transient potential of the CDW distortion could transiently modify the orbital overlap in the network of Te orbitals and lead to the observed modification of the electronic dispersions and the improved FS nesting concluded above. As such, the enhanced nesting stabilizes the CDW in the optically excited state, however, in severe competition with fluctuations due to the excess energy, which destabilizes the CDW. The dynamics of the transient CDW state are therefore governed by the balance of these competing contributions. The limit in reduction of 2Δ≥0.17 eV observed in Fig. 3 originates from compensation between the increased fluctuations and the enhanced nesting, with their respective destabilizing and stabilizing effects on the CDW. While increasing fluctuations are the direct consequence of pumping at larger F, the stabilizing tendency is evidenced by the increasingly improved nesting for higher F (inset of Fig. 4d).

The possibility of modifying the transient electronic dispersion by optical excitation promises future control pathways for transient order in broken symmetry ground states of quantum materials as a complementary approach to nonlinear phononics1215. Optimization of the excitation conditions may lead to desired ordered states, which requires modifications of the competing contributions towards stabilization and destabilization of the ordered state. Such a controlled approach might well be possible by adjusting pump photon energy and intensity since a reduction of photon energy results in lower excitation probability of secondary processes and related fluctuations while intensity defines the instantaneous excited state potential.

Methods

Static high-resolution ARPES measurements

Single crystals of RTe3 were grown by slow cooling of a binary melt22 and cleaved in ultrahigh vacuum (base pressure <7 × 10−11 mbar) at T=10 K. Static FS maps of HoTe3 (Tc1=285 K, Tc2=120 K) at T=180 K and T=300 K have been obtained using a laser-based ARPES setup with hν=7 eV and using a hemispherical electron analyzer (Scienta SES2000) (ref. 30). Energy and momentum resolution were better than 5 meV and 0.005 Å−1, respectively.

Time-resolved ARPES measurements

For the trARPES experiments, the output of a commercial amplified Ti:sapphire laser system (Coherent RegA 9050) operating at 300 kHz repetition rate was used. Part of its output was frequency-quadrupled to yield hν=6.0 eV and used as probe pulses, while another part of the fundamental beam was time-delayed for excitation of the DyTe3 (Tc1=305 K, Tc2=50 K) sample, held at T=30 K during the measurements. The sample was mounted on a 45° slanted sample holder and oriented by Laue diffraction before trARPES measurements, which is necessary to reach the FS of RTe3 at a typical kinetic energy of 0.9 eV and large emission angles of >50°.

Position-sensitive time-of-flight spectrometer

A self-built position-sensitive time-of-flight spectrometer24 was used for photoelectron detection, mounted inside an ultrahigh vacuum chamber with a base pressure of <7 × 10−11 mbar. It enables simultaneous access to both in-plane momentum components k and k, along with the kinetic energy of the electrons, and thus allows efficient mapping of the electronic band structure. It consists of a field-free drift tube for photoelectrons, combined with a micro channel plate and two-dimensional delay line detector (RoentDek Hexanode Hex80). From the arrival time and impact position on the micro channel plate, the electron kinetic energy Ekin and both in-plane momentum components k, k of each single photoelectron are calculated and stored in a three-dimensional (E, k, k) grid24. The Hexanode design of the detector allows for efficient detection of multiple electrons per laser pulse, which is a prerequisite for obtaining high-quality trARPES data as presented in this work at typical count rates of ∼150 kHz.

The overall temporal, spectral and momentum resolutions of the setup were 100 fs, 50 meV and 2 × 10−3 Å−1, respectively. For further details of the experimental setup see Schmitt et al.623.

Additional information

How to cite this article: Rettig, L. et al. Persistent order due to transiently enhanced nesting in an electronically excited charge density wave. Nat. Commun. 7:10459 doi: 10.1038/ncomms10459 (2016).