The multi-configuration self-consistent field method within a polarizable embedded framework.

We present a detailed derivation of Multi-Configuration Self-Consistent Field (MCSCF) optimization and linear response equations within the polarizable embedding scheme: PE-MCSCF. The MCSCF model enables a proper description of multiconfigurational effects in reaction paths, spin systems, excited states, and other properties which cannot be described adequately with current implementations of polarizable embedding in density functional or coupled cluster theories. In the PE-MCSCF scheme the environment surrounding the central quantum mechanical system is represented by distributed multipole moments and anisotropic dipole-dipole polarizabilities. The PE-MCSCF model has been implemented in DALTON. As a preliminary application, the low lying valence states of acetone and uracil in water has been calculated using Complete Active Space Self-Consistent Field (CASSCF) wave functions. The dynamics of the water environment have been simulated using a series of snapshots generated from classical Molecular Dynamics. The calculated shifts from gas-phase to water display between good and excellent correlation with experiment and previous calculations. As an illustration of another area of potential applications we present calculations of electronic transitions in the transition metal complex, [Fe(NO)(CN)5](2-) in a micro-solvated environment. This system is highly multiconfigurational and the influence of solvation is significant.

[1]  H. Ågren,et al.  Modeling the structure and absorption spectra of stilbazolium merocyanine in polar and nonpolar solvents using hybrid QM/MM techniques. , 2010, The journal of physical chemistry. B.

[2]  K. Mikkelsen,et al.  Linear response properties for solvated molecules described by a combined multiconfigurational self-consistent-field/molecular mechanics model , 2002 .

[3]  Mark S. Gordon,et al.  Solvent-induced frequency shifts: configuration interaction singles combined with the effective fragment potential method. , 2010, The journal of physical chemistry. A.

[4]  H. Gray,et al.  Electronic Structure of Nitroprusside Ion , 1965 .

[5]  Pär Söderhjelm,et al.  Protein Influence on Electronic Spectra Modeled by Multipoles and Polarizabilities. , 2009, Journal of chemical theory and computation.

[6]  Ignacio Tinoco,et al.  Vapor Spectra and Heats of Vaporization of Some Purine and Pyrimidine Bases1 , 1965 .

[7]  Jacob Kongsted,et al.  Color modeling of protein optical probes. , 2012, Physical chemistry chemical physics : PCCP.

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

[9]  Jiali Gao,et al.  Explicit polarization (X-Pol) potential using ab initio molecular orbital theory and density functional theory. , 2009, The journal of physical chemistry. A.

[10]  Jacob Kongsted,et al.  Solvatochromic Shifts in Uracil: A Combined MD-QM/MM Study. , 2010, Journal of chemical theory and computation.

[11]  J. Olsen,et al.  Linear and nonlinear response functions for an exact state and for an MCSCF state , 1985 .

[12]  Walter Thiel,et al.  Toward accurate barriers for enzymatic reactions: QM/MM case study on p-hydroxybenzoate hydroxylase. , 2008, The Journal of chemical physics.

[13]  Julien Toulouse,et al.  On the universality of the long-/short-range separation in multiconfigurational density-functional theory. , 2007, The Journal of chemical physics.

[14]  J. Olsen,et al.  Quadratic response functions for a multiconfigurational self‐consistent field wave function , 1992 .

[15]  Donald G Truhlar,et al.  The variational explicit polarization potential and analytical first derivative of energy: Towards a next generation force field. , 2008, The Journal of chemical physics.

[16]  Jacob Kongsted,et al.  Solvation Effects on Electronic Transitions: Exploring the Performance of Advanced Solvent Potentials in Polarizable Embedding Calculations. , 2011, Journal of chemical theory and computation.

[17]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

[18]  V. Buss,et al.  Origin and consequences of steric strain in the rhodopsin binding pocket. , 2006, Biochemistry.

[19]  K. Mikkelsen,et al.  Coupled cluster calculation of the n --> pi* electronic transition of acetone in aqueous solution. , 2005, The journal of physical chemistry. A.

[20]  A. Schäfer,et al.  Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr , 1994 .

[21]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[22]  G. Langley,et al.  Photolysis of the [Fe(CN)5(NO)]2– ion in water and poly(vinyl alcohol) films: evidence for cyano radical, cyanide ion and nitric oxide loss and redox pathways , 1995 .

