A&A 402, 315-319 (2003)
R. E. Mennickent1, - K. Matsumoto 2 - M. Diaz 3
1 - Departamento de Física, Facultad de Ciencias Físicas y Matemáticas, Universidad de Concepción, Casilla 160-C, Concepción, Chile
2 - Graduate School of Natural Science and Technology, Okayama University, Okayama 700-8530, Japan
3 - Instituto Astronômico e Geofísico, Universidade de São Paulo, Brazil
Received 5 November 2002 / Accepted 14 February 2003
We describe the infrared spectrum of the galactic supersoft X-ray source MR Vel (RX J0925.7-4758) in the wavelength range of 1.1-2.6 m. We find Paschen, Brackett, He I, He II and neutral and ionized oxygen lines in emission. No traces of the secondary star were found. P-Cygni profiles were observed in Paschen and Brackett , from which we conclude that these lines arise in a wind. Fitting the optical and infrared energy distribution with a supercritical disk model, incorporating white dwarf and secondary star contributions, we were able to constrain the distance to 2-5 kpc, which implies an absolute magnitude MV between -2 and 0.
Key words: accretion, accretion disks - stars: emission-line - X-rays: stars - infrared: stars - stars: individual: MR Vel
The supersoft X-ray sources (SSS) are a class of luminous point sources 10 36-1037 ergs s-1, emitting mostly in the energy range of 20-60 eV. The most promising model for binary SSS invokes a white dwarf undergoing steady nuclear burning as a result of very rapid accretion ( /yr) from a Roche lobe-filling companion star (e.g. van den Heuvel et al. 1992). Whereas the bulk of the energy should be released near the white dwarf into the X-ray region, the optical and infrared spectra should be dominated by the other components of the system, namely the accretion disk and eventually the secondary star.
Most persistent SSS have been found in the Magellanic Clouds, and only few of them in the galaxy, since presumably their soft X-rays are easily absorbed by the interstellar hydrogen near the galactic plane. MR Vel (RX J0925.7-4758) is one of the few known galactic SSS and it has been relatively well studied in X-rays and optical wavelengths. The orbital period is unusually long for an SSS ( d, Schmidtke & Cowley 2001). During the orbital cycle, the star shows cyclic variations in the V magnitude between 17.1 and 17.3, whereas the B-V and V-R colours show no change (Schmidtke et al. 2000). The great strength of the interstellar absorption lines seen in the optical region are consistent with the position of the source behind the Vela Sheet molecular cloud. The inclination has been constrained to i= 55 10(Matsumoto & Mennickent 2000). The optical emission lines show a large velocity amplitude, implying masses of 1-2 for the donor star, which must be a giant in order to fill its Roche lobe, and 0.5-1.7 for the compact star (Bearda et al. 2002). The X-ray spectrum is rather complex, and cannot be explained by the usual models of photo-ionized plasma or stellar atmospheres combined with thermal plasma in collisional ionization equilibrium (Bearda et al. 2002; Motch et al. 2002).
In spite of the recent progress made to understand the unusual properties of MR Vel, the infrared region has not yet been explored. This is also true for most SSS, which have been historically studied in X-rays and the optical region. Infrared observations of SSS could, in principle, be combined with X-rays and optical data, in order to analyse the overall spectral energy distribution, aiming to constrain models for the emitting region and get insights on stellar parameters, energetics and distances. In this paper we provide the first description of the infrared spectrum of MR Vel, and combine this information with optical data to model the optical-IR spectral energy distribution.
The infrared spectroscopic observations reported in this paper were obtained at the ESO-Paranal Observatory with VLT-Antu using the ISAAC spectrograph in service mode. Spectra were taken, under photometric conditions, during the night of May 20, 2000 (UT). The exposure times were 120, 240 and 480 s for the J-band, H-band and K-band spectra, respectively. Sufficient overlapping was assured for composing a single spectrum ranging from 1.09 to 2.57 m with a FWHM resolution of 12 (J) to 27 (K) Å. The data were reduced using IRAF by first applying the combined dark and flatfield images as supplied by the ESO service mode operation group. Median sky frames were then combined for each object and spectral window. This process made use of jittered images where the object is located at different positions along the slit. Wavelength calibration was achieved by measuring the location of atmospheric OH emission lines (Rousselot et al. 2000) in the sky background. The flux calibration was performed using observations of the A0 standard HD 216009, made with the same instrumental setup (but with a wider slit) by the ESO operation team as part of the service mode program. Finally, the telluric absorption features were corrected with the aid of the absorption template constructed by dividing the spectrum of HD 216009 with a low order polynomial fit, excluding a few stellar absorption lines. The IRAF task "telluric" was used to find the best scale and shift factors which, when applied to the normalized telluric template, provided a reasonable correction of the telluric absorptions in HD 216009 and science exposures. This procedure worked well, except for the regions between 1.35-1.44 microns and 1.80-1.94 microns, characterized by heavy telluric absorption. These regions, corresponding to the ends of the J, H and K spectra, were excluded from the analysis and are not shown in this paper. Synthetic J, H, and K photometry of our calibrated standard star observations were compared to broadband photometry by Carter & Meadows (1995) showing differences below 0.08 mag.
Synthetic magnitudes of MR Vel resulted in J = 13.05, J-H = 0.76 and H-K = 0.53, with probable uncertainties of 0.1 mag. When using a colour excess of E(B-V) = 1.94 (Matsumoto & Mennickent 2000), we find unreddened colours of J-H = 0.15and H-K = 0.17, comparable to those observed in the SSS QR And (RX J0019.8+2156, Fender & Burnell 1996). Our finding likely indicates the presence of absorbing circumstellar material or a late-type companion. To our knowledge, these are the first infrared magnitudes and colours obtained for MR Vel ever. Therefore, it would be interesting to compare these magnitudes with future data in order to search for variability.
The combined spectrum reveals a steep blue continuum and weak emission lines (Fig. 1). A straight-line fit to the continuum in this region results in a slope of ergs s-1 cm-2 Å-1/m, which is incompatible with a power law. After a close inspection of the normalized spectrum (Fig. 2), we identified lines of neutral hydrogen (Paschen and Bracket), single ionized helium and highly ionized oxygen. These lines are listed in Table 1 along with relevant spectroscopic data. No evidence of a transient emission-line jet like those observed in H by Motch (1998) was detected. No photospheric absorptions, such as CO bands, are observed in the spectra. One interesting feature is the P-Cygni profile observed in the Paschen and Bracket lines (Fig. 3). The velocities found in the corresponding blueshifted absorption wings, relative to the emission maximum, reach up to -3000 (-1500) km s-1, being the absorption minimum at -755 (-680) km s-1, for Paschen (Brackett ). After inspection of the data, we realized that the difference between the maximum absorption velocity of these lines are real, and they cannot be due to an error in the continuum normalization.
|Figure 1: Flux calibrated spectrum of MR Vel. The vertical scale has been normalized to 2.03E-15 erg cm-2 s-1 Å-1.|
|Open with DEXTER|
|Figure 2: From up to down, J-band, K-band and H-band continuum normalized spectra for MR Vel. Emission lines are labeled.|
|Open with DEXTER|
In the picture of the steady nuclear-burning white dwarf model,
the mass accretion rate is limited in the narrow region
Since the critical accretion rate of the white dwarf
|Figure 3: A zoom around the P-Cygni profiles of Paschen and Brackett .|
|Open with DEXTER|
|Line||Intensity||EW (Å)||FWHM (km s-1)|
O II 11334-38
|He II 11628-77||2.3||-2||N/A|
|O VI 12617||2.0||N/A||N/A|
|O III 12692||2.2||N/A||N/A|
|He II 16926||1.5||-1||540|
O I 19463
|He II 19548||1.2||-3||520|
|He II 21891||0.9||-4||740|
|He II 23480||0.7||N/A||N/A|
A model of an accretion disk with a supercritical accretion rate was originally proposed by Abramowicz et al. (1988), as an optically thick and geometrically thick disk (e.g., Kato et al. 1998 for a review). Abramowicz et al. (1988) called the model a slim disk because of its medium thickness between a thin disk and a thick one. We refer it as a supercritical accretion disk, based on a physical viewpoint of the supercritical accretion rate.
The standard model for SSS indicates that we should consider a nearly critical accretion in SSS systems instead of the standard accretion disk (Shakura & Sunyaev 1973). It is also expected that a very luminous white dwarf in SSS strongly irradiates the accretion disk and eventually the companion (e.g. Popham & Di Stefano 1996; Schandl et al. 1997; Meyer-Hofmeister et al. 1997; Matsumoto & Fukue 1998; Meyer-Hofmeister et al. 1998). These effects in SSS were investigated by Fukue & Matsumoto (2001). In this work, we have calculated the spectral energy distribution (SED) for MR Vel, including the supercritical accretion disk and the irradiation effect, and have compared it with the observed spectrum in the range of optical-IR region after dereddening with E(B-V)= 1.94 (i.e. AV 5.9, Matsumoto & Mennickent 2000). The concept and formulation of the numerical model are described in Fukue & Matsumoto (2001), and we basically conform the manner in calculations to Fukue & Matsumoto (2001) in this paper.
The inferred intrinsic X-ray luminosity highly depends on the
model-atmosphere and gravity used in the fitting process,
but in any case, X-ray observations
and theoretical models have suggested that the white dwarf should be
extremely massive (Shimura 2000; Ebisawa et al. 2001),
although recent NLTE models of hot white dwarf atmospheres fail to
represent the complex X-ray spectrum (Motch et al. 2002).
In the following, as a working hypothesis, we assume
a blackbody SED for the central source associated with a
This means that we choose the WD
luminosity matching thecore-mass-luminosity relationship derived by
Iben & Tutukov (1996) for cold WDs accreting hydrogen:
An extremely small opening angle of the supercritical accretion disk () corresponds to a thin standard accretion disk. According to Hanamoto et al. (2001), is required to be more than 0.3 for a mass-accretion rate of expected for SSS. On the other hand, larger brings no significant increase of the luminosity in Rayleigh-Jeans tail of the SED. This can constrain upper and lower limits for the distance.
For an inclination angle in the range from 45 to 65, we found that the calculated SED based on the supercritical model described above suggests a distance to the source between 2 and 4 kpc for a wide range of (0.3-5). Even in the unlikely case of , the distance is constrained to 5-6 kpc. A representative fitting is shown in Fig. 4, which is calculated for a case of (corresponding to , , T= 7200 K) and . Besides the approximations and assumptions of our disk models, the differences between observed and theoretical spectra can be explained by a small error in the dereddening and/or the flux calibration for the optical NTT data, by photometric orbital variability or by the use of non-simultaneous optical-infrared data. From the above considerations we conclude that the distance for the source probably lies between 2 and 5 kpc.
|Figure 4: A representative SED model of MR Vel calculated for the supercritical accretion case with an inclination angle of . Two solid lines denote the contributions by the white dwarf (right) and the companion (left). Total contribution from the white dwarf, companion, and supercritical accretion disk are shown by five dashed lines, which correspond to 5, 3, 1, 0.5, and 0.3, from top to bottom. The de-reddened observed spectrum is over-plotted on the model SEDs by filled circles considering the best fitted case of d = 3 kpc. Optical data is from Matsumoto & Mennickent (2000) and infrared data is from this paper.|
|Open with DEXTER|
The presence of weak emission lines in the infrared (and optical) spectrum of MR Vel, the detection of highly ionized elements like He II and O VI, and the fact that the secondary star remains hidden in the IR spectrum, are consistent with a very luminous accretion disk. Are the accretion disks in longer orbital period systems more luminous? This is possible, since the thermal time scale for the secondary, regulating the mass transfer, decreases with the square of the donor mass. For instance, in QR And, whose orbital period is only 15.85 h, the He II 16926 and 21891 equivalent widths (Quaintrell & Fender 1998) are a factor 3 larger than in MR Vel, consistent with a fainter disk.
Interestingly, a P-Cygni profile has been observed in MR Vel at lines with very different excitation potential. The feature has been observed in optical H I lines (Matsumoto & Mennickent 2000; Schmidtke et al. 2000) and in the infrared (this paper). P-Cygni absorptions are also observed in X-ray lines like Fe XVII and O VIII (Bearda et al. 2002) with approaching velocities reaching a few thousands km s-1. It is not clear, at present, if the visibility of the P-Cygni profiles in the spectrum depends on the orbital phase. One may speculate that the less dense upper parts of the disk can be pushed away by the radiation pressure arising from the inner disk, producing an extended wind region where the P-Cygni absorptions arise.
It is interesting to note that the slope of the infrared continuum reported in Sect. 3 differs from the slope of an infinitely large steady state accretion disk ( Lynden-Bell 1969). This is likely related to a non-neglectable contribution of the secondary star to the overall flux in the infrared.
Our distance estimate is consistent with the 4 kpc value given by Hartmann et al. (1999) and only marginally consistent with the range derived by Motch et al. (1994), viz. 1-2 kpc. The distance derived using the accretion disk model fitting suggests an absolute magnitude between -2 and 0. Due to the nature of the assumptions considered in our modeling, this value should be taken with extreme caution. This is especially important, since in our model the secondary should contribute about 40% to the total light of the system in the infrared, however their spectral features are not visible in our spectra.
We thank the referee C. Motch, for useful comments on a first version of this paper. We also acknowledge J. Fukue who gave some advises in the calculation. This work was supported by Grant Fondecyt 1000324 and DIUC 202.01.030-1.0. MD acknowledges the CNPq support under grant # 301029.