1
|
Ricardi N, González-Espinoza CE, Adam S, Church JR, Schapiro I, Wesołowski TA. Embedding Nonrigid Solutes in an Averaged Environment: A Case Study on Rhodopsins. J Chem Theory Comput 2023; 19:5289-5302. [PMID: 37441785 PMCID: PMC10413860 DOI: 10.1021/acs.jctc.3c00285] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2023]
Abstract
Many simulation methods concerning solvated molecules are based on the assumption that the solvated species and the solvent can be characterized by some representative structures of the solute and some embedding potential corresponding to this structure. While the averaging of the solvent configurations to obtain an embedding potential has been studied in great detail, this hinges on a single solute structure representation. This assumption is re-examined and generalized for conformationally flexible solutes and tested on 4 nonrigid systems. In this generalized approach, the solute is characterized by a set of representative structures and the corresponding embedding potentials. The representative structures are identified by means of subdividing the statistical ensemble, which in this work is generated by a constant-temperature molecular dynamics simulation. The embedding potential defined in the Frozen-Density Embedding Theory is used to characterize the average effect of the solvent in each subensemble. The numerical examples concern the vertical excitation energies of protonated retinal Schiff bases in protein environments. It is comprehensively shown that subensemble averaging leads to huge computational savings compared with explicit averaging of the excitation energies in the whole ensemble while introducing only minor errors in the case of the systems examined.
Collapse
Affiliation(s)
- Niccolò Ricardi
- Department of Physical Chemistry, University of Geneva, 1205 Geneva, Switzerland
| | | | - Suliman Adam
- Fritz Haber Center for Molecular Dynamics, Hebrew University of Jerusalem Israel, 91904 Jerusalem, Israel
| | - Jonathan R Church
- Fritz Haber Center for Molecular Dynamics, Hebrew University of Jerusalem Israel, 91904 Jerusalem, Israel
| | - Igor Schapiro
- Fritz Haber Center for Molecular Dynamics, Hebrew University of Jerusalem Israel, 91904 Jerusalem, Israel
| | | |
Collapse
|
2
|
Fu M, Wesolowski TA. Excitation Energies of Embedded Chromophores from Frozen-Density Embedding Theory Using State-Specific Electron Densities of the Environment. J Phys Chem A 2023; 127:535-545. [PMID: 36599107 DOI: 10.1021/acs.jpca.2c07747] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Starting from the Perdew-Levy theorem on extrema of the Hohenberg-Kohn functional, the expression for the vertical excitation energy is derived within the formal framework of Frozen-Density Embedding Theory (FDET) that makes it possible to use state-specific electron densities of the environment (ρB) of an embedded species. The derived general expression involves the embedded wave functions for ground and excited states that are orthogonal and is exact up to quadratic terms in the appropriate density expansion. It can be applied in practice using various methods differing in the treatment of the electron-electron correlation for embedded electrons, the method to evaluate different contributions to the excitation energy, the method to generate state-specific ρB, and the approximation used for the non-electrostatic component of the FDET embedding potential. The derived expression is applied for 47 local excitations in 10 embedded organic chromophores. The explicit treatment of the differential polarization of ρB improves indeed the accuracy of the excitation energy as compared to the implicit treatment in which the same ρB is used for all states of embedded chromophore. For 47 local excitations in 10 embedded organic chromophores, the average absolute errors in excitation energies drop from 0.04 to 0.03 eV and their standard deviations from 0.032 to 0.025 eV, respectively. The maximal errors show similar trends.
Collapse
Affiliation(s)
- Mingxue Fu
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 4Genève, Switzerland
| | - Tomasz A Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 4Genève, Switzerland
| |
Collapse
|
3
|
Teale AM, Helgaker T, Savin A, Adamo C, Aradi B, Arbuznikov AV, Ayers PW, Baerends EJ, Barone V, Calaminici P, Cancès E, Carter EA, Chattaraj PK, Chermette H, Ciofini I, Crawford TD, De Proft F, Dobson JF, Draxl C, Frauenheim T, Fromager E, Fuentealba P, Gagliardi L, Galli G, Gao J, Geerlings P, Gidopoulos N, Gill PMW, Gori-Giorgi P, Görling A, Gould T, Grimme S, Gritsenko O, Jensen HJA, Johnson ER, Jones RO, Kaupp M, Köster AM, Kronik L, Krylov AI, Kvaal S, Laestadius A, Levy M, Lewin M, Liu S, Loos PF, Maitra NT, Neese F, Perdew JP, Pernal K, Pernot P, Piecuch P, Rebolini E, Reining L, Romaniello P, Ruzsinszky A, Salahub DR, Scheffler M, Schwerdtfeger P, Staroverov VN, Sun J, Tellgren E, Tozer DJ, Trickey SB, Ullrich CA, Vela A, Vignale G, Wesolowski TA, Xu X, Yang W. DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science. Phys Chem Chem Phys 2022; 24:28700-28781. [PMID: 36269074 PMCID: PMC9728646 DOI: 10.1039/d2cp02827a] [Citation(s) in RCA: 62] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Accepted: 08/09/2022] [Indexed: 12/13/2022]
Abstract
In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.
Collapse
Affiliation(s)
- Andrew M. Teale
- School of Chemistry, University of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Trygve Helgaker
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Andreas Savin
- Laboratoire de Chimie Théorique, CNRS and Sorbonne University, 4 Place Jussieu, CEDEX 05, 75252 Paris, France.
| | - Carlo Adamo
- PSL University, CNRS, ChimieParisTech-PSL, Institute of Chemistry for Health and Life Sciences, i-CLeHS, 11 rue P. et M. Curie, 75005 Paris, France.
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, University of Bremen, P.O. Box 330440, D-28334 Bremen, Germany.
| | - Alexei V. Arbuznikov
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7Straße des 17. Juni 13510623Berlin
| | | | - Evert Jan Baerends
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56125 Pisa, Italy.
| | - Patrizia Calaminici
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav), CDMX, 07360, Mexico.
| | - Eric Cancès
- CERMICS, Ecole des Ponts and Inria Paris, 6 Avenue Blaise Pascal, 77455 Marne-la-Vallée, France.
| | - Emily A. Carter
- Department of Mechanical and Aerospace Engineering and the Andlinger Center for Energy and the Environment, Princeton UniversityPrincetonNJ 08544-5263USA
| | | | - Henry Chermette
- Institut Sciences Analytiques, Université Claude Bernard Lyon1, CNRS UMR 5280, 69622 Villeurbanne, France.
| | - Ilaria Ciofini
- PSL University, CNRS, ChimieParisTech-PSL, Institute of Chemistry for Health and Life Sciences, i-CLeHS, 11 rue P. et M. Curie, 75005 Paris, France.
| | - T. Daniel Crawford
- Department of Chemistry, Virginia TechBlacksburgVA 24061USA,Molecular Sciences Software InstituteBlacksburgVA 24060USA
| | - Frank De Proft
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium.
| | | | - Claudia Draxl
- Institut für Physik and IRIS Adlershof, Humboldt-Universität zu Berlin, 12489 Berlin, Germany. .,Fritz-Haber-Institut der Max-Planck-Gesellschaft, 14195 Berlin, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, University of Bremen, P.O. Box 330440, D-28334 Bremen, Germany. .,Beijing Computational Science Research Center (CSRC), 100193 Beijing, China.,Shenzhen JL Computational Science and Applied Research Institute, 518110 Shenzhen, China
| | - Emmanuel Fromager
- Laboratoire de Chimie Quantique, Institut de Chimie, CNRS/Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France.
| | - Patricio Fuentealba
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago, Chile.
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute, and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, USA.
| | - Giulia Galli
- Pritzker School of Molecular Engineering and Department of Chemistry, The University of Chicago, Chicago, IL, USA.
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China. .,Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA
| | - Paul Geerlings
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium.
| | - Nikitas Gidopoulos
- Department of Physics, Durham University, South Road, Durham DH1 3LE, UK.
| | - Peter M. W. Gill
- School of Chemistry, University of SydneyCamperdown NSW 2006Australia
| | - Paola Gori-Giorgi
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Andreas Görling
- Chair of Theoretical Chemistry, University of Erlangen-Nuremberg, Egerlandstrasse 3, 91058 Erlangen, Germany.
| | - Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University, Gold Coast, Qld 4222, Australia.
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany.
| | - Oleg Gritsenko
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Hans Jørgen Aagaard Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark.
| | - Erin R. Johnson
- Department of Chemistry, Dalhousie UniversityHalifaxNova ScotiaB3H 4R2Canada
| | - Robert O. Jones
- Peter Grünberg Institut PGI-1, Forschungszentrum Jülich52425 JülichGermany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, 10623, Berlin.
| | - Andreas M. Köster
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav)CDMX07360Mexico
| | - Leeor Kronik
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth, 76100, Israel.
| | - Anna I. Krylov
- Department of Chemistry, University of Southern CaliforniaLos AngelesCalifornia 90089USA
| | - Simen Kvaal
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Andre Laestadius
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Mel Levy
- Department of Chemistry, Tulane University, New Orleans, Louisiana, 70118, USA.
| | - Mathieu Lewin
- CNRS & CEREMADE, Université Paris-Dauphine, PSL Research University, Place de Lattre de Tassigny, 75016 Paris, France.
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC 27599-3420, USA. .,Department of Chemistry, University of North Carolina, Chapel Hill, NC 27599-3290, USA
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France.
| | - Neepa T. Maitra
- Department of Physics, Rutgers University at Newark101 Warren StreetNewarkNJ 07102USA
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany.
| | - John P. Perdew
- Departments of Physics and Chemistry, Temple UniversityPhiladelphiaPA 19122USA
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland.
| | - Pascal Pernot
- Institut de Chimie Physique, UMR8000, CNRS and Université Paris-Saclay, Bât. 349, Campus d'Orsay, 91405 Orsay, France.
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA. .,Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA
| | - Elisa Rebolini
- Institut Laue Langevin, 71 avenue des Martyrs, 38000 Grenoble, France.
| | - Lucia Reining
- Laboratoire des Solides Irradiés, CNRS, CEA/DRF/IRAMIS, École Polytechnique, Institut Polytechnique de Paris, F-91120 Palaiseau, France. .,European Theoretical Spectroscopy Facility
| | - Pina Romaniello
- Laboratoire de Physique Théorique (UMR 5152), Université de Toulouse, CNRS, UPS, France.
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA.
| | - Dennis R. Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS – Centre for Molecular Simulation, IQST – Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary2500 University Drive NWCalgaryAlbertaT2N 1N4Canada
| | - Matthias Scheffler
- The NOMAD Laboratory at the FHI of the Max-Planck-Gesellschaft and IRIS-Adlershof of the Humboldt-Universität zu Berlin, Faradayweg 4-6, D-14195, Germany.
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632 Auckland, New Zealand.
| | - Viktor N. Staroverov
- Department of Chemistry, The University of Western OntarioLondonOntario N6A 5B7Canada
| | - Jianwei Sun
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA.
| | - Erik Tellgren
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - David J. Tozer
- Department of Chemistry, Durham UniversitySouth RoadDurhamDH1 3LEUK
| | - Samuel B. Trickey
- Quantum Theory Project, Deptartment of Physics, University of FloridaGainesvilleFL 32611USA
| | - Carsten A. Ullrich
- Department of Physics and Astronomy, University of MissouriColumbiaMO 65211USA
| | - Alberto Vela
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav), CDMX, 07360, Mexico.
| | - Giovanni Vignale
- Department of Physics, University of Missouri, Columbia, MO 65203, USA.
| | - Tomasz A. Wesolowski
- Department of Physical Chemistry, Université de Genève30 Quai Ernest-Ansermet1211 GenèveSwitzerland
| | - Xin Xu
- Shanghai Key Laboratory of Molecular Catalysis and Innovation Materials, Collaborative Innovation Centre of Chemistry for Energy Materials, MOE Laboratory for Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, China.
| | - Weitao Yang
- Department of Chemistry and Physics, Duke University, Durham, NC 27516, USA.
| |
Collapse
|
4
|
Shao X, Lopez AC, Khan Musa MR, Nouri MR, Pavanello M. Adaptive Subsystem Density Functional Theory. J Chem Theory Comput 2022; 18:6646-6655. [PMID: 36179128 DOI: 10.1021/acs.jctc.2c00698] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Subsystem density functional theory (DFT) is emerging as a powerful electronic structure method for large-scale simulations of molecular condensed phases and interfaces. Key to its computational efficiency is the use of approximate nonadditive noninteracting kinetic energy functionals. Unfortunately, currently available nonadditive functionals lead to inaccurate results when the subsystems interact strongly such as when they engage in chemical reactions. This work disrupts the status quo by devising a workflow that extends subsystem DFT's applicability also to strongly interacting subsystems. This is achieved by implementing a fully automated adaptive definition of subsystems which is realized during geometry optimizations or ab initio molecular dynamics simulations. The new method prescribes subsystem merging and splitting events redistributing the resources (both for work and data) in an efficient way making use of modern parallelization strategies and object-oriented programming. We showcase the method with examples probing from moderate-to-strong inter-subsystem interactions, opening the door to using subsystem DFT for modeling chemical reactions in molecular condensed phases with a black box computational tool.
Collapse
Affiliation(s)
- Xuecheng Shao
- Department of Chemistry, Rutgers University, Newark, New Jersey07102, United States
| | | | - Md Rajib Khan Musa
- Department of Chemistry, Rutgers University, Newark, New Jersey07102, United States
| | - Mohammad Reza Nouri
- Department of Chemistry, Rutgers University, Newark, New Jersey07102, United States
| | - Michele Pavanello
- Department of Physics, Rutgers University, Newark, New Jersey07102, United States
| |
Collapse
|
5
|
Chakravarty C, Aksu H, Martinez B JA, Ramos P, Pavanello M, Dunietz BD. Role of Dielectric Screening in Calculating Excited States of Solvated Azobenzene: A Benchmark Study Comparing Quantum Embedding and Polarizable Continuum Model for Representing the Solvent. J Phys Chem Lett 2022; 13:4849-4855. [PMID: 35617015 DOI: 10.1021/acs.jpclett.2c00982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The low energy excited states of the conformational isomers of solvated azobenzene are calculated with several DFT methods accounting for the solute-solvent interaction implicitly with the polarizable continuum model or explicitly with subsystem DFT. For the latter, embedding potentials are calculated for 21 sampled snapshots of the solvent molecules. First, we find that accounting for the solvent implicitly or explicitly has little effect on the predicted cis-trans S1 excitation energy gap. Second, we find that azobenzene's S1 cis and trans energies are accurate as long as a screened range-separated hybrid exchange-correlation functional is employed. Finally, we also tested a simplified workflow whereby a single, averaged, embedding potential is used. Unfortunately, we find larger deviations against the experiment for the simplified workflow. This highlights a basic flaw in the approach, where the time scale of solvent averaging is much longer than that of the solute's electronic polarization.
Collapse
Affiliation(s)
- Chandrima Chakravarty
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, United States
| | - Huseyin Aksu
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, United States
- Computational Physics Laboratory, Department of Physics, Pamukkale University, 20010 Denizli, Turkey
| | - Jessica A Martinez B
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Pablo Ramos
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Barry D Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, United States
| |
Collapse
|
6
|
Wilson L, Krasny R, Luchko T. Accelerating the 3D reference interaction site model theory of molecular solvation with treecode summation and cut-offs. J Comput Chem 2022; 43:1251-1270. [PMID: 35567580 DOI: 10.1002/jcc.26889] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 02/25/2022] [Accepted: 04/15/2022] [Indexed: 11/07/2022]
Abstract
The 3D reference interaction site model (3D-RISM) of molecular solvation is a powerful tool for computing the equilibrium thermodynamics and density distributions of solvents, such as water and co-ions, around solute molecules. However, 3D-RISM solutions can be expensive to calculate, especially for proteins and other large molecules where calculating the potential energy between solute and solvent requires more than half the computation time. To address this problem, we have developed and implemented treecode summation for long-range interactions and analytically corrected cut-offs for short-range interactions to accelerate the potential energy and long-range asymptotics calculations in non-periodic 3D-RISM in the AmberTools molecular modeling suite. For the largest single protein considered in this work, tubulin, the total computation time was reduced by a factor of 4. In addition, parallel calculations with these new methods scale almost linearly and the iterative solver remains the largest impediment to parallel scaling. To demonstrate the utility of our approach for large systems, we used 3D-RISM to calculate the solvation thermodynamics and density distribution of 7-ring microtubule, consisting of 910 tubulin dimers, over 1.2 million atoms.
Collapse
Affiliation(s)
- Leighton Wilson
- Department of Mathematics, University of Michigan, Ann Arbor, Michigan, USA
| | - Robert Krasny
- Department of Mathematics, University of Michigan, Ann Arbor, Michigan, USA
| | - Tyler Luchko
- Department of Physics and Astronomy, California State University, Los Angeles, California, USA
| |
Collapse
|
7
|
Salahub DR. Multiscale molecular modelling: from electronic structure to dynamics of nanosystems and beyond. Phys Chem Chem Phys 2022; 24:9051-9081. [PMID: 35389399 DOI: 10.1039/d1cp05928a] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Important contemporary biological and materials problems often depend on interactions that span orders of magnitude differences in spatial and temporal dimensions. This Tutorial Review attempts to provide an introduction to such fascinating problems through a series of case studies, aimed at beginning researchers, graduate students, postdocs and more senior colleagues who are changing direction to focus on multiscale aspects of their research. The choice of specific examples is highly personal, with examples either chosen from our own work or outstanding multiscale efforts from the literature. I start with various embedding schemes, as exemplified by polarizable continuum models, 3-D RISM, molecular DFT and frozen-density embedding. Next, QM/MM (quantum mechanical/molecular mechanical) techniques are the workhorse of pm-to-nm/ps-to-ns simulations; examples are drawn from enzymes and from nanocatalysis for oil-sands upgrading. Using polarizable force-fields in the QM/MM framework represents a burgeoning subfield; with examples from ion channels and electron dynamics in molecules subject to strong external fields, probing the atto-second dynamics of the electrons with RT-TDDFT (real-time - time-dependent density functional theory) eventually coupled with nuclear motion through the Ehrenfest approximation. This is followed by a section on coarse graining, bridging dimensions from atoms to cells. The penultimate chapter gives a quick overview of multiscale approaches that extend into the meso- and macro-scales, building on atomistic and coarse-grained techniques to enter the world of materials engineering, on the one hand, and cell biology, on the other. A final chapter gives just a glimpse of the burgeoning impact of machine learning on the structure-dynamics front. I aim to capture the excitement of contemporary leading-edge breakthroughs in the description of physico-chemical systems and processes in complex environments, with only enough historical content to provide context and aid the next generation of methodological development. While I aim also for a clear description of the essence of methodological breakthroughs, equations are kept to a minimum and detailed formalism and implementation details are left to the references. My approach is very selective (case studies) rather than exhaustive. I think that these case studies should provide fodder to build as complete a reference tree on multiscale modelling as the reader may wish, through forward and backward citation analysis. I hope that my choices of cases will excite interest in newcomers and help to fuel the growth of multiscale modelling in general.
Collapse
Affiliation(s)
- Dennis R Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS-Centre for Molecular Simulation, IQST-Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, Calgary, Alberta, T2N 1N4, Canada.
| |
Collapse
|
8
|
An MP2/Molecular Dynamics study of the solvent effects on the conformational equilibrium of the glycine dipeptide. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
9
|
Tölle J, Neugebauer J. The Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory. J Phys Chem Lett 2022; 13:1003-1018. [PMID: 35061387 DOI: 10.1021/acs.jpclett.1c04023] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
The theoretical understanding of photoinduced processes in multichromophoric systems requires, as an essential ingredient, the possibility of accurately describing their electronically excited states. However, the size of these systems often prohibits the usage of conventional electronic-structure methods, so that often multiscale approaches based on phenomenologically motivated models are employed. In contrast, subsystem time-dependent density functional theory (sTDDFT) allows for a subsystem-based ab initio description of multichromophoric systems and therefore allows for, in principle, an exact description of photoinduced processes. This Perspective aims to outline the theoretical foundations and commonly used practical realizations as well as to illustrate benefits of recent developments and open issues in the field of sTDDFT. Prospective, potential future applications and possible methodological developments are discussed.
Collapse
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
10
|
González-Espinoza CE, Rumble CA, Borgis D, Wesolowski TA. Quantifying Fluctuations of Average Solvent Environments for Embedding Calculations. J Chem Theory Comput 2022; 18:1072-1088. [PMID: 35044168 DOI: 10.1021/acs.jctc.1c01108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The viability and effectiveness of replacing an ensemble of embedded solute calculations by a single calculation using an average description of the solvent environment are evaluated. This work explores the fluctuations of the average description of the system obtained in two ways: from calculations on an ensemble of geometries and from an average environment constructed from the same ensemble. To this end, classical molecular dynamics simulations of a rigid acetone solute in SPCE water are performed in order to generate an ensemble of solvent environments. From this ensemble of solvent configurations, a number of different approaches for constructing an average solvent environment are employed. We perform a thorough numerical analysis of the fluctuations of the electrostatic potential experienced by the solute, as well as the resulting fluctuations of the solute's electronic density, measured through its dipole moment and fitted atomic point charges. At the same time, we inspect the accuracy of the methods used to construct average environments. Finally, the proposed method for generating the embedding potential from an average environment density is applied to estimate the solvatochromic shift of the first excitation of acetone. In order to account for quantum confinement effects, which may be important in certain cases, the fluctuations in the shift due to the interaction with the solvent are evaluated using frozen-density-embedding theory. Our results demonstrate that, for normally distributed environments, the constructed average environment is a reasonably good representation of a fluctuating molecular solvent environment. We then provide guidance for future comparisons between these theoretical treatments of solute/solvent systems to experimental measurements.
Collapse
Affiliation(s)
| | - Christopher A Rumble
- Université de Genève, Départment de Chimie Physique 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Daniel Borgis
- Ecole Normale Supérieure, Départment de Chimie, UMR 8640, PSL University, Sorbonne University, CNRS, 75005 Paris, France
| | - Tomasz A Wesolowski
- Université de Genève, Départment de Chimie Physique 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
11
|
Bensberg M, Türtscher PL, Unsleber JP, Reiher M, Neugebauer J. Solvation Free Energies in Subsystem Density Functional Theory. J Chem Theory Comput 2022; 18:723-740. [PMID: 34985890 DOI: 10.1021/acs.jctc.1c00864] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
For many chemical processes the accurate description of solvent effects are vitally important. Here, we describe a hybrid ansatz for the explicit quantum mechanical description of solute-solvent and solvent-solvent interactions based on subsystem density functional theory and continuum solvation schemes. Since explicit solvent molecules may compromise the scalability of the model and transferability of the predicted solvent effect, we aim to retain both, for different solutes as well as for different solvents. The key for the transferability is the consistent subsystem decomposition of solute and solvent. The key for the scalability is the performance of subsystem DFT for increasing numbers of subsystems. We investigate molecular dynamics and stationary point sampling of solvent configurations and compare the resulting (Gibbs) free energies to experiment and theoretical methods. We can show that with our hybrid model reaction barriers and reaction energies are accurately reproduced compared to experimental data.
Collapse
Affiliation(s)
- Moritz Bensberg
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Paul L Türtscher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Jan P Unsleber
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
12
|
Benchmarking Free Energy Calculations in Liquid Aliphatic Ketone Solvents Using the 3D-RISM-KH Molecular Solvation Theory. J 2021. [DOI: 10.3390/j4040044] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
The three-dimensional reference interaction site model of the molecular solvation theory with the Kovalenko–Hirata closure is used to calculate the free energy of solvation of organic solutes in liquid aliphatic ketones. The ketone solvent sites were modeled using a modified united-atom force field. The successful application of these solvation models in calculating ketone–water partition coefficients of a large number of solutes supports the validation and benchmarking reported here.
Collapse
|
13
|
Mroginski MA, Adam S, Amoyal GS, Barnoy A, Bondar AN, Borin VA, Church JR, Domratcheva T, Ensing B, Fanelli F, Ferré N, Filiba O, Pedraza-González L, González R, González-Espinoza CE, Kar RK, Kemmler L, Kim SS, Kongsted J, Krylov AI, Lahav Y, Lazaratos M, NasserEddin Q, Navizet I, Nemukhin A, Olivucci M, Olsen JMH, Pérez de Alba Ortíz A, Pieri E, Rao AG, Rhee YM, Ricardi N, Sen S, Solov'yov IA, De Vico L, Wesolowski TA, Wiebeler C, Yang X, Schapiro I. Frontiers in Multiscale Modeling of Photoreceptor Proteins. Photochem Photobiol 2021; 97:243-269. [PMID: 33369749 DOI: 10.1111/php.13372] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 11/01/2020] [Indexed: 02/06/2023]
Abstract
This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.
Collapse
Affiliation(s)
| | - Suliman Adam
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Gil S Amoyal
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Avishai Barnoy
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Ana-Nicoleta Bondar
- Freie Universität Berlin, Department of Physics, Theoretical Molecular Biophysics Group, Berlin, Germany
| | - Veniamin A Borin
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Jonathan R Church
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Tatiana Domratcheva
- Department of Chemistry, Lomonosov Moscow State University, Moscow, Russia.,Department Biomolecular Mechanisms, Max Planck Institute for Medical Research, Heidelberg, Germany
| | - Bernd Ensing
- Van 't Hoff Institute for Molecular Science and Amsterdam Center for Multiscale Modeling, University of Amsterdam, Amsterdam, The Netherlands
| | - Francesca Fanelli
- Department of Life Sciences, Center for Neuroscience and Neurotechnology, Università degli Studi di Modena e Reggio Emilia, Modena, Italy
| | | | - Ofer Filiba
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Laura Pedraza-González
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Siena, Italy
| | - Ronald González
- Institut für Chemie, Technische Universität Berlin, Berlin, Germany
| | | | - Rajiv K Kar
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Lukas Kemmler
- Freie Universität Berlin, Department of Physics, Theoretical Molecular Biophysics Group, Berlin, Germany
| | - Seung Soo Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, Korea
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Odense, Denmark
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, CA, USA
| | - Yigal Lahav
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel.,MIGAL - Galilee Research Institute, S. Industrial Zone, Kiryat Shmona, Israel
| | - Michalis Lazaratos
- Freie Universität Berlin, Department of Physics, Theoretical Molecular Biophysics Group, Berlin, Germany
| | - Qays NasserEddin
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Isabelle Navizet
- MSME, Univ Gustave Eiffel, CNRS UMR 8208, Univ Paris Est Creteil, Marne-la-Vallée, France
| | - Alexander Nemukhin
- Department of Chemistry, Lomonosov Moscow State University, Moscow, Russia.,Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, Russia
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Siena, Italy.,Chemistry Department, Bowling Green State University, Bowling Green, OH, USA
| | - Jógvan Magnus Haugaard Olsen
- Department of Chemistry, Aarhus University, Aarhus, Denmark.,Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, Tromsø, Norway
| | - Alberto Pérez de Alba Ortíz
- Van 't Hoff Institute for Molecular Science and Amsterdam Center for Multiscale Modeling, University of Amsterdam, Amsterdam, The Netherlands
| | - Elisa Pieri
- Aix-Marseille Univ, CNRS, ICR, Marseille, France
| | - Aditya G Rao
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, Korea
| | - Niccolò Ricardi
- Département de Chimie Physique, Université de Genève, Genève, Switzerland
| | - Saumik Sen
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Ilia A Solov'yov
- Department of Physics, Carl von Ossietzky University Oldenburg, Oldenburg, Germany
| | - Luca De Vico
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Siena, Italy
| | | | - Christian Wiebeler
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Xuchun Yang
- Chemistry Department, Bowling Green State University, Bowling Green, OH, USA
| | - Igor Schapiro
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| |
Collapse
|
14
|
Jeanmairet G, Levesque M, Borgis D. Tackling Solvent Effects by Coupling Electronic and Molecular Density Functional Theory. J Chem Theory Comput 2020; 16:7123-7134. [PMID: 32894674 DOI: 10.1021/acs.jctc.0c00729] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Solvation effects can have a tremendous influence on chemical reactions. However, precise quantum chemistry calculations are most often done either in vacuum neglecting the role of the solvent or using continuum solvent model ignoring its molecular nature. We propose a new method coupling a quantum description of the solute using electronic density functional theory with a classical grand-canonical treatment of the solvent using molecular density functional theory. Unlike a previous work, both densities are minimized self-consistently, accounting for mutual polarization of the molecular solvent and the solute. The electrostatic interaction is accounted using the full electron density of the solute rather than fitted point charges. The introduced methodology represents a good compromise between the two main strategies to tackle solvation effects in quantum calculation. It is computationally more effective than a direct quantum mechanics/molecular mechanics coupling, requiring the exploration of many solvent configurations. Compared to continuum methods, it retains the full molecular-level description of the solvent. We validate this new framework onto two usual benchmark systems: a water solvated in water and the symmetrical nucleophilic substitution between chloromethane and chloride in water. The prediction for the free energy profiles are not yet fully quantitative compared to experimental data, but the most important features are qualitatively recovered. The method provides a detailed molecular picture of the evolution of the solvent structure along the reaction pathway.
Collapse
Affiliation(s)
- Guillaume Jeanmairet
- Sorbonne Université, CNRS, Physico-Chimie des Électrolytes et Nanosystèmes, Interfaciaux, PHENIX, F-75005 Paris, France.,Réseau sur le Stockage Électrochimique de l'Énergie (RS2E), FR CNRS 3459, 80039 Amiens Cedex, France
| | - Maximilien Levesque
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne, Université, CNRS, 75005 Paris, France.,Aqemia, 75006 Paris, France
| | - Daniel Borgis
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne, Université, CNRS, 75005 Paris, France.,Maison de la Simulation, CEA, CNRS, Université Paris-Sud, UVSQ, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| |
Collapse
|
15
|
Reimann M, Kaupp M. Evaluation of an Efficient 3D-RISM-SCF Implementation as a Tool for Computational Spectroscopy in Solution. J Phys Chem A 2020; 124:7439-7452. [PMID: 32838530 DOI: 10.1021/acs.jpca.0c06322] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The 3D-RISM-SCF solvent-model implementation of Gusarov et al. [ J. Phys. Chem. A 2006, 110, 6083-6090] in the Amsterdam density functional program has been improved and extended. In particular, an accurate yet efficient representation of the solute electrostatic potential is provided. The Coulomb-potential fitting of many DFT codes can be used advantageously in this context. The extra effort compared to a point-charge representation is small for a given SCF cycle and compensated by faster SCF convergence. This allows applications to large solutes, as demonstrated by evaluation of the solvatochromism of Reichardt's dye. In general, TDDFT applications to excitation energies in solution stand out and are highlighted. Applications to the 17O NMR chemical shifts of N-methylformamide in different solvents also demonstrate the distinct advantages of 3D-RISM over continuum solvents. Limitations are observed in this case for water solvent, where the solvent shielding is overestimated. This shortcoming applies also to the 17O gas-to-liquid shift of water, where we used localized molecular orbital analyses for a deeper understanding. For such cases of extremely strong solute-solvent interactions, couplings between solute and solvent orbitals induced by the magnetic perturbation are relevant. These clearly require a quantum-mechanical treatment of the most closely bound solvent molecules. Except for such extreme cases, 3D-RISM-SCF is very well suited to treat solvent effects on NMR parameters. More serious limitations pertain to the treatment of vibrational spectra, where the absence of the coupling between solute and solvent vibrational modes limits the accuracy of applications of 3D-RISM-SCF. The reported extended, efficient, and numerically accurate 3D-RISM-SCF implementation should provide a useful tool to study chemical and spectroscopic properties of molecules of appreciable size in a realistic solvent environment.
Collapse
Affiliation(s)
- Marc Reimann
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623, Berlin, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623, Berlin, Germany
| |
Collapse
|
16
|
Ricardi N, Ernst M, Macchi P, Wesolowski TA. Embedding-theory-based simulations using experimental electron densities for the environment. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2020; 76:571-579. [PMID: 32869754 PMCID: PMC7459768 DOI: 10.1107/s2053273320008062] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/05/2020] [Accepted: 06/16/2020] [Indexed: 11/26/2022]
Abstract
For the first time, the use of experimentally derived molecular electron densities as ρB(r) in calculations based on frozen-density embedding theory (FDET) of environment-induced shifts of electronic excitations for chromophores in clusters is demonstrated. ρB(r) was derived from X-ray restrained molecular wavefunctions of glycylglycine to obtain environment densities for simulating electronic excitations in clusters. The basic idea of frozen-density embedding theory (FDET) is the constrained minimization of the Hohenberg–Kohn density functional EHK[ρ] performed using the auxiliary functional , where ΨA is the embedded NA-electron wavefunction and ρB(r) is a non-negative function in real space integrating to a given number of electrons NB. This choice of independent variables in the total energy functional makes it possible to treat the corresponding two components of the total density using different methods in multi-level simulations. The application of FDET using ρB(r) reconstructed from X-ray diffraction data for a molecular crystal is demonstrated for the first time. For eight hydrogen-bonded clusters involving a chromophore (represented as ΨA) and the glycylglycine molecule [represented as ρB(r)], FDET is used to derive excitation energies. It is shown that experimental densities are suitable for use as ρB(r) in FDET-based simulations.
Collapse
Affiliation(s)
- Niccolò Ricardi
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Michelle Ernst
- University of Bern, Freiestraße 3, 3012 Bern, Switzerland
| | - Piero Macchi
- Department of Chemistry, Materials and Chemical Engineering, Polytechnic of Milan, via Mancinelli 7, Milano 20131, Italy
| | - Tomasz Adam Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
17
|
Schwarz K, Sundararaman R. The electrochemical interface in first-principles calculations. SURFACE SCIENCE REPORTS 2020; 75:10.1016/j.surfrep.2020.100492. [PMID: 34194128 PMCID: PMC8240516 DOI: 10.1016/j.surfrep.2020.100492] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
First-principles predictions play an important role in understanding chemistry at the electrochemical interface. Electronic structure calculations are straightforward for vacuum interfaces, but do not easily account for the interfacial fields and solvation that fundamentally change the nature of electrochemical reactions. Prevalent techniques for first-principles prediction of electrochemical processes range from expensive explicit solvation using ab initio molecular dynamics, through a hierarchy of continuum solvation techniques, to neglecting solvation and interfacial field effects entirely. Currently, no single approach reliably captures all relevant effects of the electrochemical double layer in first-principles calculations. This review systematically lays out the relation between all major approaches to first-principles electrochemistry, including the key approximations and their consequences for accuracy and computational cost. Focusing on ab initio methods for thermodynamic properties of aqueous interfaces, we first outline general considerations for modeling electrochemical interfaces, including solvent and electrolyte dynamics and electrification. We then present the specifics of various explicit and implicit models of the solvent and electrolyte. Finally, we discuss the compromise between computational efficiency and accuracy, and identify key outstanding challenges and future opportunities in the wide range of techniques for first-principles electrochemistry.
Collapse
Affiliation(s)
- Kathleen Schwarz
- Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Dr., Gaithersburg, Maryland 20899, USA
| | - Ravishankar Sundararaman
- Department of Materials Science and Engineering, Rensselaer Polytechnic Institute, 110 8th St., Troy, New York 12180, USA
| |
Collapse
|
18
|
Omelyan I, Kovalenko A. Enhanced solvation force extrapolation for speeding up molecular dynamics simulations of complex biochemical liquids. J Chem Phys 2019; 151:214102. [PMID: 31822083 DOI: 10.1063/1.5126410] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
We propose an enhanced approach to the extrapolation of mean potential forces acting on atoms of solute macromolecules due to their interactions with solvent atoms in complex biochemical liquids. It improves and extends our previous extrapolation schemes by additionally including new techniques such as an exponential scaling transformation of coordinate space with weights complemented by an automatically adjusted balancing between the least square minimization of force deviations and the norm of expansion coefficients in the approximation. The expensive mean potential forces are treated in terms of the three-dimensional reference interaction site model with Kovalenko-Hirata closure molecular theory of solvation. During the dynamics, they are calculated only after every long (outer) time interval, i.e., quite rarely to reduce the computational costs. At much shorter (inner) time steps, these forces are extrapolated on the basis of their outer values. The equations of motion are then solved using a multiple time step integration within an optimized isokinetic Nosé-Hoover chain thermostat. The new approach is applied to molecular dynamics simulations of various systems consisting of solvated organic and biomolecules of different complexity. For example, we consider hydrated alanine dipeptide, asphaltene in toluene solvent, miniprotein 1L2Y, and protein G in aqueous solution. It is shown that in all these cases, the enhanced extrapolation provides much better accuracy of the solvation force approximation than the existing approaches. As a result, it can be used with much larger outer time steps, leading to a significant speedup of the simulations.
Collapse
Affiliation(s)
- Igor Omelyan
- Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine, 1 Svientsitskii Street, Lviv 79011, Ukraine
| | - Andriy Kovalenko
- Department of Mechanical Engineering, University of Alberta, 9211-116 Street NW, Edmonton, Alberta T6G 1H9, Canada
| |
Collapse
|
19
|
Prentice JCA, Charlton RJ, Mostofi AA, Haynes PD. Combining Embedded Mean-Field Theory with Linear-Scaling Density-Functional Theory. J Chem Theory Comput 2019; 16:354-365. [DOI: 10.1021/acs.jctc.9b00956] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Joseph C. A. Prentice
- Department of Materials, Department of Physics, and the Thomas Young Centre for Theory and Simulation of Materials, Imperial College London, London SW7 2AZ, United Kingdom
| | - Robert J. Charlton
- Department of Materials, Department of Physics, and the Thomas Young Centre for Theory and Simulation of Materials, Imperial College London, London SW7 2AZ, United Kingdom
| | - Arash A. Mostofi
- Department of Materials, Department of Physics, and the Thomas Young Centre for Theory and Simulation of Materials, Imperial College London, London SW7 2AZ, United Kingdom
| | - Peter D. Haynes
- Department of Materials, Department of Physics, and the Thomas Young Centre for Theory and Simulation of Materials, Imperial College London, London SW7 2AZ, United Kingdom
| |
Collapse
|
20
|
Blinov N, Wishart DS, Kovalenko A. Solvent Composition Effects on the Structural Properties of the Aβ42 Monomer from the 3D-RISM-KH Molecular Theory of Solvation. J Phys Chem B 2019; 123:2491-2506. [PMID: 30811210 DOI: 10.1021/acs.jpcb.9b00480] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Structural characterization of amyloid (A)β peptides implicated in Alzheimer's disease is a challenging problem due to their intrinsically disordered nature and their high propensity for aggregation. Only limited information is currently available from experiments on conformational properties and aggregation pathways of the peptides in cellular environments. In silico modeling complements experimental information, providing atomistic insight into structure and dynamics of different Aβ species. All-atom explicit solvent molecular dynamics (MD) simulations with a properly selected force field can deliver reliable structural and dynamic information. In the case of intrinsically disordered Aβ peptides, enhanced sampling simulations beyond the nanosecond time scale are required to obtain statistically meaningful results even for simple solvent conditions. To overcome the challenges of conformational sampling in crowded cellular environments, alternative approaches have to be used, including postprocessing of MD data. In this study, we employ the statistical-mechanical, three-dimensional reference interaction site model with the Kovalenko-Hirata closure integral equation molecular theory of solvation to describe solvent composition effects on the conformational equilibrium in a structural ensemble of the Aβ42 (covering residues 1-42) monomer based on a statistical reweighting technique. The methodology enables a computationally efficient prediction on how different factors in the cellular environment, such as solvent composition, nonpolar solvation, and macromolecular crowding, affect the structural properties of the monomer. Similarities have been identified between changes in the structural ensemble caused by nonpolar solvation and crowded environments modeled by ionic solution with large negative ions. In particular, both solvent conditions reduce the random coil content and enhance the helical structure content of the monomer. In contrast to the previous studies, which reported increased α-helical content of peptides in crowded environments, this work attributes these structural features to the difference in solvent exposure of hydrophilic residues of the monomer for different secondary structure elements, rather than to (entropic) excluded volume effects.
Collapse
Affiliation(s)
- Nikolay Blinov
- Department of Mechanical Engineering , Edmonton , Alberta T6G 1H9 , Canada.,Nanotechnology Research Centre , Edmonton , Alberta T6G 2M9 , Canada
| | - David S Wishart
- Departments of Computing Science and Biological Sciences , University of Alberta , Edmonton , Alberta T6G 2E8 , Canada
| | - Andriy Kovalenko
- Department of Mechanical Engineering , Edmonton , Alberta T6G 1H9 , Canada.,Nanotechnology Research Centre , Edmonton , Alberta T6G 2M9 , Canada
| |
Collapse
|
21
|
Affiliation(s)
- Basile F. E. Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry and PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
22
|
Kovalenko A, Gusarov S. Multiscale methods framework: self-consistent coupling of molecular theory of solvation with quantum chemistry, molecular simulations, and dissipative particle dynamics. Phys Chem Chem Phys 2018; 20:2947-2969. [DOI: 10.1039/c7cp05585d] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
In this work, we will address different aspects of self-consistent field coupling of computational chemistry methods at different time and length scales in modern materials and biomolecular science.
Collapse
Affiliation(s)
- Andriy Kovalenko
- National Institute for Nanotechnology
- National Research Council of Canada
- Edmonton
- Canada
- Department of Mechanical Engineering
| | - Sergey Gusarov
- National Institute for Nanotechnology
- National Research Council of Canada
- Edmonton
- Canada
| |
Collapse
|
23
|
Howard JC, Womack JC, Dziedzic J, Skylaris CK, Pritchard BP, Crawford TD. Electronically Excited States in Solution via a Smooth Dielectric Model Combined with Equation-of-Motion Coupled Cluster Theory. J Chem Theory Comput 2017; 13:5572-5581. [DOI: 10.1021/acs.jctc.7b00833] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- J. Coleman Howard
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - James C. Womack
- Department
of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, U.K
| | - Jacek Dziedzic
- Department
of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, U.K
- Faculty
of Applied Physics and Mathematics, Gdańsk University of Technology, Gdańsk 80-233, Poland
| | - Chris-Kriton Skylaris
- Department
of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, U.K
| | - Benjamin P. Pritchard
- Molecular Sciences Software Institute, Virginia Tech, Blacksburg, Virginia 24060, United States
| | - T. Daniel Crawford
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
24
|
Ding L, Levesque M, Borgis D, Belloni L. Efficient molecular density functional theory using generalized spherical harmonics expansions. J Chem Phys 2017; 147:094107. [DOI: 10.1063/1.4994281] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Lu Ding
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Maximilien Levesque
- PASTEUR, Département de Chimie, École Normale Supérieure, UPMC Univ. Paris 06, CNRS, PSL Research University, 75005 Paris, France
- Sorbonne Universités, UPMC Univ. Paris 06, École Normale Supérieure, CNRS, Processus d’Activation Sélective par Transfert d’Énergie Uni-Électronique ou Radiatif (PASTEUR), 75005 Paris, France
| | - Daniel Borgis
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
- PASTEUR, Département de Chimie, École Normale Supérieure, UPMC Univ. Paris 06, CNRS, PSL Research University, 75005 Paris, France
- Sorbonne Universités, UPMC Univ. Paris 06, École Normale Supérieure, CNRS, Processus d’Activation Sélective par Transfert d’Énergie Uni-Électronique ou Radiatif (PASTEUR), 75005 Paris, France
| | - Luc Belloni
- LIONS, NIMBE, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| |
Collapse
|
25
|
Abstract
In complex systems, it is often the case that the region of interest forms only one part of a much larger system. The idea of joining two different quantum simulations-a high level calculation on the active region of interest, and a low level calculation on its environment-formally defines a quantum embedding. While any combination of techniques constitutes an embedding, several rigorous formalisms have emerged that provide for exact feedback between the embedded system and its environment. These three formulations: density functional embedding, Green's function embedding, and density matrix embedding, respectively, use the single-particle density, single-particle Green's function, and single-particle density matrix as the quantum variables of interest. Many excellent reviews exist covering these methods individually. However, a unified presentation of the different formalisms is so far lacking. Indeed, the various languages commonly used, functional equations for density functional embedding, diagrammatics for Green's function embedding, and entanglement arguments for density matrix embedding, make the three formulations appear vastly different. In this Account, we introduce the basic equations of all three formulations in such a way as to highlight their many common intellectual strands. While we focus primarily on a straightforward theoretical perspective, we also give a brief overview of recent applications and possible future developments. The first section starts with density functional embedding, where we introduce the key embedding potential via the Euler equation. We then discuss recent work concerning the treatment of the nonadditive kinetic potential, before describing mean-field density functional embedding and wave function in density functional embedding. We finish the section with extensions to time-dependence and excited states. The second section is devoted to Green's function embedding. Here, we use the Dyson equation to obtain equations that parallel as closely as possible the density functional embedding equations, with the hybridization playing the role of the embedding potential. Embedding a high-level self-energy within a low-level self-energy is treated analogously to wave function in density functional embedding. The numerical computation of the high-level self-energy allows us to briefly introduce the bath representation in the quantum impurity problem. We then consider translationally invariant systems to bring in the important dynamical mean-field theory. Recent developments to incorporate screening and long-range interactions are discussed. The third section concerns density matrix embedding. Here, we first highlight some mathematical complications associated with a simple Euler equation derivation, arising from the open nature of fragments. This motivates the density matrix embedding theory, where we use the Schmidt decomposition to represent the entanglement through bath orbitals. The resulting impurity plus bath formulation resembles that of dynamical mean-field theory. We discuss the numerical self-consistency associated with using a high-level correlated wave function with a mean-field low-level treatment, and connect the resulting numerical inversion to that used in density functional embedding. We finish with perspectives on the future of all three methods.
Collapse
Affiliation(s)
- Qiming Sun
- Division of Chemistry and
Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and
Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
26
|
Schlüns D, Franchini M, Götz AW, Neugebauer J, Jacob CR, Visscher L. Analytical gradients for subsystem density functional theory within the slater-function-based amsterdam density functional program. J Comput Chem 2016; 38:238-249. [DOI: 10.1002/jcc.24670] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2016] [Revised: 10/14/2016] [Accepted: 10/19/2016] [Indexed: 11/09/2022]
Affiliation(s)
- Danny Schlüns
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster; Corrensstraße 40 Münster 48149 Germany
| | - Mirko Franchini
- Amsterdam Center for Multiscale Modeling; Vrije Universiteit Amsterdam; De Boelelaan 1083 HV Amsterdam 1081 The Netherlands
- Scientific Computing & Modelling NV; Vrije Universiteit; Theoretical Chemistry, De Boelelaan 1083 HV Amsterdam 1081 The Netherlands
| | - Andreas W. Götz
- San Diego Supercomputer Center; University of California San Diego; 9500 Gilman Drive La Jolla, California 92093-0505 USA
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster; Corrensstraße 40 Münster 48149 Germany
| | - Christoph R. Jacob
- Institute of Physical and Theoretical Chemistry; TU Braunschweig, Hans-Sommer-Straße 10 Braunschweig 38106 Germany
| | - Lucas Visscher
- Amsterdam Center for Multiscale Modeling; Vrije Universiteit Amsterdam; De Boelelaan 1083 HV Amsterdam 1081 The Netherlands
| |
Collapse
|
27
|
Li W, Li Y, Lin R, Li S. Generalized Energy-Based Fragmentation Approach for Localized Excited States of Large Systems. J Phys Chem A 2016; 120:9667-9677. [DOI: 10.1021/acs.jpca.6b11193] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Wei Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Repubic of China
| | - Yunzhi Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Repubic of China
| | - Ruochen Lin
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Repubic of China
| | - Shuhua Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Repubic of China
| |
Collapse
|
28
|
Kobryn AE, Gusarov S, Kovalenko A. A closure relation to molecular theory of solvation for macromolecules. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:404003. [PMID: 27549008 DOI: 10.1088/0953-8984/28/40/404003] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We propose a closure to the integral equations of molecular theory of solvation, particularly suitable for polar and charged macromolecules in electrolyte solution. This includes such systems as oligomeric polyelectrolytes at a finite concentration in aqueous and various non-aqueous solutions, as well as drug-like compounds in solution. The new closure by Kobryn, Gusarov, and Kovalenko (KGK closure) imposes the mean spherical approximation (MSA) almost everywhere in the solvation shell but levels out the density distribution function to zero (with the continuity at joint boundaries) inside the repulsive core and in the spatial regions of strong density depletion emerging due to molecular associative interactions. Similarly to MSA, the KGK closure reduces the problem to a linear equation for the direct correlation function which is predefined analytically on most of the solvation shells and has to be determined numerically on a relatively small (three-dimensional) domain of strong depletion, typically within the repulsive core. The KGK closure leads to the solvation free energy in the form of the Gaussian fluctuation (GF) functional. We first test the performance of the KGK closure coupled to the reference interaction site model (RISM) integral equations on the examples of Lennard-Jones liquids, polar and nonpolar molecular solvents, including water, and aqueous solutions of simple ions. The solvation structure, solvation chemical potential, and compressibility obtained from RISM with the KGK closure favorably compare to the results of the hypernetted chain (HNC) and Kovalenko-Hirata (KH) closures, including their combination with the GF solvation free energy. We then use the KGK closure coupled to RISM to obtain the solvation structure and thermodynamics of oligomeric polyelectrolytes and drug-like compounds at a finite concentration in electrolyte solution, for which no convergence is obtained with other closures. For comparison, we calculate their solvation structure from molecular dynamics (MD) simulations. We further couple the 3D-RISM integral equation with the 3D-version of the KGK closure, and solve it for molecular mixtures as well as oligomeric polyelectrolytes and drug-like molecules in electrolyte solutions.
Collapse
Affiliation(s)
- Alexander E Kobryn
- National Institute for Nanotechnology, National Research Council Canada, 11421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada
| | | | | |
Collapse
|
29
|
Omelyan I, Kovalenko A. MTS-MD of Biomolecules Steered with 3D-RISM-KH Mean Solvation Forces Accelerated with Generalized Solvation Force Extrapolation. J Chem Theory Comput 2016; 11:1875-95. [PMID: 26574393 DOI: 10.1021/ct5010438] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
We developed a generalized solvation force extrapolation (GSFE) approach to speed up multiple time step molecular dynamics (MTS-MD) of biomolecules steered with mean solvation forces obtained from the 3D-RISM-KH molecular theory of solvation (three-dimensional reference interaction site model with the Kovalenko-Hirata closure). GSFE is based on a set of techniques including the non-Eckart-like transformation of coordinate space separately for each solute atom, extension of the force-coordinate pair basis set followed by selection of the best subset, balancing the normal equations by modified least-squares minimization of deviations, and incremental increase of outer time step in motion integration. Mean solvation forces acting on the biomolecule atoms in conformations at successive inner time steps are extrapolated using a relatively small number of best (closest) solute atomic coordinates and corresponding mean solvation forces obtained at previous outer time steps by converging the 3D-RISM-KH integral equations. The MTS-MD evolution steered with GSFE of 3D-RISM-KH mean solvation forces is efficiently stabilized with our optimized isokinetic Nosé-Hoover chain (OIN) thermostat. We validated the hybrid MTS-MD/OIN/GSFE/3D-RISM-KH integrator on solvated organic and biomolecules of different stiffness and complexity: asphaltene dimer in toluene solvent, hydrated alanine dipeptide, miniprotein 1L2Y, and protein G. The GSFE accuracy and the OIN efficiency allowed us to enlarge outer time steps up to huge values of 1-4 ps while accurately reproducing conformational properties. Quasidynamics steered with 3D-RISM-KH mean solvation forces achieves time scale compression of conformational changes coupled with solvent exchange, resulting in further significant acceleration of protein conformational sampling with respect to real time dynamics. Overall, this provided a 50- to 1000-fold effective speedup of conformational sampling for these systems, compared to conventional MD with explicit solvent. We have been able to fold the miniprotein from a fully denatured, extended state in about 60 ns of quasidynamics steered with 3D-RISM-KH mean solvation forces, compared to the average physical folding time of 4-9 μs observed in experiment.
Collapse
Affiliation(s)
- Igor Omelyan
- Department of Mechanical Engineering, University of Alberta , Mechanical Engineering Building 4-9, Edmonton, Alberta T6G 2G8, Canada.,National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada.,Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine , 1 Svientsitskii Street, Lviv 79011, Ukraine
| | - Andriy Kovalenko
- Department of Mechanical Engineering, University of Alberta , Mechanical Engineering Building 4-9, Edmonton, Alberta T6G 2G8, Canada.,National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| |
Collapse
|
30
|
SAMPL5: 3D-RISM partition coefficient calculations with partial molar volume corrections and solute conformational sampling. J Comput Aided Mol Des 2016; 30:1115-1127. [PMID: 27585474 DOI: 10.1007/s10822-016-9947-7] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2016] [Accepted: 08/18/2016] [Indexed: 12/15/2022]
Abstract
Implicit solvent methods for classical molecular modeling are frequently used to provide fast, physics-based hydration free energies of macromolecules. Less commonly considered is the transferability of these methods to other solvents. The Statistical Assessment of Modeling of Proteins and Ligands 5 (SAMPL5) distribution coefficient dataset and the accompanying explicit solvent partition coefficient reference calculations provide a direct test of solvent model transferability. Here we use the 3D reference interaction site model (3D-RISM) statistical-mechanical solvation theory, with a well tested water model and a new united atom cyclohexane model, to calculate partition coefficients for the SAMPL5 dataset. The cyclohexane model performed well in training and testing ([Formula: see text] for amino acid neutral side chain analogues) but only if a parameterized solvation free energy correction was used. In contrast, the same protocol, using single solute conformations, performed poorly on the SAMPL5 dataset, obtaining [Formula: see text] compared to the reference partition coefficients, likely due to the much larger solute sizes. Including solute conformational sampling through molecular dynamics coupled with 3D-RISM (MD/3D-RISM) improved agreement with the reference calculation to [Formula: see text]. Since our initial calculations only considered partition coefficients and not distribution coefficients, solute sampling provided little benefit comparing against experiment, where ionized and tautomer states are more important. Applying a simple [Formula: see text] correction improved agreement with experiment from [Formula: see text] to [Formula: see text], despite a small number of outliers. Better agreement is possible by accounting for tautomers and improving the ionization correction.
Collapse
|
31
|
Partial molar volume of nonionic surfactants in aqueous solution studied by the KB/3D-RISM–KH theory. J Mol Liq 2016. [DOI: 10.1016/j.molliq.2016.02.016] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
32
|
Krishtal A, Pavanello M. Revealing electronic open quantum systems with subsystem TDDFT. J Chem Phys 2016; 144:124118. [PMID: 27036438 DOI: 10.1063/1.4944526] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Open quantum systems (OQSs) are perhaps the most realistic systems one can approach through simulations. In recent years, describing OQSs with Density Functional Theory (DFT) has been a prominent avenue of research with most approaches based on a density matrix partitioning in conjunction with an ad-hoc description of system-bath interactions. We propose a different theoretical approach to OQSs based on partitioning of the electron density. Employing the machinery of subsystem DFT (and its time-dependent extension), we provide a novel way of isolating and analyzing the various terms contributing to the coupling between the system and the surrounding bath. To illustrate the theory, we provide numerical simulations on a toy system (a molecular dimer) and on a condensed phase system (solvated excimer). The simulations show that non-Markovian dynamics in the electronic system-bath interactions are important in chemical applications. For instance, we show that the superexchange mechanism of transport in donor-bridge-acceptor systems is a non-Markovian interaction between the donor-acceptor (OQS) with the bridge (bath) which is fully characterized by real-time subsystem time-dependent DFT.
Collapse
Affiliation(s)
- Alisa Krishtal
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, USA
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, USA
| |
Collapse
|
33
|
Aono S, Mori T, Sakaki S. 3D-RISM-MP2 Approach to Hydration Structure of Pt(II) and Pd(II) Complexes: Unusual H-Ahead Mode vs Usual O-Ahead One. J Chem Theory Comput 2016; 12:1189-206. [PMID: 26863511 DOI: 10.1021/acs.jctc.5b01137] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Solvation of transition metal complexes with water has been one of the fundamental topics in physical and coordination chemistry. In particular, Pt(II) complexes have recently attracted considerable interest for their relation to anticancer activity in cisplatin and its analogues, yet the interaction of the water molecule and the metal center has been obscured. The challenge from a theoretical perspective remains that both the microscopic solvation effect and the dynamical electron correlation (DEC) effect have to be treated simultaneously in a reasonable manner. In this work we derive the analytical gradient for the three-dimensional reference interaction site model Møller-Plesset second order (3D-RISM-MP2) free energy. On the basis of the three-regions 3D-RISM self-consistent field (SCF) method recently proposed by us, we apply a new layer of the Z-vector method to the CP-RISM equation as well as point-charge approximation to the derivatives with respect to the density matrix elements in the RISM-CPHF equation to remarkably reduce the computational cost. This method is applied to study the interaction of H2O with the d(8) square planar transition metal complexes in aqueous solution, trans-[Pt(II)Cl2(NH3)(glycine)] (1a), [Pt(II)(NH3)4](2+) (1b), [Pt(II)(CN)4](2-) (1c), and their Pd(II) analogues 2a, 2b, and 2c, respectively, to elucidate whether the usual H2O interaction through O atom (O-ahead mode) or unusual one through H atom (H-ahead mode) is stable in these complexes. We find that the interaction energy of the coordinating water and the transition metal complex changes little when switching from gas to aqueous phase, but the solvation free energy differs remarkably between the two interaction modes, thereby affecting the relative stability of the H-ahead and O-ahead modes. Particularly, in contrast to the expectation that the O-ahead mode is preferred due to the presence of positive charges in 1b, the H-ahead mode is also found to be more stable. The O-ahead mode is found to be more stable than the H-ahead one only in 2b. The energy decomposition analysis (EDA) at the 3D-RISM-MP2 level revealed that the O-ahead mode is stabilized by the electrostatic (ES) interaction, whereas the H-ahead one is mainly stabilized by the DEC effect. The ES interaction is also responsible for the difference between the Pd(II) and Pt(II) complexes; because the electrostatic potential is more negative along the z-axis in the Pt(II) complex than in the Pd(II) one, the O-ahead mode prefers the Pd(II) complexes, whereas the H-ahead becomes predominant in the Pt(II) complexes.
Collapse
Affiliation(s)
- Shinji Aono
- Fukui Institute for Fundamental Chemistry, Kyoto University , Nishihiraki-cho, Takano, Sakyo-ku, Kyoto 606-8103, Japan
| | - Toshifumi Mori
- Institute for Molecular Science , Okazaki, Aichi 444-8585, Japan.,School of Physical Sciences, The Graduate University for Advanced Studies , Okazaki, Aichi 444-8585, Japan
| | - Shigeyoshi Sakaki
- Fukui Institute for Fundamental Chemistry, Kyoto University , Nishihiraki-cho, Takano, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
34
|
Laktionov A, Chemineau-Chalaye E, Wesolowski TA. Frozen-density embedding theory with average solvent charge densities from explicit atomistic simulations. Phys Chem Chem Phys 2016; 18:21069-78. [DOI: 10.1039/c6cp00497k] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Besides molecular electron densities obtained within the Born–Oppenheimer approximation (ρB(r)) to represent the environment, the ensemble averaged density (〈ρB〉(r)) is also admissible in frozen-density embedding theory (FDET) [Wesolowski, Phys. Rev. A, 2008, 77, 11444].
Collapse
Affiliation(s)
- Andrey Laktionov
- Université de Genève
- Département de Chimie Physique 30
- CH-1211 Genève 4
- Switzerland
| | | | - Tomasz A. Wesolowski
- Université de Genève
- Département de Chimie Physique 30
- CH-1211 Genève 4
- Switzerland
| |
Collapse
|
35
|
Goez A, Neugebauer J. A Local Variant of the Conductor-Like Screening Model for Fragment-Based Electronic-Structure Methods. J Chem Theory Comput 2015; 11:5277-90. [DOI: 10.1021/acs.jctc.5b00832] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Albrecht Goez
- Theoretische Organische Chemie,
Organisch-Chemisches Institut and Center for Multiscale Theory and
Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße
40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie,
Organisch-Chemisches Institut and Center for Multiscale Theory and
Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße
40, 48149 Münster, Germany
| |
Collapse
|
36
|
Abstract
The calculation of electrostatic solute-solvent interactions in 3D RISM ("three-dimensional reference interaction site model") integral equation theory is recast in a form that allows for a computational treatment analogous to the "particle-mesh Ewald" formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems.
Collapse
Affiliation(s)
- Jochen Heil
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 6, 44227 Dortmund, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 6, 44227 Dortmund, Germany
| |
Collapse
|
37
|
Wesolowski TA, Shedge S, Zhou X. Frozen-Density Embedding Strategy for Multilevel Simulations of Electronic Structure. Chem Rev 2015; 115:5891-928. [DOI: 10.1021/cr500502v] [Citation(s) in RCA: 219] [Impact Index Per Article: 24.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Tomasz A. Wesolowski
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| | - Sapana Shedge
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| | - Xiuwen Zhou
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| |
Collapse
|
38
|
Śmiga S, Fabiano E, Laricchia S, Constantin LA, Della Sala F. Subsystem density functional theory with meta-generalized gradient approximation exchange-correlation functionals. J Chem Phys 2015; 142:154121. [DOI: 10.1063/1.4917257] [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)
- Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Torun, Poland
| | - Eduardo Fabiano
- Istituto Nanoscienze-CNR, Euromediterranean Center for Nanomaterial Modelling and Technology (ECMT), Via per Arnesano 16, I-73100 Lecce, Italy
- Center for Biomolecular Nanotechnologies@UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, LE, Italy
| | - Savio Laricchia
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Lucian A. Constantin
- Center for Biomolecular Nanotechnologies@UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, LE, Italy
| | - Fabio Della Sala
- Istituto Nanoscienze-CNR, Euromediterranean Center for Nanomaterial Modelling and Technology (ECMT), Via per Arnesano 16, I-73100 Lecce, Italy
- Center for Biomolecular Nanotechnologies@UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, LE, Italy
| |
Collapse
|
39
|
Huang W, Blinov N, Kovalenko A. Octanol-Water Partition Coefficient from 3D-RISM-KH Molecular Theory of Solvation with Partial Molar Volume Correction. J Phys Chem B 2015; 119:5588-97. [PMID: 25844645 DOI: 10.1021/acs.jpcb.5b01291] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The octanol-water partition coefficient is an important physical-chemical characteristic widely used to describe hydrophobic/hydrophilic properties of chemical compounds. The partition coefficient is related to the transfer free energy of a compound from water to octanol. Here, we introduce a new protocol for prediction of the partition coefficient based on the statistical-mechanical, 3D-RISM-KH molecular theory of solvation. It was shown recently that with the compound-solvent correlation functions obtained from the 3D-RISM-KH molecular theory of solvation, the free energy functional supplemented with the correction linearly related to the partial molar volume obtained from the Kirkwood-Buff/3D-RISM theory, also called the "universal correction" (UC), provides accurate prediction of the hydration free energy of small compounds, compared to explicit solvent molecular dynamics [ Palmer , D. S. ; J. Phys.: Condens. Matter 2010 , 22 , 492101 ]. Here we report that with the UC reparametrized accordingly this theory also provides an excellent agreement with the experimental data for the solvation free energy in nonpolar solvent (1-octanol) and so accurately predicts the octanol-water partition coefficient. The performance of the Kovalenko-Hirata (KH) and Gaussian fluctuation (GF) functionals of the solvation free energy, with and without UC, is tested on a large library of small compounds with diverse functional groups. The best agreement with the experimental data for octanol-water partition coefficients is obtained with the KH-UC solvation free energy functional.
Collapse
Affiliation(s)
- WenJuan Huang
- †Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta Canada.,‡National Institute for Nanotechnology, National Research Council of Canada, Edmonton, Alberta Canada
| | - Nikolay Blinov
- †Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta Canada.,‡National Institute for Nanotechnology, National Research Council of Canada, Edmonton, Alberta Canada
| | - Andriy Kovalenko
- †Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta Canada.,‡National Institute for Nanotechnology, National Research Council of Canada, Edmonton, Alberta Canada
| |
Collapse
|
40
|
Huang W, Blinov N, Wishart DS, Kovalenko A. Role of water in ligand binding to maltose-binding protein: insight from a new docking protocol based on the 3D-RISM-KH molecular theory of solvation. J Chem Inf Model 2015; 55:317-28. [PMID: 25545470 DOI: 10.1021/ci500520q] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Maltose-binding protein is a periplasmic binding protein responsible for transport of maltooligosaccarides through the periplasmic space of Gram-negative bacteria, as a part of the ABC transport system. The molecular mechanisms of the initial ligand binding and induced large scale motion of the protein's domains still remain elusive. In this study, we use a new docking protocol that combines a recently proposed explicit water placement algorithm based on the 3D-RISM-KH molecular theory of solvation and conventional docking software (AutoDock Vina) to explain the mechanisms of maltotriose binding to the apo-open state of a maltose-binding protein. We confirm the predictions of previous NMR spectroscopic experiments on binding modes of the ligand. We provide the molecular details on the binding mode that was not previously observed in the X-ray experiments. We show that this mode, which is defined by the fine balance between the protein-ligand direct interactions and solvation effects, can trigger the protein's domain motion resulting in the holo-closed structure of the maltose-binding protein with the maltotriose ligand in excellent agreement with the experimental data. We also discuss the role of water in blocking unfavorable binding sites and water-mediated interactions contributing to the stability of observable binding modes of maltotriose.
Collapse
Affiliation(s)
- WenJuan Huang
- Department of Mechanical Engineering, University of Alberta , Edmonton, AB T6G 2G8, Canada
| | | | | | | |
Collapse
|
41
|
Yesudas JP, Blinov N, Dew SK, Kovalenko A. Calculation of binding free energy of short double stranded oligonucleotides using MM/3D-RISM-KH approach. J Mol Liq 2015. [DOI: 10.1016/j.molliq.2014.11.017] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
42
|
Fabiano E, Laricchia S, Della Sala F. Frozen density embedding with non-integer subsystems' particle numbers. J Chem Phys 2014; 140:114101. [PMID: 24655166 DOI: 10.1063/1.4868033] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We extend the frozen density embedding theory to non-integer subsystems' particles numbers. Different features of this formulation are discussed, with special concern for approximate embedding calculations. In particular, we highlight the relation between the non-integer particle-number partition scheme and the resulting embedding errors. Finally, we provide a discussion of the implications of the present theory for the derivative discontinuity issue and the calculation of chemical reactivity descriptors.
Collapse
Affiliation(s)
- Eduardo Fabiano
- National Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via per Arnesano 16, I-73100 Lecce, Italy
| | - Savio Laricchia
- National Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via per Arnesano 16, I-73100 Lecce, Italy
| | - Fabio Della Sala
- National Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via per Arnesano 16, I-73100 Lecce, Italy
| |
Collapse
|
43
|
Shedge SV, Wesolowski TA. Nonuniform Continuum Model for Solvatochromism Based on Frozen-Density Embedding Theory. Chemphyschem 2014; 15:3291-300. [DOI: 10.1002/cphc.201402351] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2014] [Revised: 07/09/2014] [Indexed: 11/07/2022]
|
44
|
Omelyan I, Kovalenko A. Multiple time step molecular dynamics in the optimized isokinetic ensemble steered with the molecular theory of solvation: accelerating with advanced extrapolation of effective solvation forces. J Chem Phys 2014; 139:244106. [PMID: 24387356 DOI: 10.1063/1.4848716] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
We develop efficient handling of solvation forces in the multiscale method of multiple time step molecular dynamics (MTS-MD) of a biomolecule steered by the solvation free energy (effective solvation forces) obtained from the 3D-RISM-KH molecular theory of solvation (three-dimensional reference interaction site model complemented with the Kovalenko-Hirata closure approximation). To reduce the computational expenses, we calculate the effective solvation forces acting on the biomolecule by using advanced solvation force extrapolation (ASFE) at inner time steps while converging the 3D-RISM-KH integral equations only at large outer time steps. The idea of ASFE consists in developing a discrete non-Eckart rotational transformation of atomic coordinates that minimizes the distances between the atomic positions of the biomolecule at different time moments. The effective solvation forces for the biomolecule in a current conformation at an inner time step are then extrapolated in the transformed subspace of those at outer time steps by using a modified least square fit approach applied to a relatively small number of the best force-coordinate pairs. The latter are selected from an extended set collecting the effective solvation forces obtained from 3D-RISM-KH at outer time steps over a broad time interval. The MTS-MD integration with effective solvation forces obtained by converging 3D-RISM-KH at outer time steps and applying ASFE at inner time steps is stabilized by employing the optimized isokinetic Nosé-Hoover chain (OIN) ensemble. Compared to the previous extrapolation schemes used in combination with the Langevin thermostat, the ASFE approach substantially improves the accuracy of evaluation of effective solvation forces and in combination with the OIN thermostat enables a dramatic increase of outer time steps. We demonstrate on a fully flexible model of alanine dipeptide in aqueous solution that the MTS-MD/OIN/ASFE/3D-RISM-KH multiscale method of molecular dynamics steered by effective solvation forces allows huge outer time steps up to tens of picoseconds without affecting the equilibrium and conformational properties, and thus provides a 100- to 500-fold effective speedup in comparison to conventional MD with explicit solvent. With the statistical-mechanical 3D-RISM-KH account for effective solvation forces, the method provides efficient sampling of biomolecular processes with slow and/or rare solvation events such as conformational transitions of hydrated alanine dipeptide with the mean life times ranging from 30 ps up to 10 ns for "flip-flop" conformations, and is particularly beneficial for biomolecular systems with exchange and localization of solvent and ions, ligand binding, and molecular recognition.
Collapse
Affiliation(s)
- Igor Omelyan
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| | - Andriy Kovalenko
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| |
Collapse
|
45
|
Jacob CR, Neugebauer J. Subsystem density-functional theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014. [DOI: 10.1002/wcms.1175] [Citation(s) in RCA: 229] [Impact Index Per Article: 22.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Christoph R. Jacob
- Center for Functional Nanostructures and Institute of Physical Chemistry; Karlsruhe Institute of Technology (KIT); Karlsruhe Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut; Westfälische Wilhelms-Universität Münster; Münster Germany
| |
Collapse
|
46
|
|
47
|
State-specific solvation effect on the intramolecular charge transfer reaction in solution: A linear-response free energy TDDFT method. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.05.104] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
48
|
Minezawa N. Excited-state free energy surfaces in solution: time-dependent density functional theory∕reference interaction site model self-consistent field method. J Chem Phys 2014; 138:244101. [PMID: 23822221 DOI: 10.1063/1.4811201] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Constructing free energy surfaces for electronically excited states is a first step toward the understanding of photochemical processes in solution. For that purpose, the analytic free energy gradient is derived and implemented for the linear-response time-dependent density functional theory combined with the reference interaction site model self-consistent field method. The proposed method is applied to study (1) the fluorescence spectra of aqueous acetone and (2) the excited-state intramolecular proton transfer reaction of ortho-hydroxybenzaldehyde in an acetonitrile solution.
Collapse
Affiliation(s)
- Noriyuki Minezawa
- Fukui Institute for Fundamental Chemistry, Kyoto University, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
49
|
Innovations in Finite-Temperature Density Functionals. LECTURE NOTES IN COMPUTATIONAL SCIENCE AND ENGINEERING 2014. [DOI: 10.1007/978-3-319-04912-0_3] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
50
|
Joung IS, Luchko T, Case DA. Simple electrolyte solutions: comparison of DRISM and molecular dynamics results for alkali halide solutions. J Chem Phys 2013; 138:044103. [PMID: 23387564 DOI: 10.1063/1.4775743] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Using the dielectrically consistent reference interaction site model (DRISM) of molecular solvation, we have calculated structural and thermodynamic information of alkali-halide salts in aqueous solution, as a function of salt concentration. The impact of varying the closure relation used with DRISM is investigated using the partial series expansion of order-n (PSE-n) family of closures, which includes the commonly used hypernetted-chain equation (HNC) and Kovalenko-Hirata closures. Results are compared to explicit molecular dynamics (MD) simulations, using the same force fields, and to experiment. The mean activity coefficients of ions predicted by DRISM agree well with experimental values at concentrations below 0.5 m, especially when using the HNC closure. As individual ion activities (and the corresponding solvation free energies) are not known from experiment, only DRISM and MD results are directly compared and found to have reasonably good agreement. The activity of water directly estimated from DRISM is nearly consistent with values derived from the DRISM ion activities and the Gibbs-Duhem equation, but the changes in the computed pressure as a function of salt concentration dominate these comparisons. Good agreement with experiment is obtained if these pressure changes are ignored. Radial distribution functions of NaCl solution at three concentrations were compared between DRISM and MD simulations. DRISM shows comparable water distribution around the cation, but water structures around the anion deviate from the MD results; this may also be related to the high pressure of the system. Despite some problems, DRISM-PSE-n is an effective tool for investigating thermodynamic properties of simple electrolytes.
Collapse
Affiliation(s)
- In Suk Joung
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey 08854, USA
| | | | | |
Collapse
|