Showing 1–50 of 101 results for author: Morlino, G

Contribution of young massive star clusters to Galactic diffuse $γ$-ray emission

Authors: Stefano Menchiari, Giovanni Morlino, Elena Amato, Niccolò Bucciantini, Giada Peron, Germano Sacco

Abstract: Context: Young massive stellar clusters (YMSCs) have emerged as potential $γ$-ray sources, after the recent association of a dozen YMSCs with extended $γ$-ray emission. The large size of the detected halos, comparable to that of the wind-blown bubble expected around YMSCs, makes the $γ$-ray detection of individual YMSCs rather challenging. As a result, the emission from most of the Galactic YMSCs… ▽ More Context: Young massive stellar clusters (YMSCs) have emerged as potential $γ$-ray sources, after the recent association of a dozen YMSCs with extended $γ$-ray emission. The large size of the detected halos, comparable to that of the wind-blown bubble expected around YMSCs, makes the $γ$-ray detection of individual YMSCs rather challenging. As a result, the emission from most of the Galactic YMSCs could be unresolved, thus contributing to the diffuse $γ$-ray radiation observed along the Galactic Plane. Aims: In this study, we estimate the possible contribution to the Galactic diffuse $γ$-ray emission from a synthetic population of YMSCs, and we compare it with observations obtained with different experiments, from 1 GeV to hundreds of TeV, in two regions of the Galactic Plane. Methods: As the population of galactic YMSCs is only known locally, we evaluate the contribution of $γ$-ray emission relying on the simulation of synthetic populations of YMSCs based on the observed properties of local clusters. We compute the $γ$-ray emission from each cluster assuming that the radiation is purely hadronic in nature and produced by cosmic rays accelerated at the cluster's collective wind termination shock. Results: We find that the $γ$-ray emission from unresolved YMSCs can significantly contribute to the observed Galactic diffuse flux, especially in the inner part of the Galaxy. The result is independent of the assumed particle transport, but an important role is played by Wolf-Rayet stars. The predicted $γ$-ray flux should be considered as a lower limit, given that our calculation neglects the contribution of supernovae exploding in YMSCs. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: Submitted to A&A letter on 06/05/24

Probing Stellar Clusters from Gaia DR2 as Galactic PeVatrons: I -- Expected Gamma-ray and Neutrino Emission

Authors: Alison M. W. Mitchell, Giovanni Morlino, Silvia Celli, Stefano Menchiari, Andreas Specovius

Abstract: 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 bu… ▽ More 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) △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: Submitted to Astronomy and Astrophysics. 19 pages, 10 figures, 5 tables

Cygnus OB2 as a test case for particle acceleration in young massive star clusters

Authors: Stefano Menchiari, Giovanni Morlino, Elena Amato, Niccolò Bucciantini, Maria Teresa Beltrán

Abstract: In this paper, we focus on the scientific case of Cygnus OB2, a northern sky young massive stellar cluster (YMSC) located towards the Cygnus X star-forming complex. We consider a model that assumes cosmic ray acceleration occurring only at the termination shock of the collective wind of the YMSC and address the question of whether, and under what hypotheses, hadronic emission by the accelerated pa… ▽ More In this paper, we focus on the scientific case of Cygnus OB2, a northern sky young massive stellar cluster (YMSC) located towards the Cygnus X star-forming complex. We consider a model that assumes cosmic ray acceleration occurring only at the termination shock of the collective wind of the YMSC and address the question of whether, and under what hypotheses, hadronic emission by the accelerated particles can account for the observations of Cygnus OB2 obtained by Fermi-LAT, HAWC and LHAASO. In order to do so, we carefully review the available information on this source, also confronting different estimates of the relevant parameters with ad hoc developed simulations. Once other model parameters are fixed, the spectral and spatial properties of the emission are found to be very sensitive to the unknown properties of the turbulent magnetic field. Comparison with the data shows that our suggested scenario is incompatible with Kolmogorov turbulence. Assuming Kraichnan or Bohm type turbulence spectra, the model accounts well for the Very High Energy (VHE) data, but fails to reproduce the centrally peaked morphology observed by Fermi-LAT, suggesting that additional effects might be important for lower energy $γ$-ray emission. We discuss how additional progress can be made with a more detailed and extended knowledge of the spectral and morphological properties of the emission. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: Submitted to A&A. 15 pages, 8 figures plus appendices

On the mass and wind luminosity of young Galactic open clusters in Gaia DR2

Authors: Silvia Celli, Andreas Specovius, Stefano Menchiari, Alison Mitchell, Giovanni Morlino

Abstract: Context. Star clusters constitute a relevant 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, possibly contributing to the production of cosmic rays. Ai… ▽ More Context. Star clusters constitute a relevant 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, possibly contributing 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, which is crucial to determine the energy channeled into accelerated particles. Methods. To this extent, we develop 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 both the parent cluster age and mass, we conservatively estimate the mass and wind luminosity of 387 local clusters within the second data release of Gaia. Results. We compare the results of our computation with recent estimations of young cluster masses. With respect to these, we provide a sample three times more abundant, particularly above a few thousand solar masses, which is of the utmost relevance for predicting the gamma-ray emission resulting from the interaction of accelerated particles. In fact, the cluster wind luminosity distribution we obtain is found to extend up to 3 x 10^38 erg/s, a promising feature in terms of potential particle acceleration scenarios. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 11 pages, 16 figures, comments are welcome

Journal ref: A&A 686, A118 (2024)

Contribution to Galactic cosmic rays from young stellar clusters

Authors: G. Morlino, S. Menchiari, E. Amato, N. Bucciantini

Abstract: The origin of Galactic cosmic rays (CR) is still a matter of debate. Diffusive shock acceleration (DSA) applied to supernova remnant (SNR) shocks provides the most reliable explanation. However, within the current understanding of DSA several issues remain unsolved, like the CR maximum energy, the chemical composition and the transition region between Galactic and extra-Galactic CRs. These issues… ▽ More The origin of Galactic cosmic rays (CR) is still a matter of debate. Diffusive shock acceleration (DSA) applied to supernova remnant (SNR) shocks provides the most reliable explanation. However, within the current understanding of DSA several issues remain unsolved, like the CR maximum energy, the chemical composition and the transition region between Galactic and extra-Galactic CRs. These issues motivate the search for other possible Galactic sources. Recently, several young stellar clusters (YSC) have been detected in gamma rays, suggesting that such objects could be powerful sources of Galactic CRs. The energy input could come from winds of massive stars hosted in the clusters which is a function of the cluster total mass and initial mass function of stars. In this work we evaluate the total CR flux produced by a synthetic population of YSCs assuming that the CR acceleration occurs at the termination shock of the collective wind resulting from the sum of cluster's stellar winds. We show that the spectrum produced by YSC can significantly contribute to energies $\gtrsim 100$ TeV if the diffusion inside the wind-blown bubble is Bohm-like and the spectral slope is harder than the one produced by SNRs. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: Proceeding to the International Cosmic Ray Conference, ICRC 2023, Nagoya, Japan

Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, K. Abe, S. Abe, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, E. O. Angüner, L. A. Antonelli, C. Aramo, M. Araya, C. Arcaro, L. Arrabito, K. Asano, Y. Ascasíbar, J. Aschersleben , et al. (542 additional authors not shown)

Abstract: Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med… ▽ More Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 93 pages (including author list, appendix and references), 143 figures. Submitted to JCAP

Different spectra of cosmic ray H, He and heavier nuclei escaping compact star clusters

Authors: Pasquale Blasi, Giovanni Morlino

Abstract: Cosmic ray acceleration at the termination shock of compact star clusters has recently received much attention, mainly because of the detection of gamma ray emission from some of such astrophysical sources. Here we focus on the acceleration of nuclei at the termination shock and we investigate the role played by proton energy losses and spallation reactions of nuclei, especially downstream of the… ▽ More Cosmic ray acceleration at the termination shock of compact star clusters has recently received much attention, mainly because of the detection of gamma ray emission from some of such astrophysical sources. Here we focus on the acceleration of nuclei at the termination shock and we investigate the role played by proton energy losses and spallation reactions of nuclei, especially downstream of the shock. We show that for a rather generic choice of the mean gas density in the cavity excavated by the cluster wind, the spectrum of He nuclei is systematically harder than the spectrum of hydrogen, in a manner that appears to be qualitatively consistent with the observed and yet unexplained phenomenon of discrepant hardening. We also find that the spallation reactions of heavier nuclei are likely to be so severe that their spectra become very hard and with a low normalization, meaning that it is unlikely that heavy nuclei escaping star clusters can provide a sizeable contribution to the spectrum of cosmic rays at the Earth. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 9 Pages, Submitted to MNRAS

High energy cosmic rays and gamma rays from star clusters: the case of Cygnus OB2

Authors: Pasquale Blasi, Giovanni Morlino

Abstract: We investigate the acceleration of cosmic rays at the termination shock that results from the interaction of the collective wind of star clusters with the surrounding interstellar medium. The solution of the transport equation of accelerated particles in the wind-excavated cavity, including energy losses due to CR interactions with neutral gas in the bubble, shows several interesting properties th… ▽ More We investigate the acceleration of cosmic rays at the termination shock that results from the interaction of the collective wind of star clusters with the surrounding interstellar medium. The solution of the transport equation of accelerated particles in the wind-excavated cavity, including energy losses due to CR interactions with neutral gas in the bubble, shows several interesting properties that are discussed in detail. The issue of the maximum energy of the accelerated particles is discussed with special care, because of its implications for the origin of Galactic cosmic rays. Gamma ray emission is produced in the cavity due to inelastic pp scattering, while accelerated particles are advected downstream of the termination shock and diffuse at the same time. Both the spectrum and the morphology of such emission are discussed, with a comparison of our results with the observations of gamma ray emission from the Cygnus OB2 region. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

Authors: The Cherenkov Telescope Array Consortium, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Aloisio, N. Álvarez Crespo, R. Alves Batista, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, C. Aramo, C. Arcaro, T. Armstrong, K. Asano, Y. Ascasibar, J. Aschersleben, M. Backes, A. Baktash, C. Balazs, M. Balbo , et al. (334 additional authors not shown)

Abstract: The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The pote… ▽ More The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy $γ$-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte--Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a $γ$-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 hours of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with $\mathcal{O}(100)$ hours of exposure per source. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, Accepted for publication in Astroparticle Physics

Diffusive shock acceleration at EeV and associated multimessenger flux from ultra-fast outflows driven by Active Galactic Nuclei

Authors: Enrico Peretti, Alessandra Lamastra, Francesco Gabriele Saturni, Markus Ahlers, Pasquale Blasi, Giovanni Morlino, Pierre Cristofari

Abstract: Active galactic nuclei (AGN) can launch and sustain powerful winds featuring mildly relativistic velocity and wide opening angle. Such winds, known as ultra-fast outflows (UFOs), can develop a bubble structure characterized by a forward shock expanding in the host galaxy and a wind termination shock separating the fast cool wind from the hot shocked wind. In this work we explore whether diffusive… ▽ More Active galactic nuclei (AGN) can launch and sustain powerful winds featuring mildly relativistic velocity and wide opening angle. Such winds, known as ultra-fast outflows (UFOs), can develop a bubble structure characterized by a forward shock expanding in the host galaxy and a wind termination shock separating the fast cool wind from the hot shocked wind. In this work we explore whether diffusive shock acceleration can take place efficiently at the wind termination shock of UFOs. We calculate the spectrum of accelerated particles and find that protons can be energized up to the EeV range promoting UFOs to promising candidates for accelerating ultra-high energy cosmic rays (UHECRs). We also compute the associated gamma-ray and neutrino fluxes and compare them with available data in the literature. We observe that high-energy (HE) neutrinos are efficiently produced up to hundreds of PeV while the associated gamma rays could be efficiently absorbed beyond a few tens of GeV by the optical-ultraviolet AGN photon field. By assuming a typical source density of non-jetted AGN we expect that UFOs could play a dominant role as diffuse sources of UHECRs and HE neutrinos. We finally apply our model to the recently observed NGC1068 and we find out that under specific parametric conditions an obscured UFO could provide a sizeable contribution to the observed gamma-ray flux while only contributing up to ~10 per cent to the associated neutrino flux. △ Less

Submitted 22 September, 2023; v1 submitted 31 January, 2023; originally announced January 2023.

Comments: 12 pages, 6 figures, 2 tables; submitted to MNRAS; matching version published in MNRAS

Galactic Observatory Science with the ASTRI Mini-Array at the Observatorio del Teide

Authors: A. D'Aì, E. Amato, A. Burtovoi, A. A. Compagnino, M. Fiori, A. Giuliani, N. La Palombara, A. Paizis, G. Piano, F. G. Saturni, A. Tutone, A. Belfiore, M. Cardillo, S. Crestan, G. Cusumano, M. Della Valle, M. Del Santo, A. La Barbera, V. La Parola, S. Lombardi, S. Mereghetti, G. Morlino, F. Pintore, P. Romano, S. Vercellone , et al. (30 additional authors not shown)

Abstract: The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Mini-Array will be composed of nine imaging atmospheric Cherenkov telescopes at the Observatorio del Teide site. The array will be best suited for astrophysical observations in the 0.3-200 TeV range with an angular resolution of few arc-minutes and an energy resolution of 10-15\%. A core-science programme in the first four years… ▽ More The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Mini-Array will be composed of nine imaging atmospheric Cherenkov telescopes at the Observatorio del Teide site. The array will be best suited for astrophysical observations in the 0.3-200 TeV range with an angular resolution of few arc-minutes and an energy resolution of 10-15\%. A core-science programme in the first four years will be devoted to a limited number of key targets, addressing the most important open scientific questions in the very-high energy domain. At the same time, thanks to a wide field of view of about 10 degrees, ASTRI Mini-Array will observe many additional field sources, which will constitute the basis for the long-term observatory programme that will eventually cover all the accessible sky. In this paper, we review different astrophysical Galactic environments, e.g. pulsar wind nebulae, supernova remnants, and gamma-ray binaries, and show the results from a set of ASTRI Mini-Array simulations of some of these field sources made to highlight the expected performance of the array (even at large offset angles) and the important additional observatory science that will complement the core-science program. △ Less

Submitted 5 August, 2022; originally announced August 2022.

Comments: Accepted for publication in Journal of High-Energy Astrophysics

Journal ref: JHEAP, Volume 35, August 2022, Pages 139-175

ASTRI Mini-Array Core Science at the Observatorio del Teide

Authors: S. Vercellone, C. Bigongiari, A. Burtovoi, M. Cardillo, O. Catalano, A. Franceschini, S. Lombardi, L. Nava, F. Pintore, A. Stamerra, F. Tavecchio, L. Zampieri, R. Alves Batista, E. Amato, L. A. Antonelli, C. Arcaro, J. Becerra Gonzalez, G. Bonnoli, M. Bottcher, G. Brunetti, A. A. Compagnino, S. Crestan, A. D Ai, M. Fiori, G. Galanti , et al. (62 additional authors not shown)

Abstract: The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Project led by the Italian National Institute for Astrophysics (INAF) is developing and will deploy at the Observatorio del Teide a mini-array (ASTRI Mini-Array) composed of nine telescopes similar to the small-size dual-mirror Schwarzschild-Couder telescope (ASTRI-Horn) currently operating on the slopes of Mt. Etna in Sicily.… ▽ More The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Project led by the Italian National Institute for Astrophysics (INAF) is developing and will deploy at the Observatorio del Teide a mini-array (ASTRI Mini-Array) composed of nine telescopes similar to the small-size dual-mirror Schwarzschild-Couder telescope (ASTRI-Horn) currently operating on the slopes of Mt. Etna in Sicily. The ASTRI Mini-Array will surpass the current Cherenkov telescope array differential sensitivity above a few tera-electronvolt (TeV), extending the energy band well above hundreds of TeV. This will allow us to explore a new window of the electromagnetic spectrum, by convolving the sensitivity performance with excellent angular and energy resolution figures. In this paper we describe the Core Science that we will address during the first four years of operation, providing examples of the breakthrough results that we will obtain when dealing with current open questions, such as the acceleration of cosmic rays, cosmology and fundamental physics and the new window, for the TeV energy band, of the time-domain astrophysics. △ Less

Submitted 5 August, 2022; originally announced August 2022.

Comments: Published in Journal of High Energy Astrophysics. 46 Figures, 7 Tables

Journal ref: Journal of High Energy Astrophysics, Volume 35, August 2022, Pages 1-42

Extragalactic Observatory Science with the ASTRI Mini-Array at the Observatorio del Teide

Authors: F. G. Saturni, C. H. E. Arcaro, B. Balmaverde, J. Becerra González, A. Caccianiga, M. Capalbi, A. Lamastra, S. Lombardi, F. Lucarelli, R. Alves Batista, L. A. Antonelli, E. M. de Gouveia Dal Pino, R. Della Ceca, J. G. Green, A. Pagliaro, C. Righi, F. Tavecchio, S. Vercellone, A. Wolter, E. Amato, C. Bigongiari, M. Böttcher, G. Brunetti, P. Bruno, A. Bulgarelli , et al. (25 additional authors not shown)

Abstract: The ASTRI Mini-Array is a next-generation system of nine imaging atmospheric Cherenkov telescopes that is going to be built at the Observatorio del Teide site. After a first phase, in which the instrument will be operated as an experiment prioritizing a schedule of primary science cases, an observatory phase is foreseen in which other significant targets will be pointed. We focus on the observatio… ▽ More The ASTRI Mini-Array is a next-generation system of nine imaging atmospheric Cherenkov telescopes that is going to be built at the Observatorio del Teide site. After a first phase, in which the instrument will be operated as an experiment prioritizing a schedule of primary science cases, an observatory phase is foreseen in which other significant targets will be pointed. We focus on the observational feasibility of extragalactic sources and on astrophysical processes that best complement and expand the ASTRI Mini-Array core science, presenting the most relevant examples that are at reach of detection over long-term time scales and whose observation can provide breakthrough achievements in the very-high energy extragalactic science. Such examples cover a wide range of $γ$-ray emitters, including the study of AGN low states in the multi-TeV energy range, the possible detection of Seyfert galaxies with long exposures and the searches of dark matter lines above 10 TeV. Simulations of the presented objects show that the instrument performance will be competitive at multi-TeV energies with respect to current arrays of Cherenkov telescopes. △ Less

