Free Access
Issue
A&A
Volume 575, March 2015
Article Number A109
Number of page(s) 6
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/201425208
Published online 05 March 2015

© ESO, 2015

1. Introduction

Relativistic jets are surrounded by rich environments, and stars and dense medium inhomogeneities can enter into these jets leading to the formation of shocks, as considered for instance in Blandford & Koenigl (1979); Coleman & Bicknell (1985); Komissarov (1994); Bednarek & Protheroe (1997); Araudo et al. (2009, 2010, 2013); Barkov et al. (2010, 2012b); Perucho & Bosch-Ramon (2012), and Bosch-Ramon et al. (2012). These shocks are, only initially in some cases, almost standing shocks, at rest in the laboratory frame (LF). A similar situation can occur in the termination region of other relativistic outflows, such as pulsar winds in binary systems (Paredes-Fortuny et al. 2015).

In several previous papers that modelled the radiation from relativistic shocks at rest in the LF, the emitting region was circumscribed to the closest obstacle vicinity and, as the flow there is only mildly relativistic, the radiation Doppler boosting was not considered. However, there could be more jet energy available for non-thermal particles than assumed in these works, as the effective interaction area may be much larger than just the obstacle section. In addition to that, the flow initially moves at c/ 3 right after the shock, and then is quickly accelerated further downstream because of strong pressure gradients. This effect is seen for instance in Fig. 2 in Bosch-Ramon et al. (2012). Therefore, even relatively close to the obstacle at rest, the shocked flow velocity is high enough for Doppler boosting to significantly enhance radiation in the emitting fluid direction of motion. The motion of the shocked flow could be also considered for an obstacle that is itself accelerated by the jet ram pressure (e.g. Barkov et al. 2012a; Khangulyan et al. 2013), but this would require accounting for the complex dynamical evolution of the obstacle (see Bosch-Ramon et al. 2012, for hydrodynamical simulations predicting the obstacle disruption).

In this work, we study the kinetic-to-internal (possibly non-thermal) energy conversion, gamma-ray radiation efficiency, and Doppler boosting effects, in the shock and the shocked flow of a relativistic hydrodynamic jet interacting with an obstacle at rest. The obstacle could be a star or a cloud in the case of an extragalactic jet, or a stellar wind inhomogeneity in the case of a galactic microquasar jet. Since this is a relevant and illustrative case, we focus here on the case of a population of red giant stars (RG) interacting with the jet of an active galactic nucleus (AGN). We first characterize the non-thermal luminosity from the shocked jet region accounting for hydrodynamical information, and second we adopt a phenomenological description of the spatial distribution of the RG in the regions crossed by the jet as it goes through the galaxy. To compute the emission, we consider mainly inverse Compton (IC) scattering, assume that synchrotron radiation is a minor emission channel, and put special emphasis on the GeV–TeV range. Studies such as this are needed to characterize the non-thermal emission from interactions between stars and jets in AGN, which is important if the AGN gamma-ray radiation is to be understood (e.g. Barkov et al. 2012a; Khangulyan et al. 2013).

Concerning other types of sources, the characterization of the radiation properties from jet-obstacle interactions could also shed light on the main factors shaping jet variability at high energies in microquasars (e.g. Owocki et al. 2009), or the post-shock emission in pulsar-wind termination shocks in some binary systems (e.g. Bosch-Ramon 2013). Regarding jets with strong magnetic fields with polarity changes, it has been proposed that the interaction with obstacles may lead to magnetic reconnection and efficient non-thermal emission (see Bosch-Ramon 2012). This kind of interaction is very complex given the strong dynamical role of the magnetic field and its non-trivial topological properties and deserves a specific treatment, which is beyond the scope of this paper.

2. The jet-RG wind interaction

2.1. Physical scenario

