Simone Maldera

We present a catalog of sources detected above 50 GeV by the Fermi Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope using 16 yr of observations. This is the Fourth Catalog of Hard Fermi-LAT Sources (4FHL), containing 673 objects detected in the 50 GeV–2 TeV energy range. The sensitivity and localization are improved by factors of ∼2 relative to the previous LAT catalog at the same energies (i.e. the 2FHL catalog). The increased photon statistics, nearly doubling the number of detections relative to 2FHL, enables more detailed population studies, including the identification of fainter sources detected above 50 GeV. The majority of the sources (84.7%) are associated with extragalactic counterparts, primarily BL Lacertae objects, including 36 objects at redshift z > 1. Galactic sources account for 10% of the sample, mainly pulsar wind nebulae and supernova remnants, while ∼5% remain unassociated or of uncertain classification. Twenty sources are reported with no counterpart in previous LAT catalogs. The 4FHL catalog provides the deepest all-sky survey above 50 GeV to date. By focusing on this high-energy regime, it provides a fundamental link between space-based gamma-ray observations and the TeV range probed by current and future ground-based observatories. It serves as a useful reference for guiding future observations, including new potential targets for follow-up with Cherenkov telescopes.

We used 14 yr of Fermi-Large Area Telescope (LAT) data to derive a catalog of sources detected between 20 MeV and 200 MeV. For each source, we provide localization, spectral, and association information. The Second Catalog of Low-Energy Fermi-LAT sources (2FLE) contains 450 sources detected at TS > 25, with median energy flux of F _{20−200 MeV} = 1.9 × 10 ⁻¹¹ erg cm ⁻² s ⁻¹ . Using a Bayesian and likelihood-ratio method, we associate 377 2FLE sources with a known counterpart. Of the 325 objects associated with extragalactic sources, the dominant source class represented in the 2FLE is flat-spectrum radio quasars (FSRQs). Additionally, 52 sources are associated with Galactic sources. Three associated sources have not been previously detected in any γ -ray catalog and represent new γ -ray sources of known origin. We find that 387 2FLE sources have a positional cross-match with the 4FGL-DR4 sources. For the sources with a 4FGL counterpart, we compare the spectrum of the 2FLE source to the 4FGL spectral energy distributions. Of the 73 2FLE sources without multiwavelength associations, 13 have a 4FGL counterpart that lacks classification, while 60 are newly discovered γ -ray sources. Among these 60 sources, 51 are found at high latitudes (∣ b ∣ > 10°), and 38 are detected at >10 σ . Additionally, we find that out of the 450 2FLE sources, 45 are TeV emitters, and six are spatially coincident with a 1LHAASO source. The 2FLE helps us pave the way for reducing the “MeV gap,” and provides the most sensitive catalog between 20 and 200 MeV to date.

High-energy astrophysical neutrinos serve as crucial messengers for understanding hadronic acceleration processes and identifying the origins of cosmic rays, with blazars among the most promising neutrino sources. The KM3NeT experiment reported the detection of an ultrahigh-energy neutrino with an energy estimate of ∼220 PeV, the most energetic yet observed. The neutrino arrival direction has a 99% confidence region of 3° radius centered at R.A. 94 . ° 3, decl. −7 . ° 8 (J2000). In this work, 17 candidate blazars located within this region are identified. Comprehensive new observations and archival data analysis for these sources are presented. The study provides a complete multiwavelength coverage across radio, optical, X-ray, and gamma-ray bands, including proprietary data and dedicated follow-up observations. This study highlights flaring behavior in several candidate counterparts. One object exhibits a radio flare coinciding with the neutrino arrival time, with a pre-trial chance probability of 0.26%. Only for this object, the chance probability could be quantified in a robust way. Another candidate displays a rising trend in X-ray flux in a 1 yr window around the neutrino arrival time, while a third undergoes a gamma-ray flare during the same period. Based on the observational findings presented here, while none of these candidates can conclusively be linked to the neutrino, the implications of a possible blazar origin for the KM3NeT event are discussed.

