1
|
Rana R, Ali SM, Maity DK. Structure and dynamics of the Li + ion in water, methanol and acetonitrile solvents: ab initio molecular dynamics simulations. Phys Chem Chem Phys 2023; 25:31382-31395. [PMID: 37961866 DOI: 10.1039/d3cp04403c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2023]
Abstract
Fundamental understanding of the structure and dynamics of the Li+ ion in solution is of utmost importance in different fields of science and technology, especially in the field of ion batteries. In view of this, ab initio molecular dynamics (AIMD) simulations of the LiCl salt in water, methanol and acetonitrile were performed to elucidate structural parameters such as radial distribution function and coordination number, and dynamical properties like diffusion coefficient, limiting ion conductivity and hydrogen bond correlation function. In the present AIMD simulation, one LiCl in water is equivalent to 0.8 M, which is close to the concentration of the lithium salt used in the Li-ion battery. The first sphere of coordination number of the Li+ ion was reaffirmed to be 4. The radial distribution function for different pairs of atoms is seen to be in good agreement with the experimental results. The calculated potential of mean force indicates the stronger interaction of the Li+ ion with methanol over water followed by acetonitrile. The dynamical parameters convey quite high diffusion and limiting ionic conductivity of the Li+ ion in acetonitrile compared to that in water and methanol which has been attributed to the transport of the Li-Cl ion pair in a non-dissociated form in acetonitrile. The AIMD results were found to be in accordance with the experimental findings, i.e. the limiting ion conductivity was found to follow the order acetonitrile > methanol > water. This study shows the importance of atomistic level simulations in evaluating the structural and dynamical parameters and in implementing the results for predicting and synthesizing better next generation solvents for lithium ion batteries (LIBs).
Collapse
Affiliation(s)
- Reman Rana
- Homi Bhabha National Institute, Anushaktinagar, Mumbai 400094, India.
| | - Sk Musharaf Ali
- Homi Bhabha National Institute, Anushaktinagar, Mumbai 400094, India.
- Chemical Engineering Division, Bhabha Atomic Research Centre, Mumbai-400085, India
| | - Dilip K Maity
- Homi Bhabha National Institute, Anushaktinagar, Mumbai 400094, India.
| |
Collapse
|
2
|
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
|
3
|
Tran B, Cai Y, Janik MJ, Milner ST. Hydrogen Bond Thermodynamics in Aqueous Acid Solutions: A Combined DFT and Classical Force-Field Approach. J Phys Chem A 2022; 126:7382-7398. [PMID: 36190836 DOI: 10.1021/acs.jpca.2c04124] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The thermodynamics of hydrogen bonds in aqueous and acidic solutions significantly impacts the kinetics and thermodynamics of acid reaction chemistry. We utilize in this work a multiscale approach, combining density functional theory (DFT) with classical molecular dynamics (MD) to model hydrogen bond thermodynamics in an acidic solution. Using thermodynamic cycles, we split the solution phase free energy into its gas phase counterpart plus solvation free energies. We validate this DFT/MD approach by calculating the aqueous phase hydrogen bond free energy between two water molecules (H2O-···-H2O), the free energy to transform an H3O+ cation into an H5O2+ cation, and the hydrogen bond free energy of protonated water clusters (H3O+-···-H2O and H5O2+-···-H2O). The computed equilibrium hydrogen bond free energy of H2O-···-H2O is remarkably accurate, especially considering the large individual contributions to the thermodynamic cycle. Turning to cations, we find the ion to be more stable than H3O+ by roughly 1-2 kBT. This small free energy difference allows for thermal fluctuation between the two idealized motifs, consistent with spectroscopic and simulation studies. Lastly, hydrogen bonding free energies between either H+ cation and H2O in solution were found to be stronger than between two H2O, though much less so than in vacuum because of dielectric screening in solution. Altogether, our results suggest the DFT/MD approach is promising for application in modeling hydrogen bonding and proton transfer thermodynamics in condensed phases.
Collapse
Affiliation(s)
- Bolton Tran
- Department of Chemical Engineering, Pennsylvania State University, University Park, Pennsylvania16801, United States
| | - Yusheng Cai
- Department of Chemical & Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania19104, United States
| | - Michael J Janik
- Department of Chemical Engineering, Pennsylvania State University, University Park, Pennsylvania16801, United States
| | - Scott T Milner
- Department of Chemical Engineering, Pennsylvania State University, University Park, Pennsylvania16801, United States
| |
Collapse
|
4
|
Phan LX, Lynch CI, Crain J, Sansom MS, Tucker SJ. Influence of effective polarization on ion and water interactions within a biomimetic nanopore. Biophys J 2022; 121:2014-2026. [DOI: 10.1016/j.bpj.2022.05.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Revised: 03/25/2022] [Accepted: 05/04/2022] [Indexed: 11/16/2022] Open
|
5
|
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
|
6
|
Yang X, Zhuang Y, Zhu J, Le J, Cheng J. Recent progress on multiscale modeling of electrochemistry. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1559] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Xiao‐Hui Yang
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Yong‐Bin Zhuang
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Jia‐Xin Zhu
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Jia‐Bo Le
- Ningbo Institute of Materials Technology and Engineering Chinese Academy of Sciences Ningbo China
| | - Jun Cheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| |
Collapse
|
7
|
Gentry BM, Choi TH, Belfield WS, Keith JA. Computational predictions of metal-macrocycle stability constants require accurate treatments of local solvent and pH effects. Phys Chem Chem Phys 2021; 23:9189-9197. [PMID: 33885118 DOI: 10.1039/d1cp00611h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Rational design of molecular chelating agents requires a detailed understanding of physicochemical ligand-metal interactions in solvent phase. Computational quantum chemistry methods should be able to provide this, but computational reports have shown poor accuracy when determining absolute binding constants for many chelating molecules. To understand why, we compare and benchmark static- and dynamics-based computational procedures for a range of monovalent and divalent cations binding to a conventional cryptand molecule: 2.2.2-cryptand ([2.2.2]). The benchmarking comparison shows that dynamics simulations using standard OPLS-AA classical potentials can reasonably predict binding constants for monovalent cations, but these procedures fail for divalent cations. We also consider computationally efficient static procedure using Kohn-Sham density functional theory (DFT) and cluster-continuum modeling that accounts for local microsolvation and pH effects. This approach accurately predicts binding energies for monovalent and divalent cations with an average error of 3.2 kcal mol-1 compared to experiment. This static procedure thus should be useful for future molecular screening efforts, and high absolute errors in the literature may be due to inadequate modeling of local solvent and pH effects.
Collapse
Affiliation(s)
- Brian M Gentry
- Department of Chemical Engineering, University of Pittsburgh, 3700 O'Hara Street, Pittsburgh, PA 15261, USA.
| | | | | | | |
Collapse
|
8
|
Jinnouchi R, Karsai F, Verdi C, Kresse G. First-principles hydration free energies of oxygenated species at water-platinum interfaces. J Chem Phys 2021; 154:094107. [PMID: 33685177 DOI: 10.1063/5.0036097] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The hydration free energy of atoms and molecules adsorbed at liquid-solid interfaces strongly influences the stability and reactivity of solid surfaces. However, its evaluation is challenging in both experiments and theories. In this work, a machine learning aided molecular dynamics method is proposed and applied to oxygen atoms and hydroxyl groups adsorbed on Pt(111) and Pt(100) surfaces in water. The proposed method adopts thermodynamic integration with respect to a coupling parameter specifying a path from well-defined non-interacting species to the fully interacting ones. The atomistic interactions are described by a machine-learned inter-atomic potential trained on first-principles data. The free energy calculated by the machine-learned potential is further corrected by using thermodynamic perturbation theory to provide the first-principles free energy. The calculated hydration free energies indicate that only the hydroxyl group adsorbed on the Pt(111) surface attains a hydration stabilization. The observed trend is attributed to differences in the adsorption site and surface morphology.
Collapse
Affiliation(s)
| | - Ferenc Karsai
- VASP Software GmbH, Sensengasse 8/16, 1090 Vienna, Austria
| | - Carla Verdi
- Computational Materials Physics, Faculty of Physics, University of Vienna, Sensengasse 8/12, 1090 Vienna, Austria
| | - Georg Kresse
- VASP Software GmbH, Sensengasse 8/16, 1090 Vienna, Austria
| |
Collapse
|
9
|
Itkis D, Cavallo L, Yashina LV, Minenkov Y. Ambiguities in solvation free energies from cluster-continuum quasichemical theory: lithium cation in protic and aprotic solvents. Phys Chem Chem Phys 2021; 23:16077-16088. [PMID: 34291782 DOI: 10.1039/d1cp01454d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Gibbs free energies for Li+ solvation in water, methanol, acetonitrile, DMSO, dimethylacetamide, dimethoxyethane, dimethylformamide, gamma-butyrolactone, pyridine, and sulfolane have been calculated using the cluster-continuum quasichemical theory. With n independent solvent molecules S initial state forming the "monomer" thermodynamic cycle, Li+ solvation free energies are found to be on average 14 kcal mol-1 more positive compared to those from the "cluster" thermodynamic cycle where the initial state is the cluster Sn. We ascribe the inconsistency between the "monomer" and "cluster" cycles mainly to the incorrectly predicted solvation free energies of solvent clusters Sn from the SMD and CPCM continuum solvation models, which is in line with the earlier study of Bryantsev et al., J. Phys. Chem. B, 2008, 112, 9709-9719. When experimental-based solvation free energies of individual solvent molecules and solvent clusters are employed, the "monomer" and "cluster" cycles result in identical numbers. The best overall agreement with experimental-based "bulk" scale lithium cation solvation free energies was obtained for the "monomer" scale, and we recommend this set of values. We expect that further progress in the field is possible if (i) consensus on the accuracy of experimental reference values is achieved; (ii) the most recent continuum solvation models are properly parameterized for all solute-solvent combinations and become widely accessible for testing.
Collapse
Affiliation(s)
- Daniil Itkis
- N.N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russia.
| | | | | | | |
Collapse
|
10
|
Busato M, Melchior A, Migliorati V, Colella A, Persson I, Mancini G, Veclani D, D'Angelo P. Elusive Coordination of the Ag + Ion in Aqueous Solution: Evidence for a Linear Structure. Inorg Chem 2020; 59:17291-17302. [PMID: 33233885 DOI: 10.1021/acs.inorgchem.0c02494] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
X-ray absorption spectroscopy (XAS) has been employed to study the coordination of the Ag+ ion in aqueous solution. The conjunction of extended X-ray absorption fine structure (EXAFS) and X-ray absorption near-edge structure (XANES) data analysis provided results suggesting the preference for a first shell linear coordination with a mean Ag-O bond distance of 2.34(2) Å, different from the first generally accepted tetrahedral model with a longer mean Ag-O bond distance. Ab initio molecular dynamics simulations with the Car-Parrinello approach (CPMD) were also performed and were able to describe the coordination of the hydrated Ag+ ion in aqueous solution in very good agreement with the experimental data. The high sensitivity for the closest environment of the photoabsorber of the EXAFS and XANES techniques, together with the long-range information provided by CPMD and large-angle X-ray scattering (LAXS), allowed us to reconstruct the three-dimensional model of the coordination geometry around the Ag+ ion in aqueous solution. The obtained results from experiments and theoretical simulations provided a complex picture with a certain amount of water molecules with high configurational disorder at distances comprised between the first and second hydration spheres. This evidence may have caused the proliferation of the coordination numbers that have been proposed so far for Ag+ in water. Altogether these data show how the description of the hydration of the Ag+ ion in aqueous solution can be complex, differently from other metal species where hydration structures can be described by clusters with well-defined geometries. This diffuse hydration shell causes the Ag-O bond distance in the linear [Ag(H2O)2]+ ion to be ca. 0.2 Å longer than in isolated ions in solid state.
Collapse
Affiliation(s)
- Matteo Busato
- Dipartimento Politecnico di Ingegneria e Architettura (DPIA), Laboratori di Chimica, Università di Udine, via delle Scienze 99, 33100 Udine, Italy.,Dipartimento di Chimica, Università di Roma "La Sapienza", P.le A. Moro 5, 00185 Roma, Italy
| | - Andrea Melchior
- Dipartimento Politecnico di Ingegneria e Architettura (DPIA), Laboratori di Chimica, Università di Udine, via delle Scienze 99, 33100 Udine, Italy
| | - Valentina Migliorati
- Dipartimento di Chimica, Università di Roma "La Sapienza", P.le A. Moro 5, 00185 Roma, Italy
| | - Andrea Colella
- Dipartimento di Chimica, Università di Roma "La Sapienza", P.le A. Moro 5, 00185 Roma, Italy
| | - Ingmar Persson
- Department of Molecular Sciences, Swedish University of Agricultural Sciences (SLU), SE-750 07 Uppsala, Sweden
| | - Giordano Mancini
- Centro HPC, Scuola Normale Superiore, Piazza San Silvestro 12, 56125 Pisa, Italy
| | - Daniele Veclani
- Dipartimento Politecnico di Ingegneria e Architettura (DPIA), Laboratori di Chimica, Università di Udine, via delle Scienze 99, 33100 Udine, Italy
| | - Paola D'Angelo
- Dipartimento di Chimica, Università di Roma "La Sapienza", P.le A. Moro 5, 00185 Roma, Italy
| |
Collapse
|
11
|
Absolute ion hydration free energy scale and the surface potential of water via quantum simulation. Proc Natl Acad Sci U S A 2020; 117:30151-30158. [PMID: 33203676 DOI: 10.1073/pnas.2017214117] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
With a goal of determining an absolute free energy scale for ion hydration, quasi-chemical theory and ab initio quantum mechanical simulations are employed to obtain an accurate value for the bulk hydration free energy of the Na+ ion. The free energy is partitioned into three parts: 1) the inner-shell or chemical contribution that includes direct interactions of the ion with nearby waters, 2) the packing free energy that is the work to produce a cavity of size λ in water, and 3) the long-range contribution that involves all interactions outside the inner shell. The interfacial potential contribution to the free energy resides in the long-range term. By averaging cation and anion data for that contribution, cumulant terms of all odd orders in the electrostatic potential are removed. The computed total is then the bulk hydration free energy. Comparison with the experimentally derived real hydration free energy produces an effective surface potential of water in the range -0.4 to -0.5 V. The result is consistent with a variety of experiments concerning acid-base chemistry, ion distributions near hydrophobic interfaces, and electric fields near the surface of water droplets.
Collapse
|
12
|
Maldonado AM, Basdogan Y, Berryman JT, Rempe SB, Keith JA. First-principles modeling of chemistry in mixed solvents: Where to go from here? J Chem Phys 2020; 152:130902. [PMID: 32268733 DOI: 10.1063/1.5143207] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
Mixed solvents (i.e., binary or higher order mixtures of ionic or nonionic liquids) play crucial roles in chemical syntheses, separations, and electrochemical devices because they can be tuned for specific reactions and applications. Apart from fully explicit solvation treatments that can be difficult to parameterize or computationally expensive, there is currently no well-established first-principles regimen for reliably modeling atomic-scale chemistry in mixed solvent environments. We offer our perspective on how this process could be achieved in the near future as mixed solvent systems become more explored using theoretical and computational chemistry. We first outline what makes mixed solvent systems far more complex compared to single-component solvents. An overview of current and promising techniques for modeling mixed solvent environments is provided. We focus on so-called hybrid solvation treatments such as the conductor-like screening model for real solvents and the reference interaction site model, which are far less computationally demanding than explicit simulations. We also propose that cluster-continuum approaches rooted in physically rigorous quasi-chemical theory provide a robust, yet practical, route for studying chemical processes in mixed solvents.
Collapse
Affiliation(s)
- Alex M Maldonado
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Yasemin Basdogan
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Joshua T Berryman
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Susan B Rempe
- Center for Computational Biology and Biophysics, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - John A Keith
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| |
Collapse
|
13
|
Chaudhari MI, Vanegas JM, Pratt LR, Muralidharan A, Rempe SB. Hydration Mimicry by Membrane Ion Channels. Annu Rev Phys Chem 2020; 71:461-484. [PMID: 32155383 DOI: 10.1146/annurev-physchem-012320-015457] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Ions transiting biomembranes might pass readily from water through ion-specific membrane proteins if these protein channels provide environments similar to the aqueous solution hydration environment. Indeed, bulk aqueous solution is an important reference condition for the ion permeation process. Assessment of this hydration mimicry concept depends on understanding the hydration structure and free energies of metal ions in water in order to provide a comparison for the membrane channel environment. To refine these considerations, we review local hydration structures of ions in bulk water and the molecular quasi-chemical theory that provides hydration free energies. In doing so, we note some current views of ion binding to membrane channels and suggest new physical chemical calculations and experiments that might further clarify the hydration mimicry concept.
Collapse
Affiliation(s)
- Mangesh I Chaudhari
- Department of Computational Biology and Biophysics, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA;
| | - Juan M Vanegas
- Department of Computational Biology and Biophysics, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA; .,Current affiliation: Department of Physics, University of Vermont, Burlington, Vermont 05405, USA
| | - L R Pratt
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, USA
| | - Ajay Muralidharan
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, USA.,Current affiliation: Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| | - Susan B Rempe
- Department of Computational Biology and Biophysics, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA;
| |
Collapse
|
14
|
Farshad M, Rasaiah JC. Reverse Translocation of Nucleotides through a Carbon Nanotube. J Phys Chem B 2020; 124:937-943. [PMID: 31917564 DOI: 10.1021/acs.jpcb.9b09587] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report molecular dynamics (MD) simulations of the reverse translocation of single nucleotides through a narrow carbon nanotube (CNT), with a diameter of 1.36 nm, immersed in a 1 M KCl electrolyte solution under an applied electric field along the tube axis. We observe ion selectivity by the narrow CNT, which leads to a high flow of K+ ions, in contrast to a negligible and opposing current of Cl- ions. The K+ ions, driven by the electric field, force a negatively charged single nucleotide into the narrow CNT where it is trapped by the incoming K+ ions and water molecules, and the nucleotide is driven in the same direction as the K+ ions. This illustrates a novel mechanism of nucleotide reverse translocation that is controlled by ion selectivity. An increase in the CNT diameter to 2.71 nm or an increase in nucleotide chain length both lead to translocation in the normal direction of the applied field. The reverse translocation rate of single nucleotides is correlated to the ionic current of K+ ions in the narrow tube, unlike translocation in the normal direction in the wider tube.
Collapse
Affiliation(s)
- Mohsen Farshad
- Department of Chemistry , University of Maine , Orono , Maine 04469 , United States
| | - Jayendran C Rasaiah
- Department of Chemistry , University of Maine , Orono , Maine 04469 , United States
| |
Collapse
|
15
|
Xu Z, Chen K, Li M, Hu C, Yin P. Sustained release of Ag+ confined inside polyoxometalates for long-lasting bacterial resistance. Chem Commun (Camb) 2020; 56:5287-5290. [DOI: 10.1039/d0cc01676d] [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
The release of Ag+ confined in the cavities of nanoscale inorganic clusters can be selectively triggered by the Na+ present in solutions or biological media for long-lasting bacteriostasis.
Collapse
Affiliation(s)
- Zhewei Xu
- Packaging Engineering Institute
- Jinan University
- Zhuhai 519070
- P. R. China
- South China Advanced Institute for Soft Matter Science and Technology
| | - Kun Chen
- South China Advanced Institute for Soft Matter Science and Technology
- South China University of Technology
- Guangzhou 510641
- P. R. China
- State Key Laboratory of Luminescent Materials and Devices
| | - Mu Li
- South China Advanced Institute for Soft Matter Science and Technology
- South China University of Technology
- Guangzhou 510641
- P. R. China
| | - Changying Hu
- Packaging Engineering Institute
- Jinan University
- Zhuhai 519070
- P. R. China
- Department of Food Science and Engineering
| | - Panchao Yin
- South China Advanced Institute for Soft Matter Science and Technology
- South China University of Technology
- Guangzhou 510641
- P. R. China
- State Key Laboratory of Luminescent Materials and Devices
| |
Collapse
|
16
|
Basdogan Y, Groenenboom MC, Henderson E, De S, Rempe SB, Keith JA. Machine Learning-Guided Approach for Studying Solvation Environments. J Chem Theory Comput 2019; 16:633-642. [PMID: 31809056 DOI: 10.1021/acs.jctc.9b00605] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Molecular-level understanding and characterization of solvation environments are often needed across chemistry, biology, and engineering. Toward practical modeling of local solvation effects of any solute in any solvent, we report a static and all-quantum mechanics-based cluster-continuum approach for calculating single-ion solvation free energies. This approach uses a global optimization procedure to identify low-energy molecular clusters with different numbers of explicit solvent molecules and then employs the smooth overlap for atomic positions learning kernel to quantify the similarity between different low-energy solute environments. From these data, we use sketch maps, a nonlinear dimensionality reduction algorithm, to obtain a two-dimensional visual representation of the similarity between solute environments in differently sized microsolvated clusters. After testing this approach on different ions having charges 2+, 1+, 1-, and 2-, we find that the solvation environment around each ion can be seen to usually become more similar in hand with its calculated single-ion solvation free energy. Without needing either dynamics simulations or an a priori knowledge of local solvation structure of the ions, this approach can be used to calculate solvation free energies within 5% of experimental measurements for most cases, and it should be transferable for the study of other systems where dynamics simulations are not easily carried out.
Collapse
Affiliation(s)
- Yasemin Basdogan
- Department of Chemical and Petroleum Engineering Swanson School of Engineering , University of Pittsburgh , Pittsburgh 15261 , Pennsylvania , United States
| | - Mitchell C Groenenboom
- Department of Chemical and Petroleum Engineering Swanson School of Engineering , University of Pittsburgh , Pittsburgh 15261 , Pennsylvania , United States
| | - Ethan Henderson
- Department of Chemical and Petroleum Engineering Swanson School of Engineering , University of Pittsburgh , Pittsburgh 15261 , Pennsylvania , United States
| | - Sandip De
- Laboratory of Computational Science and Modelling, Institute of Materials , École Polytechnique Fédérale de Lausanne , Lausanne 1015 , Switzerland
| | - Susan B Rempe
- Department of Nanobiology , Sandia National Laboratories , Albuquerque 87185 , New Mexico , United States
| | - John A Keith
- Department of Chemical and Petroleum Engineering Swanson School of Engineering , University of Pittsburgh , Pittsburgh 15261 , Pennsylvania , United States
| |
Collapse
|
17
|
Yang XH, Cuesta A, Cheng J. Computational Ag/AgCl Reference Electrode from Density Functional Theory-Based Molecular Dynamics. J Phys Chem B 2019; 123:10224-10232. [PMID: 31693366 DOI: 10.1021/acs.jpcb.9b06650] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have developed a scheme to compute the standard potential of the Ag/AgCl reference electrode using density functional theory-based molecular dynamics, similar to the computational standard hydrogen electrode (SHE) developed by Cheng, Sulpizi, and Sprik [J. Chem. Phys. 2009, 131, 154504], with which our new computational reference electrode was compared. We have obtained a similar value of the potential of the Ag/AgCl electrode versus SHE to the experiment. The newly developed computational reference electrode will be extended to nonaqueous solvents in the future, where it will be used to predict standard equilibrium potentials to be compared with experimental data.
Collapse
Affiliation(s)
- Xiao-Hui Yang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering , Xiamen University , Xiamen 361005 , China.,Department of Chemistry , University of Aberdeen , Aberdeen AB24 3UE , U.K
| | - Angel Cuesta
- Department of Chemistry , University of Aberdeen , Aberdeen AB24 3UE , U.K
| | - Jun Cheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering , Xiamen University , Xiamen 361005 , China
| |
Collapse
|
18
|
von Rudorff GF, von Lilienfeld OA. Atoms in Molecules from Alchemical Perturbation Density Functional Theory. J Phys Chem B 2019; 123:10073-10082. [PMID: 31647233 DOI: 10.1021/acs.jpcb.9b07799] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Based on thermodynamic integration, we introduce atoms in molecules (AIM) using the orbital-free framework of alchemical perturbation density functional theory (APDFT). Within APDFT, atomic energies and electron densities in molecules are arbitrary because any reference system and integration path can be selected as long as it meets the boundary conditions. We choose the uniform electron gas (jellium) as a reference and linearly scale up all nuclear charges, situated at any query molecule's atomic coordinates. Within the approximations made when calculating one-particle electron densities, this universal choice affords unambiguous and exact definitions of energies and electron densities of AIMs. Numerical results are presented for neutral small molecules (CO, N2, BF, CO2), various small molecules with different electronic hybridization states of carbon (CH4, C2H6, C2H4, C2H2, HCN), and all of the possible BN-doped mutants connecting benzene to borazine (C2nB3-nN3-nH6, 0 ≤ n ≤ 3). Our results, as well as comparison to atomic energy estimates resulting from either DFT trained neural network models or atomic basis set overlap within CCSD, suggest that APDFT based AIMs enable meaningful, interesting, and counterintuitive interpretations of chemical bonding and molecular electron densities.
Collapse
Affiliation(s)
- Guido Falk von Rudorff
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry , University of Basel , Klingelbergstrasse 80 , CH-4056 Basel , Switzerland
| | - O Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry , University of Basel , Klingelbergstrasse 80 , CH-4056 Basel , Switzerland
| |
Collapse
|
19
|
Muralidharan A, Pratt L, Chaudhari M, Rempe S. Quasi-chemical theory for anion hydration and specific ion effects: Cl-(aq) vs. F-(aq). ACTA ACUST UNITED AC 2019. [DOI: 10.1016/j.cpletx.2019.100037] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
20
|
Baskin A, Prendergast D. "Ion Solvation Spectra": Free Energy Analysis of Solvation Structures of Multivalent Cations in Aprotic Solvents. J Phys Chem Lett 2019; 10:4920-4928. [PMID: 31322350 DOI: 10.1021/acs.jpclett.9b01569] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Using advanced molecular dynamics free energy sampling techniques-both classical and ab initio-we analyze the solvation structures of multivalent cations in aprotic solvents. In contrast to previous studies of mono- and bivalent ions in organic solvents, mainly performed using hybrid cluster-continuum quantum chemistry calculations that rely on the assumption of uniqueness of ion solvation free energies, here we find that monatomic bivalent cations may have multiple well-defined minima, as previously reported only for water, or plateaus of free energy with respect to the ion-solvent coordination. These observations are generalized in the concept of the "ion solvation spectrum" to highlight the rich phenomenology related to ion solvation as opposed to the normally expected free energy profiles with a single coordination minimum. Specifically, we show that a single chemical species may exhibit a multiplicity of distinctly different electrochemical properties. Using one- and two-dimensional projections of the free energy landscape, we analyze the stability of ion solvation structures and reveal minimum free energy pathways for ion (de-)solvation with low-dimensional approximations to associated kinetic barriers. Unexpectedly, we show that in some cases the process of opening the first ion solvation shell, by removing a solvent molecule, may actually drive the ion into a free energy basin with a higher coordination number. Our study highlights some deficiencies of conventional methodologies for studying ion solvation as a path to determine redox potentials and provides experimentally testable predictions.
Collapse
Affiliation(s)
- Artem Baskin
- Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - David Prendergast
- Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
21
|
Shevkunov SV. Water Vapor Nucleation on a Surface with Nanoscopic Grooves. 2. Features of Thermodynamic Behavior. COLLOID JOURNAL 2019. [DOI: 10.1134/s1061933x19030141] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
22
|
Doyle CC, Shi Y, Beck TL. The Importance of the Water Molecular Quadrupole for Estimating Interfacial Potential Shifts Acting on Ions Near the Liquid–Vapor Interface. J Phys Chem B 2019; 123:3348-3358. [DOI: 10.1021/acs.jpcb.9b01289] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Carrie C. Doyle
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Yu Shi
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Thomas L. Beck
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| |
Collapse
|
23
|
Muralidharan A, Pratt LR, Chaudhari MI, Rempe SB. Quasi-Chemical Theory with Cluster Sampling from Ab Initio Molecular Dynamics: Fluoride (F -) Anion Hydration. J Phys Chem A 2018; 122:9806-9812. [PMID: 30475612 DOI: 10.1021/acs.jpca.8b08474] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Accurate predictions of the hydration free energy for anions typically has been more challenging than that for cations. Hydrogen bond donation to the anion in hydrated clusters such as F(H2O) n - can lead to delicate structures. Consequently, the energy landscape contains many local minima, even for small clusters, and these minima present a challenge for computational optimization. Utilization of cluster experimental results for the free energies of gas-phase clusters shows that even though anharmonic effects are interesting they need not be of troublesome magnitudes for careful applications of quasi-chemical theory to ion hydration. Energy-optimized cluster structures for anions can leave the central ion highly exposed, and application of implicit solvation models to these structures can incur more serious errors than those for metal cations. Utilizing cluster structures sampled from ab initio molecular dynamics simulations substantially fixes those issues.
Collapse
Affiliation(s)
- A Muralidharan
- Department of Chemical and Biomolecular Engineering , Tulane University , New Orleans , Louisiana 70118 , United States
| | - L R Pratt
- Department of Chemical and Biomolecular Engineering , Tulane University , New Orleans , Louisiana 70118 , United States
| | - M I Chaudhari
- Center for Biological and Engineering Sciences , Sandia National Laboratories , Albuquerque , New Mexico 87185 , United States
| | - S B Rempe
- Center for Biological and Engineering Sciences , Sandia National Laboratories , Albuquerque , New Mexico 87185 , United States
| |
Collapse
|
24
|
Pethes I. The structure of aqueous lithium chloride solutions at high concentrations as revealed by a comparison of classical interatomic potential models. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.05.044] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
|
25
|
Hofer TS, Hünenberger PH. Absolute proton hydration free energy, surface potential of water, and redox potential of the hydrogen electrode from first principles: QM/MM MD free-energy simulations of sodium and potassium hydration. J Chem Phys 2018; 148:222814. [DOI: 10.1063/1.5000799] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Thomas S. Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, Centre for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | | |
Collapse
|
26
|
Chaudhari MI, Rempe SB, Pratt LR. Quasi-chemical theory of F -(aq): The "no split occupancies rule" revisited. J Chem Phys 2018; 147:161728. [PMID: 29096480 DOI: 10.1063/1.4986244] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
We use ab initio molecular dynamics (AIMD) calculations and quasi-chemical theory (QCT) to study the inner-shell structure of F-(aq) and to evaluate that single-ion free energy under standard conditions. Following the "no split occupancies" rule, QCT calculations yield a free energy value of -101 kcal/mol under these conditions, in encouraging agreement with tabulated values (-111 kcal/mol). The AIMD calculations served only to guide the definition of an effective inner-shell constraint. QCT naturally includes quantum mechanical effects that can be concerning in more primitive calculations, including electronic polarizability and induction, electron density transfer, electron correlation, molecular/atomic cooperative interactions generally, molecular flexibility, and zero-point motion. No direct assessment of the contribution of dispersion contributions to the internal energies has been attempted here, however. We anticipate that other aqueous halide ions might be treated successfully with QCT, provided that the structure of the underlying statistical mechanical theory is absorbed, i.e., that the "no split occupancies" rule is recognized.
Collapse
Affiliation(s)
- Mangesh I Chaudhari
- Center for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Susan B Rempe
- Center for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Lawrence R Pratt
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, USA
| |
Collapse
|
27
|
Ahlstrand E, Zukerman Schpector J, Friedman R. Computer simulations of alkali-acetate solutions: Accuracy of the forcefields in difference concentrations. J Chem Phys 2018; 147:194102. [PMID: 29166095 DOI: 10.1063/1.4985919] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
When proteins are solvated in electrolyte solutions that contain alkali ions, the ions interact mostly with carboxylates on the protein surface. Correctly accounting for alkali-carboxylate interactions is thus important for realistic simulations of proteins. Acetates are the simplest carboxylates that are amphipathic, and experimental data for alkali acetate solutions are available and can be compared with observables obtained from simulations. We carried out molecular dynamics simulations of alkali acetate solutions using polarizable and non-polarizable forcefields and examined the ion-acetate interactions. In particular, activity coefficients and association constants were studied in a range of concentrations (0.03, 0.1, and 1M). In addition, quantum-mechanics (QM) based energy decomposition analysis was performed in order to estimate the contribution of polarization, electrostatics, dispersion, and QM (non-classical) effects on the cation-acetate and cation-water interactions. Simulations of Li-acetate solutions in general overestimated the binding of Li+ and acetates. In lower concentrations, the activity coefficients of alkali-acetate solutions were too high, which is suggested to be due to the simulation protocol and not the forcefields. Energy decomposition analysis suggested that improvement of the forcefield parameters to enable accurate simulations of Li-acetate solutions can be achieved but may require the use of a polarizable forcefield. Importantly, simulations with some ion parameters could not reproduce the correct ion-oxygen distances, which calls for caution in the choice of ion parameters when protein simulations are performed in electrolyte solutions.
Collapse
Affiliation(s)
- Emma Ahlstrand
- Department of Chemistry and Biomedical Sciences, Linnæus University, 391 82 Kalmar, Sweden
| | - Julio Zukerman Schpector
- Universidade Federal de São Carlos, Departamento de Química, CP 676, 13565-905 São Carlos, SP, Brazil
| | - Ran Friedman
- Department of Chemistry and Biomedical Sciences, Linnæus University, 391 82 Kalmar, Sweden
| |
Collapse
|
28
|
Chaudhari MI, Muralidharan A, Pratt LR, Rempe SB. Assessment of Simple Models for Molecular Simulation of Ethylene Carbonate and Propylene Carbonate as Solvents for Electrolyte Solutions. Top Curr Chem (Cham) 2018; 376:7. [PMID: 29435669 PMCID: PMC5809610 DOI: 10.1007/s41061-018-0187-2] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2017] [Accepted: 01/23/2018] [Indexed: 01/13/2023]
Abstract
Progress in understanding liquid ethylene carbonate (EC) and propylene carbonate (PC) on the basis of molecular simulation, emphasizing simple models of interatomic forces, is reviewed. Results on the bulk liquids are examined from the perspective of anticipated applications to materials for electrical energy storage devices. Preliminary results on electrochemical double-layer capacitors based on carbon nanotube forests and on model solid-electrolyte interphase (SEI) layers of lithium ion batteries are considered as examples. The basic results discussed suggest that an empirically parameterized, non-polarizable force field can reproduce experimental structural, thermodynamic, and dielectric properties of EC and PC liquids with acceptable accuracy. More sophisticated force fields might include molecular polarizability and Buckingham-model description of inter-atomic overlap repulsions as extensions to Lennard-Jones models of van der Waals interactions. Simple approaches should be similarly successful also for applications to organic molecular ions in EC/PC solutions, but the important case of Li\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^+$$\end{document}+ deserves special attention because of the particularly strong interactions of that small ion with neighboring solvent molecules. To treat the Li\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^+$$\end{document}+ ions in liquid EC/PC solutions, we identify interaction models defined by empirically scaled partial charges for ion-solvent interactions. The empirical adjustments use more basic inputs, electronic structure calculations and ab initio molecular dynamics simulations, and also experimental results on Li\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^+$$\end{document}+ thermodynamics and transport in EC/PC solutions. Application of such models to the mechanism of Li\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^+$$\end{document}+ transport in glassy SEI models emphasizes the advantage of long time-scale molecular dynamics studies of these non-equilibrium materials.
Collapse
Affiliation(s)
- Mangesh I Chaudhari
- Center for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, NM 87185, USA
| | - Ajay Muralidharan
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, LA 70118, USA
| | - Lawrence R Pratt
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, LA 70118, USA
| | - Susan B Rempe
- Center for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, NM 87185, USA.
| |
Collapse
|
29
|
Duignan TT, Baer MD, Schenter GK, Mundy CJ. Electrostatic solvation free energies of charged hard spheres using molecular dynamics with density functional theory interactions. J Chem Phys 2017; 147:161716. [DOI: 10.1063/1.4994912] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- Timothy T. Duignan
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Marcel D. Baer
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Gregory K. Schenter
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Chistopher J. Mundy
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98185, USA
| |
Collapse
|
30
|
Duignan TT, Baer MD, Schenter GK, Mundy CJ. Real single ion solvation free energies with quantum mechanical simulation. Chem Sci 2017; 8:6131-6140. [PMID: 28989643 PMCID: PMC5625628 DOI: 10.1039/c7sc02138k] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2017] [Accepted: 05/26/2017] [Indexed: 01/11/2023] Open
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 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.
Collapse
Affiliation(s)
- Timothy T Duignan
- Physical Science Division , Pacific Northwest National Laboratory , P.O. Box 999 , Richland , Washington 99352 , USA . ; Tel: +1 509 3756940
| | - Marcel D Baer
- Physical Science Division , Pacific Northwest National Laboratory , P.O. Box 999 , Richland , Washington 99352 , USA . ; Tel: +1 509 3756940
| | - Gregory K Schenter
- Physical Science Division , Pacific Northwest National Laboratory , P.O. Box 999 , Richland , Washington 99352 , USA . ; Tel: +1 509 3756940
| | - Christopher J Mundy
- Physical Science Division , Pacific Northwest National Laboratory , P.O. Box 999 , Richland , Washington 99352 , USA . ; Tel: +1 509 3756940
| |
Collapse
|
31
|
Frenzel J, Meyer B, Marx D. Bicanonical ab Initio Molecular Dynamics for Open Systems. J Chem Theory Comput 2017. [DOI: 10.1021/acs.jctc.7b00263] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Johannes Frenzel
- Lehrstuhl für
Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Bernd Meyer
- Interdisziplinäres
Zentrum für Molekulare Materialien (ICMM) and Computer-Chemie-Centrum
(CCC), Friedrich-Alexander-Universität Erlangen-Nürnberg, 91052 Erlangen, Germany
| | - Dominik Marx
- Lehrstuhl für
Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
32
|
Chaudhari MI, Pratt LR, Rempe SB. Utility of chemical computations in predicting solution free energies of metal ions. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1342127] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Mangesh I. Chaudhari
- Center for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, NM, USA
| | - Lawrence R. Pratt
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, LA, USA
| | - Susan B. Rempe
- Center for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, NM, USA
| |
Collapse
|
33
|
Shi X, Xiao H, Chen X, Lackner KS. The Effect of Moisture on the Hydrolysis of Basic Salts. Chemistry 2016; 22:18326-18330. [DOI: 10.1002/chem.201603701] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2016] [Indexed: 11/11/2022]
Affiliation(s)
- Xiaoyang Shi
- Department of Earth and Environmental Engineering; Columbia University, 905E Mudd; 500 West 120th Street New York NY 10027 USA
| | - Hang Xiao
- Department of Earth and Environmental Engineering; Columbia University, 905E Mudd; 500 West 120th Street New York NY 10027 USA
| | - Xi Chen
- Department of Earth and Environmental Engineering; Columbia University, 905E Mudd; 500 West 120th Street New York NY 10027 USA
| | - Klaus S. Lackner
- School of Sustainable Engineering & Built Environment; Arizona State University; Tempe AZ 85287-9309 USA
| |
Collapse
|
34
|
Chaudhari MI, Nair JR, Pratt LR, Soto FA, Balbuena PB, Rempe SB. Scaling Atomic Partial Charges of Carbonate Solvents for Lithium Ion Solvation and Diffusion. J Chem Theory Comput 2016; 12:5709-5718. [DOI: 10.1021/acs.jctc.6b00824] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Mangesh I. Chaudhari
- Center
for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Jijeesh R. Nair
- Department
of Applied Science and Technology, Politecnico di Torino, Turin 10129, Italy
| | - Lawrence R. Pratt
- Department
of Chemical and Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States
| | - Fernando A. Soto
- Department of Chemical Engineering, Texas A&M University, College Station, Texas 77843, United States
| | - Perla B. Balbuena
- Department of Chemical Engineering, Texas A&M University, College Station, Texas 77843, United States
| | - Susan B. Rempe
- Center
for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| |
Collapse
|
35
|
Pham TA, Ogitsu T, Lau EY, Schwegler E. Structure and dynamics of aqueous solutions from PBE-based first-principles molecular dynamics simulations. J Chem Phys 2016; 145:154501. [DOI: 10.1063/1.4964865] [Citation(s) in RCA: 82] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Tuan Anh Pham
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Tadashi Ogitsu
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Edmond Y. Lau
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Eric Schwegler
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| |
Collapse
|
36
|
Soniat M, Rogers DM, Rempe SB. Dispersion- and Exchange-Corrected Density Functional Theory for Sodium Ion Hydration. J Chem Theory Comput 2016; 11:2958-67. [PMID: 26575733 DOI: 10.1021/acs.jctc.5b00357] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
A challenge in density functional theory is developing functionals that simultaneously describe intermolecular electron correlation and electron delocalization. Recent exchange-correlation functionals address those two issues by adding corrections important at long ranges: an atom-centered pairwise dispersion term to account for correlation and a modified long-range component of the electron exchange term to correct for delocalization. Here we investigate how those corrections influence the accuracy of binding free energy predictions for sodium-water clusters. We find that the dual-corrected ωB97X-D functional gives cluster binding energies closest to high-level ab initio methods (CCSD(T)). Binding energy decomposition shows that the ωB97X-D functional predicts the smallest ion-water (pairwise) interaction energy and larger multibody contributions for a four-water cluster than most other functionals - a trend consistent with CCSD(T) results. Also, ωB97X-D produces the smallest amounts of charge transfer and the least polarizable waters of the density functionals studied, which mimics the lower polarizability of CCSD. When compared with experimental binding free energies, however, the exchange-corrected CAM-B3LYP functional performs best (error <1 kcal/mol), possibly because of its parametrization to experimental formation enthalpies. For clusters containing more than four waters, "split-shell" coordination must be considered to obtain accurate free energies in comparison with experiment.
Collapse
Affiliation(s)
- Marielle Soniat
- Department of Chemistry, University of New Orleans , 2000 Lakeshore Drive, New Orleans, Louisiana 70148, United States
| | - David M Rogers
- Center for Biological and Engineering Sciences, Sandia National Laboratories , Albuquerque, New Mexico 87123, United States.,Department of Chemistry, University of South Florida , Tampa, Florida 33620, United States
| | - Susan B Rempe
- Center for Biological and Engineering Sciences, Sandia National Laboratories , Albuquerque, New Mexico 87123, United States
| |
Collapse
|
37
|
Ions interacting in solution: Moving from intrinsic to collective properties. Curr Opin Colloid Interface Sci 2016. [DOI: 10.1016/j.cocis.2016.05.009] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
38
|
Remsing RC, Weeks JD. Role of Local Response in Ion Solvation: Born Theory and Beyond. J Phys Chem B 2016; 120:6238-49. [DOI: 10.1021/acs.jpcb.6b02238] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Richard C. Remsing
- Institute
for Computational Molecular Science and Center for the Computational
Design of Functional Layered Materials, Temple University, Philadelphia, Pennsylvania 19122, United States
- Institute
for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | - John D. Weeks
- Institute
for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
- Department
of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
39
|
A theoretical study on intermolecular [2+2] radical cation cycloaddition reactions and the competition between concerted and stepwise mechanisms. COMPUT THEOR CHEM 2016. [DOI: 10.1016/j.comptc.2016.01.004] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
40
|
Solovyeva A, von Lilienfeld OA. Alchemical screening of ionic crystals. Phys Chem Chem Phys 2016; 18:31078-31091. [DOI: 10.1039/c6cp04258a] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We introduce alchemical perturbations as a rapid and accurate tool to estimate fundamental structural and energetic properties in pure and mixed ionic crystals.
Collapse
Affiliation(s)
- Alisa Solovyeva
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials MARVEL
- Department of Chemistry
- University of Basel
- 4056 Basel
- Switzerland
| | - O. Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials MARVEL
- Department of Chemistry
- University of Basel
- 4056 Basel
- Switzerland
| |
Collapse
|
41
|
Bankura A, Santra B, DiStasio RA, Swartz CW, Klein ML, Wu X. A systematic study of chloride ion solvation in water using van der Waals inclusive hybrid density functional theory. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1059959] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Arindam Bankura
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Biswajit Santra
- Department of Chemistry, Princeton University , Princeton, NJ, USA
| | | | - Charles W. Swartz
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Michael L. Klein
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Xifan Wu
- Department of Physics, Temple University , Philadelphia, PA, USA
| |
Collapse
|
42
|
Abstract
The hydration structure of Ba(2+) ion is important for understanding blocking mechanisms in potassium ion channels. Here, we combine statistical mechanical theory, ab initio molecular dynamics simulations, and electronic structure methods to calculate the hydration free energy and local hydration structure of Ba(2+)(aq). The predicted hydration free energy (-304 ± 1 kcal/mol) agrees with the experimental value (-303 kcal/mol) when a maximally occupied, unimodal inner solvation shell is treated. In the local environment defined by the first shell of hydrating waters, Ba(2+) is directly and stably coordinated by eight (8) waters. Octa-coordination resembles the crystal structure of Ba(2+) and K(+) bound in potassium ion channels, but differs from the local hydration structure of K(+)(aq) determined earlier.
Collapse
Affiliation(s)
- Mangesh I Chaudhari
- †Center for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Marielle Soniat
- ‡Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| | - Susan B Rempe
- †Center for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| |
Collapse
|
43
|
Mason PE, Ansell S, Neilson GW, Rempe SB. Neutron Scattering Studies of the Hydration Structure of Li+. J Phys Chem B 2015; 119:2003-9. [DOI: 10.1021/jp511508n] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- P. E. Mason
- Institute
of Organic Chemistry
and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo
nám. 2, 16610 Prague 6, Czech Republic
| | - S. Ansell
- Rutherford Appleton Laboratories, Chilton, Oxfordshire OX11 0QX, United Kingdom
| | - G. W. Neilson
- Department Physics, University of Bristol, Tyndall Ave., Bristol BS8 1TL, United Kingdom
| | - S. B. Rempe
- Center for Biological and
Materials Sciences, Sandia National Laboratories, Albuquerque, New Mexico 87185-1315, United States
| |
Collapse
|
44
|
Carvalho NF, Pliego JR. Cluster-continuum quasichemical theory calculation of the lithium ion solvation in water, acetonitrile and dimethyl sulfoxide: an absolute single-ion solvation free energy scale. Phys Chem Chem Phys 2015; 17:26745-55. [DOI: 10.1039/c5cp03798k] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A solvation free energy scale excluding the net electrostatic potential inside the solute cavity is presented.
Collapse
Affiliation(s)
- Nathalia F. Carvalho
- Departamento de Ciências Naturais
- Universidade Federal de São João del Rei
- São João del-Rei
- Brazil
| | - Josefredo R. Pliego
- Departamento de Ciências Naturais
- Universidade Federal de São João del Rei
- São João del-Rei
- Brazil
| |
Collapse
|
45
|
Remsing RC, Baer MD, Schenter GK, Mundy CJ, Weeks JD. The Role of Broken Symmetry in Solvation of a Spherical Cavity in Classical and Quantum Water Models. J Phys Chem Lett 2014; 5:2767-2774. [PMID: 26278076 DOI: 10.1021/jz501067w] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Insertion of a hard sphere cavity in liquid water breaks translational symmetry and generates an electrostatic potential difference between the region near the cavity and the bulk. Here, we clarify the physical interpretation of this potential and its calculation. We also show that the electrostatic potential in the center of small, medium, and large cavities depends very sensitively on the form of the assumed molecular interactions for different classical simple point-charge models and quantum mechanical DFT-based interaction potentials, as reflected in their description of donor and acceptor hydrogen bonds near the cavity. These differences can significantly affect the magnitude of the scalar electrostatic potential. We argue that the result of these studies will have direct consequences toward our understanding of the thermodynamics of ion solvation through the cavity charging process.
Collapse
Affiliation(s)
| | - Marcel D Baer
- §Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Gregory K Schenter
- §Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christopher J Mundy
- §Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | | |
Collapse
|
46
|
Annapureddy HVR, Dang LX. Understanding the rates and molecular mechanism of water-exchange around aqueous ions using molecular simulations. J Phys Chem B 2014; 118:8917-27. [PMID: 24911526 DOI: 10.1021/jp502922c] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Solvation processes occurring around aqueous ions are of fundamental importance in physics, chemistry, and biology. Over the past few decades, several experimental and theoretical studies were devoted to understanding ion solvation and the processes involved in it. In this article, we present a summary of our recent efforts that, through computer simulations, focused on providing a comprehensive understanding of solvent-exchange processes around aqueous ions. To accomplish these activities, we have looked at the mechanistic properties associated with the water-exchange process, such as potentials of mean force, time-dependent transmission coefficients, and the corresponding rate constants using transition state theory, the reactive flux method, and Grote-Hynes treatments of the dynamic response of the solvent.
Collapse
Affiliation(s)
- Harsha V R Annapureddy
- Physical Sciences Division, Pacific Northwest National Laboratory , Richland, Washington 93352, United States
| | | |
Collapse
|
47
|
Lin YL, Aleksandrov A, Simonson T, Roux B. An Overview of Electrostatic Free Energy Computations for Solutions and Proteins. J Chem Theory Comput 2014; 10:2690-709. [PMID: 26586504 DOI: 10.1021/ct500195p] [Citation(s) in RCA: 104] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Free energy simulations for electrostatic and charging processes in complex molecular systems encounter specific difficulties owing to the long-range, 1/r Coulomb interaction. To calculate the solvation free energy of a simple ion, it is essential to take into account the polarization of nearby solvent but also the electrostatic potential drop across the liquid-gas boundary, however distant. The latter does not exist in a simulation model based on periodic boundary conditions because there is no physical boundary to the system. An important consequence is that the reference value of the electrostatic potential is not an ion in a vacuum. Also, in an infinite system, the electrostatic potential felt by a perturbing charge is conditionally convergent and dependent on the choice of computational conventions. Furthermore, with Ewald lattice summation and tinfoil conducting boundary conditions, the charges experience a spurious shift in the potential that depends on the details of the simulation system such as the volume fraction occupied by the solvent. All these issues can be handled with established computational protocols, as reviewed here and illustrated for several small ions and three solvated proteins.
Collapse
Affiliation(s)
- Yen-Lin Lin
- Department of Biochemistry and Molecular Biology, University of Chicago , Chicago, Illinois 60637, United States
| | - Alexey Aleksandrov
- Department of Biology, Laboratoire de Biochimie (CNRS UMR7654), Ecole Polytechnique , 91128 Palaiseau, France
| | - Thomas Simonson
- Department of Biology, Laboratoire de Biochimie (CNRS UMR7654), Ecole Polytechnique , 91128 Palaiseau, France
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, University of Chicago , Chicago, Illinois 60637, United States.,Biosciences Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| |
Collapse
|
48
|
Kahlen J, Salimi L, Sulpizi M, Peter C, Donadio D. Interaction of Charged Amino-Acid Side Chains with Ions: An Optimization Strategy for Classical Force Fields. J Phys Chem B 2014; 118:3960-72. [DOI: 10.1021/jp412490c] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jens Kahlen
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Leila Salimi
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
- Johannes Gutenberg University Mainz, Staudinger Weg 7, 55099 Mainz, Germany
| | - Marialore Sulpizi
- Johannes Gutenberg University Mainz, Staudinger Weg 7, 55099 Mainz, Germany
| | - Christine Peter
- University of Konstanz, P.O. Box 718, 78547 Konstanz, Germany
| | - Davide Donadio
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
49
|
Murakami W, Eda K, Yamamoto M, Osakai T. A Non-Bornian Analysis of the Gibbs Energy of Ion Hydration. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2014. [DOI: 10.1246/bcsj.20130296] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Wataru Murakami
- Department of Chemistry, Graduate School of Science, Kobe University
| | - Kazuo Eda
- Department of Chemistry, Graduate School of Science, Kobe University
| | - Masahiro Yamamoto
- Department of Chemistry of Functional Molecules, Faculty of Science and Engineering, Konan University
| | - Toshiyuki Osakai
- Department of Chemistry, Graduate School of Science, Kobe University
| |
Collapse
|
50
|
Marsalek O, Chen PY, Dupuis R, Benoit M, Méheut M, Bačić Z, Tuckerman ME. Efficient Calculation of Free Energy Differences Associated with Isotopic Substitution Using Path-Integral Molecular Dynamics. J Chem Theory Comput 2014; 10:1440-53. [DOI: 10.1021/ct400911m] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Ondrej Marsalek
- Department
of Chemistry, New York University, New York, New York 10003, United States
| | - Pei-Yang Chen
- Department
of Chemistry, New York University, New York, New York 10003, United States
| | - Romain Dupuis
- Géosciences
Environnement Toulouse, OMP−Université Paul Sabatier, 14 avenue
Edouard Belin, 31400 Toulouse, France
| | - Magali Benoit
- CEMES−CNRS−UPR
8011, 29 rue Jeanne Marvig, 31055 Toulouse, France
| | - Merlin Méheut
- Géosciences
Environnement Toulouse, OMP−Université Paul Sabatier, 14 avenue
Edouard Belin, 31400 Toulouse, France
| | - Zlatko Bačić
- Department
of Chemistry, New York University, New York, New York 10003, United States
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Mark E. Tuckerman
- Department
of Chemistry and Courant Institute of Mathematical Sciences, New York University, New York, New York 10003, United States
- Institute for Pure and Applied Mathematics, 460 Portola Plaza, Los Angeles, California 90095, United States
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|