Gamma-ray Flares Research Articles

A Gamma-Ray Flare from TXS 1508+572: Characterizing the Jet of a z = 4.31 Blazar in the Early Universe

Blazars can be detected from very large distances due to their high luminosity. However, the detection of γ-ray emission of blazars beyond z = 3 has only been confirmed for a small number of sources. Such observations probe the growth of supermassive black holes close to the peak of star formation in the history of galaxy evolution. As a result from a continuous monitoring of a sample of 80 z > 3 blazars with the Fermi Large Area Telescope (Fermi-LAT), we present the first detection of a γ-ray flare from the z = 4.31 blazar TXS 1508+572. This source showed high γ-ray activity from 2022 February to August, reaching a peak luminosity comparable to the most luminous flares ever detected with Fermi-LAT. We conducted a multiwavelength observing campaign involving XMM-Newton, the Neil Gehrels Swift Observatory, the Effelsberg 100 m radio telescope, and the Very Long Baseline Array. In addition, we make use of the monitoring programs by the Zwicky Transient Facility and the Near-Earth Object Wide-field Infrared Survey Explorer at optical and infrared wavelengths, respectively. We find that the source is particularly variable in the infrared band on daily timescales. The spectral energy distribution collected during our campaign is well described by a one-zone leptonic model, with the γ-ray flare originating from an increase of external Compton emission as a result of a fresh injection of accelerated electrons.

Redshift constrain of BL Lac PKS 1424+240

ABSTRACT In the period between 2009 and 2015, several very high-energy (VHE $> 100$ GeV) gamma-ray flaring events from the BL Lac object PKS 1424+240 were observed by the Cerenkov telescopes VERITAS and MAGIC. It had uncertain redshift (z) and using spectroscopical measurement, Paiano et al. (2017) found it to be $z=0.604$. Using four different extragalactic background light (EBL) models and the photohadronic model, nine independently observed VHE gamma-ray spectra of PKS 1424+240 are analysed and a global $\chi ^2$ fit is performed on all observations to estimate the best-fitting value for the redshift for each EBL model. Confidence levels (CL) intervals for the redshift are also estimated using all the EBL models. This method is tested by comparing our analysis with the observed value. It is shown that the photohadronic scenario provides an excellent description of all the observed spectra. It is found that the EBL model of Dominguez et al. (2011) is the one that provides the most restrictive limits on the redshift of PKS 1424+240, but in our analysis, $z=0.604$ lies within the $3\sigma$ CL interval of the EBL model of Saldana-Lopez et al. (2021).

Separating source-intrinsic and Lorentz invariance violation induced delays in the very high-energy emission of blazar flares

The aim of the present study is to explore how to disentangle energy-dependent time delays due to a possible Lorentz invariance violation (LIV) at Planck scale from intrinsic delays expected in standard blazar flares. We first characterised the intrinsic time delays in BL Lacs and flat-spectrum radio quasars in standard one-zone time-dependent synchrotron self-Compton or external Compton models, during flares produced by particle acceleration and cooling processes. We simulated families of flares with both intrinsic and external LIV-induced energy-dependent delays. Discrimination between intrinsic and LIV delays is then investigated in two different ways. A technique based on Euclidean distance calculation between delays obtained in the synchrotron and in the inverse-Compton spectral bumps is used to assess their degree of correlation. A complementary study is performed using spectral hardness versus intensity diagrams in both energy ranges. We show that the presence of non-negligible LIV effects, which essentially act only at very high energy (VHE), can drastically reduce the strong correlation expected between the X-ray and the VHE gamma-ray emission in leptonic scenarios. The LIV phenomenon can then be hinted at measuring the Euclidean distance $d_ E $ from simultaneous X-ray and gamma-ray flare monitoring. Large values of minimum distance $d_ E,min $ would directly indicate the influence of non-intrinsic time delays possibly due to LIV in SSC flares. LIV effects can also significantly modify the VHE hysteresis patterns in hardness-intensity diagrams and even change their direction of rotation compared to the X-ray behaviour. Both observables could be used to discriminate between LIV and intrinsic delays, provided high-quality flare observations are available.

Establishing accretion flares from supermassive black holes as a source of high-energy neutrinos

ABSTRACT The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6σ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ∼10 per cent of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares.

