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.

This study examined if individuals who are at increased risk for drowsy-driving because of obstructive sleep apnea syndrome (OSAS), have impairments in driving performance in the moments during microsleep episodes as opposed to during periods of wakefulness. Twenty-four licensed drivers diagnosed with OSAS based on standard clinical and polysomnographic criteria, participated in an hour-long drive in a high-fidelity driving simulator with synchronous electroencephalographic (EEG) recordings for identification of microsleeps. The drivers showed significant deterioration in vehicle control during the microsleep episodes compared to driving performance in the absence of microsleeps on equivalent segments of roadway. The degree of performance decrement correlated with microsleep duration, particularly on curved roads. Results indicate that driving performance deteriorates during microsleep episodes. Detecting microsleeps in real-time and identifying how these episodes of transition between wakefulness and sleep impair driver performance is relevant to the design and implementation of countermeasures such as drowsy driver detection and alerting systems that use EEG technology.

A zero-point energy (ZPE) constraint model is proposed for classical trajectory simulations of unimolecular decomposition and applied to CH4* → H + CH3 decomposition. With this model trajectories are not allowed to dissociate unless they have ZPE in the CH3 product. If not, they are returned to the CH4* region of phase space and, if necessary, given additional opportunities to dissociate with ZPE. The lifetime for dissociation of an individual trajectory is the time it takes to dissociate with ZPE in CH3, including multiple possible returns to CH4*. With this ZPE constraint the dissociation of CH4* is exponential in time as expected for intrinsic RRKM dynamics and the resulting rate constant is in good agreement with the harmonic quantum value of RRKM theory. In contrast, a model that discards trajectories without ZPE in the reaction products gives a CH4* → H + CH3 rate constant that agrees with the classical and not quantum RRKM value. The rate constant for the purely classical simulation indicates that anharmonicity may be important and the rate constant from the ZPE constrained classical trajectory simulation may not represent the complete anharmonicity of the RRKM quantum dynamics. The ZPE constraint model proposed here is compared with previous models for restricting ZPE flow in intramolecular dynamics, and connecting product and reactant/product quantum energy levels in chemical dynamics simulations.

We present the molecular symmetry (MS) adapted treatment of nonadiabatic coupling terms (NACTs) for the excited electronic states (2(2)E' and 1(2)A1') of Na3 cluster, where the adiabatic potential energy surfaces (PESs) and the NACTs are calculated at the MRCI level by using an ab initio quantum chemistry package (MOLPRO). The signs of the NACTs at each point of the configuration space (CS) are determined by employing appropriate irreducible representations (IREPs) arising due to MS group, and such terms are incorporated into the adiabatic to diabatic transformation (ADT) equations to obtain the ADT angles. Since those sign corrected NACTs and the corresponding ADT angles demonstrate the validity of curl condition for the existence of three-state (2(2)E' and 1(2)A1') sub-Hilbert space, it becomes possible to construct the continuous, single-valued, symmetric, and smooth 3 × 3 diabatic Hamiltonian matrix. Finally, nuclear dynamics has been carried out on such diabatic surfaces to explore whether our MS-based treatment of diabatization can reproduce the pattern of the experimental spectrum for system B of Na3 cluster.