Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

IEEE International Conference on High Performance Computing, Data, and Analytics

Efficient maximum entropy algorithms for electronic structure

Richard N. Silver, H. Roeder, Arthur F. Voter, J. D. Kress

Published Apr 1, 1996• 0 citations
OpenAlex

Abstract

Two Chebyshev recursion methods are presented for calculations with very large sparse Hamiltonians, the kernel polynomial method (KPM) and the maximum entropy method (MEM). If limited statistical accuracy and energy resolution are acceptable, they provide linear scaling methods for the calculation of physical properties involving large numbers of eigenstates such as densities of states, spectral functions, thermodynamics, total energies for Monte Carlo simulations and forces for molecular dynamics. KPM provides a uniform approximation to a DOS, with resolution inversely proportional to the number of Chebyshev moments, while MEM can achieve significantly higher, but non-uniform, resolution at the risk of possible artifacts. This paper emphasizes efficient algorithms.

Related papers

Calculation of densities of states and spectral functions by Chebyshev recursion and maximum entropy

Richard N. Silver, Heinrich Röder

We present an efficient algorithm for calculating spectral properties of large sparse Hamiltonian matrices such as densities of states and spectral functions. The combination of Chebyshev recursion and maximum entropy achieves high-energy resolution without significant roundoff error, machine precision, or numerical instability limitations. If controlled statistical or systematic errors are acceptable, CPU and memory requirements scale linearly in the number of states. The inference of spectral properties from moments is much better conditioned for Chebyshev moments than for power moments. We adapt concepts from the kernel polynomial method, a linear Chebyshev approximation with optimized Gibbs damping, to control the accuracy of Fourier integrals of positive nonanalytic functions. We compare the performance of kernel polynomial and maximum entropy algorithms for an electronic structure example.

arXiv:cond-mat/9703229•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Oct 1, 1997•162 citations

DENSITIES OF STATES OF MEGA-DIMENSIONAL HAMILTONIAN MATRICES

Richard N. Silver, Heinrich Röder

We propose a statistical method to estimate densities of states (DOS) and thermodynamic functions of very large Hamiltonian matrices. Orthogonal polynomials are defined on the interval between lower and upper energy bounds. The DOS is represented by a kernel polynomial constructed out of polynomial moments of the DOS and modified to damp the Gibbs phenomenon. The moments are stochastically evaluated using matrixvector multiplications on Gaussian random vectors and the polynomial recurrence relations. The resulting kernel estimate is a controlled approximation to the true DOS, because it also provides estimates of statistical and systematic errors. For a given fractional energy resolution and statistical accuracy, the required cpu time and memory scale linearly in the number of states for sparse Hamiltonians. The method is demonstrated for the two-dimensional Heisenberg anti-ferromagnet with the number of states as large as 2 26 . Results are compared to exact diagonalization where available.

International Journal of Modern Physics C•Aug 1, 1994•151 citations

The kernel polynomial method

Alexander Weiße, Gerhard Wellein, Andreas Alvermann +1 more

Efficient and stable algorithms for the calculation of spectral quantities and correlation functions are some of the key tools in computational condensed-matter physics. In this paper basic properties and recent developments of Chebyshev expansion based algorithms and the kernel polynomial method are reviewed. Characterized by a resource consumption that scales linearly with the problem dimension these methods enjoyed growing popularity over the last decade and found broad application not only in physics. Representative examples from the fields of disordered systems, strongly correlated electrons, electron-phonon interaction, and quantum spin systems are discussed in detail. In addition, an illustration on how the kernel polynomial method is successfully embedded into other numerical techniques, such as cluster perturbation theory or Monte Carlo simulation, is provided.

arXiv:cond-mat/0504627•Reviews of Modern Physics•Mar 24, 2006•1,038 citations

Linear-scaling electronic structure theory: Electronic temperature in the Kernel Polynomial Method

Eunan J. McEniry, Ralf Drautz

Linear-scaling electronic structure methods based on the calculation of moments of the underlying electronic Hamiltonian offer a computationally efficient and numerically robust scheme to drive large-scale atomistic simulations, in which the quantum-mechanical nature of the electrons is explicitly taken into account. We compare the kernel polynomial method to the Fermi operator expansion method and establish a formal connection between the two approaches. We show that the convolution of the kernel polynomial method may be understood as an effective electron temperature. The results of a number of possible kernels are formally examined, and then applied to a representative tight-binding model.

arXiv:1701.01568•Jan 6, 2017•0 citations

Calculating the many-body density of states on a digital quantum computer

Alessandro Summer, Cecilia Chiaracane, Mark T. Mitchison +1 more

Quantum statistical mechanics allows us to extract thermodynamic information from a microscopic description of a many-body system. A key step is the calculation of the density of states, from which the partition function and all finite-temperature equilibrium thermodynamic quantities can be calculated. In this work, we devise and implement a quantum algorithm to perform an estimation of the density of states on a digital quantum computer which is inspired by the kernel polynomial method. Classically, the kernel polynomial method allows us to sample spectral functions via a Chebyshev polynomial expansion. Our algorithm computes moments of the expansion on quantum hardware using a combination of random-state preparation for stochastic trace evaluation and a controlled unitary operator. We use our algorithm to estimate the density of states of a nonintegrable Hamiltonian on the Quantinuum H1-1 trapped ion chip for a controlled register of 18 qubits. This not only represents a state-of-the-art calculation of thermal properties of a many-body system on quantum hardware, but also exploits the controlled unitary evolution of a many-qubit register on an unprecedented scale. Published by the American Physical Society 2024

arXiv:2303.13476•Physical Review Research•Jan 26, 2024•11 citations

Hybrid kernel polynomial method

Muhammad Irfan, Sathish R. Kuppuswamy, Dániel Varjas +4 more

