Theoretical studies of the tunneling splitting of malonaldehyde using the multiconfiguration time-dependent Hartree approach

Full dimensional multi-configuration time-dependent Hartree calculations of the zero point energy and the tunneling splitting of malonaldehyde using a recently published potential energy surface [Y. Wang, B. J. Braams, J. M. Bowman, S. Carter, and D. P. Tew, J. Chem. Phys. 128, 224314 (2008)] are reported. The potential energy surface has been approximated by a modified version of the n-mode representation and careful convergence check has been performed to ensure accurate results. The obtained value for the splitting (23.4 cm(-1)) is in acceptable agreement with the experimental value of 21.583 cm(-1). The computed zero-point-energy is 14,670 cm(-1) which is lower than previous results of Wang et al., but likely to be about 4 cm(-1) too low because of shortcomings of the n-mode representation of the potential. The energies reported in this abstract contain a correction to account for neglected vibrational angular momentum terms.     

Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcalmol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively.     

The ground state tunneling splitting and the zero point energy of malonaldehyde: a quantum Monte Carlo determination

Quantum dynamics calculations of the ground state tunneling splitting and of the zero point energy of malonaldehyde on the full dimensional potential energy surface proposed by Yagi et al. [J. Chem. Phys. 1154, 10647 (2001)] are reported. The exact diffusion Monte Carlo and the projection operator imaginary time spectral evolution methods are used to compute accurate benchmark results for this 21-dimensional ab initio potential energy surface. A tunneling splitting of 25.7+/-0.3 cm-1 is obtained, and the vibrational ground state energy is found to be 15 122+/-4 cm-1. Isotopic substitution of the tunneling hydrogen modifies the tunneling splitting down to 3.21+/-0.09 cm-1 and the vibrational ground state energy to 14 385+/-2 cm-1. The computed tunneling splittings are slightly higher than the experimental values as expected from the potential energy surface which slightly underestimates the barrier height, and they are slightly lower than the results from the instanton theory obtained using the same potential energy surface.     

Torsion-wagging tunneling and vibrational states in hydrazine determined from its ab initio potential energy surface

Geometries, anharmonic vibrations, and torsion-wagging (TW) multiplets of hydrazine and its deuterated species are studied using high-level ab initio methods employing the second-order Mo?ller-Plesset perturbation theory (MP2) as well as the coupled cluster singles and doubles model including connected triple corrections, CCSD(T), in conjunction with extended basis sets containing diffuse and core functions. To describe the splitting patterns caused by tunneling in TW states, the 3D potential energy surface (PES) for the large-amplitude TW modes is constructed. Stationary points in the 3D PES, including equivalent local minima and saddle points are characterized. Using this 3D PES, a flexible Hamiltonian is built numerically and then employed to solve the vibrational problem for TW coupled motion. The calculated ground state r(av) structure is expected to be more reliable than the experimental one that has been determined using a simplified structural model. The calculated fundamental frequencies allowed resolution of the assignment problems discussed earlier in the literature. The determined energy barriers, including the contributions from the small-amplitude vibrations, to the tunneling of the symmetric and antisymmetric wagging mode of 1997 cm(-1) and 3454 cm(-1), respectively, are in reasonable agreement with the empirical estimates of 2072 cm(-1) and 3312 cm(-1), respectively [W. ?odyga et al. J. Mol. Spectrosc. 183, 374 (1997)]. However, the empirical torsion barrier of 934 cm(-1) appears to be overestimated. The ab initio calculations yield two torsion barriers: cis and trans of 744 cm(-1) and 2706 cm(-1), respectively. The multiplets of the excited torsion states are predicted from the refined 3D PES.     

The ground state tunneling splitting of malonaldehyde: accurate full dimensional quantum dynamics calculations

