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.

Accelerating Fock build via hybrid analytical-numerical integration

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms

Efficient and Robust Ab Initio Self-Consistent Field Acceleration Algorithm Based on a Semiempirical Model Hamiltonian

Efficient all-electron hybrid density functionals for atomistic simulations beyond 10 000 atoms

Single Field HIE FDTD Solver with Nonuniform Grid Stratagem for More Fine Structure

Efficient Implementation of the CPML in 3-D Hybrid Implicit-Explicit FDTD Method

Efficient simulations of Hartree--Fock equations by an accelerated gradient descent method

Distributed Memory, GPU Accelerated Fock Construction for Hybrid, Gaussian Basis Density Functional Theory

Fast Multipole Accelerated Scattering Matrix Method for Multiple Scattering of A Large Number of Cylinders

Advanced Techniques for High-Performance Fock Matrix Construction on GPU Clusters

Speeding-up Hybrid Functional based Ab Initio Molecular Dynamics using Multiple Time-stepping and Resonance Free Thermostat

Efficient periodic resolution-of-the-identity Hartree–Fock exchange method with k -point sampling and Gaussian basis sets

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

Interpolative separable density fitting decomposition for accelerating Hartree-Fock exchange calculations within numerical atomic orbitals

Low-rank Approximations to Accelerate Hybrid Functional Enabled Real-Time Time-Dependent Density Functional Theory Within Plane Waves

Efficient Hartree-Fock Exchange Algorithm with Coulomb Range Separation and Long-Range Density Fitting

Tucker Tensor Approach for Accelerating Fock Exchange Computations in a Real-Space Finite-Element Discretization of Generalized Kohn-Sham Density Functional Theory

Extending GPU-Accelerated Gaussian Integrals in the TeraChem Software Package to f Type Orbitals: Implementation and Applications

Efficient HF exchange evaluation through Fourier convolution in Cartesian grid for orbital-dependent density functionals

DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations