Volume 14, Issue 2
Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System

Tiao Lu, Wei Cai, Jianguo Xin & Yinglong Guo

Commun. Comput. Phys., 14 (2013), pp. 276-300.

Published online: 2014-08

Preview Full PDF 105 616
Export citation
  • Abstract

In the first of a series of papers, we will study a discontinuous Galerkin (DG) framework for many electron quantum systems. The salient feature of this framework is the flexibility of using hybrid physics-based local orbitals and accuracy-guaranteed piecewise polynomial basis in representing the Hamiltonian of the many body system. Such a flexibility is made possible by using the discontinuous Galerkin method to approximate the Hamiltonian matrix elements with proper constructions of numerical DG fluxes at the finite element interfaces. In this paper, we will apply the DG method to the density matrix minimization formulation, a popular approach in the density functional theory of many body Schrödinger equations. The density matrix minimization is to find the minima of the total energy, expressed as a functional of the density matrix ρ(r,r ′ ), approximated by the proposed enriched basis, together with two constraints of idempotency and electric neutrality. The idempotency will be handled with the McWeeny’s purification while the neutrality is enforced by imposing the number of electrons with a penalty method. A conjugate gradient method (a PolakRibiere variant) is used to solve the minimization problem. Finally, the linear-scaling algorithm and the advantage of using the local orbital enriched finite element basis in the DG approximations are verified by studying examples of one dimensional lattice model systems.


  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • References
  • Hide All
    View All

  • BibTex
  • RIS
  • TXT
@Article{CiCP-14-276, author = {Tiao Lu, Wei Cai, Jianguo Xin and Yinglong Guo}, title = {Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System}, journal = {Communications in Computational Physics}, year = {2014}, volume = {14}, number = {2}, pages = {276--300}, abstract = {

In the first of a series of papers, we will study a discontinuous Galerkin (DG) framework for many electron quantum systems. The salient feature of this framework is the flexibility of using hybrid physics-based local orbitals and accuracy-guaranteed piecewise polynomial basis in representing the Hamiltonian of the many body system. Such a flexibility is made possible by using the discontinuous Galerkin method to approximate the Hamiltonian matrix elements with proper constructions of numerical DG fluxes at the finite element interfaces. In this paper, we will apply the DG method to the density matrix minimization formulation, a popular approach in the density functional theory of many body Schrödinger equations. The density matrix minimization is to find the minima of the total energy, expressed as a functional of the density matrix ρ(r,r ′ ), approximated by the proposed enriched basis, together with two constraints of idempotency and electric neutrality. The idempotency will be handled with the McWeeny’s purification while the neutrality is enforced by imposing the number of electrons with a penalty method. A conjugate gradient method (a PolakRibiere variant) is used to solve the minimization problem. Finally, the linear-scaling algorithm and the advantage of using the local orbital enriched finite element basis in the DG approximations are verified by studying examples of one dimensional lattice model systems.


}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.290212.240812a}, url = {http://global-sci.org/intro/article_detail/cicp/7160.html} }
TY - JOUR T1 - Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System AU - Tiao Lu, Wei Cai, Jianguo Xin & Yinglong Guo JO - Communications in Computational Physics VL - 2 SP - 276 EP - 300 PY - 2014 DA - 2014/08 SN - 14 DO - http://dor.org/10.4208/cicp.290212.240812a UR - https://global-sci.org/intro/cicp/7160.html KW - AB -

In the first of a series of papers, we will study a discontinuous Galerkin (DG) framework for many electron quantum systems. The salient feature of this framework is the flexibility of using hybrid physics-based local orbitals and accuracy-guaranteed piecewise polynomial basis in representing the Hamiltonian of the many body system. Such a flexibility is made possible by using the discontinuous Galerkin method to approximate the Hamiltonian matrix elements with proper constructions of numerical DG fluxes at the finite element interfaces. In this paper, we will apply the DG method to the density matrix minimization formulation, a popular approach in the density functional theory of many body Schrödinger equations. The density matrix minimization is to find the minima of the total energy, expressed as a functional of the density matrix ρ(r,r ′ ), approximated by the proposed enriched basis, together with two constraints of idempotency and electric neutrality. The idempotency will be handled with the McWeeny’s purification while the neutrality is enforced by imposing the number of electrons with a penalty method. A conjugate gradient method (a PolakRibiere variant) is used to solve the minimization problem. Finally, the linear-scaling algorithm and the advantage of using the local orbital enriched finite element basis in the DG approximations are verified by studying examples of one dimensional lattice model systems.


Tiao Lu, Wei Cai, Jianguo Xin & Yinglong Guo. (1970). Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System. Communications in Computational Physics. 14 (2). 276-300. doi:10.4208/cicp.290212.240812a
Copy to clipboard
The citation has been copied to your clipboard