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.
2023
21
3
777
803
https://arxiv.org/abs/2201.07035
Dai, Xiaoying; de Gironcoli, Stefano; Yang, Bin; Zhou, Aihui
File in questo prodotto:
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.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.11767/137591
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 0
  • ???jsp.display-item.citation.isi??? 0
social impact