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 

机构地区：[1]HEDPS&CAPT,LMAM&School of Mathematical Sciences,Peking University,Beijing 100871,P.R.China [2]INS,Shanghai Jiaotong University,Shanghai 200240,P.R.China [3]Department of Mathematics and Statistics,University of North Carolina at Charlotte,Charlotte,NC 28223-0001,USA [4]School of Mathematical Sciences,Peking University,Beijing 100871,P.R.China

出　　处：《Communications in Computational Physics》2013年第7期276-300,共25页计算物理通讯（英文）

基　　金：support of U.S.Army Research Office(grant number W911NF-11-1-0364);support of NSFC(grant number 11011130029)and of SRF for ROCS,SEM.

摘　　要：In the first of a series of papers,wewill 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 Schrodinger 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 theMcWeeny’s purification while the neutrality is enforced by imposing the number of electrons with a penalty method.A conjugate gradient method(a Polak-Ribiere 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.

关 键 词：Density functional theory density matrix minimization discontinuous Galerkin method linear scaling method 

分 类 号：O17[理学—数学]