  We develop a novel approach to understand the phases of one-dimensional
Bose-Hubbard models. We integrate the simplicity of the mean-field theory and
the numerical power of the density matrix renormalization group method to build
an effective numerical technique with moderate computational resources to
determine superfluid order parameters and correlation functions of large
one-dimensional systems. We demonstrate the applicability of this method to
directly identify superfluid, Mott insulator, and density wave phases in
Bose-Hubbard models.