Submitted 5 August, 2022; originally announced August 2022.

Comments: 27 pages, 10 figures, 7 tables, published on JHEAp

Journal ref: Journal of High Energy Astrophysics 35 (2022), 91-111

Forbidden Line Emission from Type Ia Supernova Remnants Containing Balmer-Dominated Shells

Authors: Chuan-Jui Li, You-Hua Chu, John C. Raymond, Bruno Leibundgut, Ivo R. Seitenzahl, Giovanni Morlino

Abstract: Balmer-dominated shells in supernova remnants (SNRs) are produced by collisionless shocks advancing into a partially neutral medium, and are most frequently associated with Type Ia supernovae. We have analyzed Hubble Space Telescope (HST) images and VLT/MUSE or AAT/WiFeS observations of five Type Ia SNRs containing Balmer-dominated shells in the LMC: 0509-67.5, 0519-69.0, N103B, DEM L71, and 0548-… ▽ More Balmer-dominated shells in supernova remnants (SNRs) are produced by collisionless shocks advancing into a partially neutral medium, and are most frequently associated with Type Ia supernovae. We have analyzed Hubble Space Telescope (HST) images and VLT/MUSE or AAT/WiFeS observations of five Type Ia SNRs containing Balmer-dominated shells in the LMC: 0509-67.5, 0519-69.0, N103B, DEM L71, and 0548-70.4. Contrary to expectations, we find bright forbidden line emission from small dense knots embedded in four of these SNRs. The electron densities in some knots are higher than 10$^4$ cm$^{-3}$. The size and density of these knots are not characteristic for interstellar medium (ISM) -- they most likely originate from a circumstellar medium (CSM) ejected by the SN progenitor. Physical property variations of dense knots in the SNRs appear to reflect an evolutionary effect. The recombination timescales for high densities are short, and HST images of N103B taken 3.5 yr apart already show brightness changes in some knots. VLT/MUSE observations detect [Fe XIV] line emission from reverse shocks into SN ejecta as well as forward shocks into the dense knots. Faint [O III] line emission is also detected from the Balmer shell in 0519-69.0, N103B, and DEM L71. We exclude the postshock origin because the [O III] line is narrow. For the preshock origin, we considered three possibilities: photoionization precursor, cosmic ray precursor, and neutral precursor. We conclude that the [O III] emission arises from oxygen that has been photoionized by [He II] $λ$304 photons and is then collisionally excited in a shock precursor heated mainly by cosmic rays. △ Less

Submitted 12 November, 2021; v1 submitted 29 September, 2021; originally announced October 2021.

Comments: 20 pages, 2 tables, 14 figures, accepted for publication in The Astrophysical Journal (ApJ)

Particle acceleration at the termination shock of stellar clusters' wind

Authors: Giovanni Morlino

Abstract: We investigate the process of particle acceleration at the termination shock that develops in the bubble excavated by winds of star clusters in the interstellar medium. We develop a theory of diffusive shock acceleration at such shock and we find that the maximum energy may reach the PeV region for very powerful clusters. We show how the maximum energy is limited by two different processes: the pa… ▽ More We investigate the process of particle acceleration at the termination shock that develops in the bubble excavated by winds of star clusters in the interstellar medium. We develop a theory of diffusive shock acceleration at such shock and we find that the maximum energy may reach the PeV region for very powerful clusters. We show how the maximum energy is limited by two different processes: the particle escape from the bubble boundary and the drop of energy gain for particles able to diffuse up to the center of the cluster. A crucial role in this problem is played by the dissipation of kinetic energy of the wind to magnetic perturbations which determines the diffusion regime of particles: in case the diffusion is close to Bohm than PeV energies can be reached. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: Proceeding of the 37th International Cosmic Ray Conference (ICRC2021)

Report number: https://pos.sissa.it/395/444

Feedback on the ASTRONET Science Vision and Infrastructure Roadmap from the CTA Consortium

Authors: A. Araudo, A. Carosi, W. Hofmann, F. Iocco, J. -P. Lenain, E. Lindfors, A. Lopez, M. Meyer, G. Morlino, B. Olmi, P. Romano, M. Santander, L. Tibaldo, R. Zanin

Abstract: Feedback on the ASTRONET Science Vision and Infrastructure Roadmap from the CTA Consortium. Feedback on the ASTRONET Science Vision and Infrastructure Roadmap from the CTA Consortium. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: Submitted as input to the ASTRONET Science Vision and Infrastructure roadmap on behalf of the CTA consortium

Cosmic-ray electrons released by supernova remnants

Authors: Giovanni Morlino, Silvia Celli

Abstract: 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 $Δs_{\rm ep} \sim 0.3$ for energies above $\sim 10$ GeV and by $Δs_{\rm ep} \sim 1.2$ above $\sim 1$ TeV. Assuming that… ▽ More 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 $Δs_{\rm ep} \sim 0.3$ for energies above $\sim 10$ GeV and by $Δ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. △ Less

Submitted 13 October, 2021; v1 submitted 11 June, 2021; originally announced June 2021.

Comments: 14 pages, 8 figures. Accepted for publication in MNRAS

Probing extreme environments with the Cherenkov Telescope Array

Authors: C. Boisson, A. M. Brown, A. Burtovoi, M. Cerruti, M. Chernyakova, T. Hassan, J. -P. Lenain, M. Manganaro, P. Romano, H. Sol, F. Tavecchio, S. Vercellone, L. Zampieri, R. Zanin, A. Zech, I. Agudo, R. Alves Batista, E. O. Anguner, L. A. Antonelli, M. Backes, C. Balazs, J. Becerra González, C. Bigongiari, E. Bissaldi, J. Bolmont , et al. (105 additional authors not shown)

Abstract: The physics of the non-thermal Universe provides information on the acceleration mechanisms in extreme environments, such as black holes and relativistic jets, neutron stars, supernovae or clusters of galaxies. In the presence of magnetic fields, particles can be accelerated towards relativistic energies. As a consequence, radiation along the entire electromagnetic spectrum can be observed, and ex… ▽ More The physics of the non-thermal Universe provides information on the acceleration mechanisms in extreme environments, such as black holes and relativistic jets, neutron stars, supernovae or clusters of galaxies. In the presence of magnetic fields, particles can be accelerated towards relativistic energies. As a consequence, radiation along the entire electromagnetic spectrum can be observed, and extreme environments are also the most likely sources of multi-messenger emission. The most energetic part of the electromagnetic spectrum corresponds to the very-high-energy (VHE, E>100 GeV) gamma-ray regime, which can be extensively studied with ground based Imaging Atmospheric Cherenkov Telescopes (IACTs). The results obtained by the current generation of IACTs, such as H.E.S.S., MAGIC, and VERITAS, demonstrate the crucial importance of the VHE band in understanding the non-thermal emission of extreme environments in our Universe. In some objects, the energy output in gamma rays can even outshine the rest of the broadband spectrum. The Cherenkov Telescope Array (CTA) is the next generation of IACTs, which, with cutting edge technology and a strategic configuration of ~100 telescopes distributed in two observing sites, in the northern and southern hemispheres, will reach better sensitivity, angular and energy resolution, and broader energy coverage than currently operational IACTs. With CTA we can probe the most extreme environments and considerably boost our knowledge of the non-thermal Universe. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: Submitted as input to ASTRONET Science Vision and Infrastructure roadmap on behalf of the CTA consortium

Grammage of cosmic rays in the proximity of supernova remnants embedded in a partially ionized medium

Authors: S. Recchia, D. Galli, L. Nava, M. Padovani, S. Gabici, A. Marcowith, V. Ptuskin, G. Morlino

Abstract: We investigate the damping of Alfvén waves generated by the cosmic ray resonant streaming instability in the context of the cosmic ray escape and propagation in the proximity of supernova remnants. We consider ion-neutral damping, turbulent damping and non linear Landau damping in the warm ionized and warm neutral phases of the interstellar medium. For the ion-neutral damping, up-to-date damping c… ▽ More We investigate the damping of Alfvén waves generated by the cosmic ray resonant streaming instability in the context of the cosmic ray escape and propagation in the proximity of supernova remnants. We consider ion-neutral damping, turbulent damping and non linear Landau damping in the warm ionized and warm neutral phases of the interstellar medium. For the ion-neutral damping, up-to-date damping coefficients are used. We investigate in particular whether the self-confinement of cosmic rays nearby sources can appreciably affect the grammage. We show that the ion-neutral damping and the turbulent damping effectively limit the residence time of cosmic rays in the source proximity, so that the grammage accumulated near sources is found to be negligible. Contrary to previous results, this also happens in the most extreme scenario where ion-neutral damping is less effective, namely in a medium with only neutral helium and fully ionized hydrogen. Therefore, the standard picture, in which CR secondaries are produced during the whole time spent by cosmic rays throughout the Galactic disk, need not to be deeply revisited. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 13 pages, 9 figures, submitted to MNRAS

Multi-messenger and transient astrophysics with the Cherenkov Telescope Array

Authors: Ž. Bošnjak, A. M. Brown, A. Carosi, M. Chernyakova, P. Cristofari, F. Longo, A. López-Oramas, M. Santander, K. Satalecka, F. Schüssler, O. Sergijenko, A. Stamerra, I. Agudo, R. Alves Batista, E. Amato, E. O. Anguner, L. A. Antonelli, M. Backes, Csaba Balazs, L. Baroncelli, J. Becker Tjus, C. Bigongiari, E. Bissaldi, C. Boisson, J. Bolmont , et al. (120 additional authors not shown)

Abstract: The discovery of gravitational waves, high-energy neutrinos or the very-high-energy counterpart of gamma-ray bursts has revolutionized the high-energy and transient astrophysics community. The development of new instruments and analysis techniques will allow the discovery and/or follow-up of new transient sources. We describe the prospects for the Cherenkov Telescope Array (CTA), the next-generati… ▽ More The discovery of gravitational waves, high-energy neutrinos or the very-high-energy counterpart of gamma-ray bursts has revolutionized the high-energy and transient astrophysics community. The development of new instruments and analysis techniques will allow the discovery and/or follow-up of new transient sources. We describe the prospects for the Cherenkov Telescope Array (CTA), the next-generation ground-based gamma-ray observatory, for multi-messenger and transient astrophysics in the decade ahead. CTA will explore the most extreme environments via very-high-energy observations of compact objects, stellar collapse events, mergers and cosmic-ray accelerators. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: Submitted to ASTRONET roadmap on behalf of the CTA consortium

Origin and role of relativistic cosmic particles

Authors: A. Araudo, G. Morlino, B. Olmi, F. Acero, I. Agudo, R. Adam, R. Alves Batista, E. Amato, E. O. Anguner, L. A. Antonelli, Y. Ascasibar, C. Balazs, J. Becker Tjus, C. Bigongiari, E. Bissaldi, J. Bolmont, C. Boisson, P. Bordas, Ž. Bošnjak, A. M. Brown, M. Burton, N. Bucciantini, F. Cangemi, P. Caraveo, M. Cardillo , et al. (99 additional authors not shown)

Abstract: This white paper briefly summarizes the importance of the study of relativistic cosmic rays, both as a constituent of our Universe, and through their impact on stellar and galactic evolution. The focus is on what can be learned over the coming decade through ground-based gamma-ray observations over the 20 GeV to 300 TeV range. The majority of the material is drawn directly from "Science with the C… ▽ More This white paper briefly summarizes the importance of the study of relativistic cosmic rays, both as a constituent of our Universe, and through their impact on stellar and galactic evolution. The focus is on what can be learned over the coming decade through ground-based gamma-ray observations over the 20 GeV to 300 TeV range. The majority of the material is drawn directly from "Science with the Cherenkov Telescope Array", which describes the overall science case for CTA. We request that authors wishing to cite results contained in this white paper cite the original work. △ Less

Submitted 15 June, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: Submitted as input to ASTRONET Science Vision and Infrastructure roadmap on behalf of the CTA consortium. arXiv admin note: text overlap with arXiv:1709.07997

Particle acceleration and multimessenger emission from starburst-driven galactic winds

Authors: Enrico Peretti, Giovanni Morlino, Pasquale Blasi, Pierre Cristofari

Abstract: The enhanced star forming activity, typical of starburst galaxies, powers strong galactic winds expanding on kiloparsec (kpc) scales and characterized by bubble structures. Here we discuss the possibility that particle acceleration may take place at the termination shock of such winds. We calculate the spectrum of accelerated particles and their maximum energy, that turns out to range up to a few… ▽ More The enhanced star forming activity, typical of starburst galaxies, powers strong galactic winds expanding on kiloparsec (kpc) scales and characterized by bubble structures. Here we discuss the possibility that particle acceleration may take place at the termination shock of such winds. We calculate the spectrum of accelerated particles and their maximum energy, that turns out to range up to a few hundred petaelectronvolt (PeV) for typical values of the parameters. Cosmic rays accelerated at the termination shock are advected towards the edge of the bubble excavated by the wind and eventually escape into extragalactic space. We also calculate the flux of gamma rays and neutrinos produced by hadronic interactions in the bubble as well as the diffuse flux resulting from the superposition of the contribution of starburst galaxies on cosmological scales. Finally, we compute the diffuse flux of cosmic rays from starburst bubbles and compare it with existing data. △ Less

Submitted 13 February, 2022; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 13 pages, 7 figures, 2 tables. Accepted for publication on MNRAS

Particle acceleration in winds of star clusters

Authors: Giovanni Morlino, Pasquale Blasi, Enrico Peretti, Pierre Cristofari

Abstract: The origin of cosmic rays in our Galaxy remains a subject of active debate. While supernova remnant shocks are often invoked as the sites of acceleration, it is now widely accepted that the difficulties of such sources in reaching PeV energies are daunting and it seems likely that only a subclass of rare remnants can satisfy the necessary conditions. Moreover the spectra of cosmic rays escaping th… ▽ More The origin of cosmic rays in our Galaxy remains a subject of active debate. While supernova remnant shocks are often invoked as the sites of acceleration, it is now widely accepted that the difficulties of such sources in reaching PeV energies are daunting and it seems likely that only a subclass of rare remnants can satisfy the necessary conditions. Moreover the spectra of cosmic rays escaping the remnants have a complex shape that is not obviously the same as the spectra observed at the Earth. Here we investigate the process of particle acceleration at the termination shock that develops in the bubble excavated by star clusters' winds in the interstellar medium. While the main limitation to the maximum energy in supernova remnants comes from the need for effective wave excitation upstream so as to confine particles in the near-shock region and speed up the acceleration process, at the termination shock of star clusters the confinement of particles upstream in guaranteed by the geometry of the problem. We develop a theory of diffusive shock acceleration at such shock and we find that the maximum energy may reach the PeV region for powerful clusters in the high end of the luminosity tail for these sources. A crucial role in this problem is played by the dissipation of energy in the wind to magnetic perturbations. Under reasonable conditions the spectrum of the accelerated particles has a power law shape with a slope $4÷4.3$, in agreement with what is required based upon standard models of cosmic ray transport in the Galaxy. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 11 pages, 5 figures, submitted to MNRAS, comments welcome

Study of the GeV to TeV morphology of the $γ$-Cygni SNR (G78.2+2.1) with MAGIC and Fermi-LAT

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (161 additional authors not shown)

Abstract: Context. Diffusive shock acceleration (DSA) is the most promising mechanism to accelerate Galactic cosmic rays (CRs) in the shocks of supernova remnants (SNRs). The turbulence upstream is supposedly generated by the CRs, but this process is not well understood. The dominant mechanism may depend on the evolutionary state of the shock and can be studied via the CRs escaping upstream into the interst… ▽ More Context. Diffusive shock acceleration (DSA) is the most promising mechanism to accelerate Galactic cosmic rays (CRs) in the shocks of supernova remnants (SNRs). The turbulence upstream is supposedly generated by the CRs, but this process is not well understood. The dominant mechanism may depend on the evolutionary state of the shock and can be studied via the CRs escaping upstream into the interstellar medium (ISM). Aims. Previous observations of the $γ$-Cygni SNR showed a difference in morphology between GeV and TeV energies. Since this SNR has the right age and is at the evolutionary stage for a significant fraction of CRs to escape, we aim to understand $γ$-ray emission in the vicinity of the $γ$-Cygni SNR. Methods. We observed the region of the $γ$-Cygni SNR with the MAGIC Imaging Atmospheric Cherenkov telescopes between May 2015 and September 2017 recording 87 h of good-quality data. Additionally we analysed Fermi-LAT data to study the energy dependence of the morphology as well as the energy spectrum in the GeV to TeV range. The energy spectra and morphology were compared against theoretical predictions, which include a detailed derivation of the CR escape process and their $γ$-ray generation. Results. The MAGIC and Fermi-LAT data allowed us to identify three emission regions, which can be associated with the SNR and dominate at different energies. Our hadronic emission model accounts well for the morphology and energy spectrum of all source components. It constrains the time-dependence of the maximum energy of the CRs at the shock, the time-dependence of the level of turbulence, and the diffusion coefficient immediately outside the SNR shock. While in agreement with the standard picture of DSA, the time-dependence of the maximum energy was found to be steeper than predicted and the level of turbulence was found to change over the lifetime of the SNR. △ Less

Submitted 29 October, 2020; originally announced October 2020.

Comments: 20 pages, 10 figures. Accepted for publication in A&A. Corresponding authors: M. Strzys, G. Morlino, S. Masuda, and I. Vovk

Journal ref: A&A 670, A8 (2023)

New high-frequency radio observations of the Cygnus Loop supernova remnant with the Italian radio telescopes

Authors: S. Loru, A. Pellizzoni, E. Egron, A. Ingallinera, G. Morlino, S. Celli, G. Umana, C. Trigilio, P. Leto, M. N. Iacolina, S. Righini, P. Reich, S. Mulas, M. Marongiu, M. Pilia, A. Melis, R. Concu, M. Bufano, C. Buemi, F. Cavallaro, S. Riggi, F. Schillirò

