On the evolution of the particle distribution and the cascade in a moving, expanding emission region in blazar jets

There is a large variety in the models explaining blazar flares. Here, we study the flare profile induced by a moving and expanding blob with special emphasize on the gamma-gamma pair production. We first develop a simple semi-analytical model to study the evolution of the particle distribution in the expanding blob and show the influence of the pair production. In a second step, we produce a realistic simulation using the OneHaLe code based upon parameters of PKS 1510-089. The semi-analytical model shows that the pair production significantly influences the flare evolution, while the opening angle and the expansion can prolong flares considerably. The simulation based on PKS 1510-089 indicate that flares of a moving expanding blob result in strongly wavelength dependant light curves including delayed, secondary flares. A moving, expanding blob can cause significant flaring events with a large variety in light curve profiles. High-cadence multiwavelength observations are necessary to derive the details causing the flare. Extended observations beyond the initial burst may provide important information on the opening angle and the particle content due to delayed secondary flares in some energy bands.


Introduction
Blazar emission (with blazars defined as active galaxies whose jets point toward the Earth) is typically explained on the basis of the so-called one-zone model, where a single zone is responsible for most of the source's radiative output.The spectral energy distribution (SED) is characterized by two broad humps.The lowenergy one is explained by electron-synchrotron emission, while the high-energy hump can be explained by inverse-Compton emission or hadronically induced processes, such as proton synchrotron or synchrotron emission from the leptonic cascade.An important role is played by photon fields external to the jet, such as the accretion disk (AD), the broad-line region (BLR), and the dusty torus (DT), as these fields may provide ample seed photons for particle-photon and photon-photon interactions (see, e.g., Böttcher 2019; Cerruti 2020, for detailed reviews).The one-zone model is well justified in flares, where the variability timescale restricts the size of the emission region.However, while the particle flow in this region is relativistic, the emission region itself is typically assumed to remain stationary with respect to the black hole (e.g., H.E.S.S. Collaboration 2019), which allows us to ignore certain complications, such as varying (external) photon fields, adiabatic expansion and cooling, and so on.On the other hand, the observation of stationary features in radio very-long-baseline interferometry (VLBI) observations (and other wavelengths where jets have been resolved) suggests Now at: Landessternwarte, Universität Heidelberg, Königstuhl 12, 69117 Heidelberg, Germany.standing recollimation shocks (e.g., Weaver et al. 2022), where particles may be accelerated and emit radiation.In such a situation, the emission region would indeed be stationary with respect to the black hole.
The same radio VLBI observations have further revealed moving components launched somewhere upstream with respect to the radio core.The interaction of such moving features with the standing features has been connected with multiwavelength flaring events (e.g., Ahnen et al. 2017; H.E.S.S. Collaboration 2021).Indeed, numerical simulations by Fichet de Clairfontaine et al. (2021Clairfontaine et al. ( , 2022) ) have shown that the interaction of moving and standing shocks can induce rapid flaring events.Thus, it is important to study the flaring characteristics of a moving, expanding emission region (or "blob").
Recently, Boula & Mastichiadis (2022) and Tramacere et al. (2022) discussed this model in detail and derived the expected time delays between the γ-ray and the radio band expected from the progressive optical thinning at radio frequencies of the expanding source (see also Saito et al. 2015).In the present work, we specifically consider the effect of γ−γ pair production on the evolution of the particle distribution and the photon fluxes.In particular, bright external photon fields can have a major influence, as they provide a significant amount of absorption for the γ rays of the emission region (e.g., Zacharias 2021).However, as these photon fields are only present in the direct vicinity of the black hole, the optical thickness for γ rays changes over time, as the blob moves down the jet.Additionally, the opening angle will strongly influence the escape time of both photons and particles owing to the increase in the radius of the emission region.This could potentially prolong the pair production as the (internal) photons have more time to interact before escape.This naturally competes with the thinning of the photon density due to expansion.In any case, it is possible to obtain variability simply via the motion of an expanding blob through the various radiation fields due to the varying injection of pairs.Firstly, we derive a simple semi-analytical model that describes the time-dependent evolution of the particle distribution, given in Sect. 2. In this section, we also introduce the basic assumptions about the escape timescale and the implications of the expansion.We then derive the particle distribution without pair injection, followed by a very simple linear cascade model for a few exemplary cases.In Sect.3, we use the numerical code OneHaLe (Zacharias 2021;Zacharias et al. 2022) to derive the light curves of a realistic simulation based on the parameters of PKS 1510−089, which is known for its bright external photon fields.We present our conclusions in Sect. 4.

