Search

Young & massive stellar clusters (SCs) are a potential source of galactic cosmic rays up to very high energies as a result of two possible acceleration scenarios. Collective stellar winds from massive member stars form a wind-blown bubble with a termination shock (TS) at which particle acceleration to PeV energies may occur. Furthermore, shock acceleration may occur at SNRs expanding inside the bubble. By applying a model of CR acceleration at both the wind TS and SNR shocks to catalogues of known SCs derived from Gaia DR2, we identify the most promising targets to search for evidence of PeVatron activity. Predictions for the secondary fluxes of gamma-ray and neutrino emission are derived based on particle acceleration at the collective wind TS and the subsequent hadronic interactions with the surrounding medium. Predictions from our modelling under baseline and optimistic scenarios are compared to data, finding consistent results. We estimate the detection prospects for future gamma-ray and neutrino experiments. We find that degree-scale angular sizes of the wind-blown bubbles are typical, that may pose a challenge for experimental detection. A shortlist of the most promising candidates is provided, with an anticipated flux range. Of order 10 SCs may be detectable with future facilities, and 1-5 could be currently operating as PeVatrons. Of these, three gamma-ray detected SCs have data within our predicted range. Our model can consistently describe gamma-ray measurements of SC emission. Several further as-yet-undetected SCs offer promising targets for future observations, although the flux range allowed by our model can be large (> factor 10). The large angular size of the wind-blown bubble may lead to low surface brightness emission, worsening the problem of source confusion. Nevertheless, we discuss how further work will help to constrain SCs as PeVatron candidates. (abridged), Comment: Submitted to Astronomy and Astrophysics. 19 pages, 10 figures, 5 tables

Context. The recent discovery of several ultra high-energy gamma-ray emitters in our Galaxy represents a significant advancement towards the characterisation of its most powerful accelerators. Nonetheless, in order to unambiguously locate the regions where the highest energy particles are produced and understand the responsible physical mechanisms, detailed spectral and morphological studies are required, especially given that most of the observed sources were found to be significantly extended. Aims. In these regards, pointing observations with the next-generation Imaging Atmospheric Cherenkov Telescopes, such as the Cherenkov Telescope Array (CTA) Observatory and the ASTRI Mini-Array (ASTRI), are expected to provide significant improvements. Here we aim to identify the most promising sources to target in future observations. Methods. For this purpose, we performed a comparative analysis of the expected performance of ASTRI and CTA, computing their differential sensitivities towards extended sources, and further explored their capabilities with respect to specific case studies, including follow-ups of existing gamma-ray source catalogues. Results. We find that almost all of the sources thus far detected by LHAASO-WCDA and in the H.E.S.S. Galactic Plane Survey will be in the reach of ASTRI and CTA with about 300 and 50 hours of exposure, respectively. For the highest energy emitters detected by LHAASO-KM2A, in turn, we provide a list of the most promising objects that would require further investigation. We additionally examined specific classes of sources in order to identify potentially detectable gamma-ray emitters, such as passive molecular clouds (i.e. illuminated by the cosmic-ray sea) and pulsars surrounded by a halo of runaway particles., Comment: Accepted for publication in A&A

Context. Star clusters constitute a significant part of the stellar population in our Galaxy. The feedback processes they exert on the interstellar medium impact multiple physical processes from the chemical to the dynamical evolution of the Galaxy. In addition, young and massive stellar clusters might act as efficient particle accelerators and contribute to the production of cosmic rays. Aims. We aim at evaluating the wind luminosity driven by the young (< 30 Myr) Galactic open stellar clusters observed by the Gaia space mission. This is crucial for determining the energy channeled into accelerated particles. Methods. To do this, we developed a method relying on the number, magnitude, and line-of-sight extinction of the stars observed per cluster. Assuming that the stellar mass function follows a Kroupa mass distribution and accounting for the maximum stellar mass allowed by the age and mass of the parent cluster, we conservatively estimated the mass and wind luminosity of 387 local clusters within the second data release of Gaia. Results. We compared the results of our computation with recent estimates of young cluster masses. With respect to these, our sample is three times more abundant, particularly above a few thousand solar masses. This is of the utmost relevance for predicting the gamma-ray emission resulting from the interaction of accelerated particles. The cluster wind luminosity distribution we obtained extends up to 3x10^38 erg/s. This is a promising feature in terms of potential particle acceleration scenarios., Comment: Accepted for publication in A&A

HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging (<10" FWHM) and broad spectral coverage (0.2-80 keV) with an effective area far superior to current facilities (including XMM-Newton and NuSTAR) to enable revolutionary new insights into a variety of important astrophysical problems. With the recent discoveries of over 40 ultra-high-energy gamma-ray sources (detected above 100 TeV) and neutrino emission in the Galactic Plane, we have entered a new era of multi-messenger astrophysics facing the exciting reality of Galactic PeVatrons. In the next decade, as more Galactic PeVatrons and TeV gamma-ray sources are expected to be discovered, the identification of their acceleration and emission mechanisms will be the most pressing issue in both particle and high-energy astrophysics. In this paper, along with its companion papers (Reynolds et al. 2023, Mori et al. 2023), we will present that HEX-P is uniquely suited to address important problems in various cosmic-ray accelerators, including Galactic PeVatrons, through investigating synchrotron X-ray emission of TeV-PeV electrons produced by both leptonic and hadronic processes., Comment: 31 pages, 12 figures, submitted to FrASS

The Fermi Gamma-ray Space Telescope has revealed a diffuse $\gamma$-ray background at energies from 0.1 GeV to 1 TeV, which can be separated into Galactic emission and an isotropic, extragalactic component. Previous efforts to understand the latter have been hampered by the lack of physical models capable of predicting the $\gamma$-ray emission produced by the many candidate sources, primarily active galactic nuclei and star-forming galaxies, leaving their contributions poorly constrained. Here we present a calculation of the contribution of star-forming galaxies to the $\gamma$-ray background that does not rely on empirical scalings, and is instead based on a physical model for the $\gamma$-ray emission produced when cosmic rays accelerated in supernova remnants interact with the interstellar medium. After validating the model against local observations, we apply it to the observed cosmological star-forming galaxy population and recover an excellent match to both the total intensity and the spectral slope of the $\gamma$-ray background, demonstrating that star-forming galaxies alone can explain the full diffuse, isotropic $\gamma$-ray background., Comment: 18 pages, 10 figures. This work has been published in Nature. The version deposited here is the author's pre-print and may not reflect post-acceptance corrections or formatting related changes. The published version (Version of Record) of this manuscript is available at https://www.nature.com/articles/s41586-021-03802-x

Interstellar clouds can act as target material for hadronic cosmic rays; gamma-rays produced through inelastic proton-proton collisions and spatially associated with the clouds provide a key indicator of efficient particle acceleration. However, even for PeVatron sources reaching PeV energies, the system of cloud and accelerator must fulfil several conditions in order to produce a detectable gamma-ray flux. In this contribution, we characterise the necessary properties of both cloud and accelerator. Using available Supernova Remnant (SNR) and interstellar cloud catalogues, and assuming particle acceleration to PeV energies in a nearby SNR, we produce a ranked shortlist of the most promising target systems; those for which a detectable gamma-ray flux is predicted. We discuss detection prospects for future facilities including CTA and SWGO; and compare our predictions with known gamma-ray sources, including the Ultra-High-Energy sources recently detected by LHAASO. A range of model scenarios are tested, including variation in the diffusion coefficient and particle spectrum, under which the best candidate clouds in our shortlist are consistently bright. On average, a detectable gamma-ray flux is more likely for more massive clouds; for systems with lower separation distance between the SNR and cloud; and for slightly older SNRs, due to the time required for particles to traverse the separation distance., Comment: 8 pages, 5 figures, Proceedings of the 37th International Cosmic Ray Conference (ICRC2021)

The process that allows cosmic rays to escape from their sources and be released into the Galaxy is still largely unknown. The comparison between cosmic-ray electron and proton spectra measured at Earth suggests that electrons are released with a spectrum steeper than protons by $\Delta s_{\rm ep} \sim 0.3$ for energies above $\sim 10$ GeV and by $\Delta s_{\rm ep} \sim 1.2$ above $\sim 1$ TeV. Assuming that both species are accelerated at supernova remnant shocks, we here explore two possible scenarios that can in principle justify steeper electron spectra: (i) energy losses due to synchrotron radiation in an amplified magnetic field, and (ii) time dependent acceleration efficiency. We account for magnetic field amplification produced by either cosmic-ray induced instabilities or by magneto-hydrodynamics instabilities my means of a parametric description. We show that both mechanisms are required to explain the electron spectrum. In particular synchrotron losses can only produce a significant electron steepening above $\sim 1$~TeV, while a time dependent acceleration can explain the spectrum at lower energies if the electron injection into diffusive shock acceleration is inversely proportional to the shock speed. We discuss observational and theoretical evidences supporting such a behaviour. Furthermore, we predict two additional spectral features: a spectral break below $\sim$ few GeV (as required by existing observations) due to the acceleration efficiency drop during the adiabatic phase, and a spectral hardening above $\sim 20$ TeV (where no data are available yet) resulting from electrons escaping from the shock precursor., Comment: 14 pages, 8 figures. Accepted for publication in MNRAS

Multi-messenger astrophysics is becoming a major avenue to explore the Universe, with the potential to span a vast range of redshifts. The growing synergies between different probes is opening new frontiers, which promise profound insights into several aspects of fundamental physics and cosmology. In this context, THESEUS will play a central role during the 2030s in detecting and localizing the electromagnetic counterparts of gravitational wave and neutrino sources that the unprecedented sensitivity of next generation detectors will discover at much higher rates than the present. Here, we review the most important target signals from multi-messenger sources that THESEUS will be able to detect and characterize, discussing detection rate expectations and scientific impact., Comment: Submitted to Experimental Astronomy

The strong constraints from the Fermi-LAT data on the isotropic gamma-ray background suggest that the neutrinos observed by IceCube might possibly come from sources that are hidden to gamma-ray observations. A possibility recently discussed in the literature is that neutrinos may come from jets of collapsing massive stars which fail to break out of the stellar envelope, and for this reason they are known as choked jets, or choked Gamma-Ray Bursts (GRBs). In this paper, we estimate the neutrino flux and spectrum expected from these sources, focusing on Type II SNe. We perform detailed calculations of pg interactions, accounting for all the neutrino production channels and scattering angles. We provide predictions of expected event rates for operating neutrino telescopes, such as ANTARES and IceCube, as well as for the future generation telescope KM3NeT. We find that for GRB energies channeled into protons spanning between 10^51 - 10^53 erg, choked GRBs may substantially contribute to the observed astrophysical neutrino flux, if their local rate is 80 - 1 Gpc^-3 yr^-1 respectively., Comment: 24 pages, 5 figures, published on JCAP

The origin of the radiation observed in the region of the supernova remnant RX J1713.7-3946, one of the brightest TeV emitters, has been debated since its discovery. The existence of atomic and molecular clouds in this object supports the idea that part of the GeV gamma rays in this region originate from proton-proton collisions. However, the observed column density of protons derived from gas observations cannot explain the whole emission. Yet there could be a fraction of protons contained in fainter structures that have note been detected so far. Here we search for faint objects in the line of sight of RX J1713.7-3946 using the principle of light extinction and the ESA/Gaia DR2 astrometric and photometric data. We reveal and locate with precision a number of dust clouds and note that only one appears to be in the vicinity of RX J1713.7-3946. We estimate the embedded mass to $M_{dust} = (7.0 \pm 0.6) \times 10^3 \, M_{\odot}$ which might be big enough to contain the missing protons. Finally, using the fact that the supernova remnant is expected to be located in a dusty environment and that there appears to be only one such structure in the vicinity of RX J1713.7-3946, we set a very precise constrain to the supernova remnant distance, at ($1.12 \pm 0.01$) kpc.

We analyze the energy distributions of final (stable) products - gamma rays, neutrinos, and electrons - produced in inelastic proton-proton collisions in the PeV energy regime. We also calculate the energy spectrum of synchrotron radiation from secondary electrons, assuming that these are promptly cooled in the surrounding magnetic field. We show that the synchrotron radiation has an energy distribution much shallower than that of primary protons, and hence we suggest to take advantage of such a feature in the spectral analysis of the highest energy (cut-off) emission region from particle accelerators. For a broad range of energy distributions in the parent protons, we propose simple analytical presentations for the spectra of secondaries in the cut-off region. These results can be used in the interpretation of high-energy radiation from PeVatrons - cosmic-ray factories accelerating protons to energies up to 1 PeV., Comment: 14 pages, 6 figures. Accepted by ApJ

The KM3NeT research infrastructure is currently under construction at two locations in the Mediterranean Sea. The KM3NeT/ORCA water-Cherenkov neutrino detector off the French coast will instrument several megatons of seawater with photosensors. Its main objective is the determination of the neutrino mass ordering. This work aims at demonstrating the general applicability of deep convolutional neural networks to neutrino telescopes, using simulated datasets for the KM3NeT/ORCA detector as an example. To this end, the networks are employed to achieve reconstruction and classification tasks that constitute an alternative to the analysis pipeline presented for KM3NeT/ORCA in the KM3NeT Letter of Intent. They are used to infer event reconstruction estimates for the energy, the direction, and the interaction point of incident neutrinos. The spatial distribution of Cherenkov light generated by charged particles induced in neutrino interactions is classified as shower- or track-like, and the main background processes associated with the detection of atmospheric neutrinos are recognized. Performance comparisons to machine-learning classification and maximum-likelihood reconstruction algorithms previously developed for KM3NeT/ORCA are provided. It is shown that this application of deep convolutional neural networks to simulated datasets for a large-volume neutrino telescope yields competitive reconstruction results and performance improvements with respect to classical approaches.

The gSeaGen code is a GENIE-based application developed to efficiently generate high statistics samples of events, induced by neutrino interactions, detectable in a neutrino telescope. The gSeaGen code is able to generate events induced by all neutrino flavours, considering topological differences between track-type and shower-like events. Neutrino interactions are simulated taking into account the density and the composition of the media surrounding the detector. The main features of gSeaGen are presented together with some examples of its application within the KM3NeT project., Comment: 45 pages, 10 figures

The escape process of particles accelerated at supernova remnant (SNR) shocks is one of the poorly understood aspects of the shock acceleration theory. Here we adopt a phenomenological approach to study the particle escape and its impact on the gamma-ray spectrum resulting from hadronic collisions both inside and outside of a middle-aged SNR. Under the assumption that in the spatial region immediately outside of the remnant the diffusion coefficient is suppressed with respect to the average Galactic one, we show that a significant fraction of particles are still located inside the SNR long time after their nominal release from the acceleration region. This fact results into a gamma-ray spectrum that resembles a broken power law, similar to those observed in several middle-aged SNRs. Above the break, the spectral steepening is determined by the diffusion coefficient outside of the SNR and by the time dependence of maximum energy. Consequently, the comparison between the model prediction and actual data will contribute to determining these two quantities, the former being particularly relevant within the predictions of the gamma-ray emission from the halo of escaping particles around SNRs which could be detected with future Cherenkov telescope facilities. We also calculate the spectrum of run-away particles injected into the Galaxy by an individual remnant. Assuming that the acceleration stops before the SNR enters the snowplow phase, we show that the released spectrum can be a featureless power law only if the accelerated spectrum has a slope alpha > 4., Comment: 18 pages, 10 figures, accepted for publication in MNRAS

Observations from the radio to the gamma-ray wavelengths indicate that supernova remnant (SNR) shocks are sites of effective particle acceleration. It has been proposed that the presence of dense clumps in the environment where supernovae explode might have a strong impact in shaping the hadronic gamma-ray spectrum. Here we present a detailed numerical study about the penetration of relativistic protons into clumps which are engulfed by a SNR shock, taking into account the magneto-hydrodynamical properties of the background plasma. We show that the spectrum of protons inside clumps is much harder than that in the diffuse inter-clump medium and we discuss the implications for the formation of the spectrum of hadronic gamma rays, which does not reflect anymore the acceleration spectrum of protons, resulting substantially modified inside the clumps due to propagation effects. For the Galactic SNR RX J1713.7-3946, we show that a hadronic scenario including dense clumps inside the remnant shell is able to reproduce the broadband gamma-ray spectrum from GeV to TeV energies. Moreover, we argue that small clumps crossed by the shock could provide a natural explanation to the non-thermal X-ray variability observed in some hot spots of RX J1713.7-3946. Finally we discuss the detectability of gamma-ray emission from clumps with the upcoming Cherenkov Telescope Array and the possible detection of the clumps themselves through molecular lines., Comment: 15 pages, 10 figures, accepted for publication in MNRAS

We discuss the discovery potential of extended very-high-energy (VHE) neutrino sources by the future KM3 Neutrino Telescope (KM3NeT) in the context of the constraining power of the Cherenkov Telescope Array (CTA), designed for deep surveys of the sky in VHE gamma rays. The study is based on a comparative analysis of sensitivities of KM3NeT and CTA. We show that a minimum gamma-ray energy flux of E^2{\phi}_{\gamma}(10 TeV) > 1x10^{-12} TeV cm^{-2} s^{-1} is required to identify a possible neutrino counterpart with a 3{\sigma} significance and 10 years of KM3NeT observations with upgoing muons, if the source has an angular size of R_{src} = 0.1 deg and emits gamma rays with an E^{-2} energy spectrum through a full hadronic mechanism. This minimum gamma-ray flux is increased to the level of E^2{\phi}_{\gamma}(10 TeV) > 2x10^{-11} TeV cm^{-2} s^{-1} in case of sources with radial extension of R_{src} = 2.0 deg. The analysis methods are applied to the supernova remnant RX J1713.7-3946 and the Galactic Center Ridge, as well as to the recent HAWC catalog of multi-TeV gamma-ray sources., Comment: 15 pages, 7 figures

The hypothesis of a PeVatron in the Galactic Center, emerged with the recent {\gamma}-ray measurements of H.E.S.S., motivates the search for neutrinos from this source. The effect of {\gamma}-ray absorption is studied: at the energies currently probed, the known background radiation fields lead to small effects, whereas it is not possible to exclude large effects due to new IR radiation fields near the very Center. Precise upper limits on neutrino fluxes are derived and the underlying hypotheses are discussed. The expected number of events for ANTARES, IceCube and KM3NeT, based on the H.E.S.S. measurements, are calculated. It is shown that km3-class telescopes in the Northern hemisphere have the potential of observing high-energy neutrinos from this important astronomical object and can check the existence of a hadronic PeV galactic accelerator., Comment: 11 pages, 4 figures, 2 tables. Accepted for publication in EPJC