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:2509.20504

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

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

Published Sep 24, 2025• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.2509.20504

Related papers

Quantum coherence in a maximally hot Hubbard chain

Anonymous, Ovidiu I. Pâţu, Balázs Dóra +1 more

We present a detailed study of the real-time dynamics and spectral properties of the one-dimensional fermionic Hubbard model at infinite temperature. Using tensor network simulations in Liouville space, we compute the single-particle Green's function and analyze its dynamics across a broad range of interaction strengths. To complement the time-domain approach, we develop a high-resolution Chebyshev expansion method within the density matrix formalism, enabling direct access to spectral functions in the frequency domain. In the non-interacting limit, we derive exact analytical expressions for the Green's function, providing a benchmark for our numerical methods. As interactions are introduced, we observe a transition in the spectral function from a sharp peak at the free dispersion to a broadened two-band structure associated with hole and doublon excitations. These features are well captured by a Hubbard-I mean-field approximation, even at intermediate coupling. At infinite interaction strength ($U = \infty$), we exploit a determinant representation of the Green's function to access both real-time and spectral properties. In this regime, the system retains a sharp, cosine-like momentum dispersion in frequency space, while the dynamics display nontrivial light-cone spreading with sub-ballistic scaling. Our results demonstrate that strong correlations and nontrivial quantum coherence can persist even at infinite temperature.

arXiv:2509.20498•Physical review. B./Physical review. B•Aug 27, 2026•0 citations

Doublon dynamics in the extended Fermi-Hubbard model

Felix Hofmann, Michael Potthoff

Two fermions occupying the same site of a lattice model with strongly repulsive Hubbard-type interaction $U$ form a doublon, a long-living excitation, the decay of which is suppressed because of energy conservation. By means of an exact-diagonalization approach based on the Krylov-space technique, we study the dynamics of a single doublon, of two doublons, and of a doublon in the presence of two additional fermions prepared locally in the initial state of the extended Hubbard model. The time dependence of the expectation value of the double occupancy at the different sites of a large one-dimensional lattice is analyzed by perturbative arguments. In this way the spatiotemporal evolution of the doublon can be understood. The initial decay takes place on a short time scale $1/U$, and the long-time average of the decayed fraction of the total double occupancy scales as $1/{U}^{2}$. We demonstrate how the dynamics of a doublon in the initial state is related to the spectrum of two-fermion excitations obtained from linear-response theory, we work out the difference between doublons composed of fermions vs doublons composed of bosons, and we show that despite the increase of phase space for inelastic decay processes, the stability of a doublon is enhanced by the presence of additional fermions on an intermediate time scale.

arXiv:1202.5149•Physical Review B•May 15, 2012•53 citations

Two-particle states in one-dimensional coupled Bose-Hubbard models

Yabo Li, Dominik Schneble, Tzu-Chieh Wei

We study dynamically coupled one-dimensional Bose-Hubbard models and solve for the wave functions and energies of two-particle eigenstates. Even though the wave functions do not directly follow the form of a Bethe ansatz, we describe an intuitive construction to express them as combinations of Choy-Haldane states for models with intra- and interspecies interactions. We find that the two-particle spectrum of the system with generic interactions comprises in general four different continua and three doublon dispersions. The existence of doublons depends on the coupling strength $\mathrm{\ensuremath{\Omega}}$ between two species of bosons, and their energies vary with $\mathrm{\ensuremath{\Omega}}$ and interaction strengths. We give details on one specific limit, i.e., with infinite interaction, and derive the spectrum for all types of two-particle states and their spatial and entanglement properties. We demonstrate the difference in time evolution under different coupling strengths and examine the relation between the long-time behavior of the system and the doublon dispersion. These dynamics can in principle be observed in cold atoms and might also be simulated by digital quantum computers.

arXiv:2201.05536•Physical Review A•May 27, 2022•3 citations

Subbands in the doped two-orbital Kanamori-Hubbard model

K. Hallberg, Y. Núñez-Fernández

We calculate and resolve with unprecedented detail the local density of states (DOS) and momentum-dependent spectral functions at zero temperature of one of the key models for strongly correlated electron materials, the degenerate two-orbital Kanamori-Hubbard model, by means of the dynamical mean-field theory, which uses the density matrix renormalization group as the impurity solver. When the system is hole doped and in the presence of a finite interorbital Coulomb interaction, we find the emergence of a novel holon-doublon in-gap subband which is split by the Hund's coupling. We also observe interesting features in the DOS, such as the splitting of the lower Hubbard band into a coherent narrowly dispersing peak around the Fermi energy, and another subband which evolves with the chemical potential. We characterize the main transitions giving rise to each subband by calculating the response functions of specific projected operators and by comparing with the energies in the atomic limit. The detailed results for the spectral functions found in this work pave the way to study with great precision the microscopic quantum behavior in correlated materials.

arXiv:2007.14923•Physical review. B./Physical review. B•Dec 23, 2020•6 citations

Finite-temperature dynamics of the Mott insulating Hubbard chain

Alberto Nocera, Fabian H. L. Eßler, Adrian E. Feiguin

We study the dynamical response of the half-filled one-dimensional Hubbard model for a range of interaction strengths $U$ and temperatures $T$ by a combination of numerical and analytical techniques. Using time-dependent density matrix renormalization group computations we find that the single-particle spectral function undergoes a crossover to a spin-incoherent Luttinger liquid regime at temperatures $T\ensuremath{\sim}J=4{t}^{2}/U$ for sufficiently large $U>4t$. At smaller values of $U$ and elevated temperatures the spectral function is found to exhibit two thermally broadened bands of excitations, reminiscent of what is found in the Hubbard-I approximation. The dynamical density-density response function is shown to exhibit a finite-temperature resonance at low frequencies inside the Mott gap, with a physical origin similar to the Villain mode in gapped quantum spin chains. We complement our numerical computations by developing an analytic strong-coupling approach to the low-temperature dynamics in the spin-incoherent regime.

arXiv:1710.06452•Physical review. B./Physical review. B•Jan 29, 2018•44 citations

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

Roman Rausch, Michael Potthoff

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.

arXiv:1510.01162•New Journal of Physics•Feb 5, 2016•25 citations
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 present a detailed study of the real-time dynamics and spectral properties of the one-dimensional fermionic Hubbard model at infinite temperature. Using tensor network simulations in Liouville space, we compute the single-particle Green's function and analyze its dynamics across a broad range of interaction strengths. To complement the time-domain approach, we develop a high-resolution Chebyshev expansion method within the density matrix formalism, enabling direct access to spectral functions in the frequency domain. In the non-interacting limit, we derive exact analytical expressions for the Green's function, providing a benchmark for our numerical methods. As interactions are introduced, we observe a transition in the spectral function from a sharp peak at the free dispersion to a broadened two-band structure associated with hole and doublon excitations. These features are well captured by a Hubbard-I mean-field approximation, even at intermediate coupling. At infinite interaction strength (U=∞U = \inftyU=∞), we exploit a determinant representation of the Green's function to access both real-time and spectral properties. In this regime, the system retains a sharp, cosine-like momentum dispersion in frequency space, while the dynamics display nontrivial light-cone spreading with sub-ballistic scaling. Our results demonstrate that strong correlations and nontrivial quantum coherence can persist even at infinite temperature.
Two fermions occupying the same site of a lattice model with strongly repulsive Hubbard-type interaction UUU form a doublon, a long-living excitation, the decay of which is suppressed because of energy conservation. By means of an exact-diagonalization approach based on the Krylov-space technique, we study the dynamics of a single doublon, of two doublons, and of a doublon in the presence of two additional fermions prepared locally in the initial state of the extended Hubbard model. The time dependence of the expectation value of the double occupancy at the different sites of a large one-dimensional lattice is analyzed by perturbative arguments. In this way the spatiotemporal evolution of the doublon can be understood. The initial decay takes place on a short time scale 1/U1/U1/U, and the long-time average of the decayed fraction of the total double occupancy scales as 1/U21/{U}^{2}1/U2. We demonstrate how the dynamics of a doublon in the initial state is related to the spectrum of two-fermion excitations obtained from linear-response theory, we work out the difference between doublons composed of fermions vs doublons composed of bosons, and we show that despite the increase of phase space for inelastic decay processes, the stability of a doublon is enhanced by the presence of additional fermions on an intermediate time scale.
We study dynamically coupled one-dimensional Bose-Hubbard models and solve for the wave functions and energies of two-particle eigenstates. Even though the wave functions do not directly follow the form of a Bethe ansatz, we describe an intuitive construction to express them as combinations of Choy-Haldane states for models with intra- and interspecies interactions. We find that the two-particle spectrum of the system with generic interactions comprises in general four different continua and three doublon dispersions. The existence of doublons depends on the coupling strength $\mathrm{\ensuremath{\Omega}}$ between two species of bosons, and their energies vary with $\mathrm{\ensuremath{\Omega}}$ and interaction strengths. We give details on one specific limit, i.e., with infinite interaction, and derive the spectrum for all types of two-particle states and their spatial and entanglement properties. We demonstrate the difference in time evolution under different coupling strengths and examine the relation between the long-time behavior of the system and the doublon dispersion. These dynamics can in principle be observed in cold atoms and might also be simulated by digital quantum computers.
We calculate and resolve with unprecedented detail the local density of states (DOS) and momentum-dependent spectral functions at zero temperature of one of the key models for strongly correlated electron materials, the degenerate two-orbital Kanamori-Hubbard model, by means of the dynamical mean-field theory, which uses the density matrix renormalization group as the impurity solver. When the system is hole doped and in the presence of a finite interorbital Coulomb interaction, we find the emergence of a novel holon-doublon in-gap subband which is split by the Hund's coupling. We also observe interesting features in the DOS, such as the splitting of the lower Hubbard band into a coherent narrowly dispersing peak around the Fermi energy, and another subband which evolves with the chemical potential. We characterize the main transitions giving rise to each subband by calculating the response functions of specific projected operators and by comparing with the energies in the atomic limit. The detailed results for the spectral functions found in this work pave the way to study with great precision the microscopic quantum behavior in correlated materials.
We study the dynamical response of the half-filled one-dimensional Hubbard model for a range of interaction strengths UUU and temperatures TTT by a combination of numerical and analytical techniques. Using time-dependent density matrix renormalization group computations we find that the single-particle spectral function undergoes a crossover to a spin-incoherent Luttinger liquid regime at temperatures $T\ensuremath{\sim}J=4{t}^{2}/U$ for sufficiently large U>4tU>4tU>4t. At smaller values of UUU and elevated temperatures the spectral function is found to exhibit two thermally broadened bands of excitations, reminiscent of what is found in the Hubbard-I approximation. The dynamical density-density response function is shown to exhibit a finite-temperature resonance at low frequencies inside the Mott gap, with a physical origin similar to the Villain mode in gapped quantum spin chains. We complement our numerical computations by developing an analytic strong-coupling approach to the low-temperature dynamics in the spin-incoherent regime.
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.