A simple electron evolution model
We derived simple semi-analytical models for the timedependent evolution of the electron distribution and while ignoring all energy dependencies.This is indeed a major simplification, however, it allows us to study three distinct cases with a linear cascade (at most).Using the time-dependent, one-zone radiation code OneHaLe, we reproduced two of the three cases well.We discuss the failure of the third case, as well as the influence of the choice of parameters on the solutions.

Escape time
The blob travels with constant speed, β Γ c, corresponding to the bulk Lorentz factor Γ = (1 − β 2 Γ ) −1/2 along the z-axis.Within a conical jet, the opening angle is constant.The VLBI observations suggest an opening angle of ∝α/Γ with α ∼ 0.26 (Pushkarev et al. 2017), but we used α as a free parameter.The radius, R, of the blob thus evolves as a function of comoving time 1 , t, and jet coordinate, z : where, in the last line, we present approximations for small angles and β Γ ≈ 1. Boula & Mastichiadis (2022) expressed this equation by an "expansion speed", u exp , which relates to our equation as u exp = αc.In their paper, they used α between 0.01 and 0.2.It should be noted that one of the common estimates for jet expansion, α = 1 (e.g., Ackermann et al. 2016), implies a radial expansion with the speed of light, while α > 1 implies a superluminal expansion and thus a causal disconnection of regions within the blob or jet.On average, particles escape the emission region on an energy-independent timescale, with η esc > 1.This resembles an advective motion of the plasma below the speed of light mimicking the trapping of charged particles in the magnetized blob.Calculating the average escape timescale for photons, we obtain η esc,ph = 3/4 (Böttcher & Chiang 2002).
In an expanding blob, the escape time increases.In turn, particles and photons take longer to escape and can only efficiently do so once the intrinsic time, t, since launch surpasses ∼t esc (t); namely, with t esc (0) = η esc R 0 /c.In the case of η esc α → 1, particles are effectively trapped in the blob without a chance for a meaningful escape.In this case, the escape of photons is also significantly slowed down.Hence, a significant cascade could still materialize at far distances from the black hole.It also suggests that a cascade that has begun developing closer to the black hole (say, within the external photon fields) can continue to grow even at far distances.

Particle density evolution without secondary injection
We are interested in the time-dependent evolution of the total particle density, while neglecting its detailed energy-dependent evolution, as the latter does not change the overall density.Considering only time-dependent primary injection and escape, but no secondary injections, the kinetic equation for the particle distribution n(t) becomes: with the analytical solution: The particle injection rate Q(t) is coupled to the particle injection luminosity L inj (t), taking the form: where E(γ) is a function of the injection energy spectrum, which is of no concern to us, while V(t) is the volume of the spherical blob.With the assumption and the definition q 0 := L 0 /[V(0)E(γ)], the second equality of Eq. ( 6) can readily be achieved.The power-law index, p, describes the decrease of the injection luminosity as a function of R(t).As the jet power is proportional to R(t) 2 , p = 2 implies a constant jet particle injection power.While p is a free parameter, we will mostly use p = 2 below.For constant injection and escape -that is a straight jet, α = 0 -the solution of Eq. ( 5) is n = q 0 t esc (0), as expected.
This equation is positive for all times and p.For p = (η esc α) −1 − 2, the density becomes Equation ( 7) holds for all particle species as long as there is no secondary injection or particle destruction.Below, we refer to Eq. ( 7) as the "standard solution".Figure 1 exemplifies Eq. ( 7) for two values of p and various values of η esc α.The parameters q 0 and R 0 are chosen to resemble typical one-zone parameters, while η esc = 3 implies α ∈ [0.0 3, 0.3], resembling the typical range of this parameter (Pushkarev et al. 2017;Boula & Mastichiadis 2022).Interestingly, both p and η esc α have a similar effect.A larger p implies a faster decrease in the injection luminosity and thus a faster decrease of the density.More importantly though, larger opening angles imply a slower escape of particles from the emission region, and thus a slower decline of the density (see Eq. ( 3)).The peak in each curve is attained at roughly t ∼ t esc (0) (red vertical line), however with a slightly earlier peak time for larger opening angles because of the quicker injection density decrease.For α → 0, the lines would approach a constant as the blob approaches the steady-state.The solid lines in Fig. 1 are derived from simulations using OneHaLe, confirming the analytical result.
The dash-double-dotted line in Fig. 1 indicates n ∼ t −3 , which is the expected evolution of the density if it were solely due to the increase in volume of the blob.Most model lines are harder than this line implying a continuous increase of particle number in the emission region, even for p = 2.It is only for small opening angles (black and blue cases) for p = 2 that the model lines are softer and the total number of particles decreases at large times.This can be easily understood from Eq. ( 7), as the evolution at late times is governed by the power-law index: (2 + p) or 1/η esc α, whichever is smaller.If the injection switches off entirely, the density drops ∝t −1/η esc α , which for large opening angles is the same behavior as with continuous injection.This shows the aforementioned trapping of particles in the (rapidly) expanding blob.

