A&A 491, 525-529 (2008)
V. Suleimanov1,2 - J. Poutanen3 - M. Falanga4,5 - K. Werner1
1 - Institute for Astronomy and Astrophysics, Kepler Center for Astro and Particle Physics, Eberhard Karls University, Sand 1, 72076 Tübingen, Germany
2 - Kazan State University, Kremlevskaja str. 18, Kazan 420008, Russia
3 - Astronomy Division, Department of Physical Sciences, PO Box 3000, 90014 University of Oulu, Finland
4 - CEA Saclay, DSM/DAPNIA/Service d'Astrophysique (CNRS FRE 2591), 91191 Gif-sur-Yvette, France
5 - AIM, UMR 7158, CEA - CNRS - Université Paris 7, France
Received 4 May 2008 / Accepted 9 Septembre 2008
Context. The majority of cataclysmic variables observed in the hard X-ray energy band are intermediate polars where the magnetic field is strong enough to channel the accreting matter to the magnetic poles of the white dwarf. A shock above the stellar surface heats the gas to fairly high temperatures (10-100 keV). The post-shock region cools mostly via optically thin bremsstrahlung.
Aims. We investigate the influence of Compton scattering on the structure and the emergent spectrum of the post-shock region. We also study the effect it has on the mass of the white dwarfs obtained from fitting the observed X-ray spectrum of intermediate polars.
Methods. We construct the model of the post-shock region taking Compton scattering into account. The radiation transfer equation is solved in the plane-parallel approximation. The feedback of Compton scattering on the structure of the post-shock region is also accounted for. A set of the post-shock region model spectra for various white dwarf masses is calculated.
Results. We find that Compton scattering does not change the emergent spectra significantly for low accretion rates or low white dwarf masses. However, it becomes important at high accretion rates and high white dwarf masses. The time-averaged, broad-band X-ray spectrum of intermediate polar V709 Cas obtained by the RXTE and INTEGRAL observatories is fitted using the set of computed spectral models. We obtained the white dwarf mass of and using models with Compton scattering taken into account and without it, respectively.
Key words: stars: binaries: close - stars: novae, cataclysmic variables - stars: white dwarfs - radiative transfer - scattering - stars: atmospheres
Cataclysmic variables (CVs) are close binary systems with a white dwarf (WD) as a primary (Warner 1995). A WD accretes matter from a companion star (typically a red dwarf) that fills its Roche lobe. This matter forms an accretion disc inside the Roche lobe of the WD. In intermediate polars (IPs), the significant magnetic field (106-107 G) can disrupt the disc at some distance from the WD surface. The matter is then free-falling along the magnetic field lines onto the WD and forms a strong shock above its surface (Aizu 1973). The temperature in the post-shock region (PSR) can be very high (10-100 keV) and the plasma there is cooled mainly via optically thin bremsstrahlung radiation in the X-ray band (Lamb & Masters 1979; King & Lasota 1979; Fabian et al. 1976). Cooling via Compton scattering is less significant due to a relatively low radiation energy density in PSR. This contrasts with the accretion columns of X-ray pulsars, where Comptonization is the dominant cooling mechanism (see Becker & Wolff 2007, for more details).
The theory of the PSR has been developed in a number of papers (Cropper et al. 1999; Aizu 1973; Wu et al. 1994; Woelk & Beuermann 1996). The most recent investigations by Canalle et al. (2005) and Saxton et al. (2005,2007) have studied the role of the two-temperature plasma and considered the dipole magnetic funneling (not the standard plane-parallel approximation). The plasma temperature in the PSR depends on the free-fall velocity at the WD surface; therefore, the X-ray spectra of IPs can be used to determine the WD mass (Rothschild et al. 1981). Such studies have been performed by Cropper et al. (1998), Cropper et al. (1999), Ramsay (2000), Beardmore et al. (2000), Revnivtsev et al. (2004), Suleimanov et al. (2005), and Falanga et al. (2005), using various X-ray data sets. Suleimanov et al. (2005) demonstrated that, to accurately measure the maximum post-shock temperature to infer the WD mass, the spectra at high energies (up to at least 100 keV) are needed.
The IPs are the most luminous and the hardest X-ray sources among accreting WDs. The interest in this class of objects has been growing in the past few years as IPs have been recently proposed as the dominant X-ray source population detected near the Galactic centre by Chandra observatory (Ruiter et al. 2006; Muno et al. 2004). The IPs also contribute significantly to the X-ray diffuse Galactic ridge emission (Revnivtsev et al. 2006), as well as dominate the hard X-ray sky there (Krivonos et al. 2007a). Moreover, most of the known CVs detected by INTEGRAL (Krivonos et al. 2007b; Bird et al. 2007) and Swift satellites are IPs (see e.g. Barlow et al. 2006; Bonnet-Bidaud et al. 2007; de Martino et al. 2008). In the strongly magnetized ( G) polar systems, cyclotron radiation is an important cooling mechanism, which suppresses the high temperature bremsstrahlung emission, whilst it should be negligible for the IPs. This could explain why the majority of the CVs observed in the hard X-ray band are IPs.
To better understand the hard X-ray production in IPs, one needs to model their observed spectral properties. The goal of the present paper is to investigate the influence of Compton scattering on the structure of the PSR and the emergent spectrum. We developed the physical model of the PSR and compare the computed spectral models to the broad-band X-ray spectrum of the intermediate polar V709 Cas obtained by RXTE and INTEGRAL observatories. We have also determined the WD mass in this system.
|Figure 1: A schematic view of the accretion flow onto the WD magnetic polar cap. The accreting gas falls quasi-radially towards the magnetic poles, forming a standing shock, below which the X-rays are emitted.|
|Open with DEXTER|
The structure of the stationary post-shock region (the accretion and emission geometry is illustrated schematically in Fig. 1) in a plane-parallel, one-dimensional geometry can be described (see e.g. Cropper et al. 1999; Suleimanov et al. 2005) by the mass continuity equation
The term 1+A is the amplification factor due to Comptonization. In contrast to the usually used integral amplification factor (which is just the ratio of the emergent flux to the incident flux of some seed soft photons, see Svensson 1984; Illarionov & Sunyaev 1975), we use here the local amplification factor, because we are interested in the vertical structure of the accretion column. In a homogeneous isothermal slab, the integral and local factors are nearly equal. As the Compton cooling rate depends on the local plasma temperature and density in a different way compared to the optically thin, thermal radiation, the local amplification factor (1+A) is a function of the geometrical coordinate z. The cooling due to cyclotron radiation is not considered here.
Equation (1) has the integral
A simple analytical solution of Eqs. (1)-(4)
can be obtained (see Frank et al. 2002), if one assumes a constant gas pressure throughout the PSR and considers the cooling only by bremsstrahlung.
The distributions of the velocity, temperature, and density along the PSR height z in this case are
We now describe the method of numerically solving the PSR structure. First of all, the input parameters, the WD mass, and the local accretion rate a are defined. Then we use two-loop iteration scheme. In the outer loop the amplification factor is changed, while in the inner loop the
PSR structure and the emergent spectrum are computed with the fixed
distribution of the amplification factor on the relative height z/z0, where z0 is the shock position.
First, an initial model without Comptonization (A(z)=0) is
calculated. Equations (1)-(4) are solved by the
shooting method from z=z0 to z=0 with the following boundary
conditions at z=z0 (see details in Suleimanov et al. 2005):
Then we again solve Eqs. (1)-(4) with the computed amplification factor. Here we assume that the amplification factor only depends on the relative position z/z0 in the PSR. It is important, because in the inner loop of iterations the shock position and the set of geometrical depth z change from iteration to iteration. Therefore, for the given shock position z0, the amplification factor at height z is computed by interpolating the dependence A(z/z0). This iteration procedure is repeated until the relative change in the shock position is less than 1 per cent (with 5-6 iterations being sufficient). As a result we obtain the self-consistent PSR model with Compton cooling, together with the emergent spectrum of radiation.
|Figure 2: Dependence of the Compton parameter for the PSR on the WD mass. The solid curve gives the results for the analytical model (7) of Frank et al. (2002) given by Eq. (21), while the dotted and dashed curves are for the numerical model of Suleimanov et al. (2005) for two accretion rates a = 1 and 10 g s-1 cm-2, respectively.|
|Open with DEXTER|
The radiative transfer equation for specific intensity
Compton scattering taken into account in the plane-parallel approximation
has the form
|Figure 3: Comparison of PSR models computed for with a local accretion rate a = 10 g s-1 cm-2 with (solid curves) and without (dashed curves) Compton scattering taken into account. The temperature and the density structures (panel a)), as well as the emergent spectra (panel b)), are shown. The dotted lines correspond to the analytical solution (7).|
|Open with DEXTER|
Before presenting the results of our calculations, let us first evaluate the importance of cooling by Compton scattering in the PSR. The Compton parameter
The influence of Comptonization on the temperature and the density structure of the PSR, as well as the emergent spectra, are demonstrated in Fig. 3. Here we take a 1.2 WD with a local accretion rate a = 10 g s-1 cm-2 and compute the PSR structure with and without Comptonization taken into account. The amplification factor 1+A due to Compton scattering is about 1.5 at the top of the PSR (near the shock) and decreases to 1 at the WD surface. Therefore, if Comptonization is taken into account, the cooling rate becomes higher, the height of PSR smaller, and the emergent spectrum slightly harder. But the difference between the emergent spectra with and without Comptonization is, however, rather small.
We computed the grid of PSR models for WD masses ranging from 0.3 to 1.3 with step 0.02 (51 models in total). The total accretion rate was assumed to be g s-1 with the fixed part f = 10-3 of the WD surface occupied by PSR. This gives different local accretion rates a for various WD masses, varying from 0.5 to 5 g s-1 cm-2. The emergent spectra of PSR with different local accretion rates (in the interval 0.1-10 g s-1 cm-2) are very similar (Suleimanov et al. 2005); therefore, we have only one free parameter, the WD mass, in our grid of model spectra. This means that the mass can be found by fitting the observed X-ray spectrum. Of course, the WD mass determination using spectral fitting is model dependent, because of using of the specific theoretical WD mass-radius relation and of the framework the theoretical spectra are calculated. However, the WD mass determination depends very little on the particular choice of and f.
The grid of the IP spectral models has been incorporated into the XSPEC software package under the name IPCOMP and is publicly available.
|Figure 4: The unfolded spectrum of V709 Cas and the best-fit spectral model of the PSR. The crosses represent the data from RXTE/PCA (3-20 keV) and INTEGRAL/ISGRI (20-100 keV) instruments. The total spectrum of the PSR model is shown by the solid (with Compton scattering) and the dashed (without Compton scattering) curves. The dotted curve represents the Gaussian iron line.|
|Open with DEXTER|
V709 Cas is a typical IP emitting in the hard X-rays (see e.g. Suleimanov et al. 2005). The present dataset was obtained using the RXTE/PCA (3-20 keV) (Jahoda et al. 1996) and INTEGRAL/IBIS/ISGRI (20-120 keV) (Ubertini et al. 2003; Lebrun et al. 2003) instruments with the total effective exposure times of 47 ks and 4.5 Ms, respectively. Our spectral model consists of the PSR model IPCOMP, an iron line modelled as a Gaussian, and the partial absorption model PCFABS described by the covering factor and the column density . The spectral analysis is performed with the XSPEC (Arnaud 1996) version 12.4.
The best fit with is obtained for , a partial covering fraction and cm-2 (see Fig. 4). The iron line energy is keV with an equivalent width of 508 eV. The unabsorbed flux between 0.5 and 150 keV is erg cm-2 s-1, which translates into a bolometric luminosity of erg s-1, for the source distance of 230 pc (Bonnet-Bidaud et al. 2001).
We can estimate the mass accretion rate using the simple relation . We get g s-1 for the WD radius of cm, which was computed using Eq. (10). To investigate the role of Compton scattering, we also implemented into XSPEC the grids of the PSR model without Compton scattering. The derived parameters are similar to those obtained with the Compton scattering model with a comparable reduced . The WD mass is found to be slightly lower . This seems counter-intuitive, because Comptonization increases the flux at high energies, and therefore would need lower temperature (and therefore lower WD mass) to radiate the same energy. However, the spectral slope in the 4-20 keV energy range is different for the two models, so the absorber parameters also change to and cm-2. Thus for V709 Cas, both mass estimations are consistent with each other within the errors, and Compton scattering does not have much effect, which is expected for such a WD mass and a rather low accretion rate.
For illustrative purpose, we also obtained a fit of the observed spectrum with a single temperature bremsstrahlung model. We obtained
cm-2, and the iron line energy
keV with the reduced
Using a linear approximation to the mass-radius relation (11), we can estimate the WD mass:
We have studied the influence of Compton scattering on the structure of the post-shock region and the broad-band X-ray spectra of intermediate polars. Compton scattering leads to slightly harder emergent spectra and smaller PSR height. The effect can be significant for the luminous intermediate polars with high-mass WDs. We have also constructed a grid of the emergent spectra in a wide range of WD masses and incorporated into the XSPEC package. We used this grid to fit the broad-band X-ray spectrum of V709 Cas observed by RXTE and INTEGRAL observatories and determined the WD mass in this object. Accounting for Compton scattering does not change the obtained mass significantly.
In this work we have not considered the cooling by cyclotron radiation. However, for WDs with stronger magnetic field (polars), Comptonization of the cyclotron photons can significantly affect the structure of the PSR and the emergent spectra. We plan to perform such a study in the future.
V.S. was suppoted by DFG (grant We 1312/35-1 and grant SFB/Transregio 7 ``Gravitational Wave Astronomy'') and by the President's programme for support of leading science schools (grant NSh-4224.2008.2). J.P. has been supported by the Academy of Finland grant 110792. M.F. acknowledges the French Space Agency (CNES) for financial support and thanks J. M. Bonnet-Bidaud for valuable discussions. V.S., J.P., and M.F. acknowledge the support of the International Space Science Institute (Bern), where part of this investigation was carried out.