1
|
Folkestad S, Paul AC, Paul née Matveeva R, Reinholdt P, Coriani S, Odelius M, Koch H. Quantum Mechanical Versus Polarizable Embedding Schemes: A Study of the Xray Absorption Spectra of Aqueous Ammonia and Ammonium. J Chem Theory Comput 2024; 20:4161-4169. [PMID: 38713524 PMCID: PMC11137810 DOI: 10.1021/acs.jctc.4c00105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Revised: 04/16/2024] [Accepted: 04/17/2024] [Indexed: 05/09/2024]
Abstract
The X-ray absorption spectra of aqueous ammonia and ammonium are computed using a combination of coupled cluster singles and doubles (CCSD) with different quantum mechanical and molecular mechanical embedding schemes. Specifically, we compare frozen Hartree-Fock (HF) density embedding, polarizable embedding (PE), and polarizable density embedding (PDE). Integrating CCSD with frozen HF density embedding is possible within the CC-in-HF framework, which circumvents the conventional system-size limitations of standard coupled cluster methods. We reveal similarities between PDE and frozen HF density descriptions, while PE spectra differ significantly. By including approximate triple excitations, we also investigate the effect of improving the electronic structure theory. The spectra computed using this approach show an improved intensity ratio compared to CCSD-in-HF. Charge transfer analysis of the excitations shows the local character of the pre-edge and main-edge, while the post-edge is formed by excitations delocalized over the first solvation shell and beyond.
Collapse
Affiliation(s)
- Sarai
Dery Folkestad
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| | - Alexander C. Paul
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| | - Regina Paul née Matveeva
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| | - Peter Reinholdt
- Department
of Physics, Chemistry and Pharmacy, University
of Southern Denmark, SDU, Campusvej 55, 5230 Odense, Denmark
| | - Sonia Coriani
- Department
of Chemistry, Technical University of Denmark,
DTU, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Michael Odelius
- Department
of Physics, Stockholm University, 10691 Stockholm, Sweden
| | - Henrik Koch
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| |
Collapse
|
2
|
Weng G, Pang A, Vlček V. Spatial Decay and Limits of Quantum Solute-Solvent Interactions. J Phys Chem Lett 2023; 14:2473-2480. [PMID: 36867592 DOI: 10.1021/acs.jpclett.3c00208] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Molecular excitations in the liquid-phase environment are renormalized by the surrounding solvent molecules. Herein, we employ the GW approximation to investigate the solvation effects on the ionization energy of phenol in various solvent environments. The electronic effects differ by up to 0.4 eV among the five investigated solvents. This difference depends on both the macroscopic solvent polarizability and the spatial decay of the solvation effects. The latter is probed by separating the electronic subspace and the GW correlation self-energy into fragments. The fragment correlation energy decays with increasing intermolecular distance and vanishes at ∼9 Å, and this pattern is independent of the type of solvent environment. The 9 Å cutoff defines an effective interacting volume within which the ionization energy shift per solvent molecule is proportional to the macroscopic solvent polarizability. Finally, we propose a simple model for computing the ionization energies of molecules in an arbitrary solvent environment.
Collapse
Affiliation(s)
- Guorong Weng
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, United States
| | - Amanda Pang
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, United States
| | - Vojtěch Vlček
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, United States
- Materials Department, University of California, Santa Barbara, California 93106-5050, United States
| |
Collapse
|
3
|
Jónsson EÖ, Rasti S, Galynska M, Meyer J, Jónsson H. Transferable Potential Function for Flexible H 2O Molecules Based on the Single-Center Multipole Expansion. J Chem Theory Comput 2022; 18:7528-7543. [PMID: 36395502 DOI: 10.1021/acs.jctc.2c00598] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
A potential function is presented for describing a system of flexible H2O molecules based on the single-center multipole expansion (SCME) of the electrostatic interaction. The model, referred to as SCME/f, includes the variation of the molecular quadrupole moment as well as the dipole moment with changes in bond length and angle so as to reproduce results of high-level electronic structure calculations. The multipole expansion also includes fixed octupole and hexadecapole moments, as well as anisotropic dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole polarizability tensors. The model contains five adjustable parameters related to the repulsive interaction and damping functions in the electrostatic and dispersion interactions. Their values are adjusted to reproduce the lowest energy isomers of small clusters, (H2O)n with n = 2-6, as well as measured properties of the ice Ih crystal. Subsequent calculations of the energy difference between the various isomer configurations of the clusters show that SCME/f gives good agreement with results of electronic structure calculations and represents a significant improvement over the previously presented rigid SCME potential function. Analysis of the vibrational frequencies of the clusters and structural properties of ice Ih crystal show the importance of accurately describing the variation of the quadrupole moment with molecular structures.
Collapse
Affiliation(s)
- Elvar Örn Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107Reykjavík, Iceland
| | - Soroush Rasti
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RALeiden, The Netherlands
| | - Marta Galynska
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107Reykjavík, Iceland
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RALeiden, The Netherlands
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107Reykjavík, Iceland
| |
Collapse
|
4
|
Graham DS, Wen X, Chulhai DV, Goodpaster J. Huzinaga Projection Embedding for Efficient and Accurate Energies of Systems with Localized Spin-densities. J Chem Phys 2022; 156:054112. [DOI: 10.1063/5.0076493] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Xuelan Wen
- Department of Chemistry, University of Minnesota Twin Cities, United States of America
| | | | | |
Collapse
|
5
|
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: 5.3] [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
|
6
|
Osella S, Knippenberg S. The influence of lipid membranes on fluorescent probes' optical properties. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2020; 1863:183494. [PMID: 33129783 DOI: 10.1016/j.bbamem.2020.183494] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Subscribe] [Scholar Register] [Received: 06/20/2020] [Revised: 10/13/2020] [Accepted: 10/15/2020] [Indexed: 02/03/2023]
Abstract
BACKGROUND Organic fluorophores embedded in lipid bilayers can nowadays be described by a multiscale computational approach. Combining different length and time scales, a full characterization of the probe localization and optical properties led to novel insight into the effect of the environments. SCOPE OF REVIEW Following an introduction on computational advancements, three relevant probes are reviewed that delineate how a multiscale approach can lead to novel insight into the probes' (non) linear optical properties. Attention is paid to the quality of the theoretical description of the optical techniques. MAJOR CONCLUSIONS Computation can assess a priori novel probes' optical properties and guide the analysis and interpretation of experimental data in novel studies. The properties can be used to gain information on the phase and condition of the surrounding biological environment. GENERAL SIGNIFICANCE Computation showed that a canonical view on some of the probes should be revisited and adapted.
Collapse
Affiliation(s)
- Silvio Osella
- Biological Systems Simulation Lab, Centre of New Technologies, University of Warsaw, Banacha 2C, 02-097 Warsaw, Poland.
| | - Stefan Knippenberg
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic; Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-10691 Stockholm, Sweden; Theoretical Physics, Hasselt University, Agoralaan Building D, 3590 Diepenbeek, Belgium.
| |
Collapse
|
7
|
Reinholdt P, Jørgensen FK, Kongsted J, Olsen JMH. Polarizable Density Embedding for Large Biomolecular Systems. J Chem Theory Comput 2020; 16:5999-6006. [PMID: 32991163 DOI: 10.1021/acs.jctc.0c00763] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We present an efficient and robust fragment-based quantum-classical embedding model capable of accurately capturing effects from complex environments such as proteins and nucleic acids. This is realized by combining the molecular fractionation with conjugate caps (MFCC) procedure with the polarizable density embedding (PDE) model at the level of Fock matrix construction. The PDE contributions to the Fock matrix of the core region are constructed using the local molecular basis of the individual fragments rather than the supermolecular basis of the entire system. Thereby, we avoid complications associated with the application of the MFCC procedure on environment quantities such as electronic densities and molecular-orbital energies. Moreover, the computational cost associated with solving self-consistent field (SCF) equations of the core region remains unchanged from that of purely classical polarized embedding models. We analyze the performance of the resulting model in terms of the reproduction of the electrostatic potential of an insulin monomer protein and further in the context of solving problems related to electron spill-out. Finally, we showcase the model for the calculation of one- and two-photon properties of the Nile red molecule in a protein environment. Based on our analyses, we find that the combination of the MFCC approach with the PDE model is an efficient, yet accurate approach for calculating molecular properties of molecules embedded in structured biomolecular environments.
Collapse
Affiliation(s)
- Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Frederik Kamper Jørgensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Jógvan Magnus Haugaard Olsen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway.,Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
8
|
Giovannini T, Egidi F, Cappelli C. Molecular spectroscopy of aqueous solutions: a theoretical perspective. Chem Soc Rev 2020; 49:5664-5677. [PMID: 32744278 DOI: 10.1039/c9cs00464e] [Citation(s) in RCA: 62] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/28/2024]
Abstract
Computational spectroscopy is an invaluable tool to both accurately reproduce the spectra of molecular systems and provide a rationalization for the underlying physics. However, the inherent difficulty to accurately model systems in aqueous solutions, owing to water's high polarity and ability to form hydrogen bonds, has severely hampered the development of the field. In this tutorial review we present a technique developed and tested in recent years based on a fully atomistic and polarizable classical modeling of water coupled with a quantum mechanical description of the solute. Thanks to its unparalleled accuracy and versatility, this method can change the perspective of computational and experimental chemists alike.
Collapse
Affiliation(s)
| | - Franco Egidi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| |
Collapse
|
9
|
Olsen JMH, Reine S, Vahtras O, Kjellgren E, Reinholdt P, Hjorth Dundas KO, Li X, Cukras J, Ringholm M, Hedegård ED, Di Remigio R, List NH, Faber R, Cabral Tenorio BN, Bast R, Pedersen TB, Rinkevicius Z, Sauer SPA, Mikkelsen KV, Kongsted J, Coriani S, Ruud K, Helgaker T, Jensen HJA, Norman P. Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems. J Chem Phys 2020; 152:214115. [PMID: 32505165 DOI: 10.1063/1.5144298] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.
Collapse
Affiliation(s)
- Jógvan Magnus Haugaard Olsen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Simen Reine
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, N-0315 Oslo, Norway
| | - Olav Vahtras
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Erik Kjellgren
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Karen Oda Hjorth Dundas
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Xin Li
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Janusz Cukras
- Department of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
| | - Magnus Ringholm
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Erik D Hedegård
- Division of Theoretical Chemistry, Lund University, SE-223 62 Lund, Sweden
| | - Roberto Di Remigio
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Nanna H List
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Rasmus Faber
- DTU Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | | | - Radovan Bast
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Thomas Bondo Pedersen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, N-0315 Oslo, Norway
| | - Zilvinas Rinkevicius
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark
| | - Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Kenneth Ruud
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Trygve Helgaker
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, N-0315 Oslo, Norway
| | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Patrick Norman
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| |
Collapse
|
10
|
Graham DS, Wen X, Chulhai DV, Goodpaster JD. Robust, Accurate, and Efficient: Quantum Embedding Using the Huzinaga Level-Shift Projection Operator for Complex Systems. J Chem Theory Comput 2020; 16:2284-2295. [DOI: 10.1021/acs.jctc.9b01185] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Daniel S. Graham
- Department of Chemistry, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455, United States
| | - Xuelan Wen
- Department of Chemistry, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455, United States
| | - Dhabih V. Chulhai
- Department of Chemistry, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455, United States
| | - Jason D. Goodpaster
- Department of Chemistry, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
11
|
Bondanza M, Nottoli M, Cupellini L, Lipparini F, Mennucci B. Polarizable embedding QM/MM: the future gold standard for complex (bio)systems? Phys Chem Chem Phys 2020; 22:14433-14448. [DOI: 10.1039/d0cp02119a] [Citation(s) in RCA: 59] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.
Collapse
Affiliation(s)
- Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| |
Collapse
|
12
|
Dohn AO, Jónsson EÖ, Jónsson H. Polarizable Embedding with a Transferable H 2O Potential Function II: Application to (H 2O) n Clusters and Liquid Water. J Chem Theory Comput 2019; 15:6578-6587. [PMID: 31692344 DOI: 10.1021/acs.jctc.9b00778] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The incorporation of polarization in multiscale quantum-mechanics/molecular-mechanics (QM/MM) simulations is important for a variety of applications, for example, charge-transfer reactions. A recently developed formalism based on a density functional theory description of the QM region and a potential energy function for H2O molecules that includes quadrupole as well as dipole polarizability of the MM region is used to simulate liquid water and water clusters. Analysis of the energy, atomic forces, MM polarization, and structure is presented. A quantitative assessment of the QM/MM-MM/MM interaction energy differences of all possible QM/MM configurations of (H2O)n clusters shows that the interquartile range of the distributions of the QM/MM binding energies is never more than 20 meV/molecule higher or lower than the binding energies produced with either of the single-model results. Comparing these interaction energy differences with the QM/MM induction differences show that they are not systematically caused by the induced MM moments of our polarizable embedding scheme. Optimized hexamer geometries as well as the liquid water structure are shown to be improved in comparison with results obtained using point-charge based embedding models neglecting polarization.
Collapse
Affiliation(s)
- Asmus Ougaard Dohn
- Science Institute and Faculty of Physical Sciences , University of Iceland , Reykjavík 107 , Iceland
| | - Elvar Örn Jónsson
- Science Institute and Faculty of Physical Sciences , University of Iceland , Reykjavík 107 , Iceland
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences , University of Iceland , Reykjavík 107 , Iceland
| |
Collapse
|
13
|
Jónsson EÖ, Dohn AO, Jónsson H. Polarizable Embedding with a Transferable H 2O Potential Function I: Formulation and Tests on Dimer. J Chem Theory Comput 2019; 15:6562-6577. [PMID: 31689104 DOI: 10.1021/acs.jctc.9b00777] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The incorporation of mutual polarization in multiscale simulations where different regions of the system are treated at different level of theory is important in studies of, for example, electronic excitations and charge transfer processes. We present here an energy functional for describing a quantum mechanics/molecular mechanics (QM/MM) scheme that includes reciprocal polarization between the two subsystems. The inclusion of polarization alleviates shortcomings inherent in electrostatic embedding QM/MM models based on point-charge force fields. A density functional theory (DFT) description of the QM subsystem is coupled to a single center multipole expansion (SCME) description of H2O molecules in the MM subsystem that includes anisotropic dipole and quadrupole polarizability as well as static multipoles up to and including the hexadecapole. The energy functional and the coupling scheme is general and can be extended to arbitrary order in terms of both the static and induced moments. Tests of the energy surface for the H2O dimer show that the QM/MM results lie in between the pure DFT and pure SCME values. The consistency of the many-body contributions to the energy and analytical forces is demonstrated for an H2O pentamer.
Collapse
Affiliation(s)
- Elvar Örn Jónsson
- Science Institute and Faculty of Physical Sciences, VR-III , University of Iceland , Reykjavík 107 , Iceland
| | - Asmus Ougaard Dohn
- Science Institute and Faculty of Physical Sciences, VR-III , University of Iceland , Reykjavík 107 , Iceland
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, VR-III , University of Iceland , Reykjavík 107 , Iceland
| |
Collapse
|
14
|
Wen X, Graham DS, Chulhai DV, Goodpaster JD. Absolutely Localized Projection-Based Embedding for Excited States. J Chem Theory Comput 2019; 16:385-398. [DOI: 10.1021/acs.jctc.9b00959] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Xuelan Wen
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| | - Daniel S. Graham
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| | - Dhabih V. Chulhai
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| | - Jason D. Goodpaster
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
15
|
Izsák R. Single‐reference coupled cluster methods for computing excitation energies in large molecules: The efficiency and accuracy of approximations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1445] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Róbert Izsák
- Max‐Planck‐Institut für Kohlenforschung Mülheim an der Ruhr Germany
| |
Collapse
|