## Attosecond Probing of Coherent Vibrational Dynamics in CBr4

Jen-Hao Ou,  Diptarka Hait, 1 2,3 Patrick Rupprecht, 1,4 John E. Beetar, 1,† Todd J. Martí nez, 2,3 and Stephen R. Leone 1,4,5,*

1 Department of Chemistry, University of California, Berkeley, CA 94720, USA

2 Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States

3 SLAC National Accelerator Laboratory, Menlo Park, California 94024, United States

4 Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720,

USA

5 Department of Physics, University of California, Berkeley, CA 94720, USA

†

Present address: NSF National eXtreme Ultrafast Science (NeXUS) Facility, Columbus, OH

43210, USA

E-mail: srl@berkeley.edu

## ABSTRACT

A coherent vibrational wavepacket is launched and manipulated in the symmetric stretch (a1) mode of CBr4, by impulsive stimulated Raman scattering from non-resonant 400 nm laser pump  pulses  with  various  peak  intensities  on  the  order  of  tens  of  10 12 W/cm 2 .  Extreme ultraviolet (XUV)  attosecond transient absorption spectroscopy (ATAS) records the wavepacket dynamics as temporal oscillations in XUV absorption energy at the bromine M4,5 3d3/2,5/2 edges around 70 eV. The results are augmented by nuclear time-dependent Schrödinger equation  simulations.  Slopes  of  the  (Br-3d3/2,5/2) -1 10a1 * core-excited  state  potential  energy surface  (PES)  along  the  a1  mode  are  calculated  to  be  -9.4  eV/Å   from  restricted  open-shell Kohn-Sham calculations. Using analytical relations derived for the small-displacement limit with the calculated slopes of the core-excited state PES, a deeper insight into the vibrational dynamics is obtained by retrieving the experimental excursion amplitude of the vibrational wavepacket and the amount of population transferred to the vibrational first-excited state, as a function of pump-pulse peak intensity. Experimentally, the results show that XUV ATAS is capable of easily resolving oscillations in the XUV absorption energy on the order of few to tens of meV and tens of femtosecond time precision, limited only by the averaging times in the experimental scans. This corresponds to oscillations of C-Br bond length on the order of 10 -4 to 10 -3 Å . The results and the analytic relationships offer a clear physical picture, on multiple levels of understanding, for how the pump-pulse intensity controls the vibrational dynamics launched by non-resonant ISRS in the small-displacement limit.

## 1. INTRODUCTION

Preparing and controlling vibrational coherent superposition states (referred to as vibrational wavepackets)  in  molecules  with  light  is  a  major  theme  in  light-matter  interactions. 1,2 This necessitates control over  the creation, detection, and  manipulation  of  the vibrational wavepacket. 1,2 In  this  work (Figure 1), such control is achieved in the carbon tetrabromide (CBr4) molecule by attosecond transient absorption spectroscopy (ATAS). 3-5 We employ a 400 nm 26 fs pump laser pulse to generate a vibrational wavepacket in the symmetric stretch (a1) mode of the electronic ground state of CBr4 with non-resonant impulsive stimulated Raman scattering (ISRS). 6-8 The  resulting vibrational wavepacket  dynamics  are  subsequently monitored by time-resolved absorption of attosecond extreme ultraviolet (XUV) probe laser pulses at the Br M4,5 3d3/2,5/2 edge around 70 eV, recorded as temporal oscillations in XUV absorption energy.

Similar vibrational dynamics in various molecules have been experimentally observed via transient absorption spectroscopy probed at XUV or soft X-ray regions. 5,9-28 Vibrational wavepackets were triggered by pump pulses mostly in the infrared (780, 800, 1200 or 1600 nm), 9-21,25-28 while some in visible 22 or UV 23,24 spectral regions. Among them, pump pulses with high (~10 14 W/cm 2 ) peak intensity 9,10,12-20,27,28 often lead to large vibrational amplitudes. In particular, relatively few works 26,27 systematically measured how the vibrational dynamics change with pump-pulse peak intensity. In this work, 400 nm pulses with lower peak intensities (8.1-20.2×10 12 W/cm 2 ) are employed, and the pump-pulse intensity dependence of vibrational dynamics is investigated in the small-displacement limit. 26

To gain greater understanding, it is valuable to explicitly understand how the temporal oscillations in XUV absorption energy connect the shape of the corresponding core-excited state potential energy surface (PES) with aspects of the vibrational wavepackets, such as the excursion  amplitude  of  the  wavepacket  and  populations  in  the  vibrational  eigenstates.  To clearly establish these connections, we derive analytical relations for the above quantities in the small-displacement limit with the harmonic approximation. We compute the slope of the Br 3d core-excited state PES along the symmetric stretch mode at the ground state equilibrium geometry from restricted open-shell Kohn-Sham (ROKS) calculations. 29-31 This slope is used to retrieve the excursion amplitude of the vibrational wavepacket (in terms of C-Br bond-length displacement) and the population transferred by ISRS to the vibrational first-excited state in the small-displacement limit, from the measured time-resolved XUV absorption spectrum.

Experimentally, the results show that XUV ATAS is capable of resolving small bondlength  changes  during  vibrational  wavepacket  dynamics  to  the  order  of  10 -4 Å with femtosecond  temporal  resolution.  With  analytical  formulations,  we  offer  a  clear  physical picture on multiple levels of understanding for how the pump-pulse peak intensity controls vibrational dynamics launched by non-resonant ISRS in the small-displacement limit.

Figure  1. Experimental  Scheme.  (a)  Attosecond  transient  absorption  spectroscopy  on  the symmetric stretch (a1) mode of CBr4. (b) Vibrational eigenstates are coherently excited by nonresonant impulsive stimulated Raman scattering (ISRS) with a 400 nm 26 fs pump pulse, and a vibrational wavepacket is launched on the electronic ground state potential energy surface.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-75d31422-ecfe-435f-959b-c2fcee95bec0.jpg)

(c) Vibrational wavepacket dynamics are monitored by time-resolved absorption of attosecond XUV probe pulses from the electronic ground state to the (Br-3d3/2,5/2) -1 10a1 * core-excited states around 70 eV. The wavepackets drawn in green and red represent the motion shifting toward left and right, respectively.

## 2. METHODS

## 2.1 Experiment

The scheme of attosecond transient absorption spectroscopy on CBr4, the 400 nm pump pulse, and the XUV probe pulse, are briefly outlined here and in  Figure 1. The laser source is a Coherent Astrella USP (800 nm, 35 fs FWHM, 7 mJ, 1 kHz), and its output is split by a 70:30 beamsplitter into pump and probe arms.

## 2.1.1 400 nm Pump Pulse

The pump arm starts with 70% of the Astrella output (800 nm, 35 fs FWHM, 4.9 mJ) and undergoes second harmonic generation (SHG) in a β-barium borate (BBO) crystal and the central wavelength becomes 400 nm. Then it is free-space focused through a Ne-filled (1.7 bar) tube (1.37 m long) to gently broaden the spectrum, and then compressed by chirped mirrors. The pulse duration of the main peak in the intensity envelope is measured to be 26±1 fs FWHM by homebuilt Self-Diffraction Frequency Resolved Optical Gating (SD-FROG) (Figure S6, Supporting Information).

## 2.1.2 Few-Cycle Visible-NIR Pulse

The probe arm starts with 30% of the Astrella output (800 nm, 35 fs FWHM, 2.1 mJ) and is focused into a Ne-filled (2.4 bar) hollow core fiber (HCF) (1.5 m long, 400 µm inner diameter), where it undergoes supercontinuum generation, and then the pulses are compressed by chirped mirrors.  After  the  chirped  mirrors,  the  spectrum  spans  from  around  500  to  900  nm.  The dispersion is further fine tuned by a piece of 2 mm thick ADP crystal and a pair of thin glass wedges. The output is few-cycle (~ 4 fs) visible-NIR pulses that are used as the driving field for the high harmonic generation (HHG) into XUV region.

## 2.1.3 XUV Probe Pulse

The few-cycle visible-NIR pulses (driving field) are focused into a Ne gas cell (~80 mmHg, 300 µm hole diameter) and an XUV probe attosecond pulse train (~50-85 eV photon energy) is  produced  by  high  harmonic  generation  (HHG).  The  remaining  800  nm  driving  field  is removed by an Al filter (0.1 µm thick).

## 2.1.4 Static and Time-Resolved (Transient) XUV Absorption Measurement

At the sample cell (4 mm long, 500 µm hole diameter), 400 nm pump and XUV probe beams recombine non-collinearly and interact with CBr4 vapor. The sample reservoir, delivery line, and cell are heated to 100-120 ℃ to obtain enough vapor pressure from solid CBr4. The time delay between the pump and probe pulses is controlled by changing the optical path length between the pump and probe beamlines with a translational stage. After a Zr filter (0.1 µm thick), the XUV beam is dispersed by a grating and detected on a XUV CCD camera.

Two types of XUV absorption spectra are collected experimentally. One is the static XUV absorbance XUV ( ) , A E

$$A ( E _ { \text{XUV} } ) = - \log _ { 1 0 } \frac { I _ { \text{XUV} } ( E _ { \text{XUV} } ) } { I _ { \text{XUV} } ^ { ( 0 ) } ( E _ { \text{XUV} } ) }$$

where E XUV  is  the  photon  energy  of  XUV  probe  pulses,  and XUV I and (0) XUV I are  the  XUV intensity with and without the CBr4 molecule in the cell, respectively. The other is the timeresolved (transient) XUV differential absorption XUV ( , ) A E   ,

$$\Delta A ( E _ { X U V }, \tau ) = - \log _ { 1 0 } \frac { I _ { X U V + p u p m p } ( E _ { X U V }, \tau ) } { I _ { X U V } ( E _ { X U V }, \tau ) }$$

where   is the time delay between the pump and probe pulses, and τ XUV+pump I and XUV I are the XUV spectrum of CBr4 with and without the pump pulse.

Temporal resolution is calibrated with a cross-correlation measurement of pump and probe pulses on the He 2s2p state. The fitted instrument response function is 26±6 fs and the time zero is -1±3 fs. The spectral axis is calibrated with a static XUV absorption spectrum of He 2snp (n ≥ 2) states. Details are in the Supporting Information.

## 2.1.5 Pump-Pulse-Intensity Dependent Measurement

The peak intensity of the 400 nm pump pulse is modified by cropping the beam with an iris, before the pump and probe pulses recombine. The peak intensity at the sample cell location is estimated  to  range  from  8.1  to  20.2×10 12 W/cm 2 .  Detailed  calibrations  are  given  in  the Supporting Information.

## 2.2 Calculation

Quantum chemical calculations are performed with the Q-Chem software package, 32 utilizing the  SCAN0 functional 33 and  the  decontracted  (i.e.  all  primitives  are  uncontracted)  aug-ccpVTZ basis set. 34-36 Scalar relativistic effects are included through the use of the spin-free exact two component one-electron (SF-X2C-1e) approach. 37-40 Local exchange-correlation integrals are calculated by quadrature over a radial grid with 250 points and an angular Lebedev grid with 974 points. The electronic ground state is modeled with the standard restricted KohnSham procedure, 41 yielding  an  equilibrium bond distance of 1.925 Å  (vs an experimental 42 value  of  of  1.942  Å )  and  a  symmetric  stretch  frequency  of  282  cm -1 (vs  267  cm -1 from experiment 43 ).

