CF3CH3 --> HF + CF2CH2: a non-RRKM reaction?

The experimental shock tube data recently reported by Kiefer et al. [J. Phys. Chem. A 2004, 108, 2443-2450] for the title reaction at temperatures between 1600 and 2400 K have been compared to master equation simulations using three models: (a) standard RRKM theory, (b) RRKM theory modified by local random matrix theory, which introduces dynamical corrections arising from slow intramolecular vibrational energy randomization, and (c) an ad hoc empirical non-RRKM model. Only the third model provides a good fit of the Kiefer et al. unimolecular reaction rate data. In separate simulations, all three models accurately reproduce the experimental 300 K chemical activation data of Marcoux and Setser [J. Phys. Chem. 1978, 82, 97-108] when the energy transfer parameters are freely varied to fit the data. When experimental energy transfer parameters for a geometrical isomer (1,1,2-trifluoroethane) are used, the standard RRKM model fits the chemical activation data better than the other models, but if energy transfer in the 1,1,1-trifluoroethane is significantly reduced in comparison to the 1,1,2 isomer, then the empirical ad hoc non-RRKM model also gives a good fit. While the ad hoc empirical non-RRKM model can be made to fit the data, it is not based on theory, and we argue that it is physically unrealistic. We also show that the master equation simulations can mimic the Kiefer et al. vibrational relaxation data, which was the first shock tube observation of double-exponential relaxation. We conclude that, until more data on the trifluoroethanes become available, the current evidence is insufficient to decide with confidence whether non-RRKM effects are important in this reaction, or whether the Kiefer et al. data can be explained in some other way.     

Dissociation of 1,1,1-trifluoroethane behind reflected shock waves: shock tube/time-of-flight mass spectrometry experiments

The dissociation of 1,1,1,-trifluoroethane, a potential non-RRKM reaction, has been studied at 600 and 1200 Torr and high temperatures (1500-1840 K) using a new shock tube/time-of-flight mass spectrometer (ST/TOF-MS). These data obtained by an independent method are in good agreement with the laser schlieren, LS, experiments of Kiefer et al. [J. Phys. Chem. A 2004, 108, 2443-2450] and extend the range of that experimental dataset. The data have been simulated by both standard RRKM calculations and the non-RRKM model reported by Kiefer et al. but with = 750 cm(-1). Both the RRKM and non-RRKM calculations provide equally good fits to the ST/TOF-MS data. Neither model simulates the combined ST/TOF-MS and LS datasets particularly well. However, the non-RRKM model predicts a pressure dependency closer to that observed in the experiments than the RRKM model.     

Overlapping resonances in the control of intramolecular vibrational redistribution

Coherent control of bound state processes via the interfering overlapping resonance scenario [Christopher et al., J. Chem. Phys. 123, 064313 (2006)] is developed to control intramolecular vibrational redistribution (IVR). The approach is applied to the flow of population between bonds in a model of chaotic OCS vibrational dynamics, showing the ability to significantly alter the extent and rate of IVR by varying quantum interference contributions.     

Theoretical investigation of highly excited vibrational states in DFCO: calculation of the out-of-plane bending states and simulation of the intramolecular vibrational energy redistribution

A previously developed modified Davidson scheme [C. Iung and F. Ribeiro, J. Chem. Phys. 121, 174105 (2005)] is applied to compute and analyze highly excited (nu2,nu6) eigenstates in DFCO. The present paper is also devoted to the simulations of the intramolecular vibrational energy redistribution (IVR) initiated by an excitation of the out-of-plane bending vibration (nnu6, n=2,4,6, . . . ,18, and 20). The multiconfiguration time-dependent Hartree method is exploited to propagate the corresponding six-dimensional wave packets. A comprehensive comparison with experimental data as well as with previous simulations of IVR in HFCO [G. Pasin et al. J. Chem. Phys. 124, 194304 (2006)] is presented.     

