Blazar Zone Research Articles

Radiative properties of plasmoids and plasmoid mergers in magnetic reconnection

ABSTRACT Magnetic reconnection is often considered as the primary particle acceleration mechanism in a magnetized blazar zone environment. The majority of radiation in the reconnection layer comes from plasmoids and their mergers. In particular, plasmoid mergers can produce strong multiwavelength flares and major variations in synchrotron polarization signatures. However, radiative properties of plasmoid mergers have not been well explored due to difficulties in tracking the merging processes. Here we use an image processing method that combines the magnetic vector potential and density to identify isolated and merging plasmoids. We find that this method can clearly distinguish radiation contributions from isolated plasmoids, merging plasmoids, and the primary current sheet of reconnection. This new method enables us to study the radiative properties of plasmoids and mergers statistically. Our results show that isolated plasmoids have similar emissivity regardless of their sizes, and they generally have non-zero polarization degree (PD) due to their quasi-circular shape. Flares due to plasmoid mergers have relative amplitudes that are antiproportional to the size ratio of the plasmoids participating in the mergers. Finally, only mergers between plasmoids of comparable sizes (width ratio ≲5) can lead to significant spectral hardening and polarization angle (PA) variations; the amplitude of the PA variations is between 0 and 180° and has a mean value of 90°. Our analyses on 2D simulations can pave the way for future analyses and machine learning techniques on radiative properties of 3D magnetic reconnection simulations.

A ring accelerator? Unusual jet dynamics in the IceCube candidate PKS 1502+106

ABSTRACT On 2019/07/30.86853 ut, IceCube detected a high-energy astrophysical neutrino candidate. The Flat Spectrum Radio Quasar PKS 1502+106 is located within the 50 per cent uncertainty region of the event. Our analysis of 15 GHz Very Long Baseline Array (VLBA) and astrometric 8 GHz VLBA data, in a time span prior and after the IceCube event, reveals evidence for a radio ring structure that develops with time. Several arc-structures evolve perpendicular to the jet ridge line. We find evidence for precession of a curved jet based on kinematic modelling and a periodicity analysis. An outflowing broad line region (BLR) based on the C iv line emission (Sloan Digital Sky Survey) is found. We attribute the atypical ring to an interaction of the precessing jet with the outflowing material. We discuss our findings in the context of a spine-sheath scenario where the ring reveals the sheath and its interaction with the surroundings (narrow line region, NLR, clouds). We find that the radio emission is correlated with the γ-ray emission, with radio lagging the γ-rays. Based on the γ-ray variability time-scale, we constrain the γ-ray emission zone to the BLR (30–200 rg) and within the jet launching region. We discuss that the outflowing BLR provides the external radiation field for γ-ray production via external Compton scattering. The neutrino is most likely produced by proton–proton interaction in the blazar zone (beyond the BLR), enabled by episodic encounters of the jet with dense clouds, i.e. some molecular cloud in the NLR.

Radiation signatures from striped blazar jet

ABSTRACT Relativistic jets from supermassive black holes are among the most powerful and luminous astrophysical systems in Universe. We propose that the open magnetic field lines through the black hole, which drive a strongly magnetized jet, may have their polarity reversing over time scales related to the growth of the magnetorotational dynamo in the disc, resulting in dissipative structures in the jet characterized by reversing toroidal field polarities, referred to as ‘stripes’. The magnetic reconnection between the stripes dissipates the magnetic energy and powers jet acceleration. The striped jet model can explain the jet acceleration, large-scale jet emission, and blazar emission signatures consistently in a unified physical picture. Specifically, we find that the jet accelerates to the bulk Lorentz factor Γ ≳ 10 within 1-parsec distance from the central engine. The acceleration slows down but continues at larger distances, with intrinsic acceleration rate $\dot{\Gamma }/\Gamma$ between $0.0005$ and $0.005~\rm {yr^{-1}}$ at tens of parsecs, which is in very good agreement with recent radio observations. Magnetic reconnection continuously accelerates non-thermal particles over large distances from the central engine, resulting in the core-shift effect and overall flat-to-inverted synchrotron spectrum. The large-scale spectral luminosity peak νpeak is antiproportional to the location of the peak of the dissipation, which is set by the minimal stripe width lmin. The blazar zone is approximately at the same location. At this distance, the jet is moderately magnetized, with the comoving magnetic field strength and dissipation power consistent with typical leptonic blazar model parameters.

Multi-epoch Modeling of TXS 0506+056 and Implications for Long-term High-energy Neutrino Emission

Abstract The IceCube report of a excess of 13 ± 5 neutrino events in the direction of the blazar TXS 0506+056 in 2014–2015 and the 2017 detection of a high-energy neutrino event, IceCube-170922A, during a gamma-ray flare from the same blazar, have revived the interest in scenarios for neutrino production in blazars. We perform comprehensive analyses on the long-term electromagnetic emission of TXS 0506+056 using optical, X-ray, and gamma-ray data from the All-Sky Automated Survey for Supernovae, the Neil Gehrels Swift Observatory, Monitor of All-sky X-ray Image, and the Fermi Large Area Telescope. We also perform numerical modeling of the spectral energy distributions (SEDs) in four epochs prior to 2017 with contemporaneous gamma-ray and lower-energy (optical and/or X-ray) data. We find that the multi-epoch SEDs are consistent with a hybrid leptonic scenario, where the gamma-rays are produced in the blazar zone via external inverse Compton scattering of accelerated electrons, and high-energy neutrinos are produced via the photomeson production process of co-accelerated protons. The multi-epoch SEDs can be satisfactorily explained with the same jet parameters and variable external photon density and electron luminosity. Using the maximal neutrino flux derived for each epoch, we put an upper limit of ∼0.4–2 on the muon neutrino number in 10 years of IceCube observations. Our results are consistent with the IceCube-170922A detection, which can be explained as an upper fluctuation from the average neutrino rate expected from the source, but in strong tension with the 2014–2015 neutrino flare.

Long-term optical spectroscopic variations in blazar 3C 454.3

Aims. Characterisation of the long-term variations in the broad line region in a luminous blazar, where Comptonisation of broad-line emission within a relativistic jet is the standard scenario for production of γ-ray emission that dominates the spectral energy distribution. Methods. We analysed ten years of optical spectroscopic data from the Steward Observatory for the blazar 3C 454.3, as well as γ-ray data from the Fermi Large Area Telescope (LAT). The optical spectra are dominated by a highly variable non-thermal synchrotron continuum with a prominent Mg II broad emission line. The line flux was obtained by spectral decomposition including significant contribution from the Fe II pseudo-continuum. Three methods were used to characterise variations in the line flux: (1) stacking of the continuum-subtracted spectra, (2) subtracting the running mean light curves calculated for different timescales, and (3) evaluating potential time delays via the discrete correlation function (DCF). Results. Despite very large variations in the γ-ray and optical continua, the line flux changes only moderately (<0.1 dex). The data suggest that the line flux responds to a dramatic change in the blazar activity from a very high state in 2010 to a deep low state in 2012. Two interpretations are possible: either the line flux is anti-correlated with the continuum or the increase in the line luminosity is delayed by ∼600 days. If this time delay results from the reverberation of poorly constrained accretion disc emission in both the broad-line region (BLR) and the synchrotron emitting blazar zone within a relativistic jet, we would obtain natural estimates for the BLR radius RBLR, MgII ≳ 0.28 pc and for the supermassive black hole mass MSMBH ∼ 8.5 × 108 M⊙. We did not identify additional examples of short-term “flares” of the line flux, in addition to the previously reported case observed in 2010.

BARYON LOADING EFFICIENCY AND PARTICLE ACCELERATION EFFICIENCY OF RELATIVISTIC JETS: CASES FOR LOW LUMINOSITY BL LACS

ABSTRACT Relativistic jets launched by supermassive black holes, so-called active galactic nuclei (AGNs), are known as the most energetic particle accelerators in the universe. However, the baryon loading efficiency onto the jets from the accretion flows and their particle acceleration efficiencies have been veiled in mystery. With the latest data sets, we perform multi-wavelength spectral analysis of quiescent spectra of 13 TeV gamma-ray detected high-frequency-peaked BL Lacs (HBLs) following one-zone static synchrotron self-Compton (SSC) model. We determine the minimum, cooling break, and maximum electron Lorentz factors following the diffusive shock acceleration (DSA) theory. We find that HBLs have and the radiative efficiency , where P B and P e is the Poynting and electron power, respectively. By assuming 10 leptons per one proton, the jet power relates to the black hole mass as , where and is the jet power and the Eddington luminosity, respectively. Under our model assumptions, we further find that HBLs have a jet production efficiency of and a mass loading efficiency of . We also investigate the particle acceleration efficiency in the blazar zone by including the most recent Swift/BAT data. Our samples ubiquitously have particle acceleration efficiencies of , which is inefficient to accelerate particles up to the ultra-high-energy-cosmic-ray (UHECR) regime. This implies that the UHECR acceleration sites should not be the blazar zones of quiescent low power AGN jets, if one assumes the one-zone SSC model based on the DSA theory.

