Leptonic Model for Neutrino Emission from Active Galactic Nuclei

It is often stated that the observation of high-energy neutrinos from an astrophysical source would constitute indisputable proof for the acceleration of hadronic cosmic rays. Here, we point out that there exists a purely leptonic mechanism to produce TeV-scale neutrinos in astrophysical environments. In particular, very high-energy synchrotron photons can scatter with x rays, exceeding the threshold for muon-antimuon pair production. When these muons decay, they produce neutrinos without any cosmic-ray protons or nuclei being involved. In order for this mechanism to be efficient, the source in question must produce very high-energy photons which interact in an environment that is dominated by keV-scale radiation. We find that such a source could potentially generate an observable neutrino flux through muon pair production for reasonable choices of physical parameters.

An extremely energetic cosmic ray observed by a surface detector array

Cosmic rays are energetic charged particles from extraterrestrial sources, with the highest-energy events thought to come from extragalactic sources. Their arrival is infrequent, so detection requires instruments with large collecting areas. In this work, we report the detection of an extremely energetic particle recorded by the surface detector array of the Telescope Array experiment. We calculate the particle's energy as [Formula: see text] (~40 joules). Its arrival direction points back to a void in the large-scale structure of the Universe. Possible explanations include a large deflection by the foreground magnetic field, an unidentified source in the local extragalactic neighborhood, or an incomplete knowledge of particle physics.

Observation of high-energy neutrinos from the Galactic plane

The origin of high-energy cosmic rays, atomic nuclei that continuously impact Earth's atmosphere, is unknown. Because of deflection by interstellar magnetic fields, cosmic rays produced within the Milky Way arrive at Earth from random directions. However, cosmic rays interact with matter near their sources and during propagation, which produces high-energy neutrinos. We searched for neutrino emission using machine learning techniques applied to 10 years of data from the IceCube Neutrino Observatory. By comparing diffuse emission models to a background-only hypothesis, we identified neutrino emission from the Galactic plane at the 4.5σ level of significance. The signal is consistent with diffuse emission of neutrinos from the Milky Way but could also arise from a population of unresolved point sources.

Higgsino Dark Matter Confronts 14 Years of Fermi γ-Ray Data

Thermal Higgsino dark matter (DM), with mass around 1 TeV, is a well-motivated, minimal DM scenario that arises in supersymmetric extensions of the standard model. Higgsinos may naturally be the lightest superpartners in split-supersymmetry models that decouple the scalar superpartners while keeping Higgsinos and gauginos close to the TeV scale. Higgsino DM may annihilate today to give continuum γ-ray emission at energies less than a TeV in addition to a linelike signature at energies equal to the mass. Previous searches for Higgsino DM, for example with the H.E.S.S. γ-ray telescope, have not reached the necessary sensitivity to probe the Higgsino annihilation cross section. In this work we make use of 14 years of data from the Fermi Large Area Telescope at energies above ∼10  GeV to search for the continuum emission near the Galactic Center from Higgsino annihilation. We interpret our results using DM profiles from Milky Way analog galaxies in the FIRE-2 hydrodynamic cosmological simulations. We set the strongest constraints to date on Higgsino-like DM. Our results show a mild, ∼2σ preference for Higgsino DM with a mass near the thermal Higgsino mass and, depending on the DM density profile, the expected cross section.

Weighing the Local Interstellar Medium Using Gamma Rays and Dust

Cold gas forms a significant mass fraction of the Milky Way disk, but is its most uncertain baryonic component. The density and distribution of cold gas is of critical importance for Milky Way dynamics, as well as models of stellar and galactic evolution. Previous studies have used correlations between gas and dust to obtain high-resolution measurements of cold gas, but with large normalization uncertainties. We present a novel approach that uses Fermi-LAT γ-ray data to measure the total gas density, achieving a similar precision as previous works, but with independent systematic uncertainties. Notably, our results have sufficient precision to probe the range of results obtained by current world-leading experiments.

Neutron star mass estimates from gamma-ray eclipses in spider millisecond pulsar binaries

Reliable neutron star mass measurements are key to determining the equation of state of cold nuclear matter, but such measurements are rare. Black widows and redbacks are compact binaries consisting of millisecond pulsars and semi-degenerate companion stars. Spectroscopy of the optically bright companions can determine their radial velocities, providing inclination-dependent pulsar mass estimates. Although inclinations can be inferred from subtle features in optical light curves, such estimates may be systematically biased due to incomplete heating models and poorly understood variability. Using data from the Fermi Large Area Telescope, we have searched for gamma-ray eclipses from 49 spider systems, discovering significant eclipses in 7 systems, including the prototypical black widow PSR B1957+20. Gamma-ray eclipses require direct occultation of the pulsar by the companion, and so the detection, or significant exclusion, of a gamma-ray eclipse strictly limits the binary inclination angle, providing new robust, model-independent pulsar mass constraints. For PSR B1957+20, the eclipse implies a much lighter pulsar (1.81 ± 0.07 solar masses) than inferred from optical light curve modelling.

Explaining the GeV Antiproton Excess, GeV γ-Ray Excess, and W-Boson Mass Anomaly in an Inert Two Higgs Doublet Model

For the newly discovered W-boson mass anomaly, one of the simplest dark matter (DM) models that can account for the anomaly without violating other astrophysical and experimental constraints is the inert two Higgs doublet model, in which the DM mass (m_{S}) is found to be within ∼54-74  GeV. In this model, the annihilation of DM via SS→bb[over ¯] and SS→WW^{*} would produce antiprotons and gamma rays, and may account for the excesses identified previously in both particles. Motivated by this, we reanalyze the AMS-02 antiproton and Fermi-LAT Galactic center γ-ray data. For the antiproton analysis, the novel treatment is the inclusion of the charge-sign-dependent three-dimensional solar modulation model as constrained by the time-dependent proton data. We find that the excess of antiprotons is more distinct than previous results based on the force-field solar modulation model. The interpretation of this excess as the annihilation of SS→WW^{*} (SS→bb[over ¯]) requires a DM mass of ∼40-80 (40-60) GeV and a velocity-averaged cross section of O(10^{-26})  cm^{3} s^{-1}. As for the γ-ray data analysis, besides adopting the widely used spatial template fitting, we employ an orthogonal approach with a data-driven spectral template analysis. The fitting to the GeV γ-ray excess yields DM model parameters overlapped with those to fit the antiproton excess via the WW^{*} channel. The consistency of the DM particle properties required to account for the W-boson mass anomaly, the GeV antiproton excess, and the GeV γ-ray excess suggests a common origin of them.

Gigaelectronvolt emission from a compact binary merger

An energetic γ-ray burst (GRB), GRB 211211A, was observed on 11 December 20211,2. Despite its long duration, typically associated with bursts produced by the collapse of massive stars, the observation of an optical-infrared kilonova points to a compact binary merger origin3. Here we report observations of a significant (more than five sigma) transient-like emission in the high-energy γ-rays of GRB 211211A (more than 0.1 gigaelectronvolts) starting 103 seconds after the burst. After an initial phase with a roughly constant flux (about 5 × 10-10 erg per second per square centimetre) lasting about 2 × 104 seconds, the flux started decreasing and soon went undetected. Our detailed modelling of public and dedicated multi-wavelength observations demonstrates that gigaelectronvolt emission from GRB 211211A is in excess with respect to the flux predicted by the state-of-the-art afterglow model at such late time. We explore the possibility that the gigaelectronvolt excess is inverse Compton emission owing to the interaction of a late-time, low-power jet with an external source of photons, and find that kilonova emission can provide the seed photons. Our results open perspectives for observing binary neutron star mergers.

