A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments.

The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes.

[1]  Reinhold Schneider,et al.  Daubechies wavelets as a basis set for density functional pseudopotential calculations. , 2008, The Journal of chemical physics.

[2]  Stefan Goedecker,et al.  Efficient and accurate three-dimensional Poisson solver for surface problems. , 2007, The Journal of chemical physics.

[3]  Fang Liu,et al.  Quantum Chemistry for Solvated Molecules on Graphical Processing Units Using Polarizable Continuum Models. , 2015, Journal of chemical theory and computation.

[4]  S. Goedecker.  Minima hopping: an efficient search method for the global minimum of the potential energy surface of complex molecular systems. , 2004, The Journal of chemical physics.

[5]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[6]  J. J. López-García,et al.  Poisson-Boltzmann description of the electrical double layer including ion size effects. , 2011, Langmuir : the ACS journal of surfaces and colloids.

[7]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[8]  Alfred B. Anderson,et al.  Electronic structure calculations of liquid-solid interfaces: Combination of density functional theory and modified Poisson-Boltzmann theory , 2008 .

[9]  Álvaro Vázquez-Mayagoitia,et al.  Norm-conserving pseudopotentials with chemical accuracy compared to all-electron calculations. , 2012, The Journal of chemical physics.

[10]  D. Chapman,et al.  LI. A contribution to the theory of electrocapillarity , 1913 .

[11]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

[12]  Teter,et al.  Separable dual-space Gaussian pseudopotentials. , 1996, Physical review. B, Condensed matter.

[13]  Luigi Genovese,et al.  Efficient and accurate solver of the three-dimensional screened and unscreened Poisson's equation with generic boundary conditions. , 2012, The Journal of chemical physics.

[14]  Nicola Marzari,et al.  A unified electrostatic and cavitation model for first-principles molecular dynamics in solution. , 2006, The Journal of chemical physics.

[15]  Stefan Goedecker,et al.  Daubechies wavelets for linear scaling density functional theory. , 2014, The Journal of chemical physics.

[16]  Stefan Goedecker,et al.  Minima hopping guided path search: an efficient method for finding complex chemical reaction pathways. , 2014, The Journal of chemical physics.

[17]  Minoru Otani,et al.  First-principles calculations of charged surfaces and interfaces: A plane-wave nonrepeated slab approach , 2006 .

[18]  Jean-François Méhaut,et al.  High Performance Computing / Le Calcul Intensif Daubechies wavelets for high performance electronic structure calculations: The BigDFT project , 2011 .

[19]  F. J. Luque,et al.  Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. , 2000, Chemical reviews.

[20]  F. Gygi,et al.  The solvation of Na+ in water: First-principles simulations , 2000 .

[21]  Nicola Marzari,et al.  Ab Initio Electrochemical Properties of Electrode Surfaces , 2010 .

[22]  C. Cramer,et al.  Implicit Solvation Models: Equilibria, Structure, Spectra, and Dynamics. , 1999, Chemical reviews.

[23]  Stefan Goedecker,et al.  Efficient solution of Poisson's equation with free boundary conditions. , 2006, The Journal of chemical physics.

[24]  J. J. Bikerman,et al.  XXXIX. Structure and capacity of electrical double layer , 1942 .

[25]  Jean-Luc Fattebert,et al.  Density functional theory for efficient ab initio molecular dynamics simulations in solution , 2002, J. Comput. Chem..

[26]  L. Onsager.  Electric Moments of Molecules in Liquids , 1936 .

[27]  Nicola Marzari,et al.  Revised self-consistent continuum solvation in electronic-structure calculations. , 2011, The Journal of chemical physics.

[28]  Jean-Luc Fattebert,et al.  First‐principles molecular dynamics simulations in a continuum solvent , 2003 .

[29]  Michael J. Frisch,et al.  Achieving Linear Scaling for the Electronic Quantum Coulomb Problem , 1996, Science.

[30]  Nino Russo,et al.  Solvation effects on reaction profiles by the polarizable continuum model coupled with the Gaussian density functional method , 1998 .

[31]  J. Tomasi,et al.  Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects , 1981 .

[32]  Jacopo Tomasi,et al.  Molecular Interactions in Solution: An Overview of Methods Based on Continuous Distributions of the Solvent , 1994 .

[33]  M. Gouy,et al.  Sur la constitution de la charge électrique à la surface d'un électrolyte , 1910 .

[34]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[35]  Stefan Goedecker,et al.  A customized 3D GPU Poisson solver for free boundary conditions , 2013, Comput. Phys. Commun..