1
|
Becker M, Loche P, Rezaei M, Wolde-Kidan A, Uematsu Y, Netz RR, Bonthuis DJ. Multiscale Modeling of Aqueous Electric Double Layers. Chem Rev 2024; 124:1-26. [PMID: 38118062 PMCID: PMC10785765 DOI: 10.1021/acs.chemrev.3c00307] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2023] [Revised: 11/17/2023] [Accepted: 11/30/2023] [Indexed: 12/22/2023]
Abstract
From the stability of colloidal suspensions to the charging of electrodes, electric double layers play a pivotal role in aqueous systems. The interactions between interfaces, water molecules, ions and other solutes making up the electrical double layer span length scales from Ångströms to micrometers and are notoriously complex. Therefore, explaining experimental observations in terms of the double layer's molecular structure has been a long-standing challenge in physical chemistry, yet recent advances in simulations techniques and computational power have led to tremendous progress. In particular, the past decades have seen the development of a multiscale theoretical framework based on the combination of quantum density functional theory, force-field based simulations and continuum theory. In this Review, we discuss these theoretical developments and make quantitative comparisons to experimental results from, among other techniques, sum-frequency generation, atomic-force microscopy, and electrokinetics. Starting from the vapor/water interface, we treat a range of qualitatively different types of surfaces, varying from soft to solid, from hydrophilic to hydrophobic, and from charged to uncharged.
Collapse
Affiliation(s)
| | - Philip Loche
- Fachbereich
Physik, Freie Universität Berlin, 14195 Berlin, Germany
- Laboratory
of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Majid Rezaei
- Fachbereich
Physik, Freie Universität Berlin, 14195 Berlin, Germany
- Institute
of Theoretical Chemistry, Ulm University, 89081 Ulm, Germany
| | | | - Yuki Uematsu
- Department
of Physics and Information Technology, Kyushu
Institute of Technology, 820-8502 Iizuka, Japan
- PRESTO,
Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Roland R. Netz
- Fachbereich
Physik, Freie Universität Berlin, 14195 Berlin, Germany
| | - Douwe Jan Bonthuis
- Institute
of Theoretical and Computational Physics, Graz University of Technology, 8010 Graz, Austria
| |
Collapse
|
2
|
The collective burst mechanism of angular jumps in liquid water. Nat Commun 2023; 14:1345. [PMID: 36906703 PMCID: PMC10008639 DOI: 10.1038/s41467-023-37069-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Accepted: 02/24/2023] [Indexed: 03/13/2023] Open
Abstract
Understanding the microscopic origins of collective reorientational motions in aqueous systems requires techniques that allow us to reach beyond our chemical imagination. Herein, we elucidate a mechanism using a protocol that automatically detects abrupt motions in reorientational dynamics, showing that large angular jumps in liquid water involve highly cooperative orchestrated motions. Our automatized detection of angular fluctuations, unravels a heterogeneity in the type of angular jumps occurring concertedly in the system. We show that large orientational motions require a highly collective dynamical process involving correlated motion of many water molecules in the hydrogen-bond network that form spatially connected clusters going beyond the local angular jump mechanism. This phenomenon is rooted in the collective fluctuations of the network topology which results in the creation of defects in waves on the THz timescale. The mechanism we propose involves a cascade of hydrogen-bond fluctuations underlying angular jumps and provides new insights into the current localized picture of angular jumps, and its wide use in the interpretations of numerous spectroscopies as well in reorientational dynamics of water near biological and inorganic systems. The role of finite size effects, as well as of the chosen water model, on the collective reorientation is also elucidated.
Collapse
|
3
|
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
|
4
|
Becker MR, Loche P, Netz RR. Electrokinetic, electrochemical, and electrostatic surface potentials of the pristine water liquid-vapor interface. J Chem Phys 2022; 157:240902. [PMID: 36586978 DOI: 10.1063/5.0127869] [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/05/2022] Open
Abstract
Although conceptually simple, the air-water interface displays rich behavior and is subject to intense experimental and theoretical investigations. Different definitions of the electrostatic surface potential as well as different calculation methods, each relevant for distinct experimental scenarios, lead to widely varying potential magnitudes and sometimes even different signs. Based on quantum-chemical density-functional-theory molecular dynamics (DFT-MD) simulations, different surface potentials are evaluated and compared to force-field (FF) MD simulations. As well explained in the literature, the laterally averaged electrostatic surface potential, accessible to electron holography, is dominated by the trace of the water molecular quadrupole moment, and using DFT-MD amounts to +4.35 V inside the water phase, very different from results obtained with FF water models which yield negative values of the order of -0.4 to -0.6 V. Thus, when predicting potentials within water molecules, as relevant for photoelectron spectroscopy and non-linear interface-specific spectroscopy, DFT simulations should be used. The electrochemical surface potential, relevant for ion transfer reactions and ion surface adsorption, is much smaller, less than 200 mV in magnitude, and depends specifically on the ion radius. Charge transfer between interfacial water molecules leads to a sizable surface potential as well. However, when probing electrokinetics by explicitly applying a lateral electric field in DFT-MD simulations, the electrokinetic ζ-potential turns out to be negligible, in agreement with predictions using continuous hydrodynamic models. Thus, interfacial polarization charges from intermolecular charge transfer do not lead to significant electrokinetic mobility at the pristine vapor-liquid water interface, even assuming these transfer charges are mobile in an external electric field.
Collapse
Affiliation(s)
| | - Philip Loche
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany
| | - Roland R Netz
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
5
|
Witek J, Heindel JP, Guan X, Leven I, Hao H, Naullage P, LaCour A, Sami S, Menger MFSJ, Cofer-Shabica DV, Berquist E, Faraji S, Epifanovsky E, Head-Gordon T. M-Chem: a Modular Software Package for Molecular Simulation that Spans Scientific Domains. Mol Phys 2022; 121:e2129500. [PMID: 37470065 PMCID: PMC10353727 DOI: 10.1080/00268976.2022.2129500] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2022] [Accepted: 09/06/2022] [Indexed: 10/10/2022]
Abstract
We present a new software package called M-Chem that is designed from scratch in C++ and parallelized on shared-memory multi-core architectures to facilitate efficient molecular simulations. Currently, M-Chem is a fast molecular dynamics (MD) engine that supports the evaluation of energies and forces from two-body to many-body all-atom potentials, reactive force fields, coarse-grained models, combined quantum mechanics molecular mechanics (QM/MM) models, and external force drivers from machine learning, augmented by algorithms that are focused on gains in computational simulation times. M-Chem also includes a range of standard simulation capabilities including thermostats, barostats, multi-timestepping, and periodic cells, as well as newer methods such as fast extended Lagrangians and high quality electrostatic potential generation. At present M-Chem is a developer friendly environment in which we encourage new software contributors from diverse fields to build their algorithms, models, and methods in our modular framework. The long-term objective of M-Chem is to create an interdisciplinary platform for computational methods with applications ranging from biomolecular simulations, reactive chemistry, to materials research.
Collapse
Affiliation(s)
- Jagna Witek
- Kenneth S. Pitzer Theory Center and Department of Chemistry
| | - Joseph P Heindel
- Kenneth S. Pitzer Theory Center and Department of Chemistry
- Chemical Sciences Division, Lawrence Berkeley National Laboratory
| | - Xingyi Guan
- Kenneth S. Pitzer Theory Center and Department of Chemistry
- Chemical Sciences Division, Lawrence Berkeley National Laboratory
| | - Itai Leven
- Kenneth S. Pitzer Theory Center and Department of Chemistry
| | - Hongxia Hao
- Kenneth S. Pitzer Theory Center and Department of Chemistry
| | | | - Allen LaCour
- Kenneth S. Pitzer Theory Center and Department of Chemistry
- Chemical Sciences Division, Lawrence Berkeley National Laboratory
| | - Selim Sami
- Kenneth S. Pitzer Theory Center and Department of Chemistry
| | - M F S J Menger
- Stratingh Institute for Chemistry, University of Groningen, 9747 AG Groningen, The Netherlands
| | - D Vale Cofer-Shabica
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA, 19128 USA
| | - Eric Berquist
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Shirin Faraji
- Stratingh Institute for Chemistry, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry
- Chemical Sciences Division, Lawrence Berkeley National Laboratory
- Department of Bioengineering and Chemical and Biomolecular Engineering University of California, Berkeley, CA, USA
| |
Collapse
|
6
|
Samanta T, Matyushov DV. Ionic mobility driven by correlated van der Waals and electrostatic forces. J Chem Phys 2022; 156:204501. [DOI: 10.1063/5.0088835] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Classical theories of dielectric friction make two critical assumptions: (i) friction due to van der Waals (vdW) forces is described by hydrodynamic drag and is independent of the ionic charge and (ii) vdW and electrostatic forces are statistically independent. Both assumptions turn out to be incorrect when tested against simulations of anions and cations with varying charge magnitude dissolved in water. Both the vdW and electrostatic components of the force variance scale linearly with the ionic charge squared. The two components are strongly anticorrelated producing simple relations for the total force variance in terms of self-variances. The inverse diffusion constant scales linearly with the charge squared. Solvation asymmetry between cations and anions extends to linear transport coefficients.
Collapse
Affiliation(s)
- Tuhin Samanta
- School of Molecular Sciences and Department of Physics, Arizona State University, P.O. Box 871604, Tempe, Arizona 85287-1604, USA
| | - Dmitry V. Matyushov
- School of Molecular Sciences and Department of Physics, Arizona State University, P.O. Box 871604, Tempe, Arizona 85287-1604, USA
| |
Collapse
|
7
|
Foucaud Y, Siboulet B, Duvail M, Jonchere A, Diat O, Vuilleumier R, Dufrêche JF. Deciphering second harmonic generation signals. Chem Sci 2021; 12:15134-15142. [PMID: 34909155 PMCID: PMC8612378 DOI: 10.1039/d1sc03960a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Accepted: 10/25/2021] [Indexed: 11/21/2022] Open
Abstract
Second harmonic generation (SHG) has emerged as one of the most powerful techniques used to selectively monitor surface dynamics and reactions for all types of interfaces as well as for imaging non-centrosymmetric structures, although the molecular origin of the SHG signal is still poorly understood. Here, we present a breakthrough approach to predict and interpret the SHG signal at the atomic level, which is freed from the hyperpolarisability concept and self-consistently considers the non-locality and the coupling with the environment. The direct ab initio method developed here shows that a bulk quadrupole contribution significantly overwhelms the interface dipole term in the purely interfacial induced second-order polarisation for water/air interfaces. The obtained simulated SHG responses are in unprecedented agreement with the experimental signal. This work not only paves the road for the prediction of SHG response from more complex interfaces of all types, but also suggests new insights in the interpretation of the SHG signal at a molecular level. In particular, it highlights the modest influence of the molecular orientation and the high significance of the bulk quadrupole contribution, which does not depend on the interface, in the total experimental response. Second harmonic generation is one of the most powerful techniques used to selectively probe interfaces of all types. The direct ab initio method developed here allows predicting the signal and highlights the importance of local and non-local effects.![]()
Collapse
Affiliation(s)
- Yann Foucaud
- ICSM, Univ Montpellier, CEA, CNRS, ENSCM Marcoule France
| | | | - Magali Duvail
- ICSM, Univ Montpellier, CEA, CNRS, ENSCM Marcoule France
| | - Alban Jonchere
- ICSM, Univ Montpellier, CEA, CNRS, ENSCM Marcoule France
| | - Olivier Diat
- ICSM, Univ Montpellier, CEA, CNRS, ENSCM Marcoule France
| | - Rodolphe Vuilleumier
- PASTEUR, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS 75005 Paris France
| | | |
Collapse
|
8
|
Shi Y, Doyle CC, Beck TL. Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data. J Phys Chem Lett 2021; 12:10310-10317. [PMID: 34662132 DOI: 10.1021/acs.jpclett.1c02328] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion-solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. Here we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.
Collapse
Affiliation(s)
- Yu Shi
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221-0172, United States
| | - Carrie C Doyle
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221-0172, United States
| | - Thomas L Beck
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221-0172, United States
| |
Collapse
|
9
|
Kathmann SM. Electric fields and potentials in condensed phases. Phys Chem Chem Phys 2021; 23:23836-23849. [PMID: 34647950 DOI: 10.1039/d1cp03571a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The electric fields and potentials inside and at the interface of matter are relevant to many branches of physics, chemistry, and biology. Accurate quantification of these fields and/or potentials is essential to control and exploit chemical and physical transformations. Before we understand the response of matter to external fields, it is first important to understand the intrinsic interior and interfacial fields and potentials, both classically and quantum mechanically, as well as how they are probed experimentally. Here we compare and contrast, beginning with the hydrogen atom in vacuum and ending with concentrated aqueous NaCl electrolyte, both classical and quantum mechanical electric potentials and fields. We make contact with experimental vibrational Stark, electrochemical, X-ray, and electron spectroscopic probes of these potentials and fields, outline relevant conceptual difficulties, and underscore the advantage of electron holography as a basis to better understand electrostatics in matter.
Collapse
Affiliation(s)
- Shawn M Kathmann
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, WA, 99354, USA.
| |
Collapse
|
10
|
Thakur AC, Remsing RC. Distributed charge models of liquid methane and ethane for dielectric effects and solvation. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1933228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Atul C. Thakur
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ, USA
| | - Richard C. Remsing
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ, USA
| |
Collapse
|
11
|
Chuev GN, Fedotova MV, Valiev M. Renormalized site density functional theory for models of ion hydration. J Chem Phys 2021; 155:064501. [PMID: 34391371 DOI: 10.1063/5.0060249] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The development of accurate statistical mechanics models of molecular liquid systems is a problem of great practical and fundamental importance. Site-density functional theory (SDFT) is one of the promising directions in this area, but its success hinges upon the ability to efficiently reconcile the co-existence of two distinct intra- and inter-molecular interaction regimes in a molecular liquid. The renormalized formulation of SDFT (RSDFT), which we have recently developed, resolves this problem by introducing an additional potential field variable that decouples two interaction scales and maps the molecular liquid problem onto the effective simple liquid mixture. This work provides a critical assessment of RSDFT for the hydrated ion system-a problem that historically has always been one of the most difficult cases for SDFT applications. Using a two-site model of water, we perform a comprehensive analysis of hydrated alkali metal and halogen ions, including both structural and free energy based characteristics. The results indicate that RSDFT provides a significant improvement over conventional three-dimensional reference interaction site model implementations and may prove useful in coarse grained simulations based on two-site solvent models.
Collapse
Affiliation(s)
- Gennady N Chuev
- Institute of Theoretical and Experimental Biophysics, Russian Academy of Science, Pushchino, Moscow Region 142290, Russia
| | - Marina V Fedotova
- G. A. Krestov Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
| | - Marat Valiev
- Molecular Sciences Software Group, Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| |
Collapse
|
12
|
Cox SJ, Mandadapu KK, Geissler PL. Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water. J Chem Phys 2021; 154:244502. [PMID: 34241373 DOI: 10.1063/5.0051399] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. We propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.
Collapse
Affiliation(s)
- Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Kranthi K Mandadapu
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Phillip L Geissler
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
13
|
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
|
14
|
Cox SJ, Thorpe DG, Shaffer PR, Geissler PL. Assessing long-range contributions to the charge asymmetry of ion adsorption at the air-water interface. Chem Sci 2020; 11:11791-11800. [PMID: 34094413 PMCID: PMC8162909 DOI: 10.1039/d0sc01947j] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Anions generally associate more favorably with the air–water interface than cations. In addition to solute size and polarizability, the intrinsic structure of the unperturbed interface has been discussed as an important contributor to this bias. Here we assess quantitatively the role that intrinsic charge asymmetry of water's surface plays in ion adsorption, using computer simulations to compare model solutes of various size and charge. In doing so, we also evaluate the degree to which linear response theory for solvent polarization is a reasonable approach for comparing the thermodynamics of bulk and interfacial ion solvation. Consistent with previous works on bulk ion solvation, we find that the average electrostatic potential at the center of a neutral, sub-nanometer solute at the air–water interface depends sensitively on its radius, and that this potential changes quite nonlinearly as the solute's charge is introduced. The nonlinear response closely resembles that of the bulk. As a result, the net nonlinearity of ion adsorption is weaker than in bulk, but still substantial, comparable to the apparent magnitude of macroscopically nonlocal contributions from the undisturbed interface. For the simple-point-charge model of water we study, these results argue distinctly against rationalizing ion adsorption in terms of surface potentials inherent to molecular structure of the liquid's boundary. Cations and anions have different affinities for the air-water interface. The intrinsic orientation of surface molecules suggests such an asymmetry, but the bias is dominated by solvent response that is spatially local and significantly nonlinear.![]()
Collapse
Affiliation(s)
- Stephen J Cox
- Department of Chemistry, University of Cambridge Lensfield Road Cambridge CB2 1EW UK
| | - Dayton G Thorpe
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA.,Department of Physics, University of California Berkeley CA 94720 USA
| | - Patrick R Shaffer
- Department of Chemistry, University of California Berkeley CA 94720 USA
| | - Phillip L Geissler
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA.,Department of Chemistry, University of California Berkeley CA 94720 USA
| |
Collapse
|
15
|
Wineman-Fisher V, Delgado JM, Nagy PR, Jakobsson E, Pandit SA, Varma S. Transferable interactions of Li + and Mg 2+ ions in polarizable models. J Chem Phys 2020; 153:104113. [PMID: 32933310 DOI: 10.1063/5.0022060] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Therapeutic implications of Li+, in many cases, stem from its ability to inhibit certain Mg2+-dependent enzymes, where it interacts with or substitutes for Mg2+. The underlying details of its action are, however, unknown. Molecular simulations can provide insights, but their reliability depends on how well they describe relative interactions of Li+ and Mg2+ with water and other biochemical groups. Here, we explore, benchmark, and recommend improvements to two simulation approaches: the one that employs an all-atom polarizable molecular mechanics (MM) model and the other that uses a hybrid quantum and MM implementation of the quasi-chemical theory (QCT). The strength of the former is that it describes thermal motions explicitly and that of the latter is that it derives local contributions from electron densities. Reference data are taken from the experiment, and also obtained systematically from CCSD(T) theory, followed by a benchmarked vdW-inclusive density functional theory. We find that the QCT model predicts relative hydration energies and structures in agreement with the experiment and without the need for additional parameterization. This implies that accurate descriptions of local interactions are essential. Consistent with this observation, recalibration of local interactions in the MM model, which reduces errors from 10.0 kcal/mol to 1.4 kcal/mol, also fixes aqueous phase properties. Finally, we show that ion-ligand transferability errors in the MM model can be reduced significantly from 10.3 kcal/mol to 1.2 kcal/mol by correcting the ligand's polarization term and by introducing Lennard-Jones cross-terms. In general, this work sets up systematic approaches to evaluate and improve molecular models of ions binding to proteins.
Collapse
Affiliation(s)
- Vered Wineman-Fisher
- Department of Cell Biology, Microbiology and Molecular Biology, University of South Florida, Tampa, Florida 33620, USA
| | - Julián Meléndez Delgado
- Department of Cell Biology, Microbiology and Molecular Biology, University of South Florida, Tampa, Florida 33620, USA
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Eric Jakobsson
- National Center for Supercomputing Applications, Center for Biophysics and Computational Biology, Department of Molecular and Integrative Physiology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Sagar A Pandit
- Department of Physics, University of South Florida, Tampa, Florida 33620, USA
| | - Sameer Varma
- Department of Cell Biology, Microbiology and Molecular Biology, University of South Florida, Tampa, Florida 33620, USA
| |
Collapse
|
16
|
Ghaani MR, Kusalik PG, English NJ. Massive generation of metastable bulk nanobubbles in water by external electric fields. SCIENCE ADVANCES 2020; 6:eaaz0094. [PMID: 32284977 PMCID: PMC7124953 DOI: 10.1126/sciadv.aaz0094] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/05/2019] [Accepted: 12/26/2019] [Indexed: 05/21/2023]
Abstract
Nanobubbles (NBs) are nanoscopic gaseous domains than can exist on solid surfaces or in bulk liquids. They have attracted substantial attention due to their long-time (meta)stability and a high potential for real-world applications. Using an approach not previously investigated, we exploit surface-electrostatic NB formation and stabilization via application of external electric fields in gas-liquid systems, with the marked result of massively increased gas uptake into the liquid in NB form. The de facto gas solubility enhancement (over many months) ranges from 2.5-fold for oxygen to 30-fold for methane vis-à-vis respective Henry's law values for gas solubility; the more hydrophobic the gas, the more spectacular the increase. Molecular dynamics simulations reveal that the origin of NBs' movement lies in dielectrophoresis, while substantial NB stabilization arises from a surface-polarization interaction.
Collapse
Affiliation(s)
- Mohammad Reza Ghaani
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
- Corresponding author. (M.R.G.); (N.J.E.)
| | - Peter G. Kusalik
- Department of Chemistry, University of Calgary, Calgary, Alberta, Canada
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
- Corresponding author. (M.R.G.); (N.J.E.)
| |
Collapse
|
17
|
Yesibolati MN, Laganà S, Sun H, Beleggia M, Kathmann SM, Kasama T, Mølhave K. Mean Inner Potential of Liquid Water. PHYSICAL REVIEW LETTERS 2020; 124:065502. [PMID: 32109081 DOI: 10.1103/physrevlett.124.065502] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/29/2019] [Revised: 11/20/2019] [Accepted: 12/23/2019] [Indexed: 06/10/2023]
Abstract
Improving our experimental and theoretical knowledge of electric potentials at liquid-solid boundaries is essential to achieve a deeper understanding of the driving forces behind interfacial processes. Electron holography has proved successful in probing solid-solid interfaces but requires knowledge of the materials' mean inner potential (MIP, V_{0}), which is a fundamental bulk material property. Combining off-axis electron holography with liquid phase transmission electron microscopy (LPTEM), we provide the first quantitative MIP determination of liquid water V_{0}=+4.48±0.19 V. This value is larger than most theoretical predictions, and to explain the disagreement we assess the dominant factors needed in quantum simulations of liquid water. A precise MIP lays the foundations for nanoscale holographic potential measurements in liquids, and provides a benchmark to improve quantum mechanical descriptions of aqueous systems and their interfaces in, e.g., electrochemistry, solvation processes, and spectroscopy.
Collapse
Affiliation(s)
- Murat Nulati Yesibolati
- DTU Nanolab, National Centre for Nano Fabrication and Characterization, Technical University of Denmark, Building 307, 2800 Kgs. Lyngby, Denmark
| | - Simone Laganà
- DTU Nanolab, National Centre for Nano Fabrication and Characterization, Technical University of Denmark, Building 307, 2800 Kgs. Lyngby, Denmark
| | - Hongyu Sun
- DTU Nanolab, National Centre for Nano Fabrication and Characterization, Technical University of Denmark, Building 307, 2800 Kgs. Lyngby, Denmark
| | - Marco Beleggia
- DTU Nanolab, National Centre for Nano Fabrication and Characterization, Technical University of Denmark, Building 307, 2800 Kgs. Lyngby, Denmark
| | - Shawn M Kathmann
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, USA
| | - Takeshi Kasama
- DTU Nanolab, National Centre for Nano Fabrication and Characterization, Technical University of Denmark, Building 307, 2800 Kgs. Lyngby, Denmark
| | - Kristian Mølhave
- DTU Nanolab, National Centre for Nano Fabrication and Characterization, Technical University of Denmark, Building 307, 2800 Kgs. Lyngby, Denmark
| |
Collapse
|
18
|
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
|
19
|
Houriez C, Réal F, Vallet V, Mautner M, Masella M. Ion hydration free energies and water surface potential in water nano drops: The cluster pair approximation and the proton hydration Gibbs free energy in solution. J Chem Phys 2019; 151:174504. [PMID: 31703526 DOI: 10.1063/1.5109777] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We estimate both single ion hydration Gibbs free energies in water droplets, comprising from 50 to 1000 molecules, and water/vacuum surface potentials in pure water droplets comprising up to 10 000 molecules. We consider four ions, namely, Li+, NH4 +, F-, and Cl-, and we model their hydration process and water/water interactions using polarizable force fields based on an induced point dipole approach. We show both ion hydration Gibbs free energies and water surface potentials to obey linear functions of the droplet radius as soon as droplets comprising a few hundred water molecules. Moreover, we also show that the differences in anion/cation hydration Gibbs free energies in droplets obey a different regime in large droplets than in small clusters comprising no more than six water molecules, in line with the earlier results computed from standard additive point charge force fields. Hence, both point charge and more sophisticated induced point dipole molecular modeling approaches suggest that methods considering only the thermodynamical properties of small ion/water clusters to estimate the absolute proton hydration Gibbs free energy in solution are questionable. In particular, taking into account the data of large ion/water droplets may yield a proton hydration Gibbs free energy in solution value to be shifted by several kBT units compared to small clusters-based approaches.
Collapse
Affiliation(s)
- Céline Houriez
- MINES ParisTech, PSL Research University, CTP - Centre Thermodynamique des Procédés, 35 rue Saint-Honoré, 77300 Fontainebleau, France
| | - Florent Réal
- Univ. Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - Valérie Vallet
- Univ. Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - Michael Mautner
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284-2006, USA and Department of Chemistry, University of Canterbury, Christchurch 8001, New Zealand
| | - Michel Masella
- Laboratoire de Biologie Structurale et Radiobiologie, Service de Bioénergétique, Biologie Structurale et Mécanismes, Institut Joliot, CEA Saclay, F-91191 Gif sur Yvette Cedex, France
| |
Collapse
|
20
|
Seyedi S, Martin DR, Matyushov DV. Screening of Coulomb interactions in liquid dielectrics. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:325101. [PMID: 31042681 DOI: 10.1088/1361-648x/ab1e6f] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The interaction of charges in dielectric materials is screened by the dielectric constant of the bulk dielectric. In dielectric theories, screening is assigned to surface charges appearing from preferential orientations of dipoles along the local field in the interface. For liquid dielectrics, such interfacial orientations are affected by the interfacial structure characterized by a separate interfacial dipolar susceptibility. We argue that dielectric properties of polar liquids should be characterized by two distinct susceptibilities responsible for local response (solvation) and long-range response (dielectric screening). A microscopic model for dipolar screening in polar liquids is developed here. It shows that the standard bulk dielectric constant is responsible for screening at large distances. The potential of mean force between ions in polar liquids becomes oscillatory at short distances. Oscillations arise from the coupling of collective longitudinal excitations (dipolarons) of the polar liquid with its interfacial density profile. Interfacial electrostatics demonstrates a crossover beyond the solute radius of ∼1 nm from the scaling laws roughly consistent with continuum expectations for small solutes to new scaling trends characterizing the much softer nano-scale solute-solvent interface. This crossover also marks the transition to a continuum-type electrostatic screening between ions, when short-distance oscillation of the potential of mean force become strongly damped. Screening oscillations are enhanced for more structured interfaces.
Collapse
Affiliation(s)
- Salman Seyedi
- Department of Physics, Arizona State University, PO Box 871504, Tempe, AZ 85287-1504, United States of America
| | | | | |
Collapse
|
21
|
Das AK, Urban L, Leven I, Loipersberger M, Aldossary A, Head-Gordon M, Head-Gordon T. Development of an Advanced Force Field for Water Using Variational Energy Decomposition Analysis. J Chem Theory Comput 2019; 15:5001-5013. [DOI: 10.1021/acs.jctc.9b00478] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
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
|
Albaugh A, Tuckerman ME, Head-Gordon T. Combining Iteration-Free Polarization with Large Time Step Stochastic-Isokinetic Integration. J Chem Theory Comput 2019; 15:2195-2205. [PMID: 30830768 DOI: 10.1021/acs.jctc.9b00072] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In order to accelerate molecular dynamics simulations using polarizable force fields, we combine a new extended Lagrangian approach that eliminates the self-consistent field step (iEL/0-SCF) with a stochastic integration scheme that allows for a long time step using a multiple time scale algorithm (SIN(R)). We consider different algorithms for the combined scheme that places different components of the nonbonded forces into different time scales, as well as splitting individual nonbonded forces across time scales, to demonstrate that the combined method works well for bulk water as well as for a concentrated salt solution, aqueous peptide, and solvated protein. Depending on system and desired accuracy, the iEL/0-SCF and SIN(R) combination yields lower bound computational speed-ups of ∼6-8 relative to a molecular dynamics Verlet integration using a standard SCF solver implemented in the reference program TINKER 8.1. The combined approach embodies a significant advance for equilibrium simulations in the canonical ensemble of many-body potential energy surfaces for condensed phase systems with speed-ups that exceed what is possible by either method alone.
Collapse
Affiliation(s)
| | - Mark E Tuckerman
- NYU-ECNU , Center for Computational Chemistry at NYU, Shanghai , Shanghai 200062 , China
| | | |
Collapse
|
24
|
Zdrali E, Baer MD, Okur HI, Mundy CJ, Roke S. The Diverse Nature of Ion Speciation at the Nanoscale Hydrophobic/Water Interface. J Phys Chem B 2019; 123:2414-2423. [DOI: 10.1021/acs.jpcb.8b10207] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Evangelia Zdrali
- Laboratory for fundamental BioPhotonics, Institutes of Bioengineering and Materials Science and Engineering, School of Engineering, and Lausanne Centre for Ultrafast Science, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Marcel D. Baer
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Halil I. Okur
- Laboratory for fundamental BioPhotonics, Institutes of Bioengineering and Materials Science and Engineering, School of Engineering, and Lausanne Centre for Ultrafast Science, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Christopher J. Mundy
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| | - Sylvie Roke
- Laboratory for fundamental BioPhotonics, Institutes of Bioengineering and Materials Science and Engineering, School of Engineering, and Lausanne Centre for Ultrafast Science, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
25
|
Remsing RC, Klein ML. Exponential Scaling of Water Exchange Rates with Ion Interaction Strength from the Perspective of Dynamic Facilitation Theory. J Phys Chem A 2019; 123:1077-1084. [PMID: 30609371 DOI: 10.1021/acs.jpca.8b09667] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Richard C. Remsing
- Institute for Computational Molecular Science and Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Michael L. Klein
- Institute for Computational Molecular Science and Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
26
|
Prasetyo N, Hünenberger PH, Hofer TS. Single-Ion Thermodynamics from First Principles: Calculation of the Absolute Hydration Free Energy and Single-Electrode Potential of Aqueous Li + Using ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations. J Chem Theory Comput 2018; 14:6443-6459. [PMID: 30284829 DOI: 10.1021/acs.jctc.8b00729] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
A recently proposed thermodynamic integration (TI) approach formulated in the framework of quantum mechanical/molecular mechanical molecular dynamics (QM/MM MD) simulations is applied to study the structure, dynamics, and absolute intrinsic hydration free energy Δs GM+,wat◦ of the Li+ ion at a correlated ab initio level of theory. Based on the results, standard values (298.15 K, ideal gas at 1 bar, ideal solute at 1 molal) for the absolute intrinsic hydration free energy [Formula: see text] of the proton, the surface electric potential jump χwat◦ upon entering bulk water, and the absolute single-electrode potential [Formula: see text] of the reference hydrogen electrode are calculated to be -1099.9 ± 4.2 kJ·mol-1, 0.13 ± 0.08 V, and 4.28 ± 0.04 V, respectively, in excellent agreement with the standard values recommended by Hünenberger and Reif on the basis of an extensive evaluation of the available experimental data (-1100 ± 5 kJ·mol-1, 0.13 ± 0.10 V, and 4.28 ± 0.13 V). The simulation results for Li+ are also compared to those for Na+ and K+ from a previous study in terms of relative hydration free energies ΔΔs GM+,wat◦ and relative electrode potentials [Formula: see text]. The calculated values are found to agree extremely well with the experimental differences in standard conventional hydration free energies ΔΔs GM+,wat• and redox potentials [Formula: see text]. The level of agreement between simulation and experiment, which is quantitative within error bars, underlines the substantial accuracy improvement achieved by applying a highly demanding QM/MM approach at the resolution-of-identity second-order Møller-Plesset perturbation (RIMP2) level over calculations relying on purely molecular mechanical or density functional theory (DFT) descriptions. A detailed analysis of the structural and dynamical properties of the Li+ hydrate indicates that a correct description of the solvation structure and dynamics is achieved as well at this level of theory. Consideration of the QM/MM potential-energy components also shows that the partitioning into QM and MM zones does not induce any significant energetic artifact for the system considered.
Collapse
Affiliation(s)
- Niko Prasetyo
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry , University of Innsbruck , Innrain 80-82 , A-6020 Innsbruck , Austria.,Austria-Indonesia Centre (AIC) for Computational Chemistry , Universitas Gadjah Mada , Sekip Utara , Yogyakarta 55281 , Indonesia.,Department of Chemistry, Faculty of Mathematics and Natural Sciences , Universitas Gadjah Mada , Sekip Utara , Yogyakarta 55281 , Indonesia
| | - Philippe H Hünenberger
- Laboratorium für Physikalische Chemie , ETH Zürich, ETH-Hönggerberg , HCI Building , CH-8093 Zürich , Switzerland
| | - Thomas S Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry , University of Innsbruck , Innrain 80-82 , A-6020 Innsbruck , Austria
| |
Collapse
|
27
|
Ruiz Pestana L, Marsalek O, Markland TE, Head-Gordon T. The Quest for Accurate Liquid Water Properties from First Principles. J Phys Chem Lett 2018; 9:5009-5016. [PMID: 30118601 DOI: 10.1021/acs.jpclett.8b02400] [Citation(s) in RCA: 53] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Developing accurate ab initio molecular dynamics (AIMD) models that capture both electronic reorganization and nuclear quantum effects associated with hydrogen bonding is key to quantitative understanding of bulk water and its anomalies as well as its role as a universal solvent. For condensed phase simulations, AIMD has typically relied on the generalized gradient approximation (GGA) of density functional theory (DFT) as the underlying model chemistry for the potential energy surface, with nuclear quantum effects (NQEs) sometimes modeled by performing classical molecular dynamics simulations at elevated temperatures. Here we show that the properties of liquid water obtained from the meta-GGA B97M-rV functional, when evaluated using accelerated path integral molecular dynamics simulations, display accuracy comparable to a computationally expensive dispersion-corrected hybrid functional, revPBE0-D3. We show that the meta-GGA DFT functional reproduces bulk water properties including radial distribution functions, self-diffusion coefficients, and infrared spectra with comparable accuracy of a much more expensive hybrid functional. This work demonstrates that the underlying quality of a good DFT functional requires evaluation with quantum nuclei and that high-temperature simulations are a poor proxy for properly treating NQEs.
Collapse
Affiliation(s)
| | - Ondrej Marsalek
- Charles University , Faculty of Mathematics and Physics , Ke Karlovu 3 , 121 16 Prague 2 , Czech Republic
| | - Thomas E Markland
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | | |
Collapse
|
28
|
Duignan TT, Baer MD, Mundy CJ. Understanding the scale of the single ion free energy: A critical test of the tetra-phenyl arsonium and tetra-phenyl borate assumption. J Chem Phys 2018; 148:222819. [PMID: 29907030 DOI: 10.1063/1.5020171] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The tetra-phenyl arsonium and tetra-phenyl borate (TATB) assumption is a commonly used extra-thermodynamic assumption that allows single ion free energies to be split into cationic and anionic contributions. The assumption is that the values for the TATB salt can be divided equally. This is justified by arguing that these large hydrophobic ions will cause a symmetric response in water. Experimental and classical simulation work has raised potential flaws with this assumption, indicating that hydrogen bonding with the phenyl ring may favor the solvation of the TB- anion. Here, we perform ab initio molecular dynamics simulations of these ions in bulk water demonstrating that there are significant structural differences. We quantify our findings by reproducing the experimentally observed vibrational shift for the TB- anion and confirm that this is associated with hydrogen bonding with the phenyl rings. Finally, we demonstrate that this results in a substantial energetic preference of the water to solvate the anion. Our results suggest that the validity of the TATB assumption, which is still widely used today, should be reconsidered experimentally in order to properly reference single ion solvation free energy, enthalpy, and entropy.
Collapse
Affiliation(s)
- Timothy T Duignan
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99354, USA
| | - Marcel D Baer
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99354, USA
| | - Christopher J Mundy
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99354, USA
| |
Collapse
|
29
|
Pollard TP, Beck TL. Re-examining the tetraphenyl-arsonium/tetraphenyl-borate (TATB) hypothesis for single-ion solvation free energies. J Chem Phys 2018; 148:222830. [DOI: 10.1063/1.5024209] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Travis P. Pollard
- Electrochemistry Branch, US Army Research Laboratory, Adelphi, Maryland 20852, USA
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, USA
| | - Thomas L. Beck
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, USA
- Department of Physics, University of Cincinnati, Cincinnati, Ohio 45221, USA
| |
Collapse
|
30
|
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
|
31
|
Boateng HA. Mesh-free hierarchical clustering methods for fast evaluation of electrostatic interactions of point multipoles. J Chem Phys 2018; 147:164104. [PMID: 29096477 DOI: 10.1063/1.4990552] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Electrostatic interactions involving point multipoles are being increasingly implemented to achieve higher accuracy in molecular simulations. A major drawback of multipolar electrostatics is the increased computational cost. Here we develop and compare two Cartesian tree algorithms which employ Taylor approximations and hierarchical clustering to speed up the evaluation of point multipole interactions. We present results from applying the algorithms to compute the free space Coulomb potential and forces of different sets of interacting point multipoles with different densities. The methods achieve high accuracy and speedup of more than an order of magnitude over direct sum calculations and scale well in parallel.
Collapse
Affiliation(s)
- H A Boateng
- Department of Mathematics, Bates College, 2 Andrews Rd., Lewiston, Maine 04240, USA
| |
Collapse
|
32
|
Xi E, Marks SM, Fialoke S, Patel AJ. Sparse sampling of water density fluctuations near liquid-vapor coexistence. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1457218] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Erte Xi
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania , Philadelphia, PA, USA
| | - Sean M. Marks
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania , Philadelphia, PA, USA
| | - Suruchi Fialoke
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania , Philadelphia, PA, USA
| | - Amish J. Patel
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania , Philadelphia, PA, USA
| |
Collapse
|
33
|
Albaugh A, Boateng HA, Bradshaw RT, Demerdash ON, Dziedzic J, Mao Y, Margul DT, Swails J, Zeng Q, Case DA, Eastman P, Wang LP, Essex JW, Head-Gordon M, Pande VS, Ponder JW, Shao Y, Skylaris CK, Todorov IT, Tuckerman ME, Head-Gordon T. Advanced Potential Energy Surfaces for Molecular Simulation. J Phys Chem B 2018; 120:9811-32. [PMID: 27513316 DOI: 10.1021/acs.jpcb.6b06414] [Citation(s) in RCA: 66] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Advanced potential energy surfaces are defined as theoretical models that explicitly include many-body effects that transcend the standard fixed-charge, pairwise-additive paradigm typically used in molecular simulation. However, several factors relating to their software implementation have precluded their widespread use in condensed-phase simulations: the computational cost of the theoretical models, a paucity of approximate models and algorithmic improvements that can ameliorate their cost, underdeveloped interfaces and limited dissemination in computational code bases that are widely used in the computational chemistry community, and software implementations that have not kept pace with modern high-performance computing (HPC) architectures, such as multicore CPUs and modern graphics processing units (GPUs). In this Feature Article we review recent progress made in these areas, including well-defined polarization approximations and new multipole electrostatic formulations, novel methods for solving the mutual polarization equations and increasing the MD time step, combining linear-scaling electronic structure methods with new QM/MM methods that account for mutual polarization between the two regions, and the greatly improved software deployment of these models and methods onto GPU and CPU hardware platforms. We have now approached an era where multipole-based polarizable force fields can be routinely used to obtain computational results comparable to state-of-the-art density functional theory while reaching sampling statistics that are acceptable when compared to that obtained from simpler fixed partial charge force fields.
Collapse
Affiliation(s)
| | - Henry A Boateng
- Department of Mathematics, Bates College , 2 Andrews Road, Lewiston, Maine 04240, United States
| | - Richard T Bradshaw
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, United Kingdom
| | | | - Jacek Dziedzic
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, United Kingdom.,Faculty of Applied Physics and Mathematics, Gdansk University of Technology , 80-223 Gdansk, Poland
| | | | | | - Jason Swails
- Department of Chemistry and Chemical Biology, Rutgers University , Piscataway, New Jersey 08854-8066, United States
| | - Qiao Zeng
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| | - David A Case
- Department of Chemistry and Chemical Biology, Rutgers University , Piscataway, New Jersey 08854-8066, United States
| | - Peter Eastman
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| | - Lee-Ping Wang
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| | - Jonathan W Essex
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, United Kingdom
| | | | - Vijay S Pande
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| | - Jay W Ponder
- Department of Chemistry, Washington University in St. Louis , St. Louis, Missouri 63130, United States
| | - Yihan Shao
- Q-Chem Inc. , 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, United Kingdom
| | - Ilian T Todorov
- STFC Daresbury Laboratory , Keckwick Lane, Daresbury, Warrington WA4 4AD, United Kingdom
| | - Mark E Tuckerman
- NYU-ECNU, Center for Computational Chemistry at NYU, Shanghai , Shanghai 200062, China
| | | |
Collapse
|
34
|
Wise PK, Ben-Amotz D. Interfacial Adsorption of Neutral and Ionic Solutes in a Water Droplet. J Phys Chem B 2018; 122:3447-3453. [PMID: 29244951 DOI: 10.1021/acs.jpcb.7b10488] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Direct (solute-water) and indirect (water-water) contributions to adsorption at an air-water interface are identified using the Widom potential distribution theorem and quantified using molecular dynamics simulations of a liquid water droplet containing either neopentane or iodide-like solutes with charges of 0 or ±1. The results are used to quantitatively compare direct and indirect energetic and entropic contributions to adsorption, as well as to critically test surface capillary wave, linear response (LR), and mean field (MF) predictions. The negative signs of the total adsorption energies and entropies of both the anionic and cationic solutes are found to result from indirect adsorption induced changes in water-water interactions, rather than from surface capillary wave perturbations, which are found to be asymmetric with respect to solute charge. The LR and MF approximations both accurately describe the adsorption of neutral (hydrophobic) solutes, while for ionic solutes the MF approximation is entirely inappropriate and LR predictions are qualitatively (but not quantitatively) accurate.
Collapse
Affiliation(s)
- Patrick K Wise
- Department of Chemistry , Purdue University , West Lafayette , Indiana 47907 United States
| | - Dor Ben-Amotz
- Department of Chemistry , Purdue University , West Lafayette , Indiana 47907 United States
| |
Collapse
|
35
|
Remsing RC, Xi E, Patel AJ. Protein Hydration Thermodynamics: The Influence of Flexibility and Salt on Hydrophobin II Hydration. J Phys Chem B 2018; 122:3635-3646. [DOI: 10.1021/acs.jpcb.7b12060] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Richard C. Remsing
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Erte Xi
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Amish J. Patel
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
36
|
Remsing RC, Duignan TT, Baer MD, Schenter GK, Mundy CJ, Weeks JD. Water Lone Pair Delocalization in Classical and Quantum Descriptions of the Hydration of Model Ions. J Phys Chem B 2018; 122:3519-3527. [DOI: 10.1021/acs.jpcb.7b10722] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Richard C. Remsing
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Timothy T. Duignan
- Chemical and Materials Science Division, Pacific Northwest National Laboratory, Richland, Washington, United States
| | - Marcel D. Baer
- Chemical and Materials Science Division, Pacific Northwest National Laboratory, Richland, Washington, United States
| | - Gregory K. Schenter
- Chemical and Materials Science Division, Pacific Northwest National Laboratory, Richland, Washington, United States
| | - Christopher J. Mundy
- Chemical and Materials Science Division, Pacific Northwest National Laboratory, Richland, Washington, United States
- Affiliate Professor, Department of Chemical Engineering, University of Washington, Seattle, Washington, 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
|
37
|
Zhong J, Zhu C, Li L, Richmond GL, Francisco JS, Zeng XC. Interaction of SO2 with the Surface of a Water Nanodroplet. J Am Chem Soc 2017; 139:17168-17174. [DOI: 10.1021/jacs.7b09900] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jie Zhong
- Department
of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
| | - Chongqin Zhu
- Department
of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
| | - Lei Li
- Department
of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
| | | | - Joseph S. Francisco
- Department
of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
| | - Xiao Cheng Zeng
- Department
of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
| |
Collapse
|
38
|
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
|
39
|
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
|
40
|
Galib M, Duignan TT, Misteli Y, Baer MD, Schenter GK, Hutter J, Mundy CJ. Mass density fluctuations in quantum and classical descriptions of liquid water. J Chem Phys 2017; 146:244501. [DOI: 10.1063/1.4986284] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Mirza Galib
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Timothy T. Duignan
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Yannick Misteli
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Marcel D. Baer
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Gregory K. Schenter
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Jürg Hutter
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Christopher J. Mundy
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| |
Collapse
|
41
|
Molavi Tabrizi A, Goossens S, Mehdizadeh Rahimi A, Cooper CD, Knepley MG, Bardhan JP. Extending the Solvation-Layer Interface Condition Continum Electrostatic Model to a Linearized Poisson–Boltzmann Solvent. J Chem Theory Comput 2017; 13:2897-2914. [DOI: 10.1021/acs.jctc.6b00832] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Amirhossein Molavi Tabrizi
- Department
of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Spencer Goossens
- Department
of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Ali Mehdizadeh Rahimi
- Department
of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Christopher D. Cooper
- Departamento
de Ingeniería Mecánica and Centro Científico
Tecnológico de Valparaíso (CCTVal), Universidad Técnica Federico Santa María, Valparaiso, Chile
| | - Matthew G. Knepley
- Department
of Computational and Applied Mathematics, Rice University, Houston, Texas 77005, United States
| | - Jaydeep P. Bardhan
- Department
of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| |
Collapse
|
42
|
Marsalek O, Markland TE. Quantum Dynamics and Spectroscopy of Ab Initio Liquid Water: The Interplay of Nuclear and Electronic Quantum Effects. J Phys Chem Lett 2017; 8:1545-1551. [PMID: 28296422 DOI: 10.1021/acs.jpclett.7b00391] [Citation(s) in RCA: 137] [Impact Index Per Article: 19.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Understanding the reactivity and spectroscopy of aqueous solutions at the atomistic level is crucial for the elucidation and design of chemical processes. However, the simulation of these systems requires addressing the formidable challenges of treating the quantum nature of both the electrons and nuclei. Exploiting our recently developed methods that provide acceleration by up to 2 orders of magnitude, we combine path integral simulations with on-the-fly evaluation of the electronic structure at the hybrid density functional theory level to capture the interplay between nuclear quantum effects and the electronic surface. Here we show that this combination provides accurate structure and dynamics, including the full infrared and Raman spectra of liquid water. This allows us to demonstrate and explain the failings of lower-level density functionals for dynamics and vibrational spectroscopy when the nuclei are treated quantum mechanically. These insights thus provide a foundation for the reliable investigation of spectroscopy and reactivity in aqueous environments.
Collapse
Affiliation(s)
- Ondrej Marsalek
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| |
Collapse
|
43
|
Demerdash O, Head-Gordon T. Parallel implementation of approximate atomistic models of the AMOEBA polarizable model. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.10.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
44
|
Hoehn RD, Carignano MA, Kais S, Zhu C, Zhong J, Zeng XC, Francisco JS, Gladich I. Hydrogen bonding and orientation effects on the accommodation of methylamine at the air-water interface. J Chem Phys 2016; 144:214701. [DOI: 10.1063/1.4950951] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
|
45
|
Reif MM, Hünenberger PH. Origin of Asymmetric Solvation Effects for Ions in Water and Organic Solvents Investigated Using Molecular Dynamics Simulations: The Swain Acity-Basity Scale Revisited. J Phys Chem B 2016; 120:8485-517. [PMID: 27173101 DOI: 10.1021/acs.jpcb.6b02156] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The asymmetric solvation of ions can be defined as the tendency of a solvent to preferentially solvate anions over cations or cations over anions, at identical ionic charge magnitudes and effective sizes. Taking water as a reference, these effects are quantified experimentally for many solvents by the relative acity (A) and basity (B) parameters of the Swain scale. The goal of the present study is to investigate the asymmetric solvation of ions using molecular dynamics simulations, and to connect the results to this empirical scale. To this purpose, the charging free energies of alkali and halide ions, and of their hypothetical oppositely charged counterparts, are calculated in a variety of solvents. In a first set of calculations, artificial solvent models are considered that present either a charge or a shape asymmetry at the molecular level. The solvation asymmetry, probed by the difference in charging free energy between the two oppositely charged ions, is found to encompass a term quadratic in the ion charge, related to the different solvation structures around the anion and cation, and a term linear in the ion charge, related to the solvation structure around the uncharged ion-sized cavity. For these simple solvent models, the two terms are systematically counteracting each other, and it is argued that only the quadratic term should be retained when comparing the results of simulations involving physical solvents to experimental data. In a second set of calculations, 16 physical solvents are considered. The theoretical estimates for the acity A are found to correlate very well with the Swain parameters, whereas the correlation for B is very poor. Based on this observation, the Swain scale is reformulated into a new scale involving an asymmetry parameter Σ, positive for acitic solvents and negative for basitic ones, and a polarity parameter Π. This revised scale has the same predictive power as the original scale, but it characterizes asymmetry in an absolute sense, the atomistic simulations playing the role of an extra-thermodynamic assumption, and is optimally compatible with the simulation results. Considering the 55 solvents in the Swain set, it is observed that a moderate basity (Σ between -0.9 and -0.3, related to electronic polarization) represents the baseline for most solvents, while a highly variable acity (Σ between 0.0 and 3.0, related to hydrogen-bond donor capacity modulated by inductive effects) represents a landmark of protic solvents.
Collapse
Affiliation(s)
- Maria M Reif
- Physics Department (T38), Technische Universität München , D-85748 Garching, Germany
| | | |
Collapse
|
46
|
Toward a quantitative theory of Hofmeister phenomena: From quantum effects to thermodynamics. Curr Opin Colloid Interface Sci 2016. [DOI: 10.1016/j.cocis.2016.06.015] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
47
|
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]
|
48
|
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
|
49
|
Baer MD, Mundy CJ. Local Aqueous Solvation Structure Around Ca2+ During Ca2+···Cl– Pair Formation. J Phys Chem B 2016; 120:1885-93. [DOI: 10.1021/acs.jpcb.5b09579] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marcel D. Baer
- Physical Science Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christopher J. Mundy
- Physical Science Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
50
|
Misin M, Fedorov MV, Palmer DS. Hydration Free Energies of Molecular Ions from Theory and Simulation. J Phys Chem B 2016; 120:975-83. [PMID: 26756333 DOI: 10.1021/acs.jpcb.5b10809] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a theoretical/computational framework for accurate calculation of hydration free energies of ionized molecular species. The method is based on a molecular theory, 3D-RISM, combined with a recently developed pressure correction (PC+). The 3D-RISM/PC+ model can provide ∼3 kcal/mol hydration free energy accuracy for a large variety of ionic compounds, provided that the Galvani potential of water is taken into account. The results are compared with direct atomistic simulations. Several methodological aspects of hydration free energy calculations for charged species are discussed.
Collapse
Affiliation(s)
| | | | - David S Palmer
- Department of Pure and Applied Chemistry, University of Strathclyde , 295 Cathedral Street, Glasgow, G1 1XL, United Kingdom
| |
Collapse
|