The kernel polynomial method allows to sample overall spectral properties of a quantum system, while sparse diagonalization provides accurate information about a few important states. We present a method combining these two approaches without loss of performance or accuracy. We apply this hybrid kernel polynomial method to improve the computation of thermodynamic quantities and the construction of perturbative effective models, in a regime where neither of the methods is sufficient on its own. To achieve this we develop a perturbative kernel polynomial method to compute arbitrary order series expansions of expectation values. We demonstrate the efficiency of our approach on three examples: the calculation of supercurrent and inductance in a Josephson junction, the interaction of spin qubits defined in a two dimensional electron gas, and the calculation of the effective band structure in a realistic model of a semiconductor nanowire.

arXiv:1909.09649•Sep 20, 2019•1 citation
Two Chebyshev recursion methods are presented for calculations with very large sparse Hamiltonians, the kernel polynomial method (KPM) and the maximum entropy method (MEM). If limited statistical accuracy and energy resolution are acceptable, they provide linear scaling methods for the calculation of physical properties involving large numbers of eigenstates such as densities of states, spectral functions, thermodynamics, total energies for Monte Carlo simulations and forces for molecular dynamics. KPM provides a uniform approximation to a DOS, with resolution inversely proportional to the number of Chebyshev moments, while MEM can achieve significantly higher, but non-uniform, resolution at the risk of possible artifacts. This paper emphasizes efficient algorithms.
We present an efficient algorithm for calculating spectral properties of large sparse Hamiltonian matrices such as densities of states and spectral functions. The combination of Chebyshev recursion and maximum entropy achieves high-energy resolution without significant roundoff error, machine precision, or numerical instability limitations. If controlled statistical or systematic errors are acceptable, CPU and memory requirements scale linearly in the number of states. The inference of spectral properties from moments is much better conditioned for Chebyshev moments than for power moments. We adapt concepts from the kernel polynomial method, a linear Chebyshev approximation with optimized Gibbs damping, to control the accuracy of Fourier integrals of positive nonanalytic functions. We compare the performance of kernel polynomial and maximum entropy algorithms for an electronic structure example.
We propose a statistical method to estimate densities of states (DOS) and thermodynamic functions of very large Hamiltonian matrices. Orthogonal polynomials are defined on the interval between lower and upper energy bounds. The DOS is represented by a kernel polynomial constructed out of polynomial moments of the DOS and modified to damp the Gibbs phenomenon. The moments are stochastically evaluated using matrixvector multiplications on Gaussian random vectors and the polynomial recurrence relations. The resulting kernel estimate is a controlled approximation to the true DOS, because it also provides estimates of statistical and systematic errors. For a given fractional energy resolution and statistical accuracy, the required cpu time and memory scale linearly in the number of states for sparse Hamiltonians. The method is demonstrated for the two-dimensional Heisenberg anti-ferromagnet with the number of states as large as 2 26 . Results are compared to exact diagonalization where available.
Efficient and stable algorithms for the calculation of spectral quantities and correlation functions are some of the key tools in computational condensed-matter physics. In this paper basic properties and recent developments of Chebyshev expansion based algorithms and the kernel polynomial method are reviewed. Characterized by a resource consumption that scales linearly with the problem dimension these methods enjoyed growing popularity over the last decade and found broad application not only in physics. Representative examples from the fields of disordered systems, strongly correlated electrons, electron-phonon interaction, and quantum spin systems are discussed in detail. In addition, an illustration on how the kernel polynomial method is successfully embedded into other numerical techniques, such as cluster perturbation theory or Monte Carlo simulation, is provided.
Linear-scaling electronic structure methods based on the calculation of moments of the underlying electronic Hamiltonian offer a computationally efficient and numerically robust scheme to drive large-scale atomistic simulations, in which the quantum-mechanical nature of the electrons is explicitly taken into account. We compare the kernel polynomial method to the Fermi operator expansion method and establish a formal connection between the two approaches. We show that the convolution of the kernel polynomial method may be understood as an effective electron temperature. The results of a number of possible kernels are formally examined, and then applied to a representative tight-binding model.
Quantum statistical mechanics allows us to extract thermodynamic information from a microscopic description of a many-body system. A key step is the calculation of the density of states, from which the partition function and all finite-temperature equilibrium thermodynamic quantities can be calculated. In this work, we devise and implement a quantum algorithm to perform an estimation of the density of states on a digital quantum computer which is inspired by the kernel polynomial method. Classically, the kernel polynomial method allows us to sample spectral functions via a Chebyshev polynomial expansion. Our algorithm computes moments of the expansion on quantum hardware using a combination of random-state preparation for stochastic trace evaluation and a controlled unitary operator. We use our algorithm to estimate the density of states of a nonintegrable Hamiltonian on the Quantinuum H1-1 trapped ion chip for a controlled register of 18 qubits. This not only represents a state-of-the-art calculation of thermal properties of a many-body system on quantum hardware, but also exploits the controlled unitary evolution of a many-qubit register on an unprecedented scale. Published by the American Physical Society 2024
The kernel polynomial method allows to sample overall spectral properties of a quantum system, while sparse diagonalization provides accurate information about a few important states. We present a method combining these two approaches without loss of performance or accuracy. We apply this hybrid kernel polynomial method to improve the computation of thermodynamic quantities and the construction of perturbative effective models, in a regime where neither of the methods is sufficient on its own. To achieve this we develop a perturbative kernel polynomial method to compute arbitrary order series expansions of expectation values. We demonstrate the efficiency of our approach on three examples: the calculation of supercurrent and inductance in a Josephson junction, the interaction of spin qubits defined in a two dimensional electron gas, and the calculation of the effective band structure in a realistic model of a semiconductor nanowire.