We study the photometric behavior of the recurrent nova RS Oph by 58 monitoring light curves (LCs), taken by 5 telescopes. All LCs show repeating time structures with some quasi-periods (QPs) in time scales from minutes to hours. In our previous work 97 QPs were detected in the LCs by local minimums of structure functions and local maximums of auto-correlation functions. The distribution of the QPs shows modes at 8, 13, 21, 30, 48 and 73 min, where the mode at 8 min is poorly unveiled. These modes follow a power function with base $1.55\approx 3/2$ with standard deviation 4.7%. This function predicts modes also at 5.3 and 3.5 min, which are not detected in the full MLCs. In the present work we analyze simple small parts from high resolution LCs. We confirm the QPs modes at 8.0, 5.3 and 3.5 min. Generally, we found 8 QP modes with regular logarithmic distribution in the time interval 3.5-73 min. We also show typical intra-night evolutions of QP modes in the minute scale -- sharp or gradual transitions from one QP mode to other. In the end we find that the parts of the LCs carry out the properties of the whole LCs at short time scale. This lead to two well pronounces dependences - between the range deviation and standard deviation of the LC, as well as between the quasi-period and the relevant level of the density function of the LC.

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 analyzed 29 pairs of time series in B and V bands of the recurrent nova RS Oph. The observations were carried out in 2008-2017 with duration 0.6 - 3.6 hours, with time resolution 0.5 - 3.3 min. We scanned digitally each series by data windows with various sizes Theta and derived two of the simplest fractal parameters for every Theta - standard deviation D and structural deviation S. Using the local minima of the structural function log S = f_S(log Theta) we unveiled 80 time structures, 42 in B band and 38 in V band, with time sizes 10-120 min. About 3/4 of the time sizes belong to the interval 10-40 min and about 1/4 lie in the interval 60-120 min. The respective cycles per day are 144-36 c/d and 24-15 c/d. On logarithmic scale, the distribution of the time sizes shows maximums at about 10, 21, 36 and 74 min. The 10 min flickering is poorly detectable in our series and we found the most widespread time structures (in about 1/5 of the cases) have time sizes about 21 min (about 69 c/d). Using the deviation function log D = f_D(log Theta) we estimated the relative cumulative energy (including the energy of the shorter structures in it), associated with the detected structure sizes, to be in the interval of the relative fluxes 2-11 %. The energies correlate weakly with the logarithms of the structure sizes, with correlation coefficients 0.60 and 0.57, under slope coefficients 0.04 and 0.03 in B and V band, respectively. The distributions of the energies occur bimodal, with maximums about 4% and 6% in B band, as well as about 3% and 5% in V band. The left and right modes of the distributions may be associated with the structure sizes 10 - 21 min and 37 - 74 min, respectively.

We performed 48.6 hours (in 28 nights) of simultaneous B and V band observations of the flickering variability of the recurrent nova RS Oph in quiescence. During the time of our observations the brightness of the system varied between 13.2 > B > 11.1 and the colour in the range 0.86 < B-V < 1.33 . We find that RS~Oph becomes more blue, as it becomes brighter, however the hot component becomes more red as it becomes brighter (assuming that the red giant is non-variable). During all the runs RS Oph exhibits flickering with amplitude 0.16 - 0.59 mag in B band. For the flickering source we find that it has colour -0.14 < B-V < 0.40, temperature in the range 7200 < T_fl < 18900, and average radius 1.1 < R_fl < 6.7 R_sun. We do not find a correlation between the temperature of the flickering and the brightness. However, we do find a strong correlation (correlation coefficient 0.81, significance 1.1x10^{-7} ) between B band magnitude and the average radius of the flickering source - as the brightness of the system increases the size of the flickering source also increases. The estimated temperature is similar to that of the bright spot of cataclysmic variables. The persistent presence of flickering indicates that the white dwarf is actively accreting material for the next outburst.

Here we present deep (16 mumJy), very high (40 mas) angular resolution 1.14 mm, polarimetric, Atacama Large Millimeter/submillimeter Array (ALMA) observations towards the massive protostar driving the HH 80-81 radio jet. The observations clearly resolve the disk oriented perpendicular to the radio jet, with a radius of ~0.171 arcsec (~291 au at 1.7 kpc distance). The continuum brightness temperature, the intensity profile, and the polarization properties clearly indicate that the disk is optically thick for a radius of R<170 au. The linear polarization of the dust emission is detected almost all along the disk and its properties suggest that dust polarization is produced mainly by self-scattering. However, the polarization pattern presents a clear differentiation between the inner (optically thick) part of the disk and the outer (optically thin) region of the disk, with a sharp transition that occurs at a radius of 0.1 arcsec (~170 au). The polarization characteristics of the inner disk suggest that dust settling has not occurred yet with a maximum dust grain size between 50 and 500 mum. The outer part of the disk has a clear azimuthal pattern but with a significantly higher polarization fraction compared to the inner disk. This pattern is broadly consistent with self-scattering of a radiation field that is beamed radially outward, as expected in the optically thin outer region, although contribution from non-spherical grains aligned with respect to the radiative flux cannot be excluded.

Here we present quasi-simultaneous observations of the flickering of the symbiotic binary star CH Cyg in U, B and V bands. We calculate the flickering source parameters and discuss the possible reason for the flickering cessation in the period 2010-2013.

We report optical CCD photometry of the recently identified symbiotic star EF Aql. Our observations in Johnson V and B bands clearly show the presence of stochastic light variations with an amplitude of about 0.2 mag on a time scale of minutes. The observations point toward a white dwarf (WD) as the hot component in the system. It is the 11-th object among more than 200 symbiotic stars known with detected optical flickering. Estimates of the mass accretion rate onto the WD and the mass loss rate in the wind of the Mira secondary star lead to the conclusion that less than 1 per cent of the wind is captured by the WD. Eight further candidates for the detection of flickering in similar systems are suggested.

We estimate the mass of the compact object in the gamma-ray binaries LSI+61303 and MWC 148, using the latest data for the inclination, orbital motion and assuming that the orbital plane coincides with the equatorial plane of the Be star. For LSI+61303 we find the mass of the compact object to be most likely in the range 1.3 M_sun < M_2 < 2.0 M_sun, which means that it is probably a neutron star. For MWC 148, we find the mass of the compact object in a higher range, 2.1 M_sum < M_2 < 7.3 M_sun, which increases the chances for this system to host a black hole companion.

We report optical spectroscopic observations of the Be/gamma-ray binaries LSI+61303, MWC 148 and MWC 656. The peak separation and equivalent widths of prominent emission lines (H-alpha, H-beta, H-gamma, HeI, and FeII) are measured. We estimated the circumstellar disc size, compared it with separation between the components, and discussed the disc truncation. We find that in LSI+61303 the compact object comes into contact with the outer parts of the circumstellar disc at periastron, in MWC 148 the compact object goes deeply into the disc during the periastron passage, and in MWC 656 the black hole is accreting from the outer parts of the circumstellar disc along the entire orbit. The interstellar extinction was estimated using interstellar lines. The rotation of the mass donors appears to be similar to the rotation of the mass donors in Be/X-ray binaries. We suggest that X-ray/optical periodicity of about 1 day deserves to be searched for.

Context. 1E 1740.7-2942 is believed to be one of the two prototypical microquasars towards the Galactic center region whose X-ray states strongly resemble those of Cygnus X-1. Yet, the bipolar radio jets of 1E 1740.7-2942 are very reminiscent of a radio galaxy. The true nature of the object has thus remained an open question for nearly a quarter of a century. Aims. Our main goal here is to confirm the Galactic membership of 1E 1740.7-2942 by searching for morphological changes of its extended radio jets in human timescales. This work was triggered as a result of recent positive detection of fast structural changes in the large-scale jets of its 'twin' GRS 1758-258. Methods. We carried out an in-depth exploration of the VLA public archives and fully recalibrated all 1E 1740.7-2942 extended data sets in the C configuration of the array. We obtained and analyzed matching beam radio maps for five epochs, covering years 1992, 1993, 1994, 1997 and 2000, with an angular resolution of a few arcseconds. Results. We clearly detected structural changes in the arcminute jets of 1E 1740.7-2942 on timescales of roughly a year, which set a firm distance upper limit of 12 kpc. Moreover, a simple precessing model was simultaneously fitted to the five observing epochs available. The observed changes in the jet flow are strongly suggestive of a precession period of ~1.3 years. Conclusions. The fitting of the precession model to the data yields a distance of ~5 kpc. This value, and the observed changes, rule out any remaining doubts about the 1E 1740.7-2942 Galactic nature. To our knowledge, this microquasar is the second whose jet precession ephemeris become available after SS433. This kind of information is relevant to the physics of compact objects, since the genesis of the precession phenomenon occurs very close to the interplay region between the accretion disk and the compact object in the system.

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.

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.

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

We report 138 spectral observations of the H-alpha emission line of the radio- and gamma-ray-emitting Be/X-ray binary LSI+61303 obtained during the period of September 1998 -- January 2013. From measuring various H-alpha parameters, we found that the orbital modulation of the H-alpha is best visible in the equivalent width ratio EW(B)/EW(R), the equivalent width of the blue hump, and in the radial velocity of the central dip. The periodogram analysis confirmed that the H-alpha emission is modulated with the orbital and superorbital periods. For the past 20 years the radius of the circumstellar disk is similar to the Roche lobe size at the periastron. It is probably truncated by a 6:1 resonance. The orbital maximum of the equivalent width of H-alpha emission peaks after the periastron and coincides on average with the X-ray and gamma-ray maxima. All the spectra are available upon request from the authors and through the CDS.

RCW49 and its ionizing cluster form an extensive, complex region, widely studied at IR and optical wavelengths. Molonglo 843 MHz and ATCA data at 1.4 and 2.4 GHz showed two shells. Recent high-resolution IR imaging revealed a complex dust structure and ongoing star formation. New high-bandwidth and high-resolution data of the RCW49 field were obtained to survey the radio emission at arcsec scale and investigate the small-scale features and nature of the HII region.- Data were collected with the new 2-GHz bandwidth receivers and the CABB correlator of the ATCA, at 5.5 and 9.0 GHz. In addition, archival observations at 1.4 and 2.4 GHz have been re-reduced and re-analyzed, together with optical, IR, X-ray, and gamma-ray observations.- The CABB data result in the most detailed radio continuum images of RCW49 to date. The radio emission closely mimics the near-IR emission observed by Spitzer, showing pillars and filaments. The brightest continuum emission comes from the region known as the bridge. The overall flattish spectral index is consistent with a free-free emission mechanism. However, hints of nonthermal components are also present in the bridge. A jet-like structure surrounded by a bubble feature whose nature is still unclear has been discovered close to the Westerlund~2 core. Two apparent bow shocks and a number of discrete sources have been detected as well in the surroundings of RCW49. We also report on and discuss the possible detection of a H RRL.- The radio results support an association between the cm continuum and molecular emission. The detection of the RRL kinematically favors a RCW49 distance of 6--7 kpc. If the negative spectral indices measured at the bridge are caused by synchrotron emission, we propose a scenario where high-energy emission could be produced. Finally, the newly discovered jet-like structure appears to deserve a detailed study by itself.

The existence of multiple planetary systems involved in mean motion conmensurabilities has increased significantly since the Kepler mission. Although most correspond to 2-planet resonances, multiple resonances have also been found. The Laplace resonance is a particular case of a three-body resonance where the period ratio between consecutive pairs is n_1/n_2 near to n_2/n_3 near to 2/1. It is not clear how this triple resonance can act in order to stabilize (or not) the systems. The most reliable extrasolar system located in a Laplace resonance is GJ876 because it has two independent confirmations. However best-fit parameters were obtained without previous knowledge of resonance structure and no exploration of all the possible stable solutions for the system where done. In the present work we explored the different configurations allowed by the Laplace resonance in the GJ876 system by varying the planetary parameters of the third outer planet. We find that in this case the Laplace resonance is a stabilization mechanism in itself, defined by a tiny island of regular motion surrounded by (unstable) highly chaotic orbits. Low eccentric orbits and mutual inclinations from -20 to 20 degrees are compatible with the observations. A definite range of mass ratio must be assumed to maintain orbital stability. Finally we give constrains for argument of pericenters and mean anomalies in order to assure stability for this kind of systems.

The CORNISH project is the highest resolution radio continuum survey of the Galactic plane to date. It is the 5 GHz radio continuum part of a series of multi-wavelength surveys that focus on the northern GLIMPSE region (10 deg < l < 65 deg), observed by the Spitzer satellite in the mid-infrared. Observations with the Very Large Array in B and BnA configurations have yielded a 1.5" resolution Stokes I map with a root-mean-squared noise level better than 0.4 mJy/beam. Here we describe the data-processing methods and data characteristics, and present a new, uniform catalogue of compact radio-emission. This includes an implementation of automatic deconvolution that provides much more reliable imaging than standard CLEANing. A rigorous investigation of the noise characteristics and reliability of source detection has been carried out. We show that the survey is optimised to detect emission on size scales up to 14" and for unresolved sources the catalogue is more than 90 percent complete at a flux density of 3.9 mJy. We have detected 3,062 sources above a 7-sigma detection limit and present their ensemble properties. The catalogue is highly reliable away from regions containing poorly-sampled extended emission, which comprise less than two percent of the survey area. Imaging problems have been mitigated by down-weighting the shortest spacings and potential artefacts flagged via a rigorous manual inspection with reference to the Spitzer infrared data. We present images of the most common source types found: regions, planetary nebulae and radio-galaxies. The CORNISH data and catalogue are available online at http://cornish.leeds.ac.uk

In the last years several exoplanets have been discovered that orbit one component of a compact binary system (separation < 50 astronomical units), the probably best-known case is gamma-Cephei. So far, all attempts to explain the in-situ formation of these planets has been unsuccessful, in part because of the strong gravitational perturbations of the secondary star on any initial planetesimal swarm. Here we test whether planetary bodies in compact binaries, in particular gamma-Cephei, could have originated from a close encounter with a passing star, assuming initial configurations for the stellar system suitable for planetary formation. In other words, we analyze whether the orbital configuration of the current binary system might have been generated after the formation of the planet, and as a consequence of a close encounter with a third star in hyperbolic orbit. We performed a series of time-reverse N-body simulations of stellar scattering events in which the present-day configuration of gamma-Cephei was used as the initial condition plus a hypothetical third star as an impactor. We analyzed which configurations and system parameters could have given birth to the current system. Depending on the maximum impact velocity allowed for accretional collisions, we find that between 1% and 5% of stellar encounters correspond to an "original" system in which planetary formation around the primary star is not inhibited by the secondary, but is acceptable within the classical core-accretion scenario. Thus, although not highly probable, it is plausible that stellar encounters may have played a significant role in shaping these types of exoplanetary systems.

Synchrotron emission is commonly found in relativistic jets from active galactic nuclei (AGNs) and microquasars, but so far its presence in jets from young stellar objects (YSOs) has not been proved. Here, we present evidence of polarized synchrotron emission arising from the jet of a YSO. The apparent magnetic field, with strength of ~0.2 milligauss, is parallel to the jet axis, and the polarization degree increases towards the jet edges, as expected for a confining helical magnetic field configuration. These characteristics are similar to those found in AGN jets, hinting at a common origin of all astrophysical jets.

We present a study of the radio emission from a massive runaway star. The star forms a bow shock that is clearly observed in the infrared. We have performed VLA observations under the assumption that the reverse shock in the stellar wind might accelerate charged particles up to relativistic energies. Non-thermal radio emission of synchrotron origin has been detected, confirming the hypothesis. We have then modeled the system and we predict a spectral energy distribution that extends up to gamma-rays. Under some simplifying assumptions, we find that the intensity at high energies is too low to be detected by current instruments, but the future Cherenkov Telescope Array might detect the source.

Aims. We report the first detailed X-ray study of the supernova remnant (SNR) G304.6+0.1, achieved with the XMM-Newton mission. Methods. The powerful imaging capability of XMM-Newton was used to study the X-ray characteristics of the remnant at different energy ranges. The X-ray morphology and spectral properties were analyzed. In addittion, radio and mid-infrared data obtained with the Molonglo Observatory Synthesis Telescope and the Spitzer Space Telescope were used to study the association with the detected X-ray emission and to understand the structure of the SNR at differents wavelengths. Results. The SNR shows an extended and arc-like internal structure in the X-ray band with out a compact point-like source inside the remnant. We find a high column density of NH in the range 2.5-3.5x1022 cm-2, which supports a relatively distant location (d $\geq$ 9.7 kpc). The X-ray spectrum exhibits at least three emission lines, indicating that the X-ray emission has a thin thermal plasma origin, although a non-thermal contribution cannot be discarded. The spectra of three different regions (north, center and south) are well represented by a combination of a non-equilibrium ionization (PSHOCK) and a power-law (PL) model. The mid-infrared observations show a bright filamentary structure along the north-south direction coincident with the NW radio shell. This suggests that Kes 17 is propagating in a non-uniform environment with high density and that the shock front is interacting with several adjacent massive molecular clouds. The good correspondence of radio and mid-infrared emissions suggests that the filamentary features are caused by shock compression. The X-ray characteristics and well-known radio parameters indicate that G304.6+0.1 is a middle-aged SNR (2.8-6.4)x104 yr old and a new member of the recently proposed group of mixed-morphology SNRs.

Aims. We report results of an X-ray study of the supernova remnant (SNR) G344.7-0.1 and the point-like X-ray source located at the geometrical center of the SNR radio structure. Methods. The morphology and spectral properties of the remnant and the central X-ray point-like source were studied using data from the XMM-Newton and Chandra satellites. Archival radio data and infrared Spitzer observations at 8 and 24 $\mu$m were used to compare and study its multi-band properties at different wavelengths. Results. The XMM-Newton and Chandra observations reveal that the overall X-ray emission of G344.7-0.1 is extended and correlates very well with regions of bright radio and infrared emission. The X-ray spectrum is dominated by prominent atomic emission lines. These characteristics suggest that the X-ray emission originated in a thin thermal plasma, whose radiation is represented well by a plane-parallel shock plasma model (PSHOCK). Our study favors the scenario in which G344.7-0.1 is a 6 x 10^3 year old SNR expanding in a medium with a high density gradient and is most likely encountering a molecular cloud on the western side. In addition, we report the discovery of a soft point-like X-ray source located at the geometrical center of the radio SNR structure. The object presents some characteristics of the so-called compact central objects (CCO). However, its neutral hydrogen absorption column (N_{H}) is inconsistent with that of the SNR. Coincident with the position of the source, we found infrared and optical objects with typical early-K star characteristics. The X-ray source may be a foreground star or the CCO associated with the SNR. If this latter possibility were confirmed, the point-like source would be the farthest CCO detected so far and the eighth member of the new population of isolated and weakly magnetized neutron stars.

The environs of massive, early-type stars have been inspected in recent years in the search for sites where particles can be accelerated up to relativistic energies. Wind regions of massive binaries that collide have already been established as sources of high-energy emission; however, there is a different scenario for massive stars where strong shocks can also be produced: the bow-shaped region of matter piled up by the action of the stellar strong wind of a runaway star interacting with the interstellar medium. We study the bow-shock region produced by a very massive runaway star, BD+43 3654, to look for nonthermal radio emission as evidence of a relativistic particle population. We observed the field of BD+43 3654 at two frequencies, 1.42 and 4.86 GHz, with the Very Large Array (VLA), and obtained a spectral index map of the radio emission. We have detected, for the first time, nonthermal radio emission from the bow shock of a massive runaway star. After analyzing the radiative mechanisms that can be at work, we conclude that the region under study could produce enough relativistic particles whose radiation might be detectable by forthcoming gamma-ray instruments, like CTA North.

We present a catalog of cross-correlated radio, infrared and X-ray sources using a very restrictive selection criteria with an IDL-based code developed by us. The significance of the observed coincidences was evaluated through Monte Carlo simulations of synthetic sources following a well-tested protocol. We found 3320 coincident radio/X-ray sources with a high statistical significance characterized by the sum of error-weighted coordinate differences. For 997 of them, 2MASS counterparts were found. The percentage of chance coincidences is less than 1%. X-ray hardness ratios of well-known populations of objects were used to provide a crude representation of their X-ray spectrum and to make a preliminary diagnosis of the possible nature of unidentified X-ray sources. The results support the fact that the X-ray sky is largely dominated by Active Galactic Nuclei at high galactic latitudes (|b| >= 10^\circ). At low galactic latitudes (|b| <= 10^\circ) most of unidentified X-ray sources (~94%) lie at |b| <= 2^\circ. This result suggests that most of the unidentified sources found toward the Milky Way plane are galactic objects. Well-known and unidentified sources were classified in different tables with their corresponding radio/infrared and X-ray properties. These tables are intended as a useful tool for researchers interested in particular identifications.

The stable propagation of jets in FRII sources is remarkable if one takes into account that large-scale jets are subjected to potentially highly disruptive three-dimensional (3D) Kelvin-Helmholtz instabilities. Numerical simulations can address this problem and help clarify the causes of this remarkable stability. Following previous studies of the stability of relativistic flows in two dimensions (2D), it is our aim to test and extend the conclusions of such works to three dimensions. We present numerical simulations for the study of the stability properties of 3D, sheared, relativistic flows. This work uses a fully parallelized code Ratpenat that solves equations of relativistic hydrodynamics in 3D. The results of the present simulations confirm those in 2D. We conclude that the growth of resonant modes in sheared relativistic flows could be important in explaining the long-term collimation of extragalactic jets.

We obtain renormalized sets of right and left eigenvectors of the flux vector Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. The renormalization procedure relies on the characterization of the degeneracy types in terms of the normal and tangential components of the magnetic field to the wavefront in the fluid rest frame. Proper expressions of the renormalized eigenvectors in conserved variables are obtained through the corresponding matrix transformations. Our work completes previous analysis that present different sets of right eigenvectors for non-degenerate and degenerate states, and can be seen as a relativistic generalization of earlier work performed in classical MHD. Based on the full wave decomposition (FWD) provided by the the renormalized set of eigenvectors in conserved variables, we have also developed a linearized (Roe-type) Riemann solver. Extensive testing against one- and two-dimensional standard numerical problems allows us to conclude that our solver is very robust. When compared with a family of simpler solvers that avoid the knowledge of the full characteristic structure of the equations in the computation of the numerical fluxes, our solver turns out to be less diffusive than HLL and HLLC, and comparable in accuracy to the HLLD solver. The amount of operations needed by the FWD solver makes it less efficient computationally than those of the HLL family in one-dimensional problems. However its relative efficiency increases in multidimensional simulations.

Aims. We aim here to contribute to the identification of unassociated bright sources of gamma-rays in the recently released catalogue obtained by the Fermi collaboration. Methods. Our work is based on a extensive cross-identification of sources from different wavelength catalogues and databases. Results. As a first result, we report the finding of a few counterpart candidates inside the 95% confidence error box of the Fermi LAT unidentified gamma-ray source 0FGL J1848.6$-$0138. The globular cluster GLIMPSE-C01 remarkably stands out among the most peculiar objects consistent with the position uncertainty of the gamma-ray source and with a conceivable physical scenario for gamma-ray production. The Fermi observed spectrum is compared against theoretical predictions in the literature making the association plausible but not yet certain due to its low X-ray to gamma-ray luminosity ratio. Other competing counterparts are also discussed. In particular, we pay a special attention to a possible Pulsar Wind Nebula inside the Fermi error box whose nature is yet to be confirmed. Conclusions.Both a globular cluster and an infrared source resembling a Pulsar Wind Nebula have been found in positional agreement with 0FGL J1848.6$-$0138. In addition, other interesting objects in the field are also reported. Future gamma-ray observations will narrow the position uncertainty and we hope to eventually confirm one of the counterpart candidates reported here. If GLIMPSE-C01 is confirmed, together with the Fermi possible detection of the well known globular cluster 47 Tuc, then it would provide strong support to theoretical predictions of globular clusters as gamma-ray sources.

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

Context. A short duration burst reminiscent of a soft gamma-ray repeater/anomalous X-ray pulsar behaviour was detected in the direction of LS I +61 303 by the Swift satellite. While the association with this well known gamma-ray binary is likely, a different origin cannot be excluded. Aims. We explore the error box of this unexpected flaring event and establish the radio, near-infrared and X-ray sources in our search for any peculiar alternative counterpart. Methods. We carried out a combined analysis of archive Very Large Array radio data of LS I +61 303 sensitive to both compact and extended emission. We also reanalysed previous near infrared observations with the 3.5 m telescope of the Centro Astronomico Hispano Aleman and X-ray observations with the Chandra satellite. Results. Our deep radio maps of the LS I +61 303 environment represent a significant advancement on previous work and 16 compact radio sources in the LS I +61 303 vicinity are detected. For some detections, we also identify near infrared and X-ray counterparts. Extended emission features in the field are also detected and confirmed. The possible connection of some of these sources with the observed flaring event is considered. Based on these data, we are unable to claim a clear association between the Swift-BAT flare and any of the sources reported here. However, this study represents the most sophisticated attempt to determine possible alternative counterparts other than LS I +61 303.

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.

(Abridged) We present numerical simulations of the spectral evolution and emission of radio components in relativistic jets. We have developed an algorithm (SPEV) for the transport of a population of non-thermal electrons including radiative losses. For large values of the ratio of gas pressure to magnetic field energy density, \ab \sim 6\times 10^4, quiescent jet models show substantial spectral evolution, with observational consequences only above radio frequencies. Larger values of the magnetic field (\ab \sim 6\times 10^2), such that synchrotron losses are moderately important at radio frequencies, present a larger ratio of shocked-to-unshocked regions brightness than the models without radiative losses, despite the fact that they correspond to the same underlying hydrodynamic structure. We also show that jets with a positive photon spectral index result if the lower limit \gamma_min of the non-thermal particle energy distribution is large enough. A temporary increase of the Lorentz factor at the jet inlet produces a traveling perturbation that appears in the synthetic maps as a superluminal component. We show that trailing components can be originated not only in pressure matched jets, but also in over-pressured ones, where the existence of recollimation shocks does not allow for a direct identification of such features as Kelvin-Helmholtz modes, and its observational imprint depends on the observing frequency. If the magnetic field is large (\ab \sim 6\times 10^2), the spectral index in the rarefaction trailing the traveling perturbation does not change much with respect to the same model without any hydrodynamic perturbation. If the synchrotron losses are considered the spectral index displays a smaller value than in the corresponding region of the quiescent jet model.

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.

We report here on the first detection at X-ray wavelengths of the Supernova Remnant (SNR) G337.8-0.1, carried out with the XMM-Newton Observatory. Using the X-ray observations, we studied the X-ray morphology of the remnant at different energy ranges, analysed the spectral properties and investigated a possible variable behavior. The SNR shows a diffuse filled-center structure in the X-ray region with an absence of a compact source in its center. We find a high column density of N_H > 6.9 * 10^{22} cm^{-2}, which supports a relatively distant location (d > 7 kpc). The X-ray spectrum exhibits emission lines, indicating that the X-ray emission has a thin thermal plasma origin, and is well represented by a non-equilibrium ionization (NEI) plasma model. The X-ray characteristics and well-known radio parameters show that G337.8-0.1 belongs to the emerging class of mixed-morphology (MM) SNRs.

