A&A 394, 983-991 (2002)
DOI: 10.1051/0004-6361:20021238
M. Ribó1 - E. Ros2 - J. M. Paredes1 - M. Massi2 - J. Martí3
1 - Departament d'Astronomia i Meteorologia, Universitat de Barcelona, Av. Diagonal 647, 08028 Barcelona, Spain
2 - Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
3 - Departamento de Física, Escuela Politécnica Superior, Universidad de Jaén, Virgen de la Cabeza 2, 23071 Jaén, Spain
Received 2 August 2002 / Accepted 27 August 2002
Abstract
In an attempt to increase the number of known microquasars, Paredes et al.
(2002) have presented a long-term project focused on the search
for new objects of this type. They performed a cross-identification between
X-ray and radio catalogs under very restrictive selection criteria for sources
with
,
and obtained a sample of 13 radio-emitting X-ray sources.
Follow-up observations of 6 of these sources with the VLA provided accurate
coordinates, which were used to discover optical counterparts for all of them.
We have observed these six sources with the EVN and MERLIN at 5 GHz. Five of
the six objects have been detected and imaged, presenting different
morphologies: one source has a two-sided jet, three sources have one-sided
jets, and one source is compact. With all the presently available information,
we conclude that two of the sources are promising microquasar candidates in
our Galaxy.
Key words: X-rays: binaries - radio continuum: stars
Microquasars are stellar-mass black holes or neutron stars that mimic, on
smaller scales, many of the phenomena seen in AGN and quasars. Microquasars
have been found in X-ray binary systems, where a compact object accretes
matter from a companion star. Radio Emitting X-ray Binaries (REXBs) with
relativistic radio jets, like SS 433, GRS 1915+105,
GRO J1655-40 or Cyg X-3, are good examples of microquasars
(see Mirabel & Rodríguez 1999 for a detailed review), while
other well known sources like LS I +61 303 could turn out to be new
relativistic jet sources (Massi et al. 2002). With the recent
addition of LS 5039 (Paredes et al. 2000),
Cyg X-1 (Stirling et al. 2001) and
XTE J1550-564 (Hannikainen et al. 2001) to the
microquasar group, the current known number of this kind of sources is 14,
among 50 REXBs within the
280 known X-ray binaries (Liu et al.
2000; Liu et al. 2001). Recent studies of microquasars can be
found in Castro-Tirado et al. (2001).
As pointed out in Sect. 1 of Paredes et al. (2002, hereafter
Paper I), the number of known microquasars remains still small, especially
when trying to study them from a statistical point of view. An interesting
aspect of microquasars is their possibility of being related to unidentified
high-energy -ray sources, as suggested by Paredes et al.
(2000). Moreover, some parameters like the jet velocity, seem to
be related to the mass of the compact object (i.e., with its potential well).
Nevertheless, the lack of meaningful statistical studies because of the small
population of microquasars with known jet velocities, prevents any definitive
statements of this kind being made. Therefore, it is worth searching for new
microquasars in order to increase the known population.
![]()
|
In this context, the sources to search for are REXBs, which are in fact
microquasar candidates. To this end, a cross-identification between the X-ray
ROSAT all sky Bright Source Catalog (RBSC) (Voges et al. 1999) and
the NRAO VLA Sky Survey (NVSS) (Condon et al. 1998) was made for
sources with
under very restrictive selection criteria, and the
obtained results have been presented in Paper I. A sample containing 13
sources was obtained, and 6 of them were observed at radio wavelengths in
order to obtain radio spectra and their variability, as well as accurate radio
positions, which allowed the authors to discover the corresponding optical
counterparts for all of them. At the end of this study, two of the sources,
namely 1RXS J001442.2+580201 and 1RXS J013106.4+612035 were
classified as promising microquasar candidates. The remaining four sources
were classified as weaker candidates for a number of reasons.
1RXS J042201.0+485610 shows a highly inverted spectrum at high radio
frequencies and an optical counterpart slightly extended.
1RXS J062148.1+174736 is an object extended in the optical.
1RXS J072259.5-073131 presents a one-sided radio jet at arcsecond
scales, supporting the possibility of being extragalactic. And finally,
1RXS J072418.3-071508 is a known quasar.
In this paper we present Very Long Baseline Interferometry (VLBI) observations of these six sources, aimed at revealing possible jet-like features at milliarcsecond scales. We describe the observations and the data reduction in Sect. 2, present the results and a discussion in Sect. 3, and we summarize our findings in Sect. 4.
We observed the six sources studied in Paper I simultaneously with the Multi-Element Radio-Linked Interferometer Network (MERLIN) and the European VLBI Network (EVN) on February 29th/March 1st 2000 (23:30-23:05 UT) at 5 GHz. Since some of the target radio sources were faint, we scheduled the observations introducing phase-reference calibrators with cycle times of around 7 min (compatible with the expected coherence times). Apart from the calibrators presented in Table 1, 4C 61.02 was used as calibrator for 1RXS J013106.4+612035, and TXS 0422+496 and NRAO 150 for 1RXS J042201.0+485610. We also observed the fringe-finder DA 193 and the MERLIN flux density calibrator 3C 286. Single dish flux density measurements were carried out with the MPIfR 100 m antenna in Effelsberg, Germany.
MERLIN is a connected radio interferometer across England, with baselines reaching up to 217 km length. This array observed with 2-bit sampling at dual polarisation with two blocks of 16 channels, each channel of 1 MHz bandwidth. We analysed the left hand circular polarisation data excluding one channel at both edges of the band, yielding a final bandwidth of 14 MHz. The correlator integration time was of 4 s.
The MERLIN data reduction was carried out at Jodrell Bank Observatory, using standard procedures within the Astronomical Image Processing System ( AIPS, developed and maintained by the US National Radio Astronomy Observatory). We did not detect 1RXS J042201.0+485610. All other sources were detected, and accurate positions for the target radio sources were obtained via phase-referencing. These positions, presented in Table 1, were used later as a priori information for the VLBI correlation. The position given in Table 1 for the quasar 1RXS J072418.3-071508 was obtained from the VLA observations presented in Paper I. The position of 1RXS J072259.5-073131 is deduced from the phase-reference offset relative to 1RXS J072418.3-071508 provided by MERLIN.
To image the sources, we averaged the data in frequency and exported them to be processed into the difference mapping software DIFMAP (Shepherd et al. 1994), where we time-averaged the data in 32 s bins after careful editing.
The EVN observations were performed with the following
array (name, code, location, diameter):
Effelsberg, EB, Germany, 100 m;
Jodrell Bank, JB, UK, 25 m;
Cambridge, CM, UK, 32 m;
Westerbork, WB, The Netherlands,
m;
Medicina, MC, Italy, 32 m;
Noto, NT, Italy, 32 m;
Shanghai, SH, China, 25 m;
Torun, TR, Poland, 32 m;
and
Onsala85, ON, Sweden, 25 m.
Data were recorded in MkIV mode with 2-bit sampling at 256 Mbps with left hand
circular polarization. A bandwidth of 64 MHz was used, divided into 8
intermediate frequency (IF) bands.
The data were processed at the EVN MkIV correlator at the Joint Institute for VLBI in Europe (JIVE), in Dwingeloo, The Netherlands. The correlator integration time was of 4 s. A first post-processing analysis was also carried out at JIVE. The data were processed using AIPS. A first a priori visibility amplitude calibration was performed using antenna gains and system temperatures measured at each antenna. The fringe fitting ( FRING) of the residual delays and fringe rates was performed for all the radio sources. No fringes were found for 1RXS J042201.0+485610 and 4C 61.02. Fringes for many baselines were missing for 1RXS J001442.2+580201, 1RXS J062148.1+174736 and TXS 0422+496.
To improve the fringe detection on all baselines for 1RXS J001442.2+580201 and 1RXS J062148.1+174736 we
used the delay, rate, and phase solutions from their corresponding
phase-reference calibrators (J0007+5706, 118
separation,
and J0630+1738, 1
59
separation)
and interpolated them to the target sources using the AIPS task
CLCAL. We fringe-fitted the target sources again using narrower search
windows and obtained solutions for all baselines. A similar attempt on 1RXS J042201.0+485610 (with respect to TXS 0422+486 and NRAO 150) was unfruitful.
Effelsberg was used as reference antenna throughout the AIPS data
reduction process.
We then averaged the data in frequency and exported them to be imaged and self-calibrated in DIFMAP. The a priori visibility amplitude calibration was not sufficient to reliably image the weakest radio sources. We improved that by first imaging in DIFMAP the calibrator sources J0007+5706, J0147+5840, and J0630+1738, with appropriate amplitude self-calibration. We deduced correction factors for each antenna, these being consistent for the three radio sources within 2%. We corrected the amplitude calibration back in AIPS and exported the data again into DIFMAP, where the final imaging was performed after editing and averaging of the visibilities in 32 s blocks.
The EVN and MERLIN arrays have one common baseline, between JB and CM, which allows to combine both data sets and map them together. We processed the data within AIPS in order to combine both arrays. The B1950.0 (u,v) coordinates of the MERLIN data had to be corrected to the ones of the EVN for the same reference system (J2000.0) with UVFIX. Then, the MERLIN data were self-calibrated with the EVN images (see below), and the phase solutions were limited to the longest MERLIN baselines. The MERLIN data were imaged and the peak-of-brightness of both data sets were checked to be similar. As a next step, the EVN data were averaged in frequency to correspond to the MERLIN data. The AIPS headers of both data sets were modified conveniently to match together, and the weighting of both data sets was also modified to be equal with WTMOD. Finally, both data sets were concatenated (using DBCON) and exported to DIFMAP to be imaged with different data weighting in (u,v) distance (tapering) after time averaging in 32 s bins.
We interleaved cross-scans (in azimuth and elevation) with the 100 m
Effelsberg antenna to measure the radio source flux densities (A. Kraus,
private communication). We fitted a Gaussian function to the flux-density
response for every cross-scan, and we averaged the different Gaussians. We
linked the flux density scale by observing primary calibrators such as
3C 286, 3C 48, or NGC 7027 (see e.g. Kraus
1997; Peng et al. 2000). We list the single dish flux
density measurements in the second column of Table 2. The flux
density values for the main VLBI calibrators were of
Jy for 3C 286,
and
Jy for DA 193.
We present all the imaging results in Figs. 1-5,
and the image parameters in Table 2. The total flux density
values for the different images diverge from each other and from the single
dish measurements, due to the amplitude self-calibration process in all cases.
Therefore, those values should be considered with care. The minimum contours
in the images are those listed as
in Table 2,
while consecutive higher contours scale with 31/2. Here follows a
detailed discussion on each source.
![]()
|
As can be seen in our images, shown in Fig. 1, this source
appears point-like at MERLIN resolution, partially resolved in the tapered
EVN+MERLIN image and clearly resolved at EVN scales. In this last case, it
shows a two-sided jet-like structure roughly in the north-south direction,
with brighter components towards the south. The trends are visible in the
closure phases, giving us confidence that the structure observed is not a
consequence of sidelobes or imaging artifacts. In the tapered EVN+MERLIN
images, the structure extends up to 20-30 mas outside of the core, and more
clearly towards the north. This discrepancy could be due to calibration
problems.
![]() |
Figure 1:
Images of 1RXS J001442.2+580201 using the different arrays and with the parameters given in Table 2. The axes units are in mas. A (u,v)-tapering with a FWHM at 10 M![]() |
Open with DEXTER |
Model fitting of the EVN visibilities with circular Gaussians provides a
parametrization of the inner structure. Five components reproduce the
visibilities. The central one has 7.8 mJy, with a FWHM of 0.4 mas. Towards the
north, one component (N2) of 0.6 mJy at 3.6 mas (PA
,
FWHM
0.7 mas) and another one (N1) of 0.5 mJy at 8.6 mas (PA
,
extended
over 3 mas) are needed. Brighter components are present southwards, one (S2)
of 1.1 mJy at 5.0 mas (PA
,
FWHM of 0.8 mas) and the other one
(S1) of 0.4 mJy at 13.7 mas (PA
,
FWHM below 0.3 mas).
If we assume that components S1 and N1 correspond to a pair of plasma clouds
ejected at the same epoch near the compact object and perpendicularly to the
accretion disk, we can estimate some parameters of the jets by using the
following equation:
Comp. | Distance | PA |
![]() |
![]() |
[mas] | [![]() |
[![]() |
||
N1 | ![]() |
2 |
![]() |
![]() |
N2 | ![]() |
-10 |
![]() |
![]() |
S2 | ![]() |
177 |
![]() |
![]() |
S1 |
![]() |
177 |
![]() |
![]() |
![]() |
Figure 2:
Images of 1RXS J013106.4+612035 using the different arrays and with the parameters given in Table 2. The axes units are in mas. A (u,v)-tapering with a FWHM at 15 M![]() |
Open with DEXTER |
A similar approach to obtain the jet parameters of the source can be performed
thanks to the brightness asymmetry of the components using the following
equation (Mirabel & Rodríguez 1999):
In order to use this approach we will consider k=3, because the components
seem to be discrete condensations, and
,
according to the
overall spectral index reported in Paper I, since we do not have spectral
index information of the components. The use of the flux densities obtained
after model fitting gives
for the S2-N2 pair
and
for the S1-N1 pair. This last value is
certainly surprising, although it can be explained by the fact that the S1
flux density obtained after model fitting does not account for the total flux
of this plasma cloud. In fact, better estimates of the flux densities can be
obtained summing together the flux densities of the CLEAN components
obtained within each one of the four plasma clouds. This yields to
for the S2-N2 pair and
for the S1-N1 pair, in good agreement with the
values computed using the distances from the components to the core, shown in
Table 3.
If we compare the VLA position reported in Paper I with the MERLIN position in
Table 1 we can see that they are different. In fact,
taking into account the errors, the MERLIN-VLA position offsets can be
expressed as:
mas and
mas. Assuming that the difference in position is due to
intrinsic proper motions and considering the time span between both
observations, 224 days, we obtain:
mas yr-1 and
mas yr-1. Although a proper motion of the source
would clearly indicate its microquasar nature, because no proper motion would
be detected in an extragalactic source, we must be cautious with this result.
First of all, the phase-reference sources where different in the VLA and
MERLIN observations, as well as the observing frequencies from which positions
were estimated. On the other hand, none of these results exceeds the
value, preventing to state that a proper motion has been detected.
Summarizing, our results indicate that this source exhibits relativistic radio
jets with
and, therefore, together with the results reported in
Paper I, we consider 1RXS J001442.2+580201 as a very promising microquasar candidate.
We show in Fig. 2 the images obtained after our observations.
Although the source appears compact at MERLIN and EVN+MERLIN scales, there is
a weak one-sided radio jet towards the northwest in the EVN image. Model
fitting with circular Gaussian components can reproduce the observed
visibilities as follows: a central 15.4 mJy component with 0.64 mas FWHM, and
a northwest 2.1 mJy component with a FWHM of 0.74 mas, located at 1.8 mas in
PA
.
As can be seen, this last component is located at a distance
of
2 times the beam size from the core.
Using the fact that we do not detect a counter-jet, we can use
Eq. (2) replacing
with the
level value. This,
of course, will only provide a lower limit to
,
expressed as
follows:
![]() |
Figure 3:
Images of 1RXS J062148.1+174736 using the different arrays and with the parameters given in Table 2. The axes units are in mas. A (u,v)-tapering with a FWHM at 8 M![]() |
Open with DEXTER |
As done for the previous source, we can compare the VLA position with the
MERLIN one, and find that they differ in
mas
and
mas. If the offsets are real, this would imply
mas yr-1 and
mas yr-1. Hence, it seems possible that we have
detected a proper motion in right ascension at a
level. However, we
must be cautious since, as in the previous case, the phase-reference sources
and observing frequencies were different in each observation.
Overall, these results are indicative of relativistic radio jets, and hence, together with the results reported in Paper I, allow us to classify this source as a promising microquasar candidate.
This radio source was marginally seen in the cross-scans carried out in
Effelsberg. In fact, this is compatible with the low flux density of
mJy at 1.4 GHz listed in the NVSS. The source was not detected
with MERLIN or with the EVN, mainly due to problems with the phase-reference
sources, as pointed out in Sect. 2.2. In fact, as discussed in
Paper I, VLA A configuration observations showed a flux density of 0.4 mJy at
5 GHz and an inverted spectrum with a spectral index up to
,
suggesting thermal radio-emission resolved at higher angular resolutions.
However, we cannot exclude the possibility of having a highly variable source.
As can be seen in our images, shown in Fig. 3, the radio source is compact on all scales. In fact, model fitting of the EVN visibilities converges to a point-like radio source (the FWHM of a circular Gaussian tends to zero). We must note that when this source was observed it was below the horizon in SH. Therefore, the obtained beam size for the EVN image is larger than the ones obtained for the other sources, as can be seen in Table 2, hence providing lower angular resolution than in the other cases. A comparison between the VLA and MERLIN positions reveals that they are perfectly compatible within the errors, suggesting an extragalactic nature or a small proper motion if it turns out to be a galactic microquasar. Although the compactness of the source is not indicative of a galactic or an extragalactic origin, the extended optical counterpart, reported in Paper I, suggests an extragalactic origin for this source. In fact, the detected radio variability from 8 to 11 mJy within a week reported in Paper I, is compatible with the compactness in an extragalactic object (see IDV phenomenon, Wagner & Witzel 1995).
The images obtained at different resolutions are plotted in
Fig. 4. The MERLIN image presents a compact structure with some
elongation eastwards, while the EVN and combined EVN+MERLIN images show a
clear one-sided jet towards the east, with a slight bent towards the south at
larger core separations. The closure phases clearly show that the source
departs from symmetry, with preferred emission to the east, both in the MERLIN
and the EVN data sets. Two distinct components are present in the EVN image,
at 9 and 17 mas from the compact core (PA of 89 and 113
,
respectively). Those components can be model fitted with elliptical Gaussians,
yielding flux densities of 3.7 and 2.8 mJy, respectively, for a core of
40.3 mJy.
![]() |
Figure 4:
Images of 1RXS J072259.5-073131 using the different arrays and with the parameters given in Table 2. The axes units are in mas. A (u,v)-tapering with a FWHM at 15 M![]() |
Open with DEXTER |
![]() |
Figure 5:
Images of the already identified quasar 1RXS J072418.3-071508 using the different arrays and with the parameters given in Table 2. The axes units are in mas. A (u,v)-tapering with a FWHM at 10 M![]() |
Open with DEXTER |
Using again the fact that we do not detect a counter-jet, we can use
Eq. (3) with
mJy (the closest component to the core
in the EVN image),
mJy (as previously done, the
value
has been taken as the root mean square noise in the image),
(see Paper I), and k=3 to be consistent with the
discrete nature of the components, to obtain
,
and hence
and
.
Therefore, these
results point towards relativistic radio jets as the origin of the elongated
radio emission present in the images.
As pointed out for 1RXS J013106.4+612035, the one-sided jet morphology does not rule out a microquasar nature for this object. However, the bending of the jet at such small angular scales resembles the ones seen in blazars. In fact, as reported in Paper I, a one-sided arcsecond scale jet is also present in VLA A configuration observations at 1.4 GHz, an unusual feature in the already known microquasars.
A comparison between the VLA and MERLIN positions reveals that they agree within the errors, suggesting an extragalactic nature or a small proper motion if it turns out to be a galactic microquasar.
Overall, these results are indicative of relativistic radio jets, although this source shows characteristics more similar to blazars than to microquasars.
This radio source has recently been (March 2002) classified as a quasar in the
SIMBAD database, and it is not any more a microquasar candidate. It is
catalogued in the Parkes-MIT-NRAO (Griffith et al. 1994) and the
Texas Survey (Douglas et al. 1996). It is listed as
PMN J0724-0715 in the NED database, and is the source
WGA J0724.3-0715 in Perlman et al. (1998), who reported
a faint and quite broad H
emission line (rest-frame
Å, FWHM=4000 km s-1), and classified it as a Flat
Spectrum Radio Quasar (FSRQ) with z=0.270. Nevertheless, we have reported
here our observational results for this source, since it was a candidate when
we performed the observations. It presents (Fig. 5) a one-sided
pc-scale jet oriented towards the northeast, changing from a PA of
at 5 mas from the core (EVN image) to
up to 200 mas,
at MERLIN scales.
Model fitting of the EVN visibilities with circular Gaussians reveals a
compact core (0.7 mas FWHM) with 245.4 mJy, and two distinct components, one
with 11.9 mJy at 6.2 mas (PA
,
1.9 mas FWHM, present in the right
panel image of Fig. 5) and another one with the size of the
beam, a flux density of 1.7 mJy at a distance of 29.6 mas in PA
(visible in the middle panel image of Fig. 5).
![]()
|
Using again the fact that we do not detect a counter-jet, we can use
Eq. (3) with
mJy (the closest component to the
core),
mJy,
(see Paper I), and k=3 to be
consistent with the discrete nature of the components, to obtain
.
Hence, an upper limit of
and a lower limit of
is obtained, pointing towards
relativistic radio jets as the origin of the elongated radio emission present
in the images of this already identified quasar.
We have presented EVN and MERLIN observations of the six sources studied by
Paredes et al. (2002) in their search for microquasar candidates at low
galactic latitudes. The first one, namely 1RXS J001442.2+580201, displays a two-sided radio
jet, which after analysis implies
and
.
1RXS J013106.4+612035, displays a one-sided radio jet, requiring
and
.
The third one, namely 1RXS J042201.0+485610,
was not detected due to its low flux density and/or to phase-referencing
problems. 1RXS J062148.1+174736 appeared compact at all scales. The fifth one, namely
1RXS J072259.5-073131, displays a bent one-sided radio jet, implying
and
.
Finally, 1RXS J072418.3-071508 shows also a bent one-sided jet,
requiring
and
.
After a detailed analysis of our data, we show in Table 4 a summary of the results obtained after the VLA and optical observations (Paper I), and the EVN+MERLIN observations reported here. As can be seen, the first two sources, 1RXS J001442.2+580201 and 1RXS J013106.4+612035, are promising microquasar candidates. 1RXS J042201.0+485610 is probably of thermal nature due to the highly inverted spectrum at high radio frequencies, while 1RXS J062148.1+174736 is probably an extragalactic object due to the extended nature of the optical counterpart. 1RXS J072259.5-073131 shows properties common to blazars, while 1RXS J072418.3-071508 is an already identified quasar. We note that 1RXS J072259.5-073131 is bright enough at radio wavelengths to attempt an Hi absorption experiment, that could allow to determine if this source is galactic or not. In any case, optical spectroscopic observations of the first five sources are in progress, to clearly unveil their galactic or extragalactic nature.
Acknowledgements
We acknowledge R. Porcas and W. Alef for their useful comments and suggestions after reading through a draft version of this paper. We acknowledge useful comments from L. F. Rodríguez, the referee of this paper. We are very grateful to S. T. Garrington, M. A. Garrett, D. C. Gabuzda, and C. Reynolds for their valuable help in the data reduction process. This paper is based on observations with the 100-m telescope of the MPIfR (Max-Planck-Institut für Radioastronomie) at Effelsberg. We thank the staff of the JIVE correlator and of the observing telescopes, especially A. Kraus for the single dish flux density measurements at the 100 m antenna in Effelsberg. The European VLBI Network is a joint facility of European, Chinese and other radio astronomy institutes funded by their national research councils. MERLIN is operated as a National Facility by the University of Manchester at Jodrell Bank Observatory on behalf of the UK Particle Physics & Astronomy Research Council. The EVN observations were carried out thanks to the TMR Access to Large-scale Facilities programme under contract No. ERBFMGECT950012. Part of the data reduction was done at JIVE with the support of the European Community - Access to Research Infrastructure action of the Improving Human Potential Programme under contract No. HPRI-CT-1999-00045. M.R., J.M.P. and J.M. acknowledge partial support by DGI of the Ministerio de Ciencia y Tecnología (Spain) under grant AYA2001-3092, as well as partial support by the European Regional Development Fund (ERDF/FEDER). During this work, M.R. has been supported by two fellowships from CIRIT (Generalitat de Catalunya, ref. 1998 BEAI 200293 and 1999 FI 00199). J.M. has been aided in this work by an Henri Chrétien International Research Grant administered by the American Astronomical Society, and has been partially supported by the Junta de Andalucía. This research has made use of the NASA's Astrophysics Data System Abstract Service, of the SIMBAD database, operated at CDS, Strasbourg, France, and of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.