Skip to main content Accessibility help
×
Home
Hostname: page-component-857cdc78dc-zmt5p Total loading time: 0.456 Render date: 2022-05-16T06:55:17.400Z Has data issue: true Feature Flags: { "shouldUseShareProductTool": true, "shouldUseHypothesis": true, "isUnsiloEnabled": true, "useRatesEcommerce": false, "useNewApi": true }

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

Published online by Cambridge University Press:  03 June 2015

Tiao Lu*
Affiliation:
HEDPS & CAPT, LMAM & School of Mathematical Sciences, Peking University, Beijing 100871, P.R. China
Wei Cai*
Affiliation:
INS, Shanghai Jiaotong University, Shanghai 200240, P.R. China Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28223-0001, USA
Jianguo Xin*
Affiliation:
Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28223-0001, USA
Yinglong Guo*
Affiliation:
School of Mathematical Sciences, Peking University, Beijing 100871, P.R. China
*
Get access

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 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.

Type
Research Article
Copyright
Copyright © Global Science Press Limited 2013

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Abgrall, R. and Shu, C.-W.Development of residual distribution schemes for the discontinuous Galerkin method: The scalar case with linear elements. Commun. Comput. Phys., 5(2-4):376390, 2009.Google Scholar
[2]Akin, J.E.The generation of elements with singularities. Int. J. Numer. Meth. Engng., 10:12491259, 1976.CrossRefGoogle Scholar
[3]Babusška, I. and Rosenzweig, B.A finite element scheme for domains with corners. Numer. Math., 20:121, 1972.CrossRefGoogle Scholar
[4]Chelikowsky, J.R., Troullier, N., and Saad, Y.Finite difference-pseudopotential method: Electronic structure calculations without a basis. Phys. Rev. Lett., 72:12401243, 1994.CrossRefGoogle ScholarPubMed
[5]Cockburn, B. and Shu, C.-W.The local discontinuous galerkin method for time-dependent convection-diffusion systems. SIAM Journal on Numerical Analysis, 35(6):24402463, 1998.CrossRefGoogle Scholar
[6]Galli, Giulia. Linear scaling methods for electronic structure calculations and quantum molecular dynamics simulations. Current Opinion in Solid State and Materials Science, 1(6):864874, 1996.CrossRefGoogle Scholar
[7]Galli, Giulia and Parrinello, Michele. Large scale electronic structure calculations. Phys. Rev. Lett., 69(24):35473550, Dec 1992.Google ScholarPubMed
[8]García-Cervera, C. J.A remark on an efficient real space method for orbital-free density-functional theory. Commun. Comput. Phys., 1(1):15, 2006.Google Scholar
[9]García-Cervera, C. J.An efficient real space method for orbital-free density-functional theory. Commun. Comput. Phys., 2(2):334357, 2007.Google Scholar
[10]García-Cervera, C.J., Lu, J., Xuan, Y., and E, W.Linear-scaling subspace-iteration algorithm with optimally localized nonorthogonal wave functions for kohn-sham density functional theory. Phys. Rev. B, 79:115110, 2009.CrossRefGoogle Scholar
[11]Goedecker, S.Low complexity algorithms for electronic structure calculations. J. Comput. Phys., 118:261268, 1995.CrossRefGoogle Scholar
[12]Goedecker, S.Linear scaling electronic structure methods. Rev. Mod. Phys., 71(4):10851123, Jul 1999.CrossRefGoogle Scholar
[13]Goedecker, S. and Colombo, L.Efficient linear scaling algorithm for tight-binding molecular dynamics. Phys. Rev. Lett., 73:122125, Jul 1994.CrossRefGoogle ScholarPubMed
[14]Hernaández, E. and Gillan, M. J.Self-consistent first-principles technique with linear scaling. Phys. Rev. B, 51(15):1015710160, Apr 1995.CrossRefGoogle Scholar
[15]Kohn, W. and Sham, L. J.Self-consistent equations including exchange and correlation effects. Phys. Rev., 140:A1133A1138, Nov 1965.CrossRefGoogle Scholar
[16]Li, X.-P., Nunes, R. W., and Vanderbilt, David. Density-matrix electronic-structure method with linear system-size scaling. Phys. Rev. B, 47(16):1089110894, Apr 1993.CrossRefGoogle ScholarPubMed
[17]Lin, L., Lu, J., Ying, L., and E, W.Adaptive local basis set for kohn-sham density functional theory in a discontinuous galerkin framework I: Total energy calculation. submitted to J. Comput. Phys., 2011.Google Scholar
[18]Lu, T. and Cai, W.A fourier spectral-discontinuous galerkin method for time-dependent 3-d schrödinger-poisson equations with discontinuous potentials. J. Comput. Appl. Math., 220:588614, October 2008.CrossRefGoogle Scholar
[19]Luenberger, D. and Ye, Y.Linear and Nonlinear Programming. Springer, third edition, 2008.Google Scholar
[20]Martin, R.Electronic Structure: Basic Theory and Practical Methods. Cambridge Univ. Press, Cambridge, 2004.CrossRefGoogle Scholar
[21]Mauri, F., Galli, G., and Car, R.Orbital formulation for electronic-structure calculations with linear system-size scaling. Phys. Rev. B, 47(15):99739976, Apr 1993.CrossRefGoogle ScholarPubMed
[22]Weeny, R. Mc. The density matrix in self-conssitent field theory. I. iterative constructon of the density matrix. Proc. Roy. Soc. Lond. A, 235:496509, 1956.Google Scholar
[23]Slater, J.C.Wave functions in a periodic potential. Phys. Rev., 51:846851, 1937.CrossRefGoogle Scholar
[24]Troullier, N. and Martins, José Luriaas. Efficient pseudopotentials for plane-wave calculations. Phys. Rev. B, 43:19932006, Jan 1991.CrossRefGoogle ScholarPubMed
[25]Wang, W. and Shu, C.-W.The wkb local discontinuous Galerkin method for the simulation of Schrödinger equation in a resonant tunneling diode. J. Scientific Computing, 40:360374, 2009.CrossRefGoogle Scholar
[26]Yang, W.Direct calculation of electron density in density-functional theory. Phys. Rev. Lett., 66:14381441, Mar 1991.CrossRefGoogle ScholarPubMed
2
Cited by

Save article to Kindle

To save this article to your Kindle, first ensure coreplatform@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about saving to your Kindle.

Note you can select to save to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

Find out more about the Kindle Personal Document Service.

Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System
Available formats
×

Save article to Dropbox

To save this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you used this feature, you will be asked to authorise Cambridge Core to connect with your Dropbox account. Find out more about saving content to Dropbox.

Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System
Available formats
×

Save article to Google Drive

To save this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you used this feature, you will be asked to authorise Cambridge Core to connect with your Google Drive account. Find out more about saving content to Google Drive.

Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System
Available formats
×
×

Reply to: Submit a response

Please enter your response.

Your details

Please enter a valid email address.

Conflicting interests

Do you have any conflicting interests? *