A, Barbano A, Barwick SW, Bastian B, Basu V, Baur S, Bay R, Beatty JJ, Becker KH, Becker Tjus J, Bellenghi C, BenZvi S, Berley D, Bernardini E, Besson DZ, Binder G, Bindig D, Blaufuss E, Blot S, Boddenberg M, Bontempo F, Borowka J, Böser S, Botner O, Böttcher J, Bourbeau E, Bradascio F, Braun J, Brinson B, Bron S, Brostean-Kaiser J, Browne S, Burgman A, Burley RT, Busse RS, Campana MA, Carnie-Bronca EG, Chen C, Chen Z, Chirkin D, Choi K, Clark BA, Clark K, Classen L, Coleman A, Collin GH, Conrad JM, Coppin P, Correa P, Cowen DF, Cross R, Dappen C, Dave P, De Clercq C, DeLaunay JJ, Delgado López D, Dembinski H, Deoskar K, Desai A, Desiati P, de Vries KD, de Wasseige G, de With M, DeYoung T, Diaz A, Díaz-Vélez JC, Dittmer M, Dujmovic H, Dunkman M, DuVernois MA, Dvorak E, Ehrhardt T, Eller P, Engel R, Erpenbeck H, Evans J, Evenson PA, Fan KL, Fazely AR, Fedynitch A, Feigl N, Fiedlschuster S, Fienberg AT, Filimonov K, Finley C, Fischer L, Fox D, Franckowiak A, Friedman E, Fritz A, Fürst P, Gaisser TK, Gallagher J, Ganster E, Garcia A, Garrappa S, Gerhardt L, Ghadimi A, Glaser C, Glauch T, Glüsenkamp T, Goldschmidt A, Gonzalez JG, Goswami S, Grant D, Grégoire T, Griswold S, Günther C, Gutjahr P, Haack C, Hallgren A, Halliday R, Halve L, Halzen F, Ha Minh M, Hanson K, Hardin J, Harnisch AA, Haungs A, Hebecker D, Helbing K, Henningsen F, Hettinger EC, Hickford S, Hignight J, Hill C, Hill GC, Hoffman KD, Hoffmann R, Hokanson-Fasig B, Hoshina K, Huang F, Huber M, Huber T, Hultqvist K, Hünnefeld M, Hussain R, Hymon K, In S, Iovine N, Ishihara A, Jansson M, Japaridze GS, Jeong M, Jin M, Jones BJP, Kang D, Kang W, Kang X, Kappes A, Kappesser D, Kardum L, Karg T, Karl M, Karle A, Katz U, Kauer M, Kellermann M, Kelley JL, Kheirandish A, Kin K, Kintscher T, Kiryluk J, Klein SR, Koirala R, Kolanoski H, Kontrimas T, Köpke L, Kopper C, Kopper S, Koskinen DJ, Koundal P, Kovacevich M, Kowalski M, Kozynets T, Kun E, Kurahashi N, Lad N, Lagunas Gualda C, Lanfranchi JL, Larson MJ, Lauber F, Lazar JP, Lee JW, Leonard K, Leszczyńska A, Li Y, Lincetto M, Liu QR, Liubarska M, Lohfink E, Lozano Mariscal CJ, Lu L, Lucarelli F, Ludwig A, Luszczak W, Lyu Y, Ma WY, Madsen J, Mahn KBM, Makino Y, Mancina S, Mariş IC, Martinez-Soler I, Maruyama R, Mase K, McElroy T, McNally F, Mead JV, Meagher K, Mechbal S, Medina A, Meier M, Meighen-Berger S, Micallef J, Mockler D, Montaruli T, Moore RW, Morse R, Moulai M, Naab R, Nagai R, Nahnhauer R, Naumann U, Necker J, Nguyen LV, Niederhausen H, Nisa MU, Nowicki SC, Nygren D, Obertacke Pollmann A, Oehler M, Oeyen B, Olivas A, O’Sullivan E, Pandya H, Pankova DV, Park N, Parker GK, Paudel EN, Paul L, Pérez de los Heros C, Peters L, Peterson J, Philippen S, Pieper S, Pittermann M, Pizzuto A, Plum M, Popovych Y, Porcelli A, Prado Rodriguez M, Price PB, Pries B, Przybylski GT, Raab C, Rack-Helleis J, Raissi A, Rameez M, Rawlins K, Rea IC, Rehman A, Reichherzer P, Reimann R, Renzi G, Resconi E, Reusch S, Rhode W, Richman M, Riedel B, Roberts EJ, Robertson S, Roellinghoff G, Rongen M, Rott C, Ruhe T, Ryckbosch D, Rysewyk Cantu D, Safa I, Saffer J, Sanchez Herrera SE, Sandrock A, Sandroos J, Santander M, Sarkar S, Sarkar S, Satalecka K, Schaufel M, Schieler H, Schindler S, Schmidt T, Schneider A, Schneider J, Schröder FG, Schumacher L, Schwefer G, Sclafani S, Seckel D, Seunarine S, Sharma A, Shefali S, Silva M, Skrzypek B, Smithers B, Snihur R, Soedingrekso J, Soldin D, Spannfellner C, Spiczak GM, Spiering C, Stachurska J, Stamatikos M, Stanev T, Stein R, Stettner J, Steuer A, Stezelberger T, Stokstad R, Stürwald T, Stuttard T, Sullivan GW, Taboada I, Ter-Antonyan S, Tilav S, Tischbein F, Tollefson K, Tönnis C, Toscano S, Tosi D, Trettin A, Tselengidou M, Tung CF, Turcati A, Turcotte R, Turley CF, Twagirayezu JP, Ty B, Unland Elorrieta MA, Valtonen-Mattila N, Vandenbroucke J, van Eijndhoven N, Vannerom D, van Santen J, Verpoest S, Walck C, Watson TB, Weaver C, Weigel P, Weindl A, Weiss MJ, Weldert J, Wendt C, Werthebach J, Weyrauch M, Whitehorn N, Wiebusch CH, Williams DR, Wolf M, Woschnagg K, Wrede G, Wulff J, Xu XW, Yanez JP, Yoshida S, Yu S, Yuan T, Zhang Z, Zhelnin P. Evidence for neutrino emission from the nearby active galaxy NGC 1068

A supermassive black hole, obscured by cosmic dust, powers the nearby active galaxy NGC 1068. Neutrinos, which rarely interact with matter, could provide information on the galaxy’s active core. We searched for neutrino emission from astrophysical objects using data recorded with the IceCube neutrino detector between 2011 and 2020. The positions of 110 known gamma-ray sources were individually searched for neutrino detections above atmospheric and cosmic backgrounds. We found that NGC 1068 has an excess of 79 − 20 + 22 neutrinos at tera–electron volt energies, with a global significance of 4.2σ, which we interpret as associated with the active galaxy. The flux of high-energy neutrinos that we measured from NGC 1068 is more than an order of magnitude higher than the upper limit on emissions of tera–electron volt gamma rays from this source.

Rapid quasi-periodic oscillations in the relativistic jet of BL Lacertae

Blazars are active galactic nuclei (AGN) with relativistic jets whose non-thermal radiation is extremely variable on various timescales^1-3. This variability seems mostly random, although some quasi-periodic oscillations (QPOs), implying systematic processes, have been reported in blazars and other AGN. QPOs with timescales of days or hours are especially rare⁴ in AGN and their nature is highly debated, explained by emitting plasma moving helically inside the jet⁵, plasma instabilities^6,7 or orbital motion in an accretion disc^7,8. Here we report results of intense optical and γ-ray flux monitoring of BL Lacertae (BL Lac) during a dramatic outburst in 2020 (ref. ⁹). BL Lac, the prototype of a subclass of blazars¹⁰, is powered by a 1.7 × 10⁸ M_Sun (ref. ¹¹) black hole in an elliptical galaxy (distance = 313 megaparsecs (ref. ¹²)). Our observations show QPOs of optical flux and linear polarization, and γ-ray flux, with cycles as short as approximately 13 h during the highest state of the outburst. The QPO properties match the expectations of current-driven kink instabilities⁶ near a recollimation shock about 5 parsecs (pc) from the black hole in the wake of an apparent superluminal feature moving down the jet. Such a kink is apparent in a microwave Very Long Baseline Array (VLBA) image.

