Rotational State Distributions Research Articles

Rotational state-to-state transition rate coefficients for H2O + H2O collisions at nonequilibrium conditions

The goal is to develop a database of rate coefficients for rotational state-to-state transitions in H$_ $O + H$_ $O collisions that is suitable for the modeling of energy transfer in nonequilibrium conditions, in which the distribution of rotational states of H$_ $O deviates from local thermodynamic equilibrium. A two-temperature model was employed that assumed that although there is no equilibrium between all possible degrees of freedom in the system, the translational and rotational degrees of freedom can be expected to achieve their own equilibria independently, and that they can be approximately characterized by Boltzmann distributions at two different temperatures, $T_ kin $ and $T_ rot Upon introducing our new parameterization of the collisional rates, taking into account their dependence on both $T_ kin $ and $T_ rot $, we find a change of up to 20<!PCT!> in the H$_ $O rotational level populations for both ortho and para-H$_ $O for the part of the cometary coma where the nonequilibrium regime occurs.

A state-selected continuous wave laser excitation method for determining CO2's rotational state distribution in a supersonic molecular beam.

State-resolved experiments can provide fundamental insight into the mechanisms behind chemical reactions. Here, we describe our methods for characterizing state-resolved experiments probing the outcome of the collision between CO2 molecules and surfaces. We create a molecular beam from a supersonic expansion that passes through an ultra-high vacuum system. The CO2 is vibrationally excited by a continuous wave infrared (IR) laser using rapid adiabatic passage. We attenuate the fractional excitation using a CO2 absorption cell in the IR beam path. We combine Monte Carlo simulations and molecular beam energy measurements to find the initial rotational state distribution of the molecular beam. We find that our pure CO2 beam from a 300K source has a rotational temperature of ∼26K.

Isotope Effects on State-to-State Photodissociation Dynamics of D2S in Its First Absorption Band.

State-to-state photodissociation dynamics of D2S in its first absorption band were explored by utilizing recently developed diabatic potential energy surfaces (PESs). Quantum dynamics calculations, involving the first two strongly coupled 1A″ states, were executed employing a Chebyshev real wavepacket method. The nonadiabatic channel via the conical intersection (CI) is facile, direct, and fast, leading to the production of rotationally and vibrationally cold SD(X̃2Π). The calculated absorption spectrum, product state distributions, and angular distributions are in reasonable agreement with the experimental results, although some discrepancies exist at 193.3 nm. Compared with H2S, there are obvious isotope effects on rotational state distributions for D2S photodissociation in its first absorption band. Moreover, we scrutinize the variation of product state distributions as a function of photon energy and the vibrational mediated photodissociation of the parent molecule. Due to the diverse shapes of the three fundamental vibrational wave functions, photoexcited wavepackets access distinct segments of the upper-state PES, resulting in a disparate absorption spectrum and ro-vibrational distributions via the nonadiabatic transition. This study provides a comprehensive figure of the isotopic effect and wavelength dependence on the photofragmentation behaviors from D2S photodissociation, which should attract more experimental and theoretical attention to this prototypical system.

Time evolution of natural orbitals in abinitio molecular dynamics.

This work combines for the first time abinitio molecular dynamics (AIMD) within the Born-Oppenheimer approximation with a global natural orbital functional (GNOF), an approximate functional of the one-particle reduced density matrix. The most prominent feature of GNOF-AIMD is its ability to display the real-time evolution of natural orbitals, providing detailed information on the time-dependent electronic structure of complex systems and processes, including reactive collisions. The quartet ground-state reaction N(4S) + H2(1Σ) → NH(3Σ) + H(2S) is taken as a validation test. Collision energy influences on integral cross sections for different initial rovibrational states of H2 and rotational-state distributions of the NH product are discussed, showing a good agreement with previous high-quality theoretical results.

Automated potential energy surface development and comprehensive dynamics for the F + CH3NH2 reaction.

This work is an extensive investigation of the F + CH3NH2 reaction dynamics using a newly-developed potential energy surface (PES). The full-dimensional spin-orbit (SO) corrected (MRCI+Q/aug-cc-pwCVDZ) PES is developed by the Robosurfer program package and the ManyHF method is used in order to fix the Hartree-Fock (HF) convergence issues in the entrance channel. On the surface, retrieved by the fitting of the iteratively extended set of the ManyHF-CCSD(T)-F12a/triple-zeta-quality and SO-corrected energy points, quasi-classical trajectory (QCT) simulations are run. By analyzing the opacity functions and integral cross sections (ICSs) of six reaction channels, the dynamics of the two most reactive hydrogen-abstraction reactions resulting in HF + CH2NH2/CH3NH products are selected for a thorough examination. Despite the statistically and thermodynamically expected results, the kinetically preferred amino hydrogen-abstraction is the dominant mechanism at low collision energies. The initial attack angle and scattering angle distributions are investigated as well. The post-reaction energy distributions show that the collision energy mostly converts into the translational energy of the products, while the reaction energy excites the vibration of the products. The computed vibrationally resolved rotational distributions and vibrational state distributions of the HF product are compared to experimental data, and the theory and experiment are found to be in good agreement.