The study presented here is based on the interaction of a jet with an RG surrounded by its dense and slow wind. The basic hydrodynamical ingredients of this study are: a highly relativistic, supersonic flow with parallel fluid lines on the scales of interest; a cold obstacle at rest; and a fraction of the internal energy downstream of the standing shock being in the form of non-thermal particles. The initial speed of the RG is taken to be c if it is to be approximated as at rest, so the problem can be considered axisymmetric around an axis crossing the star and with the jet flow direction. In addition, the RG wind must have enough inertia to stand the jet impact while it penetrates into the jet, which can be expressed as ρrgvrg2>Lj/πRj2c\hbox{$\rho_{\rm rg}v_{\rm rg}^2>L_{\rm j}/\pi R_{\rm j}^2c$}, where ρrg and vrg are the RG wind density and velocity perpendicular to the jet, and Lj and Rj are the jet power and radius, respectively. In addition, in the standing shock scenario the velocity acquired by the RG wind due to jet impact must be c. If these conditions are fulfilled, the problem becomes two dimensional, and the steady state of the physical system can be taken as the steady state of the shocked jet flow, which significantly lightens the hydrodynamical calculations presented in what follows.

2.2. Simulations

The jet-RG wind interaction was simulated in 2D assuming axisymmetry. We solved the equations of relativistic hydrodynamics using a Marquina Riemann solver (Donat & Marquina 1996), second-order spatial reconstruction scheme (minbee slope limiter, plus an additional dissipation term from Colella 1985), and an adiabatic gas with index \hbox{$\hat\gamma=4/3$}. The resolution is 600 cells in the vertical direction, the z-axis, and 300 in the radial direction, the r-axis, and the physical size is zmax = 2 × 1015 cm in the z-direction, and rmax = 1015 cm in the r-direction. The jet power injected in the grid is 1037 erg s-1, with a Lorentz factor γj = 10. We fix the RG wind radius to Rrg = 1014 cm. A snapshot of the density distribution in the steady state is presented in Fig. 1. The shocked jet material is clearly seen in light blue, forming a cometary tail pointing upwards and surrounding the RG wind, red with green borders.

For simplicity, we modelled the RG wind region as a sphere with a density 108 times that of the jet in the FF. As the shocked RG wind is much slower than the shocked jet flow, we considered it static, and thus the steady state of the simulation is quickly reached. In a more realistic setup, there would be more wind momentum flux perpendicular to the jet direction, slightly widening the shocked flow structure, but this setup would have required a simulation ~103 times longer.

We chose the grid resolution such that numerical viscosity played a negligible role in the generation of internal energy in the shocked flow. There is some mixing in the RG wind-shocked jet flow boundary, but it has a minor impact on the radiation estimates presented below. This mixing slows down the shocked jet flow close to the axis, reducing the effect of Doppler boosting and thus weakening the emission from those (small) regions.

thumbnail Fig. 1

Distribution of density in the flow frame in cgs units.

For illustrative purposes, the simulation result is interpreted as an AGN jet of Lj ≈ 1044 erg s-1 and Rj = 1 pc, much larger than the grid, interacting with an RG wind region whose size is derived assuming jet-stellar wind pressure balance Rrgwvwc/4Lj×Rj1014cm,\begin{equation} R_{\rm rg}\approx \sqrt{\dot{M}_{\rm w}v_{\rm w}c/4L_{\rm j}}\times R_{\rm j}\approx 10^{14}\,{\rm cm}, \label{ro} \end{equation}(1)for a stellar mass loss w = 10-8M/yr and a wind velocity vw ≈ 2 × 107 cm s-1, typical of an RG. Note that for a real RG wind the actual size of the obstacle would be slightly larger than for the case of an obstacle of very large inertia, as the shocked wind region should be added in Eq. (1).

Relativistic jets usually present half semi-opening angles of several degrees, αj ~ 0.1 rad, which for a supersonic jet without lateral causal contact yields γj = 1/αj ~ 10 (see Khangulyan et al. 2013, and references therein). This condition implies an interaction height of zj = γjRj ~ 10 pc.

3. Non-thermal emission

