Efficient Parallel Linear Scaling Method to get the Response Density Matrix in All-Electron Real-Space Density-Functional Perturbation Theory

Abstract

The real-space density-functional perturbation theory (DFPT) for the computations of the response properties with respect to the atomic displacement and homogeneous electric field perturbation has been recently developed and implemented into the all-electron, numeric atom-centered orbitals electronic structure package FHI-aims. It is found that the bottleneck for large scale applications is the computation of the response density matrix, which scales as $O(N^3)$. Here for the response properties with respect to the homogeneous electric field, we present an efficient parallel linear scaling algorithm for the response density matrix calculation. Our scheme is based on the second-order trace-correcting purification and the parallel sparse matrix-matrix multiplication algorithms. The new scheme reduces the formal scaling from $O(N^3)$ to $O(N)$, and shows good parallel scalability over tens of thousands of cores. As demonstrated by extensive validation, we achieve a rapid computation of accurate polarizabilities using DFPT. Finally, the computational efficiency of this scheme has been illustrated by making the scaling tests and scalability tests on massively parallel computer systems.

References in corpus (15)

• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• Phonons and related properties of extended systems from density-functional perturbation theory
• Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions
• O(N) methods in electronic structure calculations
• Order-N implementation of exact exchange in extended systems
• Calculations on millions of atoms with DFT: Linear scaling shows its potential
• Expansion algorithm for the density matrix
• ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers
• Exploiting Multiple Levels of Parallelism in Sparse Matrix-Matrix Multiplication
• Density Matrix Perturbation Theory
• Lattice Dynamics Calculations based on Density-functional Perturbation Theory in Real Space
• Ab initio linear scaling response theory: Electric polarizability by perturbed projection
• All-Electron, Real-Space Perturbation Theory for Homogeneous Electric Fields: Theory, Implementation, and Application within DFT
• Increasing the Efficiency of Sparse Matrix-Matrix Multiplication with a 2.5D Algorithm and One-Sided MPI

Cited by in corpus (1)

• Graph-based Quantum Response Theory and Shadow Born-Oppenheimer Molecular Dynamics