Dynamical Mean Field Theory with the Density Matrix Renormalization Group

A new numerical method for the solution of the Dynamical Mean Field Theory's self-consistent equations is introduced. The method uses the Density Matrix Renormalization Group technique to solve the associated impurity problem. The new algorithm makes no a priori approximations and is only limited by the number of sites that can be considered. We obtain accurate estimates of the critical values of the metal-insulator transitions and provide evidence of substructure in the Hubbard bands of the correlated metal. With this algorithm, more complex models having a larger number of degrees of freedom can be considered and finite-size effects can be minimized.Comment: 5 pages, 4 figure

Orbital selective crossover and Mott transitions in an asymmetric Hubbard model of cold atoms in optical lattices

We study the asymmetric Hubbard model at half-filling as a generic model to describe the physics of two species of repulsively interacting fermionic cold atoms in optical lattices. We use Dynamical Mean Field Theory to obtain the paramagnetic phase diagram of the model as function of temperature, interaction strength and hopping asymmetry. A Mott transition with a region of two coexistent solutions is found for all nonzero values of the hopping asymmetry. At low temperatures the metallic phase is a heavy Fermi-liquid, qualitatively analogous to the Fermi liquid state of the symmetric Hubbard model. Above a coherence temperature, an orbital-selective crossover takes place, wherein one fermionic species effectively localizes, and the resulting bad metallic state resembles the non-Fermi liquid state of the Falicov-Kimball model. We compute observables relevant to cold atom systems such as the double occupation, the specific heat and entropy and characterize their behavior in the different phases

Quantum Monte Carlo calculation of the finite temperature Mott-Hubbard transition

We present clear numerical evidence for the coexistence of metallic and insulating dynamical mean field theory(DMFT) solutions in a half-filled single-band Hubbard model with bare semicircular density of states at finite temperatures. Quantum Monte Carlo(QMC) method is used to solve the DMFT equations. We discuss important technical aspects of the DMFT-QMC which need to be taken into account in order to obtain the reliable results near the coexistence region. Among them are the critical slowing down of the iterative solutions near phase boundaries, the convergence criteria for the DMFT iterations, the interpolation of the discretized Green's function and the reduction of QMC statistical and systematic errors. Comparison of our results with those of other numerical methods is presented in a phase diagram.Comment: 4 pages, 5 figure