We use the hydrodynamical information obtained in Sect. 2 to derive some important physical properties that allow for the characterization of the non-thermal emission. These properties are the radiative efficiency, the distribution of the non-thermal energy, the Doppler boosting, and the effective section of the RG wind as an obstacle for the jet.

3.1. Radiative efficiency and non-thermal energy distribution

To estimate the non-thermal emission, the radiative efficiency of the shocked flow is needed. This efficiency can be quantified as the radiative (1 /trad) to total (1 /tc) cooling rate ratio (frad): frad=trad-1/tc-1=trad-1/(trad-1+tnrad-1)=(1+trad/tnrad)-1,\begin{equation} f_{\rm rad}=t_{\rm rad}^{-1}/t_{\rm c}^{-1}=t_{\rm rad}^{-1}/(t_{\rm rad}^{-1}+t_{\rm nrad}^{-1})=(1+t_{\rm rad}/t_{\rm nrad})^{-1}, \end{equation}(2)where 1 /tc also takes non-radiative cooling, 1 /tnrad, into account. The timescale for non-radiative cooling, i.e. particle escape and adiabatic losses, is estimated in the flow frame (FF) as the distance between each point to the RG location (drg) divided by the local velocity (vsh) times the shocked flow Lorentz factor (γsh): tnrad = drg/vshγsh. For the radiative timescale we adopt: trad = CICu-1E-1, where CIC ≈ 26 in cgs units. This timescale is derived assuming IC in Thomson, and taking a reference electron energy of EIC=me2c4/ϵ0\hbox{$E_{\rm IC}=m_{\rm e}^2c^4/\epsilon_0$}, where EIC is roughly the energy of the IC cross-section maximum, and ϵ0 is the target photon energy and u the energy density of a photon field coming from the RG, all in the FF. An effective stellar temperature of 3000 K has been taken, typical for an RG. The star luminosity, needed to obtain u, has been determined through L=vwwc\hbox{$L_*=\dot{M}_{\rm w}v_{\rm w}c$} (LF). This is a rather conservative assumption, as in fact only a fraction of all the star radiation momentum goes to the stellar wind. Note that protons are not taken into account as emitting particles.

Figure 2 shows the distribution of frad in the shocked flow. As seen in the figure, for the adopted jet and RG wind properties, frad is everywhere below one for the case simulated, i.e. the flow is adiabatic. In this case, the thermal pressure (P) of the shocked flow is a good proxy for the non-thermal energy distribution, which can be taken to be proportional to P all through the shocked flow. In Fig. 3, the distribution of P in the shocked flow is shown in a colour scale, from red to light green.

The characterization of the radiative and the non-radiative timescales, and of the non-thermal energy distribution, allows for the computation of an approximate total luminosity for IC in the FF, LIC, simplifying the IC emissivity as that at EIC. For moderately magnetized shocked flows, the synchrotron luminosity Lsync may be higher than LIC, as LsyncEmax3α\hbox{$L_{\rm sync}\propto E_{\rm max}^{3-\alpha}$}, where α is the power-law index of the particle energy distribution, and Emax is the maximum electron energy, which could be EIC. Under efficient particle acceleration and γj ≳ 10, the synchrotron spectrum may contribute to or dominate in the gamma-ray band. However, as the magnetic field is not well constrained, we will assume here that the synchrotron component can be neglected.

The adopted radiation cooling model is simple, although it accounts for the basic IC cross-section properties, and the target density behaviour under relativistic motion of the emitter. Therefore, the model should yield reasonable luminosity predictions for IC gamma rays, if gamma-ray pair-creation absorption is negligible.

thumbnail Fig. 2

Distribution of radiation efficiency for the shocked jet region.

3.2. Doppler boosting and obstacle effective radius

The radiation enhancement due to Doppler boosting is δsh4=1/γsh4(1cosθshβsh)4,\begin{equation} \delta_{\rm sh}^4=1/\gamma_{\rm sh}^4(1-\cos \theta_{\rm sh}\beta_{\rm sh})^4, \end{equation}(3)which can be approximated as 16γsh4\hbox{$ 16\gamma_{\rm sh}^4$} when the angle between the jet and the observer directions is θobs ≪ 1 /γsh.

