Showing 1–50 of 60 results for author: Cristofari, P

SPIRou monitoring of the protostar V347 Aur: binarity, magnetic fields, pulsed dynamo and accretion

Authors: J. -F. Donati, P. I. Cristofari, A. Carmona, K. Grankin

Abstract: We present in this paper an analysis of near-infrared observations of the 0.3 Msun protostar V347 Aur collected with the SPIRou high-resolution spectropolarimeter and velocimeter at the 3.6-m Canada-France-Hawaii Telescope from October 2019 to April 2023. From a set of 79 unpolarized and circularly polarized spectra of V347 Aur to which we applied Least-Squares Deconvolution (LSD), we derived radi… ▽ More We present in this paper an analysis of near-infrared observations of the 0.3 Msun protostar V347 Aur collected with the SPIRou high-resolution spectropolarimeter and velocimeter at the 3.6-m Canada-France-Hawaii Telescope from October 2019 to April 2023. From a set of 79 unpolarized and circularly polarized spectra of V347 Aur to which we applied Least-Squares Deconvolution (LSD), we derived radial velocities and longitudinal fields, along with their temporal variations over our monitoring campaign of 1258 d. Our data show that V347 Aur is an eccentric binary system with an orbital period 154.6$\pm$0.7 d, experiencing strong to extreme accretion events near periastron. The companion is a 29.0$\pm$1.6 Mjup brown dwarf, a rare member in the brown dwarf desert of close companions around M dwarfs. We detect weak longitudinal fields (<100~G) at the surface of V347 Aur, significantly weaker than those of more evolved prototypical T Tauri stars. These fields show small-amplitude rotational modulation, indicating a mainly axisymmetric parent large-scale magnetic topology, and larger fluctuations at half the orbital period, suggesting that a pulsed dynamo triggered by orbital motion operates in V347 Aur. Applying Zeeman-Doppler imaging to our circularly polarised LSD profiles, we find that the large-scale field of V347 Aur is mainly toroidal for most of our observations, with the toroidal component switching sign near periastron and apoastron. The weak large-scale dipole (~30 G) is not able to disrupt the disc beyond 1.3 Rstar even outside accretion episodes, implying longitudinally distributed (rather than localized) accretion at the surface of the protostar. △ Less

Submitted 6 July, 2024; originally announced July 2024.

Comments: submitted to MNRAS (18 pages, 12 figures, 5 tables)

Rotation and H$α$ emission in a young SMC cluster: a spectroscopic view of NGC 330

Authors: Paul I. Cristofari, Andrea K. Dupree, Antonino P. Milone, Matthew G. Walker, Mario Mateo, Aaron Dotter, John I. Bailey III

Abstract: We present an analysis of high-resolution optical spectra recorded for 30 stars of the split extended main-sequence turnoff (eMSTO) of the young ($\sim$ 40 Myr) Small Magellanic Cloud (SMC) globular cluster NGC 330. Spectra were obtained with the M2FS and MIKE spectrographs located on the Magellan-Clay 6.5m telescope. These spectra revealed the presence of Be stars, occupying primarily the cool si… ▽ More We present an analysis of high-resolution optical spectra recorded for 30 stars of the split extended main-sequence turnoff (eMSTO) of the young ($\sim$ 40 Myr) Small Magellanic Cloud (SMC) globular cluster NGC 330. Spectra were obtained with the M2FS and MIKE spectrographs located on the Magellan-Clay 6.5m telescope. These spectra revealed the presence of Be stars, occupying primarily the cool side of the split main sequence (MS). Rotational velocity ($v\sin{i}$) measurements for most of the targets are consistent with the presence of two populations of stars in the cluster: one made up of rapidly rotating Be stars ($<v\sin{i}> \approx 200$ $\rm km\,s^{-1}$), and {the other} consisting of warmer stars with slower rotation ($<\!v\sin{i}\!>\approx50$ $\rm km\,s^{-1}$). Core emission in the H$δ$ photospheric lines was observed for most of the H$α$ emitters. The shell parameter computed for the targets in our sample indicate that most of the observed stars should have inclinations below 75$^{\circ}$. These results confirm the detection of Be stars obtained through photometry, but also reveal the presence of narrow H$α$ and H$δ$ features for some targets that cannot be detected with low-resolution spectroscopy or photometry. Asymmetry variability of H$α$ line profiles on the timescales of a few years is also observed, and could provide information on the geometry of the decretion disks. Observations revealed the presence of nebular H$α$ emission, strong enough in faint targets to compromise the extraction of spectra and to impact narrow band photometry used to assess the presence of H$α$ emission. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: 15 pages, 13 figures, accepted for publication in The Astrophysical Journal

SPIRou spectropolarimetry of the T Tauri star TW Hydrae: magnetic fields, accretion and planets

Authors: J. -F. Donati, P. I. Cristofari, L. T. Lehmann, C. Moutou, S. H. P. Alencar, J. Bouvier, L. Arnold, X. Delfosse, E. Artigau, N. Cook, Á. Kóspál, F. Ménard, C. Baruteau, M. Takami, S. Cabrit, G. Hébrard, R. Doyon, the SPIRou science team