Benchmark calculations of the tunneling splitting in malonaldehyde using the full dimensional potential proposed by Yagi et al. are reported. Two exact quantum dynamics methods are used: the multiconfigurational time-dependent Hartree (MCTDH) approach and the diffusion Monte Carlo based projection operator imaginary time spectral evolution (POITSE) method. A ground state tunneling splitting of 25.7+/-0.3 cm(-1) is calculated using POITSE. The MCTDH computation yields 25 cm(-1) converged to about 10% accuracy. These rigorous results are used to evaluate the accuracy of approximate dynamical approaches, e.g., the instanton theory.     

Large-amplitude dynamics in vinyl radical: the role of quantum tunneling as an isomerization mechanism

We report tunneling splittings associated with the large amplitude 1,2 H-atom migration to the global minima in the vinyl radical. These are obtained using a recent full-dimensional ab initio potential energy surface (PES) [A. R. Sharma, B. J. Braams, S. Carter, B. C. Shepler, and J. M. Bowman, J. Chem. Phys. 130(17), 174301 (2009)] and independently, directly calculated "reaction paths." The PES is a multidimensional fit to coupled cluster single and double and perturbative treatment of triple excitations coupled-cluster single double triple (CCSD(T)) with the augmented correlation consistent triple zeta basis set (aug-cc-pVTZ). The reaction path potentials are obtained from a series of CCSD(T)/aug-cc-pVnTZ calculations extrapolated to the complete basis set limit. Approximate 1D calculations of the tunneling splitting for these 1,2-H atom migrations are obtained using each of these potentials as well as quite different 1D Hamiltonians. The splittings are calculated over a large energy ranges, with results from the two sets of calculations in excellent agreement. Though negligibly slow (>1 s) for the vibrational ground state, this work predicts tunneling-promoted 1,2 hydride shift dynamics in vinyl to exhibit exponential growth with internal vibrational excitation, specifically achieving rates on the sub-μs time scale at energies above E ≈ 7500 cm(-1). Most importantly, these results begin to elucidate the possible role of quantum isomerization through barriers without dissociation, in competition with the more conventional picture of classical roaming permitted over a much narrower window of energies immediately below the bond dissociation limit. Furthermore, when integrated over a Boltzmann distribution of thermal energies, these microcanonical tunneling rates are consistent with sub-μs time scales for 1,2 hydride shift dynamics at T > 1400 K. These results have potential relevance for combustion modeling of low-pressure flames, as well as recent observations of nuclear spin statistical mixing from high-resolution IR/microwave spectroscopy on vinyl radical.     

The lowest-lying electronic singlet and triplet potential energy surfaces for the HNO-NOH system: energetics, unimolecular rate constants, tunneling and kinetic isotope effects for the isomerization and dissociation reactions

The lowest-lying electronic singlet and triplet potential energy surfaces (PES) for the HNO-NOH system have been investigated employing high level ab initio quantum chemical methods. The reaction energies and barriers have been predicted for two isomerization and four dissociation reactions. Total energies are extrapolated to the complete basis set limit applying focal point analyses. Anharmonic zero-point vibrational energies, diagonal Born-Oppenheimer corrections, relativistic effects, and core correlation corrections are also taken into account. On the singlet PES, the (1)HNO → (1)NOH endothermicity including all corrections is predicted to be 42.23 ± 0.2 kcal mol(-1). For the barrierless decomposition of (1)HNO to H + NO, the dissociation energy is estimated to be 47.48 ± 0.2 kcal mol(-1). For (1)NOH → H + NO, the reaction endothermicity and barrier are 5.25 ± 0.2 and 7.88 ± 0.2 kcal mol(-1). On the triplet PES the reaction energy and barrier including all corrections are predicted to be 7.73 ± 0.2 and 39.31 ± 0.2 kcal mol(-1) for the isomerization reaction (3)HNO → (3)NOH. For the triplet dissociation reaction (to H + NO) the corresponding results are 29.03 ± 0.2 and 32.41 ± 0.2 kcal mol(-1). Analogous results are 21.30 ± 0.2 and 33.67 ± 0.2 kcal mol(-1) for the dissociation reaction of (3)NOH (to H + NO). Unimolecular rate constants for the isomerization and dissociation reactions were obtained utilizing kinetic modeling methods. The tunneling and kinetic isotope effects are also investigated for these reactions. The adiabatic singlet-triplet energy splittings are predicted to be 18.45 ± 0.2 and 16.05 ± 0.2 kcal mol(-1) for HNO and NOH, respectively. Kinetic analyses based on solution of simultaneous first-order ordinary-differential rate equations demonstrate that the singlet NOH molecule will be difficult to prepare at room temperature, while the triplet NOH molecule is viable with respect to isomerization and dissociation reactions up to 400 K. Hence, our theoretical findings clearly explain why (1)NOH has not yet been observed experimentally.     

