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:2507.22269 · Physical review. B./Physical review. B

Scan calculation of the density of states: Real-space cluster perturbation theory applied to the inhomogeneous Hubbard model in one dimension

Kaito Matsuki, Chisa Hotta, Kenichi Asano

Published Jul 14, 2025• 0 citations
arXivHTMLPDF

Abstract

We present the spectral analysis of a one-dimensional Hubbard model with a parabolic potential, using a real-space cluster perturbation theory (rCPT) designed to study spatially inhomogeneous electron systems with strong correlation. It is a natural extension of the conventional CPT to inhomogeneous cases by computing local Green's functions while averaging over multiple cluster boundaries. We find that the local density of states of at each site mirrors that of a homogeneous system with the same local filling. This insight offers a perspective on spectral evolution in inhomogeneous systems used to scan the occupancy-dependent features of the homogeneous system, making this setup a practical one-shot spectroscopy tool for ultracold atoms in harmonic traps.

DOI
10.1103/q6vk-cjg8

Related papers

A Cluster Approach Applied to the One-Dimensional Anderson-Hubbard Model

Donovan James McMurrich Allum

S. Johri and R. Bhatt developed a real-space renormalization group approach aimed at extracting the localized single-particle eigenstates of the Anderson model from a large system by identifying clusters of resonant site potentials. E. Campbell generalized this real-space renormalization group approach using standard perturbation theory. Both approaches were intended to approximate the single-particle density of states of the Anderson model. In this thesis, we aimed to test the potential of applying a similar real-space renormalization group approach to calculate the density of states of the interacting Anderson-Hubbard model. Our interest in the density of states of this model is due to a V-shaped zero-bias anomaly in two-dimensional systems. A real-space renormalization group approach is best applied to a one-dimensional system. We found that the zero-bias anomaly is not V-shaped in one-dimension. To test the potential of a real-space renormalization group approach, we used the cluster approach which is the same as the non-interacting renormalization group approach but without the perturbation theory and found that for strong disorder this technique could accurately calculate the density of states over a wide range of energies but deviated from exact results at the band edge, at $omega=pm U$ and near $omega=0$. The first two inaccuracies will be reduced with a proper real-space renormalization group approach. We suspect that the last inaccuracy is associated with long range physics and may be difficult to recover. We also developed a technique that adjusts the identification of clusters in the cluster approach to improve the computation time of the density of states with minimal loss of accuracy in a tunable range around the Fermi level. We found that this technique significantly reduced the computation time and was able to preserve the density of states near the Fermi level, except at the smallest energies near $omega=0$. Author Keywords: Anderson-Hubbard model, renormalization group, Strong electron correlations, Zero-bias anomaly

•Jan 1, 2019•0 citations

Inhomogeneous spectral moment sum rules for the retarded Green function and self-energy of strongly correlated electrons or ultracold fermionic atoms in optical lattices

James K. Freericks, Volodymyr Turkowski

Spectral moment sum rules are presented for the inhomogeneous many-body problem described by the fermionic Falicov-Kimball or Hubbard models. These local sum rules allow for arbitrary hoppings, site energies, and interactions. They can be employed to quantify the accuracy of numerical solutions to the inhomogeneous many-body problem such as strongly correlated multilayered devices, ultracold atoms in an optical lattice with a trap potential, strongly correlated systems that are disordered, or systems with nontrivial spatial ordering such as a charge-density wave or a spin-density wave. We also show how the spectral moment sum rules determine the asymptotic behavior of the Green function, self-energy, and dynamical mean field when applied to the dynamical mean-field theory solution of the many-body problem. In particular, we illustrate in detail how one can dramatically reduce the number of Matsubara frequencies needed to solve the Falicov-Kimball model while still retaining high precision, and we sketch how one can incorporate these results into Hirsch-Fye quantum Monte Carlo solvers for the Hubbard (or more complicated) models. Since the solution of inhomogeneous problems is significantly more time consuming than periodic systems, efficient use of these sum rules can provide a dramatic speed up in the computational time required to solve the many-body problem. We also discuss how these sum rules behave in nonequilibrium situations as well, where the Hamiltonian has explicit time dependence due to a driving field or due to the time-dependent change in a parameter such as the interaction strength or the origin of the trap potential.

arXiv:0907.1284•Physical Review B•Sep 22, 2009•13 citations

Applicability and limitations of cluster perturbation theory for Hubbard models

Nicklas Enenkel, Markus Garst, Peter Schmitteckert

We present important use cases and limitations when considering results obtained from cluster perturbation theory (CPT). CPT combines the solutions of small individual clusters of an infinite lattice system with the Bloch theory of conventional band theory to provide an approximation for the Green’s function in the thermodynamic limit. To this end, we are investigating single-band and multi-band Hubbard models in 1D and 2D systems. A special interest is taken in the supposed pseudogap regime of the 2D square lattice at half-filling and intermediate interaction strength $(U≤3t)$ as well as the metal–insulator transition. We point out that the finite-size level spacing of the cluster limits the resolution of spectral features within CPT. This restricts the investigation of asymptotic properties of the metal–insulator transition, as it would require much larger cluster sizes that are beyond computational capabilities.

