A multilevel correction adaptive finite element method for Kohn-Sham equation

In this paper, an adaptive finite element method is proposed for solving Kohn–Sham equation with the multilevel correction technique. In the method, the Kohn–Sham equation is solved on a fixed and appropriately coarse mesh with the finite element method in which the finite element space is kept improving by solving the derived boundary value problems on a series of adaptively and successively refined meshes. A main feature of the method is that solving large scale Kohn–Sham system is avoided effectively, and solving the derived boundary value problems can be handled efficiently by classical methods such as the multigrid method. Hence, the significant acceleration can be obtained on solving Kohn–Sham equation with the proposed multilevel correction technique. The performance of the method is examined by a variety of numerical experiments.

[1]  Xavier Andrade,et al.  Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems. , 2015, Physical chemistry chemical physics : PCCP.

[2]  P. Pulay Improved SCF convergence acceleration , 1982 .

[3]  Ricardo H. Nochetto,et al.  Data Oscillation and Convergence of Adaptive FEM , 2000, SIAM J. Numer. Anal..

[4]  Aihui Zhou,et al.  Numerical analysis of finite dimensional approximations of Kohn–Sham models , 2011, Adv. Comput. Math..

[5]  Gang Bao,et al.  Numerical Solution of the Kohn-Sham Equation by Finite Element Methods with an Adaptive Mesh Redistribution Technique , 2013, J. Sci. Comput..

[6]  Rob P. Stevenson,et al.  Optimality of a Standard Adaptive Finite Element Method , 2007, Found. Comput. Math..

[7]  Johnson,et al.  Modified Broyden's method for accelerating convergence in self-consistent calculations. , 1988, Physical review. B, Condensed matter.

[8]  M. Ortiz,et al.  Non-periodic finite-element formulation of Kohn–Sham density functional theory , 2010 .

[9]  Hehu Xie,et al.  Computable Error Estimates for Ground State Solution of Bose-Einstein Condensates , 2016 .

[10]  ROB STEVENSON,et al.  The completion of locally refined simplicial partitions created by bisection , 2008, Math. Comput..

[11]  Christian Kreuzer,et al.  Quasi-Optimal Convergence Rate for an Adaptive Finite Element Method , 2008, SIAM J. Numer. Anal..

[12]  Yvon Maday,et al.  Numerical analysis of the planewave discretization of some orbital-free and Kohn-Sham models , 2010, 1003.1612.

[13]  Hehu Xie,et al.  A multi-level correction scheme for eigenvalue problems , 2011, Math. Comput..

[14]  Ricardo H. Nochetto,et al.  Convergence of Adaptive Finite Element Methods , 2002, SIAM Rev..

[15]  Hehu Xie,et al.  A full multigrid method for eigenvalue problems , 2016, J. Comput. Phys..

[16]  Michael Vogelius,et al.  Feedback and adaptive finite element solution of one-dimensional boundary value problems , 1984 .

[17]  Gang Bao,et al.  An h-adaptive finite element solver for the calculations of the electronic structures , 2012, J. Comput. Phys..

[18]  Jinchao Xu,et al.  Lower bounds of the discretization error for piecewise polynomials , 2013, Math. Comput..

[19]  Jean-François Méhaut,et al.  Wavelet-Based Density Functional Theory on Massively Parallel Hybrid Architectures , 2016 .

[20]  Hehu Xie,et al.  A multigrid method for the ground state solution of Bose–Einstein condensates based on Newton iteration , 2014, BIT Numerical Mathematics.

[21]  Mark S. C. Reed,et al.  Method of Modern Mathematical Physics , 1972 .

[22]  Aihui Zhou,et al.  Adaptive Finite Element Approximations for Kohn-Sham Models , 2013, Multiscale Model. Simul..

[23]  Juan C. Meza,et al.  KSSOLV—a MATLAB toolbox for solving the Kohn-Sham equations , 2009, TOMS.

[24]  W. Rheinboldt,et al.  Error Estimates for Adaptive Finite Element Computations , 1978 .

[25]  Yvon Maday,et al.  Numerical Analysis of Nonlinear Eigenvalue Problems , 2009, J. Sci. Comput..

[26]  I. Babuska,et al.  The finite element method and its reliability , 2001 .

[27]  George Bosilca,et al.  Open MPI: Goals, Concept, and Design of a Next Generation MPI Implementation , 2004, PVM/MPI.

[28]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[29]  Donald G. M. Anderson Iterative Procedures for Nonlinear Integral Equations , 1965, JACM.

[30]  I. Babuska,et al.  A feedback element method with a posteriori error estimation: Part I. The finite element method and some basic properties of the a posteriori error estimator , 1987 .

[31]  Xie He A multigrid method for nonlinear eigenvalue problems , 2015 .

[32]  E Weinan,et al.  Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation , 2011, J. Comput. Phys..

[33]  Hansen,et al.  Supershells in metal clusters. , 1990, Physical review. B, Condensed matter.

[34]  Vikram Gavini,et al.  Higher-order adaptive finite-element methods for Kohn-Sham density functional theory , 2012, J. Comput. Phys..

[35]  Hehu Xie,et al.  A full multigrid method for eigenvalue problems , 2014, J. Comput. Phys..

[36]  Ricardo H. Nochetto,et al.  Convergence of Adaptive Finite Element Methods for General Second Order Linear Elliptic PDEs , 2005, SIAM J. Numer. Anal..

[37]  Hehu Xie,et al.  A full multigrid method for nonlinear eigenvalue problems , 2015, 1502.04657.

[38]  J. Pask,et al.  Finite element methods in ab initio electronic structure calculations , 2005 .

[39]  Bin Yang,et al.  Eigenfunction Behavior and Adaptive Finite Element Approximations of Nonlinear Eigenvalue Problems in Quantum Physics , 2011, ESAIM: Mathematical Modelling and Numerical Analysis.

[40]  W. Dörfler A convergent adaptive algorithm for Poisson's equation , 1996 .

[41]  Aihui Zhou An analysis of finite-dimensional approximations for the ground state solution of Bose?Einstein condensates , 2004 .

[42]  P. Pulay Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .