Abstract:A parallel implementation of an eigensolver designed for electronic structure calculations is presented. The method is applicable to computational tasks that solve a sequence of eigenvalue problems where the solution for a particular iteration is similar but not identical to the solution from the previous iteration. Such problems occur frequently when performing electronic structure calculations in which the eigenvectors are solutions to the Kohn-Sham equations. The eigenvectors are represented in some type of basis but the problem sizes are normally too large for direct diagonalization in that basis. Instead a subspace diagonalization procedure is employed in which matrix elements of the Hamiltonian operator are generated and the eigenvalues and eigenvectors of the resulting reduced matrix are obtained using a standard eigensolver from a package such as LAPACK or SCALAPACK. While this method works well and is widely used, the standard eigensolvers scale poorly on massively parallel computer systems for the matrix sizes typical of electronic structure calculations. We present a new method that utilizes a partitioned folded spectrum scheme (PFSM) that takes into account the iterative nature of the problem and performs well on massively parallel systems. Test results for a range of problems are presented that demonstrate an equivalent level of accuracy when compared to the standard eigensolvers, while also executing up to an order of magnitude faster. Unlike O(N) methods, the technique works equally well for metals and systems with unoccupied orbitals as for insulators and semiconductors. Timing and accuracy results are presented for a range of systems, including a 512 atom diamond cell, a cluster of 13 C60 molecules, bulk copper, a 216 atom silicon cell with a vacancy, using 40 unoccupied states/atom, and a 4000 atom aluminum supercell.

A Parallel Orbital-Updating Approach for Electronic Structure Calculations

A parallel orbital-updating based plane-wave basis method for electronic structure calculations

Numerical Analysis of the Parallel Orbital-Updating Approach for Eigenvalue Problems

Parallel Mesh Refinement Of Higher Order Finite Elements For Electronic Structure Calculations

AN ORTHOGONALIZATION-FREE PARALLELIZABLE FRAMEWORK FOR ALL-ELECTRON CALCULATIONS IN DENSITY FUNCTIONAL THEORY

Electronic Structure Interpolation Via Atomic Orbitals

Efficient single-grid and multi-grid solvers for real-space orbital-free density functional theory

Linear Scaling Electronic Structure Calculations With Numerical Atomic Basis Set

Adaptive Algorithm for Electronic Structure Calculations Using Reduction of Gaussian Mixtures

Direct determination of optimal real-space orbitals for correlated electronic structure of molecules

Adaptive Finite Differencing in High Accuracy Electronic Structure Calculations

Parallel Implementation of Large-Scale Linear Scaling Density Functional Theory Calculations with Numerical Atomic Orbitals in HONPAS.

Developing Orbital-Dependent Corrections for the Non-Additive Kinetic Energy in Subsystem Density Functional Theory

Parallel implementation of electronic structure eigensolver using a partitioned folded spectrum method

Orbital Expansion Variational Quantum Eigensolver

Alternative Approach to the Optimized Effective Potential Method

Enhancing initial state overlap through orbital optimization for faster molecular electronic ground-state energy estimation

Computationally Efficient Molecular Integrals of Solid Harmonic Gaussian Orbitals Using Quantum Entanglement of Angular Momentum

Simulating Periodic Systems on a Quantum Computer Using Molecular Orbitals.

Enhancing Reduced Density Matrix Functional Theory Calculations by Coupling Orbital and Occupation Optimizations

Advancing Natural Orbital Functional Calculations Through Deep Learning-Inspired Techniques for Large-Scale Strongly Correlated Electron Systems