Γ-ray Flux Research Articles

Swift Observations of Mrk 421 in Selected Epochs. IV. Physical Implications of X-Ray Flaring Activity and Features of Relativistic Magnetic Reconnection in 2018 April–2023 December

We present the spectral and timing results obtained during the intense observations of Mrk 421 by the Swift-based ultraviolet to X-ray instruments during 2018 April–2023 December. The source showed various strengths of X-ray flaring activity, exceeding a level of 3.5 × 10−9 erg cm−2 s−1 during the strongest 0.3–10 keV flares. Our study identifies a number of intraday brightness variability, including 61 instances that occurred within 1 ks exposures that are consistent with the shock-in-jet scenario and accompanied by significant, fast spectral changes. The source exhibited extreme spectral properties with dominance of the log-parabolic distributions of photons with energy and the frequent occurrence of hard photon indices in the 0.3–10 keV and 0.3–300 GeV bands, with the peak of synchrotron spectral energy distribution E p detected at the energies beyond 29 keV for the first time. The source showed very fast transitions of log-parabolic-to-power-law spectra, most plausibly caused by turbulence-driven relativistic magnetic reconnection. Our spectral results also demonstrate the importance of the first-order Fermi mechanism within the magnetic field of different confinement efficiencies, stochastic acceleration, transitions in the turbulence spectrum, and hadronic cascades. The X-ray, UV, and γ-ray fluxes showed a lognormal variability, which hints at the imprint of accretion disk instabilities on the blazar jet, as well as the possible presence of hadronic cascades. The UV and γ-ray variabilities demonstrated weak correlations with the X-ray flaring activity, which is not consistent with simple synchrotron self-Compton models and requires more complex particle acceleration and emission scenarios.

The Cosmological Optical Convergence: Extragalactic Background Light from TeV Gamma Rays

The intensity of the extragalactic background (EBL), the accumulated optical and infrared emissions since the first stars, is the subject of a decades-long tension in the optical band. These photons form a target field that attenuates the γ-ray flux from extragalactic sources. This Letter reports the first γ-ray measurement of the EBL spectrum at z = 0 that is purely parametric and independent of EBL evolution with redshift over a wavelength range from 0.18 to 120 μm. Our method extracts the EBL absorption imprint on more than 260 archival TeV spectra from the STeVECat catalog by marginalizing nuisance parameters describing the intrinsic emission and instrumental uncertainties. We report an intensity at 600 nm of 6.9 ± 1.9 nW m−2 sr−1 × h 70, which is indistinguishable from the intensity derived from integrated galaxy light (IGL) and compatible with direct measurements taken beyond Pluto’s orbit. We exclude with 95% confidence diffuse contributions to the EBL with an intensity relative to the IGL, f diff, greater than 20% and provide a measurement of the expansion rate of the Universe at z = 0, H0=67−6+7 km s−1 Mpc −1 × (1 + f diff), which is EBL-model independent. IGL, direct, and γ-ray measurements agree on the EBL intensity in the optical band, finally reaching a cosmological optical convergence.

Fermi Blazars in the Zwicky Transient Facility Survey: Properties of Large Optical Variations

We analyze the optical light curve data, obtained with the Zwicky Transient Facility (ZTF) survey, for 47 γ-ray blazars monitored by the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope (Fermi). These 47 sources are selected because they are among the Fermi blazars with the largest optical variations in the ZTF data. Two color–magnitude variation patterns are seen in them, with one being redder-to-stable-when-brighter (RSWB; in 31 sources) and the other being stable when brighter (in 16 sources). The patterns fit with the results recently reported in several similar studies with different data. Moreover, we find that the colors in the stable state of the sources share similar values, for which (after being corrected for the Galactic extinction) most sources are in a range of 0.4–0.55. This feature could be intrinsic and may be applied in, for example, study of the intragalactic medium. We also determine the turning points for the sources showing the RSWB pattern, after which the color changes saturate and become stable. We find a correlation between optical fluxes and γ-ray fluxes at the turning points. The physical implications of the correlation remain to be investigated, probably better with a sample of high-quality γ-ray flux measurements.

The Spectra of IceCube Neutrino (SIN) candidate sources

