Absolute ion hydration free energy scale and the surface potential of water via quantum simulation

Deep neural network based quantum simulations and quasichemical theory for accurate modeling of molten salt thermodynamics

With dual goals of efficient and accurate modeling of solvation thermodynamics in molten salt liquids, we employ ab initio molecular dynamics (AIMD) simulations, deep neural network interatomic potentials (NNIP), and quasichemical theory (QCT) to calculate the excess chemical potentials for the solute ions Na+ and Cl- in the molten NaCl liquid. NNIP-based molecular dynamics simulations accelerate the calculations by 3 orders of magnitude and reduce the uncertainty to 1 kcal mol-1. Using the Density Functional Theory (DFT) level of theory, the predicted excess chemical potential for the solute ion pair is -178.5 ± 1.1 kcal mol-1. A quantum correction of 13.7 ± 1.9 kcal mol-1 is estimated via higher-level quantum chemistry calculations, leading to a final predicted ion pair excess chemical potential of -164.8 ± 2.2 kcal mol-1. The result is in good agreement with a value of -163.5 kcal mol-1 obtained from thermo-chemical tables. This study validates the application of QCT and NNIP simulations to the molten salt liquids, allowing for significant insights into the solvation thermodynamics crucial for numerous molten salt applications.

The enhanced dissociation and associated surface structure of the anesthetic propofol at the water interface: vibrational sum frequency generation study

Propofol, the most administered drug for general anesthesia, affects the acid-base equilibrium at the interfacial region of arterial blood. Hence, the structure of propofol at the water interface under different pH conditions has been measured using the surface-selective vibrational sum frequency generation (VSFG) technique to understand the hydration as well as the dissociation of propofol at the water interface. Propofol remains in its neutral form at pH ≤ 5.8 in which the OH group of propofol forms a hydrogen bond with interfacial water molecules, where a few interfacial water molecules also interact with the π electron density of propofol. By contrast, propofol prefers to be in the deprotonated state at pH ≥ 7, due to which the surface of water becomes negatively charged and hence the interfacial water becomes oriented and the intensity of the OH stretch of water is enhanced. The pK_a of propofol at the water interface is ∼three units lower than in the bulk medium indicating that the dissociation of propofol is notably enhanced at the water interface. These VSFG studies suggest that, unlike the bulk, propofol prefers to be in the charged state at the water interface under physiological conditions, which may be important in understanding its diffusion and acid-base equilibrium in the interfacial arterial blood region.

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion-solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. Here we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

Electric fields and potentials in condensed phases

The electric fields and potentials inside and at the interface of matter are relevant to many branches of physics, chemistry, and biology. Accurate quantification of these fields and/or potentials is essential to control and exploit chemical and physical transformations. Before we understand the response of matter to external fields, it is first important to understand the intrinsic interior and interfacial fields and potentials, both classically and quantum mechanically, as well as how they are probed experimentally. Here we compare and contrast, beginning with the hydrogen atom in vacuum and ending with concentrated aqueous NaCl electrolyte, both classical and quantum mechanical electric potentials and fields. We make contact with experimental vibrational Stark, electrochemical, X-ray, and electron spectroscopic probes of these potentials and fields, outline relevant conceptual difficulties, and underscore the advantage of electron holography as a basis to better understand electrostatics in matter.

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined.In this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion-ion, ion-water, and water-water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions.DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air-water interface, ions in bulk, and at the air-water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Advanced Electrostatic Model for Monovalent Ions Based on Ab Initio Energy Decomposition

Ions play important roles in the structures and functions of biomolecules. In biomolecular simulations, ions either directly interact with biomolecules or provide an ionic environment that influences electrostatic interactions of solutes. The AMOEBA+ water model has demonstrated significant advancement of the classical force field for describing molecular interactions due to its improvements on the functional forms to account for essential physics. This work expands the applicability of the AMOEBA+ model toward alkali metal (Li, Na, K, Rb, and Cs) and halogen (F, Cl, Br, and I) ions. Various quantum chemical data on ion-ion and ion-water interactions, experimental ion hydration free energies, and lattice energies of salt crystals are used in the parametrization. The final parameters are verified with other properties outside of the parametrization data, including lattice energies of additional salt crystals and ionic activity coefficients in solution. The new model captures a wide range of ion properties from the gas phase to solution phase and crystals. More importantly, AMOEBA+ provides energy components that are consistent with ab initio energy decomposition. Thus, we expect AMOEBA+ to be more general, transferable, and valuable for the interpretation of intermolecular forces in efficient classical simulations.