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:0804.0617 · Physical Review Letters

Excitons in the One-Dimensional Hubbard Model: A Real-Time Study

K. A. Al-Hassanieh, Fernando Agustín Reboredo, Adrian E. Feiguin, Iván González, Elbio R. Dagotto

Published Apr 25, 2008• 66 citations
arXivHTMLPDF

Abstract

We study the real-time dynamics of a hole and doubly occupied site pair, namely, a holon and a doublon, in a 1D Hubbard insulator with on-site and nearest-neighbor Coulomb repulsion. Our analysis shows that the pair is long-lived and the expected decay mechanism to underlying spin excitations is actually inefficient. For a nonzero intersite Coulomb repulsion, we observe that part of the wave function remains in a bound state. Our study also provides insight on the holon-doublon propagation in real space. Because of the one-dimensional nature of the problem, these particles move in opposite directions even in the absence of an applied electric field. The potential relevance of our results to solar cell applications is discussed.

DOI
10.1103/physrevlett.100.166403

Related papers

Charge-spin-orbital dynamics of one-dimensional two-orbital Hubbard model

Hiroaki Onishi

We study the real-time evolution of a charge-excited state in a one-dimensional $e_{\rm g}$-orbital degenerate Hubbard model, by a time-dependent density-matrix renormalization group method. Considering a chain along the $z$ direction, electrons hop between adjacent $3z^2-r^2$ orbitals, while $x^2-y^2$ orbitals are localized. For the charge-excited state, a holon-doublon pair is introduced into the ground state at quarter filling. At initial time, there is no electron in a holon site, while a pair of electrons occupies $3z^2-r^2$ orbital in a doublon site. As the time evolves, the holon motion is governed by the nearest-neighbor hopping, but the electron pair can transfer between $3z^2-r^2$ orbital and $x^2-y^2$ orbital through the pair hopping in addition to the nearest-neighbor hopping. Thus holon and doublon propagate at different speed due to the pair hopping that is characteristic of multi-orbital systems.

arXiv:1101.2054•Journal of Physics Conference Series•Jan 1, 2010•2 citations

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

Excitons in one-dimensional Mott insulators

Fabian H. L. Eßler, Florian Gebhard, Eric Jeckelmann

We employ dynamical density-matrix renormalization-group (DDMRG) and field-theory methods to determine the frequency-dependent optical conductivity in one-dimensional extended, half-filled Hubbard models. The field-theory approach is applicable to the regime of ``small'' Mott gaps which is the most difficult to access by DDMRG. For very large Mott gaps the DDMRG recovers analytical results obtained previously by means of strong-coupling techniques. We focus on exciton formation at energies below the onset of the absorption continuum. As a consequence of spin-charge separation, these Mott-Hubbard excitons are bound states of spinless, charged excitations (``holon-antiholon'' pairs). We also determine exciton binding energies and sizes. In contrast to simple band insulators, we observe that excitons exist in the Mott-insulating phase only for a sufficiently strong intersite Coulomb repulsion. Furthermore, our results show that the exciton binding energy and size are not related in a simple way to the strength of the Coulomb interaction.

arXiv:cond-mat/0103406•Physical review. B, Condensed matter•Sep 11, 2001•100 citations

Charge creation via quantum tunneling in one-dimensional Mott insulators: A numerical study of the extended Hubbard model

Thomas Hansen, Lars Bojer Madsen, Yuta Murakami

Charge creation via quantum tunneling, i.e. dielectric breakdown, is one of the most fundamental and significant phenomena arising from strong light(field)-matter coupling. In this work, we conduct a systematic numerical analysis of quantum tunneling in one-dimensional Mott insulators described by the extended ($U$-$V$) Hubbard model. We discuss the applicability of the analytical formula for doublon-holon (DH) pair production, previously derived for the one-dimensional Hubbard model, which highlights the relationship between the tunneling threshold, the charge gap, and the correlation length. We test the formulas ability to predict both DH pair production and energy increase rate. Using tensor-network-based approaches, we demonstrate that the formula provides accurate predictions in the absence of excitonic states facilitated by the nearest-neighbor interaction $V$. However, when excitonic states emerge, the formula more accurately describes the rate of energy increase than the DH pair creation rate and in both cases gets improved by incorporating the exciton energy as the effective gap.

arXiv:2503.22481•Physical review. B./Physical review. B•Jun 26, 2025•1 citation

Excitons in one-dimensional Mott insulators

No author provided

