Gamma-ray binaries are systems that emit non-thermal radiation peaking at energies above 1 MeV. One proposed scenario to explain their emission consists of a pulsar orbiting a massive star, with particle acceleration taking place in shocks produced by the interaction of the stellar and pulsar winds. We develop a semi-analytical model of the non-thermal emission of the colliding-wind structure including the dynamical effects of orbital motion. We apply the model to a general case and to LS 5039. The model consists of a one-dimensional emitter the geometry of which is affected by Coriolis forces owing to orbital motion. Two particle accelerators are considered: one at the two-wind standoff location, and the other one at the turnover produced by the Coriolis force. Synchrotron and inverse Compton emission is studied, accounting for Doppler boosting and absorption processes associated to the massive star. If both accelerators are provided with the same energy budget, most of the radiation comes from the region of the Coriolis turnover and beyond, up to a few orbital separations from the binary system. The model allows us to reproduce some of the LS 5039 emission features, but not all of them. In particular, the MeV radiation is probably too high to be explained by our model alone, the GeV flux is recovered but not its modulation, and the radio emission beyond the Coriolis turnover is too low. The predicted system inclination is consistent with the presence of a pulsar in the binary. The model is quite successful in reproducing the overall non-thermal behavior of LS 5039. Some improvements are suggested to better explain the phenomenology observed in this source, like accounting for particle reacceleration beyond the Coriolis turnover, unshocked pulsar wind emission, and the three-dimensional extension of the emitter.

Dark matter may consist, at least partially, of primordial black holes formed during the radiation-dominated era. The radiation produced by accretion onto primordial black holes leaves characteristic signatures on the properties of the medium at high redshift, before and after Hydrogen recombination. Therefore, reliable modelling of accretion onto these objects is required to obtain robust constraints on their abundance. We investigate the effect of mechanical feedback, i.e. the impact of outflows (winds and/or jets) on the medium, on primordial black hole accretion, and thereby on the associated radiation. Using analytical and numerical calculations, we study for the first time whether outflows can reduce the accretion rate of primordial black holes with masses similar to those detected by the LIGO-Virgo collaboration. Despite the complexity of the accretion rate evolution, mechanical feedback is able to significantly reduce the primordial black hole accretion rate, at least by an order of magnitude, when outflows are aligned with the motion of the compact object. If the outflow is perpendicular to the direction of motion, the effect is less important but still non-negligible. Outflows from primordial black holes, even rather weak ones, can significantly decrease the accretion rate, effectively weakening abundance constraints on these objects. Our results motivate further numerical simulations with a more realistic setup, which would yield more precise quantitative predictions.

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $\gamma$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $\gamma$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries.

Cosmic reionization put an end to the dark ages that came after the recombination era. Observations seem to favor the scenario where massive stars generating photons in low-mass galaxies were responsible for the bulk of reionization. Even though a possible contribution from accretion disks of active galactic nuclei (AGN) has been widely considered, they are currently thought to have had a minor role in reionization. Our aim is to study the possibility that AGN contributed to reionization not only through their accretion disks, but also through ionizing photons coming from the AGN jets interacting with the IGM. We adopt an empirically derived AGN luminosity function at $z\simeq6$, use X-ray observations to correct it for the presence of obscured sources, and estimate the density of jetted AGN. We then use analytical calculations to derive the fraction of jet energy that goes into ionizing photons. Finally, we compute the contribution of AGN jets to the H II volume filling factor at redshifts $z\simeq15-5$. We show that the contribution of the AGN jet lobes to the reionization of the Universe at $z\sim6$ might have been as high as $\gtrsim 10$\% of that of star-forming galaxies, under the most favorable conditions of jetted and obscuration fraction. The contribution of AGN to the reionization, while most likely not dominant, could have been higher than previously assumed, thanks to the radiation originated in the jet lobes.

The persistent gamma-ray source 4FGL J0647.7-4418 is tentatively associated in the latest Fermi catalogue with the sub-dwarf O-type X-ray binary HD 49798. However, an AGN candidate is also mentioned as an alternative identification in updated versions of the catalogue accompanying paper. If the first association were correct, this would add HD 49798 to the handful of currently known gamma-ray binaries and, therefore, represent a significant breakthrough not only because of a new member addition, but also because of the apparent white dwarf companion in this system. Despite these perspectives, here we show that the stellar association is likely wrong and that the proposed AGN object, well inside the Fermi 95% confidence ellipse, is a more conceivable counterpart candidate to the Fermi source due to its strong blazar similarities.

High-mass microquasars (HMMQs) are systems from which relativistic jets are launched. At the scales of several times the binary system size, the jets are expected to follow a helical path caused by the interaction with a strong stellar wind and orbital motion. Such a trajectory has its influence on the non-thermal emission of the jets, which also depends strongly on the observing angle due to Doppler boosting effects. We explore how the expected non-thermal emission of HMMQ jets at small scales is affected by the impact of the stellar wind and the orbital motion on the jet propagation. We studied the broadband non-thermal emission, from radio to gamma rays, produced in HMMQ jets up to a distance of several orbital separations, taking into account a realistic jet trajectory, different model parameters, and orbital modulation. The jet trajectory is computed by considering momentum transfer with the stellar wind. Electrons are injected at the position where a recollimation shock in the jets is expected due to the wind impact. Their distribution along the jet path is obtained assuming local acceleration at the recollimation shock, and cooling via adiabatic, synchrotron, and inverse Compton processes. The synchrotron and inverse Compton emission is calculated taking into account synchrotron self-absorption within the jet, free-free absorption with the stellar wind, and absorption by stellar photons via pair production. Asymmetric light curves are obtained owing to the helical trajectory of the jets. The presence of helical shaped jets could be inferred from asymmetries in the light curves, which become noticeable only for large jet Lorentz factors and low magnetic fields. Model parameters could be constrained if accurate phase-resolved light curves from GeV to TeV energies were available.

The innermost parts of powerful jets in active galactic nuclei are surrounded by dense, high-velocity clouds from the broad-line region, which may penetrate into the jet and lead to the formation of a strong shock. Such jet-cloud interactions are expected to have measurable effects on the $\gamma$-ray emission from blazars. We characterise the dynamics of a typical cloud-jet interaction scenario, and the evolution of its radiative output in the 0.1-30~GeV energy range, to assess to what extent these interactions can contribute to the $\gamma$-ray emission in blazars. We use semi-analytical descriptions of the jet-cloud dynamics, taking into account the expansion of the cloud inside the jet and its acceleration. Assuming that electrons are accelerated in the interaction and making use of the hydrodynamical information, we then compute the high-energy radiation from the cloud, including the absorption of $\gamma$-rays in the ambient photon field through pair creation. Jet-cloud interactions can lead to significant $\gamma$-ray fluxes in blazars with a broad-line region, in particular when the cloud expansion and acceleration inside the jet are taken into account. This is caused by 1) the increased shocked area in the jet, which leads to an increase in the energy budget for the non-thermal emission; 2) a more efficient inverse Compton cooling with the boosted photon field of the broad-line region; and 3) an increased observer luminosity due to Doppler boosting effects. For typical broad-line region parameters, either (i) jet-cloud interactions contribute significantly to the persistent $\gamma$-ray emission from blazars or (ii) the broad-line region is far from spherical or the fraction of energy deposited in non-thermal electrons is small.

Blazars in elliptical hosts have a population of red giants surrounding their jet. These stars can carry large wind-blown bubbles into the jets, leading to gamma-ray emission through bubble-jet interactions. We study the interaction dynamics and the gamma-ray emission produced when the bubbles formed by red giant winds penetrate the jet of a blazar in an elliptical galaxy. First, we characterized the masses and penetration rates of the red giant wind bubbles that enter the jet. Then, the dynamical evolution of these bubbles under the jet impact was analysed analytically and numerically, and the radiation losses of the particles accelerated in the interaction were characterised. Finally, the synchrotron and the inverse Compton contributions above 100 MeV were estimated under different jet magnetic fields, powers, and Lorentz factors. We find that an analytical dynamical model is a reasonable approximation for the red giant wind bubble-jet interaction. The radiation produced by these wind bubbles interacting with a jet can have a duty cycle of up to 1. For realistic magnetic fields, gamma rays could be detectable from sources within the local universe, preferentially from those with high Lorentz factors (of about 10), and this could be a relatively common phenomenon for these sources. For magnetic fields in equipartition with the jet power, and high acceleration rates, synchrotron gamma rays may be detectable even for modest Lorentz factors (of about 3), but with a much lower duty cycle.

Core-collapse supernovae are found in galaxies with ongoing star-formation. In a starburst galaxy hosting an active galactic nucleus with a relativistic jet, supernovae can take place inside the jet. The collision of the supernova ejecta with the jet flow is expected to lead to the formation of an interaction region, in which particles can be accelerated and produce high-energy emission. We study the non-thermal radiation produced by electrons accelerated as a result of a supernova explosion inside the jet of an active galactic nucleus within a star-forming galaxy. We first analyzed the dynamical evolution of the supernova ejecta impacted by the jet. Then, we explored the parameter space using simple prescriptions for the observed gamma-ray lightcurve. Finally, the synchrotron and the inverse Compton spectral energy distributions for two types of sources, a radio galaxy and a powerful blazar, are computed. For a radio galaxy, the interaction between a supernova and a jet of power $\sim 10^{43}-10^{44}$~erg~s$^{-1}$ can produce apparent gamma-ray luminosities of $\sim 10^{42}-10^{43}$~erg~s$^{-1}$, with an event duty cycle of supernova remnant (SNR) interacting with the jet close to one for one galaxy. For a blazar with a powerful jet of $\sim 10^{46}$~erg~s$^{-1}$, the jet-supernova ejecta interaction could produce apparent gamma-ray luminosities of $\sim 10^{43}-10^{44}$~erg~s$^{-1}$, but with a much lower duty cycle. The interaction of supernovae with misaligned jets of moderate power can be relatively frequent, and can result in steady gamma-ray emission potentially detectable for sources in the local universe. For powerful blazars much farther away, the emission would be steady as well, and it might be detectable under very efficient acceleration, but the events would be rather infrequent.

The stellar wind in high-mass microquasars should interact with the jet. This interaction, coupled with orbital motion, is expected to make the jet follow a helical, nonballistic trajectory. The jet energy dissipated by this interaction, through shocks for example, could lead to nonthermal activity on scales significantly larger than the system size. We calculate the broadband emission from a jet affected by the impact of the stellar wind and orbital motion in a high-mass microquasar. We employ a prescription for the helical trajectory of a jet in a system with a circular orbit. Subsequently, assuming electron acceleration at the onset of the helical jet region, we compute the spatial and energy distribution of these electrons, and their synchrotron and inverse Compton emission including gamma-ray absorption effects. For typical source parameters, significant radio, X- and gamma-ray luminosities are predicted. The scales on which the emission is produced may reduce, but not erase, orbital variability of the inverse Compton emission. The wind and orbital effects on the radio emission morphology could be studied using very long baseline interferometric techniques. We predict significant broadband emission, modulated by orbital motion, from a helical jet in a high-mass microquasar. This emission may be hard to disentangle from radiation of the binary itself, although the light curve features, extended radio emission, and a moderate opacity to very high-energy gamma rays, could help to identify the contribution from an extended (helical) jet region.

Microquasars are binary systems consisting of a compact object and a star that present relativistic jets. When the companion is a high-mass star, significant interaction between the stellar wind and the jets is expected. In this work, an overview of the different effects that a strong stellar wind may have in microquasar jets is given. Both analytical estimates and numerical simulations performed in the last years are reviewed. The results of a model for the non-thermal emission of such jets at large scales are also shown and discussed. Observational studies of two high-mass microquasars, Cygnus X-1 and Cygnus X-3, are compared to the model predictions.