Isotope Effect and Heavy-Light-Heavy Reactivity Oscillation in the Cl + CHD3/CHT3 Reaction.

Previous experiments and theories have shown the existence of heavy-light-heavy (HLH) reactivity oscillation in the Cl + CH4 reaction and anticipated that similar oscillations should exist in many HLH reactions involving polyatomic reagents. However, the total reaction probabilities for the Cl + CHD3 → HCl + CD3 reaction exhibit only a step-like feature, and the total reaction probabilities for Cl + CHT3 → HCl + CT3 do not show any structure at all. Here, we report seven-dimensional state-to-state quantum dynamics studies for this reaction on the FI-NN PES, and we demonstrate that HLH reactivity oscillations also exist in these two reactions, manifesting as peaks in the reaction probabilities for low product rotational states. These oscillations, however, are obscured in the total reaction probability because of the higher excitation of j ≥ 2 product rotational states. Furthermore, the isotope replacement of nonreactive hydrogen with deuterium and tritium significantly enhances reactivity at collision energies above 0.112 eV, indicating an inverse secondary isotope effect on the probabilities, which is proved to be also caused by HLH mass combination. We also demonstrate that the highly rotational excitation of CHD3 substantially enhances reactivity and the HLH oscillations, similar to HLH triatomic reactions. These observations are completely different from those in the H + CHD3 reaction, which is also a late-barrier reaction. Therefore, the HLH mass combination is very important, which affects not only the reactivity oscillation but also the amplitude and product rotational state distribution and makes the initial rotation excitation play a pivotal role in the reaction.

Six-dimensional quantum dynamics of an Eley-Rideal reaction between gaseous and adsorbed hydrogen atoms on Cu(111).

In the form of direct abstraction of a surface adsorbate by a gaseous projectile, the Eley-Rideal (ER) reaction at the gas-surface interface manifests interesting dynamics. Unfortunately, high-dimensional quantum dynamical (QD) studies for ER reactions remain very challenging, which demands a large configuration space and the coordinate transformation of wavefunctions. Here, we report the first six-dimensional (6D) fully coupled quantum scattering method for studying the ER reaction between gas phase H(D) atoms and adsorbed D(H) atoms on a rigid Cu(111) surface. Reaction probabilities and product rovibrational state distributions obtained by this 6D model are found to be quite different from those obtained by reduced-dimensional QD models, demonstrating the high-dimensional nature of the ER reaction. Using two distinct potential energy surfaces (PESs), we further discuss the influence of the PES on the calculated product vibrational and rotational state distributions, in comparison with experimental results. The lateral corrugation and the exothermicity of the PES are found to play a critical role in controlling the energy disposal in the ER reaction.

Stereodynamic control of nonadiabatic processes in low-energy Be+(2P) + H2 (v = 0, j = 2) collisions.

Controlling the relative arrangement of colliding molecules is crucial for determining the dynamical outcomes of chemical processes and has emerged as a hot spot of experimental research. Here, the quantum scattering calculations are conducted to investigate the stereodynamic control in collisions between Be+(2P) and H2 (v = 0, j = 2), which undergo nonadiabatic transitions to the electronic ground state. Stereodynamic preparation is achieved by controlling the initial alignment of the H2 bond axis relative to the scattering frame. For product BeH+ in the reactive process, the differential cross sections (DCSs) are significantly enhanced in the forward and sideways hemispheres when the alignment angle β is 60°. For the product H2 in the quenching channel, the β = 0° preparation can result in a more than one-fold increase in the DCS at a polar scattering angle of 0°. Furthermore, varying the alignment angle β also has noteworthy effects on the rotational-state distributions of BeH+ products. Specifically, β = 0° preparation can induce the disappearance of the bimodal distribution of rotational states at a collision energy of 0.05 eV.

Molecular Simulation for Atmospheric Reactions: Non-Equilibrium Dynamics, Roaming, and Glycolaldehyde Formation following Photoinduced Decomposition of syn-Acetaldehyde Oxide.

