Adaptive Finite Element Method for Solving the Exact Kohn-Sham Equation of Density Functional Theory.

Results of the application of an adaptive finite element (FE) based solution using the FETK library of M. Holst to Density Functional Theory (DFT) approximation to the electronic structure of atoms and molecules are reported. The severe problem associated with the rapid variation of the electronic wave functions in the near singular regions of the atomic centers is treated by implementing completely unstructured simplex meshes that resolve these features around atomic nuclei. This concentrates the computational work in the regions in which the shortest length scales are necessary and provides for low resolution in regions for which there is no electron density. The accuracy of the solutions significantly improved when adaptive mesh refinement was applied, and it was found that the essential difficulties of the Kohn-Sham eigenvalues equation were the result of the singular behavior of the atomic potentials. Even though the matrix representations of the discrete Hamiltonian operator in the adaptive finite element basis are always sparse with a linear complexity in the number of discretization points, the overall memory and computational requirements for the solver implemented were found to be quite high. The number of mesh vertices per atom as a function of the atomic number Z andmore » the required accuracy e (in atomic units) was esitmated to be v (e;Z) = 122:37 * Z2:2346 /1:1173 , and the number of floating point operations per minimization step for a system of NA atoms was found to be 0(N3A*v(e,Z0) (e.g. Z=26, e=0.0015 au, and NA=100, the memory requirement and computational cost would be ~0.2 terabytes and ~25 petaflops). It was found that the high cost of the method could be reduced somewhat by using a geometric based refinement strategy to fix the error near the singularities.« less

[1]  G. Beylkin,et al.  Multiresolution quantum chemistry in multiwavelet bases: Analytic derivatives for Hartree-Fock and density functional theory. , 2004, The Journal of chemical physics.

[2]  Gregory Beylkin,et al.  Multiresolution quantum chemistry: basic theory and initial applications. , 2004, The Journal of chemical physics.

[3]  G. Strang.  Piecewise polynomials and the finite element method , 1973 .

[4]  J. Pask,et al.  Finite-element methods in electronic-structure theory , 2001 .

[5]  P. J. Hasnip,et al.  Electronic energy minimisation with ultrasoft pseudopotentials , 2006, Comput. Phys. Commun..

[6]  M. Tsukada,et al.  Electronic-structure calculations based on the finite-element method. , 1995, Physical review. B, Condensed matter.

[7]  T. Arias.  Multiresolution analysis of electronic structure: semicardinal and wavelet bases , 1998, cond-mat/9805262.

[8]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[9]  P. Dirac.  Note on Exchange Phenomena in the Thomas Atom , 1930, Mathematical Proceedings of the Cambridge Philosophical Society.

[10]  Alan Edelman,et al.  Multiscale Computation with Interpolating Wavelets , 1998 .

[11]  Gregory Beylkin,et al.  Multiresolution quantum chemistry in multiwavelet bases: Hartree-Fock exchange. , 2004, The Journal of chemical physics.

[12]  Y. Saad,et al.  Finite-difference-pseudopotential method: Electronic structure calculations without a basis. , 1994, Physical review letters.

[13]  Richard D. Hornung,et al.  Finite element approach for density functional theory calculations on locally-refined meshes , 2006, J. Comput. Phys..

[14]  Batcho Pf.  Computational method for general multicenter electronic structure calculations. , 2000 .

[15]  John R. Sabin,et al.  On some approximations in applications of Xα theory , 1979 .

[16]  Cho,et al.  Wavelets in electronic structure calculations. , 1993, Physical review letters.

[17]  Michael J. Holst,et al.  Adaptive Numerical Treatment of Elliptic Systems on Manifolds , 2001, Adv. Comput. Math..

[18]  Eric J. Bylaska,et al.  NWChem for Materials Science , 2003 .

[19]  Sullivan,et al.  Large-scale electronic-structure calculations with multigrid acceleration. , 1995, Physical review. B, Condensed matter.

[20]  Christopher Roland,et al.  Real-space multigrid methods for large-scale electronic structure problems , 1997 .

[21]  Douglas N. Arnold,et al.  Locally Adapted Tetrahedral Meshes Using Bisection , 2000, SIAM J. Sci. Comput..

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

[23]  J. Bernholc,et al.  Recent developments and applications of the real-space multigrid method , 2008 .

[24]  Wu,et al.  Higher-order finite-difference pseudopotential method: An application to diatomic molecules. , 1994, Physical review. B, Condensed matter.

[25]  Efthimios Kaxiras,et al.  Adaptive-coordinate real-space electronic-structure calculations for atoms, molecules, and solids , 1997 .

[26]  Hamann.  Generalized norm-conserving pseudopotentials. , 1989, Physical review. B, Condensed matter.

[27]  Michael J. Holst,et al.  A New Paradigm for Parallel Adaptive Meshing Algorithms , 2000, SIAM J. Sci. Comput..