The reionization of the Universe ends the dark ages that started after the recombination era. In the case of H, reionization finishes around $z\sim 6$. Faint star-forming galaxies are the best candidate sources of the H-ionizing radiation, although active galactic nuclei may have also contributed. We have explored whether the termination regions of the jets from active galactic nuclei may have contributed significantly to the ionization of H in the late reionization epoch, around $z\sim 6-7$. We assumed that, as it has been proposed, active galactic nuclei at $z\sim 6$ may have presented a high jet fraction, accretion rate, and duty cycle, and that non-thermal electrons contribute significantly to the pressure of jet termination regions. Empirical black-hole mass functions were adopted to characterize the population of active galactic nuclei. From all this, estimates were derived for the isotropic H-ionizing radiation produced in the jet termination regions, at $z\sim 6$, through inverse Compton scattering off CMB photons. We find that the termination regions of the jets of active galactic nuclei may have radiated most of their energy in the form of H-ionizing radiation at $z\sim 6$. For typical black-hole mass functions at that redshift, under the considered conditions (long-lasting, common, and very active galactic nuclei with jets), the contribution of these jets to maintain (and possibly enhance) the ionization of H may have been non-negligible. We conclude that the termination regions of jets from active galactic nuclei could have had a significant role in the reionization of the Universe at $z\gtrsim 6$.

Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous observations of FRB 121102 with the Arecibo and MAGIC telescopes during several epochs in 2016--2017. This allowed searches for millisecond-timescale burst emission in very-high-energy (VHE) gamma rays as well as the optical band. While a total of five FRBs were detected during these observations, no VHE emission was detected, neither of a persistent nature nor burst-like associated with the FRBs. The average integral flux upper limits above 100 GeV at 95% confidence level are $6.6 \times 10^{-12}~\mathrm{photons\ cm^{-2}\ s^{-1}}$ (corresponding to luminosity $L_{\rm VHE} \lesssim 10^{45}~\mathrm{erg\ s^{-1}}$) over the entire observation period, and $1.2 \times 10^{-7}~ \mathrm{photons\ cm^{-2}\ s^{-1}}$ ($L_{\rm VHE} \lesssim 10^{49}~\mathrm{erg\ s^{-1}}$) over the total duration of the five FRBs. We constrain the optical U-band flux to be below 8.6 mJy at 5-$\sigma$ level for 1-ms intervals around the FRB arrival times. A bright burst with U-band flux $29~\mathrm{mJy}$ and duration $\sim 12$ ms was detected 4.3 s before the arrival of one FRB. However, the probability of spuriously detecting such a signal within the sampled time space is 1.5% (2.2 $\sigma$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models.

The gamma-ray emission detected from several microquasars can be produced by relativistic electrons emitting through inverse Compton scattering. In particular, the GeV emission detected from Cygnus X-3, and its orbital phase dependence, strongly suggest that the emitting electrons are accelerated in a relativistic jet, and that the optical companion provides the dominant target. Here, we study the effects related to particle transport in the framework of the relativistic jet scenario. We find that even in the most compact binary systems, with parameters similar to Cygnus X-3, particle transport can have a substantial influence on the GeV lightcurve unless the jet is slow, $\beta < 0.7$. In more extended binary systems, strong impact of particle transport is nearly unavoidable. Thus, even for a very compact system such as Cygnus X-3, particle transport significantly affects the ability of one-zone models to infer the properties of the gamma-ray production site based on the shape on the GeV lightcurve. We conclude that a detailed study of the gamma-ray spectrum can further constrain the structure and other properties of the gamma-ray emitter in Cygnus X-3, although such a study should account for gamma-gamma attenuation, since it may strongly affect the spectrum above $5\rm\,GeV$.

Context: Runaway stars produce bowshocks that are usually observed at infrared (IR) wavelengths. Non-thermal radio emission has been detected so far only from the bowshock of BD+43{\deg}3654, whereas the detection of non-thermal radiation from these bowshocks at high energies remains elusive. Aims: We aim at characterising in detail the radio, X-ray, and gamma-ray emission from stellar bowshocks accounting for the structure of the region of interaction between the stellar wind and its environment. Methods: We develop a broadband-radiative, multi-zone model for stellar bowshocks that takes into account the spatial structure of the emitting region and the observational constraints. The model predicts the evolution and the emission of the relativistic particles accelerated and streaming together with the shocked flow. Results: We present broadband non-thermal spectral energy distributions for different scenarios, synthetic radio-cm synchrotron maps that reproduce the morphology of BD+43{\deg}3654, and updated predictions in X-ray and gamma-ray energy ranges. We also compare the results of the multi-zone model applied in this work with those of a refined one-zone model. Conclusions: A multi-zone model provides better constraints than a one-zone model on the relevant parameters, namely the magnetic field intensity and the amount of energy deposited in non-thermal particles. However, one-zone models can be improved by carefully characterising the intensity of the IR dust photon field and the escape rate of the plasma from the shocked region. Finally, comparing observed radio maps with those obtained from a multi-zone model enables constraints to be obtained on the direction of stellar motion with respect to the observer.

Relativistic hydrodynamical simulations of the eccentric gamma-ray binary HESS J0632$+$057, show that the energy of a putative pulsar wind should accumulate in the binary surroundings between periastron and apastron, being released by fast advection close to apastron. To assess whether this could lead to a maximum of the non-thermal emission before apastron, we derive simple prescriptions for the non-thermal energy content, the radiation efficiency, and the impact of energy losses on non-thermal particles, in the simulated hydrodynamical flow. These prescriptions are used to estimate the non-thermal emission in radio, X-rays, GeV, and TeV, from the shocked pulsar wind in a binary system simulated using a simplified 3-dimensional scheme for several orbital cycles. Lightcurves at different wavelengths are derived, together with synthetic radio images for different orbital phases. The dominant peak in the computed lightcurves is broad and appears close to, but before, apastron. This peak is followed by a quasi-plateau shape, and a minor peak only in gamma rays right after periastron. The radio maps show ejection of radio blobs before apastron in the periastron-apastron direction. The results show that a scenario with a highly eccentric high-mass binary hosting a young pulsar can explain the general phenomenology of HESS J0632$+$057: despite its simplicity, the adopted approach yields predictions that are robust at a semi-quantitative level and consistent with multiwavelength observations.

The family links between radio galaxies and microquasars have been strongly strengthened thanks to a new common phenomenon: the presence of extended winged features. The first detection of such structures in a Galactic microquasar, recently reported in Nature Communications (http://rdcu.be/zgX8), widens the already known analogy between both kinds of outflow sources (Mart\'i et al. 2017). This observational result also has potential implications affecting the black hole merger scenarios that contribute to the gravitational wave background.

Gamma-ray observations of microquasars at high and very-high energies can provide valuable information of the acceleration processes inside the jets, the jet-environment interaction and the disk-jet coupling. Two high-mass microquasars have been deeply studied to shed light on these aspects: Cygnus X-1 and Cygnus X-3. Both systems display the canonical hard and soft X-ray spectral states of black hole transients, where the radiation is dominated by non-thermal emission from the corona and jets and by thermal emission from the disk, respectively. Here, we report on the detection of Cygnus X-1 above 60 MeV using 7.5 yr of Pass8 Fermi-LAT data, correlated with the hard X-ray state. A hint of orbital flux modulation was also found, as the source is only detected in phases around the compact object superior conjunction. We conclude that the high-energy gamma-ray emission from Cygnus X-1 is most likely associated with jets and its detection allow us to constrain the production site. Moreover, we include in the discussion the final results of a MAGIC long-term campaign on Cygnus X-1 that reaches almost 100 hr of observations at different X-ray states. On the other hand, during summer 2016, Cygnus X-3 underwent a flaring activity period in radio and high-energy gamma rays, similar to the one that led to its detection in the high-energy regime in 2009. MAGIC performed comprehensive follow-up observations for a total of about 70 hr. We discuss our results in a multi-wavelength context.

The microquasar Cygnus X-1 displays the two typical soft and hard X-ray states of a black-hole transient. During the latter, Cygnus X-1 shows a one-sided relativistic radio-jet. Recent detection of the system in the high energy (HE; $E\gtrsim60$ MeV) gamma-ray range with \textit{Fermi}-LAT associates this emission with the outflow. Former MAGIC observations revealed a hint of flaring activity in the very high-energy (VHE; $E\gtrsim100$ GeV) regime during this X-ray state. We analyze $\sim97$ hr of Cygnus X-1 data taken with the MAGIC telescopes between July 2007 and October 2014. To shed light on the correlation between hard X-ray and VHE gamma rays as previously suggested, we study each main X-ray state separately. We perform an orbital phase-folded analysis to look for variability in the VHE band. Additionally, to place this variability behavior in a multiwavelength context, we compare our results with \textit{Fermi}-LAT, \textit{AGILE}, \textit{Swift}-BAT, \textit{MAXI}, \textit{RXTE}-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the VHE regime. We establish upper limits for each X-ray state, assuming a power-law distribution with photon index $\Gamma=3.2$. For steady emission in the hard and soft X-ray states, we set integral upper limits at 95\% confidence level for energies above 200 GeV at $2.6\times10^{-12}$~photons cm$^{-2}$s$^{-1}$ and $1.0\times10^{-11}$~photons cm$^{-2}$s$^{-1}$, respectively. We rule out steady VHE gamma-ray emission above this energy range, at the level of the MAGIC sensitivity, originating in the interaction between the relativistic jet and the surrounding medium, while the emission above this flux level produced inside the binary still remains a valid possibility.

HESS J0632+057 is an eccentric gamma-ray Be binary that produces non-thermal radio, X-rays, GeV, and very high-energy gamma rays. The non-thermal emission of HESS J0632+057 is modulated with the orbital period, with a dominant maximum before apastron passage. The nature of the compact object in HESS J0632+057 is not known, although it has been proposed to be a young pulsar as in PSR B1259-63, the only gamma-ray emitting high-mass binary known to host a non-accreting pulsar. In this Letter, we present hydrodynamical simulations of HESS J0632+057 in the context of a pulsar and a stellar wind interacting in an eccentric binary, and propose a scenario for the non-thermal phenomenology of the source. In this scenario, the non-thermal activity before and around apastron is linked to the accumulation of non-thermal particles in the vicinity of the binary, and the sudden drop of the emission before apastron is produced by the disruption of the two-wind interaction structure, allowing these particles to efficiently escape. In addition to providing a framework to explain the non-thermal phenomenology of the source, this scenario predicts extended, moving X-ray emitting structures similar to those observed in PSR B1259-63.

Dense stellar winds may mass-load the jets of active galactic nuclei, although it is unclear what are the time and spatial scales in which the mixing takes place. We study the first steps of the interaction between jets and stellar winds, and also the scales at which the stellar wind may mix with the jet and mass-load it. We present a detailed two-dimensional simulation, including thermal cooling, of a bubble formed by the wind of a star. We also study the first interaction of the wind bubble with the jet using a three-dimensional simulation in which the star enters the jet. Stability analysis is carried out for the shocked wind structure, to evaluate the distances over which the jet-dragged wind, which forms a tail, can propagate without mixing with the jet flow. The two-dimensional simulations point at quick wind bubble expansion and fragmentation after about one bubble shock crossing time. Three-dimensional simulations and stability analysis point at local mixing in the case of strong perturbations and relatively small density ratios between the jet and the jet dragged-wind, and to a possibly more stable shocked wind structure at the phase of maximum tail mass flux. Analytical estimates also indicate that very early stages of the star jet-penetration time may be also relevant for mass loading. The combination of these and previous results from the literature suggest highly unstable interaction structures and efficient wind-jet flow mixing on the scale of the jet interaction height, possibly producing strong inhomogeneities within the jet. In addition, the initial wind bubble shocked by the jet leads to a transient, large interaction surface. The interaction structure can be a source of significant non-thermal emission.

Fast radio bursts, or FRBs, are transient sources of unknown origin. Recent radio and optical observations have provided strong evidence for an extragalactic origin of the phenomenon and the precise localization of the repeating FRB 121102. Observations using the Karl G. Jansky Very Large Array (VLA) and very-long-baseline interferometry (VLBI) have revealed the existence of a continuum non-thermal radio source consistent with the location of the bursts in a dwarf galaxy. All these new data rule out several models that were previously proposed, and impose stringent constraints to new models. We aim to model FRB 121102 in light of the new observational results in the active galactic nucleus (AGN) scenario. We propose a model for repeating FRBs in which a non-steady relativistic $e^\pm$-beam, accelerated by an impulsive magnetohydrodynamic (MHD)-driven mechanism, interacts with a cloud at the centre of a star-forming dwarf galaxy. The interaction generates regions of high electrostatic field called cavitons in the plasma cloud. Turbulence is also produced in the beam. These processes, plus particle isotropization, the interaction scale, and light retardation effects, provide the necessary ingredients for short-lived, bright coherent radiation bursts. The mechanism studied in this work explains the general properties of FRB 121102, and may also be applied to other repetitive FRBs. Coherent emission from electrons and positrons accelerated in cavitons provides a plausible explanation of FRBs.

The central regions of galaxies are complex environments, rich in evolved and/or massive stars. For galaxies hosting an active galactic nucleus (AGN) with jets, the interaction of the jets with the winds of the stars within can lead to particle acceleration, and to extended high-energy emitting regions. We compute the non-thermal emission produced by the jet flow shocked by stellar winds on the jet scale, far from the jet-star direct interaction region. First, prescriptions for the winds of the relevant stellar populations in different types of galaxies are obtained. The scenarios adopted include galaxies with their central regions dominated by old or young stellar populations, and with jets of different power. Then, we estimate the available energy to accelerate particles in the jet shock, and compute the transport and energy evolution of the accelerated electrons, plus their synchrotron and inverse Compton emission, in the shocked flow along the jet. A significant fraction of the jet energy, $\sim 0.1-10$\%, can potentially be available for the particles accelerated in jet-wind shocks in the studied cases. The non-thermal particles can produce most of the high-energy radiation on jet scales, far from the jet shock region. This high-energy emission will be strongly enhanced in jets aligned with the line of sight due to Doppler boosting effects. The interaction of relativistic jets with stellar winds may contribute significantly to the persistent high-energy emission in some AGNs with jets. However, in the particular case of M87, this component seems too low to explain the observed gamma-ray fluxes.

The stellar winds of the massive stars in high-mass microquasars are thought to be inhomogeneous. The interaction of these inhomogeneities, or clumps, with the jets of these objects may be a major factor in gamma-ray production. Our goal is to characterize a typical scenario of clump-jet interaction, and calculate the contribution of these interactions to the gamma-ray emission from these systems. We use axisymmetric, RHD simulations to model the emitting flow in a typical clump-jet interaction. Using the simulation results we perform a numerical calculation of the high-energy emission from one of these interactions. The radiative calculations are done for relativistic electrons locally accelerated at the jet shock, and the synchrotron and IC radiation spectra are computed for different stages of the shocked clump evolution. We also explore different parameter values, such as viewing angle and magnetic field strength. We generalize phenomenologically the results derived from one clump-jet interaction to multiple interactions under different wind models, estimating the clump-jet interaction rates, and the resulting luminosities in the GeV range. If particles are efficiently accelerated in clump-jet interactions, the apparent gamma-ray luminosity through IC scattering with the stellar photons can be significant even for rather strong magnetic fields and thus efficient synchrotron cooling. Moreover, despite the standing nature or slow motion of the jet shocks for most of the interaction stage, Doppler boosting in the postshock flow is relevant even for mildly relativistic jets. For clump-to-average wind density contrasts $\gtrsim 10$, clump-jet interactions could be bright enough to match the observed GeV luminosity in Cyg~X-1 and Cyg~X-3 when a jet is present in these sources, with required non-thermal-to-total available power fractions $\sim 0.01$ and $0.1$, respectively.

Stars and their winds can contribute to the non-thermal (NT) emission in extragalactic jets. Given the complexity of jet-star interactions, the properties of the resulting emission are strongly linked to those of the emitting flows. We simulate the interaction between a stellar wind and a relativistic extragalactic jet and use the hydrodynamic results to compute the NT emission under different conditions. We perform relativistic axisymmetric hydrodynamical simulations of a relativistic jet interacting with a supersonic, non-relativistic stellar wind. We compute the corresponding streamlines out of the simulation results, and calculate the injection, evolution, and emission of NT particles accelerated in the jet shock, focusing on electrons or $e^\pm$-pairs. Several cases are explored, considering different jet-star interaction locations, magnetic fields and observer lines of sight. The jet luminosity and star properties are fixed, but the results are easily scalable under changes of these parameters. Individual jet-star interactions produce synchrotron and inverse Compton (IC) emission that peaks from X-rays to MeV energies (depending on the magnetic field), and at $\sim 100-1000$ GeV (depending on the stellar type), respectively. The radiation spectrum is hard in the scenarios explored here due to non-radiative cooling dominance, as low-energy electrons are efficiently advected even under relatively high magnetic fields. Interactions of jets with cold stars lead to an even harder IC spectrum because of the Klein-Nishina effect in the cross-section. Doppler boosting has a strong impact on the observer luminosity. The emission levels for individual interactions found here are in the line of previous, more approximate, estimates, strengthening the hypothesis that collective jet-star interactions could significantly contribute at high energies under efficient particle acceleration.

Binaries hosting a massive star and a non-accreting pulsar are powerful non-thermal emitters due to the interaction of the pulsar and the stellar wind. The winds of massive stars are thought to be inhomogeneous, which could have an impact on the non-thermal emission. We study numerically the impact of the presence of inhomogeneities or clumps in the stellar wind on the high-energy non-thermal radiation of high-mass binaries hosting a non-accreting pulsar. We compute the trajectories and physical properties of the streamlines in the shocked pulsar wind without clumps, with a small clump, and with a large one. This information is used to compute the synchrotron and inverse Compton emission from the non-thermal populations, accounting also for the effect of gamma-ray absorption through pair creation. A specific study is done for PSR B1259-63/LS2883. When stellar wind clumps perturb the two-wind interaction region, the associated non-thermal radiation in the X-ray band,of synchrotron origin, and in the GeV-TeV band, of inverse Compton origin, is affected by several effects: (i) strong changes in the the plasma velocity direction that result in Doppler boosting factor variations, (ii) strengthening of the magnetic field that mainly enhances the synchrotron radiation, (iii) strengthening of the pulsar wind kinetic energy dissipation at the shock, potentially available for particle acceleration, and (iv) changes in the rate of adiabatic losses that affect the lower energy part of the non-thermal particle population. The radiation above 100 GeV detected, presumably, during the post-periastron crossing of the Be star disc in PSR B1259-63/LS2883, can be roughly reproduced assuming that the crossing of the disc is modeled as the encounter with a large inhomogeneity.

Aims: Probe the high-energy ($>$60 MeV) emission from the black hole X-ray binary system, Cygnus X-1, and investigate its origin. Methods: We analysed 7.5 yr of data by Fermi/LAT with the latest PASS8 software version. Results: We report the detection of a signal at $\sim$8 $\sigma$ statistical significance spatially coincident with Cygnus X-1 and a luminosity above 60 MeV of 5.5$\times$10$^{33}$ erg s$^{-1}$. The signal is correlated with the hard X-ray flux: the source is observed at high energies only during the hard X-ray spectral state, when the source is known to display persistent, relativistic radio emitting jets. The energy spectrum, extending up to $\sim$20 GeV without any sign of spectral break, is well fitted by a power-law function with a photon index of 2.3$\pm$0.2. There is a hint of orbital flux variability, with high-energy emission mostly coming around the superior conjunction. Conclusions: We detected GeV emission from Cygnus X-1 and probed that the emission is most likely associated with the relativistic jets. The evidence of flux orbital variability points to the anisotropic inverse Compton on stellar photons as the mechanism at work, thus constraining the emission region to a distance $10^{11}-10^{13}$ cm from the black hole.

The binary stellar system HD 93129A is one of the most massive known binaries in our Galaxy. This system presents non-thermal emission in the radio band, which can be used to infer its physical conditions and predict its emission in the high-energy band. We intend to constrain some of the unknown parameters of HD 93129A through modelling the non-thermal emitter, and also to analyse the detectability of this source in hard X-rays and $\gamma$-rays. We develop a broadband radiative model for the wind-collision region taking into account the evolution of the accelerated particles streaming along the shocked region, the emission by different radiative processes, and the attenuation of the emission propagating through the local matter and radiation fields. From the analysis of the radio emission, we find that the binary HD~93129A is more likely to have a low inclination and a high eccentricity. The minimum energy of the non-thermal electrons seems to be between $\sim 20 - 100$MeV, depending on the intensity of the magnetic field in the wind-collision region. The latter can be in the range $\sim 20 - 1500$ mG. Our model is able to reproduce the observed radio emission, and predicts that the non-thermal radiation from HD~93129A will increase in the near future. With instruments such as \textit{NuSTAR}, \textit{Fermi}, and CTA, it will be possible to constrain the relativistic particle content of the source, and other parameters such as the magnetic field strength in the wind collision zone, which in turn can be used to obtain upper-limits of the magnetic field on the surface of the very massive stars, thereby inferring whether magnetic field amplification is taking place in the particle acceleration region.

High-mass microquasar jets propagate under the effect of the wind from the companion star, and the orbital motion of the binary system. The stellar wind and the orbit may be dominant factors determining the jet properties beyond the binary scales.} An analytical study is performed to characterize the effects of the stellar wind and the orbital motion on the jet properties. Accounting for the wind thrust transferred to the jet, we derive analytical estimates to characterize the jet evolution under the impact of the stellar wind. We include the Coriolis force effect, induced by orbital motion and enhanced by the wind presence. Large-scale evolution of the jet is sketched accounting for wind-to-jet thrust transfer, total energy conservation, and wind-jet flow mixing. If the angle of the wind-induced jet bending is larger than its half-opening angle, the following is expected: (i) a strong recollimation shock; (ii) bending against orbital motion, caused by Coriolis forces and enhanced by the wind presence; and (iii) non-ballistic helical propagation further away. Even if disrupted, the jet can reaccelerate due to ambient pressure gradients, but wind entrainment can weaken this acceleration. On large scales, the opening angle of the helical structure is determined by the wind-jet thrust relation, and the wind-loaded jet flow can be rather slow. The impact of stellar winds on high-mass microquasar jets can yield non-ballistic helical jet trajectories, jet partial disruption and wind mixing, shocks, and possibly non-thermal emission. Among several observational diagnostics at different bands, the radio morphology on milliarcsecond scales can be particularly insightful regarding the wind-jet interaction.

A mysterious X-ray-emitting object has been detected moving away from the high-mass gamma-ray binary PSR B1259-63, which contains a non-accreting pulsar and a Be star whose winds collide forming a complex interaction structure. Given the strong eccentricity of this binary, the interaction structure should be strongly anisotropic, which together with the complex evolution of the shocked winds, could explain the origin of the observed moving X-ray feature. We propose here that a fast outflow made of a pulsar-stellar wind mixture is always present moving away from the binary in the apastron direction, with the injection of stellar wind occurring at orbital phases close to periastron passage. This outflow periodically loaded with stellar wind would move with a high speed, and likely host non-thermal activity due to shocks, on scales similar to those of the observed moving X-ray object. Such an outflow is thus a very good candidate to explain this X-ray feature. This, if confirmed, would imply pulsar-to-stellar wind thrust ratios of $\sim 0.1$, and the presence of a jet-like structure on the larger scales, up to its termination in the ISM.

(abridged)The structure formed by the shocked winds of a massive star and a non-accreting pulsar in a binary suffers periodic and random variations of orbital and non-linear dynamical origin. For the 1st time, we simulate in 3 D the interaction of isotropic stellar and relativistic pulsar winds along 1 full orbit, on scales well beyond the binary size. We also investigate the impact of grid resolution and size, and of different EoOs: a gamma-constant ideal gas, and an ideal gas with gamma dependent on temperature. We carry out, with the code PLUTO, relativistic HD simulations in 2 and 3 D of the interaction of a slow wind and a relativistic wind with Gamma=2 along 1 full orbit up to ~100 x the binary size. The different 2-D simulations are carried out with equal and larger grid resolution and size, and 1 of them is done with a more realistic equation of state, than in 3 D. The simulations in 3 D confirm previous results in 2 D. The shocked flows are subject to a faster instabilities growth in 3 D, which enhances the presence of shocks, mixing, and large-scale disruption. In 2 D, higher resolution simulations confirm lower resolution results, simulations with larger grid sizes strengthen the case for the loss of global coherence of the shocked-wind structure, and simulations with 2 different EoOs yield very similar results. In addition to the KHI, we find that the Richtmyer-Meshkov and the RTI are likely acting together in the shocked flow evolution. Simulations in 3 D confirm that the interaction of stellar and pulsar winds yields structures that evolve non-linearly and get strongly entangled. The evolution is accompanied by strong kinetic energy dissipation, rapid changes in flow orientation and speed, and turbulent motion. The results strengthen the case for the loss of global coherence of the shocked structure on large scales, although higher pulsar wind speed simulations are needed.

The gamma-ray binary LS I +61 303 is composed of a Be star and a compact companion orbiting in an eccentric orbit. Variable flux modulated with the orbital period of ~26.5 d has been detected from radio to very high-energy gamma rays. In addition, the system presents a superorbital variability of the phase and amplitude of the radio outburst with a period of ~4.6 yr. We present optical photometric observations of LS I +61 303 spanning ~1.5 yr and contemporaneous Halpha equivalent width (EW Halpha) data. The optical photometry shows, for the first time, that the known orbital modulation suffers a positive orbital phase shift and an increase in flux for data obtained 1-yr apart. This behavior is similar to that already known at radio wavelengths, indicating that the optical flux follows the superorbital variability as well. The orbital modulation of the EW Halpha presents the already known superorbital flux variability but shows, also for the first time, a positive orbital phase shift. In addition, the optical photometry exhibits a lag of ~0.1-0.2 in orbital phase with respect to the EW Halpha measurements at similar superorbital phases, and presents a lag of ~0.1 and ~0.3 orbital phases with respect noncontemperaneous radio and X-ray outbursts, respectively. The phase shifts detected in the orbital modulation of thermal indicators, such as the optical flux and the EW Halpha, are in line with the observed behavior for nonthermal indicators, such as X-ray or radio emission. This shows that there is a strong coupling between the thermal and nonthermal emission processes in the gamma-ray binary LS I +61 303. The orbital phase lag between the optical flux and the EW Halpha is naturally explained considering different emitting regions in the circumstellar disk, whereas the secular evolution might be caused by the presence of a moving one-armed spiral density wave in the disk.

Binary systems containing a massive star and a non-accreting pulsar present strong interaction between the stellar and the pulsar winds. The properties of this interaction, which largely determine the non-thermal radiation in these systems, strongly depend on the structure of the stellar wind, which can be clumpy or strongly anisotropic, as in Be stars. We study numerically the influence of inhomogeneities in the stellar wind on the structure of the two-wind interaction region. We carried out for the first time axisymmetric, relativistic hydrodynamical simulations, with Lorentz factors of ~6 and accounting for the impact of instabilities, to study the impact in the two-wind interaction structure of an over-dense region of the stellar wind. We also followed the evolution of this over-dense region or clump as it faces the impact of the pulsar wind. For typical system parameters, and adopting a stellar wind inhomogeneity with a density contrast >~10, clumps with radii of a few percent of the binary size can significantly perturb the two-wind interaction region, pushing the two-wind interface to <~40% of the initial distance to the pulsar. After it is shocked, the inhomogeneity quickly expands and is disrupted when it reaches the smallest distance to the pulsar. It eventually fragments, being advected away from the binary system. The whole interaction region is quite unstable, and the shocked pulsar wind can strongly change under small perturbations. We confirm the sensitive nature of the two-wind interaction structure to perturbations, in particular when the stellar wind is inhomogeneous. For realistic over-dense regions of the stellar wind, the interaction region may shrink by a factor of a few, with the shocked flow presenting a complex spatial and temporal pattern. This can lead to strong variations in the non-thermal radiation.

We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar \c454. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelop lost by the star within the jet, and calculations of the related high-energy radiation. The model readily explains the day-long, variable on timescales of hours, GeV gamma-ray flare from \c454, observed during November 2010 on top of a weeks-long plateau. In the proposed scenario, the plateau state is caused by a strong wind generated by the heating of the star atmosphere by nonthermal particles accelerated at the jet-star interaction region. The flare itself could be produced by a few clouds of matter lost by the red giant after the initial impact of the jet. In the framework of the proposed scenario, the observations constrain the key model parameters of the source, including the mass of the central black hole: $M_{\rm BH}\simeq 10^9 M_{\odot}$, the total jet power: $L_{\rm j}\simeq 10^{48}\,\rm erg\,s^{-1}$, and the Doppler factor of the gamma-ray emitting clouds, $\delta\simeq 20$. Whereas we do not specify the particle acceleration mechanisms, the potential gamma-ray production processes are discussed and compared in the context of the proposed model. We argue that synchrotron radiation of protons has certain advantages compared to other radiation channels of directly accelerated electrons.

Galactic and extragalactic relativistic jets have rich environments that are full of moving objects, such as stars and dense clumps. These objects can enter into the jets and generate shocks and non-thermal emission. We characterize the emitting properties of the downstream region of a standing shock formed due to the interaction of a relativistic jet with an obstacle. We focus on the case of red giants interacting with an extragalactic jet. We perform relativistic axisymmetric hydrodynamical simulations of a relativistic jet meeting an obstacle of very large inertia. The results are interpreted in the framework of a red giant whose dense and slow wind interacts with the jet of an active galactic nucleus. Assuming that particles are accelerated in the standing shock generated in the jet as it impacts the red giant wind, we compute the non-thermal particle distribution, the Doppler boosting enhancement, and the non-thermal luminosity in gamma rays. The available non-thermal energy from jet-obstacle interactions is potentially enhanced by a factor of $\sim 100$ when accounting for the whole surface of the shock induced by the obstacle, instead of just the obstacle section. The observer gamma-ray luminosity, including the flow velocity and Doppler boosting effects, can be ~300(g/10)^2 times higher than when the emitting flow is assumed at rest and only the obstacle section is considered, where g is the jet Lorentz factor. For a whole population of red giants inside the jet of an AGN, the predicted persistent gamma-ray luminosities may be potentially detectable for a jet pointing to the observer. Obstacles interacting with relativistic outflows, for instance clouds and populations of stars for extragalactic jets, or stellar wind inhomogeneities in microquasar jets and in winds of pulsars in binaries, should be taken into account when investigating the non-thermal emission from these sources.

Context: Several binary systems hosting massive stars present gamma-ray emission. In most of these systems, despite detailed observational information being available, the nature and the structure of the emitter are still poorly known. Aims: We investigate the validity of the so-called one-zone approximation for the high-energy emitter in binary systems hosting a massive star. In particular, the case of LS 5039 is considered. Methods: Assuming a point-like emitter at rest, the presence of a nearby massive star, and the observed MeV and GeV fluxes as a reference, a non-thermal leptonic model is systematically applied for different locations, magnetic fields, and non-radiative losses. This allows us to identify both the emitter configurations that are most compatible with observations and inconsistencies between model predictions and the available data. Results: In the case of LS 5039, the best parameter combination is fast non-radiative cooling and a low magnetic field. However, discrepancies appear when comparing the model results at the MeV and GeV energy ranges with the observed fluxes. Predictions fail when the orbital motion is included in the analysis, because emitters and energy budgets that are too large are required. Values of X-ray and TeV fluxes that are too high are predicted in such a case, along half of the orbit. Conclusions: We show that the radiation in LS 5039 does not come from only one electron population, and the emitter is likely extended and inhomogeneous with a low magnetic field. We suggest that the emitter moves at relativistic velocities with Doppler boosting playing a significant role. Key words. Stars: high-mass, X-rays: binaries, gamma-rays: binaries

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

Relativistic stellar jets, produced in binary systems called microquasars, propagate through media with different spatial scales releasing their energy in the form of work and radiation from radio to gamma rays. There are several medium-interaction scenarios that these jets can face. In particular, in relativistic stellar jets the presence of a star is an unavoidable element whose importance deserves to be studied. In the case of high-mass stars, their powerful winds are likely to interact dynamically with the jet, but also low-mass stars in the post-main sequence phase can present dense winds that will act as an obstacle for the jet propagation. In this work, we present a semi-qualitative discussion on the importance of the star for the evolution of relativistic stellar jets.

The Fermi bubbles are part of a complex region of the Milky Way. This region presents broadband extended non-thermal radiation, apparently coming from a physical structure rooted in the Galactic Centre and with a partly-ordered magnetic field threading it. We explore the possibility of an explosive origin for the Fermi bubble region to explain its morphology, in particular that of the large-scale magnetic fields, and provide context for the broadband non-thermal radiation. We perform 3D magnetohydrodynamical simulations of an explosion from a few million years ago that pushed and sheared a surrounding magnetic loop, anchored in the molecular torus around the Galactic Centre. Our results can explain the formation of the large-scale magnetic structure in the Fermi bubble region. Consecutive explosive events may match better the morphology of the region. Faster velocities at the top of the shocks than at their sides may explain the hardening with distance from the Galactic Plane found in the GeV emission. In the framework of our scenario, we estimate the lifetime of the Fermi bubbles as $2\times10^6$ yr, with a total energy injected in the explosion(s) $> 10^{55}$ ergs. The broadband non-thermal radiation from the region may be explained by leptonic emission, more extended in radio and X-rays, and confined to the Fermi bubbles in gamma rays.

High-mass binaries hosting young pulsars can be powerful gamma-ray emitters. The stellar wind of the massive star in the system is expected to be clumpy. Since the high-energy emission comes from the pulsar-star wind interaction, the presence of clumps can affect the spectrum and variability of this radiation. We look for the main effects of the clumps on the two-wind interaction region and on the non-thermal radiation. A simple analytical model for the two-wind interaction dynamics was developed accounting for the lifetime of clumps under the pulsar-wind impact. This time plays a very important role with regard to the evolution of the clump, the magnetic field in the clump-pulsar wind interaction region, and the non-radiative and radiative cooling of the non-thermal particles. We also computed the high-energy emission produced at the interaction of long-living clumps with the pulsar wind. For reasonable parameters, the clumps will induce small variability on the X-ray and gamma-ray radiation. Sporadically, large clumps can reach closer to the pulsar increasing the magnetic field, triggering synchrotron X-ray flares and weakening other emission components like inverse Compton. The reduction of the emitter size induced by clumps also makes non-radiative losses faster. Stellar wind clumps can also enhance instability development and matter entrainment in the shocked pulsar wind when it leaves the binary. Growth limitations of the clumps from the wind acceleration region may imply that a different origin for the largest clumps is required. The large-scale wind structures behind the observed discrete absorption components in the UV may be the source of these large clumps. The presence of structure in the stellar wind can produce substantial energy-dependent variability and should not be neglected when studying the broadband emission from high-mass binaries hosting young pulsars.

Dense populations of stars surround the nuclear regions of galaxies. In active galactic nuclei, these stars can interact with the relativistic jets launched by the supermasive black hole. In this work, we study the interaction of early-type stars with relativistic jets in active galactic nuclei. A bow-shaped double-shock structure is formed as a consequence of the interaction of the jet and the stellar wind of each early-type star. Particles can be accelerated up to relativistic energies in these shocks and emit high-energy radiation. We compute, considering different stellar densities of the galactic core, the gamma-ray emission produced by non-thermal radiative processes. This radiation may be significant in some cases, and its detection might yield valuable information on the properties of the stellar population in the galaxy nucleus, as well as on the relativistic jet. This emission is expected to be particularly relevant for nearby non-blazar sources.

In the recent years massive protostars have been suggested to be high-energy emitters. Among the best candidates is IRAS 16547-4247, a protostar that presents a powerful outflow with clear signatures of interaction with its environment. This source has been revealed to be a potential high-energy source because it displays non-thermal radio emission of synchrotron origin, which is evidence of relativistic particles. To improve our understanding of IRAS 16547-4247 as a high-energy source, we analyzed XMM-Newton archival data and found that IRAS 16547-4247 is a hard X-ray source. We discuss these results in the context of a refined one-zone model and previous radio observations. From our study we find that it may be difficult to explain the X-ray emission as non-thermal radiation coming from the interaction region, but it might be produced by thermal Bremsstrahlung (plus photo-electric absorption) by a fast shock at the jet end. In the high-energy range, the source might be detectable by the present generation of Cherenkov telescopes, and may eventually be detected by Fermi in the GeV range.

We study the interaction of early-type stars with the jets of active galactic nuclei. A bow-shock will form as a consequence of the interaction of the jet with the winds of stars and particles can be accelerated up to relativistic energies in these shocks. We compute the non-thermal radiation produced by relativistic electrons from radio to gamma-rays. This radiation may be significant, and its detection might yield information on the properties of the stellar population in the galaxy nucleus, as well as on the relativistic jet. This emission is expected to be relevant for nearby non-blazar sources.

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

We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar 3C454.3. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelop lost by the star within the jet, and calculations of the related high-energy radiation. The model readily explains the day-long, variable on timescales of hours, GeV gamma-ray flare from 3C454.3, observed during November 2010 on top of a weeks-long plateau. In the proposed scenario, the plateau state is caused by a strong wind generated by the heating of the star atmosphere by nonthermal particles accelerated at the jet-star interaction region. The flare itself could be produced by a few clouds of matter lost by the red giant after the initial impact of the jet. In the framework of the proposed scenario, the observations constrain the key model parameters of the source, including the mass of the central black hole: $M_{\rm BH}\simeq 10^9 M_{\odot}$, the total jet power: $L_{\rm j}\simeq 10^{48}\,\rm erg\,s^{-1}$, and the Doppler factor of the gamma-ray emitting clouds, $\delta\simeq 20$. Whereas we do not specify the particle acceleration mechanisms, the potential gamma-ray production processes are discussed and compared in the context of the proposed model. We argue that synchrotron radiation of protons has certain advantages compared to other radiation channels of directly accelerated electrons.

The high and very-high energy spectrum of gamma-ray binaries has become a challenge for all theoretical explanations since the detection of powerful, persistent GeV emission from LS 5039 and LS I +61 303 by Fermi/LAT. The spectral cutoff at a few GeV indicates that the GeV component and the fainter, hard TeV emission above 100 GeV are not directly related. We explore the possible origins of these two emission components in the framework of a young, non-accreting pulsar orbiting the massive star, and initiating the non-thermal emission through the interaction of the stellar and pulsar winds. The pulsar/stellar wind interaction in a compact orbit binary gives rise to two potential locations for particle acceleration: the shocks at the head-on collision of the winds and the termination shock caused by Coriolis forces on scales larger than the binary separation. We explore the suitability of these two locations to host the GeV and TeV emitters, respectively, through the study of their non-thermal emission along the orbit. We focus on the application of this model to LS 5039 given its well determined stellar wind with respect to other gamma-ray binaries. The application of the proposed model to LS 5039 indicates that these two potential emitter locations provide the necessary conditions for reproduction of the two-component high-energy gamma-ray spectrum of LS 5039. In addition, the ambient postshock conditions required at each of the locations are consistent with recent hydrodynamical simulations. The scenario based on the interaction of the stellar and pulsar winds is compatible with the GeV and TeV emission observed from gamma-ray binaries with unknown compact objects, such as LS 5039 and LS I +61 303.

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

Dense populations of stars surround the nuclear regions of galaxies. In this work, we study the interaction of a WR star with relativistic jets in active galactic nuclei. A bow-shaped double-shock structure will form as a consequence of the interaction of the jet and the wind of the star. Particles can be accelerated up to relativistic energies in these shocks and emit high-energy radiation. We compute the produced gamma-ray emission obtaining that this radiation may be significant. This emission is expected to be particularly relevant for nearby non-blazar sources.

The winds from a non-accreting pulsar and a massive star in a binary system collide forming a bow-shaped shock structure. The Coriolis force induced by orbital motion deflects the shocked flows, strongly affecting their dynamics. We study the evolution of the shocked stellar and pulsar winds on scales in which the orbital motion is important. Potential sites of non-thermal activity are investigated. Relativistic hydrodynamical simulations in two dimensions, performed with the code PLUTO and using the adaptive mesh refinement technique, are used to model interacting stellar and pulsar winds on scales ~80 times the distance between the stars. The hydrodynamical results suggest the suitable locations of sites for particle acceleration and non-thermal emission. In addition to the shock formed towards the star, the shocked and unshocked components of the pulsar wind flowing away from the star terminate by means of additional strong shocks produced by the orbital motion. Strong instabilities lead to the development of turbulence and an effective two-wind mixing in both the leading and trailing sides of the interaction structure, which starts to merge with itself after one orbit. The adopted moderate pulsar-wind Lorentz factor already provides a good qualitative description of the phenomena involved in high-mass binaries with pulsars, and can capture important physical effects that would not appear in non-relativistic treatments. Simulations show that shocks, instabilities, and mass-loading yield efficient mass, momentum, and energy exchanges between the pulsar and the stellar winds. This renders a rapid increase in the entropy of the shocked structure, which will likely be disrupted on scales beyond the simulated ones. Several sites of particle acceleration and low- and high-energy emission can be identified. Doppler boosting will have significant and complex effects on radiation.

Active galactic nuclei with misaligned jets have been recently established as a class of high-energy gamma-ray sources. M87, a nearby representative of this class, shows fast TeV variability on timescales less than one day. We present calculations performed in the framework of the scenario in which gamma-ray flares in non-blazar active galactic nuclei are produced by a red giant or a gas cloud interacting with the jet. We show that both the light curve and energy spectrum of the spectacular April 2010 flare can be reproduced by this model, assuming that a relatively massive cloud of approx 1.e29 g penetrates into the jet at few tens of Schwarzschild radii from the super-massive black hole.

Context. Converging flows with strong magnetic fields of different polarity can accelerate particles through magnetic reconnection. If the particle mean free path is longer than the reconnection layer is thick, but much shorter than the entire reconnection structure, the particle will mostly interact with the incoming flows potentially with a very low escape probability. Aims. We explore, in general and also in some specific scenarios, the possibility of particles to be accelerated in a magnetic reconnection layer by interacting only with the incoming flows. Methods. We characterize converging flows that undergo magnetic reconnection, and derive analytical estimates for the particle energy distribution, acceleration rate, and maximum energies achievable in these flows. We also discuss a scenario, based on jets dominated by magnetic fields of changing polarity, in which this mechanism may operate. Results. The proposed acceleration mechanism operates if the reconnection layer is much thinner than its transversal characteristic size, and the magnetic field has a disordered component. Synchrotron losses may prevent electrons from entering in this acceleration regime. The acceleration rate should be faster, and the energy distribution of particles harder than in standard diffusive shock acceleration. The interaction of obstacles with the innermost region of jets in active galactic nuclei and microquasars may be suitable sites for particle acceleration in converging flows.

The very high energy (VHE) \gamma-ray source HESS J0632+057 has recently been confirmed to be a \gamma-ray binary. The optical counterpart is the Be star MWC 148, and a compact object of unknown nature orbits it every ~321 d with a high eccentricity of ~0.8. We monitored HESS J0632+057 with the stereoscopic MAGIC telescopes from 2010 October to 2011 March and detected significant VHE \gamma-ray emission during 2011 February, when the system exhibited an X-ray outburst. We find no \gamma-ray signal in the other observation periods when the system did not show increased X-ray flux. Thus HESS J0632+057 exhibits \gamma-ray variability on timescales of the order of one to two months possibly linked to the X-ray outburst that takes place about 100 days after the periastron passage. Furthermore our measurements provide for the first time the \gamma-ray spectrum down to about 140 GeV and indicate no turnover of the spectrum at low energies. We compare the properties of HESS J0632+057 with the similar \gamma-ray binary LS I +61 303, and discuss on the possible origin of the multi-wavelength emission of the source

Extragalactic jets are formed close to supermassive black-holes in the center of galaxies. Large amounts of gas, dust, and stars cluster in the galaxy nucleus, and interactions between this ambient material and the jet base should be frequent, having dynamical as well as radiative consequences. This work studies the dynamical interaction of an obstacle, a clump of matter or the atmosphere of an evolved star, with the innermost region of an extragalactic jet. Jet mass-loading and the high-energy outcome of this interaction are briefly discussed. Relativistic hydrodynamical simulations with axial symmetry have been carried out for homogeneous and inhomogeneous obstacles inside a relativistic jet. These obstacles may represent a medium inhomogeneity or the disrupted atmosphere of a red giant star. Once inside the jet, an homogeneous obstacle expands and gets disrupted after few dynamical timescales, whereas in the inhomogeneous case, a solid core can smoothen the process, with the obstacle mass-loss dominated by a dense and narrow tail pointing in the direction of the jet. In either case, matter is expected to accelerate and eventually get incorporated to the jet. Particles can be accelerated in the interaction region, and produce variable gamma-rays in the ambient matter, magnetic and photon fields. The presence of matter clumps or red giants into the base of an extragalactic jet likely implies significant jet mass-loading and slowing down. Fast flare-like gamma-ray events, and some level of persistent emission, are expected due to these interactions.

The termination structures of the jets of Fanaroff & Riley galaxies are observed to produce extended non-thermal emission in a wide frequency range. The study of these structures can provide valuable insights on the conditions for particle acceleration and radiation at the shock fronts.We have studied the thermal and non-thermal emission that can be expected from the jet termination regions of Fanaroff & Riley type I sources. The broadband emission from these galaxies has been recently extended to include the high-energy gamma-ray domain, owing to the Fermi detection of Cen A lobes. Exploring the physics behind the jet/medium interactions in FRI can provide valuable insights on the conditions for particle acceleration and radiation in the jet termination shocks. Making use of the results of a fully relativistic numerical simulation code of the evolution of a FRI jet we model the expected radiative output and predict spectra and lightcurves of both thermal and non-thermal emission at different source ages.

High-mass microquasars consist of a massive star and a compact object, the latter producing jets that will interact with the stellar wind. The evolution of the jets, and ultimately their radiative outcome, could depend strongly on the inhomogeneity of the wind, which calls for a detailed study. The hydrodynamics of the interaction between a jet and a clumpy wind is studied, focusing on the global wind and single clump-jet interplay. We have performed, using the code \textit{Ratpenat}, three-dimensional numerical simulations of a clumpy wind interacting with a mildly relativistic jet, and of individual clumps penetrating into a jet. For typical wind and jet velocities, filling factors of about > 0.1 are already enough for the wind to be considered as clumpy. An inhomogeneous wind makes the jet more unstable when crossing the system. Kinetic luminosities of the order 1.e37 erg/s allow the jet to reach the borders of a compact binary with an O star, as in the smooth wind case, although with a substantially higher degree of disruption. When able to enter into the jet, clumps are compressed and heated during a time of about their size divided by the sound speed in the shocked clump. Then, clumps quickly disrupt, mass-loading and slowing down the jet. We conclude that moderate wind clumpiness makes already a strong difference with the homogeneous wind case, enhancing jet disruption, mass-loading, bending, and likely energy dissipation in the form of emission. All this can have observational consequences at high-energies and also in the large scale radio jets.

We present very high energy (VHE, E > 100 GeV) {\gamma}-ray observations of the {\gamma}-ray binary system LS I+61 303 obtained with the MAGIC stereo system between 2009 October and 2010 January. We detect a 6.3{\sigma} {\gamma}-ray signal above 400 GeV in the combined data set. The integral flux above an energy of 300 GeV is F(E>300 GeV)=(1.4 +- 0.3stat +- 0.4syst) * 10^{-12} cm^{-2} s^{-1}, which corresponds to about 1.3% of the Crab Nebula flux in the same energy range. The orbit-averaged flux of LS I +61 303 in the orbital phase interval 0.6--0.7, where a maximum of the TeV flux is expected, is lower by almost an order of magnitude compared to our previous measurements between 2005 September and 2008 January. This provides evidence for a new low emission state in LS I +61 303. We find that the change to the low flux state cannot be solely explained by an increase of photon-photon absorption around the compact star.

Gamma-ray loud X-ray binaries are binary systems that show non-thermal broadband emission from radio to gamma rays. If the system comprises a massive star and a young non-accreting pulsar, their winds will collide producing broadband non-thermal emission, most likely originated in the shocked pulsar wind. Thermal X-ray emission is expected from the shocked stellar wind, but until now it has neither been detected nor studied in the context of gamma-ray binaries. We present a semi-analytic model of the thermal X-ray emission from the shocked stellar wind in pulsar gamma-ray binaries, and find that the thermal X-ray emission increases monotonically with the pulsar spin-down luminosity, reaching luminosities of the order of 10^33 erg/s. The lack of thermal features in the X-ray spectrum of gamma-ray binaries can then be used to constrain the properties of the pulsar and stellar winds. By fitting the observed X-ray spectra of gamma-ray binaries with a source model composed of an absorbed non-thermal power law and the computed thermal X-ray emission, we are able to derive upper limits on the spin-down luminosity of the putative pulsar. We applied this method to LS 5039, the only gamma-ray binary with a radial, powerful wind, and obtain an upper limit on the pulsar spin-down luminosity of ~6x10^36 erg/s. Given the energetic constraints from its high-energy gamma-ray emission, a non-thermal to spin-down luminosity ratio very close to unity may be required.

Several gamma-ray binaries show extended X-ray emission that may be associated to interactions of an outflow with the medium. Some of these systems are, or may be, high-mass binaries harboring young nonaccreting pulsars, in which the stellar and the pulsar winds collide, generating a powerful outflow that should terminate at some point in the ambient medium. This work studies the evolution and termination, as well as the related radiation, of the shocked-wind flow generated in high-mass binaries hosting powerful pulsars. A characterization, based on previous numerical work, is given for the stellar/pulsar wind interaction. Then, an analytical study of the further evolution of the shocked flow and its dynamical impact on the surrounding medium is carried out. Finally, the expected nonthermal emission from the flow termination shock, likely the dominant emitting region, is calculated. The shocked wind structure, initially strongly asymmetric, becomes a quasi-spherical, supersonically expanding bubble, with its energy coming from the pulsar and mass from the stellar wind. This bubble eventually interacts with the environment on ~pc scales, producing a reverse and, sometimes, a forward shock. Nonthermal leptonic radiation can be efficient in the reverse shock. Radio emission is expected to be faint, whereas X-rays can easily reach detectable fluxes. Under very low magnetic fields and large nonthermal luminosities, gamma rays may also be significant. We conclude that the complexity of the stellar/pulsar wind interaction is likely to be smoothed out outside the binary system, where the wind-mixed flow accelerates and eventually terminates in a strong reverse shock. This shock may be behind the extended X-rays observed in some binary systems. For very powerful pulsars, part of the unshocked pulsar wind may directly interact with the large-scale environment.

Binary systems can be powerful sources of non-thermal emission from radio to gamma rays. When the latter are detected, then these objects are known as gamma-ray binaries. In this work, we explore, in the context of gamma-ray binaries, different acceleration processes to estimate their efficiency: Fermi I, Fermi II, shear acceleration, the converter mechanism, and magnetic reconnection. We find that Fermi I acceleration in a mildly relativistic shock can provide, although marginally, the multi-10 TeV particles required to explain observations. Shear acceleration may be a complementary mechanism, giving particles the final boost to reach such a high energies. Fermi II acceleration may be too slow to account for the observed very high energy photons, but may be suitable to explain extended low-energy emission. The converter mechanism seems to require rather high Lorentz factors but cannot be discarded a priori. Standard relativistic shock acceleration requires a highly turbulent, weakly magnetized downstream medium; magnetic reconnection, by itself possibly insufficient to reach very high energies, could perhaps facilitate such a conditions. Further theoretical developments, and a better source characterization, are needed to pinpoint the dominant acceleration mechanism, which need not be one and the same in all sources.

Presently there are several classes of detected gamma-ray extragalatic sources. They are mostly associated to active galactic nuclei (AGN) and (at soft gamma rays) to gamma-ray bursts (GRB), but not only. Active galactic nuclei consist of accreting supermassive black holes hosted by a galaxy that present in some cases powerful relativistic jet activity. These sources, which have been studied in gamma rays for several decades, are probably the most energetic astrophysical objects, and their appearance depends much on whether their jets point to us. Gamma-ray bursts, thought to be associated to collapsing or merging stellar-mass objects at cosmological distances, are also accreting highly relativistic jet sources that shine strongly at high energies. These are very short-duration events, but they are also the most luminous. Recently, star formation galaxies have turned out to be also gamma-ray emitters. On the other hand, clusters of galaxies have not been detected beyond X-rays yet. These are the largest known structures in the Universe; in their formation through accretion and merging, shocks and turbulence are generated, which may lead to gamma-ray production. In this work, the gamma-ray physics of AGNs is briefly presented, as well as that of starburst galaxies, GRBs and clusters of galaxies. Afterwards, we consider some particular cases of gamma-ray production in non-blazar AGN jets interacting with their medium at different scales.

Jets are ubiquitous in the Universe. They are collimated outflows whose origin is associated to an accretion disc and a central object, and can be very powerful non-thermal emitters. Jets form in active galactic nuclei, gamma-ray bursts, microquasars, and young stellar objects. Galactic jets emitting non-thermal emission are typically associated to microquasars, although the jets of massive young stellar objects are also non-thermal sources. The production of non-thermal radiation, in particular radio synchrotron emission, is a clear indication that particle acceleration is taking place in the source, which hints to the generation of photons even at high energies. In this work, we will discuss the emitting sites in, or related to, microquasar jets, and briefly comment on the possibility of high-energy emission in jets from young stellar objects.

High energy emission can be produced in the interaction sites of both galactic and extragalactic jets with the surrounding medium. We have developed a radiative model that accounts for the continuous injection of relativistic electrons in the forward, reverse and recollimation shocks developed in the shell, cocoon and reconfinement interaction regions, respectively. We also performed hydrodynamical simulations to establish the physical properties in both galactic and extragalactic systems. The resulting non-thermal emission is studied assuming different values for the jet power, the external mass density and the source age for both FR I galaxies and galactic microquasars. The obtained fluxes are compared to current instrument sensitivities at radio, X-ray and gamma-ray bands.

Microquasars are X-ray binaries that present non-thermal radio jets. Efficient particle acceleration can take place in different regions of the jets of microquasars. The accelerated particles can emit gamma-rays via leptonic or hadronic processes, with a complex spectral and temporal behavior. The jet termination region can be also an efficient non-thermal emitter, as well as, in high-mass microquasars, the region of the binary system outside the jet. In this work, I briefly describe the physics behind the non-thermal emission observed in microquasars at different scales, focusing in the GeV and TeV bands.

The acceleration of particles up to GeV or higher energies in microquasars has been the subject of considerable theoretical and observational efforts in the past few years. Sco X-1 is a microquasar from which evidence of highly energetic particles in the jet has been found when it is in the so-called Horizontal Branch (HB), a state when the radio and hard X-ray fluxes are higher and a powerful relativistic jet is present. Here we present the first very high energy gamma-ray observations of Sco X-1 obtained with the MAGIC telescopes. An analysis of the whole dataset does not yield a significant signal, with 95% CL flux upper limits above 300 GeV at the level of 2.4x10^{-12} ph/cm^2/s. Simultaneous RXTE observations were conducted to search for TeV emission during particular X-ray states of the source. A selection of the gamma-ray data obtained during the HB based on the X-ray colors did not yield a signal either, with an upper limit of 3.4x10^{-12} ph/cm^2/s. These upper limits place a constraint on the maximum TeV luminosity to non-thermal X-ray luminosity of L_{VHE}/L_{ntX}<0.02, that can be related to a maximum TeV luminosity to jet power ratio of L_{VHE}/L_{j}<10^{-3}. Our upper limits indicate that the underlying high-energy emission physics in Sco X-1 must be inherently different from that of the hitherto detected gamma-ray binaries.

Several binary systems that contain a massive star have been detected in both the radio band and at very high energies. In the dense stellar photon field of these sources, gamma-ray absorption and pair creation are expected to occur, and the radiation from these pairs may contribute significantly to the observed radio emission. We aim at going deeper in the study of the properties, and in particular the morphology, of the pair radio emission in gamma-ray binaries. We apply a Monte-Carlo code that computes the creation location, the spatial trajectory and the energy evolution of the pairs produced in the binary system and its surroundings. The radio emission produced by these pairs, with its spectral, variability and spatial characteristics, is calculated as it would be seen from a certain direction. A generic case is studied first, and then the specific case of LS 5039 is also considered. We find that, confirming previous results, the secondary radio emission should appear as an extended radio structure of a few milliarcseconds size. This radiation would be relatively hard, with fluxes up to ~ 10 mJy. Modulation is expected depending on the gamma-ray production luminosity, system eccentricity, and wind ionization fraction, and to a lesser extent on the magnetic field structure. In gamma-ray binaries in general, the pairs created due to photon-photon interactions can contribute significantly to the core, and generate an extended structure. In the case of LS 5039, the secondary radio emission is likely to be a significant fraction of the detected core flux, with a marginal extension.

The high-mass microquasar Cygnus X-3 has been recently detected in a flaring state by the gamma-ray satellites Fermi and Agile. In the present contribution, we study the high-energy emission from Cygnus X-3 through a model based on the interaction of clumps from the Wolf-Rayet wind with the jet. The clumps inside the jet act as obstacles in which shocks are formed leading to particle acceleration and non-thermal emission. We model the high energy emission produced by the interaction of one clump with the jet and briefly discus the possibility of many clumps interacting with the jet. From the characteristics of the considered scenario, the produced emission could be flare-like due to discontinuous clump penetration, with the GeV long-term activity explained by changes in the wind properties.

We report on the observation of the region around supernova remnant G65.1+0.6 with the stand-alone MAGIC-I telescope. This region hosts the two bright GeV gamma-ray sources 1FGL J1954.3+2836 and 1FGL J1958.6+2845. They are identified as GeV pulsars and both have a possible counterpart detected at about 35 TeV by the Milagro observatory. MAGIC collected 25.5 hours of good quality data, and found no significant emission in the range around 1 TeV. We therefore report differential flux upper limits, assuming the emission to be point-like (<0.1 deg) or within a radius of 0.3 deg. In the point-like scenario, the flux limits around 1 TeV are at the level of 3 % and 2 % of the Crab Nebula flux, for the two sources respectively. This implies that the Milagro emission is either extended over a much larger area than our point spread function, or it must be peaked at energies beyond 1 TeV, resulting in a photon index harder than 2.2 in the TeV band.

Binary systems that harbor a non-accreting pulsar are efficient non-thermal emitters, from radio to gamma rays. This broadband emission is thought to come from the region where the companion star and pulsar winds collide. A paradigmatic example of this source type is PSR B1259$-$63. Whereas the high-energy radiation probes the shock structure at the binary scale, the radio emission is produced well outside the system, in regions where the shocked stellar and pulsar winds are likely mixing due to hydrodynamical instabilities. Understanding the evolution of the shocked flow depends strongly on a proper characterization of the low-energy radiation. We have performed numerical calculations of the radio emission produced in a high-mass binary hosting a young pulsar. Adopting a prescription for the shocked flow in the system and the non-thermal particle injection, we have generated synthetic radio emission maps that can be compared with observations. Preliminary results suggest that the emitting bulk flow should be rather slow, with a pulsar wind magnetization at the shock higher than in isolated pulsars.

The environment of high-mass X-ray binaries can be characterized either by the SNR that forms these systems, or by the wind from the companion massive star. These regions should be tenuous but very hot, and surrounded by a dense and cold shocked ISM shell. The interaction between the jet and such a complex medium, also affected by the system proper motion, can lead to very different jet termination structures. The evolution of the jet termination regions during the life of a high-mass microquasar is simulated to improve the present understanding of these structures. Also, the evolving emission characteristics are modeled to inform potential observational campaigns for this class of object. We have performed 2D numerical simulations of jets propagating in different scenarios, corresponding to different epochs after the formation of the high-mass X-ray binary, using the code Ratpenat. We have also made simple estimates of the non-thermal emission that could be produced in the jet termination regions. We find that, in the way through the hot and tenuous medium of the shocked wind/SNR ejecta, the jet suffers recollimation shocks in which it loses part of its thrust and ends in a strong shock inflating a hot cocoon. The jet head propagates with a speed similar to the medium sound speed, until it eventually reaches the denser and colder shocked ISM and the unperturbed ISM later on. In these last stages of evolution, the jet is significantly slowed down and can be disrupted. For relatively old sources, the microquasar peculiar velocity becomes important, leading to complete jet destruction. Extended non-thermal radiation can be generated in the jet termination regions, being hard X-rays and TeV photons the best suited wavelengths to observe these structures.

The MAGIC collaboration has recently reported correlated X-ray and very high-energy gamma-ray emission from the gamma-ray binary LS I +61 303 during ~60% of one orbit. These observations suggest that the emission in these two bands has its origin in a single particle population. We aim at improving our understanding of the source behaviour by explaining the simultaneous X-ray and VHE data through a radiation model. We use a model based on a one zone population of relativistic leptonic particles assuming dominant adiabatic losses located at the position of the compact object. The adiabatic cooling timescale is inferred from the X-ray fluxes. The model can reproduce the spectra and lightcurves in the X-ray and VHE bands. Adiabatic losses could be the key ingredient to explain the X-ray and partially the VHE lightcurves. From the best fit result, we obtain a magnetic field of B=0.2 G, a minimum luminosity budget of ~2x10^35 erg/s and a relatively high acceleration efficiency. In addition, our results seem to confirm that the GeV emission detected by Fermi does not come from the same parent particle population as the X-ray and VHE emission and the Fermi spectrum poses a constraint on the hardness of the particle spectrum at lower energies. In the context of our scenario, more sensitive observations would allow to constrain the inclination angle, which could determine the nature of the compact object.

Between 2004 and 2009 a sample of 28 X-ray selected high- and intermediate-frequency peaked blazars with a X-ray flux larger than 2 uJy at 1 keV in the redshift range from 0.018 to 0.361 was observed with the MAGIC telescope at energies above 100 GeV. Seven among them were detected and the results of these observations are discussed elsewhere. Here we concentrate on the remaining 21 blazars which were not detected during this observation campaign and present the 3 sigma (99.7 %) confidence upper limits on their flux. The individual flux upper limits lie between 1.6 % and 13.6 % of the integral flux from the Crab Nebula. Applying a stacking method to the sample of non-detections with a total of 394.1 hours exposure time, we find evidence for an excess with a cumulative significance of 4.9 standard deviations. It is not dominated by individual objects or flares, but increases linearly with the observation time as for a constant source with an integral flux level of ~1.5 % of that observed from the Crab Nebula above 150 GeV.

We report new observations of the intermediate-frequency peaked BL Lacertae object 3C 66A with the MAGIC telescopes. The data sample we use were taken in 2009 December and 2010 January, and comprises 2.3 hr of good quality data in stereoscopic mode. In this period, we find a significant signal from the direction of the blazar 3C 66A. The new MAGIC stereoscopic system is shown to play an essential role for the separation between 3C 66A and the nearby radio galaxy 3C 66B, which is at a distance of only $6^\prime$. The derived integral flux above $100\eh{GeV}$ is 8.3% of Crab Nebula flux and the energy spectrum is reproduced by a power law of photon index $3.64 \pm 0.39_{\rm stat} \pm 0.25_{\rm sys}$. Within errors, this is compatible with the one derived by VERITAS in 2009. From the spectra corrected for absorption by the extragalactic background light, we only find small differences between the four models that we applied, and constrain the redshift of the blazar to $z < 0.68$.

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied history of extreme spectral variability and is an excellent laboratory for studying the physical processes within the jets of active galactic nuclei. However, there are few detailed multiwavelength studies of Markarian 501 during its quiescent state, due to its low luminosity. A short-term multiwavelength study of Markarian 501 was coordinated in March 2009, focusing around a multi-day observation with the Suzaku X-ray satellite and including {\gamma}-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray Space Telescope with the goal of providing a well-sampled multiwavelength baseline measurement of Markarian 501 in the quiescent state. The results of these quiescent-state observations are compared to the historically extreme outburst of April 16, 1997, with the goal of examining variability of the spectral energy distribution between the two states. The derived broadband spectral energy distribution shows the characteristic double-peaked profile. We find that the X-ray peak shifts by over two orders of magnitude in photon energy between the two flux states while the VHE peak varies little. The limited shift in the VHE peak can be explained by the transition to the Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and the implied Klein-Nishina effects are explored.

Jets are found in a variety of astrophysical sources, from young stellar objects to active galactic nuclei. In all the cases the jet propagates with a supersonic velocity through the external medium, which can be inhomogeneous, and inhomogeneities could penetrate into the jet. The interaction of the jet material with an obstacle produces a bow shock in the jet in which particles can be accelerated up to relativistic energies and emit high-energy photons. In this work, we explore the active galactic nuclei scenario, focusing on the dynamical and radiative consequences of the interaction at different jet heights. We find that the produced high-energy emission could be detectable by the current gamma-ray telescopes. In general, the jet-clump interactions are a possible mechanism to produce (steady or flaring) high-energy emission in many astrophysical sources in which jets are present.

We investigate the production of gamma-rays in the inner regions of SS 433 and in its interaction between its jets and the W50 nebula. We estimate the VHE emission that can be generated within the jets at distances above 1e+13 cm from the compact object. We also apply a theoretical model of the jet/medium interaction to SS 433/W50. We compare the predicted fluxes to observations of SS 433 at TeV energies, and derive new constraints of the physical properties of this system.

Recent observations in X-rays and gamma-rays of nearby FRI radio galaxies have raised the question of the origin of the emission detected in the termination structures of their jets. The study of these structures can give information on the conditions for particle acceleration and radiation at the front shocks. In addition, an evolutionary scenario can help to disentangle the origin of the detected X-ray emission in young FRI sources, like some Gigahertz Peaked Spectrum AGNs. This work focuses on the nature and detectability of the radiation seen from the termination regions of evolving FRI jets. We use the results of a relativistic, two-dimensional numerical simulation of the propagation of an FRI jet, coupled with a radiation model, to make predictions for the spectra and lightcurves of the thermal and non-thermal emission at different stages of the FRI evolution. Our results show that under moderate magnetic fields, the synchrotron radiation would be the dominant non-thermal channel, appearing extended in radio and more compact in X-rays, with relatively small flux variations with time. The shocked jet synchrotron emission would dominate the X-ray band, although the shocked ISM/ICM thermal component alone may be significant in old sources. Inverse Compton scattering of CMB photons could yield significant fluxes in the GeV and TeV bands, with a non-negligible X-ray contribution. The IC radiation would present a bigger angular size in X-rays and GeV than in TeV, with fluxes increasing with time. We conclude that the thermal and non-thermal broadband emission from the termination regions of FRI jets could be detectable for sources located up to distances of a few 100 Mpc.

We present the first three-dimensional simulations of the evolution of a microquasar jet inside the binary-star system. The aim is to study the interaction of these jets with the stellar wind from a massive companion and the possible locations of high-energy emission sites. We have simulated two jets with different injection power in order to give a hint on the minimum power required for the jet to escape the system and become visible in larger scales. In the setup, we include a massive star wind filling the grid through which the jet evolves. We show that jets should have powers of the order of $10^{37}\rm{erg/s}$ or more in order not to be destroyed by the stellar wind. The jet-wind interaction results in regions in which high energy emission could be produced. These results imply the possible existence of a population of X-ray binaries not detected in the radio band due to jet disruption inside the region dominated by the stellar wind.

We report on the detection with the MAGIC telescopes of very high energy gamma-rays from IC 310, a head-tail radio galaxy in the Perseus galaxy cluster, observed during the interval November 2008 to February 2010. The Fermi satellite has also detected this galaxy. The source is detected by MAGIC at a high statistical significance of 7.6sigma in 20.6 hr of stereo data. The observed spectral energy distribution is flat with a differential spectral index of -2.00 \pm 0.14. The mean flux above 300 GeV, between October 2009 and February 2010, (3.1 \pm 0.5)x10^{-12} cm^{-2} s^{-1}, corresponds to (2.5 \pm 0.4)% of Crab Nebula units. Only an upper limit, of 1.9% of Crab Nebula units above 300 GeV, was obtained with the 2008 data. This, together with strong hints (>3sigma) of flares in the middle of October and November 2009, implies that the emission is variable. The MAGIC results favour a scenario with the very high energy emission originating from the inner jet close to the central engine. More complicated models than a simple one-zone SSC scenario, e.g. multi-zone SSC, external Compton or hadronic, may be required to explain the very flat spectrum and its extension over more than three orders of magnitude in energy.

Non-blazar AGN have been recently established as a class of gamma-ray sources. M87, a nearby representative of this class, show fast TeV variability on timescales of a few days. We suggest a scenario of flare gamma-ray emission in non-blazar AGN based on a red giant interacting with the jet at the base. We solve the hydrodynamical equations that describe the evolution of the envelope of a red giant blown by the impact of the jet. If the red giant is at least slightly tidally disrupted by the supermassive black hole, enough stellar material will be blown by the jet, expanding quickly until a significant part of the jet is shocked. This process can render suitable conditions for energy dissipation and proton acceleration, which could explain the detected day-scale TeV flares from M87 via proton-proton collisions. Since the produced radiation would be unbeamed, such an events should be mostly detected from non-blazar AGN. They may be frequent phenomena, detectable in the GeV-TeV range even up to distances of $\sim 1$ Gpc for the most powerful jets. The counterparts at lower energies are expected to be not too bright.} {M87, and nearby non-blazar AGN in general, can be fast variable sources of gamma-rays through red giant/jet interactions.

Cygnus X-3 is a microquasar consisting of an accreting compact object orbiting around a Wolf-Rayet star. It has been detected at radio frequencies and up to high-energy gamma rays (above 100 MeV). However, many models also predict a very high energy (VHE) emission (above hundreds of GeV) when the source displays relativistic persistent jets or transient ejections. Therefore, detecting such emission would improve the understanding of the jet physics. The imaging atmospheric Cherenkov telescope MAGIC observed Cygnus X-3 for about 70 hours between 2006 March and 2009 August in different X-ray/radio spectral states and also during a period of enhanced gamma-ray emission. MAGIC found no evidence for a VHE signal from the direction of the microquasar. An upper limit to the integral flux for energies higher than 250 GeV has been set to 2.2 x 10-12 photons cm-2 s-1 (95% confidence level). This is the best limit so far to the VHE emission from this source. The non-detection of a VHE signal during the period of activity in the high-energy band sheds light on the location of the possible VHE radiation favoring the emission from the innermost region of the jets, where absorption is significant. The current and future generations of Cherenkov telescopes may detect a signal under precise spectral conditions.

Gamma-ray binaries are suitable sources to study high-energy processes in jets and outflows in general. In the last years, there has been a lot of activity in the field of gamma-ray binaries to identify the different factors that shape their non-thermal spectra, which ranges from radio to very high energies, as well as their lightcurves. In this work, I discuss the main aspects of the non-thermal emission in this class of objects, which presently includes high-mass microquasars, high-mass binaries hosting a non-accreting pulsar and, probably, massive star binaries; few potential candidates to be gamma-ray binaries are also presented. Finally, the importance of gamma-ray absorption is discussed, and the main physical ingredients, which are likely involved in the non-thermal radiation in gamma-ray binaries, are briefly considered.

The Major Atmospheric Gamma Imaging Cerenkov (MAGIC) telescope participated in three multiwavelength (MWL) campaigns, observing the blazar Markarian (Mkn) 421 during the nights of 2006 April 28, 29, and 2006 June 14. We analyzed the corresponding MAGIC very-high energy observations during 9 nights from 2006 April 22 to 30 and on 2006 June 14. We inferred light curves with sub-day resolution and night-by-night energy spectra. A strong gamma-ray signal was detected from Mkn 421 on all observation nights. The flux (E > 250 GeV) varied on night-by-night basis between (0.92+-0.11)10^-10 cm^-2 s^-1 (0.57 Crab units) and (3.21+-0.15)10^-10 cm^-2 s^-1 (2.0 Crab units) in 2006 April. There is a clear indication for intra-night variability with a doubling time of 36+-10(stat) minutes on the night of 2006 April 29, establishing once more rapid flux variability for this object. For all individual nights gamma-ray spectra could be inferred, with power-law indices ranging from 1.66 to 2.47. We did not find statistically significant correlations between the spectral index and the flux state for individual nights. During the June 2006 campaign, a flux substantially lower than the one measured by the Whipple 10-m telescope four days later was found. Using a log-parabolic power law fit we deduced for some data sets the location of the spectral peak in the very-high energy regime. Our results confirm the indications of rising peak energy with increasing flux, as expected in leptonic acceleration models.

Context: Part of the very high energy $\gamma$-ray radiation coming from extragalactic sources is absorbed through the pair production process on the extragalactic background light photons. Extragalactic magnetic fields alter the trajectories of these cascade pairs and, in turn, convert cosmic background photons to gamma-ray energies by inverse Compton scattering. These secondary photons can form an extended halo around bright VHE sources. Aims: We searched for an extended emission around the bright blazars Mrk 421 and Mrk 501 using the MAGIC telescope data. Methods: If extended emission is present, the angular distribution of reconstructed gamma-ray arrival directions around the source is broader than for a point-like source. In the analysis of a few tens of hours of observational data taken from Mrk 421 and Mrk 501 we used a newly developed method that provides better angular resolution. This method is based on the usage of multidimensional decision trees. Comparing the measured shapes of angular distributions with those expected from a point-like source one can detect or constrain possible extended emission around the source. We also studied the influence of different types of systematic errors on the shape of the distribution of reconstructed gamma-ray arrival directions for a point source. Results: We present upper limits for an extended emission calculated for both sources for various source extensions and emission profiles. We discuss possible constraints on the extragalactic magnetic fields strength. We obtain upper limits on the extended emission around the Mrk~421 (Mrk~501) on the level of < 5% (< 4%) of the Crab Nebula flux above the energy threshold of 300 GeV. Using these results we discuss possible constraints on the extragalactic magnetic fields strength around a few times $10^{-15}$~G.

Dense and cold clouds seem to populate the broad line region surrounding the central black hole in AGNs. These clouds could interact with the AGN jet base and this could have observational consequences. We want to study the gamma-ray emission produced by these jet-cloud interactions, and explore under which conditions this radiation would be detectable. We investigate the hydrodynamical properties of jet-cloud interactions and the resulting shocks, and develop a model to compute the spectral energy distribution of the emission generated by the particles accelerated in these shocks. We discuss our model in the context of radio-loud AGNs, with applications to two representative cases, the low-luminous Centaurus A, and the powerful 3C 273. Some fraction of the jet power can be channelled to gamma-rays, which would be likely dominated by synchrotron self-Compton radiation, and show typical variability timescales similar to the cloud lifetime within the jet, which is longer than several hours. Many clouds can interact with the jet simultaneously leading to fluxes significantly higher than in one interaction, but then variability will be smoothed out. Jet-cloud interactions may produce detectable gamma-rays in non-blazar AGNs, of transient nature in nearby low-luminous sources like Cen A, and steady in the case of powerful objects of FR II type.

The blazar PG 1553+113 is a well known TeV gamma-ray emitter. In this paper, we determine its spectral energy distribution using simultaneous multi-frequency data in order to study its emission processes. An extensive campaign was carried out between March and April 2008, where optical, X-ray, high-energy (HE) gamma-ray, and very-high-energy (VHE) gamma-ray data were obtained with the KVA, Abastumani, REM, RossiXTE/ASM, AGILE and MAGIC telescopes, respectively. This is the first simultaneous broad-band (i.e., HE+VHE) gamma-ray observation, though AGILE did not detect the source. We combine data to derive source's spectral energy distribution and interpret its double peaked shape within the framework of a synchrotron self compton model

Context: Gamma-ray bursts are cosmological sources emitting radiation from the gamma-rays to the radio band. Substantial observational efforts have been devoted to the study of gamma-ray bursts during the prompt phase, i.e. the initial burst of high-energy radiation, and during the long-lasting afterglows. In spite of many successes in interpreting these phenomena, there are still several open key questions about the fundamental emission processes, their energetics and the environment. Aim: Independently of specific gamma-ray burst theoretical recipes, spectra in the GeV/TeV range are predicted to be remarkably simple, being satisfactorily modeled with power-laws, and therefore offer a very valuable tool to probe the extragalactic background light distribution. Furthermore, the simple detection of a component at very-high energies, i.e. at $\sim 100$\,GeV, would solve the ambiguity about the importance of various possible emission processes, which provide barely distinguishable scenarios at lower energies. Methods: We used the results of the MAGIC telescope observation of the moderate resdhift ($z\sim0.76$) \object{GRB\,080430} at energies above about 80\,GeV, to evaluate the perspective for late-afterglow observations with ground based GeV/TeV telescopes. Results: We obtained an upper limit of $F_{\rm 95%\,CL} = 5.5 \times 10^{-11}$\,erg\,cm$^{-2}$\,s$^{-1}$ for the very-high energy emission of \object{GRB\,080430}, which cannot set further constraints on the theoretical scenarios proposed for this object also due to the difficulties in modeling the low-energy afterglow. Nonetheless, our observations show that Cherenkov telescopes have already reached the required sensitivity to detect the GeV/TeV emission of GRBs at moderate redshift ($z \lesssim 0.8$), provided the observations are carried out at early times, close to the onset of their afterglow phase.

Active galactic nuclei present continuum and line emission. The emission lines are originated by gas located close to the central super-massive black hole. Some of these lines are broad, and would be produced in a small region called broad-line region. This region could be formed by clouds surrounding the central black hole. In this work, we study the interaction of such clouds with the base of the jets in active galactic nuclei, and we compute the produced high-energy emission. We focus on sources with low luminosities in the inner jet regions, to avoid strong gamma-ray absorption. We find that the resulting high-energy radiation may be significant in Centaurus A. Also, this phenomenon might be behind the variable gamma-ray emission detected in M87, if very large dark clouds are present. The detection of jet-cloud interactions in active galactic nuclei would give information on the properties of the jet base and the very central regions.

Several binary systems consisting of a massive star and a compact object have been detected above 100 GeV in the Galaxy. In most of these sources, gamma-rays show a modulation associated to the orbital motion, which means that the emitter should not be too far from the bright primary star. This implies that gamma-ray absorption will be non negligible, and large amounts of secondary electron-positron pairs will be created in the stellar surroundings. In this work, we show that the radio emission from these pairs should be accounted for when interpreting the radio spectrum, variability, and morphology found in gamma-ray binaries. Relevant features of the secondary radio emission are the relatively hard spectrum, the orbital motion of the radio peak center, and the extended radio structure following a spiral-like trajectory. The impact of the stellar wind free-free absorption should not be neglected.

High-mass microquasars may produce jets that will strongly interact with surrounding stellar winds on binary system spatial scales. We study the dynamics of the collision between a mildly relativistic hydrodynamical jet of supersonic nature and the wind of an OB star. We performed numerical 3D simulations of jets that cross the stellar wind with the code Ratpenat. The jet head generates a strong shock in the wind, and strong recollimation shocks occur due to the initial overpressure of the jet with its environment. These shocks can accelerate particles up to TeV energies and produce gamma-rays. The recollimation shock also strengthens jet asymmetric Kelvin-Helmholtz instabilities produced in the wind/jet contact discontinuity. This can lead to jet disruption even for jet powers of several times $10^{36}$ erg s$^{-1}$. High-mass microquasar jets likely suffer a strong recollimation shock that can be a site of particle acceleration up to very high energies, but also eventually lead to the disruption of the jet.

The Perseus galaxy cluster was observed by the MAGIC Cherenkov telescope for a total effective time of 24.4 hr during 2008 November and December. The resulting upper limits on the gamma-ray emission above 100 GeV are in the range of 4.6 to 7.5 x 10^{-12} cm^{-2} s^{-1} for spectral indices from -1.5 to -2.5, thereby constraining the emission produced by cosmic rays, dark matter annihilations, and the central radio galaxy NGC1275. Results are compatible with cosmological cluster simulations for the cosmic-ray-induced gamma-ray emission, constraining the average cosmic ray-to-thermal pressure to <4% for the cluster core region (<8% for the entire cluster). Using simplified assumptions adopted in earlier work (a power-law spectrum with an index of -2.1, constant cosmic ray-to-thermal pressure for the peripheral cluster regions while accounting for the adiabatic contraction during the cooling flow formation), we would limit the ratio of cosmic ray-to-thermal energy to E_CR/E_th<3%. The upper limit also translates into a level of gamma-ray emission from possible annihilations of the cluster dark matter (the dominant mass component) that is consistent with boost factors of ~10^4 for the typically expected dark matter annihilation-induced emission. Finally, the upper limits obtained for the gamma-ray emission of the central radio galaxy NGC1275 are consistent with the recent detection by the Fermi-LAT satellite. Due to the extremely large Doppler factors required for the jet, a one-zone synchrotron self-Compton model is implausible in this case. We reproduce the observed spectral energy density by using the structured jet (spine-layer) model which has previously been adopted to explain the high-energy emission of radio galaxies.

Massive protostars have associated bipolar outflows with velocities of hundreds of km s$^{-1}$. Such outflows can produce strong shocks when interact with the ambient medium leading to regions of non-thermal radio emission. We aim at exploring under which conditions relativistic particles are accelerated at the terminal shocks of the protostellar jets and can produce significant gamma-ray emission. We estimate the conditions necessary for particle acceleration up to very high energies and gamma-ray production in the non-thermal hot spots of jets associated with massive protostars embedded in dense molecular clouds. We show that relativistic Bremsstrahlung and proton-proton collisions can make molecular clouds with massive young stellar objects detectable by the {\it Fermi}{} satellite at MeV-GeV energies and by Cherenkov telescope arrays in the GeV-TeV range. Gamma-ray astronomy can be used to probe the physical conditions in star forming regions and particle acceleration processes in the complex environment of massive molecular clouds.

Powerful non-thermal emission has been detected coming from relativistic collimated outflows launched in the vicinity of black holes of a very wide range of masses, from few to $\sim 10^{10}$ M$_{\odot}$. At different scales along the outflows, i.e. $\sim 10-10^{10} R_{\rm Sch}$ from the black hole, the local conditions can lead to the generation of non-thermal populations of particles via, e.g., magnetic reconnection, magneto-centrifugal mechanisms, diffusive processes, or the so-called converter mechanism. These non-thermal populations of particles, interacting with dense matter, magnetic, and radiation fields, could yield radio-to-gamma-ray emission via synchrotron process, inverse Compton scattering, relativistic Bremsstrahlung, proton-proton and photo-hadron colissions, and even heavy nuclei photo-disintegration. Other processes, like pair creation or the development of electromagnetic cascades, could be also relevant in black-hole jets and their surroundings. Black holes of different masses, accretion rates and environments show different phenomenologies, as can be observed in AGNs, GRBs or microquasars. Nonetheless, these sources basically share the same fundamental physics: accretion, black-hole rotation, plus an environment, but they are individualized due to their own specific conditions. In this paper, we qualitatively review the main characteristics of the non-thermal emission produced in jets from black holes, giving also a brief overview on the physical properties of black hole/jet systems. We comment as well on some important differences and similarities between classes of sources, and on the prospects for the study of the non-thermal emission from astrophysical sources powered by black holes.

The discovery of very high energy (VHE) gamma-ray emitting X-ray binaries has triggered an intense effort to better understand the particle acceleration, absorption, and emission mechanisms in compact binary systems, which provide variable conditions along eccentric orbits. Despite this, the nature of some of these systems, and of the accelerated particles producing the VHE emission, is unclear. To answer some of these open questions, we conducted a multiwavelength campaign of the VHE gamma-ray emitting X-ray binary LS I +61 303 including the MAGIC telescope, XMM-Newton, and Swift during 60% of an orbit in 2007 September. We detect a simultaneous outburst at X-ray and VHE bands, with the peak at phase 0.62 and a similar shape at both wavelengths. A linear fit to the simultaneous X-ray/VHE pairs obtained during the outburst yields a correlation coefficient of r=0.97, while a linear fit to all simultaneous pairs provides r=0.81. Since a variable absorption of the VHE emission towards the observer is not expected for the data reported here, the correlation found indicates a simultaneity in the emission processes. Assuming that they are dominated by a single particle population, either hadronic or leptonic, the X-ray/VHE flux ratio favors leptonic models. This fact, together with the detected photon indices, suggests that in LS I +61 303 the X-rays are the result of synchrotron radiation of the same electrons that produce VHE emission as a result of inverse Compton scattering of stellar photons.

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

Massive protostars have associated bipolar outflows with velocities of hundreds of km/s. Such outflows produce strong shocks when interact with the ambient medium leading to regions of non-thermal radio emission. Under certain conditions, the population of relativistic particles accelerated at the terminal shocks of the protostellar jets can produce significant gamma-ray emission. We estimate the conditions necessary for high-energy emission in the non-thermal hot spots of jets associated with massive protostars embedded in dense molecular clouds. Our results show that particle-matter interactions can lead to the detection of molecular clouds hosting massive young stellar objects by the Fermi satellite at MeV-GeV energies and even by Cherenkov telescope arrays in the GeV-TeV range. Astronomy at gamma-rays can be used to probe the physical conditions in star forming regions and particle acceleration processes in the complex environment of massive molecular clouds.

High-mass microquasars are binary systems composed by a massive star and a compact object from which relativistic jets are launched. Regarding the companion star, observational evidence supports the idea that winds of hot stars are formed by clumps. Then, these inhomogeneities may interact with the jets producing a flaring activity. In the present contribution we study the interaction between a jet and a clump of the stellar wind in a high-mass microquasar. This interaction produces a shock in the jet, where particles may be accelerated up to relativistic energies. We calculate the spectral energy distributions of the dominant non-thermal processes: synchrotron radiation, inverse Compton scattering, and proton-proton collisions. Significant levels of X- and gamma-ray emission are predicted, with luminosities in the different domains up to ~ 10^{34} - 10^{35} erg/s on a timescale of about ~ 1 h. Finally, jet-clump interactions in high-mass microquasars could be detectable at high energies. These phenomena may be behind the fast TeV variability found in some high-mass X-ray binary systems, such as Cygnus X-1, LS 5039 and LS I+61 303. In addition, our model can help to derive information on the properties of jets and clumpy winds.

Active galactic nuclei present continuum and line emission. The former is produced by the accretion disk and the jets, whereas the latter is originated by gas located close to the super-massive black hole. The small region where the broad lines are emitted is called the broad-line region. The structure of this region is not well known, although it has been proposed that it may be formed by small and dense ionized clouds surrounding the supermassive black-hole. In this work, we study the interaction of one cloud from the broad line region with the jet of the active galactic nuclei. We explore the high-energy emission produced by this interaction close to the base of the jet. The resulting radiation may be detectable for nearby non-blazar sources as well as for powerful quasars, and its detection could give important information on the broad line region and the jet itself.

Microquasars are X-ray binaries with relativistic jets. These jets are powerful energy carriers, thought to be fed by accretion, which produce non-thermal emission at different energy bands. To date, several Galactic sources showing extended radio emission, among them at least one confirmed microquasar, Cygnus X-1, have been detected in the TeV range. All of them show complex patterns of spectral and temporal behavior. In this work, we discuss the physics behind the very high-energy emission in microquasars. In concrete, we focus on the microquasar Cygnus X-1, and also in the other two TeV binaries with detected extended outflows, LS 5039 and LS I +61 303, pointing out relevant aspects of the complex phenomena occurring in them. We conclude that the TeV emission is likely of leptonic origin, although hadrons cannot be discarded. In addition, efficient electromagnetic cascades can hardly develop since even relatively low magnetic fields suppress them. Also, the modeling of the radiation from some of the detected sources points to them as either extremely efficient accelerators, and/or having the TeV emitter at a distance from the compact object of about 10e12 cm. Finally, we point out that the role of a massive and hot stellar companion, due to its strong photon field and wind, cannot be neglected when trying to understand the behavior of microquasars at high and very high energies. The complexity of microquasars precludes straightforward generalizations to a whole population, and are better studied presently in a source by source base.

Index entry listing the contributed papers of the MAGIC collaboration to the 31th International Cosmic Ray Conference (ICRC 2009), July 7-15 2009, \L\'od\'z, Poland. The individual papers are sorted by subject: Overview and Highlight Papers; MAGIC-II Status and Components; Software and Analysis Techniques; Technical Developments; Scientific Results. This HTML document includes clickable links to the papers that exist on the astro-ph arXiv. We hope that this will make it easy to access the MAGIC contributions in a systematic way.

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

The discovery of emission of TeV gamma rays from X-ray binaries has triggered an intense effort to better understand the particle acceleration, absorption, and emission mechanisms in compact binary systems. Here we present the pioneering effort of the MAGIC collaboration to understand the very high energy emission of the prototype system LS I +61 303. We report on the variable nature of the emission from LS I +61 303 and show that this emission is indeed periodic. The system shows regular outburst at TeV energies in phase phi=0.6-0.7 and detect no signal at periastron (phi~ 0.275). Furthermore we find no indication of spectral variation along the orbit of the compact object and the spectral energy distribution is compatible with a simple power law with index Gamma=2.6+-0.2_(stat)+-0.2_(sys). To answer some of the open questions concerning the emission process of the TeV radiation we conducted a multiwavelength campaign with the MAGIC telescope, XMM-Newton, and Swift in September 2007. We detect a simultaneous outburst at X-ray and TeV energies, with the peak at phase 0.62 and a similar shape at both wavelengths. A linear fit to the strictly simultaneous X-ray/TeV flux pairs provides r=0.81 -0.21 +0.06. Here we present the observations and discuss the implications of the obtained results to the emission processes in the system.

High-mass microquasars are binary systems consisting of a massive star and an accreting compact object from which relativistic jets are launched. There is considerable observational evidence that winds of massive stars are clumpy. Individual clumps may interact with the jets in high-mass microquasars to produce outbursts of high-energy emission. Gamma-ray flares have been detected in some high-mass X-ray binaries, such as Cygnus X-1, and probably in LS 5039 and LS I+61 303. We predict the high-energy emission produced by the interaction between a jet and a clump of the stellar wind in a high-mass microquasar. Assuming a hydrodynamic scenario for the jet-clump interaction, we calculate the spectral energy distributions produced by the dominant non-thermal processes: relativistic bremsstrahlung, synchrotron and inverse Compton radiation, for leptons, and for hadrons, proton-proton collisions. Significant levels of emission in X-rays (synchrotron), high-energy gamma rays (inverse Compton), and very high-energy gamma rays (from the decay of neutral pions) are predicted, with luminosities in the different domains in the range ~ 10^{32}-10^{35} erg/s. The spectral energy distributions vary strongly depending on the specific conditions. Jet-clump interactions may be detectable at high and very high energies, and provide an explanation for the fast TeV variability found in some high-mass X-ray binary systems. Our model can help to infer information about the properties of jets and clumpy winds by means of high-sensitivity gamma-ray astronomy.

We report on the results from Suzaku broadband X-ray observations of the galactic binary source LS5039. The Suzaku data, which have continuous coverage of more than one orbital period, show strong modulation of the X-ray emission at the orbital period of this TeV gamma-ray emitting system.The X-ray emission shows a minimum at orbital phase ~ 0.1, close to the so-called superior conjunction of the compact object, and a maximum at phase ~0.7, very close to the inferior conjunction of the compact object. The X-ray spectral data up to 70 keV are described by a hard power-law with a phase-dependent photon index which varies within Gamma ~1.45 - 1.61. The amplitude of the flux variation is a factor of 2.5, but is significantly less than that of the factor ~8 variation in the TeV flux. Otherwise the two light curves are similar, but not identical. Although periodic X-ray emission has been found from many galactic binary systems, the Suzaku result implies a phenomenon different from the "standard" origin of X-rays related to the emission of the hot accretion plasma formed around the compact companion object. The X-ray radiation of LS5039is likely to be linked to very-high-energy electrons which are also responsible for the TeV gamma-ray emission. While the gamma-rays are the result of inverse Compton scattering by electrons on optical stellar photons, X-rays are produced via synchrotron radiation. Yet, while the modulation of the TeV gamma-ray signal can be naturally explained by the photon-photon pair production and anisotropic inverse Compton scattering, the observed modulation of synchrotron X-rays requires an additional process, the most natural one being adiabatic expansion in the radiation production region.

The high kinetic energy outflowing in the jets of microquasars is delivered to the surrounding interstellar medium. This energy input can cause the formation of bow shocks and cocoons that may be detectable from radio to gamma-ray energies. Evidences or hints of emission from jet/medium interactions have been reported for some sources, but little has been done regarding the theoretical modeling of the resulting non-thermal emission. We have developed an analytical model based on those successfully applied to extragalactic sources or the interaction of AGN jets with their surroundings. Focusing on the adiabatic phase of the growing structures, we give estimations of the expected luminosities through synchrotron, relativistic Bremsstrahlung and inverse Compton processes. We conclude that the interaction structures may be detectable at radio wavelengths, while extreme values for the jet kinetic power, the source age and the medium density are required to make the emission at high and very high energies detectable.

The binary LS 5039 is a non-thermal X-ray emitter that presents jet-like radio structures, and is also one of the most misterious TeV sources in our Galaxy. The presence of an O-type star in LS 5039 implies that the non-thermal emitter must be embedded in a strong stellar wind, and the role of the latter could be relevant for the understanding of the high-energy behavior of the source. In this work, we show that the lack of absorption features in the soft X-ray spectrum of LS 5039 can constrain strongly the parameters that describe the wind, and ultimately the location of the non-thermal emitter.

The interaction of microquasar jets with their environment can produce non-thermal radiation as is the case for extragalactic outflows impacting on their surroundings. We have developed an analytical model based on those successfully applied to extragalactic sources. The jet is taken to be a supersonic and mildly relativistic hydrodynamical outflow. We focus on the jet/shocked medium structure when being in its adiabatic phase, and assume that it grows in a self-similar way. We calculate the fluxes and spectra of the radiation produced via synchrotron, Inverse Compton and relativistic Bremsstrahlung processes by electrons accelerated in strong shocks. A hydrodynamical simulation is also performed to further investigate the jet interaction with the environment and check the physical parameters used in the analytical model. We conclude that microquasar jet termination regions could be detectable at radio wavelengths for current instruments sensitive to arcminute scales while at X-rays the expected luminosities are moderate, although the emitter is more compact than the radio one. The radiation at gamma-ray energies may be within the detection limits of the next generation of satellite and ground-based instruments.

LS 5039 is among the most interesting VHE sources in the Galaxy. Two scenarios have been put forward to explain the observed TeV radiation: jets vs pulsar winds. The source has been detected during the superior conjunction of the compact object, when very large gamma-ray opacities are expected. In addition, electromagnetic cascades, which may make the system more transparent to gamma-rays, are hardly efficient for realistic magnetic fields in massive star surroundings. All this makes unlikely the standard pulsar scenario for LS 5039, in which the emitter is the region located between the star and the compact object, where the opacities are the largest. Otherwise, a jet-like flow can transport energy to regions where the photon-photon absorption is much lower and the TeV radiation is not so severely absorbed.

LS 5039 is an X-ray binary that presents non-thermal radio emission. The radiation at $\sim 5$ GHz is quite steady and optically thin, consisting on a dominant core plus an extended jet-like structure. There is a spectral turnover around 1 GHz, and evidence of variability at timescales of 1 yr at 234 MHz. We investigate the radio emitter properties using the available broadband radio data, and assuming two possible scenarios to explain the turnover: free-free absorption in the stellar wind, or synchrotron self-absorption. We use the relationships between the turnover frequency, the stellar wind density, the emitter location, size and magnetic field, and the Lorentz factor of the emitting electrons, as well as a reasonable assumption on the energy budget, to infer the properties of the low-frequency radio emitter. Also, we put this information in context with the broadband radio data. The location and size of the low-frequency radio emitter can be restricted to $\ga$ few AU from the primary star, its magnetic field to $\sim 3\times 10^{-3}-1$ G, and the electron Lorentz factors to $\sim 10-100$. The observed variability of the extended structures seen with VLBA would point to electron bulk velocities $\ga 3\times 10^8$ cm s$^{-1}$, whereas much less variable radiation at 5 GHz would indicate velocities for the VLBA core $\la 10^8$ cm s$^{-1}$. The emission at 234 MHz in the high state would mostly come from a region larger than the dominant broadband radio emitter. We suggest a scenario in which secondary pairs, created via gamma-ray absorption and moving in the stellar wind, are behind the steady broadband radio core, whereas the resolved jet-like radio emission would come from a collimated, faster, outflow.

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.

Microquasars are binary systems that harbor a normal star and a compact object (black-hole or neutron star), and show relativistic outflows (or jets). The matter that forms these jets is of likely stellar origin, previously expelled from the star and trapped in the potential well of the compact object. This matter is accreted by the compact object, forming a disk due to its angular momentum, and is eventually ejected in the form of a bipolar outflow (the jets), which generates radio emission and could also be a very high-energy emitter. To study and understand the radiation from microquasars, there is a set of elements that can play a major role and are to be taken into account: the photons and the expelled matter from the star in the case of high-mass systems; the accreted matter radiation; the jet; the magnetic field carried by the jet or filling the binary system; and the medium surrounding the microquasar at large scales (~pc). In this lecture, we consider these elements of the microquasar scenario and briefly describe the physical conditions and processes involved in the production of non-thermal radiation from radio to gamma-rays. The required energetics, particle acceleration and transport, several radiative mechanisms, and the impact of different photon absorption processes, are discussed.

LS 5039 is an X-ray binary detected at very high energies. Along the orbit, there is a significant detection even during the superior conjunction of the compact object, when very large gamma-ray opacities are expected. Electromagnetic cascades, which may make the system more transparent to gamma-rays, are hardly efficient for reasonable magnetic fields in the massive star surroundings. A jet-like flow could transport energy to regions where the photon-photon absorption is much lower and the TeV radiation is not so severely absorbed. Otherwise, in the standard pulsar scenario for LS 5039, the emitter would be located between the star and the compact object, which would imply the violation of the observational constraints at X-rays.

Cygnus X-1 and LS 5039 are two X-ray binaries observed at TeV energies. Both sources are compact systems, contain jet-like (radio) structures, and harbor very luminous O stars. A TeV signal has been found around the superior conjunction of the compact object in both objects, when the highest gamma-ray opacities are expected. We investigate the implications of finding TeV emission from Cygnus X-1 and LS 5039 around the superior conjunction, since this can give information on the system magnetic field and the location of the TeV emitter. Using the very high-energy spectra and fluxes observed around the superior conjunction in Cygnus X-1 and LS 5039, we compute the absorbed luminosity that is caused by pair creation in the stellar photon field for different emitter positions with respect to the star and the observer line of sight. The role of the magnetic field and electromagnetic cascading are discussed. For the case of inefficient electromagnetic cascading, the expected secondary synchrotron fluxes are compared with the observed ones at X-ray energies. We find that, in Cygnus X-1 and LS 5039, either the magnetic field in the star surroundings is much smaller than the one expected for O stars or the TeV emitter is located at a distance >10^12 cm from the compact object. Our results strongly suggest that the TeV emitters in Cygnus X-1 and LS 5039 are located at the borders of the binary system and well above the orbital plane. This would not agree with those models for which the emitter is well inside the system, like the innermost-jet region (Cygnus X-1 and LS 5039; microquasar scenario), or the region between the pulsar and the primary star (LS 5039; standard pulsar scenario).

Thermal radio and X-ray emission has been traditionally associated with the formation of stars. However, in recent years, non-thermal radiation from massive star forming regions has been detected. Synchrotron radio emission and non-thermal X-rays from the outflows of massive young stellar objects (YSOs) provide evidence of the presence of relativistic particles in these sources. In YSOs, the acceleration of particles is likely produced where the thermal jet impacts on the surrounding medium and a shock wave is formed. Thus, particles might be accelerated up to relativistic energies through a Fermi-I type mechanism. Relativistic electrons and protons can interact with thermal particles and photons, producing then gamma-rays. These energetic photons could be detected by the new generation of instruments, making massive YSOs a new population of gamma-ray surces. In the present contribution we briefly describe some massive star forming regions from which non-thermal radio emission has been detected. In addition, we present a general model for high-energy radiation from the massive YSOs embedded in these regions. We take into account both leptonic and hadronic interactions of particles accelerated at the termination points of the collimated outflows ejected from the protostar.

Recent radio observations support a picture for star formation where there is accretion of matter onto a central protostar with the ejection of molecular outflows that can affect the surrounding medium. The impact of a supersonic outflow on the ambient gas can produce a strong shock that could accelerate particles up to relativistic energies. A strong evidence of this has been the detection of non-thermal radio emission coming from the jet termination region of some young massive stars. In the present contribution, we study the possible high-energy emission due to the interaction of relativistic particles, electrons and protons, with the magnetic, photon and matter fields inside a giant molecular cloud. Electrons lose energy via relativistic Bremsstrahlung, synchrotron radiation and inverse Compton interactions, and protons cool mainly through inelastic collisions with atoms in the cloud. We conclude that some massive young stellar objects might be detectable at gamma-rays by next generation instruments, both satellite-borne and ground based.

There are three Galactic jet sources, from which TeV emission has been detected: LS 5039, LS I +61 303 and Cygnus X-1. These three sources show power-law tails at X-rays and soft gamma-rays that could indicate a non-thermal origin for this radiation. In addition, all three sources apparently show correlated and complex behavior at X-ray and TeV energies. In some cases, this complex behavior is related to the orbital motion (e.g. LS 5039, LS I +61 303), and in some others it is related to some transient event occurring in the system (e.g. Cygnus X-1, and likely also LS I +61 303 and LS 5039). Based on modeling or energetic grounds, it seems difficult to explain the emission in the X-/soft gamma-ray and the TeV bands as coming from the same region (i.e. one-zone). We also point out the importance of the pair creation phenomena in these systems, which harbor a massive and hot star, for the radio and the X-ray emission, since a secondary pair radiation component may be significant in these energy ranges. Finally, we discuss that in fact the presence of the star can indeed have strong impact on, beside the non-thermal radiation production, the jet dynamics.

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.

Strong interactions between jets and stellar winds at binary system spatial scales could occur in high-mass microquasars. We study here, mainly from a dynamical but also a radiative point of view, the collision between a dense stellar wind and a mildly relativistic hydrodynamical jet of supersonic nature.}{We have performed numerical 2-dimensional simulations of jets, with cylindrical and planar (slab) symmetry, crossing the stellar wind material. From the results of the simulations, we derive estimates of the particle acceleration efficiency, using first order Fermi acceleration theory, and give some insight on the possible radiative outcomes. We find that, during jet launching, the jet head generates a strong shock in the wind. During and after this process, strong recollimation shocks can occur due to the initial overpressure of the jet with its environment. The conditions in all these shocks are convenient to accelerate particles up to $\sim$ TeV energies, which can lead to leptonic (synchrotron and inverse Compton) and hadronic (proton-proton) radiation. In principle, the cylindrical jet simulations show that the jet is stable, and can escape from the system even for relatively low power. However, when accounting for the wind ram pressure, the jet can be bent and disrupted for power $\la 10^{36}$ erg s$^{-1}$.

Massive hot stars produce dense ultraviolet (UV) photon fields in their surroundings. If a very high-energy (VHE) gamma-ray emitter is located close to the star, then gamma-rays are absorbed in the stellar photon field, creating secondary (electron-positron) pairs. We study the broadband emission of these secondary pairs in the stellar photon and magnetic fields. Under certain assumptions on the stellar wind and the magnetic field in the surroundings of a massive hot star, we calculate the steady state energy distribution of secondary pairs created in the system and its radiation from radio to gamma-rays. Under the ambient magnetic field, possibly high enough to suppress electromagnetic (EM) cascading, the energy of secondary pairs is radiated via synchrotron and single IC scattering producing radio-to-gamma-ray radiation. The synchrotron spectral energy distribution (SED) is hard, peaks around X-ray energies, and becomes softer. The IC SED is hard as well and peaks around 10 GeV, becoming also softer at higher energies due to synchrotron loss dominance. The radio emission from secondary pairs is moderate and detectable as a point-like and/or extended source. In X-rays, the secondary pair synchrotron component may be dominant. At energies <10 GeV, the secondary pair IC radiation may be the dominant primary gamma-ray emission and possibly detectable by the next generation of instruments.

The MAGIC Project: Contributions to ICRC 2007, Merida, Mexico. Contents pages for the Contribution on behalf of the MAGIC Collaboration to the 30th ICRC that took place in July 2007 in Merida, Mexico. The contents are in html form with clickable links to the papers that exist on the Astrophysics archive. We hope that this will make it easier to access the output of the conference in a systematic way. Comments on how useful this is/ how it could be improved should be sent to michela.demaria@iuav.it.

In high-mass microquasars (HMMQ), strong interactions between jets and stellar winds at binary system scales could occur. In order to explore this possibility, we have performed numerical 2-dimensional simulations of jets crossing the dense stellar material to study how the jet will be affected by these interactions. We find that the jet head generates strong shocks in the wind. These shocks reduce the jet advance speed, and compress and heat up jet and wind material. In addition, strong recollimation shocks can occur where pressure balance between the jet side and the surrounding medium is reached. All this, altogether with jet bending, could lead to the destruction of jets with power $<10^{36} \rm{erg/s}$. The conditions around the outflow shocks would be convenient for accelerating particles up to $\sim $TeV energies. These accelerated particles could emit via synchrotron and inverse Compton (IC) scattering if they were leptons, and via hadronic processes in case they were hadrons.

The recent detections of TeV gamma-rays from compact binary systems show that relativistic outflows (jets or winds) are sites of effective acceleration of particles up to multi-TeV energies. In this paper, we discuss the conditions of acceleration and radiation of ultra-relativistic electrons in LS 5039, the gamma-ray emitting binary system for which the highest quality TeV data are available. Assuming that the gamma-ray emitter is a jet-like structure, we performed detailed numerical calculations of the energy spectrum and lightcurves accounting for the acceleration efficiency, the location of the accelerator, the speed of the emitting flow, the inclination angle of the system, as well as specific features related to anisotropic inverse Compton scattering and pair production. We conclude that the accelerator should not be deep inside the binary system unless we assume a very efficient acceleration rate. We show that within the IC scenario both the gamma-ray spectrum and flux are strongly orbital phase dependent. Formally, our model can reproduce, for specific sets of parameter values, the energy spectrum of gamma-rays reported by HESS for wide orbital phase intervals. However, the physical properties of the source can be constrained only by observations capable of providing detailed energy spectra for narrow orbital phase intervals ($\Delta\phi\ll 0.1$).

Context: It has been proposed that the origin of the very high-energy photons emitted from high-mass X-ray binaries with jet-like features, so-called microquasars (MQs), is related to hadronic interactions between relativistic protons in the jet and cold protons of the stellar wind. Leptonic secondary emission should be calculated in a complete hadronic model that include the effects of pairs from charged pion decays inside the jets and the emission from pairs generated by gamma-ray absorption in the photosphere of the system. Aims: We aim at predicting the broadband spectrum from a general hadronic microquasar model, taking into account the emission from secondaries created by charged pion decay inside the jet. Methods: The particle energy distribution for secondary leptons injected along the jets is consistently derived taking the energy losses into account. We also compute the spectral energy distribution resulting from these leptons is calculated after assuming different values of the magnetic field inside the jets. The spectrum of the gamma-rays produced by neutral pion-decay and processed by electromagnetic cascades under the stellar photon field. Results: We show that the secondary emission can dominate the spectral energy distribution at low energies (~1 MeV). At high energies, the production spectrum can be significantly distorted by the effect of electromagnetic cascades. These effects are phase-dependent, and some variability modulated by the orbital period is predicted.

We present a model to estimate the synchrotron radio emission generated in microquasar (MQ) jets due to secondary pairs created via decay of charged pions produced in proton-proton collisions between stellar wind ions and jet relativistic protons. Signatures of electrons/positrons are obtained from consistent particle energy distributions that take into account energy losses due to synchrotron and inverse Compton (IC) processes, as well as adiabatic expansion. The space parameter for the model is explored and the corresponding spectral energy distributions (SEDs) are presented. We conclude that secondary leptonic emission represents a significant though hardly dominant contribution to the total radio emission in MQs, with observational consequences that can be used to test some still unknown processes occurring in these objects as well as the nature of the matter outflowing in their jets.

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

Massive stars form in dense and massive molecular cores. The exact formation mechanism is unclear, but it is possible that some massive stars are formed by processes similar to those that produce the low-mass stars, with accretion/ejection phenomena occurring at some point of the evolution of the protostar. This picture seems to be supported by the detection of a collimated stellar wind emanating from the massive protostar IRAS 16547-4247. A triple radio source is associated with the protostar: a compact core and two radio lobes. The emission of the southern lobe is clearly non-thermal. Such emission is interpreted as synchrotron radiation produced by relativistic electrons locally accelerated at the termination point of a thermal jet. Since the ambient medium is determined by the properties of the molecular cloud in which the whole system is embedded, we can expect high densities of particles and infrared photons. Because of the confirmed presence of relativistic electrons, inverse Compton and relativistic Bremsstrahlung interactions are unavoidable. Proton-proton collision should also occur, producing an injection of neutral pions. In this paper we aim at making quantitative predictions of the spectral energy distribution of the non-thermal spots generated by massive young stellar objects, with emphasis on the particular case of IRAS 16547-4247. We present spectral energy distributions for the southern lobe of this source, for a variety of conditions. We show that high-energy emission might be detectable from this object in the gamma-ray domain (MeV to TeV). The source may also be detectable at X-rays through long exposures with current X-ray instruments.

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

The fast repositioning system of the MAGIC Telescope has allowed during its first data cycle, between 2005 and the beginning of year 2006, observing nine different GRBs as possible sources of very high energy gammas. These observations were triggered by alerts from Swift, HETE-II, and Integral; they started as fast as possible after the alerts and lasted for several minutes, with an energy threshold varying between 80 and 200 GeV, depending upon the zenith angle of the burst. No evidence for gamma signals was found, and upper limits for the flux were derived for all events, using the standard analysis chain of MAGIC. For the bursts with measured redshift, the upper limits are compatible with a power law extrapolation, when the intrinsic fluxes are evaluated taking into account the attenuation due to the scattering in the Metagalactic Radiation Field (MRF).

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.

Microquasar (MQ) jets are sites of particle acceleration and synchrotron emission. Such synchrotron radiation has been detected coming from jet regions of different spatial scales, which for the instruments at work nowadays appear as compact radio cores, slightly resolved radio jets, or (very) extended structures. Because of the presence of relativistic particles and dense photon, magnetic and matter fields, these outflows are also the best candidates to generate the very high-energy (VHE) gamma-rays detected coming from two of these objects, LS 5039 and LS I +61 303, and may be contributing significantly to the X-rays emitted from the MQ core. In addition, beside electromagnetic radiation, jets at different scales are producing some amount of leptonic and hadronic cosmic rays (CR), and evidences of neutrino production in these objects may be eventually found. In this work, we review on the different physical processes that may be at work in or related to MQ jets. The jet regions capable to produce significant amounts of emission at different wavelengths have been reduced to the jet base, the jet at scales of the order of the size of the system orbital semi-major axis, the jet middle scales (the resolved radio jets), and the jet termination point. The surroundings of the jet could be sites of multiwavelegnth emission as well, deserving also an insight. We focus on those scenarios, either hadronic or leptonic, in which it seems more plausible to generate both photons from radio to VHE and high-energy neutrinos. We briefly comment as well on the relevance of MQ as possible contributors to the galactic CR in the GeV-PeV range.

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.

The COMPTEL unidentified source GRO J 1411-64 was observed by INTEGRAL and XMM-Newton in 2005. The Circinus Galaxy is the only source detected within the 4$\sigma$ location error of GRO J1411-64, but in here excluded as the possible counterpart. At soft X-rays, 22 reliable and statistically significant sources (likelihood $> 10$) were extracted and analyzed from XMM-Newton data. Only one of these sources, XMMU J141255.6 -635932, is spectrally compatible with GRO J1411-64 although the fact the soft X-ray observations do not cover the full extent of the COMPTEL source position uncertainty make an association hard to quantify and thus risky. At the best location of the source, detections at hard X-rays show only upper limits, which, together with MeV results obtained by COMPTEL suggest the existence of a peak in power output located somewhere between 300-700 keV for the so-called low state. Such a spectrum resembles those in blazars or microquasars, and might suggest at work by a similar scenario. However, an analysis using a microquasar model consisting on a magnetized conical jet filled with relativistic electrons, shows that it is hard to comply with all observational constrains. This fact and the non-detection at hard X-rays introduce an a-posteriori question mark upon the physical reality of this source, what is discussed here.

We consider the acceleration of energetic particles by Fermi processes (i.e., diffusive shock acceleration, second order Fermi acceleration, and gradual shear acceleration) in relativistic astrophysical jets, with particular attention given to recent progress in the field of viscous shear acceleration. We analyze the associated acceleration timescales and the resulting particle distributions, and discuss the relevance of these processes for the acceleration of charged particles in the jets of AGNs, GRBs and microquasars, showing that multi-component powerlaw-type particle distributions are likely to occur.

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

Context. The microquasar 1E 1740.7-2942 is a source located in the direction of the Galactic Center. It has been detected at X-rays, soft gamma-rays, and in the radio band, showing an extended radio component in the form of a double-sided jet. Although no optical counterpart has been found so far for 1E 1740.7-2942, its X-ray activity strongly points to a galactic nature. Aims. We aim to improve our understanding of the hard X-ray and gamma-ray production in the system, exploring whether the jet can emit significantly at high energies under the light of the present knowledge. Methods. We have modeled the source emission, from radio to gamma-rays, with a cold-matter dominated jet model. INTEGRAL data combined with radio and RXTE data, as well as EGRET and HESS upper-limits, are used to compare the computed and the observed spectra. Results. From our modeling, we find out that jet emission cannot explain the high fluxes observed at hard X-rays without violating at the same time the constraints from the radio data, favoring the corona origin of the hard X-rays. Also, 1E 1740.7-2942 might be detected by GLAST or AGILE at GeV energies, and by HESS and HESS-II beyond 100 GeV, with the spectral shape likely affected by photon-photon absorption in the disk and corona photon fields.

The COMPTEL unidentified source GRO J1411-64 was observed by INTEGRAL, and its central part, also by XMM-Newton. The data analysis shows no hint for new detections at hard X-rays. The upper limits in flux herein presented constrain the energy spectrum of whatever was producing GRO J1411-64, imposing, in the framework of earlier COMPTEL observations, the existence of a peak in power output located somewhere between 300-700 keV for the so-called low state. The Circinus Galaxy is the only source detected within the 4$\sigma$ location error of GRO J1411-64, but can be safely excluded as the possible counterpart: the extrapolation of the energy spectrum is well below the one for GRO J1411-64 at MeV energies. 22 significant sources (likelihood $> 10$) were extracted and analyzed from XMM-Newton data. Only one of these sources, XMMU J141255.6-635932, is spectrally compatible with GRO J1411-64 although the fact the soft X-ray observations do not cover the full extent of the COMPTEL source position uncertainty make an association hard to quantify and thus risky. The unique peak of the power output at high energies (hard X-rays and gamma-rays) resembles that found in the SED seen in blazars or microquasars. However, an analysis using a microquasar model consisting on a magnetized conical jet filled with relativistic electrons which radiate through synchrotron and inverse Compton scattering with star, disk, corona and synchrotron photons shows that it is hard to comply with all observational constrains. This and the non-detection at hard X-rays introduce an a-posteriori question mark upon the physical reality of this source, which is discussed in some detail.

Microquasars (MQs) 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 extended emission detectable at different energy bands through several processes: neutral pion-decay produce high-energy and very high-energy gamma-rays, secondary electrons produced by charged pion-decay generate synchrotron and bremsstrahlung emission. In addition, the jets of MQs themselves are likely sources of gamma-rays. We discuss about the association between the intrinsic and the indirect emission coming from these objects.

Microquasars are X-ray binaries with relativistic jets. The microquasar LS 5039 turned out to be the first high-energy gamma-ray microquasar candidate due to its likely association with the EGRET source 3EG J1824-1514. Further theoretical studies supported this association, which could be extended to other EGRET sources. Very recently, Aharonian et al. (2005) have communicated the detection of the microquasar LS 5039 at TeV energies. This fact confirms the EGRET source association and leaves no doubt about the gamma-ray emitting nature of this object. The aim of the present work is to show that, applying a cold-matter dominated jet model to LS 5039, we can reproduce many of the spectral and variability features observed in this source. Jet physics is explored, and some physical quantities are estimated as a by-product of the performed modeling. Although at the moment only LS 5039 has been detected on the entire electromagnetic spectrum, it does not seem unlikely that other microquasars will show similar spectral properties. Therefore, an in-depth study of the first gamma-ray microquasar, on theoretical grounds supported by observations, can render a useful knowledge applicable elsewhere.

Microquasars are sources of very high-energy gamma-rays and, very probably, high-energy gamma-ray emitters. We propose a model for a jet that can allow to give accurate observational predictions for jet emission at different energies and provide with physical information of the object using multiwavelength data.

Microquasars (MQs) present emission over the whole spectrum, from radio wavelengths to gamma-rays. The microquasar spectral energy distribution is very complex, being a signature of the different physical processes that generate the radiation emitted by these objects. In this work, we estimate the amount of broad-band emission produced by relativistic protons, released from the jet of a MQ, interacting with high density regions of the ISM. We show that a two components source, the microquasar itself and the region of interaction between the jets and the ISM, could be unveiled by the new instruments at high-energy and very high-energy gamma-rays.

Among unidentified gamma-ray sources in the galactic plane, there are some that present significant variability and have been proposed to be high-mass microquasars. To deepen the study of the possible association between variable low galactic latitude gamma-ray sources and microquasars, we have applied a leptonic jet model based on the microquasar scenario that reproduces the gamma-ray spectrum of three unidentified gamma-ray sources, 3EG J1735-1500, 3EG J1828+0142 and GRO J1411-64, and is consistent with the observational constraints at lower energies. We conclude that if these sources were generated by microquasars, the particle acceleration processes could not be as efficient as in other objects of this type that present harder gamma-ray spectra. Moreover, the dominant mechanism of high-energy emission should be synchrotron self-Compton (SSC) scattering, and the radio jets may only be observed at low frequencies. For each particular case, further predictions of jet physical conditions and variability generation mechanisms have been made in the context of the model. Although there might be other candidates able to explain the emission coming from these sources, microquasars cannot be excluded as counterparts. Observations performed by the next generation of gamma-ray instruments, like GLAST, are required to test the proposed model.

The microquasar LS 5039 has recently been detected as a source of very high energy (VHE) $\gamma$-rays. This detection, that confirms the previously proposed association of LS 5039 with the EGRET source 3EG~J1824$-$1514, makes of LS 5039 a special system with observational data covering nearly all the electromagnetic spectrum. In order to reproduce the observed spectrum of LS 5039, from radio to VHE $\gamma$-rays, we have applied a cold matter dominated jet model that takes into account accretion variability, the jet magnetic field, particle acceleration, adiabatic and radiative losses, microscopic energy conservation in the jet, and pair creation and absorption due to the external photon fields, as well as the emission from the first generation of secondaries. The radiative processes taken into account are synchrotron, relativistic Bremsstrahlung and inverse Compton (IC). The model is based on a scenario that has been characterized with recent observational results, concerning the orbital parameters, the orbital variability at X-rays and the nature of the compact object. The computed spectral energy distribution (SED) shows a good agreement with the available observational data.

Observational and theoretical studies point to microquasars (MQs) as possible counterparts of a significant fraction of the unidentified gamma-ray sources detected so far. At present, a proper scenario to explain the emission beyond soft X-rays from these objects is not known, nor what the precise connection is between the radio and the high-energy radiation. We develop a new model where the MQ jet is dynamically dominated by cold protons and radiatively dominated by relativistic leptons. The matter content and power of the jet are both related with the accretion process. The magnetic field is assumed to be close to equipartition, although it is attached to and dominated by the jet matter. For the relativistic particles in the jet, their maximum energy depends on both the acceleration efficiency and the energy losses. The model takes into account the interaction of the relativistic jet particles with the magnetic field and all the photon and matter fields. Such interaction produces significant amounts of radiation from radio to very high energies through synchrotron, relativistic Bremsstrahlung, and inverse Compton (IC) processes. Variability of the emission produced by changes in the accretion process (e.g. via orbital eccentricity) is also expected. The effects of the gamma-ray absorption by the external photon fields on the gamma-ray spectrum have been taken into account, revealing clear spectral features that might be observed. This model is consistent to the accretion scenario, energy conservation laws, and current observational knowledge, and can provide deeper physical information of the source when tested against multiwavelength data.

Population studies of EGRET gamma-ray sources indicate that there is a distinctive population of bright sources at low galactic latitudes. The sources have a distribution consistent with that of young galactic objects, with a concentration toward the inner spiral arms. There is a subgroup that displays strong variability with timescales from days to months. Following an earlier suggestion by Kaufman Bernado et al. (2002), we explore the possibility that these sources could be high-mass microquasars. Detailed models for the gamma-ray emission that include inverse Compton interactions of electrons in the relativistic jets and photons from all local fields (stellar UV photons, synchrotron photons, soft X-ray photons from the accretion disk, and hard X-ray photons from a corona) are presented. We conclude that microquasars are excellent candidates for the parent population of the subgroup of variable low-latitude EGRET sources.

Microquasars (MQs) are galactic objects with relativistic jets that constitute a source population which can be responsible for production of a non-negligible fraction of the observed galactic cosmic rays. These relativistic protons, associated with the termination of the jet, interact with the interstellar medium and, at certain surrounding conditions, may lead to production of detectable fluxes of high-energy and very high-energy gamma-rays. This radiation is accompanied by the broad-band emission of secondary electrons from decays of $\pi^\pm$-mesons produced through synchrotron, bremstrahlung and inverse Compton process. The features of broad-band emission initiated by proton-proton (pp) interactions in such a scenario is discussed in the context of the strategy of search for counterparts of high-energy and very high-energy gamma-ray sources in the galactic plane.

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

The possible associations between the microquasars LS 5039 and LS I +61 303 and the EGRET sources 3EG J1824-1514 and 3EG J0241+6103 suggest that microquasars could also be sources of high-energy gamma-rays. In this work, we present a detailed numerical inverse Compton (IC) model, based on a microquasar scenario, that reproduces the high-energy gamma-ray spectra and variability observed by EGRET for the mentioned sources. Our model considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet that interact through IC scattering with both the radiation and the magnetic fields.

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.

Population studies of unidentified EGRET sources suggest that there exist at least three different populations of galactic gamma-ray sources. One of these populations is formed by young objects distributed along the galactic plane with a strong concentration toward the inner spiral arms of the Galaxy. Variability, spectral and correlation analysis indicate that this population is not homogeneous. In particular, there is a subgroup of sources that display clear variability in their gamma-ray fluxes on timescales from days to months. Following the proposal by Kaufman Bernad\'o et al. (2002), we suggest that this group of sources might be high-mass microquasars, i.e. accreting black holes or neutron stars with relativistic jets and early-type stellar companions. We present detailed inhomogeneous models for the gamma-ray emission of these systems that include both external and synchrotron self-Compton interactions. We have included effects of interactions between the jet and all external photon fields to which it is exposed: companion star, accretion disk, and hot corona. We make broadband calculations to predict the spectral energy distribution of these objects from radio up to GeV energies. The results and predictions can be tested by present and future gamma-ray instruments like INTEGRAL, AGILE, and GLAST.

Microquasars are promising candidates to emit high-energy gamma-rays. Moreover, statistical studies show that variable EGRET sources at low galactic latitudes could be associated with the inner spiral arms. The variable nature and the location in the Galaxy of the high-mass microquasars, concentrated in the galactic plane and within 55 degrees from the galactic center, give to these objects the status of likely counterparts of the variable low-latitude EGRET sources. We consider in this work the two most variable EGRET sources at low-latitudes: 3EG J1828+0142 and 3EG J1735-1500, proposing a microquasar model to explain the EGRET data in consistency with the observations at lower energies (from radio frequencies to soft gamma-rays) within the EGRET error box.

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.

The possible association between the microquasar LS 5039 and the EGRET source 3EG J1824-1514 suggests that microquasars could also be sources of high energy gamma-rays. In this paper, we explore, with a detailed numerical model, if this system can produce the emission detected by EGRET (>100 MeV) through inverse Compton (IC) scattering. Our numerical approach considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet, interacting with both the radiation and the magnetic fields, taking into account the Thomson and Klein-Nishina regimes of interaction. The computed spectrum reproduces the observed spectral characteristics at very high energy.

We explore the possible association between the microquasar LSI+61 303 and the EGRET source 2CG 135+01/3EG J0241$+$6103 by studying, with a detailed numerical model, whether this system can produce the emission and the variability detected by EGRET ($>$100 MeV) through inverse Compton (IC) scattering. Our numerical approach considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet, interacting with both the radiation and the magnetic fields, taking into account the Thomson and Klein-Nishina regimes of interaction. Our results reproduce the observed spectral characteristics and variability at $\gamma$-rays, thus strengthening the identification of LSI+61 303 as a high-energy $\gamma$-ray source.