The isotropic polarizability of the molecule is computed via finite differences of the energy  with  a  five-point  stencil  (using  electric  field  strengths  of  0.005  and  0.01  a.u.,  see Supporting Information) over bond lengths ranging from 1.8246-2.046 Å  (in increments of

0.001 Å ). These polarizabilities are subsequently utilized to carry out nuclear time-dependent Schrödinger equation simulations of the ISRS process induced by the pump laser pulse and the resulting nuclear dynamics along the symmetric stretching mode, as described in detail in the Supporting Information.

The Br 3d excitation energies are computed with the ROKS approach, 29,30 which has been found to be quite effective at modeling core-level spectroscopies. 31,44 Excited state orbital optimization is carried out with the square gradient minimization algorithm. 45 Spin-orbit effects are  subsequently  incorporated  in  the  manner  described  previously 31 .  This  quasi-degenerate perturbation-theory-based  approach  constructs  a  zeroth-order  effective  Hamiltonian  of  the spin-orbit  free  3d  excitation  energies  at  a  given  geometry,  to  which  a  spin-orbit  coupling operator -𝐽 𝐿   ⋅ 𝑆 ⃗ is  added  (where 𝑆 is  the  electron  spin  operator  and 𝐿 ⃗ the  orbital  angular momentum  for  d  electrons).  The  coupling  constant 𝐽 is  assumed  to  be  constant  for  all geometries and chosen to be 0.40 eV since this value reproduces the experimentally observed multiplet  splitting  of  1.03  eV  between  the  3d5/2  and  3d3/2  absorption  energies.    Core-level excitation  energies  are  computed  over  bond  lengths  ranging  from  1.9046-1.9446  Å ,  in increments of 0.001 Å , and are fitted to a straight line against the bond length (see Figure 3 ) in order to determine the slope of the core-excited PES vs bond elongation.

## 3. RESULTS AND DISCUSSION

## 3.1  Detection  of  the  Coherent  Vibrational  Dynamics:  Static  and  Time-Resolved  XUV Absorption Spectrum of CBr4 at Br M4,5 3d3/2,5/2 Edge

Figure 2 (left) shows the static XUV absorption spectrum of CBr4 from the electronic ground state  to  the  (Br-3d3/2,5/2) -1 10a1 * core-excited  states  (10a1 * is  the  lowest  energy  antibonding orbital 46,47 ), and the spin-orbit doublet is resolved by fitting to a sum of two Voigt functions. The extracted peak energies are 70.932±0.011 and 69.903±0.008 eV for the Br M4 3d3/2 and

M5 3d5/2 edges, respectively (three decimal places are included here to compare with the very small (few to tens of meV) oscillation amplitude of XUV absorption energy measured and discussed below). The spin-orbit splitting is 1.03±0.01 eV (Table 1). They agree reasonably with  previous  experimental  results 46 (70.5  and  69.6  eV).  ROKS  calculations  yield  slightly lower excitation energy values of 70.05 and 69.02 eV, indicating a somewhat larger difference between theory and the current experiment (~0.9 eV) than the ~0.3 eV deviations previously observed for the K- and L- edges of lighter elements. 31,44

Next, we consider changes to the XUV absorption due to the vibrational wavepacket launched  by  400  nm  pump  pulse  (peak  intensity  20.2×10 12 W/cm 2 ).  As  the  vibrational wavepacket oscillates on the electronic ground state PES, the XUV absorption peak oscillates correspondingly with the vibrational period, as shown in the time-resolved XUV absorption spectrum in Figure 2 (right). The fitted periods are 123±1 fs for both the Br M4 3d3/2 edge and M5 3d5/2 edge, which agrees with the literature value of 126 fs for the symmetric stretch mode (267 cm -1 ). 43 The oscillations in the peak position for the two edges are almost temporally in phase with a 0.1±0.4 rad phase difference. Oscillatory behaviors corresponding to the other normal modes of the molecule are not observed (despite such modes being Raman active) because there is no net displacement along such asymmetric modes, as has been observed and discussed in a previous work for CCl4. 28

If the launching mechanism of the vibrational wavepacket is a pure non-resonant ISRS from  a  Gaussian  pump  pulse  centered  at  time  zero,  the  expectation  value  of  wavepacket position  will  be  a  sine  function  with  a  zero  initial  phase. 6,8 Nuclear  TDSE  simulations reproduce this expectation by retrieving around 0.01 π initial phase of wavepacket position for all  intensities,  using  the  temporal  intensity  envelope  of the pump pulse  retrieved from SDFROG (Figure S6). Experimentally, over the range of intensities measured, the initial phase of XUV absorption  energy  varies  from  -0.1  to  -0.6 π  ( Figure  S12,  Supporting  Information), suggesting a possible deviation from the pure non-resonant ISRS. Potential causes may be the pre-pulse  (Figure  S6)  or  multiphoton  absorption  observed  in  the  time-resolved  XUV differential absorption spectrum (Figure S11).

Figure 2. (Left) Static XUV absorption spectrum of CBr4 at Br M4,5 3d3/2,5/2 edges from the electronic ground state to the (Br-3d3/2,5/2) -1 10a1 * core-excited states. The spin-orbit doublet is resolved  by  fitting  to  a  sum  of  two  Voigt  functions.  The  extracted  peak  energies  are 70.932±0.011 and 69.903±0.008 eV for the Br M4 3d3/2 and M5 3d5/2 edges respectively, with a  1.03±0.01  eV  spin-orbit  splitting  resulting  from  the  fit.  (Right)  Time-resolved  XUV absorption spectrum of CBr4 as a function of time delay between the 400 nm pump and XUV probe pulses. The peak intensity of pump pulse is 20.2×10 12 W/cm 2 . Color shade: extracted peak energy with error bars at each time delay. Solid line: a cosine wave fitted to the extracted peak energies. The fitted period is 123±1 fs, amplitude 10.2±1.4 meV, center 70.934±0.001 eV (black dashed line) for the M4 3d3/2 edge, while the M5 3d5/2 edge yields 123±1 fs, 11.3±1.4 meV, and 69.903±0.001 eV (black dashed line). The oscillations for the two edges are almost temporally in phase with a 0.1±0.4 rad phase difference.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-f5fa2f1b-0e07-4c22-85de-118bffe19017.jpg)

Table 1. Excitation energies (eV) of (Br-3d3/2,5/2) -1 10a1 * core-excited states in CBr4.

| Core-Excited State      | Experiment (This work)   |   ROKS Calculations (This work) |   Previous Experiment 46 |
|-------------------------|--------------------------|---------------------------------|--------------------------|
| (Br-3d 3/2 ) -1 10a 1 * | 70.932±0.011             |                           70.05 |                     70.5 |
| (Br-3d 5/2 ) -1 10a 1 * | 69.903±0.008             |                           69.02 |                     69.6 |
| Spin-Orbit Splitting    | 1.03±0.01                |                            1.03 |                      0.9 |

## 3.2 Information Encoded in the Time-Resolved XUV Absorption Spectrum: Oscillation Amplitude of XUV Absorption Energy, Excursion Amplitude of Vibrational Wavepacket, and Slope of Core-Excited State Potential Energy Surface

When the vibrational  wavepacket  motion  results  in  a  normal  mode  displacement q on  the electronic ground state PES E g( q ),  the  XUV absorption energy E XUV( q )  from the electronic ground state to the core-excited state PES E ce ( q )  is ( ) = ( ) -XUV ce g ( ) .  If  the  pumpE q E q E q pulse intensity is not too strong to induce a large displacement q , the harmonic approximation holds and the XUV absorption energy can be approximated by a Taylor series in q

$$E _ { x u v } ( q ) & \approx \left [ E _ { c e } ( 0 ) - E _ { g } ( 0 ) \right ] + \left [ \left ( \frac { d E _ { c e } } { d q } \right ) _ { 0 } - \left ( \frac { d E _ { g } } { \not d q } \right ) _ { 0 } \right ] q \\ & = E _ { x u v } ( 0 ) + \left ( \frac { d E _ { c e } } { d q } \right ) _ { 0 } q$$

where  the  slope  of  the  electronic  ground  state  PES E g( q )  vanishes  at  the  ground  state equilibrium position q = 0, while the slope of the electronic core-excited state PES E ce ( q ) will generally  be  nonzero  at  the  Frank-Condon  region,  as  the  equilibrium  position  of  the  coreexcited PES is quite unlikely to coincide with that of the ground state.

As the vibrational wavepacket travels between the outer turning point q ot to the inner turning point q it on the electronic ground state PES, the range of peak XUV absorption energies that can be accessed is

$$| E _ { \text{XXV} } ( q _ { \text{ot} } ) - E _ { \text{XXV} } ( q _ { \text{it} } ) | = \left | \left ( \frac { d E _ { \text{ce} } } { d q } \right ) _ { 0 } \right | ( q _ { \text{ot} } - q _ { \text{it} } )$$

where XUV XUV ot XUV it ( ) ( ) / 2 E E q E q   -is  defined  as  the oscillation  amplitude  of  the  XUV absorption  energy , and amp ot it ( ) / 2 q q q  -, the excursion  amplitude  of  the  vibrational wavepacket , or physically the maximum displacement from equilibrium. Therefore,

$$\left | \Delta E _ { \text{XUV} } = \left | \left ( \frac { d E _ { \text{ce} } } { d q } \right ) _ { 0 } \right | q _ { \text{amp} } \right |$$

and amp q can thus be retrieved from XUV E  in the experimentally measured time-resolved XUV absorption spectrum in Figure 2 (right), if the slope of the core-excited state PES along the normal mode at the electronic ground state equilibrium geometry is known.

In this work, the slope of the (Br-3d3/2,5/2) -1 10a1 * core-excited states PES is calculated from  ROKS  to  be  around  -9.4  eV/Å   along  the  symmetric  stretch  mode  (by  fitting  coreexcitation energies to a straight line, as shown in Figure 3 ). In other words, this slope suggests that  the  Br-3d3/2,5/2→10a1 * excitation  energy  decreases  by  9.4  eV  as  the  C-Br  bond  length increases by 1 Å . This is consistent with the antibonding character of the 10a1 * orbital, as the energy of this state should significantly decrease with bond elongation. This 9.4 eV/Å   value for the slope is in the same order of magnitude (~10 eV/Å ) as previous estimates for core-to- antibonding excitations in CCl4 28 and SF6 26,27 (typical range 5 - 30 eV/Å  depending on which atom).

The ROKS results and the resulting linear fits are shown in Figure 3 , which indicates that Equation 3 holds quite well for bond stretches as high as ±0.05 Å  relative to equilibrium and higher order effects do not appear to be relevant in this regime. The ~10 meV or smaller experimental oscillation amplitudes for XUV absorption energies (Figure 2) and the computed slope -9.4 eV/Å  collectively indicate that the bond stretches remain on the order of 10 -4 Å , and the linear behavior reported in Equation 3 therefore is sufficient for analyzing the experimental results.

Figure 3. Computed XUV excitation energies from ROKS for various symmetric C-Br bond stretches (open circles) around the equilibrium geometry (1.925 Å ). The linear fits obtained from fitting to ROKS results for 41 equally spaced bond lengths between 1.9046-1.9446 Å  are also shown (solid lines), which align quite well with ROKS results for twelve bond lengths calculated outside the fit interval (six equally spaced points each between 1.8746-1.8996 Å  and 1.9496-1.9746 Å ). The slopes are about -9.4 eV/Å  for the Br-3d3/2,5/2→10a1 * excitations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-189e4d4c-c99d-4f9b-99e0-a33cdd697058.jpg)

3.3 Relative Slope between Multiple Core-Excited State PES and Insight into Property of Participating Orbitals

If multiple electronic core-excited states simultaneously probe the vibrational wavepacket on the electronic ground state, the relative slope between two or more core-excited state PES 26 -28 can be determined experimentally from the ratio of measured oscillation amplitude of XUV absorption energy E  as

XUV

$$\Delta E _ { \text{XUV} } ^ { \text{ce2} } / \Delta E _ { \text{XUV} } ^ { \text{ce1} } = \left | \left ( \frac { d E _ { \text{ce2} } } { d q } \right ) _ { 0 } \right | / \left | \left ( \frac { d E _ { \text{ce1} } } { d q } \right ) _ { 0 } \right |$$

because here the excursion amplitude of the vibrational wavepacket q amp is the same. In this work,  the  relative  slope  between  the  two  (Br-3d 3/2,5/2 ) -1 10a1 * core-excited  states  PES  is determined experimentally to be 1.1±0.13 as  the  ratio  of  measured oscillation amplitude of XUV  absorption energy XUV E  (Br-3d5/2)/ XUV E  (Br-3d3/2) (Table S4, Supporting Information). As a result, the shape of the two core-excited state PESs are quite similar around the equilibrium geometry of the electronic ground state.

Such  relative  slopes  not  only  offer  experimental  information  on  the  shape  of  coreexcited state PES, but also shed light on the nature of participating orbitals. For the discussion below, a core-excited state is approximated as a singly-excited electronic configuration, where one electron is excited from an initial core orbital to a final orbital that is higher in energy.

One scenario compares cases when the electron is excited from the same core orbital to  different  final  orbitals.  The  relative  slopes  of  the  corresponding  core-excited  state  PESs belonging to the same absorption edge gives a hint about the bonding properties of the final orbitals. For instance, if the final orbital is non-bonding, then the resulting core-excited state will  be  generally  less  sensitive  to  changes  in  the  molecular  geometry  and  its  PES  will  be relatively flat. As a result, its slope will be smaller than the slope of other core-excited states PES where the final orbital is bonding or anti-bonding. For instance, in the SF6 molecule 26,27 , an electron can be excited from S-2p core orbital to three different final orbitals a1g, t2g and eg.

Because the a1g and eg orbitals are anti-bonding while the t2g orbital is non-bonding, the slope of (S-2p) -1 a1g and (S-2p) -1 eg core-excited states is found to be larger than that of (S-2p) -1 t 2g state.

The other scenario is when the electron is excited from different core orbitals to the same  final  orbital,  or  the  corresponding  core-excited  states  originating  from  different absorption edges. This offers a potential way to compare different core orbitals on the same atom, or core orbitals on different atoms in a molecule. In this work, the electron is excited to the  same  10a1 * anti-bonding  orbital,  but  from  different  Br  3d3/2  or  Br  3d5/2  core  orbitals, corresponding to the Br M4 or M5 absorption edges. The ratio of slopes of these two coreexcited states is close to unity as mentioned above, indicating that the both Br 3d3/2 or Br 3d5/2 core orbitals respond similarly to the bond-length change. This is reasonable as these two core orbitals come from the same set of Br 3d orbitals, split by spin-orbit coupling. On the other hand, a previous study on CCl4 molecule 28 found the slope for the (Cl-2p3/2) -1 8t2 * core-excited state to be much larger than that of the (C-1s) -1 8t2 * core-excited state, despite the same final orbital (8t2 * ) for both excitations.

## 3.4  Connection  between  the  Excursion  Amplitude  of  Vibrational  Wavepacket  and Vibrational Eigenstate Populations, in the Small-Displacement Limit

In the small-displacement limit, only the vibrational ground ψ 0 and first-excited ψ 1 states of a harmonic oscillator need to be considered. 26 Long after the pump pulse, the coherent vibrational wavepacket Ψ( q , t ) can be approximated by the following linear superposition,

$$\Psi ( q, t ) = c _ { 0 } \psi _ { 0 } ( q ) e ^ { - i E _ { 0 } t / h } + c _ { 1 } \psi _ { 1 } ( q ) e ^ { - i E _ { t } t / h }$$

where q is the normal mode displacement (position), and E 0,1 is the vibrational energy of the vibrational  ground  and  first-excited  states,  respectively.  The  expectation  value  of  position ( ) q t of this vibrational wavepacket is:

$$& \langle q ( t ) \rangle \\ & = \langle \Psi ( q, t ) | q | \Psi ( q, t ) \rangle \\ & = p _ { 0 } \langle \psi _ { 0 } | q | \psi _ { 0 } \rangle \\ & \, + p _ { 1 } \langle \psi _ { 1 } | q | \psi _ { 1 } \rangle \\ & \, + 2 \sqrt { p _ { 0 } p _ { 1 } } \langle \psi _ { 0 } | q | \psi _ { 1 } \rangle \cos [ ( E _ { 1 } - E _ { 0 } ) t / \hbar { + } \phi ]$$

where 2 0,1 0,1 p c = is the population  of  the  vibrational  ground  and  first-excited  states, respectively, and ϕ is  the  initial  phase  between  the  two  coefficients 𝑐0,1 .  (If  the  launching mechanism of the  coherent  vibrational  wavepacket  is  non-resonant  ISRS  from  a  Gaussian pump  pulse  centered  at  time  zero,  the  expectation  value  of  position ( ) q t will  be  a  sine function with a zero initial phase, 6,8 or equivalently a cosine function with a π/2 initial phase in Equation 8.)

$$\text{Next,    because    } \left \langle \psi _ { 0 } | q | \psi _ { 0 } \right \rangle = \left \langle \psi _ { 1 } | q | \psi _ { 1 } \right \rangle = 0 \text{    and   } \left \langle \psi _ { 0 } | q | \psi _ { 1 } \right \rangle = \sqrt { \frac { h } { 2 \mu \omega _ { 0 } } } \text{    for                                                                                                                                                                                                        
                                                                                                                                                                                                        }$$

oscillator (where μ is the mass associated with the normal mode, which is four times the mass of Br, and 0 1 0 ( ) / E E  = -is the fundamental vibrational frequency of the normal mode), the position expectation value simplifies to:

$$\langle q ( t ) \rangle & = \left ( 2 \sqrt { p _ { 0 } p _ { 1 } } \sqrt { \frac { \hbar } { 2 \mu \omega _ { 0 } } } \right ) \cos ( \omega _ { 0 } t + \phi ) \\ & = \left ( 2 \sqrt { ( 1 - p _ { 1 } ) p _ { 1 } } \sqrt { \frac { \hbar } { 2 \mu \omega _ { 0 } } } \right ) \cos ( \omega _ { 0 } t + \phi ) \\ & = q _ { \text{amp} } \cos ( \omega _ { 0 } t + \phi )$$

as 𝑝1 +𝑝0 = 1. This shows that the excursion amplitude of the vibrational wavepacket is

$$q _ { \text{amp} } = 2 \sqrt { \left ( 1 - p _ { 1 } \right ) p _ { 1 } } \, \sqrt { \frac { h } { 2 \mu \omega _ { 0 } } }$$

Furthermore, when 𝑝1 ≪1 (as expected in the small-displacement limit),

$$\boxed { q _ { a m p } \approx 2 \sqrt { p _ { 1 } } \sqrt { \frac { \hbar } { 2 \mu \omega _ { 0 } } } }$$

so q amp  is  proportional  to 1 p .  Therefore,  if  the  excursion  amplitude  of  the  vibrational wavepacket is known, the population in the vibrational first-excited state can be determined from this relation.

The validity of the above formulation is also verified against numerical simulations of the nuclear time-dependent Schrödinger equation (TDSE) along the symmetric stretch mode alone. As shown in Figure 4, the calculated population of the vibrational ground state is close to one (on the order of 0.99), and the calculated population of the vibrational first-excited state is indeed small (around the order of 10 -3 ) for the pump pulse intensities used in this work. As a  result,  the  pump  light-matter  interaction  studied  in  this  work  is  clearly  in  the  smalldisplacement limit.

Figure 4. Calculated populations of the vibrational eigenstates (ν = 0, 1, 2) of the symmetric stretch a1 mode from nuclear TDSE simulations, as a function of the experimental 400 nm pump-pulse  intensity.  The  population  of  the  vibrational  first-excited  state  remains  small (around the order of 10 -3 ). This shows that the pump light-matter interaction is in the small-

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-e212a73c-3dfa-40a8-8eb9-b68dc77deda2.jpg)

displacement limit for the range of pump-pulse intensities used in this work. The population in the second vibrational excited state is even smaller, indicating that the two-state approximation used for the small-displacement limit is reasonable.

## 3.5  Physical  Basis  for  the  Phenomenological  Equations  in  the  Literature  to  Fit  the Oscillations in Time-Resolved XUV Absorption Spectrum due to Coherent Vibrational Dynamics

Based on the formulation above, a physical basis can be provided for the phenomenological equations that have often been used in the literature 9,12,14,23,24,27 to fit the oscillations in timeresolved  XUV  absorption  spectrum  due  to  coherent  vibrational  dynamics.  If  only  one vibrational  mode  is  involved,  the  time-resolved  XUV  absorption  energy  is  often  described empirically as

$$E _ { \text{XUV} } ( t ) = E _ { \text{XUV} } ^ { 0 } + \Delta E _ { \text{XUV} } \cos ( \omega _ { 0 } t + \phi )$$

where 0 XUV E is  the XUV absorption energy at the equilibrium geometry of electronic ground state, and XUV E  is the oscillation amplitude of XUV absorption energy. According to Equation 5,

$$\Delta E _ { \text{XUV} } = \left | \left ( \frac { d E _ { \text{ce} } } { d q } \right ) _ { 0 } \right | q _ { \text{amp} }$$

it further reveals that oscillation amplitude of XUV absorption energy is proportional to the slope  of  core-excited  state  PES  along  the  normal  mode  at  the  electronic  ground  state equilibrium  geometry  and  the  excursion  amplitude  of  vibrational  wavepacket,  within  the harmonic approximation. Larger amplitudes thus can result from steeper slopes of the coreexcited PES (which can be attained through excitation to an antibonding level, as in this work), or from higher pump-pulse peak intensity to access greater excursion amplitudes of vibrational wavepacket q amp, by adding a greater fraction of the vibrational first-excited state in the smalldisplacement limit.

## 3.6 Controlling the Motion of Vibrational Wavepacket with Different Pump-Pulse Peak Intensities in the Small-Displacement Regime

