Michael E.
Ritter
,
Simone A.
DeSouza
,
Hannah M.
Ogden
,
Tara J.
Michael
and
Amy S.
Mullin
*
Department of Chemistry and Biochemistry, University of Maryland College Park, College Park, Maryland 20742, USA. E-mail: mullin@umd.edu
First published on 5th April 2024
Collisions of optically centrifuged CO2 molecules with J = 244–282 (Erot = 22800–30300 cm−1) are investigated with high-resolution transient IR absorption spectroscopy to reveal collisional and orientational phenomena of molecules with hyper-thermal rotational energies. The optical centrifuge is a non-resonant optical excitation technique that uses ultrafast, 800 nm chirped pulses to drive molecules to extreme rotational states through sequential Raman transitions. The extent of rotational excitation is controlled by tuning the optical bandwidth of the excitation pulses. Frequencies of 30 R-branch ν3 fundamental IR probe transitions are measured for the J = 186–282 states of CO2, expanding beyond previously reported IR transitions up to J = 128. The optically centrifuged molecules have oriented angular momentum and unidirectional rotation. Polarization-sensitive transient IR absorption of individual rotational states of optically centrifuged molecules and their collision products reveals information about collisional energy transfer, relaxation kinetics, and dynamics of rotation-to-translation energy transfer. The transient IR probe also measures the extent of polarization anisotropy. Rotational energy transfer for lower energy molecules is discussed in terms of statistical models and a comparison highlights the role of increasing energy gap with J and angular momentum of the optically centrifuged molecules.
High energy rotational states of molecules cannot be populated using traditional optical approaches because selection rules for absorption and Raman processes limit changes in the rotational angular momentum to small ΔJ values. Recent advances in optical control of molecular rotation have led to new methods for controlling alignment, orientation, energy and directionality of molecular rotors.10–50 These techniques include using a pair of orthogonally polarized pulses of light with a time delay, a sequence of pulse trains with varying pulse lengths, and a pair of oppositely chirped pulses. The first two techniques can control a component of molecular motion but are limited to relatively low-J states.10–15 The third technique is the optical centrifuge which can prepare molecules in extreme rotational states and has been used to investigate the properties and dynamics for a number of systems.16–50
In the work reported here, we use a tunable optical centrifuge to prepare rotationally excited CO2 molecules and investigate their collision dynamics using high-resolution transient IR absorption spectroscopy. The collision studies focus on CO2 with J = 244–282 and rotational energies of Erot = 22800–30300 cm−1. To perform these studies, we first used high-resolution transient IR absorption spectroscopy to identify line-center frequencies of 30 R-branch transitions with J = 186–282 of the CO2ν3 fundamental absorption band and then used the IR transitions to investigate the collision dynamics of the optically centrifuged CO2 molecules.
Previously, high-resolution IR transitions for CO2 have been reported up to J = 128 in the HITRAN database.51,52 Transitions up to J = 146 have also been observed in dispersed IR emission from a flame, but populating high-J states by thermal heating leads to spectral broadening and congestion, making identification of state-resolved transition frequencies difficult.1 The HITEMP database reports transitions up to J = 200 based on semi-empirical calculations.53 The Carbon Dioxide Spectral Database at 4000 K (CDSD4000) reports IR absorption transitions for CO2 at 4000 K based on effective Hamiltonian calculations that include 41 vibrational states and 300 rotational states.54 For the IR transitions reported here, we used a combination of a third-order expansion of the rigid rotor model and the CDSD4000 as a guide for identifying the CO2 transition frequencies.
Spectral lines of high rotational states of CO2 are used to study combustion processes, and the atmospheres of Venus and exoplanets.5–8,55–57 These studies have relied on calculated transition frequencies to identify the high-J CO2 states. The new high-J transition frequencies reported here will serve as benchmark for future spectral calculations and databases.54,58
Previously, we have investigated collisions of optically centrifuged CO2 by probing P-branch transitions for states with J < 104.18,20,21,23,24 In the first of these studies in 2011, single-pass IR detection was used to investigate the collision dynamics of centrifuged CO2 molecules with 62 ≤ J ≤ 88.18 With the single-pass IR detection, a polarization-sensitive study of the J = 76 state of CO2 showed that even after 1000 collisions, there was a preference for population to remain in the original plane of rotation.20 In 2015, a J-dependent survey of CO2 states 0 ≤ J ≤ 100 showed that collisions of centrifuged molecules preferentially involve a translationally cool subset of molecules from the low-J states.21 This study measured large translational energies for states with J = 36–100 that result from collisions with optically centrifuged molecules. In 2017 and 2018, increased signal intensities from a multi-pass detection setup enabled studies of CO2J = 76–100. The first study highlighted the anisotropic kinetic energy release from collisions and orientational anisotropy values and decay rates that increase as a function of J.23 The second study identified the importance of rotational adiabaticity in collisions of centrifuged CO2 with buffer gases.24
The work presented here builds on the previous studies of optically centrifuged CO2 by extending the range of J-states up to J = 282. For comparison, the rotational energy of the J = 100 state is near 3930 cm−1 with a downward energy gap of ΔE = 150 cm−1, while the energy of the J = 282 state is 30300 cm−1 with a comparable energy gap of ΔE = 415 cm−1. Our dynamics studies use polarization-sensitive transient IR absorption to characterize the inverted rotational distribution of the centrifuged molecules, their polarization anisotropies, the kinetics of the state-resolved collisional relaxation, and the magnitude of rotation-to-translation energy transfer by measuring Doppler-broadened line profiles.
The optical centrifuge pulse uses an amplified ultrafast Ti:Sapphire laser system, which generates pulses of light centered at λ = 806 nm. The pulses are stretched in time from 40 fs to 100 ps and amplified in a regenerative amplifier. The temporal stretching induces a positive chirp. A pulse shaper spectrally splits the original pulses into pairs of pulses, one with λ > 806 nm denoted ω1 and the other with λ < 806 nm denoted ω2. A negative chirp is induced into ω2 such that the pair of pulses have an initial wavelength of λ0 = 806 nm and their optical frequency difference increases over the time of the pulse. The time-dependent frequency difference is the source of the sequential Raman transitions that lead to optically centrifuged molecules.
The oppositely chirped pulses are sent through a multi-pass amplifier, which is pumped by a 10 Hz Nd:YAG laser at λ = 532 nm, and are amplified to approximately 20 mJ per pulse each. Once amplified, the pair of pulses is recombined in space with a polarizing beam cube and in time with a delay stage. The recombined pairs of pulses are given opposite circular polarization with a ¼-wave plate, resulting in an optical field that angularly accelerates over the pulse duration. Gas-phase CO2 molecules that are trapped in the optical field through an induced-dipole interaction are angularly accelerated into extreme rotational states.
Fig. 2a shows the spectra of the pair of optical centrifuge pulses. The instantaneous angular frequency of the optical centrifuge trap is given by . Fig. 2b shows the intensity profile of the optical trap as a function of ΩOC, based on the intensity profile of ω1(t). The extent of rotational excitation in the optically centrifuged molecules is controlled by the spectral bandwidth of the optical centrifuge pulse. Optical centrifuge pulses with the full spectral bandwidth can populate CO2 rotational states near J = 360. Here we used reduced bandwidth pulses to preferentially populate rotational states with J ≤ 282. A tunable optical centrifuge trap is realized by removing a selected short wavelength portion of ω1 using a micrometer-controlled beam block.
Fig. 2c shows the geometry of the interaction region for the centrifuged molecules. The optical centrifuge propagation vector OC lies along the z-axis. Optically centrifuged molecules have unidirectional rotational motion within the x–y plane and orientated angular momenta along the z-axis. The IR probe beam is propagated along the x-axis with vector IR. A ½-wave plate controls the polarization of the IR probe, with ∥-polarization parallel to the z-axis (E∥) and ⊥-polarization (E⊥) perpendicular to the z-axis. The optical centrifuge beam is focused to a beam waist of ω0 = 52 μm. The IR beam is focused to a beam waist of 230 μm.
Individual rotational states of the centrifuged CO2 molecules are measured with high-resolution transient IR absorption spectroscopy from a λ = 4.5 μm quantum cascade laser (Daylight Solutions). The quantum cascade laser has a resolution of Δν ≤ 0.0002 cm−1. The IR output frequency is controlled via active feedback using a scanning Fabry–Perot etalon and a lock-in amplifier to tune over the CO2 transitions, while the wavelength is monitored on a wavemeter (Bristol Instruments) with a resolution of Δν ≤ 0.0004 cm−1. Inside the gas cell, the IR beam is overlapped with the optical centrifuge 11 times using an IR multipass (shown in the Fig. 1 inset) to increase signal-to-noise. Spectral features were measured with a 200 ns rise time InSb detector, while dynamics measurements used an InSb detector with a 70 ns rise time. The CO2 pressure used in this experiment is 2.95 Torr based on an equilibrium IR absorption measurement, corresponding to a gas kinetic collision time of ∼30 ns. Transient IR measurements at t = 120 ns represent CO2 molecules that have undergone approximately four gas-kinetic collisions.
Erot = BJ(J + 1) − D(J(J + 1))2 + H(J(J + 1))3 | (1) |
The optical centrifuge was used to populate CO2 rotational states with J values well above a 300 K distribution and transient IR absorption spectroscopy was used to identify the line center IR transition frequencies. CO2 in the ground vibrational state has only even J values because of nuclear spin statistics. Line-center transition frequencies (νobs) were determined by collecting transient absorption signals at discrete IR frequencies in steps of δν ≤ 0.003 cm−1. Fig. 3a shows the set of transient absorption signals used to identify the CO2 R282 transition near 2315 cm−1. Fig. 3b shows the fractional absorption intensities at t = 220 ns after the optical centrifuge pulse as a function of IR frequency. The data are fitted with a Gaussian function and the fit is used to identify the line-center transition frequency. For this transition, line-center frequency is νobs = 2314.9167 ± 0.003 cm−1. Such scans were collected for the 30 R-branch transitions reported here.
The observed IR R-branch transition frequencies for J ≤ 218 are within 0.01 cm−1 of the frequencies predicted by the 3rd order polynomial expansion in J(J + 1). The predicted values are denoted νpred. However, starting at R220 and going to R256, the observed transition frequencies deviate exponentially from the predicted values. The observed deviation for the R256 transition is 3.2 cm−1. An exponential fit of the spectral deviations predicted a frequency deviation for the R258 transition of 4.43 cm−1 but no transition was observed in this spectral region.
To identify the observed deviations for additional transitions, the predicted values from the 3rd order polynomial expansion were compared to transition frequencies from the effective Hamiltonian calculations CDSD4000.54 The transition frequency differences of νcalc and νpred are presented in Fig. 4a, showing that two spectral perturbations are present in the R-branch transitions near R258 and R292. Using the calculated frequencies as a guide, we observed the R268 to R282 transitions. The R258 to R266 transitions remain elusive, in part because of interference with strongly absorbing ambient CO2 lines. A comparison of the observed transition frequencies relative to the predicted values is shown in Fig. 4b. There is excellent agreement between the observed and calculated transition frequencies. Table 1 lists the transition frequencies measured in this study, along with calculated frequencies, deviations, and rotational energies Erot from the effective Hamiltonian calculations.
J state | ν obs (±0.003 cm−1) | ν calc (cm−1) | Δ(νobs − νcalc) (cm−1) | E rot (cm−1) |
---|---|---|---|---|
a From Heff calculations by Tashkun and Perevalov in the CDSD4000.54 | ||||
186 | 2383.921 | 2383.914 | 0.008 | 13411.817 |
188 | 2383.070 | 2383.065 | 0.005 | 13697.486 |
194 | 2380.360 | 2380.365 | −0.005 | 14571.852 |
196 | 2379.421 | 2379.413 | 0.008 | 14869.078 |
198 | 2378.443 | 2378.435 | 0.008 | 15169.181 |
200 | 2377.441 | 2377.432 | 0.009 | 15472.156 |
202 | 2376.413 | 2376.402 | 0.011 | 15777.998 |
204 | 2375.359 | 2375.347 | 0.012 | 16086.702 |
212 | 2370.885 | 2370.87 | 0.015 | 17350.035 |
214 | 2369.704 | 2369.688 | 0.016 | 17672.972 |
218 | 2367.267 | 2367.247 | 0.020 | 18327.332 |
222 | 2364.730 | 2364.707 | 0.023 | 18992.972 |
226 | 2362.095 | 2362.071 | 0.024 | 19669.847 |
228 | 2360.743 | 2360.717 | 0.026 | 20012.485 |
234 | 2356.552 | 2356.522 | 0.030 | 21057.125 |
236 | 2355.113 | 2355.083 | 0.030 | 21410.896 |
238 | 2353.658 | 2353.625 | 0.033 | 21767.435 |
240 | 2352.186 | 2352.151 | 0.035 | 22126.738 |
244 | 2349.208 | 2349.172 | 0.036 | 22853.609 |
246 | 2347.712 | 2347.682 | 0.030 | 23221.164 |
250 | 2344.796 | 2344.78 | 0.016 | 23964.486 |
256 | 2341.879 | 2341.775 | 0.104 | 25099.903 |
268 | 2327.177 | 2327.168 | 0.009 | 27443.674 |
270 | 2325.475 | 2325.466 | 0.009 | 27843.681 |
272 | 2323.728 | 2323.73 | −0.002 | 28246.351 |
274 | 2321.958 | 2321.979 | −0.021 | 28651.677 |
276 | 2320.175 | 2320.227 | −0.052 | 29059.651 |
278 | 2318.398 | 2318.494 | −0.096 | 29470.268 |
280 | 2316.637 | 2316.803 | −0.166 | 29883.520 |
282 | 2314.920 | 2315.188 | −0.268 | 30299.400 |
To understand the vibrational states responsible for the perturbations, the CDSD4000 (0000) and (0001) energies were fit to a 3rd order polynomial expansion in J(J + 1) with J ≤ 300.54 The results were used to predict unperturbed high-J energies for both vibrational states. The deviations from the calculated and predicted energies are shown in Fig. 5. The (0000) state of CO2 shows no J-dependent structure, but the (0001) state shows distinct energy deviations near J = 259 and 293. This comparison shows that the (0001) state is coupled to other CO2 rovibrational states with energies near (0001) J = 259 and 293.
Fig. 5 Differences of calculated and predicted energies for the (a) (0000) and (b) (0001) states of CO2. |
The Heff model accounts for three vibrational states of CO2 that interact with the (0001) state.54 The interacting vibrational states are (0330), (1110(1)), and (1110(2)) where (1) and (2) identify Fermi-mixed states. Based on this model, energy level crossings are present near the J = 259 rotational states of the (0001) and (0330) vibrational states and near the J = 291 rotational state of the (0001) and (1110(1)) vibrational states. The (1110(2)) state interacts indirectly through the (0330) state. Based on the effective Hamiltonian, the authors identify Fermi, Fermi + L-type, and Coriolis coupling as being responsible for the perturbations. While the perturbations affect the IR transition frequencies, they do not affect the rotational energies of the CO2 (0000) vibrational state, which is the topic of the dynamics measurements described in the next section.
Fig. 7a shows transient absorption signals I⊥ and I∥ in black for CO2J = 244 and 268. The rise times of the transient signals are detector limited, showing that these states are prepared directly by the optical centrifuge pulse. Approximately 3 collisions having occurred at t = 100 ns, based on the average collision time of 29 ns at a pressure of 2.95 Torr. The ⊥-polarized data from t = 120 ns to 1 μs were fit to a 3-parameter exponential decay function, and the ∥-polarized data over the same time range were fit to a linear function. The results are shown in red in Fig. 7a. The fit parameters were used to determine the time-dependent population and alignment values, shown in Fig. 7b and c, based on the expressions above. The optical centrifuge pulses were adjusted to optimize population in the J = 268 state. Evidence of the population inversion for the centrifuged molecules is evident with the initial population for J = 268 twice that of J = 244. The ⊥-polarized signals are an order of magnitude larger than the ∥-polarized signals, indicating a significant polarization in both states. These features are representative of the transient absorption signals for the CO2J = 244–282 states.
We now consider the J-dependence of the alignment in the x–y plane of optically centrifuged CO2 molecules with J = 244–282. Fig. 8a shows that the alignment moments at t = 120 ns are close to A(2)0 = 1.5, indicating a high degree of alignment relative to the excitation polarization of the optical centrifuge. We performed additional measurements for the CO2J = 76–100 states to check for evidence of alignment in the collision products. The J = 76–100 states are not populated directly in the optical centrifuge pulse; instead, they are populated through collisions of optically centrifuged molecules with the thermal bath. Fig. 8b shows the fractional absorption intensities I⊥ and I∥ for J = 76. These signals have rise times that are 6-fold longer than those for the J = 244–282 states, consistent with collision-induced population growth. The I⊥ and I∥ signals for J = 76–100 have similar intensities and decay rates. The alignment moment values for these states are close to zero, as shown in Fig. 8c, indicating that the collisions involve a broad range of scattering angles.
(2) |
Relaxation rate constants were determined from decay lifetimes for individual rotational states with J = 244–282 by fitting the population to exponential decay functions. Fig. 9a shows the population decay (in red) for J = 270 along with the fit results in black. The J-dependent rate constants for these states are shown in Fig. 9b, with values of kJout = (4.3 ± 0.4) × 10−11 cm3 molecule−1 s−1 for J = 280 and kJout = (3.0 ± 0.3) × 10−11 cm3 molecule−1 s−1 for J = 244. These rate constants are an order of magnitude smaller than the Lennard-Jones collision rate constant kLJ = 3.7 × 10−10 cm3 molecule−1 s−1, showing that population loss does not occur on every collision for the J = 244–282 states.
Relaxation kinetics were also measured for several states with J = 76–100 to better characterize the overall collisional relaxation. These states are populated by the collisional cascade of the optically centrifuged molecules. The population decay for J = 76 is shown (in blue) in Fig. 9c, along with an exponential fit. Fig. 9d shows that the rate constants for J = 76, 84, 92 and 100 are significantly smaller than the values for J = 244–282, consistent with the overall relaxation process.
Without state-to-state cross sections, either from experiments or calculations, it is not possible to assess how the observed relaxation rate constants compare to predicted values from statistical models that have been used for lower J states. At present, the reported rate constants for CO2 relaxation are benchmarks for future studies. In a related study on CO(J) relaxation, however, we used master equation simulations to calculate the relaxation kinetics and rotational distributions of optically centrifuged CO with J ≤ 80, based on extrapolated rate constants from state-to-state measurements for J = 0–29.27,28 The simulations successfully reproduced the qualitative features of the experimental results, with the key exception that the calculated relaxation rate constants were 5 times larger than those measured in the experiments. The statistical model for the CO system considered the increasing energy gaps for high J states, but the effects of large angular momentum and short rotational periods were not included.
Here we consider the time-dependence of the molecular alignment in the x–y plane. The A(2)0 values for CO2J = 270 are shown in Fig. 9e as a function of time. Loss of alignment in the centrifuged molecules occurs through collisions and pseudo-first order kinetics were used to describe the time-evolution of A(2)0, such that (A(2)0)t/(A(2)0)t=0 = e−kA[M]t. The decay profiles are linear for the J = 244–282 states and rate constants were determined using (A(2)0)t/(A(2)0)t=0 = −kA[M]t.
The resulting kA values are near kA = 2 × 10−11 cm3 molecule−1 s−1 for J = 244–280. However, kA for J = 282 is four times smaller, with kA = 4.7 × 10−12 cm3 molecule−1 s−1. This result shows that while collisions reduce the alignment of the optically centrifuged molecules, the effect is less apparent for the highest populated state. The alignment decay for the J = 282 state results from collisions that reorient a subset of molecules originally prepared by the optical centrifuge, but does not include incoming signal from the relaxation of higher J states. There are no J states populated above J = 282 and up-collisions from lower J are unlikely given the relatively large energy gaps of 420 cm−1 or more. In this way, the J = 282 state is unique among the set of states investigated here.
We also find that the alignment decay rates are a factor of two times smaller than the population decay rates. This set of rate constants for population and alignment decay show that collisions of optically centrifuged CO2 are more likely to reduce the J quantum number of the centrifuged molecules than to induce large changes in the mJ projection quantum number (without accompanying changes in J) that lead to randomization of their spatial orientations.
A(2)0 = 〈2P2(cosθ)〉 = (3〈cos2θ〉 − 1), | (3) |
The Doppler profiles show that collisions of optically centrifuged CO2 molecules with thermal bath molecules are impulsive, leading to rotational excitation of the bath molecules and translational energy in both collision partners. The high degree of alignment of the centrifuged molecules in the x–y plane results in collision products that are scattered predominantly outward from the z-axis with cylindrical symmetry and the most effective collisions for energy transfer will have relative velocities that are in the x–y plane. Doppler profiles measured along the x-axis therefore characterize the maximum recoil velocities of the scattered molecules.
The difference in the Doppler broadening for J = 92 and J = 280 is an indication that the optically centrifuged molecules have lower initial velocities than the thermal bath molecules. Fig. 11c considers the velocity vectors before and after a head-on collision with a b = 0 impact parameter, in the lab frame and the center of mass (cm) frame. The lab frame velocities before the collision are v1 for the centrifuged molecule and v2 for the thermal bath molecule; the corresponding cm velocities are w1 and w2. After the collision, the cm velocities w3 (for the centrifuged molecule) and w4 (for the bath molecule) are equal and opposite. However, in the lab frame, v3 < v4 and the high-J collision products have narrower Doppler profiles.
Previous studies from our laboratory provide additional support that the centrifuged molecules have sub-thermal velocity distributions.21 In depletion measurements of CO2 following the optical centrifuge pulse, the J = 0 and J = 36 states have narrow Doppler profiles that correspond to translational temperatures of 130 and 140 K, respectively, showing that a translationally cool set of molecules are centrifuged to higher states. Further translational cooling may also take place during the optical centrifuge pulse through electrostriction in the focused laser beam. This phenomenon will be the focus of future investigations.
The nearly nascent distribution of the centrifuged molecules with J = 244–282 peaks at J = 268, and the centrifuged molecules have a high degree of alignment with the initial excitation polarization. The alignment moments for this set of optically centrifuged states are near A(2)0 = 1.5, which is near the maximum value of A(2)0 = 2. The centrifuged molecules are prepared with average angular momentum projection mJ quantum numbers that are within 90% of the J vector based on measured values of A(2)0. Through collisions, the centrifuged molecules gradually lose their alignment, but the anisotropy decay is slower than the population decay. Collisions lead to rotational excitation of thermal bath molecules, and transient measurements on the J = 76–100 states show that these collision products have a low degree of alignment.
The results presented here provide a clear picture of the rotation-to-translation energy transfer mechanism for molecules with large amounts of angular momentum. Our results reveal that collisions of centrifuged molecules impart rotational energy to thermal bath molecules, along with translational energy to both collision partners. A comparison of J-dependent Doppler profiles indicates that the centrifuged molecules initially have translational energies that are less than the thermal average, showing that the centrifuged molecules are prepared with relatively low kinetic energies. Collisional relaxation rate constants for the J = 244–282 states are near one-tenth the gas kinetic collision rate constant. Reduced collisional energy transfer rates are associated with increasing rotational energy gaps. Without state-to-state rate information for CO2–CO2 collisions, it is not possible to compare the rate results to statistical models. However, statistical models that work well for low-J CO–CO collisions overestimate the observed relaxation rate constants for high-J molecules, suggesting that angular momentum effects must also be considered. Future studies will take advantage of the tunability in the optical centrifuge to investigate CO2 collisional phenomena over a broader range of rotational states in order to elucidate J-dependent behavior in more detail.
These studies demonstrate the power of modern optical methods for generating and investigating molecules with large amount of controllable angular momentum. New IR spectral lines of high-J molecules will contribute to more accurate description of high energy molecules, and enable future investigations of high energy environments and remote atmospheres. The results reported here contribute to a growing body of work that highlight the impact of large amounts of angular momentum on collisional energy transfer, and the underlying principles that control molecular collisions.
This journal is © The Royal Society of Chemistry 2024 |