arxiv: v1 [astro-ph.he] 23 Jun 2016

Transcription

1 Mon. Not. R. Astron. Soc.,?? (6) Printed 7 June 6 (MN LATEX style ile v.) Blazar lares powered by plasmoids in relativistic reconnection Maria Petropoulou, Dimitrios Giannios and Lorenzo Sironi Department o Physics and Astronomy, Purdue University, 55 Northwestern Avenue, West Laayette, IN 4797, USA Department o Astronomy, Columbia University, 55 W th St, New York, NY 7, USA Received / Accepted arxiv: v [astro-ph.he] 3 Jun 6 INTRODUCTION Blazars are a small subclass o Active Galactic Nuclei (AGN), yet they attract an ever growing interest as they are ound in increasingly large numbers by surveys at microwaves and γ-ray energies (e.g. Giommi et al. 9; Abdo et al. ; Giommi & et al ; Ackermann et al. 5). Blazars also represent the most abundant population o extragalactic sources at TeV energies (e.g. Holder 4; de Naurois 5). The extreme observational properties o blazars, such as continuum emission over the entire electromagnetic spectrum, rapid and large-amplitude variability, make them to stand out among other AGN. The blazar broadband emission, rom radio up to very high energy (VHE) γ-rays (> GeV), is believed to originate rom a relativistic jet that is nearly aligned with the observer s line o sight and emerges rom the central supermassive black hole (Blandord & Rees 978; Urry & Padovani 995). The blazar spectral energy distribution (SED) is also very distinctive due to its double-hump appearance. A typical blazar SED is composed o two broad components: a low-energy and a high-energy Einstein Post Doctoral Fellow mpetropo@purdue.edu ABSTRACT Powerul lares rom blazars with short ( min) variability timescales are challenging or current models o blazar emission. Here, we present a physically motivated ab initio model or blazar lares based on the results o recent particle-in-cell (PIC) simulations o relativistic magnetic reconnection. PIC simulations demonstrate that quasi-spherical plasmoids illed with high-energy particles and magnetic ields are a sel-consistent by-product o the reconnection process. By coupling our PIC-based results (i.e., plasmoid growth, acceleration proile, particle and magnetic content) with a kinetic equation or the evolution o the electron distribution unction we demonstrate that relativistic reconnection in blazar jets can produce powerul lares whose temporal and spectral properties are consistent with the observations. In particular, our model predicts correlated synchrotron and synchrotron sel-compton lares o duration o several hours days powered by the largest and slowest moving plasmoids that orm in the reconnection layer. Smaller and aster plasmoids produce lares o sub-hour duration with higher peak luminosities than those powered by the largest plasmoids. Yet, the observed luence in both types o lares is similar. Multiple lares with a range o lux-doubling timescales (minutes to several hours) observed over a longer period o laring activity (days or longer) may be used as a probe o the reconnection layer s orientation and the jet s magnetization. Our model shows that blazar lares are naturally expected as a result o magnetic reconnection in a magnetically-dominated jet. Key words: acceleration o particles galaxies: active magnetic reconnection radiation mechanisms: non-thermal component extending, respectively, rom radio to UV/X-rays and rom X-rays to γ-rays. Blobs, or quasi-spherical emission regions containing relativistic particles and magnetic ields, have been oten invoked to explain the broadband variable emission o blazar jets (e.g. Bloom & Marscher 996; Mastichiadis & Kirk 997; Kirk, Rieger & Mastichiadis 998; Chiaberge & Ghisellini 999; Böttcher & Chiang ). Yet, their physical origin is still not understood. Here, we present a physically motivated model or the emitting blobs in a Poynting-lux dominated blazar jet (or details see Giannios, Uzdensky & Begelman 9; Giannios 3), in the context o relativistic magnetic reconnection; the latter reers to the regime where the magnetic energy per particle exceeds its rest mass energy, or the plasma magnetization σ exceeds unity 3. MHD instabilities o a Poynting-lux dominated low lead to the ormation o current There are also certain models that aim at explaining the low-energy (Marscher & Gear 985; Marscher & Travis 996) or/and the high-energy (e.g. Maraschi, Ghisellini & Celotti 99; Reynoso, Medina & Romero ; Potter & Cotter ; Reynoso, Romero & Medina ) blazar emission in terms o a jet model where the particle population dynamically evolves along the jet rom its base to pc scales. 3 The plasma magnetization is deined as σ = B /4πρ c, where B and 6 RAS