Context. A correlation has been reported between the arrival directions of high-energy IceCube events and γ-ray blazars classified as intermediate- and high-synchrotron-peaked BL Lacs. Subsequent studies have investigated the optical properties of these sources, compiled and analyzed public multiwavelength data, and constrained their individual neutrino emission based on public IceCube point-source data. Aims. We provide a theoretical interpretation of public multiwavelength and neutrino point source data for the 32 BL Lac objects in the sample previously associated with an IceCube alert event. We combined the individual source results to draw conclusions regarding the multimesssenger properties of the sample and the required power in relativistic protons. Methods. We performed particle interaction modeling using open-source numerical simulation software. We constrained the model parameters using a novel and unique approach that simultaneously describes the host galaxy contribution, the observed synchrotron peak properties, the average multiwavelength fluxes, and, where possible, the IceCube point source constraints. Results. We show that a single-zone leptohadronic model can describe the multiwavelength broadband fluxes from all 32 IceCube candidates. In some cases, the model suggests that hadronic emission may contribute a considerable fraction of the γ-ray flux. The required power in relativistic protons ranges from a few percent to a factor of ten of the Eddington luminosity, which is energetically less demanding compared to other leptohadronic blazar models in recent literature. The model can describe the 68% confidence level IceCube flux for a large fraction of the masquerading BL Lacs in the sample, including TXS 0506+056; whereas, for true BL Lacs, the model predicts a low neutrino flux in the IceCube sensitivity range. Physically, this distinction is due to the presence of photons from broad line emission in masquerading BL Lacs, which increase the efficiency of hadronic interactions. The predicted neutrino flux peaks between a few petaelectronvolt and 100 PeV and scales positively with the flux in the gigaelectronvolt, megaelectronvolt, X-ray, and optical bands. Based on these results, we provide a list of the brightest neutrino emitters, which can be used for future searches targeting the 10–100 PeV regime.

Two Models for the Orbital Modulation of Gamma Rays in Cyg X-3

We model the currently available γ-ray data on Cyg X-3 from the Fermi Large Area Telescope. Thanks to Cyg X-3’s very strong γ-ray activity during 2018–2021, the data quality has significantly improved. We study the strong orbital modulation of the γ-rays observed at high γ-ray fluxes. The modulation, as found earlier, is well modeled by anisotropic Compton scattering of the donor blackbody emission by relativistic electrons in a jet strongly misaligned with respect to the orbital axis. We confirm that this model fits well both the average γ-ray modulation light curve and the spectrum. However, we find that if the jet were aligned with the spin axis of a rotating black hole, it would undergo geodetic precession with a period of ∼50 yr. However, its presence is ruled out by both the γ-ray and radio data. Therefore, we consider an alternative model in which the average jet direction is aligned, but it is bent outside the orbit owing to the thrust of the donor stellar wind, and thus precesses at the orbital period. The γ-ray modulation then appears as a result of the variable Doppler boosting of synchrotron self-Compton jet emission. This model also fits the data well. However, the fitted bending angle is much larger than the theoretical one based on the binary and wind parameters as currently known. Thus, both models disagree with important aspects of our current theoretical understanding of the system. We discuss possible ways to find the correct model.

Hadronuclear Interactions in Active Galactic Nuclei Jets as the Origin of the Diffuse High-energy Neutrino Background

The origin of diffuse high-energy neutrinos from TeV to PeV energies detected by the IceCube Observatory remains a mystery. In our previous work, we have shown that hadronuclear (p − p) interactions in active galactic nuclei (AGNs) jets could be important and generate detectable very-high-energy emissions. Here, we further explore these interactions in the AGN jets based on their luminosity function. The diffuse neutrino flux and corresponding γ-ray flux have been calculated and compared with observational data. In our modeling, two beaming patterns are considered separately. To make sure that the corresponding γ-ray flux does not overshoot the diffuse γ-ray background, we find that if the neutrino production region in a jet is opaque to γ-rays, p − p interactions in AGN jets with a small viewing angle (the blazar case) are able to interpret the PeV neutrino background. Similarly, AGN jets with a large viewing angle (the radio galaxy case) may interpret the TeV neutrino background. While, if the neutrino production region is transparent to γ-rays, only blazars have the potential to interpret the diffuse neutrino background around the PeV band. Some caveats are also discussed.

Detection of Increasing γ-Ray Emissions from 4FGL J1718.5+4237 with Fermi-LAT

We report a gradual brightening γ-ray source, 4FGL J1718.5+4237, in 0.1–500.0 GeV, which may be associated with a blazar NVSS J171822+423948 with a redshift of ∼2.7. We analyzed 15.25 yr of γ-ray data recorded by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope and detected significant γ-ray emissions in the direction of the blazar with a test statistic (TS) of ∼135. Based on timing analysis using a 1 yr time bin, we have observed a gradual brightening in γ-ray emissions from the target. In our analysis, we categorize them into two states: Quiet (TS ∼ 0) and Loud (TS ∼ 226) states, with the distinction occurring in 2016 August (MJD 57602.69). From the Quiet state to the brightest period (the last data point), the γ-ray flux in 0.1–500.0 GeV increased by more than 12-fold from <0.2 × 10−8 photons cm−1 s−1 to 2.6 × 10−8 photons cm−1 s−1. Additionally, we studied the spectral properties in detail for the Loud state and the overall data. While no significant variation was detected, both exhibited a spectral index Γ of ∼2.6. The origin of the brightening γ-ray emissions from the target is not yet clear. Future long-term multi-wavelength observations and studies will provide insight into the astrophysical mechanisms of the target.

Improved constraints on dark matter annihilations around primordial black holes

Cosmology may give rise to appreciable populations of both particle dark matter and primordial black holes (PBH) with the combined mass density providing the observationally inferred value ΩDM ≈ 0.26. Early studies highlighted that scenarios with both particle dark matter and PBH are strongly excluded by γ-ray limits for particle dark matter with a velocity independent thermal cross section 〈σν〉 ~ 3 × 10−26cm3/s, as is the case for classic WIMP dark matter. Here we examine the limits from di use γ-rays on velocity-dependent, including annihilations which are p-wave with 〈σν〉 ∝ v2 or d-wave 〈σν〉 ∝ v4, which we find to be considerably less constraining. This work also utilizes a refined treatment of the PBH dark matter density profile. Importantly, we highlight that even if the freeze-out process is p-wave it is typical for (loop/phase-space) suppressed s-wave processes to actually provide the leading contributions to the experimentally constrained γ-ray flux from the PBH halo.

Constraining Possible γ-Ray Burst Emission from GW230529 Using Swift-BAT and Fermi-GBM

GW230529 is the first compact binary coalescence detected by the LIGO–Virgo–KAGRA collaboration with at least one component mass confidently in the lower mass gap, corresponding to the range 3–5 M ⊙. If interpreted as a neutron star–black hole merger, this event has the most symmetric mass ratio detected so far and therefore has a relatively high probability of producing electromagnetic (EM) emission. However, no EM counterpart has been reported. At the merger time t 0, Swift-BAT and Fermi-GBM together covered 100% of the sky. Performing a targeted search in a time window [t 0 − 20 s, t 0 + 20 s], we report no detection by the Swift-BAT and Fermi-GBM instruments. Combining the position-dependent γ-ray flux upper limits and the gravitational-wave posterior distribution of luminosity distance, sky localization, and inclination angle of the binary, we derive constraints on the characteristic luminosity and structure of the jet possibly launched during the merger. Assuming a top-hat jet structure, we exclude at 90% credibility the presence of a jet that has at the same time an on-axis isotropic luminosity ≳1048 erg s−1 in the bolometric band 1 keV–10 MeV and a jet opening angle ≳15°. Similar constraints are derived by testing other assumptions about the jet structure profile. Excluding GRB 170817A, the luminosity upper limits derived here are below the luminosity of any GRB observed so far.

Search for Neutrino Emission from the Cygnus Bubble Based on LHAASO γ-Ray Observations

The Cygnus region, which contains massive molecular and atomic clouds and young stars, is a promising Galactic neutrino source candidate. Cosmic-ray transport in the region can produce neutrinos and γ-rays. Recently, the Large High Altitude Air Shower Observatory (LHAASO) detected an ultrahigh-energy γ-ray bubble (Cygnus Bubble) in this region. Using publicly available track events detected by the IceCube Neutrino Observatory in 7 yr of full detector operation, we conduct searches for correlated neutrino signals from the Cygnus Bubble with neutrino emission templates based on LHAASO γ-ray observations. No significant signals were found for any employed templates. With the 7 TeV γ-ray flux template, we set a flux upper limit of 90% confidence level for the neutrino emission from the Cygnus Bubble to be 5.7 × 10−13 TeV−1 cm−2 s−1 at 5 TeV.

X-ray emission from large scale jets of AGNs at high redshifts

The limited angular resolution of the existing telescopes, operating in radio, optical and X-ray wavebands, hinder the study on multi-spectral component (MSC) emissions from the large scale jet of AGNs. This is particularly strenuous for the sources located at high redshifts. Though a few of them were resolved in radio-optical wavebands, their X-ray counterparts have been rarely discovered. This poses a bigger challenge in understanding their emission mechanism. At high redshifts, the X-ray emission from the large scale jets are generally interpreted as inverse-Compton scattering of cosmic microwave background photons (IC/CMB), due to the increase in CMB photon density. However, recent Fermi γ-ray flux upper limit estimates on two high redshift (z>3.5) AGN jets (J1510+5702 and J1421-0643), question the validity of the IC/CMB model. Here, we propose a model by considering the advection of electrons from the sites of particle acceleration (Rahman et al., 2022), to interpret the MSC emission from these sources. We show that observed radio/optical/X-ray emissions can be explained by the synchrotron emission from the advected and accelerated electron distribution. We also show that the IC/CMB spectrum of this composite electron distribution satisfies the constraints drawn from Fermi observations.

Limits on dark matter annihilation in prompt cusps from the isotropic gamma-ray background

Recent studies indicate that thermally produced dark matter will form highly concentrated, low-mass cusps in the early universe that often survive until the present. While these cusps contain a small fraction of the dark matter, their high density significantly increases the expected γ-ray flux from dark matter annihilation, particularly in searches of large angular regions. We utilize 14 years of Fermi-LAT data to set strong constraints on dark matter annihilation through a detailed study of the isotropic γ-ray background, excluding with 95% confidence dark matter annihilation to bb¯ final states for dark matter masses below 120 GeV. Published by the American Physical Society 2024

Gamma-ray variability and multi-wavelength insights into the unprecedented outburst from 4C 31.03

The blazar 4C31.03 recently underwent a major γ-ray outburst at the beginning of 2023 after a prolonged quiescent phase. Fermi-LAT reported a daily average flux of 5×10−6 phs cm−2 s−1, which is about 60 times its average value. We investigated this extraordinary outbreak through temporal and multi-wavelength analysis. From the statistical analysis of the γ-ray lightcurves using Bayesian blocks, we identified 3 epochs of prominent flares. The fastest flux decay during this major outburst was observed within 5.5±0.7 hours. The highest energy of γ-ray photons found from the source during the active phase is ∼82GeV. Using the transparency of γ-rays against pair production and light crossing time argument, we could obtain the minimum jet Doppler factor as 17 corresponding to the flaring state. The broadband spectral energy distribution study performed using synchrotron, SSC and EC emission processes supports the external Compton scattering of IR photons as the likely mechanism for the γ-ray emission from the source. The results of this study support the scenario of the emission region in 4C31.03, being located beyond the Broad line region from the central blackhole. The long-term γ-ray flux distribution depicts a double log-normal variability, indicating that two distinct flux states are active in this energy band. The index distribution also reveals a two distinct variability pattern and hints that the γ-ray spectrum can be more precisely described by two photon indices.

Fermi Large Area Telescope Detection of Gamma Rays from the NGC 6251 Radio Lobe

We report the detection of extended γ-ray emission from lobes in the radio galaxy NGC 6251 using observation data from the Fermi Large Area Telescope. The maximum likelihood analysis results show that a radio morphology template provides a better fit than a pointlike source description for the observational data at a confidence level of 8.1σ, and the contribution of lobes constitutes more than 50% of the total γ-ray flux. Furthermore, the γ-ray energy spectra show a significant disparity in shape between the core and lobe regions, with a curved log-parabola shape observed in the core region and a power-law form observed in the lobes. Neither the core region nor the northwest (NW) lobe displays significant flux variations in the long-term γ-ray light curves. The broadband spectral energy distributions of the core region and the NW lobe can be explained with a single-zone leptonic model. The γ-rays of the core region are due to the synchrotron-self-Compton process, while the γ-rays from the NW lobe are interpreted as inverse Compton emission of the cosmic microwave background.

