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

TL;DR

This paper investigates the dynamics and decay of doublons in the one-dimensional Hubbard model at various fillings using DMRG and Bethe ansatz, revealing filling-dependent decay mechanisms and metastable states relevant for ultracold atom experiments.

Contribution

It provides a comprehensive analysis of doublon decay and propagation in the 1D Hubbard model across all fillings, highlighting the role of singly occupied sites and Bethe ansatz eigenstates.

Findings

Doublon decay is suppressed at high fillings due to kinematic constraints.

Partial decay occurs on short time scales, influenced by filling and interaction strength.

Decay products are metastable and better described by Bethe ansatz eigenstates.

Abstract

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.