Fast real-time time-dependent hybrid functional calculations with the parallel transport gauge and the adaptively compressed exchange formulation

Abstract

We present a new method to accelerate real time-time dependent density functional theory (rt-TDDFT) calculations with hybrid exchange-correlation functionals. For large basis set, the computational bottleneck for large scale calculations is the application of the Fock exchange operator to the time-dependent orbitals. Our main goal is to reduce the frequency of applying the Fock exchange operator, without loss of accuracy. We achieve this by combining the recently developed parallel transport (PT) gauge formalism and the adaptively compressed exchange operator (ACE) formalism. The PT gauge yields the slowest possible dynamics among all choices of gauge. When coupled with implicit time integrators such as the Crank-Nicolson (CN) scheme, the resulting PT-CN scheme can significantly increase the time step from sub-attoseconds to 10-100 attoseconds. At each time step $t_n$, PT-CN requires the self-consistent solution of the orbitals at time $t_{n+1}$. We use ACE to delay the update of the Fock exchange operator in this nonlinear system, while maintaining the same self-consistent solution. We verify the performance of the resulting PT-CN-ACE method by computing the absorption spectrum of a benzene molecule and the response of bulk silicon systems to an ultrafast laser pulse, using the planewave basis set and the HSE functional. We report the strong and weak scaling of the PT-CN-ACE method for silicon systems ranging from 32 to 1024 atoms, with up to 2048 computational cores. Compared to standard explicit time integrators such as the 4th order Runge-Kutta method (RK4), the PT-CN-ACE can reduce the Fock exchange operator application by nearly 70 times, thus reduce the overall wall clock time time by 46 times for the system with 1024 atoms. Hence our work enables hybrid functional rt-TDDFT calculations to be routinely performed with a large basis set for the first time.

References in corpus (12)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Maximally-localized generalized Wannier functions for composite energy bands
• Optimized norm-conserving Vanderbilt pseudopotentials
• Maximally localized Wannier functions: Theory and applications
• Optimization Algorithm for the Generation of ONCV Pseudopotentials
• The shortcomings of semi-local and hybrid functionals: what we can learn from surface science studies
• Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation
• Adaptively Compressed Exchange Operator
• Efficient formalism for large scale ab initio molecular dynamics based on time-dependent density functional theory
• DGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations
• Nonlinear electronic excitations in crystalline solids using meta-generalized gradient approximation and hybrid functional in time-dependent density functional theory
• Quantum Dynamics with the Parallel Transport Gauge

Cited by in corpus (8)

• Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance
• Quantum simulation of exact electron dynamics can be more efficient than classical mean-field methods
• Propagation of maximally localized Wannier functions in real-time TDDFT
• A Snapshot of Time-Dependent Density-Functional Theory
• Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory
• INQ, a modern GPU-accelerated computational framework for (time-dependent) density functional theory
• Theory of Moment Propagation for Quantum Dynamics in Single-Particle Description
• Parallel Transport Time-Dependent Density Functional Theory Calculations with Hybrid Functional on Summit