The decomposition dynamics of vibrationally excited syn-CH3CHOO to form vinoxy + hydroxyl (CH2CHO + OH) radicals or to recombine to form glycolaldehyde (CH2OHCHO) are characterized using statistically significant numbers of molecular dynamics simulations using a full-dimensional neural-network-based potential energy surface at the CASPT2 level of theory. The computed final OH-translational and rotational state distributions agree well with experiments and probe the still unknown O-O bond strength DeOO for which best values from 22 to 25 kcal/mol are found. OH-elimination rates are consistent with experiments and do not vary appreciably with DeOO due to the non-equilibrium nature of the process. In addition to the OH-elimination pathway, OH roaming is observed following O-O scission, which leads to glycolaldehyde formation on the picosecond time scale. Together with recent work involving the methyl-ethyl-substituted Criegee intermediate, we conclude that OH roaming is a general pathway to be included in molecular-level modeling of atmospheric processes. This work demonstrates that atomistic simulations with machine-learned energy functions provide a viable route for exploring the chemistry and reaction dynamics of atmospheric reactions.

Enhancing reactivity of SiO+ ions by controlled excitation to extreme rotational states.

Optical pumping of molecules provides unique opportunities for control of chemical reactions at a wide range of rotational energies. This work reports a chemical reaction with extreme rotational excitation of a reactant and its kinetic characterization. We investigate the chemical reactivity for the hydrogen abstraction reaction SiO+ + H2 → SiOH+ + H in an ion trap. The SiO+ cations are prepared in a narrow rotational state distribution, including super-rotor states with rotational quantum number (j) as high as 170, using a broad-band optical pumping method. We show that the super-rotor states of SiO+ substantially enhance the reaction rate, a trend reproduced by complementary theoretical studies. We reveal the mechanism for the rotational enhancement of the reactivity to be a strong coupling of the SiO+ rotational mode with the reaction coordinate at the transition state on the dominant dynamical pathway.

Rotationally Inelastic Scattering Dynamics of NO from Ag(111): Influence of Interaction Potentials

State-to-state molecular scattering from surfaces represents a sensitive probe of the molecule–surface interaction potential. For NO scattering from Ag(111), a benchmark system for studying gas–surface energy transfer, abundant experimental data have provided us a great testing ground to validate the accuracy of first-principles description of molecule–surface interactions. Here, we use an accurate neural network fitting method to construct two high-dimensional potential energy surfaces (PESs) based on density functional theory (DFT) single points using PBE and revPBE functionals, respectively, for describing the rotationally inelastic scattering of NO from Ag(111). Quasi-classical trajectory calculations based on the two PESs have been compared with a diversity of experimental data and the influence of scattering angle, initial molecular orientation, and incidence energy on the rotational state distribution is discussed. It is found in general that revPBE-based results agree better than PBE-based ones with the measured scattering angle and final rotational state distributions, capturing reasonably the rotational rainbow feature at a small incidence energy and the steric effect. But neither PES reproduces all experimental data quantitatively. The different performance of the two PESs can be attributed to their difference in the attractiveness, anisotropy, and corrugation in the entrance channel. Our results suggest the need of further improvement of the accuracy of DFT with commonly used density functionals toward a chemically accurate description of molecule–surface interactions.

Photodissociation dynamics of H2S+ near 325 nm

We study the photodissociation dynamics of the hydrogen sulfide cations (H2S+) using the time-sliced velocity map ion imaging (VMI) technique and high-accuracy calculations. High-resolution ion images of the S+(4S) products were measured at four photolysis wavelengths of 325.158, 325.200, 325.243, 325.307 nm, which correspond to the excitation to the A2A1(0,13,0) K=1 state of H2S+. Rotational state-resolved total kinetic energy releases and angular distributions have been derived as a function of the photolysis wavelengths. Notably, photolysis wavelength dependent product rotational state and anisotropy parameter distributions have been clearly observed. Full-dimensional potential energy surface characterization suggests that nonadiabatic coupling between A2A1 and B2B2 states at C2v configurations, as well as relaxation of the symmetry to Cs in the conical intersection region between the two states, plays a key role in the photodissociation process.

Velocity-Selected Rotational State Distributions of Nitric Oxide Scattered off Graphene Revealed by Surface-Velocity Map Imaging.

