On the calculation of general response properties in subsystem density functional theory.

A recently developed subsystem formalism within time-dependent density functional theory for excitation energies [J. Chem. Phys. 126, 134116 (2007)] is extended to the calculation of general response properties such as frequency-dependent polarizabilities or optical rotatory dispersion. Furthermore, explicit expressions for oscillator and rotatory strengths are presented. Polarizabilities and optical rotation tensors of the combined system are obtained as the sum of subsystem properties, so that a natural partitioning of the polarizability or the optical rotation is facilitated. This is of great advantage for the analysis of environmental effects on molecular response properties as it allows one to distinguish changes in the ground-state electronic structure of the subsystems from cooperative excited-state effects. A computational scheme that implements the formalism is presented together with applications highlighting the differences between embedding effects included in the subsystem potential and environmental response effects. The method avoids basis-set superposition errors in the calculation of response properties of large molecules, which occur in supermolecular calculations. Possible approximations to the full response algorithm are discussed.

[1]  P. Pulay.  Improved SCF convergence acceleration , 1982 .

[2]  Johannes Neugebauer,et al.  Modeling solvent effects on electron-spin-resonance hyperfine couplings by frozen-density embedding. , 2005, The Journal of chemical physics.

[3]  T. Wesołowski.  Hydrogen-bonding-induced shifts of the excitation energies in nucleic acid bases: an interplay between electrostatic and electron density overlap effects. , 2004, Journal of the American Chemical Society.

[4]  L. H. Thomas.  The calculation of atomic fields , 1927, Mathematical Proceedings of the Cambridge Philosophical Society.

[5]  Jacques Weber,et al.  Accuracy of approximate kinetic energy functionals in the model of Kohn–Sham equations with constrained electron density: The FH⋅⋅⋅NCH complex as a test case , 1996 .

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

[7]  Piet Th. van Duijnen,et al.  A discrete solvent reaction field model for calculating molecular linear response properties in solution , 2003 .

[8]  J. Autschbach,et al.  Calculating molecular electric and magnetic properties from time-dependent density functional response theory , 2002 .

[9]  L. Jensen,et al.  Microscopic and macroscopic polarization within a combined quantum mechanics and molecular mechanics model. , 2005, The Journal of chemical physics.

[10]  H. Senn,et al.  QM/MM Methods for Biological Systems , 2006 .

[11]  L. Visscher,et al.  Calculation of nuclear magnetic resonance shieldings using frozen-density embedding. , 2006, The Journal of chemical physics.

[12]  A. Warshel,et al.  Structure/function correlations of proteins using MM, QM/MM, and related approaches: methods, concepts, pitfalls, and current progress. , 2003, Advances in protein chemistry.

[13]  M. Dulak,et al.  On the electron leak problem in orbital-free embedding calculations. , 2006, The Journal of chemical physics.

[14]  E. Baerends,et al.  Effects of complex formation on vibrational circular dichroism spectra. , 2008, The journal of physical chemistry. A.

[15]  L. Seijo,et al.  The ab initio model potential representation of the crystalline environment. Theoretical study of the local distortion on NaCl:Cu+ , 1988 .

[16]  Robert Eugene Blankenship.  Molecular mechanisms of photosynthesis , 2002 .

[17]  Alessandro Laio,et al.  A Hamiltonian electrostatic coupling scheme for hybrid Car-Parrinello molecular dynamics simulations , 2002 .

[18]  Johannes Neugebauer,et al.  Comparison of frozen-density embedding and discrete reaction field solvent models for molecular properties. , 2006, Physical chemistry chemical physics : PCCP.

[19]  M. Frisch,et al.  Hartree−Fock and Density Functional Theory ab Initio Calculation of Optical Rotation Using GIAOs: Basis Set Dependence , 2000 .

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

[21]  J. Alvarellos,et al.  Fully nonlocal kinetic energy density functionals: a proposal and a general assessment for atomic systems. , 2008, The Journal of chemical physics.

[22]  Cortona,et al.  Self-consistently determined properties of solids without band-structure calculations. , 1991, Physical review. B, Condensed matter.

[23]  P. Atkins,et al.  Molecular Quantum Mechanics , 1970 .

[24]  J. Autschbach,et al.  Calculation of optical rotation with time-periodic magnetic-field-dependent basis functions in approximate time-dependent density-functional theory. , 2005, The Journal of chemical physics.

[25]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

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

[27]  J. Autschbach,et al.  Calculation of origin-independent optical rotation tensor components in approximate time-dependent density functional theory. , 2006, The Journal of chemical physics.

[28]  J. A. Aramburu,et al.  Optical and vibrational properties of MnF64− complexes in cubic fluoroperovskites: insight through embedding calculations using Kohn–Sham equations with constrained electron density , 2006 .

[29]  Johannes Neugebauer,et al.  The merits of the frozen-density embedding scheme to model solvatochromic shifts. , 2005, The Journal of chemical physics.

[30]  T. Wesołowski.  Density Functional Theory with approximate kinetic energy functionals applied to hydrogen bonds , 1997 .

[31]  A. Klamt,et al.  Refinement and Parametrization of COSMO-RS , 1998 .

[32]  J. Neugebauer.  Induced chirality in achiral media-how theory unravels mysterious solvent effects. , 2007, Angewandte Chemie.

[33]  F. Matthias Bickelhaupt,et al.  Chemistry with ADF , 2001, J. Comput. Chem..

[34]  M. Reiher,et al.  Raman optical activity spectra of chiral transition metal complexes , 2008 .

[35]  Arieh Warshel,et al.  Simulation of enzyme reactions using valence bond force fields and other hybrid quantum/classical approaches , 1993 .

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

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

[38]  K. Morokuma,et al.  ONIOM: A Multilayered Integrated MO + MM Method for Geometry Optimizations and Single Point Energy Predictions. A Test for Diels−Alder Reactions and Pt(P(t-Bu)3)2 + H2 Oxidative Addition , 1996 .

[39]  E. Baerends,et al.  Ensuring proper short-range and asymptotic behavior of the exchange-correlation Kohn-Sham potential by modeling with a statistical average of different orbital model potentials , 2000 .

[40]  H. Eschrig,et al.  Electronic structure of solids '91 : proceedings of the 75. WE-Heraeus-Seminar and 21st Annual International Symposium on Electronic Structure of Solids, held in Gaussig (Germany), March 11-15, 1991 , 1991 .

[41]  Joachim Sauer,et al.  Combining quantum mechanics and interatomic potential functions in ab initio studies of extended systems , 2000 .

[42]  Trygve Helgaker,et al.  Optical rotation studied by density-functional and coupled-cluster methods , 2002 .

[43]  Keiji Morokuma,et al.  The IMOMO method: Integration of different levels of molecular orbital approximations for geometry optimization of large systems: Test for n‐butane conformation and SN2 reaction: RCl+Cl− , 1996 .

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

[45]  Evert Jan Baerends,et al.  A density-functional theory study of frequency-dependent polarizabilities and Van der Waals dispersion coefficients for polyatomic molecules , 1995 .

[46]  D. Truhlar,et al.  Quantum mechanical methods for enzyme kinetics. , 2003, Annual review of physical chemistry.

[47]  B. Champagne,et al.  Electrostatic interaction schemes for evaluating the polarizability of silicon clusters. , 2009, The Journal of chemical physics.

[48]  Jacques Weber,et al.  Kohn-Sham equations with constrained electron density: an iterative evaluation of the ground-state electron density of interacting molecules , 1996 .

[49]  Photophysical properties of natural light-harvesting complexes studied by subsystem density functional theory. , 2008, The journal of physical chemistry. B.

[50]  Johannes Neugebauer,et al.  An explicit quantum chemical method for modeling large solvation shells applied to aminocoumarin C151. , 2005, The journal of physical chemistry. A.

[51]  T. Wesołowski.  Application of the DFT-based embedding scheme using an explicit functional of the kinetic energy to determine the spin density of Mg+ embedded in Ne and Ar matrices , 1999 .

[52]  E. Baerends,et al.  Chiroptical properties from time-dependent density functional theory. II. Optical rotations of small to medium sized organic molecules , 2002 .

[53]  L. Visscher,et al.  Exact functional derivative of the nonadditive kinetic-energy bifunctional in the long-distance limit. , 2007, The Journal of chemical physics.

[54]  J. Perdew,et al.  Density-functional approximation for the correlation energy of the inhomogeneous electron gas. , 1986, Physical review. B, Condensed matter.

[55]  J. Tomasi,et al.  Electronic excitation energies of molecules in solution within continuum solvation models: investigating the discrepancy between state-specific and linear-response methods. , 2005, The Journal of chemical physics.

[56]  William W. Parson.  Modern Optical Spectroscopy , 2007 .

[57]  E. J. Baerends,et al.  Approximation of the exchange-correlation Kohn-Sham potential with a statistical average of different orbital model potentials. , 1999 .

