Vibration-rotation energy pattern in acetylene: 13CH12CH up to 6750 cm-1

All known vibration-rotation absorption lines of 13CH12CH accessing levels up to 6750 cm-1 were gathered from the literature. They were fitted simultaneously to J-dependent Hamiltonian matrices exploiting the well known vibrational polyad or cluster block diagonalization, in terms of the pseudo-quantum-numbers Ns=v1+v2+v3 and Nr=5v1+3v2+5v3+v4+v5, and accounting also for l parity and ef symmetry properties. The anharmonic interaction coupling terms known to occur from a pure vibrational fit in this acetylene isotopologue [Robert et al., J. Chem. Phys. 123, 174302 (2005)] were included in the model. A total of 12 703 transitions accessing 158 different (v1v2v3v4v5,l4l5) vibrational states was fitted with a dimensionless standard deviation of 0.99, leading to the determination of 216 vibration-rotation parameters. The experimental data included very weak vibration-rotation transitions accessing 18 previously unreported states, some of them forming Q branches with very irregular patterns.     

Quasi-classical trajectory study of the dynamics of the Cl + CH4→ HCl + CH3 reaction

We present an on-the-fly classical trajectory study of the Cl + CH(4)→ HCl + CH(3) reaction using a specific reaction parameter (SRP) AM1 Hamiltonian that was previously optimized for the Cl + ethane reaction [S. J. Greaves et al., J. Phys Chem A, 2008, 112, 9387]. The SRP-AM1 Hamiltonian is shown to be a good model for the potential energy surface of the title reaction. Calculated differential cross sections, obtained from trajectories propagated with the SRP-AM1 Hamiltonian compare favourably with experimental results for this system. Analysis of the vibrational modes of the methyl radical shows different scattering distributions for ground and vibrationally excited products.     

Simple minimum principle to derive a quantum-mechanical/molecular-mechanical method

We propose a minimum principle to derive a QM/MM (quantum-mechanical/molecular-mechanical) method from the first principle. We approximate the Hamiltonian of a spectator substituent as the structure-dependent effective Hamiltonian in a least-squares sense. This effective Hamiltonian is expanded with the orthogonal operator set called the normal-ordered product. We determine the structure-dependent energy that corresponds to the classical MM energy and the extra one-electron potential that takes account of the interface effects. This QM/MM method is free from the double-counting problem and the artificial truncation of the localized molecular orbitals. As a numerical example we determine the one-electron effective Hamiltonian of the methyl group. This effective Hamiltonian is applied to the ethane and CH(3)CH(2)X molecules (X=CH(3), NH(2), OH, F, COOH, NH(3) (+), OH(2) (+), and COO(-)). It reproduced the relative energies, potential energy curves, and the Mulliken populations of the all-electron calculations fairly well.     

Jahn-Teller effect in tetrahedral symmetry: large-amplitude tunneling motion and rovibronic structure of CH4+ and CD4+

The energy level structures of the ground vibronic states of 12CH4+, 13CH4+, and 12CD4+ have been measured by pulsed-field-ionization zero-kinetic-energy photoelectron spectroscopy. The nuclear spin symmetries of the tunneling-rotational levels have been determined in double-resonance experiments via selected rotational levels of the v3=1 and v3=2 vibrational levels of the X 1A1 ground state of CH4. The energy level structures of 12CH4+, 13CH4+, and 12CD4+ have been analyzed with an effective tunneling-rotational Hamiltonian. The analysis together with a group theoretical treatment of the Tx(e+t2) Jahn-Teller effect in the Td(M) group prove that the equilibrium geometry of 12CH4+, 13CH4+, and 12CD4+ has C2v symmetry and characterize the pseudorotational dynamics in these fluxional cations. The tunneling behavior is discussed in terms of the relevant properties of the potential energy surface, some of which have been recalculated at the CCSD(T)/cc-pVTZ level of ab initio theory.     

Quasiclassical trajectory study of the O(3P) + CH4 --> OH + CH3 reaction with a specific reaction parameters semiempirical Hamiltonian