The dynamics of the vibrational wavepackets are initiated with five different pump-pulse peak intensities ranging from 8.1 to 20.2×10 12 W/cm 2 (detailed calibrations provided in Table S1). Figure 5(a) shows that the experimentally measured oscillation amplitudes of XUV absorption energy increase approximately linearly with pump-pulse intensity, which is also observed for the SF6 molecule 26,27 .

In  Figure  5(b),  the  excursion  amplitudes  of  the  vibrational  wavepacket q amp  are retrieved from the corresponding oscillation amplitudes of the XUV absorption energy XUV E  by Equation 5, given the -9.4 eV/Å  slope of the (Br-3d3/2,5/2) -1 10a1 * core-excited state PES from ROKS calculations. In Figure 5(c), the square root of population in the vibrational first-excited state 1 p is further estimated from the excursion amplitude of the vibrational wavepacket q amp via Equation 11. All three quantities plotted in Figure 5 are found to increase approximately linearly with the pump-pulse intensity for both the Br M4  3d3/2 and M5  3d5/2 edges.

Figure 5. Pump-pulse peak-intensity dependence of (a) experimentally measured oscillation amplitude of the XUV absorption energy XUV E  , (b) excursion amplitude of the vibrational wavepacket q amp  (retrieved  from  (a)  by  Equation  5,  given  the  -9.4  eV/Å   slope  of  the  (Br3d3/2,5/2) -1 10a1 * core-excited state PES from ROKS calculations), and (c) square root of the population  transferred  to  the  vibrational  first-excited  state 1 p (estimated  from  (b)  by Equation  11).  All  three  quantities  increase  approximately  linearly  with  the  pump-pulse intensity for both the Br M4 3d3/2 (orange circle) and M5 3d5/2 (blue square) edges. Estimates of all three quantities from nuclear TDSE simulations (black open square) for the ISRS process are also provided (scaled by 1/3, see text for details), and show clear linear relationships with the pump-pulse intensities used in this work. The outlier at 10.3×10 12 W/cm 2 may result from the fluctuation in average power measurement, as the average power is close to the lower limit of the working range for the power meter head (Coherent PowerMax PM10).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-9aa52aa6-312b-444f-9cd6-a86f538767d4.jpg)

The  experimentally  measured  and  retrieved  values  are  also  compared  with  nuclear TDSE  simulations.  The  theoretical  population  of  the  vibrational  first-excited  state p 1  and excursion amplitude of the vibrational wavepacket q amp come from numerically solving the nuclear TDSE for the non-resonant ISRS process, and the theoretical oscillation amplitude of XUV absorption energy XUV E  is estimated from Equation 5. The TDSE simulation results have  to  be  scaled  by  1/3  to  match  the  experiment.  A  few  possible  reasons  for  this  large discrepancy between experiment and theory are considered below. We note that the TDSE simulations account for rotational averaging effects on the interaction between the molecule and the pump pulse field through the use of isotropic polarizability.

From an experimental perspective, the differences with theory may arise from a lower effective  pump-pulse  peak  intensity  at  the  sample  location  than  the  estimation.  A  detailed calibration of the experimental peak intensity for the pump pulse is given in the Supporting Information, which is taken as our best estimation. However, it is uncertain whether the pump and probe pulses spatially overlap right at the beam center of each other. If the spatial overlap is off-center, the effective pump-pulse peak intensity within the probe beam area will become smaller. Before the measurements on CBr4 were made, the spatial overlap was confirmed by overlapping the recombined pump beam and the HHG driving-field beam on a beam profiler at a location with equivalent distance to the sample cell, as the XUV probe beam most likely follows the same path of the HHG driving-field beam. However, if there is an angular mismatch, the overlap could be diminished and therefore be a cause of the discrepancy. On the other hand, the outlier at 10.3×10 12 W/cm 2 may result from the fluctuations in average power measurement, as the average power is close to the lower limit of the working range for the power meter head (Coherent PowerMax PM10).

Potential computational origins of this discrepancy are errors in the Br 3d core-excited PES slope value of -9.4 eV/Å , or inadequacies in modeling the non-resonant ISRS process through TDSE (such as from errors in the molecular polarizabilities). We consider these factors to  be  somewhat  less  likely  as  the  slope  is  quite  consistent  with  previous  results  for  other molecules, and additional effort has been invested into benchmarking molecular polarizabilities used for the TDSE simulations (see Supporting Information).

## 4. CONCLUSIONS

In this work, a coherent vibrational wavepacket is launched in the symmetric stretching (a1) mode of electronic ground state CBr4 molecule by non-resonant impulsive stimulated Raman scattering  (ISRS)  with  400  nm  pump  pulses.  Dynamics  of  this  vibrational  wavepacket  are studied with time-resolved absorption of attosecond XUV pulses at the bromine M4,5 3d3/2,5/2 edges around 70 eV, with excitations to the lowest energy antibonding (10a1 * )  orbital.  The measured  XUV  absorption  energies  oscillate  in  time  with  the  period  of  the  symmetric stretching mode (126 fs), and the oscillation amplitude of XUV absorption energy equals the product of the excursion amplitude of the vibrational wavepacket and the slope of the Br (Br3d3/2,5/2) -1 10a1 * core-excited state potential energy surface along this mode. Since the measured oscillation amplitude of XUV absorption energy is on the order of few to tens of meV and the slope  is  -9.4  eV/Å   from  ROKS  calculations,  the  excursion  amplitude  of  the  vibrational wavepacket (C-Br bond-length displacement) is retrieved to be on the order of 10 -4 to 10 -3 Å . Harmonic approximation therefore holds quite well in this small-displacement limit.

Experiments with various pump-pulse peak intensities (8.1-20.2×10 12 W/cm 2 ) show an approximately linear relationship between the intensity and the measured oscillation amplitude of XUV absorption energy, and between the intensity and the retrieved excursion amplitude of vibrational wavepacket  as well. Furthermore, the excursion amplitude of vibrational wavepacket is proportional to the square root of the population in the vibrational first-excited state  in  the  small-displacement  limit.  This  relation  enables  an  indirect  estimate  of  this population from the measured oscillation amplitude of XUV absorption energy. These linear relationships are supported by numerical simulations of the nuclear time-dependent Schrodinger equation along the symmetric stretching mode for the non-resonant ISRS.

The  present  analysis  and  approach  to  extract  information  encoded  in  time-resolved XUV absorption spectra is general and can be readily applied to ISRS combined with timeresolved absorption spectroscopy probed in other spectral regions. Experimentally, the results show that XUV attosecond transient absorption spectroscopy is capable of resolving the small bond-length changes during vibration with extraordinary 10 -4 to  10 -3 Angstrom and tens of femtosecond precision. 26 Moreover, the analytical expressions provide a concrete and thorough understanding of how the pump-pulse intensity controls vibrational dynamics in the smalldisplacement limit.

## ACKNOWLEDGMENT

J.H.O., J.E.B., and S.R.L. thank the National Science Foundation under grant CHE-2243756 for  support  of  the  experimental  work.  J.H.O.  thanks  the  Government  Scholarship  to  Study Abroad, Ministry of Education, Taiwan for supplemental support, and Stefan Schippers and Christian Schroeder for discussions on fitting Fano-Gaussian lineshapes. D.H. is a Stanford Science Fellow. is. P.R. acknowledges funding by the Alexander von Humboldt Foundation (Feodor-Lynen Fellowship) and support of the AMOS program of the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. T.J.M and D.H. were supported by the AMOS program of the Office of Science of the U. S. Department of Energy. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231 using NERSC  award BESERCAP0027716.

## SUPPORTING INFORMATION AVAILABLE

- 1.  Calibration  of  Pump  Pulse  (Beam  diameter,  average  power,  peak  intensity,  and  SelfDiffraction Frequency Resolved Optical Gating (SD-FROG) measurement)
- 2. Calibration of Time Zero and Temporal Resolution for Time-Resolved XUV Absorption

## Measurements

- 3. Calibration of Spectral Axis
- 4. Measured Relative Slope of Br-3d Core-Excited State Potential Energy Surface (PES)
- 5. Comment about Initial Phase of Wavepacket Trajectory
- 6. Nuclear Time-Dependent Schrodinger Equation Simulations
- 7. Comment about Polarizability Calculations

## REFERENCES

