Real-Time TD-DFT with Classical Ion Dynamics: Methodology and Applications.

We present a method for real-time propagation of electronic wave functions, within time-dependent density functional theory (RT-TDDFT), coupled to ionic motion through mean-field classical dynamics. The goal of our method is to treat large systems and complex processes, in particular photocatalytic reactions and electron transfer events on surfaces and thin films. Due to the complexity of these processes, computational approaches are needed to provide insight into the underlying physical mechanisms and are therefore crucial for the rational design of new materials. Because of the short time step required for electron propagation (of order ∼10 attoseconds), these simulations are computationally very demanding. Our methodology is based on numerical atomic-orbital-basis sets for computational efficiency. In the computational package, to which we refer as TDAP-2.0 (Time-evolving Deterministic Atom Propagator), we have implemented a number of important features and analysis tools for more accurate and efficient treatment of large, complex systems and time scales that reach into a fraction of a picosecond. We showcase the capabilities of our method using four different examples: (i) photodissociation into radicals of opposite spin, (ii) hydrogen adsorption on aluminum surfaces, (iii) optical absorption of spin-polarized organic molecule containing a metal ion, and (iv) electron transfer in a prototypical dye-sensitized solar cell.

[1]  T. F. Heinz,et al.  Electron transport in TiO 2 probed by THz time-domain spectroscopy , 2004 .

[2]  R. Schmidt,et al.  Non-adiabatic quantum molecular dynamics: basic formalism and case study , 1996 .

[3]  Muratahan Aykol,et al.  Local environment dependent GGA + U method for accurate thermochemistry of transition metal compounds , 2014 .

[4]  Michel Dupuis,et al.  Localized Electronic States from Surface Hydroxyls and Polarons in TiO2(110) , 2009 .

[5]  Richard M. Martin,et al.  Calculation of the optical response of atomic clusters using time-dependent density functional theory and local orbitals , 2002 .

[6]  Seifert,et al.  Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon. , 1995, Physical review. B, Condensed matter.

[7]  M. Frisch,et al.  Ab initio Ehrenfest dynamics. , 2005, The Journal of chemical physics.

[8]  Rainer Eichberger,et al.  Role of molecular anchor groups in molecule-to-semiconductor electron transfer. , 2006, The journal of physical chemistry. B.

[9]  J. Tully.  Mixed quantum–classical dynamics , 1998 .

[10]  Peter Pulay,et al.  Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules , 1969 .

[11]  C. Humphreys,et al.  Electron-energy-loss spectra and the structural stability of nickel oxide: An LSDA+U study , 1998 .

[12]  Sawatzky,et al.  Local-density functional and on-site correlations: The electronic structure of La2CuO4 and LaCuO3. , 1994, Physical review. B, Condensed matter.

[13]  Jussi Enkovaara,et al.  Localized surface plasmon resonance in silver nanoparticles: Atomistic first-principles time-dependent density-functional theory calculations , 2015, 1503.07234.

[14]  I. I. Mazin,et al.  Correlated metals and the LDA+U method , 2002, cond-mat/0206548.

[15]  A. Lichtenstein,et al.  First-principles calculations of electronic structure and spectra of strongly correlated systems: the LDA+U method , 1997 .

[16]  E. Pehlke,et al.  Molecular dynamics of nonadiabatic processes at surfaces : Chemisorption of H/Al(111) , 2006 .

[17]  Mykhaylo Krykunov,et al.  The formulation of a self-consistent constricted variational density functional theory for the description of excited states , 2011 .

[18]  Nikos L. Doltsinis,et al.  FIRST PRINCIPLES MOLECULAR DYNAMICS INVOLVING EXCITED STATES AND NONADIABATIC TRANSITIONS , 2002 .

[19]  Z. Lan,et al.  An On-the-Fly Surface-Hopping Program JADE for Nonadiabatic Molecular Dynamics of Polyatomic Systems: Implementation and Applications. , 2015, Journal of chemical theory and computation.

[20]  Peter M. Oppeneer,et al.  GGA + U modeling of structural, electronic, and magnetic properties of iron porphyrin-type molecules , 2008 .

[21]  Shane R. Yost,et al.  Assessment of the ΔSCF density functional theory approach for electronic excitations in organic dyes. , 2011, The Journal of chemical physics.

[22]  D. M. Newns.  Self-Consistent Model of Hydrogen Chemisorption , 1969 .

[23]  Ralph Ernstorfer,et al.  Ultrafast interfacial electron transfer from the excited state of anchored molecules into a semiconductor , 2007 .

