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

arXiv:2208.02236 · New Journal of Physics

Many body density of states of a system of non interacting spinless fermions

R. Lefèvre, Krissia Zawadzki, Grégoire Ithier

Published May 25, 2023• 3 citations
arXivHTMLPDF

Abstract

Abstract The modeling of out-of-equilibrium many-body quantum systems requires to go beyond low-energy physics and single or few bodies densities of states. Many-body localization, presence or lack of thermalization and quantum chaos are examples of phenomena in which states at different energy scales, including highly excited ones, contribute to dynamics and therefore affect the system’s properties. Quantifying these contributions requires the many-body density of states (MBDoS), a function whose calculation becomes challenging even for non-interacting identical particles due to the difficulty to enumerate accessible states while enforcing the exchange symmetry. In the present work, we introduce a new approach to evaluate the MBDoS in the general case of non-interacting systems of identical quantum particles. The starting point of our method is the principal component analysis of a filling matrix F describing how N particles can be distributed into L single-particle energy levels. We show that the many body spectrum can be expanded as a weighted sum of singular vectors of the filling matrix. The weighting coefficients only involve renormalized energies obtained from the single body spectrum. We illustrate our method in two classes of problems that are mapped into spinless fermions : (i) non-interacting electrons in a homogeneous tight-binding model in 1D and 2D, and (ii) interacting spins in a chain under a transverse field.

DOI
10.1088/1367-2630/acd8e5

Related papers

Dynamic correlation functions in quantum many-body systems at zero temperature

Eduardo R. Gagliano, C. A. Balseiro

We present a simple numerical approach to compute dynamic correlation functions in quantum many-body systems at zero temperature. The method is based on a projective technique for the memory-function formalism, and requires as a starting point an independent evaluation of the ground-state energy and wave function. We illustrate the method with a one-dimensional model of spinless fermions and calculate the density-density, current-current, and single-particle correlation functions. The obtained spectral functions show different types of elementary excitations: particle-hole, soliton-antisoliton, exciton, etc. We compare our results with analytical results, and with those obtained through different approaches available in the literature.

Physical review. B, Condensed matter•Dec 1, 1988•52 citations

Calculation of Canonical Properties and Excited States by Path Integral Numerical Methods

P. N. Vorontsov‐Velyaminov, Mikhail Voznesenskiy, Evgeny Polyakov +1 more

