Do Specific Ion Effects on Collective Relaxation Arise from Perturbation of Hydrogen-Bonding Network Structure?

The change in the transport properties (i.e., water diffusivity, shear viscosity, etc.) when adding salts to water has been used to classify ions as either being chaotropic or kosmotropic, a terminology based on the presumption that this phenomenon arises from respective breakdown or enhancement of the hydrogen-bonding network structure. Recent quasi-elastic neutron scattering measurements of the collective structural relaxation time, τC, in aqueous salt solutions were interpreted as confirming this proposed origin of ion effects on the dynamics of water. However, we find similar changes in τC in the same salt solutions based on molecular dynamics (MD) simulations using a coarse-grained water model in which no hydrogen bonding exists, challenging this conventional interpretation of mobility change resulting from the addition of salts to water. A thorough understanding of specific ion effects should be useful in diverse material manufacturing and biomedical applications, where these effects are prevalent, but poorly understood.

Percolation transition and bimodal density distribution in hydrogen fluoride

Hydrogen-bond networks in associating fluids can be extremely robust and characterize the topological properties of the liquid phase, as in the case of water, over its whole domain of stability and beyond. Here, we report on molecular dynamics simulations of hydrogen fluoride (HF), one of the strongest hydrogen-bonding molecules. HF has more limited connectivity than water but can still create long, dynamic chains, setting it apart from most other small molecular liquids. Our simulation results provide robust evidence of a second-order percolation transition of HF's hydrogen bond network occurring below the critical point. This behavior is remarkable as it underlines the presence of two different cohesive mechanisms in liquid HF, one at low temperatures characterized by a spanning network of long, entangled hydrogen-bonded polymers, as opposed to short oligomers bound by the dispersion interaction above the percolation threshold. This second-order phase transition underlines the presence of marked structural heterogeneity in the fluid, which we found in the form of two liquid populations with distinct local densities.

A variational formulation of the Harris functional as a correction to approximate Kohn-Sham density functional theory

Accurate descriptions of intermolecular interactions are of great importance in simulations of molecular liquids. We present an electronic structure method that combines the accuracy of the Harris functional approach with the computational efficiency of approximately linear-scaling density functional theory (DFT). The non-variational nature of the Harris functional has been addressed by constructing a Lagrangian energy functional, which restores the variational condition by imposing stationarity with respect to the reference density. The associated linear response equations may be treated with linear-scaling efficiency in an atomic orbital based scheme. Key ingredients to describe the structural and dynamical properties of molecular systems are the forces acting on the atoms and the stress tensor. These first-order derivatives of the Harris Lagrangian have been derived and implemented in consistence with the energy correction. The proposed method allows for simulations with accuracies close to the Kohn-Sham DFT reference. Embedded in the CP2K program package, the method is designed to enable ab initio molecular dynamics simulations of molecular solutions for system sizes of several thousand atoms. Available subsystem DFT methods may be used to provide the reference density required for the energy correction at near linear-scaling efficiency. As an example of production applications, we applied the method to molecular dynamics simulations of the binary mixtures cyclohexane-methanol and toluene-methanol, performed within the isobaric-isothermal ensemble, to investigate the hydrogen bonding network in these non-ideal mixtures.

The influence of intermolecular correlations on the infrared spectrum of liquid dimethyl sulfoxide

Dimethyl sulfoxide (DMSO) is routinely applied as an excellent, water-miscible solvent and chemical reagent. Some of the most important data concerning its liquid structure were obtained using infrared (IR) spectroscopy. However, the actual extent of intermolecular correlations that connect the isolated monomer spectrum to the IR response of the bulk liquid is poorly studied thus far. Using ab initio molecular dynamics (AIMD) simulations, IR spectra of liquid DMSO are obtained here from first principles and further analyzed using an array of sophisticated spectral decomposition techniques. The calculated spectra when unfolded in space reveal non-trivial spatial correlations underlying the IR response of liquid DMSO. It is unequivocally demonstrated that some of the fundamental vibrations visible in the intramolecular limit are effectively suppressed by the solvation environment due to symmetry reasons and thus disappear in the bulk limit, escaping experimental detection. Overall, DMSO as an aprotic solvent with dominant dipole-dipole interactions displays strong intermolecular correlations that contribute significantly to the IR spectra, on par with the situation observed in strongly associated liquids, such as water.

