Simulating Electron Dynamics in Polarizable Environments

We propose a methodology for simulating attosecond electron dynamics in large molecular systems. Our approach is based on the combination of real time time-dependent-density-functional theory (RT-TDDFT) and polarizable Molecular Mechanics (MMpol) with the point-charge-dipole model of electrostatic i...

Full description

Saved in:
Bibliographic Details
Published in:Journal of chemical theory and computation 2017-09, Vol.13 (9), p.3985-4002
Main Authors: Wu, Xiaojing, Teuler, Jean-Marie, Cailliez, Fabien, Clavaguéra, Carine, Salahub, Dennis R, de la Lande, Aurélien
Format: Article
Language:English
Subjects:
Absorption spectra
Chemical Sciences
Computer simulation
Electric fields
Electron density
Energy absorption
Methodology
Molecular chains
Radiation effects
Time dependence
Variation
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We propose a methodology for simulating attosecond electron dynamics in large molecular systems. Our approach is based on the combination of real time time-dependent-density-functional theory (RT-TDDFT) and polarizable Molecular Mechanics (MMpol) with the point-charge-dipole model of electrostatic induction. We implemented this methodology in the software deMon2k that relies heavily on auxiliary fitted densities. In the context of RT-TDDFT/MMpol simulations, fitted densities allow the cost of the calculations to be reduced drastically on three fronts: (i) the Kohn–Sham potential, (ii) the electric field created by the (fluctuating) electron cloud which is needed in the QM/MM interaction, and (iii) the analysis of the fluctuating electron density on-the-fly. We determine conditions under which fitted densities can be used without jeopardizing the reliability of the simulations. Very encouraging results are found both for stationary and time-dependent calculations. We report absorption spectra of a dye molecule in the gas phase, in nonpolarizable water, and in polarizable water. Finally, we use the method to analyze the distance-dependent response of the environment of a peptide perturbed by an electric field. Different response mechanisms are identified. It is shown that the induction on MM sites allows excess energy to dissipate from the QM region to the environment. In this regard, the first hydration shell plays an essential role in absorbing energy. The methodology presented herein opens the possibility of simulating radiation-induced electronic phenomena in complex and extended molecular systems.
ISSN:1549-9618
1549-9626
DOI:10.1021/acs.jctc.7b00251