In this paper, we investigate the energy minimization model arising in the ensemble Kohn–Sham density functional theory for metallic systems, in which a pseudo-eigenvalue matrix and a general smearing approach are involved. We study the invariance of the energy functional and the existence of the minimizer of the ensemble Kohn–Sham model. We propose an adaptive two-parameter step size strategy and the corresponding preconditioned conjugate gradient methods to solve the energy minimization model. Under some mild but reasonable assumptions, we prove the global convergence for the gradients of the energy functional produced by our algorithms. Numerical experiments show that our algorithms are efficient, especially for large scale metallic systems. In particular, our algorithms produce convergent numerical approximations for some metallic systems, for which the traditional self-consistent field iterations fail to converge.
Mathematical Analysis and Numerical Approximations of Density Functional Theory Models for Metallic Systems / Dai, Xiaoying; de Gironcoli, Stefano; Yang, Bin; Zhou, Aihui. - In: MULTISCALE MODELING & SIMULATION. - ISSN 1540-3459. - 21:3(2023), pp. 777-803. [10.1137/22m1472103]
Mathematical Analysis and Numerical Approximations of Density Functional Theory Models for Metallic Systems
de Gironcoli, Stefano;Yang, Bin;
2023-01-01
Abstract
In this paper, we investigate the energy minimization model arising in the ensemble Kohn–Sham density functional theory for metallic systems, in which a pseudo-eigenvalue matrix and a general smearing approach are involved. We study the invariance of the energy functional and the existence of the minimizer of the ensemble Kohn–Sham model. We propose an adaptive two-parameter step size strategy and the corresponding preconditioned conjugate gradient methods to solve the energy minimization model. Under some mild but reasonable assumptions, we prove the global convergence for the gradients of the energy functional produced by our algorithms. Numerical experiments show that our algorithms are efficient, especially for large scale metallic systems. In particular, our algorithms produce convergent numerical approximations for some metallic systems, for which the traditional self-consistent field iterations fail to converge.File | Dimensione | Formato | |
---|---|---|---|
Dai-2023-Mathematical-analysis-and-numerical.pdf
accesso aperto
Tipologia:
Documento in Pre-print
Licenza:
Non specificato
Dimensione
1.34 MB
Formato
Adobe PDF
|
1.34 MB | Adobe PDF | Visualizza/Apri |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.