The gamma-ray emitting microquasar LSI+61303

LS I +61 303 is one of the most studied X-ray binary systems because of its two peculiarities: On the one hand being the probable counterpart of the variable gamma ray source 2CG 135+01 (Gregory and Taylor 1978; Tavani et al. 1998) and on the other hand being a periodic radio source (Taylor and Gregory 1982). The recent discovery of a radio emitting jet extending ca. 200 AU at both sides of a central core (Massi et al. 2004) in all evidence has shown the occurrence of accretion/ejection processes in this system. However, the radio outbursts do not occur at periastron passage, where the accretion is at its maximum, but several days later. In addition, when the gamma-ray emission of 2CG 135+01 is examined along the orbital phase of LS I +61 303 one sees that this emission seems to peak at periastron passage (Massi 2004). Here in detail we analyse the trend of the gamma-ray data versus orbital phase and discuss the delay between peaks at gamma-rays and in the radio band within the framework of a two-peak accretion/ejection model proposed by Taylor et al. (1992) and further developed by Marti' and Paredes (1995).Comment: To be published in the proceedings of the Symposium on High-Energy Gamma-Ray Astronomy, Heidelberg, July 26-30, 2004 (AIP Proceedings Series

One-sided jet at milliarcsecond scales in LSI+61303

We present Very Long Baseline Interferometry (VLBI) observations of the high mass X-ray binary LSI+61303, carried out with the European VLBI Network (EVN). Over the 11 hour observing run, performed 10 days after a radio outburst, the radio source showed a constant flux density, which allowed sensitive imaging of the emission distribution. The structure in the map shows a clear extension to the southeast. Comparing our data with previous VLBI observations we interpret the extension as a collimated radio jet as found in several other X-ray binaries. Assuming that the structure is the result of an expansion that started at the onset of the outburst, we derive an apparent expansion velocity of 0.003 c, which, in the context of Doppler boosting, corresponds to an intrinsic velocity of at least 0.4 c for an ejection close to the line of sight. From the apparent velocity in all available epochs we are able to establish variations in the ejection angle which imply a precessing accretion disk. Finally we point out that LSI+61303, like SS433 and Cygnus X-1, shows evidence for an emission region almost orthogonal to the relativistic jet.Comment: 7 pages, 4 figures, LaTeX, uses aa.cls. Accepted for publication in A&

Chandra Observations of the Gamma-ray Binary LSI+61303: Extended X-ray Structure?

We present a 50 ks observation of the gamma-ray binary LSI+61303 carried out with the ACIS-I array aboard the Chandra X-ray Observatory. This is the highest resolution X-ray observation of the source conducted so far. Possible evidence of an extended structure at a distance between 5 and 12 arcsec towards the North of LSI+61303 have been found at a significance level of 3.2 sigma. The asymmetry of the extended emission excludes an interpretation in the context of a dust-scattered halo, suggesting an intrinsic nature. On the other hand, while the obtained source flux, of F_{0.3-10 keV}=7.1^{+1.8}_{-1.4} x 10^{-12} ergs/cm^2/s, and hydrogen column density, N_{H}=0.70+/-0.06 x 10^{22} cm^{-2}, are compatible with previous results, the photon index Gamma=1.25+/-0.09 is the hardest ever found. In light of these new results, we briefly discuss the physics behind the X-ray emission, the location of the emitter, and the possible origin of the extended emission ~0.1 pc away from LSI+61303.Comment: 4 pages, 3 figures. Accepted for publication in ApJ Letter

VLBA images of the precessing jet of LSI+61303