[23]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[24]  Kurt V. Mikkelsen,et al.  A multiconfiguration self‐consistent reaction field response method , 1994 .

[25]  Mati Karelson,et al.  Theoretical treatment of solvent effects on electronic spectroscopy , 1992 .

[26]  W. Thiel,et al.  Influence of the LOV domain on low-lying excited states of flavin: a combined quantum-mechanics/molecular-mechanics investigation. , 2009, The journal of physical chemistry. B.

[27]  Mark S. Gordon,et al.  The Effective Fragment Potential Method: A QM-Based MM Approach to Modeling Environmental Effects in Chemistry , 2001 .

[28]  L. Slipchenko,et al.  Solvation of the excited states of chromophores in polarizable environment: orbital relaxation versus polarization. , 2010, The journal of physical chemistry. A.

[29]  K. Ruud,et al.  Multiconfigurational self-consistent field linear response for the polarizable continuum model: Theory and application to ground and excited-state polarizabilities of para-nitroaniline in solution , 2003 .

[30]  Haibo Yu,et al.  Accounting for polarization in molecular simulation , 2005, Comput. Phys. Commun..

[31]  Kurt V. Mikkelsen,et al.  The combined multiconfigurational self-consistent-field/molecular mechanics wave function approach , 2001 .

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

[33]  Jógvan Magnus Haugaard Olsen,et al.  PERI-CC2: A Polarizable Embedded RI-CC2 Method. , 2012, Journal of chemical theory and computation.

[34]  Walter Thiel,et al.  QM/MM study of the absorption spectra of DsRed.M1 chromophores , 2010, J. Comput. Chem..

[35]  A. Warshel,et al.  Frozen density functional approach for ab initio calculations of solvated molecules , 1993 .

[36]  Jacopo Tomasi,et al.  A second-order, quadratically convergent multiconfigurational self-consistent field polarizable continuum model for equilibrium and nonequilibrium solvation , 2002 .

[37]  J. Kongsted,et al.  How to model solvent effects on molecular properties using quantum chemistry? Insights from polarizable discrete or continuum solvation models. , 2007, The journal of physical chemistry. A.

[38]  Masaaki Fujii,et al.  Electronic spectra of uracil in a supersonic jet , 1986 .

[39]  Mark S. Gordon,et al.  An effective fragment method for modeling solvent effects in quantum mechanical calculations , 1996 .

[40]  J. Olsen,et al.  TIME-DEPENDENT RESPONSE THEORY WITH APPLICATIONS TO SELF-CONSISTENT FIELD AND MULTICONFIGURATIONAL SELF-CONSISTENT FIELD WAVE FUNCTIONS , 1995 .

[41]  Roland Lindh,et al.  Local properties of quantum chemical systems: the LoProp approach. , 2004, The Journal of chemical physics.

[42]  William L Jorgensen,et al.  Special Issue on Polarization. , 2007, Journal of chemical theory and computation.

[43]  E. Dalgaard.  Time‐dependent multiconfigurational Hartree–Fock theory , 1980 .

[44]  Kazuo Kitaura,et al.  Multiconfiguration self-consistent-field theory based upon the fragment molecular orbital method. , 2005, The Journal of chemical physics.

[45]  J. Kussmann,et al.  Linear‐Scaling Methods in Quantum Chemistry , 2007 .

[46]  J. R. Carl,et al.  Atom dipole interaction model for molecular polarizability. Application to polyatomic molecules and determination of atom polarizabilities , 1972 .

[47]  H. Ågren,et al.  Modeling Solvatochromism of Nile Red in Water , 2011 .

[48]  Hao Hu,et al.  Free energies of chemical reactions in solution and in enzymes with ab initio quantum mechanics/molecular mechanics methods. , 2008, Annual review of physical chemistry.

[49]  A. Buckingham.  Permanent and Induced Molecular Moments and Long‐Range Intermolecular Forces , 2007 .

[50]  K. Morokuma,et al.  Investigation of the S0S1 excitation in bacteriorhodopsin with the ONIOM(MO:MM) hybrid method , 2003 .

[51]  Jógvan Magnus Haugaard Olsen,et al.  Amyloid Fibril-Induced Structural and Spectral Modifications in the Thioflavin-T Optical Probe. , 2013, The journal of physical chemistry letters.

