Damped Response Theory in Combination with Polarizable Environments: The Polarizable Embedding Complex Polarization Propagator Method.

We present a combination of the polarizable embedding (PE) scheme with the complex polarization propagator (CPP) method with the aim of calculating response properties including relaxation for large and complex systems. This new approach, termed PE-CPP, will benefit from the highly advanced description of the environmental electrostatic potential and polarization in the PE method as well as the treatment of near-resonant effects in the CPP approach. The PE-CPP model has been implemented in a Kohn-Sham density functional theory approach, and we present pilot calculations exemplifying the implementation for the UV/vis and carbon K-edge X-ray absorption spectra of the protein plastocyanin. Furthermore, technical details associated with a PE-CPP calculation are discussed.

[1]  Frank Neese,et al.  QM/MM calculations with DFT for taking into account protein effects on the EPR and optical spectra of metalloproteins. Plastocyanin as a case study , 2006, J. Comput. Chem..

[2]  Robert K Szilagyi,et al.  Electronic structures of metal sites in proteins and models: contributions to function in blue copper proteins. , 2004, Chemical reviews.

[3]  Jacob Kongsted,et al.  The multi-configuration self-consistent field method within a polarizable embedded framework. , 2013, The Journal of chemical physics.

[4]  J. Rivail,et al.  A QM/MM study on the spinach plastocyanin: Redox properties and absorption spectra , 2012 .

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

[6]  Trygve Helgaker,et al.  Recent advances in wave function-based methods of molecular-property calculations. , 2012, Chemical reviews.

[7]  A. Gewirth,et al.  Electronic structure of plastocyanin: excited state spectral features , 1988 .

[8]  Luca Frediani,et al.  The Dalton quantum chemistry program system , 2013, Wiley interdisciplinary reviews. Computational molecular science.

[9]  Poul Jørgensen,et al.  Response functions from Fourier component variational perturbation theory applied to a time-averaged quasienergy , 1998 .

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

[11]  D. M. Bishop,et al.  Near-resonant absorption in the time-dependent self-consistent field and multiconfigurational self-consistent field approximations , 2001 .

[12]  Roland Lindh,et al.  Accuracy of distributed multipoles and polarizabilities: Comparison between the LoProp and MpProp models , 2007, J. Comput. Chem..

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

[14]  N. Besley,et al.  Modelling the spectroscopy and dynamics of plastocyanin. , 2010, Physical chemistry chemical physics : PCCP.

[15]  H. Bartunik,et al.  Accuracy and precision in protein structure analysis: restrained least-squares refinement of the structure of poplar plastocyanin at 1.33 A resolution. , 1992, Acta crystallographica. Section B, Structural science.

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

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

[18]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[19]  Ulf Ekström,et al.  X-ray absorption spectra from the resonant-convergent first-order polarization propagator approach , 2006 .

[20]  Patrick Norman,et al.  A perspective on nonresonant and resonant electronic response theory for time-dependent molecular properties. , 2011, Physical chemistry chemical physics : PCCP.

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

[22]  Jan H. Jensen,et al.  PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. , 2011, Journal of chemical theory and computation.

[23]  Michael Dolg,et al.  Energy‐adjusted ab initio pseudopotentials for the first row transition elements , 1987 .

[24]  Paweł Sałek,et al.  Density-functional theory of linear and nonlinear time-dependent molecular properties , 2002 .

[25]  Edward I. Solomon,et al.  Recent advances in understanding blue copper proteins , 2011 .

[26]  Poul Jørgensen,et al.  Quasienergy formulation of damped response theory. , 2009, The Journal of chemical physics.

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

[28]  N. Handy,et al.  A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP) , 2004 .

[29]  Pär Söderhjelm,et al.  How accurate can a force field become? A polarizable multipole model combined with fragment-wise quantum-mechanical calculations. , 2009, The journal of physical chemistry. A.

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

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

[32]  Timothy Clark,et al.  Efficient diffuse function‐augmented basis sets for anion calculations. III. The 3‐21+G basis set for first‐row elements, Li–F , 1983 .

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

[34]  Ulf Ekström,et al.  Relativistic four-component static-exchange approximation for core-excitation processes in molecules , 2006 .

[35]  N. Melikechi,et al.  Determinations of EF Σg+1 X Σg+1 transition frequencies in H2, D2, and HD , 2006 .

[36]  H. Ågren,et al.  Self-interaction-corrected time-dependent density-functional-theory calculations of x-ray-absorption spectra , 2007 .

[37]  Poul Jørgensen,et al.  On the Efficiency of Algorithms for Solving Hartree-Fock and Kohn-Sham Response Equations. , 2011, Journal of chemical theory and computation.

[38]  FRANCESCO AQUILANTE,et al.  MOLCAS 7: The Next Generation , 2010, J. Comput. Chem..

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

[40]  Lasse Jensen,et al.  A discrete interaction model/quantum mechanical method to describe the interaction of metal nanoparticles and molecular absorption. , 2011, The Journal of chemical physics.

[41]  Richard A. Friesner,et al.  Integrated Modeling Program, Applied Chemical Theory (IMPACT) , 2005, J. Comput. Chem..

[42]  B. Roos,et al.  Theoretical study of the electronic spectrum of plastocyanin , 1997 .

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

[44]  Trygve Helgaker,et al.  Molecular Electronic-Structure Theory: Helgaker/Molecular Electronic-Structure Theory , 2000 .

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

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

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

[48]  John Z. H. Zhang,et al.  Molecular fractionation with conjugate caps for full quantum mechanical calculation of protein-molecule interaction energy , 2003 .