Source author record

V. Bosch-Ramon

V. Bosch-Ramon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

52works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

52 published item(s)

preprint2022arXiv

3D RMHD simulations of jet-wind interactions in High Mass X-ray Binaries

The interaction of jets in High-Mass X-ray Binaries (HMXBs) with the strong winds driven by the hot companion star in the vicinity of the compact object is fundamental to understand the jet dynamics, non-thermal emission and long-term stability. However, the role of the jet magnetic field in this process is unclear. We study the dynamical role of weak and moderate-to-strong toroidal magnetic fields during the first hundreds of seconds of jet propagation, focusing on the magnetized flow dynamics and the mechanisms of energy conversion. We have developed the code Lóstrego v1.0, a new 3D RMHD code to simulate astrophysical plasmas in Cartesian coordinates. Using this tool, we performed the first 3D RMHD numerical simulations of relativistic magnetized jets propagating through the clumpy stellar wind in a HMXB. The overall morphology and dynamics of weakly magnetized jet models is similar to previous hydrodynamical simulations, where the jet head generates a strong shock in the ambient medium and the initial over-pressure with respect to the stellar wind drives one or more recollimation shocks. In the time scales of our simulations, these jets are ballistic and seem to be more stable against internal instabilities than jets with the same power in the absence of fields. However, moderate-to-strong toroidal magnetic fields favour the development of current-driven instabilities and the disruption of the jet within the binary. A detailed analysis of the energy distribution in the relativistic outflow and the ambient medium reveals that both magnetic and internal energies can contribute to the effective acceleration of the jet. We certify that the jet feedback into the ambient medium is highly dependent on the jet energy distribution at injection, where hotter, more dilute and/or more magnetized jets are more efficient, as anticipated by feedback studies in the case of jets in active galaxies.

preprint2022arXiv

Non-thermal emission in hyper-velocity and semi-relativistic stars

Context. There is a population of runaway stars that move at extremely high speeds with respect to their surroundings. The fast motion and the stellar wind of these stars, plus the wind-medium interaction, can lead to particle acceleration and non-thermal radiation. Aims. We characterise the interaction between the winds of fast runaway stars and their environment, in particular to establish their potential as cosmic-ray accelerators and non-thermal emitters. Methods. We model the hydrodynamics of the interaction between the stellar wind and the surrounding material. We self-consistently calculate the injection and transport of relativistic particles in the bow shock using a multi-zone code, and compute their broadband emission from radio to $γ$-rays. Results. Both the forward and reverse shocks are favourable sites for particle acceleration, although the radiative efficiency of particles is low and therefore the expected fluxes are in general rather faint. Conclusions. We show that high-sensitivity observations in the radio band can be used to detect the non-thermal radiation associated with bow shocks from hypervelocity and semi-relativistic stars. Hypervelocity stars are expected to be modest sources of sub-TeV cosmic rays, accounting perhaps for a $\sim 0.1$% of that of galactic cosmic rays.

preprint2020arXiv

AGN jets versus accretion as reionization sources

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.

preprint2020arXiv

The high-energy emission from HD~93129A near periastron

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 $γ$-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 $γ$-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.

preprint2018arXiv

Constraining very-high-energy and optical emission from FRB 121102 with the MAGIC telescopes

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-$σ$ 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 $σ$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models.

preprint2017arXiv

Search for very-high-energy gamma-ray emission from the microquasar Cygnus X-1 with the MAGIC telescopes

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 $Γ=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.

preprint2016arXiv

The effects of the stellar wind and orbital motion on the jets of high-mass microquasars

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.

preprint2015arXiv

Evidence of coupling between the thermal and nonthermal emission in the gamma-ray binary LS I +61 303

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.

preprint2015arXiv

Non-thermal emission from standing relativistic shocks: an application to red giant winds interacting with AGN jets

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.

preprint2015arXiv

Orbital evolution of colliding star and pulsar winds in 2D and 3D; effects of: dimensionality, EoS, resolution, and grid size

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

preprint2015arXiv

Star-Jet Interactions and Gamma-Ray Outbursts from 3c454.3

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, $δ\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.

preprint2013arXiv

Clumpy stellar winds and high-energy emission in high-mass binaries hosting a young pulsar

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.

preprint2013arXiv

Deep GMRT radio observations and a multi-wavelength study of the region around HESS J1858+020

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.

preprint2013arXiv

Relativistic stellar jets: dynamics and non-thermal radiation

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.

preprint2013arXiv

