Parallel eigensolvers in plane-wave Density Functional Theory

We consider the problem of parallelizing electronic structure computations in plane-wave Density Functional Theory. Because of the limited scalability of Fourier transforms, parallelism has to be found at the eigensolver level. We show how a recently proposed algorithm based on Chebyshev polynomials can scale into the tens of thousands of processors, outperforming block conjugate gradient algorithms for large computations.

[1]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[2]  Franccois Bottin,et al.  Large scale ab initio calculations based on three levels of parallelization , 2007, 0707.3405.

[3]  Richard B. Lehoucq,et al.  Basis selection in LOBPCG , 2006, J. Comput. Phys..

[4]  Stefano de Gironcoli,et al.  QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[5]  Constantine Bekas,et al.  Computing charge densities with partially reorthogonalized Lanczos , 2005, Comput. Phys. Commun..

[6]  Robert A. van de Geijn,et al.  Elemental: A New Framework for Distributed Memory Dense Matrix Computations , 2013, TOMS.

[8]  Eric Polizzi,et al.  Density-Matrix-Based Algorithms for Solving Eingenvalue Problems , 2009 .

[9]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

[10]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[11]  Yousef Saad,et al.  A spectrum slicing method for the Kohn-Sham problem , 2012, Comput. Phys. Commun..

[12]  Blöchl,et al.  Projector augmented-wave method. , 1994, Physical review. B, Condensed matter.

[13]  Lukas Krämer,et al.  Parallel solution of partial symmetric eigenvalue problems from electronic structure calculations , 2011, Parallel Comput..

[14]  Jack Dongarra,et al.  ScaLAPACK: a scalable linear algebra library for distributed memory concurrent computers , 1992, [Proceedings 1992] The Fourth Symposium on the Frontiers of Massively Parallel Computation.

[15]  D. Sorensen.  Numerical methods for large eigenvalue problems , 2002, Acta Numerica.

[16]  M. Segall,et al.  Electron and vibrational spectroscopies using DFT, plane waves and pseudopotentials: CASTEP implementation , 2010 .

[17]  Franz Franchetti,et al.  Large-scale electronic structure calculations of high-Z metals on the BlueGene/L platform , 2006, SC.

[18]  Y. Saad,et al.  Parallel self-consistent-field calculations via Chebyshev-filtered subspace acceleration. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[19]  Richard D. James,et al.  A spectral scheme for Kohn-Sham density functional theory of clusters , 2014, J. Comput. Phys..

[20]  T. Sakurai,et al.  A projection method for generalized eigenvalue problems using numerical integration , 2003 .

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

[22]  H. Rutishauser.  Simultaneous iteration method for symmetric matrices , 1970 .

[23]  Eric Polizzi,et al.  A Density Matrix-based Algorithm for Solving Eigenvalue Problems , 2009, ArXiv.

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

[25]  Xavier Gonze,et al.  Implementation of the projector augmented-wave method in the ABINIT code: Application to the study of iron under pressure , 2008 .

[26]  Taisuke Boku,et al.  A massively-parallel electronic-structure calculations based on real-space density functional theory , 2010, J. Comput. Phys..

[27]  Edoardo Di Napoli,et al.  An Optimized and Scalable Iterative Eigensolver for Sequences of Dense Eigenvalue Problems , 2014 .

[28]  Edoardo Di Napoli,et al.  An optimized and scalable eigensolver for sequences of eigenvalue problems , 2015, Concurr. Comput. Pract. Exp..

[29]  Stefan Goedecker,et al.  ABINIT: First-principles approach to material and nanosystem properties , 2009, Comput. Phys. Commun..

[30]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[31]  Andrew V. Knyazev,et al.  Toward the Optimal Preconditioned Eigensolver: Locally Optimal Block Preconditioned Conjugate Gradient Method , 2001, SIAM J. Sci. Comput..

[32]  Yousef Saad,et al.  Self-consistent-field calculations using Chebyshev-filtered subspace iteration , 2006, J. Comput. Phys..

[33]  T. Arias,et al.  Iterative minimization techniques for ab initio total energy calculations: molecular dynamics and co , 1992 .