A spine–sheath model for strong-line blazars

We have developed a quasi-analytical model for the production of radiation in strong-line blazars, assuming a spine-sheath jet structure. The model allows us to study how the spine and sheath spectral components depend on parameters describing the geometrical and physical structure of "the blazar zone". We show that typical broad-band spectra of strong-line blazars can be reproduced by assuming the magnetization parameter to be of order unity and reconnection to be the dominant dissipation mechanism. Furthermore, we demonstrate that the spine-sheath model can explain why gamma-ray variations are often observed to have much larger amplitudes than the corresponding optical variations. The model is also less demanding of jet power than one-zone models, and can reproduce the basic features of extreme gamma-ray events.

Magnetization of jets in luminous blazars

Luminosities of many powerful blazars are strongly dominated by $\gamma$-rays which most likely result from Comptonization of radiation produced outside a jet. This observation sets certain constraints on composition and energetics of the jet as well as the surrounding quasar environment. We study the dependence of a Compton dominance on a jet magnetization (the magnetic-to-matter energy flux) and on the location of the 'blazar zone'. Calculations are performed for two geometries of broad-emission-line and hot-dust regions: spherical and planar. The jet magnetization corresponding to the large observed Compton dominance is found to be $\sim 0.1 (\theta_{\rm j} \Gamma)^2$ for spherical geometries and $\sim 0.01 (\theta_{\rm j} \Gamma)^2$ for planar geometries, where $\theta_{\rm j}$ is the jet half-opening angle and $\Gamma$ is the jet Lorentz factor. This implies that jets in luminous blazars are matter dominated and that this domination is particularly strong for the flattened geometry of external radiation sources.

EXPLORING THE BLAZAR ZONE IN HIGH-ENERGY FLARES OF FSRQs

The gamma-ray emission offers a powerful diagnostic tool to probe jets and their surroundings in flat spectrum radio quasars (FSRQ). In particular, sources emitting at high energies (>10 GeV) give us the strongest constraints. This motivates us to start a systematic study of flares with bright emission above 10 GeV, examining archival data of Fermi-LAT gamma-ray telescope. At the same time, we began to trigger Target of Opportunity observations to the Swift observatory at the occurrence of high-energy flares, obtaining a wide coverage of the spectral energy distributions for several FSRQs during flares. Among the others we investigate the SED of a peculiar flare of 3C 454.3, showing a remarkable hard gamma-ray spectrum, quite different from the brightest flares of this source, and a bright flare of CTA 102. We modeled the SED in the framework of the one--zone leptonic model, using also archival optical spectroscopic data to derive the luminosity of the broad lines and thus estimate the disk luminosity, from which the structural parameters of the FSRQ nucleus can be inferred. The model allowed us to evaluate the magnetic field intensity in the blazar zone, and to locate the emitting region of gamma rays in the particular case in which gamma-ray spectra show neither absorption from the BLR, nor the Klein-Nishina curvature expected in leptonic models assuming the BLR as source of seed photons for the External Compton. For FSRQs bright above 10 GeV, we where able to identify short periods lasting less than 1 day characterized by high rate of high energy gamma rays, and hard gamma-ray spectra. We discussed the observed spectra and variability timescales in terms of injection and cooling of energetic particles, arguing that these flares could be triggered by magnetic reconnections events or turbulence in the flow.

The origin and emission mechanism of VHE (>100GeV) emission from FSRQs

Flat Spectrum Radio Quasars, unlike BL Lac objects, are blazars that show prominent line-emission and strong thermal components associated with the accretion disk, the broad-line region (BLR), and/or the dusty torus. The low energy peak in the continuum is from synchrotron emission (of electrons), and the high energy peak is well explained by external-Compton emission. In these models the relativistic electrons in the jet up-scatter photons from the thermal photon fields up to GeV energies. Beyond a few tens of GeV such models predict cutoffs due to Klein-Nishina effect and internal absorption via pair production. While more than 300 FSRQs have been seen with Fermi -LAT (between 100MeV−30GeV), only three have been detected at VHE (Very High Energy, E > 100 GeV) with Cherenkov telescopes. The detection of VHE emission constrains the location of the blazar zone based on internal absorption estimates, but challenges the emission models that predict cutoffs. While a number of GeV flaring states (in various FSRQs) have been observed with Cherenkov telescopes only few have resulted in detection of a VHE signal. The broadband emission characteristics of VHE FSRQs (including the VHE-detected FSRQs) are studied and put in context to better understand their location and emission mechanism.