- (1) Ohmori, K. Wave-Packet and Coherent Control Dynamics. Annual Review of Physical Chemistry 2009 , 60 (1), 487 -511. https://doi.org/10.1146/annurev.physchem.59.032607.093818.
- (2) Townsend, D.; Sussman, B. J.; Stolow, A. A Stark Future for Quantum Control. J. Phys. Chem. A 2011 , 115 (4), 357 -373. https://doi.org/10.1021/jp109095d.
- (3) Goulielmakis,  E.;  Loh,  Z.-H.;  Wirth,  A.;  Santra,  R.;  Rohringer,  N.;  Yakovlev,  V.  S.; Zherebtsov, S.; Pfeifer, T.; Azzeer, A. M.; Kling, M. F.; Leone, S. R.; Krausz, F. RealTime  Observation  of  Valence  Electron  Motion. Nature 2010 , 466 (7307),  739 -743. https://doi.org/10.1038/nature09212.
- (4) Geneaux, R.; Marroux, H. J. B.; Guggenmos, A.; Neumark, D. M.; Leone, S. R. Transient Absorption Spectroscopy Using High Harmonic Generation: A Review of Ultrafast XRay Dynamics in Molecules and Solids. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2019 , 377 (2145), 20170463. https://doi.org/10.1098/rsta.2017.0463.
- (5) Kobayashi,  Y.;  Leone,  S.  R.  Characterizing  Coherences  in  Chemical  Dynamics  with Attosecond Time-Resolved x-Ray Absorption Spectroscopy. The Journal of Chemical Physics 2022 , 157 (18), 180901. https://doi.org/10.1063/5.0119942.
- (6) Yan,  Y.  X.;  Gamble,  E.  B.;  Nelson,  K.  A.  Impulsive  Stimulated  Scattering:  General Importance  in  Femtosecond  Laser  Pulse  Interactions  with  Matter,  and  Spectroscopic Applications. J. Chem. Phys. 1985 , 83 (11), 5391 -5399. https://doi.org/10.1063/1.449708.
- (7) Ruhman, S.; Joly, A. G.; Nelson, K. A. Coherent Molecular Vibrational Motion Observed in the Time Domain Through Impulsive Stimulated Raman Scattering. IEEE J. Quantum Electron. 1988 , 24 (2), 460 -469. https://doi.org/10.1109/3.146.
- (8) Dhar, L.; Rogers, J. A.; Nelson, K. A. Time-Resolved Vibrational Spectroscopy in the Impulsive Limit. Chem. Rev. 1994 , 94 (1), 157 -193. https://doi.org/10.1021/cr00025a006.
- (9) Hosler, E. R.; Leone, S. R. Characterization of Vibrational Wave Packets by Core-Level High-Harmonic  Transient  Absorption  Spectroscopy. Physical  Review  A 2013 , 88 (2), 023420. https://doi.org/10.1103/physreva.88.023420.
- (10)  Pertot,  Y.;  Schmidt,  C.;  Matthews,  M.;  Chauvet,  A.;  Huppert,  M.;  Svoboda,  V.;  Von Conta, A.; Tehlar, A.; Baykusheva, D.; Wolf, J.-P.; Wörner, H. J. Time-Resolved x-Ray Absorption Spectroscopy with a Water Window High-Harmonic Source. Science 2017 , 355 (6322), 264 -267. https://doi.org/10.1126/science.aah6114.
- (11)  Warrick, E. R.; Bæ khøj, J. E.; Cao, W.; Fidler, A. P.; Jensen, F.; Madsen, L. B.; Leone, S.  R.;  Neumark,  D.  M.  Attosecond  Transient  Absorption  Spectroscopy  of  Molecular Nitrogen: Vibrational Coherences in the B′ 1Σ+u State. Chemical Physics Letters 2017 , 683 , 408 -415. https://doi.org/10.1016/j.cplett.2017.02.013.
- (12)  Wei, Z.; Li, J.; Wang, L.; See, S. T.; Jhon, M. H.; Zhang, Y.; Shi, F.; Yang, M.; Loh, Z.H.  Elucidating  the  Origins of Multimode  Vibrational  Coherences  of  Polyatomic Molecules Induced by Intense Laser Fields. Nat. Commun. 2017 , 8 (1), 735. https://doi.org/10.1038/s41467-017-00848-2.
- (13)  Wei,  Z.;  Tian,  L.;  Li,  J.;  Lu,  Y.;  Yang,  M.;  Loh,  Z.  H.  Tracking  Ultrafast  Bond Dissociation  Dynamics  at  0.1  Å   Resolution  by  Femtosecond  Extreme  Ultraviolet Absorption Spectroscopy. J. Phys. Chem. Lett. 2018 , 9 (19), 5742 -5747. https://doi.org/10.1021/acs.jpclett.8b02547.
- (14)  Wei, Z.; Li, J.; Zhang, H.; Lu, Y.; Yang, M.; Loh, Z.-H. Ultrafast Dissociative Ionization and Large-Amplitude Vibrational Wave Packet Dynamics of Strong-Field-Ionized DiIodomethane. J. Chem. Phys. 2019 , 151 (21), 214308. https://doi.org/10.1063/1.5132967.
- (15) Saito, N.; Sannohe, H.; Ishii, N.; Kanai, T.; Kosugi, N.; Wu, Y.; Chew, A.; Han, S.; Chang, Z.; Itatani, J. Real-Time Observation of Electronic, Vibrational, and Rotational Dynamics in Nitric Oxide with Attosecond Soft x-Ray Pulses at 400 eV. Optica 2019 , 6 (12), 1542. https://doi.org/10.1364/OPTICA.6.001542.
- (16)  Timmers, H.; Zhu, X.; Li, Z.; Kobayashi, Y.; Sabbar, M.; Hollstein, M.; Reduzzi, M.; Martí nez, T. J.; Neumark, D. M.; Leone, S. R. Disentangling Conical Intersection and Coherent Molecular Dynamics in Methyl Bromide with Attosecond Transient Absorption Spectroscopy. Nat  Commun 2019 , 10 (1),  3133.  https://doi.org/10.1038/s41467-01910789-7.
- (17)  Kobayashi,  Y.;  Neumark,  D.  M.;  Leone,  S.  R.  Attosecond  XUV  Probing  of  Vibronic Quantum Superpositions in Br2+. Phys. Rev. A 2020 , 102 (5), 051102. https://doi.org/10.1103/PhysRevA.102.051102.
- (18)  Kobayashi, Y.; Chang, K. F.; Poullain, S. M.; Scutelnic, V.; Zeng, T.; Neumark, D. M.; Leone, S. R. Coherent Electronic-Vibrational Dynamics in Deuterium Bromide Probed via Attosecond Transient-Absorption Spectroscopy. Phys. Rev. A 2020 , 101 (6), 063414. https://doi.org/10.1103/PhysRevA.101.063414.
- (19) Zinchenko, K. S.; Ardana-Lamas, F.; Seidu, I.; Neville, S. P.; Van Der Veen, J.; Lanfaloni, V.  U.;  Schuurman,  M.  S.;  Wörner,  H.  J.  Sub-7-Femtosecond  Conical-Intersection Dynamics  Probed  at the Carbon  K-Edge. Science 2021 , 371 (6528), 489 -494. https://doi.org/10.1126/science.abf1656.
- (20)  Rott,  F.;  Reduzzi,  M.;  Schnappinger,  T.;  Kobayashi,  Y.;  Chang,  K.  F.;  Timmers,  H.; Neumark, D. M.; de Vivie-Riedle, R.; Leone, S. R. Ultrafast Strong-Field Dissociation of Vinyl Bromide: An Attosecond Transient Absorption Spectroscopy and Non-Adiabatic Molecular Dynamics Study. Struct. Dyn. 2021 , 8 (3). https://doi.org/10.1063/4.0000102.
- (21)  Loh, Z.-H. Studies of Ultrafast Molecular Dynamics by Femtosecond Extreme Ultraviolet Absorption Spectroscopy. Chem. Lett. 2021 , 50 (5), 965 -973. https://doi.org/10.1246/cl.200940.
- (22)  Poullain,  S.  M.;  Kobayashi,  Y.;  Chang,  K.  F.;  Leone,  S.  R.  Visualizing  Coherent Vibrational  Motion  in  the  Molecular  Iodine  B  3  Π0+u  State  Using  Ultrafast  XUV Transient-Absorption Spectroscopy. Phys. Rev. A 2021 , 104 (2), 1 -12. https://doi.org/10.1103/PhysRevA.104.022817.
- (23)  Chang, K. F.; Wang, H.; Poullain, S. M.; Prendergast, D.; Neumark, D. M.; Leone, S. R. Mapping Wave Packet Bifurcation at a Conical Intersection in CH3I by Attosecond XUV Transient Absorption Spectroscopy. J. Chem. Phys. 2021 , 154 (23), 1ENG. https://doi.org/10.1063/5.0056299.
- (24)  Chang, K. F.; Wang, H.; Poullain, S. M.; González-Vázquez, J.; Bañares, L.; Prendergast, D.;  Neumark,  D.  M.;  Leone,  S.  R.  Conical  Intersection  and  Coherent  Vibrational Dynamics in Alkyl Iodides Captured by Attosecond Transient Absorption Spectroscopy. J. Chem. Phys. 2022 , 156 (11), 114304. https://doi.org/10.1063/5.0086775.
- (25)  Matselyukh, D. T.; Despré, V.; Golubev, N. V.; Kuleff, A. I.; Wörner, H. J. Decoherence and Revival in Attosecond Charge Migration Driven by Non-Adiabatic Dynamics. Nat. Phys. 2022 , 18 (10), 1206 -1213. https://doi.org/10.1038/s41567-022-01690-0.
- (26)  Rupprecht, P.; Aufleger, L.; Heinze, S.; Magunia, A.; Ding, T.; Rebholz, M.; Amberg, S.; Mollov, N.; Henrich, F.; Haverkort, M. W.; Ott, C.; Pfeifer, T. Resolving Vibrations in a Polyatomic  Molecule  with  Femtometer  Precision  via  X-Ray  Spectroscopy. Physical Review A 2023 , 108 (3), 032816. https://doi.org/10.1103/physreva.108.032816.
- (27)  Barreau, L.; Ross, A. D.; Kimberg, V.; Krasnov, P.; Blinov, S.; Neumark, D. M.; Leone, S.  R.  Core-Excited  States  of  SF6  Probed  with  Soft-x-Ray  Femtosecond  Transient Absorption  of  Vibrational  Wave  Packets. Physical  Review  A 2023 , 108 (1),  012805. https://doi.org/10.1103/physreva.108.012805.
- (28)  Ross, A. D.; Hait, D.; Scutelnic, V.; Neumark, D. M.; Head-Gordon, M.; Leone, S. R. Measurement  of  Coherent  Vibrational  Dynamics  with  X-Ray  Transient  Absorption Spectroscopy Simultaneously at the Carbon K- and Chlorine L$\_{2,3}$- Edges. arXiv May 8, 2024. http://arxiv.org/abs/2405.05437 (accessed 2024-06-28).
- (29)  Frank, I.; Hutter, J.; Marx, D.; Parrinello, M. Molecular Dynamics in Low-Spin Excited States. The Journal of Chemical Physics 1998 , 108 (10), 4060 -4069. https://doi.org/10.1063/1.475804.
- (30)  Kowalczyk,  T.;  Tsuchimochi,  T.;  Chen,  P.-T.;  Top,  L.;  Van  Voorhis,  T.  Excitation Energies  and  Stokes  Shifts  from  a  Restricted  Open-Shell  Kohn-Sham  Approach. The Journal of Chemical Physics 2013 , 138 (16), 164101. https://doi.org/10.1063/1.4801790.
- (31)  Hait, D.; Head-Gordon, M. Highly Accurate Prediction of Core Spectra of Molecules at Density  Functional  Theory  Cost:  Attaining  Sub-Electronvolt  Error  from  a  Restricted Open-Shell  Kohn -Sham  Approach. J. Phys.  Chem.  Lett. 2020 , 11 (3), 775 -786. https://doi.org/10.1021/acs.jpclett.9b03661.
- (32)  Epifanovsky, E.; Gilbert, A. T. B.; Feng, X.; Lee, J.; Mao, Y.; Mardirossian, N.; Pokhilko, P.; White, A. F.; Coons, M. P.; Dempwolff, A. L.; Gan, Z.; Hait, D.; Horn, P. R.; Jacobson, L.  D.;  Kaliman,  I.;  Kussmann,  J.;  Lange,  A.  W.;  Lao,  K.  U.;  Levine,  D.  S.;  Liu,  J.; McKenzie, S. C.; Morrison, A. F.; Nanda, K. D.; Plasser, F.; Rehn, D. R.; Vidal, M. L.; You, Z.-Q.; Zhu, Y.; Alam, B.; Albrecht, B. J.; Aldossary, A.; Alguire, E.; Andersen, J. H.; Athavale, V.; Barton, D.; Begam, K.; Behn, A.; Bellonzi, N.; Bernard, Y. A.; Berquist, E. J.; Burton, H. G. A.; Carreras, A.; Carter-Fenk, K.; Chakraborty, R.; Chien, A. D.; Closser,  K.  D.;  Cofer-Shabica,  V.;  Dasgupta,  S.;  De  Wergifosse,  M.;  Deng,  J.; Diedenhofen,  M.;  Do,  H.;  Ehlert,  S.;  Fang,  P.-T.;  Fatehi,  S.;  Feng,  Q.;  Friedhoff,  T.; Gayvert, J.; Ge, Q.; Gidofalvi, G.; Goldey, M.; Gomes, J.; González-Espinoza, C. E.; Gulania, S.; Gunina, A. O.; Hanson-Heine, M. W. D.; Harbach, P. H. P.; Hauser, A.; Herbst, M. F.; Hernández Vera, M.; Hodecker, M.; Holden, Z. C.; Houck, S.; Huang, X.; Hui, K.; Huynh, B. C.; Ivanov, M.; Jász, Á .; Ji, H.; Jiang, H.; Kaduk, B.; Kähler, S.; Khistyaev, K.; Kim, J.; Kis, G.; Klunzinger, P.; Koczor-Benda, Z.; Koh, J. H.; Kosenkov, D.; Koulias, L.; Kowalczyk, T.; Krauter, C. M.; Kue, K.; Kunitsa, A.; Kus, T.; Ladjánszki, I.; Landau, A.; Lawler, K. V.; Lefrancois, D.; Lehtola, S.; Li, R. R.; Li, Y.-P.; Liang, J.; Liebenthal, M.; Lin, H.-H.; Lin, Y.-S.; Liu, F.; Liu, K.-Y.; Loipersberger, M.; Luenser,
- A.; Manjanath, A.; Manohar, P.; Mansoor, E.; Manzer, S. F.; Mao, S.-P.; Marenich, A. V.; Markovich, T.; Mason, S.; Maurer, S. A.; McLaughlin, P. F.; Menger, M. F. S. J.; Mewes, J.-M.; Mewes, S. A.; Morgante, P.; Mullinax, J. W.; Oosterbaan, K. J.; Paran, G.; Paul, A. C.; Paul, S. K.; Pavošević, F.; Pei, Z.; Prager, S.; Proynov, E. I.; Rák, Á .; RamosCordoba, E.; Rana, B.; Rask, A. E.; Rettig, A.; Richard, R. M.; Rob, F.; Rossomme, E.; Scheele, T.; Scheurer, M.; Schneider, M.; Sergueev, N.; Sharada, S. M.; Skomorowski, W.; Small, D. W.; Stein, C. J.; Su, Y.-C.; Sundstrom, E. J.; Tao, Z.; Thirman, J.; Tornai, G. J.; Tsuchimochi, T.; Tubman, N. M.; Veccham, S. P.; Vydrov, O.; Wenzel, J.; Witte, J.; Yamada, A.; Yao, K.; Yeganeh, S.; Yost, S. R.; Zech, A.; Zhang, I. Y.; Zhang, X.; Zhang, Y.; Zuev, D.; Aspuru-Guzik, A.; Bell, A. T.; Besley, N. A.; Bravaya, K. B.; Brooks, B. R.; Casanova, D.; Chai, J.-D.; Coriani, S.; Cramer, C. J.; Cserey, G.; DePrince, A. E.; DiStasio, R. A.; Dreuw, A.; Dunietz, B. D.; Furlani, T. R.; Goddard, W. A.; HammesSchiffer, S.; Head-Gordon, T.; Hehre, W. J.; Hsu, C.-P.; Jagau, T.-C.; Jung, Y.; Klamt, A.; Kong, J.; Lambrecht, D. S.; Liang, W.; Mayhall, N. J.; McCurdy, C. W.; Neaton, J. B.; Ochsenfeld, C.; Parkhill, J. A.; Peverati, R.; Rassolov, V. A.; Shao, Y.; Slipchenko, L. V.; Stauch, T.; Steele, R. P.; Subotnik, J. E.; Thom, A. J. W.; Tkatchenko, A.; Truhlar, D.  G.;  Van  Voorhis,  T.;  Wesolowski,  T.  A.;  Whaley,  K.  B.;  Woodcock,  H.  L.; Zimmerman, P. M.; Faraji, S.; Gill, P. M. W.; Head-Gordon, M.; Herbert, J. M.; Krylov, A. I. Software for the Frontiers of Quantum Chemistry: An Overview of Developments in  the  Q-Chem  5  Package. The  Journal  of  Chemical  Physics 2021 , 155 (8),  084801. https://doi.org/10.1063/5.0055522.
- (33)  Hui, K.; Chai, J.-D. SCAN-Based Hybrid and Double-Hybrid Density Functionals from Models  without  Fitted  Parameters. The  Journal  of  Chemical  Physics 2016 , 144 (4), 044114. https://doi.org/10.1063/1.4940734.
- (34)  Dunning, T. H. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron through Neon and Hydrogen. The Journal of Chemical Physics 1989 , 90 (2), 1007 -1023. https://doi.org/10.1063/1.456153.
- (35)  Kendall, R. A.; Dunning, T. H.; Harrison, R. J. Electron Affinities of the First-Row Atoms Revisited. Systematic Basis Sets and Wave Functions. The Journal of Chemical Physics 1992 , 96 (9), 6796 -6806. https://doi.org/10.1063/1.462569.
- (36)  Wilson, A. K.; Woon, D. E.; Peterson, K. A.; Dunning, T. H. Gaussian Basis Sets for Use in  Correlated  Molecular  Calculations.  IX.  The  Atoms  Gallium  through  Krypton. The Journal of Chemical Physics 1999 , 110 (16), 7667 -7676. https://doi.org/10.1063/1.478678.
- (37)  Dyall, K.  G.  Interfacing  Relativistic  and Nonrelativistic  Methods.  I.  Normalized Elimination of the Small Component in the Modified Dirac Equation. The Journal of Chemical Physics 1997 , 106 (23), 9618 -9626. https://doi.org/10.1063/1.473860.
- (38)  Kutzelnigg, W.; Liu, W. Quasirelativistic Theory Equivalent to Fully Relativistic Theory. The Journal of Chemical Physics 2005 , 123 (24), 241102. https://doi.org/10.1063/1.2137315.
- (39)  Saue, T. Relativistic Hamiltonians for Chemistry: A Primer. ChemPhysChem 2011 , 12 (17), 3077 -3094. https://doi.org/10.1002/cphc.201100682.
- (40)  Cunha,  L.  A.;  Hait,  D.;  Kang,  R.;  Mao,  Y.;  Head-Gordon,  M.  Relativistic  OrbitalOptimized Density Functional Theory for Accurate Core-Level Spectroscopy. J. Phys. Chem. Lett. 2022 , 13 (15), 3438 -3449. https://doi.org/10.1021/acs.jpclett.2c00578.
- (41)  Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965 , 140 (4A), A1133 A1138. -https://doi.org/10.1103/PhysRev.140.A1133.
- (42)  Thomassen, H.; Hedberg, K. Distances, Root-Mean-Square Amplitudes of Vibration, and Vibrational Anharmonicity  in Carbon  Tetrabromide  from  a  Gas-Phase  Electron-