We investigate the nature of the unassociated sources detected by the Fermi-LAT close (∣ b ∣ < 10°) to the Galactic plane, representing 16% of all sources in the 4FGL-DR4 catalog. The bulk of these sources (referred to as soft Galactic unassociated sources, SGUs) exhibit properties not found in known classes of gamma-ray emitters, as confirmed by a machine-learning classification approach. In particular, these properties include a steep, curved spectrum peaking below 1 GeV and a specific Galactic-latitude distribution with both a narrow and a broad component (dubbed the spike and the shoulder, respectively). Some source clusters are highlighted. New plausible source classes are explored, but only star-forming regions are found to account for a significant fraction (at most 10%) of the unassociated population. A thorough search for counterparts to the 175 brightest sources brings out a number of plausible counterparts but does not reveal clues about the nature of the whole population. We investigate the possibility that SGUs originate from mismodeled clumps of diffuse emission. Using Monte Carlo simulations, the SGU spectra can be reproduced in this scenario under an ad hoc condition concerning the clump spatial extension. The possible connection between the SGUs and gas not accounted for by the ¹² CO tracer is explored using the ¹³ CO MOPRA data but leads to inconclusive results. The origin of SGUs being related to diffuse emission remains plausible. However, a scenario whereby SGUs represent a new class of gamma-ray emitters cannot be fully excluded.

We report on the average and orbital phase-resolved polarization of Cyg X-3 in the hard state (HS) during the 2023 Imaging X-ray Polarimetry Explorer (IXPE) observational campaign. We find the polarization degree of 21.4% ± 0.4% and polarization angle of 92 . ° 2 ± 0 . ° 5 , compatible within 3 σ with the first HS IXPE observation in 2022. As the observed polarization depends on both the accretion geometry and the X-ray emission mechanism, and assuming it arises from reflection from the optically thick envelope surrounding the central source, our result indicates that both are similar in the two HS epochs on year-long timescales despite the source transitioning through an intermediate state in between. We discuss time- and energy-dependent polarization properties and their implications for the geometry and stability of the accretion funnel, finding significant orbital modulation of the polarization properties.

Context. Superluminous supernovae (SLSNe) are a rare class of transients with peak luminosities 10–100 times greater than those of standard core-collapse supernovae (SNe). The mechanisms powering their extreme brightness remain debated, with circumstellar medium (CSM) interaction, or energy injection from a central engine like a magnetar wind nebula being the most plausible scenarios. While the optical properties of SLSNe are extensively studied, their γ -ray signatures remain poorly constrained. Aims. To further constrain the underlying mechanism, we carried out a systematic search for giga-electronvolt γ -ray emission using the Fermi Large Area Telescope (LAT) from a sample of nearby hydrogen-poor (Type I) and hydrogen-rich (Type II) SLSNe over the past 16 years. Our objective is to test predictions from CSM and magnetar models, and to assess the prospects for future detections with the Cherenkov Telescope Array Observatory (CTAO). Methods. For the six targets of this sample, we studied the time variability of a putative γ -ray signal at the optical position of the SLSN on a six-month timescale, and in the case of SN 2017egm, we further investigated variability on 15-day intervals and applied a Bayesian block algorithm to characterize the time variability of the signal. We then compared the temporal evolution and spectral properties to the predictions from a magnetar and CSM interaction model. Results. Among the sample, only SN 2017egm shows significant γ -ray emission, with likelihood test statistic (TS) values of 26–33 (i.e., > 5 σ ) depending on the adopted time window. The signal arises between 50 and 160 days after explosion and is well described by a power-law spectrum with index Γ = 2.17 ± 0.23. The emission is consistent both in terms of its light curve and its spectrum, with predictions from magnetar models requiring either low nebular magnetization or faster spin-down than dipole losses. The CSM shell interaction scenario can reproduce the observed flux level but not the observed timing of the γ -ray signal. In addition, the observed ratio, L _γ / L _opt ∼ 1, is inconsistent with theoretical expectations and not in line with ratio measurements in other interacting CSM-dominated objects (e.g., novae or SNe) where this ratio is less than 10 ⁻² . Conclusions. Our study strongly suggests that a central engine like a magnetar plays a key role in this SLSN and could explain the bulk of the optical and γ -ray light curves properties. In order to explain the observed late-time bumps in the optical light curve of SN 2017egm, we require either: a hybrid picture combining magnetar and multiple CSM shells for the optical bumps or a pure magnetar model with infalling matter on an accretion disk. Finally, simulations of 50 hours of CTAO observations indicate that a SN 2017egm-like event would be detectable up to ∼140 Mpc in the magnetar model but not in the CSM model due to strong γ − γ absorption.