"On-The-Fly" Calculation of the Vibrational Sum-Frequency Generation Spectrum at the Air-Water Interface

In the present work, we provide an electronic structure based method for the "on-the-fly" determination of vibrational sum frequency generation (v-SFG) spectra. The predictive power of this scheme is demonstrated at the air-water interface. While the instantaneous fluctuations in dipole moment are obtained using the maximally localized Wannier functions, the fluctuations in polarizability are approximated to be proportional to the second moment of Wannier functions. The spectrum henceforth obtained captures the signatures of hydrogen bond stretching, bending, as well as low-frequency librational modes.

Solvation Shell of the Nitrite Ion in Water: An Ab Initio Molecular Dynamics Study

We performed ab initio molecular dynamics simulation of a nitrite ion in water to investigate the structural and dynamical properties of its hydration shell. The nitrite ion is found to exhibit strong asymmetry toward hydrogen bonding due to its two different types of hydrogen bond acceptor sites. This difference is better captured through further partitioning of the hydration shell into its proximal and distal regions. The frequency shifts of the stretch modes of hydration shell water reveal that the nitrogen site forms a stronger hydrogen bond than its oxygen sites with the latter forming hydrogen bonds, which are similar in strength to that between a pair of water molecules. The escape dynamics of water from the hydration shell is found to be rather slow, which seems to classify the nitrite ion as a structure-maker. However, the dynamics of orientational and hydrogen bond relaxation reveal a faster mobility of water molecules in the hydration shell than bulk water in spite of strong ion-water interactions. It is found that the nitrite ion can hold water molecules in its solvation shell and still make them rotate fast in its vicinity through switching of their hydrogen bonds between its nitrogen and oxygen acceptor sites. The dipole moment of the solute in water is also calculated in the present study.

New Molecular-Mechanics Model for Simulations of Hydrogen Fluoride in Chemistry and Biology

Hydrogen fluoride (HF) is the most polar diatomic molecule and one of the simplest molecules capable of hydrogen-bonding. HF deviates from ideality both in the gas phase and in solution and is thus of great interest from a fundamental standpoint. Pure and aqueous HF solutions are broadly used in chemical and industrial processes, despite their high toxicity. HF is a stable species also in some biological conditions, because it does not readily dissociate in water unlike other hydrogen halides; yet, little is known about how HF interacts with biomolecules. Here, we set out to develop a molecular-mechanics model to enable computer simulations of HF in chemical and biological applications. This model is based on a comprehensive high-level ab initio quantum chemical investigation of the structure and energetics of the HF monomer and dimer; (HF)_n clusters, for n = 3-7; various clusters of HF and H₂O; and complexes of HF with analogs of all 20 amino acids and of several commonly occurring lipids, both neutral and ionized. This systematic analysis explains the unique properties of this molecule: for example, that interacting HF molecules favor nonlinear geometries despite being diatomic and that HF is a strong H-bond donor but a poor acceptor. The ab initio data also enables us to calibrate a three-site molecular-mechanics model, with which we investigate the structure and thermodynamic properties of gaseous, liquid, and supercritical HF in a wide range of temperatures and pressures; the solvation structure of HF in water and of H₂O in liquid HF; and the free diffusion of HF across a lipid bilayer, a key process underlying the high cytotoxicity of HF. Despite its inherent simplifications, the model presented significantly improves upon previous efforts to capture the properties of pure and aqueous HF fluids by molecular-mechanics methods and to our knowledge constitutes the first parameter set calibrated for biomolecular simulations.

Effects of mixing between short-chain and branched-chain alcohols in protonated clusters

