Math @ Duke

Publications [#350485] of Jianfeng Lu
Papers Published
 Yu, VWZ; Campos, C; Dawson, W; García, A; Havu, V; Hourahine, B; Huhn, WP; Jacquelin, M; Jia, W; Keçeli, M; Laasner, R; Li, Y; Lin, L; Lu, J; Moussa, J; Roman, JE; VázquezMayagoitia, Á; Yang, C; Blum, V, ELSI — An open infrastructure for electronic structure solvers,
Computer Physics Communications, vol. 256
(November, 2020),
pp. 107459107459, Elsevier BV [doi]
(last updated on 2021/01/25)
Abstract: © 2020 Elsevier B.V. Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of nonorthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite nonperiodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an opensource software interface to facilitate the implementation and optimal use of highperformance solver libraries covering cubic scaling eigensolvers, linear scaling densitymatrixbased algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spinpolarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHIaims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of largescale benchmarks performed with Kohn–Sham densityfunctional theory and densityfunctional tightbinding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. Finally, we describe a reverse communication interface encoding matrixfree iterative solver strategies that are amenable, e.g., for use with planewave basis sets. Program summary: Program title: ELSI Interface CPC Library link to program files: http://dx.doi.org/10.17632/473mbbznrs.1 Licensing provisions: BSD 3clause Programming language: Fortran 2003, with interface to C/C++ External routines/libraries: BLACS, BLAS, BSEPACK (optional), EigenExa (optional), ELPA, FortJSON, LAPACK, libOMM, MPI, MAGMA (optional), MUMPS (optional), NTPoly, ParMETIS (optional), PETSc (optional), PEXSI, PTSCOTCH (optional), ScaLAPACK, SLEPc (optional), SuperLU_DIST Nature of problem: Solving the electronic structure from given Hamiltonian and overlap matrices in electronic structure calculations. Solution method: ELSI provides a unified software interface to facilitate the use of various electronic structure solvers including cubic scaling dense eigensolvers, linear scaling density matrix methods, and other approaches.


dept@math.duke.edu
ph: 919.660.2800
fax: 919.660.2821
 
Mathematics Department
Duke University, Box 90320
Durham, NC 277080320

