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:1510.01162 · New Journal of Physics

Multiplons in the two-hole excitation spectra of the one-dimensional Hubbard model

Roman Rausch, Michael Potthoff

Published Feb 5, 2016• 25 citations
arXivHTMLPDF

Abstract

Using the density-matrix renormalization group in combination with the Chebyshev polynomial expansion technique, we study the two-hole excitation spectrum of the one-dimensional Hubbard model in the entire filling range from the completely occupied band ( n = 2) down to half-filling ( n = 1). For strong interactions, the spectra reveal multiplon physics, i.e., relevant final states are characterized by two (doublon), three (triplon), four (quadruplon) and more holes, potentially forming stable compound objects or resonances with finite lifetime. These give rise to several satellites in the spectra with largely different spectral weights as well as to different two-hole, doublon–hole, two-doublon etc continua. The complex multiplon phenomenology is analyzed by interpreting not only local and k -resolved two-hole spectra but also three- and four-hole spectra for the Hubbard model and by referring to effective low-energy models. In addition, a filter-operator technique is presented and applied which allows to extract specific information on the final states at a given excitation energy. While multiplons composed of an odd number of holes do neither form stable compounds nor well-defined resonances unless a nearest-neighbor density interaction V is added to the Hamiltonian, the doublon and the quadruplon are well-defined resonances. The k -resolved four-hole spectrum at n = 2 represents an interesting special case where a completely stable quadruplon turns into a resonance by merging with the doublon–doublon continuum at a critical wave vector. For all fillings with , the doublon lifetime is strongly k -dependent and is even infinite at the Brillouin zone edges as demonstrated by k -resolved two-hole spectra. This can be traced back to the 'hidden' charge-SU(2) symmetry of the model which is explicitly broken off half-filling and gives rise to a massive collective excitation, even for arbitrary higher-dimensional but bipartite lattices.

DOI
10.1088/1367-2630/18/2/023033

Related papers

Filling-dependent doublon dynamics in the one-dimensional Hubbard model

Roman Rausch, Michael Potthoff

The fate of a local two-hole doublon excitation in the one-dimensional Fermi-Hubbard model is systematically studied for strong Hubbard interaction $U$ in the entire filling range using the density-matrix renormalization group (DMRG) and the Bethe ansatz. For strong $U$, two holes at the same site form a compound object whose decay is impeded by the lack of phase space. Still, a partial decay is possible on an extremely short time scale where phase-space arguments do not yet apply. We argue that the initial decay and the resulting intermediate state are relevant for experiments performed with ultracold atoms loaded into an optical lattice as well as for (time-resolved) CVV Auger-electron spectroscopy. The detailed discussion comprises the mixed ballistic-diffusive real-time propagation of the doublon through the lattice, its partial decay on the short time scale as a function of filling and interaction strength, as well as the analysis of the decay products, which are metastable on the intermediate time scale that is numerically accessible and which show up in the two-hole excitation (Auger) spectrum. The ambivalent role of singly occupied sites is key to understanding the doublon physics; for high fillings, ground-state configurations with single occupancies are recognized to strongly relax the kinematic constraints and to open up decay channels. For fillings close to half-filling, however, their presence actually blocks the doublon decay. Finally, the analysis of the continua in the two-hole spectrum excludes a picture where the doublon decays into unbound electron holes for generic fillings, different from the limiting case of the completely filled band. We demonstrate that the decay products as well as the doublon propagation should rather be understood in terms of Bethe ansatz eigenstates.

arXiv:1611.05388•Physical review. B./Physical review. B•Jan 31, 2017•20 citations

One- and two-hole excitation spectra of the Hubbard model in the non-half-filled case

Andrzej M. Oleś, Guy Tréglia, D. Spanjaard +1 more

Exact results for the one- and two-hole excitation spectra of the one- dimensional Hubbard model for filling larger than one electron per site are presented. We consider short chains of up to eight atoms and use the so-called modified periodic boundary conditions, with which we are able to extrapolate the results up to infinite chains. The one-hole excitation spectra demonstrate a gradual transition from the band picture to the spectra characterized by new satellite structures, with increasing electron-electron interaction U. The weight of the satellite part of these spectra increases with increasing U, approaching the atomic-limit value. In the ferromagnetic case, a substantial reduction of the exchange splitting is found as compared with the independent-electron approximation. The two-hole excitation spectra exhibit characteristic atomic lines at large values of U, while at intermediate values of U they are much broader in the non-completely-filled chain than the respective ones of the filled chain. The distance between the satellite and the band part is found to be larger than the value of U for small and intermediate values of the interaction. The results obtained support those found with second-order perturbation theory and the t-matrix approximation.

