  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.
