1
|
Evangelakis A, Panahian Jand S, Delle Site L. Path Integral Molecular Dynamics of Liquid Water in a Mean-Field Particle Reservoir. Chemistry 2022; 11:e202100286. [PMID: 35218165 PMCID: PMC8973269 DOI: 10.1002/open.202100286] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Revised: 02/04/2022] [Indexed: 11/17/2022]
Abstract
We present a simulation scheme for path integral simulation of molecular liquids where a small open region is embedded in a large reservoir of non interacting point‐particles. The scheme is based on the latest development of the adaptive resolution technique AdResS and allows for the space‐dependent change of molecular resolution from a path integral representation with 120 degrees of freedom to a point particle that does not interact with other molecules and vice versa. The method is applied to liquid water and implies a sizable gain regarding the request of computational resources compared to full path integral simulations. Given the role of water as universal solvent with a specific hydrogen bonding network, the path integral treatment of water molecules is important to describe the quantum effects of hydrogen atoms’ delocalization in space on the hydrogen bonding network. The method presented here implies feasible computational efforts compared to full path integral simulations of liquid water which, on large scales, are often prohibitive.
Collapse
Affiliation(s)
- Antonios Evangelakis
- Freie Universität Berlin, Institute of Mathematics, Arnimallee 6, 14195, Berlin, Germany
| | - Sara Panahian Jand
- Freie Universität Berlin, Institute of Mathematics, Arnimallee 6, 14195, Berlin, Germany
| | - Luigi Delle Site
- Freie Universität Berlin, Institute of Mathematics, Arnimallee 6, 14195, Berlin, Germany
| |
Collapse
|
2
|
Sun Z, Liu Z. BAR‐Based Multi‐Dimensional Nonequilibrium Pulling for Indirect Construction of QM/MM Free Energy Landscapes: Varying the QM Region. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202100185] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| | - Zhirong Liu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| |
Collapse
|
3
|
Cortes-Huerto R, Praprotnik M, Kremer K, Delle Site L. From adaptive resolution to molecular dynamics of open systems. THE EUROPEAN PHYSICAL JOURNAL. B 2021; 94:189. [PMID: 34720711 PMCID: PMC8547219 DOI: 10.1140/epjb/s10051-021-00193-w] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/21/2021] [Accepted: 08/31/2021] [Indexed: 05/14/2023]
Abstract
ABSTRACT We provide an overview of the Adaptive Resolution Simulation method (AdResS) based on discussing its basic principles and presenting its current numerical and theoretical developments. Examples of applications to systems of interest to soft matter, chemical physics, and condensed matter illustrate the method's advantages and limitations in its practical use and thus settle the challenge for further future numerical and theoretical developments.
Collapse
Affiliation(s)
| | - Matej Praprotnik
- Laboratory for Molecular Modeling, National Institute of Chemistry, Ljubljana, Slovenia and Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana, Slovenia
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Luigi Delle Site
- Department of Mathematics and Computer Science, Institute for Mathematics, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|
4
|
Cendagorta JR, Shen H, Bačić Z, Tuckerman ME. Enhanced Sampling Path Integral Methods Using Neural Network Potential Energy Surfaces with Application to Diffusion in Hydrogen Hydrates. ADVANCED THEORY AND SIMULATIONS 2020. [DOI: 10.1002/adts.202000258] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
| | - Hengyuan Shen
- Department of Chemistry New York University Shanghai 1555 Century Avenue Pudong Shanghai 200122 China
| | - Zlatko Bačić
- Department of Chemistry New York University New York NY 10003 USA
- NYU‐ECNU Center for Computational Chemistry at NYU Shanghai 3663 Zhongshan Road, North Shanghai 200062 China
| | - Mark E. Tuckerman
- Department of Chemistry New York University New York NY 10003 USA
- NYU‐ECNU Center for Computational Chemistry at NYU Shanghai 3663 Zhongshan Road, North Shanghai 200062 China
- Courant Institute of Mathematical Sciences New York University New York NY 10012 USA
| |
Collapse
|
5
|
Vergadou N, Theodorou DN. Molecular Modeling Investigations of Sorption and Diffusion of Small Molecules in Glassy Polymers. MEMBRANES 2019; 9:E98. [PMID: 31398889 PMCID: PMC6723301 DOI: 10.3390/membranes9080098] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/12/2019] [Revised: 07/22/2019] [Accepted: 07/23/2019] [Indexed: 11/16/2022]
Abstract
With a wide range of applications, from energy and environmental engineering, such as in gas separations and water purification, to biomedical engineering and packaging, glassy polymeric materials remain in the core of novel membrane and state-of the art barrier technologies. This review focuses on molecular simulation methodologies implemented for the study of sorption and diffusion of small molecules in dense glassy polymeric systems. Basic concepts are introduced and systematic methods for the generation of realistic polymer configurations are briefly presented. Challenges related to the long length and time scale phenomena that govern the permeation process in the glassy polymer matrix are described and molecular simulation approaches developed to address the multiscale problem at hand are discussed.
Collapse
Affiliation(s)
- Niki Vergadou
- Molecular Thermodynamics and Modelling of Materials Laboratory, Institute of Nanoscience and Nanotechnology, National Center for Scientific Research Demokritos, Aghia Paraskevi Attikis, GR-15310 Athens, Greece.
| | - Doros N Theodorou
- School of Chemical Engineering, National Technical University of Athens, GR 15780 Athens, Greece
| |
Collapse
|
6
|
Ciccotti G, Delle Site L. The physics of open systems for the simulation of complex molecular environments in soft matter. SOFT MATTER 2019; 15:2114-2124. [PMID: 30761396 DOI: 10.1039/c8sm02523a] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Molecular dynamics (MD) has become one of the most powerful tools of investigation in soft matter. Despite such success, simulations of large molecular environments are mostly run using the approximation of closed systems without the possibility of exchange of matter. Due to the molecular complexity of soft matter systems, an optimal simulation strategy would require the application of concurrent multiscale resolution approaches such that each part of a large system can be considered as an open subsystem at a high resolution embedded in a large coarser reservoir of energy and particles. This paper discusses the current capability and the future perspectives of multiscale adaptive resolution MD methods to satisfy the conceptual principles of open systems and to perform simulations of complex molecular environments in soft matter.
Collapse
Affiliation(s)
- Giovanni Ciccotti
- Instituto per le Applicazioni del Calcolo, Consiglio Nazionale delle Ricerche, and Universita' La Sapienza, Rome, Italy.
| | | |
Collapse
|
7
|
Delle Site L. Simulation of Many-Electron Systems That Exchange Matter with the Environment. ADVANCED THEORY AND SIMULATIONS 2018. [DOI: 10.1002/adts.201800056] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Luigi Delle Site
- Institute for Mathematics; Freie Universität Berlin; D-14195 Berlin Germany
| |
Collapse
|
8
|
Litman Y, Donadio D, Ceriotti M, Rossi M. Decisive role of nuclear quantum effects on surface mediated water dissociation at finite temperature. J Chem Phys 2018; 148:102320. [PMID: 29544260 DOI: 10.1063/1.5002537] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Water molecules adsorbed on inorganic substrates play an important role in several technological applications. In the presence of light atoms in adsorbates, nuclear quantum effects (NQEs) influence the structural stability and the dynamical properties of these systems. In this work, we explore the impact of NQEs on the dissociation of water wires on stepped Pt(221) surfaces. By performing ab initio molecular dynamics simulations with van der Waals corrected density functional theory, we note that several competing minima for both intact and dissociated structures are accessible at finite temperatures, making it important to assess whether harmonic estimates of the quantum free energy are sufficient to determine the relative stability of the different states. We thus perform ab initio path integral molecular dynamics (PIMD) in order to calculate these contributions taking into account the conformational entropy and anharmonicities at finite temperatures. We propose that when adsorption is weak and NQEs on the substrate are negligible, PIMD simulations can be performed through a simple partition of the system, resulting in considerable computational savings. We then calculate the full contribution of NQEs to the free energies, including also anharmonic terms. We find that they result in an increase of up to 20% of the quantum contribution to the dissociation free energy compared with the harmonic estimates. We also find that the dissociation process has a negligible contribution from tunneling but is dominated by zero point energies, which can enhance the rate of dissociation by three orders of magnitude. Finally we highlight how both temperature and NQEs indirectly impact dipoles and the redistribution of electron density, causing work function changes of up to 0.4 eV with respect to static estimates. This quantitative determination of the change in the work function provides a possible approach to determine experimentally the most stable configurations of water oligomers on the stepped surfaces.
Collapse
Affiliation(s)
- Yair Litman
- Fritz Haber Institute of the Max Planck Society, Faradayweg 4-6, 14195 Berlin, Germany
| | - Davide Donadio
- Department of Chemistry, University of California Davis, One Shields Ave., Davis, California 95616, USA
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Mariana Rossi
- Fritz Haber Institute of the Max Planck Society, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
9
|
Kreis K, Kremer K, Potestio R, Tuckerman ME. From classical to quantum and back: Hamiltonian adaptive resolution path integral, ring polymer, and centroid molecular dynamics. J Chem Phys 2018; 147:244104. [PMID: 29289131 DOI: 10.1063/1.5000701] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Path integral-based methodologies play a crucial role for the investigation of nuclear quantum effects by means of computer simulations. However, these techniques are significantly more demanding than corresponding classical simulations. To reduce this numerical effort, we recently proposed a method, based on a rigorous Hamiltonian formulation, which restricts the quantum modeling to a small but relevant spatial region within a larger reservoir where particles are treated classically. In this work, we extend this idea and show how it can be implemented along with state-of-the-art path integral simulation techniques, including path-integral molecular dynamics, which allows for the calculation of quantum statistical properties, and ring-polymer and centroid molecular dynamics, which allow the calculation of approximate quantum dynamical properties. To this end, we derive a new integration algorithm that also makes use of multiple time-stepping. The scheme is validated via adaptive classical-path-integral simulations of liquid water. Potential applications of the proposed multiresolution method are diverse and include efficient quantum simulations of interfaces as well as complex biomolecular systems such as membranes and proteins.
Collapse
Affiliation(s)
- Karsten Kreis
- Max-Planck-Institut für Polymerforschung, Ackermannweg 10, 55128 Mainz, Germany
| | - Kurt Kremer
- Max-Planck-Institut für Polymerforschung, Ackermannweg 10, 55128 Mainz, Germany
| | - Raffaello Potestio
- Max-Planck-Institut für Polymerforschung, Ackermannweg 10, 55128 Mainz, Germany
| | - Mark E Tuckerman
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
| |
Collapse
|
10
|
Wang CH, Duster AW, Aydintug BO, Zarecki MG, Lin H. Chloride Ion Transport by the E. coli CLC Cl -/H + Antiporter: A Combined Quantum-Mechanical and Molecular-Mechanical Study. Front Chem 2018; 6:62. [PMID: 29594103 PMCID: PMC5859129 DOI: 10.3389/fchem.2018.00062] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2018] [Accepted: 02/26/2018] [Indexed: 01/05/2023] Open
Abstract
We performed steered molecular dynamics (SMD) and umbrella sampling simulations of Cl- ion migration through the transmembrane domain of a prototypical E. coli CLC Cl-/H+ antiporter by employing combined quantum-mechanical (QM) and molecular-mechanical (MM) calculations. The SMD simulations revealed interesting conformational changes of the protein. While no large-amplitude motions of the protein were observed during pore opening, the side chain rotation of the protonated external gating residue Glu148 was found to be critical for full access of the channel entrance by Cl-. Moving the anion into the external binding site (Sext) induced small-amplitude shifting of the protein backbone at the N-terminal end of helix F. As Cl- traveled through the pore, rigid-body swinging motions of helix R separated it from helix D. Helix R returned to its original position once Cl- exited the channel. Population analysis based on polarized wavefunction from QM/MM calculations discovered significant (up to 20%) charge loss for Cl- along the ion translocation pathway inside the pore. The delocalized charge was redistributed onto the pore residues, especially the functional groups containing π bonds (e.g., the Tyr445 side chain), while the charges of the H atoms coordinating Cl- changed almost negligibly. Potentials of mean force computed from umbrella sampling at the QM/MM and MM levels both displayed barriers at the same locations near the pore entrance and exit. However, the QM/MM PMF showed higher barriers (~10 kcal/mol) than the MM PMF (~2 kcal/mol). Binding energy calculations indicated that the interactions between Cl- and certain pore residues were overestimated by the semi-empirical PM3 Hamiltonian and underestimated by the CHARMM36 force fields, both of which were employed in the umbrella sampling simulations. In particular, CHARMM36 underestimated binding interactions for the functional groups containing π bonds, missing the stabilizations of the Cl- ion due to electron delocalization. The results suggested that it is important to explore these quantum effects for accurate descriptions of the Cl- transport.
Collapse
Affiliation(s)
- Chun-Hung Wang
- Department of Chemistry, University of Colorado Denver, Denver, CO, United States
| | - Adam W Duster
- Department of Chemistry, University of Colorado Denver, Denver, CO, United States
| | - Baris O Aydintug
- Department of Chemistry, University of Colorado Denver, Denver, CO, United States
| | - MacKenzie G Zarecki
- Department of Chemistry, University of Colorado Denver, Denver, CO, United States
| | - Hai Lin
- Department of Chemistry, University of Colorado Denver, Denver, CO, United States
| |
Collapse
|
11
|
|
12
|
Guzman HV, Junghans C, Kremer K, Stuehn T. Scalable and fast heterogeneous molecular simulation with predictive parallelization schemes. Phys Rev E 2017; 96:053311. [PMID: 29347684 DOI: 10.1103/physreve.96.053311] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2017] [Indexed: 11/07/2022]
Abstract
Multiscale and inhomogeneous molecular systems are challenging topics in the field of molecular simulation. In particular, modeling biological systems in the context of multiscale simulations and exploring material properties are driving a permanent development of new simulation methods and optimization algorithms. In computational terms, those methods require parallelization schemes that make a productive use of computational resources for each simulation and from its genesis. Here, we introduce the heterogeneous domain decomposition approach, which is a combination of an heterogeneity-sensitive spatial domain decomposition with an a priori rearrangement of subdomain walls. Within this approach, the theoretical modeling and scaling laws for the force computation time are proposed and studied as a function of the number of particles and the spatial resolution ratio. We also show the new approach capabilities, by comparing it to both static domain decomposition algorithms and dynamic load-balancing schemes. Specifically, two representative molecular systems have been simulated and compared to the heterogeneous domain decomposition proposed in this work. These two systems comprise an adaptive resolution simulation of a biomolecule solvated in water and a phase-separated binary Lennard-Jones fluid.
Collapse
Affiliation(s)
- Horacio V Guzman
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Christoph Junghans
- Computer, Computational, and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Torsten Stuehn
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|