Linear cascade evolution
In this section, we treat the additional injection of electrons through γ−γ pair production of a γ ray colliding with a soft photon.We do so by first deriving three simple (semi-)analytical scenarios.These are compared to simulations in Sec.2.3.4,where we discuss both the success and failure of the approach.We continue to neglect the energy dependency of the process n this work, as we are only interested in the rough time dependency of the total particle density.In other words, we assume that pair production takes place.This clearly is a very strong simplification, which may not hold in many cases.In this scenario, the injection rate of pairs can be written as: with the soft photon density, n soft (t), the γ-ray photon density, n γ (t), and the correlation factor, ξ ∼ σ T c, absorbing all constants and energy dependencies.
Both photon distributions can have various underlying production processes.The soft photons are most likely electronsynchrotron photons, or thermal photons from external sources, such as the AD, the BLR, and the DT.The γ rays can be produced from electron-inverse-Compton emission or protondependent processes, such as proton-synchrotron or neutralpion decay.The internal photon densities thus depend on the underlying particle distribution, implying n phot (t) ∼ n e/p (t), with the exception of SSC radiation, which depends quadratically on the electron distribution.As we only want to treat a linear cascade evolution -that is, Q γγ shall depend at most linearly on the electron distribution -we will ignore from now on the following combinations: electron-synchrotron and electroninverse-Compton, as well as external photon fields and electron-SSC.
In fact, the simplest case is hadronically induced γ rays absorbed by external photon fields, as this does not depend on the electron distribution at all.Cases with a linear dependency on the electron distribution are electron-synchrotron absorbing hadronically induced γ rays, as well as external photons absorbing electron-inverse-Compton radiation.This leaves us with three A151, page 3 of 15 distinct cases: The external photon fields comprising AD and isotropic sources, such as the BLR, are by themselves assumed to be timeindependent.However, the motion of the blob changes the distance to these external sources, implying: The first summand describes the AD field, which is roughly constant close to the accretion disk and falls off with distancesquared once the blob has traveled a distance that roughly corresponds to the AD radius, R AD (Dermer & Schlickeiser 2002).
The second summand represents isotropic photon fields (in the black hole frame), such as the BLR and the DT, within a given distance z i from the black hole.
The solution is thus separated into the primary injection following Eq.( 7) and the secondary injection leads to: While the isotropic contribution can be easily integrated, the accretion disk contribution leads to a hypergeometric integral without a simple analytical solution.
For isotropic external photons, that is, n AD = 0, the analytical result for Eq. ( 14) becomes The sum of Eqs. ( 7) and ( 15) is shown in Fig. 2. Compared to the standard solution, Eq. ( 7), there are significantly more particles injected within the boundaries of the external field, z i .The boundary of z i is noticeable by the break in the dash-dotted and solid lines.The peak of the particle density is attained later for larger opening angles compared to the standard solution.Beyond z i , the decay of the density again depends strongly on the opening angle.The standard solution is only reached for small opening angles.The semi-analytical solutions of Eq. ( 14) for the AD, that is n i = 0, are shown in Fig. 3.This behavior is comparable to the isotropic case except for an earlier peak owing to R AD z i .In both figures (Figs. 2 and 3) the horizontal dotted lines mark the maximum level of the standard solution.
the time delay for the semi-analytical solution to drop below this density.Naturally, this happens later for the isotropic field than for the AD field and depends strongly on the opening angle.In fact, the time delay can be orders of magnitude implying a significantly more pronounced and longer high state.
The linear cascade acts as a "catastrophic" injection with time scale t γγ (t) = [ξn p (t)] −1 .We directly obtain the generalization of Eqs. ( 5) and ( 7) as: As the protons follow Eq. ( 7), the integrals in the exponentials can be performed: With this solution, the remaining integral in Eq. ( 17) can be solved numerically.The result is shown in Fig. 4. The cascade develops only after t esc (0), while the opening angle dictates the further evolution.Small opening angles result in a high number density and a quick decay, while large opening angles seemingly keep the cascade going for very long times without any significant decay.This is also indicated by the dotted lines.

