  We construct partially ionized hydrogen atmosphere models for magnetized
neutron stars in radiative equilibrium with surface fields B=10^12-5 \times
10^14 G and effective temperatures T_eff \sim a few \times 10^5-10^6 K. These
models are based on the latest equation of state and opacity results for
magnetized, partially ionized hydrogen plasmas that take into account various
magnetic and dense medium effects. The atmospheres directly determine the
characteristics of thermal emission from isolated neutron stars. For the models
with B=10^12-10^13 G, the spectral features due to neutral atoms lie at extreme
UV and very soft X-ray energy bands and therefore are difficult to observe.
However, the continuum flux is also different from the fully ionized case,
especially at lower energies. For the superstrong field models (B\ga 10^14 G),
we show that the vacuum polarization effect not only suppresses the proton
cyclotron line as shown previously, but also suppresses spectral features due
to bound species; therefore spectral lines or features in thermal radiation are
more difficult to observe when the neutron star magnetic field is \ga 10^14 G.
