Unusual Influence of Fluorinated Anions on the Stretching Vibrations of Liquid Water

Competition between Solvent···Solvent and Solvent···Solute Interactions in the Microhydration of the Tetrafluoroborate Anion, BF4-(H2O)n=1,2,3,4

This study systematically examines the interactions of the tetrafluoroborate anion (BF4-) with up to four water molecules (BF4-(H2O)n=1,2,3,4). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed using a variety of density functional theory (DFT) methods (B3LYP, B3LYP-D3BJ, and M06-2X) and the MP2 ab initio method with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for B, O, and F; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the CCSD(T) ab initio method and the haTZ basis set for the mono- and dihydrate (n = 1, 2) structures. The 2-body:Many-body (2b:Mb) technique, in which CCSD(T) computations capture the 1- and 2-body contributions to the interactions and MP2 computations recover all higher-order contributions, was used to extend these demanding computations to the tri- and tetrahydrate (n = 3, 4) systems. Four, five, and eight new stationary points have been identified for the di-, tri-, and tetrahydrate systems, respectively. The global minimum of the monohydrate adopts a symmetric double ionic hydrogen bond motif with C2v symmetry and an electronic dissociation energy of 13.17 kcal mol-1 at the CCSD(T)/haTZ level of theory. This strong solvent···solute interaction, however, competes with solute···solute interactions in the lowest-energy BF4-(H2O)n=2,3,4 minima that are not seen in the other di-, tri-, or tetrahydrate minima. The latter interactions help increase the 2b:Mb dissociation energies to more than 26, 41, and 51 kcal mol-1 for n = 2, 3, and 4, respectively. Structures that form hydrogen bonds between the solvating water molecules also exhibit the largest shifts in the harmonic OH stretching frequencies for the waters of hydration. These shifts can exceed -280 cm-1 relative to an isolated H2O molecule at the 2b:Mb/haTZ level of theory.

Etalon-Assisted Determination of the Complex Index of Refraction of a Solution for the Study of Strong Cavity-Vibrational Coupling of PF6- in Acetonitrile

A new method is established using an etalon cavity to assist in the determination of the wavelength-dependent complex index of refraction of a solution throughout the mid-infrared range. The results are used to study the cavity-vibration polaritons of PF₆^- in acetonitrile. Mixed states are formed by placing solution inside a pair of parallel plate mirrors with a wavelength-scale spacing, i.e., within an etalon, such that there are cavity states that are angle-tuned into resonance with the strong P-F vibrations. The dominant ν₃ vibrations of PF₆^- consist of nearly triply degenerate oscillations of the partial-positively charged phosphorous against antisymmetric concerted motions of different sets of fluorine atoms with partial negative charges. These vibrations are dominant even though the solute is 29 times less concentrated than the solvent on a molar basis. The first part of the paper describes the method of determining the complex index of refraction of the solution from a combination of etalon transmission maxima and the attenuated total reflection (ATR) absorption spectrum of the solution. The results are presented as an analytical function including a sum of 37 vibrational contributions. Absolute integrated isolated band intensities were determined to be 463 ± 4, 462 ± 7, and 266 ± 4 km/mol for the three ν₃ PF₆^- vibrations at 841.4, 847.4, and 854.0 cm^-1, respectively, which sum to 1191 ± 9 km/mol for the ν₃ band. Then, the results are used to simulate the measured etalon transmission using the transfer matrix (TM) method with and without the ν₃ target vibrations. The etalon transmission simulations reconstruct the position of cavity modes in the absence of target vibrations. They provide input data for the testing of simple quantum mechanical models for the interaction of vibrations with cavity modes and the interactions of vibrations with other vibrations within the molecule and between solute and solvent. The model shows that the nearly degenerate ν₃ vibrations interact with each other with a vibration-vibration coupling of 33 ± 5 cm^-1. This is comparable to the cavity-vibration coupling of 30.4 ± 2.9 cm^-1 of the two strongest vibrations of PF₆^-.

Competition between Solvent-Solvent and Solvent-Solute Interactions in the Microhydration of the Hexafluorophosphate Anion, PF6-(H2O)n=1,2

This study systematically examines the interactions of the hexafluorophosphate anion (PF₆^-) with one or two solvent water molecules (PF₆^-(H₂O)_n where n = 1, 2). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed on each stationary point using a variety of common density functional theory methods (B3LYP, B3LYP-D3, M06-2X, and ωB97XD) and the MP2 and CCSD(T) ab initio methods with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for P, O, and F; denoted as haTZ). Five new stationary points of PF₆^-(H₂O)₂ have been identified, one of which has an electronic energy of approximately 2 kcal mol^-1 lower than the only other dihydrate structure reported for this system. The CCSD(T) computations also reveal that the detailed interactions between PF₆^- and H₂O can be quite difficult to model reliably, with some methods struggling to correctly characterize stationary points for n = 1 or accurately reproduce the vibrational frequency shifts induced by the formation of the hydrated complex. Although the interactions between the solvent and ionic solute are quite strong (CCSD(T) electronic dissociation energy ≈10 kcal mol^-1 for the monohydrate minimum), the solvent-solvent interactions in the lowest-energy PF₆^-(H₂O)₂ minimum give rise to appreciable cooperative effects not observed in the other dihydrate minima. In addition, this newly identified structure exhibits the largest frequency shifts in the OH stretching vibrations for the waters of hydration (with Δω exceeding -100 cm^-1 relative to the values for an isolated H₂O molecule).

Ion hydration free energies and water surface potential in water nano drops: The cluster pair approximation and the proton hydration Gibbs free energy in solution

We estimate both single ion hydration Gibbs free energies in water droplets, comprising from 50 to 1000 molecules, and water/vacuum surface potentials in pure water droplets comprising up to 10 000 molecules. We consider four ions, namely, Li⁺, NH₄ ⁺, F^-, and Cl^-, and we model their hydration process and water/water interactions using polarizable force fields based on an induced point dipole approach. We show both ion hydration Gibbs free energies and water surface potentials to obey linear functions of the droplet radius as soon as droplets comprising a few hundred water molecules. Moreover, we also show that the differences in anion/cation hydration Gibbs free energies in droplets obey a different regime in large droplets than in small clusters comprising no more than six water molecules, in line with the earlier results computed from standard additive point charge force fields. Hence, both point charge and more sophisticated induced point dipole molecular modeling approaches suggest that methods considering only the thermodynamical properties of small ion/water clusters to estimate the absolute proton hydration Gibbs free energy in solution are questionable. In particular, taking into account the data of large ion/water droplets may yield a proton hydration Gibbs free energy in solution value to be shifted by several k_BT units compared to small clusters-based approaches.