Follow-up on the Supermassive Black Hole Binary Candidate J1048+7143: Successful Prediction of the Next Gamma-Ray Flare and Refined Binary Parameters in the Framework of the Jet Precession Model

Analyzing single-dish and very long baseline interferometry radio, as well as Fermi Large Area Telescope γ-ray observations, we explained the three major flares in the γ-ray light curve of FSRQ J1048+7143 with the spin–orbit precession of the dominant mass black hole in a supermassive black hole binary system. Here, we report on the detection of a fourth γ-ray flare from J1048+7143, appearing in the time interval that was predicted in our previous work. Including this new flare, we constrained the mass ratio into a narrow range of 0.062 < q < 0.088, and consequently we were able to further constrain the parameters of the hypothetical supermassive binary black hole at the heart of J1048+7143. We predict the occurrence of the fifth major γ-ray flare that would appear only if the jet will still lay close to our line of sight. The fourth major γ-ray flare also shows the two-subflare structure, further strengthening our scenario in which the occurrence of the subflares is the signature of the precession of a spine–sheath jet structure that quasiperiodically interacts with a proton target, e.g., clouds in the broad-line region.

Идентификация транзиентных эффектов в атмосфере Земли при помощи малых космических аппаратов

The paper considers the dynamic aspects of the possibility of carrying out a space mission to study gamma-ray flares of terrestrial origin. We propose an approach to the selection of reference orbits that allows the satellites to form equilateral triangle in the near-equatorial zone. In addition, an algorithm for phasing the satellites by latitude argument using aerodynamic drag forces is proposed.

Solar flares with sustained gamma-ray emission and some characteristics of high-energy proton fluxes

We describe the characteristics of long-term gamma fluxes with quantum energies &gt;100 MeV, obtained from Fermi/LAT data during the impulsive phase of the most energetic flare phenomena. We compare GOES data on proton fluxes with energies &gt;500 MeV with Fermi/LAT data on gamma fluxes for the period 2010–2018. The results of the analysis of all data obtained on 32 gamma-ray flares from the Fermi/LAT catalog show that the flare phenomena can be classified into three different types: type 1 — gamma fluxes accompanied by energetic proton fluxes; type 2 — gamma rays recorded in the absence of increases in proton fluxes; type 3 — gamma fluxes not recorded during observed increases in energetic proton fluxes. The burst character of energy release in the hard X-ray range was noted in some flares.

Солнечные вспышки с продолжительным гамма-излучением и характеристики потоков протонов высоких энергий

We describe the characteristics of long-term gamma fluxes with quantum energies &gt;100 MeV, obtained from Fermi/LAT data during the impulsive phase of the most energetic flare phenomena. We compare GOES data on proton fluxes with energies &gt;500 MeV with Fermi/LAT data on gamma fluxes for the period 2010–2018. The results of the analysis of all data obtained on 32 gamma-ray flares from the Fermi/LAT catalog show that the flare phenomena can be classified into three different types: type 1 — gamma fluxes accompanied by energetic proton fluxes; type 2 — gamma rays recorded in the absence of increases in proton fluxes; type 3 — gamma fluxes not recorded during observed increases in energetic proton fluxes. The burst character of energy release in the hard X-ray range was noted in some flares.

Scale Invariance in Gamma-Ray Flares of the Sun and 3C 454.3

Using the gamma-ray flare samples of the Sun and 3C 454.3 observed by the Fermi telescope, we investigate the statistical properties of sizes including fluence (energy), peak flux (luminosity), duration time, and waiting time in this work. We find that the cumulative distribution of the fluctuations of these sizes follow well the Tsallis q-Gaussian function. The obtained q values from q-Gaussian distribution remain stable around 2 without any significant change, implying that there is a scale invariance structure in gamma-ray flares of the Sun and 3C 454.3. This scale invariance characteristics of the Sun and 3C 454.3 indicated by q values are also comparable to those of earthquakes, soft gamma repeaters, fast radio burst (FRB 20121102), and X-ray flares of gamma-ray bursts. On top of that, we verify the relationship between q values and the power-law indices α from the size frequency distributions, which is expressed as q = (α + 2)/α. These statistical findings could be well explained within the physical framework of a self-organizing criticality system.

High-energy Radiation and Ion Acceleration in Three-dimensional Relativistic Magnetic Reconnection with Strong Synchrotron Cooling