Diffraction Investigation. Journal of Molecular Structure 1990 , 240 , 151 -157. https://doi.org/10.1016/0022-2860(90)80506-F.

- (43)  Takehiko Shimanouchi. Tables of Molecular Vibrational Frequencies Part 1 ;  National Standard Reference Data Series; National Bureau of Standards, 1967.
- (44)  Hait, D.; Head-Gordon, M. Orbital Optimized Density Functional Theory for Electronic Excited States. J. Phys. Chem. Lett. 2021 , 12 (19), 4517 -4529. https://doi.org/10.1021/acs.jpclett.1c00744.
- (45)  Hait, D.; Head-Gordon, M. Excited State Orbital Optimization via Minimizing the Square of the Gradient: General Approach and Application to Singly and Doubly Excited States via  Density  Functional  Theory. J.  Chem.  Theory  Comput. 2020 , 16 (3),  1699 -1710. https://doi.org/10.1021/acs.jctc.9b01127.
- (46)  Boo, B. H.; Saito, N.; Suzuki, I. H.; Koyano, I. Dissociation Processes of Core-Excited CBr4  Involving  the  Br(3d,  3p,  3s)  and  C(1s)  Inner-Shells  in  the  Range  50 -460  eV. Journal  of  Electron  Spectroscopy  and  Related  Phenomena 2002 , 123 (1),  73 -84. https://doi.org/10.1016/s0368-2048(02)00005-1.
- (47)  Xu,  Z.;  Hu,  P.;  Wang,  E.;  Xu,  S.;  Wang,  X.;  Zhao,  Y.;  Deng,  J.;  Ning,  C.;  Ren,  X. Experimental and Theoretical Study of Valence Electronic Structure of Tetrabromomethane by ( e , 2 e ) Electron Momentum Spectroscopy. Phys. Rev. A 2019 , 99 (6), 062705. https://doi.org/10.1103/PhysRevA.99.062705.

## Supporting Information for 'Attosecond Probing of

## Coherent Vibrational Dynamics in CBr4'

Jen-Hao Ou,  Diptarka Hait, 1 2,3 Patrick Rupprecht, 1,4 John E. Beetar, 1,† Todd J. Martí nez, 2,3 and Stephen R. Leone 1,4,5,*

1 Department of Chemistry, University of California, Berkeley, CA 94720, USA

2 Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States

3 SLAC National Accelerator Laboratory, Menlo Park, California 94024, United States

4 Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720,

## USA

5 Department of Physics, University of California, Berkeley, CA 94720, USA

† Present address: NSF National eXtreme Ultrafast Science (NeXUS) Facility, Columbus, OH 43210, USA

E-mail: srl@berkeley.edu

## 1. Calibration of Pump Pulse

In this work, the peak intensity of pump pulse is modified by cropping the beam with an iris before the pump and probe pulses recombine. Detailed measurements of pump-pulse properties are given below.

Figure  S1. Beam  diameter  (mean μ )  of  pump  pulse  vs  iris  size. The  iris  size  is  read  by  the numbered marks on the iris. Measured with DataRay beam profiler WinCamD-UHR, 4xsigma mode. Picked off to the location at the same distance to the sample location. The beam diameter (mean μ ) remains around 109 μ m and does not change significantly when the iris size decreases from full to 3.0. Only when iris size is smaller than 3.0, the beam starts to get significantly cropped and the beam diameter increases appreciably. The averaged value of beam diameters from size 3.0 to full is 109 μ m with 1 μ m standard deviation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-62266681-8bff-4d6d-aa8a-1772d79c5c36.jpg)

Figure S2. Beam shape with iris size full measured by DataRay beam profiler WinCamD-UHR.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-a22be171-056a-4db3-b0ca-9344de5431f9.jpg)

The mean is 107.7 μ m.

Figure S3. Beam shape with iris size 3.5 measured by DataRay beam profiler WinCamD-UHR.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-051094bb-8cc6-4684-b1ee-a93b6a8859b7.jpg)

The mean is 108.1 μ m.

Figure S4. Beam shape with iris size 1.0 measured by DataRay beam profiler WinCamD-UHR.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-10adad80-d128-422a-bc9b-870af05a7c7c.jpg)

The mean is 156.2 μ m.

Figure S5. Average power of pump pulse vs iris size. The iris size is read by the numbered marks on the iris. The average power at the sample location (red circle) is estimated as 65% of that measured right after the iris (black square), according to a previous power calibration. There are 5 mirrors and 1 thin glass window between the iris and the sample cell.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-dc964397-c69e-463e-8593-e7e6ef6081d0.jpg)

Table S1. Peak intensity of pump pulse vs iris size.

|   Scan # | Iris size a   |   Average Power (mW) b | Beam Diameter, Mean μ (μm) c   | Pulse Duration (fs) d   |   Pulse Energy ( m J) e |   Peak Intensity (×10 12 W/cm 2 ) f |   Peak Intensity (×10 -4 atomic unit) g |
|----------|---------------|------------------------|--------------------------------|-------------------------|-------------------------|-------------------------------------|-----------------------------------------|
|        2 | 3.5           |                   11.9 | 109±1                          | 25.8±0.6                |                    11.9 |                                8.15 |                                    2.32 |
|        3 | 4             |                   15   | 109±1                          | 25.8±0.6                |                    15   |                               10.3  |                                    2.93 |
|        5 | 4.5           |                   18.3 | 109±1                          | 25.8±0.6                |                    18.3 |                               12.5  |                                    3.57 |
|        4 | 6.5           |                   28   | 109±1                          | 25.8±0.6                |                    28   |                               19.2  |                                    5.47 |
|        1 | full          |                   29.5 | 109±1                          | 25.8±0.6                |                    29.5 |                               20.2  |                                    5.76 |