2 Petropoulou, Giannios, & Sironi sheets where magnetic reconnection is triggered (e.g. Eichler 993; Begelman 998; Giannios & Spruit 6). Magnetic reconnection is an inherently time-dependent, highly dynamic process as solar observations and recent numerical simulations have revealed (e.g. Lin et al. 5; Karlický & Kliem ). For the highly conducting plasma o blazar jets, the reconnection current sheets are susceptible to tearing instabilities that lead to their ragmentation in a chain o plasmoids (or, magnetic islands), i.e., regions containing magnetic ields and energetic particles (e.g. Loureiro, Schekochihin & Cowley 7; Daughton & Karimabadi 7; Bhattacharjee et al. 9; Loureiro et al. ). The plasmoids grow rapidly through mergers beore leaving the reconnection region. Occasionally, plasmoids can undergo signiicant growth to a sizable raction o the reconnection region, orming monster plasmoids (Uzdensky, Loureiro & Schekochihin ). The time-dependent aspects o reconnection may prove to be crucial in understanding blazar lares (in terms o energetics and timescales) as shown by Giannios, Uzdensky & Begelman (9) and Giannios (3). The most undamental way to capture the ormation, dynamics, particle and magnetic energy content o plasmoids in reconnection layers is by means o ully kinetic particle-in-cell (PIC) simulations. PIC simulations o reconnection have been recently extended to the relativistic regime o σ that is relevant to blazars (Zenitani & Hoshino ; Guo et al. 4; Sironi & Spitkovsky 4; Nalewajko et al. 5; Sironi, Petropoulou & Giannios 5; Kagan, Nakar & Piran 6; Werner et al. 6). In particular, Sironi, Petropoulou & Giannios (5) henceorth SPG5, showed using two-dimensional (D) PIC simulations 4 o electron-positron and electron-ion plasma that relativistic reconnection can satisy all the basic conditions or the blazar emission: eicient dissipation, extended particle distributions, and rough equipartition between particles and magnetic ield in the emitting region, thus supporting the view that the dissipated energy appears in bursts associated with individual plasmoids in the layer (Giannios 3). The statistical properties o the plasmoid chain, such as size distribution, particle and magnetic energy content, were recently presented in Sironi, Giannios & Petropoulou (6) (hereater, SGP6). By employing D PIC simulations in pair plasmas extended to unprecedentedly long time and length scales, SGP6 were able to assess the basic properties o the reconnection layer, namely the particle distributions, the geometry, and the motion o individual plasmoids, as a unction o the system size. This allows us to extrapolate these results rom the small plasma scales o PIC simulations to the macroscopic scales relevant or the blazar emission. Aim o the present study is to incorporate the physics that describes the plasmoid, e.g. growth rate o a plasmoid, magnetic ield strength, and injection rate o particles, into a model or its emission. To achieve our goal, we combine (i) recent results rom PIC simulations as presented in SGP6 with (ii) the kinetic equation or the evolution o the distribution o radiating particles and their synchrotron, synchrotron sel-compton (SSC) emission. Our approach provides physical insight on the basic properties o the plasmoid-powered lares, such as their rise timescale, and leads to several robust predictions. We show that synchrotron and SSC lares o duration o several hours days are powered by the largest and slowest plasmoids that orm in a reconnection layer. Smaller ρ are the magnetic ield and mass density o the plasma outside the reconnection layer. These are measured in the rest rame o the jet luid. 4 Sironi & Spitkovsky (4) showed that the long term evolution o relativistic reconnection, including particle acceleration, proceeds similarly in D and 3D. and aster plasmoids, on the other hand, produce lares o similar, or even higher, peak luminosity but with sub-hour duration. This paper is structured as ollows. In we summarize the basic results o the PIC simulations presented in SGP6. In 3 we present an analytical model or the electron distribution in the blob and its motion. In 4 we present the basic properties o lares produced by individual plasmoids and we apply our model to the blazar emission in 5. In 5. we present analytical expressions or several observables o lares and we continue with some indicative examples o plasmoid-powered lares in 5.. We discuss various aspects o our model in 6 and we conclude in 7 with a summary o our results. BASIC RESULTS FROM PIC SIMULATIONS. Summary o SGP6 results In SGP6 we employed a suite o large-scale D PIC simulations in electron-positron plasmas to demonstrate that relativistic magnetic reconnection can naturally account or the ormation o quasispherical plasmoids illed with high-energy particles and magnetic ields. The simulations extended to unprecedentedly long temporal and spatial scales, thus allowing us to capture the asymptotic physics o plasmoid ormation independently o the initial setup. We showed that the plasmoids are continuously generated as a selconsistent by-product o the reconnection process and their most important properties are: (i) there is rough equipartition between particle kinetic energy density and magnetic ield energy density; (ii) the comoving particle density and magnetic ield strength o each plasmoid remain approximately constant during its growth; (iii) the plasmoids grow in size at. o the speed o light (i.e., at about hal o the reconnection inlow rate), with most o the growth happening while they are still non-relativistic; (iv) their growth is suppressed once they get accelerated to relativistic speeds, up to a terminal our-velocity σ c; (v) the width o the largest (monster) plasmoids is.l, independently o the size o the reconnection layer l ; (vi) plasmoids with sizes much larger than the characteristic plasma scales 5 contain isotropic particle distributions (see also Table, Figs. 7 and in SGP6); and (vii) the typical recurrence interval o the largest plasmoids is.5l /c, while smaller plasmoids are more requent. 3 FLARES FROM PLASMOIDS: AN ANALYTICAL MODEL In this section we incorporate the results rom PIC simulations regarding the dynamical evolution o the plasmoids into an analytical model or the plasmoid emission. We irst present the basic assumptions that enter our calculations ( 3.) and continue with a description o the plasmoid and particle evolution ( ). The characteristic properties o plasmoid-powered lares are presented in 4. 5 The characteristic plasma scale is the plasma skin depth deined as c/ω p, where ω p = eb /mc σ is the plasma requency. For the astrophysical application o this paper B G and σ, which results in skindepths o 6 cm l 6 cm; see next sections.

3 Blazar lares rom plasmoids 3 3. Assumptions (i) The plasmoid is homogeneous with constant particle density and magnetic ield strength throughout its volume; PIC simulations show that plasmoids have a structure, namely the particle kinetic and magnetic energy densities peak at the core o the plasmoid and decrease towards its outer parts (see e.g. Fig. A in SGP6). Although an inhomogeneous emission model is more realistic, this will not signiicantly alter the main conclusions o this work regarding the multi-wavelength spectra and light curves, since we use the volume-averaged properties o the plasmoid as determined by the PIC simulations. However, the inhomogeneous structure o the plasmoids is crucial or calculating the polarization signatures and we plan to investigate this in the uture. (ii) The plasmoid in its rest rame is a sphere; although this is a good approximation or plasmoids studied in D simulations (see e.g. Fig. 5 in SGP6), the plasmoids seen in 3D simulations o relativistic reconnection are best described as ellipsoids elongated along the direction o the electric current. Since 3D PIC studies o the plasmoid chain ormation are still premature, we will adopt the results o D simulations. (iii) The particle distribution contained in a plasmoid is isotropic or plasmoids o all sizes; PIC simulations show that anisotropy is present in the smallest plasmoids with sizes a ew tens o the plasma scale, whereas particles conined in the largest plasmoids have approximately isotropic distributions. Given that macroscopic plasmoids responsible or the blazar lares are much larger than the plasma scales o the problem (Giannios 3), it is sae to assume that they are characterized by quasi-isotropic particle distributions. Regardless, even i the anisotropy is present in the early phases o growth, when the plasmoid is still small, we show that the emission produced at this stage is a negligible raction o the emission at the peak time o a blazar lare. Thus, our assumption would not introduce any substantial errors in our estimates o peak luminosity and lare timescales. (iv) The particles are accelerated to an extended power-law distribution; the minimum and maximum Lorentz actors o the electron distribution are estimated or an electron-proton plasma. Although PIC simulations in SGP6 were perormed or an electronpositron plasma, we argue that all their basic results will hold or electron-proton reconnection as well, since or σ the ield dissipation results in nearly equal amounts o energy transerred to protons and electrons (see SPG5). So, the mean energy per particle o the two species is nearly the same, as it is the case or an electronpositron plasma (or more detailed discussion, see Sect. 3.3). 3. Plasmoid motion and Doppler boosting We assume that at a distance z diss rom the base o the jet (see Fig. ) a current sheet o length l is ormed (see 5 or the large-scale jet model and the motivation or l, z diss parameters). This is embedded in a relativistic low with bulk Lorentz actor Γ j. Let us consider a plasmoid that orms close to the central X-point o the current sheet, i.e. at X l, with initial width.6 This moves along the current sheet with a speed β co (in units o the speed o light) as measured in the jet s rest rame and Lorentz actor Γ co ( /. βco) 6 There are three reerence rames that are o relevance in our study: (i) the rest rame o the plasmoid (double-primed quantities), (ii) the rest rame o the jet (primed quantities), and (iii) the observer s rame (unprimed quantities). z diss Phase I Reconnection layer l' w '' X ' Z axis l' w '' X'=l' θ obs θ j Phase II Figure. Sketch o a reconnection layer ormed at a distance z diss in the jet. The angle between the observer s line o sight and the jet axis is θ obs while θ j is the jet s opening angle. A plasmoid ormed close to the central X-point o the current sheet grows in size as it moves along the current sheet (Phase I) and accumulates particles. At the same time it may accelerate rom nonrelativistic to relativistic speeds, reaching a terminal velocity that is close to the Alvén speed. The injection o particles ceases when the plasmoid leaves the current sheet (Phase II). The growth o the plasmoid size in this phase is caused by expansion in, e.g., the under-pressured surrounding jet plasma. Based on the results presented in SGP6 (in particular, see Fig. and eq. () therein), the plasmoid s momentum as measured in the jet rame is related to X / as ( X β co Γ co ) σ tanh ( βacc σ X X w'' ), or, () where X is the position along the current sheet and β acc is a dimensionless number determined numerically that quantiies the acceleration rate o the plasmoid and is approximately independent o the magnetization σ (i.e., β acc =..5 or σ = 3 5). Exactly at X the plasmoids are ormed with initial momentum β coγ co that is not included in eq. () or simplicity. Based on eq. (), two asymptotic regimes o the plasmoid s motion can be identiied: β co Γ co β acc (X / ), or X / 3( σ /β acc, ) β co Γ co 3 σ, or X / 3( σ /β acc, ), where σ = σ and β acc = β acc,. As the plasmoid moves along the current layer it grows in size, mainly through mergers, with a rate β g that is a signiicant raction o the speed o light. SGP6 determined that β g.6,.8, and. or σ = 3,, and 5 respectively. These values o the growth rate are appropriate or plasmoids whose speed is not too close to the Alvén speed β A c = σ/( + σ)c. To account or the act that the plasmoid growth gets slower as β co β A (see Fig. 8 in SGP6) we replace β g by a suppressed growth rate. This is modelled by β g [ + tanh(β co /β A )], so that the asymptotic growth

