skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

This content will become publicly available on September 23, 2020

Title: 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 semidiscrete (discrete in space but continuous in time) error estimate on the time-dependent Kohn-Sham orbitals, to construct an efficient 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 special reduced order quadrature to render the overlap matrix diagonal, thereby 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, wherein 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, for both pseudopotential and all-electron calculations.more » Further, for similar level of accuracy, we obtain significant computational savings by our approach as compared to state-of-the-art finite-difference methods. Here, 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.« less

Authors:
 [1];  [1]
  1. Univ. of Michigan, Ann Arbor, MI (United States)
Publication Date:
Research Org.:
Univ. of Michigan, Ann Arbor, MI (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22)
OSTI Identifier:
1564082
Grant/Contract Number:  
SC0017380
Resource Type:
Accepted Manuscript
Journal Name:
Physical Review B
Additional Journal Information:
Journal Volume: 100; Journal Issue: 11; Journal ID: ISSN 2469-9950
Publisher:
American Physical Society (APS)
Country of Publication:
United States
Language:
English
Subject:
36 MATERIALS SCIENCE

Citation Formats

Kanungo, Bikash, and Gavini, Vikram. Real time time-dependent density functional theory using higher order finite-element methods. United States: N. p., 2019. Web. doi:10.1103/PhysRevB.100.115148.
Kanungo, Bikash, & Gavini, Vikram. Real time time-dependent density functional theory using higher order finite-element methods. United States. doi:10.1103/PhysRevB.100.115148.
Kanungo, Bikash, and Gavini, Vikram. Mon . "Real time time-dependent density functional theory using higher order finite-element methods". United States. doi:10.1103/PhysRevB.100.115148.
@article{osti_1564082,
title = {Real time time-dependent density functional theory using higher order finite-element methods},
author = {Kanungo, Bikash and Gavini, Vikram},
abstractNote = {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 semidiscrete (discrete in space but continuous in time) error estimate on the time-dependent Kohn-Sham orbitals, to construct an efficient 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 special reduced order quadrature to render the overlap matrix diagonal, thereby 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, wherein 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, for both pseudopotential and all-electron calculations. Further, for similar level of accuracy, we obtain significant computational savings by our approach as compared to state-of-the-art finite-difference methods. Here, 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.},
doi = {10.1103/PhysRevB.100.115148},
journal = {Physical Review B},
number = 11,
volume = 100,
place = {United States},
year = {2019},
month = {9}
}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on September 23, 2020
Publisher's Version of Record

Save / Share:

Works referenced in this record:

NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations
journal, September 2010

  • Valiev, M.; Bylaska, E. J.; Govind, N.
  • Computer Physics Communications, Vol. 181, Issue 9, p. 1477-1489
  • DOI: 10.1016/j.cpc.2010.04.018

Self-Consistent Equations Including Exchange and Correlation Effects
journal, November 1965


Ground State of the Electron Gas by a Stochastic Method
journal, August 1980