Self-consistent models for Coulomb heated X-ray pulsar atmospheres