Evidence for PeV Proton Acceleration from Fermi-LAT Observations of SNR G106.3+2.7

The existence of a "knee" at energy ∼1  PeV in the cosmic-ray spectrum suggests the presence of Galactic PeV proton accelerators called "PeVatrons." Supernova remnant (SNR) G106.3+2.7 is a prime candidate for one of these. The recent detection of very high energy (0.1-100 TeV) gamma rays from G106.3+2.7 may be explained either by the decay of neutral pions or inverse Compton scattering by relativistic electrons. We report an analysis of 12 years of Fermi-LAT gamma-ray data that shows that the GeV-TeV gamma-ray spectrum is much harder and requires a different total electron energy than the radio and x-ray spectra, suggesting it has a distinct, hadronic origin. The nondetection of gamma rays below 10 GeV implies additional constraints on the relativistic electron spectrum. A hadronic interpretation of the observed gamma rays is strongly supported. This observation confirms the long-sought connection between Galactic PeVatrons and SNRs. Moreover, it suggests that G106.3+2.7 could be the brightest member of a new population of SNRs whose gamma-ray energy flux peaks at TeV energies. Such a population may contribute to the cosmic-ray knee and be revealed by future very high energy gamma-ray detectors.

Probing Spacetime Foam with Extragalactic Sources of High-Energy Photons

Quantum fluctuations can endow spacetime with a foamy structure. In this review article, we discuss our various proposals to observationally constrain models of spacetime foam. One way is to examine if the light wave-front from a distant quasar or GRB can be noticeably distorted by spacetime-foam-induced phase incoherence. As the phase fluctuations are proportional to the distance to the source but inversely proportional to the wavelength, ultra-high energy photons (>1 TeV) from distant sources are particularly useful. We elaborate on several proposals, including the possibility of detecting spacetime foam by observing “seeing disks” in the images of distant quasars and active galactic nuclei. We also discuss the appropriate distance measure for calculating the expected angular broadening. In addition, we discuss our more recent work in which we investigate whether wave-front distortions on small scales (due to spacetime foam) can cause distant objects become undetectable because the phase fluctuations have accumulated to the point at which image formation is impossible. Another possibility that has recently become accessible is to use interferometers to observe cosmologically distant sources, thereby giving a large baseline perpendicular to the local wave vector over which the wave front could become corrugated and thus distorted, reducing or eliminating its fringe visibility. We argue that all these methods ultimately depend on the availability of ways (if any) to carry out proper averaging of contributions from different light paths from the source to the telescope.

Antistars or Antimatter Cores in Mirror Neutron Stars?

The oscillation of the neutron n into mirror neutron n′, its partner from the dark mirror sector, can gradually transform an ordinary neutron star into a mixed star consisting in part of mirror dark matter. The implications of the reverse process taking place in the mirror neutron stars depend on the sign of baryon asymmetry in the mirror sector. Namely, if it is negative, as predicted by certain baryogenesis scenarios, then n′¯−n¯ transitions create a core of our antimatter gravitationally trapped in the mirror star interior. The annihilation of accreted gas on such antimatter cores could explain the origin of γ-source candidates with an unusual spectrum compatible with baryon–antibaryon annihilation, recently identified in the Fermi LAT catalog. In addition, some part of this antimatter escaping after the mergers of mirror neutron stars can produce the flux of cosmic antihelium and also heavier antinuclei which are hunted in the AMS-02 experiment.

A gamma-ray pulsar timing array constrains the nanohertz gravitational wave background

After large galaxies merge, their central supermassive black holes are expected to form binary systems. Their orbital motion should generate a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95% credible limit on the GWB characteristic strain of 1.0 × 10^-14 at a frequency of 1 yr^-1. The sensitivity is expected to scale with t _obs, the observing time span, as [Formula: see text]. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models.

Time-resolved hadronic particle acceleration in the recurrent nova RS Ophiuchi

Recurrent novae are repeating thermonuclear explosions in the outer layers of white dwarfs, due to the accretion of fresh material from a binary companion. The shock generated when ejected material slams into the companion star's wind can accelerate particles. We report very-high-energy (VHE, [Formula: see text]) gamma rays from the recurrent nova RS Ophiuchi, up to a month after its 2021 outburst, observed using the High Energy Stereoscopic System. The VHE emission has a similar temporal profile to lower-energy GeV emission, indicating a common origin, with a two-day delay in peak flux. These observations constrain models of time-dependent particle energization, favoring a hadronic emission scenario over the leptonic alternative. Shocks in dense winds provide favorable environments for efficient acceleration of cosmic-rays to very high energies.

Gamma-Ray Cosmology and Tests of Fundamental Physics

The propagation of gamma-rays over cosmological distances is the subject of extensive theoretical and observational research at GeV and TeV energies. The mean free path of gamma-rays in the cosmic web is limited above 100 GeV due to the production of electrons and positrons on the cosmic optical and infrared backgrounds. Electrons and positrons cool in the intergalactic medium while gyrating in its magnetic fields, which could cause either its global heating or the production of lower-energy secondary gamma-rays. The energy distribution of gamma-rays surviving the cosmological journey carries observed absorption features that gauge the emissivity of baryonic matter over cosmic time, constrain the distance scale of ΛCDM cosmology, and limit the alterations of the interaction cross section. Competitive constraints are, in particular, placed on the cosmic star-formation history as well as on phenomena expected from quantum gravity and string theory, such as the coupling to hypothetical axion-like particles or the violation of Lorentz invariance. Recent theoretical and observational advances offer a glimpse of the multi-wavelength and multi-messenger path that the new generation of gamma-ray observatories is about to open.

