Introduction

Epitaxy can be used to impose misfit strain capable of altering the properties of materials. Notable examples include the enhancement of ferroelectric and ferromagnetic order1 and even the engineering of artificial multiferroics at room temperature2. Whereas static strain engineering is a well-established paradigm3,4,5, ultrafast strain engineering has emerged only recently as an effective method to manipulate functional properties of oxides6,7, control collective excitations8,9,10,11, induce changes in the band topology12,13, and drive optoelectronic phase switching14,15. The generation of strain pulses traditionally relies on opto-acoustic conversion processes either in the functional material itself or in opto-acoustic transducers, often involving electronic excitation16. An inherently different approach is to use ultrashort pulses of light that are tuned in resonance with an infrared-active atomic vibration of a substrate, in order to transform the structural and electronic properties of an epitaxial thin film14. This mechanism, applied extensively to insulating lanthanum aluminate (LaAlO3) substrates, governs ultrafast metal–insulator transitions14, ultrasonic magnetic dynamics17, and sonic lattice waves18 in various thin films of strongly correlated oxides. However, the nature of the nonlinear lattice dynamics initiated in the substrate material is not yet fully understood. As strain is naturally related to dynamics of the crystal lattice, here we study light-induced ultrafast lattice dynamics directly in an LaAlO3 substrate. We show that impulsive optical excitation at the photon energy tuned in resonance with a polar stretching of the Al–O bonds drives a non-polar rotational mode of oxygen octahedra via ionic Raman scattering19. The anisotropic optically induced strain also generates propagating longitudinal (LA) and transverse acoustic (TA) wave packets. Importantly, shear strain wave packets are found to be produced with extraordinary efficiency close to the phonon resonance. These results uncover an hitherto unknown microscopic feature of ultrafast strain engineering that opens wide perspectives for material control via optically tunable strain.

Results and discussion

Experiment

We investigate light-induced structural dynamics in LaAlO3, an insulating substrate utilized extensively in oxide electronics for the epitaxy of correlated materials, including high-Tc cuprate superconductors20, magnetoresistive manganites21, and nickelates22. At room temperature, LaAlO3 exhibits a distorted perovskite structure (rhombohedral space group \(R\bar{3}c\), see Fig. 1a). To resonantly drive the lattice vibrations in (001) LaAlO3 single crystals, we use ultrashort pump pulses in the mid-infrared (mid-IR) frequency range. We tune the photon energy of the pump pulses continuously across the closely lying highest-frequency Eu and A2u infrared-active phonon resonances23,24. For the schematics of atomic motion corresponding to the Eu and A2u mode, see Fig. 1b and Supplementary Fig. 3, respectively. The energy was tuned in the experimentally accessible range of 70–180 meV (17–41 THz) and allowed us to compare dynamics excited in the optical transparency window (>130 meV) with structural transient dynamics induced by pulses tuned in resonance with the lattice vibrational modes centered at 81 meV. A high efficiency of ionic Raman scattering is anticipated in LaAlO3 as mutual coupling between optical phonon modes exists and even results in a small negative Lyddane–Sachs–Teller splitting between longitudinal and transverse optical modes25. The wide bandgap of LaAlO3 (5.6 eV26) as well as the absence of electronic in-gap states27 ensures the purely structural nature of the photo-induced response.

Fig. 1: Crystal structure and experimental geometry.
figure 1

a The crystal structure of LaAlO3, with indications of the different crystallographic directions. The rhombohedral distortion from the high-temperature cubic phase is due to out-of-phase rotations of the oxygen octahedra about the [111]pc axis. b Atomic motion corresponding to the infrared-active Eu stretching mode, polarized in the (111) plane in the x-direction. c Schematic illustration of the experimental scheme. The mid-IR pulse is generated by difference frequency mixing from two near-infrared pulses in a GaSe crystal, after which the mid-IR pulses are filtered by a germanium (Ge) filter. Following the mid-IR excitation, the ensuing changes in optical properties are probed with a time-delayed near-infrared pulse. The pump-induced changes to the reflection intensity ΔR and rotation ΔθB of the polarization plane are monitored.

To track the ensuing dynamics of the lattice, time-resolved optical reflectivity and birefringence measurements are performed using near-infrared probe pulses. The two complementary experimental techniques (for details, see “Methods”) are schematically illustrated in Fig. 1c. In the first scheme we monitor the transient differential reflectivity ΔR. The structural dynamics initiated by the pump pulse modulate the sample’s dielectric function resulting in a perturbation of the refractive index n, which is imprinted on the ΔR signal. In the second scheme we track the transient optical birefringence ΔθB. The phonon modes of LaAlO3 are intrinsically highly anisotropic, meaning that coherent dynamics of these modes can also modify the off-diagonal components of the permittivity tensor, thereby resulting in a transient birefringence.

Photo-induced structural dynamics

Measurements of transient changes to both the reflectivity and birefringence, using pump pulses at a photon energy tuned in resonance with the infrared-active phonon modes,  85 meV (21 THz), reveal multiple oscillatory responses at frequencies significantly below the one of the pump (see Fig. 2a, b). The highest-frequency oscillation is centered at 1.1 THz (4.3 meV) and is assigned to the Raman-active Eg soft mode of LaAlO324,28 associated with a rhombohedral instability of the \(R\bar{3}c\) lattice structure (see Supplementary Note 1). This mode comprises rotations of the oxygen octahedra around an axis perpendicular to the [111] pseudocubic direction as shown in Fig. 2c. The longer time delay further reveals oscillatory components at two discrete frequencies fTA and fLA in the GHz frequency range. Such a pattern originates from interference between light pulses reflected at the crystal surface and reflections from an acoustic wave front propagating into the bulk (Fig. 2c). In transparent materials, the frequency of the oscillations f is related to the refractive index n of the material at the probe wavelength27, the speed of sound vs, the angle θ w.r.t. the sample normal and the wavelength λ of the probe by the relation29 \(f=2n{v}_{{\rm{s}}}\cos (\theta )/\lambda\). In our experiments we vary the angle of incidence of the probe beam (the data are not shown here) and find that, while the frequency of the Eg mode remains unchanged, the frequency of the GHz oscillations decreases in agreement with the relation shown above. We extract the corresponding propagation velocities, obtaining vLA = 6.67 km s−1 and vTA = 4.87 km s−1, and find that they match the speed of LA and TA phonons in LaAlO3 propagating in the [001] direction30. This acousto-optic conversion is strongly peaked at a phonon wavenumber ka determined by the probe wavenumber31, ka = 3.2 × 105 cm−1. Our experiments show that optical excitation with ultrashort resonant mid-IR pulses initiates coherent structural dynamics in both the acoustic and optical branches of the phonon spectrum. Although the optical excitation of a broadband longitudinal acoustic wave front is expected from electrostriction in LaAlO3 32 and/or optical absorption16, the optical generation of shear strain requires the presence of an equilibrium or light-induced structural anisotropy. We discuss this aspect below, after the analysis of the THz Raman-active mode.

Fig. 2: Photo-induced structural dynamics.
figure 2

Transient changes in the intensity of the reflected probe pulse ΔR/R (a) and the polarization rotation ΔθB of the probe pulse (b) after excitation of the LaAlO3 substrate with a mid-infrared pump pulse of 86 meV (a) and 89 meV (b) and a fluence of 10 mJ cm2. The insets show the Fourier spectra of the signals. The top axis in b shows the distance z which the longitudinal sound wave has propagated at that time. c The real-space atomic motion corresponding to the excited Eg mode indicated with arrows (left) and a schematic picture of a strain wave propagating with speed vsound leading to interference between probe light (wavelength λ) reflected at the interface and scattered at the strain wave, depending on the distance z which the strain wave has propagated (right).

Coherent oxygen octahedra rotations

In our experiments only the lowest-energy Raman-active Eg optical phonon mode was observed. The frequencies of the remaining higher-energy modes (e.g. of Ag symmetry) lie outside the bandwidth of our excitation pulses and thus these modes cannot be excited in the experiment. To unveil the mechanism of excitation of the Eg mode, we vary the pump photon energy across the phonon absorption band. Figure 3a compares time-resolved transient reflectivity ΔR/R induced by pumping close to the phonon resonance (85 meV) with off-resonant pumping (124 meV), revealing a striking selectivity of the low-energy mode excitation. The inset in Fig. 3a shows that the amplitude of the excited Raman Eg mode measured for various photon energies increases strongly and peaks at 95 meV close to the absorption peak attributed to the Eu phonon mode centered at 81 meV. Recently, ionic Raman scattering (IRS) or nonlinear phononics was proposed as a mechanism for resonant non-thermal activation of coherent low-energy Raman-active (non-polar) phonon modes upon exciting infrared-active (polar) lattice vibrations19. According to this mechanism, the anharmonicity of the lattice potential leads to a nonlinear response when large-amplitude infrared-active ionic vibrations QIR are excited. The anharmonicity causes a short-living net distortion of the lattice along a Raman coordinate QR accompanied by coherent oscillations of the corresponding Raman-active mode around this displaced metastable position. In the general case, the nonlinear coupling can be described by introducing an invariant nonlinear term \({\alpha }^{\prime}{Q}_{{\rm{IR}}}^{2}{Q}_{{\rm{R}}}\) in the lattice potential, with \({\alpha }^{\prime}\) defining the strength of the coupling and Q corresponding to the normal coordinate of a phonon mode. Despite the strong correlation between the phonon absorption and the Eg mode amplitude, we note that the largest amplitude is observed at a pump photon energy shifted from the peak of the linear phonon absorption. Moreover, the lineshape of the amplitude as a function of the pump photon energy is significantly broader than the linear aborption. These observations indicate that the excitation of the Eg mode is more efficient at the reststrahlen band, where the absorption processes are not dominant and the optical response of the medium is non-dissipative.

Fig. 3: Light-induced excitation and dynamics of oxygen octahedra rotations.
figure 3

a Time-resolved reflectivity changes of the probe pulse at the central energy of 1.5 eV after mid-IR excitation at photon energies of 83 and 124 meV. Inset: the amplitude of the coherent phonon oscillations at 1.1 THz, as extracted from fitting an exponentially damped sine function to the transient reflectivity changes for different pump photon energies at constant incident fluence. The thick blue solid line indicates the linear absorption due to the Eu phonon mode, data taken from refs. 14,50. The spectral full-width at half-maximum of the excitation pulses is indicated by the shaded band. b Time-resolved measurements of the transient polarization rotation of the probe pulse after excitation with the pump pulses having a polarization of ϕ degrees with respect to the pseudocubic [100] axis at a pump energy of 89 meV. c Amplitude of the observed Eg mode excitation (square markers) as extracted by fitting sine functions to the data in b and calculated coupling of the excited phonons to the Eg mode (blue circles) vs pump polarization angle ϕ. The errorbars account only for the uncertainty of the fit.

To verify and study this nonlinear coupling in the specific case of LaAlO3, we perform a symmetry analysis and density functional theory (DFT) calculations with the ABINIT code33 to fit a nonlinear phonon-phonon model potential of bulk \(R\bar{3}c\) LaAlO3 (see “Methods” and Supplementary Discussion). Although the DFT calculations show that the coupling between the A2u and Eg mode is negligible, they confirm a strong coupling between the Eu and Eg mode. In the case of rhombohedral LaAlO3, containing a high-symmetry three-fold rotation axis along the pseudocubic [111] direction, the IR-active Eu mode has two orthogonal components (E\({\,}_{{\rm{u}}}^{x}\), E\({\,}_{{\rm{u}}}^{y}\)) oriented in the natural rhombohedral plane (pseudocubic (111)) (see inset Fig. 1a for a schematic including the different orientations). The coupling term is given by \({\alpha }^{\prime}=\alpha \cos (2\theta )\), with α a material-dependent constant and θ the angle between the laser polarization projected onto the (111) plane and the x-component of the Eu mode.