The extended X-ray emission observed in the direction of several molecular clouds in the central molecular zone (CMZ) of our Galaxy exhibits spectral and temporal properties consistent with the X-ray echo scenario. This concept postulates that the observed signal is a light-travel-time delayed reflection of a short ( δt < 1.5 yr) and bright ( L _X > 10 ³⁹ erg s ⁻¹ ) flare that was most probably produced a few hundred years ago by Sgr A*. This scenario also predicts a distinct polarization signature for the reflected X-ray continuum, with the polarization vector being perpendicular to the direction toward the primary source and the polarization degree being determined by the scattering angle. We report the results of two deep observations of the currently brightest (in reflected emission) molecular complex Sgr A taken with the Imaging X-ray Polarimetry Explorer (IXPE) in 2022 and 2023. We confirm the previous polarization measurement for a large region encompassing the Sgr A complex with high significance. We also reveal an inconsistent polarization pattern for the brightest reflection region in its center. Specifically, the X-ray polarization from this region is almost perpendicular to the expected direction in the case of Sgr A* illumination, and it shows a smaller degree of polarization compared to the large region. Taken at face value, this could indicate the simultaneous propagation of several illumination fronts throughout the CMZ, with the origin of one of them not being Sgr A*. The primary source could be associated with the Arches stellar cluster or a currently unknown source located closer to the illuminated cloud, potentially lowering the required luminosity of the primary source. Although significantly deeper observations with IXPE would be required to unequivocally distinguish between the scenarios, a combination of high-resolution imaging and micro-calorimetric spectroscopy offers an additional promising path forward.

RCW 86 is a supernova remnant whose origin has recently been linked to an off-center explosion within a cavity created by its progenitor star. In the southwestern region, the forward shock is thought to have reached the cavity wall, encountering diverse environmental conditions. We report on the spatially resolved X-ray polarimetric observation of RCW 86 with the Imaging X-ray Polarimetry Explorer (IXPE). In the 2–4.5 keV energy band we find no significant detection of polarization. Employing a dedicated background subtraction procedure and Bayesian spectropolarimetric fitting, we derive 99% upper limits on the polarization degree of the synchrotron component: 15% in higher-statistics regions and 30%–40% in lower-statistics regions. These upper limits on the polarization degree in several regions exclude the possibility of a strongly coherent magnetic field down to the subparsec scale, and that of a moderately coherent one on the scale of the synchrotron features as resolved by IXPE. The results indicate that the shocks in the southwestern rim of RCW 86 propagate more slowly than the unshocked ejecta at their locations, yet exceed the measured proper motion speeds. This behavior is consistent with reflected shocks occurring in tenuous regions of the shocked ejecta, distinct from regions that are radio-bright.

We present the X-ray polarization observation of G21.5−0.9, a young Galactic supernova remnant (SNR), conducted with the Imaging X-ray Polarimetry Explorer (IXPE) in 2023 October, with a total livetime of approximately 837 ks. Using different analysis methods, such as a space-integrated study of the entire region of the pulsar wind nebula (PWN) and a space-resolved polarization map, we detect significant polarization from the PWN at the center of the SNR, with an average polarization degree of ∼10% oriented at ∼33° (north through east). No significant energy-dependent variation in polarization is observed across the IXPE band (2–8 keV). The polarization map, corrected for the effect of polarization leakage, reveals a consistent pattern in both degree and angle, with little change across the nebula. Our findings indicate the presence of a highly polarized central torus, suggesting low levels of turbulence at particle acceleration sites. Unlike Vela, but similar to the Crab Nebula, we observe substantial differences between radio and X-ray polarization maps. This suggests a clear separation in energy of the emitting particle populations and hints at an important, yet poorly understood, role of instabilities in the turbulence dynamics of PWNe.

The steady-state gamma-ray emission from the Sun is thought to consist of two emission components due to interactions with Galactic cosmic rays: (1) a hadronic disk component, and (2) a leptonic extended component peaking at the solar edge and extending into the heliosphere. The flux of these components is expected to vary with the 11 yr solar cycle, being highest during solar minimum and lowest during solar maximum, as it varies with the cosmic-ray flux. No study has yet analyzed the flux variation of each component over solar cycles. In this work, we measure the temporal variations of the flux of each component over 15 yr of Fermi Large Area Telescope observations and compare them with the sunspot number and Galactic cosmic-ray flux from AMS-02 near Earth. We find that the flux variation of the disk anticorrelates with the sunspot number and correlates with cosmic-ray protons, as expected, confirming its emission mechanism. In contrast, the extended component exhibits a more complex variation: despite an initial anticorrelation with the sunspot number, we find neither anticorrelation with the sunspot number nor correlation with cosmic-ray electrons over the full 15 yr period. This most likely suggests that cosmic-ray transport and modulation in the inner heliosphere are unexpectedly complex and may differ for electrons and protons or, alternatively, that there is an additional, unknown component of gamma rays or cosmic rays. These findings impact space weather research and emphasize the need for close monitoring of Cycle 25 and the ongoing polarity reversal.