# Temporal and spatial tracking of ultrafast light-induced strain and polarization modulation in a ferroelectric thin film

Ruizhe Gu, Vincent Juvé, Claire Laulhé, Houssny Bouyanfif, Gwenaëlle Vaudel, Aurélie Poirier, Brahim Dkhil, Philippe Hollander, Charles Paillard, Mads Weber, et al.

► To cite this version:

Ruizhe Gu, Vincent Juvé, Claire Laulhé, Houssny Bouyanfif, Gwenaëlle Vaudel, et al. Temporal and spatial tracking of ultrafast light-induced strain and polarization modulation in a ferroelectric thin film. Science Advances, 2023, 9 (46), 10.1126/sciadv.adi1160. hal-04298208

HAL Id: hal-04298208

https://hal.science/hal-04298208v1

Submitted on 18 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
CONDENSED MATTER PHYSICS

# Temporal and spatial tracking of ultrafast light-induced strain and polarization modulation in a ferroelectric thin film

Ruizhe Gu<fnref n="1" />⁺, Vincent Juvé<fnref n="1" />*⁺, Claire Laulhé²,³, Houssny Bouyanfit⁴, Gwenaëlle Vaudel<fnref n="1" />, Aurélie Poirier<fnref n="1" />, Brahim Dkhil⁵, Philippe Hollander², Charles Paillard⁵,⁶, Mads C. Weber<fnref n="1" />, Daniel Sando⁷,⁸*, Stéphane Fusil⁹, Vincent Garcia⁹*, Pascal Ruello<fnref n="1" />*

Ultrashort light pulses induce rapid deformations of crystalline lattices. In ferroelectrics, lattice deformations couple directly to the polarization, which opens the perspective to modulate the electric polarization on an ultrafast time scale. Here, we report on the temporal and spatial tracking of strain and polar modulation in a single-domain $BiFeO_3$ thin film by ultrashort light pulses. To map the light-induced deformation of the $BiFeO_3$ unit cell, we perform time-resolved optical reflectivity and time-resolved x-ray diffraction. We show that an optical femtosecond laser pulse generates not only longitudinal but also shear strains. The longitudinal strain peaks at a large amplitude of 0.6%. The access of both the longitudinal and shear strains enables to quantitatively reconstruct the ultrafast deformation of the unit cell and to infer the corresponding reorientation of the ferroelectric polarization direction in space and time. Our findings open new perspectives for ultrafast manipulation of strain-coupled ferroic orders.

## INTRODUCTION

In piezoelectric and ferroelectric materials, the couplings between lattice distortions and the polar order are pivotal for key applications, such as smart electromechanical microelectronic components, spintronic memories, or advanced photonics (1–8). With the demand for high-rate data processing, the control of these couplings at ultrafast time scales, i.e., in the picosecond range, is desirable. One way to reach this ultrafast regime is the use of ultrashort laser pulses. For example, ultrashort light pulses have induced mechanical strains of up to 0.1% in ferroelectrics (9–11). These strain levels compare to commercial electrically driven piezoelectric sensors (1). Therefore, the ultrafast photoexcitation of piezoelectric and ferroelectric materials has triggered a large interest to control and modulate the strain (9–16) and to generate coherent acoustic (17–20) and optical phonons (21) with short light pulses. Even the switching of a ferroelectric polarization with short light pulses has recently been reported (6–8, 22, 23).

Since the polar order and crystal lattice are intrinsically coupled by the piezoelectric effect, a full understanding of the structural distortions induced by light pulses is crucial for an insight into the

light-driven manipulation of the polar order. To gain a full picture of ultrafast photo-induced distortions in ferroelectrics, access to the modulation of both the out-of-plane and in-plane components of the elastic distortions is essential. However, to date, several limitations exist. First of all, multidomain states, which are typical in ferroelectric materials, hamper this analysis, particularly in thin films with nanosize domains. The presence of various domains induces, indeed, an averaging effect on the photo-induced response, which precludes a straightforward description of the intrinsic interplay between light and lattice distortions. To obtain the clearest picture of these fundamental interactions, model systems are required, i.e., thin films consisting of a single ferroelectric domain. Second, most studies have been limited to the measurement of light-induced out-of-plane distortions precluding a full description of the light-induced lattice deformation (9, 10, 12, 14–16). In the present work, we first of all overcome the limitation of the multidomain complexity by fabricating a high-quality epitaxial single-domain $BiFeO_3$ film through strain and domain engineering (24–27). This archetypal system is a prerequisite for the understanding of the photo-induced response beyond bulk materials (11) or multidomains states (10, 15, 28). Second, we combine time-resolved optical and time-resolved x-ray diffraction (XRD) experiments to determine the complete light-induced ultrafast deformation of the unit cell in a $BiFeO_3$ thin film (Fig. 1). It enables us to access both in-plane and out-of-plane lattice distortions of $BiFeO_3$ films in the picosecond regime. Our results show that the ultrafast light excitation engenders a bipolar strain pulse with a leading compressive front, followed by a tensile tail for the longitudinal component. Concomitantly, positive and negative shear strain components are revealed in these leading front and tail, respectively. This temporal and spatial reconstruction of the longitudinal and shear strain has never been revealed before (9, 10, 12, 14–16, 20). For clarity, compressive (negative) and tensile (positive) components of the longitudinal (shear) strain are indicated by green and orange

Copyright © 2023 The American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).

<footnote n="1">¹Institut des Molécules et Matériaux du Mans, UMR 6283 CNRS, Le Mans Université, 72085 Le Mans, France. ²Synchrotron SOLEIL, L'Orme des Merisiers, Université Paris Saclay, 91190 Saint-Aubin, France. ³Université Paris-Saclay, CNRS UMR8502, Laboratoire de Physique des Solides, 91405, Orsay, France. ⁴Laboratoire de Physique de la Matière Condensée, UR2081, Université Jules Vernes Picardie, 80000 Amiens, France. ⁵Université Paris-Saclay, CentraleSupélec, CNRS-UMR8580, Laboratoire Structures, Propriétés et Modélisation des Solides, Gif-sur-Yvette, France. ⁶University of Arkansas, Physics Department, 825 W Dickson St., Fayetteville, AR 72701, USA. ⁷School of Materials Science and Engineering, UNSW Sydney, Kensington 2052, Australia. ⁸School of Physical and Chemical Sciences, University of Canterbury, Christchurch 8410 New Zealand. ⁹Unité Mixte de Physique CNRS, Thales, Université Paris-Saclay, Palaiseau 91767, France.</footnote>

<footnote>*Corresponding author. Email: vincent.juve@univ-lemans.fr (V.J.); daniel.sando@canterbury.ac.nz (D.S.); vincent.garcia@cnrs-thales.fr (V.G.); pascal.ruello@univ-lemans.fr (P.R.)</footnote>

<footnote>†These authors contributed equally to this work.</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Light-induced deformation of a ferroelectric lattice. The central sketch shows the (100)$_m$-oriented monoclinic BiFeO$_3$ film (blue) on the (110)$_c$-oriented SrTiO$_3$ substrate (gray) at equilibrium. The purple arrow denotes the ferroelectric polarization (P). $a_m, b_m$, and $c_m$ are the monoclinic axes of BiFeO$_3$ and $a, b$, and $c$ are the cubic axes of SrTiO$_3$. The femtosecond laser pulse induces an out-of-plane longitudinal strain (sketches I and III) and an in-plane shear strain (sketches II and IV) in the ferroelectric lattice. The longitudinal and shear strains can be either positive (orange arrows, sketches I and II) or negative (green arrows, sketches III and IV). Positive (negative) longitudinal and shear strain values correspond to an out-of-plane expansion (contraction) and an increase (decrease) of the monoclinic distortion, respectively. The corresponding variation of the ferroelectric polarization direction is shown in purple for each strain component. The 12-nm buffer SrRuO$_3$ layer is omitted for clarity.</figcaption></figure>

arrows in Fig. 1. The out-of-plane parameter reaches maximal compression and expansion of 0.1 and 0.6% respectively. Simultaneously, we observe a transient shear strain corresponding to the in-plane lattice distortion and of bipolar nature with a maximum amplitude of 0.04%. Our full temporal and spatial reconstruction of the transient lattice distortions points to successive reorientations of the ferroelectric polarization. Hence, ferroelectric polarization can be modulated on the picosecond time scale by a light-induced strain wave.

## RESULTS AND ANALYSIS

### Sample preparation and characterization

Figure 1 shows a sketch of our (100)$_m$-oriented [i.e., (110)$_{pc}$], 180-nm BiFeO$_3$ film grown on (110)$_c$ SrTiO$_3$ with a 12-nm SrRuO$_3$ electrode [m, monoclinic; c, cubic; pc, pseudocubic; see Methods and (29, 30) for further information on the growth process]. From XRD measurements, we identify the (100)$_m$ film orientation (Fig. 2A) and find the lattice parameters to be $a_m = 5.662$ Å, $b_m = 5.606$ Å, $c_m = 3.905$ Å, and monoclinic angle $\beta = 89.5^\circ$. The $b_m$ axis is partially relaxed, and the $c_m$ axis is fully strained (Fig. 2, B and C), as previously observed (26).

(100)$_m$-oriented BiFeO$_3$ films can show two different ferroelastic domain states, either with $c_m$ parallel or perpendicular to the [001]$_c$ axis of SrTiO$_3$ (26). For each ferroelastic domain state, two opposite polarization directions are possible, giving rise to four polar domain states in total. XRD reveals that the film is ~99.5% a pure single ferroelastic domain (see Methods and note S1). In addition, the polarization-dependent optical birefringence shows a sinusoidal shape with transmittance maxima for light polarizations parallel to $b_m$.

(Fig. 2D). As the SrTiO$_3$ substrate is cubic and hence optically isotropic, the sinusoidal shape shows the macroscopic optical anisotropy of our BiFeO$_3$ thin film and confirms the ferroelastic single-domain state.

To investigate the polar domain states within the ferroelastic domain, we used piezoresponse force microscopy (PFM) (Fig. 2, E to H; Methods; and note S1). The topography image exhibits uniaxial features parallel to the c axis of the substrate (Fig. 2E). The out-of-plane and in-plane PFM phase images are homogeneous (Fig. 2, F and G), which indicate a polar single-domain state. The out-of-plane component of the polarization points downward, and the in-plane component is antiparallel to [001]$_c$ SrTiO$_3$. We rotated the sample by 90° (Fig. 2H) and repeated the scan. The absence of an in-plane PFM signal confirms the orientation of the polar single-domain state. This ferroelectric polarization is switchable with a dc voltage of 7 V applied between the SrRuO$_3$ electrode and the PFM tip (insets of Fig. 2, F and G).

### Time-domain Brillouin light scattering

Next, we explore the impact of an ultrashort light pulse on our single-domain BiFeO$_3$ thin film. In our experiment, the light pulse impinges at normal incidence with a pump photon energy of 3 eV, which is larger than the BiFeO$_3$ bandgap (see Methods) (31). The absorption of a light pulse leads inevitably to the excitation of a strain pulse (32, 33) that propagates through the film perpendicular to its surface (11). The crystal symmetry of the BiFeO$_3$ film (point-group *m*) allows for the generation of both longitudinal and transverse acoustic waves, i.e., out-of-plane and in-plane atomic motions, respectively.
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Single-domain BiFeO₃ thin film on (110)ₐ SrTiO₃. (A) θ-2θ XRD pattern of the single-domain monoclinic (100)ₐ BiFeO₃ on the (110)ₐ SrTiO₃ substrate. (B and C) Reciprocal space maps near the 221ₐ (B) and 310ₐ (C) substrate Bragg reflections. (D) Optical transmission measurement as a function of the probe light polarization rotation angle performed with a 400-nm laser. (E to H) PFM images revealing the presence of a single ferroelectric domain. The topography (E), the out-of-plane PFM phase (F), and in-plane PFM phase in two different directions [cantilever perpendicular (G) and parallel to (H) [001]ₐ SrTiO₃, respectively] indicate that the polarization is pointing downward with an in-plane component antiparallel to [001]ₐ SrTiO₃. The insets in (F) and (G) show the 180° switched polarization after applying a dc voltage of 7 V to the bottom electrode. r.l.u., reciprocal lattice unit; A.U., arbitrary units.</figcaption></figure>

In a first step, we probe the change of the optical reflectivity in a pump-probe measurements with picosecond resolution. This allows us to reveal the different photo-induced acoustic modes through time-domain Brillouin light scattering (11, 18, 19, 34) (see Methods). The black curve in Fig. 3A shows the raw signal of the transient optical reflectivity as a function of time delay between the pump and the probe pulse. We subtract the electronic and thermal contributions of the raw signal and extracted the Brillouin oscillatory components of the time-domain optical reflectivity

(Fig. 3A in red). To identify the oscillatory components, we conducted a standard Wavelet analysis on the oscillatory signal. This yields the signal strength in the frequency domain as a function of the time delay between the pump and the probe pulses (Fig. 3B, bottom). The data reveal three distinct modes at about 31, 43, and 66 GHz. To pinpoint the origins of these modes, we compare the mode frequencies with literature data on bulk BiFeO₃ and SrTiO₃ (see also note S2). We find that the mode frequencies match with the transverse (shear strain) and longitudinal

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Time-resolved optical reflectivity. (A) Raw signal of the transient optical reflectivity (black curve) measured by time-resolved optical pump-probe experiments. The oscillatory coherent acoustic phonons signal is extracted from the raw signal (red curve). The purple line shows abrupt change in the oscillatory signal. The simulated oscillatory signal is shown by the orange dashed line. Inset: Zoom into the optical reflectivity for longer time delays. (B) Top: Comparison between the experimental acoustic signal (in red) and the simulated three acoustic modes (longitudinal, LA and shear, TA, acoustic modes in BiFeO₃ and one longitudinal acoustic mode in SrTiO₃, dashed lines) constituting the simulated acoustic signal in (A). Bottom: Wavelet analysis of the experimental acoustic signal. The colored map represents the experimental signal in the time-frequency domain and shows the three distinct acoustics modes at 31 GHz (TA mode in BiFeO₃), 43 GHz (LA mode in BiFeO₃), and at 66 GHz (LA mode in SrTiO₃). The vertical purple dashed line in (A) and (B) indicates the time delay when the LA wave in BiFeO₃ arrives at the BiFeO₃-substrate interface.</figcaption></figure>
Brillouin modes of bulk BiFeO₃ and with the longitudinal Brillouin mode of SrTiO₃, respectively. The deduced longitudinal ($v_L$ ≈ 4500 m/s) and shear ($v_T$ ≈ 3100 m/s) sound velocities in the BiFeO₃ thin film are also consistent with the values found in bulk BiFeO₃ (17, 18). Next, we investigate the signal evolution in time. Until ~40 ps, the experimental optical signal is dominated by the longitudinal acoustic phonon signal of the BiFeO₃ layer. At ~40 ps, the signal exhibits a discontinuity (dashed purple line in Fig. 3, A and B). This abrupt signal change corresponds to the arrival of coherent longitudinal acoustic phonons at the BiFeO₃-substrate interface (35–38). Thanks to a small acoustic reflection coefficient (see note S3), this acoustic mode is largely transferred into the substrate. After this acoustic transmission, the transverse acoustic mode dominates the experimental Brillouin oscillatory signal in the BiFeO₃ layer until roughly 60 ps. At this time, the transverse mode is partly transferred into the substrate and partly reflected at the BiFeO₃-substrate interface. Note that due to our measurement configuration, we can only probe the transmitted longitudinal acoustic wave in the cubic SrTiO₃ substrate (39).

To gain further insights, we have simulated the experimental strain signal in two steps. First, we model the lattice dynamics of

the BiFeO₃ thin film after the short laser pulse excitation by semi-coupled equations of motion. This allows the evaluation of the dynamics of coherent acoustic phonons in the structure (see note S3) (40). Second, the transient optical reflectivity is calculated on the basis of the standard photoelastic model (see Methods and note S3). In this simulation, because we do not know the photoelastic coefficients (for the longitudinal and the shear waves) involved in the Brillouin light scattering detection process in BiFeO₃ and SrTiO₃, we used them as free parameters to adjust the amplitude of the longitudinal and shear waves. The simulation results are displayed as dashed orange lines in Fig. 3A. The simulations are in excellent agreement with the experimental results even at longer time scales (inset of Fig. 3A). In Fig. 3B (top), we displayed the different waves that compose the total simulated signal as dashed lines. Hence, we can identify the different time windows where one specific acoustic mode dominates the signal.

These preliminary all optical experiments demonstrate the possibility to light-induce longitudinal and transversal strain waves in our single ferroelectric domain BiFeO₃ thin film. Next, we build up a detailed understanding of the evolution of the strain waves.

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Time-resolved XRD experiments. (A) Sketch of the time-resolved XRD setup. Rocking curves (ω scans) were measured to record the evolution of the Bragg peak positions as a function of the time delay (Δt) between the optical pump pulse and the x-ray probe pulse. The SrRuO₃ buffer layer is omitted for clarity. (B) Transient evolution of the 403m Bragg peak (red dots) after ultrashort laser pulse excitation compared to the reference unexcited Bragg peak (blue dots). The lines are the results of the fitting procedure (Methods).</figcaption></figure>
Time-resolved XRD

The time-domain Brillouin light scattering can probe acoustic phonons through the photoelastic interaction but cannot provide directly the amplitude of the photo-induced strain (41). In the next step, we quantitatively evaluate the strain amplitude within the unit cells. To do so, we turn to time-resolved XRD. This technique allows us to monitor the photo-induced strain pulse through the evolution of the Bragg peak positions with time (9–12, 42). First, we have to identify suitable Bragg reflections. Here, a look on symmetry properties is insightful. With the point-group symmetry (m) of BiFeO₃, the mirror plane lies perpendicular to the bₘ axis and contains the aₘ and cₘ axes. Hence, the photoexcitation cannot lead to in-plane atomic displacements along bₘ due to symmetry principles (see note S4). In other words, the photo-induced stress must have two components (11, 18, 33, 43): an out-of-plane and in-plane components. These two stress components cause then atoms to move simultaneously along out-of-plane and in-plane directions (11, 18, 33, 43). Accordingly, photo-induced deformations comprise two strain components: the longitudinal strain η_L, which describes the expansion or contraction along aₘ, and the shear strain η_T with atomic displacement along cₘ (Fig. 1). Diffraction measurements with Bragg planes (h0l)ₘ and (h0l)ₘ type gives access simultaneously to the longitudinal and shear strain (11). In contrast, the (hk0)ₘ/(hκ0)ₘ planes are only sensitive to the longitudinal strain (see note S4). From an experimental point of view, we monitored the photo-induced transient evolution of different Bragg peak positions by so-called ω scans (i.e., rocking curves; see Fig. 4A, Methods, and note S4). In this geometry, the η_L and η_T strains

translate into diffraction angles ω as follows

$$ (\omega' - \omega_0)_{h0l} = -A(\theta, h, |l|) \left[ \eta_L \frac{h}{|l|} + \sqrt{2}\eta_T \right] \\ (\omega' - \omega_0)_{h0\bar{l}} = -A(\theta, h, |l|) \left[ \eta_L \frac{h}{|l|} - \sqrt{2}\eta_T \right] \\ (\omega' - \omega_0)_{h\bar{k}0} = (\omega' - \omega_0)_{h\bar{k}0} = -B(\theta, h, |k|)\eta_L \quad (1) $$

where ω' and ω₀ denote the time-dependent diffraction angle and the diffraction angle at equilibrium, respectively. A(θ, h, |l|) and B(θ, h, |k|) are positive crystallographic parameters. Details about the derivation of Eq. 1 can be found in the note S4. As we see in Eq. 1, measuring ω' - ω₀, with these set of Bragg peaks, allows us to reconstruct fully the light-induced strain evolution of η_L(t) and η_T(t) at the picosecond time scale.

Equipped with these preliminary considerations, we embark to investigate the photo-induced strain in our BiFeO₃ thin film. We select the Bragg reflections of the (403)ₘ and (403)ₘ planes to access the longitudinal and shear strain distortions. In addition, we measure the Bragg reflections of the (530)ₘ and (530)ₘ planes as a reference for distortion of the longitudinal strain only. Figure 4B shows representative ω scans of the 403ₘ Bragg peak for selected time delays between ~50 and 190 ps (red curves). The blue curve is the reference Bragg reflection of the unperturbed system. The transient evolution of the 403ₘ Bragg peak profile reveals a complex behavior. First, within the first ~30 ps, the Bragg peak is highly asymmetric and is split into two contributions. Second, at longer time scale, only one Bragg peak exists and remains shifted to a lower angle up to the maximum measured time delay (~200

<figure><img src="image_6.png" /><figcaption>Fig. 5. Time-resolved light-induced strain. (A and B) Maps of the Bragg peak position change (ω' - ω₀) recorded as a function of the time delay for 530m and 530m reflections (A) and 403m and 403m reflections (B). The right panels show the temporal evolution of the Bragg peaks positions change for the tensile and the compressive parts (inset) for the 530m/530m (top) and 403m/403m reflections (bottom).</figcaption></figure>
ps). This general trends is also shared by the 403<sub>m</sub>, 530<sub>m</sub>, and 530<sub>m</sub> Bragg reflections.

For an in-depth analysis of the evolution of Bragg reflections, we have extracted the diffraction angles by a standard fitting procedure (Methods). Depending on the time delay, we have two or one Bragg peak to fit. Figure 5 (A and B) shows the $\omega' - \omega_0$ evolution for the four reflections as a function of the time delay. At first glance, the diffraction patterns of all four planes show similar behaviors. Past the time delay zero, the Bragg peak splits into two distinct Bragg reflections (later called negatively and positively shifted Bragg peaks), as shown in Figs. 4B and 5 (A and B). These two Bragg peaks indicate that, in the early moments, the BiFeO<sub>3</sub> layer consists of two distinct regions: one under compressive strain [($\omega' - \omega_0$) > 0] and the other under tensile strain [($\omega' - \omega_0$) < 0], as we will describe later. One Bragg reflection [($\omega' - \omega_0$) > 0] shows a steady shift from $\omega_0$ to $\omega' - \omega_0 \approx 0.2^\circ$ and disappears at ~30 ps. In contrast, the second Bragg reflection [($\omega' - \omega_0$) < 0] appears with the time delay zero at $\omega' - \omega_0 \approx -0.4^\circ$. The difference in $\omega' - \omega_0$ decreases over the first ~30 ps to $\omega' - \omega_0 \approx -0.15^\circ$, where it stabilizes for the maximal time delay of our measurements of ~200 ps. A closer look confirms the nearly identical temporal evolutions of the (530)<sub>m</sub> and (530)<sub>m</sub> Bragg reflections (Fig. 5A, right; note the apparent difference the (530)<sub>m</sub> and (530)<sub>m</sub> reflection around time delay zero results from the sharp signal change). Yet, the 403<sub>m</sub> and 403<sub>m</sub> Bragg peak shifts exhibit distinct changes in their temporal evolutions (Fig. 5B). In the time range where two Bragg peaks are detected (0 to 30 ps), whatever the peak we consider [($\omega' - \omega_0$) > 0 or ($\omega' - \omega_0$) < 0], the shift of the 403<sub>m</sub> Bragg peak is always more pronounced (positively or negatively) than the shift of the 403<sub>m</sub> Bragg peak. For time delays of 40 to 120 ps, where we detect only one Bragg peak, being negatively shifted regarding to $\omega_0$, this behavior is reversed and the 403<sub>m</sub> Bragg peak shift becomes larger than that of the 403<sub>m</sub>. Last, the diffraction angles of the 403<sub>m</sub> and 403<sub>m</sub> planes converge toward 200 ps. Crucially, the different behavior observed between 403<sub>m</sub> and 403<sub>m</sub> in Fig. 5B (right) proves that shear strain contributes to the Bragg shift, consistently with the model presented in Eq. 1.

To extract the longitudinal and shear strain amplitudes quantitatively, we rewrite Eq. 1, expressing the longitudinal strain $\eta_L$ and shear strain $\eta_T$ in terms of a combination of the ultrafast light-induced changes of the diffraction patterns

$$ \eta_L = \frac{-|l|}{2hA(\theta, h, |l|)} [(\omega' - \omega_0)_{h0l} + (\omega' - \omega_0)_{h0\bar{l}}] = \frac{-1}{B(\theta, h, |k|)} (\omega' - \omega_0)_{h\pm k0} \\ \eta_T = \frac{-1}{2\sqrt{2}A(\theta, h, |l|)} [(\omega' - \omega_0)_{h0l} - (\omega' - \omega_0)_{h0\bar{l}}] \quad (2) $$

We use the experimental data ($\omega' - \omega_0$)<sub>403<sub>m</sub></sub> and ($\omega' - \omega_0$)<sub>403<sub>m</sub></sub> shown in Fig. 5B to calculate the longitudinal and shear strains, which are displayed in Fig. 6 (A and B, respectively). We find that tensile (label 1 in Fig. 6A) and compressive (label 2 in Fig. 6A) longitudinal strain components coexist for the first 10 to 30 ps, consistently with the coexistence of positively and negatively shifted Bragg peaks. We also note that the compressive and tensile longitudinal strains deduced from the $h\bar{k}_0<sub>m</sub>$ or from the combination of $h0l<sub>m</sub>$ and $h0l\bar{\bar{l}}<sub>m</sub>$ reflections have the same overall time dependence, which directly demonstrates that the extraction of the longitudinal

strain does not depend on the selected Bragg peak, thus confirming the validity of the model presented in Eqs. 1 and 2. From the negatively shifted Bragg peak, we determine a maximal longitudinal strain of 0.6%. In the same time range, we also determine the shear strain amplitude, which we can extract from both the negatively (label 1 in Fig. 6B) and positively shifted (label 2 in Fig. 6B) Bragg peaks. We obtain 0.04% as the largest amplitude of the shear strain. The coexistence of positive and negative longitudinal and shear strains can be explained and simulated by the standard elastic response of a solid object to a photo-induced stress (see Fig. 6, C and D), which has been previously introduced to simulate the time-domain Brillouin light scattering experiment (Fig. 3 and note S3). When the light is absorbed in the bulk of the BiFeO<sub>3</sub> layer near the surface, the light-induced stress generates a lattice strain pulse that propagates in both forward and backward directions (32). When the backward emitted strain pulse arrives at the mechanically free air-BiFeO<sub>3</sub> interface, it is reflected and the sign of the strain reverses. As a result, the propagating strain pulse becomes bipolar (note S3). For example, if the initial longitudinal stress is negative (as in the case of thermal expansion), then the front of the strain pulse is compressive and its tail is tensile. Figure 6C illustrates this scenario with the simulation of the spatial dependence of the longitudinal bipolar strain pulse as it travels through the BiFeO<sub>3</sub>/SrRuO<sub>3</sub>/SrTiO<sub>3</sub> system 30 ps after photoexcitation. In our experiment, the x-rays penetrate and probe the entire BiFeO<sub>3</sub> layer at once, so that they can detect the regions exhibiting different strain fields simultaneously. This explains the coexistence of a compressive and tensile strain in the time-domain measurements shown in Figs. 6A-B, i.e., why we have a splitting of the Bragg peak in the early time-delays. Similar light-induced bipolar longitudinal strain waves have also been observed in other materials, such as InSb (42) and SrRuO<sub>3</sub> (44). Fig. 6D illustrates the relationship between the negatively and positively shifted Bragg peaks and the positive and negative longitudinal strain present in the tail and front of the strain pulse, respectively. As previously mentioned, light also induces shear stress, which generates a bipolar shear strain pulse within the film (see note S3). This additional contribution of the shear strain components to the Bragg peak shift has a different sign in the front and the tail of the strain pulse and has a different impact on the 403<sub>m</sub> and 403<sub>m</sub> peaks as shown in Fig. 5B and summarized in Fig. 6D. This is in perfect agreement with our model presented in Eq. 1.

These results are worth emphasizing for two reasons. First, to the best of our knowledge, the longitudinal strain level achieved in this study is among the largest photo-induced longitudinal strain levels ever reported in BiFeO<sub>3</sub> film. Similarly, high strain levels have been reported recently but only in multidomain BiFeO<sub>3</sub> film (15). This efficient photo-induced strain is likely connected to efficient generation mechanisms of acoustic waves that combine light-induced thermal (ultrafast thermal expansion) and nonthermal (deformation potential mechanism, inverse piezoelectric effect) processes, as already discussed in the literature for the rhombohedral BiFeO<sub>3</sub> (9–11, 15, 18). Second, it is the first observation and evaluation of ultrafast light-induced shear strain component in ferroelectric thin films. Previous works, for example, on BiFeO<sub>3</sub> and PbTiO<sub>3</sub> thin films, targeted the out-of-plane motion only (9, 10, 12, 14, 15, 20).
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Light-induced longitudinal and shear strains. (A and B) Time evolution of the longitudinal [$\eta_L$ (A)] and shear [$\eta_T$ (B)] strains extracted from the 403$_m$/40$\bar{3}_m$ and 530$_m$/5$\bar{3}_m$ Bragg reflections. The color filling [orange and green in (A) and (B)] and the small sketches in (A) and (B) indicate the unit cell distortions. The purple dashed line shows the approximate arrival time of the out-of-plane strain pulse to the BiFeO$_3$-substrate interface. (C) Simulation of the spatial dependence of the longitudinal strain in the sample at 30 ps after ultrashort laser pulse excitation. The profile shows the compressive front (green) and the tensile tail (orange) of the strain pulse propagating at the speed of sound. These opposite strain components lead to the two diffraction peaks as shown in (D) and Fig. 5B. (D) Left: Sketch of the light-induced strain components (longitudinal-L and shear-T) within the BiFeO$_3$ film illustrating the compressive front and tensile tail of the L and T strain pulses. Right: The 403$_m$ and 40$\bar{3}_m$ Bragg peaks at 10-ps time delay. The different shifts of the two Bragg peaks are related to the L and T strain contributions. Note that the shear strain component has an opposite contribution for 403$_m$ and 40$\bar{3}_m$ consistently with Eq. 1. (E) Temporal and spatial dependence of the light-induced strain in BiFeO$_3$ thin film and its effect on the ferroelectric polarization. The longitudinal and shear strains travel at velocities $v_L$ and $v_T$, respectively, distort the unit cell and the local unit-cell ferroelectric polarization direction (small purple arrows) as a function of the time delay. The uniform macroscopic ferroelectric polarization before laser excitation is shown by the large purple arrow (sketch I).</figcaption></figure>
In the next step, we can have a look at the details of the evolution of the longitudinal and shear strains as they propagate through the sample. Let us start with the longitudinal strain pulse (Fig. 6A). The compressive strain (label 2 for the front) increases as the strain pulse propagates toward the BiFeO₃-substrate interface (see Fig. 6A and fig. S5) and reaches a maximum strain level at around 20 ps, where it vanishes. The tensile longitudinal strain region (label 1 for the tail) emerges as a steep increase and reaches its maximum amplitude within the first 10 ps. Subsequently, it decays and reaches a plateau at ~40 ps (denoted by the purple dashed line in Fig. 6A). Because of the small acoustic reflection coefficient at the interface, the longitudinal strain wave is transmitted almost entirely into the substrate after ~40 ps ($R_L \approx 2\%$; note S3). The drop of the longitudinal strain at ~40 ps is perfectly consistent with the drop in the Brillouin oscillatory signal of the longitudinal mode in BiFeO₃ in Fig. 3. This efficient strain transmission also explains why the compressive component (front) disappears and why only one Bragg peak remains after ~40 ps. The constant plateau-like longitudinal strain state after ~40 ps (labels 3 and 4 in Fig. 6A) results then from an overall expansion of the thin film caused by long-living and nonpropagating (at this time scale) thermal and/or electronic contributions (9, 11, 18). Next, we investigate the evolution of the shear strain (Fig. 6B). We remind that the shear strain is extracted by combining the nonequivalent responses of the 403 m/403 m reflections (Figs. 5B and 6D) as explained in the second term of Eq. 2. The maximal amplitude of the shear strain is an order of magnitude smaller than that of the longitudinal strain. Directly after the time delay zero, the front of the shear strain pulse (label 2) shows a positive amplitude, i.e., the monoclinic distortion increases. This shear-pulse front vanishes at ~20 ps. On the contrary, the tail of the shear strain pulse (label 1) is negative from 0 to ~50 ps. Accordingly, the monoclinic distortion decreases. At ~50 ps, this propagating shear strain changes sign and the monoclinic distortion increases again. It remains positive up to the time delay of ~130 ps, where it vanishes (labels 3 and 4). The absence of the shear strain signal beyond 130 ps, i.e., the absence of a plateau-like behavior at long time scale, indicates a rather small contribution to the shear strain of the thermal and/or the long-living electronic processes. This positive shear strain from 50 to 130 ps may have two origins. (i) The reflection coefficient of the shear acoustic wave at the BiFeO₃-substrate interface has a non-negligible negative value ($R_T \approx -10\%$; note S4). Such a negative reflection coefficient can lead to a reverse of the strain sign when the tail of the shear strain impinges on the substrate. (ii) As BiFeO₃ is elastically anisotropic, a mode conversion of the incident longitudinal wave into a shear wave may also occur (43). With a longitudinal strain 10 times larger than the shear strain, even a few percent of efficiency of this mode conversion could give rise or contribute to the shear strain between 50 and 130 ps.

## DISCUSSION

Figure 6E provides a sketch of the temporal and spatial tracking of the longitudinal and shear strain fields. To help the reader, we label the sketches I to V in Fig. 6E. We link deformations of the film in Fig. 6E to the strain states in Fig 6 (A and B) with the labels 1 to 4. Sketch I represents the unperturbed monoclinic BiFeO₃ film before the arrival of the light pulse. In sketch II, the light pulse has excited the strain pulse with longitudinal and shear components traveling through the BiFeO₃ film at their respective velocities $v_L$ and $v_T$. In

the tail (label 1), the combination of a positive longitudinal and negative shear strains leads to an expansion of the unit cell along $a_m$ together with a decrease of its monoclinicity. In the front (label 2), the longitudinal and shear strain have opposite signs. This leads then to a contraction of the unit cell along $a_m$ together with an increase of its monoclinicity. As time goes on (sketch III), the strain pulse travels further through the BiFeO₃ film and the longitudinal and shear components arrive at the interface at around 40 and 50 ps, respectively. In sketch IV, the longitudinal strain pulse has been transferred to the substrate. Only the long-living longitudinal strain component remains. Hence, the unit cell remains expanded along $a_m$ (label 3). At this moment (~50 ps, label 3), the shear strain is nearly zero. Past this moment, the shear strain experiences a sign change and remains positive up to ~130 ps along with the long-living expansion along $a_m$ (label 4 and sketch V).

In BiFeO₃ films, the electric polarization is locked to the lattice, which is a fundamental property for ferroelectrics. This so-called adiabatic approximation is relevant even at the picosecond time scale as supported by theoretical calculations (45) and also recently discussed in ferroelectric superlattices submitted to a picosecond strain pulse (28). Assuming that the ferroelectric polarization equally follows the strain wave, as sketched in Fig. 6E (purple arrows), this entails an electric charge wave propagating at the speed of sound through the BiFeO₃ film. The charge wave may even include electrostatic discontinuities because of the bipolar nature of the propagating strain wave. Considering the fast dynamics and the fundamental relation between a transient polarization and the current ($\vec{J} = \frac{\partial \vec{P}}{\partial t}$), we can expect ultrafast strain-induced current approaching the terahertz range. This is a promising property for next-generation ferroelectric-based gigahertz-terahertz technologies. The full description of this charged propagating acoustic front requires additional experiments and is planned future work.

In summary, we have used a combination of time-resolved XRD and optical ultrafast spectroscopy experiments to unravel the optically induced picosecond strain dynamics of a single-domain BiFeO₃ thin film. We report maximum longitudinal and shear strain of 0.6 and 0.04%, respectively. These strain levels are comparable to those of commercially available electrically driven transducers. However, the important point is that the strain pulse here has a typical duration of tenths of picoseconds, corresponding to a typical bandwidth of 100 GHz, while electrically driven piezotransducers typically operate in the megahertz regime. This demonstrates that efficient conversion of light energy into mechanical energy is also possible at the picosecond time scale. We have also demonstrated the complex temporal and spatial profile of both the out-of-plane and in-plane strain components. As a result, the ultrafast light-induced strain dynamics modulate the ferroelectric polarization direction in time and space within the film. This work constitutes an important milestone toward the exploration of the rich lattice dynamics in complex ferroelectric textures, ranging from mosaic to self-organized periodic structures (46–50). Moreover, ultrafast light-induced anisotropic strains might be used to create new transient phases or to generate ultrashort current pulses based on strain-induced modulation of the polarization direction. In the quest for the ultrafast manipulation of the polar and magnetic orders in multiferroic compounds such as BiFeO₃, future work will require considering the dynamics of the unit-cell strain to account for
comprehensive (inverse) piezoelectric and/or (inverse) magnetostRICTIVE effects at the picosecond time scale.

METHODS

Epitaxial thin-film growth

A single-domain epitaxial BiFeO₃ thin film was grown by pulsed laser deposition on a SrRuO₃-coated (110)ₚc SrTiO₃ substrate using conditions reported previously (29, 30). Ablation was carried out by a KrF laser with wavelength λ = 248 nm at 5 Hz. The 12-nm layer of SrRuO₃ was grown from a stoichiometric target at 660°C in 100 mtorr of oxygen with fluence of ~2 J/cm². The BiFeO₃ layer was grown from a Bi₁.₁FeO₃ ceramic target, while the substrate was held at 590°C in a background oxygen pressure of 100 mtorr. The thickness of the BiFeO₃ layer is estimated at 180 ± 10 nm, extrapolated from the growth rate, taken from a sample grown directly beforehand on which x-ray reflectivity calibration was performed.

XRD structural characterization

The film was characterized using XRD using an Xpert Pro Materials Research Diffractometer (MRD) using Cu Kα₁ radiation from a two bounce Ge (220) monochromator (λ = 1.6504 Å) and one-dimensional (1D) detector (PIXcel). Single-phase growth was confirmed by standard coupled θ-2θ scans, while the domain structure and lattice parameters were determined by reciprocal space maps (RSMs) near the 220_c, 221_c, and 310_c SrTiO₃ reflections (fig. S1). In the (110)_c SrTiO₃ orientation, only two ferroelastic domains are possible (with four total ferroelectric directions) (26). The relative domain population of ferroelastic domains was determined by fitting 2D Gaussians to the BiFeO₃ peaks in the RSM around the 221_c SrTiO₃ reflection (fig. S2). The film was thus determined to have a 99.5%/0.5% volume fraction, making, it in a practical sense, a single domain.

Piezoresponse force microscopy

PFM experiments were performed in a Nanoscope V multimode (Bruker) using an external ac source (DS360, Stanford Research) to excite the BiFeO₃ thin film at a frequency of 35 kHz (far off resonance) with typical ac voltage excitations of 2 V peak to peak and external lock-in amplifiers (SR830, Standford Research). We used Cr/Pt-coated tips with a cantilever stiffness of 40 N/m. Two different orientations of the cantilever were used for the in-plane PFM measurements (Fig. 2, G and H) to characterize the in-plane domain structure of the films (fig. S3). The single ferroelectric domain could be switched with +7 V applied to the SrRuO₃ electrode while scanning with the slow scan axis along [001]c, resulting in another single domain within the same ferroelastic domain (switching by 180°) (Fig. 2).

Optical birefringence measurement

The optical birefringence was measured in transmission geometry with a 400-nm laser [photon energy of 3.1 eV, above the bandgap of BiFeO₃ to ensure absorption (31, 51–53)] with an unfocused beam with a diameter of ~2 mm.

Time-domain Brillouin light scattering

Time-domain Brillouin light scattering was conducted at IMMM with a Ti:sapphire femtosecond laser oscillator (11, 17–19, 34).

The oscillator delivered a 120-fs beam at 830 nm, which was split into two beams, the first (pump) being doubled in frequency (λ<sub>pump</sub> = 415 nm) in a BaB<sub>2</sub>O<sub>4</sub> crystal. The second beam (probe) is used to synchronously pump an optical parametric oscillator delivering a femtosecond laser beam in the visible range (λ<sub>probe</sub> = 587 nm) below the bandgap of BiFeO<sub>3</sub>, which permits the probe beam to penetrate in the film and in the substrate and hence to detect the propagating strain pulse. The transient optical reflectivity was measured through a balanced photodiode scheme. Under this condition, the optical detection process reveals only Brillouin components of the light-induced strain field (17–19, 34). These Brillouin oscillation frequencies are given by f<sub>B</sub> = 2 nv<sub>s</sub>λ<sub>probe</sub><sup>-1</sup> (for a probe beam with a normal incidence), where n is the refractive index of the probed medium and v<sub>s</sub> is the sound velocity (different for longitudinal or transverse acoustic modes).

Time-resolved XRD

The time-resolved XRD experiment was performed at the CRISTAL beamline of SOLEIL synchrotron. Optical excitation was achieved with 400-nm laser pulses of <80-fs duration as obtained by frequency doubling the fundamental output of the laser system (Micra oscillator and Legend Elite Duo HP amplifier, Coherent). The photoinduced lattice deformations were probed using x-rays of 7.082-keV energy, selected by means of a Si(111) double-crystal monochromator. The overall time-resolution of the experiment is set by the x-ray pulse duration, which is 10-ps full width at half maximum, as the storage ring was operated in the low-α mode (54). The pump beam irradiates nearly entirely the sample surface, while the probe x-ray beam is made at least twice smaller (maximum size of the x-ray footprint on sample: 90 µm by 500 µm). The various Bragg reflections studied were imaged using a photon-counting camera that was newly developed at SOLEIL, which enables simultaneous measurements at positive and negative pump-probe delays (55). Bragg reflections with high Miller indices (such as 403<sub>m</sub> and 530<sub>m</sub>) were chosen to increase the sensitivity of the measurement to the light-induced strain. The rocking curves were fitted with either one or two pseudo-Voigt functions.

Light-induced strain calculations

Light-induced strain calculations were done using the udkm1Dsim toolbox in Python, which simulates the thermal and the structural dynamics after light excitation of a 1D layered structure (40). The complete thermal and structural dynamics of the BiFeO3/SrRuO3/SrTiO3 nanostructure were modeled, although the strain of SrRuO3 layer is not displayed in the manuscript, especially in Fig. 3, as it is rather small. Further details can be found in note S3.

**Supplementary Materials**

This PDF file includes:

Figs. S1 to S9

Tables S1 and S2

Notes S1 to S5

References

REFERENCES AND NOTES

1. M. E. Lines, A. M. Glass, Principles and Applications of Ferroelectrics and Related Materials (1977).
2. M. Fiebig, T. Lottermoser, D. Meier, M. Trassin, The evolution of multiferroics. *Nat. Rev. Mater.* **1**, 16046 (2016).

3. N. Spaldin, R. Ramesh, Advances in magnetoelectric multiferroics. *Nat. Mater.* **18**, 203–212 (2019).

4. Y. Zhang, J. Dai, X. Zhong, D. Zhang, G. Zhong, J. Li, Probing ultrafast dynamics of ferro- electrics by time-resolved pump-probe spectroscopy. *Adv. Sci.* **8**, e2102488 (2021).

5. T. Linker, K.-I. Nomura, A. Aditya, S. Fukushima, R. K. Kalia, A. Krishnamoorthy, A. Nakano, P. Rajak, K. Shimmura, F. Shimojo, P. Vashishta, Exploring far-from-equilibrium ultrafast polarization control in ferroelectric oxides with excited-state neural network quantum molecular dynamics. *Science.* **8**, eabk2625 (2022).

6. T. F. Nova, A. S. Disa, M. Fechner, A. Cavalleri, Metastable ferroelectricity in optically strained $SrTiO_3$. *Science* **364**, 1075–1079 (2019).

7. X. Li, T. Qiu, J. Zhang, E. Baldini, J. Lu, A. M. Rappe, K. A. Nelson, Terahertz field induced ferroelectricity in quantum paraelectric $SrTiO_3$. *Science* **364**, 1079–1082 (2019).

8. R. Mankowsky, A. von Hoegen, M. First, A. Cavalleri, Ultrafast reversal of the ferroelectric polarization. *Phys. Rev. Lett.* **118**, 197601 (2017).

9. H. Wen, P. Chen, M. P. Cosgriff, D. A. Walko, J. H. Lee, C. Adamo, R. D. Schaller, J. F. Ihlefeld, E. M. Dufresne, D. G. Schlom, P. G. Evans, J. W. Freeland, Y. Li, Electronic origin of ultrafast photoinduced strain in $BiFeO_3$. *Phys. Rev. Lett.* **110**, 037601 (2013).

10. D. Schick, M. Herzog, H. Wen, P. Chen, C. Adamo, P. Gaal, D. G. Schlom, P. G. Evans, Y. Li, M. Bargheer, Localized excited charge carriers generate ultrafast inhomogeneous strain in the multiferroic $BiFeO_3$. *Phys. Rev. Lett.* **112**, 097602 (2014).

11. V. Juvé, R. Gu, S. Gable, T. Maroutian, G. Vaudel, S. Matzen, N. Chigarev, S. Raetz, V. E. Gusev, M. Viret, A. Jarnac, C. Laulhé, A. A. Maznev, B. Dkhil, P. Ruello, Ultrafast light-induced shear strain probed by time-resolved x-ray diffraction: Multiferroic $BiFeO_3$ as a case study. *Phys. Rev. B* **102**, 220303 (2020).

12. D. Daranciang, M. J. Highland, H. Wen, S. M. Young, N. C. Brandt, H. Y. Hwang, M. Vattilana, M. Nicoul, F. Quirin, J. Goodfellow, T. Qi, I. Grinberg, D. M. Fritz, M. Cammarata, D. Zhu, H. T. Lemke, D. A. Walko, E. M. Dufresne, Y. Li, J. Larsson, D. A. Reis, K. Sokolowski-Tinten, K. A. Nelson, A. M. Rappe, P. H. Fuoss, G. B. Stephenson, A. M. Lindenberg, Ultrafast photo-voltaic response in ferroelectric nanolayers. *Phys. Rev. Lett.* **108**, 087601 (2012).

13. C. Paillard, B. Xu, B. Dkhil, G. Geneste, L. Bellaiche, Photostriction in ferroelectrics from density functional theory. *Phys. Rev. Lett.* **116**, 247401 (2016).

14. S. Matzen, L. Guillemot, T. Maroutian, S. K. K. Patel, H. Wen, A. D. Dihiara, G. Agnus, O. G. Shpyrko, E. E. Fullerton, D. Ravelosona, P. Lecoeur, R. Kukreja, Tuning ultrafast photoinduced strain in Ferroelectric-Based devices. *Adv. Electron. Mater.* **5**, 1800709 (2019).

15. H. J. Lee, Y. Ahn, S. D. Marks, D. S. Gyan, E. C. Landahl, J. Y. Lee, T. Y. Kim, S. Unithrattil, S. H. Chun, S. Kim, S.-Y. Park, C. I. Eom, A. H. Wen, D. G. Schlom, S. Lee, J. Y. Jo, P. G. Evans, Subpicosecond optical stress generation in multiferroic $bifeo_3$. *Nano Lett.* **22**, 4294–4300 (2022).

16. C. K. Schmising, M. Bargheer, M. Kiel, N. Zhavoronkov, M. Woerner, T. Elsaesser, I. Vrejoiu, D. Hesse, M. Alexe, Coupled ultrafast lattice and polarization dynamics in ferroelectric nanolayers. *Phys. Rev. Lett.* **98**, 257601 (2007).

17. P. Ruello, T. Pézeril, S. Avanesyan, G. Vaudel, V. Gusev, I. C. Infante, B. Dkhil, Photoexcitation of gigahertz longitudinal and shear acoustic waves in $BiFeO_3$ multiferroic single crystal. *Appl. Phys. Lett.* **100**, 212906 (2012).

18. M. Lejman, G. Vaudel, I. C. Infante, P. Gemeiner, V. E. Gusev, B. Dkhil, P. Ruello, Giant ultrafast photo-induced shear strain in ferroelectric $BiFeO_3$. *Nat. Commun.* **5**, 4301 (2014).

19. M. Lejman, G. Vaudel, I. C. Infante, I. Chaban, T. Pézeril, M. Edely, G. F. Nataf, M. Guennou, J. Kreisel, V. E. Gusev, B. Dkhil, P. Ruello, Ultrafast acousto-optic mode conversion in optically birefringent ferroelectrics. *Nat. Commun.* **7**, 12345 (2016).

20. L. Y. Chen, J. C. Yang, C. W. Luo, C. W. Laing, K. H. Wu, J.-Y. Lin, T. M. Uen, J. Y. Juang, Y. H. Chu, T. Kobayashi, Ultrafast photoinduced mechanical strain in epitaxial $BiFeO_3$ thin films. *Appl. Phys. Lett.* **101**, 041902 (2012).

21. P. Khan, M. Kanamaru, K. Matsumoto, T. Ito, T. Satoh, Ultrafast light-driven simultaneous excitation of coherent terahertz magnons and phonons in multiferroic $BiFeO_3$. *Phys. Rev. B* **101**, 134413 (2020).

22. P. Chen, C. Paillard, H. J. Zhao, J. Íñiguez, L. Bellaiche, Deterministic control of ferroelectric polarization by ultrafast laser pulses. *Nat. Commun.* **13**, 1–8 (2022).

23. A. V. Kimel, A. M. Kalashnikova, A. Pogrebna, A. K. Zvezdin, Fundamentals and perspectives of ultrafast photoferroic recording. *Phys. Rep.* **852**, 1–46 (2020).

24. J.-Y. Chaleau, T. Chirac, S. Fusil, V. Garcia, W. Akhtar, J. Tranchida, P. Thibaudeau, I. Gross, C. Blouzon, A. Finco, M. Bibes, B. Dkhil, D. D. Khatyavin, P. Manuel, V. Jacques, N. Jaouen, M. Viret, Electric and antiferromagnetic chiral textures at multiferroic domain walls. *Nat. Mater.* **19**, 386–390 (2020).

25. A. Haykal, J. Fischer, W. Akhtar, J.-Y. Chaleau, D. Sando, A. Finco, F. Godel, Y. A. Birkholzer, C. Carretero, N. Jaouen, M. Bibes, M. Viret, S. Fusil, V. Jacques, V. Garcia, Antiferromagnetic textures in $BiFeO_3$ controlled by strain and electric field. *Nat. Commun.* **11**, 1704 (2020).

26. D. Sando, F. Appert, B. Xu, O. Paull, S. R. Burns, C. Carretero, B. Dupe, V. Garcia, Y. Gallais, A. Sacuto, M. Cazayous, B. Dkhil, J. M. Le Breton, A. Barthelemy, M. Bibes, L. Bellaiche, V. Nagarajan, J. Juraszek, A magnetic phase diagram for nanoscale epitaxial $BiFeO_3$ films. *Appl. Phys. Rev.* **6**, 041404 (2019).

27. D. Sando, Strain and orientation engineering in $ABO_3$ perovskite oxide thin films. *J. Phys. Condens. Matter* **34**, 153001 (2022).

28. H. J. Lee, Y. Ahn, S. D. Marks, E. C. Landahl, S. Zhuang, M. H. Yusuf, M. Dawber, J. Y. Lee, T. Y. Kim, S. Unithrattil, S. H. Chun, S. Kim, I. Eom, S.-Y. Park, K. S. Kim, S. Lee, J. Y. Jo, J. Hu, P. G. Evans, Structural evidence for ultrafast polarization rotation in ferroelectric/dielectric superlattice nanodomains. *Phys. Rev. X* **11**, 013031 (2021).

29. Y. Zhang, Y. Tan, D. Sando, L.-Q. Chen, N. Valanoor, Y. Zhu, M.-G. Han, Controlled nucleation and stabilization of ferroelectric domain wall patterns in epitaxial (110) bismuth ferrite heterostructures. *Adv. Funct. Mater.* **30**, 2003571 (2020).

30. D. Sando, T. Young, R. Bulanadi, X. Cheng, Y. Zhou, M. Weyland, P. Munroe, V. Nagarajan, Designer defect stabilization of the super tetragonal phase in >70-nm-thick $BiFeO_3$ films on $LaAlO_3$ substrates. *Jap. J. Appl. Phys.* **57**, 0902B2 (2018).

31. S. G. Choi, H. T. Yi, S.-W. Cheong, J. N. Hilfiker, R. France, A. G. Norman, Optical anisotropy and charge-transfer transition energies in $BiFeO_3$from 1.0 to 5.5 eV. *Phys. Rev. B* **83**, 100101 (2011).

32. C. Thomsen, H. T. Grahn, H. J. Maris, J. Tauc, Surface generation and detection of phonons by picosecond light pulses. *Phys. Rev. B* **34**, 4129–4138 (1986).

33. P. Ruello, V. E. Gusev, Physical mechanisms of coherent acoustic phonons generation by ultrafast laser action. *Ultrasonics* **56**, 21–35 (2015).

34. M. Lejman, C. Paillard, V. Juvé, G. Vaudel, N. Guiblin, L. Bellaiche, M. Viret, V. E. Gusev, B. Dkhil, P. Ruello, Magnetoelastic and magnetoelectric couplings across the antiferromagnetic transition in multiferroic $BiFeO_3$. *Phys. Rev. B* **99**, 104103 (2019).

35. O. B. Wright, Thickness and sound velocity measurement in thin transparent films with laser picosecond acoustics. *J. Appl. Phys.* **71**, 1617–1629 (1992).

36. V. E. Gusev, Laser hypersonics in fundamental and applied research. *Acustic-acta acustica* **82**, S37–S45 (1996).

37. O. B. Wright, O. Matsuda, Laser picosecond acoustics with oblique probe light incidence. *Rev. Sci. Instr.* **74**, 895–897 (2003).

38. C. Mechri, P. Ruello, V. Gusev, Confined coherent acoustic modes in a tubular nanoporous alumina film probed by picosecond acoustics methods. *New J. Phys.* **14**, 023048 (2012).

39. D. Mounier, P. Picart, P. Babilotte, P. Ruello, J.-M. Breteau, T. Pézeril, G. Vaudel, M. Kouyaté, V. Gusev, Jones matrix formalism for the theory of picosecond shear acoustic pulse detection. *Optic Express* **18**, 6767–6778 (2010).

40. D. Schick, A. Bojahr, R. Shayduk, M. Herzog, C. von Korff Schmising, M. Bargheer, udkm1Dsim—A simulation toolkit for 1d ultrafast dynamics in condensed matter. *Comput. Phys. Rev.* **185**, 651–660 (2014).

41. V. E. Gusev, P. Ruello, Advances in applications of time-domain brillouin scattering for nanoscale imaging. *Appl. Phys. Rev.* **5**, 031101 (2018).

42. D. A. Reis, M. F. DeCamp, P. H. Bucksbaum, R. Clarke, E. Dufresne, M. Hertlein, R. Merlin, R. Falcone, H. Kapteyn, M. M. Murnane, J. Larsson, T. Missalla, J. S. Wark, Probing impulsive strain propagation with x-ray pulses. *Phys. Rev. Lett.* **86**, 3072–3075 (2001).

43. T. Pézeril, P. Ruello, S. Gougeon, N. Chigarev, D. Mounier, J.-M. Breteau, P. Picart, V. Gusev, Generation and detection of plane coherent shear picosecond acoustic pulses by lasers: Experiment and theory. *Phys. Rev. B* **75**, 174307 (2007).

44. D. Schick, M. Herzog, A. Bojahr, W. Leitenberger, A. Hertwig, R. Shayduk, M. Bargheer, Ultrafast lattice response of photoexcited thin films studied by x-ray diffraction. *Struct. Dyn.* **1**, 064501 (2014).

45. I. Ponomareva, L. Bellaiche, Nature of dynamical coupling between polarization and strain in nanoscale ferroelectrics from first principles. *Phys. Rev. Lett.* **101**, 107602 (2008).

46. H. Huyan, L. Li, C. Addiego, W. Gao, X. Pan, Structures and electronic properties of domain walls in $BiFeO_3$ thin films. *Natl. Sci. Rev.* **6**, 669–683 (2019).

47. A. S. Everhardt, S. Matzen, N. Domingo, G. Catalan, B. Noheda, Ferroelectric domain structures in low-strain $BaTiO_3$. *Adv. Electron. Mater.* **2**, 1500214 (2016).

48. E. Gradauskaite, P. Meisenheimer, M. Muller, J. Heron, M. Trassin, Multiferroic heterostructures for spintronics. *Phys. Sci. Rev.* **6**, 0072 (2021).

49. L. W. Martin, A. M. Rappe, Thin-film ferroelectric materials and their applications. *Nat. Rev. Mater.* **2**, 16087 (2017).

50. D. Sando, M. Han, V. Govinden, O. Paull, F. Appert, C. Carrétéro, J. Fischer, A. Barthélémy, M. Bibes, V. Garcia, S. Fusil, B. Dkhil, J. Juraszek, Y. Zhu, X. Ma, V. Nagarajan, Interfacial Strain Gradients Control Nanoscale Domain Morphology in Epitaxial $BiFeO_3$ Multiferroic Films. *Adv. Funct. Mater.* **30**, 2000343 (2020).

51. A. Kumar, R. C. Rai, N. J. Podraza, S. Denev, M. Ramirez, Y.-H. Chu, L. W. Martin, J. Ihlefeld, T. Heeg, J. Schubert, D. G. Schlom, J. Orenstein, R. Ramesh, R. W. Collins, J. L. Musfeldt, V. Gopalan, Linear and nonlinear optical properties of $BiFeO_3$. *Appl. Phys. Lett.* **92**, 121915 (2008).
52. B. Gu, Y. Wang, J. Wang, W. Ji, Femtosecond third-order optical nonlinearity of $BiFeO_3$. *Optics Expr.* **17**, 10970 (2009).

53. D. Sando, P. Hermet, J. Allibe, J. Bourderionnet, S. Fusil, C. Carrétéro, E. Jacquet, J.-C. Mage, D. Dolfi, A. Barthélémy, P. Ghosez, M. Bibes, Linear electro-optic effect in multiferroic $BiFeO_3$ thin films. *Phys. Rev. B* **89**, 195106 (2014).

54. M.-A. Tordeux, J. Barros, A. Bence, P. Brunelle, N. Hubert, M. Labat, A. Nadji, L. Nadolski, P. Lebasque, J.-P. Pollina, C. Evain, Low-alpha operation for the soleil storage ring, in *Proceedings of the 2012 International Particle Accelerator Conference* (IPAC, 2012), pp. 1608.

55. D. Bachiller-Perea, Y.-M. Abiven, J. Bisou, P. Fertey, P. Grybos, A. Jarnac, B. Kanouté, A. Koziol, F. Langlois, C. Laulhé, F. Legrand, P. Maj, C. Menneglier, A. Noureddine, F. Orsini, G. Thibaux, A. Dawiec, First pump-probe-probe hard x-ray diffraction experiments with a 2d hybrid pixel detector developed at the SOLEIL synchrotron. *J. Synchrotron Radiat.* **27**, 340–350 (2020).

56. S. L. Shang, G. Sheng, Y. Wang, L. Q. Chen, Z. K. Liu, Elastic properties of cubic and rhombohedral $BiFeO_3$ from first-principles calculations. *Phys. Rev. B* **80**, 052102 (2009).

57. E. Borissenko, M. Goffinet, A. Bosak, P. Rovillain, M. Cazayous, D. Colson, P. Ghosez, M. Krisch, Lattice dynamics of multiferroic $BiFeO_3$ studied by inelastic x-ray scattering. *J. Phys. Condens. Matter* **25**, 102201 (2013).

58. P. Hemme, P. Djemia, P. Rovillain, Y. Gallais, A. Sacuto, A. Forget, D. Colson, E. Charron, B. Perrin, L. Belliard, M. Cazayous, Elastic properties assessment in the multiferroic $BiFeO_3$ by pump and probe method. *Appl. Phys. Lett.* **118**, 062902 (2021).

59. L. Maerten, A. Bojahr, M. Gohlke, M. Rössle, M. Bargheer, Coupling of GHz phonons to ferroelastic domain walls in $SrTiO_3$. *Phys. Rev. Lett.* **114**, 047401 (2015).

60. S. Yamanaka, T. Maekawa, H. Muta, T. Matsuda, S. Kobayashi, K. Kurosaki, Thermophysical properties of $SrHfO_3$ and $SrRuO_3$. *J. Solid State Chem.* **177**, 3484–3489 (2004).

61. R. O. Bell, G. Rupprecht, Elastic constants of strontium titanate. *Phys. Rev.* **129**, 90 (1963).

Acknowledgments

**Funding:** We acknowledge SOLEIL for provision of synchrotron radiation at CRISTAL beamline (proposal number 20191634) and support from the Agence Nationale de la Recherche (ANR) through the THz-MUFINS project (grant no. ANR-21-CE42-0030). This work was supported by a public grant overseen by the ANR as part of the Investissement d'Avenir programme (LABEX NanoSaclay, reference no. ANR-10-LABX-0035). V.J., G.V., M.C.W., and P.R. are grateful for the financial support by the Région des Pays de la Loire under the FASTFERRO project (Trajectoire Nationale). M.C.W. is grateful for the financial support by the Région des Pays de la Loire under the Etoile Montante Initiative (2022_11808) and the PULSAR Academy. This research was partially supported by the Australian Research Council Centre of Excellence in Future Low-Energy Electronics Technologies (project no. CE170100039) and funded by the Australian Government. **Author contributions:** V.J., B.D., D.S., V.G., and P.R. initiated the project. R.G., V.J., C.L., H.B., G.V., B.D., and P.R. performed the time-resolved XRD experiments. P.H. operated the femtosecond laser source during the experiment at SOLEIL synchrotron. R.G., V.J., C.L., and P.R. analyzed the x-ray diffraction data. The all-optical time-resolved optical reflection experiment was performed and by R.G., V.J., G.V., and A.P. V.J. and P.R. simulated the all-optical time-resolved reflection experiment. D.S. engineered, synthesized, and characterized the single-domain $BiFeO_3$ thin-film sample. V.G. and S.F. identified the polarization direction in the $BiFeO_3$ thin film and characterized its single-domain nature using scanning probe microscopy. V.J., B.D., D.S., V. G., P.R., S. F., C.P., H.B., C.L., M.C.W., R.G., and G.V. evaluated the data. V.J., P.R., B.D., C.L., M.C.W., S. F., V.G., and D.S. wrote the manuscript. All authors contributed to the discussion and to the completion of the manuscript. **Competing interests:** The authors declare that they have no competing interests. **Data and materials availability:** All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials.

Submitted 7 April 2023

Accepted 16 October 2023

Published 15 November 2023

10.1126/sciadv.adi1160
Temporal and spatial tracking of ultrafast light-induced strain and polarization
modulation in a ferroelectric thin film

Ruizhe Gu, Vincent Juvé, Claire Laulhé, Houssny Bouyanfif, Gwenaëlle Vaudel, Aurélie Poirier, Brahim Dkhil, Philippe Hollander, Charles Paillard, Mads C. Weber, Daniel Sando, Stéphane Fusil, Vincent Garcia, and Pascal Ruello

Sci. Adv. 9 (46), eadi1160. DOI: 10.1126/sciadv.adi1160

View the article online

https://www.science.org/doi/10.1126/sciadv.adi1160

Permissions

https://www.science.org/help/reprints-and-permissions