Quantum state-to-state rate constants for the rotationally inelastic collision of CH(B2Sigma(-), nu=0, N-->N') with Ar

We have calculated the state-to-state integral cross sections and rate constants for the rotationally inelastic collision of CH(B (2)Sigma(-), nu=0, N-->N') with Ar using the quantum coupled-state and close-coupling methods on an ab initio potential-energy surface constructed by Alexander et al. [J. Chem. Phys. 101, 4547 (1994)]. Overall the calculated rate constants are in good agreements with the three available experimental results. The rate constants are comparable to the usual gas kinetic and decrease with increasing N and DeltaN. For the multiquantum transition cases, the theory underestimates the experiment. We discuss some possible causes to the discrepancies among the theory and the experiments.     

A global 12-dimensional ab initio potential energy surface and dynamical studies for the SiH4+H-->SiH3+H2 reaction

A global 12-dimensional ab initio interpolated potential energy surface (PES) for the SiH(4)+H-->SiH(3)+H(2) reaction is presented. The ab initio calculations are based on the unrestricted quadratic configuration interaction treatment with all single and double excitations together with the cc-pVTZ basis set, and the modified Shepard interpolation method of Collins and co-workers [K. C. Thompson et al., J. Chem. Phys. 108, 8302 (1998); M. A. Collins, Theor. Chem. Acc. 108, 313 (2002); R. P. A. Bettens and M. A. Collins, J. Chem. Phys. 111, 816 (1999)] is applied. Using this PES, classical trajectory and variational transition state theory calculations have been carried out, and the computed rate constants are in good agreement with the available experimental data.     

Including anharmonicity in the calculation of rate constants. II. The OH+H2-->H2O+H reaction

A recently developed method for calculating anharmonic vibrational energy levels at nonstationary points along a reaction path that is based on second-order perturbation theory in curvilinear coordinates is combined with variational transition state theory with semiclassical multidimensional tunneling approximations to calculate thermal rate constants for the title reaction. Two different potential energy surfaces were employed for these calculations, an improved version of the author's surface 5 and the WSLFH surface of Wu et al. [J. Chem. Phys. 113, 3150 (2000)]. We present detailed comparisons of rate constants computed for the two surfaces with and without anharmonicity and with various approximations for incorporating tunneling along the reaction path. The results for this system are quite sensitive to the surface employed, the choice of coordinates (curvilinear versus rectilinear), and the inclusion of anharmonicity. A comparison with experiment provides information on the accuracy of these surfaces.     

Application of the Wolf damped Coulomb method to simulations of SiC

A multibody interatomic potential is developed for bulk SiC using a modification of the Wolf et al. summation technique [D. Wolf, P. Keblinski, S. R. Phillpot, and J. Eggebrecht, J. Chem. Phys. 110, 8254 (1999)] for the electrostatic interaction. The technique is modified to account for the short-range nonpoint charge effect. The nonelectrostatic interaction is modeled by a simple Morse-stretch term. This potential is then applied to beta-SiC to calculate various bulk properties using molecular dynamics simulations. The simulated x-ray diffraction pattern, radial distribution functions, lattice constant, elastic constants, and defect energy agree well with experimental data.     

Using the thermal Gaussian approximation for the Boltzmann operator in semiclassical initial value time correlation functions

The thermal Gaussian approximation (TGA) recently developed by Frantsuzov et al. [Chem. Phys. Lett. 381, 117 (2003)] has been demonstrated to be a practical way for approximating the Boltzmann operator exp(-betaH) for multidimensional systems. In this paper the TGA is combined with semiclassical (SC) initial value representations (IVRs) for thermal time correlation functions. Specifically, it is used with the linearized SC-IVR (LSC-IVR, equivalent to the classical Wigner model), and the "forward-backward semiclassical dynamics" approximation developed by Shao and Makri [J. Phys. Chem. A 103, 7753 (1999); 103, 9749 (1999)]. Use of the TGA with both of these approximate SC-IVRs allows the oscillatory part of the IVR to be integrated out explicitly, providing an extremely simple result that is readily applicable to large molecular systems. Calculation of the force-force autocorrelation for a strongly anharmonic oscillator demonstrates its accuracy, and calculation of the velocity autocorrelation function (and thus the diffusion coefficient) of liquid neon demonstrates its applicability.     

Vibrational energy relaxation rates of H2 and D2 in liquid argon via the linearized semiclassical method

The vibrational energy relaxation (VER) rates for H2 and D2 in liquid argon (T=152 K, rho=1.45x1022 cm-3) are calculated using the linearized semiclassical (LSC) method (J. Phys. Chem. 2003, 107, 9059, 9070). The calculation is based on Fermi's golden rule. The VER rate constant is expressed in terms of the quantum-mechanical force-force correlation function, which is then estimated using the LSC method. A local harmonic approximation (LHA) is employed in order to compute the multidimensional Wigner integrals underlying the LSC approximation. The H2-Ar and D2-Ar interactions are described by the three-body potential of Bissonette et al. (J. Phys. Chem. A 1996, 105, 2639). The LHA-LSC-based VER rate constants for both D2 and H2 are found to be about 2-3 orders of magnitude slower than those obtained experimentally. However, their ratio agrees quantitatively with the corresponding experimental result. In contrast, the classical VER rate constants are found to be 8-9 orders of magnitude slower than those obtained experimentally, and their ratio is found to be qualitatively different from the corresponding experimental result.     

Comparison of different methods for calculating the paramagnetic relaxation enhancement of nuclear spins as a function of the magnetic field

The enhancement of the spin-lattice relaxation rate for nuclear spins in a ligand bound to a paramagnetic metal ion [known as the paramagnetic relaxation enhancement (PRE)] arises primarily through the dipole-dipole (DD) interaction between the nuclear spins and the electron spins. In solution, the DD interaction is modulated mostly by reorientation of the nuclear spin-electron spin axis and by electron spin relaxation. Calculations of the PRE are in general complicated, mainly because the electron spin interacts so strongly with the other degrees of freedom that its relaxation cannot be described by second-order perturbation theory or the Redfield theory. Three approaches to resolve this problem exist in the literature: The so-called slow-motion theory, originating from Swedish groups [Benetis et al., Mol. Phys. 48, 329 (1983); Kowalewski et al., Adv. Inorg. Chem. 57, (2005); Larsson et al., J. Chem. Phys. 101, 1116 (1994); T. Nilsson et al., J. Magn. Reson. 154, 269 (2002)] and two different methods based on simulations of the dynamics of electron spin in time domain, developed in Grenoble [Fries and Belorizky, J. Chem. Phys. 126, 204503 (2007); Rast et al., ibid. 115, 7554 (2001)] and Ann Arbor [Abernathy and Sharp, J. Chem. Phys. 106, 9032 (1997); Schaefle and Sharp, ibid. 121, 5387 (2004); Schaefle and Sharp, J. Magn. Reson. 176, 160 (2005)], respectively. In this paper, we report a numerical comparison of the three methods for a large variety of parameter sets, meant to correspond to large and small complexes of gadolinium(III) and of nickel(II). It is found that the agreement between the Swedish and the Grenoble approaches is very good for practically all parameter sets, while the predictions of the Ann Arbor model are similar in a number of the calculations but deviate significantly in others, reflecting in part differences in the treatment of electron spin relaxation. The origins of the discrepancies are discussed briefly.     

Semiclassical calculation of nonadiabatic thermal rate constants: application to condensed phase reactions

The nonadiabatic transition state theory proposed recently by Zhao et al. [J. Chem. Phys. 121, 8854 (2004)] is extended to calculate rate constants of complex systems by using the Monte Carlo and umbrella sampling methods. Surface hopping molecular dynamics technique is incorporated to take into account the dynamic recrossing effect. A nontrivial benchmark model of the nonadiabatic reaction in the condensed phase is used for the numerical test. It is found that our semiclassical results agree well with those produced by the rigorous quantum mechanical method. Comparing with available analytical approaches, we find that the simple statistical theory proposed by Straub and Berne [J. Chem. Phys. 87, 6111 (1987)] is applicable for a wide friction region although their formula is obtained using Landau-Zener [Phys. Z. Sowjetunion 2, 46 (1932); Proc. R. Soc. London, Ser. A 137, 696 (1932)] nonadiabatic transition probability along a one-dimensional diffusive coordinate. We also investigate how the nuclear tunneling events affect the dependence of the rate constant on the friction.     

Local response dispersion method. II. Generalized multicenter interactions

Recently introduced local response dispersion method [T. Sato and H. Nakai, J. Chem. Phys. 131, 224104 (2009)], which is a first-principles alternative to empirical dispersion corrections in density functional theory, is implemented with generalized multicenter interactions involving both atomic and atomic pair polarizabilities. The generalization improves the asymptote of intermolecular interactions, reducing the mean absolute percentage error from about 30% to 6% in the molecular C(6) coefficients of more than 1000 dimers, compared to experimental values. The method is also applied to calculations of potential energy curves of molecules in the S22 database [P. Jure?ka et al., Phys. Chem. Chem. Phys. 8, 1985 (2006)]. The calculated potential energy curves are in a good agreement with reliable benchmarks recently published by Molnar et al. [J. Chem. Phys. 131, 065102 (2009)]. These improvements are achieved at the price of increasing complexity in the implementation, but without losing the computational efficiency of the previous two-center (atom-atom) formulation. A set of different truncations of two-center and three- or four-center interactions is shown to be optimal in the cost-performance balance.     

Accurate quantum wave packet calculations for the F + HCl → Cl + HF reaction on the ground 1(2)A' potential energy surface

We present converged exact quantum wave packet calculations of reaction probabilities, integral cross sections, and thermal rate coefficients for the title reaction. Calculations have been carried out on the ground 1(2)A' global adiabatic potential energy surface of Deskevich et al. [J. Chem. Phys. 124, 224303 (2006)]. Converged wave packet reaction probabilities at selected values of the total angular momentum up to a partial wave of J = 140 with the HCl reagent initially selected in the v = 0, j = 0-16 rovibrational states have been obtained for the collision energy range from threshold up to 0.8 eV. The present calculations confirm an important enhancement of reactivity with rotational excitation of the HCl molecule. First, accurate integral cross sections and rate constants have been calculated and compared with the available experimental data.     

Quantifying the non-RRKM effect in the H + O2 ⇄ OH + O reaction

In the present investigation the non-RRKM behavior in the title reaction is quantified in two different ways: (1) Quasiclassical trajectory calculations of the thermal rate coefficient are compared with results from a microcanonical variational transition-state theory/RRKM model. Results on both the Varandas DMBE IV and Melius-Blint potentials indicate that the non-RRKM behavior acts to reduce the thermal rate coefficient by about a factor of two, independent of temperature from 250 K to 5500 K. The QCT thermal rate coefficients on the two potentials are in remarkably good agreement with experiment and with each other over the entire temperature range. (2) The non-RRKM behavior as a classical phenomenon is demonstrated and quantified on both potentials by a direct test of the fundamental assumption. Complex-forming classical trajectories, started as either O + OH or H + O₂, are shown preferentially to return to the region of configuration space from which they were started. This test is discussed in detail in the text. The transition of the non-RRKM behavior from classical to quantum mechanics is also discussed. © 1997 John Wiley & Sons, Inc. Int J Chem Kinet 29: 275–287, 1997.     

Communication: Thermodynamic scaling of the Debye process in primary alcohols

The molecular dynamics of hydrogen-bonded liquids usually does not satisfy the thermodynamic scaling. However, very recently, two opposite conclusions about validity of thermodynamical scaling in monohydroxy alcohol, 2-ethyl-1-hexanol, were presented by Reiser et al. [J. Chem. Phys. 132, 181101 (2010)] and Fragiadakis et al. [J. Chem. Phys. 132, 144505 (2010)]. In this communication we present new experimental results that can explain this ostensible contradiction.     

Quantum wave packet dynamics of N(2D)+H2 reaction

The quantum wave packet dynamics of the title reaction within the coupled state approximation is examined here and initial state-selected reaction probabilities, integral reaction cross sections, and thermal rate constants are reported. The ab initio potential energy surface of the electronic ground state (1(2)A(")) of the system recently reported by Ho et al. [J. Chem. Phys., 119, 3063 (2003)] is employed in this investigation. All partial wave contributions up to the total angular momentum J=55 were necessary to obtain converged integral reaction cross sections up to a collision energy of 1.0 eV. Thermal rate constants are calculated from the reaction cross sections and compared with the available theoretical and experimental results. Typical resonances formed during the course of the reaction and elucidating the insertion type mechanism for the product formation are calculated. Vibrational energy levels supported by the deep well (approximately 5.5 eV) of the 1(2)A(") potential energy surface of NH(2) are also calculated for the total angular momentum J=0. A statistical analysis of the spacing between the adjacent levels of this energy spectrum is performed and the extent of irregularity in the spectral sequence is assessed.     

Influence of rotation and isotope effects on the dynamics of the N((2)D)+H(2) reactive system and of its deuterated variants

Integral cross sections and thermal rate constants have been calculated for the N((2)D)+H(2) reaction and its isotopic variants N((2)D)+D(2) and the two-channel N((2)D)+HD by means of quasiclassical trajectory and statistical quantum-mechanical model methods on the latest ab initio potential-energy surface [T.-S. Ho et al., J. Chem. Phys. 119, 3063 (2003)]. The effect of rotational excitation of the diatom on the dynamics of these reactions has been investigated and interesting discrepancies between the classical and statistical model calculations have been found. Whereas a net effect of reagent rotation on reactivity is always observed in the classical calculations, only a very slight effect is observed in the case of the asymmetric N((2)D)+HD reaction for the statistical quantum-mechanical method. The thermal rate constants calculated on this Potential-Energy Surface using quasiclassical trajectory and statistical model methods are in good agreement with the experimental determinations, although the latter are somewhat larger. A reevaluation of the collinear barrier of the potential surface used in the present study seems timely. Further theoretical and experimental studies are needed for a full understanding of the dynamics of the title reaction.     

A quasiclassical trajectory study for the N(4S)+O2(X 3Sigmag-)-->NO(X 2Pi)+O(3P) reaction on the new 2A' and 4A' potential-energy surfaces

A quasiclassical trajectory study with the sixth-order explicit symplectic algorithm of the N(4S)+O2(X 3Sigmag-)-->NO(X 2Pi)+O(3P) atmospheric reaction has been performed by employing the new 2A' and 4A' potential-energy surfaces reported by Sayos et al. [J. Chem. Phys. 117, 670 (2002)]. For the translational temperature considered up to 10,000 K, the larger relative translational energy and the higher rovibrational levels of O2 molecule with respect to the previous works have been taken into account, and a clearer database about the character of the total reaction cross section has been presented in this work. The dependence of microscopic rate constants on the vibrational level of O2 molecule at T=3000, 5000, and 10,000 K has been exhibited, and we can see that the values of log10 k(T,v,J) vary almost linearly with the vibrational level of O2 molecule. The thermal rate constants at the translational temperature between 300 and 10,000 K have been evaluated and compared with the experimental and previous theoretical results. It is found that the thermal rate constants determined in this work have a better agreement with the experimental data and can provide a more valid theoretical reference at the translational temperature considered for the title reaction.     

Path integral calculation of thermal rate constants within the quantum instanton approximation: application to the H + CH4 --> H2 + CH3 hydrogen abstraction reaction in full Cartesian space

The quantum instanton approximation for thermal rate constants of chemical reactions [Miller, Zhao, Ceotto, and Yang, J. Chem. Phys. 119, 1329 (2003)], which is modeled after the earlier semiclassical instanton approach, is applied to the hydrogen abstraction reaction from methane by a hydrogen atom, H + CH4 --> H2 + CH3, using a modified and recalibrated version of the Jordan-Gilbert potential surface. The quantum instanton rate is evaluated using path integral Monte Carlo approaches based on the recently proposed implementation schemes [Yamamoto and Miller, J. Chem. Phys. 120, 3086 (2004)]. The calculations were carried out using the Cartesian coordinates of all the atoms (thus involving 18 degrees of freedom), thereby taking explicit account of rotational effects of the whole system and also allowing the equivalent treatment of the four methane hydrogens. To achieve such a treatment, we present extended forms of the path integral estimators for relevant quantities that may be used for general N-atom systems with any generalized reaction coordinates. The quantum instanton rates thus obtained for the temperature range T = 200-2000 K show good agreement with available experimental data, which gives support to the accuracy of the underlying potential surface used.