Abstract: Supernova remnants (SNRs) represent a powerful laboratory to study the Cosmic-Ray acceleration processes at the shocks, and their relation to the properties of the circumstellar medium. With the aim of studying the high-frequency radio emission and investigating the energy distribution of accelerated electrons and the magnetic field conditions, we performed single-dish observations of the large an… ▽ More Supernova remnants (SNRs) represent a powerful laboratory to study the Cosmic-Ray acceleration processes at the shocks, and their relation to the properties of the circumstellar medium. With the aim of studying the high-frequency radio emission and investigating the energy distribution of accelerated electrons and the magnetic field conditions, we performed single-dish observations of the large and complex Cygnus Loop SNR from 7.0 to 24.8 GHz with the Medicina and the Sardinia Radio Telescope, focusing on the northern filament (NGC 6992) and the southern shell. Both regions show a spectrum well fitted by a power-law function ($S\proptoν^{-α}$), with spectral index $α=0.45\pm0.05$ for NGC 6992 and $α=0.49\pm0.01$ for the southern shell and without any indication of a spectral break. The spectra are significantly flatter than the whole Cygnus Loop spectrum ($α=0.54\pm0.01$), suggesting a departure from the plain shock acceleration mechanisms, which for NGC 6992 could be related to the ongoing transition towards a radiative shock. We model the integrated spectrum of the whole SNR considering the evolution of the maximum energy and magnetic field amplification. Through the radio spectral parameters, we infer a magnetic field at the shock of 10 $μ$G. This value is compatible with a pure adiabatic compression of the interstellar magnetic field, suggesting that the amplification process is currently inefficient. △ Less

Submitted 21 September, 2020; originally announced September 2020.

Comments: 19 pages, 13 figures

Impact of transport modelling on the $^{60}$Fe abundance inside Galactic cosmic ray sources

Authors: Giovanni Morlino, Elena Amato

Abstract: The ACE-CRIS collaboration has recently released the measurement of radioactive $^{60}$Fe nuclei abundance in Galactic Cosmic Rays, in the energy range $\sim 195-500$ MeV per nucleon. We model Cosmic Ray propagation and derive from this measurement the $^{60}$Fe/$^{56}$Fe ratio that is expected in the sources of Galactic Cosmic Rays. We describe Cosmic Ray origin and transport within the framework… ▽ More The ACE-CRIS collaboration has recently released the measurement of radioactive $^{60}$Fe nuclei abundance in Galactic Cosmic Rays, in the energy range $\sim 195-500$ MeV per nucleon. We model Cosmic Ray propagation and derive from this measurement the $^{60}$Fe/$^{56}$Fe ratio that is expected in the sources of Galactic Cosmic Rays. We describe Cosmic Ray origin and transport within the framework of the disk/halo diffusion model, namely a scenario in which the matter and the Cosmic Ray sources in our Galaxy are confined to a thin disk, while Cosmic Ray propagation occurs in a much larger halo with negligible matter density. We solve the Cosmic Ray transport equation accounting for spallation reactions, decay and ionization losses as well as advection. We find that the $^{60}$Fe/$^{56}$Fe ratio at the source must be very close to the value detected in the local Cosmic Ray spectrum at Earth, due to the fact that spallation reactions are more effective for $^{56}$Fe than for $^{60}$Fe. Such a result could help identify the sources of Galactic Cosmic Rays. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: 11 pages, 5 figures. Accepted for publication in PRD

Cosmic-ray anisotropies in right ascension measured by the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker , et al. (351 additional authors not shown)

Abstract: We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic e… ▽ More We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$σ$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $α_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin. △ Less

Submitted 14 February, 2020; originally announced February 2020.

Comments: Accepted by ApJ

Report number: FERMILAB-PUB-20-076

Journal ref: Published in Astrophys. J 891 (2020)142

Contribution of starburst nuclei to the diffuse gamma-ray and neutrino flux

Authors: Enrico Peretti, Pasquale Blasi, Felix Aharonian, Giovanni Morlino, Pierre Cristofari

Abstract: In nuclei of starburst galaxies, the combination of an enhanced rate of supernova explosions and a high gas density suggests that cosmic rays can be efficiently produced, and that most of them lose their energy before escaping these regions, resulting in a large flux of secondary products, including neutrinos. Although the flux inferred from an individual starburst region is expected to be well be… ▽ More In nuclei of starburst galaxies, the combination of an enhanced rate of supernova explosions and a high gas density suggests that cosmic rays can be efficiently produced, and that most of them lose their energy before escaping these regions, resulting in a large flux of secondary products, including neutrinos. Although the flux inferred from an individual starburst region is expected to be well below the sensitivity of current neutrino telescopes, such sources may provide a substantial contribution to the diffuse neutrino flux measured by IceCube. Here we compute the gamma-ray and neutrino flux due to starburst galaxies based on a physical model of cosmic ray transport in a starburst nucleus, and accounting for the redshift evolution of the number density of starburst sources as inferred from recent measurements of the star formation rate. The model accounts for gamma-ray absorption both inside the sources and in the intergalactic medium. The latter process is responsible for electromagnetic cascades, which also contribute to the diffuse gamma-ray background at lower energies. The conditions for acceleration of cosmic ray protons up to energies exceeding $ \sim 10 \, \rm PeV$ in starburst regions, necessary for the production of PeV neutrinos, are investigated in a critical way. We show that starburst nuclei can account for the diffuse neutrino flux above $\sim 200 \, \rm TeV$, thereby producing $\lesssim 40 \%$ of the extragalactic diffuse gamma-ray background. Below $\sim 200 \, \rm TeV$, the flux from starburst appears to be somewhat lower than the observed one, where both the Galactic contribution and the flux of atmospheric neutrinos may account for the difference. △ Less

Submitted 25 March, 2020; v1 submitted 14 November, 2019; originally announced November 2019.

Comments: Version accepted for publication in MNRAS, 13 pages, 5 figures

Constraining the cosmic ray spectrum in the vicinity of the supernova remnant W28: from sub-GeV to multi-TeV energies

Authors: V. H. M. Phan, S. Gabici, G. Morlino, R. Terrier, J. Vink, J. Krause, M. Menu

Abstract: Supernova remnants interacting with molecular clouds are ideal laboratories to study the acceleration of particles at shock waves and their transport and interactions in the surrounding interstellar medium. In this paper, we focus on the supernova remnant W28, which over the years has been observed in all energy domains from radio waves to very-high-energy gamma rays. The bright gamma-ray emission… ▽ More Supernova remnants interacting with molecular clouds are ideal laboratories to study the acceleration of particles at shock waves and their transport and interactions in the surrounding interstellar medium. In this paper, we focus on the supernova remnant W28, which over the years has been observed in all energy domains from radio waves to very-high-energy gamma rays. The bright gamma-ray emission detected from molecular clouds located in its vicinity revealed the presence of accelerated GeV and TeV particles in the region. An enhanced ionization rate has also been measured by means of millimetre observations, but such observations alone cannot tell us whether the enhancement is due to low energy (MeV) cosmic rays (either protons or electrons) or the X-ray photons emitted by the shocked gas. The goal of this study is to determine the origin of the enhanced ionization rate and to infer from multiwavelength observations the spectrum of cosmic rays accelerated at the supernova remnant shock in the unprecedented range spanning from MeV to multi-TeV particle energies. We developed a model to describe the transport of X-ray photons into the molecular cloud, and we fitted the radio, millimeter, and gamma-ray data to derive the spectrum of the radiating particles. The contribution from X-ray photons to the enhanced ionization rate is negligible, and therefore the ionization must be due to cosmic rays. Even though we cannot exclude a contribution to the ionization rate coming from cosmic ray electrons, we show that a scenario where cosmic ray protons explain both the gamma-ray flux and the enhanced ionization rate provides the most natural fit to multiwavelength data. This strongly suggests that the intensity of CR protons is enhanced in the region for particle energies in a very broad range covering almost 6 orders of magnitude: from $\lesssim 100$ MeV up to several tens of TeV. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 8 pages, 6 figures, submitted

Journal ref: A&A 635, A40 (2020)

AMS-02 beryllium data and its implication for cosmic ray transport

Authors: Carmelo Evoli, Giovanni Morlino, Pasquale Blasi, Roberto Aloisio

Abstract: The flux of unstable secondary cosmic ray nuclei, produced by spallation processes in the interstellar medium, can be used to constrain the residence time of cosmic rays inside the Galaxy. Among them, $^{10}$Be is especially useful because of its relatively long half-life of 1.39 Myr. In the framework of the diffusive halo model we describe cosmic ray transport taking into account all relevant int… ▽ More The flux of unstable secondary cosmic ray nuclei, produced by spallation processes in the interstellar medium, can be used to constrain the residence time of cosmic rays inside the Galaxy. Among them, $^{10}$Be is especially useful because of its relatively long half-life of 1.39 Myr. In the framework of the diffusive halo model we describe cosmic ray transport taking into account all relevant interaction channels and accounting for the decay of unstable secondary nuclei. We then compare our results with the data collected by the Alpha Magnetic Spectrometer (AMS-02) on board the International Space Station for the flux ratios Be/C, B/C, Be/O, B/O, C/O and Be/B as well as C, N and O absolute fluxes. These measurements, and especially the Be/B ratio, allow us to single out the flux of $^{10}$Be and infer a best fit propagation time of CRs in the Galaxy. Our results show that, if the cross sections for the production of secondary elements through spallation are taken at face value, AMS-02 measurements are compatible with the standard picture based on CR diffusion in a halo of size $H \gtrsim 5$ kpc. Taking into account the uncertainties in the cross sections, this conclusion becomes less reliable, although still compatible with the standard picture. Implications of our findings for alternative models of CR transport are discussed. △ Less

