Date: June 9, 2016.
The package can be downloaded at the url http://www.isdc.unige.ch/~ferrigno/images/Documents/BW_distribution/bwmod_dist.tgz
This cookbook describes the procedure to install the BW (Becker & Wolff) model, whose 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. 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.
In order to run BW, it is required to have:
To compile the BW model, you must first choose the directory in which you want to put the source code of the model and the model description file (model_bw.dat). Let us call this directory, e.g., /home/username/BW. IMPORTANT: you must have write permissions to the directory, as the local model library file libbwmod.so and some additional object files will be produced there.
The model needs the linking of GSL libraries. This can be achieved for the full Xspec installation by following the instructions at the section “Third-Party Libraries In Local Models Build” at the page https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XSappendixLocal.html.
After these instructions, it should be possible to compile and use the model in a standard way. To compile the model:
If the compilation works successfully, you should find the shared library file libbwmod.so in the directory /home/username/BW. Remember, however, that even if the compilation of the model worked successfully, you have to reload the model library every time you run xspec. In order to avoid that, go into $HOME/.xspec directory, which is automatically created for each user at first xspec running, edit (or newly create if it does not exist) the file xspec.rc, and add the line lmod bwmod /home/username/BW.
It is also worth pointing out that when you run xspec in batch mode, reading of the
$HOME/.xspec/xspec.rc file is skipped. Thus if for example you want to load some spectra
and the best-fit model (with BW) you cannot just give the command
> xspec - data.xcm model.xcm
but rather the following command
> xspec - load_bwmod.xcm data.xcm model.xcm
where the file load_bwmod.xcm, as xspec.rc, contains the line lmod compmag
/home/username/BW.
If you did not recompile XSPEC with support for the GSL, you should edit the Makefile of the local model directory.
The BW model has a large set of parameters and not all of them can be constrained in the fit. These parameters are in Table 1.
RNS | Neutron star radius in km (to be fixed) |
MNS | Neutron star mass in M⊙ (to be fixed) |
ξ | a parameter linked to the photon escape time (order of some unities) |
δ | the ratio between bulk and thermal Comptonization importances |
B | magnetic field in units of 1012 Gauss |
Ṁ | Mass accretion rate in units of 1017gs-1 |
Te | Electron temperature in units of keV |
r0 | column radius in units of m |
D | source distance in units of kpc (to be fixed) |
BBnorm | Normalization of the Black body seed photon component (fix it to zero at first) |
CYCnorm | Normalization of the Cyclotron emission seed photon component (fix it to one) |
FFnorm | Normalization of the Breemstrahlung emission seed photon component (fix it to one) |
We suggest here a procedure to use the model. It is mandatory to:
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, we strongly suggest to set the parameter BBnorm to zero and then to fix it to one for the final runs. The parameters FFnorm and CYCnorm should be fixed to one.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. For 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 filed of the continuum model to the one derived by the cyclotron scattering absorption feature(s).
For particular combinations of the parameters ξ and δ, the special functions used in the GSL libraries do not converge and an error is issued. For examples if δ is very low and ξ very high, or vice-versa. The model returns NAN values to XSPEC in these cases and a message is prompted, but the program is not exited. We noticed that this happened mostly for large values of r0. It is therefore 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. We notice that equaling the accretion and the X-ray luminosities is not granted to yield meaningful results for all sources.
For any problem and suggestions, please contact
Carlo Ferrigno email: carlo.ferrigno@unige.ch