Dynamical mean field theory of polarons and bipolarons in the half-filled Holstein model

The evolution of the properties of a finite density electronic system as a function of the electron-phonon coupling is investigated in the Holstein model using the dynamical mean-field theory (DMFT) that becomes exact in infinite dimensions. We compare the spinless fermion case, in which only isolated polarons can be formed, with the spinful model in which the polarons can bind and form bipolarons. In the latter case, the bipolaronic binding occurs within DMFT as a metal-insulator transition. In the adiabatic regime in which the phonon energy is small with respect to the electron hopping we compare numerically exact DMFT results with an analytical scheme inspired by the Born-Oppenheimer procedure. Within the latter approach, a truncation of the phononic Hilbert space leads to a mapping of the original model onto an Anderson spin-fermion model. In the anti-adiabatic regime (where the phonon energy exceeds the electronic scales) the standard treatment based on Lang-Firsov canonical transformation allows one to map the original model on to an attractive Hubbard model in the spinful case. The separate analysis of the two regimes supports the numerical evidence that the presence of well-defined polaronic lattice displacements is not necessarily associated to a metal-insulator transition, which is instead due to pairing between the carriers. The finite-dimensionality effects neglected in DMFT may lead to a finite conductivity in the bipolaronic state which is, however, not always associated with polaronic distortions. At the polaron crossover the Born-Oppenheimer approximation is shown to break down due to the entanglement of the electron-phonon state.