Concerted hydrogen exchange tunneling in formic acid dimer

The effect of conformational relaxation on the quantum dynamics of the hydrogen exchange tunneling is studied in the D2h subspace of formic acid dimer. The fully coupled quantum dynamics in up to six dimensions are derived for potential energy hypersurfaces interpolated directly from hybrid density functional calculations with and without geometry relaxation. For a calculated electronic barrier height of 35.0 kJ/mol the vibrational ground state shows a tunneling splitting of 0.0013 cm(-1). The results support the vibrational assignment of Madeja and Havenith [J. Chem. Phys. 2002, 117, 7162-7168]. Fully coupled ro-vibrational calculations demonstrate the compatibility of experimentally observed inertia defects with in-plane hydrogen exchange tunneling dynamics in formic acid dimer.     

Experimental and theoretical differential cross sections for the N(2D) + H2 reaction

In this paper, we report a combined experimental and theoretical study on the dynamics of the N(2D) + H2 insertion reaction at a collision energy of 15.9 kJ mol(-1). Product angular and velocity distributions have been obtained in crossed beam experiments and simulated by using the results of quantum mechanical (QM) scattering calculations on the accurate ab initio potential energy surface (PES) of Pederson et al. (J. Chem. Phys. 1999, 110, 9091). Since the QM calculations indicate that there is a significant coupling between the product angular and translational energy distributions, such a coupling has been explicitly included in the simulation of the experimental results. The very good agreement between experiment and QM calculations sustains the accuracy of the NH2 ab initio ground state PES. We also take the opportunity to compare the accurate QM differential cross sections with those obtained by two approximate methods, namely, the widely used quasiclassical trajectory calculations and a rigorous statistical method based on the coupled-channel theory.     

Potential energy surface of HDO up to 25,000 cm-1

A new spectroscopically determined potential energy surface (PES) for HD(16)O is presented. This surface is constructed by adjusting the high accuracy ab initio PES of Polyansky et al. [Science 299, 539 (2003)] by fitting to both published experimental data and our still unpublished data. This refinement used experimentally derived term values up to 25,000 cm(-1) and with J< or =8: a data set of 3478 energy levels once some levels with ambiguous assignment is excluded. To improve the extrapolation properties of the empirical PES, the restraint that the resulting PESs remain close to the ab initio surface was imposed. The new HDO_07 PES reproduces the experimental data, including high J levels not included in the fit, with a root mean square error of 0.035 cm(-1). Predictions for rotation-vibration term values up to J=12 are made.     

Variational calculation of specific, highly excited vibrational states in DFCO: comparison with experimental data

A stimulated emission pumping spectra of jet-cooled DFCO performed by Crane et al. (J. Mol. Spectrosc. 1997, 183, 273) has provided a great number of ro-vibrational lines up to 9000 cm(-1) of excitation energy. By combining a Jacobi-Wilson (JW) approach with a Davidson scheme, we calculate the lines provided by the experiment up to 9000 cm(-1) using an ab initio global potential energy surface (PES) developed by Kato et al. (J. Chem. Phys. 1997, 107, 6114). Comparisons between experimental and calculated data provide a critical test of the quality of the PES used. We show that the variational calculated energies can be efficiently corrected by taking into account the error observed for the A' fundamental transitions nu(i) (i = 1, ..., 5) and the first overtone 2nu(6). A detailed analysis of the eigenstates obtained by the calculation allows one to quantify the coupling between the different modes. Such an information is essential to understand and predict the energy flow through a DFCO molecule that is initially excited.     

