
==== Front
Struct Dyn
Struct Dyn
SDTYAE
Structural Dynamics
2329-7778
American Crystallographic Association

38026579
10.1063/4.0000199
4.0000199
SDY23-AR-00047
ARTICLES
Theory and Modelling
Simulation of ultrafast electron diffraction intensity under coherent acoustic phonons
https://orcid.org/0000-0003-4036-0900
Zhang Yongzhao 1,2,3
https://orcid.org/0000-0003-1180-5798
Li Jun 1
https://orcid.org/0009-0002-2071-2390
Wang Wentao 1,3
https://orcid.org/0000-0002-9473-4050
Tian Huanfang 1
https://orcid.org/0009-0005-4835-9121
Gao Wenli 1,4
https://orcid.org/0000-0002-4808-6347
Li Jianqi 1,3,5a)

https://orcid.org/0000-0002-9363-3901
Sun Shuaishuai 1a)

https://orcid.org/0000-0002-7857-2442
Yang Huaixin 1,3,6a)

1 Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
2 Institute of Quantum Materials and Physics, Henan Academy of Science, Zhengzhou 450046, China People's Republic of China
3 School of Physical Sciences, University of Chinese Academy of Science, Beijing 100190, People's Republic of China
4 School of Physics, Northwest University, Xi'an 710069, People's Republic of China
5 Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, People's Republic of China
6 Yangtze River Delta Physics Research Center Co., Ltd., Liyang, Jiangsu 213300, People's Republic of China
a) Authors to whom correspondence should be addressed: ljq@iphy.ac.cn; sss@iphy.ac.cn; and hxyang@iphy.ac.cn
13 11 2023
11 2023
13 11 2023
10 6 06410206 6 2023
22 10 2023
© 2023 Author(s).
2023
Author(s)
https://creativecommons.org/licenses/by/4.0/ All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Ultrafast electron diffraction has been proven to be a powerful tool for the study of coherent acoustic phonons owing to its high sensitivity to crystal structures. However, this sensitivity leads to complicated behavior of the diffraction intensity, which complicates the analysis process of phonons, especially higher harmonics. Here, we theoretically analyze the effects of photoinduced coherent transverse and longitudinal acoustic phonons on electron diffraction to provide a guide for the exploitation and modulation of coherent phonons. The simulation of the electron diffraction was performed in 30-nm films with different optical penetration depths based on the atomic displacements obtained by solving the wave equation. The simulation results exhibit a complex relationship between the frequencies of the phonons and diffraction signals, which highly depends on the laser penetration depth, sample thickness, and temporal stress distribution. In addition, an intensity decomposition method is proposed to account for the in-phase oscillation and high harmonics caused by inhomogeneous excitation. These results can provide new perspectives and insights for a comprehensive and accurate understanding of the lattice response under coherent phonons.

National Key Research and Development Program of China 10.13039/501100012166 2021YFA1301502 National Natural Science Foundation of China 10.13039/501100001809 U22A6005 12074408 Guangdong Major Scientific Research Project 2018KZDXM061 Youth Innovation Promotion Association of the Chinese Academy of Sciences 10.13039/501100004739 2021009 Scientific instrument developing Project of the Chinese Academy of Sciences YJKYYQ20200055 ZDKYYQ2017000 22017BA10 Strategic Priority Research Programof the Chinese Academy of Sciences XDB25000000 XDB33010100 Beijing Municipal Science and Technology Major Project Z201100001820006 IOP Hundred Talents Program Y9K5051 Postdoctoral Support Program of China 2020M670501 crossmark
==== Body
pmcINTRODUCTION

The generation and propagation of coherent acoustic phonons can produce periodic modulation of the crystal structure, which exists as collective atomic motion, and induce phase transition.1–4 Upon femtosecond-laser excitation, coherent longitudinal acoustic (LA) waves usually appear in the material due to thermoelastic effects or deformation potential generated by the carriers in the sample.5,6 The excitation of the coherent transverse acoustic (TA) waves is noted when the instability of the crystal structure or disoriented lattice is considered.7–10 In addition, several works have investigated the mechanism of piezoelectrically induced coherent phonons.11,12

Laser-induced coherent phonons typically contain multiple frequencies. In thick films, where the optical penetration depth is considerably smaller than the sample thickness, a spatially non-uniform stress distribution in longitudinal direction is excited. Both even and odd high harmonics can be observed in the coherent oscillation, corresponding to the propagating strain waves bouncing back and forth in the film.13,14 When the optical penetration depth is much larger than the thickness, uniform stress is generated in the film, resulting in a standing wave with odd high harmonics only.13,15 The modulation effect of coherent acoustic phonons on the crystal structure makes them easily detectable by ultrafast electron diffraction (UED). However, common samples used in UED measurements have an optical penetration depth comparable to the film thickness,16 where the traveling and standing waves can coexist in the system, thereby complicating the analysis. In particular, LA and TA phonons can be simultaneously excited in a given system, generating complex multifrequency oscillation signals.8,17 Structural fluctuations induced by coherent phonons can be manifested in the diffraction intensity. A crucial yet rarely discussed issue is whether the experimentally observed oscillation frequencies in the diffraction can directly reflect the intrinsic phonon frequencies.8,9,18 Recently, the complex behavior of electron diffraction under high harmonics and lattice distortion in graphite nanofilms has been discussed in detail, raising the possibility of the selective excitation of higher harmonics.19 However, several issues are yet to be discussed in detail, such as the effect of higher harmonics in TA phonons on the diffraction and mechanisms of in-phase and out-of-phase oscillation of the Friedel pairs.20

In previous several UED works8,9,17,20 reporting diffraction intensity oscillation in the tens of GHz range, multiple-frequency is usually attributed to the coexistence of TA and LA phonons, and in-phase oscillation is considered to be caused by initial deviation parameter. Here, we propose that the inhomogeneous of laser excitation has a strong influence on the diffraction intensity (both oscillation phase and frequency). In this work, we theoretically analyzed the complex correspondence between the frequencies of laser-induced coherent acoustic phonons, including the TA and LA modes, and electron diffraction intensity in two hypothetical 30-nm films with laser penetration depths of 10 and 600 nm. In addition, several excitation profiles with the rise time of 0–60 ps were used to investigate their influence on the atomic motion to give a comprehensive understanding of diffraction signal, such as the phase shift of oscillations observed in previous works.21–24 Our results could give a new interpretation of the high harmonic diffraction signals, in-phase and out-of-phase oscillations, and provide further understanding of the uniformity of the laser distribution.

EXCITATION CONDITIONS

We specified that the laser is first irradiated on the upper surface of the sample, and the film thickness was set to 30 nm. Considering the region of interest (center of the Gauss laser) and timescale (<10 ps), a uniform excitation profile in the lateral direction is assumed, and the diffusion of the carrier and heat along the lateral direction are ignored. However, the longitudinal excitation profile is distinct due to the difference in the optical penetration depth. The lattice evolution under homogeneous and inhomogeneous excitation is obtained by assuming different optical absorption coefficients α1 and α2. As shown in Fig. 1(a), the distribution of the laser energy in the film is approximately uniform when the penetration depth δ1 ( 1/α1) is 600 nm. Meanwhile, when δ2 ( 1/α2) is 10 nm, a significant inhomogeneity is created with the laser confined to a limited film depth. The carrier and heat diffusion along the longitudinal direction will smooth the inhomogeneity of laser deposited energy in different timescale. The carrier diffusion mainly occurs before the rapid relaxation through electron-phonon coupling (∼ps), then the heat diffusion dominates the thermalization within the sample, which may take several nanoseconds for a thickness of 30 nm, depending on thermal diffusivity.25,26 To simplify, we ignore the carrier and thermal diffusion along the longitudinal direction.

FIG. 1. Atomic motion under coherent acoustic phonons. (a) Laser condition used in the calculation. (b) Temporal excitation profile of stress σext in the calculation. (c) Schematic of the atomic displacements in shear mode. The total oscillation is the superposition of the equilibrium position, fundamental, and higher harmonics.

Compared with the LA breathing mode routinely excited in thin films caused by thermal-elastic stress, the photo-excitation of coherent TA shear modes requires sophisticated experimental configurations. There have been several articles reporting ultrafast coherent shear modes driven by different mechanism (transverse inhomogeneous excitation in Si;15 inverse piezoelectric effect in BiFeO3;11 structural instability in VTe2;8 miscut along a low symmetry face in TaAs;27 demagnetization in FePS328). The common feature for them is the symmetry breaking in the surface plane. In order to make the excitation of TA phonons reasonable, we assume that the film has a monoclinic crystal structure with lattice constants a = 3.52Å, b = 3.52 Å, c = 7 Å, α = 90°, β = 90.3°, and γ = 90°. The structure may appear to shear along the b-direction under laser excitation due to the structural instability. In laser-induced breathing mode, the atomic displacements are typically proportional to the laser fluence. To more conveniently study the behavior of electron diffraction under different shear displacements, we assume that the amplitude of the shear mode and the laser are also positively proportional, just like the breathing mode. In this case, the amplitudes of both shear and breathing of the lattice can be wrote as βQ, where Q is the absorbed laser fluence, and β is the linear expansion/shear ratio with respect to the lattice constant induced by the unit absorbed laser energy. The stress induced by the laser excitation can be denoted by σext as a function of time t and film depth z.29 The temporal response of the atomic motion would exhibit sine-like (same as impulsive excitation30) or cosine-like (same as displacive excitation31) oscillation depending on the excitation mechanisms. In our calculation, the excitation temporal response used is step-function like, as shown in Fig. 1(b), which is described with exponential functions with different rise times τ, and indicates the absorption of laser energy by the sample is finally stored in the lattice. In this case, all atoms are driven to a new equilibrium position, corresponding to the cosine-like response. Furthermore, given the exponential decay of laser fluence inside the sample, the stress can be expressed as σextz,t=−CβQ ftgz=−CβQ1−exp−tτexp(−αz), (1)

where C is the elastic modulus of the crystal along different directions, which is Young's modulus E and shear modulus G for the breathing and shear mode, respectively. The values of these parameters used are listed in Table I, which are set according to the typical values of two dimensional van der Waals materials.32

TABLE I. Parameters used in the calculation of atomic motion.

Definition	Symbol	Value	
Elastic modulus	E, G	50, 22 GPa	
Linear expansion/shear rate	β	2.5 × 10−3 (mJ/cm2)−1	
Absorbed laser fluence	Q	2.0 mJ/cm2	
Density	ρ	5.5 g/cm3	
Velocity of strain wave	vTA, vLA	2.0, 3.0 km/s	
Absorption coefficient	α1, α2	1.7 × 103, 1.0 × 106 cm−1	

CALCULATION OF THE ATOMIC DISPLACEMENTS

The generation and propagation of the coherent acoustic phonons can be described by the atomic displacement μz,t calculated from the wave equation in continuous medium.5 When the laser beam diameter is considerably larger than the laser penetration depth, it can be assumed that only the coherent phonons propagating along the out-of-plane direction are triggered in the region far from the sample edge. Thus, the wave equation can be directly expressed in one-dimensional form, as follows: ∂2μ∂t2=v2∂2μ∂z2+λρ∂μ∂t+∂σextρ∂z, (2)

where λ is the damping factor that determines the lifetime of the coherent phonons inside the film, which is set as zero for explicitly. The velocity of the out-of-plane TA and LA wave is set as vTA=2.0 and vLA=3.0 km/s, respectively. All the parameters used for our numerical calculation are listed in Table I. For a freestanding film, the free boundary conditions is ∂μ∂zz=0,t=∂μ∂zz=d,t=0, (3)

where d is the thickness of the film. Steady-state solutions to the wave Eq. (2) under free boundary conditions and the initial condition ∂μ∂tz,0=0 can be described as a superposition of standing waves, μz,t|fj≠0=Aj·cosjπdzsin(2πfj·t+φj). (4)

The value of amplitude Aj depends on the temporal and spatial distribution of stress σext, which will be discussed in following text. The frequency of the jth mode has the form fj≠0=jv/2d. In a cosine-like response, the phase φj has the value of π/2. Here, we defined the new atomic equilibrium position after laser excitation as μz,t|fj=0, in which fj=0=0. The actual displacement of the atoms μz,t at depth z is the sum of μz,t|fj≠0 and μz,t|fj=0 with the form μz,t|fj=0=μz,0−∑jAj·cosjπdz, (5)

where μz,0 is the initial displacement of each atoms at depth z. In a perfect crystal, μz,0=0, whereas in an actual crystal, μz,0 has nonzero values for a disordered structure. Substituting the stress σext [Eq. (1)] into the wave Eq. (2), the atomic displacements in the shear or breathing mode can be numerical solved with the initial condition μz,0=∂μ∂tz,0=0 and the free boundary conditions (3). Figure 1(c) shows the calculated equilibrium position μz,t|fj=0 and oscillation mode Aj·cosjπdz in shear mode (j = 1, 2, and 3, corresponding to the fundamental, second harmonic, and third harmonic mode, respectively).

Dynamic simulation of depth-dependent atomic displacement in the two films with optical penetration depth of 10 and 600 nm is shown in Figs. 2(a) and 2(b), respectively (The rise time τ = 0 ps). The homogeneous absorption of laser leads to the excitation of a typical standing wave. However, for inhomogeneous absorption, in addition to the standing wave, a traveling wave propagating along the surface normal becomes evident. Figure 2(c) plots the calculated shear displacements at the top layer of the film (z = 0) with the optical penetration depth δ1= 600 nm, showing an oscillation period T = 30 ps. The position of the first peak and amplitude of the oscillations under different rise times exhibits distinct features. The longer rise time τ could attenuate the amplitude and push oscillation toward a positive direction. We used Δt to denote the shift of the first peak position relative to the position, where delay = T/2 = 15 ps. The phase shift 2πΔt/T will increase with the rise time τ, but eventually converge to π/2 which can be analyzed by Green function. The corresponding information on the film with the optical penetration depth δ2= 10 nm is shown in Fig. 2(d), in which the phase shift and amplitude attenuation exhibit similar behaviors as in the homogeneous condition.

FIG. 2. Atomic motions with different rise times τ of the stress σext. (a, b) Dynamic simulation of depth-dependent atomic displacement in the two films with optical penetration depth of 10 and 600 nm, respectively. The rise time τ used is 0 ps. (c) Atomic motions at the top surface with the rise time τ of 0, 10, and 40 ps in the film with the optical penetration depth δ1= 600 nm. (d) Corresponding information of the film with the optical penetration depth δ2= 10 nm.

To quantify the effect of the rise time τ on the oscillation amplitude with the specific frequency, fast Fourier transform (FFT) was conducted on the oscillation profile [black curves in Figs. 2(c) and 2(d)] with the rise time τ = 0 ps, as shown in Fig. 3(a) and 3(c). By performing FFT of atomic motions at the top surface calculated with different rise time τ, the amplitude of each frequency against τ/T can be extracted as shown in Fig. 3(b) and 3(d), in which the positions where the amplitude drops by one order of magnitude are marked with the vertical black lines. The results shown in Fig. 3 indicate the suppressed oscillation in both homogeneous and inhomogeneous excitation as the rise time increases. In addition, the higher harmonics are more sensitive to the rise time. It is worth noting that there is only odd harmonics under homogeneous excitation, where the longitudinal excitation profile (and also the boundary condition) has central inversion symmetry about the center of the sample; however, the even harmonics appears under inhomogeneous excitation where the inversion symmetry is broken.

FIG. 3. Amplitude suppression effect related to the rise time τ of the stress σext. (a) FFT of the atomic displacements at the top surface under homogeneous excitation with τ/T = 0. (b) Amplitude of the frequency in (a) against the rise time τ/T. (c, d) Corresponding information under inhomogeneous excitation.

SIMULATION AND DECOMPOSITION OF THE DIFFRACTION INTENSITY

To investigate the effect of the coherent phonon modes on the diffraction, kinematical diffraction simulations were performed in the proposed crystal ( a = 3.52 Å, b = 3.52 Å, c = 7 Å, α = 90°, β = 90.3°, γ = 90°). The sample surface is parallel to the plane where lattice vectors a and b lie in, similar to the layered films obtained by mechanical exfoliation. In the calculation of the shear mode, the direction of the atomic displacement on the top surface is set along the a direction [Fig. 4(a)]. For the breathing mode, the displacement is along the c* direction. For the diffraction simulation and analyses in this and next section, the rise time τ in stress σext and damping factor λ are set as zero for simplicity. Subsequently, Eqs. (1) and (2) can be totally determined, and the atomic displacements μz,t can be numerically solved, which are used for the simulation of the diffraction in the following text.

FIG. 4. Decomposition of the diffraction intensity. (a) Schematic of the shear mode. (b) Influence of the shear mode on the reciprocal rod. (c) Distinction between the experimentally detected and theoretically calculated intensity distribution against the deviation parameter sz in the 27-nm MoTe2 film. The red line shows the Gaussian fit of the reciprocal rod. (d) Intensity simulation and decomposition of the Friedel pairs (200) and (2¯00) after laser excitation with the fluence of 2 mJ/cm2 and penetration depth of 10 nm. Total intensity change ΔI/I≈ΔA+ΔIΔsz. The red and blue curves represent the intensity change of the Bragg peaks (200) and (2¯00), respectively. Green solid line: intensity change ΔIΔsz caused by the shift Δsz of the reciprocal rod. Black dash line: amplitude change ΔA of the reciprocal rod. (e) Corresponding information of the film with the optical penetration depth δ= 600 nm. (f) The respective oscillation amplitude AmpΔsz and AmpΔA of ΔIΔsz and ΔA in the Bragg peaks with different scattering vectors. (g) The values of AmpΔsz and AmpΔA under different laser fluences.

In the crystal containing Nx, Ny, and Nz unit cells along the three fundamental vectors a, b, and c, the reciprocal rod [right of Fig. 4(b)] of Bragg peak with the reciprocal vector of G can be obtained by introducing the calculated atomic displacements μz,t into the kinematical diffraction equation33 Ikins,t∝A2=sin2(πNxG·a)sin2(πG·a)sin2(πNyG·b)sin2(πG·b)|∑nz=1NzFnzexp[2πiG·(rn+μz,t)]|2, (6)

where rn is the lattice position with rn=nxa +nyb +nzc. Fnz is the structure factor of the unit cell at nz-th layer. For the diffraction simulation, we propose a decomposition method to offer a new perspective on the diffraction intensity in ultrafast experiments. As shown in Figs. 4(a) and 4(b), when lattice vectors c changes to c′=c+r with the shear displacement r=k1a, reciprocal vector a′*=b×(c+r)/V=a*−k1c*. Similarly, if r=k2b, reciprocal vector b′*=b*−k2c*. As a result, all Bragg peaks (hkl) with h ≠ 0 will undergo a shift Δsz =  −hk1c* along the c* direction in reciprocal space, resulting in the intensity redistribution Ikin′sz=Ikinsz−Δsz and the relative intensity change ΔIΔsz|sz=[Ikin′(sz)−Ikinsz]/Ikinsz, where sz is deviation parameter along the c* direction. Due to the different atomic displacements of each layer, Δsz should have the form Δsz∝−hk1¯c*=−hc*∑nz=1Nz(k1,nz/Nz) for total crystal. In the kinematical diffraction theory [Eq. (6)], the diffraction intensity is the coherent superposition of each lattice layer. Meanwhile, the different vibration amplitudes of each layer reduce this coherence and introduce decrease in diffraction amplitude for all peaks (hkl) with h ≠ 0 (in the case r=k1a). We use ΔA=[Ikin′Δsz−Ikin(0)]/Ikin(0) to describe the amplitude change of the reciprocal rod.

In the actual experiments, the flatness of the sample surface dominates the width of reciprocal rod,9,20 which can be one order of magnitude larger than the theoretical value determined by Eq. (6). We have experimentally measured the reciprocal rod length of 2H-MoTe2 thin films with a thickness of 27 nm and selected area of 5 μm as shown in Fig. 4(c). The samples used for the measurement were obtained by mechanical exfoliation from a bulk crystal, and the thicknesses were determined by electron energy loss spectroscopy. The reciprocal rod was performed through the (110) peak with the electron beam incident along the [001] zone axis of the crystal. We use the Gaussian function to fit the experimentally measured results and substitute back into the simulation, which has the following form: Iexp(sz)=exp−2sz2σ2, (7)

where σ = 0.035 Å−1. The intensity change ΔIΔsz|sz induced by Δsz should be calculated as ΔIΔsz|sz=[Iexp(sz−Δsz)−Iexpsz]/Iexpsz, while the amplitude change ΔA still has the form ΔA=[Ikin′Δsz−Ikin(0)]/Ikin(0). The total intensity change at a specific deviation parameter sz can be decomposed as ΔII=(1+ΔA)·Iexp(sz−Δsz)Iexpsz−1=ΔA[1+ΔIΔsz]+ΔIΔsz≈ΔA+ΔIΔsz. (8)

Under inhomogeneous excitation, the simulated temporal intensity evolutions of the Bragg peak (200) and (2¯00) with the same deviation parameter sz=0.01Å−1 are indicated by the red and blue curves in Fig. 4(d), respectively. The corresponding intensity decomposition is shown by black dashed line ( ΔIΔsz) and green solid line ( ΔA). ΔIΔsz mainly shows fundamental mode features (fundamental with j = 1 is dominant, with the oscillation period T = 30 ps), whereas high harmonic features are dominant in ΔA (second harmonic with j = 2 is dominant, with period of T/2). Another key point is that the oscillation of ΔIΔsz in (200) and (2¯00) exhibits an out-of-phase feature [ phase(200)=π, phase(2¯00)=0] whereas ΔA shows an in-phase oscillation (both phases are π) in the Friedel pairs. Typically, the Friedel pairs have the same diffraction intensity at the same initial deviation parameter ( s0,G=s0,−G), which is determined by the central symmetry of the structural factors. However, the ΔIΔsz and ΔA would cause the unequal change of the diffraction in Friedel pairs G and –G. Oscillations caused by ΔA within a Friedel pairs are always in-phase, while whether oscillations caused by ΔIΔsz within a Friedel pairs are out-of-phase or in-phase sensitively depends on the initial deviation parameter, which is related to the tilt angle. Corresponding information of the film with the optical penetration depth δ= 600 nm is shown in Fig. 4(e). In homogeneous excitation, intensity change ΔIΔsz is dominant. In addition, ΔIΔA mainly shows the oscillation with period of T/4 and a very small amplitude.

In order to provide the dependence of ΔIΔsz and ΔIΔA on scattering vector q, the intensity oscillation of peak (h00) (h = ±1, h = ±2, and h = ±3) at the same deviation parameter sz=0.01Å−1 is calculated. We denote the oscillation amplitude of ΔA and ΔIΔsz by AmpΔA and AmpΔsz, repectively, and specify Amp <0 when phase =π, Amp >0 when phase =0, by which the phase distribution can be directly observed from the sign of the amplitude in Figs. 4(f) and 4(g). The AmpΔsz and AmpΔA are plotted against the scattering vector q in Fig. 4(f). For peak (h00), the shift distance Δsz∝−hk1¯c* is proportional to h. In addition, the intensity change can be approximately considered as linear with Δsz, considering the shape of the reciprocal rod [Eq. (6) and Fig. 4(c)]. Thus, the AmpΔsz satisfies a linear relationship with the scattering vector q. While the AmpΔA satisfies the linear relationship with q2, which can be attributed to the disorder of the entire crystal. When treating the film as a supercell, ΔA can be analogous to the Debye–Waller effect, which shows a linear relationship with q2, suggesting it is a collective intensity decrease due to the overall structural disorder in the sample. The opposite signs of AmpΔsz in (h00) and (h¯00) indicate the out-of-phase feature [ Phase(h00)=π, Phase(h¯00)=0], whereas the oscillation of ΔA shows a consistent phase of π in all eight Bragg peaks. In addition, we calculated the diffraction intensity of peak (200) at different laser fluence (Q) and performed a similar intensity decomposition [Fig. 4(g)]. As mentioned above, Δsz∝−hk1¯c*. Given the linear shear rate β, the mean expansion rate k1¯∝Q. Thus, AmpΔsz satisfies an approximately linear relationship with the energy density Q. However, AmpΔA depends on the coherence of each lattice layer, which is dramatically affected by the lattice disorder. We used the variance ( S) of lattice constant c′ under the modulation of coherent phonons to describe the disorder, expressed as S=∑nz(cnz′−c′¯)2∝∑nz(k1, nz−k1¯)2∝k1, nz2∝Q2, (9)

which is based on the relation cnz′=c−k1,nza and c′¯=∑nz(cnz′/Nz). Equation (9) achieved S∝Q2. The linear relationship between AmpΔA and Q2 shown in Fig. 4(g) indicates the amplitude of the reciprocal rod is proportional to the variance S.

We use the mean change of lattice constant Δc′¯ ( Δc′¯=∑nz(cnz′−c′¯)/Nz) to denote the mean lattice expansion or shearing. Figure 5 plots the temporal evolution of Δc′¯ and S with the modulation of shear mode in the crystal. It is clear that the oscillation of Δc′¯ shows an period of 30 ps, while S show high harmonic features, which is consistent with the results shown in Figs. 4(d) and 4(e). These results also support the conclusion that the intensity change ΔIΔsz is mainly caused by the mean change of lattice, while ΔIΔA is stemming from the lattice disorder.

FIG. 5. Temporal evolution of mean change [ Δc′¯=∑nz(cnz′−c′¯)/Nz] and the variance [S=∑nz(cnz′−c′¯)2] of lattice constant with the modulation of shear mode. (a) and (b) show the corresponding results in inhomogeneous ( δ = 10 nm) and homogeneous excitation ( δ = 600 nm), respectively.

FREQUENCY ANALYSIS

After decomposing the intensity, the frequency distribution of the diffraction intensity and its relationship to the coherent phonon mode can be clarified. The shear modes can induce the intensity change of the in-plane peaks, whereas the breathing modes exert influence on the out-of-plane peaks. Given that the strain wave propagates along the reciprocal vector c*, the FFT of the diffraction intensity of the peak (200) and (02¯2) is performed to analyze the influence of the TA phonons (shear mode) and LA phonons (breathing mode), respectively. The FFT results of peak (200) under shear mode induced by homogeneous ( δ1= 600 nm) and inhomogeneous ( δ2= 10 nm) excitation is shown in Fig. 6(a), where the upper and lower curves are obtained based on the temporal oscillation data shown in Figs. 4(f) and 4(d), respectively. Figure 6(b) shows the corresponding information of peak (02¯2) under the breathing mode. It should be noted that, in order to view peak (02¯2), the crystal should be tilted to [011] zone axis. In the actual experiment, we used the sine function y=Aisin(2πfit+φi) to describe the oscillation with the specific frequency fi. All the oscillations shown in Fig. 6 have a phase of 0 or π due to the used step-function-like temporal excitation profile ( σext with τ=0). We specify the amplitude of FFT as >0 when φi=π, while <0 when φi=0, resulting in the negative values shown in Fig. 6.

FIG. 6. Frequency domain information of the diffraction intensity oscillation. (a) and (b) FFT of the (200) and (02¯2) decomposed diffraction intensity in shear and breathing modes under homogeneous and inhomogeneous excitation. (c) and (d) FFT of the (200) and (02¯2) diffraction intensity in shear and breathing modes with the specific phonon frequency based on the atomic equilibrium position in the homogeneous condition. (e) and (f) FFT of the (200) and (02¯2) diffraction intensity in shear and breathing modes with the specific phonon frequency based on the equilibrium position in inhomogeneous condition.

Figures 6(a) and 6(b) indicate that the homogeneous excitation mainly causes the shift of the reciprocal rod, which is manifested by the intensity oscillations ΔIΔsz with frequency f1. Meanwhile, the oscillation of ΔA at the frequency of 4f1 indicates slight changes in the amplitude of the reciprocal rod. In contrast, inhomogeneous excitation causes more drastic reductions in the intensity of the reciprocal rod, resulting in a strong oscillation of ΔA with frequency of 2f1. These findings apply to both the breathing and shear modes; however, the behavior in the breathing mode is more complex. As shown in Fig. 6(b), the breathing mode causes additional oscillation of ΔA with frequency f1, which is not related to the uniformity of the optical excitation.

The frequency of atomic displacements under coherent phonons contains multiple components, mainly includes the fundamental ( f1), second harmonic ( f2), and third harmonic ( f3) components. To analyze the origin of the various frequencies of the diffraction in Figs. 6(a) and 6(b), atomic displacements under each frequency component based on the equilibrium position are introduced into the diffraction simulation. The results are shown in Figs. 6(c)–6(f), in which μ homo and μinhomo denote the atomic equilibrium position after laser excitation. While, μfi represents the atomic displacements induced by the vibration with frequency fi. Corresponding schematic can be seen in Fig. 1(c).

For the shear mode under homogeneous excitation [Fig. 6(c)], the fundamental mode (frequency: f1) generates the oscillation of diffraction intensity ΔIΔsz with the same frequency f1, which behaves as an out-of-phase oscillation in the Friedel pairs. The third harmonic mode produces the oscillation of ΔA and ΔIΔsz with frequencies of 6f1 and 3f1, respectively. Notably, the simultaneous presence of the fundamental and third harmonic modes exhibit a weak 4f1 ( ΔA) signal in the diffraction, indicating that the influence of each mode is not completely independent. Thus, the oscillation of ΔA exhibits 4f1 signals for both the shear and breathing modes in the film with the penetration depth of 600 nm. These conclusions remain valid for the breathing mode under homogeneous excitation [Fig. 6(d)], except the influence of the fundamental mode with the atomic displacements of μhomo+μf1, which causes additional oscillations of ΔA with frequency of f1.

The comparison of Figs. 6(c) and 6(e) [or Figs. 6(d) and 6(f)] indicates that, for the inhomogeneous excitation, the oscillation of ΔA shows obvious frequency i·f1 under the i-th order harmonic (i > 1), which is absent under the homogeneous condition. In addition, the amplitude of the oscillation with frequency 2i·f1 exhibits obvious attenuation, which is attributed to the asymmetric distribution of the atomic equilibrium position with respect to the center of the film. Diffraction signals in inhomogeneous excitation exhibit a 2f1 oscillation not only because of the displacement μf2 of atoms at the second harmonic, but also by the equilibrium position of the atoms μinhomo. The equilibrium position of atoms in the oscillation has great effect on diffraction, that is the reason why μhomo+ μf3 and μinhomo+ μf3 show significant difference in diffraction [Figs. 6(c) and 6(e), or Figs. 6(d) and 6(f)], and why the second harmonics manifested in the diffraction intensity oscillation is much stronger than the true phonon oscillation amplitude [Figs. 3(c), 6(c), and 6(e)]. The asymmetric atomic distribution also depends on the initial position μz,0 of each lattice layer [Eq. (5)], indicating that μz,0 can significantly influence the oscillation of ΔA with frequency i·f1 under the ith order harmonics (i > 1). The change in the initial position induced by the lattice distortion was discussed in previous works.34,35 However, this was not considered in this work.

The relationship between the frequencies of the diffraction intensity and the coherent phonons now can be clarified. In the breathing or shear modes generated by the homogeneous excitation, the oscillation of the diffraction intensity in (02¯2) or (200) with frequency f1 is attributed to the fundamental mode, whereas the frequency 4f1 stems from the coupling of the fundamental and third harmonic modes. For the inhomogeneous excitation, additional 2f1 signal is ascribed to the second harmonic mode, whereas the weak 3f1 signal is stemming from the third harmonic mode. It should be mentioned that the finite electron–phonon coupling constant in the actual material would result in a non-zero value of τ. According to the results in Fig. 3(b), the high-frequency signals will be significantly attenuated, and therefore, it is difficult to experimentally observe the 4f1 signal in homogeneous excitation.

It should be noted that the sample is assumed to be absolutely smooth and there are no wrinkles in the film. The LA and TA phonons are simply described by one dimensional wave equation and exist independently in the sample. However, the real material may consist of many wrinkles or micro-patches separated by defects, and each patch may be different in size and orientation, resulting in its own oscillation direction. In addition, TA and LA phonons could be excited simultaneously in the sample, leading to more complicated behavior in oscillation frequency and phase.

CONCLUSION

In this study, the phase shift and amplitude suppression in oscillation due to the rise time of the stress were analyzed. The results of this work are useful for understanding the generation of coherent phonons under different mechanisms. For example, long electron–phonon coupling time or small thickness can lead to a large relative rise time τ/T. This significantly weakens the amplitude of thermoelasticity induced coherent phonons, especially the higher harmonics. In addition, we demonstrated the complex behaviors of the diffraction intensity under structural modulation of the coherent LA and TA phonons generated by homogeneous and inhomogeneous distribution. The correspondence between the frequencies of the diffraction intensity oscillation and high harmonics of the strain waves was analyzed by decomposing the diffraction intensity. Thus, the phonon-induced in-phase and out-of-phase intensity oscillation in the Friedel pairs and unequal change were elaborated. These analyses deepened the understanding of the structural behaviors under coherent acoustic phonons and their corresponding influence on electron diffraction, which facilitated further structural modulations by the coherent phonons in the UED device.

ACKNOWLEDGMENTS

This work was supported by the National Key Research and Development Program of China (Grant No. 2021YFA1301502), the National Natural Science Foundation of China (Grant Nos. U22A6005 and 12074408), the Guangdong Major Scientific Research Project (No. 2018KZDXM061), the Youth Innovation Promotion Association of CAS (No. 2021009), the Scientific Instrument Developing Project of the Chinese Academy of Sciences (Grant Nos. YJKYYQ20200055, ZDKYYQ2017000, and 22017BA10), the Strategic Priority Research Program (B) of the Chinese Academy of Sciences (Grant Nos. XDB25000000 and XDB33010100), Beijing Municipal Science and Technology Major Project No. Z201100001820006, the IOP Hundred Talents Program (No. Y9K5051), the Postdoctoral Support Program of China (No. 2020M670501), and the Synergetic Extreme Condition User Facility (SECUF).

AUTHOR DECLARATIONS

Conflict of Interest

The authors have no conflicts to disclose.

Author Contributions

Yongzhao Zhang: Conceptualization (equal); Data curation (lead); Investigation (lead); Methodology (equal); Visualization (lead); Writing – original draft (lead). Jun Li: Funding acquisition (equal); Methodology (equal). Wentao Wang: Writing – review & editing (equal). Huanfang Tian: Writing – review & editing (equal). Wenli Gao: Writing – review & editing (equal). Jianqi Li: Funding acquisition (equal); Writing – review & editing (equal). Shuaishuai Sun: Conceptualization (equal); Formal analysis (equal); Methodology (equal); Writing – review & editing (equal). Huaixin Yang: Conceptualization (equal); Funding acquisition (equal); Writing – review & editing (equal).

DATA AVAILABILITY

The data that support the findings of this study are available from the corresponding authors upon reasonable request.
==== Refs
1. G. Huitric , M. Rodriguez-Fano , L. Gournay , N. Godin , M. Herve , G. Privault , J. Tranchant , Z. Khaldi , M. Cammarata , E. Collet , E. Janod , and C. Odin , “ Impact of the terahertz and optical pump penetration depths on generated strain waves temporal profiles in a V2O3 thin film,” Faraday Discuss. 237 , 389–405 (2022).10.1039/D2FD00013J 35665795
2. C. Mariette , M. Lorenc , H. Cailleau , E. Collet , L. Guerin , A. Volte , E. Trzop , R. Bertoni , X. Dong , B. Lepine , O. Hernandez , E. Janod , L. Cario , V. T. Phuoc , S. Ohkoshi , H. Tokoro , L. Patthey , A. Babic , I. Usov , D. Ozerov , L. Sala , S. Ebner , P. Bohler , A. Keller , A. Oggenfuss , T. Zmofing , S. Redford , S. Vetter , R. Follath , P. Juranic , A. Schreiber , P. Beaud , V. Esposito , Y. Deng , G. Ingold , M. Chergui , G. F. Mancini , R. Mankowsky , C. Svetina , S. Zerdane , A. Mozzanica , A. Bosak , M. Wulff , M. Levantino , H. Lemke , and M. Cammarata , “ Strain wave pathway to semiconductor-to-metal transition revealed by time-resolved x-ray powder diffraction,” Nat. Commun. 12 , 1239 (2021).10.1038/s41467-021-21316-y 33623010
3. M. Zhang , G. L. Cao , H. F. Tian , S. S. Sun , Z. W. Li , X. Y. Li , C. Guo , Z. Li , H. X. Yang , and J. Q. Li , “ Picosecond view of a martensitic transition and nucleation in the shape memory alloy Mn50Ni40Sn10 by four-dimensional transmission electron microscopy,” Phys. Rev. B 96 (17 ), 174203 (2017).10.1103/PhysRevB.96.174203
4. I. A. Mogunov , F. Fernandez , S. Lysenko , A. J. Kent , A. V. Scherbakov , A. M. Kalashnikova , and A. V. Akimov , “ Ultrafast insulator-metal transition in VO2 nanostructures assisted by picosecond strain pulses,” Phys. Rev. Appl. 11 , 014054 (2019).10.1103/PhysRevApplied.11.014054
5. P. Ruello and V. E. Gusev , “ Physical mechanisms of coherent acoustic phonons generation by ultrafast laser action,” Ultrasonics 56 , 21–35 (2015).10.1016/j.ultras.2014.06.004 25038958
6. T. Saito , O. Matsuda , and O. B. Wright , “ Picosecond acoustic phonon pulse generation in nickel and chromium,” Phys. Rev. B 67 (20 ), 205421 (2003).10.1103/PhysRevB.67.205421
7. Y. P. Qi , M. X. Guan , D. Zahn , T. Vasileiadis , H. Seiler , Y. W. Windsor , H. Zhao , S. Meng , and R. Ernstorfer , “ Traversing double-well potential energy surfaces: Photoinduced concurrent intralayer and interlayer structural transitions in XTe2 (X = Mo, W),” ACS Nano 16 (7 ), 11124–11135 (2022).10.1021/acsnano.2c03809 35793703
8. A. Nakamura , T. Shimojima , Y. Chiashi , M. Kamitani , H. Sakai , S. Ishiwata , H. Li , and K. Ishizaka , “ Nanoscale imaging of unusual photoacoustic waves in thin flake VTe2,” Nano Lett. 20 (7 ), 4932–4938 (2020).10.1021/acs.nanolett.0c01006 32463678
9. A. Nakamura , T. Shimojima , M. Matsuura , Y. Chiashi , M. Kamitani , H. Sakai , S. Ishiwata , H. Li , A. Oshiyama , and K. Ishizaka , “ Evaluation of photo-induced shear strain in monoclinic VTe2 by ultrafast electron diffraction,” Appl. Phys. Express 11 (9 ), 092601 (2018).10.7567/APEX.11.092601
10. O. Matsuda , O. B. Wright , D. H. Hurley , V. E. Gusev , and K. Shimizu , “ Coherent shear phonon generation and detection with ultrashort optical pulses,” Phys. Rev. Lett. 93 (9 ), 095501 (2004).10.1103/PhysRevLett.93.095501 15447110
11. M. Lejman , G. Vaudel , I. C. Infante , P. Gemeiner , V. E. Gusev , B. Dkhil , and P. Ruello , “ Giant ultrafast photo-induced shear strain in ferroelectric BiFeO3,” Nat. Commun. 5 , 4301 (2014).10.1038/ncomms5301 24980954
12. Y. C. Wen , T.-S. Ko , T.-C. Lu , H.-C. Kuo , J.-I. Chyi , and C.-K. Sun , “ Photogeneration of coherent shear phonons in orientated wurtzite semiconductors by piezoelectric coupling,” Phys. Rev. B 80 (19 ), 195201 (2009).10.1103/PhysRevB.80.195201
13. A. Nakamura , T. Shimojima , and K. Ishizaka , “ Finite-element simulation of photoinduced strain dynamics in silicon thin plates,” Struct. Dyn. 8 (2 ), 024103 (2021).10.1063/4.0000059 33907699
14. P. Soubelet , A. A. Reynoso , A. Fainstein , K. Nogajewski , M. Potemski , C. Faugeras , and A. E. Bruchhausen , “ The lifetime of interlayer breathing modes of few-layer 2H-MoSe2 membranes,” Nanoscale 11 (21 ), 10446–10453 (2019).10.1039/C9NR02447F 31112191
15. F. Hudert , A. Bruchhausen , D. Issenmann , O. Schecker , R. Waitz , A. Erbe , E. Scheer , T. Dekorsy , A. Mlayah , and J. R. Huntzinger , “ Confined longitudinal acoustic phonon modes in free-standing Si membranes coherently excited by femtosecond laser pulses,” Phys. Rev. B 79 (20 ), 201307 (2009).10.1103/PhysRevB.79.201307
16. X. Y. Sun , S. S. Sun , and C. Y. Ruan , “ Toward nonthermal control of excited quantum materials: Framework and investigations by ultrafast electron scattering and imaging,” C. R. Phys. 22 , 15–73 (2021).10.5802/crphys.86
17. M. Harb , W. Peng , G. Sciaini , C. T. Hebeisen , R. Ernstorfer , M. A. Eriksson , M. G. Lagally , S. G. Kruglik , and R. J. D. Miller , “ Excitation of longitudinal and transverse coherent acoustic phonons in nanometer free-standing films of (001) Si,” Phys. Rev. B 79 (9 ), 094301 (2009).10.1103/PhysRevB.79.094301
18. L. L. Wei , S. S. Sun , C. Guo , Z. W. Li , K. Sun , Y. Liu , W. J. Lu , Y. P. Sun , H. F. Tian , H. X. Yang , and J. Q. Li , “ Dynamic diffraction effects and coherent breathing oscillations in ultrafast electron diffraction in layered 1T-TaSeTe,” Struct. Dyn. 4 (4 ), 044012 (2017).10.1063/1.4979643 28470025
19. A. Ungeheuer , A. S. Hassanien , M. T. Mir , L. Noding , T. Baumert , and A. Senftleben , “ Selective excitation of higher harmonic coherent acoustic phonons in a graphite nanofilm,” J. Phys. Chem. C 126 (46 ), 19822–19833 (2022).10.1021/acs.jpcc.2c06401
20. Q. K. Qian , X. Z. Shen , D. Luo , L. X. Jia , M. Kozina , R. K. Li , M. F. Lin , A. H. Reid , S. Weathersby , S. Park , J. Yang , Y. Zhou , K. Y. Zhang , X. J. Wang , and S. X. Huang , “ Coherent lattice wobbling and out-of-phase intensity oscillations of Friedel pairs observed by ultrafast electron diffraction,” ACS Nano 14 (7 ), 8449–8458 (2020).10.1021/acsnano.0c02643 32538617
21. E. M. Mannebach , C. Nyby , F. Ernst , Y. Zhou , J. Tolsma , Y. Li , M.-J. Sher , I. C. Tung , H. Zhou , Q. Zhang , K. L. Seyler , G. Clark , Y. Lin , D. Zhu , J. M. Glownia , M. E. Kozina , S. Song , S. Nelson , A. Mehta , Y. Yu , A. Pant , B. Aslan , A. Raja , Y. Guo , A. DiChiara , W. Mao , L. Cao , S. Tongay , J. Sun , D. J. Singh , T. F. Heinz , X. Xu , A. H. MacDonald , E. Reed , H. Wen , and A. M. Lindenberg , “ Dynamic optical tuning of interlayer interactions in the transition metal dichalcogenides,” Nano Lett. 17 (12 ), 7761–7766 (2017).10.1021/acs.nanolett.7b03955 29119791
22. X. Wang , J. Li , and J. Cao , “ Coherent phonon generation in laser-heated gold nanofilm,” J. Chem. Phys. 152 (12 ), 124704 (2020).10.1063/1.5137818 32241121
23. M. Chebl , X. He , and D.-S. Yang , “ Ultrafast carrier-coupled interlayer contraction, coherent intralayer motions, and phonon thermalization dynamics of black phosphorus,” Nano Lett. 22 (13 ), 5230–5235 (2022).10.1021/acs.nanolett.2c01019 35763556
24. S. K. Sundaram and E. Mazur , “ Inducing and probing non-thermal transitions in semiconductors using femtosecond laser pulses,” Nat. Mater. 1 (14 ), 217–224 (2002).10.1038/nmat767 12618781
25. A. Block , M. Liebel , R. Yu , M. Spector , and N. F. V. Hulst , “ Tracking ultrafast hot-electron diffusion in space and time by ultrafast thermo-modulation microscopy,” Sci. Adv. 5 (5 ), eaav8965 (2018).10.1126/sciadv.aav8965
26. L. Stojchevska , I. Vaskivskyi , T. Mertelj , P. Kusar , and D. J. Svetin , “ Ultrafast switching to a stable hidden quantum state in an electronic crystal,” Science 344 (6180 ), 177–180 (2014).10.1126/science.1241591 24723607
27. M. C. Lee , N. Sirica , S. W. Teitelbaum , A. Maznev , T. Pezeril , R. Tutchton , V. Krapivin , G. A. de la Pena , Y. Huang , L. X. Zhao , G. F. Chen , B. Xu , R. Yang , J. Shi , J. X. Zhu , D. A. Yarotski , X. G. Qiu , K. A. Nelson , M. Trigo , D. A Reis , and R. P. Prasankumar , “ Direct observation of coherent longitudinal and shear acoustic phonons in TaAs using ultrafast x-ray diffraction,” Phys. Rev. Lett. 128 (15 ), 155301 (2022).10.1103/PhysRevLett.128.155301 35499894
28. A. Zong , Q. Zhang , F. Zhou , Y. Su , K. Hwangbo , X. Shen , Q. Jiang , H. Liu , T. E. Gage , D. A. Walko , M. E. Kozina , D. Luo , A. H. Reid , J. Yang , S. Park , S. H. Lapidus , J.-H. Chu , I. Arslan , X. Wang , D. Xiao , X. Xu , N. Gedik , and H. Wen , “ Spin-mediated shear oscillators in a van der Waals antiferromagnet,” Nature 620 , 988–993 (2023).10.1038/s41586-023-06279-y 37532936
29. R. Merlin , “ Generating coherent THz phonons with light pulses,” Solid State Commun. 102 (2 ), 207–220 (1997).10.1016/S0038-1098(96)00721-1
30. Y. X. Yan and K. A. Nelson , “ Impulsive stimulated light-scattering.1. General-theory,” J. Chem. Phys. 87 (11 ), 6240–6256 (1987).10.1063/1.453733
31. H. J. Zeiger , J. Vidal , T. K. Cheng , E. P. Ippen , G. Dresselhaus , and M. S. Dresselhaus , “ Theory for displacive excitation of conherent phonons,” Phys. Rev. B 45 (2 ), 768–778 (1992).10.1103/PhysRevB.45.768
32. F. Vialla and N. Del Fatti , “ Time-domain investigations of coherent phonons in van der Waals thin films,” Nanomaterials 10 (12 ), 2543 (2020).10.3390/nano10122543 33348750
33. D. B. Williams and C. B. Carter , Transmission Electron Microscopy ( Springer, New York, NY, 2009).
34. Y. Su , A. Zong , A. Kogar , D. Lu , S. S. Hong , B. Freelon , T. Rohwer , H. Hwang , and N. Gedik , “ Delamination-assisted ultrafast wrinkle formation in a freestanding film,” arXiv:2212.12082 (2022).
35. D. Du , S. Manzo , C. Zhang , V. Saraswat , K. T. Genser , K. M. Rabe , P. M. Voyles , M. S. Arnold , and J. K. Kawasaki , “ Epitaxy, exfoliation, and strain-induced magnetism in rippled Heusler membranes,” Nat. Commun. 12 (1 ), 2494 (2021).10.1038/s41467-021-22784-y 33941781