Submitted 10 January, 2020; v1 submitted 9 October, 2019; originally announced October 2019.

Comments: 11 pages, 6 figures. Accepted for publication on Physical Review D

Journal ref: Phys. Rev. D 101, 023013 (2020)

The Pierre Auger Observatory: Contributions to the 36th International Cosmic Ray Conference (ICRC 2019)

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (361 additional authors not shown)

Abstract: Contributions of the Pierre Auger Collaboration to the 36th International Cosmic Ray Conference (ICRC 2019), 24 July - 1 August 2019, Madison, Wisconsin, USA. Contributions of the Pierre Auger Collaboration to the 36th International Cosmic Ray Conference (ICRC 2019), 24 July - 1 August 2019, Madison, Wisconsin, USA. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Comments: compendium of 29 proceedings for the ICRC 2019

Exploring particle escape in supernova remnants through gamma rays

Authors: Silvia Celli, Giovanni Morlino, Stefano Gabici, Felix Aharonian

Abstract: 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 immedia… ▽ More 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. △ Less

Submitted 10 October, 2019; v1 submitted 22 June, 2019; originally announced June 2019.

Comments: 18 pages, 10 figures, accepted for publication in MNRAS

Probing the origin of ultra-high-energy cosmic rays with neutrinos in the EeV energy range using the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (367 additional authors not shown)

Abstract: Neutrinos with energies above $10^{17}$ eV are detectable with the Surface Detector Array of the Pierre Auger Observatory. The identification is efficiently performed for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for Earth-skimming $τ$ neutrinos with nearly tangential trajectories relative to the earth. No neutrino candidates were found in… ▽ More Neutrinos with energies above $10^{17}$ eV are detectable with the Surface Detector Array of the Pierre Auger Observatory. The identification is efficiently performed for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for Earth-skimming $τ$ neutrinos with nearly tangential trajectories relative to the earth. No neutrino candidates were found in $\sim\,14.7$ years of data taken up to 31 August 2018. This leads to restrictive upper bounds on their flux. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high-energy neutrinos with an $E_ν^{-2}$ spectrum in the energy range $1.0 \times 10^{17}~{\rm eV} - 2.5 \times 10^{19}~{\rm eV}$ is $E^2 {\rm d}N_ν/{\rm d}E_ν< 4.4 \times 10^{-9}~{\rm GeV~cm^{-2}~s^{-1}~sr^{-1}}$, placing strong constraints on several models of neutrino production at EeV energies and on the properties of the sources of ultra-high-energy cosmic rays. △ Less

Submitted 16 October, 2019; v1 submitted 18 June, 2019; originally announced June 2019.

Comments: Published version, 27 pages, 7 figures, 2 tables

Report number: FERMILAB-PUB-19-280

Journal ref: JCAP 10 (2019) 022

Limits on point-like sources of ultra-high-energy neutrinos with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (367 additional authors not shown)

Abstract: With the Surface Detector array (SD) of the Pierre Auger Observatory we can detect neutrinos with energy between $10^{17}\,$eV and $10^{20}\,$eV from point-like sources across the sky, from close to the Southern Celestial Pole up to $60^\circ$ in declination, with peak sensitivities at declinations around $\sim -53^\circ$ and $\sim+55^\circ$, and an unmatched sensitivity for arrival directions in… ▽ More With the Surface Detector array (SD) of the Pierre Auger Observatory we can detect neutrinos with energy between $10^{17}\,$eV and $10^{20}\,$eV from point-like sources across the sky, from close to the Southern Celestial Pole up to $60^\circ$ in declination, with peak sensitivities at declinations around $\sim -53^\circ$ and $\sim+55^\circ$, and an unmatched sensitivity for arrival directions in the Northern hemisphere. A search has been performed for highly-inclined air showers induced by neutrinos of all flavours with no candidate events found in data taken between 1 Jan 2004 and 31 Aug 2018. Upper limits on the neutrino flux from point-like steady sources have been derived as a function of source declination. An unrivaled sensitivity is achieved in searches for transient sources with emission lasting over an hour or less, if they occur within the field of view corresponding to the zenith angle range between $60^\circ$ and $~95^\circ$ where the SD of the Pierre Auger Observatory is most sensitive to neutrinos. △ Less

Submitted 8 November, 2019; v1 submitted 18 June, 2019; originally announced June 2019.

Comments: 25 pages, 9 figures

Report number: FERMILAB-PUB-19-279

Journal ref: JCAP 11 (2019) 004

Multi-Messenger Physics with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (368 additional authors not shown)

Abstract: An overview of the multi-messenger capabilities of the Pierre Auger Observatory is presented. The techniques and performance of searching for Ultra-High Energy neutrinos, photons and neutrons are described. Some of the most relevant results are reviewed, such as stringent upper bounds that were placed to a flux of diffuse cosmogenic neutrinos and photons, bounds placed on neutrinos emitted from co… ▽ More An overview of the multi-messenger capabilities of the Pierre Auger Observatory is presented. The techniques and performance of searching for Ultra-High Energy neutrinos, photons and neutrons are described. Some of the most relevant results are reviewed, such as stringent upper bounds that were placed to a flux of diffuse cosmogenic neutrinos and photons, bounds placed on neutrinos emitted from compact binary mergers that were detected by LIGO and Virgo during their first and second observing runs, as well as searches for high energy photons and neutrons from the Galactic center that constrain the properties of the putative Galactic PeVatron, observed by the H.E.S.S.\ collaboration. The observation of directional correlations between ultra-high energy cosmic rays and either high energy astrophysical neutrinos or specific source populations, weighted by their electromagnetic radiation, are also discussed. They constitute additional multi-messenger approaches aimed at identifying the sources of high energy cosmic rays. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: 28 pages, 7 figures, Review published in Frontiers in Astronomy and Space Sciences

Report number: FERMILAB-PUB-19-135

Journal ref: April 2019 | Volume 6 | Article 24

Effects of reacceleration and source grammage on secondary cosmic rays spectra

Authors: Virginia Bresci, Elena Amato, Pasquale Blasi, Giovanni Morlino

Abstract: The ratio between secondary and primary cosmic ray particles is the main source of information about cosmic ray propagation in the Galaxy. Primary cosmic rays are thought to be accelerated mainly in Supernova Remnant (SNR) shocks and then released in the interstellar medium (ISM). Here they produce secondary particles by occasional collisions with interstellar matter. As a result, the ratio betwee… ▽ More The ratio between secondary and primary cosmic ray particles is the main source of information about cosmic ray propagation in the Galaxy. Primary cosmic rays are thought to be accelerated mainly in Supernova Remnant (SNR) shocks and then released in the interstellar medium (ISM). Here they produce secondary particles by occasional collisions with interstellar matter. As a result, the ratio between the fluxes of secondary and primary particles carries information about the amount of matter cosmic rays have encountered during their journey from their sources to Earth. Recent measurements by AMS-02 revealed an unexpected behaviour of two main secondary-to-primary ratios, the Boron-to-Carbon ratio and the anti-proton-to-proton ratio. In this work we discuss how such anomalies may reflect the action of two phenomena that are usually overlooked, namely the fact that some fraction of secondary particles can be produced within the acceleration region, and the non-negligible probability that secondary particles encounter an accelerator (and are reaccelerated) during propagation. Both effects must be taken into account in order to correctly extract information about CR transport from secondary-to-primary ratios. △ Less

Submitted 23 April, 2019; originally announced April 2019.

Comments: 11 pages, 9 figures, submitted to MNRAS

Monte Carlo studies for the optimisation of the Cherenkov Telescope Array layout