Evidence of the γ-ray counterpart of nova FM Cir from Fermi–LAT

ABSTRACT We report the results of an analysis of X-ray and γ-ray data from the nova FM Cir taken by Swift and Fermi–LAT. The γ-ray emission from FM Cir can be identified with a significance level of ∼3σ within ∼40 d after the nova eruption (2018 January 19) when we bin the light curve per day. The significance can exceed the 4σ confidence level if we accumulate a longer time (i.e. 20 d) to bin the light curve. The γ-ray counterpart could be identified with a Test Statistic (TS) above 4 until ∼180 d after the eruption. The duration of the γ-ray detection is longer than those reported in previous studies of other novae detected in the GeV range. Significant X-ray emission was observed after the γ-ray flux level fell below the sensitivity limit of Fermi–LAT. The hardness ratio of the X-ray emission decreased rapidly with time, and the spectra were dominated by blackbody radiation from the hot white dwarf. Except for the longer duration of the γ-ray emission, the multiwavelength properties of FM Cir closely resemble those of other novae detected in the GeV range.

TeV flaring activity of the AGN PKS 0625–354 in November 2018

Most γ-ray detected active galactic nuclei are blazars with one of their relativistic jets pointing towards the Earth. Only a few objects belong to the class of radio galaxies or misaligned blazars. Here, we investigate the nature of the object PKS 0625−354, its γ-ray flux and spectral variability and its broad-band spectral emission with observations from H.E.S.S., Fermi-LAT, Swift-XRT, and UVOT taken in November 2018. The H.E.S.S. light curve above 200 GeV shows an outburst in the first night of observations followed by a declining flux with a halving time scale of 5.9 h. The γγ-opacity constrains the upper limit of the angle between the jet and the line of sight to ∼10°. The broad-band spectral energy distribution shows two humps and can be well fitted with a single-zone synchrotron self Compton emission model. We conclude that PKS 0625−354, as an object showing clear features of both blazars and radio galaxies, can be classified as an intermediate active galactic nuclei. Multi-wavelength studies of such intermediate objects exhibiting features of both blazars and radio galaxies are sparse but crucial for the understanding of the broad-band emission of γ-ray detected active galactic nuclei in general.

Searches for Neutrinos in the Direction of Radio-bright Blazars with the ANTARES Telescope

Active galaxies, especially blazars, are among the most promising extragalactic candidates for high-energy neutrino sources. To date, ANTARES searches included these objects and used GeV–TeV γ-ray flux to select blazars. Here, a statistically complete blazar sample selected by their bright radio emission is used as the target for searches of origins of neutrinos collected by the ANTARES neutrino telescope over 13 yr of operation. The hypothesis of a neutrino–blazar directional correlation is tested by pair counting and a complementary likelihood-based approach. The resulting posttrial p-value is 3.0% (2.2σ in the two-sided convention). Additionally, a time-dependent analysis is performed to search for temporal clustering of neutrino candidates as a means of detecting neutrino flares in blazars. None of the investigated sources alone reaches a significant flare detection level. However, the presence of 18 sources with a pretrial significance above 3σ indicates a p = 1.4% (2.5σ in the two-sided convention) detection of a time-variable neutrino flux. An a posteriori investigation reveals an intriguing temporal coincidence of neutrino, radio, and γ-ray flares of the J0242+1101 blazar at a p = 0.5% (2.9σ in the two-sided convention) level. Altogether, the results presented here suggest a possible connection of neutrino candidates detected by the ANTARES telescope with radio-bright blazars.

The Remarkable Predictive Power of Infrared Data of Blazars