Abstract Path integral Monte Carlo method in the expanded ensemble is used for calculation of the ratio of partition functions for different classes of permutations treating the problem of several interacting identical particles (fermions) in an external field. Wang‐Landau algorithm is used for adjustment of balancing factors. For systems consisting of greater than two number of particles we propose an advanced variant of our approach which implies calculation of the ratio of positive and negative contributions to the partition function. Densities and energies of the sequence of excited states starting from the ground state for a system of non interacting quantum particles are calculated in turn, one by one, by means of considering systems with artificially excluded lowest energy levels and further obtaining of the ”ground state” of each next system constructed in this way. The idea of evaluation of densities of excited states for a quantum system with interparticle interaction by evaluating density difference between systems of different number of noninteracting Fermi‐copies of the system of interest is realized in terms of cyclic expansion formalism for a simple 2D system of two spinless fermions interacting via Coulomb repulsion (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

Contributions to Plasma Physics•May 5, 2010•1 citation

Many-body density matrices on a two-dimensional square lattice: Noninteracting and strongly interacting spinless fermions

Siew Ann Cheong, Christopher L. Henley

The reduced density matrix of an interacting system can be used as the basis for a truncation scheme, or in an unbiased method to discover the strongest kind of correlation in the ground state. In this paper, we investigate the structure of the many-body fermion density matrix of a small cluster in a square lattice. The cluster density matrix is evaluated numerically over a set of finite systems, subject to nonsquare periodic boundary conditions given by the lattice vectors ${\mathbf{R}}_{1}\ensuremath{\equiv}({R}_{1x},{R}_{1y})$ and ${\mathbf{R}}_{2}\ensuremath{\equiv}({R}_{2x},{R}_{2y})$. We then approximate the infinite-system cluster density-matrix spectrum by averaging the finite-system cluster density matrix (i) over degeneracies in the ground state, and orientations of the system relative to the cluster, to ensure it has the proper point-group symmetry; and (ii) over various twist boundary conditions to reduce finite size effects. We then compare the eigenvalue structure of the averaged cluster density matrix for noninteracting and strongly interacting spinless fermions, as a function of the filling fraction $\overline{n}$, and discuss whether it can be approximated as being built up from a truncated set of single-particle operators.

arXiv:cond-mat/0508750•Physical Review B•Oct 31, 2006•4 citations

Momentum distributions and numerical methods for strongly interacting one-dimensional spinor gases

F. Deuretzbacher, D. Becker, Luis Santos

One-dimensional spinor gases with strong $\ensuremath{\delta}$ interaction fermionize and form a spin chain. The spatial degrees of freedom of this atom chain can be described by a mapping to spinless noninteracting fermions and the spin degrees of freedom are described by a spin-chain model with nearest-neighbor interactions. Here, we compute momentum and occupation-number distributions of up to 16 strongly interacting spinor fermions and bosons as a function of their spin imbalance, the strength of an externally applied magnetic field gradient, the length of their spin, and for different excited states of the multiplet. We show that the ground-state momentum distributions resemble those of the corresponding noninteracting systems, apart from flat background distributions, which extend to high momenta. Moreover, we show that the spin order of the spin chain---in particular antiferromagnetic spin order---may be deduced from the momentum and occupation-number distributions of the system. Finally, we present efficient numerical methods for the calculation of the single-particle densities and one-body density matrix elements and of the local exchange coefficients of the spin chain for large systems containing more than 20 strongly interacting particles in arbitrary confining potentials.

arXiv:1602.06816•Physical Review A•Aug 4, 2016•42 citations

Ensemble-averaged mean-field many-body level density: An indicator of integrable versus chaotic single-particle dynamics

Georg Maier, Carolyn Echter, Juan Diego Urbina +2 more

According to the quantum chaos paradigm, the nature of a system's classical dynamics, whether integrable or chaotic, is universally reflected in the fluctuations of its quantum spectrum. However, since many-body spectra in the mean-field limit are composed of independent single-particle energy levels, their spectral fluctuations always display Poissonian behavior and hence cannot be used to distinguish underlying chaotic from integrable single-particle dynamics. We demonstrate that this distinction can, instead, be revealed from the mean many-body level density (at fixed energy) and its variance after averaging over ensembles representing different types of single-particle dynamics. This is in strong contrast to the energy-averaged mean level density (of a given system) that is assumed not to carry such information and is routinely removed to focus on universal signatures. To support our claim we systematically analyze the role of single-particle level correlations, that enter through Poisson and random matrix statistics (of various symmetry classes) into the ensemble-averaged density of states and its variance, contrasting bosonic and fermionic many-body systems. Our analytical study, together with extensive numerical simulations for systems with N≥5 particles consistently reveal significant differences (up to an order of magnitude for fermions and even larger for bosons) in the mean many-body level densities, depending on the nature of the underlying dynamics. Notably, in the fermionic case Poisson-type single-particle level fluctuations precisely cancel contributions from indistinguishability, such that the average many-body spectral density equals the (Thomas-Fermi) volume term. We further highlight the difference between the mean level density and its variance as functions of the total energy E and the excitation energy Q.

arXiv:2501.02918•Physical review. E•May 5, 2025•2 citations

Universal properties of single-particle excitations across the many-body localization transition

Atanu Jana, V. Ravi Chandra, Arti Garg

Understanding the nature of the transition from the delocalized to the many-body localized (MBL) phase is an important unresolved issue. To probe the nature of the MBL transition, we investigate the universal properties of single-particle excitations produced in highly excited many-body eigenstates of a disordered interacting quantum many-body system. In a class of one-dimensional spinless fermionic models with random disorder, we study the finite-size scaling of the ratio of typical to average values of the single-particle local density of states and the scattering rates across the MBL transition. Our results indicate that the MBL transition in this class of one-dimensional models of spinless fermions is continuous in nature. For various ranges of interactions in the system, the critical exponent $\ensuremath{\nu}$ with which the correlation length $\ensuremath{\xi}$ diverges at the transition point ${W}_{c}, \ensuremath{\xi}\ensuremath{\sim}|W\ensuremath{-}{W}_{c}{|}^{\ensuremath{-}\ensuremath{\nu}}$, satisfies the Chayes-Chayes-Fisher-Spencer (CCFS) bound $\ensuremath{\nu}\ensuremath{\ge}2/d$ where $d$ is the physical dimension of the system. We also discuss why the critical exponent obtained from finite-size scaling of the conventional diagnostic of many-body localization, the level spacing ratio, strongly violates the CCFS bound while the single-particle density of states and scattering rates are consistent with the CCFS criterion.

arXiv:2212.08480•Physical review. B./Physical review. B•Jun 26, 2024•5 citations
Abstract The modeling of out-of-equilibrium many-body quantum systems requires to go beyond low-energy physics and single or few bodies densities of states. Many-body localization, presence or lack of thermalization and quantum chaos are examples of phenomena in which states at different energy scales, including highly excited ones, contribute to dynamics and therefore affect the system’s properties. Quantifying these contributions requires the many-body density of states (MBDoS), a function whose calculation becomes challenging even for non-interacting identical particles due to the difficulty to enumerate accessible states while enforcing the exchange symmetry. In the present work, we introduce a new approach to evaluate the MBDoS in the general case of non-interacting systems of identical quantum particles. The starting point of our method is the principal component analysis of a filling matrix F describing how N particles can be distributed into L single-particle energy levels. We show that the many body spectrum can be expanded as a weighted sum of singular vectors of the filling matrix. The weighting coefficients only involve renormalized energies obtained from the single body spectrum. We illustrate our method in two classes of problems that are mapped into spinless fermions : (i) non-interacting electrons in a homogeneous tight-binding model in 1D and 2D, and (ii) interacting spins in a chain under a transverse field.
We present a simple numerical approach to compute dynamic correlation functions in quantum many-body systems at zero temperature. The method is based on a projective technique for the memory-function formalism, and requires as a starting point an independent evaluation of the ground-state energy and wave function. We illustrate the method with a one-dimensional model of spinless fermions and calculate the density-density, current-current, and single-particle correlation functions. The obtained spectral functions show different types of elementary excitations: particle-hole, soliton-antisoliton, exciton, etc. We compare our results with analytical results, and with those obtained through different approaches available in the literature.
Abstract Path integral Monte Carlo method in the expanded ensemble is used for calculation of the ratio of partition functions for different classes of permutations treating the problem of several interacting identical particles (fermions) in an external field. Wang‐Landau algorithm is used for adjustment of balancing factors. For systems consisting of greater than two number of particles we propose an advanced variant of our approach which implies calculation of the ratio of positive and negative contributions to the partition function. Densities and energies of the sequence of excited states starting from the ground state for a system of non interacting quantum particles are calculated in turn, one by one, by means of considering systems with artificially excluded lowest energy levels and further obtaining of the ”ground state” of each next system constructed in this way. The idea of evaluation of densities of excited states for a quantum system with interparticle interaction by evaluating density difference between systems of different number of noninteracting Fermi‐copies of the system of interest is realized in terms of cyclic expansion formalism for a simple 2D system of two spinless fermions interacting via Coulomb repulsion (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
The reduced density matrix of an interacting system can be used as the basis for a truncation scheme, or in an unbiased method to discover the strongest kind of correlation in the ground state. In this paper, we investigate the structure of the many-body fermion density matrix of a small cluster in a square lattice. The cluster density matrix is evaluated numerically over a set of finite systems, subject to nonsquare periodic boundary conditions given by the lattice vectors ${\mathbf{R}}_{1}\ensuremath{\equiv}({R}_{1x},{R}_{1y})$ and ${\mathbf{R}}_{2}\ensuremath{\equiv}({R}_{2x},{R}_{2y})$. We then approximate the infinite-system cluster density-matrix spectrum by averaging the finite-system cluster density matrix (i) over degeneracies in the ground state, and orientations of the system relative to the cluster, to ensure it has the proper point-group symmetry; and (ii) over various twist boundary conditions to reduce finite size effects. We then compare the eigenvalue structure of the averaged cluster density matrix for noninteracting and strongly interacting spinless fermions, as a function of the filling fraction n‾\overline{n}n, and discuss whether it can be approximated as being built up from a truncated set of single-particle operators.
One-dimensional spinor gases with strong $\ensuremath{\delta}$ interaction fermionize and form a spin chain. The spatial degrees of freedom of this atom chain can be described by a mapping to spinless noninteracting fermions and the spin degrees of freedom are described by a spin-chain model with nearest-neighbor interactions. Here, we compute momentum and occupation-number distributions of up to 16 strongly interacting spinor fermions and bosons as a function of their spin imbalance, the strength of an externally applied magnetic field gradient, the length of their spin, and for different excited states of the multiplet. We show that the ground-state momentum distributions resemble those of the corresponding noninteracting systems, apart from flat background distributions, which extend to high momenta. Moreover, we show that the spin order of the spin chain---in particular antiferromagnetic spin order---may be deduced from the momentum and occupation-number distributions of the system. Finally, we present efficient numerical methods for the calculation of the single-particle densities and one-body density matrix elements and of the local exchange coefficients of the spin chain for large systems containing more than 20 strongly interacting particles in arbitrary confining potentials.
According to the quantum chaos paradigm, the nature of a system's classical dynamics, whether integrable or chaotic, is universally reflected in the fluctuations of its quantum spectrum. However, since many-body spectra in the mean-field limit are composed of independent single-particle energy levels, their spectral fluctuations always display Poissonian behavior and hence cannot be used to distinguish underlying chaotic from integrable single-particle dynamics. We demonstrate that this distinction can, instead, be revealed from the mean many-body level density (at fixed energy) and its variance after averaging over ensembles representing different types of single-particle dynamics. This is in strong contrast to the energy-averaged mean level density (of a given system) that is assumed not to carry such information and is routinely removed to focus on universal signatures. To support our claim we systematically analyze the role of single-particle level correlations, that enter through Poisson and random matrix statistics (of various symmetry classes) into the ensemble-averaged density of states and its variance, contrasting bosonic and fermionic many-body systems. Our analytical study, together with extensive numerical simulations for systems with N≥5 particles consistently reveal significant differences (up to an order of magnitude for fermions and even larger for bosons) in the mean many-body level densities, depending on the nature of the underlying dynamics. Notably, in the fermionic case Poisson-type single-particle level fluctuations precisely cancel contributions from indistinguishability, such that the average many-body spectral density equals the (Thomas-Fermi) volume term. We further highlight the difference between the mean level density and its variance as functions of the total energy E and the excitation energy Q.
Understanding the nature of the transition from the delocalized to the many-body localized (MBL) phase is an important unresolved issue. To probe the nature of the MBL transition, we investigate the universal properties of single-particle excitations produced in highly excited many-body eigenstates of a disordered interacting quantum many-body system. In a class of one-dimensional spinless fermionic models with random disorder, we study the finite-size scaling of the ratio of typical to average values of the single-particle local density of states and the scattering rates across the MBL transition. Our results indicate that the MBL transition in this class of one-dimensional models of spinless fermions is continuous in nature. For various ranges of interactions in the system, the critical exponent $\ensuremath{\nu}$ with which the correlation length $\ensuremath{\xi}$ diverges at the transition point ${W}_{c}, \ensuremath{\xi}\ensuremath{\sim}|W\ensuremath{-}{W}_{c}{|}^{\ensuremath{-}\ensuremath{\nu}}$, satisfies the Chayes-Chayes-Fisher-Spencer (CCFS) bound $\ensuremath{\nu}\ensuremath{\ge}2/d$ where ddd is the physical dimension of the system. We also discuss why the critical exponent obtained from finite-size scaling of the conventional diagnostic of many-body localization, the level spacing ratio, strongly violates the CCFS bound while the single-particle density of states and scattering rates are consistent with the CCFS criterion.