Star-Jet Interactions and Gamma-Ray Outbursts from 3C454.3

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, $δ\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.

preprint2013arXiv

Studying the non-thermal lobes of IRAS 16547-4247 through a multi-wavelength approach

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.

preprint2012arXiv

Fermi~I particle acceleration in converging flows mediated by magnetic reconnection

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.

preprint2012arXiv

Simulations of stellar/pulsar wind interaction along one full orbit

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.

preprint2011arXiv

A search for Very High Energy gamma-ray emission from Scorpius X-1 with the MAGIC telescopes

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.

preprint2011arXiv

Detection of the gamma-ray binary LS I +61 303 in a low flux state at Very High Energy gamma-rays with the MAGIC Telescopes in 2009

We present very high energy (VHE, E > 100 GeV) γ-ray observations of the γ-ray binary system LS I+61 303 obtained with the MAGIC stereo system between 2009 October and 2010 January. We detect a 6.3σ γ-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.

preprint2011arXiv

Gamma rays from extragalactic astrophysical 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.

preprint2011arXiv

Gamma-ray excess from a stacked sample of high- and intermediate-frequency peaked blazars observed with the MAGIC telescope

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.

preprint2011arXiv

Jet/medium interactions at large scales

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.

preprint2011arXiv

Large-scale flow dynamics and radiation in pulsar gamma-ray binaries

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.

preprint2011arXiv

MAGIC Upper Limits for two Milagro-detected, Bright Fermi Sources in the Region of SNR G65.1+0.6

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.

preprint2011arXiv

Multifrequency Behavior of Microquasars in the GeV--TeV era: A review

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.

preprint2011arXiv

Non-Thermal Emission from Galactic Jets

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.

preprint2011arXiv

On the origin of correlated X-ray/VHE emission from LS I +61 303

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.

preprint2011arXiv

Radio emission from high-mass binaries with non-accreting pulsars

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.

preprint2011arXiv

The termination region of high-mass microquasar jets

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.

preprint2011arXiv

Thermal X-ray Emission from the Shocked Stellar Wind of Pulsar Gamma-ray Binaries

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.

preprint2010arXiv

3D simulations of wind-jet interaction in massive X-ray binaries

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.

preprint2010arXiv

Detection of very high energy gamma-ray emission from the Perseus cluster head-tail galaxy IC 310 by the MAGIC telescopes

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.

preprint2010arXiv

Magic constraints on Gamma-ray emission from Cygnus X-3

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.

preprint2010arXiv

MAGIC Gamma-ray Telescope Observation of the Perseus Cluster of Galaxies: implications for cosmic rays, dark matter, and NGC1275

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.

preprint2010arXiv

MAGIC TeV Gamma-Ray Observations of Markarian 421 during Multiwavelength Campaigns in 2006

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.

preprint2010arXiv

Observations of the Blazar 3C 66A with the MAGIC Telescopes in Stereoscopic Mode

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

preprint2010arXiv

On the interaction of jets with stellar winds in XRBs

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.

preprint2010arXiv

Search for an extended VHE gamma-ray emission from Mrk 421 and Mrk 501 with the MAGIC Telescope

Context: Part of the very high energy $γ$-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.

preprint2010arXiv

Secondary emission behind the radio outflows in gamma-ray binaries?

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.

preprint2010arXiv

Simultaneous multi-frequency observation of the unknown redshift blazar PG 1553+113 in March-April 2008

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

preprint2010arXiv

Spectral energy distribution of gamma-ray binaries: Sources and Processes

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.

preprint2010arXiv

Spectral Energy Distribution of Markarian 501: Quiescent State vs. Extreme Outburst

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

preprint2009arXiv

Massive protostars as gamma-ray sources

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.

preprint2009arXiv

Microquasar interaction with the surrounding medium

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.

preprint2009arXiv

Radio continuum and near-infrared study of the MGRO J2019+37 region

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

preprint2009arXiv

Study of the Spectral and Temporal Characteristics of X-Ray Emission of the Gamma-Ray Binary LS 5039 with Suzaku

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.

preprint2008arXiv

Gamma-ray emission from massive star forming regions

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.

preprint2008arXiv

On the X-ray/TeV connection in Galactic jet sources

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.

preprint2006arXiv

Spectral and variability properties of LS 5039 from radio to very high-energy gamma-rays

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.

preprint2004arXiv

High-energy gamma-ray Emission from Microquasars: LS 5039 and LS I +61 303

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.

preprint2004arXiv

Microquasar models for 3EG J1828+0142 and 3EG J1735-1500

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.