61,324 research outputs found
Decay properties of spectral projectors with applications to electronic structure
Motivated by applications in quantum chemistry and solid state physics, we
apply general results from approximation theory and matrix analysis to the
study of the decay properties of spectral projectors associated with large and
sparse Hermitian matrices. Our theory leads to a rigorous proof of the
exponential off-diagonal decay ("nearsightedness") for the density matrix of
gapped systems at zero electronic temperature in both orthogonal and
non-orthogonal representations, thus providing a firm theoretical basis for the
possibility of linear scaling methods in electronic structure calculations for
non-metallic systems. We further discuss the case of density matrices for
metallic systems at positive electronic temperature. A few other possible
applications are also discussed.Comment: 63 pages, 13 figure
ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers
Solving the electronic structure from a generalized or standard eigenproblem
is often the bottleneck in large scale calculations based on Kohn-Sham
density-functional theory. This problem must be addressed by essentially all
current electronic structure codes, based on similar matrix expressions, and by
high-performance computation. We here present a unified software interface,
ELSI, to access different strategies that address the Kohn-Sham eigenvalue
problem. Currently supported algorithms include the dense generalized
eigensolver library ELPA, the orbital minimization method implemented in
libOMM, and the pole expansion and selected inversion (PEXSI) approach with
lower computational complexity for semilocal density functionals. The ELSI
interface aims to simplify the implementation and optimal use of the different
strategies, by offering (a) a unified software framework designed for the
electronic structure solvers in Kohn-Sham density-functional theory; (b)
reasonable default parameters for a chosen solver; (c) automatic conversion
between input and internal working matrix formats, and in the future (d)
recommendation of the optimal solver depending on the specific problem.
Comparative benchmarks are shown for system sizes up to 11,520 atoms (172,800
basis functions) on distributed memory supercomputing architectures.Comment: 55 pages, 14 figures, 2 table
Wavelet-Based Linear-Response Time-Dependent Density-Functional Theory
Linear-response time-dependent (TD) density-functional theory (DFT) has been
implemented in the pseudopotential wavelet-based electronic structure program
BigDFT and results are compared against those obtained with the all-electron
Gaussian-type orbital program deMon2k for the calculation of electronic
absorption spectra of N2 using the TD local density approximation (LDA). The
two programs give comparable excitation energies and absorption spectra once
suitably extensive basis sets are used. Convergence of LDA density orbitals and
orbital energies to the basis-set limit is significantly faster for BigDFT than
for deMon2k. However the number of virtual orbitals used in TD-DFT calculations
is a parameter in BigDFT, while all virtual orbitals are included in TD-DFT
calculations in deMon2k. As a reality check, we report the x-ray crystal
structure and the measured and calculated absorption spectrum (excitation
energies and oscillator strengths) of the small organic molecule
N-cyclohexyl-2-(4-methoxyphenyl)imidazo[1,2-a]pyridin-3-amine
QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion Quantum Monte Carlo
We review recent advances in the capabilities of the open source ab initio
Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for
greater efficiency and reproducibility. The auxiliary field QMC (AFQMC)
implementation has been greatly expanded to include k-point symmetries,
tensor-hypercontraction, and accelerated graphical processing unit (GPU)
support. These scaling and memory reductions greatly increase the number of
orbitals that can practically be included in AFQMC calculations, increasing
accuracy. Advances in real space methods include techniques for accurate
computation of band gaps and for systematically improving the nodal surface of
ground state wavefunctions. Results of these calculations can be used to
validate application of more approximate electronic structure methods including
GW and density functional based techniques. To provide an improved foundation
for these calculations we utilize a new set of correlation-consistent effective
core potentials (pseudopotentials) that are more accurate than previous sets;
these can also be applied in quantum-chemical and other many-body applications,
not only QMC. These advances increase the efficiency, accuracy, and range of
properties that can be studied in both molecules and materials with QMC and
QMCPACK
Two-level Chebyshev filter based complementary subspace method: pushing the envelope of large-scale electronic structure calculations
We describe a novel iterative strategy for Kohn-Sham density functional
theory calculations aimed at large systems (> 1000 electrons), applicable to
metals and insulators alike. In lieu of explicit diagonalization of the
Kohn-Sham Hamiltonian on every self-consistent field (SCF) iteration, we employ
a two-level Chebyshev polynomial filter based complementary subspace strategy
to: 1) compute a set of vectors that span the occupied subspace of the
Hamiltonian; 2) reduce subspace diagonalization to just partially occupied
states; and 3) obtain those states in an efficient, scalable manner via an
inner Chebyshev-filter iteration. By reducing the necessary computation to just
partially occupied states, and obtaining these through an inner Chebyshev
iteration, our approach reduces the cost of large metallic calculations
significantly, while eliminating subspace diagonalization for insulating
systems altogether. We describe the implementation of the method within the
framework of the Discontinuous Galerkin (DG) electronic structure method and
show that this results in a computational scheme that can effectively tackle
bulk and nano systems containing tens of thousands of electrons, with chemical
accuracy, within a few minutes or less of wall clock time per SCF iteration on
large-scale computing platforms. We anticipate that our method will be
instrumental in pushing the envelope of large-scale ab initio molecular
dynamics. As a demonstration of this, we simulate a bulk silicon system
containing 8,000 atoms at finite temperature, and obtain an average SCF step
wall time of 51 seconds on 34,560 processors; thus allowing us to carry out 1.0
ps of ab initio molecular dynamics in approximately 28 hours (of wall time).Comment: Resubmitted version (version 2
- …