State-resolved infrared spectrum of the protonated water dimer: revisiting the characteristic proton transfer doublet peak

The infrared (IR) spectra of protonated water clusters encode precise information on the dynamics and structure of the hydrated proton. However, the strong anharmonic coupling and quantum effects of these elusive species remain puzzling up to the present day. Here, we report unequivocal evidence that the interplay between the proton transfer and the water wagging motions in the protonated water dimer (Zundel ion) giving rise to the characteristic doublet peak is both more complex and more sensitive to subtle energetic changes than previously thought. In particular, hitherto overlooked low-intensity satellite peaks in the experimental spectrum are now unveiled and mechanistically assigned. Our findings rely on the comparison of IR spectra obtained using two highly accurate potential energy surfaces in conjunction with highly accurate state-resolved quantum simulations. We demonstrate that these high-accuracy simulations are important for providing definite assignments of the complex IR signals of fluxional molecules.

We reveal the intricate dynamics of the proton shuttling motion in the Zundel ion by computing 900 high-accuracy vibrational eigenstates. We show how very subtle energetic changes in the vibrational modes lead to vastly different infrared spectra.     

Eigen and Zundel forms of small protonated water clusters: structures and infrared spectra

The spectral properties of protonated water clusters, especially the difference between Eigen (H3O+) and Zundel (H5O2+) conformers and the difference between their unhydrated and dominant hydrated forms are investigated with the first principles molecular dynamics simulations as well as with the high level ab initio calculations. The vibrational modes of the excess proton in H3O+ are sensitive to the hydration, while those in H5O2+ are sensitive to the messenger atom such as Ar (which was assumed to be weakly bound to the water cluster during acquisitions of experimental spectra). The spectral feature around approximately 2700 cm-1 (experimental value: 2665 cm-1) for the Eigen moiety appears when H3O+ is hydrated. This feature corresponds to the hydrating water interacting with H3O+, so it cannot appear in the Eigen core. Thus, H3O+ alone would be somewhat different from the Eigen forms in water. For the Zundel form (in particular, H5O2+), there have been some differences in spectral features among different experiments as well as between experiments and theory. When an Ar messenger atom is introduced at a specific temperature corresponding to the experimental condition, the calculated vibrational spectra for H5O2+.Ar are in good agreement with the experimental infrared spectra showing the characteristic Zundel frequency at approximately 1770 cm-1. Thus, the effect of hydration, messenger atom Ar, and temperature are crucial to elucidating the nature of vibrational spectra of Eigen and Zundel forms and to assigning the vibrational modes of small protonated water clusters.     

Gas phase infrared multiple-photon dissociation spectra of methanol, ethanol and propanol proton-bound dimers, protonated propanol and the propanol/water proton-bound dimer

The infrared multiphoton dissociation (IRMPD) spectra of three homogenous proton-bound dimers are presented and the major features are assigned based on comparisons with the neutral alcohol and with density functional theory calculations. As well, the IRMPD spectra of protonated propanol and the propanol/water proton-bound dimer (or singly hydrated protonated propanol) are presented and analysed. Two primary IRMPD photoproducts were observed for each of the alcohol proton bound dimers and were found to vary with the frequency of the radiation impinging upon the ions. For example, when the proton-bound dimer absorbs weakly a larger amount of S(N)2 product, protonated ether and water, are observed. When the proton-bound dimer absorbs more strongly, an increase in the simple dissociation product, protonated alcohol and neutral alcohol, is observed. With the aid of RRKM calculations this frequency dependence of the branching ratio is explained by assuming that photon absorption is faster than dissociation for these species and that only a few photons extra are necessary to make the higher-energy dissociation channel (simple cleavage) competitive with the lower energy (S(N)2) reaction channel.     

IR spectrum and structure of protonated ethanol dimer: implications for the mobility of excess protons in solution

This Communication reports IR spectra and density functional calculations for the isolated protonated ethanol dimer and its N2-microsolvated complexes, (EtOH)2H+-(N2)n (n = 0-2) to investigate the degree of delocalization of the excess proton in this fundamental building block of an alcohol proton wire. The first spectroscopic characterization of isolated and microsolvated (EtOH)2H+ suggests that the excess proton is (nearly) equally shared between both EtOH units under symmetric solvation conditions (Zundel-type ion, n = 0 and 2), whereas it is largely localized on a single EtOH molecule for asymmetric solvation (Eigen-type ion, n = 1).     

Calculation of vibrational transition frequencies and intensities in water dimer: comparison of different vibrational approaches

We have calculated frequencies and intensities of fundamental and overtone vibrational transitions in water and water dimer with use of different vibrational methods. We have compared results obtained with correlation-corrected vibrational self-consistent-field theory and vibrational second-order perturbation theory both using normal modes and finally with a harmonically coupled anharmonic oscillator local mode model including OH-stretching and HOH-bending local modes. The coupled cluster with singles, doubles, and perturbative triples ab initio method with augmented correlation-consistent triple-zeta Dunning and atomic natural orbital basis sets has been used to obtain the necessary potential energy and dipole moment surfaces. We identify the strengths and weaknesses of these different vibrational approaches and compare our results to the available experimental results.     

Quantum chemical calculation of infrared spectra of acidic groups in chabazite in the presence of water

The changes in the spectra of the acidic group in chabazite are studied by quantum chemical calculations. The zeolite is modeled by two clusters consisting of eight tetrahedral atoms arranged in a ring and seven tetrahedral atoms coordinated around the zeolite OH group. The potential energy and dipole surfaces were constructed from the zeolite OH stretch, in-plane and out-of-plane bending coordinates, and the intermolecular stretch coordinate that corresponds to a movement of the water molecule as a whole. Both the anharmonicities of the potential energy and dipole were taken into account by calculation of the frequencies and intensities. The matrix elements of the vibrational Hamiltonian were calculated within the discrete variable representation basis set. We have assigned the experimentally observed frequencies at approximately 2900, approximately 2400, and approximately 1700 cm(-1) to the strongly perturbed zeolite OH vibrations caused by the hydrogen bonding with the water molecule. The ABC triplet is a Fermi resonance of the zeolite OH stretch mode with the overtone of the in-plane bending (the A band) and the overtone of the out-of-plane bending (the C band). In the B band the stretch is also coupled with the second overtone of the out-of-plane bending. The frequencies at approximately 3700 and approximately 3550 cm(-1) we have assigned to the OH stretch frequencies of a slightly perturbed water molecule.     

The structure of protonated acetone and its dimer: infrared photodissociation spectroscopy from 800 to 4,000 cm-1

The infrared spectra of protonated acetone and the proton bound acetone dimer are obtained revealing vibrational resonances associated with the shared proton motions, which are in agreement with the predictions from ab initio, MP2, harmonic frequency calculations.     

Infrared multiphoton dissociation spectra as a probe of ion molecule reaction mechanism: the formation of the protonated water dimer via sequential bimolecular reactions with 1,1,3,3-tetrafluorodimethyl ether

The gas-phase ion-molecule reactions of 1,1,3,3-tetrafluorodimethyl ether and water have been examined using Fourier transform ion cyclotron resonance mass spectrometry, infrared multiphoton dissociation (IRMPD) spectroscopy, and ab initio molecular orbital calculations. This reaction sequence leads to the efficient bimolecular production of the proton-bound dimer of water (H5O2+). Evidence for the dominant mechanistic pathway involving the reaction of CF2H-O=CHF+, an ion of m/z 99, with water is presented. The primary channel occurs via nucleophilic attack of water on the ion of m/z 99 (CF2H-O=CHF+), to lose formyl fluoride and yield-protonated difluoromethanol (m/z 69). Association of a second water molecule with protonated difluoromethanol generates a reactive intermediate that decomposes via a 1,4-elimination to release hydrogen fluoride and yield the proton-bound dimer of water and formyl fluoride (m/z 67). Last, the elimination of formyl fluoride occurs by the association of a third water molecule to produce H5O2+ (m/z 37). The most probable isomeric forms of the ions with m/z 99 and 69 were found using IRMPD spectroscopy and electronic structure theory calculations. Thermochemical information for reactant, transition state, and product species was obtained using MP2(full)/6-311+G**//6-31G* level of theory.     

Contribution of water dimer absorption to the millimeter and far infrared atmospheric water continuum

We present a rigorous calculation of the contribution of water dimers to the absorption coefficient alpha(nu,T) in the millimeter and far infrared domains, over a wide range (276-310 K) of temperatures. This calculation relies on the explicit consideration of all possible transitions within the entire rovibrational bound state manifold of the dimer. The water dimer is described by the flexible 12-dimensional potential energy surface previously fitted to far IR transitions [C. Leforestier et al., J. Chem. Phys. 117, 8710 (2002)], and which was recently further validated by the good agreement obtained for the calculated equilibrium constant Kp(T) with experimental data [Y. Scribano et al., J. Phys. Chem. A. 110, 5411 (2006)]. Transition dipole matrix elements were computed between all rovibrational states up to an excitation energy of 750 cm(-1), and J=K=5 rotational quantum numbers. It was shown by explicit calculations that these matrix elements could be extrapolated to much higher J values (J=30). Transitions to vibrational states located higher in energy were obtained from interpolation of computed matrix elements between a set of initial states spanning the 0-750 cm(-1) range and all vibrational states up to the dissociation limit (approximately 1200 cm(-1)). We compare our calculations with available experimental measurements of the water continuum absorption in the considered range. It appears that water dimers account for an important fraction of the observed continuum absorption in the millimeter region (0-10 cm(-1)). As frequency increases, their relative contribution decreases, becoming small (approximately 3%) at the highest frequency considered nu=944 cm(-1).     

Matrix isolation infrared spectra of the complexes between methylacetate and water or hydrochloric acid

The hydrogen-bonded complexes between methylacetate (MeAc) and water or hydrochloric acid have been studied by infrared spectrometry in a low temperature Argon matrix. The ν_CO, ν_CC and ν_CC ional modes of MeAc show a splitting to the low and to the high frequency side of the free molecule. This suggests that water and HCl interact with the keto and ether oxygens. This conclusions is supported by the appearance of two main absorptions in the ν_HCl region. These results are discussed as a function of the gas-phase basicity of the two oxygen atoms, derived from the O(1s) core electron binding energies and from the ionization potentials.     

Density functional potential energy hypersurface of protonated ozone: A comparison between different gradient-corrected nonlocal functionals

The more stable isomers obtained by protonation of both cyclic and open-chain forms of ozone (O₃) were studied by the linear combination of Gaussian-type orbitals-density functional (LCGTO-DF) method. Features of the first-order saddle points connecting them were also investigated. Nonlocal corrections to the exchange-correlation energy were added using different kinds of functionals. The calculated proton affinity (PA) values at 298 K compare favorably with the recent experimental determination. © 1997 John Wiley & Sons, Inc.     

Structure, stability, thermodynamic properties, and infrared spectra of the protonated water octamer H(+)(H2O)8

We investigated various two-dimensional (2D) and three-dimensional (3D) structures of H (+)(H 2O) 8, using density functional theory (DFT), Moller-Plesset second-order perturbation theory (MP2), and coupled cluster theory with single, double, and perturbative triple excitations (CCSD(T)). The 3D structure is more stable than the 2D structure at all levels of theory on the Born-Oppenheimer surface. With the zero-point energy (ZPE) correction, the predicted structure varies depending on the level of theory. The DFT employing Becke's three parameters with Lee-Yang-Parr functionals (B3LYP) favors the 2D structure. At the complete basis set (CBS) limit, the MP2 calculation favors the 3D structure by 0.29 kcal/mol, and the CCSD(T) calculation favors the 3D structure by 0.27 kcal/mol. It is thus expected that both 2D and 3D structures are nearly isoenergetic near 0 K. At 100 K, all the calculations show that the 2D structure is much more stable in free binding energy than the 3D structure. The DFT and MP2 vibrational spectra of the 2D structure are consistent with the experimental spectra. First-principles Car-Parrinello molecular dynamics (CPMD) simulations show that the 2D Zundel-type vibrational spectra are in good agreement with the experiment.     

A sodium atom in a large water cluster: electron delocalization and infrared spectra

Ab initio molecular dynamics simulations modeling low-energy collisions of a sodium atom with a cluster with more than 30 water molecules are presented. We follow the dynamics of the atom-cluster interaction and the delocalization of the valence electron of sodium together with the changes in the electron binding energy. This electron tends to be shared by the nascent sodium cation and the water cluster. IR spectra of the sodium-water cluster are both computationally and experimentally obtained, with a good agreement between the two approaches.     

Infrared spectra of the protonated neurotransmitter histamine: competition between imidazolium and ammonium isomers in the gas phase

The infrared (IR) spectrum of protonated histamine (histamineH(+)) was recorded in the 575-1900 cm(-1) fingerprint range by means of IR multiple photon dissociation (IRMPD) spectroscopy. The IRMPD spectrum of mass-selected histamineH(+) ions was obtained in a Fourier transform ion cyclotron resonance mass spectrometer coupled to an electrospray ionization source and an IR free electron laser. A variety of isomers were identified and characterized by quantum chemical calculations at the B3LYP and MP2 levels of theory using the cc-pVDZ basis set. The low-energy isomers are derived from various favourable protonation sites--all of which are N atoms--and different orientations of the ethylamine side chain with respect to the heterocyclic imidazole ring. The measured IRMPD spectrum was monitored in the NH(3) loss channel and exhibits 14 bands in the investigated spectral range, which were assigned to vibrational transitions of the most stable isomer, denoted A. This imidazolium-type isomer A with protonation at the imidazole ring and gauche conformation of the ethylamine side chain is significantly stabilized by an intramolecular ionic Nπ-H(+)···Nα hydrogen bond to the ethylamino group. The slightly less stable ammonium-type isomer B with protonation at the ethylamino group is only a few kJ mol(-1) higher in energy and may also provide a minor contribution to the observed IRMPD spectrum. Isomer B is derived from A by simple proton transfer from imidazole to the ethylamino group along the intramolecular Nπ-H(+)···Nα hydrogen bond via a low barrier, which is calculated to be of the order of 5-15 kJ mol(-1). Significantly, the most stable structure of isolated histamineH(+) differs from that in the condensed phase by both the protonation site and the conformation of the side chain, emphasizing the important effects of solvation on the structure and function of this neurotransmitter. The effects of protonation on the geometric and electronic structure of histamine are evaluated by comparing the calculated properties of isomer A with those of the most stable structure of neutral histamine A(n).     

Effect of substituents on the chemical shift of the ring protons in the PMR spectra of monosubstituted sym-triazines

A correlation equation that links the relative chemical shift of the protons of the triazine ring with the F and R constants of the substituents was obtained on the basis of data from the PMR spectra of solutions of sym-triazine and its monosubstituted derivatives [OCH₃, N(CH₃)₂, CH₃, C₆H₅, COOC₂H₅, and CN] in dimethyl sulfoxide (DMSO). The equation was analyzed by comparison with the corresponding equations for monosubstituted benzenes and pyrimidines.Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 3, pp. 410–414, March, 1982.