Polarizable continuum model with the fragment molecular orbital‐based time‐dependent density functional theory

The polarizable continuum model (PCM) for describing the solvent effect was combined with the fragment molecular orbital‐based time‐dependent density functional theory (TDDFT). Several levels of the many‐body expansion were implemented, and the importance of the many‐body contributions to the singlet‐excited states was discussed. To calibrate the accuracy, we performed a number of the model calculations using our method and the regular TDDFT in solution, applying them to phenol and polypeptides at the long‐range corrected BLYP/6‐31G* level. It was found that for systems up to 192 atoms the largest error in the excitation energy was 0.006 eV (vs. the regular TDDFT/PCM of the full system). The solvent shifts and the conformer effects were discussed, and the scaling was found to be nearly linear. Finally, we applied our method to the lowest singlet excitation of the photoactive yellow protein (PYP) in aqueous solution and determined the excitation energy to be in reasonable agreement with experiment. The excitation energy analysis provided the contributions of individual residues, and the main factors as well as their solvent shifts were determined. © 2008 Wiley Periodicals, Inc. J Comput Chem 2008

[1]  T. Nakano,et al.  Parallelized integral-direct CIS(D) calculations with multilayer fragment molecular orbital scheme , 2007 .

[2]  E. Herington,et al.  The ultra-violet absorption spectra and dissociation constants of certain phenols in aqueous solution , 1957 .

[3]  Kazuo Kitaura,et al.  Second order Møller-Plesset perturbation theory based upon the fragment molecular orbital method. , 2004, The Journal of chemical physics.

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

[5]  G. Borgstahl,et al.  1.4 A structure of photoactive yellow protein, a cytosolic photoreceptor: unusual fold, active site, and chromophore. , 1995, Biochemistry.

[6]  Kazuo Kitaura,et al.  Coupled-cluster theory based upon the fragment molecular-orbital method. , 2005, The Journal of chemical physics.

[7]  Yuto Komeiji,et al.  Change in a protein's electronic structure induced by an explicit solvent: An ab initio fragment molecular orbital study of ubiquitin , 2007, J. Comput. Chem..

[8]  Tomoo Miyahara,et al.  Symmetry-adapted-cluster/symmetry-adapted-cluster configuration interaction methodology extended to giant molecular systems: ring molecular crystals. , 2007, The Journal of chemical physics.

[9]  K. Hirao,et al.  A long-range correction scheme for generalized-gradient-approximation exchange functionals , 2001 .

[10]  Mark Earl Casida,et al.  In Recent Advances in Density-Functional Methods , 1995 .

[11]  R. Ahlrichs,et al.  Erratum: “Time-dependent density functional methods for excited state properties” [J. Chem. Phys. 117, 7433 (2002)] , 2004 .

[12]  K. Kitaura,et al.  Multilayer formulation of the fragment molecular orbital method (FMO). , 2005, The journal of physical chemistry. A.

[13]  Filipp Furche,et al.  Adiabatic time-dependent density functional methods for excited state properties , 2002 .

[14]  K. Kitaura,et al.  Time-dependent density functional theory with the multilayer fragment molecular orbital method , 2007 .

[15]  Kazuo Kitaura,et al.  CHAPTER 1 – Theoretical development of the fragment molecular orbital (FMO) method , 2006 .

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

[17]  Hui Li,et al.  The polarizable continuum model (PCM) interfaced with the fragment molecular orbital method (FMO) , 2006, J. Comput. Chem..

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

[19]  Andreas Savin,et al.  Combining long-range configuration interaction with short-range density functionals , 1997 .

[20]  N. Nakatsuji,et al.  Cluster expansion of the wavefunction. Excited states , 1978 .

[21]  Kazuo Kitaura,et al.  Pair interaction energy decomposition analysis , 2007, J. Comput. Chem..

[22]  Hui Li,et al.  Improving the efficiency and convergence of geometry optimization with the polarizable continuum model: New energy gradients and molecular surface tessellation , 2004, J. Comput. Chem..

