Ultrafast vibrational spectroscopy of water and aqueous N-methylacetamide: Comparison of different electronic structure/molecular dynamics approaches

Kwac and Cho [J. Chem. Phys. 119, 2247 (2003)] have recently developed a combined electronic structure/molecular dynamics approach to vibrational spectroscopy in liquids. The method involves fitting ab initio vibrational frequencies for a solute in a cluster of solvent molecules to a linear combination of the electrostatic potentials on the solute atoms due to the charges on the solvent molecules. These authors applied their method to the N-methylacetamide-D/D(2)O system. We (S. A. Corcelli, C. P. Lawrence, and J. L. Skinner, [J. Chem. Phys. 120, 8107 (2004)]) have recently explored a closely related method, where instead of the electrostatic potential, the solute vibrational frequencies are fit to the components of the electric fields on the solute atoms due to the solvent molecules. We applied our method to the HOD/D(2)O and HOD/H(2)O systems. In order to make a direct comparison of these two approaches, in this paper we apply their method to the water system, and our method to the N-methylacetamide system. For the water system we find that the electric field method is superior to the potential approach, as judged by comparison with experiments for the absorption line shape. For the N-methylacetamide system the two methods are comparable.     

Infrared diode laser spectroscopy of the Ne-D2O van der Waals complex: strong Coriolis and angular-radial coupling

Four internal-rotation/vibration bands of the Ne-D(2)O complex have been measured in the v(2) bend region of D(2)O using a tunable infrared diode laser spectrometer to probe a slit supersonic expansion. Three ortho bands are excited from the ground state Σ(0(00)) to the Σ and Π(1(11), υ(2) = 1) internal rotor states and the n = 1, Σ(0(00), υ(2) = 1) stretching-internal rotor combination state. Strong perturbations between the excited vibrational states are evident. The observed spectra are analyzed separately with a three-state J-dependent Coriolis plus J-independent angular-radial coupling model [M. J. Weida and D. J. Nesbitt, J. Chem. Phys. 106, 3078 (1997)] and a three-state Coriolis coupling model [R. C. Cohen and R. J. Saykally, J. Chem. Phys. 95, 7891 (1991)]. The former model works more successfully than the latter. Molecular constants for the ground and excited vibrational states of ortho (20)Ne-D(2)O isotopomer as well as the Coriolis and angular-radial coupling constants are determined accurately. The van der Waals stretching frequency is estimated to be ν(s) = 24.85 cm(-1) in the ground state and decreases to about 20.8 cm(-1) upon vibrational excitation of the D(2)O bend.     

Quantum Monte Carlo calculations of the dissociation energy of the water dimer

We report diffusion quantum Monte Carlo (DMC) calculations of the equilibrium dissociation energy D(e) of the water dimer. The dissociation energy measured experimentally, D(0), can be estimated from D(e) by adding a correction for vibrational effects. Using the measured dissociation energy and the modern value of the vibrational energy Mas et al., [J. Chem. Phys. 113, 6687 (2000)] leads to D(e)=5.00+/-0.7 kcal mol(-1), although the result Curtiss et al., [J. Chem. Phys. 71, 2703 (1979)] D(e)=5.44+/-0.7 kcal mol(-1), which uses an earlier estimate of the vibrational energy, has been widely quoted. High-level coupled cluster calculations Klopper et al., [Phys. Chem. Chem. Phys. 2, 2227 (2000)] have yielded D(e)=5.02+/-0.05 kcal mol(-1). In an attempt to shed new light on this old problem, we have performed all-electron DMC calculations on the water monomer and dimer using Slater-Jastrow wave functions with both Hartree-Fock approximation (HF) and B3LYP density functional theory single-particle orbitals. We obtain equilibrium dissociation energies for the dimer of 5.02+/-0.18 kcal mol(-1) (HF orbitals) and 5.21+/-0.18 kcal mol(-1) (B3LYP orbitals), in good agreement with the coupled cluster results.     

The relaxation of OH (v=1) and OD (v=1) by H2O and D2O at temperatures from 251 to 390 K

