Revised self-consistent continuum solvation in electronic-structure calculations

Abstract

The solvation model proposed by Fattebert and Gygi [Journal of Computational Chemistry 23, 662 (2002)] and Scherlis et al. [Journal of Chemical Physics 124, 074103 (2006)] is reformulated, overcoming some of the numerical limitations encountered and extending its range of applicability. We first recast the problem in terms of induced polarization charges that act as a direct mapping of the self-consistent continuum dielectric; this allows to define a functional form for the dielectric that is well behaved both in the high-density region of the nuclear charges and in the low-density region where the electronic wavefunctions decay into the solvent. Second, we outline an iterative procedure to solve the Poisson equation for the quantum fragment embedded in the solvent that does not require multi-grid algorithms, is trivially parallel, and can be applied to any Bravais crystallographic system. Last, we capture some of the non-electrostatic or cavitation terms via a combined use of the quantum volume and quantum surface [Physical Review Letters 94, 145501 (2005)] of the solute. The resulting self-consistent continuum solvation (SCCS) model provides a very effective and compact fit of computational and experimental data, whereby the static dielectric constant of the solvent and one parameter allow to fit the electrostatic energy provided by the PCM model with a mean absolute error of 0.3 kcal/mol on a set of 240 neutral solutes. Two parameters allow to fit experimental solvation energies on the same set with a mean absolute error of 1.3 kcal/mol. A detailed analysis of these results, broken down along different classes of chemical compounds, shows that several classes of organic compounds display very high accuracy, with solvation energies in error of 0.3-0.4 kcal/mol, whereby larger discrepancies are mostly limited to self-dissociating species and strong hydrogen-bond forming compounds.

The following article has been accepted by The Journal of Chemical Physics. After it is published, it will be found at http://link.aip.org/link/?jcp/

References in corpus (6)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Nuclear quantum effects in water
• Electrostatics in Periodic Boundary Conditions and Real-space Corrections
• Minimal parameter implicit solvent model for ab initio electronic structure calculations
• Electronic-enthalpy functional for finite systems under pressure
• Electronic levels and electrical response of periodic molecular structures from plane-wave orbital-dependent calculations

Cited by in corpus (55)

• Advanced capabilities for materials modelling with Quantum ESPRESSO
• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways
• Grand canonical electronic density-functional theory: algorithms and applications to electrochemistry
• JDFTx: software for joint density-functional theory
• An In Situ Surface-Enhanced Infrared Absorption Spectroscopy Study of Electrochemical CO2 Reduction: Selectivity Dependence on Surface C-Bound and O-Bound Reaction Intermediates
• The charge-asymmetric nonlocally-determined local-electric (CANDLE) solvation model
• The importance of nonlinear fluid response in joint density-functional theory studies of battery systems
• Dielectric continuum methods for quantum chemistry
• First-principles theory of field-effect doping in transition-metal dichalcogenides: Structural properties, electronic structure, Hall coefficient, and electrical conductivity
• A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments
• Evaluating continuum solvation models for the electrode-electrolyte interface: challenges and strategies for improvement
• Electrostatics of solvated systems in periodic boundary conditions
• Quantum-continuum simulation of underpotential deposition at electrified metal-solution interfaces
• Continuum models of the electrochemical diffuse layer in electronic-structure calculations
• Function-Space Based Solution Scheme for the Size-Modified Poisson-Boltzmann Equation in Full-Potential DFT
• Towards First Principles prediction of Voltage Dependences of Electrolyte/Electrolyte Interfacial Processes in Lithium Ion Batteries
• Roadmap on Electronic Structure Codes in the Exascale Era
• Advancement of the Homogeneous Background Method for the Computational Simulation of Electrochemical Interfaces
• Spicing up continuum solvation models with SaLSA: the spherically-averaged liquid susceptibility ansatz
• Improving Accuracy of Electrochemical Capacitance and Solvation Energetics in First-Principles Calculations
• Janus $\beta$-Te${}_2$X (X = S, Se) Monolayers for Efficient Excitonic Solar Cells and Photocatalytic Water Splitting
• Two-Dimensional $\beta$-PdX${}_2$ (X = S, Te) Monolayers for Efficient Solar Energy Conversion Applications
• The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients
• Quantum-continuum simulation of the electrochemical response of pseudocapacitor electrodes from realistic conditions
• PCMSolver: an Open-Source Library for Solvation Modeling
• Weighted-density functionals for cavity formation and dispersion energies in continuum solvation models
• A recipe for free-energy functionals of polarizable molecular fluids
• Operando XANES from first-principles and its application to iridium oxide
• A first-principle calculation of the XANES spectrum of Cu${}^{2+}$ in water
• Existence of pearled patterns in the planar Functionalized Cahn-Hilliard equation
• Voltage-dependent cluster expansion for electrified solid-liquid interfaces: Application to the electrochemical deposition of transition metals
• Optimization of constrained density functional theory
• Substrate effects on charged defects in two-dimensional materials
• Accurate non-empirical range-separated hybrid van der Waals density functional for complex molecular problems, solids, and surfaces
• Multimodel Approach to the Optical Properties of Molecular Dyes in Solution
• Voltage-dependent reconstruction of layered Bi${}_2$WO${}_6$ and Bi${}_2$MoO${}_6$ photocatalysts and its influence on charge separation for water splitting
• Quantum-continuum calculation of the surface states and electrical response of silicon in solution
• Electrochemical oxidation of CO on Cu single crystals under alkaline conditions
• Competition and Complexity in Amphiphilic Polymer Morphology
• First principles approaches and concepts for electrochemical systems
• Functional extrapolations to tame unbound anions in density-functional theory calculations
• Multi-scale theoretical approach to X-ray absorption spectra in disordered systems: an application to the study of Zn(II) in water
• Modeling solvation effects in real-space and real-time within Density Functional Approaches
• Voltage effects on the stability of Pd ensembles in Pd-Au/Au(111) surface alloys
• High resistance of superconducting TiN thin films against environmental attacks
• Gaussian-Based Periodic Grand Canonical Density Functional Theory with Implicit Solvation for Computational Electrochemistry
• Voltage-dependent first-principles simulation of insertion of chloride ions into Al/Al${}_2$O${}_3$ interfaces using the Quantum Continuum Approximation
• Resilience of gas-phase anharmonicity in the vibrational response of adsorbed carbon monoxide and breakdown under electrical conditions
• Influence of surface restructuring on the activity of SrTiO3 photoelectrodes for photocatalytic hydrogen reduction
• Quantifying multi-point ordering in alloys
• Phase engineering of MoS${}_2$ monolayers: A pathway to enhanced lithium-polysulfide battery performance
• Pressure Correction for Solvation Theories
• Electrochemical stability and light-harvesting ability of silicon photoelectrodes in aqueous environments
• Uncoupling System and Environment Simulation Cells for Fast-Scaling Modeling of Complex Continuum Embeddings