[23]  M. E. Casida.  Time-Dependent Density Functional Response Theory for Molecules , 1995 .

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

[25]  Vincenzo Barone,et al.  A Theoretical Investigation of the Ground and Excited States of Selected Ru and Os Polypyridyl Molecular Dyes , 2002 .

[26]  M. Head‐Gordon,et al.  Long-range charge-transfer excited states in time-dependent density functional theory require non-local exchange , 2003 .

[27]  Vincenzo Barone,et al.  Time-dependent density functional theory for molecules in liquid solutions , 2001 .

[28]  Kazuo Kitaura,et al.  Molecular recognition mechanism of FK506 binding protein: An all‐electron fragment molecular orbital study , 2007, Proteins.

[29]  Kazuo Kitaura,et al.  Time-dependent density functional theory based upon the fragment molecular orbital method. , 2007, The Journal of chemical physics.

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

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

[32]  Mark S. Gordon,et al.  A new hierarchical parallelization scheme: Generalized distributed data interface (GDDI), and an application to the fragment molecular orbital method (FMO) , 2004, J. Comput. Chem..

[33]  Kazuo Kitaura,et al.  Extending the power of quantum chemistry to large systems with the fragment molecular orbital method. , 2007, The journal of physical chemistry. A.

[34]  Kazuo Kitaura,et al.  On the accuracy of the 3-body fragment molecular orbital method (FMO) applied to density functional theory , 2004 .

[35]  Takeshi Ishikawa,et al.  A fully quantum mechanical simulation study on the lowest n-π* state of hydrated formaldehyde , 2007 .

[36]  G. Scuseria,et al.  An efficient implementation of time-dependent density-functional theory for the calculation of excitation energies of large molecules , 1998 .

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

[38]  K. Hirao,et al.  A long-range-corrected time-dependent density functional theory. , 2004, The Journal of chemical physics.

[39]  K. Hirao,et al.  Long-range corrected time-dependent density functional study on fluorescence of 4,4'-dimethylaminobenzonitrile. , 2007, The Journal of chemical physics.

[40]  M. Kataoka,et al.  Photoreaction cycle of photoactive yellow protein from Ectothiorhodospira halophila studied by low-temperature spectroscopy. , 1996, Biochemistry.

[41]  Yuji Mochizuki,et al.  Configuration interaction singles method with multilayer fragment molecular orbital scheme , 2005 .

[42]  E. Gross,et al.  Density-Functional Theory for Time-Dependent Systems , 1984 .

[43]  Ying Xue,et al.  A theoretical study of solvent effects on tautomerism and electronic absorption spectra of 3‐hydroxy‐2‐mercaptopyridine and 2,3‐dihydroxypyridine , 2004, J. Comput. Chem..

[44]  K. Kitaura,et al.  Fragment molecular orbital method: an approximate computational method for large molecules , 1999 .

[45]  E. V. Gromov,et al.  Electronic structure of the PYP chromophore in its native protein environment. , 2007, Journal of the American Chemical Society.

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

[47]  Hiroshi Nakatsuji,et al.  Mechanism of color tuning in retinal protein: SAC-CI and QM/MM study , 2005 .

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

[49]  Marat Valiev,et al.  Fast electron correlation methods for molecular clusters in the ground and excited states , 2005 .

[50]  R. Ahlrichs,et al.  Treatment of electronic excitations within the adiabatic approximation of time dependent density functional theory , 1996 .

[51]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[52]  Jacopo Tomasi,et al.  Evaluation of Solvent Effects in Isotropic and Anisotropic Dielectrics and in Ionic Solutions with a Unified Integral Equation Method: Theoretical Bases, Computational Implementation, and Numerical Applications , 1997 .

[53]  Gerrit Groenhof,et al.  Signal transduction in the photoactive yellow protein. I. Photon absorption and the isomerization of the chromophore , 2002, Proteins.

[54]  Kimihiko Hirao,et al.  Excited state geometry optimizations by analytical energy gradient of long-range corrected time-dependent density functional theory. , 2006, The Journal of chemical physics.