We report rate coefficients for the relaxation of OH(v=1) and OD(v=1) by H2O and D2O as a function of temperature between 251 and 390 K. All four rate coefficients exhibit a negative dependence on temperature. In Arrhenius form, the rate coefficients for relaxation (in units of 10(-12) cm3 molecule-1 s-1) can be expressed as: for OH(v=1)+H2O between 263 and 390 K: k=(2.4+/-0.9) exp((460+/-115)/T); for OH(v=1)+D2O between 256 and 371 K: k=(0.49+/-0.16) exp((610+/-90)/T); for OD(v=1)+H2O between 251 and 371 K: k=(0.92+/-0.16) exp((485+/-48)/T); for OD(v=1)+D2O between 253 and 366 K: k=(2.57+/-0.09) exp((342+/-10)/T). Rate coefficients at (297+/-1 K) are also reported for the relaxation of OH(v=2) by D2O and the relaxation of OD(v=2) by H2O and D2O. The results are discussed in terms of a mechanism involving the formation of hydrogen-bonded complexes in which intramolecular vibrational energy redistribution can occur at rates competitive with re-dissociation to the initial collision partners in their original vibrational states. New ab initio calculations on the H2O-HO system have been performed which, inter alia, yield vibrational frequencies for all four complexes: H2O-HO, D2O-HO, H2O-DO and D2O-DO. These data are then employed, adapting a formalism due to Troe (J. Troe, J. Chem. Phys., 1977, 66, 4758), in order to estimate the rates of intramolecular energy transfer from the OH (OD) vibration to other modes in the complexes in order to explain the measured relaxation rates-assuming that relaxation proceeds via the hydrogen-bonded complexes.     

Spectral diffusion in a fluctuating charge model of water

We apply the combined electronic structure/molecular dynamics approach of Corcelli, Lawrence, and Skinner [J. Chem. Phys. 120, 8107 (2004)] to the fluctuating charge (SPC-FQ) model of liquid water developed by Rick, Stuart, and Berne [J. Chem. Phys. 101, 6141 (1994)]. For HOD in H(2)O the time scale for the long-time decay of the OD stretch frequency time-correlation function, which corresponds to the time scale for hydrogen-bond rearrangement in the liquid, is about 1.5 ps. This result is significantly longer than the 0.9 ps decay previously calculated for the nonpolarizable SPC/E water model. Our results for the SPC-FQ model are in better agreement with recent vibrational echo experiments.     

Vibrational spectra of Na, K, Mn2+, Ni2+ and Zn2+ salts of 1,2,4,5-benzenetetracarboxylic (pyromellitic) acid--a short hydrogen bond evidence

Five salts of 1,2,4,5-benzenetetracarboxylic acid (pyromellitic acid), [C6H2(COO)4H4], have been synthesized and investigated by infrared and Raman spectroscopy and by single crystal X-ray diffraction methods: sodium salt [Na2(H2O)2][C6H2(COO)4H2], potassium salt [K(H2O)3][C6H2(COO)4H3] and transition metal salts [M(H2O)6][C6H2(COO)4H2], which M = Mn, Ni and Zn. Crystal structures of all five compounds show short intramolecular asymmetric hydrogen bonds (SHB) between adjacent carboxyl groups with O...O distance average of 2.40 A. The Raman and infrared spectra reported indicate the presence of short hydrogen bonds in all salts, in agreement with the X-ray data. The O-H stretching mode [nu(OH)] had been observed at about 2500 cm(-1). Deuterated analogues were synthesized and their Raman spectra show that nu(OH)/nu(OD) ratio average is about unit. The symmetric [nu(sym)(O..H..O)] and asymmetric [nu(asym)(O..H..O)] stretching modes have been attributed about 300 and 870 cm(-1), respectively, in all salts, and for deuterated analogues, the ratio nu(OH)/nu(OD) to nu(sym)(O..H..O, O..D..O) is close to unit like it occurs in nu(OH). The vibrational modes, mainly SHB modes, are tentatively assigned by molecular orbital ab initio calculations of pyromellitic acid and anions [C6H2(COO)4H3]- and [C6H2(COO)4H2]2-. Geometry optimizations showed a good agreement with experimental data. Frequency calculation confirms the assignment of specific vibrational modes. Ab initio calculations show that nu(C=O) and nu(sym)(COO) are strongly coupled with in plane OH bending [delta(OH)]. In Raman spectra of deuterated analogues is observed a frequency shift of these bands.     

The I2 dissociation mechanisms in the chemical oxygen-iodine laser revisited

The recently suggested mechanism of I(2) dissociation in the chemical oxygen-iodine laser (COIL) [K. Waichman, B. D. Barmashenko, and S. Rosenwaks, J. Appl. Phys. 106, 063108 (2009); and J. Chem. Phys. 133, 084301 (2010)] was largely based on the suggestion of V. N. Azyazov, S. Yu. Pichugin, and M. C. Heaven [J. Chem. Phys. 130, 104306 (2009)] that the vibrational population of O(2)(a) produced in the chemical generator is high enough to play an essential role in the dissociation. The results of model calculations based on this mechanism agreed very well with measurements of the small signal gain g, I(2) dissociation fraction F, and temperature T in the COIL. This mechanism is here revisited, following the recent experiments of M. V. Zagidullin [Quantum Electron. 40, 794 (2010)] where the observed low population of O(2)(b, v = 1) led to the conclusion that the vibrational population of O(2)(a) at the outlet of the generator is close to thermal equilibrium value. This value corresponds to a very small probability, ～0.05, of O(2)(a) energy pooling to the states O(2)(X,a,b, v > 0). We show that the dissociation mechanism can reproduce the experimentally observed values of g, F, and T in the COIL only if most of the energy released in the processes of O(2)(a) energy pooling and O(2)(b) quenching by H(2)O ends up as vibrational energy of the products, O(2)(X,a,b), where the vibrational states v = 2 and 3 are significantly populated. We discuss possible reasons for the differences in the suggested vibrational population and explain how these differences can be reconciled.     

The vibrational predissociation spectra of the H5O2 +RGn(RG = Ar,Ne) clusters: correlation of the solvent perturbations in the free OH and shared proton transitions of the Zundel ion

Predissociation spectra of the H(5)O(2) (+)RG(n)(RG = Ar,Ne) cluster ions are reported in energy regions corresponding to both the OH stretching (3350-3850 cm(-1)) and shared proton (850-1950 cm(-1)) vibrations. The two free OH stretching bands displayed by the Ne complex are quite close to the band origins identified earlier in bare H(5)O(2) (+) [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)], indicating that the symmetrical H(5)O(2) (+) "Zundel" ion remains largely intact in H(5)O(2) (+)Ne. The low-energy spectrum of the Ne complex is simpler than that observed previously for H(5)O(2) (+)Ar, and is dominated by two sharp transitions at 928 and 1047 cm(-1), with a weaker feature at 1763 cm(-1). The H(5)O(2) (+)Ar(n),n = 1-5 spectra generally exhibit complex band structures reflecting solvent-induced symmetry breaking of the Zundel core ion. The extent of solvent perturbation is evaluated with electronic structure calculations, which predict that the rare gas atoms should attach to the spectator OH groups of H(5)O(2) (+) rather than to the shared proton. In the asymmetric complexes, the shared proton resides closer to the more heavily solvated water molecule, leading to redshifts in the rare gas atom-solvated OH stretches and to blueshifts in the shared proton vibrations. The experimental spectra are compared with recent full-dimensional vibrational calculations (diffusion Monte Carlo and multimode/vibrational configuration interaction) on H(5)O(2) (+). These results are consistent with assignment of the strong low-energy bands in the H(5)O(2) (+)Ne spectrum to the vibration of the shared proton mostly along the O-O axis, with the 1763 cm(-1) band traced primarily to the out-of-phase, intramolecular bending vibrations of the two water molecules.     

Complexes of type C6H7(+)·L (L = N2 and CO2) studied by explicitly correlated coupled cluster theory

Complexes of the benzenium ion (C(6)H(7)(+)) with N(2) or CO(2) have been studied by explicitly correlated coupled cluster theory at the CCSD(T)-F12x (x = a, b) level [T. B. Adler et al., J. Chem. Phys. 127, 221106 (2007)] and the double-hybrid density functional B2PLYP-D [T. Schwabe and S. Grimme, Phys. Chem. Chem. Phys. 9, 3397 (2007)]. Improved harmonic vibrational wavenumbers for C(6)H(7)(+) have been obtained by CCSD(T?)-F12a calculations with the VTZ-F12 basis set. Combining them with previous B2PLYP-D anharmonic contributions we arrive at anharmonic wavenumbers which are in excellent agreement with recent experimental data from p-H(2) matrix isolation IR spectroscopy [M. Bahou et al., J. Chem. Phys. 136, 154304 (2012)]. The energetically most favourable conformer of C(6)H(7)(+)·N(2) shows a π-bonded structure similar to C(6)H(7)(+)·Rg (Rg = Ne, Ar) [P. Botschwina and R. Oswald, J. Phys. Chem. A 115, 13664 (2011)] with D(e) ≈ 870 cm(-1). For C(6)H(7)(+)·CO(2), a slightly lower energy is calculated for a conformer with the CO(2) ligand lying in the ring-plane of the C(6)H(7)(+) moiety (D(e) ≈ 1508 cm(-1)). It may be discriminated from other conformers through a strong band predicted at 1218 cm(-1), red-shifted by 21 cm(-1) from the corresponding band of free C(6)H(7)(+).     

Communication: determination of the bond dissociation energy (D0) of the water dimer, (H2O)2, by velocity map imaging

The bond dissociation energy (D(0)) of the water dimer is determined by using state-to-state vibrational predissociation measurements following excitation of the bound OH stretch fundamental of the donor unit of the dimer. Velocity map imaging and resonance-enhanced multiphoton ionization (REMPI) are used to determine pair-correlated product velocity and translational energy distributions. H(2)O fragments are detected in the ground vibrational (000) and the first excited bending (010) states by 2 + 1 REMPI via the C? (1)B(1) (000) ← X? (1)A(1) (000 and 010) transitions. The fragments' velocity and center-of-mass translational energy distributions are determined from images of selected rovibrational levels of H(2)O. An accurate value for D(0) is obtained by fitting both the structure in the images and the maximum velocity of the fragments. This value, D(0) = 1105 ± 10 cm(-1) (13.2 ± 0.12 kJ/mol), is in excellent agreement with the recent theoretical value of D(0) = 1103 ± 4 cm(-1) (13.2 ± 0.05 kJ∕mol) suggested as a benchmark by Shank et al. [J. Chem. Phys. 130, 144314 (2009)].     

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.     

A short hydrogen bond investigation by polarized Raman spectra of Co2+ and Zn2+ salts of pyromellitic acid

Cobalt and zinc salts of 1,2,4,5-benzenetetracarboxylic acid (pyromellitic acid), [C(6)H(2)(COO)(4)H(4)], have been synthesized and investigate by polarized Raman spectroscopy. These compounds present short intramolecular hydrogen bonds (SHB) between adjacent carboxyl groups. Raman spectra indicate the presence of this interaction in both salts. Three specific vibrational of SHB modes have been investigated: O-H-O symmetric [nu(sym)(OHO)] and asymmetric [nu(asym)(OHO)] stretching modes and O-H stretching mode [nu(O-H)], which they were observed around 300, 850 and 2500 cm(-1), respectively. In crystallographic point of view, the cobalt salt presents a symmetric SHB while the zinc salt presents an asymmetric SHB. In cobalt salt all three vibrational modes of O-H-O groups in polarized Raman spectra occur in A(g) orientation although in zinc salts two of them are observed in A(g) orientation and one in B(g). Spectra analysis indicate that nu(sym)(OHO) mode is observed as A(g) to cobalt salt and B(g) to zinc salt. This mode occurs in a crowded spectral region and its identification was made by deconvolution techniques. Comparing spectra of the two salts, it is observed a small difference in relative intensity and wavenumber shift of nu(sym)(OHO) (deviance of 43 cm(-1)) and nu(OH) (deviance of 21 cm(-1)) modes due probably to differences in O...O distance between salts and in orientation of pyromellitate anion in unit cell. The nu(asym)(OHO) mode does not present significant wavenumber shift due difference in SHB. The nu(OH) band presents a great potential for hydrogen bond studies due to the fact that in its vibrational region (around 2500 cm(-1)) it is not observed other vibrational modes of these compounds.     

Structure of the aqueous solvated electron from resonance Raman spectroscopy: lessons from isotopic mixtures

The structure and thermodynamics of the hydrated electron are probed with resonance Raman spectroscopy of isotopic mixtures of H(2)O and D(2)O. The strongly enhanced intramolecular bends of e(-)(H(2)O) and e(-)(D(2)O) produce single downshifted bands, whereas the e(-)(HOD) bend consists of two components: one slightly upshifted from the 1,446 cm(-1) bulk frequency to 1,457 cm(-1) and the other strongly downshifted to approximately 1,396 cm(-1). This 60 cm(-1) split and the 200 (120) cm(-1) downshifts of the OH (OD) stretch frequencies relative to bulk water reveal that the water molecules that are Franck-Condon coupled to the electron are in an asymmetric environment, with one proton forming a strong hydrogen bond to the electron. The downshifted bend and librational frequencies also indicate significantly weakened torsional restoring forces on the water molecules of e(-)(aq), which suggests that the outlying proton is a poor hydrogen bond donor to the surrounding solvent. A 1.6-fold thermodynamic preference of the electron for H(2)O is observed based on the relative intensities of the e(-)(H(2)O) and e(-)(D(2)O) bands in a 50:50 isotopic mixture. This equilibrium isotope effect is consistent with the downshifted vibrational frequencies and a relative reduction of the zero-point energy of H(2)O bound to the electron. Our results enhance the cavity model of the solvated electron and support only those models that contain water monomers as opposed to other molecular species.     

The method of Gaussian weighted trajectories. V. On the 1GB procedure for polyatomic processes

In recent years, many chemical reactions have been studied by means of the quasiclassical trajectory (QCT) method within the Gaussian binning (GB) procedure. The latter consists of "quantizing" the final vibrational actions in Bohr spirit by putting strong emphasis on the trajectories reaching the products with vibrational actions close to integer values. A major drawback of this procedure is that if N is the number of product vibrational modes, the amount of trajectories necessary to converge the calculations is ～10(N)×larger than with the standard QCT method. Applying it to polyatomic processes is thus problematic. In a recent paper, however, Czako? and Bowman propose to quantize the total vibrational energy instead of the vibrational actions [G. Czako? and J. M. Bowman, J. Chem. Phys. 131, 244302 (2009)], a procedure called 1GB here. The calculations are then only ～10 times more time consuming than with the standard QCT method, allowing thereby for considerable numerical saving. In this paper, we propose some theoretical arguments supporting the 1GB procedure and check its validity on model test cases as well as the prototype four-atom reaction OH+D(2)→HOD+D.     

Conformational dependence of intramolecular vibrational redistribution in methanol

Previous state-selected spectra of methanol in the 5nu(1) OH stretch overtone region [O. V. Boyarkin, T. R. Rizzo, and D. S. Perry, J. Chem. Phys. 110, 11346 (1999)] revealed a structure indicating an intramolecular vibrational redistribution on three time scales. Whereas in that work, methanol in the 5nu(1) bright state was prepared close to the staggered conformation, methanol in the "partially eclipsed" conformation is prepared here by double resonance excitation through a torsionally excited intermediate state. The excited molecules are detected by infrared laser assisted photofragment spectroscopy. In partially eclipsed methanol, the strong coupling of the nu(1) OH stretch to the nu(2) CH stretch becomes weaker, but the coupling responsible for the widths of the narrowest features becomes stronger.     