Figure 4 shows the distribution of the factor of enhancement of the emission from the shocked flow because of Doppler boosting. As seen in Fig. 4, values of δsh4\hbox{$\delta_{\rm sh}^4$} of a few are found right after the shock, but further downstream δsh4\hbox{$\delta_{\rm sh}^4$} achieves values as high as several thousand.

thumbnail Fig. 3

Distribution of pressure in cgs units for the shocked jet region.

thumbnail Fig. 4

Distribution of the Doppler boosting enhancement (δsh4\hbox{$\delta_{\rm sh}^4$}) for the shocked jet region.

Another important result from our simulation is derived when computing \hbox{$\dot{U}_{\rm sh}$}, the rate of the total thermal energy generated in the interaction region in the LF, where the total thermal energy is the sum of internal energy plus the pressure component. The value of \hbox{$\dot{U}_{\rm sh}$} has been computed as the fraction of Lj that escapes the grid in the form of thermal energy. From this, we can derive the effective radius (Reff) of the RG wind when acting as an obstacle for the jet sh=(Reff/Rj)2Lj4×1036ergs-1,\begin{equation} \dot{U}_{\rm sh}=(R_{\rm eff}/R_{\rm j})^2L_{\rm j}\approx 4\times 10^{36}\,{\rm erg~s}^{-1}, \end{equation}(4)which turns out to be Reff ≈ 6 Rrg because of effective kinetic to internal energy conversion even in regions far from the RG. This effective conversion manifests itself locally as well, as shown by Fig. 5.

The value of Reff is actually dependent on Rrg and zmax. To check this dependence, two additional simulations were carried out, both with twice the cell size of the presented simulation, one with the same grid size, and the other with a grid that is twice as large, and both with the same remaining parameters. These simulations show that Reff grows slowly with rmax. This is expected to happen as the oblique shock becomes weaker farther away from the RG, but the actual weakening of the shock is sensitive to the specific value of the unshocked flow internal energy, which is unknown. Note also that to derive Reff all the grid is taken into account; in reality, the efficiency of energy conversion is larger if the grid outermost region without shock presence is neglected. Therefore, Reff = 6 Rrg may be considered a lower limit.

thumbnail Fig. 5

Distribution of the ratio of fluxes of thermal energy to kinetic plus thermal energy.

4. The IC luminosity from a jet-obstacle interaction

4.1. The IC luminosity to the observer

Using the hydrodynamics results and under the simplified emission model described, we can derive an approximate observer IC luminosity of the shocked flow within the grid from LIC,obs~shδsh4(1cosθobs)ηNTγ̂P(γ̂1)-1trad-1dvsh,\begin{equation} L_{\rm IC,obs}\sim\int_{\rm sh}\delta_{\rm sh}^4(1-\cos{\theta_{\rm *obs}})\eta_{\rm NT}\hat\gamma P(\hat\gamma-1)^{-1}t_{\rm rad}^{-1}{\rm d}v_{\rm sh}, \label{icobs} \end{equation}(5)where dvsh is the shocked flow differential volume, and θ∗ obs the angle between the stellar photons and the observer direction, both in the FF. The term (1 − cosθ∗ obs) is a term related to the IC scattering probability-dependence on the stellar photon and observer directions. The quantity ηNT is the ratio of non-thermal to total thermal energy density; it is a free parameter of value between 0 and 1. For a typical accelerated electron population, say with a particle energy distribution harder than E-3, given the value of ϵ0, and accounting for Doppler boosting, the IC radiation will be mostly radiated in the GeV–TeV range, as seen by the observer.

From Eq. (5), the total observer IC luminosity is computed, obtaining LIC,obs ~ 1035ηNT erg s-1, which is ≈ 1 ηNT% of the jet luminosity injected in the grid, and a 3 ηNT% of \hbox{$\dot{U}_{\rm sh}$}. On the other hand, neglecting relativistic effects and confining the emitter to the RG wind vicinity, one would get LIC,obs ~ 3 × 1032ηNT erg s-1, which is ~300 times smaller than the more realistic value. This shows that the radiation can be underestimated if the interaction region is not globally taken into account.