The Blazar Sequence and Its Physical Understanding

Introduced in 1998 to attempt a first unified view of the broad-band emission properties of blazars, the blazar sequence has been extensively used in the past 25 years to guide observations as well as the physical interpretation of the overall emission from these galaxies. In this review, we describe the evolution of the sequence along with the tremendous advances in the observational field, in particular in the gamma-ray band. A new version of the sequence built on TeV-detected objects is also presented. Two extreme classes of objects (MeV and hard-TeV blazars) are included in the discussion, given their relevance for future observatories. Finally, the current physical understanding at the base of the sequence is presented along with the major criticisms to the blazar sequence.

The Quasar CTD 135 Is Not a Compact Symmetric Object

The radio-loud quasar CTD 135 (2234+282, J2236+2828) has been proposed as a candidate compact symmetric object (CSO), based on its symmetric radio structure revealed by multi-frequency very long baseline interferometry (VLBI) imaging observations on milliarcsec angular scales. CSOs are known as young jetted active galactic nuclei (AGN) whose relativistic plasma jets are misaligned with respect to the line of sight. The peculiarity of CTD 135 as a CSO candidate was its detection in γ-rays, while the vast majority of known γ-ray emitting AGN are blazars with jets pointing close to our viewing direction. Since only a handful of CSOs are known as γ-ray sources, the unambiguous identification of a single candidate is important for studying this rare class of objects. By collecting and interpreting observational data from the recent literature, we revisit the classification of CTD 135. We present evidence that the object, based on its flat-spectrum radio core with high brightness temperature, variability at multiple wavebands, and infrared colours should be classified as a blazar rather than a CSO.

Nonthermal Radiation of the Extreme TeV Blazar 1ES 0229+200 from Electromagnetic Cascades on Infrared Photon Field

Extreme TeV blazars (ETBs) are active galactic nuclei with jets presumably pointing towards the observer having their intrinsic (compensated for the effect of γ-ray absorption on extragalactic background light photons) spectral energy distributions (SEDs) peaked at an energy in excess of 1 TeV. These sources typically reveal relatively weak and slow variability as well as higher frequency of the low-energy SED peak compared to other classes of blazars. It proved to be exceedingly hard to incorporate all these peculiar properties of ETBs into the framework of conventional γ-ray emission models. ETB physics have recently attracted great attention in the astrophysical community, underlying the importance of the development of self-consistent ETB emission model(s). We propose a new scenario for the formation of X-ray and γ-ray spectra of ETBs assuming that electromagnetic cascades develop in the infrared photon field surrounding the central blazar engine. This scenario does not invoke compact fast-moving sources of radiation (so-called “blobs”), in agreement with the apparent absence of fast and strong variability of ETBs. For the case of the extreme TeV blazar 1ES 0229+200 we propose a specific emission model in the framework of the considered scenario. We demonstrate that this model allows to obtain a good fit to the measured SED of 1ES 0229+200.

γ-ray Flux and Spectral Variability of Blazar Ton 599 during Its 2021 Flare

Blazars are known to emit exceptionally variable non-thermal emission over the wide range (from radio to γ-rays) of electromagnetic spectrum. We present here the results of our γ-ray flux and spectral variability study of the blazar Ton 599, which has been recently observed in the γ-ray flaring state. Using 0.1–300 GeV γ-ray data from the Fermi Gamma-ray Space Telescope (hereinafter Fermi), we generated one-day binned light curve of Ton 599 for a period of about one-year from MJD 59,093 to MJD 59,457. During this one year period, the maximum γ-ray flux detected was 2.24 ± 0.25 × 10−6 ph cm−2 s−1 at MJD 59,399.50. We identified three different flux states, namely, epoch A (quiescent), epoch B (pre-flare) and epoch C (main-flare). For each epoch, we calculated the γ-ray flux variability amplitude (Fvar) and found that the source showed largest flux variations in epoch C with Fvar∼ 35%. We modelled the γ-ray spectra for each epoch and found that the Log-parabola model adequately describes the γ-ray spectra for all the three epochs. We estimated the size of the γ-ray emitting region as 1.03 × 1016 cm and determined that the origin of γ-ray radiation, during the main-flare, could be outside of the broad line region.

Astrophysical Neutrinos and Blazars

We review and discuss recent results on the search for correlations between astrophysical neutrinos and γ-ray-detected sources, with many extragalactic studies reporting potential associations with different types of blazars. We investigate possible dependencies on blazar sub-classes by using the largest catalogues and all the multi-frequency data available. Through the study of similarities and differences in these sources we conclude that blazars come in two distinct flavours: LBLs and IHBLs (low-energy-peaked and intermediate-high-energy-peaked objects). These are distinguished by widely different properties such as the overall spectral energy distribution shape, jet speed, cosmological evolution, broad-band spectral variability, and optical polarisation properties. Although blazars of all types have been proposed as neutrino sources, evidence is accumulating in favour of IHBLs being the counterparts of astrophysical neutrinos. If this is indeed the case, we argue that the peculiar observational properties of IHBLs may be indirectly related to proton acceleration to very high energies.

Fermi-LAT Observation of PSR B1259-63 during Its 2021 Periastron Passage

PSR B1259-63 is a γ-ray binary system, where the compact object is a pulsar. The system has an orbital period of 1236.7 days and shows peculiar γ-ray flares (in 100 MeV–300 GeV) after its periastron time. We analyzed the Fermi-LAT observation of PSR B1259-63 during its latest periastron passage, as well as its previous three periastrons. The bright GeV flares started about 60 days after the periastron epoch in 2021. This delay is larger than that around the 2017 periastron and much larger than earlier periastrons. The delay of the GeV flux peak time in each periastron passage is apparent in our results. We discussed the possible origin of this delay and made a prediction of the GeV flux peak time in next periastron passage, based on observation of the previous delays.

A GeV-TeV particle component and the barrier of cosmic-ray sea in the Central Molecular Zone

Cosmic rays are important probe of a number of fundamental physical problems such as the acceleration of high and very high energy particles in extreme astrophysical environments. The Galactic center is widely anticipated to be an important cosmic-ray source and the observations of some Imaging Atmospheric Cherenkov Telescopes did successfully reveal a component of TeV-PeV cosmic rays in the vicinity of the Galactic center. Here we report the identification of GeV-TeV cosmic rays in the central molecular zone with the γ-ray observations of the Fermi Large Area Telescope, whose spectrum and spatial gradient are consistent with that measured by the Imaging Atmospheric Cherenkov Telescopes but the corresponding cosmic-ray energy density is substantially lower than the so-called cosmic-ray sea component, suggesting the presence of a high energy particle accelerator at the Galactic center and the existence of a barrier that can effectively suppress the penetration of the particles from the cosmic-ray sea to the central molecular zone.

Dissecting the Inner Galaxy with γ-Ray Pixel Count Statistics

We combine adaptive template fitting and pixel count statistics in order to assess the nature of the Galactic Center excess in Fermi-LAT data. We reconstruct the flux distribution of point sources well below the Fermi-LAT detection threshold, and measure their radial and longitudinal profiles in the inner Galaxy. We find that all point sources and the bulge-correlated diffuse emission each contributes O(10%) of the total inner Galaxy emission, and disclose a potential subthreshold point-source contribution to the Galactic Center excess.

A New Sample of Gamma-Ray Emitting Jetted Active Galactic Nuclei—Preliminary Results

We are compiling a new list of gamma-ray jetted active galactic nuclei (AGN), starting from the fourth catalog of point sources of the Fermi Large Area Telescope (LAT). Our aim is to prepare a list of jetted AGN with known redshifts and classifications to be used to calibrate jet power. We searched in the available literature for all the published optical spectra and multiwavelength studies useful to characterize the sources. We found new, missed, or even forgotten information leading to a substantial change in the redshift values and classification of many sources. We present here the preliminary results of this analysis and some statistics based on the gamma-ray sources with right ascension within the interval 0h--12h (J2000). Although flat-spectrum radio quasars and BL Lac objects are still the dominant populations, there is a significant increase in the number of other objects, such as misaligned AGN, narrow-line Seyfert 1 galaxies, and Seyfert galaxies. We also introduced two new classes of objects: changing-look AGN and ambiguous sources. About one third of the sources remain unclassified.

Advanced Virgo: Status of the Detector, Latest Results and Future Prospects

The Virgo detector, based at the EGO (European Gravitational Observatory) and located in Cascina (Pisa), played a significant role in the development of the gravitational-wave astronomy. From its first scientific run in 2007, the Virgo detector has constantly been upgraded over the years; since 2017, with the Advanced Virgo project, the detector reached a high sensitivity that allowed the detection of several classes of sources and to investigate new physics. This work reports the main hardware upgrades of the detector and the main astrophysical results from the latest five years; future prospects for the Virgo detector are also presented.

Signatures of Recent Cosmic-Ray Acceleration in the High-latitude Gamma-Ray Sky

Cosmic-ray (CR) sources temporarily enhance the relativistic particle density in their vicinity over the background distribution accumulated from the Galaxy-wide past injection activity and propagation. If individual sources are close enough to the solar system, their localized enhancements may present as features in the measured spectra of the CRs and in the associated secondary electromagnetic emissions. Large-scale loop-like structures visible in the radio sky are possible signatures of such nearby CR sources. If so, these loops may also have counterparts in the high-latitude γ-ray sky. Using ~10 yr of data from the Fermi Large Area Telescope, applying Bayesian analysis including Gaussian Processes, we search for extended enhanced emission associated with putative nearby CR sources in the energy range from 1 GeV to 1 TeV for the sky region |b| > 30°. We carefully control the systematic uncertainty due to imperfect knowledge of the interstellar gas distribution. Radio Loop IV is identified for the first time as a γ-ray emitter, and we also find significant emission from Loop I. Strong evidence is found for asymmetric features about the Galactic l = 0° meridian that may be associated with parts of the so-called "Fermi Bubbles," and some evidence is also found for γ-ray emission from other radio loops. Implications for the CRs producing the features and possible locations of the sources of the emissions are discussed.

Gamma-Ray Observation of the Cygnus Region in the 100-TeV Energy Region