We present a theoretical study of the O(3P) + CH4 --> OH + CH3 reaction using electronic structure, kinetics, and dynamics calculations. We calculate a grid of ab initio points at the PMP2/AUG-cc-pVDZ level to characterize the potential energy surface in regions of up to 1.3 eV above reagents. This grid of ab initio points is used to derive a set of specific reaction parameters (SRP) for the MSINDO semiempirical Hamiltonian. The resulting SRP-MSINDO Hamiltonian improves the quality of the standard Hamiltonian, particularly in regions of the potential energy surface beyond the minimum-energy reaction path. Quasiclassical-trajectory calculations are used to study the reaction dynamics with the original and the improved MSINDO semiempirical Hamiltonians, and a prior surface. The SRP-MSINDO semiempirical Hamiltonian yields OH rotational distributions in agreement with experimental results, improving over the results of the other surfaces. Thermal rate constants estimated with Variational Transition State Theory using the SRP-MSINDO Hamiltonian are also in agreement with experiments. Our results indicate that reparametrized semiempirical Hamiltonians are a good alternative to generating potential energy surfaces for accurate dynamics studies of polyatomic reactions.     

Fermi resonance in CHX3: a hamiltonian in symmetrized curvilinear internal coordinates

The Fermi resonance between the CH stretching and CH bending vibrations in CHX₃ molecules is investigated using a curvilinear internal coordinate Hamiltonian.The eigenvalues are calculated variationally using a separable basis set consisting of functions which are products of a Morse oscillator eigenfunction for the stretch and a two-dimensional harmonic-oscillator wavefunction for the bend. The potential energy parameters are optimized by the least-squares method. The model is applied to the observed infrared spectrum of CHF₃. The energy levels can be reproduced within the uncertainty of the experimental results.     

Vibrational energy levels with arbitrary potentials using the Eckart-Watson Hamiltonians and the discrete variable representation

An effective and general algorithm is suggested for variational vibrational calculations of N-atomic molecules using orthogonal, rectilinear internal coordinates. The protocol has three essential parts. First, it advocates the use of the Eckart-Watson Hamiltonians of nonlinear or linear reference configuration. Second, with the help of an exact expression of curvilinear internal coordinates (e.g., valence coordinates) in terms of orthogonal, rectilinear internal coordinates (e.g., normal coordinates), any high-accuracy potential or force field expressed in curvilinear internal coordinates can be used in the calculations. Third, the matrix representation of the appropriate Eckart-Watson Hamiltonian is constructed in a discrete variable representation, in which the matrix of the potential energy operator is always diagonal, whatever complicated form the potential function assumes, and the matrix of the kinetic energy operator is a sparse matrix of special structure. Details of the suggested algorithm as well as results obtained for linear and nonlinear test cases including H(2)O, H(3) (+), CO(2), HCNHNC, and CH(4) are presented.     

Anharmonic force field and vibrational dynamics of CH2F2 up to 5000 cm(-1) studied by Fourier transform infrared spectroscopy and state-of-the-art ab initio calculations

Difluoromethane (CH(2)F(2), HFC-32) is a molecule used in refrigerant mixtures as a replacement of the more environmentally hazardous, ozone depleting, chlorofluorocarbons. On the other hand, presenting strong vibration-rotation bands in the 9 μm atmospheric window, it is a greenhouse gas which contributes to global warming. In the present work, the vibrational and ro-vibrational properties of CH(2)F(2), providing basic data for its atmospheric modeling, are studied in detail by coupling medium resolution Fourier transform infrared spectroscopy to high-level electronic structure ab initio calculations. Experimentally a full quantum assignment and accurate integrated absorption cross sections are obtained up to 5000 cm(-1). Ab initio calculations are carried out by using CCSD(T) theory and large basis sets of either the correlation consistent or atomic natural orbital hierarchies. By using vibrational perturbation theory to second order a complete set of vibrational and ro-vibrational parameters is derived from the ab initio quartic anharmonic force fields, which well compares with the spectroscopic constants retrieved experimentally. An excellent agreement between theory and experiment is achieved for vibrational energy levels and integrated absorption cross sections: transition frequencies up to four quanta of vibrational excitation are reproduced with a root mean square deviation (RMSD) of 7 cm(-1) while intensities are predicted within few km mol(-1) from the experiment. Basis set performances and core correlation effects are discussed throughout the paper. Particular attention is focused in the understanding of the anharmonic couplings which rule the vibrational dynamics of the |ν(1)>, |2ν(8)>, |2ν(2)> three levels interacting system. The reliability of the potential energy and dipole moment surfaces in reproducing the vibrational eigenvalues and intensities as well as in modeling the vibrational and ro-vibrational mixings over the whole 400-5000 cm(-1) region is also demonstrated by spectacular spectral simulations carried out by using the ro-vibrational Hamiltonian constants, and the relevant coupling terms, obtained from the perturbation treatment of the ab initio anharmonic force field. The present results suggest CH(2)F(2) as a prototype molecule to test ab initio calculations and theoretical models.     

