A&A 445, 1037-1040 (2006)
DOI: 10.1051/0004-6361:20053303
J. M. Bonnet-Bidaud1 - M. Mouchet2,3 - D. de Martino4 - R. Silvotti4 - C. Motch5
1 - Service d'Astrophysique, DSM/DAPNIA/SAp, CE Saclay, 91191 Gif-sur-Yvette
Cedex, France
2 -
LUTH (CNRS-UMR 8102), Observatoire de Paris, Section de Meudon, 92195 Meudon Cedex, France
3 -
APC (CNRS-UMR 7164), Université Denis Diderot, 2 place Jussieu, 75005 Paris, France
4 -
INAF - Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131 Napoli, Italy
5 -
Observatoire Astronomique de Strasbourg, 67000 Strasbourg, France
Received 25 April 2005 / Accepted 20 July 2005
Abstract
We report the first time-resolved photometric and spectroscopic optical observations of the
X-ray source RX J2133.7+5107, identified in the ROSAT survey.
A clear persistent optical light pulsation is discovered with fast photometry
at a period
of
s which we associate
with the spin period of an accreting white dwarf.
Radial velocity curves of the strong emission lines show modulation with
a period of
h, identified as the orbital period.
These observations establish that the source is a member of the intermediate
polar class (IPs) of magnetic cataclysmic variables. With only 4 IPs with longer
orbital periods, RX J2133.7+5107 is among the widest systems. It is a unique IP with an orbital period in the middle of the
so-called (6-10) h IP gap and it shows a significant degree
of asynchronism with a ratio
/
of 0.02.
When attributed to the motion of the white dwarf, the emission lines orbital
modulation yields a mass function of
which, for a probable
inclination
and a white dwarf mass
,
corresponds to a
secondary mass
.
Key words: stars: binaries: general - stars: individual: RX J2133.7+5107 - stars: novae, cataclysmic variables - accretion, accretion disks - stars: variables: general - white dwarfs
The X-ray source RX J2133.7+5107 was discovered from the ROSAT Galactic Plane Survey and identified as a relatively bright (
)
cataclysmic variable (CV) (Motch et al. 1998, hereafter M98). Its X-ray hardness ratios indicate a hard X-ray source, similar to the intermediate polar (IPs) class among cataclysmic variables (CVs).
IPs are Magnetic Cataclysmic Variables (MCVs), where a white dwarf (WD) with a relatively strong (
MG) magnetic field accretes material from a late-type secondary companion overflowing its Roche-lobe.
The accretion flow is channelled towards the polar regions where a
strong shock develops above the WD surface (see Patterson 1994, for a review), which means that the signature of channelled accretion in these systems is a strong pulsation at the WD rotational period.
In an effort to elucidate the nature of RX J2133.7+5107, we obtained optical observations to search for a signature of the WD rotation and its orbital motion.
![]() |
Figure 1:
RX J2133.7+5107 photometric light curves in relative white light amplitudes with a 30 s resolution. A ![]() |
Open with DEXTER |
The photometric time series were analyzed performing a Discrete
Fourier Transform (DFT) of the relative amplitudes, after detrending from the low frequency
underlying variations. A significant power was
present in two groups of frequencies inside the overall range (140-160) day-1 (see Fig. 2). This structure
closely resembles the one observed in the optical power spectra of many IPs with the presence of the spin frequency ()
and its orbital sideband (
), where
is the orbital frequency.
The same structure was also found, though with much lower intensity, at twice these
frequencies, which are identified as the first harmonics (see Fig. 2, top panels).
Given the strong influence of the data window, a number of aliases were also introduced
that prevent any clear identification of the
true peaks. To remove this effect, we used the CLEAN algorithm (Roberts et al. 1987),
and results are also shown in Fig. 2 (right panels). The most significant peaks were found at
day-1
with its orbital sideband at
day-1and harmonics at 2
day-1, where the error bars were determined from the peaks'half widths.
The fundamental and the first harmonic of the spin
match strictly, while for the beat frequency, the weak peak observed at
day-1,
is slightly offset from the expected beat harmonic (
)
frequency. This discrepancy is most likely
due to the very small amplitude (full amplitude 0.8
)
and imperfect removal of the data window effect.
![]() |
Figure 2:
Power spectrum of the photometric data (pre-whitened from
the low frequency trend) around the spin fundamental ( bottom)
and first harmonic ( top) frequencies. Note the different power
scale for the two ranges of frequencies. Left panels: direct
Fourier transform of the time series, central panels: synthetic
power spectrum using the best fit periods (see text).
Right panels: time series "clean'' power spectra after removing the
effect of the data window. The location of the fundamental (![]() ![]() ![]() |
Open with DEXTER |
![]() |
(1) |
We checked the consistency by producing a synthetic
DFT of the composite sinusoidal function, including the four most significant
frequencies above.
Results are shown in Fig. 2 (central panels).
The good agreement between the observed and synthetic spectrum
confirms that the spin and beat modulations have been determined accurately.
The data, pre-whitened from the
low frequency trend, were folded at the derived spin period and the pulse shape
is shown in Fig. 4. The pulse has a full amplitude of 8
and a quasi-sinusoidal shape with a less steep rise than the decay, indicative
of the presence of the first harmonic.
From the photometric data, the measured beat frequency yields an
orbital frequency of
day-1 corresponding to
h, in agreement with the orbital
period found from spectroscopy (see below).
The limited number of nights and the probable
contamination from a likely non-periodic variability do not allow a more precise
determination of the orbital period.
Standard reduction was performed in the ESO-MIDAS package, including cosmic removal, bias subtraction, flat-field correction, and wavelength calibration. The wavelength calibration was checked on sky lines that were found within 1.6 Å and 0.7 Å from their expected wavelengths in 1998 and 2004, respectively. All radial velocity measurements were corrected from the Earth motion and from the small instrumental shifts measured on the OI 5577 Å line and times were converted in the heliocentric system. Flux calibration was performed using the standards BD+28 4211 and BD+33 2642.
The mean optical spectrum of RX J2133.7+5107 (Fig. 5) is typical of magnetic CVs with strong emission lines
of the Balmer series, HeII (4686 Å and 5411 Å) and CIII-NIII (4655 Å), superposed on a relatively blue continuum.
In 2004, when the flux calibration was reliable, the mean B flux was estimated at
,
which corresponds to a B magnitude of
15.8, similar to what was reported in the identification paper by M 98. The characteristics of the main emission lines are given in Table 3.
Similar equivalent widths were observed in the two observations in 1998 and 2004, but the intensities of the lines change by a factor 2.
The radial velocities were measured for the main lines. They showed a significant modulation of (30-70)
amplitudes that were analysed to constrain the orbital period. The exposure of the individual spectra (30 min.) were close to a multiple of the spin period so that no strong influence from the spin modulation was expected, but the period search was significantly affected by the 1-day sampling of the observation. The best period was determined from an
periodogram where the radial velocity measurements were folded at trial periods and tested against a sine modulation. The most accurate determination was obtained for the strongest line (H
) in the longest observation of 2004 (Fig. 6). The minimum
corresponds to a period of
h, where the error bar is at a 99% level, computed for two independent parameters. A secondary minimum, also observed at
10.3 h, is a 1-day alias that can be excluded at a level better than
from an F-test.
A similar value was obtained independently in 1998 by making use of the velocities of the strongest H
line, but the accuracy on the period value does not allow us to keep track of the phase between the two observations, so the orbital period cannot be refined. The orbital ephemeris for RX J2133.7+5107 is determined as:
![]() |
(2) |
![]() |
Figure 3:
Radial velocities of the main lines folded with the orbital period
![]() |
Open with DEXTER |
Table 1: The parameters of the orbital modulation of the emission lines' radial velocities.
With a longer period than 7 h, a secondary of a relatively early type
(up to K0) can be expected.
Such stars usually show a significant sodium line Na D (5889-5896 Å) in absorption.
This line is indeed detected in the spectrum of RX J2133.7+5107, partly blended
with the HeI (5875 Å) line, at the resolution used
in our observations (see Fig. 5). However the Na D line observed
in RX J2133.7+5107 is at least partly interstellar,
indicative of both significant distance and absorption. This is confirmed
by the presence in the spectrum of the diffuse bands at DIB(4430)
and DIB(
5780) similar to what is observed in some other CVs
(see Thorstensen & Taylor 2001).
The equivalent widths of the (
4430) and (
5780) DIBs are measured at 0.7 Å and 0.3 Å respectively,
and are both consistent with a reddenning of
(Herbig 1975; Snow et al. 2002). If purely interstellar, the
Na D1 (
5890) component would have indicated a much higher reddenning
of
(Munari & Zwitter 1997), inconsistent with
the above value and indicating that the secondary also contributes to the line.
From the measured equivalent widths, the star contribution can be estimated to be less than 50%.
The Na D equivalent widths were measured on representative stars
from the Jacoby et al. (1984) catalogue with spectral types from K0 to M0
and compared to the EW upper limit measured in RX J2133.7+5107.
If the secondary contributes for less than half of the observed Na D,
the contribution of a K0, K5, or M0 star to the optical flux is respectively less than 39%, 37%, and 9%, which translates into a lower limit to their
distance of 1200, 650, and 600 pc respectively, assuming standard
absolute magnitudes (Allen 1970).
X-ray observations of the source by ROSAT do not give better indications
of absorption and distance.
From the observed ROSAT hardness ratios and the
measured value
(M 98),
using the
distribution for CVs in Fig. 4 (M 98) gives a
rough estimation of
with, however, a large uncertainty.
The measured ROSAT flux and
the count conversion given in M 98 yield a [0.1-2.4 keV] luminosity of
(d/500 pc)2
.
A complementary constraint on the system can also be obtained from the
velocity amplitude of the emission lines, if they are assumed
to be associated with the white dwarf. The mean velocity of the Balmer
lines yields a mass function of
.
For an inclination
of
,
suggested by the absence of double-peaked lines,
and a white dwarf mass
,
the secondary mass is
,
consistent with a star at least more massive than M2V.
The relatively long orbital period of RX J2133.7+5107 makes it fall into what has been called the "IP gap'' (Schenker et al. 2004), a lack of IPs into the (6-10) h orbital period interval. From the list of 39 IPs in Norton et al. 2004, all are either shorter than 6.14 h or longer than 9.87 h, this last value being the one newly determined for 1RXS J1548-4528 by de Martino et al. (2005). Though it could still be due to an observational bias, this may indicate a particular state of evolution.
With a spin-to-orbit period ratio of 0.02, RX J2133.7+5107 is also a
relatively highly asynchronous system with only 8 IPs showing lower values.
This suggests a young system still far from equilibrium.
For magnetic CVs, Norton et al. (2004) have computed the
equilibrium rate when the spin-up accretion torque is balanced by
the braking effect of the magnetic torque, using a detailed
plasma-magnetic field interaction.
For a value of
and an
orbital period of 7.2 h, the equilibrium is reached for a
magnetic moment of
G cm3, here assuming a mass ratio q=0.5,
but the variation is less than 10% for q in the range
.
The magnetic moment is probably much lower in RX J2133.7+5107, since accretion will only
take place if the Alven radius (
)
is lower than the corotation radius (
)
(see Warner 1995).
For RX J2133.7+5107,
cm (for
in
the range (0.6-1.0)
). For
G cm3,
the (
)
condition will require an accretion rate that is higher than
g s-1, so much too high for the observed
X-ray luminosity.
Assuming a more canonical value of
g s-1, the magnetic
moment then has to be lower than
G cm3; therefore RX J2133.7+5107
is most probably an IP with a weak magnetic field.
A definite confirmation of its IP nature will be possible with
the detection of strong X-ray pulses expected at the 571 s period.
Note added in proofs. After the completion of this paper, another long period IP, 1RXSJ173021.5-055933, was found by Gaensicke et al. (2005), MNRAS, in press, [arXiv:astro-ph/0504460]. With an orbital period of 15.42 h and a spin of 128 s, it also shows a very high degree of asynchronism.
Table 2: Log of photometric and spectroscopic observations.
Table 3: Equivalent widths, FWHM and fluxes of the strongest lines.
![]() |
Figure 4: Pulse shape of the 571 s spin modulation. A significant distortion from a pure sine shape is visible due to the presence of a first harmonic. |
Open with DEXTER |
![]() |
Figure 5: The mean optical spectrum of RX J2133.7+5107 obtained in July 2004, showing typical strong Balmer and Helium emission lines. Note the significant CIII-NIII and the presence of diffuse interstellar bands (DIB) and Na absorption, blended with the HeI (5875 Å)emission. |
Open with DEXTER |
![]() |
Figure 6:
![]() ![]() ![]() ![]() |
Open with DEXTER |