-
Scalable molecular simulation of electrolyte solutions with quantum chemical accuracy
Authors:
Junji Zhang,
Joshua Pagotto,
Tim Gould,
Timothy T. Duignan
Abstract:
Unleashing the predictive power of molecular dynamics (MD), Neural Network Potentials (NNPs) trained on Density Functional Theory (DFT) calculations are revolutionizing our ability to simulate chemical systems with unprecedented accuracy and efficiency. Electrolyte solutions are a natural initial system to apply this tool to because they are critically important for a wide range of applications an…
▽ More
Unleashing the predictive power of molecular dynamics (MD), Neural Network Potentials (NNPs) trained on Density Functional Theory (DFT) calculations are revolutionizing our ability to simulate chemical systems with unprecedented accuracy and efficiency. Electrolyte solutions are a natural initial system to apply this tool to because they are critically important for a wide range of applications and their properties cannot currently be predicted. Unfortunately, however most DFT approximations are not sufficiently accurate to predict many practically relevant properties of electrolytes. Additionally, tracking the position of every atom in a system during molecular simulations is inherently limited, even with NNP-MD. Here, we use a state-of-the-art DFT approximation to demonstrate highly accurate all-atom NNPs with minimal training data. We demonstrate that NNPs can reliably be recursively trained on a subset of their own output to enable coarse-grained continuum solvent molecular simulations that can access much longer timescales. We apply our technique to simulate lithium chloride, potassium chloride, and lithium bromide in water. We reproduce key structural, thermodynamic, and kinetic properties of these solutions in agreement with experimental data. The formation of a previously unknown Li cation dimer is observed, along with identical anion-anion interactions of chloride and bromide. Finally, the coarse-grained model is capable of reproducing crystal phase behavior and infinite dilution pairing free energies despite being trained solely on moderate concentration solutions, disproving the notion that NNPs are only useful for interpolation. This approach should be scalable to determine the properties of electrolyte solutions over a much wider range of conditions and compositions than is possible experimentally.
△ Less
Submitted 2 June, 2024; v1 submitted 19 October, 2023;
originally announced October 2023.
-
The surface potential explains ion specific bubble coalescence inhibition
Authors:
Timothy T. Duignan
Abstract:
Hypothesis: Some ions can prevent bubbles from coalescing in water. The Gibbs-Marangoni pressure has been proposed as an explanation of this phenomenon. This repulsive pressure occurs during thin film drainage whenever surface enhanced or surface depleted solutes are present. However, bubble coalescence inhibition is known to depend on which particular combination of ions are present in a peculiar…
▽ More
Hypothesis: Some ions can prevent bubbles from coalescing in water. The Gibbs-Marangoni pressure has been proposed as an explanation of this phenomenon. This repulsive pressure occurs during thin film drainage whenever surface enhanced or surface depleted solutes are present. However, bubble coalescence inhibition is known to depend on which particular combination of ions are present in a peculiar and unexplained way. This dependence may be explained by the electrostatic surface potential created by the distribution of ions at the interface, which will alters the natural surface propensity of the ions and hence the Gibbs-Marangoni pressure.
Calculations: A generalised form of the Gibbs-Marangoni pressure is derived for a mixture of solutes and a modified Poisson-Boltzmann equation model is used to calculate this pressure for five different electrolyte solutions made up of four different ions.
Findings: Combining ions with differing surface propensities, i.e., one enhanced and one depleted, creates a significant electrostatic surface potential which dampens the natural surface propensity of these ions, resulting in a reduced Gibbs-Marangoni pressure, which allows bubble coalescence. This mechanism explains why the ability of electrolytes to inhibit bubble coalescence is correlated with surface tension for pure electrolytes but not for mixed electrolytes.
△ Less
Submitted 12 May, 2021; v1 submitted 2 September, 2020;
originally announced September 2020.
-
Mass Density Fluctuations in Quantum and Classical descriptions of Liquid Water
Authors:
Mirza Galib,
Timothy T. Duignan,
Yannick Misteli,
Marcel D. Baer,
Gregory K. Schenter,
Jurg Hutter,
Christopher J. Mundy
Abstract:
First principles molecular dynamics simulation protocol is established using revised functional of Perdew-Burke-Ernzerhof (revPBE) in conjunction with Grimme's third generation of dispersion (D3) correction to describe properties of water at ambient conditions. This study also demonstrates the consistency of the structure of water across both isobaric (NpT) and isothermal (NVT) ensembles. Going be…
▽ More
First principles molecular dynamics simulation protocol is established using revised functional of Perdew-Burke-Ernzerhof (revPBE) in conjunction with Grimme's third generation of dispersion (D3) correction to describe properties of water at ambient conditions. This study also demonstrates the consistency of the structure of water across both isobaric (NpT) and isothermal (NVT) ensembles. Going beyond the standard structural benchmarks for liquid water, we compute properties that are connected to both local structure and mass density uctuations that are related to concepts of solvation and hydrophobicity. We directly compare our revPBE results to the Becke-Lee-Yang-Parr (BLYP) plus Grimme dispersion corrections (D2) and both the empirical fixed charged model (SPC/E) and many body interaction potential model (MB-pol) to further our understanding of how the computed properties herein depend on the form of the interaction potential.
△ Less
Submitted 20 March, 2017;
originally announced March 2017.
-
Real single ion solvation free energies with quantum mechanical simulation
Authors:
Timothy T. Duignan,
Marcel D. Baer,
Gregory K. Schenter,
Christopher J. Mundy
Abstract:
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 l…
▽ More
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 functional is developed and very good agreement with the experimental value for the lithium fluoride pair is obtained. Moreover, this method partitions the energies into physically intuitive terms such as surface potential, cavity and charging energies which are amenable to descriptions with reduced models. Our research suggests that lithium's solvation free energy is dominated by the free energetics of a charged hard sphere, whereas fluoride exhibits significant quantum mechanical behavior that cannot be simply described with a reduced model.
△ Less
Submitted 18 July, 2017; v1 submitted 10 March, 2017;
originally announced March 2017.
-
Electrostatic solvation free energies of charged hard spheres using molecular dynamics with density functional theory interactions
Authors:
Timothy T. Duignan,
Marcel D. Baer,
Gregory K. Schenter,
Christopher J. Mundy
Abstract:
Determining the solvation free energies of single ions in water is one of the most fundamental problems in physical chemistry and yet many unresolved questions remain. In particular, the ability to decompose the solvation free energy into simple and intuitive contributions will have important implications for models of electrolyte solution. Here, we provide definitions of the various types of sing…
▽ More
Determining the solvation free energies of single ions in water is one of the most fundamental problems in physical chemistry and yet many unresolved questions remain. In particular, the ability to decompose the solvation free energy into simple and intuitive contributions will have important implications for models of electrolyte solution. Here, we provide definitions of the various types of single ion solvation free energies based on different simulation protocols. We calculate solvation free energies of charged hard spheres using density functional theory interaction potentials with molecular dynamics simulation (DFT-MD) and isolate the effects of charge and cavitation, comparing to the Born (linear response) model. We show that using uncorrected Ewald summation leads to unphysical values for the single ion solvation free energy and that charging free energies for cations are approximately linear as a function of charge but that there is a small non-linearity for small anions. The charge hydration asymmetry (CHA) for hard spheres, determined with quantum mechanics, is much larger than for the analogous real ions. This suggests that real ions, particularly anions, are significantly more complex than simple charged hard spheres, a commonly employed representation.
△ Less
Submitted 18 July, 2017; v1 submitted 16 February, 2017;
originally announced February 2017.