The previous analysis of the neat protonated branched-chain alcohol clusters revealed the impact of steric repulsion and dispersion of the bulky alkyl group on the hydrogen-bonded (H-bonded) structures and their temperature-dependence. To further understand the influence of the alkyl groups in H-bonded clusters, we studied the mixing of the two extremes of alcohols, methanol (MeOH) and tert-butyl alcohol (t-BuOH), with an excess proton. Infrared spectroscopy and a structural search of first principles calculations on the size-selected clusters H⁺(MeOH)_m(t-BuOH)_t (m + t = 4 and 5) were conducted. Temperature-dependence of the dominant H-bonded structures was explored by the Ar-tagging technique and quantum harmonic superposition approach. By introducing the dispersion-corrected density functional theory methods, it was shown that the effects of dispersion due to the bulky alkyl groups in the mixed clusters cannot be ignored for t≥ 2. The computational results qualitatively depicted the characteristics of the observed IR spectra, but overestimation of the temperature-dependence with dispersion correction was clearly seen due to the unbalanced correction between linear H-bonded structures and compact cyclic ones. These results demonstrate the importance of extensive investigation and benchmarks on different levels of theory, and that a properly sampled structure database is crucial to evaluate theoretical models.

An ab initio molecular dynamics study of benzene in water at supercritical conditions: Structure, dynamics, and polarity of hydration shell water and the solute

Anisotropic structure and dynamics of the hydration shell of a benzene solute in supercritical water are investigated by means of ab initio molecular dynamics simulations. The polarity and structural distortion of the benzene solute in supercritical water are also investigated in this study. Calculations are done at 673 K for three different densities of the solvent. The simulations are carried out using the Becke-Lee-Yang-Parr (BLYP) and also the Becke-Lee-Yang-Parr functional including dispersion corrections of Grimme (BYLP-D). The structural anisotropy is found to exist even at supercritical conditions as elucidated by the radial distribution functions of different conical regions and also by angular and spatial distribution functions. The benzene-water πH-bond and also the water-water hydrogen bonds are found to exist even at the supercritical temperature of 673 K. However, the numbers of these hydrogen bonds are reduced substantially with a decrease in water density. The water molecules in the axial region of benzene are found to be preferably oriented with one OH vector pointing toward the benzene ring, whereas the water molecules located in the equatorial region are found to orient their dipoles mostly parallel to the ring plane. The orientational distributions, however, are found to be rather broad at the supercritical temperature due to thermal fluctuations. Although the water molecules have faster dynamics at these supercritical conditions, a slight difference is observed in the dynamics of the solvation shell and bulk molecules. The conformational flexibility of the ring is found to be enhanced which causes an increase in polarity of the benzene solute in water under supercritical conditions.

Ab initio molecular dynamics studies of hydroxide coordination of alkaline earth metals and uranyl

Ab initio molecular dynamics (AIMD) simulations of the Mg2+, Ca2+, Sr2+ and UO22+ ions in either a pure aqueous environment or an environment containing two hydroxide ions have been carried out at the density functional level of theory, employing the generalised gradient approximation via the PBE exchange-correlation functional. Calculated mean M-O bond lengths in the first solvation shell of the aquo systems compared very well to existing experimental and computational literature, with bond lengths well within values measured previously and coordination numbers in line with previously calculated values. When applied to systems containing additional hydroxide ions, the methodology revealed increased bond lengths in all systems. Proton transfer events (PTEs) were recorded and were found to be most prevalent in the strontium hydroxide systems, likely due to the low charge density of the ion and the consequent lack of hydroxide coordination. For all alkaline earths, intrashell PTEs which occurred outside of the first solvation shell were most prevalent. Only three PTEs were identified in the entire simulation data of the uranium dihydroxide system, indicating the clear impact of the increased charge density of the hexavalent uranium ion on the strength of metal-oxygen bonds in aqueous solution. Broadly, systems containing more charge dense ions were found to exhibit fewer PTEs than those containing ions of lower charge density.

Temperature dependence of the ultrafast vibrational echo spectroscopy of OD modes in liquid water from first principles simulations

The temperature dependence of the vibrational spectral diffusion of OD modes in liquid water is investigated through calculations of vibrational echo spectral observables from first principles molecular dynamics.

Hydrogen bond network structures of protonated short-chain alcohol clusters

Protonated alcohol clusters enable extraction of the physical essence of the nature of hydrogen bond networks.

Static vs dynamic DFT prediction of IR spectra of flexible molecules in the condensed phase: The (ClCF2CF(CF3)OCF2CH3) liquid as a test case

