49 research outputs found
Accelerating Atomic Orbital-based Electronic Structure Calculation via Pole Expansion and Selected Inversion
We describe how to apply the recently developed pole expansion and selected
inversion (PEXSI) technique to Kohn-Sham density function theory (DFT)
electronic structure calculations that are based on atomic orbital
discretization. We give analytic expressions for evaluating the charge density,
the total energy, the Helmholtz free energy and the atomic forces (including
both the Hellman-Feynman force and the Pulay force) without using the
eigenvalues and eigenvectors of the Kohn-Sham Hamiltonian. We also show how to
update the chemical potential without using Kohn-Sham eigenvalues. The
advantage of using PEXSI is that it has a much lower computational complexity
than that associated with the matrix diagonalization procedure. We demonstrate
the performance gain by comparing the timing of PEXSI with that of
diagonalization on insulating and metallic nanotubes. For these quasi-1D
systems, the complexity of PEXSI is linear with respect to the number of atoms.
This linear scaling can be observed in our computational experiments when the
number of atoms in a nanotube is larger than a few hundreds. Both the wall
clock time and the memory requirement of PEXSI is modest. This makes it even
possible to perform Kohn-Sham DFT calculations for 10,000-atom nanotubes with a
sequential implementation of the selected inversion algorithm. We also perform
an accurate geometry optimization calculation on a truncated (8,0)
boron-nitride nanotube system containing 1024 atoms. Numerical results indicate
that the use of PEXSI does not lead to loss of accuracy required in a practical
DFT calculation
An Approximate Framework for Quantum Transport Calculation with Model Order Reduction
A new approximate computational framework is proposed for computing the
non-equilibrium charge density in the context of the non-equilibrium Green's
function (NEGF) method for quantum mechanical transport problems. The framework
consists of a new formulation, called the X-formulation, for single-energy
density calculation based on the solution of sparse linear systems, and a
projection-based nonlinear model order reduction (MOR) approach to address the
large number of energy points required for large applied biases. The advantages
of the new methods are confirmed by numerical experiments
Efficient Computation of Sparse Matrix Functions for Large-Scale Electronic Structure Calculations: The CheSS Library
We present CheSS, the “Chebyshev Sparse Solvers” library, which has been designed to solve typical problems arising in large-scale electronic structure calculations using localized basis sets. The library is based on a flexible and efficient expansion in terms of Chebyshev polynomials and presently features the calculation of the density matrix, the calculation of matrix powers for arbitrary powers, and the extraction of eigenvalues in a selected interval. CheSS is able to exploit the sparsity of the matrices and scales linearly with respect to the number of nonzero entries, making it well-suited for large-scale calculations. The approach is particularly adapted for setups leading to small spectral widths of the involved matrices and outperforms alternative methods in this regime. By coupling CheSS to the DFT code BigDFT, we show that such a favorable setup is indeed possible in practice. In addition, the approach based on Chebyshev polynomials can be massively parallelized, and CheSS exhibits excellent scaling up to thousands of cores even for relatively small matrix sizes.We gratefully acknowledge the support of the MaX (SM) and POP (MW) projects, which have received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement
No. 676598 and 676553, respectively. This work was also supported by the Energy oriented Centre of Excellence (EoCoE), grant agreement number 676629, funded within the Horizon2020 framework of the European
Union, as well as by the Next-Generation Supercomputer project (the K computer project) and the FLAGSHIP2020 within the priority study5 (Development of new fundamental technologies for high-efficiency
energy creation, conversion/storage and use) from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan. We (LG, DC, WD, TN) gratefully acknowledge the joint CEA-RIKEN collaboration action.Peer ReviewedPostprint (author's final draft