Evidence for a Peierls phase-transition in a three-dimensional multiple charge-density waves solid.

Evidence for a Peierls phase-transition in a three-dimensional multiple charge-density waves solid.

Barbara Mansart Laboratory for Ultrafast Electrons, ICMP, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH Laboratory of Ultrafast Spectroscopy, ISIC, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH    Mathieu J.G. Cottet Laboratory for Ultrafast Electrons, ICMP, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH    Thomas J. Penfold Laboratory of Ultrafast Spectroscopy, ISIC, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH Laboratory of Computational Chemistry and Biochemistry, ISIC, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH SwissFEL, PSI, CH-5232 Villigen, CH    Stephen B. Dugdale H. H. Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol BS8 1TL, UK    Riccardo Tediosi Département de Physique de la Matière Condensée, Université de Genève, CH-1211 Genève 4, CH    Majed Chergui Laboratory of Ultrafast Spectroscopy, ISIC, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH    Fabrizio Carbone Laboratory for Ultrafast Electrons, ICMP, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH Laboratory for Ultrafast Electrons, ICMP, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH Laboratory of Ultrafast Spectroscopy, ISIC, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH [ Laboratory for Ultrafast Electrons, ICMP, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH Laboratory of Ultrafast Spectroscopy, ISIC, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH Laboratory of Computational Chemistry and Biochemistry, ISIC, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH SwissFEL, PSI, CH-5232 Villigen, CH H. H. Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol BS8 1TL, UK Département de Physique de la Matière Condensée, Université de Genève, CH-1211 Genève 4, CH Laboratory of Ultrafast Spectroscopy, ISIC, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH Laboratory for Ultrafast Electrons, ICMP, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, CH
Abstract

The effect of dimensionality on materials properties has become strikingly evident with the recent discovery of graphene. Charge ordering phenomena can be induced in one dimension by periodic distortions of a material’s crystal structure, termed Peierls ordering transition. Charge-density waves can also be induced in solids by strong Coulomb repulsion between carriers, and at the extreme limit, Wigner predicted that crystallization itself can be induced in an electrons gas in free space close to the absolute zero of temperature. Similar phenomena are observed also in higher dimensions, but the microscopic description of the corresponding phase transition is often controversial, and remains an open field of research for fundamental physics. Here, we photoinduce the melting of the charge ordering in a complex three-dimensional solid and monitor the consequent charge redistribution by probing the optical response over a broad spectral range with ultrashort laser pulses. Although the photoinduced electronic temperature far exceeds the critical value, the charge-density wave is preserved until the lattice is sufficiently distorted to induce the phase transition. Combining this result with ab initio electronic structure calculations, we identified the Peierls origin of multiple charge-density waves in a three-dimensional system for the first time.

Accepted for publication in Proc. Natl. Acad. Sci. USA, doi:10.1073/pnas.1117028109

Charge ordering phenomena occurring upon symmetry breaking are important in solids as they give rise to new current and spin flow patterns in promising materials such as organic conductors Lorenz2002 (), multilayered graphene Ho2010 () and transition metal oxides Qazilbash2011 (). The possibility to investigate the microscopic steps through which such ordering transition occurs gives also the opportunity to speculate on more general aspects of critical phenomena. Charge-density waves (CDWs) Middleton1992 (); Narayan1994 (), sandpile automata tadic () and Josephson arrays tinkham () have been investigated in relation to the scale invariance of self-organized critical phenomena Bak1987 (), of which avalanches are dramatic manifestations vajont (). In one dimension, Peierls demonstrated that at low temperature an instability can be induced by the coupling between carriers and a periodic lattice distortion. Such an instability triggers a charge ordering phenomenon and a metal-insulator phase transition, called Peierls transition, occurs peierls (). Like for Bardeen-Cooper-Schrieffer (BCS) superconductors, such an electron-phonon interaction driven transition is expected to be second order peierls (). While this situation is fairly established in mono-dimensional organic materials tcnq (), increased hybridization leading to higher dimensionality of a solid perturbs this scenario and makes the assessment of the microscopic origin of charge localization phenomena more difficult nbse3 (); tas2 (); mcmillan ().

Contrary to other low-dimensional CDW gruner () systems studied so far by time-resolved spectroscopies yusupov (); Eichberger2010 (); Schmitt2008 (); Tomeljak2009 (), LuIrSi presents a complex three-dimensional structure with several substructures such as one-dimensional Lu chains and three-dimensional cages in which a variety of many-body effects (including superconductivity) originate smaalen () (Fig. 6A). Although a CDW occurs below = 83 K, its microscopic origin is still debated since this transition is first order and isotropic without hysteresis Becker1999 (); Tediosi (), in contrast to the standard Peierls paradigm. Also, owing to their complexity, a detailed microscopic description of the properties of these solids is still lacking and their band structure has never been reported.

The charge redistribution induced by a phase transition in a solid can be obtained via the optical frequency-sum (f-sum) rule Tediosi (), which states that the integral over all frequencies of the optical conductivity (), termed spectral weight (), is constant and corresponds to the number of electrons per unit cell:

 SW=πe22meVNeff=∫∞0σ1(ω)dω (1)

where is the free electron mass, the unit cell volume and the total number of carriers. Performing the integral up to a frequency cut-off (partial sum-rule) gives the number of carriers contained in the electronic levels included in this energy range SI (). Optical weight redistributions are known to accompany every BCS-like phase transition, as described by the Ferrel-Glover-Tinkham sum-rule fgt (). Recently, it has been shown that information on the can be obtained in a model-independent fashion based on the analytical continuation of holomorphic functions devin ().

In this study we aim at obtaining the temporal evolution of the partial sum rule during the photoinduced annihilation of the CDW order parameter to identify the different states through which the electronic and the lattice structure evolve. To achieve this, the changes in the optical constants are monitored over a broad energy range with a temporal resolution better than the electron-phonon coupling time. We performed pump-probe reflectivity using a broad (1.5-3 eV) supercontinuum of 50 fs polarized pulses as a probe SI (). The temperature of the sample was controlled between 10 K and room temperature. The 1.55 eV pump fluence was varied between 0.8 and 3.1 mJ/cm, corresponding to an absorbed fluence between 0.4 and 1.55 mJ/cm; in these conditions and for an initial temperature of 10 K, the electronic temperature reaches values between 302 and 1100 K (see the transient temperature analysis below), corresponding to energies in the range 26-94 meV. Therefore, for the lowest pump fluences, the absorbed energy lies below the CDW gap estimated by static optical spectroscopy (80 meV) Tediosi (). Our ultrafast reflectivity data were combined with static optical spectroscopy Tediosi () and new ab initio band structure calculations, allowing the identification of the spectroscopic signature of the CDW melting, the description of the orbitals and atomic motions involved and their temporal evolution.

The static optical reflectivity of LuIrSi for light polarized along the one-dimensional chains (-axis), and perpendicular to them (-axis) is displayed in Fig. 6B. Fig. 6C shows the calculated reflectivity which agrees with the experiments; some discrepancies are observed at low energy along the -axis where correlation effects are expected due to the presence of the monodimensional chains. Figs 6(D-G) show the transient optical reflectivity for different temperatures and light polarizations. The low temperature dynamic reflectivity decays faster than at room temperature (200 fs compared with 300 fs, for the same pump fluence and energy of 1.9 eV). Its spectral behavior is also rather different: at 10 K, a large region is observed in which the photoinduced change in reflectivity switches sign (see white traces). The timescales and the absolute magnitude of the changes also depend on the polarization of the probe.

The evolution of is obtained from the transient reflectivity both via a Drude-Lorentz model (3D color maps in Fig. 8) and the model-independent approach (white lines with error bars in Fig. 8A and BSI (). In these plots, the integral in Eq. 1 is displayed as a function of the frequency cut-off, emphasizing the spectral regions of positive and negative change of the effective number of carriers. When the sample temperature is below and the probe light is polarized along the -axis (Fig. 8B), the change, defined as ((t,)-(t)), becomes negative after 200 fs. Instead, it is always positive when the polarization is perpendicular to the -axis (Fig. 8A), or at room temperature (Fig. 8C). This sign change at low temperature along the -axis signals a transfer from energies higher than 1.5 eV, typical of strongly correlated solids Rozenberg1996 (). In Fig. 8D and E, the quantities ((,)-(, ))/(,) are displayed for = 60 K and = 120 K, with light polarized along the -axis (red line) or the -axis (black line), and = 120 K and = 290 K for light polarized along the -axis (light blue line) (data taken from Tediosi ()). The transfer across the CDW temperature along the -axis becomes negative above 1.5 eV, as observed in the transient data; therefore, the spectroscopic signature of the CDW melting is the redistribution of carriers above 1.5 eV where bands composed of Lu and Ir states are found. Remarkably, we observe that this feature appears several tens of fs after laser excitation, suggesting that before the CDW melts, an unconventional state of matter is produced in which some charge order survives despite the large energy deposited in the electronic subsystem.

To verify that indeed the behavior of the optical constants in the 2-3 eV region is mainly sensitive to the CDW melting, we measured the transient reflectivity spectra at different fluences (0.8 to 3.1 mJ/cm) and temperatures (10 K, 80 K and 300 K). In Fig. 7B and C, the transient reflectivity at 1.77 and 2.21 eV are shown for the different pump fluences. At 1.77 eV, the reflectivity changes are always positive, while at 2.21 eV a negative dip develops at fluences higher than 2.1 mJ/cm. In Fig. 7D, the temperature dependence of the transient reflectivity is also shown, at a probe energy of 2.21 eV and for a fluence of 3.1 mJ/cm. Also in this case, the change of reflectivity has a negative part, associated to the CDW melting, at 10 K, while it is always positive at temperatures higher than .

The time scales involved in the melting of the CDW are obtained through a transient temperature analysis using the three-temperature model Kaganov1957 (); allen (), appropriate for anisotropic materials with selective electron-phonon coupling perfetti (); Mansart2010 (); Carbone2008 (). This model describes the energy transfer from excited carriers to the lattice as an effective heat exchange between them. We obtained the temporal evolution of electronic and lattice temperatures SI () and display them in Fig. 7A together with the transient at 10 K and polarization along the -axis.

The CDW melting, indicated by becoming negative at high energy, occurs after the lattice reaches a temperature above . While the electronic temperature crosses and reaches values as high as 1100 K, the system exists in a transient state where the signature of the molten CDW state is still absent. These simulations provide an estimate of the Eliashberg electron-phonon coupling parameter 0.6 and a fraction of strongly coupled modes 20, in agreement with the value of = 0.66 determined by combined heat capacity and magnetic susceptibility measurements under pressure Shelton1986 (). The rise time of our signal is in our temporal resolution range, i.e. 40-50 fs. This, together with the three-temperature model simulations, provides an estimate of the excited electron thermalization time around 50 fs, much faster than the CDW melting time. This rules out any possible electronic melting as in this case a clear signature of the molten CDW state should be present as early as the electron thermalization time, like recently observed in TaS cavalleri2011 ().

To speculate on the microscopic nature of the intermediate state, we consider that during the first 200 fs the electronic structure is relaxing after being put abruptly out-of-equilibrium by the pump pulse. By comparing the spectra of the photoinduced reflectivity changes at different time delays to the difference spectra between the static reflectivity recorded at different temperatures, the impact of the out-of-equilibrium electronic structure can be estimated. In Fig. 4A, the transient reflectivity spectra at 60, 120, 200 and 500 fs after excitation are displayed together with R(150 K) - R(25 K). When the electronic structure and the lattice are at equilibrium, the temperature difference spectrum exhibits a positive peak around 2 eV and a minimum around 2.7 eV across the phase transition (red line). This behavior is quite similar to that of the transient reflectivity after 500 fs, at which time the electronic structure is likely at equilibrium with the lattice. shows also a positive peak around 2 eV and a negative minimun around 2.7 eV; this similarity indicates that photoexcitation and temperature induce the same carriers redistribution in the solid. At early times, while the shape of the transient spectra is still similar to the static one showing a positive peak around 2 eV and a minimum around 2.7 eV, a quantitative discrepancy is more pronounced and could originate from the very large electronic temperature jump induced by laser excitation ( 1000 K) compared to the static case in which the system is heated by a hundred degrees. These spectra suggest that, consistently with the three-temperature model analysis, a thermal electronic distribution of carriers is established within few tens of fs.

It is important to note at this point that while a negative reflectivity change can be associated to the phase transition, the complete melting of the CDW state can only be claimed when the optical SW, which is the cumulative effect of the positive and negative changes, becomes negative. At a fluence of 3.1 mJ/cm, around 1 carrier per unit cell is photoexcited. Photoexcited carriers create ”hot areas” in which the CDW can be perturbed. This electronic perturbation diffuses in space via its interaction with the lattice and after some time, a thermally CDW molten state is reached. The delay of this process is governed by the microscopic details of the coupling between carriers and phonons in the CDW state. This scenario is pictorially represented in the cartoon in Fig. 4B-D.

The presented situation is very distinct from an electronically driven melting of the charge order, since in that case the negative optical should be observable without any delay; therefore, even though an instantaneous ( 60 fs) CDW melting would locally happen upon light excitation, it should propagate and melt the whole charge order within the electronic thermalization time, which we found to be around 60 fs. Such a situation is observed for example in TaS cavalleri2011 (), where a non-Peierls phase transition is claimed to take place.

Further microscopic details are obtained through ab initio electronic structure calculations SI (), performed via both an all-electron full-potential linearized augmented plane-wave method elk () and via pseudopotentials qe (). Eleven bands having mixed character (color coded in Fig. 9B) cross the Fermi level. The density of states (DOS) (Fig. 9B) shows that the Fermi surface has mostly Ir character. Some Lu character is found, coming from the 5 orbitals strongly hybridized with both Si and Ir, producing delocalized electrons as observed in transition metal-silicides CarboneXAS (). The Fermi surface is presented in Fig. 9A, and from calculations of the generalized susceptibility, we identified two bands having nesting vectors of 2/7 and 3/7 of , as observed by X-ray diffraction below  Becker1999 (). This suggests the presence of multiple CDWs occurring at the same transition temperature. These two nested bands account for approximately 30 of the Fermi energy DOS, in agreement with estimates based on magnetic susceptibility measurements Shelton1986 (). While one of them comes mostly from Ir and Lu1 (Lu1 atoms being along the chains Becker1999 (), Fig. 9C), the second one is mainly formed by Lu2 - and all the orbitals of both Lu3 and Ir. This complex and three-dimensional nature of the orbitals reveals the possibility of a 3D Peierls transition. We also show in Fig. 9C the motion of two stronglycoupled optical phonons, yellow and orange arrows respectively, having frequencies (1.4 THz and 2.8 THz) close to the coherent oscillations observed in the same material Tomeljak2006 (). Ab initio calculations also estimate the electron-phonon coupling parameter of each individual lattice mode, yielding a total 0.6 for the 18 of the most coupled modes, in agreement with the three-temperature model simulations.¨

The orbital occupation-number change induced by the CDW melting is signalled by the optical transfer around a pivot energy close to 1 eV, both in static and time-resolved data. Upon excitation, acting as a photodoping process, carriers are transferred from the Lu and Ir states forming the CDW SI (). This transfer depopulates the nested bands and redistributes the excited carriers above the Fermi level, increasing the low-energy as well as the effective electronic temperature. The observation of a delay between photoexcitation and the melting of the CDW strongly suggests the Peierls nature of the phase transition in these compounds, rather than a pure electronic effects as in different systems vo2 (). In fact, while CDW are known to couple to specific lattice distortions, the intimate microscopic origin of the mechanism described in peierls () is that a periodic modulation of the lattice induces the Fermi surface instability causing the consequent charge ordering. In other words, to understand correctly the origin of such a phase transition, a notion of causality must be invoked, which is directly obtainable via time-resolved experiments. Indeed, one must wait until the atoms have moved sufficiently to suppress the distortion promoting the CDW state Schmitt2008 (); Cavalleri2004 (); the latency time can be of the order of one quarter of the corresponding phonon period. In our case, 200 fs gives a frequency of 1.4 THz, corresponding to the mode represented in orange arrows in Fig. 9C. While the direct involvement of the lattice in the CDW formation reveals the Peierls nature of the metal-insulator transition in this material, its weakly one-dimensional electronic structure, together with the signatures of a first order CDW transition and the suggested presence of multiple CDWs, are at odds with the conventional scenario describing Peierls transitions gruner (). Recently, the notion that Peierls transitions can happen in higher dimensions has been challenged by new ultrafast experiments in TaS cavalleri2011 (), which has a 2D charge ordering transition. Also in this case, the delay between excitation and melting has been used as the key-argument for determining the microscopic nature of the phase transition. Our experiments contribute to this exciting and current debate by showing evidence for a Peierls transition in a three-dimensional system.

The proper description of phase transitions is important for fundamental physics; in fact, phase transitions in complex systems often evolve through several unconventional intermediate states baum (). Critical phenomena involving abrupt perturbation of ordered systems are observed in giant land-slides vajont () and sandpile automata tadic (). In these cases and in CDWs, an external perturbation may induce a transient state in which the order is preserved while undergoing a collective reorganization Middleton1992 (); Narayan1994 (). The speed of rock-slides, CDW slides, and flux jumps in the vortex state of type II superconductors, are self-organized critical phenomena and obey similar critical laws Middleton1992 (); tinkham (); Bak1987 (); richard (). In our experiments, the disappearance of the charge order manifests itself after a significant delay from the abrupt perturbation, and in the intermediate state the order seems preserved despite the large energy deposited.

I Acknowledgements

The authors acknowledge useful discussions with D. van der Marel, J. Demsar, C. Giannetti and A. Pasquarello. The LuIrSi single crystal have been provided by J. A. Mydosh. We thank U. Roethlisberger (EPFL, Lausanne) and ACRC (Bristol) for computer time. This work was supported by the ERC grant N 258697 USED and by the Swiss NSF via contract .

References

• (1) Lorenz T, et al. (2002) Evidence for spin–charge separation in quasi-one-dimensional organic conductors. Nature 418:614-617.
• (2) Ho CH, et al. (2010) Magneto-electronic properties of rhombohedral trilayer graphene: Peierls tight-binding model. Annals of Physics 326:721-739.
• (3) Qazilbash MM, et al. (2011) Nanoscale imaging of the electronic and structural transitions in vanadium dioxide. Phys. Rev. B 83:165108.
• (4) Middleton AA, Biham O, Littlewood PB, Sibani P (1992) Complete mode locking in models of charge-density waves. Phys. Rev. Lett. 68:1586-1589.
• (5) Narayan O, Middleton AA (1994) Avalanches and the renormalization group for pinned charge-density waves. Phys. Rev. B 49:244-256.
• (6) Tadić B, Nowak U, Usadel KD, Ramaswamy R, Padlewski S (1992) Scaling behavior in disordered sandpile automata. Phys. Rev. A 45:8536-8545.
• (7) Benz SP, Rzchowski MS, Tinkham M, Lobb CJ (1990) Fractional giant Shapiro steps and spatially correlated phase motion in 2D Josephson arrays. Phys. Rev. Lett. 64:693-696.
• (8) Bak P, Tang C, Wiesenfeld K (1987) Self-organized criticality: An explanation of the 1/ noise. Phys. Rev. Lett. 59:381-384.
• (9) Kilburn CRJ, Petley DN (2003) Forecasting giant, catastrophic slope collapse: lessons from Vajont, Northern Italy. Geomorphology 54:21-32.
• (10) Peierls RE (1955) Quantum Theory of Solids Oxford University Press, London p. 108.
• (11) Coleman LB et al. (1973) Superconducting fluctuations and the peierls instability in an organic solid. Solid State Comm. 12:1125-1132.
• (12) Ong NP, Monceau P (1977) Anomalous transport properties of a linear-chain metal: NbSe. Phys. Rev. B 16:3443-3455.
• (13) Moncton DE, Axe JD, DiSalvo FJ (1975) Study of Superlattice Formation in 2-NbSe and 2-TaSe by Neutron Scattering. Phys. Rev. Lett. 34:734-737.
• (14) McMillan WL (1977) Microscopic model of charge-density waves in 2-TaSe Phys. Rev. B 16:643-650.
• (15) Grüner G (1988) The dynamics of charge-density waves. Rev. Mod. Phys. 60:1129-1181.
• (16) Yusupov R et al. (2010) Coherent dynamics of macroscopic electronic order through a symmetry breaking transition. Nature Physics 6:681-684.
• (17) Eichberger M et al. (2010) Snapshots of cooperative atomic motions in the optical suppression of charge density waves. Nature 468:799-802.
• (18) Schmitt F et al. (2008) Transient Electronic Structure and Melting of a Charge Density Wave in TbTe. Science 321:1649-1652.
• (19) Tomeljak A et al. (2009) Dynamics of Photoinduced Charge-Density-Wave to Metal Phase Transition in KMoO. Phys. Rev. Lett. 102:066404.
• (20) van Smaalen S et al. (2004) Multiple charge-density waves in RIrSi (R=Ho, Er, Tm, and Lu). Phys. Rev. B 69:014103.
• (21) Becker B et al. (1999) Strongly coupled charge-density wave transition in single-crystal LuIrSi. Phys. Rev. B 59:7266-7269.
• (22) Tediosi R et al. (2009) Evidence for strongly coupled charge-density-wave ordering in three-dimensional IrSi compounds from optical measurements. Phys. Rev. B 80:035107.
• (23) see Supporting Information.
• (24) Tinkham M, Ferrell RA (1959) Determination of the Superconducting Skin Depth from the Energy Gap and Sum Rule. Phys. Rev. Lett. 2:331-333.
• (25) Kuzmenko AB, van der Marel D, Carbone F, Marsiglio F (2007) Model-independent sum rule analysis based on limited-range spectral data. New J. Phys. 9:229.
• (26) Rozenberg MJ, Kotliar G, Kajueter H (1996) Transfer of spectral weight in spectroscopies of correlated electron systems. Phys. Rev. B 54:8452-8468.
• (27) Kaganov MI, Lifshitz IM, Tanatarov LV (1957) Relaxation between electrons and the crystalline lattice. Sov. Phys. JETP 4:173.
• (28) Allen PB (1987) Theory of thermal relaxation of electrons in metals. Phys. Rev. Lett. 59:1460-1463.
• (29) Perfetti L et al. (2007) Ultrafast Electron Relaxation in Superconducting BiSrCaCuO by Time-Resolved Photoelectron Spectroscopy. Phys. Rev. Lett. 99:197001.
• (30) Mansart B et al. (2010) Ultrafast transient response and electron-phonon coupling in the iron-pnictide superconductor Ba(FeCo)As. Phys. Rev. B 82:024513.
• (31) Carbone F, Yang DS, Giannini E, Zewail AH (2008) Direct role of structural dynamics in electron-lattice coupling of superconducting cuprates. Proc Natl Acad Sci USA 105:20161-20166.
• (32) Shelton RN et al. (1986) Electronic phase transition and partially gapped Fermi surface in superconducting LuIrSi. Phys. Rev. B 34:4590-4594.
• (33) Petersen JC et al. (2011) Clocking the Melting Transition of Charge and Lattice Order in 1-TaS with Ultrafast Extreme-Ultraviolet Angle-Resolved Photoemission Spectroscopy. Phys. Rev. Lett. 107:177402.
• (34)
• (35) Giannozzi P et al. (2009) QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys.: Condens. Matter 21:395502.
• (36) Carbone F et al. (2006) Electronic structure of MnSi: The role of electron-electron interactions. Phys. Rev. B 73:085114.
• (37) Tomeljak A et al. (2006) Femtosecond real-time studies of LuIrSi charge density wave compound. Workshop ”recent developments in low dimensional charge density wave conductors”, Croatia.
• (38) Wentzcovitch RM, Schulz WW, Allen PB (1994) VO: Peierls or Mott-Hubbard? A view from band theory. Phys. Rev. Lett. 72:3389-3392.
• (39) Cavalleri A, Dekorsy Th, Chong HHW, Kieffer JC, Schoenlein RW (2004) Evidence for a structurally-driven insulator-to-metal transition in VO : A view from the ultrafast timescale. Phys. Rev. B 70:161102(R).
• (40) Baum P, Yang D-S, Zewail AH (2007) 4D Visualization of Transitional Structures in Phase Transformations by Electron Diffraction. Science 318:788-792.
• (41) Richard J, Monceau P, Renard M (1982) Charge-density-wave motion in NbSe. II. Dynamical properties. Phys. Rev. B 25:948-970.

Supporting Information: Evidence for a Peierls phase-transition in a three-dimensional multiple charge-density waves solid.

B. Mansart ]barbara.mansart@epfl.ch M. Cottet T. J. Penfold S. B. Dugdale R. Tediosi M. Chergui F. Carbone

Appendix A Experimental setup

We performed an experimental study of the transient optical properties of LuIrSi using pump-probe broadband reflectivity. Our experimental setup, located at the Laboratory of Ultrafast Spectroscopy, EPFL, Switzerland, is presented in Fig. S1. An amplified Ti:sapphire laser system provides 1.55 eV pulses of 40 fs duration at a repetition rate of 1kHz. This laser beam is separated into two paths by a beamsplitter. The pump beam, whose energy is fixed at 1.55 eV, and the probe beam which is converted into a supercontinuum by a CaF nonlinear crystal. The probe spectrum covers a range between 1.5 and 3.2 eV. The fluence of the white light probe is several orders of magnitude lower than the pump beam and has no effect on the measurements. The pump pulse is used to excite the system at the time , and the consequent variations of reflectivity are recorded with a delay by changing the optical path difference between pump and probe (which can be achieved by setting a variable delay in the probe arm).

Small reflectivity changes (down to the order of 10) are detectable owing to the signal-to-noise ratio of 10, obtained by the use of an optical chopper on the pump beam which reduces its repetition rate by a factor of two. This allows the recording, for each time delay, of a single-shot spectrum with and without the pumping pulse. Furthermore, the variations of the white light intensity are taken into account using a reference signal, measured just before the sample position.

Two identical spectrometers are used to detect the transient reflectivity signal and the reference one. By dispersing the incident light with optical gratings and sending it to an array detector, we obtain the complete spectrum for each single pulse with an energy resolution better than 1 meV. Subsequently, the group velocity dispersion introduced by the chromatic aberration of the optics is corrected before the transient reflectivity analysis.

The pump-beam diameter on the sample surface was 250 m, while that of the probe was set to 100 m in order to probe a uniformly excited area. Both spot sizes and spatial overlap were rigorously verified using a beam profiler and a charge-coupled device camera. During this experiment, the pump fluence was varied between 0.8 and 3.1 mJ/cm.

Prior to these experiments, the sample was oriented using X-ray diffraction and polished to obtain good optical quality. Its surface contained both (001) and (100) directions. Therefore, pump and probe polarizations were changed independently with quarter-waveplates and were oriented either parallel or orthogonal to the -axis (containing the Lu chains). Finally, in order to cool down to 10 K and precisely monitor the sample temperature, the latter was placed in a liquid He cryostat.

Such kind of experimental setup has been successfully employed to investigate photoexcited states in strongly correlated systems, such as the charge-transfer insulator CuGeO [S1], as well as in bulk graphite and graphene paper [S2].

Appendix B Spectral Weight analysis

The definition of the optical spectral weight () in a solid, given in Eq. (1) of the main article, makes clear that the evaluation of requires the knowledge of the real part of the optical conductivity down to zero frequency. Since this is not experimentally feasible, common procedures (including extrapolations or modeling based on the Drude-Lorentz formalism) are invoked to obtain the whole spectrum. Recently, it has been shown mathematically that according to the theorems on analytical continuation of holomorphic functions, and taking advantage of the Kramers-Kronig relation between the real and the imaginary part of the complex optical conductivity (), the partial integral can be re-written as a function of a limited spectral interval in the form of Carleman-Goluzin-Krylov equations [S3]:

 SW(Ωc)=Re {limn→∞∫ωmaxωminun(ω′,Ωc)σ(ω′)dω′} (2)

where are properly chosen kernel functions. In the case of real experimental data (with finite energy resolution and statistical noise), Eq. (2) cannot be strictly applied, and a numerical procedure based on the idea of the analytical continuation has to be invoked in order to obtain the integral in a model-independent fashion [S4]. A code called “Devin” has been developed which allows the integral to be extracted from the complex dielectric function in a limited data range in a model-independent fashion [S4].

We performed a Drude-Lorentz analysis of our experimental data, and verified the results via the model-independent routine “Devin”. The approach is described below in more detail.

To obtain the complex dielectric function, we fitted our transient reflectivity data to a highly flexible Drude-Lorentz function [S5]; the underlying Drude-Lorentz model of the optical data is based on the knowledge of the static optical spectra between a few meV and 4.6 eV, data taken from Ref. [S6].

For the fitting of the transient data, the Lorentz oscillators in our experimental range are allowed to change in order to reproduce the dynamical reflectivity. This procedure yields the transient and , displayed together with the static spectra in Fig. S2. In our experimental frequency range, the real and imaginary part of the dielectric constant have rather similar absolute values (Fig. S2A and E). For this reason, the optical reflectivity, defined as , is equally sensitive to the reactive and the absorptive components of the dielectric function. As a consequence, the complex obtained through the Drude-Lorentz fitting of the optical reflectivity matches very well the one obtained via spectroscopic ellipsometry (see Fig. S2A and E and Ref. [S6]). The transient dielectric function is displayed in Fig. S2(B-D, F-H) for light polarized along the -axis and the -axis at 10 K, and for light polarized along the -axis at room temperature.

Such a Drude-Lorentz modeling of the data also yields the spectral weight; the overall static spectrum for both light polarizations is displayed in Fig. S3. In the same plot, we show the derivative of as a function of the number of carriers which exhibits an inflection point between 1 and 1.5 eV for light polarized along the -axis; this inflection point is found to be between 10 and 15 carriers. Since Lu has a electronic configuration, and we found a strong hybridization with at least one or two Ir bands, we expect between 10 and 20 electrons to form the material conduction band and therefore the low energy -axis conductivity should be mainly sensitive to these orbitals.

By integrating the transient conductivity, we calculated at each time step and obtained the time evolution of via Drude-Lorentz model, presented as 3D color maps in Fig. 2(A-C) of the main article.

In order to verify that the transient changes of the we obtained by extrapolating the Drude-Lorentz model down to low frequency are robust and do not depend on the latter extrapolation, we fed into the “Devin” routine the real and imaginary parts of the dielectric function as an input in order to generate the up to a given frequency cut-off. Through this procedure it is also possible to obtain an estimate of the uncertainty related to the knowledge of the spectrum being in a limited range with a limited resolution and with a given noise level. Because we are interested in the relative changes of the induced by pump pulses, we use as an error our estimated accuracy in the determination of the transient dielectric function, which is of the order of , while the resolution of our spectrometer is approximately 1 meV. The transient dielectric function is fed into “Devin” and the model-independent changes are shown in Fig. 2(A-B) of the main article as white lines with errorbars. The fact that both these analyses give consistent results validate our description of the transient optical reflectivity.

The reason why we can successfully obtain information on the CDW-induced spectral weight transfer from a limited spectrum in the visible light region is that, in this sample, the optical anisotropy, related to the Lu one-dimensional chains, strongly manifests itself around 2 eV. This is visible in the static reflectivity spectrum (Fig. 1 of the main article), and also in the temperature-dependent dielectric constants obtained via ellipsometry and reported in Ref. [S6], which in fact show a particularly large CDW-induced change between 1.5 and 2.5 eV.

Also, in strongly correlated solids, a phase transition can induce spectral weight transfers across a very broad energy range, extending to the visible light region. In these cases, it has been shown that the spectral region in the proximity of the screened plasma frequency of the material is very sensitive to these shifts [S7, S8]. Here, the screened plasma frequency along the -axis (direction of the charge-density-wave) lies within our experimental range (zero crossing of in Fig. S3A).

Appendix C Three-temperature model

The two-temperature model was first introduced by Kaganov et al. [S9] in order to describe the energy transfer between the electron and lattice subsystems, in the case where the electrons are driven to a higher temperature than the phonons. Its validity relies on the electronic temperature being of the order of the Debye temperature (366 K in LuIrSi [S10]) and on the isotropy of the electron-phonon (e-ph) coupling function. For anisotropic materials like cuprates [S11] or iron-pnictides [S12], a selective coupling between electrons and a subset of the total phonon modes may be taken into account using the three-temperature model (3TM), governed by the following equations:

 2Ce∂Te∂t=2(1−R)lsI(t)+∂∂z(κe∂Te∂z)−g(Te−Th)αCL∂Th∂t=g(Te−Th)−gc(Th−Tc)(1−α)CL∂Tc∂t=gc(Th−Tc) (3)

where , and are the temperatures of electrons, efficiently coupled phonons (“hot phonons”) and remaining modes (“cold phonons”), respectively. = is the electronic specific heat, is that of the lattice. is the fraction of efficiently coupled modes, is the static reflectivity and is the penetration depth (both at the pump energy). the pumping intensity, and is the electronic thermal conductivity (which is ignored in our case). Indeed, the electronic diffusion term changes neither the relaxation rates nor the e-ph thermalization time. The constant governs the energy transfer rate from electrons to hot phonons, and is related to the second moment of the Eliashberg function through  [S13], being the dimensionless e-ph coupling constant. governs the energy relaxation from hot phonons to the remaining lattice.

To perform the 3TM simulations, we used an iteration procedure, calculating at each time step the depth-dependent temperature profiles. At each depth and time step, we iterate the electronic and lattice part of the specific heat. The electronic part is taken as , with =23.42  [S10], and the lattice part is taken from Ref. [S14].

The e-ph coupling constant was obtained via these simulations, given a logarithmic average of the phonon energies of 30 meV, obtained by static optical spectroscopy in Ref. [S6]. The e-ph coupling constant should in any case depend on the probing energy; as a check, we simulated the experimental data for several probe wavelengths in order to show the result constancy in the experimental energy range. The constant obtained through these simulations is given in Fig. S4C, from which we may notice its energy-independence in the range 1.6 –2 eV.

In fact, as this energy range is far from the electronic inter-band transitions (except very close to 2 eV, see the Lorentz oscillator positions in Fig. 1B of the article), the probe is mainly sensitive to intra-band transitions.

Considering , both experimental and theoretical determinations show that the carriers involved in the optical conductivity at 1.5 eV are mainly the Lu 5 and some Ir 5, which mainly belong to the conduction band of the material; in fact, the plasma edge of the solid is close to 1.7 eV along the -axis. Therefore, even though the numerous possible inter-band transitions should give a non-negligible contribution to the optical absorption, screening effects prevent the existence of long lived excited states at those energies (1.5 eV) in such a metallic system, and the typical electron-hole lifetime should not exceed 10-20 fs. Therefore, on the time-scales of our experiments, 50 fs, most if not all of the light excitation creates electron-hole pairs which relax instantaneously back in the conduction band with an excess energy that can be described by an effective temperature. Consequently, applying the 3TM in this energy range allows the measurement of the coupling between electrons at the Fermi level and lattice.

The linear dependence of the fast relaxation time as a function of the maximum electronic temperature reached (see Fig. S4A), as expected considering this relaxation as due to e-ph coupling [S13], reflects also the applicability of the 3TM in our case.

We may notice the excellent agreement between the parameters obtained from this effective temperature model and the calculated partial e-ph coupling constants presented in the following section. Indeed, both the fraction of efficiently coupled phonons (Fig. S4D) and the sum of their partial e-ph coupling (Fig. S4C) may be retrieved by 3TM simulations, enhancing our confidence in the reliability of this model in describing our transient data.

Appendix D Electronic structure calculations

Two separate sets of calculations were performed, the first using the highly accurate all-electron full-potential linearized augmented plane-wave method (often referred to as FP-LAPW), as implemented in the ELK code [S15], and the second using pseudopotentials, as implemented within the QUANTUM ESPRESSO code [S16], for the calculation of the phonons and e-ph coupling.

The lattice constants used were the experimental ones (Å, Å), as were the atomic positions [S17], giving a unit cell of two formula units (38 atoms). The calculations made using the ELK code [S15] had a cutoff (for planewaves in the interstitial region) determined by , where is the smallest muffin-tin radius. These muffin-tin radii were Lu (2.3328 a.u.), Ir (2.3654 a.u.) and Si (1.8369 a.u.). The self-consistent cycle was carried out on a mesh of 64 -points in the irreducible (one sixteenth) of the Brillouin zone, and the Fermi surface evaluated from the converged potential on a denser mesh of 1936 -points within one eighth of the Brillouin zone.

The calculations presented are scalar-relativistic and omit spin-orbit coupling (after the effect of its inclusion was investigated and found to be negligibly small at the Fermi level.) The exchange-correlation functional was that of Perdew-Wang/Ceperley-Alder, Ref. [S18].

The calculated band structure in the energy range [-8 eV, +2 eV] is given in Fig. S5A. In this range, the latter is composed of strongly hybridized Lu and , Ir and Si , the Lu giving localized density of states around -4.5 eV, while electrons at the Fermi surface come mainly from Lu , Ir and Si as shown in Fig. S5B and Fig. 5 of the main article. A large amount of “interstitial” electronic density is also found at the Fermi level, indicating their delocalized nature as well as the strong hybridization degree in this system.

The e-ph coupling resulting from the mass enhancement may be calculated from the density of states (DOS) at the Fermi level, following the formula , where is the renormalized DOS and the bare one, both at the Fermi level ; is the free electron mass and the effective one. The renormalized DOS obtained by specific heat measurement being 0.26/eV/atom/spin, (i.e. 19.76/eV), and the bare DOS calculated here of 11.02/eV if one considers that 30 of the Fermi surface is gapped, one obtains giving an e-ph coupling constant .

We also present in Fig. S6 the real and imaginary parts of the calculated susceptibility for two bands, corresponding of the Fermi surfaces indicated by arrows in Fig. 5A of the main article. One may notice the peak at two different wavevector, in both real and imaginary parts. These peaks appear at nesting wavevectors which are different between the two presented bands, indicating the formation of several charge-density waves.

For the phonon and e-ph coupling calculations, linear response within the pseudopotential approach provided by the QUANTUM ESPRESSO code was used [S16]. The Lu pseudopotential had the electrons as core states, and the resulting electronic structure was verified against that produced by the all-electron ELK calculation. Both ultrasoft and norm-conserving pseudopotentials were used, the latter needed for the calculation of the dielectric function with the EPSILON package within QUANTUM ESPRESSO. The cut-offs for the energy and charge density were 40 Ry and 400 Ry, respectively. Phonons at the point were calculated from self-consistent calculations performed on a (6,6,16) -point mesh and e-ph parameters were evaluated on a denser mesh of (12,12,32). For the calculations of dielectric function, the routine EPSILON within the QUANTUM ESPRESSO package was used. Here, Brillouin zone integrations were carried out on a mesh of over 6000 -points with a Gaussian smearing of 0.25 eV applied to the evaluation of the inter-band contribution.

We present the partial e-ph coupling constants for each of the 114 modes of the system as a function of the mode energy at in Fig. S7B. Two strongly coupled modes are found having respectively a partial e-ph coupling constant of 0.166 and 0.131, involving mainly Si atom motion since the latter are the lighest elements of the crystal structure. One of these two modes is represented with yellow arrows in Fig. 5C of the main paper. The sum over all modes gives , so in good agreement with the value obtained from the DOS renormalization at the Fermi level (when one considers that the electron-phonon calculation has only been made at the point).

The histogram of partial e-ph coupling constants is shown in Fig. S7A. From such a representation, one is able to extract the partial e-ph coupling of some percentage of the most efficiently coupled modes, as represented with a red-dashed line in Fig. S7A. This area contains 21 modes, (or 18 of the phonons), and summing up their constant one obtain 0.62. It indicates that the 18 of the phonon modes having the largest e-ph coupling have, all together, an e-ph coupling constant of 0.62. This observation is in excellent agreement with the results of 3TM simulations, also indicating the reliability of these electronic structure and phonon calculations. This value of gives a superconducting transition temperature of 9 K using the McMillan formula [S19], as experimentally observed in these system above the critical pressure, and suggests that conventional e-ph mediated superconductivity takes place in LuIrSi, obeying Bardeen-Cooper-Shrieffer theory. Indeed, since upon photoexcitation the CDW is molten, and therefore does not compete anymore with superconductivity (as in the system under pressure), obtaining a similar value may reflect effectively the order parameter giving rise to superconductivity and not the one related to Peierls ordering.

References

• (42)
• (43) Giannetti C et al. (2009) Disentangling thermal and nonthermal excited states in a charge-transfer insulator by time- and frequency-resolved pump-probe spectroscopy. Phys. Rev. B 80:235129.
• (44) Carbone F et al. (2011) Femtosecond carrier dynamics in bulk graphite and graphene paper. Chem. Phys. Lett. 504:37-40.
• (45) Aizenberg LA (1993) Carleman’s Formulas in Complex Analysis. (Dordrecht: Kluwer).
• (46) Kuzmenko AB, van der Marel D, Carbone F, Marsiglio F (2007) Model-independent sum rule analysis based on limited-range spectral data. New J. Phys. 9:229.
• (47) Kuzmenko AB (2005) Kramers-Kronig constrained variational analysis of optical spectra. Rev. Sci. Instrum. 76:083108.
• (48) Tediosi R et al. (2009) Evidence for strongly coupled charge-density-wave ordering in three-dimensional IrSi compounds from optical measurements. Phys. Rev. B 80:035107.
• (49) Carbone F et al. (2006) Doping dependence of the redistribution of optical spectral weight in BiSrCaCuO. Phys. Rev. B 74:064510.
• (50) Molegraaf HJA, Presura C, van der Marel D, Kes PH, Li M (2002) Superconductivity-Induced Transfer of In-Plane Spectral Weight in BiSrCaCuO. Science 295:2239-2241.
• (51) Kaganov MI, Lifshitz IM, Tanatarov LV (1957) Relaxation between electrons and the crystalline lattice. Sov. Phys. JETP 4:173.
• (52) Shelton RN et al. (1986) Electronic phase transition and partially gapped Fermi surface in superconducting LuIrSi. Phys. Rev. B 34:4590-4594.
• (53) Perfetti L et al. (2007) Ultrafast Electron Relaxation in Superconducting BiSrCaCuO by Time-Resolved Photoelectron Spectroscopy. Phys. Rev. Lett. 99:197001.
• (54) Mansart B et al. (2010) Ultrafast transient response and electron-phonon coupling in the iron-pnictide superconductor Ba(FeCo)As. Phys. Rev. B 82:024513.
• (55) Allen PB (1987) Theory of thermal relaxation of electrons in metals. Phys. Rev. Lett. 59:1460-1463.
• (56) Kuo Y-K et al. (2003) Ionic size and atomic disorder effects on the charge-density-wave transitions in RIrSi (R=Dy-Lu). Phys. Rev. B 67:195101.
• (57)
• (58) Giannozzi P et al. (2009) QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys.: Condens. Matter 21:395502.
• (59) Becker B et al. (1999) Strongly coupled charge-density wave transition in single-crystal LuIrSi. Phys. Rev. B 59:7266-7269.
• (60) Perdew JP, Wang Y (1992) Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B 45:13244-13249.
• (61) McMillan WL (1968) Transition Temperature of Strong-Coupled Superconductors. Phys.Rev. 167:331-344.
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters