Vibrationally Induced Proton Transfer in F−(H2O) and F−(D2O)

Evaluation of Ar tagging toward the vibrational spectra and zero point energy of X-HOH, X-DOH, and X-HOD, for X = F, Cl, Br

In this study, we theoretically evaluated the effect of argon tagging toward the binding energy and vibrational spectra of water halide anion complexes Ar.X-HOH, Ar.X-HOD, and Ar.X-DOH (X = F, Cl, Br). The ionic hydrogen bond (IHB) OH stretching mode was calculated to have a strong peak in the vibrational spectra, and coupling to intermolecular modes as well as bending modes was observed as combination bands and Fermi resonances. We found that the argon tagging affected the IHB OH stretching peak position in Ar.F-H2O, but not in Ar.Cl-H2O and Ar.Br-H2O. Furthermore, D-binding is favored for Cl and Br based on zero point energies, but for F our calculated zero point energies did not show a preference between H- and D-binding. We show that the competition of the energy lowering in the zero point energy of the anharmonic IHB OH (OD) stretching mode versus the intermolecular out-of-plane IHB OH (OD) wagging mode is important for determining the preference between H- and D-binding for these monohydrated halide clusters. We also found that for X-HOD the HOD bending fundamental peak is blue shifted compared to bare HOD, but is redshifted for F-DOH.

Electrostatics, Charge Transfer, and the Nature of the Halide-Water Hydrogen Bond

Binary halide-water complexes X^-(H₂O) are examined by means of symmetry-adapted perturbation theory, using charge-constrained promolecular reference densities to extract a meaningful charge-transfer component from the induction energy. As is known, the X^-(H₂O) potential energy surface (for X = F, Cl, Br, or I) is characterized by symmetric left and right hydrogen bonds separated by a C_2v-symmetric saddle point, with a tunneling barrier height that is <2 kcal/mol except in the case of F^-(H₂O). Our analysis demonstrates that the charge-transfer energy is correspondingly small (<2 kcal/mol except for X = F), considerably smaller than the electrostatic interaction energy. Nevertheless, charge transfer plays a crucial role determining the conformational preferences of X^-(H₂O) and provides a driving force for the formation of quasi-linear X··· H-O hydrogen bonds. Charge-transfer energies correlate well with measured O-H vibrational redshifts for the halide-water complexes and also for OH^-(H₂O) and NO₂^-(H₂O), providing some indication of a general mechanism.

Water Radical Cations in the Gas Phase: Methods and Mechanisms of Formation, Structure and Chemical Properties

Water radical cations, (H₂O)_n^+•, are of great research interest in both fundamental and applied sciences. Fundamental studies of water radical reactions are important to better understand the mechanisms of natural processes, such as proton transfer in aqueous solutions, the formation of hydrogen bonds and DNA damage, as well as for the discovery of new gas-phase reactions and products. In applied science, the interest in water radicals is prompted by their potential in radiobiology and as a source of primary ions for selective and sensitive chemical ionization. However, in contrast to protonated water clusters, (H₂O)_nH⁺, which are relatively easy to generate and isolate in experiments, the generation and isolation of radical water clusters, (H₂O)_n^+•, is tremendously difficult due to their ultra-high reactivity. This review focuses on the current knowledge and unknowns regarding (H₂O)_n^+• species, including the methods and mechanisms of their formation, structure and chemical properties.

Near-Infrared Spectroscopy and Anharmonic Theory of Protonated Water Clusters: Higher Elevations in the Hydrogen Bonding Landscape

Near-infrared spectroscopy measurements are presented for protonated water clusters, H⁺(H₂O) _n, in the size range of n = 1-8. Clusters are produced in a pulsed-discharge supersonic expansion, mass selected, and studied with infrared laser photodissociation spectroscopy in the regions of 3600-4550 and 4850-7350 cm^-1. Although there is some variation with cluster size, the main features of these spectra are a broad absorption near 5300 cm^-1, a sharp doublet near 7200 cm^-1, as well as a structured absorption near 4100 cm^-1 for n ≥ 2. The vibrational patterns measured for the hydronium, Zundel, and Eigen ions are compared to those predicted by different forms of anharmonic theory. Second-order vibrational perturbation theory (VPT2) and a local mode treatment of the OH stretches both capture key aspects of the spectra but suffer understandable deficiencies in the quantitative description of band positions and intensities.

Vibrational Characterization of Radical Ion Adducts between Imidazole and CO2

We address the molecular level origins of the dramatic difference in the catalytic mechanisms of CO₂ activation by the seemingly similar molecules pyridine (Py) and imidazole (Im). This is accomplished by comparing the fundamental interactions of CO₂ radical anions with Py and Im in the isolated, gas phase PyCO₂^- and ImCO₂^- complexes. These species are prepared by condensation of the neutral compounds onto a (CO₂) _n^- cluster ion beam by entrainment in a supersonic jet ion source. The structures of the anionic complexes are determined by theoretical analysis of their vibrational spectra, obtained by IR photodissociation of weakly bound CO₂ molecules in a photofragmentation mass spectrometer. Although the radical PyCO₂^- system adopts a carbamate-like configuration corresponding to formation of an N-C covalent bond, the ImCO₂^- species is revealed to be best described as an ion-molecule complex in which an oxygen atom in the CO₂^- radical anion is H-bonded to the NH group. Species that feature a covalent N-C interaction in ImCO₂^- are calculated to be locally stable structures, but are much higher in energy than the largely electrostatically bound ion-molecule complex. These results support the suggestion from solution phase electrochemical studies (Bocarsly et al. ACS Catal. 2012, 2, 1684-1692) that the N atoms are not directly involved in the catalytic activation of CO₂ by Im.

Structures and Spectroscopic Properties of F-(H2O) n with n = 1-10 Clusters from a Global Search Based On Density Functional Theory

Using a genetic algorithm incorporated in density functional theory, we explore the ground state structures of fluoride anion-water clusters F^-(H₂O) _n with n = 1-10. The F^-(H₂O) _n clusters prefer structures in which the F^- anion remains at the surface of the structure and coordinates with four water molecules, as the F^-(H₂O) _n clusters have strong F^--H₂O interactions as well as strong hydrogen bonds between H₂O molecules. The strong interaction between the F^- anion and adjacent H₂O molecule leads to a longer O-H distance in the adjacent molecule than in an individual water molecule. The simulated infrared (IR) spectra of the F^-(H₂O)_1-5 clusters obtained via second-order vibrational perturbation theory (VPT2) and including anharmonic effects reproduce the experimental results quite well. The strong interaction between the F^- anion and water molecules results in a large redshift (600-2300 cm^-1) of the adjacent O-H stretching mode. Natural bond orbital (NBO) analysis of the lowest-energy structures of the F^-(H₂O)_1-10 clusters illustrates that charge transfer from the lone pair electron orbital of F^- to the antibonding orbital of the adjacent O-H is mainly responsible for the strong interaction between the F^- anion and water molecules, which leads to distinctly different geometric and vibrational properties compared with neutral water clusters.

Effects of vibrational excitation on the F + H2O → HF + OH reaction: dissociative photodetachment of overtone-excited [F-H-OH]. Chem Sci 2017;8:7821-7833. [PMID: 29163919 PMCID: PMC5674243 DOI: 10.1039/c7sc03364h]

Photodetaching vibrationally excited FH₂O^– channels energy into the reaction coordinate of the F + H₂O reaction, as shown in this joint experimental-theoretical study.

The reaction F + H₂O → HF + OH is a four-atom system that provides an important benchmark for reaction dynamics. Hydrogen atom transfer at the transition state for this reaction is expected to exhibit a strong dependence on reactant vibrational excitation. In the present study, the vibrational effects are examined by photodetachment of vibrationally excited F^–(H₂O) precursor anions using photoelectron-photofragment coincidence (PPC) spectroscopy and compared with full six-dimensional quantum dynamical calculations on ab initio potential energy surfaces. Prior to photodetachment at hν_UV = 4.80 eV, the overtone of the ionic hydrogen bond mode in the precursor F^–(H₂O), 2ν_IHB at 2885 cm^–1, was excited using a tunable IR laser. Experiment and theory show that vibrational energy in the anion can be effectively carried away by the photoelectron upon a Franck–Condon photodetachment, and also show evidence for an increase of branching into the F + H₂O reactant channel. The experimental results suggest a greater role for product rotational excitation than theory. Improved potential energy surfaces and longer wavepacket propagation times would be helpful to further examine the nature of the discrepancy.

Disentangling the Complex Vibrational Spectrum of the Protonated Water Trimer, H+(H2O)3, with Two-Color IR-IR Photodissociation of the Bare Ion and Anharmonic VSCF/VCI Theory

Vibrational spectroscopy of the protonated water trimer provides a stringent constraint on the details of the potential energy surface (PES) and vibrational dynamics governing excess proton motion far from equilibrium. Here we report the linear spectrum of the cold, bare H⁺(H₂O)₃ ion using a two-color, IR-IR photofragmentation technique and follow the evolution of the bands with increasing ion trap temperature. The key low-energy features are insensitive to both D₂ tagging and internal energy. The D₂-tagged D⁺(D₂O)₃ spectrum is reported for the first time, and the isotope dependence of the band pattern is surprisingly complex. These spectra are reproduced by large-scale vibrational configuration interaction calculations based on a new full-dimensional PES, which treat the anharmonic effects arising from large amplitude motion. The results indicate such extensive mode mixing in both isotopologues that one should be cautious about assigning even the strongest features to particular motions, especially for the absorptions that occur close to the intramolecular bending mode of the water molecule.

Theoretical calculation of the vibrational state dependent photodetachment spectra of X-H2O, X = F, Cl, Br

Vibrational spectra of X^-H₂O (X = F, Cl, Br) were simulated using full dimensional vibrational calculations using quantum chemistry potential energy surfaces. Furthermore, utilizing the reflection approximation, we simulated the photodetachment spectra obtained from different vibrational excited states. From these spectra, we can observe changes in the hydrogen bond interaction between the anion and the neutral XH₂O system. Notably, for F^-H₂O, the excitation of the ionic hydrogen bonded (IHB) OH stretching vibration generates a large tail on the low energy side of the photodetachment spectra compared to the detachment from the zero-point vibration state. This shows that the IHB OH stretching vibration of F^-H₂O causes charge delocalization from F^- to the oxygen atom in H₂O, and that the photodetachment from FHOH^- occurs at lower energies.

An investigation on the structure, spectroscopy and thermodynamic aspects of Br2((-))(H2O)n clusters using a conjunction of stochastic and quantum chemical methods

In this work we obtained global as well as local structures of Br2((-))(H2O)n clusters for n = 2 to 6 followed by the study of IR-spectral features and thermochemistry for the structures. The way adopted by us to obtain structures is not the conventional one used in most cases. Here we at first generated excellent quality pre-optimized structures by exploring the suitable empirical potential energy surface using stochastic optimizer simulated annealing. These structures are then further refined using quantum chemical calculations to obtain the final structures, and spectral and thermodynamic features. We clearly showed that our approach results in very quick and better convergence which reduces the computational cost and obviously using the strategy we are able to get one [i.e. global] or more than one [i.e. global and local(s)] energetically lower structures than those which are already reported for a given cluster size. Moreover, IR-spectral results and the evolutionary trends in interaction energy, solvation energy and vertical detachment energy for global structures of each size have also been presented to establish the utility of the procedure employed.

Rovibrational energy levels of the F−(H2O) and F−(D2O) complexes

A new accurate 6D PES is determined obtained from CCSD(T)-F12 calculations including two dissociation channels (HF + OH⁻ and F⁻ + H₂O). A novel way is developed to use complex coordinates in variational nuclear motion computations. The rovibrational energies of F⁻(H₂O) (the complete set up to 3700 cm⁻¹) and F⁻(D₂O) have been computed. The tunneling splittings describing the two complexes are obtained.

Structural and spectroscopic studies of iodine dimer radical anion hydrated clusters: an approach using a combination of stochastic and quantum chemical methods

Structures of I₂⁽⁻⁾(H₂O)₅ clusters after evaluation by simulated annealing and subsequent DFT calculation respectively.

O-H anharmonic vibrational motions in Cl(-)···(CH3OH)(1-2) ionic clusters. Combined IRPD experiments and AIMD simulations

The structures of Cl(-)-(Methanol)1,2 clusters have been unraveled combining Infrared Predissociation (IR-PD) experiments and DFT-based molecular dynamics simulations (DFT-MD) at 100 K. The dynamical IR spectra extracted from DFT-MD provide the initial 600 cm(-1) large anharmonic red-shift of the O-H stretch from uncomplexed methanol (3682 cm(-1)) to Cl(-)-(Methanol)1 complex (3085 cm(-1)) as observed in the IR-PD experiment, as well as the subtle supplementary blue- and red-shifts of the O-H stretch in Cl(-)-(Methanol)2 depending on the structure. The anharmonic vibrational calculations remarkably provide the 100 cm(-1) O-H blue-shift when the two methanol molecules are simultaneously organized in the anion first hydration shell (conformer 2A), while they provide the 240 cm(-1) O-H red-shift when the second methanol is in the second hydration shell of Cl(-) (conformer 2B). RRKM calculations have also shown that 2A/2B conformers interconvert on a nanosecond time-scale at the estimated 100 K temperature of the clusters formed by evaporative cooling of argon prior to the IR-PD process.

Ab initio potential energy and dipole moment surfaces of the F(-)(H2O) complex