Figure 6 shows the structure of the observer IC emitter for θobs ≪ 1 /γsh, produced in ring-like cells of volume 2πrjΔrjΔzj, where Δrj = rmax/ 300 and Δzj = zmax/ 600. As seen in the figure, most of the emission comes from the closest regions to the RG, but also from farther away, close to the jet shock.

thumbnail Fig. 6

Distribution of the observer IC luminosity in cgs units for the shocked jet region.

4.2. Generalization of the results

For the adopted approximation zj ~ γjRj, our results correspond to an emitter located at relatively high heights (~10 pc) in a jet of moderate power (~1044 erg s-1). These results can nevertheless be generalized through LICtnrad/trad~Ljzj-1,\begin{equation} L_{\rm IC}\propto t_{\rm nrad}/t_{\rm rad}\sim L_{\rm j}z_{\rm j}^{-1}, \label{coolz} \end{equation}(6)so LIC,obs will grow closer to the jet base and/or with a more powerful jet. In saturation, Eq. (6) is not valid anymore as tradtnrad, and thus Eq. (5) becomes LIC,obsmax~102ηNTsh4×1038ηNTergs-1,\begin{equation} L_{\rm IC,obs}^{\rm max}\sim 10^2\eta_{\rm NT}\dot{U}_{\rm sh}\approx 4\times 10^{38}\eta_{\rm NT}\,{\rm erg~s}^{-1}, \end{equation}(7)i.e. it gets constant, for the same remaining parameters: γj and vw.

The value of LIC,obs is also γj2\hbox{$ \gamma_{\rm j}^2$}, for θobs ≪ 1 /γsh and γj,sh ≫ 1. This relation comes from the relation γshγj, strictly valid in the highly relativistic case, and the link between the IC luminosity in the LF (LIC,lab) and LIC,obs: LIC,obs/LIC,labδsh4/γsh216γsh2γj2.\begin{equation} L_{\rm IC,obs}/L_{\rm IC,lab}\propto\ \delta_{\rm sh}^4/\gamma_{\rm sh}^2\approx 16\,\gamma_{\rm sh}^2\propto \gamma_{\rm j}^2. \label{dopsh} \end{equation}(8)Therefore, the computed estimates can be analytically generalized through LIC,obsγj2Rrg2Ljzj-1γj2vwLjzj-1,\begin{equation} L_{\rm IC,obs}\propto \gamma_{\rm j}^2R_{\rm rg}^2L_{\rm j}z_{\rm j}^{-1}\propto \gamma_{\rm j}^2\dot{M}v_{\rm w}L_{\rm j}z_{\rm j}^{-1}, \end{equation}(9)or in saturation, through LIC,obsγj2vw,\begin{equation} L_{\rm IC,obs}\propto \gamma_{\rm j}^2\dot{M}v_{\rm w}, \end{equation}(10)relations that also account for different types of winds via vw.

For a microquasar jet, given the typical clump sizes (e.g. Owocki et al. 2009), clumps with Reff ~ Rj may reach the jet. Also, IC emission may take place in saturation given the proximity of the stellar companion, yielding luminosities as high asLIC,obs~16γsh2ηNTLj,\begin{equation} L_{\rm IC,obs}\sim 16\gamma_{\rm sh}^2\eta_{\rm NT}L_{\rm j}, \end{equation}(11)which could be Lj for some relatively exceptional, and brief, events.

5. A population of red giants within an AGN jet

We explore in what follows whether collectively RG may be relevant for the non-thermal emission from AGN jets. About a 5% of the stars in the AGN vicinity are expected to be RG (Syer & Ulmer 1999). In addition, there are observational hints of stars interacting with AGN jets (e.g. Hardcastle et al. 2003; Müller et al. 2014). Therefore, the scenario explored in this work, and already postulated as the origin of high-energy emission by several authors (see Sect. 1), could be actually a rather common phenomenon.

We take the nearby, gamma-ray emitting AGN M87 as an example. With ~1011 stars in the inner ~1 kpc of the galaxy (see Fig. 7 in Gebhardt & Thomas 2009), the number of RG within one jet should be about NRG~0.05×1011×3αj2/44×107,\begin{equation} N_{\rm RG}\sim 0.05\times 10^{11}\times 3\alpha_{\rm j}^2/4\approx 4\times 10^7, \end{equation}(12)where we take αj ~ 0.1, and 3αj2/4\hbox{$3\alpha_{\rm j}^2/4$} is the jet to total volume ratio.

Assuming an RG density nRG1/zjξ\hbox{$n_{\rm RG}\propto 1/z_{\rm j}^\xi$}, an average vww~1025\hbox{$\dot{M}_{\rm w}v_{\rm w}\sim 10^{25}$} g cm s-2, the relation from Eq. (6), negligible synchrotron radiation, a constant γj = 10, and Lj = 1044 erg s-1, each interaction would produce in gamma rays LIC,obs(zj)~LIC,obs(zmin)×(zmin/zj),\begin{equation} L_{\rm IC,obs}(z_{\rm j})\sim L_{\rm IC,obs}(z_{\rm min})\times (z_{\rm min}/z_{\rm j}), \end{equation}(13)where zmin ≈ 8 × 1015 cm is found through LIC,obs(10 pc) = 1035ηNT erg s-1, and fixing the zmin-luminosity to the saturation value: LIC,obs(zmin)=LIC,obsmax\hbox{$L_{\rm IC,obs}(z_{\rm min})=L_{\rm IC,obs}^{\rm max}$}. All this allows us to derive LIC,obstot==NRG(3ξ2ξ)(zmax2ξzmin2ξzmax3ξzmin3ξ)LIC,obs(zmin)×zmin,\begin{eqnarray} L_{\rm IC,obs}^{\rm tot}&=&\int_{\rm z_{\rm min}}^{\rm z_{\rm max}} n_{\rm RG}(z_{\rm j})L_{\rm IC,obs}(z_{\rm j})\,{\rm d}v_{\rm j} \label{lmaxic}\notag \\[2mm] &=&N_{\rm RG}\left(\frac{3-\xi}{2-\xi}\right)\left(\frac{z_{\rm max}^{2-\xi}-z_{\rm min}^{2-\xi}}{z_{\rm max}^{3-\xi}-z_{\rm min}^{3-\xi}}\right)L_{\rm IC,obs}(z_{\rm min})\times z_{\rm min}, \end{eqnarray}(14)where dvj is the jet differential volume.

The dominant contribution to LIC,obstot\hbox{$L_{\rm IC,obs}^{\rm tot}$} will come from the inner-jet regions, or from galactic scales, depending on whether ξ is > or < than 2. For ξ ≈ 2, the zmin,max-dependence of LIC,obstot\hbox{$L_{\rm IC,obs}^{\rm tot}$} is weak, but the external jet regions dominate. For this ξ-value1, Eq. (14) gives LIC,obsmax~5×1041ηNT\hbox{$L_{\rm IC,obs}^{\rm max}\sim 5\times 10^{41}\eta_{\rm NT}$} erg s-1, implying that the contribution from all RG within 1 kpc could be up to a factor of a few above the persistent GeV–TeV luminosities in M87 (e.g. Abdo et al. 2009; Aleksić et al. 2012). It seems unlikely however that the acceleration efficiency will be ηNT ~ 1, and we may be overestimating the average value of vw. On the other hand, we do not take into account the asymptotic giant branch (AGB) and massive star contribution, both of which should be added to that of RG. Non-stellar objects such as clouds could also be relevant. Finally, the value of L could be significantly higher, at least in some stars. In conclusion, LIC,obsmax\hbox{$L_{\rm IC,obs}^{\rm max}$} may explain the persistent GeV–TeV emission in M87, and could be significant for AGN jets in general.

Note that the viewing angle with respect to the M87 jet is not zero, but rather ~20° (e.g. Acciari et al. 2009). However, given the obtained values of γsh (4), an angle of 20° is small enough for the condition θobs ≪ 1 /γsh to marginally apply. Here, the gamma-ray luminosity for M87 has been computed using the proper viewing angle, although taking θobs = 0 yields very similar results.

6. Discussion

The emission coming from the jet in Blazar AGN will be more boosted than that from the downstream region of a standing shock. Thus, strong jet activity could mask the putative jet-obstacle emission. However, for slightly misaligned, relatively nearby AGN, or those without significant non-thermal activity intrinsic of the jet, the contribution from populations of stars interacting with the jet could be important, ~0.11% of Lj (but strongly dependent on ηNT, γj, and θobs). This collective emission should be constant unless ξ were significantly softer than 2, which does not seem realistic. There could also be rare flaring events when a dense cloud, an AGB, a Wolf-Rayet or a luminous blue variable (WR/LBV) star, entered into the jet close to its base (e.g. Araudo et al. 2010, 2013; Barkov et al. 2010; Hubbard & Blackman 2006). In the case of microquasar jets, the ratio Reff/Rj can be larger than in the extragalactic case and fluctuations may be more common, with the arrival of larger clumps leading to brighter events (Owocki et al. 2009; Araudo et al. 2009). This may also happen in binary systems hosting a non-accreting pulsar (see Bosch-Ramon 2013; Paredes-Fortuny et al. 2015, and references therein).

