Accuracy and computational efficiency of real-time subspace propagation schemes for the time-dependent density functional theory.

Time-dependent Density Functional Theory (TDDFT) has become successful for its balance of economy and accuracy. However, the application of TDDFT to large systems or long time scales remains computationally prohibitively expensive. In this paper, we investigate the numerical stability and accuracy of two subspace propagation methods to solve the time-dependent Kohn-Sham equations with finite and periodic boundary conditions. The bases considered are the Lánczos basis and the adiabatic eigenbasis. The results are compared to a benchmark fourth-order Taylor expansion of the time propagator. Our results show that it is possible to use larger time steps with the subspace methods, leading to computational speedups by a factor of 2-3 over Taylor propagation. Accuracy is found to be maintained for certain energy regimes and small time scales.

[1]  André D. Bandrauk,et al.  EXPONENTIAL PROPAGATORS (INTEGRATORS) FOR THE TIME-DEPENDENT SCHRÖDINGER EQUATION , 2013 .

[2]  R. Levis,et al.  Coulomb explosion of large polyatomic molecules assisted by nonadiabatic charge localization. , 2004, Physical review letters.

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

[4]  Jonathan S. Parker,et al.  Numerical integration of the time-dependent Schrödinger equation for laser-driven helium , 1998 .

[5]  C. Lanczos.  An iteration method for the solution of the eigenvalue problem of linear differential and integral operators , 1950 .

[6]  Carsten A. Ullrich,et al.  A Brief Compendium of Time-Dependent Density Functional Theory , 2013, 1305.1388.

[7]  R. Baer,et al.  Time-dependent density-functional studies of the D2 coulomb explosion. , 2006, The journal of physical chemistry. A.

[8]  J. A. Alonso,et al.  Interaction of the Charged Deuterium Cluster D3+ with Femtosecond Laser Pulses† , 2007 .

[9]  K. Varga,et al.  Electron-ion dynamics in laser-assisted desorption of hydrogen atoms from H-Si(111) surface , 2011 .

[10]  G. Schatz,et al.  From Discrete Electronic States to Plasmons: TDDFT Optical Absorption Properties of Agn(n= 10, 20, 35, 56, 84, 120) Tetrahedral Clusters , 2008 .

[11]  Christine M Aikens,et al.  Electronic structure and TDDFT optical absorption spectra of silver nanorods. , 2009, The journal of physical chemistry. A.

[12]  Uwe Manthe,et al.  Quantum calculations of thermal rate constants and reaction probabilities: H2+CN→H+HCN , 1998 .

[13]  K. Varga,et al.  Time-dependent density-functional study of field emission from tipped carbon nanotubes , 2009 .

[14]  Lin-Wang Wang,et al.  Efficient real-time time-dependent density functional theory method and its application to a collision of an ion with a 2D material. , 2015, Physical review letters.

[15]  U. Manthe,et al.  Dissociation and predissociation on coupled electronic potential energy surfaces: A three‐dimensional wave packet dynamical study , 1991 .

[16]  Sullivan,et al.  Real-space multigrid-based approach to large-scale electronic structure calculations. , 1996, Physical review. B, Condensed matter.

[17]  K. Varga,et al.  Time-dependent density functional study of field emission from nanotubes composed of C, BN, SiC, Si, and GaN , 2011, Nanotechnology.

[18]  T. Akama,et al.  Development of efficient time-evolution method based on three-term recurrence relation. , 2015, The Journal of chemical physics.

[19]  Eric Polizzi,et al.  Spectral-based propagation schemes for time-dependent quantum systems with application to carbon nanotubes , 2010 .

[20]  M. Feit,et al.  Solution of the Schrödinger equation by a spectral method , 1982 .

[21]  S. X. Hu,et al.  Parallel solver for the time-dependent linear and nonlinear Schrödinger equation. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[22]  H. Zeng,et al.  Strong-field dissociative double ionization of acetylene. , 2014, Physical review letters.

[23]  S. Gray,et al.  Quantum dynamics with real wave packets, including application to three-dimensional (J=0) D + H2 → HD + H reactive scattering , 1998 .

[24]  R. Kosloff,et al.  A fourier method solution for the time dependent Schrödinger equation as a tool in molecular dynamics , 1983 .

[25]  S. Chu,et al.  Time-dependent density-functional theory for molecular processes in strong fields: Study of multiphoton processes and dynamical response of individual valence electrons of N2 in intense laser fields , 2001 .

[26]  Zhigang Sun,et al.  Communication: An exact short-time solver for the time-dependent Schrödinger equation. , 2011, The Journal of chemical physics.

[27]  H. Rabitz,et al.  Control of quantum phenomena: past, present and future , 2009, 0912.5121.

[28]  Zhigang Sun,et al.  Comparison of second-order split operator and Chebyshev propagator in wave packet based state-to-state reactive scattering calculations. , 2009, The Journal of chemical physics.

[29]  Sixth-order factorization of the evolution operator for time-dependent potentials. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  B. Piraux,et al.  Explicit schemes for time propagating many-body wave functions , 2014, 1401.6318.

[31]  K. Varga,et al.  Time-dependent density-functional study of the ionization and fragmentation ofC2H2andH2by strong circularly polarized laser pulses , 2015 .

[32]  S. Pantelides,et al.  Simulation of high-energy ion collisions with graphene fragments , 2012 .

[33]  J. Feist,et al.  Attosecond streaking of correlated two-electron transitions in helium. , 2011, Physical Review Letters.

[34]  André D. Bandrauk,et al.  Complex integration steps in decomposition of quantum exponential evolution operators , 2006 .

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

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

[37]  A. Zunger,et al.  Self-interaction correction to density-functional approximations for many-electron systems , 1981 .

[38]  A. Baltuska,et al.  High energy proton ejection from hydrocarbon molecules driven by highly efficient field ionization. , 2011 .

[39]  A. Krasheninnikov,et al.  Role of electronic excitations in ion collisions with carbon nanostructures. , 2006, Physical review letters.

[40]  Martins,et al.  Efficient pseudopotentials for plane-wave calculations. , 1991, Physical review. B, Condensed matter.

[41]  George F. Bertsch,et al.  Time-dependent local-density approximation in real time , 1996 .

[42]  Y. Saad,et al.  Finite-difference-pseudopotential method: Electronic structure calculations without a basis. , 1994, Physical review letters.

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

[44]  A. Baltuska,et al.  Selective control over fragmentation reactions in polyatomic molecules using impulsive laser alignment. , 2014, Physical review letters.

[45]  M. E. Casida,et al.  Progress in time-dependent density-functional theory. , 2011, Annual review of physical chemistry.

[46]  Spectral difference Lanczos method for efficient time propagation in quantum control theory. , 2004, The Journal of chemical physics.

[47]  A. Baltuska,et al.  Strong laser-pulse-driven ionization and Coulomb explosion of hydrocarbon molecules , 2012 .

[48]  Claude Leforestier,et al.  A comparison of different propagation schemes for the time dependent Schro¨dinger equation , 1991 .

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

[50]  E. Kaxiras,et al.  A TDDFT study of the optical response of DNA bases, base pairs, and their tautomers in the gas phase. , 2005, The journal of physical chemistry. A.

[51]  Hongbing Jiang,et al.  Selective steering of molecular multiple dissociative channels with strong few-cycle laser pulses. , 2011, Physical review letters.

[52]  P. Ehrenfest.  Bemerkung über die angenäherte Gültigkeit der klassischen Mechanik innerhalb der Quantenmechanik , 1927 .

[53]  H. Tal-Ezer,et al.  An accurate and efficient scheme for propagating the time dependent Schrödinger equation , 1984 .

[54]  Xavier Andrade,et al.  Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems. , 2015, Physical chemistry chemical physics : PCCP.

[55]  A. Baltuska,et al.  Attosecond-recollision-controlled selective fragmentation of polyatomic molecules. , 2012, Physical review letters.

[56]  Lin-wang Wang,et al.  Nonadiabatic molecular dynamics simulation for carrier transport in a pentathiophene butyric acid monolayer , 2013 .

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

[58]  M. Beck,et al.  Dynamical effects in the interaction of ion beams with solids. , 2008, Physical review letters.

[59]  First-principles study of field emission from carbon nanotubes and graphene nanoribbons , 2011 .

[60]  M. Jungen,et al.  Lanczos wave packet propagation on coupled potential energy surfaces: the three body predissociation of rotating D3 and H3 3 2 A ′ ?> ( 2 sa 1 ′ ?> ) , 2015 .

[61]  K. Varga,et al.  First-principles time-dependent simulation of laser assisted desorption of hydrogen atoms from H–Si(111) surface , 2011 .

[62]  Hua Guo,et al.  A single Lanczos propagation method for calculating transition amplitudes , 1999 .

[63]  Angel Rubio,et al.  Real-space, real-time method for the dielectric function , 2000 .