Ab initio potential energy and dipole moment surfaces for H5O2+

Experimental and modeling study of the ion-molecule association reaction H3O++H2O(+M)→H5O2+(+M)

Experimental results for the rate of the association reaction H3O+ + H2O (+M) --> H5O2(+) (+M) obtained with the Cinetique de Reactions en Ecoulement Supersonique Uniforme flow technique are reported. The reaction was studied in the bath gases M=He and N2, over the temperature range of 23-170 K, and at pressures between 0.16 and 3.1 mbar. At the highest temperatures, the reaction was found to be close to the limiting low-pressure termolecular range, whereas the limiting high-pressure bimolecular range was approached at the lowest temperatures. Whereas the low-pressure rate coefficients can satisfactorily be reproduced by standard unimolecular rate theory, the derived high-pressure rate coefficients in the bath gas He at the lowest temperatures are found to be markedly smaller than given by simple ion-dipole capture theory. This result differs from previous observations on the related reaction NH4(+) + NH3 (+M) --> N2H7(+) (+M). This observation is tentatively attributed to more pronounced contributions of the valence part of the potential-energy surface to the reaction in H5O2(+) than in N2H7(+). Falloff curves of the reaction H3O+ + H2O (+M) --> H5O2(+) (+M) are constructed over wide ranges of conditions and represented in compact analytical form.

Hydrogen-Bonded Networks in Ethanol Proton Wires: IR Spectra of (EtOH)qH+−Ln Clusters (L = Ar/N2, q ≤ 4, n ≤ 5)

Isolated and microsolvated protonated ethanol clusters, (EtOH)qH+-Ln with L = Ar and N2, are characterized by infrared photodissociation (IRPD) spectroscopy in the 3 microm range and quantum chemical calculations. For comparison, also the spectrum of the protonated methanol dimer, (MeOH)2H+, is presented. The IRPD spectra carry the signature of H-bonded (EtOH)qH+ chain structures, in which the excess proton is either strongly localized on one or (nearly) equally shared between two EtOH molecules, corresponding to Eigen-type ion cores (EtOH2+ for q = 1, 3) or Zundel-type ion cores (EtOH-H+-HOEt for q = 2, 4), respectively. In contrast to neutral (EtOH)q clusters, no cyclic (EtOH)qH+ isomers are detected in the size range investigated (q < or = 4), indicative of the substantial impact of the excess proton on the properties of the H-bonded ethanol network. The acidity of the two terminal OH groups in the (EtOH)qH+ chains decreases with the length of the chain (q). Comparison between (ROH)qH+ with R = CH3 and C2H5 shows that the acidity of the terminal O-H groups increases with the length of the aliphatic rest (R). The most stable (EtOH)qH+-Ln clusters with n < or = 2 feature intermolecular H-bonds between the inert ligands and the two available terminal OH groups of the (EtOH)qH+ chain. Asymmetric microsolvation of (EtOH)qH+ with q = 2 and 4 promotes a switch from Zundel-type to Eigen-type cores, demonstrating that the fundamental structural motif of the (EtOH)qH+ proton wire sensitively depends on the environment. The strength of the H-bonds between L and (EtOH)qH+ is shown to provide a rather sensitive probe of the acidity of the terminal OH groups.

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.

Ab initio global potential-energy surface for H5+→H3++H2

An accurate global potential-energy surface (PES) is reported for H5(+) based on more than 100,000 CCSD(T)/aug-cc-pVTZ ab initio energies. This PES has full permutational symmetry with respect to interchange of H atoms and dissociates to H3(+) and H2. Ten known stationary points of H5(+) are characterized and compared to previous ab initio calculations. Quantum diffusion Monte Carlo calculations are performed on the PES to obtain the zero-point energy of H5(+) and the anharmonic dissociation energy (D0) of H5(+) --> H3(+) + H2. The rigorous zero-point state of H4D+ is also calculated and discussed within the context of a strictly classical approach to obtain the branching ratio of the reaction H4D+ --> H3(+) + HD and H2D+ + H2. Such an approach is taken using the PES and critiqued based on the properties of the quantum zero-point state. Finally, a simple procedure for adding the long range-interaction energy is described.

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H(5)O(2) (+) using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh et al. [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.