Trajectory dynamics study of collision-induced dissociation of the Ar + CH4 reaction at hyperthermal conditions: vibrational excitation and isotope substitution

We investigate the role of vibrational energy excitation of methane and two deuterated species (CD(4) and CH(2)D(2)) in the collision-induced dissociation (CID) process with argon at hyperthermal energies. The quasi-classical trajectory method has been applied, and the reactive Ar + CH(4) system has been modeled by using a modified version of the CH(4) potential energy surface of Duchovic et al. (J. Phys. Chem. 1984, 88, 1339) and the Ar-CH(4) intermolecular potential function obtained by Troya (J. Phys. Chem. A 2005, 109, 5814). This study clearly shows that CID is markedly enhanced with vibrational excitation and, to a lesser degree, with collision energy. In general, CID increases by exciting stretch vibrational modes of the reactant molecule. For the direct dissociation of CH(4), however, the CID cross sections appear to be essentially independent of which vibrational mode is initially excited. In all situations studied, the CID cross sections are always greater for the Ar + CD(4) reaction than for the Ar + CH(4) one, the Ar + CH(2)D(2) being an intermediate situation. A detailed analysis of the energy transfer processes, including their relation with CID, is also presented.     

Correlation problems in atomic and molecular systems. VII. Application of the open-shell coupled-cluster approach to simple π-electron model systems

The coupled-cluster variational-like direct approach to the calculation of ionization and electron attachment energies and of excitation energies is applied to several π-electron model systems, using the PPP Hamiltonian with various parametrizations. A simple approximation, which represents the triexcited clusters in terms of disconnected W₁T₂ terms, is employed. All the necessary diagrams for both excitation energy and ionization potential (electron affinity) calculations are given in the compact Hugenholtz nonoriented form. The results of the calculations for benzene, trans-butadiene, all-trans-hexatriene, and fulvene are compared with the corresponding full CI results, and the conclusions about the validity and efficiency of this approach are drawn.     

Calculation of the tau dependence of the vibration-internal rotation-overall rotation interactions in CH3OH from molecular structure and molecular dynamics

The Guan and Quade theory for vibration-large-amplitude internal-motion-rotation interactions has been applied to the internal rotation problem in CH(3)OH. Through the molecular dynamics, the cos 3tau and sin 3tau dependence of the torsional-rotational coefficients in the effective Hamiltonian have been calculated from molecular structure. The internal rotation coordinate tau(') for the vibrationally distorted molecule is shown to have the necessary threefold symmetry for all values of tau('). For the methyl deformation modes, the vibrational dependence of the internal rotation potential energy is shown to have a threefold symmetry. The S(t) and S(t)S(t) dependence of the inertia tensor and Coriolis coupling coefficients has been developed in terms of curvilinear internal coordinates. The T transformation separating rotation from vibrations in zeroth order is then applied, the kinetic-energy tensor inverted to momentum space, and finally the effective torsion-rotation coefficients are calculated by Van Vleck perturbation theory. When compared to the empirical results, the kinetic-energy contributions to the cos 3tau and sin 3tau dependence of the coefficients are as follows: 54% of P(a)(2) is accounted for, 28% of P(a)P(b), 16% of P(a)P(c), and 91% of the asymmetry. The calculation is inadequate to account for the P(b)(2),P(c)(2), and P(b)P(c) coefficients, ranging from factors of 20-70, even with the incorrect sign for some of the terms. Anharmonic force contributions from the vibrations have not been used in the calculation since these forces are not known at this time.     

Rotational spectra of rare isotopic species of bromofluoromethane: determination of the equilibrium structure from ab initio calculations and microwave spectroscopy

Guided by theoretical predictions, the rotational spectra of the mono- and bideuterated species of bromofluoromethane, CDH(79)BrF, CDH(81)BrF, CD(2) (79)BrF, and CD(2) (81)BrF, have been recorded for the first time. Assignment of a few hundred rotational transitions led to the accurate determination of the ground-state rotational constants, all of the quartic and most of the sextic centrifugal distortion constants, as well as the full bromine quadrupole-coupling tensor for both (79)Br and (81)Br, in good agreement with corresponding theoretical predictions based on high-level coupled-cluster calculations. The rotational spectra of the (13)C containing species (13)CH(2) (79)BrF and (13)CH(2) (81)BrF have been observed in natural abundance and have been assigned, thus allowing the determination of the rotational and centrifugal distortion constants as well as the bromine quadrupole-coupling tensor. Furthermore, empirical equilibrium structures have been obtained within a least-squares fit procedure using the available experimental ground-state rotational constants for various isotopic species. Vibrational effects have been accounted for in the analysis using vibration-rotation interaction constants derived from anharmonic force fields computed at the second-order Moller-Plesset perturbation theory as well as coupled-cluster (CC) levels. The empirical equilibrium geometries obtained in this way agree well with the corresponding theoretical predictions obtained from CC calculations [at the CCSD(T) level] after extrapolation to the complete basis set limit and inclusion of core-valence correlation corrections and relativistic effects.     

Torsion-vibration coupling in methanol: diabatic behavior in the CH overtone region

Through a fit to methanol CH overtone data, a previously developed 4-dimensional torsion-vibration Hamiltonian is extended to high CH stretch excitation as well as to high torsional excitation. The strength of the torsion-vibration coupling is found to increase with CH stretch excitation. Systematic patterns of near degeneracy (3-, 4-, and 6-fold) are found in different regions of quantum number space. In the region of the CH fundamentals, an approximate a diabatic separation of the torsion (slow degree of freedom) from the CH stretches (fast degrees of freedom) accounts for the pattern of the energy levels and for the signs of the torsional tunneling splittings. For the higher CH overtones (v(CH) > or = 4), a diabatic representation accounts for the torsional structure obtained from the fully coupled calculation and for certain trends found in the pattern of the energy levels.     

CH3＋HNCO反应机理的理论研究

在6-311＋＋G**基组水平上,采用UMP2方法对自由基CH3与HNCO反应机理进行了研究，全参数优化了反应通道上各驻点的几何构型.结果表明, 自由基CH3与HNCO分子间反应有三条反应通道,第一为CH3与HNCO分子间经过生成一个稳定化能为4.56 kJ•mol－1的含氢键的分子复合物M后，经过渡态TS生成另一个产物复合物M′，然后分解为甲烷和NCO自由基；第二是CH3与HNCO分子间通过生成稳定反式中间体trans-int，其经过渡态trans-ts分解成产物CH3NH和CO；第三是CH3与HNCO分子间通过生成稳定顺式中间体cis-int，其经过渡态cis-ts分解成产物CH3NH和CO.比较三条反应通道的反应活化能，表明CH3与HNCO反应较易生成CH4＋NCO.     

The Jahn-Teller effect in the triply degenerate electronic state of methane radical cation

A quantum dynamics study is performed to examine the complex nuclear motion underlying the first photoelectron band of methane. The broad and highly overlapping structures of the latter are found to originate from transitions to the ground electronic state, X(2)T(2), of the methane radical cation. Ab initio calculations have also been carried out to establish the potential energy surfaces for the triply degenerate electronic manifold of CH(4)(+). A suitable diabatic vibronic Hamiltonian has been devised and the nonadiabatic effects due to Jahn-Teller conical intersections on the vibronic dynamics investigated in detail. The theoretical results show fair accord with experiment.     

