Chatura
Perera
a,
Ethan
Ross
a,
Yang
Liu
b,
Hua
Guo
b and
Arthur G.
Suits
*a
aDepartment of Chemistry, University of Missouri, Columbia, MO 65211, USA. E-mail: suitsa@missouri.edu
bDepartment of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, NM 87131, USA
First published on 17th December 2024
We present state-to-state differential cross sections for rotationally inelastic collisions of vibrationally excited NO X2Π (ν = 9) with Ar using a near-counterpropagating molecular beam geometry. These were obtained using the stimulated emission pumping technique coupled with velocity map imaging. Collision energies well over ∼1 eV were achieved and rotational excitations up to ∼Δj = 60 recorded for the first time for inelastic collisions. This allowed us to investigate scattering of a diatomic molecule in a 2Π state which is initially well described by Hund's case (a) into final states well described by Hund's case (b) as the rotational level splitting becomes larger than the spin–orbit splitting. Differential cross sections for both parity-changing and parity-conserving collisions exhibit very similar structures at the high collision energies. Quantum scattering calculations have been carried out to obtain approximate integral cross sections, which confirm the high rotational excitation. These studies will take the arena of rotationally inelastic collisions to a new regime while providing insight into dynamics under extreme non-equilibrium conditions. Furthermore, these present a unique challenge to both quantum and quasi-classical scattering calculations to validate the methods and the potential energy surfaces used to assess their applicability under extreme conditions.
NO in its ground electronic state has a 1σ22σ23σ24σ21π45σ22π1 electronic configuration which gives rise to the 2Π term. The lowest rotational states of NO(X2Π) are well described by Hund's case (a) where the projection of the electronic orbital momentum and the spin angular momentum on the internuclear axis couple together. This results in two spin–orbit states with a lower energy 2Π1/2 state indicated as Ω = 1/2 and the higher energy 2Π3/2 state with Ω = 3/2. Furthermore, the rotational manifolds that arise from states 2Π1/2 and 2Π3/2 are denoted as F1 and F2 respectively. Each rotational level further splits into two closely spaced Λ-doublet components that differ in their total parities. This is due to interaction between electronic orbital and rotational angular momenta. The symmetric and antisymmetric combination of the +Ω and −Ω wavefunctions are indicated by e Λ-doublet and f Λ-doublet sublevels respectively. The total parity of the Λ-doublet sublevels is then given by ε(−1)j−1/2 where ε is +1 and −1 for e and f, respectively. The f levels are slightly higher in energy than e levels. When a closed shell target approaches, the degeneracy of NO breaks yielding two potential energy surfaces (PESs) of A′ or A′′ symmetry.15 The sum potential (A′ + A′′)/2 is said to govern the spin–orbit conserving collisions while the difference potential (A′ − A′′)/2 induces the spin–orbit changing collisions. Studies involving spin–orbit changing collisions are less common due to their low collisional cross sections compared to spin–orbit conserving collisions.16–18 Furthermore, collisions considering initial and final state parities have been very useful in testing and validating the accuracy of calculated PESs. Parity conserving collisions show multiple peaks in the differential cross sections (DCSs) due to quantum interference between trajectories at two different NO orientation angles which scatters into the same final state.19 Generally only a single peak is observed in the DCS for parity changing collisions. It is also shown that the inelastic cross sections are unchanged with respect to the e/f labelling of initial and final states for a number of molecules well described in the Hund's case (a) limit.20,21 In other words, cross sections will be identical for (je → j′e) and (jf → j′f) and similarly for (je → j′f) and (jf → j′e). Depending on the symmetric and antisymmetric nature of the electronic wave function with respect to the plane of rotation, Λ-doublet states are also labeled as A′ or A′′. When rotational level splitting becomes comparable or larger than the spin–orbit coupling constant an intermediate case or Hund's case (b) is more suitable in describing the NO in 2Π state.20 The unpaired electron of NO lies in the plane of rotation in the A′ states and perpendicular to that plane in A′′ states.22,23 Dagdigian and co-workers have discussed in detail how e/f symmetry disappears in the high j limit and the observed cross section asymmetry is a quantum mechanical interference effect.24 This further demonstrates that two potential energy surfaces are necessary in interpreting collisions of molecules in the 2Π state with closed shell partners.10,15,22 Studies have been carried out for species with both singly filled π1 orbital and triply filled π3 orbital giving rise to a 2Π state.22,25
State-to-state studies offer the deepest insight into these and although final state detection is usually straightforward with laser-based probes, initial state preparation is not straightforward. Preparing NO in single quantum states is difficult due to its open-shell character and low rotational constant but essential for a detailed comparison with theories. This is very important in revealing fine quantum effects including resonances and stereodynamics. Using supersonically expanded beams and hexapoles, the Stark effect of NO has been used to achieve greater detail in scattering experiments by preparing NO in its low-field-seeking j = 0.5, −f state for stereodynamics investigations.16,17,19,26,27 Using this, fully quantum state-resolved differential cross sections (DCSs) have been observed for NO–Ar inelastic collisions in a crossed molecular beam experiment over a modest range of collision energies and orientation dependence of the scattering revealed. These results were accompanied by quantum mechanical close-coupling (QMCC) scattering calculations employing coupled-cluster PESs.19,28 Experiments and QMCC calculations were able to capture the oscillations in the DCSs that appeared in both parity-conserving and parity-changing collisions validating the accuracy of calculated PESs. Stark deceleration has also been further used to reveal rich details in scattering including diffraction oscillations in the DCSs by achieving very narrow velocity spreads for state-selected reactants.29–31 These oscillations are extremely sensitive to the effective size of the target.
However, all these studies have been focused on inelastic collisions of the vibrational ground level of NO. Even though QMCC calculations successfully reproduced the experimental results, the accuracy of the NO–Ar PESs was largely untested for the NO molecule stretched in high vibrational levels and for low energy collisions that are dominated by the attractive potential. We recently extended such studies to include DCSs for rotationally inelastic and spin–orbit changing collisions of highly vibrationally excited NO (ν = 10) molecules with closed shell partners Ar and Ne, and were able to test the accuracy of both attractive and repulsive parts of both MRCI (multi-reference configuration interaction) and CCSD(T) (coupled cluster singles, doubles, and perturbative triples) PESs using QMCC scattering calculations.32–37 We utilized the stimulated emission pumping (SEP) technique which was pioneered by Field and coworkers and used by Wodtke and co-workers to prepare NO in highly excited vibrational levels to investigate rovibrational relaxation rates and surface scattering dynamics.38–40 We applied SEP to prepare NO (ν = 10) in a single rovibrational level in both spin-manifolds coupled with velocity map imaging (VMI) to measure DCSs for inelastic collisions of NO in a near-copropagating (4° crossing angle) crossed beam apparatus. Fully quantum state-resolved state-to-state DCSs were measured for rotationally inelastic collisions of vibrationally excited NO (ν = 10) at a broad range of collision energies well above room temperature down to the cold limit (∼2 K), and for spin–orbit conserving and spin–orbit changing collisions.
The utilization of counterpropagating beam geometries in crossed beam experiments is relatively unconventional. In the early 1990s, Henning Meyer successfully employed a counterpropagating pulsed beam scattering technique in pioneering stereodynamics studies to enable measurements of angle-resolved state-specific product distributions using time-of-flight detection of scattered molecules.41–44 State-to-state inelastic DCSs for NO–Ne and NO–Ar collisions, obtained by combining the Stark deceleration and VMI techniques were later studied using counter-propagating beams by van de Meerakker and co-workers.45 All the collisions were carried out at collision energies up to ∼1300 cm−1. In the current work we present state-to-state DCSs obtained for rotationally inelastic collisions between vibrationally excited NO (v = 9) and Ar at two very high collision energies ∼1 eV utilizing a near-counterpropagating beam geometry. We have recently reported use of this experimental geometry in a study of spin–orbit changing collisions over a very broad range of collision energies, although that investigation was confined to low j product states.46 Here we applied the same SEP technique as mentioned above in our recent work in preparing vibrationally excited NO (ν = 9). In this study we obtained DCSs for very high rotational excitation, up to ∼j = 60, which has not been reported before for state-to-state scattering.
A mixture consisting of 15% NO seeded in H2 with a backing pressure maintained at 3.5 bar was used to prepare vibrationally excited NO molecules. The supersonically expanded NO beam was introduced into the chamber through a 0.5 mm nozzle in the new pulsed valve. A pure Ar beam with a backing pressure of 9 bar was expanded through the side beam of the dual pulsed valve positioned off axis on the opposite side of the chamber giving an intersection angle of 176° in the crossing region. Both beams were collimated by skimmers with 1 mm aperture which are 2 cm downstream from valve nozzles. All lasers and both valves were operated at 10 Hz and nominal pulse widths of 25 μs and 20 μs were used respectively for NO and Ar beams. SEP lasers which were loosely focused intersected the molecular beam perpendicularly 1 cm upstream from the interaction region to prepare vibrationally excited NO molecules.
Three tunable dye lasers were used in this experiment. SEP was used to prepare NO molecules in single ro-vibrational states: NO molecules in the beam were electronically excited via the R11 + Q21(0.5) transition of the 2–0 band from the NO X2Π (ν = 0, Ω = 0.5, j = 0.5, +e) to the A2Σ+ (ν = 2, j = 1.5, F1) level by the pump laser. The dump laser was then used to prepare vibrationally excited NO molecules in the spin–orbit ground state via the R11 + Q21(0.5) transition of the 2–9 band by dumping electronically excited NO molecules into the NO X2Π (v = 9, Ω = 0.5, j = 0.5, +e) level. Scattered NO products in the NO X2Π (ν = 9, Ω′ = 0.5, N) level due to rotationally inelastic collisions were detected using a (1 + 1) REMPI scheme around 278 nm via the 3–9 band of the A2Σ+ ← X2Π transition. The ions produced were then focused onto the microchannel plate (MCP) detector coupled to a phosphor screen from our sliced velocity map imaging setup, and the flashes in the phosphor screen were detected by a CCD camera and the positions of the ion spots were centroided and recorded utilizing our NuAcq acquisition software.48,49 The measured state-to-state differential cross sections (DCSs) were converted from the LAB frame to the COM frame using a Monte Carlo (MC) forward convolution simulation.33,50
An accurate comparison with the measured DCSs is challenging for two reasons. First, an appropriate description of highly excited rotational states of NO needs to follow Hund's case (b), which necessitates a transformation of the basis for NO.24 Second, the high rotational excitation in the NO product requires a very large basis. We defer the full quantum scattering calculation of the DCS to a future work.
Fig. 3 shows the parity conserving collisions studied at a collision energy of 6800 cm−1 (0.84 eV). NO molecules prepared in the X2Π (ν = 9, Ω = 0.5, j = 0.5, +e) level following collisions with Ar populate a range of final rotational levels detected using a (1 + 1) REMPI scheme. This initial state of NO is still very much within the description of Hund's case (a). For parity conserving collisions, products were detected in two different rotational levels which are NO X2Π (ν = 9, Ω′ = 0.5, N = 41, +f) and NO X2Π (ν = 9, Ω′ = 0.5, N = 53, +f). These will be referred to as 41f and 53f hereafter. Due to the high velocity of the NO beam resulting in high collision energies, relatively high voltages were applied for the electrostatic lenses in the DC slice VMI step up. Voltages of 3000 V, 2648 V, and 2328 V were applied to the repeller, first extractor lens and the second extractor lens. The third extractor lens was grounded. Sufficient statistics for a scattering image were obtained with 300000 laser shots for a particular rotational level. Images were recorded with 60
000 laser shots and summed together to prevent blurring effects from the valve heating or other variation in experimental conditions. The “raw” experimental center-of-mass DCS was obtained by directly integrating the edge of the scattering image.
As we observed before with near-copropagating beams, a clear symmetry in the scattering image is seen here for near-counterpropagating beams due to the smaller density-to-flux correction compared to a conventional 90° crossing geometry. We have scanned the probe laser delay relative to SEP lasers to capture all the scattered products uniformly during image acquisition. This was necessary to account for the fact that our scattering was not confined to the probe region. The underlying center of mass DCSs were extracted using our previously used Monte-Carlo (MC) simulation code by modifying the code to accommodate the near-counterpropagating geometry. We incorporated all the aspects of the current experiment in this forward-convolution fitting program including measured beam velocity spreads and delays. More details about the fitting program can be found elsewhere.33 As reported in many previous studies, multiple peaks were observed in the DCSs for both 41f and 53f rotational levels where the parity is conserved.
The parity dependence of these collisions was also studied by probing parity changing collisions at 6800 cm−1 as shown in Fig. 4. Here the rotationally inelastic products were detected at NO X2Π (ν = 9, Ω′ = 0.5, N = 41, −e) and NO X2Π (ν = 9, Ω′ = 0.5, N = 51, −e). These will be referred to as 41e and 51e hereafter. We have observed multiple peaks in the DCSs for both parity conserving (Δj = odd, e → f) and parity changing (Δj = odd, e → e) collisions in these high collision energy studies involving very high rotational excitations. In contrast, at a collision energy of 530 cm−1 we observed only a single main peak for parity changing collisions involving vibrationally excited NO, and this is consistent with the results of Brouard and coworkers for low energy collisions of the vibrationless molecule.55 In the present high collision energy results, strong sideways scattering peaks are observed at both 41f and 41e. Sideways scattering peaks were shifted to longer scattering angles all the way to backscattered peaks with the increase in rotational excitation up to 53f and 51e at the same collision energy. The reason to study two different final rotational states of 53f and 51e for parity conserving and changing collisions was to avoid overlapping transitions and background due to the amount of available energy at these collision energies and the presence of other vibrational levels populated by spontaneous emission.
A second experiment was carried out at a higher collision energy of 11200 cm−1 (1.4 eV). For this, 15% NO seeded in H2 collided with a secondary beam which was 15% Ar seeded in H2. This new gas composition yielded an Ar beam velocity of 1485 ms−1. A new mechanical backing pump had to be used in line with the existing pump to facilitate H2 pumping and maintain single collision conditions in carrying out these new studies. Here we observed rotational excitations up ∼Δj = 60 with the large available energy as shown in Fig. 5.
In this case the rotationally inelastic products were detected at NO X2Π (ν = 9, Ω′ = 0.5, N = 59, +f) and NO X2Π (ν = 9, Ω′ = 0.5, N = 55, −e). These will be referred to as 59f and 55e hereafter. Energy taken up by the rotational excitation at Δj = 60 is well over 6000 cm−1. At both collision energies we find both parity-conserving and -changing collisions exhibit multiple peaks in the DCSs.
An intermediate Hund's case or Hund's case (b) provides a better description outside the above-mentioned limits. Ω = 1/2 and Ω = 3/2 wave functions are mixed by the spin-rotation terms when A ≪ BJ. In this case the true molecular wave functions are a linear combination of case (a) wave functions. As mentioned before, in case (a), cross sections for (je → j′e) and (jf → j′f) will be identical and similarly for (je → j′f) and (jf → j′e) independent from the A′ or A′′ PESs. Therefore, equal populations in the Λ doublet levels for all final states will be observed if an experiment is carried out with equal populations in both Λ doublet levels (e and f) of an initial j. The deviation seen from Hund's case (a) to (b) resulting in unequal final state Λ doublets levels at high j limit even with equal initial Λ doublet levels was studied by Dagdigian and co-workers.24 They investigated the quantum formalism for scattering of a diatomic molecule in the 2Π state which is well described by Hund's case (b). They showed that the differing Λ doublet propensities were due to an interference between the terms in the expansion of the two sum and difference potentials of A′ or A′′ reflection symmetry. This describes the interaction between a molecule in 2Π state and a closed-shell partner. They calculated the cross sections for the CH(X2Π)–He system at a collision energy of 1000 cm−1 and the results were compared with crossed molecular beam studies done by Liu and Macdonald.56 Furthermore it is also shown that for interaction of a closed-shell partner with a molecule in the 2Π electronic state resulting from a singly filled π orbital (π1) such as CH(X2Π) and NO(X2Π), the repulsive interaction will be larger for the A′ state. This is because the π1 orbital lies in the triatomic plane for the A′ state and perpendicular to the triatomic plane for the A′′ state. In conclusion the preferentially populated final states are F1f and F2e of A′′ symmetry whose electronic wave function is antisymmetric.
In this study it is clear that the highly rotationally excited products of vibrationally excited NO with Ar fall closer to the Hund's case (b) limit. By detecting the scattered products in both 41f and 41e at 6800 cm−1, also in 59f and 55e at 11200 cm−1 we have looked at the parity dependence on the DCSs of these collisions involving high rotational excitations. All the studies done so far including our recent work of vibrationally excited NO (ν = 10) have shown multiple peaks in the DCSs for parity conserving (Δj = odd, e → f) and a single peak for parity changing (Δj = odd, e → e) collisions. This was also a confirmation that our initial state preparation by the SEP technique does not populate more than one single rovibrational state. Otherwise, this oscillatory structure would not be resolved if more than a single rovibrational level is populated. Brouard and co-workers observed this product parity sensitivity towards the state-to-state DCS for the NO–Ar system.19 They utilized a simple four-path hard shell model originally introduced by McCurdy and Miller and extended by Korsch and Schinke to predict quantum interference structures in DCSs and ICSs.57,58 This model only involves the repulsive part of the PESs which largely governs rotationally inelastic scattering conducted at collision energies corresponding to room temperature or higher. The conclusion of that study was that the parity-conserving DCS is sensitive to both the heteronuclear and homonuclear terms in the interaction potential, whereas the parity-changing collisions are only sensitive to heteronuclear terms. Collisions at two different NO orientation angles interfere to give rise to these oscillations in the DCS as result of homonuclear terms dominating the NO molecule. On the other hand, parity-changing collisions are sensitive only to interferences from different paths involving collisions at the two different ends of the molecule. As NO is a near homonuclear molecule, the phase shifts associated with these different trajectories do not give rise to structures in the DCS at modest collision energies as they are varying slowly. However, in the current study we observe multiple peaks in the DCS for both parity-conserving and parity-changing collisions. This could well be explained by the quantum interferences explained by Dagdigian and co-workers in Hund's case (b).24 Preferentially populating F1f could play a critical role as is seen in almost identical scattering signal intensity and features in the DCSs for parity changing collisions. To gain some insight into these processes, we have adapted the four-path model discussed above. Fig. 6 shows the DCSs obtained from the model for both parity-conserving and parity-changing collisions of NO X2Π (ν = 9, Ω = 0.5, j = 0.5, +e) at a collision energy of 6800 cm−1 resulting in 41f and 41e. This model also predicts multiple peaks in the DCSs for both parity-conserving and parity-changing collisions, suggesting one aspect of this is more rapidly varying phase shifts at this higher collision energy.
Integral cross sections (ICSs) for spin–orbit changing collisions are generally significantly lower compared to spin–orbit conserving collisions.2,4,16–18 This is mainly due to the fact that the difference potential is smaller than the sum potential. Previous experimental and theoretical work involving rotationally inelastic collisions of spin–orbit changing collisions have shown an increase in the ICS with the increase in rotational excitation (jfinal − jinital) for low Δj transitions. This behavior is in complete contrast with the spin–orbit conserving manifold. Some initial quantum calculations carried out with the CS approximation show the same trend for the studies we have carried out here which involve vibrationally excited NO molecules. Fig. 7 shows the calculated ICSs for both spin–orbit conserving and spin–orbit changing collisions arising from the NO X2Π (ν = 9, Ω′′ = 0.5, j′′ = 0.5, +e) initial state. These results clearly show the high rotational excitation of the scattered NO. Furthermore, with the increase in rotational excitations for low Δj transitions, the spin–orbit changing collisions show a clear increase in the ICS while spin–orbit conserving collisions show a gradual decrease with the increasing rotational excitation. Fig. 8 shows the ICSs for vibrationally inelastic scattering into v′ = 8 and 10. These are lower than the vibrationally elastic cross sections by three to four orders of magnitude. As expected, the exoergic process is favored, in this case by roughly a factor of two.
![]() | ||
Fig. 7 Calculated ICSs for spin–orbit conserving (red) and spin–orbit changing (blue) collisions arising from the NO X2Π (v = 9, Ω = 0.5, j = 0.5, +e) initial state. |
This journal is © the Owner Societies 2025 |