Abstract:Spiral galaxies M31 and M33 are among the gamma -ray sources detected by the Fermi Large Area Telescope (LAT). We aim to model the broadband non-thermal emission of the central region of M31 (a LAT point source) and of the disk of M33 (a LAT extended source) as part of our continued survey of non-thermal properties of local galaxies that includes the Magellanic Clouds (Persic Rephaeli 2022). We analysed the observed emission from the central region of M31 ($R < 5.5$ kpc) and the disk-sized emission from M33 ($R 9$ kpc). For each galaxy, we self-consistently modelled the broadband spectral energy distribution of the diffuse non-thermal emission based on published radio and gamma -ray data. All relevant radiative processes involving relativistic and thermal electrons (synchrotron, Compton scattering, bremsstrahlung, and free-free emission and absorption), along with relativistic protons ($ decay following interaction with thermal protons), were considered, using exact emissivity formulae. We also used the Fermi/LAT-validated gamma -ray emissivities for pulsars. Joint spectral analyses of the emission from the central region of M31 and the extended disk of M33 indicate that the radio emission is composed of both primary and secondary electron synchrotron and thermal bremsstrahlung, whereas the gamma -ray emission may be explained as a combination of diffuse pionic, pulsar, and nuclear-BH--related emissions in M31 and plain diffuse pionic emission (with an average proton energy density of 0.5 eV cm$^ $) in M33. The observed gamma -ray emission from M33 appears to be mainly hadronic. This situation is similar to other local galaxies, namely, the Magellanic Clouds. In contrast, we have found suggestions of a more complex situation in the central region of M31, whose emission could be an admixture of pulsar emission and hadronic emission, with the latter possibly originating from both the disk and the vicinity of the nuclear black hole. The alternative modelling of the spectra of M31 and M33 is motivated by the different hydrogen distribution in the two galaxies: The hydrogen deficiency in the central region of M31 partially unveils emissions from the nuclear BH and the pulsar population in the bulge and inner disk. If this were to be the case in M33 as well, these emissions would be outshined by diffuse pionic emission originating within the flat central-peak gas distribution in M33.

Searching for GeV gamma-ray emission from the bulge of M31

Non-thermal emission in M31 and M33

Discovery of GeV $\Gamma$-Ray Emission from M33 and Arp 299 with Fermi-LAT

On the Gamma-Ray Emission of the Andromeda Galaxy M31

Fermi-LAT Observations of γ-Ray Emission toward the Outer Halo of M31

Chandra Detection of Diffuse Hot Gas in and Around the M31 Bulge

GeV γ -Ray Emission from M33 and Arp 299

Identifying the Gamma-ray Emission of the Nearby Galaxy M83

Detection of Gamma-ray Halos around Nearby Late-type Galaxies

Diffuse Γ-Ray Emission from the Vicinity of Young Massive Star Cluster RSGC 1

DISCOVERY OF GeV EMISSION FROM THE DIRECTION OF THE LUMINOUS INFRARED GALAXY NGC 2146

Millisecond Pulsar Interpretation of the Galactic Center Gamma-Ray Excess

XMM-Newton RGS Spectroscopy of the M31 Bulge. I. Evidence for a Past AGN Half a Million Years Ago

Search for gamma-ray emission from the 12 nearby dwarf spheroidal galaxies with 12¬†years of Fermi-LAT data

The Fermi-LAT GeV Excess Traces Stellar Mass in the Galactic Bulge

Robust inference of the Galactic centre gamma-ray excess spatial properties

A New Determination of the Millisecond Pulsar Gamma-Ray Luminosity Function and Implications for the Galactic Center Gamma-Ray Excess

The Gamma-Ray Pulsar Population of Globular Clusters: Implications for the GeV Excess

Search for dark matter signals with Fermi-LAT observation of globular clusters NGC 6388 and M 15