[58]  Kurt V. Mikkelsen,et al.  Linear Response Properties of Liquid Water Calculated Using CC2 and CCSD within Different Molecular Mechanics Methods , 2004 .

[59]  A. Lembarki,et al.  Obtaining a gradient-corrected kinetic-energy functional from the Perdew-Wang exchange functional. , 1994, Physical review. A, Atomic, molecular, and optical physics.

[60]  Christoph R. Jacob,et al.  A flexible implementation of frozen‐density embedding for use in multilevel simulations , 2008, J. Comput. Chem..

[61]  E. Fermi.  Eine statistische Methode zur Bestimmung einiger Eigenschaften des Atoms und ihre Anwendung auf die Theorie des periodischen Systems der Elemente , 1928 .

[62]  T. Truong,et al.  Embedded density functional approach for calculations of adsorption on ionic crystals , 1996 .

[63]  J. G. Snijders,et al.  A discrete solvent reaction field model within density functional theory , 2003 .

[64]  T. Wesołowski,et al.  Orbital-free embedding applied to the calculation of induced dipole moments in CO2...X (X = He, Ne, Ar, Kr, Xe, Hg) van der Waals complexes. , 2005, The Journal of chemical physics.

[65]  Tomasz Adam Wesolowski,et al.  Generalization of the Kohn–Sham equations with constrained electron density formalism and its time‐dependent response theory formulation , 2004 .

[66]  Lucas Visscher,et al.  NMR solvent shifts of acetonitrile from frozen density embedding calculations. , 2008, The journal of physical chemistry. A.

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

[68]  Jeffrey A. Nichols,et al.  Direct atomic‐orbital‐based time‐dependent Hartree–Fock calculations of frequency‐dependent polarizabilities , 1992 .

[69]  Evert Jan Baerends,et al.  Asymptotic correction of the exchange-correlation kernel of time-dependent density functional theory for long-range charge-transfer excitations. , 2004, The Journal of chemical physics.

[70]  P. Wipf,et al.  Contribution of a solute's chiral solvent imprint to optical rotation. , 2007, Angewandte Chemie.

[71]  Kurt V. Mikkelsen,et al.  Linear response functions for coupled cluster/molecular mechanics including polarization interactions , 2003 .

[72]  S. Clima,et al.  Embedding Fragment ab Initio Model Potentials in CASSCF/CASPT2 Calculations of Doped Solids: Implementation and Applications. , 2008, Journal of chemical theory and computation.

[73]  M. Head‐Gordon,et al.  Failure of time-dependent density functional theory for long-range charge-transfer excited states: the zincbacteriochlorin-bacteriochlorin and bacteriochlorophyll-spheroidene complexes. , 2004, Journal of the American Chemical Society.

[74]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[75]  Samuel Fux,et al.  Analysis of electron density distributions from subsystem density functional theory applied to coordination bonds , 2008 .

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

[77]  M. Reiher,et al.  Topological analysis of electron densities from Kohn-Sham and subsystem density functional theory. , 2008, The Journal of chemical physics.

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

[79]  Marvin L. Cohen,et al.  Electronic structure of solids , 1984 .

[80]  T. Wesołowski,et al.  Properties of CO adsorbed in ZSM5 zeolite: Density functional theory study using the embedding scheme based on electron density partitioning , 2001 .

[81]  David J. Tozer,et al.  Relationship between long-range charge-transfer excitation energy error and integer discontinuity in Kohn–Sham theory , 2003 .

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

[83]  G. Schaftenaar,et al.  Molden: a pre- and post-processing program for molecular and electronic structures* , 2000, J. Comput. Aided Mol. Des..

[84]  Johannes Neugebauer,et al.  Exploring the ability of frozen-density embedding to model induced circular dichroism. , 2006, The journal of physical chemistry. A.

[85]  Evert Jan Baerends,et al.  Molecular calculations of excitation energies and (hyper)polarizabilities with a statistical average of orbital model exchange-correlation potentials , 2000 .

[86]  Jerzy Leszczynski,et al.  COMPUTATIONAL CHEMISTRY: Reviews of Current Trends , 2006 .

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

[88]  Johannes Neugebauer,et al.  Assessment of a simple correction for the long-range charge-transfer problem in time-dependent density-functional theory. , 2006, The Journal of chemical physics.

[89]  Johannes Neugebauer,et al.  Couplings between electronic transitions in a subsystem formulation of time-dependent density functional theory. , 2007, The Journal of chemical physics.

[90]  D. Chong.  Recent Advances in Density Functional Methods Part III , 2002 .