We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

We simulate the spectral behaviour of the non-transient gamma-ray binaries using archival observations as a reference. With this we test the CTA capability to measure the sources' spectral parameters and detect variability on various time scales. We review the known properties of gamma-ray binaries and the theoretical models that have been used to describe their spectral and timing characteristics. We show that CTA is capable of studying these sources on time scales comparable to their characteristic variability time scales. For most of the binaries, the unprecedented sensitivity of CTA will allow the spectral evolution to be studied on a time scale as short as 30 min. This will enable a direct comparison of the TeV and lower energy (radio to GeV) properties of these sources from simultaneous observations. We also review the source-specific questions that can be addressed with such high-accuracy CTA measurements.

We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-flaring source at the time such neutrinos are recorded. The use of neutrino-triggered alerts thus aims at increasing the availability of simultaneous multi-messenger data during potential neutrino flaring activity, which can increase the discovery potential and constrain the phenomenological interpretation of the high-energy emission of selected source classes (e.g. blazars). The requirements of a fast and stable online analysis of potential neutrino signals and its operation are presented, along with first results of the program operating between 14 March 2012 and 31 December 2015.

The BL Lac object 1ES 1011+496 was discovered at Very High Energy gamma-rays by MAGIC in spring 2007. Before that the source was little studied in different wavelengths. Therefore a multi-wavelength (MWL) campaign was organized in spring 2008. Along MAGIC, the MWL campaign included the Metsahovi radio observatory, Bell and KVA optical telescopes and the Swift and AGILE satellites. MAGIC observations span from March to May, 2008 for a total of 27.9 hours, of which 19.4 hours remained after quality cuts. The light curve showed no significant variability. The differential VHE spectrum could be described with a power-law function. Both results were similar to those obtained during the discovery. Swift XRT observations revealed an X-ray flare, characterized by a harder when brighter trend, as is typical for high synchrotron peak BL Lac objects (HBL). Strong optical variability was found during the campaign, but no conclusion on the connection between the optical and VHE gamma-ray bands could be drawn. The contemporaneous SED shows a synchrotron dominated source, unlike concluded in previous work based on nonsimultaneous data, and is well described by a standard one zone synchrotron self Compton model. We also performed a study on the source classification. While the optical and X-ray data taken during our campaign show typical characteristics of an HBL, we suggest, based on archival data, that 1ES 1011+496 is actually a borderline case between intermediate and high synchrotron peak frequency BL Lac objects.

The Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula have been claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 hours of good-quality data, and searched for emission from the pulsar and pulsar wind nebula. We did not find any significant detection, and derived 95% confidence level upper limits. The resulting upper limits of 5.3 x 10^{-13} TeV cm^{-2} s^{-1} for the Geminga pulsar and 3.5 x 10^{-12} TeV cm^{-2} s^{-1} for the surrounding nebula at 50 GeV are the most constraining ones obtained so far at VHE. To complement the VHE observations, we also analyzed 5 years of Fermi-LAT data from Geminga, finding that the sub-exponential cut-off is preferred over the exponential cut-off that has been typically used in the literature. We also find that, above 10 GeV, the gamma-ray spectra from Geminga can be described with a power law with index softer than 5. The extrapolation of the power-law Fermi-LAT pulsed spectra to VHE goes well below the MAGIC upper limits, indicating that the detection of pulsed emission from Geminga with the current generation of Cherenkov telescopes is very difficult.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

The Cherenkov Telescope Array (CTA), the new generation very high-energy gamma-ray observatory, will improve the flux sensitivity of the current Cherenkov telescopes by an order of magnitude over a continuous range from about 10 GeV to above 100 TeV. With tens of telescopes distributed in the Northern and Southern hemispheres, the large effective area and field of view coupled with the fast pointing capability make CTA a crucial instrument for the detection and understanding of the physics of transient, short-timescale variability phenomena (e.g. Gamma-Ray Bursts, Active Galactic Nuclei, gamma-ray binaries, serendipitous sources). The key CTA system for the fast identification of flaring events is the Real-Time Analysis (RTA) pipeline, a science alert system that will automatically detect and generate science alerts with a maximum latency of 30 seconds with respect to the triggering event collection and ensure fast communication to/from the astrophysics community. According to the CTA design requirements, the RTA search for a true transient event should be performed on multiple time scales (from minutes to hours) with a sensitivity not worse than three times the nominal CTA sensitivity. Given the CTA requirement constraints on the RTA efficiency and the fast response ability demanded by the transient science, we perform a preliminary evaluation of the RTA sensitivity as a function of the CTA high-level technical performance (e.g. effective area, point spread function) and the observing time. This preliminary approach allows the exploration of the complex parameter space defined by the scientific and technological requirements, with the aim of defining the feasibility range of the input parameters and the minimum background rejection capability of the RTA pipeline.

The Cherenkov Telescope Array (CTA) observatory will be one of the largest ground-based very high-energy gamma-ray observatories. The On-Site Analysis will be the first CTA scientific analysis of data acquired from the array of telescopes, in both northern and southern sites. The On-Site Analysis will have two pipelines: the Level-A pipeline (also known as Real-Time Analysis, RTA) and the level-B one. The RTA performs data quality monitoring and must be able to issue automated alerts on variable and transient astrophysical sources within 30 seconds from the last acquired Cherenkov event that contributes to the alert, with a sensitivity not worse than the one achieved by the final pipeline by more than a factor of 3. The Level-B Analysis has a better sensitivity (not be worse than the final one by a factor of 2) and the results should be available within 10 hours from the acquisition of the data: for this reason this analysis could be performed at the end of an observation or next morning. The latency (in particular for the RTA) and the sensitivity requirements are challenging because of the large data rate, a few GByte/s. The remote connection to the CTA candidate site with a rather limited network bandwidth makes the issue of the exported data size extremely critical and prevents any kind of processing in real-time of the data outside the site of the telescopes. For these reasons the analysis will be performed on-site with infrastructures co-located with the telescopes, with limited electrical power availability and with a reduced possibility of human intervention. This means, for example, that the on-site hardware infrastructure should have low-power consumption. A substantial effort towards the optimization of high-throughput computing service is envisioned to provide hardware and software solutions with high-throughput, low-power consumption at a low-cost.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