Authors: A. Acharyya, I. Agudo, E. O. Angüner, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, J. -P. Amans, L. Amati, E. Amato, G. Ambrosi, L. A. Antonelli, C. Aramo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, H. Ashkar, C. Balazs, M. Balbo, B. Balmaverde, P. Barai, A. Barbano, M. Barkov , et al. (445 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possi… ▽ More The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Comments: 48 pages, 16 figures, accepted for publication in Astroparticle Physics

Data-driven estimation of the invisible energy of cosmic ray showers with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (367 additional authors not shown)

Abstract: The determination of the primary energy of extensive air showers using the fluorescence detection technique requires an estimation of the energy carried away by particles that do not deposit all their energy in the atmosphere. This estimation is typically made using Monte Carlo simulations and thus depends on the assumed primary particle mass and on model predictions for neutrino and muon producti… ▽ More The determination of the primary energy of extensive air showers using the fluorescence detection technique requires an estimation of the energy carried away by particles that do not deposit all their energy in the atmosphere. This estimation is typically made using Monte Carlo simulations and thus depends on the assumed primary particle mass and on model predictions for neutrino and muon production. In this work we present a new method to obtain the invisible energy from events detected by the Pierre Auger Observatory. The method uses measurements of the muon number at ground level, and it allows us to reduce significantly the systematic uncertainties related to the mass composition and the high energy hadronic interaction models, and consequently to improve the estimation of the energy scale of the Observatory. △ Less

Submitted 6 November, 2019; v1 submitted 23 January, 2019; originally announced January 2019.

Comments: Published version, 18 pages, 10 figures, 4 tables

Report number: FERMILAB-PUB-19-040

Journal ref: Phys. Rev. D 100, 082003 (2019)

Cosmic ray transport and radiative processes in nuclei of starburst galaxies

Authors: Enrico Peretti, Pasquale Blasi, Felix Aharonian, Giovanni Morlino

Abstract: The high rate of star formation and supernova explosions of starburst galaxies make them interesting sources of high energy radiation. Depending upon the level of turbulence present in their interstellar medium, the bulk of cosmic rays produced inside starburst galaxies may lose most of their energy before escaping, thereby making these sources behave as calorimeters, at least up to some maximum e… ▽ More The high rate of star formation and supernova explosions of starburst galaxies make them interesting sources of high energy radiation. Depending upon the level of turbulence present in their interstellar medium, the bulk of cosmic rays produced inside starburst galaxies may lose most of their energy before escaping, thereby making these sources behave as calorimeters, at least up to some maximum energy. Contrary to previous studies, here we investigate in detail the conditions under which cosmic ray confinement may be effective for electrons and nuclei and we study the implications of cosmic ray confinement in terms of multifrequency emission from starburst nuclei and production of high energy neutrinos. The general predictions are then specialized to three cases of active starbursts, namely M82, NGC253 and Arp220. Both primary and secondary electrons, as well as electron-positron pairs produced by gamma ray absorption inside starburst galaxies are taken into account. Electrons and positrons produced as secondary products of hadronic interactions are found to be responsible for most of the emission of leptonic origin. In particular, synchrotron emission of very high energy secondary electrons produces an extended emission of hard X-rays that represent a very interesting signature of hadronic process in starburst galaxies, potentially accessible to current and future observations in the X-ray band. A careful understanding of both the production and absorption of gamma rays in starburst galaxies is instrumental to the assessment of the role of these astrophysical sources as sources of high energy astrophysical neutrinos. △ Less

Submitted 26 April, 2019; v1 submitted 5 December, 2018; originally announced December 2018.

Comments: Version accepted for publication in MNRAS

Interplay between Physics and Geometry in Balmer filaments: the Case of SN 1006

Authors: Rino Bandiera, Giovanni Morlino, Sladjana Knežević, John C. Raymond

Abstract: The analysis of Balmer-dominated emission in supernova remnants is potentially a very powerful way to derive information on the shock structure, on the physical conditions of the ambient medium and on the cosmic-ray acceleration efficiency. However, the outcome of models developed in plane-parallel geometry is usually not easily comparable with the data, since they often come from regions with rat… ▽ More The analysis of Balmer-dominated emission in supernova remnants is potentially a very powerful way to derive information on the shock structure, on the physical conditions of the ambient medium and on the cosmic-ray acceleration efficiency. However, the outcome of models developed in plane-parallel geometry is usually not easily comparable with the data, since they often come from regions with rather a complex geometry. We present here a general scheme to disentangle physical and geometrical effects in the data interpretation, which is especially powerful when the transition zone of the shock is spatially resolved and the spectral resolution is high enough to allow a detailed investigation of spatial changes of the line profile. We then apply this technique to re-analyze very high quality data of a region along the northwestern limb of the remnant of SN~1006. We show how some observed features, previously interpreted only in terms of spatial variations of physical quantities, naturally arise from geometrical effects. With these effects under control, we derive new constraints on physical quantities in the analyzed region, like the ambient density (in the range 0.03-$0.1{\,\rm cm^{-3}}$), the upstream neutral fraction (more likely in the range 0.01-0.1), the level of face-on surface brightness variations (with factors up to $\sim 3$) and the typical scale lengths related to such variations ($\ge 0.1{\,\rm pc}$, corresponding to angular scales $\ge 10{\,\rm arcsec}$). △ Less

Submitted 28 November, 2018; originally announced November 2018.

Comments: 22 pages, 24 figures, 3 tables, accepted on 2018 November 21 for publication on MNRAS

Measurement of the average shape of longitudinal profiles of cosmic-ray air showers at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvar ez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (363 additional authors not shown)

Abstract: The profile of the longitudinal development of showers produced by ultra-high energy cosmic rays carries information related to the interaction properties of the primary particles with atmospheric nuclei. In this work, we present the first measurement of the average shower profile in traversed atmospheric depth at the Pierre Auger Observatory. The shapes of profiles are well reproduced by the Gais… ▽ More The profile of the longitudinal development of showers produced by ultra-high energy cosmic rays carries information related to the interaction properties of the primary particles with atmospheric nuclei. In this work, we present the first measurement of the average shower profile in traversed atmospheric depth at the Pierre Auger Observatory. The shapes of profiles are well reproduced by the Gaisser-Hillas parametrization within the range studied, for E > 10^{17.8} eV. A detailed analysis of the systematic uncertainties is performed using 10 years of data and a full detector simulation. The average shape is quantified using two variables related to the width and asymmetry of the profile, and the results are compared with predictions of hadronic interaction models for different primary particles. △ Less

Submitted 16 May, 2019; v1 submitted 12 November, 2018; originally announced November 2018.

Comments: Published version, 20 pages, 6 figures, 2 tables

Report number: FERMILAB-PUB-18-633

Journal ref: JCAP03(2019)018

Numerical simulations of mass loading in the tails of Bow Shock Pulsar Wind Nebulae

Authors: Barbara Olmi, Niccolò Bucciantini, Giovanni Morlino

Abstract: When a pulsar is moving through a partially ionized medium, a fraction of neutral Hydrogen atoms penetrate inside the pulsar wind and can be photo-ionized by the nebula UV radiation. The resulting protons remains attached to the magnetic field of the light leptonic pulsar wind enhancing its inertia and changing the flow dynamics of the wind. We present here the first numerical simulations of such… ▽ More When a pulsar is moving through a partially ionized medium, a fraction of neutral Hydrogen atoms penetrate inside the pulsar wind and can be photo-ionized by the nebula UV radiation. The resulting protons remains attached to the magnetic field of the light leptonic pulsar wind enhancing its inertia and changing the flow dynamics of the wind. We present here the first numerical simulations of such effect in the tails of bow shock nebulae. We produce a set of different models representative of pulsars moving in the interstellar medium with different velocities, from highly subsonic to supersonic, by means of 2D hydrodynamic relativistic simulations. We compare the different tail morphologies with results from theoretical models of mass loading in bow shocks. As predicted by analytical models we observe a fast sideways expansion of the tail with the formation of secondary shocks in the ISM. This effect could be at the origin of the head-and-shoulder morphology observed in many BSPWNe. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 7 pages, 5 figures, 1 table

Large-scale cosmic-ray anisotropies above 4 EeV measured by the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker, J. A. Bellido , et al. (352 additional authors not shown)

Abstract: We present a detailed study of the large-scale anisotropies of cosmic rays with energies above 4 EeV measured using the Pierre Auger Observatory. For the energy bins [4,8] EeV and $E\geq 8$ EeV, the most significant signal is a dipolar modulation in right ascension at energies above 8 EeV, as previously reported. In this paper we further scrutinize the highest-energy bin by splitting it into three… ▽ More We present a detailed study of the large-scale anisotropies of cosmic rays with energies above 4 EeV measured using the Pierre Auger Observatory. For the energy bins [4,8] EeV and $E\geq 8$ EeV, the most significant signal is a dipolar modulation in right ascension at energies above 8 EeV, as previously reported. In this paper we further scrutinize the highest-energy bin by splitting it into three energy ranges. We find that the amplitude of the dipole increases with energy above 4 EeV. The growth can be fitted with a power law with index $β=0.79\pm 0.19$. The directions of the dipoles are consistent with an extragalactic origin of these anisotropies at all the energies considered. Additionally we have estimated the quadrupolar components of the anisotropy: they are not statistically significant. We discuss the results in the context of the predictions from different models for the distribution of ultrahigh-energy sources and cosmic magnetic fields. △ Less

Submitted 5 December, 2018; v1 submitted 10 August, 2018; originally announced August 2018.

Comments: Published version, 16 pages, 7 figures, 8 tables

Report number: FERMILAB-PUB-18-393

Journal ref: ApJ 868 (2018) 4

Self-Generated Cosmic-Ray Confinement in TeV Halos: Implications for TeV Gamma-Ray Emission and the Positron Excess

Authors: Carmelo Evoli, Tim Linden, Giovanni Morlino

Abstract: Recent observations have detected extended TeV gamma-ray emission surrounding young and middle-aged pulsars. The morphology of these "TeV halos" requires cosmic-ray diffusion to be locally suppressed by a factor of ~100-1000 compared to the typical interstellar medium. No model currently explains this suppression. We show that cosmic-ray self-confinement can significantly inhibit diffusion near pu… ▽ More Recent observations have detected extended TeV gamma-ray emission surrounding young and middle-aged pulsars. The morphology of these "TeV halos" requires cosmic-ray diffusion to be locally suppressed by a factor of ~100-1000 compared to the typical interstellar medium. No model currently explains this suppression. We show that cosmic-ray self-confinement can significantly inhibit diffusion near pulsars. The steep cosmic-ray gradient generates Alfven waves that resonantly scatter the same cosmic-ray population, suppressing diffusion within ~20 pc of pulsars younger than ~100 kyr. In this model, TeV halos evolve through two phases, a growth phase where Alfven waves are resonantly generated and cosmic-ray diffusion becomes increasingly suppressed, and a subsequent relaxation phase where the diffusion coefficient returns to the standard interstellar value. Intriguingly, cosmic-rays are not strongly confined early in the TeV halo evolution, allowing a significant fraction of injected e+e- to escape. If these e+e- also escape from the surrounding supernova remnant, they would provide a natural explanation for the positron excess observed by PAMELA and AMS-02. Recently created TeV cosmic-rays are confined in the TeV halo, matching observations by HAWC and H.E.S.S. While our default model relaxes too rapidly to explain the confinement of TeV cosmic rays around mature pulsars, such as Geminga, models utilizing a Kraichnan turbulence spectrum experience much slower relaxation. Thus, observations of TeV halos around mature pulsars may provide a probe into our understanding of interstellar turbulence. △ Less

Submitted 26 September, 2018; v1 submitted 24 July, 2018; originally announced July 2018.

Comments: 10 pages, 3 figures, published on PRD

Journal ref: Phys. Rev. D 98, 063017 (2018)

Observation of inclined EeV air showers with the radio detector of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, N. Arsene, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (370 additional authors not shown)

Abstract: With the Auger Engineering Radio Array (AERA) of the Pierre Auger Observatory, we have observed the radio emission from 561 extensive air showers with zenith angles between 60$^\circ$ and 84$^\circ$. In contrast to air showers with more vertical incidence, these inclined air showers illuminate large ground areas of several km$^2$ with radio signals detectable in the 30 to 80\,MHz band. A compariso… ▽ More With the Auger Engineering Radio Array (AERA) of the Pierre Auger Observatory, we have observed the radio emission from 561 extensive air showers with zenith angles between 60$^\circ$ and 84$^\circ$. In contrast to air showers with more vertical incidence, these inclined air showers illuminate large ground areas of several km$^2$ with radio signals detectable in the 30 to 80\,MHz band. A comparison of the measured radio-signal amplitudes with Monte Carlo simulations of a subset of 50 events for which we reconstruct the energy using the Auger surface detector shows agreement within the uncertainties of the current analysis. As expected for forward-beamed radio emission undergoing no significant absorption or scattering in the atmosphere, the area illuminated by radio signals grows with the zenith angle of the air shower. Inclined air showers with EeV energies are thus measurable with sparse radio-antenna arrays with grid sizes of a km or more. This is particularly attractive as radio detection provides direct access to the energy in the electromagnetic cascade of an air shower, which in case of inclined air showers is not accessible by arrays of particle detectors on the ground. △ Less

Submitted 24 October, 2018; v1 submitted 14 June, 2018; originally announced June 2018.

Comments: Published version, 14 pages, 12 figures

Report number: FERMILAB-PUB-18-259

Journal ref: JCAP10(2018)026

Origin of the Cosmic Ray Galactic Halo Driven by Advected Turbulence and Self-Generated Waves

Authors: Carmelo Evoli, Pasquale Blasi, Giovanni Morlino, Roberto Aloisio

Abstract: The diffusive paradigm for the transport of Galactic cosmic rays is central to our understanding of the origin of these high energy particles. However, it is worth recalling that the normalization, energy dependence, and spatial extent of the diffusion coefficient in the interstellar medium are fitted to the data and typically are not derived from more basic principles. Here, we discuss a scenario… ▽ More The diffusive paradigm for the transport of Galactic cosmic rays is central to our understanding of the origin of these high energy particles. However, it is worth recalling that the normalization, energy dependence, and spatial extent of the diffusion coefficient in the interstellar medium are fitted to the data and typically are not derived from more basic principles. Here, we discuss a scenario in which the diffusion properties of cosmic rays are derived from a combination of wave self-generation and advection from the Galactic disc, where the sources of cosmic rays are assumed to be located. We show for the first time that a halo naturally arises from these phenomena, with a size of a few kiloparsecs, compatible with the value that typically best fits observations in simple parametric approaches to cosmic ray diffusion. We also show that transport in such a halo results in a hardening in the spectra of primary cosmic rays at ~300 GV. △ Less

Submitted 10 July, 2018; v1 submitted 11 June, 2018; originally announced June 2018.

Comments: 4 pages, 3 figures, accepted for publication in PRL

Journal ref: Phys. Rev. Lett. (2018) 121, 021102

Supernova remnants in clumpy media: particle propagation and gamma-ray emission

Authors: Silvia Celli, Giovanni Morlino, Stefano Gabici, Felix Aharonian

Abstract: 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 p… ▽ More 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. △ Less

Submitted 21 May, 2019; v1 submitted 27 April, 2018; originally announced April 2018.

Comments: 15 pages, 10 figures, accepted for publication in MNRAS

What causes the ionization rates observed in molecular clouds? The role of cosmic ray protons and electrons

Authors: V. H. M. Phan, G. Morlino, S. Gabici

Abstract: Cosmic rays are usually assumed to be the main ionization agent for the interior of molecular clouds where UV and X-ray photons cannot penetrate. Here we test this hypothesis by limiting ourselves to the case of diffuse clouds and assuming that the average cosmic ray spectrum inside the Galaxy is equal to the one at the position of the Sun as measured by Voyager 1 and AMS-02. To calculate the cosm… ▽ More Cosmic rays are usually assumed to be the main ionization agent for the interior of molecular clouds where UV and X-ray photons cannot penetrate. Here we test this hypothesis by limiting ourselves to the case of diffuse clouds and assuming that the average cosmic ray spectrum inside the Galaxy is equal to the one at the position of the Sun as measured by Voyager 1 and AMS-02. To calculate the cosmic ray spectrum inside the clouds, we solve the one-dimensional transport equation taking into account advection, diffusion and energy losses. While outside the cloud particles diffuse, in its interior they are assumed to gyrate along magnetic field lines because ion-neutral friction is effective in damping all the magnetic turbulence. We show that ionization losses effectively reduce the CR flux in the cloud interior for energies below $\approx 100$ MeV, especially for electrons, in such a way that the ionization rate decreases by roughly 2 order of magnitude with respect to the case where losses are neglected. As a consequence, the predicted ionization rate is more than 10 times smaller than the one inferred from the detection of molecular lines. We discuss the implication of our finding in terms of spatial fluctuation of the Galactic cosmic ray spectra and possible additional sources of low energy cosmic rays. △ Less

Submitted 26 April, 2018; originally announced April 2018.

Comments: 8 pages, 7 figures, submitted

Indication of anisotropy in arrival directions of ultra-high-energy cosmic rays through comparison to the flux pattern of extragalactic gamma-ray sources

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, N. Arsene, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Balaceanu, F. Barbato, R. J. Barreira Luz, J. J. Beatty, K. H. Becker, J. A. Bellido , et al. (368 additional authors not shown)

Abstract: A new analysis of the dataset from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above 20 EeV with zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been created for t… ▽ More A new analysis of the dataset from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above 20 EeV with zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL) and starburst galaxies from a sample that was examined with Fermi-LAT. Flux-limited samples, which include all types of galaxies from the Swift-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0 sigma, the highest value of the test statistic being for energies above 39 EeV. The three alternative models are favored against isotropy with 2.7-3.2 sigma significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed. △ Less