Structure and dynamics of water dangling OH bonds in hydrophobic hydration shells. Comparison of simulation and experiment

Molecular dynamics and electric field strength simulations are performed in order to quantify the structural, dynamic, and vibrational properties of non-H-bonded (dangling) OH groups in the hydration shell of neopentane, as well as in bulk water. The results are found to be in good agreement with the experimentally observed high-frequency (～3660 cm(-1)) OH band arising from the hydration shell of neopentanol dissolved in HOD/D(2)O, obtained by analyzing variable concentration Raman spectra using multivariate curve resolution (Raman-MCR). The simulation results further indicate that hydration shell dangling OH groups preferentially point toward the central carbon atom of neopentane to a degree that increases with the lifetime of the dangling OH.     

Full-dimensional quantum dynamics study of exchange processes for the D + H2O and D + HOD reactions

The exchange processes of D + H(2)O and D + HOD reactions are studied using initial state-selected time-dependent wave packet approach in full dimension. The total reaction probabilities for different partial waves, together with the integral cross sections, are obtained both by the centrifugal sudden (CS) approximation and exact coupled-channel (CC) calculations, for the H(2)O(HOD) reactant initially in the ground rovibrational state. In the CC calculations, small resonance peaks in the reaction probabilities and quick diminishing of the resonance peaks with the increase of total angular momenta J do not lead to clear step-like features just above the threshold in the cross sections for the title reactions, which are different in other isotopically substituted reactions where the hydrogen atom was included as the reactant instead of the deuterium atom [B. Fu, Y. Zhou, and D. H. Zhang, Chem. Sci. 3, 270 (2012); B. Fu and D. H. Zhang, J. Phys. Chem. A 116, 820 (2012)]. It is interesting that the shape resonance-induced features resulting from the reaction tunneling are significantly diminished accordingly in the reactions of the deuterium atom and H(2)O or HOD, owing to the weaker tunneling capability of the reagent deuterium atom in the title reactions than the reagent hydrogen atom in other reactions. In the CS calculations, the resonance peaks persist in many partial waves but cannot survive the partial-wave summations. The cross sections for the D(') + H(2)O → D(')OH + H and D(') + HOD → D(')OD + H reactions are substantially larger than those for the D(') + HOD → HOD(') + D reaction, indicating that the D(')/H exchange reactions are much more favored than the D(')/D exchange.     

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.     

Communication: quasiclassical trajectory calculations of correlated product-state distributions for the dissociation of (H2O)2 and (D2O)2

Stimulated by recent experiments [B. E. Rocher-Casterline, L. C. Ch'ng, A. K. Mollner, and H. Reisler, J. Chem. Phys. 134, 211101 (2011)], we report quasiclassical trajectory calculations of the dissociation dynamics of the water dimer, (H(2)O)(2) (and also (D(2)O)(2)) using a full-dimensional ab initio potential energy surface. The dissociation is initiated by exciting the H-bonded OH(OD)-stretch, as done experimentally for (H(2)O)(2). Normal mode analysis of the fragment pairs is done and the correlated vibrational populations are obtained by (a) standard histogram binning (HB), (b) harmonic normal-mode energy-based Gaussian binning (GB), and (c) a modified version of (b) using accurate vibrational energies obtained in the Cartesian space. We show that HB allows opening quantum mechanically closed states, whereas GB, especially via (c), gives physically correct results. Dissociation of both (H(2)O)(2) and (D(2)O)(2) mainly produces either fragment in the bending excited (010) state. The H(2)O(J) and D(2)O(J) rotational distributions are similar, peaking at J = 3-5. The computations do not show significant difference between the ro-vibrational distributions of the donor and acceptor fragments. Diffusion Monte Carlo computations are performed for (D(2)O)(2) providing an accurate zero-point energy of 7247 cm(-1), and thus, a benchmark D(0) of 1244 ± 5 cm(-1).