We present full-dimensional, ab initio potential energy and dipole moment surfaces for the F(-)(H2O) complex. The potential surface is a permutationally invariant fit to 16,114 coupled-cluster single double (triple)/aVTZ energies, while the dipole surface is a covariant fit to 11,395 CCSD(T)/aVTZ dipole moments. Vibrational self-consistent field/vibrational configuration interaction (VSCF/VCI) calculations of energies and the IR-spectrum are presented both for F(-)(H2O) and for the deuterated analog, F(-)(D2O). A one-dimensional calculation of the splitting of the ground state, due to equivalent double-well global minima, is also reported.

Correlated ab initio investigations on the intermolecular and intramolecular potential energy surfaces in the ground electronic state of the O2(-)(X2Πg)-HF(X1Σ+) complex

This work reports the first highly correlated ab initio study of the intermolecular and intramolecular potential energy surfaces in the ground electronic state of the O(2)(-)(X(2)Π(g))-HF(X(1)Σ(+)) complex. Accurate electronic structure calculations were performed using the coupled cluster method including single and double excitations with addition of the perturbative triples correction [CCSD(T)] with the Dunning's correlation consistent basis sets aug-cc-pVnZ, n = 2-5. Also, the explicitly correlated CCSD(T)-F12a level of theory was employed with the AVnZ basis as well as the Peterson and co-workers VnZ-F12 basis sets with n = 2 and 3. Results of all levels of calculations predicted two equivalent minimum energy structures of planar geometry and C(s) symmetry along the A" surface of the complex, whereas the A' surface is repulsive. Values of the geometrical parameters and the counterpoise corrected dissociation energies (Cp-D(e)) that were calculated using the CCSD(T)-F12a/VnZ-F12 level of theory are in excellent agreement with those obtained from the CCSD(T)/aug-cc-pV5Z calculations. The minimum energy structure is characterized by a very short hydrogen bond of length of 1.328 Å, with elongation of the HF bond distance in the complex by 0.133 Å, and D(e) value of 32.313 Kcal/mol. Mulliken atomic charges showed that 65% of the negative charge is localized on the hydrogen bonded end of the superoxide radical and the HF unit becomes considerably polarized in the complex. These results suggest that the hydrogen bond is an incipient ionic bond. Exploration of the potential energy surface confirmed the identified minimum and provided support for vibrationally induced intramolecular proton transfer within the complex. The T-shaped geometry that possesses C(2v) symmetry presents a saddle point on the top of the barrier to the in-plane bending of the hydrogen above and below the axis that connects centers of masses of the monomers. The height of this barrier is 7.257 Kcal/mol, which is higher in energy than the hydrogen bending frequency by 909.2 cm(-1). The calculated harmonic oscillator vibrational frequencies showed that the H-F stretch vibrational transition in the complex is redshifted by 2564 cm(-1) and gained significant intensity (by at least a factor of 30) with respect to the transition in the HF monomer. These results make the O(2)(-)-HF complex an excellent prototype for infrared spectroscopic investigations on open-shell complexes with vibrationally induced proton transfer.

A local entropic signature of specific ion hydration

Monovalent ion hydration entropies are analyzed via energetic partitioning of the potential distribution theorem free energy. Extensive molecular dynamics simulations and free energy calculations are performed over a range of temperatures to determine the electrostatic and van der Waals components of the entropy. The far-field electrostatic contribution is negative and small in magnitude, and it does not vary significantly as a function of ion size, consistent with dielectric models. The local electrostatic contribution, however, varies widely as a function of ion size; the sign yields a direct indication of the kosmotropic (strongly hydrated) or chaotropic (weakly hydrated) nature of the ion hydration. The results provide a thermodynamic signature for specific ion effects in hydration and are consistent with experiments that suggest minimal perturbations of water structure outside the first hydration shell. The hydration entropies are also examined in relation to the corresponding entropies for the isoelectronic rare gas pairs; an inverse correlation is observed, as expected from thermodynamic hydration data.

Infrared spectra of SF6(-) x HCOOH x Ar(n) (n = 0-2): infrared triggered reaction and Ar-induced reactive inhibition

We present the infrared spectra of SF(6)(-) x HCOOH x Ar(m) (m=0-2) complexes. We find that the binding motif involves a single hydrogen bond between the SF(6)(-) anion and the OH group of the formic acid, with the CH group weakly tethered to a neighboring F atom. Similar to the case of hydrated SF(6)(-), the SF bond involved in the (OH-F) bond is significantly stretched and weakened by the attachment of the HCOOH ligand. The bare complex undergoes reaction upon infrared absorption in the CH/OH stretching region of the formic acid moiety, leading predominantly to the formation of SF(4)(-) + 2HF + CO(2). The reaction can be inhibited by attachment of two Ar atoms. We discuss a likely reaction mechanism in the framework of ab initio calculations, suggesting that reaction proceeds via tunneling through the potential barrier.