Case 3: External and electron-IC photons
This case is similar to the case 2, except that we have to replace n p (t) with n ext (t) and the γ rays stem from IC scattering of the external photon fields.Then the electron density becomes n e (t) = q 0,e 1 + η esc α t esc (0) t with the modified Heaviside function H For isotropic external photon fields, that is, n AD = 0, Eq. ( 19) can be analytically integrated: ) with the incomplete Gamma-function Γ(q, x).The result is shown in Fig. 5, where compared to Fig. 2, a strong pile-up is visible, which results from ξn i t > 1 close to z i .This may be due to an extreme choice of parameters, however, it does display the importance of the immediate feedback of the cascade on the γ rays in this case compared to case 1.We note that the pile-up for η esc α = 0.1 in Fig. 5 might be influenced by some approximation inaccuracies for small values of η esc α and so, a smaller pile-up is to be expected in this case.
photons before the γ rays have fully developed, resulting in a reduced cascade compared to the isotropic case.
In both cases, the dotted lines are roughly as long as in case 1.In the isotropic case, this is significantly influenced by the pile-up close to z i .With a slightly smaller size for the isotropic field and a resulting smaller pile-up for case 3, the duration of the high-state in case 3 would be shorter than in case 1.However, no such strong influence is expected in the AD case.

Simulations
We use the OneHaLe code (Zacharias 2021;Zacharias et al. 2022, see the appendix for a brief description) to model the timedependent evolution of the blob.In order to comply with some of the assumptions above, we have to slightly tweak the code.The γ−γ pair production injection rate is given by Aharonian et al. (1983), Cerruti et al. ( 2021) with the photon distribution, n ph , photon energies, and ˜ , normalized to the electron rest energy, and the electron Lorentz factor, γ.Equation ( 22) requires 1/˜ , ensuring, together with ˜ < 1, an absorber with low photon energies for the γ rays.
In order to suppress severe cooling (due to the extreme parameter settings; see Table 1), we switch off electron IC cooling.While this is a massive change, it is necessary to create sufficient synchrotron or IC photons.Furthermore, the magnetic field is kept constant throughout these simulations.Lastly, as we are only interested in the additional injection of pairs from γ−γ processes, the secondary injections from muon decay and Bethe-Heitler pair production have also been disabled.The input parameters for the simulations are given in Table 1, while the results are shown as solid lines in Figs.2-6.We emphasize that these parameters are not chosen for realism but merely to reproduce the semi-analytical models as closely as possible.
In this regard, the external fields labeled "BLR" and "DT" should just be considered as two separate external fields without any intended resemblance to actual BLR and DT photon fields.In case 1, the external fields act as absorbers.In case 3, they are needed for the IC process, while the absorption is done by electron-synchrotron emission.
Cases 1 and 3 can be reproduced reasonably well.In case 1, the isotropic photon field is made of two thermal fields with different temperature and luminosity, while in the corresponding case in case 3, only one thermal field is needed.In the latter case, the pile-up close to the edge of the external field is present, but not as strong as in the analytical model2 .The accretion disk field is represented by a Shakura-Sunyaev disk (Shakura & Sunyaev 1973), with luminosity L AD = η Edd L Edd .In case 3, the AD field absorbs IC emission scattering both isotropic fields (but not the AD field).
We have found no parameter set that could reproduce the semi-analytical result in case 2. The simulations with OneHaLe with the parameters given in Table 1 show a pronounced maximum shortly after t esc (0) followed by a decrease depending on the opening angle.This seems reasonable, as the dashed lines roughly correspond to the evolution of the γ rays (psynchrotron).In turn, the number of "absorbees" (i.e., γ rays) is quickly reduced.The long evolution suggested by the semianalytical model would only be achievable if the absorber were to increase in such a way that a higher fraction of the γ rays would be progressively absorbed.Then, at late times, practically all γ rays would be absorbed to keep the cascade going.This is implausible and marks the limit of the simple analytical model that ignores all the energy-dependencies of the cross-section.
The influence of the parameters on the simulations is shown in Fig. 7, where three variations for each case have been derived.
Notes.The first five parameters are the same in all simulations.The electron distribution parameters of case 1 correspond to the simulation used to reproduce the standard solution in Fig. 1.
It is clear that reproducing the analytical model requires a specific set of parameters.This also shows the limits of the approach, however, it is nonetheless reassuring that we have managed to reproduce two out of the three cases overall.

