1
|
O’Neill N, Shi BX, Fong K, Michaelides A, Schran C. To Pair or not to Pair? Machine-Learned Explicitly-Correlated Electronic Structure for NaCl in Water. J Phys Chem Lett 2024; 15:6081-6091. [PMID: 38820256 PMCID: PMC11181334 DOI: 10.1021/acs.jpclett.4c01030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Revised: 05/23/2024] [Accepted: 05/24/2024] [Indexed: 06/02/2024]
Abstract
The extent of ion pairing in solution is an important phenomenon to rationalize transport and thermodynamic properties of electrolytes. A fundamental measure of this pairing is the potential of mean force (PMF) between solvated ions. The relative stabilities of the paired and solvent shared states in the PMF and the barrier between them are highly sensitive to the underlying potential energy surface. However, direct application of accurate electronic structure methods is challenging, since long simulations are required. We develop wave function based machine learning potentials with the random phase approximation (RPA) and second order Møller-Plesset (MP2) perturbation theory for the prototypical system of Na and Cl ions in water. We show both methods in agreement, predicting the paired and solvent shared states to have similar energies (within 0.2 kcal/mol). We also provide the same benchmarks for different DFT functionals as well as insight into the PMF based on simple analyses of the interactions in the system.
Collapse
Affiliation(s)
- Niamh O’Neill
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Cavendish
Laboratory, Department of Physics, University
of Cambridge, Cambridge CB3 0HE, United
Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Benjamin X. Shi
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Kara Fong
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Angelos Michaelides
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Christoph Schran
- Cavendish
Laboratory, Department of Physics, University
of Cambridge, Cambridge CB3 0HE, United
Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| |
Collapse
|
2
|
Kastinen T, Batys P, Tolmachev D, Laasonen K, Sammalkorpi M. Ion-Specific Effects on Ion and Polyelectrolyte Solvation. Chemphyschem 2024:e202400244. [PMID: 38712639 DOI: 10.1002/cphc.202400244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Revised: 05/02/2024] [Accepted: 05/02/2024] [Indexed: 05/08/2024]
Abstract
Ion-specific effects on aqueous solvation of monovalent counter ions, Na+ ${^+ }$ , K+ ${^+ }$ , Cl- ${^- }$ , and Br- ${^- }$ , and two model polyelectrolytes (PEs), poly(styrene sulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) were here studied with ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations based on the OPLS-aa force-field which is an empirical fixed point-charge force-field. Ion-specific binding to the PE charge groups was also characterized. Both computational methods predict similar response for the solvation of the PEs but differ notably in description of ion solvation. Notably, AIMD captures the experimentally observed differences in Cl- ${^- }$ and Br- ${^- }$ anion solvation and binding with the PEs, while the classical MD simulations fail to differentiate the ion species response. Furthermore, the findings show that combining AIMD with the computationally less costly classical MD simulations allows benefiting from both the increased accuracy and statistics reach.
Collapse
Affiliation(s)
- Tuuva Kastinen
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076, Aalto, Finland
- Academy of Finland Center of Excellence in Life-Inspired Hybrid Materials (LIBER), Aalto University, P.O. Box 16100, 00076, Aalto, Finland
- Faculty of Engineering and Natural Sciences, Tampere University, P.O. Box 541, 33014, Tampere University, Finland
| | - Piotr Batys
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, PL-30239, Krakow, Poland
| | - Dmitry Tolmachev
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076, Aalto, Finland
- Academy of Finland Center of Excellence in Life-Inspired Hybrid Materials (LIBER), Aalto University, P.O. Box 16100, 00076, Aalto, Finland
| | - Kari Laasonen
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076, Aalto, Finland
| | - Maria Sammalkorpi
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076, Aalto, Finland
- Academy of Finland Center of Excellence in Life-Inspired Hybrid Materials (LIBER), Aalto University, P.O. Box 16100, 00076, Aalto, Finland
- Department of Bioproducts and Biosystems, Aalto University, P.O. Box 16100, 00076, Aalto, Finland
| |
Collapse
|
3
|
Savoj R, Agnew H, Zhou R, Paesani F. Molecular Insights into the Influence of Ions on the Water Structure. I. Alkali Metal Ions in Solution. J Phys Chem B 2024; 128:1953-1962. [PMID: 38373140 DOI: 10.1021/acs.jpcb.3c08150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2024]
Abstract
In this study, we explore the impact of alkali metal ions (Li+, Na+, K+, Rb+, and Cs+) on the hydration structure of water using molecular dynamics simulations carried out with MB-nrg potential energy functions (PEFs). Our analyses include radial distribution functions, coordination numbers, dipole moments, and infrared spectra of water molecules, calculated as a function of solvation shells. The results collectively indicate a highly local influence of all of the alkali metal ions on the hydrogen-bond network established by the surrounding water molecules, with the smallest and most densely charged Li+ ion exerting the most pronounced effect. Remarkably, the MB-nrg PEFs demonstrate excellent agreement with available experimental data for the position and size of the first solvation shells, underscoring their potential as predictive models for realistic simulations of ionic aqueous solutions across various thermodynamic conditions and environments.
Collapse
Affiliation(s)
- Roya Savoj
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Henry Agnew
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Ruihan Zhou
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- Halicioğlu Data Science Institute, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
4
|
Dasgupta S, Palos E, Pan Y, Paesani F. Balance between Physical Interpretability and Energetic Predictability in Widely Used Dispersion-Corrected Density Functionals. J Chem Theory Comput 2024; 20:49-67. [PMID: 38150541 DOI: 10.1021/acs.jctc.3c00903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2023]
Abstract
We assess the performance of different dispersion models for several popular density functionals across a diverse set of noncovalent systems, ranging from the benzene dimer to molecular crystals. By analyzing the interaction energies and their individual components, we demonstrate that there exists variability across different systems for empirical dispersion models, which are calibrated for reproducing the interaction energies of specific systems. Thus, parameter fitting may undermine the underlying physics, as dispersion models rely on error compensation among the different components of the interaction energy. Energy decomposition analyses reveal that, the accuracy of revPBE-D3 for some aqueous systems originates from significant compensation between dispersion and charge transfer energies. However, revPBE-D3 is less accurate in describing systems where error compensation is incomplete, such as the benzene dimer. Such cases highlight the propensity for unpredictable behavior in various dispersion-corrected density functionals across a wide range of molecular systems, akin to the behavior of force fields. On the other hand, we find that SCAN-rVV10, a targeted-dispersion approach, affords significant reductions in errors associated with the lattice energies of molecular crystals, while it has limited accuracy in reproducing structural properties. Given the ubiquitous nature of noncovalent interactions and the key role of density functional theory in computational sciences, the future development of dispersion models should prioritize the faithful description of the dispersion energy, a shift that promises greater accuracy in capturing the underlying physics across diverse molecular and extended systems.
Collapse
Affiliation(s)
- Saswata Dasgupta
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Etienne Palos
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Yuanhui Pan
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
5
|
Palos E, Caruso A, Paesani F. Consistent density functional theory-based description of ion hydration through density-corrected many-body representations. J Chem Phys 2023; 159:181101. [PMID: 37947509 DOI: 10.1063/5.0174577] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 10/23/2023] [Indexed: 11/12/2023] Open
Abstract
Delocalization error constrains the accuracy of density functional theory in describing molecular interactions in ion-water systems. Using Na+ and Cl- in water as model systems, we calculate the effects of delocalization error in the SCAN functional for describing ion-water and water-water interactions in hydrated ions, and demonstrate that density-corrected SCAN (DC-SCAN) predicts n-body and interaction energies with an accuracy approaching coupled cluster theory. The performance of DC-SCAN is size-consistent, maintaining an accurate description of molecular interactions well beyond the first solvation shell. Molecular dynamics simulations at ambient conditions with many-body MB-SCAN(DC) potentials, derived from the many-body expansion, predict the solvation structure of Na+ and Cl- in quantitative agreement with reference data, while simultaneously reproducing the structure of liquid water. Beyond rationalizing the accuracy of density-corrected models of ion hydration, our findings suggest that our unified density-corrected MB formalism holds great promise for efficient DFT-based simulations of condensed-phase systems with chemical accuracy.
Collapse
Affiliation(s)
- Etienne Palos
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Alessandro Caruso
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
6
|
Stevens MJ, Rempe SLB. Binding of carboxylate and water to monovalent cations. Phys Chem Chem Phys 2023; 25:29881-29893. [PMID: 37889481 DOI: 10.1039/d3cp04200f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2023]
Abstract
The interactions of carboxylate anions with water and cations are important for a wide variety of systems, both biological and synthetic. To gain insight on properties of the local complexes, we apply density functional theory, to treat the complex electrostatic interactions, and investigate mixtures with varied numbers of carboxylate anions (acetate) and waters binding to monovalent cations, Li+, Na+ and K+. The optimal structure with overall lowest free energy contains two acetates and two waters such that the cation is four-fold coordinated, similar to structures found earlier for pure water or pure carboxylate ligands. More generally, the complexes with two acetates have the lowest free energy. In transitioning from the overall optimal state, exchanging an acetate for water has a lower free energy barrier than exchanging water for an acetate. In most cases, the carboxylates are monodentate and in the first solvation shell. As water is added to the system, hydrogen bonding between waters and carboxylate O atoms further stabilizes monodentate structures. These structures, which have strong electrostatic interactions that involve hydrogen bonds of varying strength, are significantly polarized, with ChelpG partial charges that vary substantially as the bonding geometry varies. Overall, these results emphasize the increasing importance of water as a component of binding sites as the number of ligands increases, thus affecting the preferential solvation of specific metal ions and clarifying Hofmeister effects. Finally, structural analysis correlated with free energy analysis supports the idea that binding to more than the preferred number of carboxylates under architectural constraints are a key to ion transport.
Collapse
Affiliation(s)
- Mark J Stevens
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, NM 87185, USA.
| | - Susan L B Rempe
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, NM 87185, USA.
| |
Collapse
|
7
|
Baker S, Pagotto J, Duignan TT, Page AJ. High-Throughput Aqueous Electrolyte Structure Prediction Using IonSolvR and Equivariant Graph Neural Network Potentials. J Phys Chem Lett 2023; 14:9508-9515. [PMID: 37845640 DOI: 10.1021/acs.jpclett.3c01783] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
Neural network potentials have recently emerged as an efficient and accurate tool for accelerating ab initio molecular dynamics (AIMD) in order to simulate complex condensed phases such as electrolyte solutions. Their principal limitation, however, is their requirement for sufficiently large and accurate training sets, which are often composed of Kohn-Sham density functional theory (DFT) calculations. Here we examine the feasibility of using existing density functional tight-binding (DFTB) molecular dynamics trajectory data available in the IonSolvR database in order to accelerate the training of E(3)-equivariant graph neural network potentials. We show that the solvation structure of Na+ and Cl- in aqueous NaCl solutions can be accurately reproduced with remarkably small amounts of data (i.e., 100 MD frames). We further show that these predictions can be systematically improved further via an embarrassingly parallel resampling approach.
Collapse
Affiliation(s)
- Sophie Baker
- Discipline of Chemistry, College of Engineering, Science and Environment, University of Newcastle, Callaghan, Newcastle, NSW 2308, Australia
| | - Joshua Pagotto
- School of Chemical Engineering, The University of Queensland, St Lucia, Brisbane, QLD 4072, Australia
| | - Timothy T Duignan
- School of Chemical Engineering, The University of Queensland, St Lucia, Brisbane, QLD 4072, Australia
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan, Brisbane, QLD 4111, Australia
| | - Alister J Page
- Discipline of Chemistry, College of Engineering, Science and Environment, University of Newcastle, Callaghan, Newcastle, NSW 2308, Australia
| |
Collapse
|
8
|
Belleflamme F, Hutter J. Radicals in aqueous solution: assessment of density-corrected SCAN functional. Phys Chem Chem Phys 2023; 25:20817-20836. [PMID: 37497572 DOI: 10.1039/d3cp02517a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
We study self-interaction effects in solvated and strongly-correlated cationic molecular clusters, with a focus on the solvated hydroxyl radical. To address the self-interaction issue, we apply the DC-r2SCAN method, with the auxiliary density matrix approach. Validating our method through simulations of bulk liquid water, we demonstrate that DC-r2SCAN maintains the structural accuracy of r2SCAN while effectively addressing spin density localization issues. Extending our analysis to solvated cationic molecular clusters, we find that the hemibonded motif in the [CH3S∴CH3SH]+ cluster is disrupted in the DC-r2SCAN simulation, in contrast to r2SCAN that preserves the (three-electron-two-center)-bonded motif. Similarly, for the [SH∴SH2]+ cluster, r2SCAN restores the hemibonded motif through spin leakage, while DC-r2SCAN predicts a weaker hemibond formation influenced by solvent-solute interactions. Our findings demonstrate the potential of DC-r2SCAN combined with the auxiliary density matrix method to improve electronic structure calculations, providing insights into the properties of solvated cationic molecular clusters. This work contributes to the advancement of self-interaction corrected electronic structure theory and offers a computational framework for modeling condensed phase systems with intricate correlation effects.
Collapse
Affiliation(s)
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Zurich, Switzerland.
| |
Collapse
|
9
|
Jing Z, Zhou Y, Yamaguchi T, Yoshida K, Ikeda K, Ohara K, Wang G. Hydration of Alkali Metal and Halide Ions from Static and Dynamic Viewpoints. J Phys Chem Lett 2023:6270-6277. [PMID: 37399074 DOI: 10.1021/acs.jpclett.3c01302] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2023]
Abstract
Ion hydration in aqueous solutions plays a paramount role in many fields. Despite many studies on ion hydration, the nature of ion hydration is not consistently understood at the molecular level. Combining neutron scattering (NS), wide-angle X-ray scattering (WAXS), and molecular dynamics (MD), we quantify the ionic hydration degree (hydration ability) systematically for a series of alkali metal and halide ions based on static and dynamic hydration numbers. The former is based on the orientational correlation of water molecules bound to an ion derived from the positional information from NS and WAXS. The latter is defined as the mean number of water molecules remaining in the first coordination shell of an ion over a residence time of bound water molecules around the ion from MD. The static and dynamic hydration numbers distinguish hydration from coordination and quantify the ionic hydration degree, which provides a valuable reference for understanding various phenomena in nature.
Collapse
Affiliation(s)
- Zhuanfang Jing
- Key Laboratory of Comprehensive and Highly Efficient Utilization of Salt Lake Resources, Qinghai Institute of Salt Lakes, Chinese Academy of Sciences, Xining 810008, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yongquan Zhou
- Key Laboratory of Comprehensive and Highly Efficient Utilization of Salt Lake Resources, Qinghai Institute of Salt Lakes, Chinese Academy of Sciences, Xining 810008, China
| | - Toshio Yamaguchi
- Key Laboratory of Comprehensive and Highly Efficient Utilization of Salt Lake Resources, Qinghai Institute of Salt Lakes, Chinese Academy of Sciences, Xining 810008, China
| | - Koji Yoshida
- Department of Chemistry, Faculty of Science, Fukuoka University, 8-19-1 Nanakuma, Jonan, Fukuoka 814-0180, Japan
| | - Kazutaka Ikeda
- High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan
| | - Koji Ohara
- Diffraction and Scattering Division, Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo, Hyogo 679-5198, Japan
| | - Guangguo Wang
- Key Laboratory of Comprehensive and Highly Efficient Utilization of Salt Lake Resources, Qinghai Institute of Salt Lakes, Chinese Academy of Sciences, Xining 810008, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
10
|
Kostal V, Mason PE, Martinez-Seara H, Jungwirth P. Common Cations Are Not Polarizable: Effects of Dispersion Correction on Hydration Structures from Ab Initio Molecular Dynamics. J Phys Chem Lett 2023; 14:4403-4408. [PMID: 37140439 DOI: 10.1021/acs.jpclett.3c00856] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
We employed density functional theory-based ab initio molecular dynamics simulations to examine the hydration structure of several common alkali and alkali earth metal cations. We found that the commonly used atom pairwise dispersion correction scheme D3, which assigns dispersion coefficients based on the neutral form of the atom rather than its actual oxidation state, leads to inaccuracies in the hydration structures of these cations. We evaluated this effect for lithium, sodium, potassium, and calcium and found that the inaccuracies are particularly pronounced for sodium and potassium compared to the experiment. To remedy this issue, we propose disabling the D3 correction specifically for all cation-including pairs, which leads to a much better agreement with experimental data.
Collapse
Affiliation(s)
- Vojtech Kostal
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, 166 10 Prague 6, Czech Republic
| | - Philip E Mason
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, 166 10 Prague 6, Czech Republic
| | - Hector Martinez-Seara
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, 166 10 Prague 6, Czech Republic
| | - Pavel Jungwirth
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 2, 166 10 Prague 6, Czech Republic
| |
Collapse
|
11
|
McElhany SJ, Summers TJ, Shiery RC, Cantu DC. Analysis of the First Ion Coordination Sphere: A Toolkit to Analyze the Coordination Sphere of Ions. J Chem Inf Model 2023; 63:2699-2706. [PMID: 37083437 DOI: 10.1021/acs.jcim.3c00294] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/22/2023]
Abstract
Rapid and accurate approaches to characterizing the coordination structure of an ion are important for designing ligands and quantifying structure-property trends. Here, we introduce AFICS (Analysis of the First Ion Coordination Sphere), a tool written in Python 3 for analyzing the structural and geometric features of the first coordination sphere of an ion over the course of molecular dynamics simulations. The principal feature of AFICS is its ability to quantify the distortion a coordination geometry undergoes compared to uniform polyhedra. This work applies the toolkit to analyze molecular dynamics simulations of the well-defined coordination structure of aqueous Cr3+ along with the more ambiguous structure of aqueous Eu3+ chelated to ethylenediaminetetraacetic acid. The tool is targeted for analyzing ions with fluxional or irregular coordination structures (e.g., solution structures of f-block elements) but is generalized such that it may be applied to other systems.
Collapse
Affiliation(s)
- Stuart J McElhany
- Department of Chemical and Materials Engineering, University of Nevada, Reno, Reno, Nevada 89557, United States
| | - Thomas J Summers
- Department of Chemical and Materials Engineering, University of Nevada, Reno, Reno, Nevada 89557, United States
| | - Richard C Shiery
- Department of Chemical and Materials Engineering, University of Nevada, Reno, Reno, Nevada 89557, United States
| | - David C Cantu
- Department of Chemical and Materials Engineering, University of Nevada, Reno, Reno, Nevada 89557, United States
| |
Collapse
|
12
|
Gao A, Remsing RC, Weeks JD. Local Molecular Field Theory for Coulomb Interactions in Aqueous Solutions. J Phys Chem B 2023; 127:809-821. [PMID: 36669139 DOI: 10.1021/acs.jpcb.2c06988] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Coulomb interactions play a crucial role in a wide array of processes in aqueous solutions but present conceptual and computational challenges to both theory and simulations. We review recent developments in an approach addressing these challenges─local molecular field (LMF) theory. LMF theory exploits an exact and physically suggestive separation of intermolecular Coulomb interactions into strong short-range and uniformly slowly varying long-range components. This allows us to accurately determine the averaged effects of the long-range components on the short-range structure using effective single particle fields and analytical corrections, greatly reducing the need for complex lattice summation techniques used in most standard approaches. The simplest use of these ideas in aqueous solutions leads to the short solvent (SS) model, where both solvent-solvent and solute-solvent Coulomb interactions have only short-range components. Here we use the SS model to give a simple description of pairing of nucleobases and biologically relevant ions in water.
Collapse
Affiliation(s)
- Ang Gao
- Department of Physics, Beijing University of Posts and Telecommunications, Beijing, China 100876
| | - Richard C Remsing
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - John D Weeks
- Institute for Physical Science and Technology and Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
13
|
Zhou K, Qian C, Liu Y. Quantifying the Structure of Water and Hydrated Monovalent Ions by Density Functional Theory-Based Molecular Dynamics. J Phys Chem B 2022; 126:10471-10480. [PMID: 36451081 DOI: 10.1021/acs.jpcb.2c05330] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
The accurate description of the structures of water and hydrated ions is important in electrochemical desalination, ion separation, and supercapacitors. In this work, we present an ab initio atomistic simulation-based study to explore the structure of water and hydrated monovalent ions (Li+, Na+, K+, Rb+, F-, and Cl-) at ambient conditions using generalized gradient approximation (GGA)-based methods with and without van der Waals correction (PBE, PBE + D3, and revPBE + D3) and recently developed strongly constrained and appropriately normed (SCAN) meta-GGA. We find that both revPBE + D3 and SCAN can well capture the structure of bulk water with +30 K artificial high temperature in contrast to overstructuring water using PBE and PBE + D3. However, being the same as PBE + D3, revPBE + D3 overestimates the structure of the hydration shell, especially for monovalent cations. Surprisingly, SCAN can well match the experimental results of hydrated monovalent ions. Detailed structure analyzes of entropy reveal that the hydration shell under the level of PBE + D3 and revPBE + D3 is more disordered and looser than SCAN. The successful prediction of the flexible SCAN functional could facilitate the exploration of complex ionic processes in the aqueous phase, the interactions of hydrated ions with surfaces, and solvation states in nanopores at an accurate, efficient, predictive, and ab initio level.
Collapse
Affiliation(s)
- Ke Zhou
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou215006, China.,Laboratory for Multiscale Mechanics and Medical Science, SV LAB, School of Aerospace, Xi'an Jiaotong University, Xi'an710049, China
| | - Chen Qian
- Department of Mechanical Engineering, Zhejiang University, Hangzhou310058, China
| | - Yilun Liu
- Laboratory for Multiscale Mechanics and Medical Science, SV LAB, School of Aerospace, Xi'an Jiaotong University, Xi'an710049, China
| |
Collapse
|
14
|
Zhuang D, Riera M, Zhou R, Deary A, Paesani F. Hydration Structure of Na + and K + Ions in Solution Predicted by Data-Driven Many-Body Potentials. J Phys Chem B 2022; 126:9349-9360. [PMID: 36326071 DOI: 10.1021/acs.jpcb.2c05674] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
The hydration structure of Na+ and K+ ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na+ and K+ in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion-water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.
Collapse
Affiliation(s)
- Debbie Zhuang
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Ruihan Zhou
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Alexander Deary
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States.,Materials Science and Engineering, University of California San Diego, La Jolla, California92093, United States.,San Diego Supercomputer Center, University of California San Diego, La Jolla, California92093, United States
| |
Collapse
|
15
|
Xi C, Zheng F, Gao G, Song Z, Zhang B, Dong C, Du XW, Wang LW. Ion Solvation Free Energy Calculation Based on Ab Initio Molecular Dynamics Using a Hybrid Solvent Model. J Chem Theory Comput 2022; 18:6878-6891. [PMID: 36253911 DOI: 10.1021/acs.jctc.1c01298] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Free energy calculation of small molecules or ion species in aqueous solvent is one of the most important problems in electrochemistry study. Although there are many previous approaches to calculate such free energies, they are based on ab initio methods and suffer from various limitations and approximations. In the current work, we developed a hybrid approach based on ab initio molecular dynamics (AIMD) simulations to calculate the ion solvation energy. In this approach, a small water cluster surrounding the central ion is used, and implicit solvent model is applied outside the water cluster. A dynamic potential well is used during AIMD to keep the water cluster together. Quasi-harmonic approximation is used to calculate the entropy contribution, while the total energy average is used to calculate the enthalpy term. The obtained solvation voltages of the bulk metal agree with experiments within 0.3 eV, and the simulation results for the solvation energies of gaseous ions are close to the experimental observations. Besides the free energies, radial pair distribution functions and coordination numbers of hydrated cations are also obtained. The remaining challenges of this method are also discussed.
Collapse
Affiliation(s)
- Cong Xi
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States.,Institute of New Energy Materials, School of Materials Science and Engineering, Tianjin University, Tianjin30072, People's Republic of China
| | - Fan Zheng
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Guoping Gao
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Zhigang Song
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Buyu Zhang
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Cunku Dong
- Institute of New Energy Materials, School of Materials Science and Engineering, Tianjin University, Tianjin30072, People's Republic of China
| | - Xi-Wen Du
- Institute of New Energy Materials, School of Materials Science and Engineering, Tianjin University, Tianjin30072, People's Republic of China
| | - Lin-Wang Wang
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| |
Collapse
|
16
|
Rehl B, Ma E, Parshotam S, DeWalt-Kerian EL, Liu T, Geiger FM, Gibbs JM. Water Structure in the Electrical Double Layer and the Contributions to the Total Interfacial Potential at Different Surface Charge Densities. J Am Chem Soc 2022; 144:16338-16349. [PMID: 36042195 DOI: 10.1021/jacs.2c01830] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The electric double layer governs the processes of all charged surfaces in aqueous solutions; however, elucidating the structure of the water molecules is challenging for even the most advanced spectroscopic techniques. Here, we present the individual Stern layer and diffuse layer OH stretching spectra at the silica/water interface in the presence of NaCl over a wide pH range using a combination of vibrational sum frequency generation spectroscopy, heterodyned second harmonic generation, and streaming potential measurements. We find that the Stern layer water molecules and diffuse layer water molecules respond differently to pH changes: unlike the diffuse layer, whose water molecules remain net-oriented in one direction, water molecules in the Stern layer flip their net orientation as the solution pH is reduced from basic to acidic. We obtain an experimental estimate of the non-Gouy-Chapman (Stern) potential contribution to the total potential drop across the insulator/electrolyte interface and discuss it in the context of dipolar, quadrupolar, and higher order potential contributions that vary with the observed changes in the net orientation of water in the Stern layer. Our findings show that a purely Gouy-Chapman (Stern) view is insufficient to accurately describe the electrical double layer of aqueous interfaces.
Collapse
Affiliation(s)
- Benjamin Rehl
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| | - Emily Ma
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
| | - Shyam Parshotam
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| | - Emma L DeWalt-Kerian
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| | - Tianli Liu
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| | - Franz M Geiger
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
| | - Julianne M Gibbs
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| |
Collapse
|
17
|
Yamaguchi A, Nagata K, Kobayashi K, Tanaka K, Kobayashi T, Tanida H, Shimojo K, Sekiguchi T, Kaneta Y, Matsuda S, Yokoyama K, Yaita T, Yoshimura T, Okumura M, Takahashi Y. EXAFS spectroscopy measurements and ab initio molecular dynamics simulations reveal the hydration structure of the radium(II) ion. iScience 2022; 25:104763. [PMID: 35992079 PMCID: PMC9386089 DOI: 10.1016/j.isci.2022.104763] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Revised: 06/27/2022] [Accepted: 07/09/2022] [Indexed: 12/02/2022] Open
Abstract
Radium is refocused from the viewpoint of an environmental pollutant and cancer therapy using alpha particles, where it mainly exists as a hydrated ion. We investigated the radium hydration structure and the dynamics of water molecules by extended X-ray absorption fine structure (EXAFS) spectroscopy and ab initio molecular dynamics (AIMD) simulation. The EXAFS experiment showed that the coordination number and average distance between radium ion and the oxygen atoms in the first hydration shell are 9.2 ± 1.9 and 2.87 ± 0.06 Å, respectively. They are consistent with those obtained from the AIMD simulations, 8.4 and 2.88 Å. The AIMD simulations also revealed that the water molecules in the first hydration shell of radium are less structured and more mobile than those of barium, which is an analogous element of radium. Our results indicate that radium can be more labile than barium in terms of interactions with water. Extended X-ray absorption fine structure (EXAFS) measurement revealed the hydration structure of radium ion Ab initio molecular dynamics (AIMD) simulation brought consistent results AIMD revealed the structural and dynamic properties of the water molecules The hydration structure of radium ion is more labile than that of barium ion
Collapse
|
18
|
Herrero C, Pauletti M, Tocci G, Iannuzzi M, Joly L. Connection between water's dynamical and structural properties: Insights from ab initio simulations. Proc Natl Acad Sci U S A 2022; 119:e2121641119. [PMID: 35588447 PMCID: PMC9173753 DOI: 10.1073/pnas.2121641119] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 04/12/2022] [Indexed: 01/25/2023] Open
Abstract
SignificanceFirst-principles calculations, which explicitly account for the electronic structure of matter, can shed light on the molecular structure and dynamics of water in its supercooled state. In this work, we use density functional theory, which relies on a functional to describe electronic exchange and correlations, to evaluate which functional best describes the temperature evolution of bulk water transport coefficients. We also assess the validity of the Stokes-Einstein relation for all the functionals in the temperature range studied, and explore the link between structure and dynamics. Based on these results, we show how transport coefficients can be computed from structural descriptors, which require shorter simulation times to converge, and we point toward strategies to develop better functionals.
Collapse
Affiliation(s)
- Cecilia Herrero
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Michela Pauletti
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Gabriele Tocci
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Marcella Iannuzzi
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Laurent Joly
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
- Institut Universitaire de France (IUF), 75005 Paris, France
| |
Collapse
|
19
|
Beverley KM, Shahi PK, Kabra M, Zhao Q, Heyrman J, Steffen J, Pattnaik BR. Kir7.1 disease mutant T153I within the inner pore affects K+ conduction. Am J Physiol Cell Physiol 2022; 323:C56-C68. [PMID: 35584325 DOI: 10.1152/ajpcell.00093.2022] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Inward-rectifier potassium channel 7.1 (Kir7.1) is present in polarized epithelium, including the RPE. A single amino acid change at position 153 in the KCNJ13 gene, a substitution of threonine to isoleucine in Kir7.1 protein, causes blindness. We hypothesized that the disease caused by this single amino acid substitution within the transmembrane protein domain could alter the translation, localization, or ion transport properties. We assessed the effects of amino acid side-chain length, arrangement, and polarity on channel structure and function. We showed that the T153I mutation yielded a full-length protein localized to the cell membrane. Whole-cell patch-clamp recordings and chord conductance analyses revealed that the T153I mutant channel had negligible K+ conductance and failed to hyperpolarize the membrane potential. However, the mutant channel exhibited enhanced inward current when Rb+ was used as a charge carrier, suggesting that an inner pore had formed, and the channel was dysfunctional. Substituting with a polar, non-polar, or short side-chain amino acid did not affect the localization of the protein. Still, it had an altered channel function due to differences in pore radius. Polar side chains (cysteine and serine) with inner pore radii comparable to wildtype exhibited normal inward K+ conductance. Short side-chains (glycine and alanine) produced a channel with wider than expected inner pore size and lacked the biophysical characteristics of the wildtype channel. Leucine substitution produced results similar to the T153I mutant channel. This study provides direct electrophysiological evidence for the structure and function of the Kir7.1 channel's narrow inner pore in regulating conductance.
Collapse
Affiliation(s)
- Katie M Beverley
- Endocrinology and Reproductive Physiology Graduate Program, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States.,Department of Pediatrics, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States.,McPherson Eye Research Institute, University of Wisconsin, Madison, WI, United States
| | - Pawan K Shahi
- Department of Pediatrics, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States.,McPherson Eye Research Institute, University of Wisconsin, Madison, WI, United States
| | - Meha Kabra
- Department of Pediatrics, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States.,McPherson Eye Research Institute, University of Wisconsin, Madison, WI, United States
| | - Qianqian Zhao
- Department of Biostatistics and Medical Informatics, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States
| | - Joseph Heyrman
- Department of Pediatrics, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States
| | - Jack Steffen
- Department of Pediatrics, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States
| | - Bikash R Pattnaik
- Endocrinology and Reproductive Physiology Graduate Program, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States.,Department of Pediatrics, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States.,McPherson Eye Research Institute, University of Wisconsin, Madison, WI, United States.,Department of Ophthalmology and Visual Sciences, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States
| |
Collapse
|
20
|
Khanal R, Irle S. Quantum chemical investigation of the effect of alkali metal ions on the dynamic structure of water in aqueous solutions. RSC Adv 2022; 12:25500-25510. [PMID: 36275866 PMCID: PMC9480497 DOI: 10.1039/d2ra04563j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2022] [Accepted: 08/26/2022] [Indexed: 11/21/2022] Open
Abstract
Cation size primarily governs dynamic correlations in aqueous solutions. However, drawing a clear line between chaotropic and cosmotrope ions is complicated by the fact that all cations delay dynamic correlation decay.
Collapse
Affiliation(s)
- Rabi Khanal
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee, 37831, USA
| | - Stephan Irle
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee, 37831, USA
| |
Collapse
|
21
|
Confalonieri G, Vezzalini G, Quattrini F, Quartieri S, Dejoie C, Arletti R. Ce-exchange capacity of zeolite L in different cationic forms: a structural investigation. J Appl Crystallogr 2021. [DOI: 10.1107/s1600576721010827] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
Cerium exchange by microporous materials, such as zeolites, has important applications in different fields, for example, rare earth element recovery from waste or catalytic processes. This work investigated the Ce-exchange capacity of zeolite L in three different cationic forms (the as-synthesized K form and Na- and NH4-exchanged ones) from a highly concentrated solution. Chemical analyses and structural investigations allowed determination of the mechanisms involved in the exchanges and give new insights into the interactions occurring between the cations and the zeolite framework. Different cation sites are involved: (i) K present in the original LTL in the cancrinite cage (site KB) cannot be exchanged; (ii) the cations in KD (in the 12-membered ring channel) are always exchanged; while (iii) site KC (in the eight-membered ring channel) is involved only when K+ is substituted by NH4
+, thus promoting a higher exchange rate for NH4
+ → K+ than for Na+ → K+. In the Ce-exchanged samples, a new site occupied by Ce appears in the centre of the main channel, accompanied by an increase in the number of and a rearrangement of H2O molecules. In terms of Ce exchange, the three cationic forms behave similarly, from both the chemical and structural point of view (exchanged Ce ranges from 38 to 42% of the pristine cation amount). Beyond the intrinsic structural properties of the zeolite L framework, the Ce exchange seems thus also governed by the water coordination sphere of the cation. Complete Ce recovery from zeolite pores was achieved.
Collapse
|
22
|
Duignan TT, Zhao XS. Prediction of the Osmotic/Activity Coefficients of Alkali Hydroxide Electrolytes. Ind Eng Chem Res 2021. [DOI: 10.1021/acs.iecr.1c02950] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Timothy T. Duignan
- School of Chemical Engineering, The University of Queensland, St Lucia, Brisbane 4072, Australia
| | - X. S. Zhao
- School of Chemical Engineering, The University of Queensland, St Lucia, Brisbane 4072, Australia
| |
Collapse
|
23
|
Yamaguchi A, Kobayashi K, Takahashi Y, Machida M, Okumura M. Hydration structures of barium ions: Ab initio molecular dynamics simulations using the SCAN meta-GGA density functional and EXAFS spectroscopy studies. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138945] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
24
|
Duignan TT, Kathmann SM, Schenter GK, Mundy CJ. Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes. Acc Chem Res 2021; 54:2833-2843. [PMID: 34137593 DOI: 10.1021/acs.accounts.1c00107] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
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.
Collapse
Affiliation(s)
- Timothy T. Duignan
- School of Chemical Engineering, The University of Queensland, St Lucia, Brisbane 4072, Australia
| | - Shawn M. Kathmann
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Gregory K. Schenter
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Christopher J. Mundy
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
- Affiliate Professor, Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
25
|
Osti NC, Thapaliya BP, Dai S, Tyagi M, Mamontov E. Strong Enhancement of Nanoconfined Water Mobility by a Structure Breaking Salt. J Phys Chem Lett 2021; 12:4038-4044. [PMID: 33881871 DOI: 10.1021/acs.jpclett.1c00461] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
For the majority of the water present on earth, the two most important factors influencing its behavior are confinement, in either inorganic or organic matrixes, and the presence of solutes. Here, we investigate the effect of confinement in 3 nm pores on water diffusivity in aqueous solutions with archetypical solutes, a structure making (kosmotrope) NaCl and a structure breaking (chaotrope) KCl, up to 1.0 M in concentration. The water diffusivity in bulk aqueous solutions in such a concentration range is known to decrease very slightly in the presence of NaCl and increase very slightly in the presence of KCl. However, here we observe the water diffusivity in confined H2O-KCl increases by a factor of 2 compared to the pure water diffusivity in the same confinement. This unusually strong cumulative effect of confinement and a structure breaking additive may have profound implications for the mobility and transport of aqueous species in nature.
Collapse
Affiliation(s)
- Naresh C Osti
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Bishnu Prasad Thapaliya
- Chemical Science Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37916, United States
| | - Sheng Dai
- Chemical Science Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37916, United States
| | - Madhusudan Tyagi
- NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
- Department of Materials Science, University of Maryland, College Park, Maryland 20742, United States
| | - Eugene Mamontov
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| |
Collapse
|
26
|
Zhan C, Sun Y, Aydin F, Wang YM, Pham TA. Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel: A first-principles study. J Chem Phys 2021; 154:164706. [PMID: 33940836 DOI: 10.1063/5.0047554] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS2 nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS2 channel can be significantly altered. We found that larger cations with weaker hydration energy (K+, Rb+, and Cs+) exhibited a distinctive CN decrease under confinement, while smaller cations (Li+ and Na+) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K+, Rb+, and Cs+) in MoS2 showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS2 channel.
Collapse
Affiliation(s)
- Cheng Zhan
- Material Science Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Yangyunli Sun
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Fikret Aydin
- Material Science Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Y Morris Wang
- Department of Materials Science and Engineering, University of California, Log Angeles, California 90095, USA
| | - Tuan Anh Pham
- Material Science Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| |
Collapse
|
27
|
Heindel JP, Xantheas SS. The Many-Body Expansion for Aqueous Systems Revisited: II. Alkali Metal and Halide Ion-Water Interactions. J Chem Theory Comput 2021; 17:2200-2216. [PMID: 33709708 DOI: 10.1021/acs.jctc.0c01309] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
We present a detailed study of the many-body expansion (MBE) for alkali metal and halide ion-water interactions and quantify the effect of these ions on the strength of the surrounding aqueous hydrogen bonding environment. Building on our previous work on neutral water clusters [J. P. Heindel and S. S. Xantheas, J. Chem. Theor. Comput. 16 (11), 6843-6855 (2020)], we carry out the MBE for the alkali metal and halide ion-water clusters, Z+/-(H2O)9, where Z = Li+, Na+, K+, Rb+, Cs+, F-, Cl-, Br-, and I- and compare them with the results for a pure water cluster with the same number of "bodies", viz., (H2O)10. The 2-B ion-water (I-W) interaction accounts for a larger percentage of the total cluster binding energy compared to a pure water cluster of the same size, with the total 3-B term being smaller and of opposite sign (repulsive), whereas higher order terms are essentially negligible. The same oscillating behavior around zero for the MBE terms higher than the 5-B with a basis set that was reported for water clusters is also observed for the ion-water clusters considered here, with the basis set superposition error (BSSE) corrections amending this as in the water cluster case. A remarkable, linear anticorrelation between the total 2-B (I-W), the total 2-B (W-W), and also the 3-B (W-W-W) interactions is found, quantifying the effect of the different ions in disrupting and altering (weakening) the neighboring hydrogen bonded water network: stronger (I-W) interactions result in weaker (W-W) interactions. Additional linear correlations across the two series of alkali metals and halide ions were found between the 3-B (I-W-W) and the 2-B (I-W) as well as between the 3-B (I-W-W) and the 3-B (W-W-W) interactions, suggesting the existence of previously unrealized underlying physics governing these 2-B intermolecular and 3-B collective interactions. Our results further suggest a universal behavior of the two different families of ions (alkali metals and halides) for both the correlations of the various components of the total binding energies and the estimate of the 2-B BSSE correction, which is reported to follow a common profile for ion-water and water-water interactions when cast in terms of reduced distances and energies of the respective dimers. We expect the current results that quantify the interplay between ion-water and water-water interactions in aqueous clusters to impact the development of classical, ab initio-based force fields for monatomic ion solvation, whereas the insights into the nature of the BSSE to be critical in future ab initio-based, many-body molecular dynamics studies.
Collapse
Affiliation(s)
- Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States.,Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| |
Collapse
|
28
|
Fuentes-Monteverde JC, Becker S, Rezaei-Ghaleh N. Biomolecular phase separation through the lens of sodium-23 NMR. Protein Sci 2020; 30:1315-1325. [PMID: 33314347 PMCID: PMC8197435 DOI: 10.1002/pro.4010] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2020] [Revised: 12/04/2020] [Accepted: 12/04/2020] [Indexed: 01/10/2023]
Abstract
Phase separation is a fundamental physicochemical process underlying the spatial arrangement and coordination of cellular events. Detailed characterization of biomolecular phase separation requires experimental access to the internal environment of dilute and especially condensed phases at high resolution. In this study, we take advantage from the ubiquitous presence of sodium ions in biomolecular samples and present the potentials of 23Na NMR as a proxy to report the internal fluidity of biomolecular condensed phases. After establishing the temperature and viscosity dependence of 23Na NMR relaxation rates and translational diffusion coefficient, we demonstrate that 23Na NMR probes of rotational and translational mobility of sodium ions are capable of capturing the increasing levels of confinement in agarose gels in dependence of agarose concentration. The 23Na NMR approach is then applied to a gel‐forming phenylalanine‐glycine (FG)‐containing peptide, part of the nuclear pore complex involved in controlling the traffic between cytoplasm and cell nucleus. It is shown that the 23Na NMR together with the 17O NMR provide a detailed picture of the sodium ion and water mobility within the interior of the FG peptide hydrogel. As another example, we study phase separation in water‐triethylamine (TEA) mixture and provide evidence for the presence of multiple microscopic environments within the TEA‐rich phase. Our results highlight the potentials of 23Na NMR in combination with 17O NMR in studying biological phase separation, in particular with regards to the molecular properties of biomolecular condensates and their regulation through various physico‐ and biochemical factors.
Collapse
Affiliation(s)
| | - Stefan Becker
- Department for NMR-based Structural Biology, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Nasrollah Rezaei-Ghaleh
- Department for NMR-based Structural Biology, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany.,Department of Neurology, University Medical Center Göttingen, Göttingen, Germany
| |
Collapse
|
29
|
A chemometric investigation on the influence of the nature and concentration of supporting electrolyte on charging currents in electrochemistry. J Electroanal Chem (Lausanne) 2020. [DOI: 10.1016/j.jelechem.2020.114296] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
30
|
Wang X, Toroz D, Kim S, Clegg SL, Park GS, Di Tommaso D. Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions. Phys Chem Chem Phys 2020; 22:16301-16313. [PMID: 32647838 DOI: 10.1039/d0cp01957g] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell.
Collapse
Affiliation(s)
- Xiangwen Wang
- School of Biological and Chemical Sciences, Materials Research Institute, Thomas Young Centre, Queen Mary University of London, Mile End Road, E1 4NS, London, UK.
| | - Dimitrios Toroz
- School of Biological and Chemical Sciences, Materials Research Institute, Thomas Young Centre, Queen Mary University of London, Mile End Road, E1 4NS, London, UK.
| | - Seonmyeong Kim
- Center for THz-driven Biological Systems, Department of Physics and Astronomy, Seoul National University, Seoul, 08826, Republic of Korea and Advanced Institutes of Convergence Technology, Seoul National University, Suwon-Si, Gyeonggi-do 16229, Republic of Korea
| | - Simon L Clegg
- School of Environmental Sciences, University of East Anglia, Norwich, NR4 7TJ, UK
| | - Gun-Sik Park
- Center for THz-driven Biological Systems, Department of Physics and Astronomy, Seoul National University, Seoul, 08826, Republic of Korea and Advanced Institutes of Convergence Technology, Seoul National University, Suwon-Si, Gyeonggi-do 16229, Republic of Korea
| | - Devis Di Tommaso
- School of Biological and Chemical Sciences, Materials Research Institute, Thomas Young Centre, Queen Mary University of London, Mile End Road, E1 4NS, London, UK.
| |
Collapse
|
31
|
Duignan TT, Mundy CJ, Schenter GK, Zhao XS. Method for Accurately Predicting Solvation Structure. J Chem Theory Comput 2020; 16:5401-5409. [DOI: 10.1021/acs.jctc.0c00300] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Timothy T. Duignan
- School of Chemical Engineering, The University of Queensland, St Lucia, Brisbane 4072, Australia
| | - Christopher J. Mundy
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
- Affiliate Professor, Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| | - Gregory K. Schenter
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - X. S. Zhao
- School of Chemical Engineering, The University of Queensland, St Lucia, Brisbane 4072, Australia
| |
Collapse
|
32
|
Duignan TT, Zhao XS. The Born model can accurately describe electrostatic ion solvation. Phys Chem Chem Phys 2020; 22:25126-25135. [DOI: 10.1039/d0cp04148c] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The solvation free energies of ions in water are consistent with the Born linear response model if the centre on which the ion–water repulsion force acts is moved from the oxygen atom towards the hydrogens.
Collapse
Affiliation(s)
- Timothy T. Duignan
- School of Chemical Engineering
- The University of Queensland
- St Lucia
- Australia
| | - X. S. Zhao
- School of Chemical Engineering
- The University of Queensland
- St Lucia
- Australia
| |
Collapse
|
33
|
Aydin F, Zhan C, Ritt C, Epsztein R, Elimelech M, Schwegler E, Pham TA. Similarities and differences between potassium and ammonium ions in liquid water: a first-principles study. Phys Chem Chem Phys 2020; 22:2540-2548. [DOI: 10.1039/c9cp06163k] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Understanding ion solvation in liquid water is critical in optimizing materials for a wide variety of emerging technologies, including water desalination and purification.
Collapse
Affiliation(s)
- Fikret Aydin
- Lawrence Livermore National Laboratory
- Livermore
- USA
| | - Cheng Zhan
- Lawrence Livermore National Laboratory
- Livermore
- USA
| | - Cody Ritt
- Department of Chemical and Environmental Engineering
- Yale University
- New Haven
- USA
| | - Razi Epsztein
- Department of Chemical and Environmental Engineering
- Yale University
- New Haven
- USA
- Faculty of Civil and Environmental Engineering
| | - Menachem Elimelech
- Department of Chemical and Environmental Engineering
- Yale University
- New Haven
- USA
| | | | | |
Collapse
|