We present the results of 3D particle-in-cell simulations that explore relativistic magnetic reconnection in pair plasma with strong synchrotron cooling and a small mass fraction of nonradiating ions. Our results demonstrate that the structure of the current sheet is highly sensitive to the dynamic efficiency of radiative cooling. Specifically, stronger cooling leads to more significant compression of the plasma and magnetic field within the plasmoids. We demonstrate that ions can be efficiently accelerated to energies exceeding the plasma magnetization parameter, ≫σ, and form a hard power-law energy distribution, f i ∝ γ −1. This conclusion implies a highly efficient proton acceleration in the magnetospheres of young pulsars. Conversely, the energies of pairs are limited to either σ in the strong cooling regime or the radiation burnoff limit, γ syn, when cooling is weak. We find that the high-energy radiation from pairs above the synchrotron burnoff limit, ε c ≈ 16 MeV, is only efficiently produced in the strong cooling regime, γ syn < σ. In this regime, we find that the spectral cutoff scales as ε cut ≈ ε c (σ/γ syn) and the highest energy photons are beamed along the direction of the upstream magnetic field, consistent with the phenomenological models of gamma-ray emission from young pulsars. Furthermore, our results place constraints on the reconnection-driven models of gamma-ray flares in the Crab Nebula.

High-energy neutrino-induced cascade from the direction of the flaring radio blazar TXS 0506 + 056 observed by Baikal-GVD in 2021

ABSTRACT The existence of high-energy astrophysical neutrinos has been unambiguously demonstrated, but their sources remain elusive. IceCube reported an association of a 290-TeV neutrino with a gamma-ray flare of TXS 0506 + 056, an active galactic nucleus with a compact radio jet pointing to us. Later, radio-bright blazars were shown to be associated with IceCube neutrino events with high statistical significance. These associations remained unconfirmed with the data of independent experiments. Here, we report on the detection of a rare neutrino event with the estimated energy of 224 ± 75 TeV from the direction of TXS 0506 + 056 by the new Baikal Gigaton Volume Detector (Baikal-GVD) in April 2021. This event is the highest energy cascade detected so far by the Baikal-GVD neutrino telescope from a direction below horizon. The result supports previous suggestions that radio blazars in general, and TXS 0506 + 056 in particular, are the sources of high-energy neutrinos, and opens up the cascade channel for the neutrino astronomy.

Quantitative comparisons of very-high-energy gamma-ray blazar flares with relativistic reconnection models

The origin of extremely fast variability is one of the long-standing questions in the gamma-ray astronomy of blazars. While many models explain the slower, lower energy variability, they cannot easily account for such fast flares reaching hour-to-minute timescales. Magnetic reconnection, a process where magnetic energy is converted to the acceleration of relativistic particles in the reconnection layer, is a candidate solution to this problem. In this work, we employ state-of-the-art particle-in-cell simulations in a statistical comparison with observations of a flaring episode of a well-known blazar, Mrk 421, at a very high energy (VHE, E &gt; 100 GeV). We tested the predictions of our model by generating simulated VHE light curves that we compared quantitatively with methods that we have developed for a precise evaluation of theoretical and observed data. With our analysis, we can constrain the parameter space of the model, such as the magnetic field strength of the unreconnected plasma, viewing angle and the reconnection layer orientation in the blazar jet. Our analysis favours parameter spaces with magnetic field strength 0.1 G, rather large viewing angles (6 − 8°), and misaligned layer angles, offering a strong candidate explanation for the Doppler crisis often observed in the jets of high synchrotron peaking blazars.

Repeated patterns of gamma-ray flares suggest structured jets of blazars as likely neutrino sources