Interlude
The sections above have demonstrated the influence of the expansion of the moving blob on the development of the pair cascade.The effects are two-fold.First, a significantly higher density can be achieved through the pair production, as compared to no secondaries.Second, a larger opening angle can drag out the peak of the density evolution.While the densities for the standard solution peak around t esc (0), the curves with pair injection typically peak later than this timescale.The exception to this rule are secondary injections involving AD photons, because of the timescale, namely, R AD /Γ/c < t esc (0).For small values of z i , this could become the case as well for the isotropic photon field.For internal photon fields (case 2), however, the peak is always attained after t esc (0), as the internal fields only reach their maximum at this time and the developing cascade reinforces the interaction.Non-linear cascades, which we have not treated here, might increase this effect, as both absorber and absorbee are increased at later times.The dotted lines in Figs.2-6 mark the peak density of the standard solution and indicate the time at which the models with pair injection drop below this level.This typically occurs at least an order of magnitude after t esc (0).Of course, this also strongly depends on the amount of pairs produced, the parameters, and the opening angle -highlighting the importance of pair production and the geometry of the jet for the development and duration of a flux high state.
Naturally, the inclusion of the energy dependency on the particle distribution, as well as the pair production process, can have a significant influence.In the absence of efficient reacceleration of particles, the particle cooling will lead to a drop in the pair production process, as less energetic γ-ray photons will be produced.The adiabatic expansion of the blob will also cause a reduction in the energy density of the particles and the photon fields, similarly reducing the pair production.Hence, the time evolution might be shorter than envisaged in the present semianalytical model.

Modeling based on PKS 1510−089
The flat spectrum radio quasar (FSRQ) named PKS 1510−089, located at a redshift of 0.361, is one of the most famous blazars known to date.It was detected at very-high-energy γ rays with H.E.S.S. in 2009 (H.E.S.S. Collaboration 2013) and has since received significant attention (e.g., Ahnen et al. 2017;H.E.S.S. Collaboration 2021).Its modeling typically requires both the BLR and DT photon fields (Barnacka et al. 2014) or multiple zones (Nalewajko et al. 2012;Prince et al. 2019).
While no modeling of the source or a specific light curve is approached in this work, we use this source to study the impact of the pair cascade in a moving, expanding blob in terms of enhancing flaring events in a realistic setting.The modeling parameters are given in Table 2 and are based upon the works by Nalewajko et al. (2012) and Barnacka et al. (2014).We derived multiple modelings: (1) a steady-state model similar to the aforementioned works to provide a baseline flux that may describe PKS 1510−089 on average; (2) a moving blob using only leptonic processes similar to Saito et al. (2015); and (3) a moving blob including hadronic interactions.For the moving A151, page 7 of 15 ] esc =0.5 se=2.5, Le = 2.5 × 10 43 erg/s Edd=3.0LDT=0 erg/s Fig. 7. Small parameter study using OneHaLe.The black lines correspond in style to those in Figs.2-6, but for the case η esc α = 0.5.The solid colored lines show simulations varying the input parameter as labeled.In the bottom panels, the variation in the electron parameters (magenta lines) preserves the initial injection rate of ∼1 cm −3 s −1 .
Notes.The steady-state simulation uses a stationary emission region.Parameters below the horizontal line are used in each simulation.
blobs we again derive models for various opening angles as before.We also derive specifically models without γ−γ pair production to show the difference.For the moving-blob models 2 and 3, the launching position is chosen relatively close to the black hole, where the jet has not expanded much (Zacharias et al. 2022).At this position the blob is still small in order to fit into the jet, and its magnetic field is high compared to the steady-state one-zone model, which is located further down in the jet.One important addition to the simulations of the previous section is the consideration of magnetic flux conservation, that is, B(z) = B 0 R 0 /R(z), assuming a dominating toroidal guide field (Kaiser 2006).This also implies that the larger the opening angle, the faster the strength of the magnetic field drops.This may have a significant influence on the synchrotron emission, compared to the expectation from the simple model above.