We find that the formation of MWC 656 (the first Be binary containing a black hole) involves a common envelope phase and a supernova explosion. This result supports the idea that a rapidly rotating Be star can emerge out of a common envelope phase, which is very intriguing because this evolutionary stage is thought to be too fast to lead to significant accretion and spin up of the B star. We predict $\sim 10-100$ of B BH binaries to currently reside in the Galactic disk, among which around $1/3$ contain a Be star, but there is only a small chance to observe a system with parameters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be BH binary population, it may indicate that standard evolutionary theory needs to be revised. This would pose another evolutionary problem in understanding BH binaries, with BH X-ray Novae formation issue being the prime example. The future evolution of MWC 656 with a $\sim 5$ M$_{\odot}$ black hole and with a $\sim 13$ M$_{\odot}$ main sequence companion on a $\sim 60$ day orbit may lead to the formation of a coalescing BH-NS system. The estimated Advanced LIGO/Virgo detection rate of such systems is up to $\sim 0.2$ yr$^{-1}$. This empirical estimate is a lower limit as it is obtained with only one particular evolutionary scenario, the MWC 656 binary. This is only a third such estimate available (after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so far undetected BH--NS binaries.

The MAGIC telescopes are two Imaging Atmospheric Cherenkov Telescopes (IACTs) located on the Canary island of La Palma. The telescopes are designed to measure Cherenkov light from air showers initiated by gamma rays in the energy regime from around 50 GeV to more than 50 TeV. The two telescopes were built in 2004 and 2009, respectively, with different cameras, triggers and readout systems. In the years 2011-2012 the MAGIC collaboration undertook a major upgrade to make the stereoscopic system uniform, improving its overall performance and easing its maintenance. In particular, the camera, the receivers and the trigger of the first telescope were replaced and the readout of the two telescopes was upgraded. This paper (Part I) describes the details of the upgrade as well as the basic performance parameters of MAGIC such as raw data treatment, dead time of the system, linearity in the electronic chain and sources of noise. In Part II, we describe the physics performance of the upgraded system.

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma, Spain. During summer 2011 and 2012 it underwent a series of upgrades, involving the exchange of the MAGIC-I camera and its trigger system, as well as the upgrade of the readout system of both telescopes. We use observations of the Crab Nebula taken at low and medium zenith angles to assess the key performance parameters of the MAGIC stereo system. For low zenith angle observations, the standard trigger threshold of the MAGIC telescopes is ~50GeV. The integral sensitivity for point-like sources with Crab Nebula-like spectrum above 220GeV is (0.66+/-0.03)% of Crab Nebula flux in 50 h of observations. The angular resolution, defined as the sigma of a 2-dimensional Gaussian distribution, at those energies is < 0.07 degree, while the energy resolution is 16%. We also re-evaluate the effect of the systematic uncertainty on the data taken with the MAGIC telescopes after the upgrade. We estimate that the systematic uncertainties can be divided in the following components: < 15% in energy scale, 11-18% in flux normalization and +/-0.15 for the energy spectrum power-law slope.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

Stellar-mass black holes have all been discovered through X-ray emission, which arises from the accretion of gas from their binary companions (this gas is either stripped from low-mass stars or supplied as winds from massive ones). Binary evolution models also predict the existence of black holes accreting from the equatorial envelope of rapidly spinning Be-type stars (stars of the Be type are hot blue irregular variables showing characteristic spectral emission lines of hydrogen). Of the ~80 Be X-ray binaries known in the Galaxy, however, only pulsating neutron stars have been found as companions. A black hole was formally allowed as a solution for the companion to the Be star MWC 656 (also known as HD 215227), although that was based on a single radial velocity curve of the Be star, a mistaken spectral classification and rough estimates of the inclination angle. Here we report observations of an accretion disk line mirroring the orbit of the Be star. This, together with an improved radial velocity curve of the Be star through fitting sharp Fe II profiles from the equatorial disk, and a refined Be classification (to that of a B1.5-B2 III star), reveals a black hole of 3.8 to 6.9 solar masses orbiting MWC 656, the candidate counterpart of the gamma-ray source AGL J2241+4454. The black hole is X-ray quiescent and fed by a radiatively inefficient accretion flow giving a luminosity less than 1.6 x 10-7 times the Eddington luminosity. This implies that Be binaries with black-hole companions are difficult to detect by conventional X-ray surveys.

There are a number of very high energy sources in the Galaxy that remain unidentified. Multi-wavelength and variability studies, and catalogue searches, are powerful tools to identify the physical counterpart, given the uncertainty in the source location and extension. This work carries out a thorough multi-wavelength study of the unidentified, very high energy source HESS J1858+020 and its environs. Giant Metrewave Radio Telescope observations at 610 MHz and 1.4 GHz have been done to obtain a deep, low-frequency radio image of the region surrounding HESS J1858+020. Archival radio, infrared, and X-ray data have been analysed as well. This observational information, combined with molecular data, catalogue sources, and a nearby Fermi gamma-ray detection of unidentified origin, are combined to explore possible counterparts to the very high energy source. We provide with a deep radio image of a supernova remnant that might be related to the GeV and TeV emission in the region. We confirm the presence of an H II region next to the supernova remnant and coincident with molecular emission. A potential region of star formation is also identified. We identify several radio and X-ray sources in the surroundings. Some of these sources are known planetary nebulae, whereas others may be non-thermal extended emitters and embedded young stellar objects. Three old, background Galactic pulsars also neighbour HESS J1858+020 along the line of sight. The region surrounding HESS J1858+020 is rich in molecular structures and non-thermal objects that may potentially be linked to this unidentified very high energy source. In particular, a supernova remnant interacting with nearby molecular clouds may be a good candidate, but a star forming region, or a non-thermal radio source of yet unclear nature, may also be behind the gamma-ray source. Further observational studies are needed.

Several types of binary systems have been detected up to now at high and very high gamma-ray energies, including microquasars, young pulsars around massive stars and colliding wind binaries. The study of the sources already known, and of the new types of sources expected to be discovered with the unprecedented sensitivity of CTA, will allow us to qualitatively improve our knowledge on particle acceleration, emission and radiation reprocessing, and on the dynamics of flows and their magnetic fields. Here we present some examples of the capabilities of CTA to probe the flux and spectral changes that typically occur in these astrophysical sources, as well as to search for delays in correlated X-ray/TeV variability with CTA and satellites of the CTA era. Our results show that our knowledge of the high-energy physics in binary systems will significantly deepen with CTA.

The binary systems that have been detected in gamma rays have proven very useful to study high-energy processes, in particular particle acceleration, emission and radiation reprocessing, and the dynamics of the underlying magnetized flows. Binary systems, either detected or potential gamma-ray emitters, can be grouped in different subclasses depending on the nature of the binary components or the origin of the particle acceleration: the interaction of the winds of either a pulsar and a massive star or two massive stars; accretion onto a compact object and jet formation; and interaction of a relativistic outflow with the external medium. We evaluate the potentialities of an instrument like the Cherenkov telescope array (CTA) to study the non-thermal physics of gamma-ray binaries, which requires the observation of high-energy phenomena at different time and spatial scales. We analyze the capability of CTA, under different configurations, to probe the spectral, temporal and spatial behavior of gamma-ray binaries in the context of the known or expected physics of these sources. CTA will be able to probe with high spectral, temporal and spatial resolution the physical processes behind the gamma-ray emission in binaries, significantly increasing as well the number of known sources. This will allow the derivation of information on the particle acceleration and emission sites qualitatively better than what is currently available.

Context: The blazar Markarian 421 is one of the brightest TeV gamma-ray sources of the northern sky. From December 2007 until June 2008 it was intensively observed in the very high energy (VHE, E > 100 GeV) band by the single-dish Major Atmospheric Gamma-ray Imaging Cherenkov telescope (MAGIC-I). Aims: We aimed to measure the physical parameters of the emitting region of the blazar jet during active states. Methods: We performed a dense monitoring of the source in VHE with MAGIC-I, and also collected complementary data in soft X-rays and optical-UV bands; then, we modeled the spectral energy distributions (SED) derived from simultaneous multi-wavelength data within the synchrotron self--compton (SSC) framework. Results: The source showed intense and prolonged gamma-ray activity during the whole period, with integral fluxes (E > 200 GeV) seldom below the level of the Crab Nebula, and up to 3.6 times this value. Eight datasets of simultaneous optical-UV (KVA, Swift/UVOT), soft X-ray (Swift/XRT) and MAGIC-I VHE data were obtained during different outburst phases. The data constrain the physical parameters of the jet, once the spectral energy distributions obtained are interpreted within the framework of a single-zone SSC leptonic model. Conclusions: The main outcome of the study is that within the homogeneous model high Doppler factors (40 <= delta <= 80) are needed to reproduce the observed SED; but this model cannot explain the observed short time-scale variability, while it can be argued that inhomogeneous models could allow for less extreme Doppler factors, more intense magnetic fields and shorter electron cooling times compatible with hour or sub-hour scale variability.

We argue that the bright flare of the binary pulsar \object{PSR B1259$-$63/LS2883} detected by the {\it Fermi} Large Area Telescope (LAT), is due to the inverse Compton (IC) scattering of the unshocked electron-positron pulsar wind with a Lorentz factor $\Gamma_0 \approx 10^4$. The combination of two effects both linked to the circumstellar disk (CD), is a key element in the proposed model. The first effect is related to the impact of the surrounding medium on the termination of the pulsar wind. Inside the disk, the "early" termination of the wind results in suppression of its gamma-ray luminosity. When the pulsar escapes the disk, the conditions for termination of the wind undergo significant changes. This would lead to a dramatic increase of the pulsar wind zone, and thus to the proportional increase of the gamma-ray flux. On the other hand, if the parts of the CD disturbed by the pulsar can supply infrared photons of density high enough for efficient Comptonization of the wind, almost the entire kinetic energy of the pulsar wind would be converted to radiation, thus the gamma-ray luminosity of the wind could approach to the level of the pulsar's spin-down luminosity as reported by the {\it Fermi} collaboration.

We present optical spectroscopy of MWC 656 and MWC 148, the proposed optical counterparts of the gamma-ray sources AGL J2241+4454 and HESS J0632+0 57, respectively. The main parameters of the Halpha emission line (EW, FWHM and centroid velocity) in these stars are modulated on the proposed orbital periods of 60.37 and 321 days, respectively. These modulations are likely produced by the resonant interaction of the Be discs with compact stars in eccentric orbits. We also present radial velocity curves of the optical stars folded on the above periods and obtain the first orbital elements of the two gamma-ray sources thus confirming their binary nature. Our orbital solution support eccentricities e~0.4 and 0.83+-0.08 for MWC 656 and MWC 148, respectively. Further, our orbital elements imply that the X-ray outbursts in HESS J0632+057/MWC 148 are delayed ~0.3 orbital phases after periastron passage, similarly to the case of LS I +61 303. In addition, the optical photometric light curve maxima in AGL J2241+4454/MWC 656 occur ~0.25 phases passed periastron, similar to what is seen in LS I +61 303. We also find that the orbital eccentricity is correlated with orbital period for the known gamma-ray binaries. This is explained by the fact that small stellar separations are required for the efficient triggering of VHE radiation. Another correlation between the EW of Halpha and orbital period is also observed, similarly to the case of Be/X-ray binaries. These correlations are useful to provide estimates of the key orbital parameters Porb and e from the Halpha line in future Be gamma-ray binary candidates.

PSR B1259-63 is a 48 ms pulsar in a highly eccentric 3.4 year orbit around the young massive star LS 2883. During the periastron passage the system displays transient non-thermal unpulsed emission from radio to very high energy gamma rays. It is one of the three galactic binary systems clearly detected at TeV energies, together with LS 5039 and LS I +61 303. We observed PSR B1259-63 after the 2007 periastron passage with the Australian Long Baseline Array at 2.3 GHz to trace the milliarcsecond (mas) structure of the source at three different epochs. We have discovered extended and variable radio structure. The peak of the radio emission is detected outside the binary system near periastron, at projected distances of 10-20 mas (25-45 AU assuming a distance of 2.3 kpc). The total extent of the emission is ~50 mas (~120 AU). This is the first observational evidence that non-accreting pulsars orbiting massive stars can produce variable extended radio emission at AU scales. Similar structures are also seen in LS 5039 and LS I +61 303, in which the nature of the compact object is unknown. The discovery presented here for the young non-accreting pulsar PSR B1259-63 reinforces the link with these two sources and supports the presence of pulsars in these systems as well. A simple kinematical model considering only a spherical stellar wind can approximately trace the extended structures if the binary system orbit has a longitude of the ascending node of omega ~ -40 deg and a magnetization parameter of sigma~0.005.

LS 5039 is one of the few X-ray binaries detected at VHE, and potentially contains a young non-accreting pulsar. The outflow of accelerated particles emitting synchrotron emission can be directly mapped with high resolution radio observations. The morphology of the radio emission strongly depends on the properties of the compact object and on the orbital parameters of the binary system. We present VLBA observations of LS 5039 covering an orbital cycle, which show morphological and astrometric variability at mas scales. On the other hand, we discuss the possible association of LS 5039 with the supernova remnant SNR G016.8-01.1.

Context. With the latest infrared surveys, the number of massive protostellar candidates has increased significantly. New studies have posed additional questions on important issues about the formation, evolution, and other phenomena related to them. Complementary to infrared data, radio observations are a good tool to study the nature of these objects, and to diagnose the formation stage. Aims. Here we study the far-infrared source IRAS 16353-4636 with the aim of understanding its nature and origin. In particular, we search for young stellar objects (YSOs), possible outflow structure, and the presence of non-thermal emission. Methods. Using high-resolution, multi-wavelength radio continuum data obtained with the Australia Telescope Compact Array, we image IRAS 16353-4636 and its environment from 1.4 to 19.6 GHz, and derive the distribution of the spectral index at maximum angular resolution. We also present new JHKs photometry and spectroscopy data obtained at ESO NTT. 13 CO and archival HI line data, and infrared databases (MSX, GLIMPSE, MIPSGal) are also inspected. Results. The radio continuum emission associated with IRAS 16353-4636 was found to be extended (~10 arcsec), with a bow-shaped morphology above 4.8 GHz, and a strong peak persistent at all frequencies. The NIR photometry led us to identify ten near-IR sources and classify them according to their color. We used the HI line data to derive the source distance, and analyzed the kinematical information from the CO and NIR lines detected. Conclusions. We have identified the source IRAS 16353-4636 as a new protostellar cluster. In this cluster we recognized three distinct sources: a low-mass YSO, a high-mass YSOs, and a mildly confined region of intense and non-thermal radio emission. We propose the latter corresponds to the terminal part of an outflow.

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

Microquasars, X-ray binaries displaying relativistic jets driven by accretion onto a compact object, are some of the most efficient accelerators in the Galaxy. Theoretical models predict Very High Energy (VHE) emission at the base of the jet where particles are accelerated to multi-TeV energies. This emission could be detected by present IACTs. %Moreover, gamma-ray fluxes should increase during flaring events when accretion rates are enhanced. The MAGIC telescope observed the microquasars GRS 1915+105, Cyg X-3, Cyg X-1 and SS433 for ~ 150 hours in total from 2005 to 2008. We triggered our observations by using multi wavelength information through radio flaring alerts with the RATAN telescope as well as by ensuring the low/hard state of the source through RXTE/ASM and Swift/BAT monitoring data. We report on the upper limits on steady and variable emission from these sources over this long period.

XTE J1818-245 is an X-ray nova that experienced an outburst in 2005, first seen by the RXTE satellite. The source was observed simultaneously at various wavelengths up to soft gamma-rays with the INTEGRAL satellite, from 2005 February to September. X-ray novae are extreme systems that often harbor a black hole, and are known to emit throughout the electromagnetic spectrum when in outburst. We analyzed radio, (N)IR, optical, X-ray and soft gamma-ray observations and constructed simultaneous broad-band X-ray spectra. Analyzing both the light curves in various energy ranges and the hardness-intensity diagram enabled us to study the long-term behavior of the source. Spectral parameters were typical of the Soft Intermediate States and the High Soft States of a black hole candidate. The source showed relatively small spectral variations in X-rays with considerable flux variation in radio. Spectral studies showed that the accretion disc cooled down from 0.64 to 0.27 keV in 100 days and that the total flux decreased while the relative flux of the hot medium increased. Radio emission was detected several times, and, interestingly, five days after entering the HSS. Modeling the spectral energy distribution from the radio to the soft gamma-rays reveals that the radio flares arise from several ejection events. XTE J1818-245 is likely a black hole candidate transient source that might be closer than the Galactic Bulge. The results from the data analysis trace the physical changes that took place in the system at a maximum bolometric luminosity of (0.4-0.9)e38 erg/s (assuming a distance between 2.8-4.3 kpc) and they are discussed within the context of disc and jet models.

LS 5039 is one of the four TeV emitting X-ray binaries detected up to now. The powering source of its multi-wavelength emission can be accretion in a microquasar scenario or wind interaction in a young non-accreting pulsar scenario. These two scenarios predict different morphologic and peak position changes along the orbital cycle of 3.9 days, which can be tested at milliarcsecond scales using VLBI techniques. Here we present a campaign of 5 GHz VLBA observations conducted in June 2000 (2 runs five days apart). The results show a core component with a constant flux density, and a fast change in the morphology and the position angle of the elongated extended emission, but maintaining a stable flux density. These results are difficult to fit comfortably within a microquasar scenario, whereas they appear to be compatible with the predicted behavior for a non-accreting pulsar.

The source 4U 2206+54 is one of the most enigmatic high-mass X-ray binaries. In spite of intensive searches, X-ray pulsations have not been detected in the time range 0.001-1000 s. A cyclotron line at ~30 keV has been suggested by various authors but never detected with significance. The stellar wind of the optical companion is abnormally slow. The orbital period, initially reported to be 9.6 days, disappeared and a new periodicity of 19.25 days emerged. Our new long and uninterrupted RXTE observations allow us to search for long (~1 hr) pulsations for the first time. We have discovered 5560-s pulsations in the light curve of 4U 2206+54. Initially detected in RXTE data, these pulsations are also present in INTEGRAL and EXOSAT observations. The average X-ray luminosity in the energy range 2-10 keV is 1.5 x 10^{35} erg s^{-1} with a ratio Fmax/Fmin ~ 5. This ratio implies an eccentricity of ~0.4, somewhat higher than previously suggested. The source also shows a soft excess at low energies. If the soft excess is modelled with a blackbody component, then the size and temperature of the emitting region agrees with its interpretation in terms of a hot spot on the neutron star surface. The source displays variability on time scales of days, presumably due to changes in the mass accretion rate as the neutron star moves around the optical companion in a moderately eccentric orbit.

LS 5039 is one of the four TeV emitting X-ray binaries detected up to now. The powering source of its multi-wavelength emission can be accretion in a microquasar scenario or wind interaction in a young non-accreting pulsar scenario. These two scenarios predict different morphologic and peak position changes along the orbital cycle of 3.9 days, which can be tested at milliarcsecond scales using VLBI techniques. Here we present a campaign of 5 GHz VLBA observations conducted in June 2000 (2 runs five days apart). The results show a core component with a constant flux density, and a fast change in the morphology and the position angle of the elongated extended emission, but maintaining a stable flux density. These results are difficult to fit comfortably within a microquasar scenario, whereas they appear to be compatible with the predicted behavior for a non-accreting pulsar.

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.

Context. LS 5039 is one of the few TeV emitting X-ray binaries detected so far. The powering source of its multiwavelength emission can be accretion in a microquasar scenario or wind interaction in a young nonaccreting pulsar scenario. Aims. To present new high-resolution radio images and compare them with the expected behavior in the different scenarios. Methods. We analyze Very Long Baseline Array (VLBA) radio observations that provide morphological and astrometric information at milliarcsecond scales. Results. We detect a changing morphology between two images obtained five days apart. In both runs there is a core component with a constant flux density, and an elongated emission with a position angle (PA) that changes by 12+/-3 degrees between both runs. The source is nearly symmetric in the first run and asymmetric in the second one. The astrometric results are not conclusive. Conclusions. A simple and shockless microquasar scenario cannot easily explain the observed changes in morphology. An interpretation within the young nonaccreting pulsar scenario requires the inclination of the binary system to be very close to the upper limit imposed by the absence of X-ray eclipses.

Supergiant Fast X-ray Transients are obviously related to persistent Supergiant X-ray Binaries. Any convincing explanation for their behaviour must consistently take into account all types of X-ray sources powered by wind accretion. Here we present a common framework for wind accreting sources, within the context of clumpy wind models, that allows a coherent interpretation of their different behaviours as an immediate consequence of diverse orbital geometries.

The window of TeV Gamma-Ray Astrophysics was opened less than two decades ago, when the Crab Nebula was detected for the first time. After several years of development, the technique used by imaging atmospheric Cherenkov telescopes like HESS, MAGIC or VERITAS, is now allowing to conduct sensitive observations in the TeV regime. Water Cherenkov instruments like Milagro are also providing the first results after years of integration time. Different types of extragalactic and galactic sources have been detected, showing a variety of interesting phenomena that are boosting theory in very high energy gamma-ray astrophysics. Here I review some of the most interesting results obtained up to now, making special emphasis in the field of X-ray/gamma-ray binaries.

We report on the results from the observations in very high energy band (VHE, E_gamma > 100GeV) of the black hole X-ray binary (BHXB) Cygnus X-1. The observations were performed with the MAGIC telescope, for a total of 40 hours during 26 nights, spanning the period between June and November 2006. We report on the results of the searches for steady and variable gamma-ray signals, including the first experimental evidence for an intense flare, of duration between 1.5 and 24 hours.

LS 5039 has been observed with several X-ray instruments so far. The source presents X-ray variability at orbital timescales in flux and photon index. The system harbors an O-type main sequence star with moderate mass loss. At present, the link between the X-rays and the stellar wind is unclear. New XMM-Newton observations have been performed around periastron and apastron passages in September 2005, when the stellar wind activity was apparently higher. April 2005 Chandra observations on LS 5039 are also revisited. Moreover, a compilation of Halpha EW data obtained since 1992, from which the stellar mass loss evolution can be approximately inferred, is carried out. XMM-Newton observations show higher and harder emission around apastron than around periastron. No signatures of thermal emission or a reflection iron line indicating the presence of an accretion disk are found in the spectrum, and the hydrogen column density (NH) is compatible with being the same in both observations and consistent with the interstellar value. 2005 Chandra observations show a hard X-ray spectrum, and possibly high fluxes. The Halpha EW shows yearly variations of a ~ 10%, and does not seem to be correlated with X-ray fluxes obtained at similar phases, unlike it would be expected in the wind accretion scenario. 2005 XMM-Newton and Chandra observations are consistent with 2003 RXTE/PCA results. The constancy of the NH seems to imply that either the X-ray emitter is located at >~ 1012 cm from the compact object, or the real N_H is 3-27 times smaller than the one predicted by a spherical symmetric wind model. We suggest that the multiwavelength non-thermal emission of LS 5039 is related to the radio jets and unlikely produced inside the binary system.

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.

The space velocity of a stellar black hole encodes the history of its formation and evolution. Here we measure the 3-dimensional motion of the microquasar GRS 1915+105, using a decade of astrometry with the NRAO Very Long Baseline Array, together with the published radial velocity. The velocity in the Galactic Plane deviates from circular rotation by 53-80 +_ 8 km/s, where the range covers any specific distance from 6-12 kpc. Perpendicular to the plane, the velocity is only 10 +_ 4 km/s. The peculiar velocity is minimized at a distance 9-10 kpc, and is then nearly in the radial direction towards the Galactic Center. We discuss mechanisms for the origin of the peculiar velocity, and conclude that it is most likely a consequence of Galactic velocity diffusion on this old binary, rather than the result of a supernova kick during the formation of the 14 Mo black hole. Finally, a brief comparison is made with 4 other BH binaries whose kinematics are well determined.

1WGA J1346.5-6255 is a ROSAT X-ray source found within the radio lobes of the supernova remnant (SNR) G309.2-00.6. This source also appears to coincide with the bright and early-type star HD 119682, which is in the middle of the galactic open cluster NGC 5281. The radio morphology of the remnant, consisting of two brightened and distorted arcs of emission on opposite sides of the 1WGA J1346.5-6255 source and of a jet-like feature and break in the shell, led to the suggestion that 1WGA J1346.5-6255/G309.2-00.6 is a young analog of the microquasar SS 433 powering the W50 nebula. This motivated us to study this source at X-ray and optical wavelengths. We here present new Chandra observations of 1WGA J1346.5-6255, archival XMM-Newton observations of G309.2-00.6, and optical spectroscopic observations of HD 119682, in order to search for X-ray jets from 1WGA J1346.5-6255, study its association with the SNR, and test for whether HD 119682 represents its optical counterpart. We do not find evidence for jets from 1WGA J1346.5-6255 down to an unabsorbed flux of 2.6E-13 ergs/cm2/s (0.5-7.5 keV), we rule out its association with G309.2-00.6, and we confirm that HD 119682 is its optical counterpart. We derive a distance of 1.2+/-0.3 kpc, which is consistent with the distance estimate to NGC 5281 (1.3+/-0.3 kpc), and much smaller than the distance derived to the SNR G309.2-00.6. We discuss the nature of the source, unveil that HD 119682 is a Be star and suggest it is a new member of the recently proposed group of $\gamma$-Cas analogs. The Chandra and XMM X-ray lightcurves show variability on timescales of hundreds of seconds, and the presence of a possible period of about 1500 s that could be the rotational period of an accreting neutron star or white dwarf in this $\gamma$-Cas analog.

We report the results of simultaneous multiwavelength observations of the X-ray transient source SWIFT J1753.5-0127 performed with INTEGRAL, RXTE, NTT, REM and VLA on 2005 August 10-12. The source, which underwent an X-ray outburst since 2005 May 30, was observed during the INTEGRAL Target of Opportunity program dedicated to new X-ray novae located in the Galactic Halo. Broad-band spectra and fast timing variability properties of SWIFT J1753.5-0127 are analyzed together with the optical, near infra-red and radio data. We show that the source was significantly detected up to 600 keV with Comptonization parameters and timing properties typical of the so-called Low/Hard State of black hole candidates. We build a spectral energy distribution and we show that SWIFT J1753.5-0127 does not follow the usual radio/X-ray correlation of X-ray binaries in the Low/Hard State. We give estimates of distance and mass. We conclude that SWIFT J1753.5-0127 belongs to the X-ray nova class and that it is likely a black hole candidate transient source of the Galactic Halo which remained in the Low/Hard State during its main outburst. We discuss our results within the context of Comptonization and jet models.

We report high-energy results obtained with INTEGRAL and Rossi-XTE on two microquasars: the persistent high-mass system Cygnus X-1 and the transient low-mass binary SWIFT J1753.5-0127. INTEGRAL observed Cygnus X-1 from 2002 to 2004: the spectral (5-1000 keV) properties of the source, seen at least in three distinct spectral states, show disc and corona changes. In 2003 June, a high-energy tail at several hundred keV in excess of the thermal Comptonization model was observed, suggesting the presence of an additional non-thermal component. At that time, we detected an unusual correlation between radio data and high-energy hardness. We also report and compare the results obtained with simultaneous observations of the transient source SWIFT J1753.5-0127 performed with Rossi-XTE, INTEGRAL, VLA, REM and NTT on 2005 August 10-12 near its hard X-ray outburst. Broad-band spectra and fast time-variability properties are derived on this source (probably located in the galactic halo) together with radio, IR and optical data. We build a spectral energy distribution of the source and derive interesting multiwavelength constraints. Significantly detected up to 600 keV in a typical Low/Hard State, the transient does not seem to follow the usual radio/X-ray correlation.

LS 5039 is the only X-ray binary persistently detected at TeV energies by the Cherenkov HESS telescope. It is moreover a gamma-ray emitter in the GeV and possibly MeV energy ranges. To understand important aspects of jet physics, like the magnetic field content or particle acceleration, and emission processes, such as synchrotron and inverse Compton (IC), a complete modeling of the multiwavelength data is necessary. LS 5039 has been detected along almost all the electromagnetic spectrum thanks to several radio, infrared, optical and soft X-ray detections. However, hard X-ray detections above 20 keV have been so far elusive and/or doubtful, partly due to source confusion for the poor spatial resolution of hard X-ray instruments. We report here on deep (300 ksec) serendipitous INTEGRAL hard X-ray observations of LS 5039, coupled with simultaneous VLA radio observations. We obtain a 20-40 keV flux of 1.1 +/- 0.3 mCrab (5.9 (+/-1.6) X 10^{-12} erg cm^{-2} s^{-1}), a 40-100 keV upper limit of 1.5 mCrab (9.5 x 10^{-12} erg cm^{-2}s^{-1}), and typical radio flux densities of about 25 mJy at 5GHz. These hard X-ray fluxes are significantly lower than previous estimates obtained with BATSE in the same energy range but, in the lower interval, agree with extrapolation of previous RXTE measurements. The INTEGRAL observations also hint to a break in the spectral behavior at hard X-rays. A more sensitive characterization of the hard X-ray spectrum of LS 5039 from 20 to 100 keV could therefore constrain key aspects of the jet physics, like the relativistic particle spectrum and the magnetic field strength. Future multiwavelength observations would allow to establish whether such hard X-ray synchrotron emission is produced by the same population of relativistic electrons as those presumably producing TeV emission through IC.

A new sample of hard X-ray sources in the Galactic Plane is being revealed by the regular observations performed by the INTEGRAL satellite. The full characterization of these sources is mandatory to understand the hard X-ray sky. Here we report new multifrequency radio, infrared and optical observations of the source IGR J18027-1455, as well as a multi-wavelength study from radio to hard X-rays. The radio counterpart of IGR J18027-1455 is not resolved at any observing frequency. The radio flux density is well fitted by a simple power law with a spectral index alpha=-0.75+/-0.02. This value is typical of optically thin non-thermal synchrotron emission originated in a jet. The NIR and optical spectra show redshifted emission lines with z=0.034, and a broad Halpha line profile with FWHM ~3400 km/s. This suggests an Active Galactic Nucleus (AGN) of type 1 as the optical counterpart of IGR J18027-1455. We confirm the Seyfert 1 nature of the source, which is intrinsically bright at high energies both in absolute terms and when scaled to a normalized 6 cm luminosity. Finally, comparing its X-ray luminosity with isotropic indicators, we find that the source is Compton thin and AGN dominated. This indicates that INTEGRAL might have just seen the tip of the iceberg, and several tens of such sources should be unveiled during the course of its lifetime.

Massive X-ray binaries are usually classified depending on the properties of the donor star in classical, supergiant and Be X-ray binaries. The massive X-ray binary 4U 2206+54 does not fit in any of these groups, and deserves a detailed study to understand how the transfer of matter and the accretion on to the compact object take place. To this end we study an IUE spectrum of the donor and obtain a wind terminal velocity (v_inf) of ~350 km/s, which is abnormally slow for its spectral type. We also analyse here more than 9 years of available RXTE/ASM data. We study the long-term X-ray variability of the source and find it to be similar to that observed in the wind-fed supergiant system Vela X-1, reinforcing the idea that 4U 2206+54 is also a wind-fed system. We find a quasi-period decreasing from ~270 to ~130 d, noticed in previous works but never studied in detail. We discuss possible scenarios and conclude that long-term quasi-periodic variations in the mass-loss rate of the primary are probably driving such variability in the measured X-ray flux. We obtain an improved orbital period of 9.5591 d with maximum X-ray flux at MJD 51856.6. Our study of the orbital X-ray variability in the context of wind accretion suggests a moderate eccentricity around 0.15. Moreover, the low value of v_inf solves the long-standing problem of the relatively high X-ray luminosity for the unevolved nature of the donor, BD +53 2790, which is probably an O9.5 V star. We note that changes in v_inf and/or the mass-loss rate of the primary alone cannot explain the diferent patterns displayed by the orbital X-ray variability. We finally emphasize that 4U 2206+54, together with LS 5039, could be part of a new population of wind-fed HMXBs with main sequence donors, the natural progenitors of supergiant X-ray binaries. (Abridged)

The population of high energy and very high energy gamma-ray sources, detected with EGRET and the new generation of ground-based Cherenkov telescopes, conforms a reduced but physically important sample. Most of these sources are extragalactic (e.g., blazars), while among the galactic ones there are pulsars and SN remnants. The microquasar LS 5039, previously proposed to be associated with an EGRET source by Paredes et al. (2000), has recently been detected at TeV energies, confirming that microquasars should be regarded as a class of high energy gamma-ray sources. To model and understand how the energetic photons are produced and escape from LS 5039 it is crucial to unveil the nature of the compact object, which remains unknown. Here we present new intermediate-dispersion spectroscopy of this source which, combined with values reported in the literature, provides an orbital period of 3.90603+/-0.00017 d, a mass function f(M)=0.0053+/-0.0009 M_sun, and an eccentricity e=0.35+/-0.04. Atmosphere model fitting to the spectrum of the optical companion, together with our new distance estimate of d=2.5+/-0.1 kpc, yields R_opt=9.3+0.7-0.6 R_sun, log (L_opt/L_sun)=5.26+/-0.06, and M_opt=22.9+3.4-2.9 M_sun. These, combined with our dynamical solution and the assumption of pseudo-synchronization, yield an inclination i=24.9+/-2.8 degree and a compact object mass M_X=3.7+1.3-1.0 M_sun. This is above neutron star masses for most of the standard equations of state and, therefore, we propose that the compact object in LS 5039 is a black hole. We finally discuss about the implications of our orbital solution and new parameters of the binary system on the CNO products, the accretion/ejection energetic balance, the SN explosion scenario, and the behaviour of the TeV emission with the new orbital period.

The properties of the orbit and the donor star in the high mass X-ray binary microquasar LS 5039 indicate that accretion processes should mainly occur via a radiatively driven wind. In such a scenario, significant X-ray variability would be expected due to the eccentricity of the orbit. The source has been observed at X-rays by several missions, although with a poor coverage that prevents to reach any conclusion about orbital variability. Therefore, we conducted RossiXTE observations of the microquasar system LS 5039 covering a full orbital period of 4 days. Individual observations are well fitted with an absorbed power-law plus a Gaussian at 6.7 keV, to account for iron line emission that is probably a diffuse background feature. In addition, we have taken into account that the continuum is also affected by significant diffuse background contamination. Our results show moderate power-law flux variations on timescales of days, as well as the presence of miniflares on shorter timescales. The new orbital ephemeris of the system recently obtained by Casares et al. have allowed us to show, for the first time, that an increase of emission is seen close to the periastron passage, as expected in an accretion scenario. Moreover, the detected orbital variability is a factor of ~4 smaller than the one expected by using a simple wind accretion model, and we suggest that an accretion disk around the compact object could be responsible for this discrepancy. On the other hand, significant changes in the photon index are also observed clearly anti-correlated with the flux variations. We interpret the overall X-ray spectral characteristics of LS 5039 in the context of X-ray radiation produced by inverse Compton and/or synchrotron processes in the jet of this microquasar.

The majority of High Mass X-ray Binaries (HMXBs) behave as X-ray pulsars, revealing that they contain a magnetised neutron star. Among the four HMXBs not showing pulsations, and that do not show the characteristics of accreting black holes, there is the unusual HMXB 4U 2206+54. Here we present contemporaneous high-energy and radio observations of this system conducted with INTEGRAL and the VLA in order to unveil its nature. The high-energy spectra show clear indications of the presence of an absorption feature at ~32 keV. This is the third high-energy observatory which reveals marginal evidence of this feature, giving strong support to the existence of a cyclotron resonance scattering feature, which implies a magnetic field of 3.6 x 10^12 G. On the other hand, the source is not detected at centimetre radio wavelengths with a 3 sigma upper limit of 0.039 mJy. The expected radio emission for an accreting black hole in the low/hard state, inferred from X-ray flux measurements, would be at least 60 times greater than the measured upper limit. Both results firmly indicate that, in spite of the absence of pulsations, 4U 2206+54 hosts a magnetic accreting neutron star, the first one not to be observed as an X-ray pulsar.

We present multiepoch VLBA observations of the compact jet of GRS 1915+105 conducted at 15.0 and 8.4 GHz during a {\it plateau} state of the source in 2003 March-April. These observations show that the compact jet is clearly asymmetric. Assuming an intrinsically symmetric continuous jet flow, using Doppler boosting arguments and an angle to the line of sight of $\theta=70\degr$, we obtain values for the velocity of the flow in the range 0.3--0.5$c$. These values are much higher than in previous observations of such compact jet, although much lower than the highly relativistic values found during individual ejection events. These preliminary results are compatible with current ideas on the jet flow velocity for black holes in the low/hard state.

We present the first optical polarimetric observations of the runaway microquasar LS 5039. Our results reveal the presence of a large amount (~5%) of polarized emission towards this binary system. By combining polarimetric and spectroscopic observations of some stars in the field together with available statistical information on the galactic interstellar polarization of the region, we have estimated and subtracted the contribution of the interstellar polarization in this direction. As a result, we obtain an intrinsic polarization of ~3% for the object, much higher than what would be expected from jet emission in the optical domain. We suggest that the polarized light originates by electron Thomson scattering in the stellar envelope of the companion star. This allows us to constrain the size of the region where the polarized emission originates, as well as estimating the scattering electronic density and the wind velocity at such distance.

The discovery of the microquasar LS 5039 well within the 95% conficence contour of the Unidentified EGRET Source (UES) 3EG J1824-1514 was a major step towards the possible association between microquasars (MQs) and UESs. The recent discovery of precessing relativistic radio jets in LS I +61 303, a source associated for long time with 2CG 135+01 and with the UES 3EG J0241+6103, has given further support to this idea. Finally, the very recently proposed association between the microquasar candidate AX J1639.0-4642 and the UES 3EG J1639-4702 points towards a population of High Mass X-ray Binary (HMXB)/Neutron Star (NS) microquasars as counterparts of low-latitude unidentified EGRET sources.

We present a multiwavelength study of the environment of the unidentified X-ray/gamma-ray sources IGR J18027-1455 and IGR J21247+5058, recently discovered by the IBIS/ISGRI instrument, onboard the INTEGRAL satellite. The main properties of the sources found inside their position error circles, give us clues about the nature of these high-energy sources.

We present a multiwavelength study of the field containing the unidentified X-ray source AX J1639.0-4642, discovered with the ASCA observatory and recently detected with the IBIS telescope, onboard the INTEGRAL satellite, dubbed IGR J16393-4643. The huge hydrogen column density towards the source, the hard spectral index in the 0.7-10 keV band and its flux variability suggest that the source is a High Mass X-ray Binary (HMXB) enshrouded by dust. Our search reveals the presence of a non-thermal radio counterpart within the X-ray error box. After a study of the broadband emission from X-rays to the radio domain, we propose that AX J1639.0-4642 is a dust-enshrouded Microquasar (MQ) candidate. In addition, the X-ray source is well within the 95% location contour of the unidentified gamma-ray source 3EG J1639-4702. The main properties of AX J1639.0-4642/3EG J1639-4702 are consistent with those of two other MQs previously proposed to display high-energy gamma-ray emission.

We present the international collaboration MINE (Multi-lambda Integral NEtwork) aimed at conducting multi-wavelength observations of X-ray binaries and microquasars simultaneously with the INTEGRAL gamma-ray satellite. We will focus on the 2003 March-April campaign of observations of the peculiar microquasar GRS 1915+105 gathering radio, IR and X-ray data. The source was observed 3 times in the plateau state, before and after a major radio and X-ray flare. It showed strong steady optically thick radio emission corresponding to powerful compact jets resolved in the radio images, bright near-infrared emission, a strong QPO at 2.5 Hz in the X-rays and a power law dominated spectrum without cutoff in the 3-300 keV range. We compare the different observations, their multi-wavelength light curves, including JEM-X, ISGRI and SPI, and the parameters deduced from fitting the spectra obtained with these instruments on board INTEGRAL.

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.

I partially review our current knowledge on microquasars, making special emphasis on radio interferometric observations, and I comment on emerging trends in this field of astrophysics.

We present here the results of a search for new microquasars at low galactic latitudes, based on a cross-identification between the ROSAT all sky Bright Source Catalog (RBSC) and the NRAO VLA Sky Survey (NVSS) and follow-up observations. The results obtained up to now suggest that persistent/silent microquasars such as LS 5039 are rare objects in our Galaxy, and indicate that future deeper surveys, and harder than the RBSC in X-rays, will play a fundamental role in order to discover them.

Recent studies of relativistic jet sources in the Galaxy, also known as microquasars, have been very useful in trying to understand the accretion/ejection processes that take place near compact objects. However, the number of sources involved in such studies is still small. In an attempt to increase the number of known microquasars we have carried out a search for new Radio Emitting X-ray Binaries (REXBs). These sources are the ones to be observed later with VLBI techniques to unveil their possible microquasar nature. To this end, we have performed a cross-identification between the X-ray ROSAT all sky survey Bright Source Catalog (RBSC) and the radio NRAO VLA Sky Survey (NVSS) catalogs under very restrictive selection criteria for sources with |b|<5 degrees. We have also conducted a deep observational radio and optical study for six of the selected candidates. At the end of this process two of the candidates appear to be promising, and deserve additional observations aimed to confirm their proposed microquasar nature.

We report optical spectroscopic observations of a sample of 6 low-galactic latitude microquasar candidates selected by cross-identification of X-ray and radio point source catalogs for |b|<5 degrees. Two objects resulted to be of clear extragalactic origin, as an obvious cosmologic redshift has been measured from their emission lines. For the rest, none exhibits a clear stellar-like spectrum as would be expected for genuine Galactic microquasars. Their featureless spectra are consistent with being extragalactic in origin although two of them could be also highly reddened stars. The apparent non-confirmation of our candidates suggests that the population of persistent microquasar systems in the Galaxy is more rare than previously believed. If none of them is galactic, the upper limit to the space density of new Cygnus X-3-like microquasars within 15 kpc would be 1.1\times10^{-12} per cubic pc. A similar upper limit for new LS 5039-like systems within 4 kpc is estimated to be 5.6\times10^{-11} per cubic pc.

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 |b|<5 degrees, 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.

We present the international collaboration MINE (Multi-lambda INTEGRAL NEtwork) aimed at conducting multi-wavelength observations of microquasars simultaneously with the INTEGRAL satellite. The first results on GRS 1915+105 are encouraging and those to come should help us to understand the physics of the accretion and ejection phenomena around a compact object.

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.

We present the international collaboration MINE (Multi-lambda Integral NEtwork) aimed at conducting multi-wavelength observations of X-ray binaries and microquasars simultaneously with INTEGRAL. The first results on GRS 1915+105 are encouraging and those to come should help us to understand the physics of the accretion and ejection phenomena around a compact object. A collaboration such as MINE could be very useful for observing quasars and microquasars simultaneously with HESS and later with ANTARES.

We present the result of multi-wavelength observations of the microquasar GRS 1915+105 in a plateau state with a luminosity of ~7.5x10^{38) erg s-1 (~40% L_Edd), conducted simultaneously with the INTEGRAL and RXTE satellites, the ESO/NTT, the Ryle Telescope, the NRAO VLA and VLBA, in 2003 April 2-3. For the first time were observed concurrently in GRS 1915+105 all of the following properties: a strong steady optically thick radio emission corresponding to a powerful compact jet resolved with the VLBA, bright near-IR emission, a strong QPO at 2.5 Hz in the X-rays and a power law dominated spectrum without any cutoff in the 3-400 keV range.

We report on the results of the spectral and timing analysis of a BeppoSAX observation of the microquasar system LS 5039/RX J1826.2-1450. The source was found in a low-flux state with Fx(1-10 keV)= 4.7 x 10^{-12} erg cm^{-2} s^{-1}, which represents almost one order of magnitude lower than a previous RXTE observation 2.5 years before. The 0.1--10 keV spectrum is described by an absorbed power-law continuum with photon-number spectral index Gamma=1.8+-0.2 and hydrogen column density of NH=1.0^{+0.4}_{-0.3} x 10^{22} cm^{-2}. According to the orbital parameters of the system the BeppoSAX observation covers the time of an X-ray eclipse should one occur. However, the 1.6-10 keV light curve does not show evidence for such an event, which allows us to give an upper limit to the inclination of the system. The low X-ray flux detected during this observation is interpreted as a decrease in the mass accretion rate onto the compact object due to a decrease in the mass-loss rate from the primary.

We present new XMM-Newton X-ray observations of the z=1.246 narrow-line radioquasar RX J1011.2+5545 serendipitously discovered by ROSAT. The flat X-ray spectrum previously measured by ROSAT and ASCA is shown to be the result of a steep Gamma~1.8 power law spectrum seen through a moderate intrinsic absorbing column NH~4E21 cm^-2. The position of the X-ray source is entirely coincident with the nucleus of the radio source that we have resolved in new sensitive VLA observations at 3.6 and 6 cm, implying that scattering in the radio lobes is not responsible for the bulk of X-ray emission. In the EPIC pn image, a faint patch of X-ray emission is apparent 14'' to the NE of the main X-ray source. The former is positionally coincident with an apparently extended optical object with R~21.9, but there is no associated radio emission, thus ruling out the possibility that this represents a hotspot in a jet emanating from the primary X-ray source. No reflection features are detected in the X-ray spectrum of the narrow-line radioquasar, although an Fe line with equivalent width of up to 600 eV cannot be ruled out.

We present multiepoch Very Large Array (VLA) observations at 1.4 GHz, 4.9 GHz, 8.5 GHz and 14.9 GHz for a sample of eight RS CVn binary systems. Circular polarization measurements of these systems are also reported. Most of the fluxes observed are consistent with incoherent emission from mildly relativistic electrons. Several systems show an increase of the degree of circular polarization with increasing frequency in the optically thin regime, in conflict with predictions by gyrosynchrotron models. We observed a reversal in the sense of circular polarization with increasing frequency in three non-eclipsing systems: EI Eri, DM Uma and HD 8358. We find clear evidence for coherent plasma emission at 1.4 GHz in the quiescent spectrum of HD 8358 during the helicity reversal. The degrees of polarization of the other two systems could also be accounted for by a coherent emission process. The observations of ER Vul revealed two U-shaped flux spectra at the highest frequencies. The U-shape of the spectra may be accounted for by an optically thin gyrosynchrotron source for the low frequency part whereas the high frequency part is dominated by a thermal emission component.

We report the discovery of a new radio galaxy within the location error box of the gamma-ray source 3EG J1735-1500. The galaxy is a double-sided jet source forming a large angle with the line of sight. Optical observations reveal a V ~ 18 magnitude galaxy at the position of the radio core. Although the association with the EGRET source is not confirmed at the present stage, because there is a competing, alternative gamma-ray candidate within the location error contours which is also studied here, the case deserves further attention. The new radio galaxy can be used to test the recently proposed possibility of gamma-ray emitting radio galaxies beyond the already known case of Centaurus A.

We present here new observations conducted with the EVN and MERLIN of the persistent microquasar LS 5039 discovered by Paredes et al. (2000) with the VLBA. The new observations confirm the presence of an asymmetric two-sided jet reaching up to 1000 AU on the longest jet arm. The results suggest a bending of the jets with increasing distance from the core and/or precession. The origin and location of the high-energy gamma-ray emission associated with the system is discussed and an estimate of the magnetic field at the base of the jet given. Our results suggest a well collimated radio jet. We also comment on new observing strategies to be used with satellites and forthcoming detectors, since this persistent source appears to be a rather good laboratory to explore the accretion/ejection processes taking place near compact objects.

The unidentified EGRET source 3EG J1828+0142, located at (l,b) ~ (31.9,5.8) (Hartman et al. 1999), is particularly enigmatic since it displays very high variability over timescales of months (Torres et al. 2000, Tompkins 1999) yet no blazar has been found within the 95% confidence contour of the gamma-ray location error-box (Punsly et al 2000). It also has an unusually steep gamma-ray spectral index of Gamma=2.76 +/- 0.39. It's relatively low latitude may be indicative of a galactic origin of the variable GeV flux. A multifrequency study of this source has been initiated and may reveal the presence of a novel type of compact object: a galactic 'microblazar'. X-ray data, especially those of high spatial resolution from CHANDRA, would be particularly useful in establishing such a connection.

The microquasar nature of LS5039 was revealed by May 1999 VLBA+VLA observations showing a two-sided jet at milliarcsecond scales. Here we present follow-up interferometric observations carried out with the EVN and MERLIN at 5 GHz in March 2000. The obtained maps with both the EVN and MERLIN show a two-sided jet with a similar position angle to the previous VLBA+VLA map. The total length of the jet arms is ~60 mas in the EVN map and ~300 mas in the MERLIN map. A brightness and length asymmetry of the jets, compatible with the earlier observations, is also present in the maps. Overall, these observations confirm the existence of a two-sided jet structure in LS 5039 and seem to indicate their persistent nature.

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 have compiled optical and radio astrometric data of the microquasar LS 5039 and derived its proper motion. This, together with the distance and radial velocity of the system, allows us to state that this source is escaping from its own regional standard of rest, with a total systemic velocity of about 150 km/s and a component perpendicular to the galactic plane larger than 100 km/s. This is probably the result of an acceleration obtained during the supernova event that created the compact object in this binary system. We have computed the trajectory of LS 5039 in the past, and searched for OB associations and supernova remnants in its path. In particular, we have studied the possible association between LS 5039 and the supernova remnant G016.8-01.1, which, despite our efforts, remains dubious. We have also discovered and studied an HI cavity in the ISM, which could have been created by the stellar wind of LS 5039 or by the progenitor of the compact object in the system. Finally, in the symmetric supernova explosion scenario, we estimate that at least 17 solar masses were lost in order to produce the high eccentricity observed. Such a mass loss could also explain the observed runaway velocity of the microquasar.

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.

Microquasars are stellar x-ray binaries that behave as a scaled down version of extragalactic quasars. The star LS 5039 is a new microquasar system with apparent persistent ejection of relativistic plasma at a 3 kiloparsec distance from the sun. It may also be associated with a gamma-ray source discovered by the Energetic Gamma Ray Experiment Telescope (EGRET) on board the COMPTON-Gamma Ray Observatory satellite. Before the discovery of LS 5039, merely a handful of microquasars had been identified in the Galaxy, and none of them was detected in high-energy gamma-rays.

We present the results of analyzing H$\alpha$ spectra of the radio emitting X-ray binary LS I+61303. For the first time, the same 26.5 d radio period is clearly detected in the H$\alpha$ emission line. Moreover, the equivalent width and the peak separation of the H$\alpha$ emission line seem also to vary over a time scale of 1600 days. This points towards the $\sim4$ yr modulation, detected in the radio outburst amplitude, being probably a result of variations in the mass loss rate of the Be star and/or density variability in the circumstellar disk. In addition, the dependence of the peak separation from the equivalent width informs us that the LS I+61303 circumstellar disk is among the densest of Be-stars.

RX J1826.2-1450/LS 5039 has been recently proposed to be a radio emitting high mass X-ray binary. In this paper, we present an analysis of its X-ray timing and spectroscopic properties using different instruments on board the RXTE satellite. The timing analysis indicates the absence of pulsed or periodic emission on time scales of 0.02-2000 s and 2-200 d, respectively. The source spectrum is well represented by a power-law model, plus a Gaussian component describing a strong iron line at 6.6 keV. Significant emission is seen up to 30 keV, and no exponential cut-off at high energy is required. We also study the radio properties of the system according to the GBI-NASA Monitoring Program. RX J1826.2-1450/LS 5039 continues to display moderate radio variability with a clearly non-thermal spectral index. No strong radio outbursts have been detected after several months.