We employ dynamical density-matrix renormalization-group (DDMRG) and field-theory methods to determine the frequency-dependent optical conductivity in one-dimensional extended. half-filled Hubbard models. The field-theory approach is applicable to the regime of "small" Mott gaps which is the most difficult to access by DDMRG. For very large Mott gaps the DDMRG recovers analytical results obtained previously by means of strong-coupling techniques. We focus on exciton formation at energies below the onset of the absorption continuum. As a consequence of spin-charge separation, these Mott-Hubbard excitons are bound states of spinless, charged excitations ("holon-antiholon" pairs). We also determine exciton binding energies and sizes. In contrast to simple band insulators, we observe that excitons exist in the Mott-insulating phase only for a sufficiently strong intersite Coulomb repulsion. Furthermore, our results show that the exciton binding energy and size are not related in a simple way to the strength of the Coulomb interaction.

•Sep 10, 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
We study the real-time dynamics of a hole and doubly occupied site pair, namely, a holon and a doublon, in a 1D Hubbard insulator with on-site and nearest-neighbor Coulomb repulsion. Our analysis shows that the pair is long-lived and the expected decay mechanism to underlying spin excitations is actually inefficient. For a nonzero intersite Coulomb repulsion, we observe that part of the wave function remains in a bound state. Our study also provides insight on the holon-doublon propagation in real space. Because of the one-dimensional nature of the problem, these particles move in opposite directions even in the absence of an applied electric field. The potential relevance of our results to solar cell applications is discussed.
We study the real-time evolution of a charge-excited state in a one-dimensional ege_{\rm g}eg​-orbital degenerate Hubbard model, by a time-dependent density-matrix renormalization group method. Considering a chain along the zzz direction, electrons hop between adjacent 3z2−r23z^2-r^23z2−r2 orbitals, while x2−y2x^2-y^2x2−y2 orbitals are localized. For the charge-excited state, a holon-doublon pair is introduced into the ground state at quarter filling. At initial time, there is no electron in a holon site, while a pair of electrons occupies 3z2−r23z^2-r^23z2−r2 orbital in a doublon site. As the time evolves, the holon motion is governed by the nearest-neighbor hopping, but the electron pair can transfer between 3z2−r23z^2-r^23z2−r2 orbital and x2−y2x^2-y^2x2−y2 orbital through the pair hopping in addition to the nearest-neighbor hopping. Thus holon and doublon propagate at different speed due to the pair hopping that is characteristic of multi-orbital systems.
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.
We employ dynamical density-matrix renormalization-group (DDMRG) and field-theory methods to determine the frequency-dependent optical conductivity in one-dimensional extended, half-filled Hubbard models. The field-theory approach is applicable to the regime of ``small'' Mott gaps which is the most difficult to access by DDMRG. For very large Mott gaps the DDMRG recovers analytical results obtained previously by means of strong-coupling techniques. We focus on exciton formation at energies below the onset of the absorption continuum. As a consequence of spin-charge separation, these Mott-Hubbard excitons are bound states of spinless, charged excitations (``holon-antiholon'' pairs). We also determine exciton binding energies and sizes. In contrast to simple band insulators, we observe that excitons exist in the Mott-insulating phase only for a sufficiently strong intersite Coulomb repulsion. Furthermore, our results show that the exciton binding energy and size are not related in a simple way to the strength of the Coulomb interaction.
Charge creation via quantum tunneling, i.e. dielectric breakdown, is one of the most fundamental and significant phenomena arising from strong light(field)-matter coupling. In this work, we conduct a systematic numerical analysis of quantum tunneling in one-dimensional Mott insulators described by the extended (UUU-VVV) Hubbard model. We discuss the applicability of the analytical formula for doublon-holon (DH) pair production, previously derived for the one-dimensional Hubbard model, which highlights the relationship between the tunneling threshold, the charge gap, and the correlation length. We test the formulas ability to predict both DH pair production and energy increase rate. Using tensor-network-based approaches, we demonstrate that the formula provides accurate predictions in the absence of excitonic states facilitated by the nearest-neighbor interaction VVV. However, when excitonic states emerge, the formula more accurately describes the rate of energy increase than the DH pair creation rate and in both cases gets improved by incorporating the exciton energy as the effective gap.
We employ dynamical density-matrix renormalization-group (DDMRG) and field-theory methods to determine the frequency-dependent optical conductivity in one-dimensional extended. half-filled Hubbard models. The field-theory approach is applicable to the regime of "small" Mott gaps which is the most difficult to access by DDMRG. For very large Mott gaps the DDMRG recovers analytical results obtained previously by means of strong-coupling techniques. We focus on exciton formation at energies below the onset of the absorption continuum. As a consequence of spin-charge separation, these Mott-Hubbard excitons are bound states of spinless, charged excitations ("holon-antiholon" pairs). We also determine exciton binding energies and sizes. In contrast to simple band insulators, we observe that excitons exist in the Mott-insulating phase only for a sufficiently strong intersite Coulomb repulsion. Furthermore, our results show that the exciton binding energy and size are not related in a simple way to the strength of the Coulomb interaction.
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.