Blazars are the brightest and most abundant persistent sources in the extragalactic γ-ray sky. Due to their significance, they are often observed across various energy bands, where the data of which can be used to explore potential correlations between emission at different energies, yielding valuable insights into the emission processes of their powerful jets. In this study we utilized IR data at 3.4 and 4.6 μm from the Near-Earth Object Wide-field Infrared Survey Explorer Reactivation Mission, spanning 8 yr of observations, X-ray data from the Neil Gehrels Swift Observatory collected throughout the satellite’s lifetime, and 12 years of γ-ray measurements from the Fermi Large Area Telescope’s all-sky survey. Our analysis reveals that the IR spectral slope reliably predicts the peak frequency and maximum intensity of the synchrotron component of blazar spectral energy distributions, provided it is uncontaminated by radiation unrelated to the jet. A notable correlation between the IR and γ-ray fluxes was observed, with the BL Lacertae subclass of blazars displaying a strong correlation coefficient of r = 0.80. IR band variability is more pronounced in flat spectrum radio quasars than in BL Lacertae objects, with mean fractional variability values of 0.65 and 0.35, respectively. We also observed that the synchrotron peak intensity of intermediate-high-energy-peaked objects can forecast their detectability at very high γ-ray energies. We used this predicting power to identify objects in current catalogs that could meet the detection threshold of the Cerenkov Telescope Array extragalactic survey, which should encompass approximately 180 blazars.

Broad-band spectral and temporal study of Ton 599 during the brightest 2023 January flare

ABSTRACT In this work, we provide a detailed analysis of the broad-band temporal and spectral properties of the blazar Ton 599 by using observations from the Fermi Large Area Telescope (LAT) and Swift X-Ray Telescope (XRT)/Ultraviolet–Optical Telescope (UVOT), during its brightest γ-ray flaring. The one-day bin γ-ray light curve exhibits multiple substructures with asymmetric and symmetric profiles. Notably, the γ-ray light curve shows a maximum flux of $\rm 3.63 \times 10^{-6}\, photon\, cm^{-2}\, s^{-1}$ on MJD 59954.50, which is the highest flux ever observed from this source. The correlation between the γ-ray flux and γ-ray spectral indices suggests a moderate ‘harder when brighter’ trend. Taking the γ-ray light curve as the reference, a strong correlation is observed with X-ray, optical, and UV energies. Additionally, the γ-rays and optical/UV emission exhibit higher variability compared with X-rays. To understand the parameter variation during the active state of the source, we conducted a statistical broad-band spectral modelling of the source in 10 flux intervals of equal duration. A one-zone leptonic model involving synchrotron, synchrotron-self-Compton, and external Compton processes successfully reproduces the broad-band spectral energy distribution (SED) in each of these flux intervals. We observed that flux variation during the active state is associated mainly with variation in the magnetic field and particle spectral indices.

Broad-band study of gamma-ray blazars at redshifts z = 2.0–2.5

ABSTRACT High redshift blazars are among the most powerful non-explosive sources in the Universe and play a crucial role in understanding the evolution of relativistic jets. To understand these bright objects, we performed a detailed investigation of the multiwavelength properties of 79 γ-ray blazars with redshifts ranging from z = 2.0 to 2.5, using data from Fermi LAT, Swift XRT/UVOT, and NuSTAR observations. In the γ-ray band, the spectral analysis revealed a wide range of flux and photon indices, from 5.32 × 10−10 to 3.40 × 10−7 photon cm−2 s−1 and from 1.66 to 3.15, respectively, highlighting the diverse nature of these sources. The detailed temporal analysis showed that flaring activities were observed in 31 sources. Sources such as 4C+71.07, PKS 1329-049, and 4C + 01.02, demonstrated significant increase in the γ-ray luminosity and flux variations, reaching peak luminosity exceeding 1050 erg s−1. The temporal analysis extended to X-ray and optical/ultraviolet (UV) bands, showed clear flux changes in some sources in different observations. The time-averaged properties of high redshift blazars were derived through modeling the spectral energy distributions with a one-zone leptonic scenario, assuming the emission region is within the broad-line region (BLR) and the X-ray and γ-ray emissions are due to inverse Compton scattering of synchrotron and BLR-reflected photons. This modeling allowed us to constrain the emitting particle distribution, estimate the magnetic field inside the jet, and evaluate the jet luminosity, which is discussed in comparison with the disc luminosity derived from fitting the excess in the UV band.