Intermolecular potential energy surface and spectra of He-HCl with generalization to other rare gas-hydrogen halide complexes

A two-dimensional (rigid monomer) intermolecular potential energy surface (PES) of the He-HCl complex has been obtained from ab initio calculations utilizing the symmetry-adapted perturbation theory (SAPT) and an spdfg basis set including midbond functions. The bond length in HCl was chosen to be equal to the expectation value in the ground vibrational state of isolated HCl. The rigid-monomer potential should be a very good approximation to the complete (three-dimensional) potential for H-Cl distances corresponding to the lowest vibrational levels of the monomer since the He-HCl interaction energy was found to be only weakly dependent on the HCl bond length in this region, at least as compared to systems such as Ar-HF. The calculated points were fitted using an analytic function with ab initio computed asymptotic coefficients. As expected, the complex is loosely bound, with the dispersion energy providing the majority of the attraction. Our SAPT PES agrees with the semiempirical PES of Willey et al. [J. Chem. Phys. 96, 898 (1992)], in finding that, atypically for rare gas-hydrogen halide complexes including the lighter halide atoms, the global minimum is on the Cl side (with intermonomer separation 3.35 A and depth of 32.8 cm(-1)), rather than on the H side, where there is only a local minimum (3.85 A, 30.8 cm(-1)). The ordering of the minima was confirmed by single-point calculations in larger basis sets and complete basis set extrapolations, and also using higher levels of theory. We show that the opposite findings in the recent calculations of Zhang and Shi [J. Mol. Struct: THEOCHEM 589, 89 (2002)] are due to the lack of midbond functions in their basis set. Despite the closeness in depth of the two linear minima, the existence of a relatively high barrier between them invalidates the assumption of isotropy, a feature of some literature potentials. The trends concerning the locations of minima within the family of rare gas-hydrogen halide complexes are rationalized in terms of the physical components of the intermolecular forces and related to monomer properties. The accuracy of the SAPT PES was tested by performing calculations of rovibrational levels. The transition frequencies obtained were found to be in excellent agreement (to within 0.02 cm(-1)) with the measurements of Lovejoy and Nesbitt [J. Chem. Phys. 93, 5387 (1990)]. The SAPT PES predicts a dissociation energy for the complex of 7.74 cm(-1) which is probably more accurate than the experimental value of 10.1+/-1.2 cm(-1). Our analysis of the ground-state rovibrational wave function shows that the He-HCl configuration is favored over the He-ClH configuration despite the ordering of minima. This is due to the greater volume of the well in the former case. We have also determined positions and widths of three low-lying resonance states through scattering calculations. These predictions are expected to be more accurate than values derived from experiment.     

A reaction surface Hamiltonian study of malonaldehyde

We report calculations using a reaction surface Hamiltonian for which the vibrations of a molecule are represented by 3N-8 normal coordinates, Q, and two large amplitude motions, s(1) and s(2). The exact form of the kinetic energy operator is derived in these coordinates. The potential surface is first represented as a quadratic in Q, the coefficients of which depend upon the values of s(1),s(2) and then extended to include up to Q(6) diagonal anharmonic terms. The vibrational energy levels are evaluated by solving the variational secular equations, using a basis of products of Hermite polynomials and appropriate functions of s(1),s(2). Our selected example is malonaldehyde (N=9) and we choose as surface parameters two OH distances of the migrating H in the internal hydrogen transfer. The reaction surface Hamiltonian is ideally suited to the study of the kind of tunneling dynamics present in malonaldehyde. Our results are in good agreement with previous calculations of the zero point tunneling splitting and in general agreement with observed data. Interpretation of our two-dimensional reaction surface states suggests that the OH stretching fundamental is incorrectly assigned in the infrared spectrum. This mode appears at a much lower frequency in our calculations due to substantial transition state character.     

