Ion Pair Potentials-of-Mean-Force in Water

In this paper we compare different force fields that are widely used (Gromacs, Charmm-22/x-Plor, Charmm-27, Amber-1999, OPLS-AA) in biophysical simulations containing aqueous NaCl. We show that the uncertainties of the microscopic parameters of, in particular, sodium and, to a lesser extent, chloride translate into large differences in the computed radial-distribution functions. This uncertainty reflects the incomplete experimental knowledge of the structural properties of io...

We present a method to derive implicit solvent models of electrolyte solutions from all-atom descriptions; providing analytical expressions of the thermodynamic and structural properties of the ions consistent with the underlying explicit solvent representation. Effective potentials between ions in solution are calculated to perform perturbation theory calculations, in order to derive the best possible description in terms of charged hard spheres. Applying this method to NaCl...

Long-range ion induced water-water correlations were recently observed in femtosecond elastic second harmonic scattering experiments of electrolyte solutions. To further the qualitative understanding of these correlations, we derive an analytical expression that quantifies ion induced dipole-dipole correlations in a non-interacting gas of dipoles. This model is a logical extension of Debye-H\"uckel theory that can be used to qualitatively understand how the combined electric ...

Single ion solvation free energies are one of the most important properties of electrolyte solutions and yet there is ongoing debate about what these values are. Only the values for neutral ion pairs are known. Here, we use DFT interaction potentials with molecular dynamics simulation (DFT-MD) combined with a modified version of the quasi-chemical theory (QCT) to calculate these energies for the lithium and fluoride ions. A method to correct for the error in the DFT functiona...

The solvation structures and ion dynamics of CaCl$_2$ aqueous electrolytes have been investigated using ab initio molecular dynamics simulations and molecular dynamics simulations with deep learning potentials. We found multiple solvation structures around the Ca$^{2+}$ ion, including fully hydrated single Ca$^{2+}$ ion, Ca-Cl contact ion pair, and Ca-2Cl bridged ion pair, could coexist. The ion-pairing condition plays an important role in the translational and orientational ...

Despite the ubiquity of nanoconfined aqueous electrolytes, a theoretical framework that accounts for the nonlinear coupling of water and ion polarization is still missing. We introduce a nonlocal and nonlinear field theory for the nanoscale polarization of ions and water and derive the electrolyte dielectric properties as a function of salt concentration to first order in a loop expansion. Classical molecular dynamics simulations are favorably compared with the calculated die...

The propensity for ion-pairing can often dictate the thermodynamic and kinetic properties of electrolyte solutions. Fast and accurate estimates of ion-pairing can thus be extremely valuable for supplementing design and screening efforts for novel electrolytes. Here, we introduce an efficient cluster model to estimate the local ion-pair potential-of-mean-force (PMF) between ionic solutes in electrolytes. The model incorporates an ion-pair and a few layers of explicit solvent i...

In this paper we extend the familiar continuum electrostatic model with a perturbation to the usual macroscopic boundary condition. The perturbation is based on the mean spherical approximation (MSA), to derive a multiscale hydration-shell boundary condition (HSBC). We show that the HSBC/MSA model reproduces MSA predictions for Born ions in a variety of polar solvents, including both protic and aprotic solvents. Importantly, the HSBC/MSA model predicts not only solvation free...

The Poisson-Boltzmann mean-field description of ionic solutions has been successfully used in predicting charge distributions and interactions between charged macromolecules. While the electrostatic model of charged fluids, on which the Poisson-Boltzmann description rests, and its statistical mechanical consequences have been scrutinized in great detail, much less is understood about its probable shortcomings when dealing with various aspects of real physical, chemical and bi...

Ions are critical to the structure and stability of polyelectrolytes such as nucleic acids. In this work, we systematically calculated the potentials of mean force between two like-charged nanoparticles in salt solutions by Monte Carlo simulations. The pseudo-spring method is employed to calculate the potential of mean force and compared systematically with the inversed-Boltzmann method. An effective attraction is predicted between two like-charged nanoparticles in divalent/t...