Abstract: In this paper we report near-infrared observations of the classical T Tauri star TW Hya with the SPIRou high-resolution spectropolarimeter and velocimeter at the 3.6-m Canada-France-Hawaii Telescope in 2019, 2020, 2021 and 2022. By applying Least-Squares Deconvolution (LSD) to our circularly polarized spectra, we derived longitudinal fields that vary from year to year from -200 to +100 G, and exhi… ▽ More In this paper we report near-infrared observations of the classical T Tauri star TW Hya with the SPIRou high-resolution spectropolarimeter and velocimeter at the 3.6-m Canada-France-Hawaii Telescope in 2019, 2020, 2021 and 2022. By applying Least-Squares Deconvolution (LSD) to our circularly polarized spectra, we derived longitudinal fields that vary from year to year from -200 to +100 G, and exhibit low-level modulation on the 3.6 d rotation period of TW Hya, despite the star being viewed almost pole-on. We then used Zeeman-Doppler Imaging to invert our sets of unpolarized and circularly-polarized LSD profiles into brightness and magnetic maps of TW Hya in all 4 seasons, and obtain that the large-scale field of this T Tauri star mainly consists of a 1.0-1.2 kG dipole tilted at about 20° to the rotation axis, whereas the small-scale field reaches strengths of up to 3-4 kG. We find that the large-scale field is strong enough to allow TW Hya to accrete material from the disc on the polar regions at the stellar surface in a more or less geometrically stable accretion pattern, but not to succeed in spinning down the star. We also report the discovery of a radial velocity signal of semi-amplitude $11.1^{+3.3}_{-2.6}$ m/s (detected at 4.3$σ$ at a period of 8.3 d in the spectrum of TW Hya, whose origin may be attributed to either a non-axisymmetric density structure in the inner accretion disc, or to a $0.55^{+0.17}_{-0.13}$ Jupiter mass candidate close-in planet (if orbiting in the disc plane), at an orbital distance of $0.075\pm0.001$ au. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: MNRAS, in press (23 pages, 16 figures, 6 tables)

The population of Galactic supernova remnants in the TeV range

Authors: Rowan Batzofin, Pierre Cristofari, Kathrin Egberts, Constantin Steppa, Dominique Meyer

Abstract: Supernova remnants (SNRs) are likely to be significant sources of cosmic rays up to the knee of the local cosmic-ray (CR) spectrum. They produce gamma-rays in the very-high-energy (VHE) ($E>0.1$ TeV) range via: hadronic interactions with the interstellar medium and leptonic interactions with soft photons. Current observations have lead to the detection of about a dozen of VHE SNRs and future instr… ▽ More Supernova remnants (SNRs) are likely to be significant sources of cosmic rays up to the knee of the local cosmic-ray (CR) spectrum. They produce gamma-rays in the very-high-energy (VHE) ($E>0.1$ TeV) range via: hadronic interactions with the interstellar medium and leptonic interactions with soft photons. Current observations have lead to the detection of about a dozen of VHE SNRs and future instruments should increase this number. The details of particle acceleration at SNRs, and of the mechanisms producing VHE gamma-rays at SNRs are poorly understood. We aim to study the population of SNRs detected in the TeV range and its properties, and to address fundamental questions of particle acceleration at SNR shocks: What is the spectrum of accelerated particles? What is the efficiency of acceleration? Is the VHE emission dominated by hadronic or leptonic interactions? By means of Monte Carlo methods, we simulate the population of SNRs in the VHE domain and confront our simulations to H.E.S.S. Galactic Plane Survey (HGPS). We explore the parameter space: the slope of accelerated particles $α$, the electron-to-proton ratio $K_{\rm ep}$, and the efficiency of particle acceleration $ξ$. We found sets of parameters for which $\gtrsim 90$\% of realisations are found in agreements with the HGPS data. These parameters are found $ 4.2 \gtrsim α\gtrsim 4.1 $, $10^{-5} \lesssim K_{\rm ep} \lesssim 10^{-4.5}$, and $0.03 \lesssim ξ\lesssim 0.1 $ . We were able to strongly argue against some regions of the parameter space: $α\lesssim 4.05$, $α\gtrsim 4.35$, or $K_{\rm ep} \gtrsim 10^{-3}$. Our model is so far able to explain the SNR population of the HGPS. Our approach, confronted to the results of future systematic surveys, will help remove degeneracy in the solutions, and to better understand particle acceleration at SNR shocks. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 15 pages, 12 figures. Accepted for publication in A&A

Spectrum and extension of the inverse-Compton emission of the Crab Nebula from a combined Fermi-LAT and H.E.S.S. analysis

Authors: F. Aharonian, F. Ait Benkhali, J. Aschersleben, H. Ashkar, M. Backes, A. Baktash, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, J. Borowska, F. Bradascio, M. Breuhaus, R. Brose, A. Brown, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin , et al. (137 additional authors not shown)

Abstract: The Crab Nebula is a unique laboratory for studying the acceleration of electrons and positrons through their non-thermal radiation. Observations of very-high-energy $γ$ rays from the Crab Nebula have provided important constraints for modelling its broadband emission. We present the first fully self-consistent analysis of the Crab Nebula's $γ$-ray emission between 1 GeV and $\sim$100 TeV, that is… ▽ More The Crab Nebula is a unique laboratory for studying the acceleration of electrons and positrons through their non-thermal radiation. Observations of very-high-energy $γ$ rays from the Crab Nebula have provided important constraints for modelling its broadband emission. We present the first fully self-consistent analysis of the Crab Nebula's $γ$-ray emission between 1 GeV and $\sim$100 TeV, that is, over five orders of magnitude in energy. Using the open-source software package Gammapy, we combined 11.4 yr of data from the Fermi Large Area Telescope and 80 h of High Energy Stereoscopic System (H.E.S.S.) data at the event level and provide a measurement of the spatial extension of the nebula and its energy spectrum. We find evidence for a shrinking of the nebula with increasing $γ$-ray energy. Furthermore, we fitted several phenomenological models to the measured data, finding that none of them can fully describe the spatial extension and the spectral energy distribution at the same time. Especially the extension measured at TeV energies appears too large when compared to the X-ray emission. Our measurements probe the structure of the magnetic field between the pulsar wind termination shock and the dust torus, and we conclude that the magnetic field strength decreases with increasing distance from the pulsar. We complement our study with a careful assessment of systematic uncertainties. △ Less

Submitted 21 March, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: 18+6 pages, 15+2 figures. Accepted for publication in A&A. Corresponding authors: M. Meyer, L. Mohrmann, T. Unbehaun. v2: after A&A language editing

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

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, 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, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 43 pages JCAP style (excluding author list and references), 19 figures

The classical T Tauri star CI Tau observed with SPIRou: magnetospheric accretion and planetary formation

Authors: J. -F. Donati, B. Finociety, P. I. Cristofari, S. H. P. Alencar, C. Moutou, X. Delfosse, P Fouqué, L. Arnold, C. Baruteau, Á. Kóspál, F. Ménard, A. Carmona, K. Grankin, M. Takami, E. Artigau, R. Doyon, G. Hébrard, the SLS collaboration

Abstract: We report new observations of the classical T~Tauri star CI~Tau with the SPIRou near-infrared spectropolarimeter and velocimeter at the Canada-France-Hawaii Telescope (CFHT) in late 2019, 2020 and 2022, complemented with observations obtained with the ESPaDOnS optical spectropolarimeter at CFHT in late 2020. From our SPIRou and ESPaDOnS spectra, to which we applied Least-Squares Deconvolution, we… ▽ More We report new observations of the classical T~Tauri star CI~Tau with the SPIRou near-infrared spectropolarimeter and velocimeter at the Canada-France-Hawaii Telescope (CFHT) in late 2019, 2020 and 2022, complemented with observations obtained with the ESPaDOnS optical spectropolarimeter at CFHT in late 2020. From our SPIRou and ESPaDOnS spectra, to which we applied Least-Squares Deconvolution, we infer longitudinal fields clearly modulated with the 9-d rotation period of CI~Tau. Using Zeeman-Doppler imaging, we reconstruct the large-scale magnetic topology, first from SPIRou data only in all three seasons, then from our 2020 SPIRou and ESPaDOnS data simultaneously. We find that CI~Tau hosts a mainly axisymmetric poloidal field, with a 1~kG dipole slightly tilted to the rotation axis and dark spots close to the pole that coincide with the footpoints of accretion funnels linking the star to the inner disc. Our results also suggest that CI~Tau accretes mass from the disc in a stable fashion. We further find that radial velocities (RV) derived from atomic and CO lines in SPIRou spectra are both rotationally modulated, but with a much lower amplitude than that expected from the putative candidate planet CI~Tau~b. We confirm the presence of a RV signal at a period of 23.86~d reported in a separate analysis, but detect it clearly in CO lines only and not in atomic lines, suggesting that it likely traces a non-axisymmetric structure in the inner disc of CI~Tau rather than a massive close-in planet. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: MNRAS, in press (22 pages, 11 figures, 6 tables)

Cosmic-ray induced ionization rates and non-thermal emissions from nuclei of starburst galaxies

Authors: Vo Hong Minh Phan, Enrico Peretti, Pierre Cristofari, Antoine Gusdorf, Philipp Mertsch

Abstract: Cosmic rays are the only agent capable of ionizing the interior of dense molecular clouds and, thus, they are believed to play an essential role in determining the physical and chemical evolution of star-forming regions. In this work, we aim to study cosmic-ray induced ionization rates in starburst environments using non-thermal emissions of cosmic rays from starburst nuclei. To this end, we first… ▽ More Cosmic rays are the only agent capable of ionizing the interior of dense molecular clouds and, thus, they are believed to play an essential role in determining the physical and chemical evolution of star-forming regions. In this work, we aim to study cosmic-ray induced ionization rates in starburst environments using non-thermal emissions of cosmic rays from starburst nuclei. To this end, we first revisit cosmic-ray models which could explain data of non-thermal emissions from radio to X-ray and gamma-ray from nuclei of three prototypical starburst galaxies NGC 253, M82, and Arp 220. These models are then applied to predict ionization rates in starburst environments which gives values around $10^{-14}$ s$^{-1}$. Such a high value of the ionization rate, which is 2 to 3 orders of magnitude higher than the typical values found in the Milky Way, is probably due to relatively high rates of supernova explosions occurring within the nuclei of these starburst galaxies. We also discuss in more details the case of NGC 253 where our predicted ionization rate is found to be, in most cases, a few times smaller than the values inferred from molecular line observations of clouds in the starburst nucleus. The general framework provided in this work illustrates how the use of non-thermal emission data could help to provide more insights into ionization rates or, more generally, cosmic-ray impact in starburst environments. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 12 pages and 5 figures

The galactic bubbles of starburst galaxies The influence of galactic large-scale magnetic fields

Authors: Z. Meliani, P. Cristofari, A. Rodriguez-Gonzalez, G. Fichet de Clairfontaine, E. Proust, E. Peretti

Abstract: Context. The galactic winds of starburst galaxies (SBGs) give rise to remarkable structures on kiloparsec scales. However, the evolution and shape of these giant wind bubbles, as well as the properties of the shocks they develop, are not yet fully understood. Aims. We aim to understand what shapes the galactic winds of SBGs, with a particular focus on the role of large-scale magnetic fields in the… ▽ More Context. The galactic winds of starburst galaxies (SBGs) give rise to remarkable structures on kiloparsec scales. However, the evolution and shape of these giant wind bubbles, as well as the properties of the shocks they develop, are not yet fully understood. Aims. We aim to understand what shapes the galactic winds of SBGs, with a particular focus on the role of large-scale magnetic fields in the dynamical evolution of galactic wind-inflated bubbles. In addition, we aim to explore where the conditions for efficient particle acceleration are met in these systems. Methods. We performed magnetohydrodynamic simulations with the AMRVAC code (Adaptive Mesh Refinement Versatile Advection Code) with various configurations of the galactic medium density profile and magnetization. Results. We observe that the large-scale magnetic field, in which galactic winds expand, can impact the structure and evolution of inflated bubbles. However, the typical structures observed in starburst galaxies, such as M82, cannot be solely explained by the magnetic field structures that have been considered. This highlights the importance of other factors, such as the galactic disk, in shaping the galactic bubble. Furthermore, in all the magnetized cases we investigated, the forward wave resulting from the expanding bubbles only results in compression waves, whereas the wind termination shock features high Mach numbers, making it a promising site for diffusive shock acceleration up to $\sim 10^{2}$ PeV. The synthetic X-ray images generated from our models reveal an envelope surrounding the bubbles that extends up to 2 kpc, which could correspond to the polarized emission observed from planar geometry in M82, as well as a large structure inside the bubble corresponding to the shocked galactic wind. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: 17 pages

Origin of $^{60}$Fe nuclei in cosmic rays: the contribution of local OB associations

Authors: Nicolas de Séréville, Vincent Tatischeff, Pierre Cristofari, Stefano Gabici, Roland Diehl

Abstract: The presence of live $^{60}$Fe nuclei (lifetime of 3.8~Myr) in cosmic rays detected by the ACE/CRIS instrument suggests a nearby nucleosynthesis source. $^{60}$Fe is primarily produced in core-collapse supernovae, and we aim to clarify whether the detected $^{60}$Fe nuclei can be associated with a particular local supernova. We consider 25 OB associations and sub-groups located within 1 kpc of the… ▽ More The presence of live $^{60}$Fe nuclei (lifetime of 3.8~Myr) in cosmic rays detected by the ACE/CRIS instrument suggests a nearby nucleosynthesis source. $^{60}$Fe is primarily produced in core-collapse supernovae, and we aim to clarify whether the detected $^{60}$Fe nuclei can be associated with a particular local supernova. We consider 25 OB associations and sub-groups located within 1 kpc of the solar system based on recent $Gaia$ census. A model is developed that combines stellar population synthesis within these OB associations, cosmic-ray acceleration within associated superbubbles, and cosmic-ray transport to the solar system. The most critical model parameter impacting $^{60}$Fe cosmic-ray production is the explodability criterion, which determines if a massive star ends its life as a supernova. Our study points to the Sco-Cen OB association as the most probable origin of the observed $^{60}$Fe nuclei, particularly suggesting they were accelerated in the Sco-Cen superbubble by a young supernova aged $\leq500$ kyr with a progenitor mass of approximately $13-20~M_\odot$. A less likely source is the supernova at the origin of the Geminga pulsar 342 kyr ago, if the progenitor originated in the Orion OB1 association. The contribution of local OB associations to the cosmic-ray density of stable $^{56}$Fe is estimated to be around 20\%, with some sensitivity to cosmic ray acceleration efficiency and diffusion coefficient. These findings shed light on the origins of cosmic-ray nuclei, connecting them to nucleosynthesis events within our local cosmic neighborhood. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Comprehensive High-resolution Chemical Spectroscopy of Barnard's Star with SPIRou

Authors: Farbod Jahandar, René Doyon, Étienne Artigau, Neil J. Cook, Charles Cadieux, David Lafrenière, Thierry Forveille, Jean-François Donati, Pascal Fouqué, Andrés Carmona, Ryan Cloutier, Paul Cristofari, Eric Gaidos, João Gomes da Silva, Lison Malo, Eder Martioli, J. -D. do Nascimento Jr., Stefan Pelletier, Thomas Vandal, Kim Venn

Abstract: Determination of fundamental parameters of stars impacts all fields of astrophysics, from galaxy evolution to constraining the internal structure of exoplanets. This paper presents a detailed spectroscopic analysis of Barnard's star that compares an exceptionally high-quality (an average signal-to-noise ratio of $\sim$1000 in the entire domain), high-resolution NIR spectrum taken with CFHT/SPIRou… ▽ More Determination of fundamental parameters of stars impacts all fields of astrophysics, from galaxy evolution to constraining the internal structure of exoplanets. This paper presents a detailed spectroscopic analysis of Barnard's star that compares an exceptionally high-quality (an average signal-to-noise ratio of $\sim$1000 in the entire domain), high-resolution NIR spectrum taken with CFHT/SPIRou to PHOENIX-ACES stellar atmosphere models. The observed spectrum shows thousands of lines not identified in the models with a similar large number of lines present in the model but not in the observed data. We also identify several other caveats such as continuum mismatch, unresolved contamination and spectral lines significantly shifted from their expected wavelengths, all of these can be a source of bias for abundance determination. Out of $>10^4$ observed lines in the NIR that could be used for chemical spectroscopy, we identify a short list of a few hundred lines that are reliable. We present a novel method for determining the effective temperature and overall metallicity of slowly-rotating M dwarfs that uses several groups of lines as opposed to bulk spectral fitting methods. With this method, we infer $T_{eff}$ = 3231 $\pm$ 21 K for Barnard's star, consistent with the value of 3238 $\pm$ 11 K inferred from the interferometric method. We also provide abundance measurements of 15 different elements for Barnard's star, including the abundances of four elements (K, O, Y, Th) never reported before for this star. This work emphasizes the need to improve current atmosphere models to fully exploit the NIR domain for chemical spectroscopy analysis. △ Less

Submitted 1 April, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: 29 pages, 23 figures, accepted for publication in ApJ

Measuring small-scale magnetic fields of 44 M dwarfs from SPIRou spectra with ZeeTurbo

Authors: P. I. Cristofari, J. -F. Donati, C. Moutou, L. T. Lehmann, P. Charpentier, P. Fouqué, C. P. Folsom, T. Masseron, A. Carmona, X. Delfosse, P. Petit, E. Artigau, N. J. Cook, the SLS consortium

Abstract: We present the results of an analysis aimed at probing the small-scale magnetic fields of M dwarfs observed with SPIRou, the nIR high-resolution spectro-polarimeter installed at the Canada-France-Hawaii Telescope, in the context of the SPIRou Legacy Survey. Our analysis relies on high-resolution median spectra built from several tens of spectra recorded between 2019 and 2022, and on synthetic spec… ▽ More We present the results of an analysis aimed at probing the small-scale magnetic fields of M dwarfs observed with SPIRou, the nIR high-resolution spectro-polarimeter installed at the Canada-France-Hawaii Telescope, in the context of the SPIRou Legacy Survey. Our analysis relies on high-resolution median spectra built from several tens of spectra recorded between 2019 and 2022, and on synthetic spectra computed with the ZeeTurbo code for various combination of atmospheric parameters and magnetic field strengths. We pursue the efforts undertaken in a previous study and focus on 44 weakly to moderately active M dwarfs. We derive average magnetic field strengths (<$B$>) ranging from 0.05 to 1.15 kG, in good agreement with activity estimates and rotation periods. We found that including magnetic fields in our models has virtually no impact on our derived atmospheric parameters, and that a priori assumptions on the stellar surface gravity can affect our estimated <$B$>. Our results suggest that small-scale magnetic fields account for more than 70% of the overall average magnetic field for most targets whose large-scale fields were previously measured. We derived low magnetic fluxes for several targets in our sample, and found no clear evidence that <$B$> decreases with increasing Rossby number in the unsaturated dynamo regime. We even identified counterexamples (GJ 1289 and GJ 1286) where the small-scale field is unusually strong despite the long rotation period. Along with similar results on the large-scale fields, our findings further suggest that dynamo processes may operate in a non-conventional mode in these strongly magnetic, slowly-rotating stars. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 27 pages, 25 figures, including appendix. Accepted for publication in MNRAS

Chasing Gravitational Waves with the Cherenkov Telescope Array

Authors: Jarred Gershon Green, Alessandro Carosi, Lara Nava, Barbara Patricelli, Fabian Schüssler, Monica Seglar-Arroyo, Cta Consortium, :, Kazuki Abe, Shotaro Abe, Atreya Acharyya, Remi Adam, Arnau Aguasca-Cabot, Ivan Agudo, Jorge Alfaro, Nuria Alvarez-Crespo, Rafael Alves Batista, Jean-Philippe Amans, Elena Amato, Filippo Ambrosino, Ekrem Oguzhan Angüner, Lucio Angelo Antonelli, Carla Aramo, Cornelia Arcaro, Luisa Arrabito , et al. (545 additional authors not shown)

Abstract: The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very… ▽ More The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA. △ Less

Submitted 5 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/30

Monitoring the young planet host V1298 Tau with SPIRou: planetary system and evolving large-scale magnetic field

Authors: B. Finociety, J. -F. Donati, P. I. Cristofari, C. Moutou, C. Cadieux, N. J. Cook, E. Artigau, C. Baruteau, F. Debras, P. Fouqué, J. Bouvier, S. H. P Alencar, X. Delfosse, K. Grankin, A. Carmona, P. Petit, Á. Kóspál, the SLS/SPICE consortium

Abstract: We report results of a spectropolarimetric monitoring of the young Sun-like star V1298~Tau based on data collected with the near-infrared spectropolarimeter SPIRou at the Canada-France-Hawaii Telescope between late 2019 and early 2023. Using Zeeman-Doppler Imaging and the Time-dependent Imaging of Magnetic Stars methods on circularly polarized spectra, we reconstructed the large-scale magnetic top… ▽ More We report results of a spectropolarimetric monitoring of the young Sun-like star V1298~Tau based on data collected with the near-infrared spectropolarimeter SPIRou at the Canada-France-Hawaii Telescope between late 2019 and early 2023. Using Zeeman-Doppler Imaging and the Time-dependent Imaging of Magnetic Stars methods on circularly polarized spectra, we reconstructed the large-scale magnetic topology of the star (and its temporal evolution), found to be mainly poloidal and axisymmetric with an average strength varying from 90 to 170 G over the ~3.5 years of monitoring. The magnetic field features a dipole whose strength evolves from 85 to 245 G, and whose inclination with respect to the stellar rotation axis remains stable until 2023 where we observe a sudden change, suggesting that the field may undergo a polarity reversal, potentially similar to those periodically experienced by the Sun. Our data suggest that the differential rotation shearing the surface of V1298 Tau is about 1.5 times stronger than that of the Sun. When coupling our data with previous photometric results from K2 and TESS and assuming circular orbits for all four planets, we report a $3.9σ$ detection of the radial velocity signature of the outermost planet (e), associated with a most probable mass, density and orbital period of $M_e=0.95^{+0.33}_{-0.24} \ \rm M_{\rm jup}$, $ρ_e=1.66^{+0.61}_{-0.48}$ $\rm g\,cm^{-3}$ and $P_e=53.0039\pm0.0001 \ \rm d$, respectively. For the 3 inner planets, we only derive 99\% confidence upper limits on their mass of $0.44\ \rm M_{\rm jup}$, $0.22\ \rm M_{\rm jup}$ and $0.25\ \rm M_{\rm jup}$, for b, c and d, respectively. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 45 pages, 28 figures, accepted for publication in MNRAS

First light of VLT/HiRISE: High-resolution spectroscopy of young giant exoplanets

Authors: A. Vigan, M. El Morsy, M. Lopez, G. P. P. L. Otten, J. Garcia, J. Costes, E. Muslimov, A. Viret, Y. Charles, G. Zins, G. Murray, A. Costille, J. Paufique, U. Seemann, M. Houllé, H. Anwand-Heerwart, M. Phillips, A. Abinanti, P. Balard, I. Baraffe, J. -A. Benedetti, P. Blanchard, L. Blanco, J. -L. Beuzit, E. Choquet , et al. (24 additional authors not shown)

Abstract: A major endeavor of this decade is the direct characterization of young giant exoplanets at high spectral resolution to determine the composition of their atmosphere and infer their formation processes and evolution. Such a goal represents a major challenge owing to their small angular separation and luminosity contrast with respect to their parent stars. Instead of designing and implementing comp… ▽ More A major endeavor of this decade is the direct characterization of young giant exoplanets at high spectral resolution to determine the composition of their atmosphere and infer their formation processes and evolution. Such a goal represents a major challenge owing to their small angular separation and luminosity contrast with respect to their parent stars. Instead of designing and implementing completely new facilities, it has been proposed to leverage the capabilities of existing instruments that offer either high contrast imaging or high dispersion spectroscopy, by coupling them using optical fibers. In this work we present the implementation and first on-sky results of the HiRISE instrument at the very large telescope (VLT), which combines the exoplanet imager SPHERE with the recently upgraded high resolution spectrograph CRIRES using single-mode fibers. The goal of HiRISE is to enable the characterization of known companions in the $H$ band, at a spectral resolution of the order of $R = λ/Δλ= 100\,000$, in a few hours of observing time. We present the main design choices and the technical implementation of the system, which is constituted of three major parts: the fiber injection module inside of SPHERE, the fiber bundle around the telescope, and the fiber extraction module at the entrance of CRIRES. We also detail the specific calibrations required for HiRISE and the operations of the instrument for science observations. Finally, we detail the performance of the system in terms of astrometry, temporal stability, optical aberrations, and transmission, for which we report a peak value of $\sim$3.9% based on sky measurements in median observing conditions. Finally, we report on the first astrophysical detection of HiRISE to illustrate its potential. △ Less

Submitted 22 November, 2023; v1 submitted 21 September, 2023; originally announced September 2023.

Comments: 17 pages, 15 figures, 3 tables. Accepted in A&A on 24 October 2023

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

The population of Galactic supernova remnants in the TeV range

Authors: Rowan Batzofin, Pierre Cristofari, Kathrin Egberts, Constantin Steppa

Abstract: SNRs are likely to be significant sources of Galactic cosmic rays up to the knee. They produce gamma rays in the very-high-energy (E>100 GeV) range mainly via two mechanisms: hadronic interactions of accelerated protons with the interstellar medium and leptonic interactions of accelerated electrons with soft photons. Observations with current instruments have lead to the detection of about a dozen… ▽ More SNRs are likely to be significant sources of Galactic cosmic rays up to the knee. They produce gamma rays in the very-high-energy (E>100 GeV) range mainly via two mechanisms: hadronic interactions of accelerated protons with the interstellar medium and leptonic interactions of accelerated electrons with soft photons. Observations with current instruments have lead to the detection of about a dozen SNRs in VHE gamma rays and future instruments will help significantly increase this number. Yet, the details of particle acceleration at SNRs, and of the mechanisms producing VHE gamma-ray at SNRs remain poorly understood: What is the spectrum of accelerated particles? What is the efficiency of particle acceleration? Is the gamma-ray emission dominated by hadronic or leptonic origin? To address these questions, we simulate the population of SNRs in the gamma-ray domain, and confront it to the current population of TeV SNRs. This method allows us to investigate several crucial aspects of particle acceleration at SNRs, such as the level of magnetic field around SNR shocks or scanning the parameter space of the accelerated particles (spectral index, electron to proton ratio and the acceleration efficiency of the shock) with the possibility to constrain some of the parameters. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 8 pages, 5 figures, Proceedings for the 38th International Cosmic Ray Conference (ICRC2023)

Magnetic fields & rotation periods of M dwarfs from SPIRou spectra

Authors: J. -F. Donati, L. T. Lehmann, P. I. Cristofari, P. Fouqué, C. Moutou, P. Charpentier, M. Ould-Elhkim, A. Carmona, X. Delfosse, E. Artigau, S. H. P. Alencar, C. Cadieux, L. Arnold, P. Petit, J. Morin, T. Forveille, R. Cloutier, R. Doyon, G. Hébrard, the SLS collaboration

Abstract: We present near-infrared spectropolarimetric observations of a sample of 43 weakly- to moderately-active M dwarfs, carried with SPIRou at the Canada-France-Hawaii Telescope in the framework of the SPIRou Legacy Survey from early 2019 to mid 2022. We use the 6700 circularly polarised spectra collected for this sample to investigate the longitudinal magnetic field and its temporal variations for all… ▽ More We present near-infrared spectropolarimetric observations of a sample of 43 weakly- to moderately-active M dwarfs, carried with SPIRou at the Canada-France-Hawaii Telescope in the framework of the SPIRou Legacy Survey from early 2019 to mid 2022. We use the 6700 circularly polarised spectra collected for this sample to investigate the longitudinal magnetic field and its temporal variations for all sample stars, from which we diagnose, through quasi-periodic Gaussian process regression, the periodic modulation and longer-term fluctuations of the longitudinal field. We detect the large-scale field for 40 of our 43 sample stars, and infer a reliable or tentative rotation period for 38 of them, using a Bayesian framework to diagnose the confidence level at which each rotation period is detected. We find rotation periods ranging from 14 to over 60d for the early-M dwarfs, and from 70 to 200d for most mid- and late-M dwarfs (potentially up to 430d for one of them). We also find that the strength of the detected large-scale fields does not decrease with increasing period or Rossby number for the slowly rotating dwarfs of our sample as it does for higher-mass, more active stars, suggesting that these magnetic fields may be generated through a different dynamo regime than those of more rapidly rotating stars. We also show that the large-scale fields of most sample stars evolve on long timescales, with some of them globally switching sign as stars progress on their putative magnetic cycles. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: MNRAS, in press (25 pages, 15 figures, 3 tables)

$\texttt{Wapiti}$: a data-driven approach to correct for systematics in RV data -- Application to SPIRou data of the planet-hosting M dwarf GJ 251

Authors: M. Ould-Elhkim, C. Moutou, J-F. Donati, É. Artigau, P. Fouqué, N. J. Cook, A. Carmona, P. I. Cristofari, E. Martioli, F. Debras, X. Dumusque, J. H. C. Martins, G. Hébrard, C. Cadieux, X. Delfosse, R. Doyon, B. Klein, J. Gomes da Silva, T. Forveille, T. Hood, P. Charpentier

Abstract: Context: Recent advances in the development of precise radial velocity (RV) instruments in the near-infrared (nIR) domain, such as SPIRou, have facilitated the study of M-type stars to more effectively characterize planetary systems. However, the nIR presents unique challenges in exoplanet detection due to various sources of planet-independent signals which can result in systematic errors in the R… ▽ More Context: Recent advances in the development of precise radial velocity (RV) instruments in the near-infrared (nIR) domain, such as SPIRou, have facilitated the study of M-type stars to more effectively characterize planetary systems. However, the nIR presents unique challenges in exoplanet detection due to various sources of planet-independent signals which can result in systematic errors in the RV data. Aims: In order to address the challenges posed by the detection of exoplanetary systems around M-type stars using nIR observations, we introduce a new data-driven approach for correcting systematic errors in RV data. The effectiveness of this method is demonstrated through its application to the star GJ 251. Methods: Our proposed method, referred to as $\texttt{Wapiti}$ (Weighted principAl comPonent analysIs reconsTructIon), uses a dataset of per-line RV time-series generated by the line-by-line (LBL) algorithm and employs a weighted principal component analysis (wPCA) to reconstruct the original RV time-series. A multi-step process is employed to determine the appropriate number of components, with the ultimate goal of subtracting the wPCA reconstruction of the per-line RV time-series from the original data in order to correct systematic errors. Results: The application of $\texttt{Wapiti}$ to GJ 251 successfully eliminates spurious signals from the RV time-series and enables the first detection in the nIR of GJ 251b, a known temperate super-Earth with an orbital period of 14.2 days. This demonstrates that, even when systematics in SPIRou data are unidentified, it is still possible to effectively address them and fully realize the instrument's capability for exoplanet detection. Additionally, in contrast to the use of optical RVs, this detection did not require to filter out stellar activity, highlighting a key advantage of nIR RV measurements. △ Less

Submitted 4 May, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: Submitted to A&A. For the publicly available Wapiti code, see https://github.com/HkmMerwan/wapiti

Journal ref: A&A 675, A187 (2023)

The magnetic field and multiple planets of the young dwarf AU~Mic

Authors: J. -F. Donati, P. I. Cristofari, B. Finociety, B. Klein, C. Moutou, E. Gaidos, C. Cadieux, E. Artigau, A. C. M. Correia, G. Boué, N. J. Cook, A. Carmona, L. T. Lehmann, J. Bouvier, E. Martioli, J. Morin, P. Fouqué, X. Delfosse, R. Royon, G. Hébrard, S. H. P. Alencar, J. Laskar, L. Arnold, P. Petit, A. Kospal , et al. (3 additional authors not shown)

Abstract: In this paper we present an analysis of near-infrared spectropolarimetric and velocimetric data of the young M dwarf AU Mic, collected with SPIRou at the Canada-France-Hawaii telescope from 2019 to 2022, mostly within the SPIRou Legacy Survey. With these data, we study the large- and small-scale magnetic field of AU Mic, detected through the unpolarized and circularly-polarized Zeeman signatures o… ▽ More In this paper we present an analysis of near-infrared spectropolarimetric and velocimetric data of the young M dwarf AU Mic, collected with SPIRou at the Canada-France-Hawaii telescope from 2019 to 2022, mostly within the SPIRou Legacy Survey. With these data, we study the large- and small-scale magnetic field of AU Mic, detected through the unpolarized and circularly-polarized Zeeman signatures of spectral lines. We find that both are modulated with the stellar rotation period (4.86 d), and evolve on a timescale of months under differential rotation and intrinsic variability. The small-scale field, estimated from the broadening of spectral lines, reaches $2.61\pm0.05$ kG. The large-scale field, inferred with Zeeman-Doppler imaging from Least-Squares Deconvolved profiles of circularly-polarized and unpolarized spectral lines, is mostly poloidal and axisymmetric, with an average intensity of $550\pm30$ G. We also find that surface differential rotation, as derived from the large-scale field, is $\simeq$30% weaker than that of the Sun. We detect the radial velocity (RV) signatures of transiting planets b and c, although dwarfed by activity, and put an upper limit on that of candidate planet d, putatively causing the transit-timing variations of b and c. We also report the detection of the RV signature of a new candidate planet (e) orbiting further out with a period of $33.39\pm0.10$ d, i.e., near the 4:1 resonance with b. The RV signature of e is detected at 6.5$σ$ while those of b and c show up at $\simeq$4$σ$, yielding masses of $10.2^{+3.9}_{-2.7}$ and $14.2^{+4.8}_{-3.5}$ Earth masses for b and c, and a minimum mass of $35.2^{+6.7}_{-5.4}$ Earth masses for e. △ Less

Submitted 24 April, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: MNRAS, in press (20 pages and 12 figures + 9 pages of supplementary material)

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

Constraining atmospheric parameters and surface magnetic fields with $\texttt{ZeeTurbo}$: an application to SPIRou spectra

Authors: P. I. Cristofari, J. -F. Donati, C. P. Folsom, T. Masseron, P. Fouqué, C. Moutou, E. Artigau, A. Carmona, P. Petit, X. Delfosse, E. Martioli

Abstract: We report first results on a method aimed at simultaneously characterising atmospheric parameters and magnetic properties of M dwarfs from high-resolution nIR spectra recorded with SPIRou in the framework of the SPIRou Legacy Survey. Our analysis relies on fitting synthetic spectra computed from MARCS model atmospheres to selected spectral lines, both sensitive and insensitive to magnetic fields.… ▽ More We report first results on a method aimed at simultaneously characterising atmospheric parameters and magnetic properties of M dwarfs from high-resolution nIR spectra recorded with SPIRou in the framework of the SPIRou Legacy Survey. Our analysis relies on fitting synthetic spectra computed from MARCS model atmospheres to selected spectral lines, both sensitive and insensitive to magnetic fields. We introduce a new code, $\texttt{ZeeTurbo}$, obtained by including the Zeeman effect and polarised radiative transfer capabilities to $\texttt{Turbospectrum}$. We compute a grid of synthetic spectra with $\texttt{ZeeTurbo}$ for different magnetic field strengths and develop a process to simultaneously constrain $T_{\rm eff}$, $\log{g}$, [M/H], [$α$/Fe] and the average surface magnetic flux. In this paper, we present our approach and assess its performance using simulations, before applying it to six targets observed in the context of the SPIRou Legacy Survey (SLS), namely AU Mic, EV Lac, AD Leo, CN Leo, PM J18482+0741, and DS Leo. Our method allows us to retrieve atmospheric parameters in good agreement with the literature, and simultaneously yields surface magnetic fluxes in the range 2-4 kG with a typical precision of 0.05 kG, in agreement with literature estimates, and consistent with the saturated dynamo regime in which most of these stars are. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 17 pages plus supplementary material. Accepted for publication in MNRAS

The SPIRou Legacy Survey Rotation period of quiet M dwarfs from circular polarization in near-infrared spectral lines: I. The SPIRou APERO analysis

Authors: P. Fouqué, E. Martioli, J. -F. Donati, L. T. Lehmann, B. Zaire, S. Bellotti, E. Gaidos, J. Morin, C. Moutou, P. Petit, S. H. P. Alencar, L. Arnold, É. Artigau, T. -Q. Cang, A. Carmona, N. J. Cook, P. Cortés-Zuleta, P. I. Cristofari, X. Delfosse, R. Doyon, G. Hébrard, L. Malo, C. Reylé, C. Usher

Abstract: Context. The rotation period of stars is an important parameter along with mass, radius, effective temperature. It is an essential parameter for any radial velocity monitoring, as stellar activity can mimic the presence of a planet at the stellar rotation period. Several methods exist to measure it, including long sequences of photometric measurements or temporal series of stellar activity indicat… ▽ More Context. The rotation period of stars is an important parameter along with mass, radius, effective temperature. It is an essential parameter for any radial velocity monitoring, as stellar activity can mimic the presence of a planet at the stellar rotation period. Several methods exist to measure it, including long sequences of photometric measurements or temporal series of stellar activity indicators. Aims. Here, we use the circular polarization in near-infrared spectral lines for a sample of 43 quiet M dwarfs and compare the measured rotation periods to those obtained with other methods. Methods. From Stokes V spectropolarimetric sequences observed with SPIRou at CFHT and the data processed with the APERO pipeline, we compute the least squares deconvolution profiles using different masks of atomic stellar lines with known Landé factor appropriate to the effective temperature of the star. We derive the longitudinal magnetic field to examine its possible variation along the 50 to 200 observations of each star. For determining the stellar rotation period, we apply a Gaussian process regression enabling us to determine the rotation period of stars with evolving longitudinal field. Results. Among the 43 stars of our sample, we were able to measure a rotation period for 27 stars. For 8 stars, the rotation period was previously unknown. We find a good agreement of our rotation periods with periods found in the literature based on photometry and activity indicators and confirm that near-infrared spectropolarimetry is an important tool to measure rotation periods, even for magnetically quiet stars. Furthermore, we compute ages for 20 stars of our sample using gyrochronology. △ Less

Submitted 8 February, 2023; v1 submitted 7 February, 2023; originally announced February 2023.

Journal ref: A&A 672, A52 (2023)

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

Optical and near-infrared stellar activity characterization of the early M dwarf Gl~205 with SOPHIE and SPIRou

Authors: P. Cortes-Zuleta, I. Boisse, B. Klein, E. Martioli, P. I. Cristofari, A. Antoniadis-Karnavas, J-F. Donati, X. Delfosse, C. Cadieux, N. Heidari, E. Artigau, S. Bellotti, X. Bonfils, A. Carmona, N. J. Cook, R. F. Diaz, R. Doyon, P. Fouque, C. Moutou, P. Petit, T. Vandal, L. Acuña, L. Arnold, N. Astudillo-Defru, V. Bourrier , et al. (19 additional authors not shown)

Abstract: The stellar activity of M dwarfs is the main limitation for discovering and characterizing exoplanets orbiting them since it induces quasi-periodic RV variations. We aim to characterize the magnetic field and stellar activity of the early, moderately active, M dwarf Gl205 in the optical and nIR domains. We obtained high-precision quasi-simultaneous spectra in the optical and nIR with the SOPHIE sp… ▽ More The stellar activity of M dwarfs is the main limitation for discovering and characterizing exoplanets orbiting them since it induces quasi-periodic RV variations. We aim to characterize the magnetic field and stellar activity of the early, moderately active, M dwarf Gl205 in the optical and nIR domains. We obtained high-precision quasi-simultaneous spectra in the optical and nIR with the SOPHIE spectrograph and SPIRou spectropolarimeter between 2019 and 2022. We computed the RVs from both instruments and the SPIRou Stokes V profiles. We used ZDI to map the large-scale magnetic field over the time span of the observations. We studied the temporal behavior of optical and nIR RVs and activity indicators with the Lomb-Scargle periodogram and a quasi-periodic GP regression. In the nIR, we studied the equivalent width of Al I, Ti I, K I, Fe I, and He I. We modeled the activity-induced RV jitter using a multi-dimensional GP regression with activity indicators as ancillary time series. The optical and nIR RVs have similar scatter but nIR shows a more complex temporal evolution. We observe an evolution of the magnetic field topology from a poloidal dipolar field in 2019 to a dominantly toroidal field in 2022. We measured a stellar rotation period of Prot=34.4$\pm$0.5 d in the longitudinal magnetic field. Using ZDI we measure the amount of latitudinal differential rotation (DR) shearing the stellar surface yielding rotation periods of Peq=32.0$\pm$1.8 d at the stellar equator and Ppol=45.5$\pm$0.3 d at the poles. We observed inconsistencies in the activity indicators' periodicities that could be explained by these DR values. The multi-dimensional GP modeling yields an RMS of the RV residuals down to the noise level of 3 m/s for both instruments, using as ancillary time series H$α$ and the BIS in the optical, and the FWHM in the nIR. △ Less

Submitted 22 February, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: 41 pages, 24 figures. Accepted for publication in A&A. Improved quality of figures and reduced size of Appendix

Journal ref: A&A 673, A14 (2023)

A sub-Neptune planet around TOI-1695 discovered and characterized with SPIRou and TESS

Authors: F. Kiefer, G. Hébrard, E. Martioli, E. Artigau, R. Doyon, J. -F. Donati, C. Cadieux, A. Carmona, D. R. Ciardi, P. I. Cristofari, L. de Almeida, P. Figueira, E. Gaidos, E. Gonzales, A. Lecavelier, K. G. Stassun, L. Arnold, B. Benneke, I. Boisse, X. Bonfils, N. J. Cook, P. Cortés-Zuleta, X. Delfosse, J. Dias do Nascimento, M. Fausnaugh , et al. (17 additional authors not shown)

Abstract: TOI-1695 is a V-mag=13 M-dwarf star from the northern hemisphere at 45$\,$pc from the Sun, around which a 3.134-day periodic transit signal from a super-Earth candidate was identified in TESS photometry. With a transit depth of 1.3$\,$mmag, the radius of candidate TOI-1695.01 was estimated by the TESS pipeline to be 1.82$\,$R$_\oplus$ with an equilibrium temperature of $\sim 620\,$K. We successful… ▽ More TOI-1695 is a V-mag=13 M-dwarf star from the northern hemisphere at 45$\,$pc from the Sun, around which a 3.134-day periodic transit signal from a super-Earth candidate was identified in TESS photometry. With a transit depth of 1.3$\,$mmag, the radius of candidate TOI-1695.01 was estimated by the TESS pipeline to be 1.82$\,$R$_\oplus$ with an equilibrium temperature of $\sim 620\,$K. We successfully detect a reflex motion of the star and establish it is due to a planetary companion at an orbital period consistent with the photometric transit period thanks to a year-long radial-velocity monitoring of TOI-1695 by the SPIRou infrared spectropolarimeter. We use and compare different methods to reduce and analyse those data. We report a 5.5-$σ$ detection of the planetary signal, giving a mass of $5.5 \pm 1.0\,$M$_\oplus$ and a radius of $2.03 \pm 0.18\,$R$_\oplus$. We derive a mean equilibrium planet temperature of $590 \pm 90\,$K. The mean density of this small planet of $3.6 \pm 1.1\,$g$\,$cm$^{-3}$ is similar (1.7-$σ$ lower) than that of the Earth. It leads to a non-negligible fraction of volatiles in its atmosphere with $f_{H,He}=0.28^{+0.46}_{-0.23}$% or $f_\text{water}=23 \pm 12$%. TOI-1695$\,$b is a new sub-Neptune planet at the border of the M-dwarf radius valley that can help test formation scenarios for super-Earth/sub-Neptune-like planets. △ Less

Submitted 14 November, 2022; v1 submitted 11 November, 2022; originally announced November 2022.

Comments: 15 pages (+9 appendix pages), 19 figures, 8 tables, accepted for publication in A&A

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

Gamma-ray observations of MAXI J1820+070 during the 2018 outburst

Authors: H. Abe, S. Abe, V. A. Acciari, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari , et al. (418 additional authors not shown)

Abstract: MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS expe… ▽ More MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi-LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi-LAT upper limits, and closer to them for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ~ 500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential high-energy and very-high-energy gamma-ray emitting region should be located at a distance from the black hole ranging between 10^11 and 10^13 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA. △ Less

Submitted 6 October, 2022; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 15 pages, 5 figures. Accepted for publication in MNRAS

Estimating the atmospheric properties of 44 M dwarfs from SPIRou spectra

Authors: P. I. Cristofari, J. -F. Donati, T. Masseron, P. Fouqué, C. Moutou, A. Carmona, E. Artigau, E. Martioli, G. Hébrard, E. Gaidos, X. Delfosse

Abstract: We describe advances on a method designed to derive accurate parameters of M dwarfs. Our analysis consists in comparing high-resolution infrared spectra acquired with the near-infrared spectro-polarimeter SPIRou to synthetic spectra computed from MARCS model atmospheres, in order to derive the effective temperature ($T_{\rm eff}$), surface gravity ($\rm \log{g}$), metallicity ([M/H]) and alpha-enh… ▽ More We describe advances on a method designed to derive accurate parameters of M dwarfs. Our analysis consists in comparing high-resolution infrared spectra acquired with the near-infrared spectro-polarimeter SPIRou to synthetic spectra computed from MARCS model atmospheres, in order to derive the effective temperature ($T_{\rm eff}$), surface gravity ($\rm \log{g}$), metallicity ([M/H]) and alpha-enhancement ($\rm [α/Fe]$) of 44 M dwarfs monitored within the SPIRou Legacy Survey (SLS). Relying on 12 of these stars, we calibrated our method by refining our selection of well modelled stellar lines, and adjusted the line list parameters to improve the fit when necessary. Our retrieved $T_{\rm eff}$, $\rm \log{g}$ and [M/H] are in good agreement with literature values, with dispersions of the order of 50 K in $T_{\rm eff}$ and 0.1 dex in $\rm \log{g}$ and [M/H]. We report that fitting $\rm [α/Fe]$ has an impact on the derivation of the other stellar parameters, motivating us to extend our fitting procedure to this additional parameter. We find that our retrieved $\rm [α/Fe]$ are compatible with those expected from empirical relations derived in other studies. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: 13 pages, 9 figures + appendix and supplementary material; accepted for publication in MNRAS

A deep spectromorphological study of the $γ$-ray emission surrounding the young massive stellar cluster Westerlund 1

Authors: F. Aharonian, H. Ashkar, M. Backes, V. Barbosa Martins, Y. Becherini, D. Berge, B. Bi, M. Böttcher, M. de Bony de Lavergne, F. Bradascio, R. Brose, F. Brun, T. Bulik, C. Burger-Scheidlin, F. Cangemi, S. Caroff, S. Casanova, M. Cerruti, T. Chand, S. Chandra, A. Chen, O. Chibueze, P. Cristofari, J. Damascene Mbarubucyeye, A. Djannati-Ataï , et al. (134 additional authors not shown)

Abstract: Young massive stellar clusters are extreme environments and potentially provide the means for efficient particle acceleration. Indeed, they are increasingly considered as being responsible for a significant fraction of cosmic rays (CRs) accelerated within the Milky Way. Westerlund 1, the most massive known young stellar cluster in our Galaxy is a prime candidate for studying this hypothesis. While… ▽ More Young massive stellar clusters are extreme environments and potentially provide the means for efficient particle acceleration. Indeed, they are increasingly considered as being responsible for a significant fraction of cosmic rays (CRs) accelerated within the Milky Way. Westerlund 1, the most massive known young stellar cluster in our Galaxy is a prime candidate for studying this hypothesis. While the very-high-energy $γ$-ray source HESS J1646-458 has been detected in the vicinity of Westerlund 1 in the past, its association could not be firmly identified. We aim to identify the physical processes responsible for the $γ$-ray emission around Westerlund 1 and thus to better understand the role of massive stellar clusters in the acceleration of Galactic CRs. Using 164 hours of data recorded with the High Energy Stereoscopic System (H.E.S.S.), we carried out a deep spectromorphological study of the $γ$-ray emission of HESS J1646-458. We furthermore employed H I and CO observations of the region to infer the presence of gas that could serve as target material for interactions of accelerated CRs. We detected large-scale ($\sim 2^\circ$ diameter) $γ$-ray emission with a complex morphology, exhibiting a shell-like structure and showing no significant variation with $γ$-ray energy. The combined energy spectrum of the emission extends to several tens of TeV, and is uniform across the entire source region. We did not find a clear correlation of the $γ$-ray emission with gas clouds as identified through H I and CO observations. We conclude that, of the known objects within the region, only Westerlund 1 can explain the bulk of the $γ$-ray emission. Several CR acceleration sites and mechanisms are conceivable, and discussed in detail. (abridged) △ Less

Submitted 10 November, 2022; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: 15 pages, 9 figures. Corresponding authors: L. Mohrmann, S. Ohm, R. Rauth, A. Specovius. v2: corrected affiliation of M. Vecchi (still incorrect in journal version)

Journal ref: A&A 666, A124 (2022)

Microphysics of diffusive shock acceleration: impact on the spectrum of accelerated particles

Authors: Pierre Cristofari, Pasquale Blasi, Damiano Caprioli

Abstract: Diffusive shock acceleration at collisionless shocks remains the most likely process for accelerating particles in a variety of astrophysical sources. While the standard prediction for strong shocks is that the spectrum of accelerated particles is universal, $f(p)\propto p^{-4}$, numerous phenomena affect this simple conclusion. In general, the non-linear dynamical reaction of accelerated particle… ▽ More Diffusive shock acceleration at collisionless shocks remains the most likely process for accelerating particles in a variety of astrophysical sources. While the standard prediction for strong shocks is that the spectrum of accelerated particles is universal, $f(p)\propto p^{-4}$, numerous phenomena affect this simple conclusion. In general, the non-linear dynamical reaction of accelerated particles leads to a concave spectrum, steeper than $p^{-4}$ at momenta below a few tens of GeV/c and harder than the standard prediction at high energies. However, the non-linear effects become important in the presence of magnetic field amplification, which in turn leads to higher values of the maximum momentum $p_{max}$. It was recently discovered that the self-generated perturbations that enhance particle scattering, when advected downstream, move in the same direction as the background plasma, so that the effective compression factor at the shock decreases and the spectrum becomes steeper. We investigate the implications of the excitation of the non-resonant streaming instability on these spectral deformations, the dependence of the spectral steepening on the shock velocity and the role played by the injection momentum. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: 11 pages, 6 figures, accepted for publication in ApJ

TOI-1759 b: a transiting sub-Neptune around a low mass star characterized with SPIRou and TESS

Authors: Eder Martioli, Guillaume Hébrard, Pascal Fouqué, Étienne Artigau, Jean-François Donati, Charles Cadieux, Stefano Bellotti, Alain Lecavelier des Etangs, Réne Doyon, J. -D. do Nascimento Jr., L. Arnold, A. Carmona, N. J. Cook, P. Cortes-Zuleta, L. de Almeida, X. Delfosse, C. P. Folsom, P. -C. König, C. Moutou, M. Ould-Elhkim, P. Petit, K. G. Stassun, A. A. Vidotto, T. Vandal, B. Benneke , et al. (35 additional authors not shown)

Abstract: We report the detection and characterization of the transiting sub-Neptune TOI-1759 b, using photometric time-series from TESS and near infrared spectropolarimetric data from SPIRou on the CFHT. TOI-1759 b orbits a moderately active M0V star with an orbital period of $18.849975\pm0.000006$ d, and we measure a planetary radius and mass of $3.06\pm0.22$ R$_\oplus$ and $6.8\pm2.0$ M$_\oplus$. Radial… ▽ More We report the detection and characterization of the transiting sub-Neptune TOI-1759 b, using photometric time-series from TESS and near infrared spectropolarimetric data from SPIRou on the CFHT. TOI-1759 b orbits a moderately active M0V star with an orbital period of $18.849975\pm0.000006$ d, and we measure a planetary radius and mass of $3.06\pm0.22$ R$_\oplus$ and $6.8\pm2.0$ M$_\oplus$. Radial velocities were extracted from the SPIRou spectra using both the CCF and the LBL methods, optimizing the velocity measurements in the near infrared domain. We analyzed the broadband SED of the star and the high-resolution SPIRou spectra to constrain the stellar parameters and thus improve the accuracy of the derived planet parameters. A LSD analysis of the SPIRou Stokes $V$ polarized spectra detects Zeeman signatures in TOI-1759. We model the rotational modulation of the magnetic stellar activity using a GP regression with a quasi-periodic covariance function, and find a rotation period of $35.65^{+0.17}_{-0.15}$ d. We reconstruct the large-scale surface magnetic field of the star using ZDI, which gives a predominantly poloidal field with a mean strength of $18\pm4$ G. Finally, we perform a joint Bayesian MCMC analysis of the TESS photometry and SPIRou RVs to optimally constrain the system parameters. At $0.1176\pm0.0013$ au from the star, the planet receives $6.4$ times the bolometric flux incident on Earth, and its equilibrium temperature is estimated at $433\pm14$ K. TOI-1759 b is a likely gas-dominated sub-Neptune with an expected high rate of photoevaporation. Therefore, it is an interesting target to search for neutral hydrogen escape, which may provide important constraints on the planetary formation mechanisms responsible for the observed sub-Neptune radius desert. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: Accepted for publication in the 10. Planets and planetary systems section of Astronomy & Astrophysics

Journal ref: A&A 660, A86 (2022)

The first days of type II-P core collapse supernovae in the gamma-ray range

Authors: P. Cristofari, A. Marcowith, M. Renaud, V. V. Dwarkadas, V. Tatischeff, G. Giacinti, E. Peretti, H. Sol

Abstract: Type II-P supernovæ~(SNe), the most common core-collapse SNe type, result from the explosions of red supergiant stars. Their detection in the radio domain testifies of the presence of relativistic electrons, and shows that they are potentially efficient energetic particle accelerators. If hadrons can also be accelerated, these energetic particles are expected to interact with the surrounding mediu… ▽ More Type II-P supernovæ~(SNe), the most common core-collapse SNe type, result from the explosions of red supergiant stars. Their detection in the radio domain testifies of the presence of relativistic electrons, and shows that they are potentially efficient energetic particle accelerators. If hadrons can also be accelerated, these energetic particles are expected to interact with the surrounding medium to produce a gamma-ray signal even in the multi--TeV range. The intensity of this signal depends on various factors, but an essential one is the density of the circumstellar medium. Such a signal should however be limited by electron-positron pair production arising from the interaction of the gamma-ray photons with optical photons emitted by the supernova photosphere, which can potentially degrade the gamma-ray signal by over ten orders of magnitude in the first days/weeks following the explosion. We calculate the gamma-gamma opacity from a detailed modelling of the time evolution of the forward shock and supernova photosphere, taking a full account of the non-isotropy of the photon interactions. We discuss the time-dependent gamma-ray TeV emission from type II-P SNe as a function of the stellar progenitor radius and mass-loss rate, as well as the explosion energy and mass of the ejected material. We evaluate the detectability of the SNe with the next generation of Cherenkov telescopes. We find that, while most extragalactic events may be undetectable, type II-P SNe exploding in our Galaxy or in the Magellanic Clouds should be detected by gamma-ray observatories such as the upcoming Cherenkov Telescope Array. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: 9 pages, 8 figures, Accepted for publication in MNRAS

Evidence for gamma-ray emission from the remnant of Kepler's supernova based on deep H.E.S.S. observations

Authors: H. E. S. S. Collaboration, F. Aharonian, F. Ait Benkhali, E. O. Anguner, H. Ashkar, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernloehr, M. Boettcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun, T. Bulik, T. Bylund, F. Cangemi, S. Caroff, S. Casanova, M. Cerruti, T. Chand , et al. (136 additional authors not shown)

Abstract: Observations with imaging atmospheric Cherenkov telescopes (IACTs) have enhanced our knowledge of nearby supernova (SN) remnants with ages younger than 500 years by establishing Cassiopeia A and the remnant of Tycho's SN as very-high-energy (VHE) gamma-ray sources. The remnant of Kepler's SN, which is the product of the most recent naked-eye supernova in our Galaxy, is comparable in age to the oth… ▽ More Observations with imaging atmospheric Cherenkov telescopes (IACTs) have enhanced our knowledge of nearby supernova (SN) remnants with ages younger than 500 years by establishing Cassiopeia A and the remnant of Tycho's SN as very-high-energy (VHE) gamma-ray sources. The remnant of Kepler's SN, which is the product of the most recent naked-eye supernova in our Galaxy, is comparable in age to the other two, but is significantly more distant. If the gamma-ray luminosities of the remnants of Tycho's and Kepler's SNe are similar, then the latter is expected to be one of the faintest gamma-ray sources within reach of the current generation IACT arrays. Here we report evidence at a statistical level of 4.6 sigma for a VHE signal from the remnant of Kepler's SN based on deep observations by the High Energy Stereoscopic System (H.E.S.S.) with an exposure of 152 hours. The measured integral flux above an energy of 226 GeV is ~0.3% of the flux of the Crab Nebula. The spectral energy distribution (SED) reveals a gamma-ray emitting component connecting the VHE emission observed with H.E.S.S. to the emission observed at GeV energies with Fermi-LAT. The overall SED is similar to that of the remnant of Tycho's SN, possibly indicating the same non-thermal emission processes acting in both these young remnants of thermonuclear SNe. △ Less

Submitted 23 March, 2024; v1 submitted 15 January, 2022; originally announced January 2022.

Comments: Published in Astronomy and Astrophysics. A typo in the normalization constant in Sect. 4 was corrected

Journal ref: A&A 662, A65 (2022)

A MeerKAT, e-MERLIN, H.E.S.S. and Swift search for persistent and transient emission associated with three localised FRBs

Authors: James O. Chibueze, M. Caleb, L. Spitler, H. Ashkar, F. Schussler, B. W. Stappers, C. Venter, I. Heywood, A. M. S. Richards, D. R. A. Williams, M. Kramer, R. Beswick, M. C. Bezuidenhout, R. P. Breton, L. N. Driessen, F. Jankowski, E. F. Keane, M. Malenta, M. Mickaliger, V. Morello, H. Qiu, K. Rajwade, S. Sanidas, M. Surnis, T. W. Scragg , et al. (134 additional authors not shown)

Abstract: We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, an… ▽ More We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst. △ Less

Submitted 31 December, 2021; originally announced January 2022.

Comments: 15 pages, 9 figures

Estimating fundamental parameters of nearby M dwarfs from SPIRou spectra

Authors: P. I. Cristofari, J. -F. Donati, T. Masseron, P. Fouqué, C. Moutou, X. Delfosse, E. Artigau, C. P. Folsom, A. Carmona, E. Gaidos, J. -D. do Nascimento Jr., F. Jahandar, G. Hébrard

Abstract: We present the results of a study aiming at retrieving the fundamental parameters of M dwarfs from spectra secured with SPIRou, the near-infrared high-resolution spectropolarimeter installed at the Canada-France-Hawaii Telescope (CFHT), in the framework of the SPIRou Legacy Survey (SLS). Our study relies on comparing observed spectra with two grids of synthetic spectra, respectively computed from… ▽ More We present the results of a study aiming at retrieving the fundamental parameters of M dwarfs from spectra secured with SPIRou, the near-infrared high-resolution spectropolarimeter installed at the Canada-France-Hawaii Telescope (CFHT), in the framework of the SPIRou Legacy Survey (SLS). Our study relies on comparing observed spectra with two grids of synthetic spectra, respectively computed from PHOENIX and MARCS model atmospheres, with the ultimate goal of optimizing the precision at which fundamental parameters can be determined. In this first step, we applied our technique to 12 inactive M dwarfs with effective temperatures ($T_{\rm eff}$) ranging from 3000 to 4000 K. We implemented a benchmark to carry out a comparison of the two models used in this study. We report that the choice of model has a significant impact on the results and may lead to discrepancies in the derived parameters of 30 K in $T_{\rm eff}$ and 0.05 dex to 0.10 dex in surface gravity ($\log{g}$) and metallicity ([M/H]), as well as systematic shifts of up to 50 K in $T_{\rm eff}$ and 0.4 dex $\log{g}$ and [M/H]. The analysis is performed on high signal-to-noise ratio template SPIRou spectra, averaged over multiple observations corrected from telluric absorption features and sky lines, using both a synthetic telluric transmission model and principal component analysis. With both models, we retrieve $T_{\rm eff}$ , $\log{g}$ and [M/H] estimates in good agreement with reference literature studies, with internal error bars of about 30 K, 0.05 dex and 0.1 dex, respectively. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 20 pages, plus supplementary material, accepted for publication in MNRAS

The Hunt for Pevatrons: The Case of Supernova Remnants

Authors: P. Cristofari

Abstract: The search for Galactic pevatrons is now a well-identified key science project of all instruments operating in the very-high-energy domain. Indeed, in this energy range, the detection of gamma rays clearly indicates that efficient particle acceleration is taking place, and observations can thus help identify which astrophysical sources can energize particles up to the $\sim$PeV range, thus being… ▽ More The search for Galactic pevatrons is now a well-identified key science project of all instruments operating in the very-high-energy domain. Indeed, in this energy range, the detection of gamma rays clearly indicates that efficient particle acceleration is taking place, and observations can thus help identify which astrophysical sources can energize particles up to the $\sim$PeV range, thus being $pevatrons$. In the search for the origin of Galactic cosmic rays (CRs), the PeV range is an important milestone, since the sources of Galactic CRs are expected to accelerate PeV particles. This is how the central scientific goal that is 'solving the mystery of the origin of CRs' has often been distorted into 'finding (a) pevatron(s)'. Since supernova remnants (SNRs) are often cited as the most likely candidates for the origin of CRs, 'finding (a) pevatron(s)' has often become 'confirming that SNRs are pevatrons'. Pleasingly, the first detection(s) of pevatron(s) were not associated to SNRs. Moreover, all clearly detected SNRs have yet revealed to not be pevatrons, and the detection from VHE gamma rays from regions unassociated with SNRs, are reminding us that other astrophysical sites might well be pevatrons. This short review aims at highlighting a few important results on the search for Galactic pevatrons. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: Accepted for publication in the Special Issue "High-Energy Gamma-Ray Astronomy: Results on Fundamental Questions after 30 Years of Ground-Based Observations" Universe MDPI

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

Gamma rays and neutrinos from RX J1713.7-3946 in a lepto-hadronic scenario

Authors: Pierre Cristofari, Viviana Niro, Stefano Gabici

Abstract: The gamma-ray emission of RX J1713.7-3946, despite being extensively studied in the GeV and TeV domain, remains poorly understood. This is mostly because in this range, two competing mechanisms can efficiently produce gamma rays: the inverse Compton scattering of accelerated electrons, and interactions of accelerated protons with nuclei of the ISM. In addition to the acceleration of particles from… ▽ More The gamma-ray emission of RX J1713.7-3946, despite being extensively studied in the GeV and TeV domain, remains poorly understood. This is mostly because in this range, two competing mechanisms can efficiently produce gamma rays: the inverse Compton scattering of accelerated electrons, and interactions of accelerated protons with nuclei of the ISM. In addition to the acceleration of particles from the thermal pool, the reacceleration of pre-existing CRs is often overlooked, and shall in fact also been taken into account. Especially, because of the distance to the SNR ($\sim 1$ kpc), and the low density in which the shock is currently expanding ($\sim 10^{-2}$ cm$^{-3}$), the re-acceleration of CR electrons pre-existing in the ISM, can account for a significant fraction of the observed gamma-ray emission, and contribute to the shaping of the spectrum in the GeV-TeV range. Remarkably, this emission of leptonic origin is found to be close to the level of the gamma-ray signal in the TeV range, provided that the spectrum of pre-exisiting cosmic ray electrons is similar to the one observed in the local interstellar medium. The overall gamma-ray spectrum of RX J1713.7-3946 is naturally produced as the sum of a leptonic emission from reaccelerated CR electrons, and a subdominant hadronic emission from accelerated protons. We also argue that neutrino observations with next-generation detectors might lead to a detection even in the case of a lepto-hadronic origin of the gamma-ray emission. △ Less

Submitted 1 September, 2021; v1 submitted 26 May, 2021; originally announced May 2021.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS

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

Cosmic ray protons and electrons from supernova remnants

Authors: P. Cristofari, P. Blasi, D. Caprioli

Abstract: The spectrum of cosmic ray protons and electrons released by supernova remnants throughout their evolution is poorly known, because of the difficulty in accounting for particle escape and confinement in the downstream of a shock front, where both adiabatic and radiative losses are present. Here we calculate the spectrum of cosmic ray protons released during the evolution of supernovae of different… ▽ More The spectrum of cosmic ray protons and electrons released by supernova remnants throughout their evolution is poorly known, because of the difficulty in accounting for particle escape and confinement in the downstream of a shock front, where both adiabatic and radiative losses are present. Here we calculate the spectrum of cosmic ray protons released during the evolution of supernovae of different types, accounting for the escape from upstream and for adiabatic losses of particles advected downstream of the shock and liberated at later times. The same calculation is carried out for electrons. The magnetic field in the post-shock region is calculated by using an analytic treatment of the magnetic field amplification due to non--resonant and resonant streaming instability and their saturation. We find that when the field is the result of the growth of the cosmic-ray--driven non--resonant instability alone, the spectrum of electrons and protons released by a supernova remnant are indeed different, but such a difference becomes appreciable only at energies $\gtrsim 100-1000$ GeV, while observations of the electron spectrum require such a difference to be present at energies as low as $\sim 10$ GeV. An effect at such low energies requires substantial magnetic field amplification in the late stages of the supernova remnant evolution (shock velocity $\ll 1000$ km/s), perhaps not due to streaming instability but hydrodynamical processes. We comment on the feasibility of such conditions and speculate on the possibility that the difference in spectral shape between electrons and protons may reflect either some unknown acceleration effect, or additional energy losses in cocoons around the sources. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: 11 pages, 5 figures, Accepted for publication in Astronomy \& Astrophysics

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

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

Authors: The Cherenkov Telescope Array Consortium, :, H. Abdalla, H. Abe, F. Acero, A. Acharyya, R. Adam, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves B, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, A. Araudo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, Y. Ascasíbar, M. Ashley , et al. (474 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nucle… ▽ More The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology. △ Less

Submitted 26 February, 2021; v1 submitted 3 October, 2020; originally announced October 2020.

Comments: 71 pages (including affiliations and references), 13 figures, 6 tables. Accepted in JCAP; matches published version. Corresponding authors: Jonathan Biteau, Julien Lefaucheur, Humberto Martinez-Huerta, Manuel Meyer, Santiago Pita, Ievgen Vovk

Journal ref: JCAP 02 (2021) 048

Direct characterization of young giant exoplanets at high spectral resolution by coupling SPHERE and CRIRES+

Authors: G. P. P. L. Otten, A. Vigan, E. Muslimov, M. N'Diaye, E. Choquet, U. Seemann, K. Dohlen, M. Houllé, P. Cristofari, M. W. Phillips, Y. Charles, I. Baraffe, J. -L. Beuzit, A. Costille, R. Dorn, M. El Morsy, M. Kasper, M. Lopez, C. Mordasini, R. Pourcelot, A. Reiners, J. -F. Sauvage

Abstract: Studies of atmospheres of directly imaged exoplanets with high-resolution spectrographs have shown that their characterization is predominantly limited by noise on the stellar halo at the location of the studied exoplanet. An instrumental combination of high-contrast imaging and high spectral resolution that suppresses this noise and resolves the spectral lines can therefore yield higher quality s… ▽ More Studies of atmospheres of directly imaged exoplanets with high-resolution spectrographs have shown that their characterization is predominantly limited by noise on the stellar halo at the location of the studied exoplanet. An instrumental combination of high-contrast imaging and high spectral resolution that suppresses this noise and resolves the spectral lines can therefore yield higher quality spectra. We study the performance of the proposed HiRISE fiber coupling between the SPHERE and CRIRES+ at the VLT for spectral characterization of directly imaged planets. Using end-to-end simulations of HiRISE we determine the S/N of the detection of molecular species for known exoplanets in $H$ and $K$ bands, and compare them to CRIRES+. We investigate the ultimate detection limits of HiRISE as a function of stellar magnitude, and we quantify the impact of different coronagraphs and of the system transmission. We find that HiRISE largely outperforms CRIRES+ for companions around bright hosts like $β$ Pic or 51 Eri. For an $H=3.5$ host, we observe a gain of a factor of up to 16 in observing time with HiRISE to reach the same S/N on a companion at 200 mas. More generally, HiRISE provides better performance than CRIRES+ in two-hour integration times between 50-350 mas for hosts with $H<8.5$ and between 50-700 mas for $H<7$. For fainter hosts like PDS 70 and HIP 65426, no significant improvements are observed. We find that using no coronagraph yields the best S/N when characterizing known exoplanets due to higher transmission and fiber-based starlight suppression. We demonstrate that the overall transmission of the system is in fact the main driver of performance. Finally, we show that HiRISE outperforms the best detection limits of SPHERE for bright stars, opening major possibilities for the characterization of future planetary companions detected by other techniques. △ Less

Submitted 10 December, 2020; v1 submitted 3 September, 2020; originally announced September 2020.

Comments: 21 pages, 17 figures. 2 tables. Originally accepted for publication in A&A. Revised version after error in sky background was discovered (too high) and corrected in proofing stage. Revisions made in Abstract, Sections 3-6 and Figures 6-12,C1,C2 that address the changes in performance. Resubmitted to A&A

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

Authors: The Cherenkov Telescope Array Consortium, :, A. Acharyya, R. Adam, C. Adams, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, L. Amati, G. Ambrosi, E. O. Angüner, L. A. Antonelli, C. Aramo, A. Araudo, T. Armstrong, F. Arqueros, K. Asano, Y. Ascasíbar, M. Ashley, C. Balazs, O. Ballester , et al. (427 additional authors not shown)

Abstract: We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models giv… ▽ More We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]" △ Less

Submitted 30 January, 2021; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 68 pages (including references) and 26 figures; text identical to the version published in JCAP

Journal ref: JCAP01(2021)057

The low rate of Galactic pevatrons

Authors: P. Cristofari, P. Blasi, E. Amato

Abstract: Although supernova remnants remain the main suspects as sources of Galactic cosmic rays up to the knee, the supernova paradigm still has many loose ends. The weakest point in this construction is the possibility that individual supernova remnants can accelerate particles to the rigidity of the knee, $\sim 10^{6}$ GV. This scenario heavily relies upon the possibility to excite current driven non-re… ▽ More Although supernova remnants remain the main suspects as sources of Galactic cosmic rays up to the knee, the supernova paradigm still has many loose ends. The weakest point in this construction is the possibility that individual supernova remnants can accelerate particles to the rigidity of the knee, $\sim 10^{6}$ GV. This scenario heavily relies upon the possibility to excite current driven non-resonant hybrid modes while the shock is still at the beginning of the Sedov phase. These modes can enhance the rate of particle scattering thereby leading to potentially very-high maximum energies. Here we calculate the spectrum of particles released into the interstellar medium from the remnants of different types of supernovae. We find that only the remnants of very powerful, rare core-collapse supernova explosions can accelerate light elements such as hydrogen and helium nuclei, to the knee rigidity, and that the local spectrum of cosmic rays directly constrains the rate of such events, if they are also source of PeV cosmic rays. On the other hand, for remnants of typical core-collapse supernova explosions, the Sedov phase is reached at late times, when the maximum energy is too low and the spectrum at very-high energies is very steep, being mostly produced during the ejecta dominated phase. For typical thermonuclear explosions, resulting in type Ia supernovae, we confirm previous findings that these objects can only produce cosmic rays up to $\lesssim 10^{5}$ GeV. The implications for the overall cosmic ray spectrum observed at the Earth and for the detection of PeVatrons by future gamma-ray observatories are discussed. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: 17 pages, 6 Figures, accepted to Astroparticle Physics

Time-dependent high-energy gamma-ray signal from accelerated particles in core-collapse supernovae: the case of SN 1993J

Authors: Pierre Cristofari, Matthieu Renaud, Alexandre Marcowith, Vikram V. Dwarkadas, Vincent Tatischeff

Abstract: Some core-collapse supernovae are likely to be efficient cosmic-ray accelerators up to the PeV range, and therefore, to potentially play an important role in the overall Galactic cosmic-ray population. The TeV gamma-ray domain can be used to study particle acceleration in the multi-TeV and PeV range. This motivates the study of the detectability of such supernovae by current and future gamma-ray f… ▽ More Some core-collapse supernovae are likely to be efficient cosmic-ray accelerators up to the PeV range, and therefore, to potentially play an important role in the overall Galactic cosmic-ray population. The TeV gamma-ray domain can be used to study particle acceleration in the multi-TeV and PeV range. This motivates the study of the detectability of such supernovae by current and future gamma-ray facilities. The gamma-ray emission of core-collapse supernovae strongly depends on the level of the two-photon annihilation process: high-energy gamma-ray photons emitted at the expanding shock wave following the supernova explosion can interact with soft photons from the supernova photosphere through the pair production channel, thereby strongly suppressing the flux of gamma rays leaving the system. In the case of SN 1993J, whose photospheric and shock-related parameters are well measured, we calculate the temporal evolution of the expected gamma-ray attenuation by accounting for the temporal and geometrical effects. We find the attenuation to be of about $10$ orders of magnitude in the first few days after the SN explosion. The probability of detection of a supernova similar to SN 1993J with the Cherenkov Telescope Array is highest if observations are performed either earlier than 1 day, or later than 10 days after the explosion, when the gamma-ray attenuation decreases to about $2$ orders of magnitude. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: Accepted to MNRAS, 7 pages, 7 figures

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

Gamma rays from reaccelerated particles at supernova remnant shocks

Authors: Pierre Cristofari, Pasquale Blasi

Abstract: Diffusive shock acceleration is considered as the main mechanism for particle energization in supernova remnants, as well as in other classes of sources. The existence of some remnants that show a bilateral morphology in the X-rays and gamma rays suggests that this process occurs with an efficiency that depends upon the inclination angle between the shock normal and the large scale magnetic field… ▽ More Diffusive shock acceleration is considered as the main mechanism for particle energization in supernova remnants, as well as in other classes of sources. The existence of some remnants that show a bilateral morphology in the X-rays and gamma rays suggests that this process occurs with an efficiency that depends upon the inclination angle between the shock normal and the large scale magnetic field in which the shock propagates. This interpretation is additionally supported by recent particle-in-cell simulations that show how ions are not injected if the shock is more oblique than $\sim 45^{o}$. These shocks provide an excellent test bench for the process of reacceleration at the same shock: non-thermal seed particles that are reached by the shock front are automatically injected and accelerated. This process was recently discussed as a possible reason for some anomalous behaviour of the spectra of secondary cosmic ray nuclei. Here we discuss how gamma--ray observations of selected supernova remnants can provide us with precious information about this process and lead us to a better assessment of particle diffusive shock reacceleration for other observables in cosmic ray physics. △ Less

Submitted 31 July, 2019; originally announced July 2019.

Comments: 9 pages, 6 figures, Accepted for publication in MNRAS

The Southern Wide-Field Gamma-Ray Observatory (SWGO): A Next-Generation Ground-Based Survey Instrument for VHE Gamma-Ray Astronomy

Authors: P. Abreu, A. Albert, R. Alfaro, C. Alvarez, R. Arceo, P. Assis, F. Barao, J. Bazo, J. F. Beacom, J. Bellido, S. BenZvi, T. Bretz, C. Brisbois, A. M. Brown, F. Brun, M. Buscemi, K. S. Caballero Mora, P. Camarri, A. Carramiñana, S. Casanova, A. Chiavassa, R. Conceição, G. Cotter, P. Cristofari, S. Dasso , et al. (73 additional authors not shown)

Abstract: We describe plans for the development of the Southern Wide-field Gamma-ray Observatory (SWGO), a next-generation instrument with sensitivity to the very-high-energy (VHE) band to be constructed in the Southern Hemisphere. SWGO will provide wide-field coverage of a large portion of the southern sky, effectively complementing current and future instruments in the global multi-messenger effort to und… ▽ More We describe plans for the development of the Southern Wide-field Gamma-ray Observatory (SWGO), a next-generation instrument with sensitivity to the very-high-energy (VHE) band to be constructed in the Southern Hemisphere. SWGO will provide wide-field coverage of a large portion of the southern sky, effectively complementing current and future instruments in the global multi-messenger effort to understand extreme astrophysical phenomena throughout the universe. A detailed description of science topics addressed by SWGO is available in the science case white paper [1]. The development of SWGO will draw on extensive experience within the community in designing, constructing, and successfully operating wide-field instruments using observations of extensive air showers. The detector will consist of a compact inner array of particle detection units surrounded by a sparser outer array. A key advantage of the design of SWGO is that it can be constructed using current, already proven technology. We estimate a construction cost of 54M USD and a cost of 7.5M USD for 5 years of operation, with an anticipated US contribution of 20M USD ensuring that the US will be a driving force for the SWGO effort. The recently formed SWGO collaboration will conduct site selection and detector optimization studies prior to construction, with full operations foreseen to begin in 2026. Throughout this document, references to science white papers submitted to the Astro2020 Decadal Survey with particular relevance to the key science goals of SWGO, which include unveiling Galactic particle accelerators [2-10], exploring the dynamic universe [11-21], and probing physics beyond the Standard Model [22-25], are highlighted in red boldface. △ Less

Submitted 17 July, 2019; originally announced July 2019.

Comments: Astro2020 APC White Paper

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