Quantum, classical, and hybrid QM/MM calculations in solution: general implementation of the ddCOSMO linear scaling strategy.

We present the general theory and implementation of the Conductor-like Screening Model according to the recently developed ddCOSMO paradigm. The various quantities needed to apply ddCOSMO at different levels of theory, including quantum mechanical descriptions, are discussed in detail, with a particular focus on how to compute the integrals needed to evaluate the ddCOSMO solvation energy and its derivatives. The overall computational cost of a ddCOSMO computation is then analyzed and decomposed in the various steps: the different relative weights of such contributions are then discussed for both ddCOSMO and the fastest available alternative discretization to the COSMO equations. Finally, the scaling of the cost of the various steps with respect to the size of the solute is analyzed and discussed, showing how ddCOSMO opens significantly new possibilities when cheap or hybrid molecular mechanics/quantum mechanics methods are used to describe the solute.

[1]  Roberto Cammi,et al.  Continuum Solvation Models in Chemical Physics , 2007 .

[2]  J. Kongsted,et al.  Electronic Energy Transfer in Condensed Phase Studied by a Polarizable QM/MM Model. , 2009, Journal of chemical theory and computation.

[3]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[4]  Jacopo Tomasi,et al.  A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics , 1997 .

[5]  Walter Thiel,et al.  Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 1. Implementation. , 2005, The journal of physical chemistry. B.

[6]  A. Becke.  A multicenter numerical integration scheme for polyatomic molecules , 1988 .

[7]  Martin Karplus,et al.  A Smooth Solvation Potential Based on the Conductor-Like Screening Model , 1999 .

[8]  Mark S. Gordon,et al.  Self‐consistent molecular orbital methods. XXIII. A polarization‐type basis set for second‐row elements , 1982 .

[9]  Giovanni Scalmani,et al.  Energies, structures, and electronic properties of molecules in solution with the C‐PCM solvation model , 2003, J. Comput. Chem..

[10]  Giovanni Scalmani,et al.  Non covalent interactions in RNA and DNA base pairs: a quantum-mechanical study of the coupling between solvent and electronic density. , 2009, Physical chemistry chemical physics : PCCP.

[11]  Benjamin Stamm,et al.  Quantum Calculations in Solution for Large to Very Large Molecules: A New Linear Scaling QM/Continuum Approach. , 2014, The journal of physical chemistry letters.

[12]  Giovanni Scalmani,et al.  Continuous surface charge polarizable continuum models of solvation. I. General formalism. , 2010, The Journal of chemical physics.

[13]  Giovanni Scalmani,et al.  Analytical First and Second Derivatives for a Fully Polarizable QM/Classical Hamiltonian. , 2012, Journal of chemical theory and computation.

[14]  Walter Thiel,et al.  Solvent Boundary Potentials for Hybrid QM/MM Computations Using Classical Drude Oscillators: A Fully Polarizable Model. , 2012, Journal of chemical theory and computation.

[15]  K. N. Dollman,et al.  - 1 , 1743 .

[16]  G. Scuseria,et al.  Achieving linear-scaling computational cost for the polarizable continuum model of solvation , 2004 .

[17]  Benjamin Stamm,et al.  Scalable Evaluation of Polarization Energy and Associated Forces in Polarizable Molecular Dynamics: I. Toward Massively Parallel Direct Space Computations , 2014 .

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

[19]  Giovanni Scalmani,et al.  A variational formulation of the polarizable continuum model. , 2010, The Journal of chemical physics.

[20]  Benedetta Mennucci,et al.  Polarizable continuum model , 2012 .

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

[22]  C. Cramer,et al.  Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. , 2009, The journal of physical chemistry. B.

[23]  Benedetta Mennucci,et al.  The escaped charge problem in solvation continuum models , 2001 .

[24]  D. Chipman.  Simulation of volume polarization in reaction field theory , 1999 .

[25]  Luca Frediani,et al.  Excitation energies in solution: the fully polarizable QM/MM/PCM method. , 2011, The journal of physical chemistry. B.

[26]  K. Diederichs,et al.  Structural Basis of Light Harvesting by Carotenoids: Peridinin-Chlorophyll-Protein from Amphidinium carterae , 1996, Science.

[27]  John M. Herbert,et al.  Polarizable Continuum Reaction-Field Solvation Models Affording Smooth Potential Energy Surfaces , 2010 .

[28]  John M Herbert,et al.  A smooth, nonsingular, and faithful discretization scheme for polarizable continuum models: the switching/Gaussian approach. , 2010, Journal of Chemical Physics.

[29]  Benjamin Stamm,et al.  Fast Domain Decomposition Algorithm for Continuum Solvation Models: Energy and First Derivatives. , 2013, Journal of chemical theory and computation.

[30]  Benedetta Mennucci,et al.  Toward a Unified Modeling of Environment and Bridge-Mediated Contributions to Electronic Energy Transfer: A Fully Polarizable QM/MM/PCM Approach. , 2012, Journal of chemical theory and computation.

[31]  Benjamin Stamm,et al.  Domain decomposition for implicit solvation models. , 2013, The Journal of chemical physics.

[32]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[33]  Benedetta Mennucci,et al.  Continuum Solvation Models: What Else Can We Learn from Them? , 2010 .

[34]  J. Stewart.  Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements , 2007, Journal of molecular modeling.

[35]  B. Mennucci,et al.  Comprar Continuum Solvation Models in Chemical Physics: Theory and Applications | Benedetta Mennucci | 9780470029381 | Wiley , 2007 .

[36]  R. Cammi.  Quantum cluster theory for the polarizable continuum model. I. The CCSD level with analytical first and second derivatives. , 2009, The Journal of chemical physics.

[37]  Giovanni Scalmani,et al.  Self-Consistent Field and Polarizable Continuum Model: A New Strategy of Solution for the Coupled Equations. , 2011, Journal of chemical theory and computation.

[38]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

[39]  Vincenzo Barone,et al.  A gauge invariant multiscale approach to magnetic spectroscopies in condensed phase: general three-layer model, computational implementation and pilot applications. , 2013, The Journal of chemical physics.

[40]  Vincenzo Barone,et al.  Polarizable Force Fields and Polarizable Continuum Model: A Fluctuating Charges/PCM Approach. 1. Theory and Implementation. , 2011, Journal of chemical theory and computation.

[41]  A. Klamt.  The COSMO and COSMO‐RS solvation models , 2011 .

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

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

[44]  Vincenzo Barone,et al.  Linear Response Theory and Electronic Transition Energies for a Fully Polarizable QM/Classical Hamiltonian. , 2012, Journal of chemical theory and computation.

[45]  M. Frisch,et al.  Electronic excitation energies in solution at equation of motion CCSD level within a state specific polarizable continuum model approach. , 2010, The Journal of chemical physics.

[46]  Benedetta Mennucci,et al.  New applications of integral equations methods for solvation continuum models: ionic solutions and liquid crystals , 1998 .

[47]  Michael J. Frisch,et al.  Achieving linear scaling in exchange-correlation density functional quadratures , 1996 .