Context: In 2004, changes in the radio morphology of the Be/X-ray binary system LSI+61303 suggested that it is a precessing microquasar. In 2006, a set of VLBA observations performed throughout the entire orbit of the system were not used to study its precession because the changes in radio morphology could tentatively be explained by the alternative pulsar model. However, a recent radio spectral index data analysis has confirmed the predictions of the two-peak microquasar model, which therefore does apply in LSI+61303. Aims: We revisit the set of VLBA observations performed throughout the orbit to determine the precession period and improve our understanding of the physical mechanism behind the precession. Methods: By reanalyzing the VLBA data set, we improve the dynamic range of images by a factor of four, using self-calibration. Different fitting techniques are used and compared to determine the peak positions in phase-referenced maps. Results: The improved dynamic range shows that in addition to the images with a one-sided structure, there are several images with a double-sided structure. The astrometry indicates that the peak in consecutive images for the whole set of observations describes a well-defined ellipse, 6-7 times larger than the orbit, with a period of about 28 d. Conclusions: A double-sided structure is not expected to be formed from the expanding shocked wind predicted in the pulsar scenario. In contrast, a precessing microquasar model can explain the double- and one-sided structures in terms of variable Doppler boosting. The ellipse defined by the astrometry could be the cross-section of the precession cone, at the distance of the 8.4 GHz-core of the steady jet, and 28d the precession period.Comment: 7 pages, 5 figures, Accepted for publication in Astronomy and Astrophysics, added references for sect.

Radio Spectral Index Analysis and Classes of Ejection in LS I +61 303

LS I +61303 is a gamma-ray binary with periodic radio outbursts coincident with the orbital period of P=26.5 d. The origin of the radio emission is unclear,it could be due either to a jet, as in microquasars, or to the shock boundary between the Be star and a possible pulsar wind. We here analyze the radio spectral index over 6.7 yr from Green Bank Interferometer data at 2.2 GHz and 8.3 GHz. We find two new characteristics in the radio emission. The first characteristic is that the periodic outbursts indeed consist of two consecutive outbursts; the first outburst is optically thick, whereas the second outburst is optically thin. The spectrum of LS I +61 303 is well reproduced by the shock-in-jet model commonly used in the context of microquasars and AGNs: the optically thin spectrum is due to shocks caused by relativistic plasma ("transient jet") traveling through a pre-existing much slower steady flow ("steady jet"). This steady flow is responsible for the preceding optically thick spectrum. The second characteristic we find is that the observed spectral evolution, from optically thick to optically thin emission, occurs twice during the orbital period. We observed this occurrence at the orbital phase of the main 26.5 d outburst and also at an earlier phase, shifted by $\Delta \Phi \sim$ 0.3 (i.e almost 8 days before). We show that this result qualitatively and quantitatively agrees with the two-peak accretion/ejection model proposed in the past for LS I +61303. We conclude that the radio emission in LS I +61303 originates from a jet and suggest that the variable TeV emission comes from the usual Compton losses expected as an important by-product in the shock-in-jet theory.Comment: 27 pages, 7 figures, accepted for publication in Ap

A One-sided, Highly Relativistic Jet from Cygnus X-3

Very Long Baseline Array images of the X-ray binary, Cygnus X-3, were obtained 2, 4 and 7 days after the peak of a 10 Jy flare on 4 February 1997. The first two images show a curved one-sided jet, the third a scatter-broadened disc, presumably at the position of the core. The jet curvature changes from the first to the second epoch, which strongly suggests a precessing jet. The ratio of the flux density in the approaching to that in the (undetected) receding jet is > 330; if this asymmetry is due to Doppler boosting, the implied jet speed is > 0.81c. Precessing jet model fits, together with the assumptions that the jet is intrinsically symmetric and was ejected during or after the major flare, yield the following constraints: the jet inclination to the line of sight must be < 14 degrees; the cone opening angle must be < 12 degrees; and the precession period must be > 60 days.Comment: 12 pages 7 figures, accepted by Ap

Feasibility Study of Lense-Thirring Precession in LS I +61303