Submitted 6 February, 2018; v1 submitted 18 January, 2018; originally announced January 2018.

Comments: Published version, 10 pages, 3 figures, 2 tables

Report number: FERMILAB-PUB-17-358

Journal ref: Aab. et al., The Astrophysical Journal Letters, 853:L29 (2018)

Diffuse γ-ray emission from self-confined cosmic rays around Galactic sources

Authors: M. D'Angelo, G. Morlino, E. Amato, P. Blasi

Abstract: The propagation of particles accelerated at supernova remnant shocks and escaping the parent remnants is likely to proceed in a strongly non-linear regime, due to the efficient self-generation of Alfvén waves excited through streaming instability near the sources. Depending on the amount of neutral hydrogen present in the regions around the sites of supernova explosions, cosmic rays may accumulate… ▽ More The propagation of particles accelerated at supernova remnant shocks and escaping the parent remnants is likely to proceed in a strongly non-linear regime, due to the efficient self-generation of Alfvén waves excited through streaming instability near the sources. Depending on the amount of neutral hydrogen present in the regions around the sites of supernova explosions, cosmic rays may accumulate an appreciable grammage in the same regions and get self-confined for non-negligible times, which in turn results in an enhanced rate of production of secondaries. Here we calculate the contribution to the diffuse gamma-ray background due to the overlap along lines of sight of several of these extended halos as due to pion production induced by self-confined cosmic rays. We find that if the density of neutrals is low, the halos can account for a substantial fraction of the diffuse emission observed by Fermi-LAT, depending on the orientation of the line of sight with respect to the direction of the Galactic centre. △ Less

Submitted 30 October, 2017; originally announced October 2017.

Comments: 12 pages, 8 figures. Accepted for publication by MNRAS