Physical review. B, Condensed matter•Oct 15, 1986•17 citations

Two-site correlation in analysis of the Hubbard model

Hideki Matsumoto, Ferdinando Mancini

Electronic states of a two-dimensional Hubbard model near half filling is investigated by use of the composite-operator method, which is improved by including two-site correlations in combined use with the resolvent method. The densities of states obtained for the upper and lower Hubbard bands have three-peek structures at half filling and give shapes above a certain temperature similar to the ones obtained in the numerical simulation by Dagotto, Ortolani, and Scalapino [Phys. Rev. B 46, 3183 (1992)], but show narrowing of the bands at low temperature due to the formation of spin-singlet states among the nearest-neighbor electrons. A doped case shows a coherent peak at the Fermi level and the feature of the density of states is very similar to the one from numerical simulation. The one-particle spectra in the (k,k) direction is also presented. The results show that, at half filling, the spectral weight is not evenly distributed between two bands; most of the weight of the lower Hubbard band appears between (0,0) and (\ensuremath{\pi}/2a,\ensuremath{\pi}/2a), while that of the upper Hubbard band lies in the region between (\ensuremath{\pi}/2a,\ensuremath{\pi}/2a) and (\ensuremath{\pi}/a,\ensuremath{\pi}/a), keeping roughly the original band distribution. With hole doping, a flat band develops at the Fermi level with a large Fermi surface and the lower Hubbard band extends through the whole Brillouin zone.

Physical review. B, Condensed matter•Jan 15, 1997•27 citations

Hot but Coherent: Doublons at Infinite Temperature in the Hubbard chain

Cătălin Paşcu Moca, Balázs Dóra, Gergely Zaránd

We investigate the spectral properties and the dynamics of doublons in the one-dimensional Hubbard model at infinite temperature. Using a Chebyshev expansion approach formulated in the superfermionic representation, we compute the momentum- and frequency-resolved doublon spectral function across a wide range of interaction strengths $U$ and in the presence of an external electric field. Increasing the interaction, the spectrum gradually splits into separate bands associated with two-hole and two-particle excitations. Despite the presence of strong correlations, we find that doublons retain their coherence and undergo long-lived Bloch oscillations, as well as rich quantum walk dynamics characterized by light-cone spreadings at strong coupling, which we analyze through the time-and space-resolved Green's function.

arXiv:2509.20504•Sep 24, 2025•0 citations

In-gap band in the one-dimensional two-orbital Kanamori-Hubbard model with interorbital Coulomb interaction

N. Aucar Boidi, H. Fernández García, Y. Núñez-Fernández +1 more

We study the electronic spectral properties at zero temperature of the one-dimensional (1D) version of the degenerate two-orbital Kanamori-Hubbard model, one of the well-established frameworks to study transition metal compounds, using state-of-the-art numerical techniques based on the density matrix renormalization group. While the system is Mott insulating for the half-filled case, as expected for an interacting 1D system, we find interesting and rich structures in the single-particle density of states (DOS) for the hole-doped system. In particular, we find the existence of in-gap states which are pulled down to lower energies from the upper Hubbard band with increasing the interorbital Coulomb interaction $V$. We analyze the composition of the DOS by projecting it onto different local excitations, and we observe that for large dopings these in-gap excitations are formed mainly by interorbital holon-doublon (HD) states and their energies follow approximately the HD states in the atomic limit. We observe that the Hund interaction $J$ increases the width of the in-gap band, as expected from the two-particle fluctuations in the Hamiltonian. The observation of a finite density of states within the gap between the Hubbard bands for this extended 1D model indicates that these systems present a rich excitation spectra which could help us understand the microscopic physics behind multiorbital compounds.

arXiv:2112.03794•Physical Review Research•Dec 23, 2021•8 citations

One-hole excitation spectra of the one-dimensional Hubbard model

Andrzej M. Oleś, Guy Tréglia, D. Spanjaard +1 more

Exact results for the one-hole excitation spectra of the one-dimensional half-filled Hubbard Hamiltonian are presented. We consider short chains of N=2, 4, 6 atoms and introduce modified periodic boundary conditions which allow us to make a comparison between different N and to extrapolate the results up to infinite chains. By the increasing weight of the scattering states in the one-hole spectrum, we demonstrate a gradual disappearance of the free-electron (band) picture with increasing local electron-electron interaction U. The exact spectra of the N=6 chain verify the results of second-order perturbation theory for U/W1, where W is the bandwidth. The extrapolation to the infinite system, made with the finite-cell-scaling method, gives the gap which agrees rather well with the exact result for U/W>0.2. The average energies and weights of the band and scattering states obtained with this method for an infinite chain do not differ significantly from those of the N=6 chain.

Physical review. B, Condensed matter•Aug 15, 1985•15 citations
Using the density-matrix renormalization group in combination with the Chebyshev polynomial expansion technique, we study the two-hole excitation spectrum of the one-dimensional Hubbard model in the entire filling range from the completely occupied band ( n = 2) down to half-filling ( n = 1). For strong interactions, the spectra reveal multiplon physics, i.e., relevant final states are characterized by two (doublon), three (triplon), four (quadruplon) and more holes, potentially forming stable compound objects or resonances with finite lifetime. These give rise to several satellites in the spectra with largely different spectral weights as well as to different two-hole, doublon–hole, two-doublon etc continua. The complex multiplon phenomenology is analyzed by interpreting not only local and k -resolved two-hole spectra but also three- and four-hole spectra for the Hubbard model and by referring to effective low-energy models. In addition, a filter-operator technique is presented and applied which allows to extract specific information on the final states at a given excitation energy. While multiplons composed of an odd number of holes do neither form stable compounds nor well-defined resonances unless a nearest-neighbor density interaction V is added to the Hamiltonian, the doublon and the quadruplon are well-defined resonances. The k -resolved four-hole spectrum at n = 2 represents an interesting special case where a completely stable quadruplon turns into a resonance by merging with the doublon–doublon continuum at a critical wave vector. For all fillings with , the doublon lifetime is strongly k -dependent and is even infinite at the Brillouin zone edges as demonstrated by k -resolved two-hole spectra. This can be traced back to the 'hidden' charge-SU(2) symmetry of the model which is explicitly broken off half-filling and gives rise to a massive collective excitation, even for arbitrary higher-dimensional but bipartite lattices.
The fate of a local two-hole doublon excitation in the one-dimensional Fermi-Hubbard model is systematically studied for strong Hubbard interaction UUU in the entire filling range using the density-matrix renormalization group (DMRG) and the Bethe ansatz. For strong UUU, two holes at the same site form a compound object whose decay is impeded by the lack of phase space. Still, a partial decay is possible on an extremely short time scale where phase-space arguments do not yet apply. We argue that the initial decay and the resulting intermediate state are relevant for experiments performed with ultracold atoms loaded into an optical lattice as well as for (time-resolved) CVV Auger-electron spectroscopy. The detailed discussion comprises the mixed ballistic-diffusive real-time propagation of the doublon through the lattice, its partial decay on the short time scale as a function of filling and interaction strength, as well as the analysis of the decay products, which are metastable on the intermediate time scale that is numerically accessible and which show up in the two-hole excitation (Auger) spectrum. The ambivalent role of singly occupied sites is key to understanding the doublon physics; for high fillings, ground-state configurations with single occupancies are recognized to strongly relax the kinematic constraints and to open up decay channels. For fillings close to half-filling, however, their presence actually blocks the doublon decay. Finally, the analysis of the continua in the two-hole spectrum excludes a picture where the doublon decays into unbound electron holes for generic fillings, different from the limiting case of the completely filled band. We demonstrate that the decay products as well as the doublon propagation should rather be understood in terms of Bethe ansatz eigenstates.
Exact results for the one- and two-hole excitation spectra of the one- dimensional Hubbard model for filling larger than one electron per site are presented. We consider short chains of up to eight atoms and use the so-called modified periodic boundary conditions, with which we are able to extrapolate the results up to infinite chains. The one-hole excitation spectra demonstrate a gradual transition from the band picture to the spectra characterized by new satellite structures, with increasing electron-electron interaction U. The weight of the satellite part of these spectra increases with increasing U, approaching the atomic-limit value. In the ferromagnetic case, a substantial reduction of the exchange splitting is found as compared with the independent-electron approximation. The two-hole excitation spectra exhibit characteristic atomic lines at large values of U, while at intermediate values of U they are much broader in the non-completely-filled chain than the respective ones of the filled chain. The distance between the satellite and the band part is found to be larger than the value of U for small and intermediate values of the interaction. The results obtained support those found with second-order perturbation theory and the t-matrix approximation.
Electronic states of a two-dimensional Hubbard model near half filling is investigated by use of the composite-operator method, which is improved by including two-site correlations in combined use with the resolvent method. The densities of states obtained for the upper and lower Hubbard bands have three-peek structures at half filling and give shapes above a certain temperature similar to the ones obtained in the numerical simulation by Dagotto, Ortolani, and Scalapino [Phys. Rev. B 46, 3183 (1992)], but show narrowing of the bands at low temperature due to the formation of spin-singlet states among the nearest-neighbor electrons. A doped case shows a coherent peak at the Fermi level and the feature of the density of states is very similar to the one from numerical simulation. The one-particle spectra in the (k,k) direction is also presented. The results show that, at half filling, the spectral weight is not evenly distributed between two bands; most of the weight of the lower Hubbard band appears between (0,0) and (\ensuremath{\pi}/2a,\ensuremath{\pi}/2a), while that of the upper Hubbard band lies in the region between (\ensuremath{\pi}/2a,\ensuremath{\pi}/2a) and (\ensuremath{\pi}/a,\ensuremath{\pi}/a), keeping roughly the original band distribution. With hole doping, a flat band develops at the Fermi level with a large Fermi surface and the lower Hubbard band extends through the whole Brillouin zone.
We investigate the spectral properties and the dynamics of doublons in the one-dimensional Hubbard model at infinite temperature. Using a Chebyshev expansion approach formulated in the superfermionic representation, we compute the momentum- and frequency-resolved doublon spectral function across a wide range of interaction strengths UUU and in the presence of an external electric field. Increasing the interaction, the spectrum gradually splits into separate bands associated with two-hole and two-particle excitations. Despite the presence of strong correlations, we find that doublons retain their coherence and undergo long-lived Bloch oscillations, as well as rich quantum walk dynamics characterized by light-cone spreadings at strong coupling, which we analyze through the time-and space-resolved Green's function.
We study the electronic spectral properties at zero temperature of the one-dimensional (1D) version of the degenerate two-orbital Kanamori-Hubbard model, one of the well-established frameworks to study transition metal compounds, using state-of-the-art numerical techniques based on the density matrix renormalization group. While the system is Mott insulating for the half-filled case, as expected for an interacting 1D system, we find interesting and rich structures in the single-particle density of states (DOS) for the hole-doped system. In particular, we find the existence of in-gap states which are pulled down to lower energies from the upper Hubbard band with increasing the interorbital Coulomb interaction VVV. We analyze the composition of the DOS by projecting it onto different local excitations, and we observe that for large dopings these in-gap excitations are formed mainly by interorbital holon-doublon (HD) states and their energies follow approximately the HD states in the atomic limit. We observe that the Hund interaction JJJ increases the width of the in-gap band, as expected from the two-particle fluctuations in the Hamiltonian. The observation of a finite density of states within the gap between the Hubbard bands for this extended 1D model indicates that these systems present a rich excitation spectra which could help us understand the microscopic physics behind multiorbital compounds.
Exact results for the one-hole excitation spectra of the one-dimensional half-filled Hubbard Hamiltonian are presented. We consider short chains of N=2, 4, 6 atoms and introduce modified periodic boundary conditions which allow us to make a comparison between different N and to extrapolate the results up to infinite chains. By the increasing weight of the scattering states in the one-hole spectrum, we demonstrate a gradual disappearance of the free-electron (band) picture with increasing local electron-electron interaction U. The exact spectra of the N=6 chain verify the results of second-order perturbation theory for U/W1, where W is the bandwidth. The extrapolation to the infinite system, made with the finite-cell-scaling method, gives the gap which agrees rather well with the exact result for U/W>0.2. The average energies and weights of the band and scattering states obtained with this method for an infinite chain do not differ significantly from those of the N=6 chain.