1
|
Medrano Sandonas L, Van Rompaey D, Fallani A, Hilfiker M, Hahn D, Perez-Benito L, Verhoeven J, Tresadern G, Kurt Wegner J, Ceulemans H, Tkatchenko A. Dataset for quantum-mechanical exploration of conformers and solvent effects in large drug-like molecules. Sci Data 2024; 11:742. [PMID: 38972891 PMCID: PMC11228031 DOI: 10.1038/s41597-024-03521-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2024] [Accepted: 06/13/2024] [Indexed: 07/09/2024] Open
Abstract
We here introduce the Aquamarine (AQM) dataset, an extensive quantum-mechanical (QM) dataset that contains the structural and electronic information of 59,783 low-and high-energy conformers of 1,653 molecules with a total number of atoms ranging from 2 to 92 (mean: 50.9), and containing up to 54 (mean: 28.2) non-hydrogen atoms. To gain insights into the solvent effects as well as collective dispersion interactions for drug-like molecules, we have performed QM calculations supplemented with a treatment of many-body dispersion (MBD) interactions of structures and properties in the gas phase and implicit water. Thus, AQM contains over 40 global and local physicochemical properties (including ground-state and response properties) per conformer computed at the tightly converged PBE0+MBD level of theory for gas-phase molecules, whereas PBE0+MBD with the modified Poisson-Boltzmann (MPB) model of water was used for solvated molecules. By addressing both molecule-solvent and dispersion interactions, AQM dataset can serve as a challenging benchmark for state-of-the-art machine learning methods for property modeling and de novo generation of large (solvated) molecules with pharmaceutical and biological relevance.
Collapse
Affiliation(s)
- Leonardo Medrano Sandonas
- Department of Physics and Materials Science, University of Luxembourg, L-1511, Luxembourg City, Luxembourg.
- Institute for Materials Science and Max Bergmann Center of Biomaterials, TU Dresden, 01062, Dresden, Germany.
| | - Dries Van Rompaey
- Drug Discovery Data Sciences (D3S), Janssen Pharmaceutica NV, Turnhoutseweg 30, 2340, Beerse, Belgium.
| | - Alessio Fallani
- Department of Physics and Materials Science, University of Luxembourg, L-1511, Luxembourg City, Luxembourg
- Drug Discovery Data Sciences (D3S), Janssen Pharmaceutica NV, Turnhoutseweg 30, 2340, Beerse, Belgium
| | - Mathias Hilfiker
- Department of Physics and Materials Science, University of Luxembourg, L-1511, Luxembourg City, Luxembourg
| | - David Hahn
- Computational Chemistry, Janssen Pharmaceutica NV, Turnhoutseweg 30, 2340, Beerse, Belgium
| | - Laura Perez-Benito
- Computational Chemistry, Janssen Pharmaceutica NV, Turnhoutseweg 30, 2340, Beerse, Belgium
| | - Jonas Verhoeven
- Drug Discovery Data Sciences (D3S), Janssen Pharmaceutica NV, Turnhoutseweg 30, 2340, Beerse, Belgium
| | - Gary Tresadern
- Computational Chemistry, Janssen Pharmaceutica NV, Turnhoutseweg 30, 2340, Beerse, Belgium
| | - Joerg Kurt Wegner
- Drug Discovery Data Sciences (D3S), Janssen Pharmaceutica NV, Turnhoutseweg 30, 2340, Beerse, Belgium
- Drug Discovery Data Sciences (D3S), Johnson & Johnson Innovative Medicine, 301 Binney Street, MA 02142, Cambridge, USA
| | - Hugo Ceulemans
- Drug Discovery Data Sciences (D3S), Janssen Pharmaceutica NV, Turnhoutseweg 30, 2340, Beerse, Belgium
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511, Luxembourg City, Luxembourg.
| |
Collapse
|
2
|
M V, Singh S, Bononi F, Andreussi O, Karmodak N. Thermodynamic and kinetic modeling of electrocatalytic reactions using a first-principles approach. J Chem Phys 2023; 159:111001. [PMID: 37728202 DOI: 10.1063/5.0165835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Accepted: 08/28/2023] [Indexed: 09/21/2023] Open
Abstract
The computational modeling of electrochemical interfaces and their applications in electrocatalysis has attracted great attention in recent years. While tremendous progress has been made in this area, however, the accurate atomistic descriptions at the electrode/electrolyte interfaces remain a great challenge. The Computational Hydrogen Electrode (CHE) method and continuum modeling of the solvent and electrolyte interactions form the basis for most of these methodological developments. Several posterior corrections have been added to the CHE method to improve its accuracy and widen its applications. The most recently developed grand canonical potential approaches with the embedded diffuse layer models have shown considerable improvement in defining interfacial interactions at electrode/electrolyte interfaces over the state-of-the-art computational models for electrocatalysis. In this Review, we present an overview of these different computational models developed over the years to quantitatively probe the thermodynamics and kinetics of electrochemical reactions in the presence of an electrified catalyst surface under various electrochemical environments. We begin our discussion by giving a brief picture of the different continuum solvation approaches, implemented within the ab initio method to effectively model the solvent and electrolyte interactions. Next, we present the thermodynamic and kinetic modeling approaches to determine the activity and stability of the electrocatalysts. A few applications to these approaches are also discussed. We conclude by giving an outlook on the different machine learning models that have been integrated with the thermodynamic approaches to improve their efficiency and widen their applicability.
Collapse
Affiliation(s)
- Vasanthapandiyan M
- Department of Chemistry, Shiv Nadar Institution of Eminence, Dadri, Gautam Buddha Nagar, Uttar Pradesh 201314, India
| | - Shagun Singh
- Department of Chemistry, Shiv Nadar Institution of Eminence, Dadri, Gautam Buddha Nagar, Uttar Pradesh 201314, India
| | - Fernanda Bononi
- Department of Physics, University of North Texas, Denton, Texas 76203, USA
| | - Oliviero Andreussi
- Department of Chemistry and Biochemistry, Boise State University, Boise, Idaho 83725, USA
| | - Naiwrit Karmodak
- Department of Chemistry, Shiv Nadar Institution of Eminence, Dadri, Gautam Buddha Nagar, Uttar Pradesh 201314, India
| |
Collapse
|
3
|
Sharma S, Zagalskaya A, Weitzner SE, Eggart L, Cho S, Hsu T, Chen X, Varley JB, Alexandrov V, Orme CA, Pham TA, Wood BC. Metal dissolution from first principles: Potential-dependent kinetics and charge transfer. Electrochim Acta 2022. [DOI: 10.1016/j.electacta.2022.141443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
4
|
Aziz A, Carrasco J. Modelling magnesium surfaces and their dissolution in an aqueous environment using an implicit solvent model.. J Chem Phys 2022; 156:174702. [DOI: 10.1063/5.0087683] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Magnesium has attracted a growing interest for its use in various applications, primarily due to its, abundance, lightweight properties and relatively low-cost. However, one major drawback to its widespread use remains its reactivity in aqueous environments, which is poorly understood at the atomistic level. Ab initio density functional theory methods are particularly well suited to bridge this knowledge gap, but the explicit simulation of electrified water/metal interfaces is often too costly from a computational viewpoint. Here we investigate water/Mg interfaces using the computationally efficient implicit solvent model VASPsol. We show that the Mg (0001), (10-10), and (10-11) surfaces each form different electrochemical double layers due to the anisotropic smoothing of the electron density at their surfaces, following Smoluchowski rules. We highlight the dependence that the position of the diffuse cavity surrounding the interface has on the potential of zero charge and the electron double layer capacitance, and how these parameters are also affected by the addition of explicated water and adsorbed OH. Lastly, we calculate the equilibrium potential of Mg2+ / Mg0 in an aqueous environment as 2.46 V vs. standard hydrogen electrode in excellent agreement with experiment.
Collapse
Affiliation(s)
| | - Javier Carrasco
- Power Storage: Batteries and Supercaps, CIC energiGUNE, Spain
| |
Collapse
|
5
|
Ringe S, Hörmann NG, Oberhofer H, Reuter K. Implicit Solvation Methods for Catalysis at Electrified Interfaces. Chem Rev 2021; 122:10777-10820. [PMID: 34928131 PMCID: PMC9227731 DOI: 10.1021/acs.chemrev.1c00675] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Implicit solvation
is an effective, highly coarse-grained approach
in atomic-scale simulations to account for a surrounding liquid electrolyte
on the level of a continuous polarizable medium. Originating in molecular
chemistry with finite solutes, implicit solvation techniques are now
increasingly used in the context of first-principles modeling of electrochemistry
and electrocatalysis at extended (often metallic) electrodes. The
prevalent ansatz to model the latter electrodes and the reactive surface
chemistry at them through slabs in periodic boundary condition supercells
brings its specific challenges. Foremost this concerns the difficulty
of describing the entire double layer forming at the electrified solid–liquid
interface (SLI) within supercell sizes tractable by commonly employed
density functional theory (DFT). We review liquid solvation methodology
from this specific application angle, highlighting in particular its
use in the widespread ab initio thermodynamics approach
to surface catalysis. Notably, implicit solvation can be employed
to mimic a polarization of the electrode’s electronic density
under the applied potential and the concomitant capacitive charging
of the entire double layer beyond the limitations of the employed
DFT supercell. Most critical for continuing advances of this effective
methodology for the SLI context is the lack of pertinent (experimental
or high-level theoretical) reference data needed for parametrization.
Collapse
Affiliation(s)
- Stefan Ringe
- Department of Energy Science and Engineering, Daegu Institute of Science and Technology (DGIST), Daegu 42988, Republic of Korea.,Energy Science & Engineering Research Center, Daegu Institute of Science and Technology (DGIST), Daegu 42988, Republic of Korea
| | - Nicolas G Hörmann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany.,Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany.,Chair for Theoretical Physics VII and Bavarian Center for Battery Technology (BayBatt), University of Bayreuth, Universitätsstraße 30, 95447 Bayreuth, Germany
| | - Karsten Reuter
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| |
Collapse
|
6
|
Hörmann NG, Reuter K. Thermodynamic cyclic voltammograms: peak positions and shapes. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:264004. [PMID: 33848987 DOI: 10.1088/1361-648x/abf7a1] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2021] [Accepted: 04/13/2021] [Indexed: 06/12/2023]
Abstract
Based on a mean-field description of thermodynamic cyclic voltammograms (CVs), we analyze here in full generality, how CV peak positions and shapes are related to the underlying interface energetics, in particular when also including electrostatic double layer (DL) effects. We show in particular, how non-Nernstian behaviour is related to capacitive DL charging, and how this relates to common adsorbate-centered interpretations such as a changed adsorption energetics due to dipole-field interactions and the electrosorption valency - the number of exchanged electrons upon electrosorption per adsorbate. Using Ag(111) in halide-containing solutions as test case, we demonstrate that DL effects can introduce peak shifts that are already explained by rationalizing the interaction of isolated adsorbates with the interfacial fields, while alterations of the peak shape are mainly driven by the coverage-dependence of the adsorbate dipoles. In addition, we analyze in detail how changing the experimental conditions such as the ion concentrations in the solvent but also of the background electrolyte can affect the CV peaks via their impact on the potential drop in the DL and the DL capacitance, respectively. These results suggest new routes to analyze experimental CVs and use of those for a detailed assessment of the accuracy of atomistic models of electrified interfaces e.g. with and without explicitly treated interfacial solvent and/or approximate implicit solvent models.
Collapse
Affiliation(s)
- Nicolas Georg Hörmann
- Theoretical Chemistry, Technische Universitaet Muenchen, Lichtenbergstraße 4, Garching, DE 85748, Germany
- Theory, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, Berlin, DE 14195, Germany
| | - Karsten Reuter
- Theory, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, Berlin, DE 14195, Germany
| |
Collapse
|
7
|
Hörmann NG, Reuter K. Thermodynamic Cyclic Voltammograms Based on Ab Initio Calculations: Ag(111) in Halide-Containing Solutions. J Chem Theory Comput 2021; 17:1782-1794. [PMID: 33606513 PMCID: PMC8023662 DOI: 10.1021/acs.jctc.0c01166] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Cyclic voltammograms
(CVs) are a central experimental tool for
assessing the structure and activity of electrochemical interfaces.
Based on a mean-field ansatz for the interface energetics under applied
potential conditions, we here derive an ab initio thermodynamics approach to efficiently simulate thermodynamic CVs.
All unknown parameters are determined from density functional theory
(DFT) calculations coupled to an implicit solvent model. For the showcased
CVs of Ag(111) electrodes in halide-anion-containing solutions, these
simulations demonstrate the relevance of double-layer contributions
to explain experimentally observed differences in peak shapes over
the halide series. Only the appropriate account of interfacial charging
allows us to capture the differences in equilibrium coverage and total
electronic surface charge that cause the varying peak shapes. As a
case in point, this analysis demonstrates that prominent features
in CVs do not only derive from changes in adsorbate structure or coverage
but can also be related to variations of the electrosorption valency.
Such double-layer effects are proportional to adsorbate-induced changes
in the work function and/or interfacial capacitance. They are thus
especially pronounced for electronegative halides and other adsorbates
that affect these interface properties. In addition, the analysis
allows us to draw conclusions on how the possible inaccuracy of implicit
solvation models can indirectly affect the accuracy of other predicted
quantities such as CVs.
Collapse
Affiliation(s)
- Nicolas G Hörmann
- Chair of Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, 85748 Garching, Germany.,Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Karsten Reuter
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
8
|
Abidi N, Lim KRG, Seh ZW, Steinmann SN. Atomistic modeling of electrocatalysis: Are we there yet? WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1499] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Nawras Abidi
- Univ Lyon, Ens de Lyon, CNRS UMR 5182 Université Claude Bernard Lyon 1, Laboratoire de Chimie, F69342, Lyon France
| | - Kang Rui Garrick Lim
- Institute of Materials Research and Engineering, Agency for Science, Technology and Research (A*STAR) Singapore
| | - Zhi Wei Seh
- Institute of Materials Research and Engineering, Agency for Science, Technology and Research (A*STAR) Singapore
| | - Stephan N. Steinmann
- Univ Lyon, Ens de Lyon, CNRS UMR 5182 Université Claude Bernard Lyon 1, Laboratoire de Chimie, F69342, Lyon France
| |
Collapse
|
9
|
Heenen HH, Gauthier JA, Kristoffersen HH, Ludwig T, Chan K. Solvation at metal/water interfaces: An ab initio molecular dynamics benchmark of common computational approaches. J Chem Phys 2020; 152:144703. [DOI: 10.1063/1.5144912] [Citation(s) in RCA: 63] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Hendrik H. Heenen
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Joseph A. Gauthier
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, USA
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | | | - Thomas Ludwig
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, USA
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Karen Chan
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| |
Collapse
|
10
|
Sakaushi K, Kumeda T, Hammes-Schiffer S, Melander MM, Sugino O. Advances and challenges for experiment and theory for multi-electron multi-proton transfer at electrified solid–liquid interfaces. Phys Chem Chem Phys 2020; 22:19401-19442. [DOI: 10.1039/d0cp02741c] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Understanding microscopic mechanism of multi-electron multi-proton transfer reactions at complexed systems is important for advancing electrochemistry-oriented science in the 21st century.
Collapse
Affiliation(s)
- Ken Sakaushi
- Center for Green Research on Energy and Environmental Materials
- National Institute for Materials Science
- Ibaraki 305-0044
- Japan
| | - Tomoaki Kumeda
- Center for Green Research on Energy and Environmental Materials
- National Institute for Materials Science
- Ibaraki 305-0044
- Japan
| | | | - Marko M. Melander
- Nanoscience Center
- Department of Chemistry
- University of Jyväskylä
- Jyväskylä
- Finland
| | - Osamu Sugino
- The Institute of Solid State Physics
- the University of Tokyo
- Chiba 277-8581
- Japan
| |
Collapse
|
11
|
Leung K. DFT modelling of explicit solid–solid interfaces in batteries: methods and challenges. Phys Chem Chem Phys 2020; 22:10412-10425. [DOI: 10.1039/c9cp06485k] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Density Functional Theory (DFT) calculations of electrode material properties in high energy density storage devices like lithium batteries have been standard practice for decades.
Collapse
Affiliation(s)
- Kevin Leung
- Sandia National Laboratories
- MS 1415
- Albuquerque
- USA
| |
Collapse
|
12
|
Stein CJ, Herbert JM, Head-Gordon M. The Poisson–Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients. J Chem Phys 2019; 151:224111. [DOI: 10.1063/1.5131020] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Affiliation(s)
- Christopher J. Stein
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
13
|
Gauthier JA, Dickens CF, Heenen HH, Vijay S, Ringe S, Chan K. Unified Approach to Implicit and Explicit Solvent Simulations of Electrochemical Reaction Energetics. J Chem Theory Comput 2019; 15:6895-6906. [DOI: 10.1021/acs.jctc.9b00717] [Citation(s) in RCA: 61] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Joseph A. Gauthier
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Colin F. Dickens
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Hendrik H. Heenen
- Department of Physics, Technical University of Denmark, DK-2800, Kongens Lyngby, Denmark
| | - Sudarshan Vijay
- Department of Physics, Technical University of Denmark, DK-2800, Kongens Lyngby, Denmark
| | - Stefan Ringe
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Karen Chan
- Department of Physics, Technical University of Denmark, DK-2800, Kongens Lyngby, Denmark
| |
Collapse
|
14
|
Zhdanov VP. Intracellular RNA delivery by lipid nanoparticles: Diffusion, degradation, and release. Biosystems 2019; 185:104032. [PMID: 31563119 DOI: 10.1016/j.biosystems.2019.104032] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2019] [Revised: 07/15/2019] [Indexed: 01/09/2023]
Abstract
Various RNAs (siRNAs, miRNAs, or mRNAs) can be delivered into cells by lipid nanoparticles (LNPs) of 50-150 nm in diameter. The subsequent RNA release from LNPs may occur via various scenarios. Herein, two related kinetic models are proposed. The first model takes into account that LNPs are often porous so that RNA molecules diffuse in and detach from nanopores. The analysis is focused on RNA diffusion from a pore. The analytical expression obtained for the RNA escape rate constant is used to identify the difference in the release of siRNAs, miRNAs, and mRNAs. The key message here is that the mRNA diffusion from pores appears to be too slow, and accordingly the mRNA release seems to occur primarily via degradation of LNPs. The second coarse-grained model describes the diffusion-mediated release of RNA from a LNP in the situation when this process is accompanied by the LNP degradation at the lipid-solution interface. The corresponding kinetics are shown in detail at different relative rates of the RNA diffusion and LNP degradation. Potentially, this can help to interpret drug plasma levels after various dosing regimens.
Collapse
Affiliation(s)
- Vladimir P Zhdanov
- Section of Biological Physics, Department of Physics, Chalmers University of Technology, Göteborg, Sweden; Boreskov Institute of Catalysis, Russian Academy of Sciences, Novosibirsk, Russia.
| |
Collapse
|
15
|
Le Bahers T, Takanabe K. Combined theoretical and experimental characterizations of semiconductors for photoelectrocatalytic applications. JOURNAL OF PHOTOCHEMISTRY AND PHOTOBIOLOGY C-PHOTOCHEMISTRY REVIEWS 2019. [DOI: 10.1016/j.jphotochemrev.2019.01.001] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
|
16
|
Affiliation(s)
- Matthew Truscott
- Department of Physics, University of North Texas, Denton, Texas 76207, United States
| | - Oliviero Andreussi
- Department of Physics, University of North Texas, Denton, Texas 76207, United States
| |
Collapse
|
17
|
Nattino F, Truscott M, Marzari N, Andreussi O. Continuum models of the electrochemical diffuse layer in electronic-structure calculations. J Chem Phys 2019; 150:041722. [PMID: 30709273 DOI: 10.1063/1.5054588] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Continuum electrolyte models represent a practical tool to account for the presence of the diffuse layer at electrochemical interfaces. However, despite the increasing popularity of these in the field of materials science, it remains unclear which features are necessary in order to accurately describe interface-related observables such as the differential capacitance (DC) of metal electrode surfaces. We present here a critical comparison of continuum diffuse-layer models that can be coupled to an atomistic first-principles description of the charged metal surface in order to account for the electrolyte screening at electrified interfaces. By comparing computed DC values for the prototypical Ag(100) surface in an aqueous solution to experimental data, we validate the accuracy of the models considered. Results suggest that a size-modified Poisson-Boltzmann description of the electrolyte solution is sufficient to qualitatively reproduce the main experimental trends. Our findings also highlight the large effect that the dielectric cavity parameterization has on the computed DC values.
Collapse
Affiliation(s)
- Francesco Nattino
- Theory and Simulations of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Matthew Truscott
- Department of Physics, University of North Texas, Denton, Texas 76207, USA
| | - Nicola Marzari
- Theory and Simulations of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Oliviero Andreussi
- Department of Physics, University of North Texas, Denton, Texas 76207, USA
| |
Collapse
|
18
|
Hille C, Ringe S, Deimel M, Kunkel C, Acree WE, Reuter K, Oberhofer H. Generalized molecular solvation in non-aqueous solutions by a single parameter implicit solvation scheme. J Chem Phys 2019; 150:041710. [PMID: 30709294 DOI: 10.1063/1.5050938] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
In computer simulations of solvation effects on chemical reactions, continuum modeling techniques regain popularity as a way to efficiently circumvent an otherwise costly sampling of solvent degrees of freedom. As effective techniques, such implicit solvation models always depend on a number of parameters that need to be determined earlier. In the past, the focus lay mostly on an accurate parametrization of water models. Yet, non-aqueous solvents have recently attracted increasing attention, in particular, for the design of battery materials. To this end, we present a systematic parametrization protocol for the Self-Consistent Continuum Solvation (SCCS) model resulting in optimized parameters for 67 non-aqueous solvents. Our parametrization is based on a collection of ≈6000 experimentally measured partition coefficients, which we collected in the Solv@TUM database presented here. The accuracy of our optimized SCCS model is comparable to the well-known universal continuum solvation model (SMx) family of methods, while relying on only a single fit parameter and thereby largely reducing statistical noise. Furthermore, slightly modifying the non-electrostatic terms of the model, we present the SCCS-P solvation model as a more accurate alternative, in particular, for aromatic solutes. Finally, we show that SCCS parameters can, to a good degree of accuracy, also be predicted for solvents outside the database using merely the dielectric bulk permittivity of the solvent of choice.
Collapse
Affiliation(s)
- Christoph Hille
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, 85747 Garching, Germany
| | - Stefan Ringe
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, USA
| | - Martin Deimel
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, 85747 Garching, Germany
| | - Christian Kunkel
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, 85747 Garching, Germany
| | - William E Acree
- Department of Chemistry, University of North Texas, 1155 Union Circle Drive #305070, Denton, Texas 76203, USA
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, 85747 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, 85747 Garching, Germany
| |
Collapse
|
19
|
Melander MM, Kuisma MJ, Christensen TEK, Honkala K. Grand-canonical approach to density functional theory of electrocatalytic systems: Thermodynamics of solid-liquid interfaces at constant ion and electrode potentials. J Chem Phys 2019; 150:041706. [PMID: 30709274 DOI: 10.1063/1.5047829] [Citation(s) in RCA: 77] [Impact Index Per Article: 15.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Properties of solid-liquid interfaces are of immense importance for electrocatalytic and electrochemical systems, but modeling such interfaces at the atomic level presents a serious challenge and approaches beyond standard methodologies are needed. An atomistic computational scheme needs to treat at least part of the system quantum mechanically to describe adsorption and reactions, while the entire system is in thermal equilibrium. The experimentally relevant macroscopic control variables are temperature, electrode potential, and the choice of the solvent and ions, and these need to be explicitly included in the computational model as well; this calls for a thermodynamic ensemble with fixed ion and electrode potentials. In this work, a general framework within density functional theory (DFT) with fixed electron and ion chemical potentials in the grand canonical (GC) ensemble is established for modeling electrocatalytic and electrochemical interfaces. Starting from a fully quantum mechanical description of multi-component GC-DFT for nuclei and electrons, a systematic coarse-graining is employed to establish various computational schemes including (i) the combination of classical and electronic DFTs within the GC ensemble and (ii) on the simplest level a chemically and physically sound way to obtain various (modified) Poisson-Boltzmann (mPB) implicit solvent models. The detailed and rigorous derivation clearly establishes which approximations are needed for coarse-graining as well as highlights which details and interactions are omitted in vein of computational feasibility. The transparent approximations also allow removing some of the constraints and coarse-graining if needed. We implement various mPB models within a linear dielectric continuum in the GPAW code and test their capabilities to model capacitance of electrochemical interfaces as well as study different approaches for modeling partly periodic charged systems. Our rigorous and well-defined DFT coarse-graining scheme to continuum electrolytes highlights the inadequacy of current linear dielectric models for treating properties of the electrochemical interface.
Collapse
Affiliation(s)
- Marko M Melander
- Nanoscience Center, Department of Chemistry, University of Jyväskylä, P.O. Box 35 (YN), FI-40014 Jyväskylä, Finland
| | - Mikael J Kuisma
- Nanoscience Center, Department of Chemistry, University of Jyväskylä, P.O. Box 35 (YN), FI-40014 Jyväskylä, Finland
| | | | - Karoliina Honkala
- Nanoscience Center, Department of Chemistry, University of Jyväskylä, P.O. Box 35 (YN), FI-40014 Jyväskylä, Finland
| |
Collapse
|
20
|
Hörmann NG, Andreussi O, Marzari N. Grand canonical simulations of electrochemical interfaces in implicit solvation models. J Chem Phys 2019; 150:041730. [DOI: 10.1063/1.5054580] [Citation(s) in RCA: 87] [Impact Index Per Article: 17.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Nicolas G. Hörmann
- Theory and Simulation of Materials (THEOS), and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Oliviero Andreussi
- Department of Physics, University of North Texas, Denton, Texas 76207, USA
| | - Nicola Marzari
- Theory and Simulation of Materials (THEOS), and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
21
|
Gauthier JA, Ringe S, Dickens CF, Garza AJ, Bell AT, Head-Gordon M, Nørskov JK, Chan K. Challenges in Modeling Electrochemical Reaction Energetics with Polarizable Continuum Models. ACS Catal 2018. [DOI: 10.1021/acscatal.8b02793] [Citation(s) in RCA: 114] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Joseph A. Gauthier
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Stefan Ringe
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Colin F. Dickens
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Alejandro J. Garza
- The Joint Center for Artificial Photosynthesis, Lawrence Berkeley National Laboratory, Berkeley California 94720, United States
| | - Alexis T. Bell
- The Joint Center for Artificial Photosynthesis, Lawrence Berkeley National Laboratory, Berkeley California 94720, United States
- Department of Chemical and Biomolecular Engineering, University of California at Berkeley, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- The Joint Center for Artificial Photosynthesis, Lawrence Berkeley National Laboratory, Berkeley California 94720, United States
- Department of Chemistry, University of California at Berkeley, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Jens K. Nørskov
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Karen Chan
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
22
|
Womack JC, Anton L, Dziedzic J, Hasnip PJ, Probert MIJ, Skylaris CK. DL_MG: A Parallel Multigrid Poisson and Poisson-Boltzmann Solver for Electronic Structure Calculations in Vacuum and Solution. J Chem Theory Comput 2018; 14:1412-1432. [PMID: 29447447 DOI: 10.1021/acs.jctc.7b01274] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The solution of the Poisson equation is a crucial step in electronic structure calculations, yielding the electrostatic potential-a key component of the quantum mechanical Hamiltonian. In recent decades, theoretical advances and increases in computer performance have made it possible to simulate the electronic structure of extended systems in complex environments. This requires the solution of more complicated variants of the Poisson equation, featuring nonhomogeneous dielectric permittivities, ionic concentrations with nonlinear dependencies, and diverse boundary conditions. The analytic solutions generally used to solve the Poisson equation in vacuum (or with homogeneous permittivity) are not applicable in these circumstances, and numerical methods must be used. In this work, we present DL_MG, a flexible, scalable, and accurate solver library, developed specifically to tackle the challenges of solving the Poisson equation in modern large-scale electronic structure calculations on parallel computers. Our solver is based on the multigrid approach and uses an iterative high-order defect correction method to improve the accuracy of solutions. Using two chemically relevant model systems, we tested the accuracy and computational performance of DL_MG when solving the generalized Poisson and Poisson-Boltzmann equations, demonstrating excellent agreement with analytic solutions and efficient scaling to ∼109 unknowns and 100s of CPU cores. We also applied DL_MG in actual large-scale electronic structure calculations, using the ONETEP linear-scaling electronic structure package to study a 2615 atom protein-ligand complex with routinely available computational resources. In these calculations, the overall execution time with DL_MG was not significantly greater than the time required for calculations using a conventional FFT-based solver.
Collapse
Affiliation(s)
- James C Womack
- Department of Chemistry , University of Southampton , Highfield, Southampton SO17 1BJ , United Kingdom
| | - Lucian Anton
- Cray U.K. Ltd. , Broad Quay House, Prince Street , Bristol BS1 4DJ , United Kingdom
| | - Jacek Dziedzic
- Department of Chemistry , University of Southampton , Highfield, Southampton SO17 1BJ , United Kingdom.,Faculty of Applied Physics and Mathematics , Gdańsk University of Technology , Gdańsk 80-233 , Poland
| | - Phil J Hasnip
- Department of Physics , University of York , Heslington, York YO10 5DD , United Kingdom
| | - Matt I J Probert
- Department of Physics , University of York , Heslington, York YO10 5DD , United Kingdom
| | - Chris-Kriton Skylaris
- Department of Chemistry , University of Southampton , Highfield, Southampton SO17 1BJ , United Kingdom
| |
Collapse
|
23
|
Liu Z, Timmermann J, Reuter K, Scheurer C. Benchmarks and Dielectric Constants for Reparametrized OPLS and Polarizable Force Field Models of Chlorinated Hydrocarbons. J Phys Chem B 2017; 122:770-779. [PMID: 28832148 DOI: 10.1021/acs.jpcb.7b06709] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The knowledge of dielectric response properties of the environment is of paramount importance in many theoretical embedding methods and studies of solutes and of catalytic sites and processes in condensed phases. In particular, the realistic embedding of active sites into solid/liquid and liquid/liquid interfaces is a crucial point in the context of modeling energy conversion (e.g., electrochemical, photochemical, power-to-X) processes. Recently, the finding that the dielectric permeability of liquids near solid/liquid interfaces is far from being constant but deviates strongly from the bulk value within several nanometers from the interface has raised the interest in a more fundamental understanding of the response properties near interfaces. As these questions are hard to study experimentally, reliable theoretical models are required. Here we describe a careful first-principles based reparametrization of nonpolarizable molecular mechanics force fields for a class of technological relevant organic chlorinated hydrocarbon solvents which are immiscible with water. For the solvent 1,2-dichloroethane (1,2-DCE) we also present a new polarizable force field based on the Drude oscillator model. Its parametrization needs particular attention to avoid unphysical couplings between the internal torsional degree of freedom and the Drude oscillators, which could severely skew the response properties. The performance of this new set of force fields is critically assessed based on a comprehensive molecular dynamics study.
Collapse
Affiliation(s)
- Zhu Liu
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstrasse 4, 85747 Garching, Germany
| | - Jakob Timmermann
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstrasse 4, 85747 Garching, Germany
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstrasse 4, 85747 Garching, Germany
| | - Christoph Scheurer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstrasse 4, 85747 Garching, Germany
| |
Collapse
|