Theoretical studies of collisional energy transfer in highly excited molecules: The importance of intramolecular vibrational redistribution in successive collision modeling of He+CS2

Comparison of intermolecular energy transfer from vibrationally excited benzene in mixed nitrogen-benzene baths at 140 K and 300 K

Gas phase intermolecular energy transfer (IET) is a fundamental component of accurately explaining the behavior of gas phase systems in which the internal energy of particular modes of molecules is greatly out of equilibrium. In this work, chemical dynamics simulations of mixed benzene/N₂ baths with one highly vibrationally excited benzene molecule (Bz^*) are compared to experimental results at 140 K. Two mixed bath models are considered. In one, the bath consists of 190 N₂ and 10 Bz, whereas in the other bath, 396 N₂ and 4 Bz are utilized. The results are compared to results from 300 K simulations and experiments, revealing that Bz^*-Bz vibration-vibration IET efficiency increased at low temperatures consistent with longer lived "chattering" collisions at lower temperatures. In the simulations, at the Bz^* excitation energy of 150 kcal/mol, the averaged energy transferred per collision, ⟨ΔE_c⟩, for Bz^*-Bz collisions is found to be ∼2.4 times larger in 140 K than in 300 K bath, whereas this value is ∼1.3 times lower for Bz^*-N₂ collisions. The overall ⟨ΔE_c⟩, for all collisions, is found to be almost two times larger at 140 K compared to the one obtained from the 300 K bath. Such an enhancement of IET efficiency at 140 K is qualitatively consistent with the experimental observation. However, the possible reasons for not attaining a quantitative agreement are discussed. These results imply that the bath temperature and molecular composition as well as the magnitude of vibrational energy of a highly vibrationally excited molecule can shift the overall timescale of rethermalization.

Modeling Collisional Transitions in Thermal Unimolecular Reactions: Successive Trajectories and Two-Dimensional Master Equation for Trifluoromethane Decomposition in an Argon Bath

Collisional transition processes in thermal unimolecular reactions are modeled by collision frequency, Z, and probability distribution function, P(E, J; E', J'), which describes the probabilities of collisional transitions from the initial state specified by the total energy and angular momentum, (E', J'), to the final states, (E, J). The validity of the collisional transition model, consisting of Z and P(E, J; E', J'), is assessed here for the title reaction. The present model and its parameters are derived from the moments of transition probabilities calculated by classical trajectory simulations. The model explicitly accounts for coupling between the energy and angular momentum transfer and the dependence of transition probability on the initial state. The performance of the model is evaluated by comparing the rate constants calculated by solving the two-dimensional master equation with those obtained from the classical trajectory calculations of the sequence of successive collisions. The rate constants are also compared with available experimental data. The present collisional transition model is found to perform fairly well for predicting the pressure-dependent rate constants. The uncertainty in the prediction and sensitivities of the rate constants to the model parameters are discussed. A simplified version of the model is proposed, which performs as well as the full model. The simplifications and robust procedures for calculating the model parameters are described.

Pressure effects on the vibrational and rotational relaxation of vibrationally excited OH (ν, J) in an argon bath

Quasi-classical molecular dynamics simulations were used to study the energy relaxation of an initially non-rotating, vibrationally excited (ν = 4) hydroxyl radical (OH) in an Ar bath at 300 K and at high pressures from 50 atm to 400 atm. A Morse oscillator potential represented the OH, and two sets of interaction potentials were used based on whether the Ar-H potential was a Buckingham (Exp6) or a Lennard-Jones (LJ) potential. The vibrational and rotational energies were monitored for 25 000-90 000 ps for Exp6 trajectories and 5000 ps for LJ trajectories. Comparisons to measured vibrational relaxation rates show that Exp6 rates are superior. Simulated initial vibrational relaxation rates are linearly proportional to pressure, implying no effect of high-pressure breakdown in the isolated binary collision approximation. The vibrational decay curves upward from single-exponential decay. A model based on transition rates that exponentially depend on the anharmonic energy gap between vibrational levels fits the vibrational decay well at all pressures, suggesting that anharmonicity is a major cause of the curvature. Due to the competition of vibration-to-rotation energy transfer and bath gas relaxation, the rotational energy overshoots and then relaxes to its thermal value. Approximate models with adjustable rates for this competition successfully reproduced the rotational results. These models show that a large fraction of the vibrational energy loss is initially converted to rotational energy but that fraction decreases rapidly as the vibrational energy content of OH decreases. While simulated rates change dramatically between Exp6 and LJ potentials, the mechanisms remain the same.

Nonequilibrium internal energy distributions during dissociation

In this work, we propose a model for nonequilibrium vibrational and rotational energy distributions in nitrogen using surprisal analysis. The model is constructed by using data from direct molecular simulations (DMSs) of rapidly heated nitrogen gas using an ab initio potential energy surface (PES). The surprisal-based model is able to capture the overpopulation of high internal energy levels during the excitation phase and also the depletion of high internal energy levels during the quasi-steady-state (QSS) dissociation phase. Due to strong coupling between internal energy and dissociation chemistry, such non-Boltzmann effects can influence the overall dissociation rate in the gas. Conditions representative of the flow behind strong shockwaves, relevant to hypersonic flight, are analyzed. The surprisal-based model captures important molecular-level nonequilibrium physics, yet the simple functional form leads to a continuum-level expression that now accounts for the underlying energy distributions and their coupling to dissociation.

Energy transfer of highly vibrationally excited phenanthrene and diphenylacetylene

The energy transfer between Kr atoms and highly vibrationally excited, rotationally cold phenanthrene and diphenylacetylene in the triplet state was investigated using crossed-beam/time-of-flight mass spectrometer/time-sliced velocity map ion imaging techniques. Compared to the energy transfer between naphthalene and Kr, energy transfer between phenanthrene and Kr shows a larger cross-section for vibrational to translational (V → T) energy transfer, a smaller cross-section for translational to vibrational and rotational (T → VR) energy transfer, and more energy transferred from vibration to translation. These differences are further enlarged in the comparison between naphthalene and diphenylacetylene. In addition, less complex formation and significant increases in the large V → T energy transfer probabilities, termed supercollisions in diphenylacetylene and Kr collisions were observed. The differences in the energy transfer between these highly vibrationally excited molecules are attributed to the low-frequency vibrational modes, especially those vibrations with rotation-like wide-angle motions.

Quasi-Classical Trajectory Simulations of Intramolecular Vibrational Energy Redistribution in HONO2 and DONO2

By use of an analytic potential energy surface developed in this work for nitric acid, the quasi-classical trajectory method was used to simulate intramolecular vibrational energy redistribution (IVR). A method was developed for monitoring the average vibrational energy in the OH (or OD) mode that uses the mean-square displacement of the bond length calculated during the trajectories. This method is effective for both rotating and nonrotating molecules. The calculated IVR time constant for HONO(2) decreases exponentially with increasing excitation energy, is almost independent of rotational temperature, and is in excellent agreement with the experimental determination (Bingemann, D.; Gorman, M. P.; King, A. M.; Crim, F. F. J. Chem.Phys. 1997, 107, 661). In DONO(2), the IVR time constants show more complicated behavior with increasing excitation energy, apparently due to 2:1 Fermi-resonance coupling with lower frequency modes. This effect should be measurable in experiments.