4 4 Petropoulou, Giannios, & Sironi σ= σ= 8 7 σ= β co Γ co.5 β co Γ co.5 β co Γ co log (w /l ) log (w /l ) log (w /l ) Figure. Plasmoid momentum β co Γ co as a unction o its width normalized to the hal-length o the layer, /l, or dierent magnetizations marked on the plot. Each coloured line denotes the evolutionary path o a single plasmoid orming at X =.l with initial width (black dashed line, nearly horizontal at Γ co β co ) and exiting the reconnection layer (X = l ) with (black dashed-dotted line). The hatched region corresponds to plasmoids with >.3l that are very rare. The paths o plasmoids that leave the layer with β co Γ co > are shown with blue coloured lines. Cyan coloured lines correspond to plasmoids that exit the layer while being non-relativistic. The red coloured dashed line marks the asymptotic value ( σ) o the plasmoid dimensionless our-velocity. Other parameters used are: β g =.6, β acc =. (let panel), β g =.8, β acc =. (middle panel) and β g =., β acc =.5 (right panel), as inormed by the SGP6 PIC simulations. rate is β g /3. 7 Thus, the equation that governs the plasmoid growth in its rest rame is d = β gcdt g ( X ) = β g cdt Γ co ( X ) g ( X ), () where g(x) + tanh (β co (x)/β A ). In the numerical calculations that ollow the modiied β g will be used, whereas β g will be kept constant in the analytical calculations. This introduces only a small error while greatly simpliies the analytical calculations. Solving or β co = dx /c dt rom eq. () we ind ( ) X dx = cβ co dt = cdt w + ( ). (3) X The equation that relates the size o the plasmoid with its position X along the current sheet is ound by combining eqs. () and (3) ( ) ( ) dx X = β X d g g. (4) The rate o momentum change at the initial stages o a plasmoid s evolution can be written as d(β co Γ co )/dt = c(β acc β g )Γ co β co /, where eqs. () and (4) have been used and g. Thus, the plasmoid will accelerate, or β acc > β g as ound in PIC simulations. The momentum β co Γ co as a unction o /l is plotted in Fig. or iducial plasmoids and dierent magnetizations marked on the plots. Each coloured line denotes the evolutionary path o a single plasmoid orming at X =.l with initial width (black dashed line) and exiting the reconnection layer with (black dashed-dotted line). The hatched region corresponds to extremely rare, big plasmoids with >.3l. Large plasmoids (.l ) leave the reconnection layer with a mildly relativistic speed β co Γ co, whereas smaller plasmoids (.4l ) exit with the asymptotic (dimensionless) our-velocity σ (see also Lyubarsky (5)). As case studies o a large/slow plasmoid and a small/ast plasmoid we, respectively, adopt =.l and =.4l. 7 We have checked that a dierent orm o the suppression actor does not change our main conclusions. We urther assume that the current sheet and, in turn, the plasmoid s direction o motion orm an angle θ with respect to the jet axis, which or convenience we deine as the z-axis. Without loss o generality, we adopt the y z plane to be the plasmoid s plane o motion. The plasmoid s velocity in the observer s rame is deined as β p and its components parallel and normal to the jet axis are given by β p,z = β j + β co cos θ + β j β co cos θ, β β co sin θ p,y = Γ j ( + β j β co cos θ ) The corresponding Lorentz actor Γ p = ( β p) / is given by Γ p = Γ j Γ co ( + βj β co cos θ ), (6) while the angle o the plasmoid velocity with respect to the z axis as measured in the lab rame is β co sin θ tan θ = (. (7) Γ j βj + β co cos θ ) Let the observer s position also lie on the y z plane and θ obs be the angle between the observer s line o sight (see Fig. ) and jet axis and ω = θ θ obs be the angle between the plasmoid s direction o motion and the line o sight. The Doppler actor o the plasmoid is then deined as δ p = ( Γ p βp cos ω ), (8) where cos ω = cos θ cos θ obs + sin θ sin θ obs. The Doppler actor corresponding to the bulk motion o the jet alone is deined as δ j = Γ j ( βj cos θ obs ). (9) In the extreme case o perect alignment among the direction o the plasmoid s motion, the jet s bulk motion and the line o sight, namely θ = θ obs =, one derives the maximum Doppler boosting δ p 4Γ co Γ j (see also Giannios, Uzdensky & Begelman 9) or relativistic plasmoids. For non-relativistic plasmoids and perect alignment the Doppler actor simpliies to δ p δ j Γ j. The evolution o δ p as a plasmoid grows in size and accelerates is illustrated in Fig. 3 or σ =, two values o θ marked on the plot and a viewing angle θ obs =.5/Γ j. Blue and cyan (5)