a Iris size is read by the numbered marks on the iris.

b At the sample cell location. Estimated according to previous power calibration (Figure S5)

c Measured with DataRay beam profiler WinCamD-UHR, 4xsigma mode. Details in Figure S1.

d Measured by SD-FROG (Figure S6(d))

e Average power divided by 1 kHz repetition rate

f Assume Gaussian spatial profile and FROG-retrieved temporal profile (Figure S6)

g Peak intensity in W/cm  divided by the conversion factor 3.50944758 2 × 10 16 from the literature . 1

Figure S6. Self-Diffraction Frequency Resolved Optical Gating (SD-FROG) measurement for 400 nm pump pulse (a) Measured and (b) reconstructed spectrogram. (c) Measured (solid line) and reconstructed (dash line) spectrum. The central wavelength is around 400 nm. (d) Reconstructed temporal intensity envelope. The pulse duration is considered as 25.8±0.6 fs FWHM from fitting the main pulse to a Gaussian function.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-6af6b7d2-ca25-486b-84b5-c250a4fcecdb.jpg)

## 2. Calibration of Time Zero and Temporal Resolution for Time-Resolved XUV Absorption Measurements

The time zero and temporal resolution (instrument response function) is calibrated with a crosscorrelation measurement of pump and probe pulses in He 2s2p state.

Figure S7. The cross-correlation measurement of pump and probe pulses in He 2s2p state. (a)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-9fa7399d-5be6-4817-b719-c33015600fa0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-8074ae89-9ad1-4938-94c0-c316fff73283.jpg)

Time-resolved XUV differential absorption ( Δ A) spectrum of He 2s2p state. The XUV absorption is suppressed ( Δ A turns negative) by the pump pulse when pump and probe pulses overlap. (b) Circle: Negative of measured differential absorption (Δ A) at 60.16 eV. Black solid line: Fit by a convolution of two Gaussian (main and pre-pulse) and an exponential decay function in time delay. Red dash line: Deconvoluted Gaussian (main pulse). The fitted FWHM and center of this Gaussian is 26 ±6 fs and -1±3 fs, and these are taken as the instrument response function and time zero for this work, respectively. Red dash line: Deconvoluted Gaussian (pre-pulse)

## 3. Calibration of Spectral Axis

In this work, the XUV beam is dispersed by a grating onto a XUV CCD camera. As a result, the spectral axis of the raw data is recorded in camera pixels. We measure the static XUV absorption of He 2snp (n = 2 to 5) states, and establish the calibration curve for converting pixel to wavelength and photon energy.

Figure S8. Static XUV absorption spectrum of He 2snp (n ≥ 2) excited states, as a function of

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-b05288b0-293e-461f-ad50-4ddd50beff8e.jpg)

XUV CCD camera pixel.

Table S2. Calibration Table with He 2snp (n = 2 to 5) States

| Electronic State   | Resonance Location (Pixel) a   |   Resonance Energy (eV), Literature 2 |   Resonance Wavelength (nm) b |
|--------------------|--------------------------------|---------------------------------------|-------------------------------|
| He 2s5p            | 608 ±1                         |                                64.814 |                        19.129 |
| He 2s4p            | 621.2 ±0.4                     |                                64.467 |                        19.232 |
| He 2s3p            | 652.2 ±0.2                     |                                63.658 |                        19.476 |
| He 2s2p            | 793. 4±0.2                     |                                60.147 |                        20.613 |

a Fitted to a convolution of Fano lineshape and Gaussian function, to determine the resonance location in pixel. Estimated manually for 2s5p.

b Converted from resonance energy (eV) in the literature. 2

Figure S9. The calibration curve for wavelength vs pixel. Cross: Measured resonance location.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-80a57267-ad86-4a9c-9afd-b83446c1b939.jpg)

Red line: Linear fit. This calibration curve enables the conversion from pixel to wavelength and photon energy.

Figure S10. Calibrated static XUV absorption spectrum of He 2snp (n ≥ 2) excited states, as a function of photon energy.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-f7894a2e-30dc-4b53-8858-fe2f4dc4b169.jpg)

Figure  S11. Time-resolved  XUV  differential  absorption  spectrum  of  CBr4  with  pump  peak intensity  20.2 × 10 12 W/cm 2 .  The  negative  oscillating  signals  around  69  to  71  eV  encode  the coherent  vibrational  dynamics  on  the  electronic  ground  state.  (Upper)  Differential  absorption integrated from time delay of 300 to 1100 fs. Three major peaks between 64 to 66 eV correspond to dissociated atomic Br. These are likely to come from multiphoton absorption of the 400 nm pump pulse (photon energy 3.1 eV), because the lowest transition to valence-excited states (A bands) of CBr4 starts at 4.9 eV (39550 cm -1 ) 3,4

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-363348b9-c6e2-40ef-8f4b-83941e62666f.jpg)

Table S3. Energy of dissociated atomic Br 3d absorption. Our values agree reasonably with the literature. This also confirms the validity of our spectral axis calibration.

| Transition           |   This work Energy (eV) |   Previous Literature 5 Energy (eV) |
|----------------------|-------------------------|-------------------------------------|
| Br 2 P 3/2 → 2 D 3/2 |                   65.55 |                               65.58 |
| Br 2 P 1/2 → 2 D 3/2 |                   65.11 |                               65.13 |
| Br 2 P 3/2 → 2 D 5/2 |                   64.56 |                               64.54 |

## 4. Measured Relative Slope of Br-3d Core-Excited State Potential Energy Surface (PES)

Table  S4. Measured  relative  slope  of  (Br-3d3/2,5/2) -1 10a1 * core-excited  state  PES  along  the symmetric stretch mode. Relative slope is determined experimentally as the ratio of measured oscillation amplitude of XUV absorption energy XUV E  (Br-3d5/2)/ XUV E  (Br-3d3/2).

| Scan #      | Peak Intensity (×10 12 W/cm 2 )   | XUV E  (Br- 3d 3/2 ) (meV)   | XUV E  (Br- 3d 5/2 ) (meV)   | Relative Slope   |
|-------------|-----------------------------------|-------------------------------|-------------------------------|------------------|
| 2           | 8.15                              | 2.8±0.6                       | 4.1±0.9                       | 1.5±0.4          |
| 3           | 10.3                              | 9.6±1.2                       | 9.0±1.3                       | 0.9±0.2          |
| 5           | 12.5                              | 7.2±1.7                       | 7.6±1.6                       | 1.0±0.3          |
| 4           | 19.2                              | 11.9±1.9                      | 12.2±1.7                      | 1.0±0.2          |
| 1           | 20.2                              | 10.2±1.4                      | 11.3±1.4                      | 1.1± 0.2         |
| Mean ±Error |                                   |                               |                               | 1.1 ±0.13        |

## 5. Comment about Initial Phase of Wavepacket Trajectory

If the launching mechanism of the vibrational wavepacket is a pure non-resonant ISRS from a Gaussian pump pulse centered at t = 0, the expectation value of the wavepacket displacement will be a sine function with a zero initial phase. 6,7 Nuclear TDSE simulations reproduce this expectation by  retrieving  around 0.01  π initial  phase  of  wavepacket  position  for  all  intensities,  with  the experimental pump pulse intensity profile measured by SD-FROG (as shown in Figure S6(d)). Experimentally, over the range of intensities measured, the initial phase of XUV absorption energy varies from -0.1 to -0.6 π (Figure S12), suggesting a possible deviation from the pure non-resonant ISRS. Potential cause may be the pre-pulse (Figure S6) or the multiphoton absorption observed in the time-resolved XUV differential absorption spectrum (Figure S11).

Figure S12. The initial phase at time zero from fitting the measured XUV absorption energy as a function of time-delay to a sine function for different pump-pulse intensities.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00696v1-figure-6d9c338e-51fa-437f-b182-91712a96cf4a.jpg)

## 6. Nuclear Time-Dependent Schrodinger Equation Simulations

We simulated  the  dynamics  of  the  symmetric  stretch  mode  in  a  manner  similar  to  what  was described in a previous work.  In brief, under the harmonic approximation (which is valid in the 8 small-displacement limit), the normal modes are all independent of each other and can be treated separately. The symmetric stretching mode evolves under the influence of the Hamiltonian:

$$\widehat { H } ( t ) = \frac { \hat { p } ^ { 2 } } { 2 \mu } + \frac { 1 } { 2 } \mu \omega _ { 0 } ^ { 2 } \, \hat { q } ^ { 2 } - \frac { 1 } { 2 } \alpha ( \hat { q } ) | \mathcal { E } ( t ) | ^ { 2 }$$

where 𝑝̂ is the momentum operator along this mode, 𝑞̂ is the operator for bond stretching (relative to  the  equilibrium  bond length), 𝛼(𝑞 ̂) is  the  isotropic  polarizability  as  a  function  of 𝑞 ̂ , 𝜇 is  the normal mode mass, 𝜔0 is the normal mode frequency, and |ℰ(𝑡)| 2 = 𝐼(𝑡) is the temporal intensity envelope of the pump pulse. Details of the above terms are as follows.

The isotropic polarizability 𝛼(𝑞 ̂) tends  to  increase  with  bond  stretching  as  the  bonding electron density becomes 'softer' (i.e. more diffuse) in response . 9 For small displacements, it can be approximated with 𝛼(𝑞) = 𝛼0  + 𝛼1𝑞  + 𝛼2𝑞 2 .  Fitting  SCAN0/decontracted aug-cc-pVTZ polarizabilities across 201 uniformly spaced bond stretches spanning ± 0.1 Å relative to equilibrium for CBr4 yields 𝛼0   =  97.76, 𝛼 1 =  38.70, 𝛼 2 =  16.90 , in atomic units.

The normal mode mass 𝜇 = 575437.04 atomic units (i. e. 4 times the mass of a single 79 Br atom, assuming all four Br atoms in the molecule correspond to this isotope). The normal mode frequency is 0.00129 atomic units (282.32 cm -1 ) from the ground state DFT calculations.

The temporal intensity envelope of the pump pulse |ℰ(𝑡)| 2 = 𝐼(𝑡) is calculated as

$$I ( t ) = I _ { \mathrm p k } ( \mathrm a. u. ) I _ { \text{FROG} } ( t ) \cos ^ { 2 } ( \omega t )$$

where  pk (a.u.) is the calibrated peak intensity in atomic units (Table S1), I I FROG( ) is the normalized t temporal  intensity  envelope  from  SD-FROG  measurement  (Figure  S6(d)), ω is  the  angular frequency of 400 nm pump pulse and equals 0.1139 atomic units.

