A&A 462, 1091-1095 (2007)
DOI: 10.1051/0004-6361:20066354
A. Shporer1,
-
J. Hartman2 -
T. Mazeh1 -
W. Pietsch3
1 - Wise Observatory, Tel Aviv University, 69978 Tel Aviv, Israel
2 -
Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA
3 -
Max-Planck-Institut für extraterrestrische Physik, Giessenbachstrasse, 85741 Garching, Germany
Received 6 September 2006 / Accepted 26 November 2006
Abstract
Aims. Study the high-mass X-ray binary X-7 in M 33 using broad-band optical data.
Methods. We used recently published CFHT r' and i' data for variable stars in M 33 to extract the light curve of the optical counterpart of X-7. We combined these data with DIRECT B and V measurements in order to search for an independent optical modulation with the X-ray periodicity. The periodic modulation is modelled with the ellipsoidal effect. We used UBVRr'i' magnitudes of the system to constrain the temperature and radius of the optical component.
Results. The optical data revealed a periodicity of
days, which is consistent with the known X-ray period. Double modulation, which we attributed to ellipsoidal modulation, is clearly seen in four different optical bands. The absolute magnitude in six optical bands is most consistent with a stellar counterpart with
K and
.
We modelled the optical periodic modulation and derived the masses of the two components as a function of the orbital inclination and the radius of the stellar component. The resulting mass range for the compact object is
.
Conclusions. The system is probably a black hole HMXB, similar to Cyg X-1, LMC X-1 and LMC X-3.
Key words: stars: binaries: eclipsing - X-rays: individuals: M 33 X-7
M 33 X-7 (hereafter X-7) is a well known bright variable X-ray source
in the nearby M 33 galaxy (Long et al. 1981). Its variability was identified
as periodic by Peres et al. (1989), who derived a period of 1.7857 days and
were the first to suggest its XRB nature.
Using data from Einstein, ROSAT and ASCA, Larson & Schulman (1997) have
suggested that the actual period is 3.4531 days. This was confirmed
by Dubus et al. (1999), who reported a period of
days.
Pietsch et al. (2004) have analyzed X-ray data from XMM-Newton and Chandra and derived an improved period of
days. They also used optical B and V data to suggest the optical
companion is of spectral type between B0 I and O7 I, with masses of 25-35
,
making the compact binary component a black
hole. If this is true, X-7 is the first eclipsing black hole high mass
X-ray binary (HMXB).
In the course of the Chandra ACIS-I survey of M 33
(ChASeM 33, Sasaki 2005), the X-ray eclipse ingress and egress
of X-7 were observed for the first time (Sasaki et al. 2005).
Pietsch et al. (2006) combined Einstein, ROSAT and XMM-Newton data,
and derived an improved period, of
days and a
possible period decay rate of
yr-1. Pietsch et al. (2006) have also used HST WFPC2 images
to resolve the OB association HS 13 (Humphreys & Sandage 1980) in which X-7
resides, and have identified the optical counterpart to be an O6 III star whose mass is at least 20
.
Based on the optical
counterpart identification, possible orbital period decay rate, lack
of pulsations, X-ray spectrum and fits to the optical light curves of
Pietsch et al. (2004), Pietsch et al. (2006) suggested that X-7 is a black hole
HMXB.
In the previous studies the periodicity of X-7 was derived from X-ray data, and no independent analysis of optical data was performed. Optical light curves of X-7 were presented in B and V filters by Pietsch et al. (2004), who folded the optical measurements on the X-ray period. This paper presents an analysis of optical light curves of the X-7 counterpart obtained by Hartman et al. (2006) in a search for variable stars in M 33, performed at the Canada-France-Hawaii-Telescope (CFHT) with the Sloan r' and i'filters. Hartman et al. (2006) measurements are combined here with the B and V observations of Pietsch et al. (2004) to search independently for a periodic modulation. The combined data, as presented in Sect. 2 do show clear periodicity with the X-ray period, as presented in Sect. 3. Based on broad-band UBVRr'i' data from Pietsch et al. (2006), Massey et al. (2006) and Hartman et al. (2006) we apply a photometric modelling to derive the stellar temperature and radius of the optical component in Sect. 4.1, and periodic ellipsoidal model in Sect. 4.2. We briefly discuss our results in Sect. 5.
The present analysis is based on the optical measurements of X-7 in
four bands, obtained from two sets of observations carried out
4 years apart. The B and V light curves of
Pietsch et al. (2004, Table 4) were obtained from a special
photometric analysis of
DIRECT
(Mochejska et al. 2001) images, taken at the Kitt Peak National
Observatory 2.1-m telescope in 1999. The V light curve includes 70 measurements and the B light curve only 30 points. Sloan r' and i' light curves were taken from the publicly available data of the
M 33 CFHT variability
survey
(Hartman et al. 2006, object 243 718). Their data were obtained at the
CFHT 3.6-m telescope atop Mauna Kea on two consecutive observing
seasons in 2003 and 2004. An i' finding chart of X-7 is given in
Fig. 1. A comparison with the finding chart of Pietsch et al. (2006,
their Fig. 5) shows this object is at the position of X-7
optical counterpart identified by Pietsch et al. (2006). The
identification of the same period in the optical as in the X-rays
(presented hereafter) confirms the identification. The i' light
curve contains 32 measurements, where we ignored two measurements with
uncertainty larger than 0.04 mag. The r' light curve consists of 31 measurements, where we ignored three obvious outliers.
![]() |
Figure 1:
X-7 finding chart, in i', from the CFHT data. Position of the
variable object is indicated by an arrow. Image FOV is
![]() |
Open with DEXTER |
To search for periodicity we have applied a multi-band double-harmonic
fitting (Shporer & Mazeh 2006) to the four available light curves together,
which included 158 individual measurements. For each trial period we
fitted the entire data set with a double-harmonic function, with
independent zero points for each of the four light curves. The
periodogram value for each trial period was taken as the power (sum of
squared harmonic coefficients) divided by the goodness-of-fit
parameter of the fitted light curve.
Figure 2 presents the resulting periodogram. The
strongest peak corresponds to a period of
=
days. The
second strongest peak is exactly at the first harmonic of
.
The third strongest peak is at a frequency corresponding to
the lunar cycle. The optical period agrees well with the known X-rays
period,
days, of Pietsch et al. (2006).
![]() |
Figure 2: Multi-band double-harmonic periodogram. |
Open with DEXTER |
![]() |
Figure 3: Folded light curves, in magnitude, in ( top to bottom) i', r', V and B. All light curves are folded on the period identified here, of 3.4530 days, and phase zero was taken to be at HJD = 2 453 639.119, which is time of mid-eclipse in the X-rays (Pietsch et al. 2006). |
Open with DEXTER |
It appears as if the amplitude of the modulation decreases with
decreasing wavelength, from i' to B, even when ignoring the four
faintest points in i' and r'. The origin of this varying modulation
is not clear. It could result from some real wavelength dependence of the
system's modulation. On the other hand, the
different modulation could
come either from some secular changes during the 4 years
between the KPNO B and V and CFHT i' and r' observations, or
from some wavelength-dependent blending of light from a near-by star.
Table 1: Harmonic fits amplitude and minimum phase near 0.0.
In our modelling of the X-7 binary system, we preferred to find first a stellar model that fits the broad-band photometry of the system, assuming the optical brightness is coming from the optical star alone. From the stellar radius, the observed X-ray eclipse width and the periodic ellipsoidal modulation we can derive the two masses as a function of the orbital inclination.
We have used the following broad-band magnitudes: U =
18.1, B = 18.8, V = 18.9 (Pietsch et al. 2006, based on HST data),
R = 19.0 (Massey et al. 2006, LGGS),
r' = 19.1 and
i' = 19.42 mag (Hartman et al. 2006), assuming an
uncertainty of 0.1 mag on all values. We converted the observed
measurements to absolute magnitudes by using an M 33 distance modulus
of 24.62 mag (Freedman et al. 2001),
mag from
Pietsch et al. (2006) and
values from Schlegel et al. (1998).
We compared these magnitudes to the ones determined from the
integrated magnitudes from Girardi et al. (2004,2002), where we
assume
(Girardi et al. 2002) and fit for the
radius for each value of
and
.
We did this for
both [Fe/H] = 0 and [Fe/H] = -0.5 tables, since assuming X-7 is at a
galactocentric distance of 2 kpc results in [Fe/H] = -0.22, based on
the [O/H] gradient of Garnett et al. (1997) and the [O/H] to [Fe/H]
conversion of Maciel et al. (2003).
The data allow a large range of stellar models. One such model, for
[Fe/H] = 0, with
,
has
K,
and
.
For [Fe/H] = -0.5, with
,
we get
K,
and
.
The
95% uncertainty range, where
increases by up to 8.0 (for
three parameters), is:
K,
.
The value of
is unconstrained. Figure 4
brings three representative models; the solid line presents the model
with
K,
,
the dotted line
represent the model with
K,
and
the dashed line shows the model with
K,
.
The figure shows that the absolute magnitudes of the
X-7 optical counterpart can not constrain the stellar radius, and any
value in the range of 15-35
is possible. However, as we
show in the next section, a stellar radius of
is
very unlikely, since the star will be too massive. Therefore, our most
likely models are in the range of
and
K.
![]() |
Figure 4:
Photometric broad-band modelling of the X-7 stellar
component. Error bars present the assumed uncertainty of 0.1 mag
in all bands. Solid line presents a model with
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 5:
BVr'i' light curves phased at
P = 3.453014 day and
T0 = 2453639.119 from Pietsch et al. (2006), with model light curves
for
![]() ![]() ![]() |
Open with DEXTER |
We used two different codes, PHOEBE and ELC, to model the obtained periodic light curves and derive an estimate of the masses of both components. For both codes, the main effect is due to the tidally induced ellipsoidal shape of the optical component.
The first program is the "PHysics Of Eclipsing BinariEs''
program of Prsa & Zwitter (2005), a front-end code for the Wilson-Devinney
program (Wilson & Devinney 1971; Wilson 1979,1990). Following
Pietsch et al. (2006), who used PHOEBE to model the X-7 B and V light
curves, we used the eclipse half angle value of
to find the ratio of the radius of the optical component to the
semi-major axis as a function of inclination. We varied the
inclination between
and about
and fit for the
mass ratio and a luminosity scale to match each light curve, assuming
a temperature of 27 000 K and linear limb-darkening
coefficients from Claret (2000,2004), assuming a metallicity of
[Fe/H] = -0.2.
We performed the fits ignoring the four faintest points in each of the r' and i' light curves which occur near phases 0 and 0.5 - no model was able to match these points. Including those points does not substantially change the mass-ratio, while it does slightly affect the luminosity scale.
There is not enough data to constrain the stellar radius or the
inclination of the system; a good fit to the light curves could be
found for large ranges of these two parameters. One such
model is plotted in Fig. 5. The PHOEBE light curve is
drawn with a dotted line and the other model with a dashed line. For
both, values of
,
K and
were taken, although changing these values results in
essentially no difference for the model light curve.
The range of allowed inclinations and stellar radii did not enable us
to constrain the masses of the two components and they could only be
derived as a function of the assumed orbital inclination and optical
radius. We found that for the PHOEBE code the mass ratio (
)
varies from 0.29 at
to 0.15 at
.
Below
the optical component would exceed its Roche
lobe. Assuming
,
the mass of the optical component
would vary between
and
,
while that of
the compact object would vary between
and
.
Because the fits are fairly insensitive to the
temperature of the optical component, the resulting masses scale as
the radius of the star to the third power. For
,
the
mass of the optical component varies between
and
,
while the compact object mass varies between
and
.
We note, however, that assuming a
stellar radius larger than
results in an unrealistic
large stellar mass. We therefore limit the discussion to stellar radii
smaller than
.
The second program used was the ELC code of Orosz & Hauschildt (2000). This
program has been specifically tailored to model optical light curves
of X-ray binary systems, including effects such as the presence of an
accretion disk, X-ray heating of the optical component and the ability
to explicitly take one of the components to be invisible. To model the
X-ray heating we assumed an average X-ray luminosity of 1037.5 erg s-1, and to model the disk we assumed an inner disk
temperature of 107 K (Pietsch et al. 2006). We assumed a disk opening
angle of
,
and a power-law exponent for the temperature
profile of the disk of -0.75. We used the black-body mode of the
model rather than using the detailed model atmosphere since the model
atmosphere would not necessarily be a good match to the Sloan r' or
i' filters. We proceeded as above, stepping through a range of
inclination angles for two values of the optical radius and fitting
for the mass ratio. This time we also fit for the inner and outer disk
radii. One possible model, for
,
K
and
,
is plotted in Fig. 5.
The resulting mass ratio varies from 0.23 at
inclination to 0.07 at
inclination. For a radius of 20
this
corresponds to a mass range of 83 to 42
for the optical
component and a range of 19 to 3
for the compact
object. Below
the ELC model star would exceed its Roche
lobe. In units of the inner Lagrange radius of the compact object, the
inner disk radius varies from
to
while the outer disk radius varies from 0.61 to 0.81.
The masses of both components as a function of inclination, derived by both codes for two values of the optical radius, are plotted in Fig. 6.
![]() |
Figure 6:
Plots showing the mass vs. inclination for two values of
the stellar radius. The left plot is for the optical component, the
right plot is for the compact object. The solid and the dashed
lines show the results for
![]() ![]() |
Open with DEXTER |
We have conducted a broad-band optical analysis of the X-7 optical
counterpart, including period analysis, photometric modelling and light
curve modelling. We applied a multi-band double-harmonic period
analysis of X-7 optical light curves in four bands and identified a
period of
=
days. This period, derived from an
independent optical analysis, is in good agreement with the known
X-rays period (Pietsch et al. 2006).
Using photometric absolute magnitudes in six bands, we find a range of
models for the optical counterpart, with
and
K. However, stellar radius larger than
results in an unrealistic large stellar mass. We
therefore suggest that the stellar radius is smaller than
.
The possible temperature range is reduced accordingly
to
.
The classification of the optical star
as O6 III by Pietsch et al. (2006) is well within our uncertainties.
We have modelled the optical light curve using two programs. While the ELC model, incorporating both a disk and X-ray heating matches better to the data than the PHOEBE model (see Fig. 5), neither model yielded a good fit to the light curves. In particular, the minima in r' and i' appear to be underestimated by all models that we have tried. At this point it is unclear if this is due to random or systematic errors in the observations or to inadequacies in the physical models.
The simplistic models can nevertheless yield a rough estimate of the
masses of the two components as a function of the optical radius and
the orbital inclination. For the smallest likely optical radius
(
)
we get
and
from the
PHOEBE and ELC codes respectively. Therefore, the present analysis can
not exclude completely a neutron star as the compact object, although
this option is very unlikely.
On the other hand, the similarity between M 33 X-7 and the three known
black-hole HMXB - Cyg X-1, LMC X-1 and LMC X-3
(e.g., McClintock & Remillard 2003; Cowley 1992), is striking. All four systems
have orbital periods of a few days, their optical counterpart is an
early-type star and the compact object has a mass of 6-10 .
Radial velocity measurements of M 33 X-7 will enable us to
better constrain the masses of the two components. Large telescopes
and presently efficient spectrographs render such measurements
feasible. Together with the X-ray eclipse and the optical light
curves, we should be able to understand this system better than any
other black-hole HMXB.
Finally, if this system would be proven to be a black-hole HMXB, it is interesting to note that three out of the four such known systems were found in external galaxies. This might indicate that we are not very efficient in detecting such systems in our own galaxy. Maybe the dust in the Galactic plane hides from our telescopes many more Galactic black holes.
Acknowledgements
We are grateful to J. McClintock for helpful discussions and to J. Orosz for sharing his ELC program with us. We thank the anonymous referee for his helpful comments. This work was partially funded by the German-Israeli Foundation for Scientific Research and Development and by the Israeli Science Foundation. This research has made use of NASA's Astrophysics Data System Abstract Service and of the SIMBAD database, operated at CDS, Strasbourg, France.