We report observations of gamma-ray emissions with energies in the 100-TeV energy region from the Cygnus region in our Galaxy. Two sources are significantly detected in the directions of the Cygnus OB1 and OB2 associations. Based on their positional coincidences, we associate one with a pulsar PSR J2032+4127 and the other mainly with a pulsar wind nebula PWN G75.2+0.1, with the pulsar moving away from its original birthplace situated around the centroid of the observed gamma-ray emission. This work would stimulate further studies of particle acceleration mechanisms at these gamma-ray sources.

Multi-Wavelength Properties of the 2021 Periastron Passage of PSR B1259-63

PSR B1259-63 is a gamma-ray binary system hosting a radio pulsar orbiting around a O9.5Ve star, LS 2883, with a period of ∼3.4 years. The interaction of the pulsar wind with the LS 2883 outflow leads to unpulsed broadband emission in the radio, X-ray, GeV, and TeV domains. One of the most unusual features of the system is an outburst of GeV energies around the periastron, during which the energy release substantially exceeds the spin down luminosity under the assumption of the isotropic emission. In this paper, we present the first results of a recent multi-wavelength campaign (radio, optical, and X-ray bands) accompanied by the analysis of publicly available GeV Fermi/LAT data. The campaign covered a period of more than 100 days around the 2021 periastron and revealed substantial differences from previously observed passages. We report a major delay of the GeV flare, weaker X-ray flux during the peaks, which are typically attributed to the times when the pulsar crosses the disk, and the appearance of a third X-ray peak never observed before. We argue that these features are consistent with the emission cone model proposed by us previously, in the case of a sparser and clumpier disk of the Be star.

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622+3749 Observed by LHAASO-KM2A

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High-Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of 8.2σ for E>25  TeV assuming a Gaussian template. The best-fit location is (right ascension, declination) =(95.47°±0.11°,37.92°±0.09°), and the extension is 0.40°±0.07°. The energy spectrum can be described by a power-law spectrum with an index of -2.92±0.17_{stat}±0.02_{sys}. No clear extended multiwavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of ∼160  TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

Ultrahigh-energy photons up to 1.4 petaelectronvolts from 12 γ-ray Galactic sources

The extension of the cosmic-ray spectrum beyond 1 petaelectronvolt (PeV; 10¹⁵ electronvolts) indicates the existence of the so-called PeVatrons-cosmic-ray factories that accelerate particles to PeV energies. We need to locate and identify such objects to find the origin of Galactic cosmic rays¹. The principal signature of both electron and proton PeVatrons is ultrahigh-energy (exceeding 100 TeV) γ radiation. Evidence of the presence of a proton PeVatron has been found in the Galactic Centre, according to the detection of a hard-spectrum radiation extending to 0.04 PeV (ref. ²). Although γ-rays with energies slightly higher than 0.1 PeV have been reported from a few objects in the Galactic plane^3-6, unbiased identification and in-depth exploration of PeVatrons requires detection of γ-rays with energies well above 0.1 PeV. Here we report the detection of more than 530 photons at energies above 100 teraelectronvolts and up to 1.4 PeV from 12 ultrahigh-energy γ-ray sources with a statistical significance greater than seven standard deviations. Despite having several potential counterparts in their proximity, including pulsar wind nebulae, supernova remnants and star-forming regions, the PeVatrons responsible for the ultrahigh-energy γ-rays have not yet been firmly localized and identified (except for the Crab Nebula), leaving open the origin of these extreme accelerators.

Gamma Rays as Probes of Cosmic-Ray Propagation and Interactions in Galaxies

Continuum gamma-ray emission produced by interactions of cosmic rays with interstellar matter and radiation fields is a probe of non-thermal particle populations in galaxies. After decades of continuous improvements in experimental techniques and an ever-increasing sky and energy coverage, gamma-ray observations reveal in unprecedented detail the properties of galactic cosmic rays. A variety of scales and environments are now accessible to us, from the local interstellar medium near the Sun and the vicinity of cosmic-ray accelerators, out to the Milky Way at large and beyond, with a growing number of gamma-ray emitting star-forming galaxies. Gamma-ray observations have been pushing forward our understanding of the life cycle of cosmic rays in galaxies and, combined with advances in related domains, they have been challenging standard assumptions in the field and have spurred new developments in modelling approaches and data analysis methods. We provide a review of the status of the subject and discuss perspectives on future progress.

20 Years of Indian Gamma Ray Astronomy Using Imaging Cherenkov Telescopes and Road Ahead

The field of ground-based γ-ray astronomy has made very significant advances over the last three decades with the extremely successful operations of several atmospheric Cherenkov telescopes worldwide. The advent of the imaging Cherenkov technique for indirect detection of cosmic γ rays has immensely contributed to this field with the discovery of more than 220 γ-ray sources in the Universe. This has greatly improved our understanding of the various astrophysical processes involved in the non-thermal emission at energies above 100 GeV. In this paper, we summarize the important results achieved by the Indian γ-ray astronomers from the GeV-TeV observations using imaging Cherenkov telescopes over the last two decades. We mainly emphasize the results obtained from the observations of active galactic nuclei with the TACTIC (TeV Atmospheric Cherenkov Telescope with Imaging Camera) telescope, which has been operational since 1997 at Mount Abu, India. We also discuss the future plans of the Indian γ-ray astronomy program with special focus on the scientific objectives of the recently installed 21 m diameter MACE (Major Atmospheric Cherenkov Experiment) telescope at Hanle, India.

