Atmospheres and spectra of strongly magnetized neutron stars. III. Partially ionized hydrogen models

We construct partially ionized hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with surface fields B = 10(12)-5 x 10(14) G and effective temperatures T-eff similar to a few x 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 greater than or similar to 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 greater than or similar to 10(14) G.