arXiv:2309.14035•Jan 1, 2023•0 citations

Applicability and limitations of cluster perturbation theory for Hubbard models

Nicklas Enenkel, Markus Garst, Peter Schmitteckert

Abstract We present important use cases and limitations when considering results obtained from cluster perturbation theory (CPT). CPT combines the solutions of small individual clusters of an infinite lattice system with the Bloch theory of conventional band theory to provide an approximation for the Green’s function in the thermodynamic limit. To this end, we are investigating single-band and multi-band Hubbard models in 1D and 2D systems. A special interest is taken in the supposed pseudogap regime of the 2D square lattice at half-filling and intermediate interaction strength ( $$U \le 3t$$ U ≤ 3 t ) as well as the metal–insulator transition. We point out that the finite-size level spacing of the cluster limits the resolution of spectral features within CPT. This restricts the investigation of asymptotic properties of the metal–insulator transition, as it would require much larger cluster sizes that are beyond computational capabilities.

The European Physical Journal Special Topics•Sep 27, 2023•9 citations

The cumulant Green’s functions method for the Hubbard model

R.N. Lira, Peter S. Riseborough, Jereson Silva-Valencia +1 more

Abstract We use the cumulant Green’s functions method (CGFM) to study the single-band Hubbard model. The starting point of the method is to diagonalize a cluster (‘seed’) containing N correlated sites and employ the cumulants calculated from the cluster solution to obtain the full Green’s functions for the lattice. All calculations are done directly; no variational or self-consistent process is needed. We benchmark the one-dimensional results for the gap, the double occupancy, and the ground-state energy as functions of the electronic correlation at half-filling and the occupation numbers as functions of the chemical potential obtained from the CGFM against the corresponding results of the thermodynamic Bethe ansatz and the quantum transfer matrix methods. The particle-hole symmetry of the density of states is fulfilled, and the gap, occupation numbers, and ground-state energy tend systematically to the known results as the cluster size increases. We include a straightforward application of the CGFM to simulate the singles occupation of an optical lattice experiment with lithium-6 atoms in an eight-site Fermi-Hubbard chain near half-filling. The method can be applied to any parameter space for one, two, or three-dimensional Hubbard Hamiltonians and extended to other strongly correlated models, like the Anderson Hamiltonian, the t − J, Kondo, and Coqblin-Schrieffer models.

arXiv:2205.01161•Journal of Physics Condensed Matter•Mar 21, 2023•2 citations

Spectral function of electron-phonon models by cluster perturbation theory

Martin Hohenadler, Markus Aichhorn, Wolfgang von der Linden

Cluster perturbation theory in combination with the Lanczos method is used to compute the one-electron spectral function of the Holstein polaron in one and two dimensions. It is shown that the method allows reliable calculations using relatively small clusters, and at the same time significantly reduces finite-size effects. Results are compared with exact data and the relation to existing work is discussed. We also use a strong-coupling perturbation theory---equivalent to the Hubbard I approximation---to calculate the spectral function of the quarter-filled Holstein model of spinless fermions, starting from the exact atomic-limit Green function. The results agree well with previous calculations within the many-body coherent potential approximation.

