Jet cooled spectroscopy of H2DO+: Barrier heights and isotope-dependent tunneling dynamics from H3O+ to D3O+

The first high resolution spectroscopic data for jet cooled H2DO+ are reported, specifically via infrared laser direct absorption in the OH stretching region with a slit supersonic jet discharge source. Transitions sampling upper (0-) and lower (0+) tunneling states for both symmetric (nu1+ <-- 0+, nu1- <-- 0-, and nu1- <-- 0+) and antisymmetric (nu3+ <-- 0+ and nu3- <-- 0-) OH stretching bands are observed, where +/- refers to wave function reflection symmetry with respect to the planar umbrella mode transition state. The spectra can be well fitted to a Watson asymmetric top Hamiltonian, revealing band origins and rotational constants for benchmark comparison with high-level ab initio theory. Of particular importance are detection and assignment of the relatively weak band (nu1- <-- 0+) that crosses the inversion tunneling gap, which is optically forbidden in H3O+ or D3O+, but weakly allowed in H2DO+ by lowering of the tunneling transition state symmetry from D(3h) to C(2v). In conjunction with other H2DO+ bands, this permits determination of the tunneling splittings to within spectroscopic precision for each of the ground [40.518(10) cm(-1)], nu1 = 1 [32.666(6) cm(-1)], and nu3 = 1 [25.399(11) cm(-1)] states. A one-dimensional zero-point energy corrected potential along the tunneling coordinate is constructed from high-level ab initio CCSD(T) calculations (AVnZ, n = 3,4,5) and extrapolated to the complete basis set limit to extract tunneling splittings via a vibrationally adiabatic treatment. Perturbative scaling of the potential to match splittings for all four isotopomers permits an experimental estimate of DeltaV0 = 652.9(6) cm(-1) for the tunneling barrier, in good agreement with full six-dimensional ab initio results of Rajamaki, Miani, and Halonen (RMH) [J. Chem. Phys. 118, 10929 (2003)]. (DeltaV0 (RMH) = 650 cm(-1)). The 30%-50% decrease in tunneling splitting observed upon nu1 and nu3 vibrational excitations arises from an increase in OH stretch frequencies at the planar transition state, highlighting the transition between sp2 and sp3 hybridizations of the OHD bonds as a function of inversion bending angle.     

High-resolution infrared studies in slit supersonic discharges: CH2 stretch excitation of jet-cooled CH2Cl radical

First high-resolution infrared spectra are presented for jet-cooled CH2 35Cl and CH2 37Cl radicals in the symmetric (nu1) CH2 stretching mode. A detailed spectral assignment yields refined lower and upper state rotational constants, as well as fine structure spin-rotation parameters from least-squares fits to the sub-Doppler line shapes for individual transitions. The rotational constants are consistent with a nearly planar structure, but do not exclude substantial large amplitude bending motion over a small barrier to planarity accessible with zero-point excitation. High level coupled cluster (singles/doubles/triples) calculations, extrapolated to the complete basis set limit, predict a slightly nonplanar equilibrium structure (theta approximately 11 degrees), with a vibrationally adiabatic treatment of the bend coordinate yielding a v = 1<--0 anharmonic frequency (393 cm(-1)) in excellent agreement with matrix studies (nu(bend) approximately 400 cm(-1)). The antisymmetric CH2 stretch vibration is not observed despite high sensitivity detection (signal to noise ratio >20:1) in the symmetric stretch band. This is consistent with density functional theory intensity calculations indicating a >35-fold smaller antisymmetric stretch transition moment for CH2Cl, and yet contrasts dramatically with high-resolution infrared studies of CH2F radical, for which both symmetric and antisymmetric CH2 stretches are observed in a nearly 2:1 intensity ratio. A simple physical model is presented based on a competition between bond-dipole and "charge-sloshing" contributions to the transition moment, which nicely explains the trends in CH2X symmetric versus asymmetric stretch intensities as a function of electron withdrawing group (X = D,Br,Cl,F).     

