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.

Context. The nearby elliptical galaxy M87 contains one of only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. Aims. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. Methods. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high-energy (VHE) gamma-rays as well as details of the individual observations and light curves. We also conducted phenomenological modelling to investigate the basic source properties. Results. We present the first VHE gamma-ray flare from M87 detected since 2010. The flux above 350 GeV more than doubled within a period of approximate to 36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Conclusions. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE gamma-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and it emphasises the need for combined image and spectral modelling.

A 2.1 yr periodic oscillation of the gamma-ray flux from the blazar PG 1553+113 has previously been tentatively identified in ∼7 yr of data from the Fermi Large Area Telescope. After 15 yr of Fermi sky-survey observations, doubling the total time range, we report >7 cycle gamma-ray modulation with an estimated significance of 4σ against stochastic red noise. Independent determinations of oscillation period and phase in the earlier and the new data are in close agreement (chance probability <0.01). Pulse timing over the full light curve is also consistent with a coherent periodicity. Multiwavelength new data from Swift X-Ray Telescope, Burst Alert Telescope, and UVOT, and from KAIT, Catalina Sky Survey, All-Sky Automated Survey for Supernovae, and Owens Valley Radio Observatory ground-based observatories as well as archival Rossi X-Ray Timing Explorer satellite-All Sky Monitor data, published optical data of Tuorla, and optical historical Harvard plates data are included in our work. Optical and radio light curves show clear correlations with the gamma-ray modulation, possibly with a nonconstant time lag for the radio flux. We interpret the gamma-ray periodicity as possibly arising from a pulsational accretion flow in a sub-parsec binary supermassive black hole system of elevated mass ratio, with orbital modulation of the supplied material and energy in the jet. Other astrophysical scenarios introduced include instabilities, disk and jet precession, rotation or nutation, and perturbations by massive stars or intermediate-mass black holes in polar orbit.

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 use pulsar timing arrays, which perform long-term monitoring of millisecond pulsars at radio wavelengths. We used 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 x 10(-14) at a frequency of 1 year(-1). The sensitivity is expected to scale with t(obs), the observing time span, as t(obs)(-13/6). This direct measurement provides an independent probe of the GWB while offering a check on radio noise models.

Blazars can show variability on a wide range of timescales. However, the search for periodicity in the γ-ray emission of blazars remains an on-going challenge. This contribution will show the results obtained when a systematic pipeline is used to implement ten well-established methods to search for periodicity. We analyze the most promising candidates selected from our previous work, extending the Fermi-LAT light curves over three more years, for a total telescope time of twelve years. These improvements have allowed us to build the first sample of blazars that display a periodicity detected at a significance > 5σ. Finally, we will discuss the potential origins for the periodic behavior observed in blazars.

Massive black holes at the centers of galaxies can launch powerful wide-angle winds that, if sustained over time, can unbind the gas from the stellar bulges of galaxies. These winds may be responsible for the observed scaling relation between the masses of the central black holes and the velocity dispersion of stars in galactic bulges. Propagating through the galaxy, the wind should interact with the interstellar medium creating a strong shock, similar to those observed in supernovae explosions, which is able to accelerate charged particles to high energies. In this work we use data from the Fermi Large Area Telescope to search for the γ-ray emission from galaxies with an ultra-fast outflow (UFO): a fast (v ∼ 0.1c), highly ionized outflow, detected in absorption at hard X-rays in several nearby active galactic nuclei (AGN). Adopting a sensitive stacking analysis we are able to detect the average γ-ray emission from these galaxies and exclude that it is due to processes other than the UFOs. Moreover, our analysis shows that the γ-ray luminosity scales with the AGN bolometric luminosity and that these outflows transfer ∼0.04 % of their mechanical power to γ rays. Interpreting the observed γ-ray emission as produced by cosmic rays (CRs) accelerated at the shock front, we find that the γ-ray emission may attest to the onset of the wind-host interaction and that these outflows can energize charged particles up to the transition region between galactic and extragalactic CRs. A preprint of the full analysis is available on the arXiv: 2105.11469.

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi Large Area Telescope (LAT) catalog of gamma-ray sources. Based on the first 12 years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss the degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources.

Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma-rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four supernova remnants (SNRs), IC 443, W44, W49B, and W51C, with the Fermi Large Area Telescope. This detection provided direct evidence that cosmic-ray protons are (re-)accelerated in SNRs. Here, we present a comprehensive search for low-energy spectral breaks among 311 4FGL catalog sources located within 5 degrees from the Galactic plane. Using 8 yr of data from the Fermi Large Area Telescope between 50 MeV and 1 GeV, we find and present the spectral characteristics of 56 sources with a spectral break confirmed by a thorough study of systematic uncertainty. Our population of sources includes 13 SNRs for which the proton-proton interaction is enhanced by the dense target material; the high-mass gamma-ray binary LS I+61 303; the colliding wind binary eta Carinae; and the Cygnus star-forming region. This analysis better constrains the origin of the gamma-ray emission and enlarges our view to potential new cosmic-ray acceleration sites.

An incremental version of the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi Large Area Telescope is presented. This version (4LAC-DR3) derives from the third data release of the 4FGL catalog based on 12 yr of E > 50 MeV gamma-ray data, where the spectral parameters, spectral energy distributions (SEDs), yearly light curves, and associations have been updated for all sources. The new reported AGNs include 587 blazar candidates and four radio galaxies. We describe the properties of the new sample and outline changes affecting the previously published one. We also introduce two new parameters in this release, namely the peak energy of the SED high-energy component and the corresponding flux. These parameters allow an assessment of the Compton dominance, the ratio of the inverse-Compton to the synchrotron-peak luminosities, without relying on X-ray data.

We present the first Fermi Large Area Telescope (LAT) catalog of long-term gamma-ray transient sources (1FLT). This comprises sources that were detected on monthly time intervals during the first decade of Fermi-LAT operations. The monthly timescale allows us to identify transient and variable sources that were not yet reported in other Fermi-LAT catalogs. The monthly data sets were analyzed using a wavelet-based source detection algorithm that provided the candidate new transient sources. The search was limited to the extragalactic regions of the sky to avoid the dominance of the Galactic diffuse emission at low Galactic latitudes. The transient candidates were then analyzed using the standard Fermi-LAT maximum likelihood analysis method. All sources detected with a statistical significance above 4 sigma in at least one monthly bin were listed in the final catalog. The 1FLT catalog contains 142 transient gamma-ray sources that are not included in the 4FGL-DR2 catalog. Many of these sources (102) have been confidently associated with active galactic nuclei (AGNs): 24 are associated with flat-spectrum radio quasars, 1 with a BL Lac object, 70 with blazars of uncertain type, 3 with radio galaxies, 1 with a compact steep-spectrum radio source, 1 with a steep-spectrum radio quasar, and 2 with AGNs of other types. The remaining 40 sources have no candidate counterparts at other wavelengths. The median gamma-ray spectral index of the 1FLT-AGN sources is softer than that reported in the latest Fermi-LAT AGN general catalog. This result is consistent with the hypothesis that detection of the softest gamma-ray emitters is less efficient when the data are integrated over year-long intervals.