A NOVEL APPROACH IN CONSTRAINING ELECTRON SPECTRA IN BLAZAR JETS: THE CASE OF Markarian 421

We report results from the observations of the well studied TeV blazar Mrk 421 with the Swift and the Suzaku satellites in December 2008. During the observation, Mrk 421 was found in a relatively low activity state, with the corresponding 2-10 keV flux of $3 \times 10^{-10}$ erg/s/cm^2. For the purpose of robust constraining the UV-to-X-ray emission continuum we selected only the data corresponding to truly simultaneous time intervals between Swift and Suzaku, allowing us to obtain a good-quality, broad-band spectrum despite a modest length (0.6 ksec) exposure. We analyzed the spectrum with the parametric forward-fitting SYNCHROTRON model implemented in XSPEC assuming two different representations of the underlying electron energy distribution, both well motivated by the current particle acceleration models: a power-law distribution above the minimum energy $\gamma_{\rm min}$ with an exponential cutoff at the maximum energy $\gamma_{\rm max}$, and a modified ultra-relativistic Maxwellian with an equilibrium energy $\gamma_{\rm eq}$. We found that the latter implies unlikely physical conditions within the blazar zone of Mrk 421. On the other hand, the exponentially moderated power-law electron distribution gives two possible sets of the model parameters: (i) flat spectrum $dN'_e/d\gamma \propto \gamma^{-1.91}$ with low minimum electron energy $\gamma_{\rm min}<10^3$, and (ii) steep spectrum $\propto \gamma^{-2.77}$ with high minimum electron energy $\gamma_{\rm min}\simeq 2\times10^4$. We discuss different interpretations of both possibilities in the context of a diffusive acceleration of electrons at relativistic, sub- or superluminal shocks. We also comment on how exactly the gamma-ray data can be used to discriminate between the proposed different scenarios.

GeV BREAKS IN BLAZARS AS A RESULT OF GAMMA-RAY ABSORPTION WITHIN THE BROAD-LINE REGION

Spectra of the brightest blazars detected by the Fermi Gamma-ray Space Telescope Large Area Telescope cannot be described by a simple power law model. A much better description is obtained with a broken power law, with the break energies of a few GeV. We show here that the sharpness and the position of the breaks can be well reproduced by absorption of gamma-rays via photon--photon pair production on HeII Lyman recombination continuum and lines. This implies that the blazar zone lies inside the region of the highest ionization of the broad-line region (BLR) within a light-year from a super-massive black hole. The observations of gamma-ray spectral breaks open a way of studying the BLR photon field in the extreme-UV/soft X-rays, which are otherwise hidden from our view.

3C 454.3 Reveals the Structure and Physics of Its “Blazar Zone”

Recent multi-wavelength observations of 3C454.3, in particular during its giant outburst in 2005, put severe constraints on the location of the 'blazar zone', its dissipative nature, and high energy radiation mechanisms. As the optical, X-ray, and millimeter light-curves indicate, significant fraction of the jet energy must be released in the vicinity of the millimeter-photosphere, i.e. at distances where, due to the lateral expansion, the jet becomes transparent at millimeter wavelengths. We conclude that this region is located at ~10 parsecs, the distance coinciding with the location of the hot dust region. This location is consistent with the high amplitude variations observed on ~10 day time scale, provided the Lorentz factor of a jet is ~20. We argue that dissipation is driven by reconfinement shock and demonstrate that X-rays and gamma-rays are likely to be produced via inverse Compton scattering of near/mid IR photons emitted by the hot dust. We also infer that the largest gamma-to-synchrotron luminosity ratio ever recorded in this object - having taken place during its lowest luminosity states - can be simply due to weaker magnetic fields carried by a less powerful jet.

The role of kink instability in Poynting-flux dominated jets

The role of kink instability in magnetically driven jets is explored through numerical one-dimensional steady relativistic MHD calculations. The instability is shown to have enough time to grow and influence the dynamics of Poynting-flux dominated jets. In the case of AGN jets, the flow becomes kinetic flux dominated at distances because of the rapid dissipation of Poynting flux. When applied to GRB outflows, the model predicts more gradual Poynting dissipation and moderately magnetized flow at distances of ~1016 cm where the deceleration of the ejecta due to interaction with the external medium is expected. The energy released by the instability can power the compact “blazar zone” emission and the prompt emission of GRB outflows with high radiative efficiencies.