The obtained numerical results have been used to characterize the effective section of the RG wind as an obstacle for the jet, the distribution of the non-thermal particles through the constant ηNT, and the effect of Doppler boosting in the shocked flow for γj ≫ 1. Previous research focussed on a region downstream of the jet shock next to the obstacle and with a similar size (e.g. Araudo et al. 2009, 2010; Barkov et al. 2010, 2012b, and Doppler boosting effects were not considered as the flow is mildly relativistic there. As shown here, both approximations can lead to the underestimation of the non-thermal emission from jet-obstacle interactions.

Our work has some caveats. The simulation, in axisymmetry, without thermal cooling nor magnetic field, with modest resolution, and with a simple prescription of the RG wind, is not very realistic, although it captures the basic features of the shocked region of a hydrodynamic jet. The value of frad is a rough characterization of the radiative efficiency for IC, and our model of the IC emission is very simplified: it assumes IC only in Thomson, and accounts only for the electron energy at the maximum of the IC cross section. This simplification introduces an uncertainty in the predicted luminosity by a factor of a few. On the other hand, our radiation calculations account for the main features of IC scattering, and so we expect our results to be correct within an order of magnitude. Note that the value ηNT cannot be derived from first principles, and despite ηNT ≳ 0.1 are not uncommon in relativistic jet sources, ηNT ≪ 1 are also possible. It is also difficult to constrain the magnetic field strength in jets, which affects the calculation of synchrotron radiation. More accurate calculations will be presented in future work.


1

A slope ξ ~ 2 can be derived from Gebhardt & Thomas (2009) for kpc scales, becoming harder closer to the centre of M87.

Acknowledgments

V.B.-R. acknowledges support by the Spanish Ministerio de Economía y Competitividad (MINECO) under grant AYA2013-47447-C3-1-P. This research has been supported by the Marie Curie Career Integration Grant 321520. V.B-R. also acknowledges financial support from MINECO and European Social Funds through a Ramón y Cajal fellowship.

References

  1. Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, ApJ, 707, 55 [NASA ADS] [CrossRef] [Google Scholar]
  2. Acciari, V. A., Aliu, E., Arlen, T., et al. 2009, Science, 325, 444 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  3. Aleksić, J., Alvarez, E. A., Antonelli, L. A., et al. 2012, A&A, 544, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Araudo, A. T., Bosch-Ramon, V., & Romero, G. E. 2009, A&A, 503, 673 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Araudo, A. T., Bosch-Ramon, V., & Romero, G. E. 2010, A&A, 522, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Araudo, A. T., Bosch-Ramon, V., & Romero, G. E. 2013, MNRAS, 436, 3626 [NASA ADS] [CrossRef] [Google Scholar]
  7. Barkov, M. V., Aharonian, F. A., & Bosch-Ramon, V. 2010, ApJ, 724, 1517 [NASA ADS] [CrossRef] [Google Scholar]
  8. Barkov, M. V., Aharonian, F. A., Bogovalov, S. V., Kelner, S. R., & Khangulyan, D. 2012a, ApJ, 749, 119 [NASA ADS] [CrossRef] [Google Scholar]
  9. Barkov, M. V., Bosch-Ramon, V., & Aharonian, F. A. 2012b, ApJ, 755, 170 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bednarek, W., & Protheroe, R. J. 1997, MNRAS, 287, L9 [NASA ADS] [CrossRef] [Google Scholar]
  11. Blandford, R. D., & Koenigl, A. 1979, Astrophys. Lett., 20, 15 [NASA ADS] [Google Scholar]
  12. Bosch-Ramon, V. 2012, A&A, 542, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Bosch-Ramon, V. 2013, A&A, 560, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Bosch-Ramon, V., Perucho, M., & Barkov, M. V. 2012, A&A, 539, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Colella, P. 1985, SIAM J. Sci. Stat. Comput., 6, 104 [CrossRef] [Google Scholar]
  16. Coleman, C. S., & Bicknell, G. V. 1985, MNRAS, 214, 337 [NASA ADS] [CrossRef] [Google Scholar]
  17. Donat, R., & Marquina, A. 1996, J. Comput. Phys., 125, 42 [NASA ADS] [CrossRef] [Google Scholar]
  18. Gebhardt, K., & Thomas, J. 2009, ApJ, 700, 1690 [NASA ADS] [CrossRef] [Google Scholar]
  19. Hardcastle, M. J., Worrall, D. M., Kraft, R. P., et al. 2003, ApJ, 593, 169 [NASA ADS] [CrossRef] [Google Scholar]
  20. Hubbard, A., & Blackman, E. G. 2006, MNRAS, 371, 1717 [NASA ADS] [CrossRef] [Google Scholar]
  21. Khangulyan, D. V., Barkov, M. V., Bosch-Ramon, V., Aharonian, F. A., & Dorodnitsyn, A. V. 2013, ApJ, 774, 113 [NASA ADS] [CrossRef] [Google Scholar]
  22. Komissarov, S. S. 1994, MNRAS, 269, 394 [NASA ADS] [CrossRef] [Google Scholar]
  23. Müller, C., Kadler, M., Ojha, R., et al. 2014, A&A, 569, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Owocki, S. P., Romero, G. E., Townsend, R. H. D., & Araudo, A. T. 2009, ApJ, 696, 690 [NASA ADS] [CrossRef] [Google Scholar]
  25. Paredes-Fortuny, X., Bosch-Ramon, V., Perucho, M., & Ribó, M. 2015, A&A, 574, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Perucho, M., & Bosch-Ramon, V. 2012, A&A, 539, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Syer, D., & Ulmer, A. 1999, MNRAS, 306, 35 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1

Distribution of density in the flow frame in cgs units.

In the text
thumbnail Fig. 2

Distribution of radiation efficiency for the shocked jet region.

In the text
thumbnail Fig. 3

Distribution of pressure in cgs units for the shocked jet region.

In the text
thumbnail Fig. 4

Distribution of the Doppler boosting enhancement (δsh4\hbox{$\delta_{\rm sh}^4$}) for the shocked jet region.

In the text
thumbnail Fig. 5

Distribution of the ratio of fluxes of thermal energy to kinetic plus thermal energy.

In the text
thumbnail Fig. 6

Distribution of the observer IC luminosity in cgs units for the shocked jet region.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.