Abstract:The distribution of the eigenvalues of a Hermitian matrix (or of a Hermitian matrix pencil) reveals important features of the underlying problem, whether a Hamiltonian system in physics, or a social network in behavioral sciences. However, computing all the eigenvalues explicitly is prohibitively expensive for real-world applications. This paper presents two types of methods to efficiently estimate the spectral density of a matrix pencil $(A, B)$ when both $A$ and $B$ are Hermitian and, in addition, $B$ is positive definite. The first one is based on the Kernel Polynomial Method (KPM) and the second on Gaussian quadrature by the Lanczos procedure. By employing Chebyshev polynomial approximation techniques, we can avoid direct factorizations in both methods, making the resulting algorithms suitable for large matrices. Under some assumptions, we prove bounds that suggest that the Lanczos method converges twice as fast as the KPM method. Numerical examples further indicate that the Lanczos method can provide more accurate spectral densities when the eigenvalue distribution is highly non-uniform. As an application, we show how to use the computed spectral density to partition the spectrum into intervals that contain roughly the same number of eigenvalues. This procedure, which makes it possible to compute the spectrum by parts, is a key ingredient in the new breed of eigensolvers that exploit "spectrum slicing".

On the generalised eigenvalue method and its relation to Prony and generalised pencil of function methods

Beyond Generalized Eigenvalues in Lattice Quantum Field Theory

Deflation of eigenvalues for iterative methods in lattice QCD

The Energy-Diminishing Weak Galerkin Finite Element Method for the Computation of Ground State and Excited States in Bose-Einstein Condensates

The static energy of a quark-antiquark pair from Laplacian eigenmodes

Conserved energy-momentum tensor for real-time lattice simulations

On the equivalence of Prony and Lanczos methods for Euclidean correlation functions

Accurate and robust splitting methods for the generalized Langevin equation with a positive Prony series memory kernel

Newton-Based Methods for Finding the Positive Ground State of Gross-Pitaevskii Equations

Estimating energy levels from lattice QCD correlation functions using a transfer matrix formalism

Conserved energy–momentum tensor for real-time lattice simulations

Mixed finite elements for the Gross-Pitaevskii eigenvalue problem: a priori error analysis and guaranteed lower energy bound

QED on the lattice and numerical perturbative computation of $g-2$

Spectral-function determination of complex electroweak amplitudes with lattice QCD

Analytic Solutions of the DGLAP Evolution and Theoretical Uncertainties

Derivative of Electron Density in Non-Equilibrium Green'S Function Technique and Its Application to Boost Performance of Convergence

A simple high-accuracy method for solving bound-state equations with the Cornell potential in momentum space

Block Lanczos for lattice QCD spectroscopy and matrix elements

Lattice QCD Beyond Ground States

Ewald method for polytropic potentials in arbitrary dimensionality

Fast computation of spectral densities for generalized eigenvalue problems