Direct versus resonances mediated F+OH collisions on a new 3A" potential energy surface

A theoretical study of the F(2P) + OH(2Pi) --> HF(1Sigma+) + O(3P) reactive collisions is carried out on a new global potential energy surface (PES) of the ground 3A" adiabatic electronic state. The ab initio calculations are based on multireference configuration interaction calculations, using the aug-cc-pVTZ extended basis sets of Dunning et al. A functional representation of the PES shows no nominal barrier to reaction, contrary to previous results by others. Wave packet and quasiclassical trajectory calculations have been performed for this PES to study the F + OH(v = 0,j) reactive collision. The comparison was performed at fixed and constant values of the total angular momentum from 0 to 110 and relative translational energy up to 0.8 eV. The reaction presents a dynamical barrier, essentially due to the zero-point energy for the bending vibration near the saddle point. This determines two different reaction mechanisms. At energies higher than approximately 0.125 eV the reaction is direct, while below that value it is indirect and mediated by heavy-light-heavy resonances. Such resonances, also found in the simulations of the photodetachment spectrum of the triatomic anion, manifest themselves in the quasiclassical simulations, too, where they are associated to periodic orbits.     

Improved dissociation energy of the 39K85Rb molecule

The predissociation data for the 1 (1)Pi state of (39)K(85)Rb of Kasahara et al. [J. Chem. Phys. 111, 8857 (1999)] are combined with the recent determination of the long range C(6) coefficients of the predissociating 2 (3)Sigma(+) approximately 2(0(-)), 2(1) states [Wang et al., Eur. Phys. J. D31, 165 (2004) ] molecule: to infer a more precise dissociation energy of the (39)K(85)Rb molecule D(0)=4180.06+/-0.42 cm(-1) and D(e)=4217.91+/-0.42 cm(-1).     

Quantum dynamics of H2, D2, and HD in the small dodecahedral cage of clathrate hydrate: evaluating H2-water nanocage interaction potentials by comparison of theory with inelastic neutron scattering experiments

We have performed rigorous quantum five-dimensional (5D) calculations and analysis of the translation-rotation (T-R) energy levels of one H(2), D(2), and HD molecule inside the small dodecahedral (H(2)O)(20) cage of the structure II clathrate hydrate, which was treated as rigid. The H(2)- cage intermolecular potential energy surface (PES) used previously in the molecular dynamics simulations of the hydrogen hydrates [Alavi et al., J. Chem. Phys. 123, 024507 (2005)] was employed. This PES, denoted here as SPC/E, combines an effective, empirical water-water pair potential [Berendsen et al., J. Phys. Chem. 91, 6269 (1987)] and electrostatic interactions between the partial charges placed on H(2)O and H(2). The 5D T-R eigenstates of HD were calculated also on another 5D H(2)-cage PES denoted PA-D, used by us earlier to investigate the quantum T-R dynamics of H(2) and D(2) in the small cage [Xu et al., J. Phys. Chem. B 110, 24806 (2006)]. In the PA-D PES, the hydrogen-water pair potential is described by the ab initio 5D PES of the isolated H(2)-H(2)O dimer. The quality of the SPC/E and the PA-D H(2)-cage PESs was tested by direct comparison of the T-R excitation energies calculated on them to the results of two recent inelastic neutron scattering (INS) studies of H(2) and HD inside the small clathrate cage. The translational fundamental and overtone excitations, as well as the triplet splittings of the j=0-->j=1 rotational transitions, of H(2) and HD in the small cage calculated on the SPC/E PES agree very well with the INS results and represent a significant improvement over the results computed on the PA-D PES. Our calculations on the SPC/E PES also make predictions about several spectroscopic observables for the encapsulated H(2), D(2), and HD, which have not been measured yet.     

An eight-degree-of-freedom, time-dependent quantum dynamics study for the H2+C2H reaction on a new modified potential energy surface

An eight-dimensional time-dependent quantum dynamics wave packet approach is performed for the study of the H2+C2H-->H+C2H2 reaction system on a new modified potential energy surface (PES) [L.-P. Ju et al., Chem. Phys. Lett. 409, 249 (2005)]. This new potential energy surface is obtained by modifying Wang and Bowman's old PES [J. Chem. Phys. 101, 8646 (1994)] based on the new ab initio calculation. This new modified PES has a much lower transition state barrier height at 2.29 kcal/mol than Wang and Bowman's old PES at 4.3 kcal/mol. This study shows that the reactivity for this diatom-triatom reaction system is enhanced by vibrational excitations of H2, whereas the vibrational excitations of C2H only have a small effect on the reactivity. Furthermore, the bending excitations of C2H, compared to the ground state reaction probability, hinder the reactivity. The comparison of the rate constant between this calculation and experimental results agrees with each other very well. This comparison indicates that the new modified PES corrects the large barrier height problem in Wang and Bowman's old PES.     

He-ThO(1Σ+) interactions at low temperatures: elastic and inelastic collisions, transport properties, and complex formation in cold 4He gas

We present an ab initio study of cold (4)He + ThO((1)Σ(+)) collisions based on an accurate potential energy surface (PES) evaluated by the coupled cluster method with single, double, and noniterative triple excitations using an extended basis set augmented by bond functions. Variational calculations of rovibrational energy levels show that the (4)He-ThO van der Waals complex has a binding energy of 10.9 cm(-1) in its ground J = 0 rotational state. The calculated energy levels are used to obtain the temperature dependence of the chemical equilibrium constant for the formation of the He-ThO complex. We find that complex formation is thermodynamically favored at temperatures below 1 K and predict the maximum abundance of free ground-state ThO(v = 0, j = 0) molecules between 2 and 3 K. The calculated cross sections for momentum transfer in elastic He + ThO collisions display a rich resonance structure below 5 cm(-1) and decline monotonically above this collision energy. The cross sections for rotational relaxation accompanied by momentum transfer decline abruptly to zero at low collision energies (<0.1 cm(-1)). We find that Stark relaxation in He + ThO collisions can be enhanced by applying an external dc electric field of less than 100 kV∕cm. Finally, we present calculations of thermally averaged diffusion cross sections for ThO in He gas, and find these to be insensitive to small variations of the PES at temperatures above 1 K.     

Calibration-quality adiabatic potential energy surfaces for H3(+) and its isotopologues

Calibration-quality ab initio adiabatic potential energy surfaces (PES) have been determined for all isotopologues of the molecular ion H(3)(+). The underlying Born-Oppenheimer electronic structure computations used optimized explicitly correlated shifted Gaussian functions. The surfaces include diagonal Born-Oppenheimer corrections computed from the accurate electronic wave functions. A fit to the 41,655 ab initio points is presented which gives a standard deviation better than 0.1 cm(-1) when restricted to the points up to 6000 cm(-1) above the first dissociation asymptote. Nuclear motion calculations utilizing this PES, called GLH3P, and an exact kinetic energy operator given in orthogonal internal coordinates are presented. The ro-vibrational transition frequencies for H(3)(+), H(2)D(+), and HD(2)(+) are compared with high resolution measurements. The most sophisticated and complete procedure employed to compute ro-vibrational energy levels, which makes explicit allowance for the inclusion of non-adiabatic effects, reproduces all the known ro-vibrational levels of the H(3)(+) isotopologues considered to better than 0.2 cm(-1). This represents a significant (order-of-magnitude) improvement compared to previous studies of transitions in the visible. Careful treatment of linear geometries is important for high frequency transitions and leads to new assignments for some of the previously observed lines. Prospects for further investigations of non-adiabatic effects in the H(3)(+) isotopologues are discussed. In short, the paper presents (a) an extremely accurate global potential energy surface of H(3)(+) resulting from high accuracy ab initio computations and global fit, (b) very accurate nuclear motion calculations of all available experimental line data up to 16,000 cm(-1), and (c) results suggest that we can predict accurately the lines of H(3)(+) towards dissociation and thus facilitate their experimental observation.