Quantum Calculations in Solution for Large to Very Large Molecules: A New Linear Scaling QM/Continuum Approach.

We present a new implementation of continuum solvation models for semiempirical Hamiltonians that allows the description of environmental effects on very large molecular systems. In this approach based on a domain decomposition strategy of the COSMO model (ddCOSMO), the solution to the COSMO equations is no longer the computational bottleneck but becomes a negligible part of the overall computation time. In this Letter, we analyze the computational impact of COSMO on the solution of the SCF equations for large to very large molecules, using semiempirical Hamiltonians, for both the new ddCOSMO implementation and the most recent, linear scaling one, based on the fast multipole method. A further analysis is on the simulation of the UV/visible spectrum of a light-harvesting pigment-protein complex. All of the results show how the new ddCOSMO algorithm paves the way to routine computations for large molecular systems in the condensed phase.

[1]  Benedetta Mennucci,et al.  Modeling environment effects on spectroscopies through QM/classical models. , 2013, Physical chemistry chemical physics : PCCP.

[2]  Jacob Kongsted,et al.  Molecular Properties through Polarizable Embedding , 2011 .

[3]  G. Fleming,et al.  Solvation Dynamics in Protein Environments Studied by Photon Echo Spectroscopy , 1999 .

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

[5]  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.

[6]  Donald G. Truhlar,et al.  Implicit Solvation Models: Equilibria, Structure, Spectra, and Dynamics , 1999 .

[7]  Michael C. Zerner,et al.  An intermediate neglect of differential overlap technique for spectroscopy: Pyrrole and the azines , 1973 .

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

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

[10]  D. Truhlar,et al.  QM/MM: what have we learned, where are we, and where do we go from here? , 2007 .

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

[12]  V. Sundström,et al.  Carotenoid to chlorophyll energy transfer in the peridinin–chlorophyll-a–protein complex involves an intramolecular charge transfer state , 2002, Proceedings of the National Academy of Sciences of the United States of America.

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

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

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

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

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

[18]  R. Friesner,et al.  Ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic catalysis. , 2005, Annual review of physical chemistry.

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

[20]  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.

[21]  J. Stewart.  Optimization of parameters for semiempirical methods I. Method , 1989 .

[22]  Hao Hu,et al.  Development and application of ab initio QM/MM methods for mechanistic simulation of reactions in solution and in enzymes. , 2009, Theochem.

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

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

[25]  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.

[26]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[27]  Thomas R. Cundari,et al.  Reviews in Computational Chemistry, Reviews in Computational Chemistry , 2000 .

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

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

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

[31]  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.

[32]  Jacopo Tomasi,et al.  Quantum Mechanical Continuum Solvation Models , 2005 .

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

[34]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[35]  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.

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

[37]  John R. Sabin,et al.  Combining quantum mechanics and molecular mechanics : some recent progresses in QM/MM methods , 2010 .

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

[39]  Hai Lin,et al.  QM/MM: What Have We Learned, Where Are We, and Where Do We Go from Here? , 2007 .