With  these  parameters,  the time-dependent  Schrödinger  equation  (TDSE)  for  the symmetric stretching mode was numerically simulated within the interaction picture for each of the experimentally studied pump laser intensities, taking the field-free harmonic oscillator terms as the reference Hamiltonian and the - 1 2 𝛼(𝑞 ̂)|ℰ(𝑡)| 2 term as the perturbation. The simulations were done in the basis of the ten lowest energy eigenstates of the unperturbed harmonic oscillator Hamiltonian. The wavefunction at t = -4134.14 ato mic units (100 fs) was assumed to be the ground state of the unperturbed harmonic oscillator, and the wave-function was subsequently propagated in time to t = 41339.86 atomic units (1000 fs) with step sizes of 2 atomic units (0.048 fs).

## 7. Comment about Polarizability Calculations

The TDSE simulations along the symmetric stretch mode appear to predict bond stretches that are roughly  three  times  as  great  for  a  given  pump  power  intensity  as  those  estimated  from  the experimental Δ𝐸𝑋𝑈𝑉 values and the computed Br core-excited PES slope of -9.4 eV/ Å . A potential source of this discrepancy might be the use of too large polarizabilities for the TDSE simulations, as more simplified previous models  have indicated that the bond stretch amplitude is directly 7 proportional to 𝛼1 (i.e. the slope of the polarizability vs bond stretch at the equilibrium geometry).

Such an overestimation could potentially arise from either errors in the finite difference approximation used for computing polarizabilities or errors intrinsic to the SCAN0/decontracted aug-cc-pVTZ model chemistry. Both scenarios appear to be unlikely, as discussed below.

7.1 Finite difference: Q-Chem presently lacks analytic polarizabilities with the X2C relativistic method. As a result, the polarizability  was  computed  via finite differences of the energy with different electric fields. As the polarizability is the negative of the second derivative of the energy vs applied field strength, we can utilize standard finite difference formulae and obtain:

$$\alpha _ { z z } = - \left ( \frac { d ^ { 2 } } { d \mathcal { E } _ { z } ^ { 2 } } E ( \mathcal { E } _ { z } ) \right ) _ { \mathcal { E } _ { z } = 0 } = \frac { 3 0 E ( 0 ) - 1 6 E ( h ) - 1 6 E ( - h ) + E ( 2 h ) + E ( - 2 h ) } { 1 2 h ^ { 2 } } + O ( h ^ { 4 } )$$

where 𝐸(ℰ𝑧 ) is the energy of the system computed with an applied field ℰ𝑧 along the z direction and ℎ = 0.005 atomic  units  is  the  applied  electric  field  strength  for  the  finite  difference calculations. As CBr4 is tetrahedral, the isotropic polarizability 𝛼 = 𝛼𝑧𝑧 . Furthermore,the lack of a permanent dipole moment ensures 𝐸(ℎ) = 𝐸(-ℎ) . So the previous five-point stencil simplifies to the following three point formula:

$$\alpha _ { z z } = - \left ( \frac { d ^ { 2 } } { d \varepsilon _ { z } ^ { 2 } } E ( \varepsilon _ { z } ) \right ) _ { \varepsilon _ { z } = 0 } = \frac { 1 5 E ( 0 ) - 1 6 E ( h ) + E ( 2 h ) } { 6 h ^ { 2 } } + O ( h ^ { 4 } )$$

As a result, for any given bond stretch,the energies are computed for zero external fields, and  with  fields  of  0.005  a.u.  and  0.01  a.u.  along  the  z  direction  ( 𝐸(0), 𝐸(ℎ) and 𝐸(2ℎ) , respectively) to compute the polarizability. The scale of the finite difference error can be estimated by comparing the finite difference results in the nonrelativistic (i.e. X2C free) SCAN0/decontracted aug-cc-pVTZ case, for which analytic values can be computed. The resulting values for the 201 uniformly spaced bond stretches spanning ± 0.1 Å relative to the equilibrium

geometry  that  were  used  to  fit  the  quadratic  polynomial 𝛼(𝑞) = 𝛼0  + 𝛼1𝑞  + 𝛼2𝑞 2 indicate very little finite difference error. Specifically, the largest deviation between analytic and finite difference nonrelativistic SCAN0/decontracted aug-cc-pVTZ isotropic polarizabilities is 0.0008 atomic  units  for  the  geometries  utilized  and  fits  from  the  two  approaches  predict  identical 𝛼1 values of 38.18 atomic units. Finite difference errors in the X2C SCAN0/decontracted aug-ccpVTZ isotropic  polarizabilities  therefore  appear  to  be  not  a  major  concern.  We  also  take  the opportunity to note that relativistic effects appear to slightly increase the isotropic polarizability in this  molecule. All the computed polarizability values are provided in the accompanying Excel spreadsheet.

7.2  Electronic  Structure: The  other  potential  source  for  errors  in  the  polarizability  may  be inadequacies of the SCAN0/decontracted aug-cc-pVTZ model chemistry, even though it has been shown to be accurate for lighter elements 10 . We therefore carried out nonrelativistic CCSD(T)/augcc-pVTZ  calculations  (utilizing  the  aforementioned  finite  difference  formula)  for  a  few  bond stretches  to  determine  if  SCAN0/decontracted aug-cc-pVTZ was acceptable. Using 41 equally spaced geometries between 1.9046-1.946 Å  (same as those used in the main text for finding the Br core-excited PES slope from ROKS), we find that nonrelativistic CCSD(T)/aug-cc-pVTZ yields very  similar  polarizabilities  (albeit  slightly  smaller  by  0.1-0.2  atomic  units)  as  nonrelativistic SCAN0/ decontracted aug-cc-pVTZ. The fitted 𝛼1 values in this interval are: 38.02 atomic units from SCAN0/decontracted aug-cc-pVTZ and 36.64 atomic units from CCSD(T)/aug-cc-pVTZ. However, the CCSD(T) calculations utilize a frozen core (C 1s, Br 1s 2s 2p) and the decontracted aug-cc-pVTZ basis used in the SCAN0 calculations is more flexible (and hence likely to lead to larger  polarizabilities)  than  the  aug-cc-pVTZ  basis  used  for  the  CCSD(T)  calculations.  It

nonetheless appears clear that the polarizabilities obtained from finite differences on SCAN0/decontracted aug-cc-pVTZ with X2C are quite reasonable, and lead to 𝛼1 within 5% of CCSD(T). The polarizability fits are therefore extremely unlikely to contribute to the factor of 3 difference between the bond stretching amplitude predicted by TDSE simulations and experiment.

## Reference

- (1)  Krainov, V. P.; Reiss, H. R.; Smirnov, B. M. Radiative Processes in Atomic Physics , 1st ed.; Wiley, 1997. https://doi.org/10.1002/3527605606.
- (2)  Domke, M.; Schulz, K.; Remmers, G.; Kaindl, G.; Wintgen, D. High-Resolution Study of 1 P o Double-Excitation States in Helium. Phys. Rev. A 1996 , 53 (3), 1424 -1438. https://doi.org/10.1103/PhysRevA.53.1424.
- (3)  Causley, G. C.; Russell, B. R. Vacuum Ultraviolet Absorption Spectra of the Bromomethanes. The Journal of Chemical Physics 1975 , 62 (3), 848 -857. https://doi.org/10.1063/1.430535.
- (4)  Olthoff, J. K.; Moore, J. H.; Tossell, J. A. Electron Attachment by Chloro and Bromomethanes. The Journal of Chemical Physics 1986 , 85 (1), 249 -254. https://doi.org/10.1063/1.451644.
- (5)  Nahon, L.; Morin, P.; Farnoux, F. C. Relaxation of the Nd \rightarrow (N+ 1)p Resonances in Atomic Bromine and Lodine. Physica Scripta 1992 , T41 , 104 -112. https://doi.org/10.1088/0031-8949/1992/t41/017.
- (6)  Yan,  Y.  X.;  Gamble,  E.  B.;  Nelson,  K.  A.  Impulsive  Stimulated  Scattering:  General Importance  in  Femtosecond  Laser  Pulse  Interactions  with  Matter,  and  Spectroscopic Applications. J. Chem. Phys. 1985 , 83 (11), 5391 -5399. https://doi.org/10.1063/1.449708.
- (7)  Dhar,  L.;  Rogers,  J.  A.;  Nelson,  K.  A.  Time-Resolved  Vibrational  Spectroscopy  in  the Impulsive Limit. Chem. Rev. 1994 , 94 (1), 157 -193. https://doi.org/10.1021/cr00025a006.
- (8)  Ross,  A.  D.;  Hait,  D.;  Scutelnic,  V.;  Neumark,  D.  M.;  Head-Gordon,  M.;  Leone,  S.  R. Measurement of Coherent Vibrational Dynamics with X-Ray Transient Absorption Spectroscopy Simultaneously at the Carbon K- and Chlorine L$\_{2,3}$- Edges. arXiv May 8, 2024. http://arxiv.org/abs/2405.05437 (accessed 2024-06-28).

(9) Hait, D.; Head-Gordon, M. When Is a Bond Broken? The Polarizability Perspective. Angew Chem Int Ed 2023 , 62 (46), e202312078. https://doi.org/10.1002/anie.202312078.

- (10) Hait,  D.;  Head-Gordon,  M.  How  Accurate  Are  Static  Polarizability  Predictions  from Density Functional Theory? An Assessment over 132 Species at Equilibrium Geometry. Phys. Chem. Chem. Phys. 2018 , 20 (30), 19800 -19810. https://doi.org/10.1039/C8CP03569E.
Jen-Hao Ou Diptarka Hait Patrick Rupprecht John E. Beetar Todd J. Martínez Stephen R. Leone
Attosecond Probing of Coherent Vibrational Dynamics in CBr$_4$
A coherent vibrational wavepacket is launched and manipulated in the
symmetric stretch (a$_1$) mode of CBr$_4$, by impulsive stimulated Raman
scattering from non-resonant 400 nm laser pump pulses with various peak
intensities on the order of tens of 10$^{12}$ W/cm$^2$. Extreme ultraviolet
(XUV) attosecond transient absorption spectroscopy (ATAS) records the
wavepacket dynamics as temporal oscillations in XUV absorption energy at the
bromine M$_{4,5}$ 3d$_{3/2,5/2}$ edges around 70 eV. The results are augmented
by nuclear time-dependent Schr\"odinger equation simulations. Slopes of the
(Br-3d$_{3/2,5/2}$)$^{-1}$10a$_1^*$ core-excited state potential energy surface
(PES) along the a$_1$ mode are calculated to be -9.4 eV/{\AA} from restricted
open-shell Kohn-Sham calculations. Using analytical relations derived for the
small-displacement limit with the calculated slopes of the core-excited state
PES, a deeper insight into the vibrational dynamics is obtained by retrieving
the experimental excursion amplitude of the vibrational wavepacket and the
amount of population transferred to the vibrational first-excited state, as a
function of pump-pulse peak intensity. Experimentally, the results show that
XUV ATAS is capable of easily resolving oscillations in the XUV absorption
energy on the order of few to tens of meV and tens of femtosecond time
precision, limited only by the averaging times in the experimental scans. This
corresponds to oscillations of C-Br bond length on the order of 10$^{-4}$ to
10$^{-3}$ {\AA}. The results and the analytic relationships offer a clear
physical picture, on multiple levels of understanding, for how the pump-pulse
intensity controls the vibrational dynamics launched by non-resonant ISRS in
the small-displacement limit.