arXiv:cond-mat/0306740•Physical review. B, Condensed matter•Nov 19, 2003•56 citations
We present the spectral analysis of a one-dimensional Hubbard model with a parabolic potential, using a real-space cluster perturbation theory (rCPT) designed to study spatially inhomogeneous electron systems with strong correlation. It is a natural extension of the conventional CPT to inhomogeneous cases by computing local Green's functions while averaging over multiple cluster boundaries. We find that the local density of states of at each site mirrors that of a homogeneous system with the same local filling. This insight offers a perspective on spectral evolution in inhomogeneous systems used to scan the occupancy-dependent features of the homogeneous system, making this setup a practical one-shot spectroscopy tool for ultracold atoms in harmonic traps.
S. Johri and R. Bhatt developed a real-space renormalization group approach aimed at extracting the localized single-particle eigenstates of the Anderson model from a large system by identifying clusters of resonant site potentials. E. Campbell generalized this real-space renormalization group approach using standard perturbation theory. Both approaches were intended to approximate the single-particle density of states of the Anderson model. In this thesis, we aimed to test the potential of applying a similar real-space renormalization group approach to calculate the density of states of the interacting Anderson-Hubbard model. Our interest in the density of states of this model is due to a V-shaped zero-bias anomaly in two-dimensional systems. A real-space renormalization group approach is best applied to a one-dimensional system. We found that the zero-bias anomaly is not V-shaped in one-dimension. To test the potential of a real-space renormalization group approach, we used the cluster approach which is the same as the non-interacting renormalization group approach but without the perturbation theory and found that for strong disorder this technique could accurately calculate the density of states over a wide range of energies but deviated from exact results at the band edge, at omega=pmUomega=pm Uomega=pmU and near omega=0omega=0omega=0. The first two inaccuracies will be reduced with a proper real-space renormalization group approach. We suspect that the last inaccuracy is associated with long range physics and may be difficult to recover. We also developed a technique that adjusts the identification of clusters in the cluster approach to improve the computation time of the density of states with minimal loss of accuracy in a tunable range around the Fermi level. We found that this technique significantly reduced the computation time and was able to preserve the density of states near the Fermi level, except at the smallest energies near omega=0omega=0omega=0. Author Keywords: Anderson-Hubbard model, renormalization group, Strong electron correlations, Zero-bias anomaly
Spectral moment sum rules are presented for the inhomogeneous many-body problem described by the fermionic Falicov-Kimball or Hubbard models. These local sum rules allow for arbitrary hoppings, site energies, and interactions. They can be employed to quantify the accuracy of numerical solutions to the inhomogeneous many-body problem such as strongly correlated multilayered devices, ultracold atoms in an optical lattice with a trap potential, strongly correlated systems that are disordered, or systems with nontrivial spatial ordering such as a charge-density wave or a spin-density wave. We also show how the spectral moment sum rules determine the asymptotic behavior of the Green function, self-energy, and dynamical mean field when applied to the dynamical mean-field theory solution of the many-body problem. In particular, we illustrate in detail how one can dramatically reduce the number of Matsubara frequencies needed to solve the Falicov-Kimball model while still retaining high precision, and we sketch how one can incorporate these results into Hirsch-Fye quantum Monte Carlo solvers for the Hubbard (or more complicated) models. Since the solution of inhomogeneous problems is significantly more time consuming than periodic systems, efficient use of these sum rules can provide a dramatic speed up in the computational time required to solve the many-body problem. We also discuss how these sum rules behave in nonequilibrium situations as well, where the Hamiltonian has explicit time dependence due to a driving field or due to the time-dependent change in a parameter such as the interaction strength or the origin of the trap potential.
We present important use cases and limitations when considering results obtained from cluster perturbation theory (CPT). CPT combines the solutions of small individual clusters of an infinite lattice system with the Bloch theory of conventional band theory to provide an approximation for the Green’s function in the thermodynamic limit. To this end, we are investigating single-band and multi-band Hubbard models in 1D and 2D systems. A special interest is taken in the supposed pseudogap regime of the 2D square lattice at half-filling and intermediate interaction strength (U≤3t)(U≤3t)(U≤3t) as well as the metal–insulator transition. We point out that the finite-size level spacing of the cluster limits the resolution of spectral features within CPT. This restricts the investigation of asymptotic properties of the metal–insulator transition, as it would require much larger cluster sizes that are beyond computational capabilities.
Abstract We present important use cases and limitations when considering results obtained from cluster perturbation theory (CPT). CPT combines the solutions of small individual clusters of an infinite lattice system with the Bloch theory of conventional band theory to provide an approximation for the Green’s function in the thermodynamic limit. To this end, we are investigating single-band and multi-band Hubbard models in 1D and 2D systems. A special interest is taken in the supposed pseudogap regime of the 2D square lattice at half-filling and intermediate interaction strength ( U≤3tU \le 3tU≤3t U ≤ 3 t ) as well as the metal–insulator transition. We point out that the finite-size level spacing of the cluster limits the resolution of spectral features within CPT. This restricts the investigation of asymptotic properties of the metal–insulator transition, as it would require much larger cluster sizes that are beyond computational capabilities.
Abstract We use the cumulant Green’s functions method (CGFM) to study the single-band Hubbard model. The starting point of the method is to diagonalize a cluster (‘seed’) containing N correlated sites and employ the cumulants calculated from the cluster solution to obtain the full Green’s functions for the lattice. All calculations are done directly; no variational or self-consistent process is needed. We benchmark the one-dimensional results for the gap, the double occupancy, and the ground-state energy as functions of the electronic correlation at half-filling and the occupation numbers as functions of the chemical potential obtained from the CGFM against the corresponding results of the thermodynamic Bethe ansatz and the quantum transfer matrix methods. The particle-hole symmetry of the density of states is fulfilled, and the gap, occupation numbers, and ground-state energy tend systematically to the known results as the cluster size increases. We include a straightforward application of the CGFM to simulate the singles occupation of an optical lattice experiment with lithium-6 atoms in an eight-site Fermi-Hubbard chain near half-filling. The method can be applied to any parameter space for one, two, or three-dimensional Hubbard Hamiltonians and extended to other strongly correlated models, like the Anderson Hamiltonian, the t − J, Kondo, and Coqblin-Schrieffer models.
Cluster perturbation theory in combination with the Lanczos method is used to compute the one-electron spectral function of the Holstein polaron in one and two dimensions. It is shown that the method allows reliable calculations using relatively small clusters, and at the same time significantly reduces finite-size effects. Results are compared with exact data and the relation to existing work is discussed. We also use a strong-coupling perturbation theory---equivalent to the Hubbard I approximation---to calculate the spectral function of the quarter-filled Holstein model of spinless fermions, starting from the exact atomic-limit Green function. The results agree well with previous calculations within the many-body coherent potential approximation.