Time-dependent non-equilibrium dielectric response in QM/continuum approaches.

The Polarizable Continuum Models (PCMs) are some of the most inexpensive yet successful methods for including the effects of solvation in quantum-mechanical calculations of molecular systems. However, when applied to the electronic excitation process, these methods are restricted to dichotomously assuming either that the solvent has completely equilibrated with the excited solute charge density (infinite-time limit), or that it retains the configuration that was in equilibrium with the solute prior to excitation (zero-time limit). This renders the traditional PCMs inappropriate for resolving time-dependent solvent effects on non-equilibrium solute electron dynamics like those implicated in the instants following photoexcitation of a solvated molecular species. To extend the existing methods to this non-equilibrium regime, we herein derive and apply a new formalism for a general time-dependent continuum embedding method designed to be propagated alongside the solute's electronic degrees of freedom in the time domain. Given the frequency-dependent dielectric constant of the solvent, an equation of motion for the dielectric polarization is derived within the PCM framework and numerically integrated simultaneously with the time-dependent Hartree fock/density functional theory equations. Results for small molecular systems show the anticipated dipole quenching and electronic state dephasing/relaxation resulting from out-of-phase charge fluctuations in the dielectric and embedded quantum system.