[52]  Adrian J Mulholland,et al.  High-accuracy computation of reaction barriers in enzymes. , 2006, Angewandte Chemie.

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

[54]  Kenneth Ruud,et al.  Solvatochromic shift of phenol blue in water from a combined Car-Parrinello molecular dynamics hybrid quantum mechanics-molecular mechanics and ZINDO approach. , 2010, The Journal of chemical physics.

[55]  I. Renge.  Solvent dependence of n-pi* absorption in acetone. , 2009, The journal of physical chemistry. A.

[56]  Jacob Kongsted,et al.  The polarizable embedding coupled cluster method. , 2011, Journal of Chemical Physics.

[57]  Hans Ågren,et al.  A direct, restricted-step, second-order MC SCF program for large scale ab initio calculations , 1986 .

[58]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[59]  W. Thiel,et al.  Hybrid Models for Combined Quantum Mechanical and Molecular Mechanical Approaches , 1996 .

[60]  Jan H. Jensen,et al.  Chapter 10 The Effective Fragment Potential: A General Method for Predicting Intermolecular Interactions , 2007 .

[61]  Vincenzo Barone,et al.  Solvent effect on vertical electronic transitions by the polarizable continuum model , 2000 .

[62]  Jun-Ya Hasegawa,et al.  Theoretical Studies on the Color-Tuning Mechanism in Retinal Proteins. , 2007, Journal of chemical theory and computation.

[63]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[64]  K. Kistler,et al.  Solvatochromic shifts of uracil and cytosine using a combined multireference configuration interaction/molecular dynamics approach and the fragment molecular orbital method. , 2009, The journal of physical chemistry. A.

[65]  Helmut Grubmüller,et al.  Ultrafast deactivation of an excited cytosine-guanine base pair in DNA. , 2007, Journal of the American Chemical Society.

[66]  Marat Valiev,et al.  Hybrid coupled cluster and molecular dynamics approach: application to the excitation spectrum of cytosine in the native DNA environment. , 2006, The Journal of chemical physics.

[67]  W. Hauswirth,et al.  Fluorescence of the Purine and Pyrimidine Bases of the Nucleic Acids in Neutral Aqueous Solution at 300�K , 1971, Science.

[68]  K. Bravaya,et al.  An opsin shift in rhodopsin: retinal S0-S1 excitation in protein, in solution, and in the gas phase. , 2007, Journal of the American Chemical Society.

[69]  Walter Thiel,et al.  QM/MM calculation of solvent effects on absorption spectra of guanine , 2010, J. Comput. Chem..

[70]  Jacopo Tomasi,et al.  Excited states and solvatochromic shifts within a nonequilibrium solvation approach: A new formulation of the integral equation formalism method at the self-consistent field, configuration interaction, and multiconfiguration self-consistent field level , 1998 .

[71]  Kirk A Peterson,et al.  Systematically convergent basis sets for transition metals. I. All-electron correlation consistent basis sets for the 3d elements Sc-Zn. , 2005, The Journal of chemical physics.

[72]  Adrian J. Mulholland,et al.  Investigations of enzyme-catalysed reactions with combined quantum mechanics/molecular mechanics (QM/MM) methods , 2010 .

[73]  Thom Vreven,et al.  Combining Quantum Mechanics Methods with Molecular Mechanics Methods in ONIOM. , 2006, Journal of chemical theory and computation.

[74]  Trygve Helgaker,et al.  A multiconfigurational self‐consistent reaction‐field method , 1988 .

[75]  Donald G. Truhlar,et al.  Two-response-time model based on CM2/INDO/S2 electrostatic potentials for the dielectric polarization component of solvatochromic shifts on vertical excitation energies , 2000 .

[76]  Jacob Kongsted,et al.  Excited States in Solution through Polarizable Embedding , 2010 .

[77]  M. Gordon,et al.  Solvent effects on optical properties of molecules: a combined time-dependent density functional theory/effective fragment potential approach. , 2008, The Journal of chemical physics.

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

[79]  Feliu Maseras,et al.  IMOMM: A new integrated ab initio + molecular mechanics geometry optimization scheme of equilibrium structures and transition states , 1995, J. Comput. Chem..