We report velocity-dependent internal energy distributions of nitric oxide molecules, NO, scattered off graphene supported on gold to further explore the dynamics of the collision process between NO radicals and graphene. These experiments were performed by directing a molecular beam of NO onto graphene in a surface-velocity map imaging setup, which allowed us to record internal energy distributions of the NO radicals as a function of their velocity. We do not observe bond formation but (1) major contributions from direct inelastic scattering and (2) a smaller trapping-desorption component where some physisorbed NO molecules have residence times on the order of microseconds. This is in agreement with our classical molecular dynamics simulations which also observe a small proportion of two- and multi-bounce collisions events but likewise a small proportion of NO radicals trapped at the surface for the entire length of the molecular dynamics simulations (a few picoseconds). Despite a collision energy of 0.31 eV, which would be sufficient to populate NO(v = 1), we do not detect vibrationally excited nitric oxide.

Theoretical study on photodissociation dynamics of vibrational excited states of H2S in the first absorption band

The photodissociation quantum dynamics for the first absorption band of H2S in different initial vibrational states have been investigated using Chebyshev real wave packet method. Because of the difference of the wave functions for the initial vibrational states, the calculated absorption spectra and the distributions of vibrational and rotational state of the products display different dynamic characteristics. The width and peak position of the absorption spectra for initial stretching excited states (1,0,0) and (0,0,1) are different from that of the vibrational ground state, while the (0,1,0) vibrational state has two almost equally high peaks in its absorption spectrum because of the change of wave function in angular coordinate. The product vibrational state distribution for (0,1,0) initial state weakly depends on the excitation energy and is dominated by the products of v=0, but SH(v=1) fragment is dominant at lower energy for (1,0,0) and (0,0,1) vibrational states. The rotational state distributions of products are very cold with the peak at j=l for these four states and weakly depend on the total energy. Besides, the rotational state distribution from (0,1,0) vibrational state displays strong oscillation, and its anisotropic parameter with rotational quantum numbers is also different from that of the other three vibrational states.

Ozone Photodissociation in the Singlet Channel at 226 nm.

We report the rotational state distribution and vector correlations of the O2(a 1Δg, v = 0) fragments arising from the 226 nm photodissociation of jet-cooled O3. Consistent with previously reported trends, the rotational distribution is shifted to higher rotational states with decreasing wavelength. We observe highly suppressed odd rotational state populations due to a strong Λ-doublet propensity. The measured rotational distribution is in agreement with classical trajectory calculations for the v = 0 products, although the distribution is slightly narrower than predicted. The spatial anisotropy follows the previously observed trend of decreasing β with increasing photon energy with β = 0.72 ± 0.14 for v = 0, j = 38. As expected for a triatomic molecule, the v-j correlation is consistent with v perpendicular to j, but the measured correlation is nonlimiting due, in part, to rotational and translational depolarization. The j-dependent line width of the O2(a 1Δg) REMPI spectrum is also discussed in connection with the lifetime of the resonant O2(d 1Πg) state due to predissociation via the II 1Πg valence state.

NO (A) Rotational State Distributions from Photodissociation of the N2-NO Complex.

We have recorded the resonance-enhanced multiphoton ionization spectrum for NO (A) products from photodissociation of the N2-NO complex. We made measurements at excitation energies ranging from 28 to 758 cm-1 above the threshold to produce NO (A) + N2 (X) products, and the resulting spectra reveal the NO (A) rotational states formed during dissociation, allowing us to determine the rotational state distribution. At the lowest available energies, 28 and 50 cm-1 above threshold, we observed contributions from NO (A) rotational states that exceed the available energy and must originate from excitation due to hotbands of the complex. At all higher energies, we did not observe any energetically disallowed NO (A) rotational states, and for all available energies above 259 cm-1 the observed rotational transitions do not extend to the maximum allowed by energy conservation. Furthermore, the observed distributions were typically biased toward low rotational states, in contrast with expectations from vibrational predissociation. From the rotational state distributions, we determined the average fraction of energy partitioned into NO (A) rotation, fNOrot,ave, to be 0.088 at the highest available energy, and this fraction increased as the available energy decreased. By combining the average NO (A) rotational energy along with the average center-of-mass translational energy from our previous work, we determined the average rotational energy for the undetected N2 (X) photoproduct. The results showed that the N2 fragment has a higher average rotational energy relative to the NO fragment. Finally, we found that the NO (A) rotational state distribution was colder than expected for a statistical dissociation.

Rotational Cooling Effect on the Rate Constant in the CH3F + Ca+ Reaction at Low Collision Energies.

