Skip to main content

The NASA/GSFC network is currently experiencing intermittent issues that affect our services. We are working on it and apologize for the inconvenience.
The offline tool hark (also on github) can be used to find and download data from the cloud.

Xspec Home Page


nsmax: Neutron Star Magnetic Atmosphere

This model has been superseded by nsmaxg.

This model interpolates from a grid of neutron star (NS) atmosphere spectra to produce a final spectrum that depends on the parameters listed below. The atmosphere spectra are obtained using the latest equation of state and opacity results for a partially ionized, strongly magnetized hydrogen or mid-Z element plasma. The models are constructed by solving the coupled radiative transfer equations for the two photon polarization modes in a magnetized medium, and the atmosphere is in radiative and hydrostatic equilibrium. The atmosphere models mainly depend on the surface effective temperature $T_{eff}$ and magnetic field strength $B$ and inclination $\Theta_B$; there is also a dependence on the surface gravity $g = (1+z_g)GM/R^2$, where $1+z_g = \sqrt{1-2GM/R}$ is the gravitational redshift and $M$ and $R$ are the NS mass and radius, respectively.

Two sets of models are given: one set with a single surface $B$ and $T_{eff}$ and a set which is constructed with $B$ and $T_{eff}$ varying across the surface according to the magnetic dipole model (for the latter, $\theta_m$ is the angle between the direction to the observer and the magnetic axis). The effective temperatures span the range $\log T_{eff} = 5.5 - 6.8$ for hydrogen and $\log T_{eff} = 5.8
- 6.9$ for mid-Z elements. The models with single ($B$,$T_{eff}$) cover the energy range 0.05–10 keV, while the models with ($B$,$T_{eff}$)-distributions cover the range 0.09–5 keV.

par1 $\log T_{eff}$, surface (unredshifted) effective temperature
par2 $1 + z_g$, gravitational redshift
par3 switch indicating model to use
norm $A = (R_{em}/d)^2$, normalization, where $R_{em}$ is the size (in km) of the emission region and $d$ is the distance (kpc) to the object Note: A is added automatically by XSPEC.

The models available by setting par3 are:

Switch Element $B$($10^{12}$ G) $\Theta_B$ $\theta_m$ $\log
g$ (cm/s$^2$) $\log T_{eff}$ $E$ (keV)
1000 H 0.01 0 N/A 2.4 5.5 – 6.7 0.05 - 10
1060 H 0.04 0 N/A 2.4 5.5 – 6.7 0.05 - 10
1085 H 0.07 0 N/A 2.4 5.5 – 6.7 0.05 - 10
1100 H 0.1 0 N/A 2.4 5.5 – 6.7 0.05 - 10
1200 H 1.0 0 N/A 0.4 – 2.5 5.5 – 6.7 0.05 - 10
1230 H 2.0 0 N/A 2.4 5.5 – 6.8 0.05 - 10
1260 H 4.0 0 N/A 2.4 5.5 – 6.8 0.05 - 10
1280 H 7.0 0 N/A 2.4 5.5 – 6.8 0.05 - 10
1300 H 10.0 0 N/A 0.4 – 2.5 5.5 – 6.7 0.05 - 10
1330 H 20.0 0 N/A 2.4 5.6 – 6.8 0.05 - 10
1350 H 30.0 0 N/A 2.4 5.7 – 6.8 0.05 - 10
1211 H 1.26 0 N/A 1.6 5.5 – 6.8 0.05 - 10
1281 H 7.0 0 N/A 1.6 5.5 – 6.8 0.05 - 10
12006 C 1.0 0 N/A 2.4 5.8 – 6.9 0.05 - 10
13006 C 10.0 0 N/A 2.4 5.8 – 6.9 0.05 - 10
12008 O 1.0 0 N/A 2.4 5.8 – 6.9 0.05 - 10
13008 O 10.0 0 N/A 2.4 5.8 – 6.9 0.05 - 10
12010 Ne 1.0 0 N/A 2.4 5.8 – 6.9 0.05 - 10
13010 Ne 10.0 0 N/A 2.4 5.8 – 6.9 0.05 - 10
123100 H 1.0 – 1.82 0-90 0 1.6 5.5 – 6.8 0.09 - 5
123190 H 1.0 – 1.82 0-90 90 1.6 5.5 – 6.8 0.09 - 5
130100 H 5.5 – 10.0 0-90 0 1.6 5.5 – 6.8 0.09 - 5
130190 H 5.5 – 10.0 0-90 90 1.6 5.5 – 6.8 0.09 - 5

If you publish results obtained using nsmax, please reference Ho, Potekhin & Chabrier (2008) and also Mori & Ho (2007) if using the mid-Z models. See Ho (2014) for discussion of nsmaxg and nsmax.