1
|
Wang Y, Teng C, Begin E, Bussiere M, Bao JL. PW-SMD: A Plane-Wave Implicit Solvation Model Based on Electron Density for Surface Chemistry and Crystalline Systems in Aqueous Solution. J Chem Theory Comput 2024. [PMID: 39024317 DOI: 10.1021/acs.jctc.4c00594] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
Electron density-based implicit solvation models are a class of techniques for quantifying solvation effects and calculating free energies of solvation without an explicit representation of solvent molecules. Integral to the accuracy of solvation modeling is the proper definition of the solvation shell separating the solute molecule from the solvent environment, allowing for a physical partitioning of the free energies of solvation. Unlike state-of-the-art implicit solvation models for molecular quantum chemistry calculations, e.g., the solvation model based on solute electron density (SMD), solvation models for systems under periodic boundary conditions with plane-wave (PW) basis sets have been limited in their accuracy. Furthermore, a unified implicit solvation model with both homogeneous solution-phase and heterogeneous interfacial structures treated on equal footing is needed. In order to address this challenge, we developed a high-accuracy solvation model for periodic PW calculations that is applicable to molecular, ionic, interfacial, and bulk-phase chemistry. Our model, PW-SMD, is an extension of the SMD molecular solvation model to periodic systems in water. The free energy of solvation is partitioned into the electrostatic and cavity-dispersion-solvent structure (CDS) contributions. The electrostatic contributions of the solvation shell surrounding solute structures are parametrized based on their geometric and physical properties. In addition, the nonelectrostatic contribution to the solvation energy is accounted for by extending the CDS formalism of SMD to incorporate periodic boundary conditions. We validate the accuracy and robustness of our solvation model by comparing predicted solvation free energies against experimental data for molecular and ionic systems, carved-cluster composite energetic models of solvated reaction energies and barriers on surface systems, and deep-learning-accelerated ab initio molecular dynamics (AIMD). Our developed periodic implicit solvation model shows significantly improved accuracy compared to previous work (namely, solvation models in aqueous solution) and can be applied to simulate solvent effects in a wide range of surface and crystalline materials.
Collapse
Affiliation(s)
- Yang Wang
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Chong Teng
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Elijah Begin
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Mason Bussiere
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| |
Collapse
|
2
|
Labat M, Giner E, Jeanmairet G. Coupling molecular density functional theory with converged selected configuration interaction methods to study excited states in aqueous solution. J Chem Phys 2024; 161:014113. [PMID: 38958166 DOI: 10.1063/5.0213426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Accepted: 06/12/2024] [Indexed: 07/04/2024] Open
Abstract
This paper presents the first implementation of a coupling between advanced wavefunction theories and molecular density functional theory (MDFT). This method enables the modeling of solvent effect into quantum mechanical (QM) calculations by incorporating an electrostatic potential generated by solvent charges into the electronic Hamiltonian. Solvent charges are deduced from the spatially and angularly dependent solvent particle density. Such a density is obtained through the minimization of the functional associated with the molecular mechanics (MM) Hamiltonian describing the interaction between the fluid particles. The introduced QM/MDFT framework belongs to QM/MM family of methods, but its originality lies in the use of MDFT as the MM solver, offering two main advantages. First, its functional formulation makes it competitive with respect to sampling-based molecular mechanics. Second, it preserves a molecular-level description lost in macroscopic continuum approaches. The excited state properties of water and formaldehyde molecules solvated into water have been computed at the selected configuration interaction (SCI) level. The excitation energies and dipole moments have been compared with experimental data and previous theoretical work. A key finding is that using the Hartree-Fock method to describe the solute allows for predicting the solvent charge around the ground state with sufficient precision for the subsequent SCI calculations of excited states. This significantly reduces the computational cost of the described procedure, paving the way for the study of more complex molecules.
Collapse
Affiliation(s)
- Maxime Labat
- Sorbonne Université, CNRS, Physico-Chimie des électrolytes et Nanosystèmes Interfaciaux, PHENIX, F-75005 Paris, France
| | - Emmanuel Giner
- Sorbonne Université, CNRS, Laboratoire de Chimie Théorique, Sorbonne Université, F-75005 Paris, France
| | - 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
| |
Collapse
|
3
|
Si P, Jayanth A, Andreussi O. Soft-sphere continuum solvation models for nonaqueous solvents. J Comput Chem 2024; 45:719-737. [PMID: 38112395 DOI: 10.1002/jcc.27254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Revised: 10/24/2023] [Accepted: 10/30/2023] [Indexed: 12/21/2023]
Abstract
Solvation effects profoundly influence the characteristics and behavior of chemical systems in liquid solutions. The interaction between solute and solvent molecules intricately impacts solubility, reactivity, stability, and various chemical processes. Continuum solvation models gained prominence in quantum chemistry by implicitly capturing these interactions and enabling efficient investigations of diverse chemical systems in solution. In comparison, continuum solvation models in condensed matter simulation are very recent. Among these, the self-consistent continuum solvation (SCCS) and the soft-sphere continuum solvation models (SSCS) have been among the first to be successfully parameterized and extended to model periodic systems in aqueous solutions and electrolytes. As most continuum approaches, these models depend on a number of parameters that are linked to experimental or theoretical properties of the solvent, or that can be tuned based on reference data. Here, we present a systematic parameterization of the SSCS model for over 100 nonaqueous solvents. We validate the model's efficacy across diverse solvent environments by leveraging experimental solvation-free energies and partition coefficients from comprehensive databases. The average root means square error over all the solvents was calculated as 0.85 kcal/mol which is below the chemical accuracy (1 kcal/mol). Similarly to what has been reported by Hille et al. (J. Chem. Phys. 2019, 150, 041710.) for the SCCS model, a single-parameter model accurately reproduces experimental solvation energies, showcasing the transferability and predictive power of these continuum approaches. Our findings underscore the potential for a unified approach to predict solvation properties, paving the way for enhanced computational studies across various chemical environments.
Collapse
Affiliation(s)
- Pradip Si
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| | - Ajay Jayanth
- Texas Academy of Math and Science, University of North Texas, Denton, Texas, USA
| | - Oliviero Andreussi
- Department of Chemistry and Biochemistry, Boise State University, Boise, Idaho, USA
| |
Collapse
|
4
|
Luan B, McDonagh JL. Developing semi-empirical water model for efficiently simulating temperature-dependent chemisorption of CO 2 in amine solvents. Phys Chem Chem Phys 2024; 26:3540-3547. [PMID: 38214052 DOI: 10.1039/d3cp05874c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
Classical molecular dynamics (MD) simulations without bond forming/breaking cannot be used to model chemical reactions (CRs) among small molecules. Although the first-principle MD simulation can adequately describe CRs with explicit water molecules, such simulation is normally too costly for most researchers to afford. Generally, water molecules in a solvent can exert hydrophobic forces on reacting molecules, which yields a so-called caging effect that cannot be ignored when constructing a free energy landscape for reacting molecules. Many recently developed semi-empirical methods (such as DFTB, PM6 and xTB) are highly efficient for modeling CRs, however none of them can be directly used to model bulk water properly. Here, we developed a modified xTB approach that enables the simulation of CRs in explicit water. Using the chemisorption of CO2 by amines in water as an example application, we demonstrate that our approach yielded results comparable with the first-principle ones, while only using a limited computing resource. Potentially, our proposed semi-empirical water model can be utilized for the computational study of any CR in water.
Collapse
Affiliation(s)
- Binquan Luan
- IBM Thomas J. Watson Research, Yorktown Heights, NY 10598, USA.
| | - James L McDonagh
- IBM Research Europe, Hartree Centre, SciTech Daresbury, Warrington, Chesire WA4 4AD, UK
| |
Collapse
|
5
|
Eggert T, Hörmann NG, Reuter K. Cavity formation at metal-water interfaces. J Chem Phys 2023; 159:194702. [PMID: 37966001 DOI: 10.1063/5.0167406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2023] [Accepted: 10/26/2023] [Indexed: 11/16/2023] Open
Abstract
The free energy cost of forming a cavity in a solvent is a fundamental concept in rationalizing the solvation of molecules and ions. A detailed understanding of the factors governing cavity formation in bulk solutions has inter alia enabled the formulation of models that account for this contribution in coarse-grained implicit solvation methods. Here, we employ classical molecular dynamics simulations and multistate Bennett acceptance ratio free energy sampling to systematically study cavity formation at a wide range of metal-water interfaces. We demonstrate that the obtained size- and position-dependence of cavitation energies can be fully rationalized by a geometric Gibbs model, which considers that the creation of the metal-cavity interface necessarily involves the removal of interfacial solvent. This so-called competitive adsorption effect introduces a substrate dependence to the interfacial cavity formation energy that is missed in existing bulk cavitation models. Using expressions from scaled particle theory, this substrate dependence is quantitatively reproduced by the Gibbs model through simple linear relations with the adsorption energy of a single water molecule. Besides providing a better general understanding of interfacial solvation, this paves the way for the derivation and efficient parametrization of more accurate interface-aware implicit solvation models needed for reliable high-throughput calculations toward improved electrocatalysts.
Collapse
Affiliation(s)
- Thorben Eggert
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
- Chair of Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, 85748 Garching, Germany
| | - Nicolas G Hörmann
- 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
|
6
|
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
|
7
|
Unraveling the Nature of Hydrogen Bonds of "Proton Sponges" Based on Car-Parrinello and Metadynamics Approaches. Int J Mol Sci 2023; 24:ijms24021542. [PMID: 36675059 PMCID: PMC9860969 DOI: 10.3390/ijms24021542] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2022] [Revised: 12/26/2022] [Accepted: 01/05/2023] [Indexed: 01/14/2023] Open
Abstract
The nature of intra- and intermolecular non-covalent interactions was studied in four naphthalene derivatives commonly referred to as "proton sponges". Special attention was paid to an intramolecular hydrogen bond present in the protonated form of the compounds. The unsubstituted "proton sponge" served as a reference structure to study the substituent influence on the hydrogen bond (HB) properties. We selected three compounds substituted by methoxy, amino, and nitro groups. The presence of the substituents either retained the parent symmetry or rendered the compounds asymmetric. In order to reveal the non-covalent interaction properties, the Hirshfeld surface (HS) was computed for the crystal structures of the studied compounds. Next, quantum-chemical simulations were performed in vacuo and in the crystalline phase. Car-Parrinello molecular dynamics (CPMD), Path Integral Molecular Dynamics (PIMD), and metadynamics were employed to investigate the time-evolution changes of metric parameters and free energy profile in both phases. Additionally, for selected snapshots obtained from the CPMD trajectories, non-covalent interactions and electronic structure were studied. Quantum theory of atoms in molecules (QTAIM) and the Density Overlap Regions Indicator (DORI) were applied for this purpose. It was found based on Hirshfeld surfaces that, besides intramolecular hydrogen bonds, other non-covalent interactions are present and have a strong impact on the crystal structure organization. The CPMD results obtained in both phases showed frequent proton transfer phenomena. The proton was strongly delocalized in the applied time-scale and temperature, especially in the PIMD framework. The use of metadynamics allowed for tracing the free energy profiles and confirming that the hydrogen bonds present in "proton sponges" are Low-Barrier Hydrogen Bonds (LBHBs). The electronic and topological analysis quantitatively described the temperature dependence and time-evolution changes of the electronic structure. The covalency of the hydrogen bonds was estimated based on QTAIM analysis. It was found that strong hydrogen bonds show greater covalency, which is additionally determined by the proton position in the hydrogen bridge.
Collapse
|
8
|
Prentice JCA. Efficiently Computing Excitations of Complex Systems: Linear-Scaling Time-Dependent Embedded Mean-Field Theory in Implicit Solvent. J Chem Theory Comput 2022; 18:1542-1554. [PMID: 35133827 PMCID: PMC9082505 DOI: 10.1021/acs.jctc.1c01133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Quantum embedding schemes have the
potential to significantly reduce
the computational cost of first-principles calculations while maintaining
accuracy, particularly for calculations of electronic excitations
in complex systems. In this work, I combine time-dependent embedded
mean field theory (TD-EMFT) with linear-scaling density functional
theory and implicit solvation models, extending previous work within
the ONETEP code. This provides a way to perform multilevel calculations
of electronic excitations on very large systems, where long-range
environmental effects, both quantum and classical in nature, are important.
I demonstrate the power of this method by performing simulations on
a variety of systems, including a molecular dimer, a chromophore in
solution, and a doped molecular crystal. This work paves the way for
high accuracy calculations to be performed on large-scale systems
that were previously beyond the reach of quantum embedding schemes.
Collapse
Affiliation(s)
- Joseph C A Prentice
- Department of Materials, University of Oxford, Parks Road, Oxford OX1 3PH, United Kingdom
| |
Collapse
|
9
|
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
|
10
|
Oğuz IC, Vassetti D, Labat F. Assessing the performances of different continuum solvation models for the calculation of hydration energies of molecules, polymers and surfaces: a comparison between the SMD, VASPsol and FDPB models. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02799-w] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
11
|
Yang X, Zhuang Y, Zhu J, Le J, Cheng J. Recent progress on multiscale modeling of electrochemistry. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1559] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Xiao‐Hui Yang
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Yong‐Bin Zhuang
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Jia‐Xin Zhu
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Jia‐Bo Le
- Ningbo Institute of Materials Technology and Engineering Chinese Academy of Sciences Ningbo China
| | - Jun Cheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| |
Collapse
|
12
|
Hu Y, Luo P. Energy change mechanisms of HMX solute molecules in pure solvents and binary solvent mixtures. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.115898] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
13
|
Herbert JM. Dielectric continuum methods for quantum chemistry. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1519] [Citation(s) in RCA: 37] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- John M. Herbert
- Department of Chemistry and Biochemistry The Ohio State University Columbus Ohio USA
| |
Collapse
|
14
|
Mao Y, Loipersberger M, Kron KJ, Derrick JS, Chang CJ, Sharada SM, Head-Gordon M. Consistent inclusion of continuum solvation in energy decomposition analysis: theory and application to molecular CO 2 reduction catalysts. Chem Sci 2020; 12:1398-1414. [PMID: 34163903 PMCID: PMC8179122 DOI: 10.1039/d0sc05327a] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2020] [Accepted: 11/26/2020] [Indexed: 12/13/2022] Open
Abstract
To facilitate computational investigation of intermolecular interactions in the solution phase, we report the development of ALMO-EDA(solv), a scheme that allows the application of continuum solvent models within the framework of energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs). In this scheme, all the quantum mechanical states involved in the variational EDA procedure are computed with the presence of solvent environment so that solvation effects are incorporated in the evaluation of all its energy components. After validation on several model complexes, we employ ALMO-EDA(solv) to investigate substituent effects on two classes of complexes that are related to molecular CO2 reduction catalysis. For [FeTPP(CO2-κC)]2- (TPP = tetraphenylporphyrin), we reveal that two ortho substituents which yield most favorable CO2 binding, -N(CH3)3 + (TMA) and -OH, stabilize the complex via through-structure and through-space mechanisms, respectively. The coulombic interaction between the positively charged TMA group and activated CO2 is found to be largely attenuated by the polar solvent. Furthermore, we also provide computational support for the design strategy of utilizing bulky, flexible ligands to stabilize activated CO2 via long-range Coulomb interactions, which creates biomimetic solvent-inaccessible "pockets" in that electrostatics is unscreened. For the reactant and product complexes associated with the electron transfer from the p-terphenyl radical anion to CO2, we demonstrate that the double terminal substitution of p-terphenyl by electron-withdrawing groups considerably strengthens the binding in the product state while moderately weakens that in the reactant state, which are both dominated by the substituent tuning of the electrostatics component. These applications illustrate that this new extension of ALMO-EDA provides a valuable means to unravel the nature of intermolecular interactions and quantify their impacts on chemical reactivity in solution.
Collapse
Affiliation(s)
- Yuezhi Mao
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA
| | | | - Kareesa J Kron
- Mork Family Department of Chemical Engineering and Material Science, University of Southern California Los Angeles CA 90089 USA
| | - Jeffrey S Derrick
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA
| | - Christopher J Chang
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
- Department of Molecular and Cell Biology, University of California Berkeley Berkeley CA 94720 USA
| | - Shaama Mallikarjun Sharada
- Mork Family Department of Chemical Engineering and Material Science, University of Southern California Los Angeles CA 90089 USA
- Department of Chemistry, University of Southern California Los Angeles CA 90089 USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
| |
Collapse
|
15
|
Bononi FC, Chen Z, Rocca D, Andreussi O, Hullar T, Anastasio C, Donadio D. Bathochromic Shift in the UV–Visible Absorption Spectra of Phenols at Ice Surfaces: Insights from First-Principles Calculations. J Phys Chem A 2020; 124:9288-9298. [DOI: 10.1021/acs.jpca.0c07038] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Fernanda C. Bononi
- Department of Chemistry, University of California Davis, Davis, California 95616-5270, United States
| | - Zekun Chen
- Department of Chemistry, University of California Davis, Davis, California 95616-5270, United States
| | - Dario Rocca
- Université de Lorraine, CNRS, LPTC, F-54000 Nancy, France
| | - Oliviero Andreussi
- Department of Physics, University of North Texas Denton, Texas 76203, United States
| | - Ted Hullar
- Department of Land, Air and Water Resources, University of California Davis Davis, California 95616-8627, United States
| | - Cort Anastasio
- Department of Land, Air and Water Resources, University of California Davis Davis, California 95616-8627, United States
| | - Davide Donadio
- Department of Chemistry, University of California Davis, Davis, California 95616-5270, United States
| |
Collapse
|
16
|
Collinge G, Yuk SF, Nguyen MT, Lee MS, Glezakou VA, Rousseau R. Effect of Collective Dynamics and Anharmonicity on Entropy in Heterogenous Catalysis: Building the Case for Advanced Molecular Simulations. ACS Catal 2020. [DOI: 10.1021/acscatal.0c01501] [Citation(s) in RCA: 40] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Greg Collinge
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Simuck F. Yuk
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Manh-Thuong Nguyen
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Mal-Soon Lee
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Vassiliki-Alexandra Glezakou
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Roger Rousseau
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
17
|
Kühne TD, Iannuzzi M, Del Ben M, Rybkin VV, Seewald P, Stein F, Laino T, Khaliullin RZ, Schütt O, Schiffmann F, Golze D, Wilhelm J, Chulkov S, Bani-Hashemian MH, Weber V, Borštnik U, Taillefumier M, Jakobovits AS, Lazzaro A, Pabst H, Müller T, Schade R, Guidon M, Andermatt S, Holmberg N, Schenter GK, Hehn A, Bussy A, Belleflamme F, Tabacchi G, Glöß A, Lass M, Bethune I, Mundy CJ, Plessl C, Watkins M, VandeVondele J, Krack M, Hutter J. CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations. J Chem Phys 2020; 152:194103. [PMID: 33687235 DOI: 10.1063/5.0007045] [Citation(s) in RCA: 924] [Impact Index Per Article: 231.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.
Collapse
Affiliation(s)
- Thomas D Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Mauro Del Ben
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Vladimir V Rybkin
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Patrick Seewald
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Frederick Stein
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Teodoro Laino
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | - Rustam Z Khaliullin
- Department of Chemistry, McGill University, CH-801 Sherbrooke St. West, Montreal, Quebec H3A 0B8, Canada
| | - Ole Schütt
- Department of Materials, ETH Zürich, CH-8092 Zürich, Switzerland
| | | | - Dorothea Golze
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| | - Jan Wilhelm
- Institute of Theoretical Physics, University of Regensburg, Universitätsstraße 31, D-93053 Regensburg, Germany
| | - Sergey Chulkov
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | | | - Valéry Weber
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | | | | | | | | | - Hans Pabst
- Intel Extreme Computing, Software and Systems, Zürich, Switzerland
| | - Tiziano Müller
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Robert Schade
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Manuel Guidon
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Samuel Andermatt
- Integrated Systems Laboratory, ETH Zürich, CH-8092 Zürich, Switzerland
| | - Nico Holmberg
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076 Aalto, Finland
| | - Gregory K Schenter
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Anna Hehn
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Augustin Bussy
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Fabian Belleflamme
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Gloria Tabacchi
- Department of Science and High Technology, University of Insubria and INSTM, via Valleggio 9, I-22100 Como, Italy
| | - Andreas Glöß
- BASF SE, Carl-Bosch-Straße 38, D-67056 Ludwigshafen am Rhein, Germany
| | - Michael Lass
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Iain Bethune
- Hartree Centre, Science and Technology Facilities Council, Sci-Tech Daresbury, Warrington WA4 4AD, United Kingdom
| | - Christopher J Mundy
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Christian Plessl
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Matt Watkins
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | - Joost VandeVondele
- Swiss National Supercomputing Centre (CSCS), ETH Zürich, Zürich, Switzerland
| | - Matthias Krack
- Laboratory for Scientific Computing and Modelling, Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
18
|
Prentice JCA, Aarons J, Womack JC, Allen AEA, Andrinopoulos L, Anton L, Bell RA, Bhandari A, Bramley GA, Charlton RJ, Clements RJ, Cole DJ, Constantinescu G, Corsetti F, Dubois SMM, Duff KKB, Escartín JM, Greco A, Hill Q, Lee LP, Linscott E, O'Regan DD, Phipps MJS, Ratcliff LE, Serrano ÁR, Tait EW, Teobaldi G, Vitale V, Yeung N, Zuehlsdorff TJ, Dziedzic J, Haynes PD, Hine NDM, Mostofi AA, Payne MC, Skylaris CK. The ONETEP linear-scaling density functional theory program. J Chem Phys 2020; 152:174111. [PMID: 32384832 DOI: 10.1063/5.0004445] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.
Collapse
Affiliation(s)
- Joseph C A Prentice
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Jolyon Aarons
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - James C Womack
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Alice E A Allen
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Lampros Andrinopoulos
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Lucian Anton
- UKAEA, Culham Science Centre, Abingdon OX14 3DB, United Kingdom
| | - Robert A Bell
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Arihant Bhandari
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Gabriel A Bramley
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Robert J Charlton
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Rebecca J Clements
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Daniel J Cole
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| | - Gabriel Constantinescu
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Fabiano Corsetti
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Simon M-M Dubois
- Institute of Condensed Matter and Nanosciences, Université Catholique de Louvain, Louvain-la-Neuve, Belgium
| | - Kevin K B Duff
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - José María Escartín
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Andrea Greco
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Quintin Hill
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Louis P Lee
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Edward Linscott
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - David D O'Regan
- School of Physics, AMBER, and CRANN Institute, Trinity College Dublin, The University of Dublin, Dublin 2, Ireland
| | - Maximillian J S Phipps
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Laura E Ratcliff
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Álvaro Ruiz Serrano
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Edward W Tait
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Gilberto Teobaldi
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Valerio Vitale
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nelson Yeung
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Peter D Haynes
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nicholas D M Hine
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Arash A Mostofi
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Mike C Payne
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
19
|
Bramley G, Nguyen MT, Glezakou VA, Rousseau R, Skylaris CK. Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study. J Chem Theory Comput 2020; 16:2703-2715. [DOI: 10.1021/acs.jctc.0c00034] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Gabriel Bramley
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, U.K
| | - Manh-Thuong Nguyen
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | | - Roger Rousseau
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | |
Collapse
|
20
|
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
|
21
|
Nattino F, Dupont C, Marzari N, Andreussi O. Functional Extrapolations to Tame Unbound Anions in Density-Functional Theory Calculations. J Chem Theory Comput 2019; 15:6313-6322. [DOI: 10.1021/acs.jctc.9b00552] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
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, Lausanne CH-1015, Switzerland
| | - Céline Dupont
- Laboratoire Interdisciplinaire Carnot de Bourgogne (ICB), UMR 6303, CNRS, Université de Bourgogne Franche-Comté, BP 47870, Dijon Cedex 21078, France
| | - 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, Lausanne CH-1015, Switzerland
| | - Oliviero Andreussi
- Department of Physics, University of North Texas, Denton, Texas 76207, United States
| |
Collapse
|
22
|
Lange AW, Herbert JM, Albrecht BJ, You ZQ. Intrinsically smooth discretisation of Connolly's solvent-excluded molecular surface. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1644384] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Adrian W. Lange
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH, USA
| | - Benjamin J. Albrecht
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH, USA
| | - Zhi-Qiang You
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH, USA
| |
Collapse
|
23
|
Abstract
Bilayer graphene consists of two stacked graphene layers bound together by van der Waals interaction. As the molecular analog of bilayer graphene, molecular bilayer graphene (MBLG) can offer useful insights into the structural and functional properties of bilayer graphene. However, synthesis of MBLG, which requires discrete assembly of two graphene fragments, has proved to be challenging. Here, we show the synthesis and characterization of two structurally well-defined MBLGs, both consisting of two π-π stacked nanographene sheets. We find they have excellent stability against variation of concentration, temperature and solvents. The MBLGs show sharp absorption and emission peaks, and further time-resolved spectroscopic studies reveal drastically different lifetimes for the bright and dark Davydov states in these MBLGs.
Collapse
|
24
|
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
|
25
|
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
|
26
|
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
|
27
|
Andreussi O, Hörmann NG, Nattino F, Fisicaro G, Goedecker S, Marzari N. Solvent-Aware Interfaces in Continuum Solvation. J Chem Theory Comput 2019; 15:1996-2009. [DOI: 10.1021/acs.jctc.8b01174] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Oliviero Andreussi
- Department of Physics, University of North Texas, Denton, Texas 76207, United States
- 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, Station 9, CH-1015 Lausanne, Switzerland
| | - Nicolas Georg Hörmann
- 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, Station 9, CH-1015 Lausanne, Switzerland
| | - 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, Station 9, CH-1015 Lausanne, Switzerland
| | - Giuseppe Fisicaro
- Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland
- Istituto per la Microelettronica e Microsistemi (CNR-IMM), VIII Strada 5, 95121 Catania, Italy
| | - Stefan Goedecker
- Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland
| | - 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, Station 9, CH-1015 Lausanne, Switzerland
| |
Collapse
|
28
|
Castañeda Medina A, Schmid R. High Order Compact Multigrid Solver for Implicit Solvation Models. J Chem Theory Comput 2019; 15:1293-1301. [DOI: 10.1021/acs.jctc.8b00774] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Arcesio Castañeda Medina
- Computational Materials Chemistry Group, Lehrstuhl für Anorganische Chemie 2, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Rochus Schmid
- Computational Materials Chemistry Group, Lehrstuhl für Anorganische Chemie 2, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
29
|
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
|
30
|
Labat F, Civalleri B, Dovesi R. Implicit Solvation Using a Generalized Finite-Difference Approach in CRYSTAL: Implementation and Results for Molecules, Polymers, and Surfaces. J Chem Theory Comput 2018; 14:5969-5983. [PMID: 30347161 DOI: 10.1021/acs.jctc.8b00762] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the implementation of an implicit solvation model in the CRYSTAL code. The solvation energy is separated into two components: the electrostatic contribution arising from a self-consistent reaction field treatment obtained within a generalized finite-difference Poisson model, augmented by a nonelectrostatic contribution proportional to the solvent-accessible surface area of the solute. A discontinuous dielectric boundary is used, along with a solvent-excluded surface built from interlocking atom-centered spheres on which apparent surface point charges are mapped. The procedure is general and can be performed at both the Hartree-Fock and density functional theory levels, with pure or hybrid functionals, for systems periodic in 0, 1, and 2 directions, that is, for isolated molecules and extended polymers and surfaces. The Poisson equation resolution and apparent surface charge formalism is first validated on model analytical test cases. The good agreement obtained on solvation free energies is further confirmed by calculations performed on a large test set of 501 neutral molecules, for which a mean unsigned error of 1.3 kcal/mol is obtained when compared to the available experimental data. Importantly, the self-consistent reaction field procedure converges well for all molecules tested. This is further verified for all polymers and surfaces considered. In particular, for periodic systems, results obtained on an infinite glycine chain and on the wettability parameters of SiO2 surfaces are in good agreement with previously published data. The size extensivity of the energetic terms involved in the electrostatic contribution to the solvation energy is also well verified. These encouraging results constitute a first step to take into account complex environments in the CRYSTAL code, potentially allowing for a more accurate modeling of complex processes for both periodic and nonperiodic systems.
Collapse
Affiliation(s)
- Frédéric Labat
- PSL Research University, Chimie Paristech-CNRS , Institut de Recherche de Chimie de Paris , 11 rue P. et M. Curie , 75005 Paris , France
| | - Bartolomeo Civalleri
- Dipartimento di Chimica IFM , Università di Torino and NIS - Nanostructured Interfaces and Surfaces - Centre of Excellence , Via P. Giuria 7 , 10125 Torino , Italy
| | - Roberto Dovesi
- Dipartimento di Chimica IFM , Università di Torino and NIS - Nanostructured Interfaces and Surfaces - Centre of Excellence , Via P. Giuria 7 , 10125 Torino , Italy
| |
Collapse
|
31
|
Howard JC, Crawford TD. Calculating Optical Rotatory Dispersion Spectra in Solution Using a Smooth Dielectric Model. J Phys Chem A 2018; 122:8557-8564. [DOI: 10.1021/acs.jpca.8b07803] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- J. Coleman Howard
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T. Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
32
|
Coons MP, Herbert JM. Quantum chemistry in arbitrary dielectric environments: Theory and implementation of nonequilibrium Poisson boundary conditions and application to compute vertical ionization energies at the air/water interface. J Chem Phys 2018; 148:222834. [DOI: 10.1063/1.5023916] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Marc P. Coons
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
33
|
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
|
34
|
Yoo SH, Todorova M, Neugebauer J. Selective Solvent-Induced Stabilization of Polar Oxide Surfaces in an Electrochemical Environment. PHYSICAL REVIEW LETTERS 2018; 120:066101. [PMID: 29481276 DOI: 10.1103/physrevlett.120.066101] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/17/2017] [Revised: 12/10/2017] [Indexed: 06/08/2023]
Abstract
The impact of an electrochemical environment on the thermodynamic stability of polar oxide surfaces is investigated for the example of ZnO(0001) surfaces immersed in water using density functional theory calculations. We show that solvation effects are highly selective: They have little effect on surfaces showing a metallic character, but largely stabilize semiconducting structures, particularly those that have a high electrostatic penalty in vacuum. The high selectivity is shown to have direct consequences for the surface phase diagram and explains, e.g., why certain surface structures could be observed only in an electrochemical environment.
Collapse
Affiliation(s)
- Su-Hyun Yoo
- Department of Computational Materials Design, Max-Planck-Institut für Eisenforschung GmbH, Max-Planck-Str. 1, D-40237 Düsseldorf, Germany
| | - Mira Todorova
- Department of Computational Materials Design, Max-Planck-Institut für Eisenforschung GmbH, Max-Planck-Str. 1, D-40237 Düsseldorf, Germany
| | - Jörg Neugebauer
- Department of Computational Materials Design, Max-Planck-Institut für Eisenforschung GmbH, Max-Planck-Str. 1, D-40237 Düsseldorf, Germany
| |
Collapse
|
35
|
Giannozzi P, Andreussi O, Brumme T, Bunau O, Buongiorno Nardelli M, Calandra M, Car R, Cavazzoni C, Ceresoli D, Cococcioni M, Colonna N, Carnimeo I, Dal Corso A, de Gironcoli S, Delugas P, DiStasio RA, Ferretti A, Floris A, Fratesi G, Fugallo G, Gebauer R, Gerstmann U, Giustino F, Gorni T, Jia J, Kawamura M, Ko HY, Kokalj A, Küçükbenli E, Lazzeri M, Marsili M, Marzari N, Mauri F, Nguyen NL, Nguyen HV, Otero-de-la-Roza A, Paulatto L, Poncé S, Rocca D, Sabatini R, Santra B, Schlipf M, Seitsonen AP, Smogunov A, Timrov I, Thonhauser T, Umari P, Vast N, Wu X, Baroni S. Advanced capabilities for materials modelling with Quantum ESPRESSO. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:465901. [PMID: 29064822 DOI: 10.1088/1361-648x/aa8f79] [Citation(s) in RCA: 1505] [Impact Index Per Article: 215.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Quantum EXPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum EXPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.
Collapse
Affiliation(s)
- P Giannozzi
- Department of Mathematics, Computer Science, and Physics, University of Udine, via delle Scienze 206, I-33100 Udine, Italy
| | - O Andreussi
- Institute of Computational Sciences, Università della Svizzera Italiana, Lugano, Switzerland
- Theory and Simulation of Materials (THEOS), and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - T Brumme
- Wilhelm-Ostwald-Institute of Physical and Theoretical Chemistry, Leipzig University, Linnéstr. 2, D-04103 Leipzig, Germany
| | - O Bunau
- IMPMC, UMR CNRS 7590, Sorbonne Universités-UPMC University Paris 06, MNHN, IRD, 4 Place Jussieu, F-75005 Paris, France
| | - M Buongiorno Nardelli
- Department of Physics and Department of Chemistry, University of North Texas, Denton, TX, United States of America
| | - M Calandra
- IMPMC, UMR CNRS 7590, Sorbonne Universités-UPMC University Paris 06, MNHN, IRD, 4 Place Jussieu, F-75005 Paris, France
| | - R Car
- Department of Chemistry, Princeton University, Princeton, NJ 08544, United States of America
| | - C Cavazzoni
- CINECA-Via Magnanelli 6/3, I-40033 Casalecchio di Reno, Bologna, Italy
| | - D Ceresoli
- Institute of Molecular Science and Technologies (ISTM), National Research Council (CNR), I-20133 Milano, Italy
| | - M Cococcioni
- Theory and Simulation of Materials (THEOS), and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - N Colonna
- Theory and Simulation of Materials (THEOS), and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - I Carnimeo
- Department of Mathematics, Computer Science, and Physics, University of Udine, via delle Scienze 206, I-33100 Udine, Italy
| | - A Dal Corso
- SISSA-Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, I-34136 Trieste, Italy
- CNR-IOM DEMOCRITOS, Istituto Officina dei Materiali, Consiglio Nazionale delle Ricerche, Italy
| | - S de Gironcoli
- SISSA-Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, I-34136 Trieste, Italy
- CNR-IOM DEMOCRITOS, Istituto Officina dei Materiali, Consiglio Nazionale delle Ricerche, Italy
| | - P Delugas
- SISSA-Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, I-34136 Trieste, Italy
| | - R A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853, United States of America
| | - A Ferretti
- CNR Istituto Nanoscienze, I-42125 Modena, Italy
| | - A Floris
- School of Mathematics and Physics, College of Science, University of Lincoln, United Kingdom
| | - G Fratesi
- Dipartimento di Fisica, Università degli Studi di Milano, via Celoria 16, I-20133 Milano, Italy
| | - G Fugallo
- ETSF, Laboratoire des Solides Irradiés, Ecole Polytechnique, F-91128 Palaiseau cedex, France
| | - R Gebauer
- The Abdus Salam International Centre for Theoretical Physics (ICTP), Strada Costiera 11, I-34151 Trieste, Italy
| | - U Gerstmann
- Department Physik, Universität Paderborn, D-33098 Paderborn, Germany
| | - F Giustino
- Department of Materials, University of Oxford, Parks Road, Oxford OX1 3PH, United Kingdom
| | - T Gorni
- IMPMC, UMR CNRS 7590, Sorbonne Universités-UPMC University Paris 06, MNHN, IRD, 4 Place Jussieu, F-75005 Paris, France
- SISSA-Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, I-34136 Trieste, Italy
| | - J Jia
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853, United States of America
| | - M Kawamura
- The Institute for Solid State Physics, Kashiwa, Japan
| | - H-Y Ko
- Department of Chemistry, Princeton University, Princeton, NJ 08544, United States of America
| | - A Kokalj
- Department of Physical and Organic Chemistry, Jožef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia
| | - E Küçükbenli
- SISSA-Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, I-34136 Trieste, Italy
| | - M Lazzeri
- IMPMC, UMR CNRS 7590, Sorbonne Universités-UPMC University Paris 06, MNHN, IRD, 4 Place Jussieu, F-75005 Paris, France
| | - M Marsili
- Dipartimento di Fisica e Astronomia, Università di Padova, via Marzolo 8, I-35131 Padova, Italy
| | - N Marzari
- Theory and Simulation of Materials (THEOS), and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - F Mauri
- Dipartimento di Fisica, Università di Roma La Sapienza, Piazzale Aldo Moro 5, I-00185 Roma, Italy
| | - N L Nguyen
- Theory and Simulation of Materials (THEOS), and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - H-V Nguyen
- Institute of Physics, Vietnam Academy of Science and Technology, 10 Dao Tan, Hanoi, Vietnam
| | - A Otero-de-la-Roza
- Department of Chemistry, University of British Columbia, Okanagan, Kelowna BC V1V 1V7, Canada
| | - L Paulatto
- IMPMC, UMR CNRS 7590, Sorbonne Universités-UPMC University Paris 06, MNHN, IRD, 4 Place Jussieu, F-75005 Paris, France
| | - S Poncé
- Department of Materials, University of Oxford, Parks Road, Oxford OX1 3PH, United Kingdom
| | - D Rocca
- Université de Lorraine, CRM2, UMR 7036, F-54506 Vandoeuvre-lès-Nancy, France
- CNRS, CRM2, UMR 7036, F-54506 Vandoeuvre-lès-Nancy, France
| | - R Sabatini
- Orionis Biosciences, Newton, MA 02466, United States of America
| | - B Santra
- Department of Chemistry, Princeton University, Princeton, NJ 08544, United States of America
| | - M Schlipf
- Department of Materials, University of Oxford, Parks Road, Oxford OX1 3PH, United Kingdom
| | - A P Seitsonen
- Institut für Chimie, Universität Zurich, CH-8057 Zürich, Switzerland
- Département de Chimie, École Normale Supérieure, F-75005 Paris, France
| | - A Smogunov
- SPEC, CEA, CNRS, Université Paris-Saclay, F-91191 Gif-Sur-Yvette, France
| | - I Timrov
- Theory and Simulation of Materials (THEOS), and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - T Thonhauser
- Department of Physics, Wake Forest University, Winston-Salem, NC 27109, United States of America
| | - P Umari
- Dipartimento di Fisica e Astronomia, Università di Padova, via Marzolo 8, I-35131 Padova, Italy
- CNR-IOM DEMOCRITOS, Istituto Officina dei Materiali, Consiglio Nazionale delle Ricerche, Italy
| | - N Vast
- Laboratoire des Solides Irradiés, École Polytechnique, CEA-DRF-IRAMIS, CNRS UMR 7642, Université Paris-Saclay, F-91120 Palaiseau, France
| | - X Wu
- Department of Physics, Temple University, Philadelphia, PA 19122-1801, United States of America
| | - S Baroni
- SISSA-Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, I-34136 Trieste, Italy
| |
Collapse
|
36
|
Provorse Long MR, Isborn CM. Combining Explicit Quantum Solvent with a Polarizable Continuum Model. J Phys Chem B 2017; 121:10105-10117. [DOI: 10.1021/acs.jpcb.7b06693] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
| | - Christine M. Isborn
- Chemistry
and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
37
|
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
|
38
|
Sinstein M, Scheurer C, Matera S, Blum V, Reuter K, Oberhofer H. Efficient Implicit Solvation Method for Full Potential DFT. J Chem Theory Comput 2017; 13:5582-5603. [PMID: 28910530 DOI: 10.1021/acs.jctc.7b00297] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
With the advent of efficient electronic structure methods, effective continuum solvation methods have emerged as a way to, at least partially, include solvent effects into simulations without the need for expensive sampling over solvent degrees of freedom. The multipole moment expansion (MPE) model, while based on ideas initially put forward almost 100 years ago, has recently been updated for the needs of modern electronic structure calculations. Indeed, for an all-electron code relying on localized basis sets and-more importantly-a multipole moment expansion of the electrostatic potential, the MPE method presents a particularly cheap way of solving the macroscopic Poisson equation to determine the electrostatic response of a medium surrounding a solute. In addition to our implementation of the MPE model in the FHI-aims electronic structure theory code [ Blum , V. ; Comput. Phys. Commun. 2009 , 180 , 2175 - 2196 , DOI: 10.1016/j.cpc.2009.06.022 ], we describe novel algorithms for determining equidistributed points on the solvation cavity-defined as a charge density isosurface-and the determination of cavity surface and volume from just this collection of points and their local density gradients. We demonstrate the efficacy of our model on an analytically solvable test case, against high-accuracy finite-element calculations for a set of ≈140000 2D test cases, and finally against experimental solvation free energies of a number of neutral and singly charged molecular test sets [ Andreussi , O. ; J. Chem. Phys. 2012 , 136 , 064102 , DOI: 10.1063/1.3676407 ; Marenich , A. V. ; Minnesota Solvation Database , Version 2012; University of Minnesota : Minneapolis, MN, USA , 2012 . ]. In all test cases we find that our MPE approach compares very well with given references at computational overheads < 20% and sometimes much smaller compared to a plain self-consistency cycle.
Collapse
Affiliation(s)
- Markus Sinstein
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Christoph Scheurer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Sebastian Matera
- Institut für Mathematik, Freie Universität Berlin , Arnimallee 9, D-14195 Berlin, Germany
| | - Volker Blum
- Department of Mechanical Engineering and Materials Science, Duke University , Durham, North Carolina 27708, United States
| | - Karsten Reuter
- 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
| |
Collapse
|
39
|
Fisicaro G, Genovese L, Andreussi O, Mandal S, Nair NN, Marzari N, Goedecker S. Soft-Sphere Continuum Solvation in Electronic-Structure Calculations. J Chem Theory Comput 2017. [DOI: 10.1021/acs.jctc.7b00375] [Citation(s) in RCA: 65] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Giuseppe Fisicaro
- Department
of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland
| | - Luigi Genovese
- Laboratoire
de simulation atomistique (L_Sim), SP2M, INAC, CEA-UJF, F-38054 Grenoble, France
| | - Oliviero Andreussi
- Institute
of Computational Science, Università della Svizzera Italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland
- 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, Station 12, CH-1015 Lausanne, Switzerland
| | - Sagarmoy Mandal
- Department
of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Nisanth N. Nair
- Department
of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - 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, Station 12, CH-1015 Lausanne, Switzerland
| | - Stefan Goedecker
- Department
of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland
| |
Collapse
|
40
|
Mao Y, Shao Y, Dziedzic J, Skylaris CK, Head-Gordon T, Head-Gordon M. Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis. J Chem Theory Comput 2017; 13:1963-1979. [PMID: 28430427 DOI: 10.1021/acs.jctc.7b00089] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The importance of incorporating solvent polarization effects into the modeling of solvation processes has been well-recognized, and therefore a new generation of hybrid quantum mechanics/molecular mechanics (QM/MM) approaches that accounts for this effect is desirable. We present a fully self-consistent, mutually polarizable QM/MM scheme using the AMOEBA force field, in which the total energy of the system is variationally minimized with respect to both the QM electronic density and the MM induced dipoles. This QM/AMOEBA model is implemented through the Q-Chem/LibEFP code interface and then applied to the evaluation of solute-solvent interaction energies for various systems ranging from the water dimer to neutral and ionic solutes (NH3, NH4+, CN-) surrounded by increasing numbers of water molecules (up to 100). In order to analyze the resulting interaction energies, we also utilize an energy decomposition analysis (EDA) scheme which identifies contributions from permanent electrostatics, polarization, and van der Waals (vdW) interaction for the interaction between the QM solute and the solvent molecules described by AMOEBA. This facilitates a component-wise comparison against full QM calculations where the corresponding energy components are obtained via a modified version of the absolutely localized molecular orbitals (ALMO)-EDA. The results show that the present QM/AMOEBA model can yield reasonable solute-solvent interaction energies for neutral and cationic species, while further scrutiny reveals that this accuracy highly relies on the delicate balance between insufficiently favorable permanent electrostatics and softened vdW interaction. For anionic solutes where the charge penetration effect becomes more pronounced, the QM/MM interface turns out to be unbalanced. These results are consistent with and further elucidate our findings in a previous study using a slightly different QM/AMOEBA model ( Dziedzic et al. J. Chem. Phys. 2016 , 145 , 124106 ). The implications of these results for further refinement of this model are also discussed.
Collapse
Affiliation(s)
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma , Norman, Oklahoma 73019, United States
| | - Jacek Dziedzic
- School 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
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, U.K
| | | | - Martin Head-Gordon
- Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
41
|
Ringe S, Oberhofer H, Reuter K. Transferable ionic parameters for first-principles Poisson-Boltzmann solvation calculations: Neutral solutes in aqueous monovalent salt solutions. J Chem Phys 2017; 146:134103. [DOI: 10.1063/1.4978850] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Affiliation(s)
- Stefan Ringe
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| |
Collapse
|
42
|
Solubility study of cefpodoxime acid antibiotic in terms of free energy of solution - Insights from polarizable continuum model (PCM) analysis. J Mol Liq 2016. [DOI: 10.1016/j.molliq.2016.10.019] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
|
43
|
Provorse MR, Peev T, Xiong C, Isborn CM. Convergence of Excitation Energies in Mixed Quantum and Classical Solvent: Comparison of Continuum and Point Charge Models. J Phys Chem B 2016; 120:12148-12159. [DOI: 10.1021/acs.jpcb.6b09176] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Makenzie R. Provorse
- Chemistry and Chemical Biology, University of California at Merced, Merced, California 95343, United States
| | - Thomas Peev
- Chemistry and Chemical Biology, University of California at Merced, Merced, California 95343, United States
| | - Chou Xiong
- Chemistry and Chemical Biology, University of California at Merced, Merced, California 95343, United States
| | - Christine M. Isborn
- Chemistry and Chemical Biology, University of California at Merced, Merced, California 95343, United States
| |
Collapse
|
44
|
Cole DJ, Hine NDM. Applications of large-scale density functional theory in biology. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:393001. [PMID: 27494095 DOI: 10.1088/0953-8984/28/39/393001] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Density functional theory (DFT) has become a routine tool for the computation of electronic structure in the physics, materials and chemistry fields. Yet the application of traditional DFT to problems in the biological sciences is hindered, to a large extent, by the unfavourable scaling of the computational effort with system size. Here, we review some of the major software and functionality advances that enable insightful electronic structure calculations to be performed on systems comprising many thousands of atoms. We describe some of the early applications of large-scale DFT to the computation of the electronic properties and structure of biomolecules, as well as to paradigmatic problems in enzymology, metalloproteins, photosynthesis and computer-aided drug design. With this review, we hope to demonstrate that first principles modelling of biological structure-function relationships are approaching a reality.
Collapse
Affiliation(s)
- Daniel J Cole
- Theory of Condensed Matter group, Cavendish Laboratory, 19 JJ Thomson Ave, Cambridge CB3 0HE, UK. School of Chemistry, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
| | | |
Collapse
|
45
|
Dziedzic J, Mao Y, Shao Y, Ponder J, Head-Gordon T, Head-Gordon M, Skylaris CK. TINKTEP: A fully self-consistent, mutually polarizable QM/MM approach based on the AMOEBA force field. J Chem Phys 2016; 145:124106. [DOI: 10.1063/1.4962909] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
- Faculty of Applied Physics and Mathematics, Gdańsk University of Technology, Gdańsk, Poland
| | - Yuezhi Mao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yihan Shao
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Jay Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, USA
| | - Teresa Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
- Department of Bioengineering, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
46
|
Montemore MM, Andreussi O, Medlin JW. Hydrocarbon adsorption in an aqueous environment: A computational study of alkyls on Cu(111). J Chem Phys 2016; 145:074702. [PMID: 27544118 DOI: 10.1063/1.4961027] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Hydrocarbon chains are important intermediates in various aqueous-phase surface processes, such as CO2 electroreduction, aqueous Fischer-Tropsch synthesis, and aqueous phase reforming of biomass-derived molecules. Further, the interaction between water and adsorbed hydrocarbons represents a difficult case for modern computational methods. Here, we explore various methods for calculating the energetics of this interaction within the framework of density functional theory and explore trade-offs between the use of low water coverages, molecular dynamics approaches, and minima hopping for identification of low energy structures. An effective methodology for simulating low temperature processes is provided by using a unit cell in which the vacuum space is filled with water, employing the minima hopping algorithm to search for low-lying minima, and including dispersion (van der Waals) interactions. Using this methodology, we show that a high coverage of adsorbed alkyls is destabilized by the presence of water, while a low coverage of alkyls is stabilized. Solvation has a small effect on the energetics of hydrocarbon chain growth, generally decreasing its favorability at low temperatures. We studied higher temperatures by running molecular dynamics simulations starting at the minima found by the minima hopping algorithm and found that increased temperatures facilitate chain growth. The self-consistent continuum solvation method effectively describes the alkyl-water interaction and is in general agreement with the explicit solvation results in most cases, but care should be taken at high alkyl coverage.
Collapse
Affiliation(s)
- Matthew M Montemore
- Department of Mechanical Engineering, University of Colorado Boulder, UCB 427, Boulder, Colorado 80309, USA
| | - Oliviero Andreussi
- Faculty of Informatics, Institute of Computational Science, Università della Svizzera Italiana, Via G. Buffi 13, 6904 Lugano, Switzerland
| | - J Will Medlin
- Department of Chemical and Biological Engineering, University of Colorado Boulder, UCB 596, Boulder, Colorado 80309, USA
| |
Collapse
|
47
|
Timrov I, Micciarelli M, Rosa M, Calzolari A, Baroni S. Multimodel Approach to the Optical Properties of Molecular Dyes in Solution. J Chem Theory Comput 2016; 12:4423-9. [DOI: 10.1021/acs.jctc.6b00417] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Iurii Timrov
- SISSA − Scuola Internazionale
Superiore di Studi Avanzati, Via Bonomea
265, 34136 Trieste, Italy
| | - Marco Micciarelli
- SISSA − Scuola Internazionale
Superiore di Studi Avanzati, Via Bonomea
265, 34136 Trieste, Italy
| | - Marta Rosa
- SISSA − Scuola Internazionale
Superiore di Studi Avanzati, Via Bonomea
265, 34136 Trieste, Italy
| | - Arrigo Calzolari
- CNR-NANO, Istituto Nanoscienze, Centro S3, Via Campi 213A, 41125 Modena, Italy
| | - Stefano Baroni
- SISSA − Scuola Internazionale
Superiore di Studi Avanzati, Via Bonomea
265, 34136 Trieste, Italy
| |
Collapse
|
48
|
Ringe S, Oberhofer H, Hille C, Matera S, Reuter K. Function-Space-Based Solution Scheme for the Size-Modified Poisson-Boltzmann Equation in Full-Potential DFT. J Chem Theory Comput 2016; 12:4052-66. [PMID: 27323006 DOI: 10.1021/acs.jctc.6b00435] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The size-modified Poisson-Boltzmann (MPB) equation is an efficient implicit solvation model which also captures electrolytic solvent effects. It combines an account of the dielectric solvent response with a mean-field description of solvated finite-sized ions. We present a general solution scheme for the MPB equation based on a fast function-space-oriented Newton method and a Green's function preconditioned iterative linear solver. In contrast to popular multigrid solvers, this approach allows us to fully exploit specialized integration grids and optimized integration schemes. We describe a corresponding numerically efficient implementation for the full-potential density-functional theory (DFT) code FHI-aims. We show that together with an additional Stern layer correction the DFT+MPB approach can describe the mean activity coefficient of a KCl aqueous solution over a wide range of concentrations. The high sensitivity of the calculated activity coefficient on the employed ionic parameters thereby suggests to use extensively tabulated experimental activity coefficients of salt solutions for a systematic parametrization protocol.
Collapse
Affiliation(s)
- Stefan Ringe
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Christoph Hille
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Sebastian Matera
- Fachbereich f. Mathematik u. Informatik, Freie Universität Berlin , Otto-von-Simson-Str. 19, D-14195 Berlin, Germany
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstr. 4, D-85747 Garching, Germany
| |
Collapse
|
49
|
Lipparini F, Mennucci B. Perspective: Polarizable continuum models for quantum-mechanical descriptions. J Chem Phys 2016; 144:160901. [DOI: 10.1063/1.4947236] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Filippo Lipparini
- Institut für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, D55128 Mainz, Germany
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
50
|
Zuehlsdorff TJ, Haynes PD, Hanke F, Payne MC, Hine NDM. Solvent Effects on Electronic Excitations of an Organic Chromophore. J Chem Theory Comput 2016; 12:1853-61. [DOI: 10.1021/acs.jctc.5b01014] [Citation(s) in RCA: 52] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- T. J. Zuehlsdorff
- Cavendish
Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, U.K
| | | | - F. Hanke
- Dassault Systèmes
BIOVIA, 334 Science Park, Cambridge CB4 0WN, U.K
| | - M. C. Payne
- Cavendish
Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, U.K
| | - N. D. M. Hine
- Department
of Physics, University of Warwick, Coventry CV4 7AL, U.K
| |
Collapse
|