[24]  Michael Grätzel,et al.  Porphyrin-Sensitized Solar Cells with Cobalt (II/III)–Based Redox Electrolyte Exceed 12 Percent Efficiency , 2011, Science.

[25]  Lars Nordstrom,et al.  Multipole decomposition of LDA+U energy and its application to actinide compounds , 2009, 0904.0742.

[26]  Michael Grätzel,et al.  Photoelectrochemical cells , 2001, Nature.

[27]  Stefano de Gironcoli,et al.  Linear response approach to the calculation of the effective interaction parameters in the LDA + U method , 2004, cond-mat/0405160.

[28]  Walter Thiel,et al.  Ground States of Molecules. 38. The MNDO Method. Approximations and Parameters , 1977 .

[29]  Tian Lan,et al.  Phonon anharmonicity of rutile TiO_2 studied by Raman spectrometry and molecular dynamics simulations , 2012 .

[30]  Efthimios Kaxiras,et al.  Mechanisms for ultrafast nonradiative relaxation in electronically excited eumelanin constituents. , 2008, Biophysical journal.

[31]  Roger B. Sidje,et al.  Expokit: a software package for computing matrix exponentials , 1998, TOMS.

[32]  Kiyoyuki Terakura,et al.  Is Hund's Second Rule Responsible for the Orbital Magnetism in Solids? , 1998 .

[33]  Alain Delgado,et al.  Charge separation dynamics and opto-electronic properties of a diaminoterephthalate- C 60 diad , 2014 .

[34]  F. Matthias Bickelhaupt,et al.  Voronoi deformation density (VDD) charges: Assessment of the Mulliken, Bader, Hirshfeld, Weinhold, and VDD methods for charge analysis , 2004, J. Comput. Chem..

[35]  Mats Persson,et al.  Dynamics of the spin transition in the adsorption of hydrogen atoms on metals , 2004 .

[36]  E. Kaxiras,et al.  Dynamics of the Photogenerated Hole at the Rutile TiO2(110)/Water Interface: A Nonadiabatic Simulation Study , 2014 .

[37]  Nicola Marzari,et al.  First-principles photoemission spectroscopy and orbital tomography in molecules from koopmans-compliant functionals. , 2015, Physical review letters.

[38]  M. Batzill,et al.  Why is anatase a better photocatalyst than rutile? - Model studies on epitaxial TiO2 films , 2014, Scientific Reports.

[39]  E. Pehlke,et al.  Ab initio simulation of the spin transition during chemisorption: H/Al(111). , 2006, Physical review letters.

[40]  T. Park,et al.  Unitary quantum time evolution by iterative Lanczos reduction , 1986 .

[41]  B. Lundqvist,et al.  Exchange and correlation in atoms, molecules, and solids by the spin-density-functional formalism , 1976 .

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

[43]  Chi-Cheng Lee,et al.  Dynamical linear response of TDDFT with LDA + U functional: Strongly hybridized Frenkel excitons in NiO , 2010 .

[44]  Andreas Görling,et al.  Density-functional theory beyond the Hohenberg-Kohn theorem , 1999 .

[45]  Peter M W Gill,et al.  Self-consistent field calculations of excited states using the maximum overlap method (MOM). , 2008, The journal of physical chemistry. A.

[46]  Yagi,et al.  Electronic conduction above 4 K of slightly reduced oxygen-deficient rutile TiO2-x. , 1996, Physical review. B, Condensed matter.

[47]  Sándor Suhai,et al.  Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties , 1998 .

[48]  Daniela Kohen,et al.  Model studies of nonadiabatic dynamics , 1998 .

[49]  A. Akimov,et al.  Advanced Capabilities of the PYXAID Program: Integration Schemes, Decoherence Effects, Multiexcitonic States, and Field-Matter Interaction. , 2014, Journal of chemical theory and computation.

[50]  Efthimios Kaxiras,et al.  Design of Dye Acceptors for Photovoltaics from First-Principles Calculations , 2011 .

[51]  D. Sánchez-Portal,et al.  The SIESTA method for ab initio order-N materials simulation , 2001, cond-mat/0111138.

[52]  Jack Ballantyne,et al.  A Blue Spectral Shift of the Hemoglobin Soret Band Correlates with the Age (Time Since Deposition) of Dried Bloodstains , 2010, PloS one.

[53]  Angel Rubio,et al.  Propagators for the time-dependent Kohn-Sham equations. , 2004, The Journal of chemical physics.

