(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.

MGRO J2019+37 is an unidentified source of very high energy gamma-rays originally reported by the MILAGRO collaboration as the brightest TeV source in the Cygnus region. Despite the poor angular resolution of MILAGRO, this object seems to be most likely an extended source or, alternatively, a superposition of point-like TeV sources. In order to contribute to the understanding of this peculiar object, we have mosaiced it with the Giant Metrewave Radio Telescope (GMRT) in Pune, India, at the 610 MHz frequency covering a field of view of about 6 square degrees down to a typical rms noise of a few tenths of mJy. We also observed the central square degree of this mosaic in the near infrared Ks-band using the 3.5 m telescope and the OMEGA2000 camera at the Calar Alto observatory (Spain). We present here a first account of our observations and results, together with a preliminar cross correlation of the radio and infrared source catalog as well as a correlation with the available X-ray observations of the region.

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).

In this paper we present a multiresolution-based method for period determination that is able to deal with unevenly sampled data. This method allows us to detect superimposed periodic signals with lower signal-to-noise ratios than in classical methods. This multiresolution-based method is a generalization of the wavelet-based method for period detection that is unable to deal with unevenly sampled data, presented by the authors in Otazu et al. (2002). This new method is a useful tool for the analysis of real data and, in particular, for the analysis of astronomical data.

Here we discuss two consecutive MERLIN observations of the X-ray binary LS I +61 303. The first observation shows a double-sided jet extending up to about 200 AU on both sides of a central source. The jet shows a bent S-shaped structure similar to the one displayed by the well-known precessing jet of SS 433. The precession suggested in the first MERLIN image becomes evident in the second one, showing a one-sided bent jet significantly rotated with respect to the jet of the day before. We conclude that the derived precession of the relativistic (beta=0.6) jet explains puzzling previous VLBI results. Moreover, the fact that the precession is fast could be the explanation of the never understood short term (days) variability of the associated gamma-ray source 2CG 135+01 / 3EG J0241+6103.

In this paper we present a wavelet based algorithm that is able to detect superimposed periodic signals in data with low signal-noise ratio. In this context, the results given by classical period determination methods highly depend on the intrinsic characteristics of each periodic signal, like amplitude or profile. It is then difficult to detect the different periods present in the data set. The results given by the wavelet based method for period determination we present here are independent of the characteristics of the signals.

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.

We report millimetre-wave continuum observations of the X-ray binaries Cygnus X-3, SS 433, LSI+61 303, Cygnus X-1 and GRS 1915+105. The observations were carried out with the IRAM 30 m-antenna at 250 GHz (1.25 mm) from 1998 March 14 to March 20. These millimetre measurements are complemented with centimetre observations from the Ryle Telescope, at 15 GHz (2.0 cm) and from the Green Bank Interferometer at 2.25 and 8.3 GHz (13 and 3.6 cm). Both Cygnus X-3 and SS 433 underwent moderate flaring events during our observations, whose main spectral evolution properties are described and interpreted. A significant spectral steepening was observed in both sources during the flare decay, that is likely to be caused by adiabatic expansion, inverse Compton and synchrotron losses. Finally, we also report 250 GHz upper limits for three additional undetected X-ray binary stars: LSI+65 010, LSI+61 235 and X Per.