@Article{EAJAM-13-2,
author = {Hu, Guanghui and Wang, Ting and Zhou, Jie},
title = {A Linearized Structure-Preserving Numerical Scheme for a Gradient Flow Model of the Kohn-Sham Density Functional Theory},
journal = {East Asian Journal on Applied Mathematics},
year = {2023},
volume = {13},
number = {2},
pages = {299--319},
abstract = {<p style="text-align: justify;">Dai et al. [Multiscale Model. Simul. 18 (2020)] proposed a gradient flow
model and a numerical scheme for ground state calculations in Kohn-Sham density functional theory. It is a feature that orthonormality of all wave functions can be preserved
automatically during the simulation which makes such a method attractive towards simulations for large scale systems. In this paper, two extensions are proposed for further
improving the efficiency of the method. The first one is a linearization of the original
nonlinear scheme. It is shown analytically that both the orthonormality of wave functions and the decay of the total energy can be preserved well by this linear scheme,
while a significant acceleration can be observed from the numerical experiments due
to the removal of an iteration process in the nonlinear scheme. The second one is the
introduction of the adaptivity in the algorithm both temporally and spatially — i.e. an $h$-adaptive mesh method is employed to control the total amount of mesh grids, and
an adaptive stop criterion in time propagation process is designed based on an observation that total energy always decays much faster at the beginning. Plenty of numerical
experiments successfully demonstrate effectiveness of our method.</p>},
issn = {2079-7370},
doi = {https://doi.org/10.4208/eajam.2022-134.081022},
url = {https://global-sci.com/article/82426/a-linearized-structure-preserving-numerical-scheme-for-a-gradient-flow-model-of-the-kohn-sham-density-functional-theory}
}