Light curves
The multi-wavelength light curves are shown in Fig. 8.The wavebands are the very-high-energy γ-ray (VHE) band (E > 30 GeV 3 , corresponding to the energy band of the forthcoming CTA observatory), the high-energy γ-ray (HE) band (100 MeV < E < 100 GeV, corresponding to the energy band of the Fermi-LAT), the X-ray band (2−10 keV, corresponding to the range of Swift-XRT), and the optical R-band.It is immediately obvious that the flare evolution depends strongly on the wavelength and the presence of relativistic protons.
Models without pair production (solid lines) do not show major differences between the leptonic and hadronic simulations.Except for the VHE band, which is fully absorbed at early times, all light curves rise quickly to a maximum at the initial escape timescale, t esc (0), followed by a subsequent decay.The switch from IC-BLR cooling to IC-DT or synchrotron cooling during the crossing of the BLR edge4 results in a minor bump in the X-rays and the R-band, which is absent in the HE band.
The VHE light curve shows a minor bump at this point, as the environment is now optically thin to VHE photons.Except for the HE band, the decay of the light curve is faster for larger opening angles, which is a consequence of the diminishing magnetic field, and hence a faster decay of the synchrotron emission.We should note here that the X-ray emission from the blob is also influenced by electron-synchrotron emission, as the (initial) high magnetic field shifts the synchrotron peak to higher energies compared to the steady-state model.As the HE band is governed by inverse-Compton emission of external photon fields, the magnetic field has basically no effect, and the decay of the light curve is fully governed by the evolution of the electron density.In the hadronic model and for large opening angles, the VHE light curve shows an additional bump around the time of the DT-edge crossing.This radiation is electron-synchrotron emission from secondary particles injected from muon decay.
The inclusion of the pair production has almost no consequence for the VHE band, except for hadronic simulations of large opening angles, where the DT-edge-crossing bump is much more pronounced and shows a plateau phase lasting for a few hours.At that time, a peak also shows up in the X-ray domain for large opening angles.It is induced by pair production and is also visible (although less pronounced) in the leptonic simulations.The peak flux in the main flare is significantly higher compared to the case without pair production.Similarly, the HE band exhibits higher peak fluxes compared to the solid lines.In both bands, HE and X-rays, the peak fluxes are much more pronounced in the hadronic simulations indicating the significant impact of relativistic protons on the pair production.The R-band shows no change between pair or no pair production, implying that pairs are produced predominantly at higher energies.While cooling will reduce the energy of the secondary pairs, their density is insufficient to add to the synchrotron emission of the primaries in the R-band.
The most important effects in the light curves shown in Fig. 8 take place within a few days in the observer's frame.Figure 9 displays the first 100 days in the evolution of the radio band at 15 GHz.The rise of the light curve begins after a few days, when the flares in the other bands have already ceased (see also Boula & Mastichiadis 2022;Tramacere et al. 2022).For small opening angles, the radio flare diminishes within 50 days, while for large opening angles, the light curve keeps rising.While it is not shown, for η esc α = 0.9 the peak is reached only after more than 1000 days.There is barely a difference between the cases with and without pair production for the same reason as in the R-band.

