We construct partially ionized hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with fixed surface fields between B = 10 ^ { 12 } and 2 \times 10 ^ { 13 } G and effective temperatures \log T _ { \mathrm { eff } } = 5.5 – 6.8 , as well as with surface B and T _ { \mathrm { eff } } distributions around these values . The models are based on the latest equation of state and opacity results for magnetized , partially ionized hydrogen plasmas . The atmospheres directly determine the characteristics of thermal emission from the surface of neutron stars . We also incorporate these model spectra into XSPEC , under the model name NSMAX , thus allowing them to be used by the community to fit X-ray observations of neutron stars .