ABSTRACT Fermi-Large Area Telescope (LAT) observations provide continuous and regularly sampled measurements of gamma-ray photon flux for hundreds of blazars. Many of these light curves, spanning almost 15 yr, have been thoroughly examined for periodicity in multiple studies. However, the possibility that blazars may exhibit irregularly repeating flaring patterns in their gamma-ray light curves has not been systematically explored. In this study, we aim to find repeating episodes of flaring activity in the 100 brightest blazars using Fermi-LAT light curves with various integration times. We use a Bayesian Blocks representation to convert the time series into strings of symbols and search for repeating sub-strings using a fuzzy search algorithm. As a result, we identify 27 repeated episodes in the gamma-ray light curves of 10 blazars. We find that the patterns are most likely produced in structured jets composed of a fast spine and a slower sheath. When individual emission features propagate in the spine, they scatter seed photons produced in the non-uniform sheath through the inverse Compton mechanism, resulting in a set of gamma-ray flares with a similar profile every such passage. Additionally, we explore the theoretically predicted possibility that the spine-sheath structure facilitates the production of high-energy neutrinos in blazar jets. Using the catalogue of track-like events detected by the IceCube neutrino telescope, we find evidence supporting this hypothesis at a 2.8σ significance level.

Flare Duty Cycle of Gamma-Ray Blazars and Implications for High-energy Neutrino Emission

Gamma-ray flares of blazars may be accompanied by high-energy neutrinos due to interactions of high-energy cosmic rays in the jet with photons, as suggested by the detection of the high-energy neutrino IceCube-170922A during a major gamma-ray flare from blazar TXS 0506+056 at the ∼3σ significance level. In this work, we present a statistical study of gamma-ray emission from blazars to constrain the contribution of gamma-ray flares to their neutrino output. We construct weekly binned light curves for 145 gamma-ray bright blazars in the Fermi Large Area Telescope Monitored Source List adding TXS 0506+056. We derive the fraction of time spent in the flaring state (flare duty cycle) and the fraction of energy released during each flare from the light curves with a Bayesian blocks algorithm. We find that blazars with lower flare duty cycles and energy fractions are more numerous among our sample. We identify a significant difference in flare duty cycles between blazar subclasses at a significance level of 5%. Then using a general scaling relation for the neutrino and gamma-ray luminosities, with a weighting exponent of γ = 1.0–2.0, normalized to the quiescent gamma-ray or X-ray flux of each blazar, we evaluate the neutrino energy flux of each gamma-ray flare. The gamma-ray flare distribution indicates that blazar neutrino emission may be dominated by flares for γ ≳ 1.5. The neutrino energy fluxes for 1 week and 10 yr bins are compared with the decl.-dependent IceCube sensitivity to constrain the standard neutrino emission models for gamma-ray flares. Finally, we present the upper-limit contribution of blazar gamma-ray flares to the isotropic diffuse neutrino flux.

Investigating Possible Correlations between Gamma-Ray and Optical Lightcurves for TeV-Detected Northern Blazars over 8 Years of Observations

Blazars are a subclass of active galactic nuclei (AGN) having relativistic jets aligned within a few degrees of our line-of-sight and form the majority of the AGN detected in the TeV regime. The Fermi-Large Area Telescope (LAT) is a pair-conversion telescope, sensitive to photons having energies between 20 MeV and 2 TeV, and is capable of scanning the entire gamma-ray sky every three hours. Despite the remarkable success of the Fermi mission, many questions still remain unanswered, such as the site of gamma-ray production and the emission mechanisms involved. The Asteroid Terrestrial-impact Last Alert System (ATLAS) is a high cadence all sky survey system optimized to be efficient for finding potentially dangerous asteroids, as well as in tracking and searching for highly variable and transient sources, such as AGN. In this study, we investigate possible correlations between the Fermi-LAT observations in the 100 MeV–300 GeV energy band and the ATLAS optical data in the R-band, centered at 679 nm, for a sample of 18 TeV-detected northern blazars over 8 years of observations between 2015 and 2022. Under the assumption that the optical and gamma-ray flares are produced by the same outburst propagating down the jet, the strong correlations found for some sources suggest a single-zone leptonic model of emission.

High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations

ABSTRACT The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with 2D and 3D particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field B and plasma density n, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the n-B space are developed and used to visualize and quantify the effects of compression. The n-B histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.

Ultrafast Variability in AGN Jets: Intermittency and Lighthouse Effect

