A search for microquasar candidates at low galactic latitudes

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.

Period-Luminosity-Colour distribution and classification of Galactic oxygen-rich LPVs I. Luminosity calibrations

The absolute K magnitudes and kinematic parameters of about 350 oxygen-rich Long-Period Variable stars are calibrated, by means of an up-to-date maximum-likelihood method, using HIPPARCOS parallaxes and proper motions together with radial velocities and, as additional data, periods and V-K colour indices. Four groups, differing by their kinematics and mean magnitudes, are found. For each of them, we also obtain the distributions of magnitude, period and de-reddened colour of the base population, as well as de-biased period-luminosity-colour relations and their two-dimensional projections. The SRa semiregulars do not seem to constitute a separate class of LPVs. The SRb appear to belong to two populations of different ages. In a PL diagram, they constitute two evolutionary sequences towards the Mira stage. The Miras of the disk appear to pulsate on a lower-order mode. The slopes of their de-biased PL and PC relations are found to be very different from the ones of the Oxygen Miras of the LMC. This suggests that a significant number of so-called Miras of the LMC are misclassified. This also suggests that the Miras of the LMC do not constitute a homogeneous group, but include a significant proportion of metal-deficient stars, suggesting a relatively smooth star formation history. As a consequence, one may not trivially transpose the LMC period-luminosity relation from one galaxy to the other.

Optical counterpart of galactic plane variable radio sources

We present I-band deep CCD exposures of the fields of galactic plane radio variables. An optical counterpart, based on positional coincidence, has been found for 15 of the 27 observed program objects. The Johnson I magnitude of the sources identified is in the range 18-21.

The Besançon Galaxy Model renewed. I. Constraints on the Galactic thin disc evolution from Tycho data

Context. The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allow to test hypotheses on the star formation history, star evolution, as well as chemical and dynamical evolution of the Galaxy. Aims. The new version of the Besanc¸on Galaxy Model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the Initial Mass Function (IMF) and Star Formation Rate (SFR) of the Galactic disc. Methods. We present a new strategy for the generation of thin disc stars which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. We keep in this new scheme the local dynamical self-consistency as in Bienayme et al (1987). We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. Results. For the first time, a full sky comparison is performed between BGM and data. This strategy allows to constrain the IMF slope at high masses which is found to be close to 3.0, excluding a shallower slope such as Salpeter"s one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 M pc−3 implying that there is no compelling evidence for significant amount of dark matter in the disc. While the model is fitted to Tycho2 data, a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. Conclusions. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

Optical spectroscopy of microquasar candidates at low Galactic latitudes

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.

LS 5039: A runaway microquasar ejected from the galactic plane

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.

EVN+MERLIN observations of microquasar candidates at low galactic latitudes

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.

The Besançon Galaxy Model renewed. I. Constraints on the Galactic thin disc evolution from Tycho data

Context. The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allow to test hypotheses on the star formation history, star evolution, as well as chemical and dynamical evolution of the Galaxy. Aims. The new version of the Besanc¸on Galaxy Model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the Initial Mass Function (IMF) and Star Formation Rate (SFR) of the Galactic disc. Methods. We present a new strategy for the generation of thin disc stars which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. We keep in this new scheme the local dynamical self-consistency as in Bienayme et al (1987). We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. Results. For the first time, a full sky comparison is performed between BGM and data. This strategy allows to constrain the IMF slope at high masses which is found to be close to 3.0, excluding a shallower slope such as Salpeter"s one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 M pc−3 implying that there is no compelling evidence for significant amount of dark matter in the disc. While the model is fitted to Tycho2 data, a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. Conclusions. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

Probing quantum gravity using photons from a flare of the active galactic nucleus Markarian 501 observed by the MAGIC telescope

We analyze the timing of photons observed by the MAGIC telescope during a flare of the active galactic nucleus Mkn 501 for a possible correlation with energy, as suggested by some models of quantum gravity (QG), which predict a vacuum refractive index similar or equal to 1 + (E/M-QGn)(n), n = 1, 2. Parametrizing the delay between gamma-rays of different energies as Delta t = +/-tau E-1 or Delta t = +/-tau E-q(2), we find tau(1) = (0.030 +/- 0.012) s/GeV at the 2.5-sigma level, and tau(q) = (3.71 +/- 2.57) x 10(-6) s/GeV2, respectively. We use these results to establish lower limits M-QG1 > 0.21 X 10(18) GeV and M-QG2 > 0.26 x 10(11) GeV at the 95% C.L. Monte Carlo studies confirm the MAGIC sensitivity to propagation effects at these levels. Thermal plasma effects in the source are negligible, but we cannot exclude the importance of some other source effect.

