|Steps in the construction of a finite-element basis. 1) The domain (unit cell volume) is partitioned into subdomains called elements. 2) Polynomial basis functions are defined within each element. 3) Polynomials on neighboring elements are pieced together to form the strictly local piecewise-polynomial basis functions of the method. Because the basis functions are polynomials, the method is completely general and systematically improvable. Because they are strictly local, the method realizes the significant advantages of real-space grid approaches for the solution of large problems.|
The finite-element method is a powerful and general approach for solving partial differential and integral equations. The solutions of the Schrödinger and Poisson equations constitute the most time-consuming steps in density-functional based electronic-structure calculations. These calculations provide a means of determining materials properties completely from quantum-mechanical first principles (ab initio), with no adjustable parameters; and so provide a robust and system-independent means for understanding and predicting a wide range of properties in diverse materials systems. Many such electronic-structure methods provide accurate results for unit cells consisting of a small number of atoms (see Methods page for examples). Few methods, however, can provide accurate results for larger systems of several hundred atoms or more, a physically important regime where complex surfaces, interfaces, and materials defects can begin to be modeled. The goal of this research is to develop and apply a finite-element based method for large-scale ab initio electronic-structure calculations, extending the range of materials systems accessible by such rigorous, quantum mechanical means.
The finite-element method combines the significant advantages of both basis-oriented and real-space grid-based approaches. The finite-element method is an expansion method which uses a strictly local, piecewise polynomial basis (See Figure 1). Because the method is basis-oriented, it is variational, allows for arbitrarily accurate integrations, and allows for increased efficiency by choosing basis functions based upon physical insight. Because the basis functions are polynomials, the method is completely general and systematically improvable. Because they are strictly local, the method achieves a number of additional advantages with respect to large-scale calculations: the method produces sparse, structured matrices which are well suited to solution by efficient iterative methods and require only O(N) storage, where N is the number of atoms; the basis functions can be concentrated where needed in real-space (where the solution varies most rapidly) in order to increase the efficiency of the representation; all computations are performed directly in real-space, eliminating the need for computationally expensive transforms such as Fourier transforms which can incur large communications costs on large-scale computational platforms; and the method is well suited to parallel implementation. The method therefore combines the advantages of basis-oriented approaches such as the planewave method with the significant advantages with respect to large-scale calculations of real-space grid approaches such as the finite-difference method.
We have developed and implemented a finite-element based approach for the solution of the equations of density functional theory. We have used this approach routinely to calculate positron distributions and lifetimes in support of LLNL's experimental research program on positron annihilation in materials . We have incorporated finite-element Schrödinger and Poisson solvers into a general self-consistent density-functional electronic-structure method applicable to materials of arbitrary composition and symmetry, and to metals and insulators alike. Recent efforts, in collaboration with N. Sukumar and W. Hu at U. C. Davis, focus on increasing the efficiency of the FE representation substantially by employing modern partition-of-unity FE (PUFE) techniques to build known atomic physics into the FE basis. Initial results show order-of-magnitude reductions in basis set size relative to established planewave based methods, and thus the potential for dramatic speedups in large-scale quantum molecular dynamics on massively parallel computers. Ongoing work, in collaboration with Z. Bai at U. C. Davis, focuses the development of efficient parallel solvers and preconditioners for the large, sparse eigenproblems produced by discretization of the density-functional equations in the FE and PUFE bases.
Most recently, in collaboration with L. Lin and C. Yang at Lawrence Berkeley National Laboratory, we have undertaken the development and application of new discontinuous Galerkin based electronic structure methods to advance understanding of the chemistry and dynamics of Li-ion batteries www.dgdft-scidac.org.
- "Finite elements in electronic structure," J.E. Pask, in Encyclopedia of Applied and Computational Mathematics, Björn Engquist (Ed.), Springer, Heidelberg, 2013. (Review)
- dftatom: A robust and general Schrödinger and Dirac solver for atomic structure calculations," O. Certik, J.E. Pask, and J. Vackar, Comput. Phys. Commun., in press
- "Efficient adaptive integration of functions with sharp gradients and cusps in n-dimensional parallelepipeds," S.E. Mousavi, J.E. Pask, and N. Sukumar, Int. J. Numer. Meth. Engng. 91, 343–357 (2012).
- "Linear scaling solution of the all-electron Coulomb problem in sol," J.E. Pask, N. Sukumar, and S.E. Mousavi, Int. J. Multiscale Comput. Engng. 10, 83–99 (2012). Preprint: arXiv:1004.1765v3
- "Origin of large moments in Mn(x)Si(1-x) at small x," M. Shaughnessy, C.Y. Fong, R. Snow, K. Liu, J.E. Pask, and L.H. Yang, Appl. Phys. Lett. 95, 022515 (2009).
- "Classical and enriched finite element formulations for Bloch-periodic boundary conditions," N. Sukumar and J.E. Pask, Int. J. Numer. Meth. Engng. 77, 1121 (2009)
- "Finite element methods in ab initio electronic structure calculations," J.E. Pask and P.A. Sterne, Modelling Simul. Mater. Sci. Eng. 13, R71 (2005). (Review)
- "Real-space formalism for the electrostatic potential and total energy of solids," J.E. Pask and P.A. Sterne, Phys. Rev. B, 71, 113101 (2005).
- "Calculation of positron observables using a finite element-based approach," P.A. Sterne, J.E. Pask, and B.M. Klein, Appl. Surf. Sci. 149, 238 (1999).
Maintained by Randolph Q. Hood
The EOS & Materials Theory (EMT) Group in CMMD performs theoretical and computational condensed-matter and materials physics research in support of major NNSA and LLNL programs.