[54]  Niranjan Govind,et al.  Modeling Fast Electron Dynamics with Real-Time Time-Dependent Density Functional Theory: Application to Small Molecules and Chromophores. , 2011, Journal of chemical theory and computation.

[55]  A. Akimov,et al.  The PYXAID Program for Non-Adiabatic Molecular Dynamics in Condensed Matter Systems. , 2013, Journal of chemical theory and computation.

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

[57]  D. Bird,et al.  Electron–hole pair creation by atoms incident on a metal surface , 2003, 0909.5501.

[58]  Heather J Kulik,et al.  Systematic study of first-row transition-metal diatomic molecules: a self-consistent DFT+U approach. , 2010, The Journal of chemical physics.

[59]  E. Kaxiras,et al.  Anatomy of the Photochemical Reaction: Excited-State Dynamics Reveals the C-H Acidity Mechanism of Methoxy Photo-oxidation on Titania. , 2015, The journal of physical chemistry letters.

[60]  Xavier Andrade,et al.  Modified Ehrenfest Formalism for Efficient Large-Scale ab initio Molecular Dynamics. , 2008, Journal of chemical theory and computation.

[61]  Rüdiger Schmidt,et al.  Non-adiabatic quantum molecular dynamics: General formalism and case study H2+ in strong laser fields , 2003 .

[62]  C. Schütte,et al.  Quantum‐classical molecular dynamics as an approximation to full quantum dynamics , 1996 .

[63]  Dominik Marx,et al.  Charge localization dynamics induced by oxygen vacancies on the TiO₂(110) surface. , 2010, Physical review letters.

[64]  Marco Buongiorno Nardelli,et al.  Reformulation of DFT + U as a pseudohybrid hubbard density functional for accelerated materials discovery , 2015 .

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

[66]  H Wende,et al.  Substrate-induced magnetic ordering and switching of iron porphyrin molecules. , 2007, Nature materials.

[67]  Ville Havu,et al.  Nonadiabatic Ehrenfest molecular dynamics within the projector augmented-wave method. , 2011, The Journal of chemical physics.

[68]  F. L. Hirshfeld.  Bonded-atom fragments for describing molecular charge densities , 1977 .

[69]  M. Dewar,et al.  Ground States of Molecules. 38. The MNDO Method. Approximations and Parameters , 1977 .

[70]  Efthimios Kaxiras,et al.  Real-time, local basis-set implementation of time-dependent density functional theory for excited state dynamics simulations. , 2008, The Journal of chemical physics.

[71]  Sampsa Juhana Riikonen.  First principles study of nanostructured surface reconstructions induced by the deposition of metals on vicinal si (111) surfaces , 2007 .

[72]  M. Engelund,et al.  Delta self-consistent field method to obtain potential energy surfaces of excited molecules on surfaces , 2008 .

[73]  Mark E. Casida,et al.  Time-dependent density-functional theory for molecules and molecular solids , 2009 .

[74]  David L. Beveridge,et al.  Approximate molecular orbital theory , 1970 .

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

[76]  L. Reining,et al.  Electronic excitations: density-functional versus many-body Green's-function approaches , 2002 .

[77]  V. Anisimov,et al.  Band theory and Mott insulators: Hubbard U instead of Stoner I. , 1991, Physical review. B, Condensed matter.

[78]  J. Zaanen,et al.  Density-functional theory and strong interactions: Orbital ordering in Mott-Hubbard insulators. , 1995, Physical review. B, Condensed matter.

[79]  Xavier Andrade,et al.  Time-dependent density-functional theory in massively parallel computer architectures: the octopus project , 2012, Journal of physics. Condensed matter : an Institute of Physics journal.

[80]  Hans Lischka,et al.  The on-the-fly surface-hopping program system Newton-X: Application to ab initio simulation of the nonadiabatic photodynamics of benchmark systems , 2007 .

[81]  P. Marquetand,et al.  SHARC: ab Initio Molecular Dynamics with Surface Hopping in the Adiabatic Representation Including Arbitrary Couplings. , 2011, Journal of chemical theory and computation.

[82]  K. Reuter,et al.  Excited-state potential-energy surfaces of metal-adsorbed organic molecules from linear expansion Δ-self-consistent field density-functional theory (ΔSCF-DFT). , 2013, The Journal of chemical physics.

[83]  J. Tully.  Molecular dynamics with electronic transitions , 1990 .

[84]  Arvi Rauk,et al.  On the calculation of multiplet energies by the hartree-fock-slater method , 1977 .

[85]  U. V. Barth,et al.  Local-density theory of multiplet structure , 1979 .