1
|
Bui AT, Cox SJ. A classical density functional theory for solvation across length scales. J Chem Phys 2024; 161:104103. [PMID: 39248237 DOI: 10.1063/5.0223750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2024] [Accepted: 08/14/2024] [Indexed: 09/10/2024] Open
Abstract
A central aim of multiscale modeling is to use results from the Schrödinger equation to predict phenomenology on length scales that far exceed those of typical molecular correlations. In this work, we present a new approach rooted in classical density functional theory (cDFT) that allows us to accurately describe the solvation of apolar solutes across length scales. Our approach builds on the Lum-Chandler-Weeks (LCW) theory of hydrophobicity [K. Lum et al., J. Phys. Chem. B 103, 4570 (1999)] by constructing a free energy functional that uses a slowly varying component of the density field as a reference. From a practical viewpoint, the theory we present is numerically simpler and generalizes to solutes with soft-core repulsion more easily than LCW theory. Furthermore, by assessing the local compressibility and its critical scaling behavior, we demonstrate that our LCW-style cDFT approach contains the physics of critical drying, which has been emphasized as an essential aspect of hydrophobicity by recent theories. As our approach is parameterized on the two-body direct correlation function of the uniform fluid and the liquid-vapor surface tension, it straightforwardly captures the temperature dependence of solvation. Moreover, we use our theory to describe solvation at a first-principles level on length scales that vastly exceed what is accessible to molecular simulations.
Collapse
Affiliation(s)
- Anna T Bui
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
2
|
Zhao R, Zou Z, Weeks JD, Tiwary P. Quantifying the Relevance of Long-Range Forces for Crystal Nucleation in Water. J Chem Theory Comput 2023; 19:9093-9101. [PMID: 38084039 DOI: 10.1021/acs.jctc.3c01120] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
Understanding nucleation from aqueous solutions is of fundamental importance in a multitude of fields, ranging from materials science to biophysics. The complex solvent-mediated interactions in aqueous solutions hamper the development of a simple physical picture, elucidating the roles of different interactions in nucleation processes. In this work, we make use of three complementary techniques to disentangle the role played by short- and long-range interactions in solvent-mediated nucleation. Specifically, the first approach we utilize is the local molecular field (LMF) theory to renormalize long-range Coulomb electrostatics. Second, we use well-tempered metadynamics to speed up rare events governed by short-range interactions. Third, the deep learning-based State Predictive Information Bottleneck approach is employed in analyzing the reaction coordinate of the nucleation processes obtained from the LMF treatment coupled with well-tempered metadynamics. We find that the two-step nucleation mechanism can largely be captured by the short-range interactions, while the long-range interactions further contribute to the stability of the primary crystal state under ambient conditions. Furthermore, by analyzing the reaction coordinate obtained from the combined LMF-metadynamics treatment, we discern the fluctuations on different time scales, highlighting the need for long-range interactions when accounting for metastability.
Collapse
Affiliation(s)
- Renjie Zhao
- Chemical Physics Program and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | - Ziyue Zou
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, 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
| | - Pratyush Tiwary
- Institute for Physical Science and Technology and Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| |
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
|
Peterson BN, Alfieri ME, Hood DJ, Hettwer CD, Costantino DV, Tabor DP, Kidwell NM. Solvent-Mediated Charge Transfer Dynamics of a Model Brown Carbon Aerosol Chromophore: Photophysics of 1-Phenylpyrrole Induced by Water Solvation. J Phys Chem A 2022; 126:4313-4325. [PMID: 35776530 DOI: 10.1021/acs.jpca.2c00585] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Nitrogen heterocycles are known to be important light-absorbing chromophores in a newly discovered class of aerosols, commonly referred to as "brown carbon" (BrC) aerosols. Due to their significant absorption and spectral overlap with the solar actinic flux, these BrC chromophores steer the physical and optical properties of aerosols. To model the local aqueous solvation environment surrounding BrC chromophores, we generated cold molecular complexes with water and a prototypical BrC chromophore, 1-phenylpyrrole (1PhPy), using supersonic jet-cooling and explored their intermolecular interactions using single-conformation spectroscopy. Herein, we utilized resonant two-photon ionization (R2PI) and UV holeburning (UV HB) double-resonance spectroscopies to obtain a molecular-level understanding of the role of water microsolvation in charge transfer upon photoexcitation of 1PhPy. Quantum chemical calculations and one-dimensional discrete variable representation simulations revealed insights into the charge transfer efficacy of 1PhPy with and without addition of a single water molecule. Taken together, our results indicate that the intermolecular interactions with water guide the geometry of 1PhPy to adopt a more twisted intramolecular charge transfer (TICT) configuration, thus facilitating charge transfer from the pyrrole donor to the phenyl ring acceptor. Furthermore, the water network surrounding 1PhPy reports on the charge transfer such that the H2O solvent primarily interacts with the pyrrole ring donor in the ground state, whereas it preferentially interacts with the phenyl ring acceptor in the excited state. Large Franck-Condon activity is evident in the 1PhPy + 1H2O excitation spectrum for the water-migration vibronic bands, supporting H2O solvent reorganization upon excitation of the 1PhPy chromophore. Fluorescence measurements with increasing H2O % volume corroborated our gas-phase studies by indicating that a polar water solvation environment stabilizes the TICT configuration of 1PhPy in the excited electronic state, from which emission is observed at a lower energy compared to the locally excited configuration.
Collapse
Affiliation(s)
- Brianna N Peterson
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Megan E Alfieri
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - David J Hood
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Christian D Hettwer
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Daniel V Costantino
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Daniel P Tabor
- Department of Chemistry, Texas A&M University, College Station, Texas 77843, United States
| | - Nathanael M Kidwell
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| |
Collapse
|
5
|
Gao A, Remsing RC. Self-consistent determination of long-range electrostatics in neural network potentials. Nat Commun 2022; 13:1572. [PMID: 35322046 PMCID: PMC8943018 DOI: 10.1038/s41467-022-29243-2] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Accepted: 03/07/2022] [Indexed: 12/19/2022] Open
Abstract
Machine learning has the potential to revolutionize the field of molecular simulation through the development of efficient and accurate models of interatomic interactions. Neural networks can model interactions with the accuracy of quantum mechanics-based calculations, but with a fraction of the cost, enabling simulations of large systems over long timescales. However, implicit in the construction of neural network potentials is an assumption of locality, wherein atomic arrangements on the nanometer-scale are used to learn interatomic interactions. Because of this assumption, the resulting neural network models cannot describe long-range interactions that play critical roles in dielectric screening and chemical reactivity. Here, we address this issue by introducing the self-consistent field neural network - a general approach for learning the long-range response of molecular systems in neural network potentials that relies on a physically meaningful separation of the interatomic interactions - and demonstrate its utility by modeling liquid water with and without applied fields.
Collapse
Affiliation(s)
- Ang Gao
- Department of Physics, Beijing University of Posts and Telecommunications, 100876, Beijing, China.
| | - Richard C Remsing
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ, 08854, USA.
| |
Collapse
|
6
|
Wang D, Zhao R, Weeks JD, Tiwary P. Influence of Long-Range Forces on the Transition States and Dynamics of NaCl Ion-Pair Dissociation in Water. J Phys Chem B 2022; 126:545-551. [PMID: 34985884 DOI: 10.1021/acs.jpcb.1c09454] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We study NaCl ion-pair dissociation in a dilute aqueous solution using computer simulations both for the full system with long-range Coulomb interactions and for a well-chosen reference system with short-range intermolecular interactions. Analyzing results using concepts from Local Molecular Field (LMF) theory and the recently proposed AI-based analysis tool "State predictive information bottleneck" (SPIB), we show that the system with short-range interactions can accurately reproduce the transition rate for the dissociation process, the dynamics for moving between the underlying metastable states, and the transition state ensemble. Contributions from long-range interactions can be largely neglected for these processes because long-range forces from the direct interionic Coulomb interactions are almost completely canceled (>90%) by those from solvent interactions over the length scale where the transition takes place. Thus, for this important monovalent ion-pair system, short-range forces alone are able to capture detailed consequences of the collective solvent motion, allowing the use of physically suggestive and computationally efficient short-range models for the dissociation event. We believe that the framework here should be applicable to disentangling mechanisms for more complex processes such as multivalent ion disassociation, where previous work has suggested that long-range contributions may be more important.
Collapse
Affiliation(s)
- Dedi Wang
- Biophysics Program and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | - Renjie Zhao
- Chemical Physics Program and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | - John D Weeks
- Department of Chemistry and Biochemistry and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | - Pratyush Tiwary
- Department of Chemistry and Biochemistry and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
7
|
Niblett SP, Galib M, Limmer DT. Learning intermolecular forces at liquid-vapor interfaces. J Chem Phys 2021; 155:164101. [PMID: 34717371 DOI: 10.1063/5.0067565] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
By adopting a perspective informed by contemporary liquid-state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. We find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces but typically fail for properties that depend on unbalanced long-ranged interactions that build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly varying long-ranged interactions and training neural networks only on the short-ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid, where the exact long-ranged interaction is known, and in an ab initio water model, where it is approximated.
Collapse
Affiliation(s)
- Samuel P Niblett
- Department of Chemistry, University of California, Berkeley California 94609, USA
| | - Mirza Galib
- Department of Chemistry, University of California, Berkeley California 94609, USA
| | - David T Limmer
- Department of Chemistry, University of California, Berkeley California 94609, USA
| |
Collapse
|
8
|
Yonetani Y. Dielectric continuum model examination of real-space electrostatic treatments. J Chem Phys 2021; 154:044103. [PMID: 33514106 DOI: 10.1063/5.0033053] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Electrostatic interaction is long ranged; thus, the accurate calculation is not an easy task in molecular dynamics or Monte Carlo simulations. Though the rigorous Ewald method based on the reciprocal space has been established, real-space treatments have recently become an attractive alternative because of the efficient calculation. However, the construction is not yet completed and is now a challenging subject. In an earlier theoretical study, Neumann and Steinhauser employed the Onsager dielectric continuum model to explain how simple real-space cutoff produces artificial dipolar orientation. In the present study, we employ this continuum model to explore the fundamental properties of the recently developed real-space treatments of three shifting schemes. The result of the distance-dependent Kirkwood function GK(R) showed that the simple bare cutoff produces a well-known hole-shaped artifact, whereas the shift treatments do not. Two-dimensional mapping of electric field well explained how these shift treatments remove the hole-shaped artifact. Still, the shift treatments are not sufficient because they do not produce a flat GK(R) profile unlike ideal no-cutoff treatment. To test the continuum model results, we also performed Monte Carlo simulations of dipolar particles. The results found that the continuum model could predict the qualitative tendency as to whether each electrostatic treatment produces the hole-shaped artifact of GK(R) or not. We expect that the present study using the continuum model offers a stringent criterion to judge whether the primitive electrostatic behavior is correctly described or not, which will be useful for future construction of electrostatic treatments.
Collapse
Affiliation(s)
- Yoshiteru Yonetani
- Quantum Beam Science Research Directorate and Institute for Quantum Life Science, National Institutes for Quantum and Radiological Science and Technology, 8-1-7 Umemidai, Kizugawa, Kyoto 619-0215, Japan
| |
Collapse
|
9
|
Abstract
The dielectric nature of polar liquids underpins much of their ability to act as useful solvents, but its description is complicated by the long-ranged nature of dipolar interactions. This is particularly pronounced under the periodic boundary conditions commonly used in molecular simulations. In this article, the dielectric properties of a water model whose intermolecular electrostatic interactions are entirely short-ranged are investigated. This is done within the framework of local molecular-field theory (LMFT), which provides a well-controlled mean-field treatment of long-ranged electrostatics. This short-ranged model gives a remarkably good performance on a number of counts, and its apparent shortcomings are readily accounted for. These results not only lend support to LMFT as an approach for understanding solvation behavior, but also are relevant to those developing interaction potentials based on local descriptions of liquid structure.
Collapse
|
10
|
Baker EB, Rodgers JM, Weeks JD. Local Molecular Field Theory for Nonequilibrium Systems. J Phys Chem B 2020; 124:5676-5684. [PMID: 32511933 DOI: 10.1021/acs.jpcb.0c03295] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
We provide a framework for extending equilibrium local molecular field (LMF) theory to a statistical ensemble evolving under a time-dependent applied field. In this context, the self-consistency of the original LMF equation is achieved dynamically, which provides an efficient method for computing the equilibrium LMF potential, in addition to providing the nonequilibrium generalization. As a concrete example, we investigate water confined between hydrophobic or charged walls, systems that are very sensitive to the treatment of long-ranged electrostatics. We then analyze confined water in the presence of a time-dependent applied electric field, generated by a sinusoidal or abrupt variation of the magnitudes of uniform charge densities on each wall. Very accurate results are found from the time-dependent LMF formalism even for strong static fields and for time-dependent systems that are driven far from equilibrium where linear response methods fail.
Collapse
Affiliation(s)
- Edward B Baker
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | | | - John D Weeks
- Department of Chemistry and Biochemistry and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
11
|
Budkov YA, Kolesnikov AL. Nonlocal density functional theory of water taking into account many-body dipole correlations: binodal and surface tension of 'liquid-vapour' interface. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:365001. [PMID: 32272457 DOI: 10.1088/1361-648x/ab884a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/05/2019] [Accepted: 04/09/2020] [Indexed: 06/11/2023]
Abstract
In this paper we formulate a nonlocal density functional theory of inhomogeneous water. We model a water molecule as a couple of oppositely charged sites. The negatively charged sites interact with each other through the Lennard-Jones potential (steric and dispersion interactions), square-well potential (short-range specific interactions due to electron charge transfer), and Coulomb potential, whereas the positively charged sites interact with all types of sites by applying the Coulomb potential only. Taking into account the nonlocal packing effects via the fundamental measure theory, dispersion and specific interactions in the mean-field approximation, and electrostatic interactions at the many-body level through the random phase approximation, we describe the liquid-vapour interface. We demonstrate that our model without explicit account of the association of water molecules due to hydrogen bonding and with explicit account of the electrostatic interactions at the many-body level is able to describe the liquid-vapour coexistence curve and the surface tension at the ambient pressures and temperatures. We obtain very good agreement with available in the literature MD simulation results for density profile of liquid-vapour interface at ambient state parameters. The formulated theory can be used as a theoretical background for describing of the capillary phenomena, occurring in micro- and mesoporous materials.
Collapse
Affiliation(s)
- Yu A Budkov
- School of Applied Mathematics, Tikhonov Moscow Institute of Electronics and Mathematics, National Research University Higher School of Economics, Tallinskaya st. 34, 123458 Moscow, Russia
- G A Krestov Institute of Solution Chemistry of the Russian Academy of Sciences, Academicheskaya st., 1, 153045 Ivanovo, Russia
| | - A L Kolesnikov
- Institut für Nichtklassische Chemie e.V., Permoserstr. 15, 04318 Leipzig, Germany
| |
Collapse
|
12
|
Abstract
Coulomb interactions play a major role in determining the thermodynamics, structure, and dynamics of condensed-phase systems, but often present significant challenges. Computer simulations usually use periodic boundary conditions to minimize corrections from finite cell boundaries but the long range of the Coulomb interactions generates significant contributions from distant periodic images of the simulation cell, usually calculated by Ewald sum techniques. This can add significant overhead to computer simulations and hampers the development of intuitive local pictures and simple analytic theory. In this paper, we present a general framework based on local molecular field theory to accurately determine the contributions from long-ranged Coulomb interactions to the potential of mean force between ionic or apolar hydrophobic solutes in dilute aqueous solutions described by standard classical point charge water models. The simplest approximation leads to a short solvent (SS) model, with truncated solvent-solvent and solute-solvent Coulomb interactions and long-ranged but screened Coulomb interactions only between charged solutes. The SS model accurately describes the interplay between strong short-ranged solute core interactions, local hydrogen-bond configurations, and long-ranged dielectric screening of distant charges, competing effects that are difficult to capture in standard implicit solvent models.
Collapse
|
13
|
Pan C, Yi S, Hu Z. Analytic theory of finite-size effects in supercell modeling of charged interfaces. Phys Chem Chem Phys 2019; 21:14858-14864. [PMID: 31232403 DOI: 10.1039/c9cp02518a] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The Ewald3D sum with the tinfoil boundary condition (e3dtf) evaluates the electrostatic energy of a finite unit cell inside an infinitely periodic supercell. Although it has been used as a de facto standard treatment of electrostatics for simulations of extended polar or charged interfaces, the finite-size effect on simulated properties has yet to be fully understood. There is, however, an intuitive way to quantify the average effect arising from the difference between the e3dtf and Coulomb potentials on the response of mobile charges to contact surfaces with fixed charges and/or to an applied external electric field. Although any charged interface formed by mobile countercharges that compensate the fixed charges fluctuates upon a change in the acting electric field, the distance between a pair of oppositely charged interfaces is found to be nearly stationary, which allows an analytic finite-size correction to the amount of countercharges. Application of the theory to solvated electric double layers (insulator/electrolyte interfaces) predicts that the state of complete charge compensation is invariant with respect to solvent permittivities, which is confirmed by a proper analysis of simulation data in the literature.
Collapse
Affiliation(s)
- Cong Pan
- State Key Laboratory of Supramolecular Structure and Materials, Jilin University, Changchun, 130012, P. R. China.
| | - Shasha Yi
- State Key Laboratory of Supramolecular Structure and Materials, Jilin University, Changchun, 130012, P. R. China.
| | - Zhonghan Hu
- State Key Laboratory of Supramolecular Structure and Materials, Jilin University, Changchun, 130012, P. R. China.
| |
Collapse
|
14
|
Theiss M, Gross J. A third and fourth order perturbation theory for dipolar hard spheres. J Chem Phys 2018; 149:044901. [DOI: 10.1063/1.5035421] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Marc Theiss
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| |
Collapse
|
15
|
Remsing RC, Weeks JD. Alchemical free energy calculations and umbrella sampling with local molecular field theory. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2018. [DOI: 10.1142/s0219633618400035] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Understanding the thermodynamic driving forces underlying any chemical process requires a description of the underlying free energy surface. However, computation of free energies is difficult, often requiring advanced sampling techniques. Moreover, these computations can be further complicated by the evaluation of any long-ranged interactions in the system of interest, such as Coulomb interactions in charged and polar media. Local molecular field theory is a promising approach to avoid many of the conceptual and computational difficulties associated with long-ranged interactions. We present frameworks for performing alchemical free energy calculations and non-Boltzmann sampling with local molecular field theory. We demonstrate that local molecular field theory can be used to perform these free energy calculations with accuracy comparable to traditional methodologies while eliminating the need for explicit treatment of long-ranged interactions in simulations.
Collapse
Affiliation(s)
- Richard C. Remsing
- Institute for Computational Molecular Science, Department of Chemistry, Temple University, Philadelphia, PA 19122, USA
| | - John D. Weeks
- Institute for Physical Science and Technology, Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742, USA
| |
Collapse
|
16
|
Chen HC, Liu YC. Creating functional water by treating excited gold nanoparticles for the applications of green chemistry, energy and medicine: A review. J IND ENG CHEM 2018. [DOI: 10.1016/j.jiec.2017.09.026] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
|
17
|
Gyawali G, Sternfield S, Kumar R, Rick SW. Coarse-Grained Models of Aqueous and Pure Liquid Alkanes. J Chem Theory Comput 2017. [DOI: 10.1021/acs.jctc.7b00389] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Gaurav Gyawali
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| | - Samuel Sternfield
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| | - Revati Kumar
- Department
of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70808, United States
| | - Steven W. Rick
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| |
Collapse
|
18
|
Wu X, Pickard FC, Brooks BR. Isotropic periodic sum for multipole interactions and a vector relation for calculation of the Cartesian multipole tensor. J Chem Phys 2017; 145:164110. [PMID: 27802614 DOI: 10.1063/1.4966019] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
Isotropic periodic sum (IPS) is a method to calculate long-range interactions based on the homogeneity of simulation systems. By using the isotropic periodic images of a local region to represent remote structures, long-range interactions become a function of the local conformation. This function is called the IPS potential; it folds long-ranged interactions into a short-ranged potential and can be calculated as efficiently as a cutoff method. It has been demonstrated that the IPS method produces consistent simulation results, including free energies, as the particle mesh Ewald (PME) method. By introducing the multipole homogeneous background approximation, this work derives multipole IPS potentials, abbreviated as IPSMm, with m being the maximum order of multipole interactions. To efficiently calculate the multipole interactions in Cartesian space, we propose a vector relation that calculates a multipole tensor as a dot product of a radial potential vector and a directional vector. Using model systems with charges, dipoles, and/or quadrupoles, with and without polarizability, we demonstrate that multipole interactions of order m can be described accurately with the multipole IPS potential of order 2 or m - 1, whichever is higher. Through simulations with the multipole IPS potentials, we examined energetic, structural, and dynamic properties of the model systems and demonstrated that the multipole IPS potentials produce very similar results as PME with a local region radius (cutoff distance) as small as 6 Å.
Collapse
Affiliation(s)
- Xiongwu Wu
- Laboratory of Computational Biology, NHLBI, NIH, Bethesda, Maryland 20892, USA
| | - Frank C Pickard
- Laboratory of Computational Biology, NHLBI, NIH, Bethesda, Maryland 20892, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, NHLBI, NIH, Bethesda, Maryland 20892, USA
| |
Collapse
|
19
|
Pan C, Yi S, Hu Z. The effect of electrostatic boundaries in molecular simulations: symmetry matters. Phys Chem Chem Phys 2017; 19:4861-4876. [PMID: 28134360 DOI: 10.1039/c6cp07406e] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Artifacts arise when the long-ranged electrostatic interaction is inappropriately treated in molecular simulations of electrolytes. When the usual Ewald3D sum method with the tinfoil boundary condition (e3dtf) is used for simulations of an interfacial liquid under an external electric field, a straightforward analysis of the liquid structure often suggests unphysical dielectric properties as a consequence of the inaccurate treatment of the electrostatics. In order to understand the underlying mechanism that leads to this apparent violation of thermodynamics, we now derive a new equation in the weak-field limit that, in a mean field view, accounts for the average effect arising from the difference between e3dtf and the sophisticated Ewald2D sum method (e2d). Numerical simulations of a water system in slab geometry confirm the validity of the weak-field limit equation for a series of parameter setup associated with e3dtf. Moreover, a similar procedure applied to a spherically confined water system suggests that corrections to the seemingly inappropriate treatment of the electrostatics in fact vanish. This cancellation of the boundary effect due to symmetry immediately sheds light on the long-lasting problem of the validity of the ad hoc application of e3dtf for bulk systems. In total, we argue that artifacts arising from e3dtf are often predictable and analytical corrections to the straightforward analysis might be applied to reveal consistent thermodynamic properties in liquid simulations.
Collapse
Affiliation(s)
- Cong Pan
- State Key Laboratory of Supramolecular Structure and Materials, College of Chemistry, Jilin University, Changchun, 130012, P. R. China.
| | - Shasha Yi
- State Key Laboratory of Supramolecular Structure and Materials, College of Chemistry, Jilin University, Changchun, 130012, P. R. China.
| | - Zhonghan Hu
- State Key Laboratory of Supramolecular Structure and Materials, College of Chemistry, Jilin University, Changchun, 130012, P. R. China.
| |
Collapse
|
20
|
Wirnsberger P, Fijan D, Šarić A, Neumann M, Dellago C, Frenkel D. Non-equilibrium simulations of thermally induced electric fields in water. J Chem Phys 2016; 144:224102. [DOI: 10.1063/1.4953036] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- P. Wirnsberger
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - D. Fijan
- Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - A. Šarić
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
- Department of Physics and Astronomy, Institute for the Physics of Living Systems, University College London, London WC1E 6BT, United Kingdom
| | - M. Neumann
- Faculty of Physics, University of Vienna, 1090 Vienna, Austria
| | - C. Dellago
- Faculty of Physics, University of Vienna, 1090 Vienna, Austria
| | - D. Frenkel
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
21
|
Bellissent-Funel MC, Hassanali A, Havenith M, Henchman R, Pohl P, Sterpone F, van der Spoel D, Xu Y, Garcia AE. Water Determines the Structure and Dynamics of Proteins. Chem Rev 2016; 116:7673-97. [PMID: 27186992 DOI: 10.1021/acs.chemrev.5b00664] [Citation(s) in RCA: 540] [Impact Index Per Article: 67.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.
Collapse
Affiliation(s)
| | - Ali Hassanali
- International Center for Theoretical Physics, Condensed Matter and Statistical Physics 34151 Trieste, Italy
| | - Martina Havenith
- Ruhr-Universität Bochum , Faculty of Chemistry and Biochemistry Universitätsstraße 150 Building NC 7/72, D-44780 Bochum, Germany
| | - Richard Henchman
- Manchester Institute of Biotechnology The University of Manchester , 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - Peter Pohl
- Johannes Kepler University , Gruberstrasse, 40 4020 Linz, Austria
| | - Fabio Sterpone
- Institut de Biologie Physico-Chimique Laboratoire de Biochimie Théorique 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - David van der Spoel
- Department of Cell and Molecular Biology, Computational and Systems Biology, Uppsala University , 751 24 Uppsala, Sweden
| | - Yao Xu
- Ruhr-Universität Bochum , Faculty of Chemistry and Biochemistry Universitätsstraße 150 Building NC 7/72, D-44780 Bochum, Germany
| | - Angel E Garcia
- Center for Non Linear Studies, Los Alamos National Laboratory , Los Alamos, New Mexico 87545, United States
| |
Collapse
|
22
|
Remsing RC, Liu S, Weeks JD. Long-ranged contributions to solvation free energies from theory and short-ranged models. Proc Natl Acad Sci U S A 2016; 113:2819-26. [PMID: 26929375 PMCID: PMC4801310 DOI: 10.1073/pnas.1521570113] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Long-standing problems associated with long-ranged electrostatic interactions have plagued theory and simulation alike. Traditional lattice sum (Ewald-like) treatments of Coulomb interactions add significant overhead to computer simulations and can produce artifacts from spurious interactions between simulation cell images. These subtle issues become particularly apparent when estimating thermodynamic quantities, such as free energies of solvation in charged and polar systems, to which long-ranged Coulomb interactions typically make a large contribution. In this paper, we develop a framework for determining very accurate solvation free energies of systems with long-ranged interactions from models that interact with purely short-ranged potentials. Our approach is generally applicable and can be combined with existing computational and theoretical techniques for estimating solvation thermodynamics. We demonstrate the utility of our approach by examining the hydration thermodynamics of hydrophobic and ionic solutes and the solvation of a large, highly charged colloid that exhibits overcharging, a complex nonlinear electrostatic phenomenon whereby counterions from the solvent effectively overscreen and locally invert the integrated charge of the solvated object.
Collapse
Affiliation(s)
- Richard C Remsing
- Institute for Physical Science and Technology and Chemical Physics Program, University of Maryland, College Park, MD 20742; Institute for Computational Molecular Science, Temple University, Philadelphia, PA 19122
| | - Shule Liu
- Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742; Department of Chemistry, James Franck Institute and Computation Institute, University of Chicago, Chicago, IL 60637
| | - John D Weeks
- Institute for Physical Science and Technology and Chemical Physics Program, University of Maryland, College Park, MD 20742; Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742;
| |
Collapse
|
23
|
Arslanargin A, Powers A, Beck TL, Rick SW. Models of Ion Solvation Thermodynamics in Ethylene Carbonate and Propylene Carbonate. J Phys Chem B 2015; 120:1497-508. [DOI: 10.1021/acs.jpcb.5b06891] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ayse Arslanargin
- Department
of Physics, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - August Powers
- Department
of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Thomas L. Beck
- Department
of Physics, University of Cincinnati, Cincinnati, Ohio 45221, United States
- Department
of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Steven W. Rick
- Department
of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| |
Collapse
|
24
|
Pan C, Hu Z. Optimized Ewald sum for electrostatics in molecular self-assembly systems at interfaces. Sci China Chem 2015. [DOI: 10.1007/s11426-014-5303-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
25
|
Remsing RC, Patel AJ. Water density fluctuations relevant to hydrophobic hydration are unaltered by attractions. J Chem Phys 2015; 142:024502. [DOI: 10.1063/1.4905009] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Richard C. Remsing
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Amish J. Patel
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
26
|
Remsing RC, Weeks JD. Hydrophobicity Scaling of Aqueous Interfaces by an Electrostatic Mapping. J Phys Chem B 2014; 119:9268-77. [DOI: 10.1021/jp509903n] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Richard C. Remsing
- Institute
for Physical Science and Technology, Department of Chemistry and Biochemistry,
and Chemical Physics Program, University of Maryland, College Park, Maryland 20742, United States
- Department
of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - John D. Weeks
- Institute
for Physical Science and Technology, Department of Chemistry and Biochemistry,
and Chemical Physics Program, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
27
|
Hu Z. Infinite Boundary Terms of Ewald Sums and Pairwise Interactions for Electrostatics in Bulk and at Interfaces. J Chem Theory Comput 2014; 10:5254-64. [DOI: 10.1021/ct500704m] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Zhonghan Hu
- State Key Laboratory
of Supramolecular Structure and Materials and Institute of Theoretical Chemistry, Jilin University, Changchun 130012, P. R. China
| |
Collapse
|
28
|
Qian Z, Wei G. Electric-Field-Induced Phase Transition of Confined Water Nanofilms between Two Graphene Sheets. J Phys Chem A 2014; 118:8922-8. [DOI: 10.1021/jp500989t] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Zhenyu Qian
- State Key
Laboratory of Surface
Physics, Key Laboratory for Computational Physical Sciences (MOE),
and Department of Physics, Fudan University, Shanghai 200433, P. R. China
| | - Guanghong Wei
- State Key
Laboratory of Surface
Physics, Key Laboratory for Computational Physical Sciences (MOE),
and Department of Physics, Fudan University, Shanghai 200433, P. R. China
| |
Collapse
|
29
|
Pan C, Hu Z. Rigorous Error Bounds for Ewald Summation of Electrostatics at Planar Interfaces. J Chem Theory Comput 2014; 10:534-42. [DOI: 10.1021/ct400839x] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Cong Pan
- State
Key Laboratory of Supramolecular Structure and Materials, Jilin University, Changchun, 130012, P. R. China
- Institute
of Theoretical Chemistry, Jilin University, Changchun, 130012, P. R. China
| | - Zhonghan Hu
- State
Key Laboratory of Supramolecular Structure and Materials, Jilin University, Changchun, 130012, P. R. China
- Institute
of Theoretical Chemistry, Jilin University, Changchun, 130012, P. R. China
- Kavli
Institute for Theoretical Physics China, CAS, Beijing 100190, P. R. China
| |
Collapse
|
30
|
Hu Z. Symmetry-preserving mean field theory for electrostatics at interfaces. Chem Commun (Camb) 2014; 50:14397-400. [DOI: 10.1039/c4cc06618a] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A novel method is developed for complex nonuniform electrostatics in computer simulations of molecular liquids at interfaces.
Collapse
Affiliation(s)
- Zhonghan Hu
- State Key Laboratory of Supramolecular Structure and Materials
- Institute of Theoretical Chemistry
- Jilin University
- Changchun, P. R. China
| |
Collapse
|
31
|
Remsing RC, Weeks JD. Dissecting Hydrophobic Hydration and Association. J Phys Chem B 2013; 117:15479-91. [DOI: 10.1021/jp4053067] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Richard C. Remsing
- Institute
for Physical Science and Technology and Chemical
Physics Program, University of Maryland, College Park, Maryland 20742, 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
|
32
|
Shi Y, Beck TL. Length scales and interfacial potentials in ion hydration. J Chem Phys 2013; 139:044504. [DOI: 10.1063/1.4814070] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
|
33
|
Beck TL. The influence of water interfacial potentials on ion hydration in bulk water and near interfaces. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2013.01.008] [Citation(s) in RCA: 79] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
34
|
Archer AJ, Evans R. Relationship between local molecular field theory and density functional theory for non-uniform liquids. J Chem Phys 2013; 138:014502. [PMID: 23298050 DOI: 10.1063/1.4771976] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The local molecular field theory (LMF) developed by Weeks and co-workers has proved successful for treating the structure and thermodynamics of a variety of non-uniform liquids. By reformulating LMF in terms of one-body direct correlation functions we recast the theory in the framework of classical density functional theory (DFT). We show that the general LMF equation for the effective reference potential φ(R)(r) follows directly from the standard mean-field DFT treatment of attractive interatomic forces. Using an accurate (fundamental measures) DFT for the non-uniform hard-sphere reference fluid we determine φ(R)(r) for a hard-core Yukawa liquid adsorbed at a planar hard wall. In the approach to bulk liquid-gas coexistence we find the effective potentials exhibit rich structure that can include damped oscillations at large distances from the wall as well as the repulsive hump near the wall required to generate the low density "gas" layer characteristic of complete drying. We argue that it would be difficult to obtain the same level of detail from other (non-DFT based) implementations of LMF. LMF emphasizes the importance of making an intelligent division of the interatomic pair potential of the full system into a reference part and a remainder that can be treated in mean-field approximation. We investigate different divisions for an exactly solvable one-dimensional model where the pair potential has a hard-core plus a linear attractive tail. Results for the structure factor and the equation of state of the uniform fluid show that including a significant portion of the attraction in the reference system can be much more accurate than treating the full attractive tail in mean-field approximation. We discuss further aspects of the relationship between LMF and DFT.
Collapse
Affiliation(s)
- A J Archer
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, United Kingdom.
| | | |
Collapse
|
35
|
Song W, Lin Y, Baumketner A, Deng S, Cai W, Jacobs DJ. Effect of the Reaction Field on Molecular Forces and Torques Revealed by an Image-Charge Solvation Model. COMMUNICATIONS IN COMPUTATIONAL PHYSICS 2013; 13:129-149. [PMID: 23833681 PMCID: PMC3702192 DOI: 10.4208/cicp.290711.180711s] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
We recently developed the Image-Charge Solvation Model (ICSM), which is an explicit/implicit hybrid model to accurately account for long-range electrostatic forces in molecular dynamics simulations [Lin et al., J. Chem. Phys., 131, 154103, 2009]. The ICSM has a productive spherical volume within the simulation cell for which key physical properties of bulk water are reproduced, such as density, radial distribution function, diffusion constants and dielectric properties. Although the reaction field (RF) is essential, it typically accounts for less than 2% of the total electrostatic force on a water molecule. This observation motivates investigating further the role of the RF within the ICSM. In this report we focus on distributions of forces and torques on water molecules as a function of distance from the origin and make extensive tests over a range of model parameters where Coulomb forces are decomposed into direct interactions from waters modeled explicitly and the RF. Molecular torques due to the RF typically account for 20% of the total torque, revealing why the RF plays an important role in the dielectric properties of simulated water. Moreover, it becomes clear that the buffer layer in the ICSM is essential to mitigate artifacts caused by the discontinuous change in dielectric constants at the explicit/implicit interface.
Collapse
Affiliation(s)
- Wei Song
- Departments of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC 28262, USA
| | - Yuchun Lin
- Departments of Physics and Optical Science, University of North Carolina at Charlotte, Charlotte, NC 28262, USA
- Departments of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28262, USA
| | - Andrij Baumketner
- Departments of Physics and Optical Science, University of North Carolina at Charlotte, Charlotte, NC 28262, USA
- On leave from Institute for Condensed Matter Physics, 1 Svientsitsky Str., Lviv 79011, Ukraine
| | - Shaozhong Deng
- Departments of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28262, USA
| | - Wei Cai
- Departments of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28262, USA
| | - Donald J. Jacobs
- Departments of Physics and Optical Science, University of North Carolina at Charlotte, Charlotte, NC 28262, USA
| |
Collapse
|
36
|
Nguyen TD, Carrillo JMY, Dobrynin AV, Brown WM. A Case Study of Truncated Electrostatics for Simulation of Polyelectrolyte Brushes on GPU Accelerators. J Chem Theory Comput 2012; 9:73-83. [DOI: 10.1021/ct300718x] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Trung Dac Nguyen
- National Center
for Computational
Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831,
United States
| | - Jan-Michael Y. Carrillo
- Polymer Program,
Institute of
Materials Science and Department of Physics, University of Connecticut,
Storrs, Connecticut 06269, United States
| | - Andrey V. Dobrynin
- Polymer Program,
Institute of
Materials Science and Department of Physics, University of Connecticut,
Storrs, Connecticut 06269, United States
| | - W. Michael Brown
- National Center
for Computational
Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831,
United States
| |
Collapse
|
37
|
Hansen JS, Schrøder TB, Dyre JC. Simplistic Coulomb Forces in Molecular Dynamics: Comparing the Wolf and Shifted-Force Approximations. J Phys Chem B 2012; 116:5738-43. [DOI: 10.1021/jp300750g] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- J. S. Hansen
- DNRF Centre “Glass and Time”, IMFUFA,
Department of Science, Systems and Models, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| | - Thomas B. Schrøder
- DNRF Centre “Glass and Time”, IMFUFA,
Department of Science, Systems and Models, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| | - Jeppe C. Dyre
- DNRF Centre “Glass and Time”, IMFUFA,
Department of Science, Systems and Models, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
38
|
Jamadagni SN, Godawat R, Garde S. Hydrophobicity of proteins and interfaces: insights from density fluctuations. Annu Rev Chem Biomol Eng 2012; 2:147-71. [PMID: 22432614 DOI: 10.1146/annurev-chembioeng-061010-114156] [Citation(s) in RCA: 170] [Impact Index Per Article: 14.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Macroscopic characterizations of hydrophobicity (e.g., contact angle measurements) do not extend to the surfaces of proteins and nanoparticles. Molecular measures of hydrophobicity of such surfaces need to account for the behavior of hydration water. Theory and state-of-the-art simulations suggest that water density fluctuations provide such a measure; fluctuations are enhanced near hydrophobic surfaces and quenched with increasing surface hydrophilicity. Fluctuations affect conformational equilibria and dynamics of molecules at interfaces. Enhanced fluctuations are reflected in enhanced cavity formation, more favorable binding of hydrophobic solutes, increased compressibility of hydration water, and enhanced water-water correlations at hydrophobic surfaces. These density fluctuation-based measures can be used to develop practical methods to map the hydrophobicity/philicity of heterogeneous surfaces including those of proteins. They highlight that the hydrophobicity of a group is context dependent and is significantly affected by its environment (e.g., chemistry and topography) and especially by confinement. The ability to include information about hydration water in mapping hydrophobicity is expected to significantly impact our understanding of protein-protein interactions as well as improve drug design and discovery methods and bioseparation processes.
Collapse
Affiliation(s)
- Sumanth N Jamadagni
- The Howard P. Isermann Department of Chemical and Biological Engineering, Rensselaer Polytechnic Institute, Troy, NY 12180, USA.
| | | | | |
Collapse
|
39
|
Rodgers JM, Sørensen J, de Meyer FJM, Schiøtt B, Smit B. Understanding the Phase Behavior of Coarse-Grained Model Lipid Bilayers through Computational Calorimetry. J Phys Chem B 2012; 116:1551-69. [DOI: 10.1021/jp207837v] [Citation(s) in RCA: 63] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jocelyn M. Rodgers
- Physical Biosciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Jesper Sørensen
- Interdisciplinary Nanoscience Center (iNANO), Aarhus University, Ny Munkegade 118, 8000 Aarhus C, Denmark
- Center for Insoluble Protein Structures (inSPIN), Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Frédérick J.-M. de Meyer
- Department of Chemical Engineering, University of California, Berkeley, 101B Gilman Hall, Berkeley, California 94720-1462, United States
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Birgit Schiøtt
- Interdisciplinary Nanoscience Center (iNANO), Aarhus University, Ny Munkegade 118, 8000 Aarhus C, Denmark
- Center for Insoluble Protein Structures (inSPIN), Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Berend Smit
- Department of Chemical Engineering, University of California, Berkeley, 101B Gilman Hall, Berkeley, California 94720-1462, United States
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, 101B Gilman Hall, Berkeley, California 94720-1462, United States
| |
Collapse
|
40
|
Hansen JS, Dyre JC, Daivis PJ, Todd BD, Bruus H. Nanoflow hydrodynamics. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:036311. [PMID: 22060496 DOI: 10.1103/physreve.84.036311] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/28/2011] [Revised: 08/04/2011] [Indexed: 05/31/2023]
Abstract
We show by nonequilibrium molecular dynamics simulations that the Navier-Stokes equation does not correctly describe water flow in a nanoscale geometry. It is argued that this failure reflects the fact that the coupling between the intrinsic rotational and translational degrees of freedom becomes important for nanoflows. The coupling is correctly accounted for by the extended Navier-Stokes equations that include the intrinsic angular momentum as an independent hydrodynamic degree of freedom.
Collapse
Affiliation(s)
- J S Hansen
- Danish National Research Foundation (DNRF) Centre Glass and Time, IMFUFA, Department of Science, Systems and Models, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark.
| | | | | | | | | |
Collapse
|
41
|
Abstract
Monovalent ion hydration entropies are analyzed via energetic partitioning of the potential distribution theorem free energy. Extensive molecular dynamics simulations and free energy calculations are performed over a range of temperatures to determine the electrostatic and van der Waals components of the entropy. The far-field electrostatic contribution is negative and small in magnitude, and it does not vary significantly as a function of ion size, consistent with dielectric models. The local electrostatic contribution, however, varies widely as a function of ion size; the sign yields a direct indication of the kosmotropic (strongly hydrated) or chaotropic (weakly hydrated) nature of the ion hydration. The results provide a thermodynamic signature for specific ion effects in hydration and are consistent with experiments that suggest minimal perturbations of water structure outside the first hydration shell. The hydration entropies are also examined in relation to the corresponding entropies for the isoelectronic rare gas pairs; an inverse correlation is observed, as expected from thermodynamic hydration data.
Collapse
Affiliation(s)
- Thomas L Beck
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221-0172, United States.
| |
Collapse
|
42
|
Rodgers JM, Hu Z, Weeks JD. On the efficient and accurate short-ranged simulations of uniform polar molecular liquids. Mol Phys 2011. [DOI: 10.1080/00268976.2011.554332] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
43
|
Hu Z, Weeks JD. Efficient solutions of self-consistent mean field equations for dewetting and electrostatics in nonuniform liquids. PHYSICAL REVIEW LETTERS 2010; 105:140602. [PMID: 21230822 DOI: 10.1103/physrevlett.105.140602] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/28/2010] [Revised: 08/09/2010] [Indexed: 05/30/2023]
Abstract
We use a new configuration-based version of linear response theory to efficiently solve self-consistent mean field equations relating an effective single particle potential to the induced density. The versatility and accuracy of the method is illustrated by applications to dewetting of a hard sphere solute in a Lennard-Jones fluid, the interplay between local hydrogen bond structure and electrostatics for water confined between two hydrophobic walls, and to ion pairing in ionic solutions. Simulation time has been reduced by more than an order of magnitude over previous methods.
Collapse
Affiliation(s)
- Zhonghan Hu
- Institute for Physical Science and Technology and Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, USA
| | | |
Collapse
|
44
|
LeBard DN, Matyushov DV. Ferroelectric Hydration Shells around Proteins: Electrostatics of the Protein−Water Interface. J Phys Chem B 2010; 114:9246-58. [DOI: 10.1021/jp1006999] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- David N. LeBard
- Center for Biological Physics, Arizona State University, PO Box 871604, Tempe, Arizona 85287-1604
| | - Dmitry V. Matyushov
- Center for Biological Physics, Arizona State University, PO Box 871604, Tempe, Arizona 85287-1604
| |
Collapse
|
45
|
Rodgers JM, Weeks JD. Accurate thermodynamics for short-ranged truncations of Coulomb interactions in site-site molecular models. J Chem Phys 2009; 131:244108. [DOI: 10.1063/1.3276729] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
46
|
Mella M, Harris KDM. Pathways for hydrogen bond switching in a tetrameric methanol cluster. Phys Chem Chem Phys 2009; 11:11340-6. [PMID: 20024403 DOI: 10.1039/b911556k] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Computational techniques (second order Møller-Plesset MP2 perturbation theory in conjunction with medium and large size basis sets) are applied to explore structural aspects of a hydrogen-bonded tetrameric cluster of methanol molecules, based geometrically on a tetrahedral arrangement of the four oxygen atoms of the cluster. The hydrogen-bonded structures that represent minima on the potential energy surface are established, and the complete set of pathways that allow interconversion between these structures through "switching" of the hydrogen bonding arrangement are elucidated. The implications of these results in terms of dynamic properties of the cluster are discussed.
Collapse
Affiliation(s)
- Massimo Mella
- School of Chemistry, Cardiff University, Park Place, Cardiff, Wales, UK CF10 3AT.
| | | |
Collapse
|
47
|
Raghunathan AV, Aluru NR. An empirical potential based quasicontinuum theory for structural prediction of water. J Chem Phys 2009; 131:184703. [DOI: 10.1063/1.3258280] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
48
|
Wu X, Brooks BR. Isotropic periodic sum of electrostatic interactions for polar systems. J Chem Phys 2009; 131:024107. [PMID: 19603970 PMCID: PMC2723922 DOI: 10.1063/1.3160730] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2009] [Accepted: 06/09/2009] [Indexed: 11/14/2022] Open
Abstract
Isotropic periodic sum (IPS) is a method to calculate long-range interactions based on homogeneity of simulation systems. Long-range interactions are represented by interactions with isotropic periodic images of a defined local region and can be reduced to short ranged IPS potentials. The original electrostatic three-dimensional (3D)-IPS potential was derived based on a nonpolar homogeneous approximation and its application is limited to nonpolar or weak polar systems. This work derived a polar electrostatic 3D-IPS potential based on polar interactions. For the convenience of application, polynomial functions with rationalized coefficients are proposed for electrostatic and Lennard-Jones 3D-IPS potentials. Model systems of various polarities and several commonly used solvent systems are simulated to evaluate the 3D-IPS potentials. It is demonstrated that for polar systems the polar electrostatic 3D-IPS potential has much improved accuracy as compared to the nonpolar 3D-IPS potential. For homogeneous systems, the polar electrostatic 3D-IPS potential with a local region radius or cutoff distance of as small as 10 A can satisfactorily reproduce energetic, structural, and dynamic properties from the particle-meshed-Ewald method. For both homogeneous and heterogeneous systems, the 3D-IPS/discrete fast Fourier transform method using either the nonpolar or the polar electrostatic 3D-IPS potentials results in very similar simulation results.
Collapse
Affiliation(s)
- Xiongwu Wu
- Laboratory of Computational Biology, NHLBI, NIH, Bethesda, Maryland 20892, USA.
| | | |
Collapse
|
49
|
González-Mozuelos P, de la Cruz MO. Asymmetric charge renormalization for nanoparticles in aqueous media. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 79:031901. [PMID: 19391965 DOI: 10.1103/physreve.79.031901] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/27/2008] [Revised: 01/19/2009] [Indexed: 05/27/2023]
Abstract
The effective renormalized charge of nanoparticles in an aqueous electrolyte is essential to determine their solubility. By using a molecular model for the supporting aqueous electrolyte, we find that the effective renormalized charge of the nanoparticles is strongly dependent on the sign of the bare charge. Negatively charged nanoparticles have a lower effective renormalized charge than positively charged nanoparticles. The degree of asymmetry is a nonmonotonic function of the bare charge of the nanoparticle. We show that the effect is due to the asymmetric charge distribution of the water molecules, which we model using a simple three-site molecular structure of point charges.
Collapse
Affiliation(s)
- P González-Mozuelos
- Department of Chemistry and Department of Materials Science and Engineering, Northwestern University, Evanston, Illinois 60208, USA
| | | |
Collapse
|
50
|
Abstract
Hydrophobicity manifests itself differently on large and small length scales. This review focuses on large-length-scale hydrophobicity, particularly on dewetting at single hydrophobic surfaces and drying in regions bounded on two or more sides by hydrophobic surfaces. We review applicable theories, simulations, and experiments pertaining to large-scale hydrophobicity in physical and biomolecular systems and clarify some of the critical issues pertaining to this subject. Given space constraints, we cannot review all the significant and interesting work in this active field.
Collapse
Affiliation(s)
- Bruce J Berne
- Department of Chemistry, Columbia University, New York, New York 10027, USA.
| | | | | |
Collapse
|