Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 12, Article number: 6621 (2022 ) Cite this article Free-Space Laser
The advent of free-electron lasers opens new routes for experimental high-pressure physics, which allows studying dynamics of condensed matter with femtosecond resolution. A rapid compression, that can be caused by laser-induced shock impact, leads to the cascade of high-pressure phase transitions. Despite many decades of study, a complete understanding of the lattice response to such a compression remains elusive. Moreover, in the dynamical case (in contrast to quasi-static loading) the thresholds of phase transitions can change significantly. Using the third harmonic pump–probe technique combined with molecular dynamics to simulate the terahertz (THz) spectrum, we revealed the dynamics of ultrafast laser-induced phase transitions in MgF2 in all-optical experiment. Tight focusing of femtosecond laser pulse into the transparent medium leads to the generation of sub-TPa shock waves and THz coherent phonons. The laser-induced shock wave propagation drastically displaces atoms in the lattice, which leads to phase transitions. We registered a cascade of ultrafast laser-induced phase transitions (P42/mnm ⇒ Pa-3 ⇒ Pnam) in magnesium fluoride as a change in the spectrum of coherent phonons. The phase transition has the characteristic time of 5–10 ps, and the lifetime of each phase is on the order of 40–60 ps. In addition, phonon density of states, simulated by molecular dynamics, together with third-harmonic time-resolved spectra prove that laser-excited phonons in a bulk of dielectrics are generated by displacive excitation (DECP) mechanism in plasma mediated conditions.
With the advent of femtosecond ultrashort laser systems, they have become widespread in fundamental science and the practical field. The femtosecond duration of the laser pulse allows to strongly localize the effect, thereby opening the way for new femtotechnologies for three-dimensional micro- and nano-processing of the bulk of various materials, competing with expensive and proven methods of electron lithography1. Furthermore, tight focusing of femtosecond low-energy (up to 10 μJ) laser radiation into the volume of a condensed medium creates the extreme states of matter that occur in stars and can thus be observed in laboratory conditions2.
The evolution of the extreme state of a matter created by laser radiation is a cascade of complex processes, ranging from the formation of a non-equilibrium laser microplasma, its thermalization with the transfer of energy from the electronic subsystem to the ion core, generation of shock waves, and the formation of residual micromodifications3. These processes occur on a wide time scale from femtoseconds to microseconds, and the final state strongly depends on the features of structural rearrangements in the substance at each of the time stages4,5. Understanding microscopic processes in materials and devices that can be switched by light requires experimental access to dynamics on nanometer length and femtosecond time scales that enable new paths to material processing2, phase transitions6,7, and material properties manipulations8,9. The extreme impact of the laser pulse on the matter induces phase transitions that change its structural and physical properties4,10. The structural properties of matter are primarily dependent on the symmetry. The symmetry change that could be produced under phase transition opens the way to manipulate material properties on ultrafast (up to sub-ps and fs) timescale. Thereby, it is essential to characterize the dynamics of phase transitions11.
Traditionally, changes in the structure of a substance can be diagnosed using X-ray methods12,13,14. However, the sub-picosecond time resolution is only possible with mega-science facilities such as synchrotron radiation sources or free-electron lasers15. The development of optical methods (such as Raman scattering16,17 or phonon spectroscopy7,18,19 for registering structural changes in matter is an essential direction in developing express methods of structural diagnosis with high spatial and temporal resolution. This method can be used and spread more widely, making the developed approach highly relevant to condensed matter physics11,20. The generation of coherent phonons upon laser excitation is an important indicator of changes in the structure of a substance11,21. The registration of phonon vibrations spectrum in time sheds light on the initiated fast phase transitions. It can provide essential information on the formation of new, stationary and non-stationary phases of matter, which significantly expands the understanding of phase diagrams of matter11. In this Letter, we have retrieved the dynamics of ultrafast phase transitions in MgF2 from the change of phonon spectrum probed by pump-probe third-harmonic generation technique. The following study will help better understand ultrafast processes during energy transfer from electron plasma to lattice under high deposited energy densities.
At ambient conditions, MgF2 is in a tetragonal rutile type structure with space group P42/mnm, which is the same as TiO2 and SiO222. Under pressure of 9.1 GPa, it transforms to the orthorombic CaCl2-type phase with a space group Pnnm23. In MgF2 crystal, the phase transition at room temperature occurs under static pressure loading22. Under further compression, MgF2 transforms to the modified fluoride structure (PdF2 type, with space group Pa-3) at 14 GPa, and then the cotunnite structure (α-PbCl2 type, with space group Pnam) at 35 GPa22. The lattice structure of MgF2 is presented in Fig. 1. It is essential to mention that under high pressure, the phase transition should occur in a cascade manner, as it is demonstrated for the shock-wave induced phase transitions in Si: α-diamond ⇒ β-Sn ⇒ Imma ⇒ simple hex24,25,26. The possible paths of atom movement during solid-to-solid phase transitions are also marked in Fig. 1.
Atoms position inside the MgF2 lattice in different phases. Arrows indicate the atom shifts during phase transitions. Mg is imaged as a blue sphere and F as a green.
The pressures achieved at the wavefront of the laser-induced shock wave could overcome the threshold pressure. We performed shadow photography measurements to retrieve the shock wave pressure (see “Methods” section). We obtained shadow photographs (see Fig. 2c) for different energies (the threshold energy of plasma formation is about 1 μJ). In the shadow photographs, the shock wave is observed as a divergent spherical wave (see Fig. 2a). We varied the energy from 1 μJ up to 14 μJ. By approximating the dependence of the shock wave diameter on time (see the inset in Fig. 2a) with an exponentially decaying function, we obtained the dependence of the shock wave velocity on energy. Then, using shock adiabat (see “Methods” section), we restored the pressures for given laser energies (see Fig. 2b). In experiments, we achieved pressure higher than 50 GPa with damping rate ~ 10 ns (on the moving shock front). Initially the shock wave velocity under such conditions is about 12 km/s (for 4 μJ pump pulse).
(a) Pressure at the shock front as a function of the laser energy normalized to the threshold plasma formation energy (1 μJ for our experimental conditions). The colored dotted lines show the pressures at which the phase transition is observed under static conditions. (b) An example of the shock wave diameter evolution registered in shadow photography experiments and the pressure on the shock wavefront. (c) Example of shock wave shadow photographs.
According to the phase diagram of MgF222, the pressures achieved in the experiments are higher than the threshold values for low-temperature phase transitions in MgF2 (Pnam, Pa-3, Pnmn) under static conditions. These phase transitions are characterized by the change of the lattice symmetry due to the decrease of interatomic distances. Under these conditions, the interatomic potential changes abruptly, determining new atomic positions. The presented in Fig. 2 threshold pressure is determined under static conditions when the new phase can exist infinitely in time. However, the shock wave propagation induces ultrafast (with a characteristic time of a picosecond) change in pressure, and the achieved pressures could not serve as a trustworthy indicator of a phase transition27.
Moreover, in the dynamical case, the phase transition is metastable and could be reversible, returning to its original state after the passage of the shock wave. If the new phase exists, it has to affect the macroscopic properties of the medium. To verify the possibility of phase transitions under laser impact, we retrieve the phonon spectrum of MgF2, which is determined by the lattice symmetry.
To retrieve the phonon spectrum, we recorded the dependence of the third harmonic (TH) signal on the time delay between the pump and probe laser pulses (see Fig. 3). The phonons, induced by the displacement of ions in lattice, drive the modulation of the third-order susceptibility χ(3)28:
where χ0(3) is the third-order susceptibility of the unperturbed sample, Q is the displacement of atoms in a lattice and \(\left( {\frac{{\partial \chi^{(3)} }}{\partial Q}} \right)\) is the polarizability tensor. The efficiency of third harmonic generation (THG) process is proportional to:
where ∆n is the refractive index difference between the third harmonic and the fundamental pulse29. Moreover, this refractive index difference between (without plasma impact) does not depend on phonon frequency \(\Delta n = \Delta n_{0} + \left[ {n_{2} \left( {\lambda_{1} } \right) - \,n_{2} \left( {\lambda_{3} } \right)} \right] \times I\) , where I is laser pulse intensity, ∆no is phase mismatch of unaffected MgF2, n2 is nonlinear refractive index λ1 = 1240 nm and λ3 = 413 nm. Thereby the third harmonic signal would be modulated on the phonon frequency, and the permitted phonon modes would be determined by \(\left( {\frac{{\partial \chi^{(3)} }}{\partial Q}} \right)\) tensor.
Dependence of the third harmonic signal on the time delay between the pump and probe pulses. Bottom black lines show the spectra of this signal in time windows, marked with a colored dotted line. The red lines show the phonon spectra (multiplied by the ω−3 function), obtained in numerical simulations for pressures of 1 bar, 30, and 50 GPa. The modes in the simulated spectra are marked if they are Raman or IR active.
Compared with the experiments on the phonon generation in Sb230, the peak fluences (~ 104 J/cm2) and intensities (~ 1013 W/cm2) in our experiments are about four orders higher, leading to the electron plasma ignition. According to the femtosecond laser-matter interaction regime, laser pulse energy cannot be transferred to the lattice directly, and phonons are generated with a delay of ~ 10 ps. As a result, the laser-induced plasma acts as a broadband excitation (close to delta function) for coherent phonons. The delay for the coherent phonons generation corresponds to the characteristic time of electron–phonon coupling28. For stationary case it could be estimated as \(\frac{1}{{\tau_{k} }} \approx \frac{{\left| \xi \right|^{2} m_{c} \theta k}}{{\pi \rho c_{s}^{2} \hbar^{3} }}\) , where ξ is deformation potential constant, mc is the effective mass of the electron, θ = kbT (T is the crystal lattice temperature, kb is Boltzmann constant), k is the wave vector of electrons in the plasma, ρ is the MgF2 density and cs is the speed of sound. In our estimation, we used the following parameters |ξ |~ 11 eV (an order of band gap), mc ~ 0.5 me, where me is electron mass, T = 300 K, ρ = 3.15 g/cm3, cs = 3.5 km/s, the wavevector of the electrons was estimated from the kinetic energy estimate of 4 eV28. Under such parameters, the energy transport time from electrons to phonons could be estimated as 6 ps.
There is also a nonzero component in the third harmonic intensity on time; it could be caused by such processes as shock wave generation (an increase of density on the shock wavefront will lead to the change in χ(3) and third harmonic generation) and micromodification formation (decrease of density in the center will lead to the change in χ(3) and third harmonic generation). The shock wave will leave the laser spot (~ 2 μm radius) in about 150 ps, which corresponds to the disappearance of coherent phonons.
In the pump-probe experiments (see “Methods” section), the energy of the pump pulse was about 2 μJ, which equals about two thresholds of plasma formation. The energy of the probe pulse was about 100 nJ, more than an order less than the pump energy. The third harmonic signal could be divided into three regions over time (see Fig. 3): 10–40 ps, 50–90 ps, and 100–140 ps. We performed the spectral analysis of these regions. In the first region, the signal at a frequency of 1.2 THz is distinguished, which, with an increase in the time delay, changes its frequency, first by 0.9 THz (second region), and then down to 0.8 THz (third region). In addition, novel spectral components start to appear in the spectrum.
The abrupt change of the spectrum results from the change in the lattice symmetry. If the distance between the atoms changes monotonically, the frequency of the phonon mode will change similarly. In contrast, we obtain a “jump” in the frequency. The jump could only result from the phase transition that leads to the change of the lattice symmetry, which determines permitted phonon modes. However, due to the limits of the method (recorded spectrum lies in the range up to 4 THz), the spectrum itself does not give direct information about phase transitions and achieved pressures.
Interestingly, the obtained modes are limited by 4 THz. That is followed by the generation mechanism, namely displacive excitation of THz phonons (DECP)31. Following the approach, the displacement of atoms Q has an inversive quadratic dependence on phonon frequency (ω) \(Q\sim 1/\omega^{2}\) . On the other hand, \(\left( {\frac{{\partial \chi^{(3)} }}{\partial Q}} \right)\) tensor has a square root dependence on phonon frequency \(\left( {\frac{{\partial \chi^{(3)} }}{\partial Q}} \right)\sim \sqrt \omega\) 32. Thereby the efficiency of the third harmonic generation has an inverse cubic dependence on frequency:
Following Eq. (3), the phonon spectrum, retrieved from the molecular dynamics (MD), was normalized on the ω−3 function. Interestingly, in the framework of the alternative approach, describing the phonon generation mechanism of impulsive stimulated Raman scattering (ISRS)33, the displacement of atoms is inversely proportional to frequency Q ~ 1/ω, but we obtained the inverse cubic dependence that corresponds to a DECP mechanism of phonon excitation. Comparing results of numerical simulation of the phonon spectrum (considering that the intensity of the third harmonic signal is inversely proportional to the cube of the vibration frequency that corresponds to the best coincidence) with experimental data allows us to estimate achieved pressures and identify the phases of MgF2 during laser impact (see Fig. 3). The first region (10–40 ps) corresponds to the P42/mnm phase and atmospheric pressure, the second region (50–90 ps) to Pa-3 phase (30 GPa), and the third region (100–140 ps) to Pnam (50 GPa), respectively. Thus, it can be argued that the change in the characteristic modes of phonon vibrations, observed experimentally under femtosecond laser irradiation is caused by a rearrangement of the crystal lattice.
On time scales up to about 40 ps, the crystal lattice still corresponds to an unexcited medium. Then, the amplitude of the phonon vibrations drops significantly. Further, the amplitude of the TH increases, and the fundamental frequency shifts from 1.2 to 0.9 THz. As a visual analogy of the process, one can imagine the two-dimensional array of nonlinear coupled oscillators. After external force is applied, they will start to oscillate on eigenfrequencies. Then the part of oscillators abruptly changes its eigenfrequency; moreover, the fraction of oscillators at the new frequency increases (due to the phase transition). Under such a process, the amplitude of oscillations will decrease because the energy is not dissipated after a short period, the oscillations on the new frequency will re-appear. The characteristic time of starting a new oscillation should correspond to phonon–phonon scattering time ~ 10 ps34. A similar process occurs when the frequency shifts from 0.9 to 0.65 THz. Thereby we can claim that after 50 ps after laser-plasma excitation, a phase transition P42/mnm ⇒ Pa-3 occurs. The atoms' shift during phase transition is presented in Fig. 3. This phase is metastable and exists during about 50 ps, after which there occurs phase transition Pa-3 ⇒ Pnam. The last phase also exists over about 50 ps, after which the crystal lattice is finally destroyed (the micromodification is formed), and the third harmonic signal disappears accordingly. An increase in the laser pulse energy (up to 4 μJ) leads to the more rapid destruction of the material and a decrease in the number of observed phases is taken place.
Remarkably, we did not observe phonon spectrum specific for the Pnmn phase. It could be caused by a relatively narrow range of pressures, where this phase is stable (see Fig. 2a). The spectrum retrieved from the MD corresponds to the phonon density of states of MgF2. This spectrum indicates all possible frequencies of atom oscillations, and some of them could not be permitted due to zero components of \(\left( {\frac{{\partial \chi^{(3)} }}{\partial Q}} \right)\) the tensor. The MD could also give information about Raman and IR activity of the obtained modes. We calculated the autocorrelation functions of dipole moments (their projections on cartesian coordinates). If a vibration changes, the molecular dipole moment (d), it is IR-active. If it leads to the change in the polarizability, it is Raman active (for example, dx2 + dy2 and z or xy components are nonzero). We marked the IR and Raman activity of the phonon in Fig. 3. In the experiment, we could only observe phonons that are Raman-active. The analysis of the diagonal components of observed modes shows that probably 1.2 THz mode (10–40 ps) is A1g (diagonal components of Raman tensor a, a, b), the mode with 0.9 THz is B1g (40–100 ps), and 0.65 THz (100–140 ps) is Eg mode. The performed numerical simulation was performed at 300 K and did not consider the possible increase of pressure caused by the laser heat of the lattice. Nevertheless, there is a good coincidence between the simulated and experimental data. That could be caused by low thermodiffusion rates (thermodiffusion time ~ 1–10 μs); thereby, the high temperature is observed only in the area of micromodification and does not influence phase transitions.
Applying Fourier analysis with a sliding time window (~ 0.1 ps) to analyze the time signal of the third harmonic, it is possible to visualize the temporal dynamics of phase transitions more clearly (see Fig. 4). The phase transition P42/mnm ⇒ Pa-3 is accompanied by a smooth increase in frequency. The smooth decrease in frequency (with time delays of more than 140 ps) precedes the formation of micromodification, which corresponds to the lattice destruction. It is also worth noting that the spectrum of phonon vibrations changes abruptly in other cases. At phase transitions, a jump-like change in the spectrum is probably caused by the fact that the system’s symmetry changes during phase transitions. With a “new” crystal lattice, “old” vibrations become impossible, leading to a redistribution of energy between phonon modes with a close frequency, similar to the process observed in Ref.35. The apparent increase in the signal amplitude in the time dependence is caused by the fact that at higher pressure, the phonon oscillations have a lower frequency, which means that the third harmonic signal will increase.
Heat map of the spectrum of phonon vibrations obtained as Fourier transform of the third harmonic signal evolution with a sliding time window. Logarithmic scale.
To sum up, with the help of the third harmonic-based pump–probe technique and molecular dynamics simulation, we reconstructed the dynamics of laser-plasma induced phase transitions in magnesium fluoride. First, the tightly focused femtosecond laser pulse generates non-equilibrium electron plasma. Once after its relaxation (transferring the energy to the lattice), the coherent phonons at the 1.2 THz frequency are excited. The phonon spectrum at the initial stage of energy transfer processes corresponds to the unperturbed MgF2 sample. Then due to a sub-TPa shock wave generation, a cascade of picosecond phase transitions P42/mnm ⇒ Pa-3 ⇒ Pnam occurs. The transition between phases has the characteristic time of 5–10 ps. The phase lifetime is on the order of 40–60 ps; after 140 ps, lattice disordering is observed, and the coherent phonons disappear. In addition, phonon density of states, simulated by molecular dynamics, together with third-harmonic time-resolved spectra prove that laser-excited phonons in a bulk of dielectrics are generated by displacive excitation (DECP) mechanism under plasma mediated conditions.
To retrieve the dynamics of the phase transitions, we used the third harmonic-based pump–probe technique. Under tight focusing of laser radiation inside the unperturbed isotropic medium, the efficiency of third-harmonic generation equals to zero. However, any inhomogeneities violated the destructive interference. The coherent phonons (atoms oscillations) act as inhomogeneities in our case. It should be noted that the intensity of the third harmonic signal, recorded in the experiment, is proportional to ω−3 (see above), where ω is the phonon frequency is 10 Hz). In the experiments, we used Cr: Forsterite laser system (central wavelength is 1240 nm, pulse duration is about 100 fs, repetition rate is 10 Hz). Using Michelson interferometer with one moving mirror, we divided the initial laser pulse onto pump and probe pulses (both pulses have the same wavelength)21.
The pump pulse (energy is varied from 1 μJ up to 14 μJ) generates plasma (electron density ~ 1018 cm−3), while the second acts as a probe (laser energy is 100 nJ). Both pulses are tightly focused into the bulk of magnesium fluoride using a Thorlabs A240 TM lens (NA = 0.5 f = 8 mm). The resulting fluence and intensity inside the sample are ~ 104 J/cm2 and 1013 W/cm2 respectively. The laser radiation is focused into a spot with a diameter of ~ 4 μm. The sample was automatically moved using motorized translation stages synchronously to the laser pulse repetition rate. Under such experimental conditions, a single-pulse interaction regime with the sample was established. The energy of the pump and probe pulses is controlled using calibrated photodetectors. A Michelson interferometer was assembled to change the time delay between the laser pulses. By varying the length of one interferometer arm, we changed the time delay between the pump and probe pulses. In one of the arms of the interferometer, a quarter-wave plate was placed. A double pass through the plate rotated the polarization by 90 degrees. A polarizer was placed in front of the detection system, adjusted in such a way as to transmit only a probe pulse.
In contrast to the traditional pump-probe schemes, in our experiments, the signal of non-synchronous third harmonic (TH), generated on coherent phonons, was recorded. The TH signal was registered by a photomultiplier tube (Hamamatsu H10722-113). As we have shown in our publications earlier, this technique is much more effective than the traditional one, within which the radiation is recorded at the fundamental wavelength29,36. A detailed description of the experimental setup can be found elsewhere29,36,37.
Time-resolved shadow photography was used to determine the achieved laser-induced pressure in the bulk of magnesium fluoride. The main aim of the technique is to determine the shock wave velocity that determines shock wave pressure. As in the THG technique, the temporal resolution is achieved by changing the optical path between the pump and probe pulses. Within the framework of this technique, the second harmonic of the Cr:Forsterite laser radiation (620 nm), generated in a type I BBO crystal, was used as a probe pulse. As in the case of probing by the third harmonic, the pump pulse is tightly focused into the sample; however, the probe pulse was preliminarily scattered on the diffusion plate and directed opposite to the pump pulse. As a result, the focusing lens collects the signal from the waist on the CCD camera. The probe pulse passing through the irradiated area undergoes refraction on the changes in the refractive index induced in the sample during the propagation of the laser-induced shock wave, which in shadow photographs looks like dark areas. We determined the shock wave diameter for each time delay and fitted the resulted function as an exponentially decaying function with a dumping parameter α38. From the fitting parameters, we determined the shock wave velocity and dumping parameter. Because the pressure and the velocity of the shock wave are uniquely related by the equation of the shock adiabat, then using the velocity of the shock wave at the initial moment (calculated from several shadow photographs taken for different time delays), it is possible to estimate the pressure of the shock wave.
where p is shock wave pressure, p0 is the pressure in the unperturbed material, ρ0 is density, Us is shock wave speed, U is particle velocity c0, and λ are empirical constants c0 = 3.5 km/s and λ = 1.56 km/s39. Then we calculated the pressure on the shock wavefront for each time delay and pump energy, as it is shown in Fig. 2.
The LAMMPS software package was used for modeling40 the phonon spectrum of MgF2. The applied potential for fluorides was presented in the paper41. This model is valid up to the pressures of 130 GPa. However, in a numerical experiment, the destruction of the grid was observed at pressures above 100 GPa. Within the simulation framework, one supercell of magnesium fluoride with a size of 100 by 100 per 100 atoms was created. External stationary pressure was applied to each boundary. The simulation was carried out in the n-p-T ensemble. Before the simulation, the system was equilibrated at room temperature and atmospheric pressure (1 bar). Since the phase transition does not occur in a system with ideal periodic boundary conditions, 100 atoms, randomly distributed along with the simulation cell, were removed. They served as seed centers for the new phase. The phonon (vibrational–rotational) spectrum was calculated as follows. First, the velocity autocorrelation function was calculated for 1,000,000 steps (time step 0.1 fs). Then the Fourier transform of the autocorrelation function was performed. The spectrum obtained in this way is the rotational-vibrational spectrum of the test substance. The simulations were performed with a 10 GPa step. The phonon spectra obtained in this way agree with the tabulated data. Figure 5 below demonstrates the comparison of the calculated phonon spectrum with the estimated one from the time-domain density functional theory. We also simulated the instantaneous value of dipole moment (d) and its protections on cartesian coordinates. For this purpose, we added the charge of each atom to the measurement of the autocorrelation function. To determine if the obtained mode is IR or Raman active, we plotted the FFT spectrum of (d, dx2 + dy2, dx × dy, and dz). If there is a nonzero component in d the mode is IR active and if dx2 + dy2 and dz, or dx × dy are nonzero than the mode is Raman active.
Phonon density of states retrieved from MD simulations (solid line) and reproduced from Ref.42.
Data is available upon reasonable request.
Schaffer, C. B., Brodeur, A. & Mazur, E. Laser-induced breakdown and damage in bulk transparent materials induced by tightly focused femtosecond laser pulses. Meas. Sci. Technol. 12, 1784–1794 (2001).
Article ADS CAS Google Scholar
Potemkin, F. et al. Controlled energy deposition and void-like modification inside transparent solids by two-color tightly focused femtosecond laser pulses. Appl. Phys. Lett. 110, 163903 (2017).
Lenzner, M. et al. Femtosecond optical breakdown in dielectrics. Phys. Rev. Lett. 80, 4076–4079 (1998).
Article ADS CAS Google Scholar
Mizeikis, V. et al. Synthesis of super-dense phase of aluminum under extreme pressure and temperature conditions created by femtosecond laser pulses in sapphire. Proc. SPIE 8249, 82490A (2012).
Juodkazis, S. et al. Femtosecond laser-induced confined microexplosion: Tool for creation high-pressure phases. MRS Adv. 1, 1149–1155 (2016).
Rini, M. et al. Control of the electronic phase of a manganite by mode-selective vibrational excitation. Nature 449, 72–74 (2007).
Article ADS CAS Google Scholar
Kim, K. W. et al. Ultrafast transient generation of spin-density-wave order in the normal state of BaFe2 As2 driven by coherent lattice vibrations. Nat. Mater. 11, 497–501 (2012).
Article ADS CAS Google Scholar
Afanasiev, D. et al. Ultrafast control of magnetic interactions via light-driven phonons. Nat. Mater. 20, 607–611 (2021).
Article ADS CAS Google Scholar
Maehrlein , S. , Paarmann , A. , Wolf , M. & Kampfrath , T. Terahertz sum-frequency excitation of a Raman-active phonon .Phys.Rev. Fr.Lett.Rev. 119, 1–6 (2017).
Gamaly, E. G. The physics of ultra-short laser interaction with solids at non-relativistic intensities. Phys. Rep. 508, 91–243 (2011).
Wall, S. et al. Ultrafast changes in lattice symmetry probed by coherent phonons. Nat. Commun. https://doi.org/10.1038/ncomms1719 (2012).
Mcbride, E. E. et al. Phase transition lowering in dynamically compressed silicon. Nat. Phys. 15, 89–94 (2019).
Article MathSciNet CAS Google Scholar
Kluge, T. et al. Nanometer-scale characterization of laser-driven compression, shocks, and phase transitions, by X-ray scattering using free electron lasers. Phys. Plasmas 24, 102709 (2017).
Kirschner, M. S. et al. Photoinduced, reversible phase transitions in all-inorganic perovskite nanocrystals. Nat. Commun. 10, 1–8 (2019).
Potemkin, F. V. et al. Hybrid X-ray laser-plasma/laser-synchrotron facility for pump-probe studies of the extreme state of matter at NRC ‘kurchatov Institute’. Rev. Sci. Instrum. 92, 1–9 (2021).
Hortensius, J. R., Afanasiev, D., Sasani, A., Bousquet, E. & Caviglia, A. D. Ultrafast strain engineering and coherent structural dynamics from resonantly driven optical phonons in LaAlO3. NPJ Quantum Mater. 5, 1–6 (2020).
Pawbake, A. et al. Pressure-induced phase transitions in Germanium telluride: Raman signatures of anharmonicity and oxidation. Phys. Rev. Lett. 122, 1–5 (2019).
Ionin, A. A., Seleznev, L. V. & Sunchugasheva, E. S. Controlling plasma channels through ultrashort laser pulse filamentation. Proc. SPIE 8898, 88980Z (2013).
Lindenberg, A. M. et al. Time-resolved X-ray diffraction from coherent phonons during a laser-induced phase transition. Phys. Rev. Lett. 84, 111–114 (2000).
Article ADS CAS Google Scholar
Siegal, Y., Glezer, E. N., Huang, L. & Mazur, E. Laser-induced phase transitions. Annu. Rev. Mater. Sci. 25, 223–247 (1995).
Article ADS CAS Google Scholar
Potemkin, F. V., Mareev, E. I., Mikheev, P. M. & Khodakovskij, N. G. Resonant laser-plasma excitation of coherent THz phonons under extreme conditions of femtosecond plasma formation in a bulk of fluorine-containing crystals. Laser Phys. Lett. 10, 076003 (2013).
Article ADS CAS Google Scholar
Zhang, T., Cheng, Y., Lv, Z.-L., Ji, G.-F. & Gong, M. Phase transition and thermodynamic properties of magnesium fluoride by first principles. Int. J. Mod. Phys. B 28, 1450026 (2014).
Haines, J. et al. X-ray diffraction and theoretical studies of the high-pressure structures and phase transitions in magnesium fluoride. Phys. Rev. B Condens. Matter Mater. Phys. 64, 1341101 (2001).
Mcbride, E. E. et al. Phase transition lowering in dynamically-compressed silicon. Nat. Phys. 15, 89 (2019).
Article MathSciNet CAS Google Scholar
Brown, S. B. et al. Direct imaging of ultrafast lattice dynamics. Sci. Adv. https://doi.org/10.1126/sciadv.aau8044 (2019).
Article PubMed PubMed Central Google Scholar
Wippermann, S., He, Y., Vörös, M. & Galli, G. Novel silicon phases and nanostructures for solar energy conversion. Appl. Phys. Rev. 3, 040807 (2016).
Mareev, E. I. & Potemkin, F. V. Dynamics of ultrafast phase transitions in (001) Si on the shock-wave front. Int. J. Mol. Sci. 23, 2115 (2022).
Potemkin, F. V. & Mikheev, P. M. Efficient generation of coherent THz phonons with a strong change in frequency excited by femtosecond laser plasma formed in a bulk of quartz. Eur. Phys. J. D 66, 248 (2012).
Gordienko, V. M., Potemkin, F. V. & Mikheev, P. M. Evolution of a femtosecond laser-induced plasma and energy transfer processes in a SiO2 microvolume detected by the third harmonic generation technique. JETP Lett. 90, 263–267 (2009).
Article ADS CAS Google Scholar
Stevens, T., Kuhl, J. & Merlin, R. Coherent phonon generation and the two stimulated Raman tensors. Phys. Rev. B 65, 3–6 (2002).
Zeiger, H. J. et al. Theory for displacive excitation of coherent phonons. Phys. Rev. B 45, 768–778 (1992).
Article ADS CAS Google Scholar
Penzkofer, A., Laubereau, A. & Kaiser, W. High intensity Raman interactions. Prog. Quantum Electron. 6, 55–140 (1979).
Merlin, R. Generating coherent Thz phonons whith light pulses. Solid State Commun. 102, 207–220 (1997).
Gibbons, T. M., Bebek, M. B., Kang, B., Stanley, C. M. & Estreicher, S. K. Phonon-phonon interactions: First principles theory. J. Appl. Phys. 118, 085103 (2015).
Zhang, L., Wang, Y., Cui, T., Ma, Y. & Zou, G. First-principles study of the pressure-induced rutile-CaCl2 phase transition in MgF2. Solid State Commun. 145, 283–287 (2008).
Article ADS CAS Google Scholar
Mareev, E. I., Migal, E. A. & Potemkin, F. V. Ultrafast third harmonic generation imaging of microplasma at the threshold of laser-induced plasma formation in solids. Appl. Phys. Lett. 114, 031106 (2019).
Potemkin, F. V. & Mikheev, P. M. Efficient generation of coherent THz phonons with a strong change in frequency excited by femtosecond laser plasma formed in a bulk of quartz. Eur. Phys. J. D. https://doi.org/10.1140/epjd/e2012-30048-x (2012).
Mareev, E. I., Rumiantsev, B. V. & Potemkin, F. V. Study of the parameters of laser-induced shock waves for laser shock peening of silicon. JETP Lett. 112, 739–744 (2020).
Article ADS CAS Google Scholar
Batsanov, S. S. Features of solid-phase transformations induced by shock compression. Russ. Chem. Rev. 75, 601–616 (2006).
Article ADS CAS Google Scholar
Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 1–42 (1995).
Article ADS CAS Google Scholar
Barreta, C. D. et al. Ionic solids at elevated temperatures and high pressures: MgF2. J. Chem. Phys. 107, 4337–4344 (1997).
Petretto, G. et al. High-throughput density-functional perturbation theory phonons for inorganic materials. Sci. Data 5, 1–12 (2018).
This work in a part of laser excitation of dielectrics was supported by Russian Science Foundation (Project No. 17-72-20130). The research in a part of coherent phonons excitation was funded by Russian Foundation for Basic Research (RFBR) (Projects No. 18-02-40018, 19-29-12037) and by Russian Foundation for Basic Research (RFBR) and Moscow city Government according to the project No. 21-32-70021.
These authors contributed equally: Evgenii Mareev and Fedor Potemkin.
Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow, Russia, 119991
Evgenii Mareev & Fedor Potemkin
Institute of Photon Technologies of Federal Scientific Research Centre “Crystallography and Photonics” of Russian Academy of Sciences, Troitsk, Russia, 108840
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
The authors contributed equally to this work.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Mareev, E., Potemkin, F. Dynamics of ultrafast phase transitions in MgF2 triggered by laser-induced THz coherent phonons. Sci Rep 12, 6621 (2022). https://doi.org/10.1038/s41598-022-09815-4
DOI: https://doi.org/10.1038/s41598-022-09815-4
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Mechanically Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.