Galactic Center Excess in a New Light: Disentangling the γ-Ray Sky with Bayesian Graph Convolutional Neural Networks

A fundamental question regarding the Galactic Center excess (GCE) is whether the underlying structure is pointlike or smooth, often framed in terms of a millisecond pulsar or annihilating dark matter (DM) origin for the emission. We show that Bayesian neural networks (NNs) have the potential to resolve this debate. In simulated data, the method is able to predict the flux fractions from inner Galaxy emission components to on average ∼0.5%. When applied to the Fermi photon-count map, the NN identifies a smooth GCE in the data, suggestive of the presence of DM, with the estimates for the background templates being consistent with existing results.

Inference of the Local Interstellar Spectra of Cosmic-Ray Nuclei Z ⩽ 28 with the GalProp-HelMod Framework

Composition and spectra of Galactic cosmic rays (CRs) are vital for studies of high-energy processes in a variety of environments and on different scales, for interpretation of γ-ray and microwave observations, for disentangling possible signatures of new phenomena, and for understanding of our local Galactic neighborhood. Since its launch, AMS-02 has delivered outstanding-quality measurements of the spectra of p ¯ , e ±, and nuclei: 1H-8O, 10Ne, 12Mg, 14Si. These measurements resulted in a number of breakthroughs; however, spectra of heavier nuclei and especially low-abundance nuclei are not expected until later in the mission. Meanwhile, a comparison of published AMS-02 results with earlier data from HEAO-3-C2 indicates that HEAO-3-C2 data may be affected by undocumented systematic errors. Utilizing such data to compensate for the lack of AMS-02 measurements could result in significant errors. In this paper we show that a fraction of HEAO-3-C2 data match available AMS-02 measurements quite well and can be used together with Voyager 1 and ACE-CRIS data to make predictions for the local interstellar spectra (LIS) of nuclei that are not yet released by AMS-02. We are also updating our already-published LIS to provide a complete set from 1H-28Ni in the energy range from 1 MeV nucleon-1 to ~100-500 TeV nucleon-1, thus covering 8-9 orders of magnitude in energy. Our calculations employ the GalProp-HelMod framework, which has proved to be a reliable tool in deriving the LIS of CR p ¯ , e -, and nuclei 1H-8O.

Spurious Point Source Signals in the Galactic Center Excess

We reexamine evidence that the Galactic Center Excess (GCE) originates primarily from point sources (PSs). We show that in our region of interest, non-Poissonian template fitting evidence for GCE PSs is an artifact of unmodeled north-south asymmetry of the GCE. This asymmetry is strongly favored by the fit (although it is unclear if this is physical), and when it is allowed, the preference for PSs becomes insignificant. We reproduce this behavior in simulations, including detailed properties of the spurious PS population. We conclude that the non-Poissonian template fitting evidence for GCE PSs is highly susceptible to certain systematic errors and should not at present be taken to robustly disfavor a dominantly smooth GCE.

Very Long Baseline Interferometry Observations of the Proposed Radio Counterpart of an EGRET Source

We present high-resolution radio interferometric imaging observations of the radio source NVSS J182659+343113 (hereafter J1826+3431), the proposed radio counterpart of the γ-ray source, 3EG J1824+3441 detected by the Energetic Gamma Ray Experiment Telescope (EGRET) on board the Compton Gamma Ray Observatory satellite. We analyzed eight epochs of archival multi-frequency very long baseline interferometry data. We imaged the asymmetric core–jet structure of the source, and detected apparent superluminal motion in the jet. At the highest observing frequency, 15.3 GHz, the core shows high brightness temperature indicating Doppler boosting. Additionally, the radio features undergo substantial flux density variability. These findings strengthen the previous claim of the association of the blazar J1826+3431 with the possible γ-ray source, 3EG J1824+3441.

What If Planet 9 Is a Primordial Black Hole?

We highlight that the anomalous orbits of trans-Neptunian objects (TNOs) and an excess in microlensing events in the 5-year Optical Gravitational Lensing Experiment data set can be simultaneously explained by a new population of astrophysical bodies with mass several times that of the Earth (M_{⊕}). We take these objects to be primordial black holes (PBHs) and point out the orbits of TNOs would be altered if one of these PBHs was captured by the Solar System, inline with the Planet 9 hypothesis. Capture of a free floating planet is a leading explanation for the origin of Planet 9, and we show that the probability of capturing a PBH instead is comparable. The observational constraints on a PBH in the outer Solar System significantly differ from the case of a new ninth planet. This scenario could be confirmed through annihilation signals from the dark matter microhalo around the PBH.

Testing the Sensitivity of the Galactic Center Excess to the Point Source Mask

The Fermi Large Area Telescope (Fermi-LAT) Collaboration has an updated point source catalog, referred to as 4FGL. We perform the first template fit using a mask based on this new catalog and find that the excess in gamma rays detected at the Galactic Center in Fermi-LAT data persists. On the other hand, we find that a search for point sources is highly sensitive to the use of the 4FGL catalog: no sizable excess of bright pixels is apparent in the inner Galaxy when we mask out 4FGL point sources. Combining these observations restricts the ability of point sources to contribute to the Galactic Center excess. After identifying which bright sources have no known counterpart, we place strong constraints on any point source luminosity function capable of explaining the smooth emission identified in the template fit.