New method for calculating bound states: the A(1) states of Li(3) on the spin-aligned Li(3)(1 (4)A(')) potential energy surface

In this paper, we present a calculation for the bound states of A(1) symmetry on the spin-aligned Li(3)(1 (4)A(')) potential energy surface. We apply a mixture of discrete variable representation and distributed approximating functional methods to discretize the Hamiltonian. We also introduce a new method that significantly reduces the computational effort needed to determine the lowest eigenvalues and eigenvectors (bound state energies and wave functions of the full Hamiltonian). In our study, we have found the lowest 150 energy bound states converged to less than 0.005% error, and most of the excited energy bound states converged to less than 2.0% error. Furthermore, we have estimated the total number of the A(1) bound states of Li(3) on the spin-aligned Li(3)(1 (4)A(')) potential surface to be 601.     

Theoretical unimolecular kinetics for CH4 + M ⇄ CH3 + H + M in eight baths, M = He, Ne, Ar, Kr, H2, N2, CO, and CH4

Ensembles of classical trajectories are used to study collisional energy transfer in highly vibrationally excited CH(4) for eight bath gases. Several simplifying assumptions for the CH(4) + M interaction potential energy surface are tested against full dimensional direct dynamics trajectory calculations for M = He, Ne, and H(2). The calculated energy transfer averages are confirmed to be sensitive to the shape of the repulsive wall of the intermolecular potential, with an exponential repulsive wall required for quantitative predictions. For the diatomic baths, the usual "separable pairwise" approximation for the interaction potential is unable to describe the orientation dependence of the interaction potential accurately, and the ambiguity in the resulting parametrizations contributes an additional uncertainty to the predicted energy transfer averages of 20-40%. On the other hand, the energy transfer averages are shown to be insensitive to the level of theory used to describe the intramolecular CH(4) potential, with a computationally efficient semiempirical tight binding potential for hydrocarbons performing equally well as an MP2 potential. The relative collisional energy transfer efficiencies of the eight bath gases are discussed and shown to be a function of temperature. The ensemble-averaged energy transferred in deactivating collisions <ΔE(d)> for each bath is used to parametrize a single-exponential-down model for collisional energy transfer in master equation calculations. The predicted decomposition rate coefficients for CH(4) agree well with available experimental rate coefficients for M = He, Ar, Kr, and CH(4). The effect of vibrational anharmonicity on the predicted rate coefficients is considered briefly.     

Ab initio and direct quasiclassical-trajectory study of the F+CH4-->HF+CH3 reaction

We present an electronic structure and dynamics study of the F+CH4-->HF+CH3 reaction. CCSD(T)/aug-cc-pVDZ geometry optimizations, harmonic-frequency, and energy calculations indicate that the potential-energy surface is remarkably isotropic near the transition state. In addition, while the saddle-point F-H-C angle is 180 degrees using MP2 methods, CCSD(T) geometry optimizations predict a bent transition state, with a 153 degrees F-H-C angle. We use these high-quality ab initio data to reparametrize the parameter-model 3 (PM3) semiempirical Hamiltonian so that calculations with the improved Hamiltonian and employing restricted open-shell wave functions agree with the higher accuracy data. Using this specific-reaction-parameter PM3 semiempirical Hamiltonian (SRP-PM3), we investigate the reaction dynamics by propagating quasiclassical trajectories. The results of our calculations using the SRP-PM3 Hamiltonian are compared with experiments and with the estimates of two recently reported potential-energy surfaces. The trajectory calculations using the SRP-PM3 Hamiltonian reproduce quantitatively the measured HF vibrational distributions. The calculations also agree with the experimental HF rotational distributions and capture the essential features of the excitation function. The results of the SRP semiempirical Hamiltonian developed here clearly improve over those using the two prior potential-energy surfaces and suggest that reparametrization of semiempirical Hamiltonians is a promising strategy to develop accurate potential-energy surfaces for reaction dynamics studies of polyatomic systems.     

High resolution spectral analysis of 13CH3OH in the excited torsional states

In this work, we have extended our previous analysis of the Hamiltonian of 13C substituted methanol to include a large number of spectral lines involving the second excited torsional state using an improved model. The data set consisted of 2529 Fourier transform and microwave transitions with the rotational angular momentum J < or = 10, K < or = 6 and n < or = 2 (with 336 MW lines). The data set was fitted with the new Hamiltonian model to derive the molecular parameters. The results indicate that the model developed for the other methanol species (CH3OH, CH3(18)OH and CH3OD) is also valid for the C-13 substituted species. The results will allow the energy levels of the molecule to be calculated for higher torsional levels above the internal rotational barrier with improved precision and allow the analysis to be carried out for more excited torsional states.     

Ab initio characterization of the Mg-HF van der Waals complex

The equilibrium structure and the three-dimensional potential energy surface of the Mg-HF van der Waals complex in its ground electronic state have been determined from accurate ab initio calculations using the coupled-cluster method, CCSD(T), in conjunction with the basis sets of triple- through quintuple-zeta quality. The core-electron correlation, high-order valence-electron correlation, and scalar relativistic effects were investigated. The Mg-HF complex was confirmed to be linear at equilibrium, with a vibrationless dissociation energy (into Mg and HF) D(e) of 280 cm(-1). The vibration-rotation energy levels of two isotopologues, (24)Mg-HF and (24)Mg-DF, were predicted using the variational method. The predicted spectroscopic constants can be useful in a further analysis of high-resolution vibration-rotation spectra of the Mg-HF complex.