The rotational cooling effect on the reaction rate constant of the gas-phase ion-polar-molecule reaction CH3F + Ca+ → CH3 + CaF+ was experimentally studied at low collision energies. Fluoromethane molecules showed higher reactivity as the rotational temperature decreased. The experimental rate constants were compared with the capture rate constants which were obtained by the Perturbed Rotational State (PRS) theory assuming the rotational level distribution corresponding to the experimental conditions. The PRS result shows a strong dependence of the capture rate constants on the rotational level distribution in accordance with the experimental findings. However, the PRS capture rate constants deviate from the measurement values as the average collision energy increases especially when the fluoromethane molecules are rotationally cooled far below room temperature. The present paper suggests that the rotational state distribution significantly affects the rate constants of ion-polar-molecule reactions and is one of the important issues to be considered in the study of molecular synthesis in the interstellar medium, where the thermal equilibrium is not necessarily established.

State-to-state integral cross sections and rate constants for the N+(3P)+HD→NH+/ND++D/H reaction: Accurate quantum dynamics studies

The reactive collisions of nitrogen ion with hydrogen and its isotopic variations have great significance in the field of astrophysics. Herein, the state-to-state quantum time-dependent wave packet calculations of N+(3P) + HD → NH+/ND+ + D/H reaction are carried out based on the recently developed potential energy surface [Phys. Chem. Chem. Phys. 21 22203 (2019)]. The integral cross sections (ICSs) and rate coefficients of both channels are precisely determined at the state-to-state level. The results of total ICSs and rate coefficients present a dramatic preference on the ND+ product over the NH+ product, conforming to the long-lived complex-forming mechanism. Product state-resolved ICSs indicate that both the product molecules are difficult to excite to higher vibrational states, and the ND+ product has a hotter rotational state distribution. Moreover, the integral cross sections and rate coefficients are precisely determined at the state-to-state level and insights are provided about the differences between the two channels. The present results would provide an important reference for the further experimental studies at the finer level for this interstellar chemical reaction. The datasets presented in this paper, including the ICSs and rate coefficients of the two products for the title reaction, are openly available at https://www.doi.org/10.57760/sciencedb.j00113.00034.

Phase curves of small bodies from the SLOAN Moving Objects Catalog

Context. Extensive photometric surveys continue to produce enormous stores of data on small bodies. These data are typically sparsely obtained at arbitrary (or unknown) rotational phases. Therefore, new methods for processing such data need to be developed to make the most of these vast catalogs. Aims. We aim to produce a method of recreating the phase curves of small bodies by considering the uncertainties introduced by the nominal errors in the magnitudes and the effect introduced by rotational variations. We use the SLOAN Moving Objects Catalog data as a benchmark to construct phase curves of all small bodies in u′, g′, r′, i′, and z′ filters. From the phase curves, we obtain the absolute magnitudes and we use them to set up the absolute colors, which are the colors of the asteroids that are not affected by changes in the phase angle. Methods. We selected objects with ≥3 observations taken in at least one filter and spanning over a minimum of 5 degrees in the phase angle. We developed a method that combines Monte Carlo simulations and Bayesian inference to estimate the absolute magnitudes using the HG12* photometric system. Results. We obtained almost 15 000 phase curves, with about 12 000 of these including all five filters. The absolute magnitudes and absolute colors are compatible with previously published data that support our method. Conclusions. The method we developed is fully automatic and well suited for a run based on large amounts of data. Moreover, it includes the nominal uncertainties in the magnitudes and the whole distribution of possible rotational states of the objects producing what are possibly less precise values, that is, larger uncertainties, but more accurate, namely, closer to the actual value. To our knowledge, this work is the first to include the effect of rotational variations in such a manner.

Long-range model of vibrational autoionization in core-nonpenetrating Rydberg states of NO.

In high orbital angular momentum (ℓ ≥ 3) Rydberg states, the centrifugal barrier hinders the close approach of the Rydberg electron to the ion-core. As a result, these core-nonpenetrating Rydberg states can be well described by a simplified model in which the Rydberg electron is only weakly perturbed by the long-range electric properties (i.e., multipole moments and polarizabilities) of the ion-core. We have used a long-range model to describe the vibrational autoionization dynamics of high-ℓ Rydberg states of nitric oxide (NO). In particular, our model explains the extensive angular momentum exchange between the ion-core and the Rydberg electron that had been previously observed in vibrational autoionization of f (ℓ = 3) Rydberg states. These results shed light on a long-standing mechanistic question around these previous observations and support a direct, vibrational mechanism of autoionization over an indirect, predissociation-mediated mechanism. In addition, our model correctly predicts newly measured total decay rates of g (ℓ = 4) Rydberg states because for ℓ ≥ 4, the non-radiative decay is dominated by autoionization rather than predissociation. We examine the predicted NO+ ion rotational state distributions generated by vibrational autoionization of g states and discuss applications of our model to achieve quantum state selection in the production of molecular ions.