compmag: Thermal and bulk Comptonization for cylindrical accretion onto the polar cap of a magnetized neutron star

To find and download HEASARC data in the cloud, you can use astroquery.heasarc or download our new tool, hark.

The HEASARC and NuSTAR teams are greatly saddened by the sudden passing of Katja Pottschmidt. Most recently Katja was the lead scientist for the NuSTAR Guest Observer Facility (GOF), a role she had supported for many years. During her science career she worked on many other high energy astrophysics missions and played an integral role in advancing our knowledge of the universe. She was a wonderful colleague and friend and will be keenly missed by all who knew her.

compmag: Thermal and bulk Comptonization for cylindrical accretion onto the polar cap of a magnetized neutron star

This model describes the spectral formation in the accretion column onto the polar cap of a magnetized neutron star, with both thermal and bulk Comptonization processes taken into account. The details for the method adopted for the numerical solution of the radiative transfer equation are reported in Farinelli et al. (2012, A&A, 538, A67).

This model can be used for spectral fitting of both accreting X-ray pulsars and Supergiant Fast X-ray Trasients.

par1	$kT_{bb}$ , temperature of the seed blackbody spectrum (keV).
par2	, electron temperature of the accretion column (keV).
par3	$\tau$ , vertical optical depth of the accretion column, with electron cross-section equal to $10^{-3}$ of the Thomson cross-section.
par4	$\eta$ , index of the velocity profile when the accretion velocity increases towards the neutron star (valid when par8=1).
par5	$\beta_0$ , terminal velocity of the accreting matter at the neutron star surface (valid when par8=1).
par6	, radius of the accretion column in units of the neutron star Schwarzschild radius.
par7	, albedo at the neutron star surface.
par8	Flag for setting the velocity profile of the accretion column (1 $\Rightarrow$ $\beta(z) = A(Z_s/Z)^{-\eta}$ , where $A = \beta_0 (Z_0/Z_s)^\beta$ , 2 $\Rightarrow$ $\beta(\tau) = -\alpha\tau$ )
norm	$(R^2_{km} /D^2_{10}$ , where $R_{km}$ and $D_{10}$ are the accretion column radius in km and the source distance in units of 10 kpc, respectively.