5 Blazar lares rom plasmoids 5 log δ p Γ j σ.6γ j σ Γ j 4Γ j / σ log (w /l ) θ = θ =π/4 Figure 3. Log-log plot o the plasmoid Doppler actor δ p as a unction o its width /l or σ =, and two choices o the inclination o the layer: θ = (dashed lines) and θ = π/4 (solid lines). Here, the viewing angle is θ obs =.5/Γ j. Blue and cyan coloured lines have the same meaning as in Fig.. The horizontal grey dashed lines indicate the limiting values or δ p in the case o perect alignment (i.e., θ = θ obs = ) assuming a nonrelativistic plasmoid (bottom dashed line) or a relativistic plasmoid moving with the maximum speed (top dashed line). The solid red lines mark the limiting Doppler actor o a relativistic plasmoid or a non-avorable orientation (i.e., θ obs =.5Γ j, θ = ; bottom line) and a avorable orientation (i.e., θ obs =.5/Γ j, θ = π/4; top line). coloured lines have the same meaning as in Fig.. The horizontal grey dashed lines indicate the limiting values or δ p in the case o perect alignment (i.e., θ = θ obs = ) assuming a non-relativistic (bottom dashed line) or a relativistic (top dashed line) plasmoid moving with the Alvén speed (β co Γ co Γ co σ). Fig. 3 shows that the Doppler actor o a plasmoid varies signiicantly during its growth phase and that it depends strongly on the combination o θ obs and θ. The curves δ p ( /l ) or the two choices o the layer s orientation considered here are not only quantitatively but also qualitatively dierent. A ew remarks on the plot ollow. For a non-avorable orientation, i.e. θ obs =.5/Γ j and θ =, the plasmoid Doppler actor is written as ( ) Γco ( + β j β co ) δ p = Γ j Γ co ( + β j β co ) +. () There are two regimes o interest in the evolution o a single plasmoid, namely (i) the non-relativistic regime, where β co Γ co. Equation () simpliies then into δ p (8/5)Γ j 6Γ j,, in agreement with with the asymptotic behavior at small sizes in Fig. (3). In act, (8/5)Γ j is the initial Doppler actor o all plasmoids, since they orm as non-relativistic structures. (ii) the asymptotic and relativistic regime, where β co Γ co Γ co σ. In this case, the Doppler actor is written as δp 4Γ j / σ Γ j, / σ. This value is indicated in Fig. 3 by the lower red solid line. For a avorable orientation, i.e. θ obs =.5/Γ j and θ = π/4, the Doppler actor in the non-relativistic regime is the same as or θ = (see point (i) above), whereas in the relativistic regime it is given by δ p Γ p, since ω = θ θ obs.. Thus, δ p Γ j Γ co ( + cos θ )/[ + (Γ j Γ co ω( + cos θ ) ], which results in δ p 3.4 σγ j /[+.3Γ j, σ ] or plasmoids moving with the asymptotic our-velocity (top red solid line in Fig. 3). Plasmoids that leave the layer with inal sizes.3l are accelerated close to their asymptotic our-velocity (see also Fig. ) and their Doppler actors depend sensitively on the orientation o the layer and the observer (compare blue solid and dashed lines). The dependence o δ p on the plasmoid size is determined by both the acceleration and growth processes, while it aects the relative timescales as measured in the observer s rame (see eq. () below and 4.). On the contrary, or a ixed viewing angle the orientation o the layer becomes irrelevant or non-relativistic plasmoids and δ p Γ j (cyan solid and dashed lines). Fig. 4 shows the ratio δ p ( )/δ j or dierent orientations o the observer (θ obs ) and the plasmoid s direction o motion (θ ) with respect to the jet axis. The results or σ = and two dierent plasmoid sizes are shown in the let ( =.4l ) and right ( =.l ) panels. For a small and ast moving plasmoid (let), the plasmoid s Doppler actor may be up to 5 times larger than that o the jet s bulk motion alone. For a bigger, yet slowly moving, plasmoid (right), the ratio δ p ( )/δ j varies at most by a actor o.4 and it is less angle-dependent. The time or a plasmoid to grow to its size as measured in the observer s rame is given by t( + z) = dt δ p ( ) = g( w) d w cβ g δ p ( w), () where X ( ) is determined by eq. (4) and g(x) is deined below eq. (). For =, the above expression determines the typical peak time or lares powered by the plasmoid. As we detail below (see next section), the calculation o the electron evolution and emission are perormed in the rest rame o the plasmoid. I L, L (ν ) are, respectively, the bolometric and dierential photon luminosities in the plasmoid s rame, their values in the observer s rame are L = δ 4 pl and L(ν) = δp 3+α L (ν ), where α is the spectral index. Here, we neglect inverse Compton scattering on external (to the jet) radiation ields which would result in a boosting actor o δ 4+α p (Dermer 995). 3.3 Plasmoid evolution The plasmoid evolution can be divided into the ollowing phases (see also Fig. ): Phase I that corresponds to the time that the plasmoid spends in the reconnection layer while growing in size, and Phase II that begins with the plasmoid leaving the reconnection layer and undergoing expansion in the jet s bulk low. The modeling o Phase I is entirely based on the results o PIC simulations that we summarize below: (i) the volume o the plasmoid increases as V 3, (ii) the magnetic ield strength in the rest rame o the plasmoid, B p, is approximately independent o its size, (iii) the electron number density n remains approximately constant with respect to the plasmoid s size (see e.g. Fig. 5 in SGP6), (iv) the total number o electrons N increases as N = n V 3, (v) in Phase I the injected particles cool because o radiative energy losses (i.e., synchrotron and inverse Compton scattering); adiabatic cooling is not relevant in this phase, since the number density remains constant, and (vi) the injected particle distribution is a power law rom γ min to γ max and slope p as determined by PIC simulations.

6 6 Petropoulou, Giannios, & Sironi Figure 4. Two-dimensional maps o the ratio o δp (w )/δj created or σ = and dierent values o θ (i.e., the angle between the jet axis and the plasmoid s direction o motion in the jet rame) and θobs (i.e., the angle between the jet axis and the observer s line o sight). Overplotted with black lines are contours o ixed ratios δp (w )/δj. For illustration purposes, the colour scale in the two plots is dierent. Let and right panels correspond to plasmoids with w =.4`, respectively. The smaller plasmoid is characterized by relativistic motion in the jet rame and its emission is strongly beamed or avorable and w =.` inclinations o the observer. Steep particle spectra with p > are ound or σ., whereas p. or σ & (Sironi & Spitkovsky 4; Sironi, Petropoulou & Giannios 5; Guo et al. 4; Werner et al. 6). For electron-proton plasmas with p > and pair multiplicity N±, the minimum Lorentz actor is given by γmin ' rec σ p mp, N± p me () where rec is the raction o dissipated magnetic energy in relativistic reconnection that is transerred to electrons (see eq. (3) in SPG5). In the derivation above, we assumed that N± mp /me and we also made use o the act that the mean energy per particle available or dissipation is σ/, where σ/ is implicitly assumed. Sironi, Petropoulou & Giannios (5) showed that rec ' 5% or σ = 3, but reaches the asymptotic value o 5% or σ > (see Figs. 3 and 4, therein). These values or rec were obtained or electron-ion reconnection with N±. I the jet is dominated by pairs in terms o number (N± ), these are expected to pick up most o the dissipated magnetic energy at the expense o the ew protons that are present in the jet. Then, the typical value o rec will be.5, as appropriate or electron-positron plasmas; nearly hal o the available energy will still remain in the magnetic ield. For σ. the exact value o γmax does not aect the results due to the steepness o the power-law spectrum (p > ). On the contrary, or σ the energy is carried by the highest energy particles. Thus, the maximum Lorentz actor γmax is limited by σ (Sironi & Spitkovsky 4; Guo et al. 5; Werner et al. 6) and is given by #/( p) " rec σ ( p)mp, (3) γmax ' N± (p )me where γmin '. These are the gross characteristics o particle distributions in plasmoid-driven reconnection. 3.4 Electron distribution All the calculations that ollow are perormed in the rest rame o the plasmoid. The evolution o the electron distribution, N(γ, t ), can be described by a partial dierential equation, which has the general orm! N dγ + N = Q(γ, t ), (4) t γ dt where dγ /dt < is the total energy loss rate and Q(γ, t ) is the injection rate o particles with Lorentz actors between γ, γ +dγ. The term that describes the escape o particles rom the plasmoid (i.e., N/tesc ) is omitted rom eq. (4), since particles with γ 6 γmax are conined in the blob (see also SGP6, or details). Phase I We express irst eq. (4) in terms o the plasmoid width w using the relation dw = βg cdt8 NI ks,i NI γ = QI (γ, w ), w γ (5) where the subscript I is used to remind o Phase I and ks,i = σt B p 6πme c βg. (6) Inverse Compton cooling on a ixed photon target ield could be easily incorporated in the above expression as long as the scatterings take place in the Thomson regime. We also remark that our 8 The growth speed is assumed to be constant in all analytical calculations, i.e, without the suppression actor g(x /w ).

7 Blazar lares rom plasmoids 7 analytical calculations presented in this section and Appendix A do not apply to cases where SSC cooling dominates the electron energy losses (see e.g. Schlickeiser 9). The injection rate or < may be written as Q I = π w n p γ p H[ ]H[γ γ min ]H[γ max γ ], (7) where p = (p )/(γ p+ min γ max p+ ) and the normalization o Q I is determined by the condition dn I /d = (π/)n. The solution to eq. (5), N I (γ, ), is derived in Appendix A. cooling break and the lower/upper cutos o the distribution can be used directly to calculate the synchrotron and SSC emission. The total number o electrons in Phases I and II can be ound by integrating the expressions (A9) and (A33), respectively, over γ. As expected, the total number o electrons increases as 3 during Phase I and remains constant ater that. This, in combination with the constant magnetic ield in Phase I and the decaying magnetic ield in Phase II, implies that the peak o the emission is expected at the end o Phase I. Phase II In contrast to Phase I, this evolutionary stage is not benchmarked with PIC simulations. In the analytical treatment (see Appendix A) we assume that the injection o particles ceases abruptly when the plasmoid leaves the layer, i.e. X = l. A gradual cessation o particles is used though in the numerical examples presented in 5.. By this time, the plasmoid has grown up to a size that can be determined rom eq. (4) given an initial condition X ( ) = X. In Phase II particles undergo synchrotron and adiabatic losses. Because o the expansion o the plasmoid the strength o the magnetic ield is expected to decrease, i.e. B q. The exact value o the exponent q depends upon the magnetic ield topology in the plasmoid. For example, i B is turbulent it can be treated as a relativistic luid with adiabatic index 4/3, leading to q =. In addition, the plasmoid s expansion is expected to depend on the properties o the jet s bulk low, such as pressure and density proiles. In order to keep our ormalism as general as possible we model the expansion rate and magnetic ield strength as power laws o the size : ( ) w a β exp = β exp,, (8) B = B p ( w ) q, q >. (9) where β exp is the expansion velocity and a is a ree parameter that could be or a constant expansion rate or a or an accelerating/decelerating expansion o the plasmoid. The kinetic equation that describes the evolution o electrons during Phase II is written as N II w + γ where (N II dγ d ) = Q II (γ, ), () dγ = γ d w k γ s,ii, () w q+a with the coeicient k s,ii deined as k s,ii = σ TB p q+a. () 6πm e c β exp, The source term in eq. () may be written as Q II (γ, ) = N I (γ, )δ ( ), (3) namely it describes an instantaneous event o particle injection when the source has width, carrying the evolution history o the distribution during Phase I. The electron distribution in Phase II is presented in Appendix A. As long as SSC cooling does not dominate the electron energy losses (see e.g. Schlickeiser 9; Zacharias & Schlickeiser ), eqs. (A9) and (A33) with the accompanying expressions or the 4 PROPERTIES OF PLASMOID-POWERED FLARES 4. Peak lare luminosity The observed peak luminosity o a lare will depend on: (i) the plasmoid s size,, when it leaves the reconnection layer, (ii) the plasmoid s Doppler actor at the end o Phase I, δ p ( ), (iii) the total number o radiating electrons, i.e., N n 3, (iv) the timescale over which the plasmoid grows (dynamical) t dyn /cβ g, and (v) the cooling timescale t cool o radiating particles. Assuming that the radiating particles are ast cooling (t cool t dyn ), which is reasonable or UV/X-ray emitting electrons, the observed peak luminosity o a lare at a given requency depends on (i)-(iv) and scales as L pk (ν) [ δ p ( )] 3+α w, (4) with α = p/. The inal size and momentum o a plasmoid carries inormation o its prior acceleration in the current sheet (see Fig. ). Although larger plasmoids contain more particles, smaller plasmoids may power as luminous, or even more powerul, lares as their emission is more strongly Doppler boosted towards the observer. We argue that the brightest lares are not necessarily powered by the biggest plasmoids. The above mentioned arguments are exempliied in Fig. 5 where the measure o the peak luminosity, [δ p ( )]3+α ( /l ), is plotted against /l or dierent σ (let to right). For each σ, the results or dierent inal plasmoid sizes are shown. For all σ values, the peak lare luminosity is not a monotonic unction o the size or those plasmoids whose radiation is beamed towards the direction o the observer (here, or θ = π/4 and π/7) and its dependence on /l can be understood as ollows. For avorable orientations between the layer and the observer the plasmoid s Doppler actor is δ p Γ p = Γ j Γ co ( + β j β co cos θ ). Plasmoids leave the layer with non-relativistic speeds i >,c, where,c. β acc, l. (5) In the above, we used eq. (), the condition β co Γ co = and the approximation tanh(x) x. 9 The Doppler actor o plasmoids leaving the layer with non-relativistic speeds is Γ j and thereore independent o the plasmoid s size. Thus, the peak lare luminosity depends quadratically on /l or large plasmoids. This is exactly illustrated in Fig. 5 (see blue and orange coloured lines or >.l ). Similarly, δ p ( ) is independent o w or plasmoids that 9 This approximation, albeit useul, underestimates the actual,c by a actor o two. We numerically ound that β co Γ co or plasmoids with given by eq. (5). Equivalently, the corresponding Doppler actor is higher than Γ j and, in particular, δ p, 4Γ j instead o Γ j. Regardless, the Doppler actor o plasmoids leaving the layer with non-relativistic speeds is independent o the plasmoid s size.

8 8 Petropoulou, Giannios, & Sironi 6 σ=3 6 σ= 6 σ=5 log[δ p 3+α () (w /l ) ] θ = θ =π/6 θ =π/4 θ =π/7 θ =π/ log [w /l ] log[δ p 3+α () (w /l ) ] θ = θ =π/6 θ =π/4 θ =π/7 θ =π/ log [w /l ] log[δ p 3+α () (w /l ) ] θ = θ =π/6 θ =π/4 θ =π/7 θ =π/ log [w /l ] Figure 5. Log-log plot o the quantity [ δ p ( )] 3+α ( /l ), which is a measure o the peak luminosity o a plasmoid-powered lare, as a unction o /l or σ = 3,, 5 and various values o θ marked on the plot. A ixed viewing angle θ obs =.5/Γ j has been adopted. The hatched region corresponds to >.3l. Here, α = p/ is adopted with p = 3, and.5 or σ = 3, and 5, respectively (Sironi & Spitkovsky 4). All other parameters are the same as in Fig.. For avorable orientations o the current sheet small plasmoids power bright lares (see dashed orange and solid blue lines). Figure 6. Two-dimensional maps o the quantity [ δ p ( )] 3+α ( /l ) or a plasmoid with =.4l (let) and.l (right). Here, σ =. Each map is created or dierent viewing angles θ obs (in the lab rame) and directions o plasmoid s motion with respect to the jet axis θ (in the jet s rest rame). The colour scale is also dierent in the two plots. All other parameters are the same as in Fig.. leave the layer with β co Γ co σ. The characteristic size can be estimated as, σ β accl σ =.3 β acc, l σ. (6) Plasmoids with inal sizes (β g /3β acc ), < σ w, have σ been accelerated to their asymptotic momentum even beore exiting the layer. The respective Doppler actor is then δ p 4Γ j σ or, σ and the peak lare luminosity scales as w (see the plasmoid s Lorentz actor also Fig. 5). For, σ < w <,c is Γ co < σ and scales as Γ co / ; here, eq. () and the linear approximation o tanh(x) were used. For avorable orientations δ p Γ co and the peak luminosity measure scales as α. The scaling is more accurate or higher σ where where the range o Lorentz actors between and σ is wider (see right panel in Fig. 5). The scalings o the peak luminosity measure or,c and, apply also to the non-avorable orientations, e.g. σ θ = (black solid lines). However, or the intermediate regime where the Doppler actor is δ p /Γ co, the peak luminosity measure is 5+α. The strong dependence on is more evident in the high σ cases (right panel in Fig. 5).,c For plasmoids with avorable orientation and, < σ w < we ind that brighter lares are powered by smaller and aster plasmoids when exiting the current sheet (see also Fig. ). The extent o this intermediate region o plasmoid sizes is larger or higher σ. For >,c, more luminous lares are powered by larger plasmoids that leave the current sheet with non-relativistic speeds (in the jet s rest rame). A similar trend is ound or lares powered by plasmoids with <,. The eect o the layer s orientation σ

9 Blazar lares rom plasmoids 9 on the peak luminosity is negligible or lares produced by monster plasmoids with.l. The results presented in Fig. 5 are obtained or a ixed viewing angle o the observer (θ obs =.5/Γ j ). The dependence o the peak luminosity on the orientation o the plasmoid relative to the observer is illustrated in Fig. 6 or σ = and two plasmoids: a small, ast-moving with =.4l (let panel) and a large, slow-moving (monster plasmoid) with =.l (right panel). A wide range o values or the peak luminosity measure is ound or smaller plasmoids due to their relativistic motion. As larger plasmoids leave the layer without becoming relativistic the angular dependence o the peak luminosity measure is less pronounced. An on-axis observer (θ obs.9 /Γ j, ) would see bright lares powered by a small plasmoid (let panel) or angles θ ranging roughly rom 3 to. In contrast, an o-axis observer (θ obs.9 /Γ j, ) would observe bright lares rom plasmoids moving at θ > 9. Relatively bright lares may thereore be expected rom AGN jets whose axis is misaligned with our line o sight. The intra-day TeV variability observed rom M87 could be such an example (Giannios, Uzdensky & Begelman ). The maximum possible peak luminosity measure or a smaller plasmoid is 3 times larger than or a rare, monster plasmoid, as indicated by the colour bars. Because o the relativistic motion o the smaller plasmoids the range o θ that avors the brightest lares is smaller compared to the larger plasmoids. Especially, or monster plasmoids and an on-axis observer almost all orientations θ -8 would lead to lares o similar peak luminosity. 4. Flux-doubling timescale Among the dierent timescales characterizing a lare, the luxdoubling timescale is one that can be directly related to the properties o the plasmoid powering the lare. This is deined as the time needed (in the observer s rame) or the luminosity (or, lux) to increase rom hal its peak value L pk / to its peak value L pk. As mentioned in the previous section, a plasmoid-powered lare reaches its peak luminosity when the plasmoid leaves the current sheet with a terminal size and Doppler actor δ p ( ). The size o the plasmoid, /, at the moment the luminosity reaches hal o its peak value is determined by the condition [ δp ( / )] 3+α w / =.5 [ δ p ( )] 3+α w, (7) where we used eq. (4). I the plasmoid s Doppler actor is constant, the above condition leads to / = w /. This is a good approximation or plasmoids with w, σ or or larger, nonrelativistically moving plasmoids, i.e. w > w,c ; in both cases, the plasmoid momentum is not signiicantly evolving over the relevant time window. The lux-doubling timescale is then calculated using eq. () as τ / ( + z) = g( w) d w cβ g δ p ( w). (8) / I the acceleration o the plasmoid and the suppression o its growth rate are ignored, the above expression results in τ /,apr ( + z) ( δ p ( )β ), (9) gc which is a naive estimate o τ / as it underestimates the actual lux-doubling timescale (see Appendix B). The lux-doubling timescale as determined by eq. (8) or σ = 3, and 5 is illustrated in Fig. 7 or a ixed viewing angle θ obs =.5/Γ j and dierent angles θ marked on the plot. A wide range o lux-doubling timescales is expected or plasmoidpowered lares. For the particular choice o l = 6 cm, we ind τ / rom min to days. For small and relativistically moving plasmoids the typical timescale may vary by two orders o magnitude among plasmoids with dierent θ. In addition, shorter τ / are obtained or the higher σ cases as a result o higher Doppler actors. For >,c (see eq. (5)), where both the Doppler actor and the growth rate are constant, the lux-doubling timescale scales linearly with respect to or all angles θ. A similar trend is also expected or suiciently small plasmoids that have been accelerated to β A c beore exiting the layer (, ). The growth rate o such plasmoids is also constant but suppressed by a actor o three or our choice o the σ suppression actor g(x / ). For the intermediate range o sizes, w σ <,c the lux-doubling timescale depends sensitively on the acceleration and growth history o the plasmoids mainly through δ p. The exact shape o the curves in this intermediate range o sizes, which extends towards smaller or higher σ, since, σ σ /, depends mainly on the plasmoid acceleration and angle θ. The dependence on the latter is evident in Fig. 7 or a given σ value. Similarly, the curves or non-avorable orientations (e.g. θ = ; black solid curves) dier among dierent magnetizations, as a result o dierences in the plasmoid acceleration. The lux-doubling timescale o lares powered by plasmoids with avorable orientation (blue and orange lines) and.3 < /l <.3 changes by at least one order o magnitude, as illustrated in Fig. 7. On the contrary, the peak luminosity o such lares varies by a actor o a ew only (see Fig. 5). For plasmoids whose radiation is beamed away rom the observer (e.g. θ = ; black solid curves) both the observed luminosity and lux-doubling timescales are strongly dependent on. The two-dimensional maps o the lux-doubling timescale or σ = and two plasmoids with =.4l (let) and.l (right) are shown in Fig. 8. We ind that the τ / o lares powered by larger plasmoids is approximately independent o θ or all viewing angles (right), whereas it strongly depends upon the combination θ and θ obs or lares produced by smaller plasmoids (let). Inspection o Figs. 6 and 8 reveals a clear anti-correlation between the peak luminosity and lux-doubling timescale. Observing ultra-ast ( min) and bright lares (as e.g. in PKS (Aharonian & al. 7) and Mrk 5 (Albert et al. 7)), requires avorable orientation o the current sheet towards the observer (see also Giannios, Uzdensky & Begelman 9). 5 APPLICATION TO THE BLAZAR EMISSION Jets are likely to be launched as Poynting-lux dominated lows (e.g. Blandord & Znajek 977; Blandord & Payne 98) and remain such ater their acceleration and collimation (e.g. Spruit 996; Vlahakis & Königl 4). In this picture reconnection is a natural candidate or the dissipation o energy and emission (or reasons why shocks are disavored, see also SPG5) provided that the magnetic ield reverses polarity in the jet. Such reversals may result rom MHD instabilities. Current driven instabilities, in particular, re-organize globally the jet s magnetic ield and may lead to current sheets occupying a large raction o the jet s cross section. However, the distance at which these instabilities may develop in the jet and the resulting size o the current sheets remain poorly understood. Alternatively, the jet could be a striped wind (e.g.

10 Petropoulou, Giannios, & Sironi log[ τ / ] (hr) θ = θ =π/6 θ =π/4 θ =π/7 θ =π/ σ= log [w /l ] log[ τ / ] (hr) σ= θ = -.8 θ =π/6 θ =π/4 -. θ =π/7 θ =π/ log [w /l ] log[ τ / ] (hr) σ=5 θ = θ =π/6 θ =π/4 θ =π/7 θ =π/ log [w /l ] Figure 7. Log-log plot o the lux-doubling timescale τ /, as determined by eq. (8), versus /l (solid lines) or the same parameters as in Fig., l = 6 cm, and dierent values o θ marked on the plot. The variability timescale does not depend much on θ or large pasmoids in contrast to smaller ones. Figure 8. Two-dimensional map o the lux-doubling timescale or σ = and two plasmoids with =.4l (let) and.l (right). Each map is created or dierent viewing angles θ obs and directions o plasmoid s motion with respect to the jet axis θ (in the jet s rest rame). Here, l = 6 cm is assumed. Lovelace, Newman & Romanova 997; Lyubarsky & Kirk ; Drenkhahn & Spruit ; McKinney & Uzdensky ). I the magnetic ield reverses polarity on a short timescale R g /c ( as, e.g., motivated by Parrey, Giannios & Beloborodov 5), where R g 3 4 M 9 cm is the gravitational radius o a M 9 = M/ 9 M black hole, the resulting stripes would have width 3 5 M 9 cm (in the black hole rame) and the scale o the reconnection layers (in the jet rame) would be l = 3 6 Γ j, M 9 cm. It can then be shown that the dissipation takes place at a distance z diss Γ j l /β rec pc, where β rec. is the reconnection speed (Giannios 3). The location o the blazar zone is hotly debated and is proposed to lie either at sub-pc scales (e.g. Ghisellini et al. b; Tavecchio et al. ; Nalewajko, Begelman & Sikora 4) or at multi-pc scales (e.g. Sikora, Moderski & Madejski 8; Marscher et al. ; Tavecchio et al. 3; Petropoulou & Dimitrakoudis 5). Regardless, the pc-scale location lies in middle o most estimates. Thus, we adopt z diss = pc and l = 6 cm as characteristic values. The absolute power o a two-sided jet is written as (e.g. Celotti & Ghisellini 8; Dermer & Menon 9) L j = πϖ cβ j Γ j (u j + P j), (3) where ϖ zθ j is the jet s cross section, β j, u j and P j u j are the energy density and pressure o the unreconnected jet low. The energy density in the jet s rame is given by u j.6 L j,46 ϖ 6 Γ j, erg cm 3. (3) and is approximately equal to the energy density in the plasmoid s rest rame, i.e. u j u j (SGP6). The blazar jet composition is uncertain but several observations (e.g. Celotti & Ghisellini 8; Ghisellini et al. a, 4) are compatible with baryon-loaded jets with several pairs per proton. Since relativistic reconnection generally results in the same energy per particle regardless o the particle species, we neglect the proton contribution to the total en-

11 Blazar lares rom plasmoids ergy density o a plasmoid, which is taken to be (moderately) dominated by leptons over the magnetic ield (SPG5). This translates to u j u e u B. 5. Observables Here we derive simple analytical expressions or two observables o lares, namely their peak luminosity L pk and lux-doubling timescale τ /. We neglect the evolution o β g with the plasmoid size (see 3), we assume perect alignment between the layer and the observer and consider lares produced by two types o plasmoids: (i) those that exit the layer with β co Γ co = and size w,c.β acc, l (see eq. (5)) and those with β co Γ co = σ and size w, σ =.3β acc, l / σ ; or the deinition see eq. (6). In general, the bolometric peak lare luminosity can be written as L pk,bol = π β gcw δ4 p, u e (3) where u e rec u j = rec L j /4πϖ cβ j Γ j. The expression or the peak bolometric luminosity is derived under the assumption that the radiating particles are ast cooling (t cool t dyn ). This is, in general, true or monster plasmoids (see 5.), but or smaller plasmoids the expressions or the peak bolometric luminosities should be treated as upper limits. (i) Plasmoids with w,c : substitution o w,c in the above expression results in L pk,bol = 3 rec L j Γ j β ( ) gβ accβ rec Γj θ j, (33) where we also used l /ϖ β rec /Γ j θ j and δ p, 4Γ j (see ootnote ollowing eq. (5)). The lux-doubling timescale can be calculated using eq. (8) and is given by τ / β accl 4Γ j β g c.3 β acc, l 6 β g, Γ j, hr. (34) The energy release during τ / can be approximated by E bol L pk,bol τ / 8 c recl j Γ j ( Γj θ j ) β 3 acc β recl. (35) Interestingly, the energy release during lares powered by monster plasmoids (w,c =.β acc, l ) does not depend strongly on the magnetization o the jet, since β acc, β rec, and rec are almost independent o σ, as shown in SPG5 and SGP6. It is the large-scale properties o the jet (e.g. Γ j, L j ) and the length o the reconnecting layer that determine the luence o such lares. (ii) Plasmoids with w, σ : substitution o w, σ in eq. (3) leads to ( ) L pk,bol = 3 rec L j Γ j σβas g β accβ rec Γj θ j, (36) where β as g β g /3 (the growth o small plasmoids is suppressed by a actor o 3, or our choice o g(x / )) and β g listed in Table. This is to be compared to eq. (33) where β g appears, since the growth rate is not suppressed or non-relativistic plasmoids. Thus, the bolometric peak luminosity o lares produced by small and ast moving plasmoids is higher by a actor o σ(β as g /β g ) σ/3 compared to L pk,bol o lares produced by the larger and non-relativistic plasmoids. The lux-doubling timescale is given by τ / l β acc 4Γ j σβ as g c l 6 4 β acc, σ Γ j, β as min, (37) g, where we used δ p, 4Γ j Γ co 4 σγ j. The luence o a lare over τ / is written as E bol L pk,bol τ / 8 c recl j Γ j ( Γj θ j ) β 3 acc β recl, (38) which is same as eq. (35). Thus, the observed luence o lares powered by either monster plasmoids or much smaller plasmoids with β co Γ co σ that orm in the same reconnecting layer is nearly the same despite o the large dierences in τ / and L pk. Interestingly, the E bol is nearly independent o the plasma magnetization, since β acc and β rec do not vary much with σ. Observations o ultra-ast ( min) and slower varying ( hr) lares o similar luence during a longer period o blazar laring activity (e.g. lasting or several hours to days) would indicate the emission rom the same current sheet. The ratio o the luxdoubling timescales o these lares is σ (see eq. (34) and (37)), i.e. it can be used to iner the magnetization on the jet. The inerred σ could be potentially cross-checked with the inerred slope o the particle distribution, which also depends on the magnetization (see 5.). For avorable orientations, both ast and slower varying lares are expected to have similar peak lare luminosities, so they are likely to be both detectable. The observed timescale required or the energy to be dissipated in an active region o the jet with (comoving) volume l 3 deines the total observed duration o a single reconnection event and is given by l /Γ j β rec c 88 hr l 6 Γ j,. The dissipated energy will be radiated away in multiple bursts over this longer period. The largest (monster) plasmoids leave the layer at nonrelativistic or mildly-relativistic speeds giving rise to lares o duration l /δ j c 8 hr l 6 δ j, (see e.g. eq. (34) and Fig. 7). Several lares with such durations may be expected rom a single reconnection event, since monster plasmoids are ejected rom the current sheet every ew dynamical times (see also SGP6). Smaller plasmoids may produce much aster lares depending on the orientation o the observer and the reconnecting layer. For avorable orientations, several bright and ultra-ast ( min) lares may be produced by smaller plasmoids that are requently ormed in the reconnection layer. Detection o multiple lares within a period o increased activity could serve as a probe o the layer s orientation with respect to the observer. This is exempliied in Fig. 9 where the peak (bolometric) luminosity o lares produced by plasmoids o dierent sizes is plotted against the respective lux-doubling timescale. The results are shown or σ =, rec =.5, θ obs =.5/Γ j, and various angles θ. The growth rate is assumed to be constant and equal to β g /( + tanh(β co, /β A ), where β co, is the inal plasmoid velocity. Here, the jet parameters are l = 6 cm, L j = 46 erg s, Γ j =, and Γ j θ j =. (Pushkarev et al. 9). Similar results are expected or other magnetizations except or the minimum τ / which will be lower or higher σ. Data rom lares with a range o observed rise timescales could be included in Fig. 9 in order to iner θ or an on-axis observer. 5. Indicative examples Henceorth, we adopt u e 3u B as inerred rom PIC simulations or all σ 3 (see also SPG5) and, as an example, we set u e =.5 erg cm 3. This translates to u B =.5 erg cm 3 or B p G and corresponds to L j = 46 erg s or rec =.5, β rec =., Γ j =, Γ j θ j =., and l = 6 cm. The adopted energy densities are within the ranges inerred rom leptonic modelling o blazar emission (e.g. Celotti & Ghisellini 8; Tavecchio & Ghisellini