Real-time time-dependent density functional theory using higher order finite element methods

Abstract

We present a computationally efficient approach to solve the time-dependent Kohn-Sham equations in real-time using higher-order finite-element spatial discretization, applicable to both pseudopotential and all-electron calculations. To this end, we develop an a priori mesh adaption technique, based on the semi-discrete (discrete in space but continuous in time) error estimate on the time-dependent Kohn-Sham orbitals, to construct a close to optimal finite-element discretization. Subsequently, we obtain the full-discrete error estimate to guide our choice of the time-step. We employ spectral finite-elements along with Gauss-Legendre-Lobatto quadrature to render the overlap matrix diagonal, simplifying the inversion of the overlap matrix that features in the evaluation of the discrete time-evolution operator. We use the second-order Magnus operator as the time-evolution operator in all our calculations. Furthermore, the action of the discrete Magnus operator, expressed as exponential of a matrix, on the Kohn-Sham orbitals is obtained efficiently through an adaptive Lanczos iteration. We observe close to optimal rates of convergence of the dipole moment with respect to spatial and temporal discretization, for both pseudopotential and all-electron calculations. We demonstrate a staggering 100-fold reduction in the computational time afforded by higher-order finite-elements over linear finite-elements. We present comparative studies, in terms of accuracy and efficiency, of our approach against finite-difference based discretization for pseudopotential calculations, and demonstrate significant computational savings when compared to the finite-difference method. We also demonstrate the competence of higher-order finite-elements for all-electron benchmark systems. Lastly, we observe good parallel scalability of the proposed method on many hundreds of processors.

References in corpus (23)

• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• A real-space grid implementation of the Projector Augmented Wave method
• The Magnus expansion and some of its applications
• Daubechies wavelets as a basis set for density functional pseudopotential calculations
• Time-dependent quantum transport: A practical scheme using density functional theory
• Localized surface plasmon resonance in silver nanoparticles: Atomistic first-principles time-dependent density-functional theory calculations
• DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations
• Excitations in time-dependent density-functional theory
• Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration
• Time-dependent electron localization function
• Higher-order adaptive finite-element methods for Kohn-Sham density functional theory
• Time-dependent quantum transport: an exact formulation based on TDDFT
• Center of mass and relative motion in time dependent density functional theory
• Real-space local polynomial basis for solid-state electronic-structure calculations: A finite-element approach
• Initial-state dependence in time-dependent density functional theory
• All-electron density functional theory and time-dependent density functional theory with high-order finite elements
• Large-scale all-electron density functional theory calculations using an enriched finite element basis
• A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization
• Higher-order adaptive finite-element methods for orbital-free density functional theory
• Configurational forces in electronic structure calculations using Kohn-Sham density functional theory
• Wavelet-Based Linear-Response Time-Dependent Density-Functional Theory
• Spectrum-splitting approach for Fermi-operator expansion in all-electron Kohn-Sham DFT calculations

Cited by in corpus (4)

• Roadmap on Electronic Structure Codes in the Exascale Era
• DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization
• Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements
• Ionic forces and stress tensor in all-electron DFT calculations using enriched finite element basis