Emissió ràdio a baixa freqüència en fenòmens eruptius i en fonts esteses de raigs gamma d'alta i molt alta energia

La tasca investigadora presentada en aquesta memòria s'ha centrat en les fonts galàctiques de raigs gamma de molt alta energia LS I +61 303, HESS J1708-410 i HESS J1858+020. La primera és una binària de raigs gamma molt estudiada, formada per una estrella massiva i un objecte compacte. S'ha proposat un escenari on l'objecte compacte seria un púlsar jove, i la interacció del seu vent amb el vent de l'estrella generaria els raigs gamma. De totes formes, no s'ha detectat polsos procedents d'aquest putatiu púlsar. L'investigador va realitzar observacions en fase a 1280 MHz amb el radiotelescopi GMRT, sense trobar-hi polsos, cosa que implica un estricte límit superior de 0,38 mJy a la densitat mitjana de flux polsat en un putatiu púlsar amb un període major que 2 mil•lisegons en el sistema binari LS I +61 303. Per altra banda, HESS J1708-410 i HESS J1858+020 són dues fonts esteses de raigs gamma de molt alta energia de les quals no es coneix cap contrapart a d'altres longituds d'ona. L'investigador les va observar amb el GMRT, quatre vegades HESS J1708-410 (dues a 610 MHz i dues a 1400 MHz) i dues vegades HESS J1858+020 (una a cada freqüència). En les imatges realitzades amb aquestes dades no hi ha emissió estesa coincident amb les regions d'emissió de raigs gamma. HESS J1858+020 se solapa parcialment amb una font estesa que podria ser un SNR. De confirmar-se la falta de contrapartida ràdio de HESS J1708-410, estaríem parlant d'un accelerador hadrònic extraordinàriament eficient, d'una classe desconeguda fins ara.

Electromagnetic radiation initiated by hadronic jets from microquasars in the ISM

Microquasars are potential candidates to produce a non-negligible fraction of the observed galactic cosmic rays. The protons accelerated at the jet termination shock interact with the interstellar medium and may produce detectable fluxes of extended emission at different energy bands: high-energy and very high-energy gamma-rays produced by neutral pion-decay, synchrotron and bremsstrahlung emission in a wide energy range generated by the secondary electrons produced by charged pion-decay. We discuss the association between this scenario and some of the unidentified EGRET sources in the galactic plane.

A new maximum likelihood method for luminosity calibrations

A new statistical parallax method using the Maximum Likelihood principle is presented, allowing the simultaneous determination of a luminosity calibration, kinematic characteristics and spatial distribution of a given sample. This method has been developed for the exploitation of the Hipparcos data and presents several improvements with respect to the previous ones: the effects of the selection of the sample, the observational errors, the galactic rotation and the interstellar absorption are taken into account as an intrinsic part of the formulation (as opposed to external corrections). Furthermore, the method is able to identify and characterize physically distinct groups in inhomogeneous samples, thus avoiding biases due to unidentified components. Moreover, the implementation used by the authors is based on the extensive use of numerical methods, so avoiding the need for simplification of the equations and thus the bias they could introduce. Several examples of application using simulated samples are presented, to be followed by applications to real samples in forthcoming articles.

Microquasars as high energy gamma-ray sources

Galactic microquasars are certainly one of the most recent additions to the field of high energy astrophysics and have attracted increasing interest over the last decade. However, the high energy part of the spectrum of microquasars is the most poorly known, mainly due the lack of sensitive instrumentation in the past. Microquasars are now primary targets for all of the observatories working in the X-ray and gamma-ray domains. They also appear as the possible counterparts for some of the unidentified sources of high-energy gamma-rays detected by the experiment EGRET on board the satellite COMPTON-GRO. This paper provides a general review of the main observational results obtained up to now as well as a summary of the scenarios for production of high-energy gamma-rays at the present moment.

CP2 stars as viewed by the uvby Hbeta system

A new statistical parallax method using the Maximum Likelihood principle is presented, allowing the simultaneous determination of a luminosity calibration, kinematic characteristics and spatial distribution of a given sample. This method has been developed for the exploitation of the Hipparcos data and presents several improvements with respect to the previous ones: the effects of the selection of the sample, the observational errors, the galactic rotation and the interstellar absorption are taken into account as an intrinsic part of the formulation (as opposed to external corrections). Furthermore, the method is able to identify and characterize physically distinct groups in inhomogeneous samples, thus avoiding biases due to unidentified components. Moreover, the implementation used by the authors is based on the extensive use of numerical methods, so avoiding the need for simplification of the equations and thus the bias they could introduce. Several examples of application using simulated samples are presented, to be followed by applications to real samples in forthcoming articles.