bwcycl: Becker-Wolff self-consistent cyclotron line model

This is the implementation by Carlo Ferrigno of the model by Becker & Wolff (2007). This implementation has been validated in Ferrigno, Becker et al. (2009)

The original model is described in P. A. Becker & M. Wolff Thermal and Bulk Comptonization in Accretion-powered X-Ray Pulsars 2007, ApJ 654, 435B.

The application details are reported in C. Ferrigno, P. A. Becker, A. Segreto, T. Mineo & A. Santangelo, Study of the accreting pulsar 4U 0115+63 using a bulk and thermal Comptonization model, 2009, A&A, 498, 825.

This package is also under version control at https://gitlab.astro.unige.ch/ferrigno/bwmodel.

It is mandatory to:

- freeze the model normalization to one;
- set the source distance in kpc;
- set and freeze the neutron star parameters (default values are a good choice);
- select which source terms should be computed.

The computation of the black-body source term is time consuming, because it involves the numerical solution of an integral. Since the contribution of this component is generally negligible, the parameter `BBnorm` should be set to zero and then fixed to one for the final runs. The parameters `FFnorm` and `CYCnorm` should be fixed to one.^{6.1}

The mass accretion rate is strongly degenerate with the accretion column radius and the parameter ; it is therefore advisable to fix to a suitable value, which can be derived by equaling the X-ray luminosity to the accretion luminosity or a fraction of it. For a source in which the magnetic field is well above the plasma temperature and the contribution by the cyclotron emission term is minor, it is suggested to link the magnetic field of the continuum model to the one derived by the cyclotron scattering absorption feature(s).

For particular combinations of the parameters, the special functions used in the GSL libraries do not provide a finite value and a “Not a number” (`NAN`) is returned to Xspec.
The following parameter constraints avoid most of `NAN` occurrences:

- ,
- ,
- ,
- for ,
- .

It is important to limit the parameter ranges in a customary way and maybe tune the mass accretion rate to a value which keeps these parameters in the range suggested by physical considerations. Equaling the accretion and the X-ray luminosities is not guaranteed to yield meaningful results for all sources.

It is possible to define a derived model as:
`mdefine newbw bwcycl(Radius,Mass,csi,csiDel/(csi+10.85)**1.63,
B,Mdot,Te,r0,D,BBnorm,CYCnorm,FFnorm)`.

in the Xspec prompt or:

csiDel/(csi+10.85)**1.63,B,Mdot,Te,r0,D,BBnorm,CYCnorm,FFnorm)')

in pyXspec, with the parameter

This complex setting could permit a safe exploration of the parameter space, while avoiding most of `NAN` in the model computation.

par1 | : Neutron star radius in km (to be fixed) |

par2 | : Neutron star mass in (to be fixed) |

par3 | : Parameter linked to the photon escape time (order of some unities) |

par4 | : Ratio between bulk and thermal Comptonization importances |

par5 | : Magnetic field in units of Gauss |

par6 | : Mass accretion rate in units of |

par7 | : Electron temperature in units of keV |

par8 | : Column radius in units of m |

par9 | : Source distance in units of kpc (to be fixed) |

par10 | BBnorm : Normalization of the blackbody seed photon component (fix it to zero at first) |

par11 | CYCnorm : Normalization of the cyclotron emission seed photon component (fix it to one) |

par12 | FFnorm : Normalization of the Bremsstrahlung emission seed photon component (fix it to one) |

Contacts

For any problem and suggestions, please contact

Carlo Ferrigno
email: carlo.ferrigno@unige.ch