Filling-dependent doublon dynamics in the one-dimensional Hubbard model
  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. 
Roman Rausch and Michael Potthoff
cond-mat.str-el cond-mat.quant-gas