First-principles molecular dynamics (FPMD) simulations in the framework of Density Functional Theory (DFT) are carried out for the prediction of the infrared spectrum of the fluorinated molecule ClCF₂CF(CF₃)OCF₂CH₃ in liquid and gas phase. This molecule is characterized by a flexible structure, allowing the co-existence of several stable conformers, that differ by values of the torsional angles. FPMD computed spectra are compared to the experimental ones, and to Boltzmann weighted IR spectra based on gas phase calculations.

Anisotropic structure and dynamics of the solvation shell of a benzene solute in liquid water from ab initio molecular dynamics simulations

The anisotropic structure and dynamics of the hydration shell of a benzene solute in liquid water have been investigated by means of ab initio molecular dynamics simulations using the BLYP (Becke–Lee–Yang–Parr) and dispersion corrected BLYP-D functionals.

Ultrafast Vibrational Echo Spectroscopy of Liquid Water from First-Principles Simulations

Vibrational echo spectroscopy has become a powerful technique to study vibrational spectral diffusion in water and aqueous solutions. The dynamics of vibrational spectral diffusion is intimately related to the hydrogen bond fluctuations in liquid water and other hydrogen bonded liquids. Earlier theoretical calculations of vibrational echo spectroscopy of aqueous systems were based on classical molecular dynamics simulations involving empirical force fields of water. In the current work, we have employed the method of ab initio molecular dynamics simulation to calculate the spectral observables of vibrational echo and two-dimensional infrared (2D-IR) spectroscopy of liquid water at room temperature under Condon and cumulant approximations. The time scales extracted from the temporal decay of the frequency-time correlation function (FTCF), short-time slope of three pulse photon echo (SP3E), dynamic line width (DLW), and the slope of nodal line of 2D-IR spectra are found to be in reasonably close agreement with each other which reinforces the assertion that signatures of FTCF can be captured using three pulse photon echo and 2D-IR spectroscopy.

Quantum mechanical force field for hydrogen fluoride with explicit electronic polarization

The explicit polarization (X-Pol) theory is a fragment-based quantum chemical method that explicitly models the internal electronic polarization and intermolecular interactions of a chemical system. X-Pol theory provides a framework to construct a quantum mechanical force field, which we have extended to liquid hydrogen fluoride (HF) in this work. The parameterization, called XPHF, is built upon the same formalism introduced for the XP3P model of liquid water, which is based on the polarized molecular orbital (PMO) semiempirical quantum chemistry method and the dipole-preserving polarization consistent point charge model. We introduce a fluorine parameter set for PMO, and find good agreement for various gas-phase results of small HF clusters compared to experiments and ab initio calculations at the M06-2X/MG3S level of theory. In addition, the XPHF model shows reasonable agreement with experiments for a variety of structural and thermodynamic properties in the liquid state, including radial distribution functions, interaction energies, diffusion coefficients, and densities at various state points.

Critical analysis of the accuracy of models predicting or extracting liquid structure information

This work aims at a critical assessment of properties predicting or extracting information on the density and structure of liquids. State-of-the-art NVT and NpT molecular dynamics (MD) simulations have been performed on five liquids: methanol, chloroform, acetonitrile, tetrahydrofuran, and ethanol. These simulations allow the computation of properties based on first principles, including the equilibrium density and radial distribution functions (RDFs), characterizing the liquid structure. Refinements have been incorporated in the MD simulations by taking into account basis set superposition errors (BSSE). An extended BSSE model for an instantaneous evaluation of the BSSE corrections has been proposed, and their impact on the liquid properties has been assessed. If available, the theoretical RDFs have been compared with the experimentally derived RDFs. For some liquids, significant discrepancies have been observed, and a profound but critical investigation is presented to unravel the origin of these deficiencies. This discussion is focused on tetrahydrofuran where the experiment reveals some prominent peaks completely missing in any MD simulation. Experiments providing information on liquid structure consist mainly of neutron diffraction measurements offering total structure factors as the primary observables. The splitting of these factors in reciprocal space into intra- and intermolecular contributions is extensively discussed, together with their sensitivity in reproducing correct RDFs in coordinate space.