Luminosity evolution
An important measure for any blazar model is the luminosity of the emission region.Power is contained in particles, magnetic field, and radiation.The evolution of the luminosities as a function of time is shown in Fig. 10.By design, the magnetic field luminosity is constant with a value of L obs B = 1.8 × 10 44 erg s −1 , which is why we do not display it in Fig. 10.The total luminosity in the bottom row of that figure is the sum of the electron, proton, magnetic field, and radiation luminosities.The individual luminosities are calculated from the respective energy densities in the comoving frame u i (t) according to L obs i (t) = u i (t)Γ 2 πR(t) 2 c deriving the bulk Lorentz factor Γ from the assumption Γ = δ.
For small opening angles (black and blue lines), the respective luminosities quickly rise to the maximum as the source is being filled with particles and radiation, and then decrease.In turn, the total luminosity reaches a bottom value that is given by the constant magnetic field luminosity.It also implies that the initial particle and radiative powers are much larger than the magnetic field power.For η esc α = 0.5, the particle powers reach a constant level implying that the injection and escape of particles is fully countered by the expansion.In the presence of protons (right column in Fig. 10), the electrons slowly continue to rise in the case where pair production is enabled (magenta dashed line).This implies continuous pair creation even at very late times.For large opening angles (orange and yellow lines), the particle luminosities keep rising, as the rapid expansion does not allow for a meaningful escape of particles.The radiative luminosity in this case continues to decrease initially, but reverses this trend for late time, as the continuing accumulation of particles also increases the amount of photons, which also remain longer and longer in the emission region.
Interestingly, the total luminosity is initially dominated by radiation and only after about 1 day in the observer's frame does it change to particle or magnetic field dominance.It implies that transfer of power into radiation is initially very efficient and that the particles are efficiently cooled.With the decrease in the magnetic field strength, and the decrease of density of internal and external photon fields due to expansion and motion, the radiation production becomes less efficient over time and its luminosity drops more rapidly than that of the particles.
A151, page 10 of 15 In all cases, the comparison between particles and magnetic fields suggests that the emission region is initially particle dominated.For large opening angles, this is maintained throughout, while for small opening angles, the magnetic field dominates the particles at late times.This is contrary to most models, where an initial magnetic field dominance is expected which changes to particle dominance later (e.g., Zacharias et al. 2022, and references therein).The reason is that we chose a relatively modest initial magnetic field of only 20 G along with a linear decrease with radius.If the magnetic field were initially mostly poloidal, one would expect a quadratic decrease with radius (Kaiser 2006), and a much faster decrease.In this case, the magnetic power would also drop as a function of time.While our initial value is in line with the measurement in M 87 (Event Horizon Telescope Collaboration 2021), theoretical works suggest much larger values along with a rapid decrease (Zdziarski et al. 2022).Testing such scenarios is left for a future work.
The black horizontal lines in Fig. 10 mark the Eddington luminosity and the AD luminosity, respectively.As the accretion process fuels the jet, the jet luminosity must be compared with these values.Apparently, all models remain initially close to or below the limit of the Eddington luminosity.In case of protons and pair production (right column, dashed lines), the Eddington luminosity is initially surpassed by a factor 5. This is, however, due to the large radiative power, which includes also the IC scattering of the external fields.While the BLR and the DT are less luminous than the AD, they are strongly beamed in the comoving frame.Thus, the short excess of the Eddington luminosity can be expected and is within reasonable bounds.The continued increase of jet power for large opening angles has already been discussed and naturally surpasses the Eddington limit by far.Clearly, the model is not realistic in these cases.

