A quantum Monte Carlo study of the one-dimensional ionic Hubbard model

Abstract

Quantum Monte Carlo methods are used to study a quantum phase transition in a 1D Hubbard model with a staggered ionic potential (D). Using recently formulated methods, the electronic polarization and localization are determined directly from the correlated ground state wavefunction and compared to results of previous work using exact diagonalization and Hartree-Fock. We find that the model undergoes a thermodynamic transition from a band insulator (BI) to a broken-symmetry bond ordered (BO) phase as the ratio of U/D is increased. Since it is known that at D = 0 the usual Hubbard model is a Mott insulator (MI) with no long-range order, we have searched for a second transition to this state by (i) increasing U at fixed ionic potential (D) and (ii) decreasing D at fixed U. We find no transition from the BO to MI state, and we propose that the MI state in 1D is unstable to bond ordering under the addition of any finite ionic potential. In real 1D systems the symmetric MI phase is never stable and the transition is from a symmetric BI phase to a dimerized BO phase, with a metallic point at the transition.