Numerical Solution of the Kohn-Sham Equation by Finite Element Methods with an Adaptive Mesh Redistribution Technique

被引:25
作者
Bao, Gang [1 ,2 ]
Hu, Guanghui [2 ]
Liu, Di [2 ]
机构
[1] Zhejiang Univ, Dept Math, Hangzhou 310027, Peoples R China
[2] Michigan State Univ, Dept Math, E Lansing, MI 48824 USA
基金
美国国家科学基金会;
关键词
Adaptive mesh redistribution; Harmonic map; Finite element method; Kohn-Sham equation; Density functional theory; DENSITY-FUNCTIONAL THEORY; HARMONIC MAPS; ENERGY; FRAMEWORK;
D O I
10.1007/s10915-012-9636-1
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
A mesh redistribution method is introduced to solve the Kohn-Sham equation. The standard linear finite element space is employed for the spatial discretization, and the self-consistent field iteration scheme is adopted for the derived nonlinear generalized eigenvalue problem. A mesh redistribution technique is used to optimize the distribution of the mesh grids according to wavefunctions obtained from the self-consistent iterations. After the mesh redistribution, important regions in the domain such as the vicinity of the nucleus, as well as the bonding between the atoms, may be resolved more effectively. Consequently, more accurate numerical results are obtained without increasing the number of mesh grids. Numerical experiments confirm the effectiveness and reliability of our method for a wide range of problems. The accuracy and efficiency of the method are also illustrated through examples.
引用
收藏
页码:372 / 391
页数:20
相关论文
共 42 条
[11]   Moving Finite Element Simulations for Reaction-Diffusion Systems [J].
Hu, Guanghui ;
Qiao, Zhonghua ;
Tang, Tao .
ADVANCES IN APPLIED MATHEMATICS AND MECHANICS, 2012, 4 (03) :365-381
[12]   Simulating finger phenomena in porous media with a moving finite element method [J].
Hu, Guanghui ;
Zegeling, Paul Andries .
JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 230 (08) :3249-3263
[13]  
Huang W., 2010, APPL MATH SCI
[14]   Preconditioners and Electron Density Optimization in Orbital-Free Density Functional Theory [J].
Hung, Linda ;
Huang, Chen ;
Carter, Emily A. .
COMMUNICATIONS IN COMPUTATIONAL PHYSICS, 2012, 12 (01) :135-161
[15]  
Johnson R.D., 2011, NIST Computational Chemistry Comparison and Benchmark Database, release, V15b
[16]   EFFICACIOUS FORM FOR MODEL PSEUDOPOTENTIALS [J].
KLEINMAN, L ;
BYLANDER, DM .
PHYSICAL REVIEW LETTERS, 1982, 48 (20) :1425-1428
[17]   Block locally optimal preconditioned eigenvalue Xolvers (BLOPEX) in hypre and PETSc [J].
Knyazev, A. V. ;
Argentati, M. E. ;
Lashuk, I. ;
Ovtchinnikov, E. E. .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2007, 29 (05) :2224-2239
[18]   SELF-CONSISTENT EQUATIONS INCLUDING EXCHANGE AND CORRELATION EFFECTS [J].
KOHN, W ;
SHAM, LJ .
PHYSICAL REVIEW, 1965, 140 (4A) :1133-&
[19]   Local-density-functional calculations of the energy of atoms [J].
Kotochigova, S ;
Levine, ZH ;
Shirley, EL ;
Stiles, MD ;
Clark, CW .
PHYSICAL REVIEW A, 1997, 55 (01) :191-199
[20]   Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set [J].
Kresse, G ;
Furthmuller, J .
PHYSICAL REVIEW B, 1996, 54 (16) :11169-11186