In Fig. 3c we show the DFT results for the evolution of the effective coupling coefficient \({\alpha }^{\prime}\) with respect to the pump polarization angle ϕ, oriented in the (001) plane, such that ϕ = 0 corresponds to the pump polarization oriented along the [100] axis. In order to calculate the coupling constant \({\alpha }^{\prime}\), the laser polarization in the pseudocubic (001) plane needs to be projected on the (111) plane. As a result \({\alpha }^{\prime}\) evolves as a non-trivial periodic function with extrema around ϕ = 0 and 112.5, see Supplementary Discussion. To verify this behavior, we measured the amplitude of the Eg oscillation for the pump polarizations oriented along several pseudocubic crystallographic directions, as shown in Fig. 3b. Figure 3c summarizes the observations, showing a good agreement with predictions of the nonlinear phonon model built from DFT. Together with the dependence on the pump photon energy, this confirms that excitation of the Eg Raman-active mode is governed by the IRS mechanism. Moreover, these selection rules are another strong indication that these lattice dynamics are not driven by the absorption, but are rather non-dissipative relying on the resonant enhancement of the scattering process19. The relevance of this mechanism is further corroborated by measurements of the fluence dependence, revealing a linear increase in the Eg mode amplitude34 (see Supplementary Discussion). These findings are therefore a clear manifestation of the efficient nonlinear phononics mechanism in a wide bandgap insulator in conditions promoting exclusively coherent phonon–phonon coupling. In this sense IRS differs substantially from regular impulsive stimulated Raman scattering in which excitation of coherent phonons is mediated by virtual electronic transitions34. Note that the observation of the net structural distortion along the Eg coordinate responsible for the oscillations is not feasible in an all-optical experiment alone and requires use of time-resolved X-ray diffraction35.

Tunable shear strain

The Fourier analysis (FFT) of the light-induced coherent strain waves as a function of the pump photon energy is summarized in Fig. 4a. We observe that acoustic waves are excited both in the optical transparency window and in the reststrahlen band. The inset of Fig. 4a shows that upon reaching the reststrahlen band, the amplitudes of both strain waves experience a pronounced growth indicating an enhancement at the phonon resonance. The dependence of the strain waves on the pump fluence and polarization are given in Supplementary Note 2. Conventionally, the generation of strain in the transparency window is described by electrostriction16,36. In principle, the enhancement of the TA and LA mode amplitudes close to the phonon resonance could simply originate from an increased photon absorption and the anisotropic elastic response of [001] LaAlO3, as this crystal cut is different from the high-symmetry [111] direction, providing a coupling between longitudinal and transverse strain. However, analyzing the strain wave amplitudes, we observe that the presence of the phonon resonance also dramatically renormalizes the ratio r between the amplitudes of the TA and LA modes (see Fig. 4b). Our experiments indicate that while their ratio r is around 0.2 in the transparency window, it shows a pronounced increase up to 2 in the reststrahlen band. As seen in the inset of Fig. 4a, the renormalization is achieved due to a strong enhancement in the generation of the TA strain wave, which seemingly occurs at the expense of the LA phonon mode. Such an evolution of the ratio with photon energy cannot be explained by direct TA mode excitation via dissipative effects due to an increased photon absorption, especially as the peak of the ratio does not coincide with the phonon absorption peak. Strikingly, the ratio peaks at a higher photon energy, comparable to the position of the maximum amplitude of the Eg mode (see inset Fig. 3a). We propose that this relative and absolute increase in TA mode generation is the result of a phonon-driven enhancement of the coupling between TA and LA modes, in combination with an ionic enhancement of the electrostriction37. In the Supplementary Discussion we quantify the elastic constants of LaAlO3 using DFT and show that, out of equilibrium, the distortion along the Eg Raman coordinate driven by the rectification of the phonon field alters the coupling between longitudinal and transverse strain. This change of elastic constants shows that the anisotropic Eg mode can transfer longitudinal strain into shear strain.

Fig. 4: Tunable shear strain.
figure 4

a Fourier transforms of the measured polarization rotation signals corresponding to strain wave propagation after excitation with different pump energies . Inset: Pump photon energy dependent amplitude of the two oscillations corresponding to LA and TA strain waves for fixed pump and probe polarizations as obtained from the Fourier transforms shown in panel a. The solid and dashed lines are guides to the eyes. The errorbars account for the uncertainty of the fit. b Linear absorption due to the infrared-active phonon modes, as taken from ref. 14 (left axis) and ratio r of the TA and LA mode amplitudes (right axis) vs photon energy at constant incident fluence. The black line serves as a guide to the eye. The spectral full-width at half-maximum of the excitation pulses is indicated by the errorbars.

Using a combination of time-resolved reflectivity and polarimetry we have studied coherent structural dynamics in LaAlO3 induced by ultrashort excitation of a selected IR-active phonon mode. Our experimental and theoretical analysis uncovers a previously unknown remarkable feature associated with the excitation of the crystal lattice. In addition to the displacement along a Raman coordinate and coherent THz atomic vibrations, expected within a nonlinear lattice excitation regime, we observe an efficient generation of shear strain wave packets. Shear strain following resonant pumping of the crystal lattice in LaAlO3 is likely to be a key element of the metal–insulator transitions, ultrasonic magnetic dynamics, and sonic lattice waves observed in recent years. Importantly, the ratio between the longitudinal and transverse strain waves can be tuned by the pump photon energy in vicinity of the phonon resonance, which hints at a close relation between nonlinear lattice dynamics, coherent lattice symmetry breaking, and (shear) strain generation. Tunable shear strain available on the ultrafast timescales via resonant lattice excitation can be exploited for material control using a wide array of perovskite wide-bandgap anisotropic substrates beyond LaAlO3. Since equilibrium shear strain is an important element for ferroelectric1, flexoelectric38, piezoelectric, and magnetoelectric effects, we envision opportunities for ultrafast manipulation of collective excitations in solids. Note: A recent complementary work by Neugebauer et al.39 compares coherent phonons generated by electronic and ionic Raman scattering in LaAlO3.

Methods

Sample and experimental setup

In our experiments we use commercially available 5 × 5 mm (001) LaAlO3 single crystals with a thickness of 0.5 mm from Crystec GMBH. The mid-infrared pulses (200 fs) are generated in a 0.35-mm-thick GaSe crystal by difference frequency mixing the output of two optical parametric amplifiers (OPAs). The mid-IR pulse is filtered from the OPA output with a germanium filter. The OPAs share the same white light, generated in a sapphire crystal, by the output of a laser amplifier (800 nm, 100 fs, 5 mJ, 1 kHz), which ensures carrier-envelope-phase-stability of the pulses40. A small part of the laser output is used to probe the structural dynamics. The transient differential reflectivity ΔR is monitored using a balanced photodetector. The optical birefringence θB is tracked using an optical polarization bridge (Wollaston prism) and a balanced photodetector. In both experimental configurations, the probe pulses were focused to a spot with a diameter of 80 μm. The spatial overlap between the pump and probe pulses is obtained by co-propagation of the beams, using an off-axis parabolic mirror, which focuses the pump beam to a spot with a diameter of about 150 μm.

DFT calculations

We simulated the \(R\bar{3}c\) phase of LaAlO3 through DFT41,42 as implemented in the ABINIT package (Ver 8.10.2)43,44. We used norm-conserving pseudopotentials45 to account for the interaction of the nuclei and the electrons. These pseudopotentials were downloaded from the Pseudodojo website46. For La we considered 5s, 5p, 5d, 6s, and 4f as valence states and for Al and O the valence states were considered to be 3s, 3p, and 2s, 2p, respectively. We used the PBEsol GGA functional for the exchange correlation interaction47 and all the calculations were done with a 5 × 5 × 5 mesh of k-points for sampling of reciprocal space and a cut-off energy on the plane wave expansion of 45 Hartree. To calculate the phonons, we used density functional perturbation theory as implemented in ABINIT.48,49