Equation of motion for the solvent polarization apparent charges in the polarizable continuum model: application to real-time TDDFT.

When a solute charge density is evolving in time, e.g., due to an external perturbation, the solvent reaction field also becomes time-dependent, in a nontrivial way due to the delayed response of the solvent polarization rooted in its frequency-dependent dielectric constant. In polarizable continuum models, the time-dependent reaction field is represented by time-dependent apparent surface charges. Here, we derive general expressions for such charges. In particular, for all the main flavors of PCM, including IEF-PCM, we show how the frequency-dependent dielectric function terms can be singled-out in diagonal matrices, most convenient for Fourier transforming. For spherical cavities such formulation highlights the relation with multipolar solvation models and, when applied to the related context of metal nanoparticles, discloses a direct connection with multipolar plasmons. Using the Debye dielectric function, we derive a simple equation of motion for the apparent charges, free from system history. Such an equation has been coupled to real time time-dependent density functional theory (RT-TDDFT), to simulate the time evolution of the solute density rigorously accounting for the delayed solvent reaction field. The presented method seamlessly encompasses previous nonequilibrium approaches limited to an instantaneous solute potential change (e.g., a sudden electronic excitation), does not require additional assumptions besides the basic PCM's, and is not limited to iterative inversion procedures. Numerical examples are given, showing the importance of accounting for the delayed solvent-response effects.

[1]  D. Chandler,et al.  Dielectric solvation dynamics of molecules of arbitrary shape and charge distribution , 1998 .

[2]  Troy Van Voorhis,et al.  Simulating molecular conductance using real-time density functional theory , 2006 .

[3]  V. Barone,et al.  Quantum Calculation of Molecular Energies and Energy Gradients in Solution by a Conductor Solvent Model , 1998 .

[4]  J. Tomasi,et al.  Quantum mechanical calculations coupled with a dynamical continuum model for the description of dielectric relaxation: Time dependent Stokes shift of coumarin C153 in polar solvents , 2003 .

[5]  A. Klamt.  Conductor-like Screening Model for Real Solvents: A New Approach to the Quantitative Calculation of Solvation Phenomena , 1995 .

[6]  J. Tomasi,et al.  A time-dependent polarizable continuum model: theory and application. , 2005, The Journal of chemical physics.

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

[8]  H. Berendsen,et al.  Inclusion of reaction fields in molecular dynamics. Application to liquid water , 1978 .

[9]  Jacopo Tomasi,et al.  Enhanced response properties of a chromophore physisorbed on a metal particle , 2001 .

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

[11]  Daniel M. Chipman,et al.  Charge penetration in dielectric models of solvation , 1997 .

[12]  J. Hynes,et al.  Time-dependent fluorescence solvent shifts, dielectric friction, and nonequilibrium solvation in polar solvents , 1985 .

[13]  Feizhi Ding,et al.  Modeling ultrafast solvated electronic dynamics using time-dependent density functional theory and polarizable continuum model. , 2012, The journal of physical chemistry. A.

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

[15]  Jacopo Tomasi,et al.  Molecular properties in solution described with a continuum solvation model , 2002 .

[16]  Dmitri A Romanov,et al.  A time-dependent Hartree-Fock approach for studying the electronic optical response of molecules in intense fields. , 2005, Physical chemistry chemical physics : PCCP.

[17]  U. Hohenester,et al.  Interaction of Single Molecules With Metallic Nanoparticles , 2008, IEEE Journal of Selected Topics in Quantum Electronics.

[18]  M. Vener,et al.  Advanced Continuum Approaches for Treating Time Correlation Functions. The Role of Solute Shape and Solvent Structure , 1999 .

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

[20]  L. Onsager.  Electric Moments of Molecules in Liquids , 1936 .

[21]  Jacopo Tomasi,et al.  On the Calculation of Local Field Factors for Microscopic Static Hyperpolarizabilities of Molecules in Solution with the Aid of Quantum-Mechanical Methods , 1998 .

[22]  S. Pipolo,et al.  The cavity electromagnetic field within the polarizable continuum model of solvation: An application to the real-time time dependent density functional theory , 2014 .

[23]  S. Pipolo,et al.  The cavity electromagnetic field within the polarizable continuum model of solvation. , 2014, The Journal of chemical physics.

[24]  W. Magnus.  On the exponential solution of differential equations for a linear operator , 1954 .

[25]  M. Vener,et al.  An advanced dielectric continuum approach for treating solvation effects: Time correlation functions. I. Local treatment , 1998 .

[26]  Giovanni Scalmani,et al.  New developments in the polarizable continuum model for quantum mechanical and classical calculations on molecules in solution , 2002 .

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

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

[29]  Mukamel,et al.  Time-dependent density-matrix functional in Liouville space and the optical response of many-electron systems. , 1995, Physical review. A, Atomic, molecular, and optical physics.

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

[31]  Fernando Casas,et al.  Improved High Order Integrators Based on the Magnus Expansion , 2000 .

[32]  Jean-Louis Rivail,et al.  A quantum chemical approach to dielectric solvent effects in molecular liquids , 1976 .

[33]  R. Fuchs,et al.  Theory of the optical properties of ionic crystal cubes , 1975 .

[34]  J. Hynes,et al.  Dynamical polar solvent effects on solution reactions: A simple continuum model , 1982 .

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

[36]  C. Cramer,et al.  A universal approach to solvation modeling. , 2008, Accounts of chemical research.

[37]  Chao-Ping Hsu,et al.  Time-Dependent Stokes Shift and Its Calculation from Solvent Dielectric Dispersion Data , 1997 .

[38]  S. Mukamel,et al.  Density-matrix representation of nonadiabatic couplings in time-dependent density functional (TDDFT) theories , 2000 .

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

[40]  Jacopo Tomasi,et al.  Formation and relaxation of excited states in solution: a new time dependent polarizable continuum model based on time dependent density functional theory. , 2006, The Journal of chemical physics.

[41]  C. Isborn,et al.  Time-dependent density functional theory Ehrenfest dynamics: collisions between atomic oxygen and graphite clusters. , 2007, The Journal of chemical physics.

[42]  P. D. Nguyen,et al.  Solvated First-Principles Excited-State Charge-Transfer Dynamics with Time-Dependent Polarizable Continuum Model and Solvent Dielectric Relaxation , 2012 .

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

[44]  J. Ponder,et al.  Calculation of the reaction field due to off-center point multipoles , 1997 .

[45]  D. Chipman.  Reaction field treatment of charge penetration , 2000 .

[46]  J. Tomasi,et al.  Analytical Hartree–Fock calculation of the dynamical polarizabilities α, β, and γ of molecules in solution , 1996 .