Gamma-ray flares from active galactic nuclei (AGNs) show substantial variability on ultrafast timescales (i.e., shorter than the light-crossing time of the AGN’s supermassive black hole). We propose that ultrafast variability is a by-product of the turbulent dissipation of the jet Poynting flux. Due to the intermittency of the turbulent cascade, the dissipation is concentrated in a set of reconnecting current sheets. Electrons energized by reconnection have a strong pitch-angle anisotropy, i.e., their velocity is nearly aligned with the guide magnetic field. Then each current sheet produces a narrow radiation beam, which dominates the emission from the whole jet when it is directed toward the observer. The ultrafast variability is set by the light-crossing time of a single current sheet, which is much shorter than the light-crossing time of the whole emission region. The predictions of our model are (i) the bolometric luminosity of ultrafast AGN flares is dominated by the inverse Compton (IC) emission as the lower-energy synchrotron emission is suppressed due to the pitch-angle anisotropy; (ii) if the observed luminosity includes a nonflaring component, the variations of the synchrotron luminosity have a small amplitude; and (iii) the synchrotron and IC emission are less variable at lower frequencies, as the cooling time of the radiating particles exceeds the light-crossing time of the current sheet. Simultaneous multiwavelength observations of ultrafast AGN flares can test these predictions.

The Game of Life on a Magnetar Crust: From Gamma-Ray Flares to FRBs

This Letter presents a model to unify the diverse range of magnetar activity, through the building and release of elastic stress from the crust. A cellular automaton drives both local and global yielding of the crust, leading to braiding of coronal loops and energy release. The model behaves like a real magnetar in many ways: giant flares and small bursts both occur, as well as periods of quiescence whose typical duration is either ≲1 yr or ∼10–30 yr. The burst energy distribution broadly follows an earthquake-like power law over the energy range 1040–1045 erg. The local nature of coronal loops allows for the possibility of high-energy and fast radio bursts from the same magnetar. Within this paradigm, magnetar observations can be used to constrain the poorly understood mechanical properties of the neutron-star crust.

Similar properties between gamma-ray emission of 3C 454.3 and solar GeV flares

ABSTRACTBased on the survey data observed by Fermi-Large Area Telescope (LAT), we elaborate the statistical characteristics of gamma-ray flares from 3C 454.3 and solar GeV flares among flare parameters, such as isotropic energy (Eγ), peak luminosity (LP), and duration time (TDuration). We find two significant correlations as: $T_{\rm Duration} \propto E_{\gamma }^{0.31\pm 0.03}$ and $L_{\rm P} \propto E_{\gamma }^{0.61\pm 0.03}$ for 3C 454.3. The exponents are in a better agreement those of the Sun from the first Fermi-LAT solar flare catalogue, namely $T_{\rm Duration,\odot } \propto E_{\gamma ,\odot }^{0.38\pm 0.08}$ and $L_{\rm P,\odot } \propto E_{\gamma ,\odot }^{0.81\pm 0.08}$. The relationship of TDuration–Eγ and LP–Eγ could be interpreted naturally as due to magnetic dissipation through reconnection. On top of that the frequency distributions of isotropic energy, peak luminosity, and duration time of gamma-ray emission of 3C 454.3 show power-law forms, and the waiting time distribution can be described by a non-stationary Poisson process. These distribution behaviours are comparable to those of the Sun, Sagittarius A*, and M87, and follow the prediction of a self-organized criticality model. All statistical properties suggest that similar to the physical process accounting for solar GeV flares and X-ray flares in supermassive black hole systems, magnetic reconnection cloud govern the energy-release, and particle acceleration process for gamma-ray flares of 3C 454.3.

Radiative Reconnection-powered TeV Flares from the Black Hole Magnetosphere in M87

Active galactic nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region near the nucleus. However, the emission mechanism for these flares is not well understood. Recent high-resolution general-relativistic magnetohydrodynamic simulations show the occurrence of episodic magnetic reconnection events that can power flares near the black hole event horizon. In this work, we analyze the radiative properties of the reconnecting current layer under the extreme plasma conditions applicable to the black hole in M87 from first principles. We show that abundant pair production is expected in the vicinity of the reconnection layer, to the extent that the produced secondary pair plasma dominates the reconnection dynamics. Using analytic estimates backed by two-dimensional particle-in-cell simulations, we demonstrate that in the presence of strong synchrotron cooling, reconnection can produce a hard power-law distribution of pair plasma imprinted in the outgoing synchrotron (up to a few tens of MeV) and the inverse-Compton signal (up to TeV). We produce synthetic radiation spectra from our simulations, which can be directly compared with the results of future multiwavelength observations of M87* flares.