Ground-state and vibrationally assisted tunneling in the formic acid dimer

The previously developed instanton theory [G. V. Mil'nikov and H. Nakamura, J. Chem. Phys. 122, 124311 (2005)] is applied to the calculation of vibrationally assisted tunneling splitting of the deuterated formic acid dimer (DCOOH)2 with all the degrees of freedom taken into account. The ground-state tunnel splitting is determined by the density-functional theory combined with coupled cluster level of quantum chemistry to be 0.0038 cm(-1) which is comparable to the experimental value of 0.0029 cm(-1). Further, the tunnel splittings of fundamental excitations are estimated for frequencies below 300 cm(-1). In this energy range it is found that the excitation modes may either enhance or suppress tunneling as compared to the ground state. For the higher-frequency modes a rapid growth of the tunnel splitting is observed. At frequencies above 1000 cm(-1) the semiclassical solution becomes unstable and no reliable tunneling splittings can be obtained. This is in vast contrast to the adiabatic approximation to the instanton theory in which the tunnel splittings can be retrieved up to 3000 cm(-1). We discuss this disparity from the viewpoint of the multidimensional character of tunneling in hydrogen bonds and the adiabatic approximation is concluded to be inaccurate.     

High-resolution infrared spectroscopy of jet-cooled vinyl radical: symmetric CH2 stretch excitation and tunneling dynamics

First high-resolution IR spectra of jet-cooled vinyl radical in the C-H stretch region are reported. Detailed spectral assignments and least squares fits to an A-reduction Watson asymmetric top Hamiltonian yield rotational constants and vibrational origins for three A-type bands, assigned to single quantum excitation of the symmetric CH(2) stretch. Two of the observed bands arise definitively from ground state vinyl radical, as rigorously confirmed by combination differences predicted from previous midinfrared CH(2) wagging studies of Kanamori et al. [J. Chem. Phys. 92, 197 (1990)] as well as millimeter wave rotation-tunneling studies of Tanaka et al. [J. Chem. Phys. 120, 3604 (2004)]. The two bands reflect transitions out of symmetric (0(+)) and antisymmetric (0(-)) tunneling levels of vinyl radical populated at 14 K slit-jet expansion temperatures. The band origins for the lower-lower (0(+)<--0(+)) and upper-upper (0(-)<--0(-)) transitions occur at 2901.8603(7) and 2901.9319(4) cm(-1), respectively, which indicates an increase in the tunneling splitting and therefore a decrease in the effective tunneling barrier upon CH(2) symmetric stretch excitation. The third A-type band with origin at 2897.2264(3) cm(-1) exhibits rotational constants quite close to (but at high-resolution distinguishable from) the vinyl radical ground state, consistent with a CH(2) symmetric stretch hot band built on one or more quanta of excitation in a low frequency vibration. The observed CH(2) symmetric stretch bands are in excellent agreement with anharmonically scaled high level density functional theory (DFT) calculations and redshifted considerably from previous low resolution assignments. Of particular dynamical interest, Boltzmann analysis indicates that the pair of 0(+) and 0(-) tunneling bands exhibits 1:1 nuclear spin statistics for K(a)=even:odd states. This differs from the expected 3:1 ratio for feasible exchange of the two methylenic H atoms but is consistent with a 4:4 ratio predicted for interchange between all three H atoms. This suggests the novel dynamical possibility of large amplitude "roaming" of all three H atoms in vinyl radical, promoted by high internal vibrational excitation arising from dissociative electron attachment in the discharge.     

Slit discharge IR spectroscopy of a jet-cooled cyclopropyl radical: structure and intramolecular tunneling dynamics

High-resolution infrared spectra of a jet-cooled cyclopropyl radical are reported for the first time, specifically sampling the in-phase antisymmetric CH2 stretch (nu7) vibration. In addition to yielding the first precise gas-phase structural information, the spectra reveal quantum level doubling into lower (+) and upper (-) states due to tunneling of the lone alpha-CH with respect to the CCC plane. The bands clearly reveal intensity alternation due to H atom nuclear spin statistics (6:10 and 10:6 for even:odd Ka+Kc in lower (+) and upper (-) tunneling levels, respectively) consistent with C2v symmetry of the cyclopropyl-tunneling transition state. The two ground-state-tunneling levels fit extremely well to a rigid asymmetric rotor Hamiltonian, but there is clear evidence for both local and global state mixing in the vibrationally excited nu7 tunneling levels. In particular, the upper (-) tunneling component of the nu7 state is split by anharmonic coupling with a nearly isoenergetic dark state, which thereby acquires oscillator strength via intensity sharing with this bright state. From thermal Boltzmann analysis of fractional populations, tunneling splittings for a cyclopropyl radical are estimated to be 3.2 +/- 0.3 cm(-1) and 4.9 +/- 0.3 cm(-1) in the ground and nu7-excited states, respectively. This analysis indicates ground-state stereoracemization of the alpha-CH radical center to be a very fast process [k approximately 2.0(4) x 10(11) s(-1)], with the increase in the tunneling rate upon CH2 in-phase asymmetric stretch excitation consistent with ab initio predictions of equilibrium vs transition-state zero-point energies. Modeling of the ground-state-tunneling splittings with high level ab initio 1D potentials indicates an improved V0 = 1115 +/- 35 cm(-1) barrier height for alpha-CH inversion through the cyclopropyl CCC plane.     

Ab initio and direct dynamics studies of the reaction of singlet methylene with acetylene and the lifetime of the cyclopropene complex

The energetics of the (1)CH(2) + C(2)H(2) --> H + C(3)H(3) reaction are accurately calculated using an extrapolated coupled-cluster/complete basis set (CBS) method based on the cc-pVDZ, cc-pVTZ, and cc-pVQZ basis sets. The reaction enthalpy (0 K) is predicted to be -20.33 kcal/mol. This reaction has no classical barrier in either the entrance or exit channel. However, there are several stable intermediates-cyclopropene (c-C(3)H(4)), allene (CH(2)CCH(2)), and propyne (CH(3)CCH)-along the minimum energy path. These intermediates with zero-point energy corrections lie below the reactants by 87.11 (c-C(3)H(4)), 109.69 (CH(2)CCH(2)), and 110.78 kcal/mol (CH(3)CCH). The vibrationally adiabatic ground-state (VAG) barrier height for c-C(3)H(4) isomerization to allene is obtained as 45.2 kcal/mol, and to propyne as 37.2 kcal/mol. In addition, the (1)CH(2) + C(2)H(2) reaction is investigated utilizing the dual-level "scaling all correlation" (SAC) ab initio method of Truhlar et al., i.e., the UCCSD(SAC)/cc-pVDZ theory. Results show that the reaction occurs via long-lived complexes. The lifetime of the cyclopropene intermediate is obtained as 3.2 +/- 0.4 ps. It is found that the intermediate propyne can be formed directly from reactants through the insertion of (1)CH(2) into a C-H bond of C(2)H(2). However, compared to the major mechanism in which the propyne is produced through a ring-opening of the cyclopropene complex, this reaction pathway is much less favorable. Finally, the theoretical thermal rate constant exhibits a negative temperature dependence, which is in excellent agreement with the previous results. The temperature dependence is consistent with the earlier RRKM results but weaker than the experimental observations at high temperatures.     

A new potential energy surface and predicted infrared spectra of He-CO2: dependence on the antisymmetric stretch of CO2

A new potential energy surface involving the antisymmetric Q(3) normal mode of CO(2) for the He-CO(2) van der Waals complex is constructed at the coupled-cluster singles and doubles with noniterative inclusion of connected triple [CCSD(T)] level with augmented correlation-consistent quadruple-zeta (aug-cc-pVQZ) basis set plus bond functions. Two vibrationally adiabatic potentials with CO(2) at both the ground and the first excited vibrational states are generated from the integration of the three-dimensional potential over the Q(3) coordinate. The potential has a T-shaped global minimum and two equivalent linear local minima. The bound rovibrational energy levels are obtained using the radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm. The observed band origin shift of the complex (0.0946 cm(-1)) is successfully reproduced by our calculation (0.1034 cm(-1)). The infrared spectra of the complex are also predicted. The fundamental band is in excellent agreement with the experiment. Most of the transitions corresponding to the observed hot band [M. J. Weida et al., J. Chem. Phys. 101, 8351 (1994)] are assigned reasonably.     

Protonated carbonyl sulfide: prospects for the spectroscopic observation of the elusive HSCO+ isomer

Spurred by the apparent conflict between ab initio predictions and infrared spectroscopic evidence regarding the relative stability of isomers of protonated carbonyl sulfide, key stationary points on the isomerization surface of HOCS(+) have been examined via systematic extrapolations of ab initio energies. Electron correlation has been accounted for using second-order M?ller-Plesset perturbation theory and coupled cluster theory through triple excitations [CCSD, CCSD(T), and CCSDT] in conjunction with the correlation consistent hierarchy of basis sets, cc-pVXZ (X=D,T,Q,5,6). HSCO(+) is predicted to lie lower in energy than HOCS(+) by 4.86 kcal mol(-1), computed using the focal point extrapolation scheme of Allen and co-workers [J. Chem. Phys. 99, 4638 (1993)] with corrections for anharmonic zero-point vibrational energy, core correlation, non-Born-Oppenheimer, and scalar relativistic effects. A transition state has been located, constituting the barrier to isomerization of HSCO(+) to HOCS(+), lying 68.9 kcal mol(-1) higher in energy than HSCO(+). This is well above predicted exothermicity [DeltaH(r) (o)(0 K)=48.1 kcal mol(-1), cc-pVQZ CCSD(T)] for the reaction considered in the experiments (HSCO(+)+H(2)-->OCS+H(3) (+)). Though proton tunneling will lead to a lower effective barrier, this prediction is consistent with the lack of HSCO(+) in electrical discharges in H(2)OCS, since the relative populations of HOCS(+) and HSCO(+) will depend on the experimental details of the protonation route rather than the relative thermodynamic stability of the isomers. Anharmonic vibrational frequencies and vibrationally corrected rotational constants from cc-pVTZ CCSD(T) cubic and quartic force constants are provided, to aid in the spectroscopic observation of the energetically favorable but apparently elusive HSCO(+) isomer.     

Quantum tunneling in the midrange vibrational fundamentals of tropolone

The Fourier transform infrared spectrum of tropolone(OH) vapor in the 1175-1700 cm(-1) region is reported at 0.0025 and 0.10 cm(-1) spectral resolutions. The 12 vibrational fundamentals in this region of rapidly rising vibrational state density are dominated by mixtures of the CC, CO, CCH, and COH internal coordinates. Estimates based on the measurement of sharp Q branch peaks are reported for 11 of the spectral doublet component separations DS(v) = |Delta(v) +/- Delta(0)|. Delta(0) = 0.974 cm(-1) is the known zero-point splitting, and three a(1) modes show tunneling splittings Delta(v) approximately Delta(0), four b(2) modes show splittings Delta(v) approximately 0.90Delta(0), and the remaining four modes show splittings Delta(v) falling 5-14% from Delta(0.) Significantly, the splitting for the nominal COH bending mode nu(8) (a(1)) is small, that is, 10% from Delta(0). Many of the vibrational excited states demonstrate strong anharmonic behavior, but there are only mild perturbations on the tautomerization mechanism driving Delta(0). The data suggest, especially for the higher frequency a(1) fundamentals, the onset of selective intramolecular vibrational energy redistribution processes that are fast on the time scale of the tautomerization process. These appear to delocalize and smooth out the topographical modifications of the zero-point potential energy surface that are anticipated to follow absorption of the nu(v) photon. Further, the spectra show the propensity for the Delta(v) splittings of b(2) and other complex vibrations to be damped relative to Delta(0).     

Bound state potential energy surface construction: ab initio zero-point energies and vibrationally averaged rotational constants

Collins' method of interpolating a potential energy surface (PES) from quantum chemical calculations for reactive systems (Jordan, M. J. T.; Thompson, K. C.; Collins, M. A. J. Chem. Phys. 1995, 102, 5647. Thompson, K. C.; Jordan, M. J. T.; Collins, M. A. J. Chem. Phys. 1998, 108, 8302. Bettens, R. P. A.; Collins, M. A. J. Chem. Phys. 1999, 111, 816) has been applied to a bound state problem. The interpolation method has been combined for the first time with quantum diffusion Monte Carlo calculations to obtain an accurate ground state zero-point energy, the vibrationally average rotational constants, and the vibrationally averaged internal coordinates. In particular, the system studied was fluoromethane using a composite method approximating the QCISD(T)/6-311++G(2df,2p) level of theory. The approach adopted in this work (a) is fully automated, (b) is fully ab initio, (c) includes all nine nuclear degrees of freedom, (d) requires no assumption of the functional form of the PES, (e) possesses the full symmetry of the system, (f) does not involve fitting any parameters of any kind, and (g) is generally applicable to any system amenable to quantum chemical calculations and Collins' interpolation method. The calculated zero-point energy agrees to within 0.2% of its current best estimate. A0 and B0 are within 0.9 and 0.3%, respectively, of experiment.     

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.     

Anharmonic vibrational frequencies and vibrationally averaged structures and nuclear magnetic resonance parameters of FHF-

The anharmonic vibrational frequencies of FHF(-) were computed by the vibrational self-consistent-field, configuration-interaction, and second-order perturbation methods with a multiresolution composite potential energy surface generated by the electronic coupled-cluster method with various basis sets. Anharmonic vibrational averaging was performed for the bond length and nuclear magnetic resonance indirect spin-spin coupling constants, where the latter computed by the equation-of-motion coupled-cluster method. The calculations placed the vibrational frequencies at 580 (nu(1)), 1292 (nu(2)), 1313 (nu(3)), 1837 (nu(1) + nu(3)), and 1864 cm(-1) (nu(1) + nu(2)), the zero-point H-F bond length (r(0)) at 1.1539 A, the zero-point one-bond spin-spin coupling constant [(1)J(0)(HF)] at 124 Hz, and the bond dissociation energy (D(0)) at 43.3 kcal/mol. They agreed excellently with the corresponding experimental values: nu(1) = 583 cm(-1), nu(2) = 1286 cm(-1), nu(3) = 1331 cm(-1), nu(1) + nu(3) = 1849 cm(-1), nu(1) + nu(2) = 1858 cm(-1), r(0) = 1.1522 A, (1)J(0)(HF) = 124+/-3 Hz, and D(0) = 44.4+/-1.6 kcal/mol. The vibrationally averaged bond lengths matched closely the experimental values of five excited vibrational states, furnishing a highly dependable basis for correct band assignments. An adiabatic separation of high- (nu(3)) and low-frequency (nu(1)) stretching modes was examined and found to explain semiquantitatively the appearance of a nu(1) progression on nu(3). Our calculations predicted a value of 186 Hz for experimentally inaccessible (2)J(0)(FF).     

Theoretical prediction of the ionization energies of the C4H7 radicals: 1-methylallyl, 2-methylallyl, cyclopropylmethyl, and cyclobutyl radicals

The ionization energies (IEs) for the 1-methylallyl, 2-methylallyl, cyclopropylmethyl, and cyclobutyl radicals have been calculated by the wave function based ab initio CCSD(T)/CBS approach, which involves the approximation to the complete basis set (CBS) limit at the coupled cluster level with single and double excitations plus quasiperturbative triple excitation [CCSD(T)]. The zero-point vibrational energy correction, the core-valence electronic correction, and the scalar relativistic effect correction are included in these calculations. The present CCSD(T)/CBS results are then compared with the IEs determined in the photoelectron experiment by Schultz et al. [J. Am. Chem. Soc. 106, 7336 (1984)] The predicted IE value (7.881 eV) of 2-methylallyl radical is found to compare very favorably with the experimental value of 7.90+/-0.02 eV. Two ionization transitions for cis-1-methylallyl and trans-1-methylallyl radicals have been considered here. The comparison between the predicted IE values and the previous measurements shows that the photoelectron peak observed by Schultz et al. likely corresponds to the adiabatic ionization transition for the trans-1-methylallyl radical to form trans-1-methylallyl cation. Although a precise IE value for the cyclopropylmethyl radical has not been directly determined, the experimental value deduced indirectly using other known energetic data is found to be in good accord with the present CCSD(T)/CBS prediction. We expect that the Franck-Condon factor for ionization transition of c-C4H7-->bicyclobutonium is much less favorable than that for ionization transition of c-C4H7-->planar-C4H7+, and the observed IE in the previous photoelectron experiment is likely due to the ionization transition for c-C4H7-->planar-C4H7+. Based on our CCSD(T)/CBS prediction, the ionization transition of c-C4H7-->bicyclobutonium with an IE value around 6.92 eV should be taken as the adiabatic ionization transition for the cyclobutyl radical. The present study provides support for the conclusion that the CCSD(T)/CBS approach with high-level energetic corrections can be used to provide reliable IE predictions for C4 hydrocarbon radicals with an uncertainty of +/-22 meV. The CCSD(T)/CBS predictions to the heats of formation for the aforementioned radicals and cations are also presented.     

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.     

Spectroscopy of the UO+2 cation and the delayed ionization of UO2

Vibronically resolved spectra for the UO+2 cation have been recorded using the pulsed field ionization zero electron kinetic energy (PFI-ZEKE) technique. For the ground state, long progressions in both the bending and symmetric stretch vibrations were observed. Bend and stretch progressions of the first electronically excited state were also observed, and the origin was found at an energy of 2678 cm(-1) above the ground state zero-point level. This observation is consistent with a recent theoretical prediction [Infante et al., J. Chem. Phys. 127, 124308 (2007)]. The ionization energy for UO2, derived from the PFI-ZEKE spectrum, namely, 6.127(1) eV, is in excellent agreement with the value obtained from an earlier photoionization efficiency measurement. Delayed ionization of UO2 in the gas phase has been reported previously [Han et al., J. Chem. Phys. 120, 5155 (2004)]. Here, we extend the characterization of the delayed ionization process by performing a quantitative study of the ionization rate as a function of the energy above the ionization threshold. The ionization rate was found to be 5 x 10(6) s(-1) at threshold, and increased linearly with increasing energy in the range investigated (0-1200 cm(-1)).     

Infrared laser spectroscopy of CH3...HF in helium nanodroplets: The exit-channel complex of the F + CH4 reaction

High-resolution infrared laser spectroscopy is used to study the CH3...HF and CD3...HF radical complexes, corresponding to the exit-channel complex in the F + CH4 --> HF + CH3 reaction. The complexes are formed in helium nanodroplets by sequential pickup of a methyl radical and a HF molecule. The rotationally resolved spectra presented here correspond to the fundamental v = 1 <-- 0 H-F vibrational band, the analysis of which reveals a complex with C(3v) symmetry. The vibrational band origin for the CH3...HF complex (3797.00 cm(-1)) is significantly redshifted from that of the HF monomer (3959.19 cm(-1)), consistent with the hydrogen-bonded structure predicted by theory [E. Ya. Misochko et al., J. Am. Chem. Soc. 117, 11997 (1995)] and suggested by previous matrix isolation experiments [M. E. Jacox, Chem. Phys. 42, 133 (1979)]. The permanent electric dipole moment of this complex is experimentally determined by Stark spectroscopy to be 2.4+/-0.3 D. The wide amplitude zero-point bending motion of this complex is revealed by the vibrational dependence of the A rotational constant. A sixfold reduction in the line broadening associated with the H-F vibrational mode is observed in going from CH3...HF to CD3...HF. The results suggest that fast relaxation in the former case results from near-resonant intermolecular vibration-vibration (V-V) energy transfer. Ab initio calculations are also reported (at the MP2 level) for the various stationary points on the F + CH4 surface, including geometry optimizations and vibrational frequency calculations for CH3...HF.     

Communication: rigorous calculation of dissociation energies (D0) of the water trimer, (H2O)3 and (D2O)3

Using a recent, full-dimensional, ab initio potential energy surface [Y. Wang, X. Huang, B. C. Shepler, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 134, 094509 (2011)] together with rigorous diffusion Monte Carlo calculations of the zero-point energy of the water trimer, we report dissociation energies, D(0), to form one monomer plus the water dimer and three monomers. The calculations make use of essentially exact zero-point energies for the water trimer, dimer, and monomer, and benchmark values of the electronic dissociation energies, D(e), of the water trimer [J. A. Anderson, K. Crager, L. Fedoroff, and G. S. Tschumper, J. Chem. Phys. 121, 11023 (2004)]. The D(0) results are 3855 and 2726 cm(-1) for the 3H(2)O and H(2)O + (H(2)O)(2) dissociation channels, respectively, and 4206 and 2947 cm(-1) for 3D(2)O and D(2)O + (D(2)O)(2) dissociation channels, respectively. The results have estimated uncertainties of 20 and 30 cm(-1) for the monomer plus dimer and three monomer of dissociation channels, respectively.     

Intramolecular proton transfer in malonaldehyde: accurate multilayer multi-configurational time-dependent Hartree calculations

Full-dimensional (multilayer) multi-configurational time-dependent Hartree calculations studying the intramolecular proton transfer in malonaldehyde based on a recent potential energy surface (PES) [Wang et al., J. Chem. Phys. 128, 224314 (2008)] are presented. The most accurate calculations yield a ground state tunneling splitting of 23.8 cm(-1) and a zero point energy of 14,678 cm(-1). Extensive convergence tests indicate an error margin of the quantum dynamics calculations for the tunneling splitting of about 0.2 cm(-1). These results are to be compared with the experimental value of the tunneling splitting of 21.58 cm(-1) and results of Monte Carlo calculations of Wang et al. on the same PES which yielded a zero point energy of 14,677.9 cm(-1) with statistical errors of 2-3 cm(-1) and a tunneling splitting of 21.6 cm(-1). The present data includes contributions resulting from the vibrational angular momenta to the tunneling splitting and the zero point energy of 0.2 cm(-1) and 2.4 cm(-1), respectively, which have been computed using a perturbative approach.     

Vibrational stimulation of the coherent tunneling transition in the cyclopentanone molecule

The tunneling interconversion of the cyclopentanone molecule, which leads to the appearance of tunneling doublets in the microwave spectrum of the system, is studied. The dynamics of interconversion is described by two generalized coordinates, one of which corresponds to bending (non-tunneling promoting mode), while the other of which corresponds to twisting of the molecular plane (tunneling coordinate). The coupling between two coordinates is symmetric. A method for quasi-classical calculation of the wave functions in the tunneling region and of the tunneling splittings of the vibrationally excited states in a two-dimensional potential with symmetric coupling is proposed. The tunneling spectrum of cyclopentanone is calculated. It agrees well with the experimental one, and the tunneling splitting increases by 140 times when the transverse quantum number goes from 0 to 6. The dynamic effect of the vibrationally assisted tunneling is shown to be due to the increase in the width of the tunneling channel with the quantum number of bending mode, as well as to the simultaneous shortening of the tunneling distance. The transition state geometry is found using the wave function at the dividing line of the potential.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 12, pp. 2098–2105, December, 1994.This work was supported by the Russian Foundation for Basic Research (Project 94-03-08863). The authors express their gratitude to W. Miller for helpful discussions and to H. Nakamura for a preprint of their work.     

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.