Conclusions
In this paper, we have discussed the particle evolution within a moving, expanding emission region within the jet of a blazar (Boula & Mastichiadis 2022;Tramacere et al. 2022).With a simple semi-analytical model, we have shown that the opening angle has an important influence on the particle distribution and the electron-positron cascade.With an increasing opening angle the high-state (including the cascade) may last much longer owing to the increased escape timescale for both the particles and photons.
With the help of the OneHaLe code (Zacharias 2021; Zacharias et al. 2022), we have shown that the simple model can be reproduced for linear cascades induced by external photon fields, while the linear cascade resulting from internal processes was not reproduced successfully.Naturally, due to the neglect of the energy dependencies of the particle evolution and the cascade, it required quite specific, and extreme parameters in the simulation.Nonetheless, the exercise has demonstrated the importance of the opening angle on the evolution of the particle distribution.
The simple model may be improved by considering deltafunction approximations to the energy dependency of the particle distribution and the cascade.This may provide a more reasonable view on the evolution providing a broader range of application.These calculations are left to the interested reader.
In order to obtain realistic simulations, we have made use of parameters from PKS 1510−089 (Nalewajko et al. 2012;Barnacka et al. 2014).The leptonic and hadro-leptonic simulations were conducted for various opening angles.From the resulting light curves in Fig. 8, we can derive the following conclusions.A moving blob induces a rapid flare.The evolution depends critically on the blob's speed and on the initial size at launch.Observations in the X-ray and VHE γ-ray are crucial to determine the opening angle and the particle composition, as for large opening angles these bands may (or may not) show a delayed flare when the blob moves out of the DT.These secondary flares are significantly influenced by the pair production process and are more pronounced if a hard and highly energetic proton distribution is present.Interestingly, the optical band is also useful to obtain information on the opening angle, as a small opening angle induces a longer decay phase due to the slower decay of the magnetic field.In the HE γ-ray band, the flare is bright and short with barely a difference between the various opening angles.Importantly, the HE band peaks slightly before the other bands (cf. the red vertical line in Fig. 8), whereas the X-ray and R-band peak at the time corresponding to t esc (0).In all bands, the (primary) flare is asymmetric with a fast rise and a slow decay (Saito et al. 2015;Boula & Mastichiadis 2022;Tramacere et al. 2022).
Truly simultaneous and high-cadence observations are paramount to obtain all the potential information.The main flare lasts just for a few hours, and the subtle differences between peak times and decay profiles might be hidden in the statistical errors or a delay in observation.It should also be noted that it is worthwhile to continue observing for a few days in the X-ray and VHE γ-ray band to search for the secondary flare.
Radio observations may also provide indications for the opening angle.Large values of η esc α result in years-long variations, while small opening angles result in flares lasting a few tens of days.Single-dish radio observations every few days for up to 100 days after the event at higher energies should provide sufficient information on the opening angle.The downside could be the higher chance of multiple events superposing each other.
The jet luminosities are mostly within reasonable bounds compared to the Eddington luminosity.However, for large opening angles, the particle luminosities keep increasing due to the lack of escape, which results in unreasonable jet powers.Additionally, the magnetic field evolution is such that its power remains constant.In turn, the initially partical-dominated emission region becomes magnetic-field-dominated at later times and small opening angles.This is contrary to standard expectations.A scenario, which is more in line with the usual jet evolution, will be discussed elsewhere.
The above mentioned evolutionary details of the flare depend on the type of object.The given results are based on parameters describing the FSRQ PKS 1510−089.In BL Lac objects with much weaker or even absent external fields, the effects induced from the pair cascade are probably much less pronounced (Zacharias et al. 2022).In this case, the solid lines in Fig. 8 might already be an adequate description and the secondary flares in the X-ray and VHE γ-ray band are absent.Similarly, the details of the magnetic field decay may have a significant influence on the evolution of the flux, especially in the R-band and needs to be discerned from the expansion profile.
Lastly, using the standard one-zone description implies that at all times an instantaneous particle spectrum was used.While shock acceleration can quickly accelerate particles (Böttcher & Baring 2019), it would be useful to properly include the particle acceleration in this model as it might have an important influence on the light curve evolution at early, but also at late times.At early times, the delayed injection of high-energetic particles would slow down the flare evolution, and might also reduce the strength of the cascade development.At late times, A151, page 12 of 15 the weakened magnetic field might not be able to accelerate particles to the highest energies, which might reduce or even inhibit the development of the secondary flare.It would also substantially reduce the particle luminosities.In order to adequately account for these effects, a two-zone model is needed (Weidinger & Spanier 2015;Chen et al. 2015;Dmytriiev et al. 2021), which is beyond the scope of this paper.
where we followed the prescription of Zdziarski et al. (2014).
As we assume Γ = const., the cooling due to external isotropic photon fields such as the BLR is constant in time, as long as the blob is within the distance of the external field: with the energy density of the external photon field in the black hole frame, u ext , and the radius of the external photon field, z ext .Again, we approximated β Γ ≈ 1.
Figure A.1 shows the evolution with time of the cooling and escape timescale.The adiabatic and acceleration (not shown) timescales behave exactly like the escape timescale, while the external-Compton cooling simply stops once the blob leaves the respective region.The synchrotron cooling timescale (with b = 1) increases rapidly for t > t esc (0).The SSC cooling timescale (including the time-dependency of Eq. ( 7)) initially decreases due to the increase in particle density.This reverses at t ∼ t esc (0), as the maximum density has been passed and particles overall leave the emission region.

Fig. 8 .Fig. 9 .
Fig.8.Model light curves based upon parameters describing PKS 1510−089 for VHE γ rays (top row), HE γ rays (2 nd row), X-rays (3 rd row), and the optical R-band (bottom row) for purely leptonic processes (left column) and including hadronic processes (right column).Displayed models are either with (thick dashed) or without (solid) γ−γ pair production.Colors are for various opening angles as labeled.The vertical lines mark t esc (0) (red), the crossing of the BLR edge (blue) and of the DT edge (magenta); all times in the observer's frame.For clarity, the launch of the blob takes place at t obs = 0.1 d.

Table 1 .
OneHaLe input parameters to simulated cases 1 to 3.