Molecular dynamics simulations were used to study relaxation of a vibrationally excited C6F6* molecule in a N2 bath. Ab initio calculations were performed to develop N2-N2 and N2-C6F6 intermolecular potentials for the simulations. Energy transfer from "hot" C6F6 is studied versus the bath density (pressure) and number of bath molecules. For the large bath limit, there is no heating of the bath. As C6F6* is relaxed, the average energy of C6F6* is determined versus time, i.e., ⟨E(t)⟩, and for each bath density ⟨E(t)⟩ is energy dependent and cannot be fit by a single exponential. In the long-time limit C6F6 is fully equilibrated with the bath. For a large bath and low pressures, the simulations are in the fixed temperature, independent collision regime and the simulation results may be compared with gas phase experiments of collisional energy transfer. The derivative d[⟨E(t)⟩]/dt divided by the collision frequency ω of the N2 bath gives the average energy transferred from C6F6* per collision ⟨ΔE(c)⟩, which is in excellent agreement with experiment. For the ~100-300 ps simulations reported here, energy transfer from C6F6* is to N2 rotation and translation in accord with the equipartition model, with no energy transfer to N2 vibration. The energy transfer dynamics from C6F6* is not statistically sensitive to fine details of the N2-C6F6 intermolecular potential. Tests, with simulation ensembles of different sizes, show that a relatively modest ensemble of only 24 trajectories gives statistically meaningful results.

Relaxation of vibrationally excited C 6 F 6 * in a thermalized bath of N 2 molecules is studied by condensed-phase chemical dynamics simulations. The average energy of C 6 F 6 as a function of time, ⟨E(t)⟩, was determined using two different models for the N 2 −C 6 F 6 intermolecular potential, and both gave statistically the same result. A simulation with a N 2 bath density of 20 kg/m 3 was performed to check the convergence and validate the results obtained previously with a higher bath density of 40 kg/m 3 . The initial ensemble of C 6 F 6 is nearly monoenergetically excited, but the ensemble acquires as energy distribution P(E) as it is collisionally relaxed. An evaluation of P(E) and the root-mean-square deviation ⟨ΔE 2 ⟩ 1/2 of P(E), versus time, shows that P(E) first broadens and then narrows. Simulations with the C 6 F 6 vibrational excitation energy of 85.8 kcal/mol, studied experimentally, show that the width of P(E) does not affect the average collisional deactivation rate. The role of the intramolecular vibrational frequencies on the energy transfer dynamics was studied by changing the mass of the F-atoms of C 6 F 6 to that of an H-atom. The resulting increase in the frequencies decreased the efficiency of collisional energy transfer. Simulations with classical and quantum microcanoical ensembles of the C 6 F 6 initial vibrational energy were compared, and the initial relaxation rate was slower for the quantum ensemble.

Chemical dynamics simulations are performed to study the unimolecular dissociation of the benzene (Bz)−hexafluorobenzene (HFB) complex at five different temperatures ranging from 1000 to 2000 K, and the results are compared with that of the Bz dimer at common simulation temperatures. Bz−HFB, in comparison with Bz dimer, possesses a much attractive intermolecular interaction, a very different equilibrium geometry, and a lower average quantum vibrational excitation energy at a given temperature. Six low-frequency modes of Bz−HFB are formed by Bz + HFB association which are weakly coupled with the vibrational modes of Bz and HFB. However, this coupling is found much stronger in Bz−HFB compared to the same in the Bz dimer. The simulations are done with very good potential energy parameters taken from the literature. Considering the canonical (TST) model, the unimolecular dissociation rate constant at each temperature is calculated and fitted to the Arrhenius equation. An activation energy of 5.0 kcal/mol and a pre-exponential factor of 2.39 × 10 12 s −1 are obtained, which are of expected magnitudes. The responsible vibrational mode for dissociation is identified by performing normal-mode analysis. Simulations with random excitations of high-frequency Bz and HFB modes and lowfrequency inter-Bz−HFB vibrational modes of the Bz−HFB complex are also performed. The intramolecular vibrational energy redistribution (IVR) time and the unimolecular dissociation rate constants are calculated from these simulations. The latter shows good agreement with the same obtained from simulation with random excitation of all vibrational modes.

Classical chemical dynamics simulations were performed to study the intramolecular and unimolecular dissociation dynamics of the benzene dimer, Bz2 → 2 Bz. The dissociation of microcanonical ensembles of Bz2 vibrational states, at energies E corresponding to temperatures T of 700-1500 K, were simulated. For the large Bz2 energies and large number of Bz2 vibrational degrees of freedom, s, the classical microcanonical (RRKM) and canonical (TST) rate constant expressions become identical. The dissociation rate constant for each T is determined from the initial rate dN(t)/dt of Bz2 dissociation, and the k(T) are well-represented by the Arrhenius eq k(T) = A exp(-E(a)/RT). The E(a) of 2.02 kcal/mol agrees well with the Bz2 dissociation energy of 2.32 kcal/mol, and the A-factor of 2.43 × 10(12) s(-1) is of the expected order-of-magnitude. The form of N(t) is nonexponential, resulting from weak coupling between the Bz2 intramolecular and intermolecular modes. With this weak coupling, large Bz2 vibrational excitation, and low Bz2 dissociation energy, most of the trajectories dissociate directly. Simulations, with only the Bz2 intramolecular modes excited at 1000 K, were also performed to study intramolecular vibrational energy redistribution (IVR) between the intramolecular and intermolecular modes. Because of restricted IVR, the initial dissociation is quite slow, but N(t) ultimately becomes exponential, suggesting an IVR time of 20.7 ps.

Chemical dynamics simulations were performed to investigate collisional energy transfer from highly vibrationally excited azulene (Az*) in a N2 bath. The intermolecular potential between Az and N2, used for the simulations, was determined from MP2/6-31+G* ab initio calculations. Az* is prepared with an 87.5 kcal/mol excitation energy by using quantum microcanonical sampling, including its 95.7 kcal/mol zero-point energy. The average energy of Az* versus time, obtained from the simulations, shows different rates of Az* deactivation depending on the N2 bath density. Using the N2 bath density and Lennard-Jones collision number, the average energy transfer per collision ⟨ΔEc⟩ was obtained for Az* as it is collisionally relaxed. By comparing ⟨ΔEc⟩ versus the bath density, the single collision limiting density was found for energy transfer. The resulting ⟨ΔEc⟩, for an 87.5 kcal/mol excitation energy, is 0.30 ± 0.01 and 0.32 ± 0.01 kcal/mol for harmonic and anharmonic Az potentials, respectively. For comparison, the experimental value is 0.57 ± 0.11 kcal/mol. During Az* relaxation there is no appreciable energy transfer to Az translation and rotation, and the energy transfer is to the N2 bath.