Context. A considerable fraction of the gamma-ray sources discovered with the Energetic Gamma-Ray Experiment Telescope (EGRET) remain unidentified. The EGRET sources that have been properly identified are either pulsars or variable sources at both radio and gamma-ray wavelengths. Most of the variable sources are strong radio blazars.However, some low galactic-latitude EGRET sources, with highly variable gamma-ray emission, lack any evident counterpart according to the radio data available until now. Aims. The primary goal of this paper is to identify and characterise the potential radio counterparts of four highly variable gamma-ray sources in the galactic plane through mapping the radio surroundings of the EGRET confidence contours and determining the variable radio sources in the field whenever possible. Methods. We have carried out a radio exploration of the fields of the selected EGRET sources using the Giant Metrewave Radio Telescope (GMRT) interferometer at 21 cm wavelength, with pointings being separated by months. Results. We detected a total of 151 radio sources. Among them, we identified a few radio sources whose flux density has apparently changed on timescales of months. Despite the limitations of our search, their possible variability makes these objects a top-priority target for multiwavelength studies of the potential counterparts of highly variable, unidentified gamma-ray sources.

We present numerical relativistic magnetohydrodynamic and emission simulations aimed to study the role played by the magnetic field in the dynamics and emission of relativistic jets in Active Galactic Nuclei. We focus our analysis on the study of the emission from recollimation shocks since they may provide an interpretation for the stationary components seen at parsec-scales in multiple sources. We show that the relative brightness of the knots associated with the recollimation shocks decreases with increasing jet magnetization, suggesting that jets presenting stationary components may have a relatively weak magnetization, with magnetic fields of the order of equipartition or below.

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.

Aims. We revisit the vicinity of the microquasar Cygnus X-3 at radio wavelengths. We aim to improve our previous search for possible associated extended radio features/hot spots in the position angle of the Cygnus X-3 relativistic jets focusing on shorter angular scales than previously explored. Methods. Our work is mostly based on analyzing modern survey and archive radio data, mainly including observations carried out with the Very Large Array and the Ryle Telescopes. We also used deep near-infrared images that we obtained in 2005. Results. We present new radio maps of the Cygnus X-3 field computed after combining multi-configuration Very Large Array archive data at 6 cm and different observing runs at 2 cm with the Ryle Telescope. These are probably among the deepest radio images of Cygnus X-3 reported to date at cm wavelengths. Both interferometers reveal an extended radio feature within a few arc-minutes of the microquasar position, thus making our detection more credible. Moreover, this extended emission is possibly non-thermal, although this point still needs confirmation. Its physical connection with the microquasar is tentatively considered under different physical scenarios. We also report on the serendipitous discovery of a likely Fanaroff-Riley type II radio galaxy only 3 arc-minute away from Cygnus X-3.

Microquasars are ideal natural laboratories for understanding accretion/ejection processes, studying the physics of relativistic jets, and testing gravitational phenomena. Nevertheless, these objects are difficult to find in our Galaxy. The main goal of this work is to increase the number of known systems of this kind. We have developed an improved search strategy based on positional cross-identification with very restrictive selection criteria to find new MQs, taking advantage of more sensitive modern radio and X-ray data. We find 86 sources with positional coincidence in the NVSS/XMM catalogs at |b|<10 deg. Among them, 24 are well-known objects and the remaining 62 sources are unidentified. For the fully coincident sources, whenever possible, we analyzed color-color and hardness ratio diagrams and found that at least 3 of them display high-mass X-ray binary characteristics, making them potential microquasar candidates.

The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. In this paper, we test the present theoretical and observational models via a relativistic numerical simulation of the jets in the radio galaxy 3C 31. We use the parameters derived from the modelling presented by \cite{lb02a,lb02b} as input parameters for the simulation of the evolution of the source, thus assuming that they have not varied over the lifetime of the source. We simulate about 10 % of the total lifetime of the jets in 3C 31. Realistic density and pressure gradients for the atmosphere are used. The simulation includes an equation of state for a two-component relativistic gas that allows a separate treatment of leptonic and baryonic matter. We compare our results with the modelling of the observational data of the source. Our results show that the bow shock evolves self-similarly at a quasi-constant speed, with slight deceleration by the end of the simulation, in agreement with recent X-ray observations that show the presence of bow shocks in FRI sources. The jet expands until it becomes underpressured with respect to the ambient medium, and then recollimates. Subsequent oscillations around pressure equilibrium and generation of standing shocks lead to the mass loading and disruption of the jet flow. We derive an estimate for the minimum age of the source of $t>1. 10^8 \rm{yrs}$, which may imply continuous activity of 3C 31 since the triggering of its activity. The simulation shows that weak CSS sources may be the young counterparts of FRIs. We conclude that the observed properties of the jets in 3C 31 are basically recovered by the standing shock scenario.

We present optical photo-polarimetric observations with high temporal resolution of the blazar AO 0235+164. Our data, the first to test the photo-polarimetric behaviour of this object at very short time-scales, show significant micro-variability in total flux, colour index, linear polarization degree, and position angle. Strong inter-night variations are also detected for these parameters. Although no correlation between colour index and total flux was found, our data seem to support the general bluer-when-brighter trend already known for this object. The polarization degree, in turn, shows no correlation with total flux, but a clear trend in the sense that colour index is redder (the spectrum is softer) when the measured polarization is higher.

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.

Context. The modern generation of Cherenkov telescopes has revealed a new population of gamma-ray sources in the Galaxy. Some of them have been identified with previously known X-ray binary systems while other remain without clear counterparts a lower energies. Our initial goal here was reporting on extensive radio observations of the first extended and yet unidentified source, namely TeV J2032+4130. This object was originally detected by the HEGRA telescope in the direction of the Cygnus OB2 region and its nature has been a matter of debate during the latest years. Aims. We aim to pursue our radio exploration of the TeV J2032+4130 position that we initiated in a previous paper but taking now into account the latest results from new Whipple and MILAGRO TeV telescopes. Methods. Our investigation is mostly based on interferometric radio observations with the Giant Metre Wave Radio Telescope (GMRT) close to Pune (India) and the Very Large Array (VLA) in New Mexico (USA). We also conducted near infrared observations with the 3.5 m telescope and the OMEGA2000 camera at the Centro Astronomico Hispano Aleman (CAHA) in Almeria (Spain). Results. We present deep radio maps centered on the TeV J2032+4130 position at different wavelengths. In particular, our 49 and 20 cm maps cover a field of view larger than half a degree that fully includes the Whipple position and the peak of MILAGRO emission. Our most important result here is a catalogue of 153 radio sources detected at 49 cm within the GMRT antennae primary beam with a full width half maximum (FWHM) of 43 arc-minute. Moreover, our multi-configuration VLA images reveal the non-thermal extended emission previously reported by us with improved angular resolution.

Qualitatively new aspects of the (linear and non-linear) stability of sheared relativistic (slab) jets are analyzed. The linear problem has been solved for a wide range of jet models well inside the ultrarelativistic domain (flow Lorentz factors up to 20; specific internal energies $\approx 60c^2$). As a distinct feature of our work, we have combined the analytical linear approach with high-resolution relativistic hydrodynamical simulations, which has allowed us i) to identify, in the linear regime, resonant modes specific to the relativistic shear layer ii) to confirm the result of the linear analysis with numerical simulations and, iii) more interestingly, to follow the instability development through the non-linear regime. We find that very high-order reflection modes with dominant growth rates can modify the global, long-term stability of the relativistic flow. We discuss the dependence of these resonant modes on the jet flow Lorentz factor and specific internal energy, and on the shear layer thickness. The results could have potential applications in the field of extragalactic relativistic jets.

In this study, we test the theoretically predicted mass-loss behaviour as a function of stellar effective temperature across the so-called `bi-stability' jump (BSJ). We gathered radio observations of 30 OB supergiants (O8-B3). We derived the radio mass-loss rates and wind efficiencies, and compared our results with Halpha mass-loss rates and predictions based on radiation-driven wind models. he wind efficiency shows the possible presence of a local maximum around an effective temperature of 21~000 K -- in qualitative agreement with predictions. We also find that the radio mass-loss rates show good agreement with empirical Halpha rates. However, the empirical mass-loss rates are larger than the predicted rates from radiation-driven wind theory for objects above the BSJ temperature, whilst they are smaller for the rest. A new wind momenta-luminosity relation for O8-B0 stars has been derived.

TeV J2032+4130 is the first extended very high energy gamma-ray source, which has remained enigmatic since its discovery, due to the lack of identification. We report here deep radio observations covering the TeV J2032+4130 field and revealing for the first time an extended and diffuse radio emission, as well as a remarkable population of compact radio sources. Some of these radio sources are in positional coincidence with X-ray and optical/IR sources. Future follow up studies of these new radio sources will likely contribute to solve the mystery of this extended unidentified TeV source.

We aim to investigate the nature of the high energy source KS 1741-293 by revisiting the radio and infrared associations proposed in the early 1990s. Our work is mostly based on the analysis of modern survey and archive data, including the NRAO, MSX, 2MASS and Chandra archives, and catalogues. We also have obtained deep CCD optical observations by ourselves. The coincidence of KS 1741-293 with an extended radio and far-infrared source, tentatively suggested in 1994, is no longer supported by modern observational data. Instead, a Chandra source is the only peculiar object found to be consistent with all high-energy error circles of KS 1741-293 and we propose it to be its most likely X-ray counterpart. We also report the existence of a non-thermal radio nebula in the vicinity of the KS 1741-293 position with the appearance of a supernova remnant. The possibility of being associated to this X-ray binary is discussed.

We present the results of a first campaign of radio continuum observations toward supergiants with spectral types in the range O8 - B3. Three targets out of twelve were detected with the ATCA and the VLA. The measured fluxes and the upper limits yielded values of the stellar mass-loss rates dM/dt and the wind efficiencies. The comparison of predicted and derived values of dM/dt shows a good agreement.

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.

Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed in the parsec-scale jet of the well-known quasar 3C 273. We present here the results from a 3D relativistic hydrodynamics numerical simulation based on the parameters given for the jet by Lobanov & Zensus (2001), and one in which the effects of jet precession and the injection of discrete components have been taken into account. We compare the model with the structures observed in 3C 273 using very long baseline interferometry and constrain the basic properties of the flow. We find growing perturbation modes in the simulation with similar wavelengths to those observed, but with a different set of wave speeds and mode identification. If the observed longest helical structure is produced by the precession of the flow, longer precession periods should be expected. Our results show that some of the observed structures could be explained by growing Kelvin-Helmholtz instabilities in a slow moving region of the jet. However, we point towards possible errors in the mode identification that show the need of more complete linear analysis in order to interpret the observations. We conclude that, with the given viewing angle, superluminal components and jet precession cannot explain the observed structures.

We search for ultra luminous X-ray source (ULXs) radio counterparts located in nearby galaxies in order to constrain their physical nature. Our work is based on a systematic cross-identification of the most recent and extensive available ULX catalogues and archival radio data. A catalogue of 70 positional coincidences is reported. Most of them are located within the galaxy nucleus. Among them, we find 11 new cases of non-nuclear ULX sources with possibly associated radio emission.

We report on the design and status of the FLAMINGOS-2 instrument - a fully-cryogenic facility near-infrared imager and multi-object spectrograph for the Gemini 8-meter telescopes. FLAMINGOS-2 has a refractive all-spherical optical system providing 0.18-arcsecond pixels and a 6.2-arcminute circular field-of-view on a 2048x2048-pixel HAWAII-2 0.9-2.4 mm detector array. A slit/decker wheel mechanism allows the selection of up to 9 multi-object laser-machined plates or 3 long slits for spectroscopy over a 6x2-arcminute field of view, and selectable grisms provide resolutions from $\sim$ 1300 to $\sim $3000 over the entire spectrograph bandpass. FLAMINGOS-2 is also compatible with the Gemini Multi-Conjugate Adaptive Optics system, providing multi-object spectroscopic capabilities over a 3x1-arcminute field with high spatial resolution (0.09-arcsec/pixel). We review the designs of optical, mechanical, electronics, software, and On-Instrument WaveFront Sensor subsystems. We also present the current status of the project, currently in final testing in mid-2006.

Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed in the parsec-scale jet of the well known quasar 3C 273. We compare the model with the radio structure observed in 3C 273 on parsec scales using very long baseline interferometry (VLBI) and constrain the basic properties of the flow. Our results show that Kelvin-Helmholtz instabilities are the most plausible mechanism to generate the observed structures.

To confirm, or rule out, the possible hot spot nature of two previously detected radio sources in the vicinity of the Cygnus X-3 microquasar. We present the results of a radio and near infrared exploration of the several arc-minute field around the well known galactic relativistic jet source Cygnus X-3 using the Very Large Array and the Calar Alto 3.5~m telescope. The data this paper is based on do not presently support the hot spot hypothesis. Instead, our new observations suggest that these sources are most likely background or foreground objects. Actually, none of them appears to be even barely extended as would be expected if they were part of a bow shock structure. Our near infrared observations also include a search for extended emission in the Bracket $\gamma$ (2.166 $\mu$m) and $H_{2}$ (2.122 $\mu$m) lines as possible tracers of shocked gas in the Cygnus X-3 surroundings. The results were similarly negative and the corresponding upper limits are reported.

The linear and non-linear stability of sheared, relativistic planar jets is studied by means of linear stability analysis and numerical hydrodynamical simulations. Our results extend the previous Kelvin-Hemlholtz stability studies for relativistic, planar jets in the vortex sheet approximation performed by Perucho et al. (2004a,b) by including a shear layer between the jet and the external medium and more general perturbations. The models considered span a wide range of Lorentz factors ($2.5-20$) and internal energies ($0.08 c^2-60 c^2$) and are classified into three classes according to the main characteristics of their long-term, non-linear evolution. We observe a clear separation of these three groups in a relativistic Mach-number Lorentz-factor plane. Jets with a low Lorentz factor and small relativistic Mach number are disrupted after saturation. Those with a large Lorentz factor and large relativistic Mach number are the stablest, due to the appearance of short wavelength resonant modes which generate local mixing and heating in the shear layer around a fast, unmixed core, giving a plausible solution for the problem of the long-term stability of relativistic jets. A third group is present between them, including jets with intermediate values of Lorentz factor and relativistic Mach number, which are disrupted by a slow process of mixing favored by an efficient and continuous conversion of kinetic into internal energy. In the long term, all the models develop a distinct transversal structure (shear/transition layers) as a consequence of KH perturbation growth, depending on the class they belong to. The properties of these shear layers are analyzed in connection with the parameters of the original jet models.

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.

We have extended the procedure to find the exact solution of the Riemann problem in relativistic hydrodynamics to a particular case of relativistic magnetohydrodynamics in which the magnetic field of the initial states is tangential to the discontinuity and orthogonal to the flow velocity. The wave pattern produced after the break up of the initial discontinuity is analogous to the non--magnetic case and we show that the problem can be understood as a purely relativistic hydrodynamical problem with a modified equation of state. The new degree of freedom introduced by the non-zero component of the magnetic field results in interesting effects consisting in the change of the wave patterns for given initial thermodynamical states, in a similar way to the effects arising from the introduction of tangential velocities. Secondly, when the magnetic field dominates the thermodynamical pressure and energy, the wave speeds approach the speed of light leading to fast shocks and fast and arbitrarily thin rarefaction waves. Our approach is the first non-trivial exact solution of a Riemann problem in relativistic magnetohydrodynamics and it can also be of great interest to test numerical codes against known analytical or exact solutions.

We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbolic character in the solution procedure. All theoretical ingredients necessary to build up high-resolution shock-capturing schemes based on the solution of local Riemann problems (i.e. Godunov-type schemes) are described. In particular, we use a renormalized set of regular eigenvectors of the flux Jacobians of the relativistic magnetohydrodynamics equations. In addition, the paper describes a procedure based on the equivalence principle of general relativity that allows the use of Riemann solvers designed for special relativistic magnetohydrodynamics in GRMHD. Our formulation and numerical methodology are assessed by performing various test simulations recently considered by different authors. These include magnetized shock tubes, spherical accretion onto a Schwarzschild black hole, equatorial accretion onto a Kerr black hole, and magnetized thick accretion disks around a black hole prone to the magnetorotational instability.

We present the results of a deep search for associated radio features in the vicinity of the microquasar Cygnus X-3. The motivation behind is to find out evidence for interaction between its relativistic jets and the surrounding interstellar medium, which could eventually allow us to perform calorimetry of the total energy released by this microquasar during its flaring lifetime. Remarkably, two radio sources with mJy emission level at centimeter wavelengths have been detected in excellent alignment with the position angle of the inner radio jets. We propose that these objects could be the hot spots where the relativitic outflow collides with the ambient gas in analogy with Fanaroff-Riley II radio galaxies. These candidate hot spots are within a few arc-minutes of Cygnus X-3 and, if physically related, the full linear extent of the jet would reach tens of parsecs. We discuss here the evidence currently available to support this hypothesis based on both archival data and our own observations.

New optical spectroscopy of the HMXB microquasar LSI +61 303 is presented. Eccentric orbital fits to our radial velocity measurements yield updated orbital parameters in good agreement with previous work. Our orbital solution indicates that the periastron passage occurs at radio phase 0.23 and the X-ray/radio outbursts are triggered 2.5--4 days after the compact star passage. The spectrum of the optical star is consistent with a B0 V spectral type and contributes ~65 percent of the total light, the remaining due to emission by a circumstellar disc. We also measure the projected rotational velocity to be v sin i = 113 km/s.

Intensive VLBI monitoring programs of jets in AGN are showing the existence of intricate emission patterns, such as upstream motions or slow moving and quasi-stationary componentes trailing superluminal features. Relativistic hydrodynamic and emission simulations of jets are in very good agreement with these observations, proving as a powerful tool for the understanding of the physical processes taking place in the jets of AGN, microquasars and GRBs. These simulations show that the variability of the jet emission is the result of a complex combination of phase motions, viewing angle selection effects, and non-linear interactions between perturbations and the underlying jet and/or ambient medium. Both observations and simulations suggest that shock-in-jet models may be an overly simplistic idealization when interpreting the emission structure observed in actual jets.

We present a hydrodynamical 3D simulation of the relativistic jet in 3C273, in comparison to previous linear perturbation analysis of Kelvin-Helmholtz instability developing in the jet. Our aim is to assess advantages and limitations of both analytical and numerical approaches and to identify spatial and temporal scales on which the linear regime of Kelvin-Helmholtz instability can be applied in studies of morphology and kinematics of parsec-scale jets.

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.

Compact Symmetric Objects (CSOs) are considered the young counterparts of large doubles according to advance speeds measured or inferred from spectral ageing. Here we present a simple power law model for the CSO/FRII evolution based on the study of sources with well defined hot-spots. The luminosity of the hot spots is estimated under minimum energy conditions. The advance of the source is considered to proceed in ram pressure equilibrium with the ambient medium. Finally, we also assume that the jets feeding the hot spots are relativistic and have a time dependent power. Comparison with observational data allows to interpret the CSO-FRII evolution in terms of decreasing jet power with time.

In this paper we continue our study of the Kelvin-Helmholtz (KH) instability in relativistic planar jets following the long-term evolution of the numerical simulations which were introduced in Paper I. The models have been classified into four classes (I to IV) with regard to their evolution in the nonlinear phase, characterized by the process of jet/ambient mixing and momentum transfer. Models undergoing qualitatively different non-linear evolution are clearly grouped in well-separated regions in a jet Lorentz factor/jet-to-ambient enthalpy diagram. Jets with a low Lorentz factor and small enthalpy ratio are disrupted by a strong shock after saturation. Those with a large Lorentz factor and enthalpy ratio are unstable although the process of mixing and momentum exchange proceeds to a longer time scale due to a steady conversion of kinetic to internal energy in the jet. In these cases, the high value of the initial Lorentz seems to prevent transversal velocity from growing far enough to generate the strong shock that breaks the slower jets. Finally, jets with either high Lorentz factors and small enthalpy ratios or low Lorentz factors and large enthalpy ratios appear as the most stable.

The effects of relativistic dynamics and thermodynamics in the development of Kelvin-Helmholtz instabilities in planar, relativistic jets along the early phases (namely linear and saturation phases) of evolution has been studied by a combination of linear stability analysis and high-resolution numerical simulations for the most unstable first reflection modes in the temporal approach. Three different values of the jet Lorentz factor (5, 10 and 20) and a few different values of specific internal energy of the jet matter (from 0.08 to $60.0 c^2$) have been considered. Figures illustrating the evolution of the perturbations are also shown.

We assess the suitability of a recent high-resolution central scheme developed by Kurganov & Tadmor (2000) for the solution of the relativistic hydrodynamics equations. The novelty of this approach relies on the absence of Riemann solvers in the solution procedure. The computations we present are performed in one and two spatial dimensions in Minkowski spacetime. Standard numerical experiments such as shock tubes and the relativistic flat-faced step test are performed. As an astrophysical application the article includes two-dimensional simulations of the propagation of relativistic jets using both Cartesian and cylindrical coordinates. The simulations reported clearly show the capabilities of the numerical scheme to yield satisfactory results, with an accuracy comparable to that obtained by the so-called high-resolution shock-capturing schemes based upon Riemann solvers (Godunov-type schemes), even well inside the ultrarelativistic regime. Such central scheme can be straightforwardly applied to hyperbolic systems of conservation laws for which the characteristic structure is not explicitly known, or in cases where the exact solution of the Riemann problem is prohibitively expensive to compute numerically. Finally, we present comparisons with results obtained using various Godunov-type schemes as well as with those obtained using other high-resolution central schemes which have recently been reported in the literature.

We present preliminary results of a campaign undertaken with different radio interferometers to observe a sample of the most variable unidentified EGRET sources. We expect to detect which of the possible counterparts of the gamma-ray sources (any of the radio emitters in the field) varies in time with similar timescales as the gamma-ray variation. If the gamma-rays are produced in a jet-like source, as we have modelled theoretically, synchrotron emission is also expected at radio wavelengths. Such radio emission should appear variable in time and correlated with the gamma-ray variability.

We present electro-photometric UBV and high-speed U-band flickering observations of the recurrent nova T CrB during a period when its U brightness varies by more than 2 mag. The V band is dominated by the ellipsoidal variability of the red giant, however, the variability of the hot component also causes 0.15 mag variations in V. We define a set of parameters which characterise the flickering. The Fourier spectra of all 27 nights are similar to each other. The power spectral density of the variations has a power law component (f^{-1.46} on average). We do not detect a dependence of the Fourier spectra and autocorrelation function on the brightness of the object. Having subtracted the contribution of the red giant, we show that the flickering amplitude correlates with the average flux of the accreting component. A comparison with CH Cyg and MWC 560 indicates that the flickering of T CrB is more stable (at least during the time of our observations), than that in the jet-ejecting symbiotic stars. The data are available in electronic form from the authors.

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.

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.

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.

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 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 photometric and spectroscopic optical and near-infrared (NIR) observations (Based on observations collected at the European Southern Observatory, Chile (ESO ID 63.H-0493 and 64.H-0382)) taken during the outburst of the microquasar V4641 Sgr = SAX J1819.3-2525 (in'tzand et al., 2000) in September 1999. We observed an increase in the J-Ks colour between 5 and 8 days after the outburst, which we interpret as likely evidence for the presence of dust around the source. We also observed an extraordinarily strong, broad and variable H_alpha line, with a velocity width of 4560 km/s suggesting the presence of a high-velocity outflow component. We constrain the distance of the system between 3 and 8 kpc, locating it further away than previously derived from radio observations (Hjellming et al., 2000), but consistent with Orosz et al. (2001). We then discuss the nature of this system, showing that the companion star is either a B3-A2 main sequence star, or a B3-A2 sub-giant crossing the Hertzsprung gap. The system is therefore an Intermediate or High Mass X-ray Binary System (IMXB or HMXB). The distance derived by these optical/NIR observations implies that the jets observed by Hjellming et al. (2000) would then exhibit apparent velocities of ~ 10 c. We finally discuss the possibility of an interaction between the jets and surroundings of the source, and also of this source being a ``microblazar''.

We present the results of a three-dimensional, relativistic, hydrodynamic simulation of a precessing jet into which a compact blob of matter is injected. A comparison of synthetic radio maps computed from the hydrodynamic model, taking into account the appropriate light travel time delays, with those obtained from observations of actual superluminal sources shows that the variability of the jet emission is the result of a complex combination of phase motions, viewing angle selection effects, and non-linear interactions between perturbations and the underlying jet and/or the external medium. These results question the hydrodynamic properties inferred from observed apparent motions and radio structures, and reveal that shock-in-jet models may be overly simplistic.

A set of simulations concerning the influence of internal energy on the stability of relativistic jets is presented. Results show that perturbations saturate when the amplitude of the velocity perturbation approaches the speed of light limit. Also, contrary to what predicted by linear stability theory, jets with higher specific internal energy appear to be more stable.

We study the influence of the matter content of extragalactic jets on their morphology, dynamics and emission properties. For this purpose we consider jets of extremely different compositions including pure leptonic and baryonic plasmas. Our work is based on two-dimensional relativistic hydrodynamic simulations of the long-term evolution of powerful extragalactic jets propagating into a homogeneous environment. The equation of state used in the simulations accounts for an arbitrary mixture of electrons, protons and electron-positron pairs. Using the hydrodynamic models we have also computed synthetic radio maps and the thermal Bremsstrahlung X-ray emission from their cavities. Although there is a difference of about three orders of magnitude in the temperatures of the cavities inflated by the simulated jets, we find that both the morphology and the dynamic behaviour are almost independent on the assumed composition of the jets. Their evolution proceeds in two distinct epochs. During the first one multidimensional effects are unimportant and the jets propagate ballistically. The second epoch starts when the first larger vortices are produced near the jet head causing the beam cross section to increase and the jet to decelerate. The evolution of the cocoon and cavity is in agreement with a simple theoretical model. The beam velocities are relativistic ($\Gamma \simeq 4$) at kiloparsec scales supporting the idea that the X-ray emission of several extragalactic jets may be due to relativistically boosted CMB photons. The radio emission of all models is dominated by the contribution of the hot spots. All models exhibit a depression in the X-rays surface brightness of the cavity interior in agreement with recent observations.

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.

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 report on several European Southern Observatory (ESO) near-infrared (NIR) observational campaigns aimed at understanding the nature of galactic black hole candidates. Our results, including NIR photometry of the sources GRO J1655-40, GRS 1739-278, GRS 1716-249, GRS 1121-68 and GX 339-4, show that all the sources but GRO J1655-40 are consistent with low-mass stars as the companion star of the binary system. By locating the counterparts on a colour-magnitude diagram (CMD), we better constrain the spectral type of the companion star of three of the systems considered here, and confirm a fourth one. The spectral types are respectively: M0-5 V for GRS 1716-249, F8-G2 III for GX 339-4 and later than F5 V for GRS 1739-278. We confirm the already known spectral type of the companion in GRS 1121-68 (K0-5 V). The location of GRO J1655-40 on the CMD is consistent with the sub-giant luminosity class and with this source crossing the Hertzsprung gap. However, a non-stellar emission seems to contribute to the NIR flux of this source.

We present a model to determine the physical parameters of jets and hot spots of a sample of CSOs under very basic assumptions like synchrotron emission and minimum energy conditions. Based on this model we propose a simple evolutionary scenario for these sources assuming that they evolve in ram pressure equilibrium with the external medium and constant jet power. The parameters of our model are constrained from fits of observational data (radio luminosity, hot spot radius and hot spot advance speed) versus projected linear size. From these plots we conclude that CSOs evolve self-similarly and that their radio luminosity increases with linear size along the first kiloparsec. Assuming that the jets feeding CSOs are relativistic from both kinematical and thermodynamical points of view, we use the values of the pressure and particle number density within the hot spots to estimate the fluxes of momentum (thrust), energy, and particles of these relativistic jets. The mean jet power obtained in this way is within an order of magnitude that inferred for FRII sources, which is consistent with CSOs being the possible precursors of large doubles. The inferred flux of particles corresponds to, for a barionic jet, about a 10% of the mass accreted by a black hole of $10^8 {\rm M_{\odot}}$ at the Eddington limit, pointing towards a very efficient conversion of accretion flow into ejection, or to a leptonic composition of jets.

We tested the ejector-propeller model of the Be/X-ray binary LSI+61303 (V 615 Cas, GT 0236+620) by using the parameters predicted by the model in the calculations of the X-ray and radio variability. The results are: (1) in terms of the Ejector-Propeller model, the X-ray maximum is due to the periastron passage; (2) the radio outburst can be really a result of the transition from the propeller to ejector regimes; (3) the radio outburst will delay with respect to the X-ray maximum every orbital period. The proposed scenario seems to be in good agreement with the observations.

Observations of the extraordinarily bright optical afterglow (OA) of GRB 991208 started 2.1 d after the event. The flux decay constant of the OA in the R-band is -2.30 +/- 0.07 up to 5 d, which is very likely due to the jet effect, and after that it is followed by a much steeper decay with constant -3.2 +/- 0.2, the fastest one ever seen in a GRB OA. A negative detection in several all-sky films taken simultaneously to the event implies either a previous additional break prior to 2 d after the occurrence of the GRB (as expected from the jet effect). The existence of a second break might indicate a steepening in the electron spectrum or the superposition of two events. Once the afterglow emission vanished, contribution of a bright underlying SN is found, but the light curve is not sufficiently well sampled to rule out a dust echo explanation. Our determination of z = 0.706 indicates that GRB 991208 is at 3.7 Gpc, implying an isotropic energy release of 1.15 x 10E53 erg which may be relaxed by beaming by a factor > 100. Precise astrometry indicates that the GRB coincides within 0.2" with the host galaxy, thus given support to a massive star origin. The absolute magnitude is M_B = -18.2, well below the knee of the galaxy luminosity function and we derive a star-forming rate of 11.5 +/- 7.1 Mo/yr. The quasi-simultaneous broad-band photometric spectral energy distribution of the afterglow is determined 3.5 day after the burst (Dec 12.0) implying a cooling frequency below the optical band, i.e. supporting a jet model with p = -2.30 as the index of the power-law electron distribution.

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 photometric and spectroscopic optical and near-infrared (NIR) observations (Based on observations collected at the European Southern Observatory, Chile (ESO ID 63.H-0493 and 64.H-0382)) taken during the outburst of the recently flaring source V4641 Sgr = SAX J1819.3-2525 (in't Zand et al., 2000), on September 1999. This source was independently detected as the RXTE source XTE J1819-254 (Markwardt et al. 1999), and afterwards identified with the variable object V4641 Sgr (Kato et al. 1999). It underwent a bright optical outburst on 1999 Sept. 15.7 UT, going from magnitude 14 to 8.8 in the V-band (Stubbings 1999), and reaching 12.2 Crab in the X-rays (Smith et al. 1999) and Ks nearly 13. This outburst was therefore bright, but very brief, with an e-fold decay time of 0.6 days. A radio source was resolved, making of V4641 Sgr a new microquasar (Hjellming et al. 2000). We discuss the nature of this system, showing that our observations suggest a distance farther than previously derived from the radio observations (Hjellming et al. 2000). The distance of the system would be between 3 and 8 kpc, the companion star being a B3-A2 main sequence star. Another possibility is that the companion star is crossing the Hertzsprung gap (type B3-A2 IV), and in this case the distance cited above would be the minimum distance of the system. The system is therefore an Intermediate or High Mass X-ray Binary System (IMXB or HMXB). The inconsistency regarding the distance between the radio and optical/NIR observations could be explained by the detection of an interaction between matter ejected before the X-ray outburst and the surrounding medium of the source. If this is confirmed, this source could be added to the short list of microquasars where such an interaction has been detected.

We present a new centimeter polarimetric VLBI image of the BL Lac object 0735+178. This source exhibits one of the most pronounced curvatures observed in jets of AGNs, with two sharp apparent bends of 90 degrees within the inner 2 milliarcseconds from the core. Through the analysis of the data gathered over the past decades we study whether this curvature is produced by precession of the jet with ejection of ballistic components, or a precession in such a way that components' velocity vector are always parallel to the jet axis. These possibilities are also studied by comparison with 3D hydrodynamic relativistic simulations of precessing jets.

We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component, instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness, and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component.

We present a comparative study of the circumstellar disks in Be/X-ray binaries and isolated Be stars based upon the H-alpha emission line. From this comparison it follows that the overall structure of the disks in the Be/X-ray binaries is similar to the disks of other Be stars, i.e. they are axisymmetric and rotationally supported. The factors for the line broadening (rotation and temperature) in the disks of the Be stars and the Be/X-ray binaries seem to be identical. However, we do detect some intriguing differences between the envelopes. On average, the circumstellar disks of the Be/X-ray binaries are twice as dense as the disks of the isolated Be stars. The different distribution of the Be/X-ray binaries and the Be stars seen in the full width at half maximum (FWHM) versus peak separation diagram indicates that the disks in Be/X-ray binaries have on average a smaller size, probably truncated by the compact object.

Compact symmetric objects conform a class of sources characterized by high luminosity radio emission located symmetrically on both sides of the active galactic nucleus on linear scales of less than 1 kpc. Given their small size, the hot spots of the jets in CSOs provide a unique laboratory for the study of the physics of relativistic jets and their environment close to the central engine. We present a simple model for the hot spots in CSOs assuming synchrotron emission, minimum energy and ram pressure equilibrium with the external medium. Further comparison of our model with observational data allows us to constrain the physical parameters in the hot spots and the jets feeding them, and the density profile of the external medium.

We present photometric and spectroscopic optical and near-infrared (NIR) observations taken during the outburst on September 1999 of the source V4641 Sgr = SAX J1819.3-2525 (in'tzand et al., 2000) = XTE J1819-254 (Markwardt et al., 1999). We show that our observations suggest a distance between 4 and 8 kpc, the spectral type of the companion star being constrained between B3 and A2 V, making the system a High Mass X-ray Binary System (HMXB). In view of the radio and optical/NIR observations, it is possible that the ejecta of the source interacted with the surrounding medium of the source.

We report the results of our search for the near infrared counterpart of the microquasar 1E 1740.7-2942 using the VLT (Based on observations collected at the European Southern Observatory, Chile (ESO No 63.H-0261)). For the first time, several counterpart candidates have been found in our Ks-band images that may be consistent with the best radio and X-ray positions available for 1E 1740.7-2942. However, the non-detection of variability between two observing epochs and the positional uncertainty still remaining at the nearly less than 1 arcsec level prevent us from identifying an unambiguous counterpart. Alternatively, the VLT images set new upper limits significantly deeper than previously reported that constrain the binary companion to be later than B8 V or earlier than G5 III.

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.

We report our attempts to detect and confirm a narrow moving component in the H-alpha emission line of the radio emitting Be/X-ray binary LSI+61 303 (V615 Cas, GT0236+610). The existence of this spectral feature was already suspected in the past. As a result, we find that this component does exist and that its radial velocity varies in agreement with the radio period of the system. We interpret it tentatively as due to a denser region, or bulge, orbiting near the outer edge of the H-alpha emitting disk.

A remarkable correlation between the H-alpha emission line and the radio behaviour of LSI+61 303 (V615 Cas, GT 0236+610) over its 4 yr modulation is discovered. The radio outburst peak is shifted by a quarter of the 4 yr modulation period (about 400 days) with respect to the equivalent width of the H-alpha emission line variability. The onset of the LSI+61 303 radio outbursts varies in phase with the changes of the H-alpha emission line, at least during the increase of H-alpha equivalent width. This is the first clear correlation between the emission associated to the compact object and the Be circumstellar disk in a Be/X-ray binary system.

We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics for arbitrary tangential flow velocities. The solution is obtained by solving the jump conditions across shocks plus an ordinary differential equation arising from the self-similarity condition along rarefaction waves, in a similar way as in purely normal flow. The dependence of the solution on the tangential velocities is analysed, and the impact of this result on the development of multidimensional relativistic hydrodynamic codes (of Godunov type) is discussed.

We present near infrared spectroscopy of the superluminal microquasar \grs\ obtained with the first unit of the VLT (Based on observations collected at the European Southern Observatory, Chile (ESO No 63.H-0261).) and the ISAAC spectro-imager. The emission features detected in the VLT data have been identified as He I, Br $\gamma$, He II and Na I. The detection of Na I is reported here for the first time, while our confirmation of weak He II emission provides support to previous marginal detections of this feature. By comparing the observed spectra with those of massive stars, we find that our results are very consistent with \grs being a high mass X-ray system with an early type primary, as previously proposed by Chaty et al. (1996) and Mirabel et al. (1997). The VLT spectra also provide evidence of P Cygni profiles, that turn into blue emission wings when the system is in outburst. This observed line profile evolution implies that \grs must be surrounded by an expanding envelope, that is partially blown out during the X-ray outbursts. The presence of such circumstellar gaseous material around \grs is more naturally understood in the context of a massive luminous star than if the system was a low-mass X-ray binary.

Using a collapsar progenitor model of MacFadyen & Woosley we have simulated the propagation of an axisymmetric jet through a collapsing rotating massive star with the GENESIS multi-dimensional relativistic hydrodynamic code. The jet forms as a consequence of an assumed (constant or variable) energy deposition in the range $10^{50}$ erg s$^{-1}$ to $10^{51}$ erg s$^{-1}$ within a $30^{\circ}$ cone around the rotation axis. The jet flow is strongly beamed ($\la$ few degrees), spatially inhomogeneous, and time dependent. The jet reaches the surface of the stellar progenitor ($R_{\ast} = 2.98 \times 10^{10} $cm) intact. At breakout the maximum Lorentz factor of the jet flow is 33. After breakout the jet accelerates into the circumstellar medium, whose density is assumed to decrease exponentially and then being constant $\rho_{\rm ext} = 10^{-5}$ gcm$^{-3}$. Outside the star the flow begins to expand also laterally ($v \sim c$), but the beam remains very well collimated. At a distance of $2.54 R_{\ast}$, where the simulation ends, the Lorentz factor has increased to 44.

We present the first radio emission simulations from high resolution three dimensional relativistic hydrodynamic jets, which allow for a study of the observational implications of the interaction between the jet and external medium. This interaction gives rise to a stratification of the jet where a fast spine is surrounded by a slow high energy shear layer. The stratification, and in particular the large specific internal energy and slow flow in the shear layer largely determines the emission from the jet. If the magnetic field in the shear layer becomes helical (e.g., resulting from an initial toroidal field and an aligned field component generated by shear) the emission shows a cross section asymmetry, in which either the top or the bottom of the jet dominates the emission. This, as well as limb or spine brightening, is a function of the viewing angle and flow velocity, and the top/bottom jet emission predominance can be reversed if the jet changes direction with respect to the observer, or presents a change in velocity. The asymmetry is more prominent in the polarized flux, because of field cancellation (or amplification) along the line of sight. Recent observations of jet cross section emission asymmetries in the blazar 1055+018 can be explained assuming the existence of a shear layer with a helical magnetic field.

Our contribution concerns with the numerical solution of the 3D general relativistic hydrodynamical system of equations within the framework of the 3+1 formalism. We summarize the theoretical ingredients which are necessary in order to build up a numerical scheme based on the solution of local Riemann problems. Hence, the full spectral decomposition of the Jacobian matrices of the system, i.e., the eigenvalues and the right and left eigenvectors, is explicitly shown. An alternative approach consists in using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows. Our proposal relies on a local change of coordinates in terms of which the spacetime metric is locally Minkowskian and permits an accurate description of numerical general relativistic hydrodynamics.

We have studied the relativistic beamed outflow proposed to occur in the collapsar model of gamma-ray bursts. A jet forms as a consequence of an assumed energy deposition of $\sim 10^{50}- 10^{51}$ erg/s within a $30^{\circ}$ cone around the rotation axis of the progenitor star. The generated jet flow is strongly beamed ($\la$ few degrees) and reaches the surface of the stellar progenitor (r $\approx 3 10^{10} $cm) intact. At break-out the maximum Lorentz factor of the jet flow is about 33. Simulations have been performed with the GENESIS multi-dimensional relativistic hydrodynamic code.

We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics for arbitrary tangential flow velocities. The solution is obtained by solving the jump conditions across shocks plus an ordinary differential equation arising from the self-similarity condition along rarefaction waves, in a similar way as in purely normal flow. The dependence of the solution on the tangential velocities is analysed. This solution has been used to build up an exact Riemann solver implemented in a multidimensional relativistic (Godunov-type) hydro-code.

The multidimensional relativistic hydrodynamical code GENESIS has been used to obtain first results of {\it 3D} simulations of relativistic jets. We have studied the influence of a slight perturbation of the injection velocity field on the morphodynamics of otherwise axisymmetric relativistic jets.

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.

This review is concerned with a discussion of numerical methods for the solution of the equations of special relativistic hydrodynamics (SRHD). Particular emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods in SRHD. Results obtained with different numerical SRHD methods are compared, and two astrophysical applications of SRHD flows are discussed. An evaluation of the various numerical methods is given and future developments are analyzed.

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.

The main features of a three dimensional, high-resolution special relativistic hydro code based on relativistic Riemann solvers are described. The capabilities and performance of the code are discussed. In particular, we present the results of extensive test calculations which demonstrate that the code can accurately and efficiently handle strong shocks in three spatial dimensions. Results of the performance of the code on single and multi-processor machines are given. Simulations (in double precision) with 7x10^6 computational cells require less than 1 Gb of RAM memory and ~7x10^-5 CPU seconds per zone and time step (on a SCI Cray-Origin 2000 with a R10000 processor). Currently, a version of the numerical code is under development, which is suited for massively parallel computers with distributed memory architecture (like, e.g., Cray T3E).

We present a general and practical procedure to solve the general relativistic hydrodynamic equations by using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows. Our proposal relies on a local change of coordinates in terms of which the spacetime metric is locally Minkowskian and permits accurate numerical calculations of general relativistic hydrodynamics problems using the numerical tools developed for the special relativistic case with negligible computational cost. The feasibility of the method has been confirmed by a number of numerical experiments.

We present the results of a multiwavelength monitoring campaign targeting the gamma-ray source 2CG 135+1 in an attempt to confirm the association of this object with the radio/Be/X-ray binary system LSI +61o 303. The campaign included simultaneous radio, optical, infrared, and hard x-ray/gamma-ray observations carried out with a variety of instruments, covering (not continously) almost three binary cycles of LSI +61o 303 during the period April-July 1994. Three separate OSSE observations of the gamma-ray source were carried out, covering different phases of the radio lightcurve. Hard X-ray/gamma-ray emission was detected from the direction of 2CG 135+1 during the first of these OSSE observations. The signal to noise ratio of the OSSE observations was insufficient to establish a spectral or intensity correlation of the high-energy emission with simultaneous radio, optical and infrared emission of LSI +61o 303. We briefly discuss the theoretical implications of our observations.

We report simultaneous observations in the X-ray, infrared, and radio wavelengths of the galactic superluminal source GRS 1915+105. During episodes of rapid disappearance and follow up replenishment of the inner accretion disk evidenced by the X-ray oscillating flux, we observe the ejection of relativistic plasma clouds in the form of synchrotron flares at infrared and radio wavelengths. The expelled clouds contain very energetic particles with Lorentz factors of about 1000, or more. These ejections can be viewed as small-scale analogs of the more massive ejecta with relativistic bulk motions that have been previously observed in GRS 1915+105.

The radio, infrared and optical counterparts of the hard X-ray source GRS 1734-292 (l{II}=358.9 deg, b{II}=+1.4 deg) are reported. The optical spectrum exhibits broad (>1000 km s^-1) H alpha emission at a redshift of 0.0214 +/- 0.0005. The radio counterpart is a double-sided synchrotron jet of 5 arcsec, which at a distance of 87 Mpc corresponds to a size of 2 kpc. The multiwavelength observations of GRS 1734-292 indicate that this X-ray source is a Seyfert 1 galaxy behind 6 +/- 1 magnitudes of visual absorption.

We present simulations of relativistic jets propagating in a nonuniform medium. Specifically, we study the bending of jets propagating obliquely to the vector of the density-gradient. Our results are applied to the NGC 4258, where such a medium is assumed to be provided by the atmosphere of the sub-parsec accretion disk tilted with respect to the original direction of the jet propagation. As a result, the jet is bent on a scale comparable to the density scaleheight of the disk atmosphere. The magnitude of the bending effect is found to be largest for light jets with low Lorentz factors. The predicted direction of bending is consistent with the observations.

In this paper we present a full general relativistic one-dimensional hydro-code which incorporates a modern high-resolution shock-capturing algorithm, with an approximate Riemann solver, for the correct modelling of formation and propagation of strong shocks. The efficiency of this code in treating strong shocks is demonstrated by some numerical experiments. The interest of this technique in several astrophysical scenarios is discussed.

We present new optical and infrared photometric observations and high resolution H$\alpha$ spectra of the periodic radio star \lsi. The optical photometric data set covers the time interval 1985-1993 and amounts to about a hundred nights. A period of $\sim$26 days is found in the V band. The infrared data also present evidence for a similar periodicity, but with higher amplitude of variation (0\rmag 2). The spectroscopic observations include 16 intermediate and high dispersion spectra of \lsi\ collected between January 1989 and February 1993. The H$\alpha$ emission line profile and its variations are analyzed. Several emission line parameters -- among them the H$\alpha$ EW and the width of the H$\alpha$ red hump -- change strongly at or close to radio maximum, and may exhibit periodic variability. We also observe a significant change in the peak separation. The H$\alpha$ profile of \lsi\ does not seem peculiar for a Be star. However, several of the observed variations of the H$\alpha$ profile can probably be associated with the presence of the compact, secondary star.