Very recent analysis of the radio spectral index and high energy observations have shown that the two-peak accretion/ejection microquasar model applies for LSI+61303. The fast variations of the position angle observed with MERLIN and confirmed by consecutive VLBA images must therefore be explained in the context of the microquasar scenario. We calculate what could be the precessional period for the accretion disk in LSI+61303 under tidal forces of the Be star (P_{tidal-forces}) or under the effect of frame dragging produced by the rotation of the compact object (P_{Lense-Thirring}). P_{tidal-forces}$is more than one year. P_{Lense-Thirring} depends on the truncated radius of the accretion disk,$R_{tr}$. We determined R_{tr}=300 r_g for observed QPO at 2 Hz. This value is much above the few$r_g$, where the Bardeen-Petterson effect should align the midplane of the disk. For this truncated radius of the accretion disk P_{Lense-Thirring} for a slow rotator results in a few days. Therefore, Lense-Thirring precession induced by a slowly rotating compact object could be compatible with the daily variations of the ejecta angle observed in LSI+61303.Comment: 6 pages, 5 figures, accepted for publication in Astronomy and Astrophysic

Radio continuum and near-infrared study of the MGRO J2019+37 region

(abridged) MGRO J2019+37 is an unidentified extended source of VHE gamma-rays originally reported by the Milagro Collaboration as the brightest TeV source in the Cygnus region. Its extended emission could be powered by either a single or several sources. The GeV pulsar AGL J2020.5+3653, discovered by AGILE and associated with PSR J2021+3651, could contribute to the emission from MGRO J2019+37, although extrapolation of the GeV spectrum does not explain the detected multi-TeV flux. Our aim is to identify radio and NIR sources in the field of the extended TeV source MGRO J2019+37, and study potential counterparts that could contribute to its emission. We surveyed a region of about 6 square degrees with the Giant Metrewave Radio Telescope (GMRT) at the frequency 610 MHz. We also observed the central square degree of this survey in the NIR Ks-band using the 3.5 m telescope in Calar Alto. Archival X-ray observations of some specific fields are included. VLBI observations of an interesting radio source were performed. We explored possible scenarios to produce the multi-TeV emission from MGRO J2019+37 and studied which of the sources could be the main particle accelerator. We present a catalogue of 362 radio sources detected with the GMRT in the field of MGRO J2019+37, and the results of a cross-correlation of this catalog with one obtained at NIR wavelengths, as well as with available X-ray observations of the region. Some peculiar sources inside the ~1 degree uncertainty region of the TeV emission from MGRO J2019+37 are discussed in detail, including the pulsar PSR J2021+3651 and its pulsar wind nebula PWN G75.2+0.1, two new radio-jet sources, the HII region Sh 2-104 containing two star clusters, and the radio source NVSS J202032+363158.Comment: 10 pages, 6 figures, 2 tables, accepted for publication in Astronomy and Astrophysic

Radio continuum and near-infrared study of the MGRO J2019+37 region

(abridged) MGRO J2019+37 is an unidentified extended source of VHE gamma-rays originally reported by the Milagro Collaboration as the brightest TeV source in the Cygnus region. Its extended emission could be powered by either a single or several sources. The GeV pulsar AGL J2020.5+3653, discovered by AGILE and associated with PSR J2021+3651, could contribute to the emission from MGRO J2019+37, although extrapolation of the GeV spectrum does not explain the detected multi-TeV flux. Our aim is to identify radio and NIR sources in the field of the extended TeV source MGRO J2019+37, and study potential counterparts that could contribute to its emission. We surveyed a region of about 6 square degrees with the Giant Metrewave Radio Telescope (GMRT) at the frequency 610 MHz. We also observed the central square degree of this survey in the NIR Ks-band using the 3.5 m telescope in Calar Alto. Archival X-ray observations of some specific fields are included. VLBI observations of an interesting radio source were performed. We explored possible scenarios to produce the multi-TeV emission from MGRO J2019+37 and studied which of the sources could be the main particle accelerator. We present a catalogue of 362 radio sources detected with the GMRT in the field of MGRO J2019+37, and the results of a cross-correlation of this catalog with one obtained at NIR wavelengths, as well as with available X-ray observations of the region. Some peculiar sources inside the ~1 degree uncertainty region of the TeV emission from MGRO J2019+37 are discussed in detail, including the pulsar PSR J2021+3651 and its pulsar wind nebula PWN G75.2+0.1, two new radio-jet sources, the HII region Sh 2-104 containing two star clusters, and the radio source NVSS J202032+363158.Comment: 10 pages, 6 figures, 2 tables, accepted for publication in Astronomy and Astrophysic