Showing 1–50 of 136 results for author: Acero, F

Best of both worlds: Fusing hyperspectral data from two generations of spectro-imagers for X-ray astrophysics

Authors: Julia Lascar, Jérôme Bobin, Fabio Acero

Abstract: With the launch of the X-Ray Imaging and Spectroscopy Mission (XRISM) and the advent of microcalorimeter detectors, X-ray astrophysics is entering in a new era of spatially resolved high resolution spectroscopy. But while this new generation of X-ray telescopes have much finer spectral resolutions than their predecessors (e.g. XMM-Newton, Chandra), they also have coarser spatial resolutions, leadi… ▽ More With the launch of the X-Ray Imaging and Spectroscopy Mission (XRISM) and the advent of microcalorimeter detectors, X-ray astrophysics is entering in a new era of spatially resolved high resolution spectroscopy. But while this new generation of X-ray telescopes have much finer spectral resolutions than their predecessors (e.g. XMM-Newton, Chandra), they also have coarser spatial resolutions, leading to problematic cross-pixel contamination. This issue is currently a critical limitation for the study of extended sources such as galaxy clusters of supernova remnants. To increase the scientific output of XRISM's hyperspectral data, we propose to fuse it with XMM-Newton data, and seek to obtain a cube with the best spatial and spectral resolution of both generations. This is the aim of hyperspectral fusion. In this article, we have implemented an algorithm that jointly deconvolves the spatial response of XRISM and the spectral response of XMM-Newton. To do so, we construct a forward model adapted for instrumental systematic degradations and Poisson noise, then tackle hyperspectral fusion as a regularized inverse problem. We test three methods of regularization: low rank approximation with Sobolev regularization; low rank approximation with 2D wavelet sparsity ; and 2D-1D wavelet sparsity. We test our method on toy models constructed from hydrodynamic simulations of supernova remnants. We find that our method reconstructs the ground truth well even when the toy model is complex. For the regularization term, we find that while the low rank approximation worked well as a spectral denoiser in models with less spectral variability, it introduced a bias in models with more spectral variability, in which case the 2D-1D wavelet sparsity regularization worked best. After demonstrating a proof of concept in this article, we aim to apply this method to real X-ray astrophysical data in the near future. △ Less

Submitted 22 July, 2024; originally announced July 2024.

Comments: This is a pre-print of an article currently in the process of peer review. The code described in this article is available at: https://github.com/JMLascar/HIFReD_fusion

Are Large Language Models Strategic Decision Makers? A Study of Performance and Bias in Two-Player Non-Zero-Sum Games

Authors: Nathan Herr, Fernando Acero, Roberta Raileanu, María Pérez-Ortiz, Zhibin Li

Abstract: Large Language Models (LLMs) have been increasingly used in real-world settings, yet their strategic abilities remain largely unexplored. Game theory provides a good framework for assessing the decision-making abilities of LLMs in interactions with other agents. Although prior studies have shown that LLMs can solve these tasks with carefully curated prompts, they fail when the problem setting or p… ▽ More Large Language Models (LLMs) have been increasingly used in real-world settings, yet their strategic abilities remain largely unexplored. Game theory provides a good framework for assessing the decision-making abilities of LLMs in interactions with other agents. Although prior studies have shown that LLMs can solve these tasks with carefully curated prompts, they fail when the problem setting or prompt changes. In this work we investigate LLMs' behaviour in strategic games, Stag Hunt and Prisoner Dilemma, analyzing performance variations under different settings and prompts. Our results show that the tested state-of-the-art LLMs exhibit at least one of the following systematic biases: (1) positional bias, (2) payoff bias, or (3) behavioural bias. Subsequently, we observed that the LLMs' performance drops when the game configuration is misaligned with the affecting biases. Performance is assessed based on the selection of the correct action, one which agrees with the prompted preferred behaviours of both players. Alignment refers to whether the LLM's bias aligns with the correct action. For example, GPT-4o's average performance drops by 34% when misaligned. Additionally, the current trend of "bigger and newer is better" does not hold for the above, where GPT-4o (the current best-performing LLM) suffers the most substantial performance drop. Lastly, we note that while chain-of-thought prompting does reduce the effect of the biases on most models, it is far from solving the problem at the fundamental level. △ Less

Submitted 16 July, 2024; v1 submitted 5 July, 2024; originally announced July 2024.

Comments: 8 pages (19 with appendix), 6 figures in the main body (4 in the appendix), 4 tables in the main body

A detailed study of the very-high-energy Crab pulsar emission with the LST-1

Authors: CTA-LST Project, :, K. Abe, S. Abe, A. Abhishek, F. Acero, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González , et al. (272 additional authors not shown)

Abstract: Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov… ▽ More Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov Telescope Array Observatory (CTAO). Its improved performance over previous IACTs makes it well suited for studying pulsars. Aims: To study the Crab pulsar emission with the LST-1, improving and complementing the results from other telescopes. These observations can also be used to characterize the potential of the LST-1 to study other pulsars and detect new ones. Methods: We analyzed a total of $\sim$103 hours of gamma-ray observations of the Crab pulsar conducted with the LST-1 in the period from September 2020 to January 2023. The observations were carried out at zenith angles less than 50 degrees. A new analysis of the Fermi-LAT data was also performed, including $\sim$14 years of observations. Results: The Crab pulsar phaseogram, long-term light-curve, and phase-resolved spectra are reconstructed with the LST-1 from 20 GeV to 450 GeV for P1 and up to 700 GeV for P2. The pulsed emission is detected with a significance of 15.2$σ$. The two characteristic emission peaks of the Crab pulsar are clearly detected (>10$σ$), as well as the so-called bridge emission (5.7$σ$). We find that both peaks are well described by power laws, with spectral indices of $\sim$3.44 and $\sim$3.03 respectively. The joint analysis of Fermi-LAT and LST-1 data shows a good agreement between both instruments in the overlapping energy range. The detailed results obtained in the first observations of the Crab pulsar with LST-1 show the potential that CTAO will have to study this type of sources. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted by A&A

Towards Generalist Robot Learning from Internet Video: A Survey

Authors: Robert McCarthy, Daniel C. H. Tan, Dominik Schmidt, Fernando Acero, Nathan Herr, Yilun Du, Thomas G. Thuruthel, Zhibin Li

Abstract: This survey presents an overview of methods for learning from video (LfV) in the context of reinforcement learning (RL) and robotics. We focus on methods capable of scaling to large internet video datasets and, in the process, extracting foundational knowledge about the world's dynamics and physical human behaviour. Such methods hold great promise for developing general-purpose robots. We open w… ▽ More This survey presents an overview of methods for learning from video (LfV) in the context of reinforcement learning (RL) and robotics. We focus on methods capable of scaling to large internet video datasets and, in the process, extracting foundational knowledge about the world's dynamics and physical human behaviour. Such methods hold great promise for developing general-purpose robots. We open with an overview of fundamental concepts relevant to the LfV-for-robotics setting. This includes a discussion of the exciting benefits LfV methods can offer (e.g., improved generalization beyond the available robot data) and commentary on key LfV challenges (e.g., missing information in video and LfV distribution shifts). Our literature review begins with an analysis of video foundation model techniques that can extract knowledge from large, heterogeneous video datasets. Next, we review methods that specifically leverage video data for robot learning. Here, we categorise work according to which RL knowledge modality (KM) benefits from the use of video data. We additionally highlight techniques for mitigating LfV challenges, including reviewing action representations that address missing action labels in video. Finally, we examine LfV datasets and benchmarks, before concluding with a discussion of challenges and opportunities in LfV. Here, we advocate for scalable foundation model approaches that can leverage the full range of internet video data, and that target the learning of the most promising RL KMs: the policy and dynamics model. Overall, we hope this survey will serve as a comprehensive reference for the emerging field of LfV, catalysing further research in the area and facilitating progress towards the development of general-purpose robots. △ Less

Submitted 7 June, 2024; v1 submitted 30 April, 2024; originally announced April 2024.

Comments: Updated formatting. Reduced paper length and made other minor improvements

Bayesian insights in Tycho supernova remnant : a detailed mapping of ejecta properties

Authors: L. Godinaud, F. Acero, A. Decourchelle, J. Ballet

Abstract: While Tycho's supernova remnant is one of the most studied type Ia Galactic supernova remnants, a global view of the physical properties of its ejecta is lacking, to understand its mysteries. In particular, the spatial distribution of the Si-rich ejecta line-of-sight velocity presents a large-scale unexplained asymmetry, with the north dominantly blueshifted and the south redshifted. To investigat… ▽ More While Tycho's supernova remnant is one of the most studied type Ia Galactic supernova remnants, a global view of the physical properties of its ejecta is lacking, to understand its mysteries. In particular, the spatial distribution of the Si-rich ejecta line-of-sight velocity presents a large-scale unexplained asymmetry, with the north dominantly blueshifted and the south redshifted. To investigate the origin of this line-of-sight velocity asymmetry in the ejecta, we carry out a detailed X-ray spatially-resolved spectral analysis of the entire shocked ejecta in Tycho's SNR to determine the physical properties of its various components. This study is based on the archival deep X-ray observations from the Chandra space telescope. The spatially-resolved spectral analysis in 211 regions over the entire SNR is based on a tesselation method applied to the line-of-sight velocity map. A Bayesian tool is used to conduct the fitting, using a nested sampling algorithm. It allows us to obtain a complete view of the statistical landscape. We provide maps of the physical parameters of the various components across the SNR ejecta. The Doppler shift map confirms spectrally the large-scale north-south asymmetry in the line-of-sight velocity. We reveal different spatial distributions of temperature and ionization time for IMEs and for iron-rich ejecta, but none of these maps shows structure associated to the large-scale north-south asymmetry in the line-of-sight velocity distribution. The abundance maps show spatial variations, depending on the element, perhaps due to an origin in different layers during the explosion. We compare these abundances with some nucleosynthesis models. In addition, we observe for the first time an emission line at 0.654 keV possibly related to oxygen. Its spatial distribution differs from the other elements, so that this line may arise in the ambient medium. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: Submitted to A&A, 20 pages, 15 figures

Early-time gamma-ray constraints on cosmic-ray acceleration in the core-collapse SN 2023ixf with the Fermi Large Area Telescope

Authors: G. Martí-Devesa, C. C. Cheung, N. Di Lalla, M. Renaud, G. Principe, N. Omodei, F. Acero

Abstract: While SNRs have been considered the most relevant Galactic CR accelerators for decades, CCSNe could accelerate particles during the earliest stages of their evolution and hence contribute to the CR energy budget in the Galaxy. Some SNRs have indeed been associated with TeV gamma-rays, yet proton acceleration efficiency during the early stages of an SN expansion remains mostly unconstrained. The mu… ▽ More While SNRs have been considered the most relevant Galactic CR accelerators for decades, CCSNe could accelerate particles during the earliest stages of their evolution and hence contribute to the CR energy budget in the Galaxy. Some SNRs have indeed been associated with TeV gamma-rays, yet proton acceleration efficiency during the early stages of an SN expansion remains mostly unconstrained. The multi-wavelength observation of SN 2023ixf, a Type II SN in the nearby galaxy M101, opens the possibility to constrain CR acceleration within a few days after the collapse of the RSG stellar progenitor. With this work, we intend to provide a phenomenological, quasi-model-independent constraint on the CR acceleration efficiency during this event at photon energies above 100 MeV. We performed a maximum-likelihood analysis of gamma-ray data from the Fermi Large Area Telescope up to one month after the SN explosion. We searched for high-energy emission from its expanding shock, and estimated the underlying hadronic CR energy reservoir assuming a power-law proton distribution consistent with standard diffusive shock acceleration. We do not find significant gamma-ray emission from SN 2023ixf. Nonetheless, our non-detection provides the first limit on the energy transferred to the population of hadronic CRs during the very early expansion of a CCSN. Under reasonable assumptions, our limits would imply a maximum efficiency on the CR acceleration of as low as 1%, which is inconsistent with the common estimate of 10% in generic SNe. However, this result is highly dependent on the assumed geometry of the circumstellar medium, and could be relaxed back to 10% by challenging spherical symmetry. A more sophisticated, inhomogeneous characterisation of the shock and the progenitor's environment is required before establishing whether or not Type II SNe are indeed efficient CR accelerators at early times. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Accepted in A&A. 13 pages, 12 figures, 4 tables

SUSHI: An algorithm for source separation of hyperspectral images with non-stationary spectral variation

Authors: Julia Lascar, Jérôme Bobin, Fabio Acero

Abstract: Hyperspectral images are data cubes with two spatial dimensions and a third spectral dimension, providing a spectrum for each pixel, and thus allow the mapping of extended sources' physical properties. In this article, we present the Semi-blind Unmixing with Sparsity for Hyperspectral Images (SUSHI), an algorithm for non-stationary unmixing of hyperspectral images with spatial regularization of sp… ▽ More Hyperspectral images are data cubes with two spatial dimensions and a third spectral dimension, providing a spectrum for each pixel, and thus allow the mapping of extended sources' physical properties. In this article, we present the Semi-blind Unmixing with Sparsity for Hyperspectral Images (SUSHI), an algorithm for non-stationary unmixing of hyperspectral images with spatial regularization of spectral parameters. The method allows for the disentangling of physical components without the assumption of a unique spectrum for each component. Thus, unlike most source separation methods used in astrophysics, all physical components obtained by SUSHI vary in spectral shape and in amplitude across the data cube. Non-stationary source separation is an ill-posed inverse problem that needs to be constrained. We achieve this by training a spectral model and applying a spatial regularization constraint on its parameters. For the spectral model, we used an Interpolatory Auto-Encoder, a generative model that can be trained with limited samples. For spatial regularization, we applied a sparsity constraint on the wavelet transform of the model parameter maps. We applied SUSHI to a toy model meant to resemble supernova remnants in X-ray astrophysics, though the method may be used on any extended source with any hyperspectral instrument. We compared this result to the one obtained by a classic 1D fit on each individual pixel. We find that SUSHI obtains more accurate results, particularly when it comes to reconstructing physical parameters. We applied SUSHI to real X-ray data from the supernova remnant Cassiopeia A and to the Crab Nebula. The results obtained are realistic and in accordance with past findings but have a much better spatial resolution. Thanks to spatial regularization, SUSHI can obtain reliable physical parameters at fine scales that are out of reach for pixel-by-pixel methods. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 17 pages, 22 figures, code available at https://github.com/JMLascar/SUSHI

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

Deep Reinforcement Learning and Mean-Variance Strategies for Responsible Portfolio Optimization

Authors: Fernando Acero, Parisa Zehtabi, Nicolas Marchesotti, Michael Cashmore, Daniele Magazzeni, Manuela Veloso

Abstract: Portfolio optimization involves determining the optimal allocation of portfolio assets in order to maximize a given investment objective. Traditionally, some form of mean-variance optimization is used with the aim of maximizing returns while minimizing risk, however, more recently, deep reinforcement learning formulations have been explored. Increasingly, investors have demonstrated an interest in… ▽ More Portfolio optimization involves determining the optimal allocation of portfolio assets in order to maximize a given investment objective. Traditionally, some form of mean-variance optimization is used with the aim of maximizing returns while minimizing risk, however, more recently, deep reinforcement learning formulations have been explored. Increasingly, investors have demonstrated an interest in incorporating ESG objectives when making investment decisions, and modifications to the classical mean-variance optimization framework have been developed. In this work, we study the use of deep reinforcement learning for responsible portfolio optimization, by incorporating ESG states and objectives, and provide comparisons against modified mean-variance approaches. Our results show that deep reinforcement learning policies can provide competitive performance against mean-variance approaches for responsible portfolio allocation across additive and multiplicative utility functions of financial and ESG responsibility objectives. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: Presented at the AAAI 2024 Workshop on AI in Finance for Social Impact

Distilling Reinforcement Learning Policies for Interpretable Robot Locomotion: Gradient Boosting Machines and Symbolic Regression

Authors: Fernando Acero, Zhibin Li

Abstract: Recent advancements in reinforcement learning (RL) have led to remarkable achievements in robot locomotion capabilities. However, the complexity and ``black-box'' nature of neural network-based RL policies hinder their interpretability and broader acceptance, particularly in applications demanding high levels of safety and reliability. This paper introduces a novel approach to distill neural RL po… ▽ More Recent advancements in reinforcement learning (RL) have led to remarkable achievements in robot locomotion capabilities. However, the complexity and ``black-box'' nature of neural network-based RL policies hinder their interpretability and broader acceptance, particularly in applications demanding high levels of safety and reliability. This paper introduces a novel approach to distill neural RL policies into more interpretable forms using Gradient Boosting Machines (GBMs), Explainable Boosting Machines (EBMs) and Symbolic Regression. By leveraging the inherent interpretability of generalized additive models, decision trees, and analytical expressions, we transform opaque neural network policies into more transparent ``glass-box'' models. We train expert neural network policies using RL and subsequently distill them into (i) GBMs, (ii) EBMs, and (iii) symbolic policies. To address the inherent distribution shift challenge of behavioral cloning, we propose to use the Dataset Aggregation (DAgger) algorithm with a curriculum of episode-dependent alternation of actions between expert and distilled policies, to enable efficient distillation of feedback control policies. We evaluate our approach on various robot locomotion gaits -- walking, trotting, bounding, and pacing -- and study the importance of different observations in joint actions for distilled policies using various methods. We train neural expert policies for 205 hours of simulated experience and distill interpretable policies with only 10 minutes of simulated interaction for each gait using the proposed method. △ Less

Submitted 21 March, 2024; originally announced March 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 23 July, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: 44 pages JCAP style (excluding author list and references), 19 figures; minor changes to match published version

Journal ref: JCAP 07 (2024) 047

Return level estimations for extreme rainfall over the Iberian Peninsula: comparing methodologies

Authors: F. J. Acero, S. Parey, J. A. García, D. Dacunha-Castelle

Abstract: Different ways to estimate future return levels for extreme rainfall are described and applied to the Iberian Peninsula (IP), based on Extreme Value Theory (EVT). This study is made for an ensemble of high quality rainfall time series observed in the Iberian Peninsula over the period 1961-2010. Both, peaks-over-threshold (POT) approach and block maxima with the Generalized Extreme Value (GEV) dist… ▽ More Different ways to estimate future return levels for extreme rainfall are described and applied to the Iberian Peninsula (IP), based on Extreme Value Theory (EVT). This study is made for an ensemble of high quality rainfall time series observed in the Iberian Peninsula over the period 1961-2010. Both, peaks-over-threshold (POT) approach and block maxima with the Generalized Extreme Value (GEV) distribution will be used and their results compared when linear trends are assumed in the parameters: threshold and scale parameter for POT and location and scale parameter for GEV. Both all-days and rainy-days-only data sets were considered, because rainfall over the IP is a special variable in that a large number of the values are 0. Another methodology is then tested, for rainy days only, considering the role of how the mean, variance, and number of rainy days evolve. The 20-year return levels (RLs) expected in 2020 were estimated using these methodologies for three seasons: fall, spring and winter. GEV is less reliable than POT because fixed blocks lead to the selection of non-extreme values. Future RLs obtained with POT are higher than those estimated with GEV, mainly for some gauges showing significant positive trend for the number of rainy days. Fall becomes the season with heaviest rainfall, rather than winter nowadays, for some regions. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 16 pages, 11 figures

Changes in heat waves characteristics over Extremadura (SW Spain)

Authors: F. J. Acero, M. I. Fernández-Fernández, V. M. S. Carrasco, S. Parey, T. T. Huong Hoang, D. Dacunha-Castelle, J. A. García

Abstract: Heat wave (HW) events are becoming more frequent, and they have important consequences because of the negative effects they can have not only on the human population in health terms, but also on biodiversity and agriculture. This motivated a study of the trends in HW events over Extremadura, a region in the southwest of Spain, with much of its area in summer devoted to the production of irrigated… ▽ More Heat wave (HW) events are becoming more frequent, and they have important consequences because of the negative effects they can have not only on the human population in health terms, but also on biodiversity and agriculture. This motivated a study of the trends in HW events over Extremadura, a region in the southwest of Spain, with much of its area in summer devoted to the production of irrigated crops such as maize and tomatoes. Heat waves were defined for the study as two consecutive days with temperatures above the 95th percentile of the summer (June-August) maximum temperature (Tmax) time series. Two datasets were used: one consisted of 13 daily temperature records uniformly distributed over the Region, and the other was the SPAIN02 gridded observational dataset, extracting just the points corresponding to Extremadura. The trends studied were in the duration, intensity, and frequency of HW events, and in other parameters such as the mean, low (25th percentile), and high (75th percentile) values. In general terms, the results showed significant positive trends in those parameters over the east, the northwest, and a small area in the south of the Region. In order to study changes in HW characteristics (duration, frequency and intensity) considering different subperiods, a stochastic model was used to generate 1000 time series equivalent to the observed ones. The results showed that there were no significant changes in HW duration in the last 10-year subperiod in comparison with the first. But the results were different for warm events (WE), defined with a lower threshold (the 75th percentile), which are also important for agriculture. For several sites, there were significant changes in WE duration, frequency, and intensity. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 36 pages, 11 figures

A limit for the values of the Dst geomagnetic index

Authors: F. J. Acero, J. M. Vaquero, M. C. Gallego, J. A. García

Abstract: The study of the extreme weather space events is important for a technological dependent society. Extreme Value Theory could be decisive to characterize those extreme events in order to have the knowledge to make decisions in technological, economic and social matters, in all fields with possible impacts. In this work, the hourly values of the Dst geomagnetic index has been studied for the period… ▽ More The study of the extreme weather space events is important for a technological dependent society. Extreme Value Theory could be decisive to characterize those extreme events in order to have the knowledge to make decisions in technological, economic and social matters, in all fields with possible impacts. In this work, the hourly values of the Dst geomagnetic index has been studied for the period 1957-2014 using the peaks-over-threshold technique. The shape parameter obtained from the fit of the generalized Pareto distribution to the extreme values of the |Dst| index leads to a negative value implying an upper bound for this time series. This result is relevant, because the estimation of this limit for the extreme values lead to 850 nT as the highest expected value for this geomagnetic index. Thus, from the previous characterization of the Carrington geomagnetic storm and our results, it could be considered the worst case scenario. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 12 pages, 5 figures

Mapping thermal emission in the synchrotron-dominated SNRs 330.2+1.0, 3C58, and RX J1713.7-3946

Authors: Adrien Picquenot, Brian J. Williams, Fabio Acero, Koji Mori

Abstract: Since the discovery of synchrotron X-ray emission from the shell of the supernova remnant (SNR) SN 1006, multiple observations from Chandra and XMM-Newton have shown that many young SNRs produce synchrotron emission in X-rays. Among those, a few peculiar SNRs have their X-ray emission largely dominated by synchrotron radiation, showing no or only faint traces of thermal emission. In this paper, we… ▽ More Since the discovery of synchrotron X-ray emission from the shell of the supernova remnant (SNR) SN 1006, multiple observations from Chandra and XMM-Newton have shown that many young SNRs produce synchrotron emission in X-rays. Among those, a few peculiar SNRs have their X-ray emission largely dominated by synchrotron radiation, showing no or only faint traces of thermal emission. In this paper, we report our mapping of the thermal emission in three emblematic synchrotron-dominated SNRs: G330.2+1.0, 3C58, and RX J1713.7-3946. We used a blind source separation method able to retrieve faint components from X-ray data in the form of Chandra and XMM-Newton observations. The thermal candidates disentangled by the algorithm were then used to select regions of extraction. We then analyzed the extracted spectra to assess their physical nature. We conclude that the components retrieved by the algorithm indeed represent the spatial distribution of the thermal emission in G330.2+1.0 and 3C58, and a likely thermal candidate in RX J1713.7-3946. Our findings confirm and expand on past studies. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Modular Neural Network Policies for Learning In-Flight Object Catching with a Robot Hand-Arm System

Authors: Wenbin Hu, Fernando Acero, Eleftherios Triantafyllidis, Zhaocheng Liu, Zhibin Li

Abstract: We present a modular framework designed to enable a robot hand-arm system to learn how to catch flying objects, a task that requires fast, reactive, and accurately-timed robot motions. Our framework consists of five core modules: (i) an object state estimator that learns object trajectory prediction, (ii) a catching pose quality network that learns to score and rank object poses for catching, (iii… ▽ More We present a modular framework designed to enable a robot hand-arm system to learn how to catch flying objects, a task that requires fast, reactive, and accurately-timed robot motions. Our framework consists of five core modules: (i) an object state estimator that learns object trajectory prediction, (ii) a catching pose quality network that learns to score and rank object poses for catching, (iii) a reaching control policy trained to move the robot hand to pre-catch poses, (iv) a grasping control policy trained to perform soft catching motions for safe and robust grasping, and (v) a gating network trained to synthesize the actions given by the reaching and grasping policy. The former two modules are trained via supervised learning and the latter three use deep reinforcement learning in a simulated environment. We conduct extensive evaluations of our framework in simulation for each module and the integrated system, to demonstrate high success rates of in-flight catching and robustness to perturbations and sensory noise. Whilst only simple cylindrical and spherical objects are used for training, the integrated system shows successful generalization to a variety of household objects that are not used in training. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 8 pages. Accepted and presented at IEEE IROS 2023

Prospects for a survey of the Galactic plane with the Cherenkov Telescope Array

Authors: CTA Consortium

Abstract: Approximately one hundred sources of very-high-energy (VHE) gamma rays are known in the Milky Way. A survey of the entire Galactic Plane in the energy range from a few tens of GeV to a few hundred TeV has been proposed as a Key Science Project for the upcoming Cherenkov Telescope Array Observatory (CTAO). This article presents the status of the studies towards the Galactic Plane Survey (GPS). We b… ▽ More Approximately one hundred sources of very-high-energy (VHE) gamma rays are known in the Milky Way. A survey of the entire Galactic Plane in the energy range from a few tens of GeV to a few hundred TeV has been proposed as a Key Science Project for the upcoming Cherenkov Telescope Array Observatory (CTAO). This article presents the status of the studies towards the Galactic Plane Survey (GPS). We build and make publicly available a sky model that combines data from observations of known gamma-ray emitters with state-of-the-art physically-driven models of synthetic populations of the main classes of established Galactic VHE sources, as well as of interstellar emission from cosmic-ray interactions in the Milky Way. We also perform an optimisation of the observation strategy. We use the improved sky model and observation strategy to simulate GPS data that are analysed using the methods and software tools under development for real data. We show that the GPS has the potential to increase the number of known Galactic VHE emitters by almost a factor of five. This corresponds to the detection of more than two hundred pulsar wind nebulae and a few tens of supernova remnants at average integral fluxes one order of magnitude lower than in the existing sample above 1 TeV, therefore opening the possibility to perform unprecedented population studies. The GPS also has the potential to provide new VHE detections of binary systems and pulsars, and to detect bright PeVatrons. Furthermore, the GPS will constitute a pathfinder for deeper follow-up observations of these source classes. Finally, we show that we can extract from GPS data an estimate of the contribution to diffuse emission from unresolved sources, and that there are good prospects of detecting interstellar emission and statistically distinguishing different scenarios. (Abridged) △ Less

Submitted 16 July, 2024; v1 submitted 4 October, 2023; originally announced October 2023.

Comments: Accepted for publication in JCAP. The sky model used in the paper is available at https://doi.org/10.5281/zenodo.8402519. Corresponding authors: F. Acero, M. Chernyakova, B. Olmi, Q. Remy, L. Tibaldo

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

A fresh perspective on the 3D dynamics of Tycho's supernova remnant: Ejecta asymmetries in the X-ray band

Authors: Leila Godinaud, Fabio Acero, Anne Decourchelle, Jean Ballet

Abstract: 450 years after the explosion of the Type Ia SN1572, the dynamics of the Tycho supernova remnant can give us keys to understand the explosion mechanism and the interaction of the remnant with the interstellar medium. To probe the asymmetries and the evolution of the SNR, we track the ejecta dynamics using new methods applied to the deep X-ray observations available in the Chandra space telescope a… ▽ More 450 years after the explosion of the Type Ia SN1572, the dynamics of the Tycho supernova remnant can give us keys to understand the explosion mechanism and the interaction of the remnant with the interstellar medium. To probe the asymmetries and the evolution of the SNR, we track the ejecta dynamics using new methods applied to the deep X-ray observations available in the Chandra space telescope archive. For the line of sight velocity measurement Vz, we use the Doppler effect focused on the bright Si line in the 1.6-2.1 keV band. Using the component separation tool General Morphological Component Analysis (GMCA), we successfully disentangle the red and blueshifted Si ejecta emission. This allows us to reconstruct a map of the peak energy of the Si line with a total coverage of the SNR at a 2'' resolution and a proxy of the velocity in the line of sight. For the proper motions in the plane of the sky Vxy, we develop a new method, named Poisson Optical Flow, to measure the displacement of 2D features between the observations of 2003 and 2009. The result is a field of 1700 velocity vectors covering the entire SNR. These exhaustive 3D velocity measurements reveal the complex and patchy dynamics of the SNR. At the large-scale, an asymmetry with the North being dominantly blueshifted and the South redshifted is observed. The proper motion vector field Vxy highlights different dynamics between the East and the West parts of the SNR. The eastern velocity field is more disturbed by external inhomogeneities and the South-East ejecta knot. In particular, a slow-down is observed in the North-East which could be due to the interaction with higher densities as seen in other wavelengths. The vector field is also used to backtrace the center of the explosion which is then compared with potential stellar progenitors distances from the latest Gaia DR3, leaving only stars B and E as possible candidates. △ Less

Submitted 29 February, 2024; v1 submitted 4 September, 2023; originally announced September 2023.

Comments: Published in Astronomy & Astrophysics, Volume 680, id.A80, 18 pp. 18 pages, 12 figures

Journal ref: A&A 680, A80 (2023)

Gammapy: A Python package for gamma-ray astronomy

Authors: Axel Donath, Régis Terrier, Quentin Remy, Atreyee Sinha, Cosimo Nigro, Fabio Pintore, Bruno Khélifi, Laura Olivera-Nieto, Jose Enrique Ruiz, Kai Brügge, Maximilian Linhoff, Jose Luis Contreras, Fabio Acero, Arnau Aguasca-Cabot, David Berge, Pooja Bhattacharjee, Johannes Buchner, Catherine Boisson, David Carreto Fidalgo, Andrew Chen, Mathieu de Bony de Lavergne, José Vinícius de Miranda Cardoso, Christoph Deil, Matthias Füßling, Stefan Funk , et al. (20 additional authors not shown)

Abstract: In this article, we present Gammapy, an open-source Python package for the analysis of astronomical $γ$-ray data, and illustrate the functionalities of its first long-term-support release, version 1.0. Built on the modern Python scientific ecosystem, Gammapy provides a uniform platform for reducing and modeling data from different $γ$-ray instruments for many analysis scenarios. Gammapy complies w… ▽ More In this article, we present Gammapy, an open-source Python package for the analysis of astronomical $γ$-ray data, and illustrate the functionalities of its first long-term-support release, version 1.0. Built on the modern Python scientific ecosystem, Gammapy provides a uniform platform for reducing and modeling data from different $γ$-ray instruments for many analysis scenarios. Gammapy complies with several well-established data conventions in high-energy astrophysics, providing serialized data products that are interoperable with other software packages. Starting from event lists and instrument response functions, Gammapy provides functionalities to reduce these data by binning them in energy and sky coordinates. Several techniques for background estimation are implemented in the package to handle the residual hadronic background affecting $γ$-ray instruments. After the data are binned, the flux and morphology of one or more $γ$-ray sources can be estimated using Poisson maximum likelihood fitting and assuming a variety of spectral, temporal, and spatial models. Estimation of flux points, likelihood profiles, and light curves is also supported. After describing the structure of the package, we show, using publicly available $γ$-ray data, the capabilities of Gammapy in multiple traditional and novel $γ$-ray analysis scenarios, such as spectral and spectro-morphological modeling and estimations of a spectral energy distribution and a light curve. Its flexibility and power are displayed in a final multi-instrument example, where datasets from different instruments, at different stages of data reduction, are simultaneously fitted with an astrophysical flux model. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: 26 pages, 16 figures

Journal ref: A&A 678, A157 (2023)

Accelerating Cutting-Plane Algorithms via Reinforcement Learning Surrogates

Authors: Kyle Mana, Fernando Acero, Stephen Mak, Parisa Zehtabi, Michael Cashmore, Daniele Magazzeni, Manuela Veloso

Abstract: Discrete optimization belongs to the set of $\mathcal{NP}$-hard problems, spanning fields such as mixed-integer programming and combinatorial optimization. A current standard approach to solving convex discrete optimization problems is the use of cutting-plane algorithms, which reach optimal solutions by iteratively adding inequalities known as \textit{cuts} to refine a feasible set. Despite the e… ▽ More Discrete optimization belongs to the set of $\mathcal{NP}$-hard problems, spanning fields such as mixed-integer programming and combinatorial optimization. A current standard approach to solving convex discrete optimization problems is the use of cutting-plane algorithms, which reach optimal solutions by iteratively adding inequalities known as \textit{cuts} to refine a feasible set. Despite the existence of a number of general-purpose cut-generating algorithms, large-scale discrete optimization problems continue to suffer from intractability. In this work, we propose a method for accelerating cutting-plane algorithms via reinforcement learning. Our approach uses learned policies as surrogates for $\mathcal{NP}$-hard elements of the cut generating procedure in a way that (i) accelerates convergence, and (ii) retains guarantees of optimality. We apply our method on two types of problems where cutting-plane algorithms are commonly used: stochastic optimization, and mixed-integer quadratic programming. We observe the benefits of our method when applied to Benders decomposition (stochastic optimization) and iterative loss approximation (quadratic programming), achieving up to $45\%$ faster average convergence when compared to modern alternative algorithms. △ Less

Submitted 27 February, 2024; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: Extended version (includes Supplementary Material). Accepted at AAAI 24 Main Track with Oral Presentation

RObotic MAnipulation Network (ROMAN) $\unicode{x2013}$ Hybrid Hierarchical Learning for Solving Complex Sequential Tasks

Authors: Eleftherios Triantafyllidis, Fernando Acero, Zhaocheng Liu, Zhibin Li

Abstract: Solving long sequential tasks poses a significant challenge in embodied artificial intelligence. Enabling a robotic system to perform diverse sequential tasks with a broad range of manipulation skills is an active area of research. In this work, we present a Hybrid Hierarchical Learning framework, the Robotic Manipulation Network (ROMAN), to address the challenge of solving multiple complex tasks… ▽ More Solving long sequential tasks poses a significant challenge in embodied artificial intelligence. Enabling a robotic system to perform diverse sequential tasks with a broad range of manipulation skills is an active area of research. In this work, we present a Hybrid Hierarchical Learning framework, the Robotic Manipulation Network (ROMAN), to address the challenge of solving multiple complex tasks over long time horizons in robotic manipulation. ROMAN achieves task versatility and robust failure recovery by integrating behavioural cloning, imitation learning, and reinforcement learning. It consists of a central manipulation network that coordinates an ensemble of various neural networks, each specialising in distinct re-combinable sub-tasks to generate their correct in-sequence actions for solving complex long-horizon manipulation tasks. Experimental results show that by orchestrating and activating these specialised manipulation experts, ROMAN generates correct sequential activations for accomplishing long sequences of sophisticated manipulation tasks and achieving adaptive behaviours beyond demonstrations, while exhibiting robustness to various sensory noises. These results demonstrate the significance and versatility of ROMAN's dynamic adaptability featuring autonomous failure recovery capabilities, and highlight its potential for various autonomous manipulation tasks that demand adaptive motor skills. △ Less

Submitted 7 July, 2023; v1 submitted 30 June, 2023; originally announced July 2023.

Comments: To appear in Nature Machine Intelligence. Includes the main and supplementary manuscript. Total of 70 pages, with a total of 9 Figures and 17 Tables

Value Functions are Control Barrier Functions: Verification of Safe Policies using Control Theory

Authors: Daniel C. H. Tan, Fernando Acero, Robert McCarthy, Dimitrios Kanoulas, Zhibin Li

Abstract: Guaranteeing safe behaviour of reinforcement learning (RL) policies poses significant challenges for safety-critical applications, despite RL's generality and scalability. To address this, we propose a new approach to apply verification methods from control theory to learned value functions. By analyzing task structures for safety preservation, we formalize original theorems that establish links b… ▽ More Guaranteeing safe behaviour of reinforcement learning (RL) policies poses significant challenges for safety-critical applications, despite RL's generality and scalability. To address this, we propose a new approach to apply verification methods from control theory to learned value functions. By analyzing task structures for safety preservation, we formalize original theorems that establish links between value functions and control barrier functions. Further, we propose novel metrics for verifying value functions in safe control tasks and practical implementation details to improve learning. Our work presents a novel method for certificate learning, which unlocks a diversity of verification techniques from control theory for RL policies, and marks a significant step towards a formal framework for the general, scalable, and verifiable design of RL-based control systems. Code and videos are available at this https url: https://rl-cbf.github.io/ △ Less

Submitted 5 December, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

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

The X-ray synchrotron rims in Cassiopeia A narrow with energy

Authors: Adrien Picquenot, Brian J. Williams, Fabio Acero, Benson T. Guest

Abstract: Some young supernova remnants exhibit thin filaments of X-ray synchrotron radiation coinciding with the forward shock due to accelerated electrons interacting with the local magnetic field. The two main models accounting for the radial brightness evolution of these filaments differ in their prediction of the narrowing (or not) of the filaments with increasing photon energy. In this paper, we repor… ▽ More Some young supernova remnants exhibit thin filaments of X-ray synchrotron radiation coinciding with the forward shock due to accelerated electrons interacting with the local magnetic field. The two main models accounting for the radial brightness evolution of these filaments differ in their prediction of the narrowing (or not) of the filaments with increasing photon energy. In this paper, we report our observation of such a narrowing of the synchrotron filaments in Cassiopeia A at X-ray energies, and how this finding could help in understanding the mechanisms at stake in their formation. We used a new blind source separation method on the 1 Ms Chandra observation of Cassiopeia A, in order to obtain detailed and unpolluted images of the synchrotron emission in three energy bands. We then extracted the profiles of several filaments at the forward shock and the reverse shock to estimate and compare their widths. We find that there is indeed a narrowing with energy of the synchrotron filaments both at the forward and at the reverse shocks in Cassiopeia A. The energy dependency of this narrowing seems stronger at high energy, which is indicative of a damping effect, confirmed by radio observations. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: Accepted by A&A

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

Fermi-LAT Gamma-ray Emission Discovered from the Composite Supernova Remnant B0453-685 in the Large Magellanic Cloud

Authors: Jordan Eagle, Daniel Castro, Peter Mahhov, Joseph Gelfand, Matthew Kerr, Patrick Slane, Jean Ballet, Fabio Acero, Samayra Straal, Marco Ajello

Abstract: We report the second extragalactic pulsar wind nebula (PWN) to be detected in the MeV-GeV band by the Fermi-LAT, located within the Large Magellanic Cloud (LMC). The only other known PWN to emit in the Fermi band outside of the Milky Way Galaxy is N 157B which lies to the west of the newly detected gamma-ray emission at an angular distance of 4 degrees. Faint, point-like gamma-ray emission is disc… ▽ More We report the second extragalactic pulsar wind nebula (PWN) to be detected in the MeV-GeV band by the Fermi-LAT, located within the Large Magellanic Cloud (LMC). The only other known PWN to emit in the Fermi band outside of the Milky Way Galaxy is N 157B which lies to the west of the newly detected gamma-ray emission at an angular distance of 4 degrees. Faint, point-like gamma-ray emission is discovered at the location of the composite supernova remnant (SNR) B0453-685 with a ~ 4 sigma significance from energies 300 MeV - 2 TeV. We present the Fermi-LAT data analysis of the new gamma-ray source, coupled with a detailed multi-wavelength investigation to understand the nature of the observed emission. Combining the observed characteristics of the SNR and the physical implications from broadband modeling, we argue it is unlikely the SNR is responsible for the gamma-ray emission. While the gamma-ray emission is too faint for a pulsation search, we try to distinguish between any pulsar and PWN component of SNR B0453-685 that would be responsible for the observed gamma-ray emission using semi-analytic models. We determine the most likely scenario is that the old PWN (t ~ 14,000 years) within B0453-685 has been impacted by the return of the SNR reverse shock with a possible substantial pulsar component below 5 GeV. △ Less

Submitted 3 February, 2023; originally announced February 2023.

Comments: 18 pages, 8 figures, accepted for publication in ApJ Feb. 1, 2023

Semi-Blind Source Separation with Learned Constraints

Authors: Rémi Carloni Gertosio, Jérôme Bobin, Fabio Acero

Abstract: Blind source separation (BSS) algorithms are unsupervised methods, which are the cornerstone of hyperspectral data analysis by allowing for physically meaningful data decompositions. BSS problems being ill-posed, the resolution requires efficient regularization schemes to better distinguish between the sources and yield interpretable solutions. For that purpose, we investigate a semi-supervised so… ▽ More Blind source separation (BSS) algorithms are unsupervised methods, which are the cornerstone of hyperspectral data analysis by allowing for physically meaningful data decompositions. BSS problems being ill-posed, the resolution requires efficient regularization schemes to better distinguish between the sources and yield interpretable solutions. For that purpose, we investigate a semi-supervised source separation approach in which we combine a projected alternating least-square algorithm with a learning-based regularization scheme. In this article, we focus on constraining the mixing matrix to belong to a learned manifold by making use of generative models. Altogether, we show that this allows for an innovative BSS algorithm, with improved accuracy, which provides physically interpretable solutions. The proposed method, coined sGMCA, is tested on realistic hyperspectral astrophysical data in challenging scenarios involving strong noise, highly correlated spectra and unbalanced sources. The results highlight the significant benefit of the learned prior to reduce the leakages between the sources, which allows an overall better disentanglement. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Journal ref: Signal Processing, Volume 202, January 2023, 108776

Constraining leptonic emission scenarios for the PeVatron candidate HESS J1702-420 with deep XMM-Newton observations

Authors: L. Giunti, F. Acero, B. Khelifi, K. Kosack, A. Lemiere, R. Terrier

Abstract: The unidentified TeV source HESS J1702-420 has recently been proposed as a new hadronic PeVatron candidate, based on the discovery of a small-scale emission sub-region with extremely hard gamma-ray spectrum up to 100 TeV (named HESS J1702-420A). Given the difficulty to discriminate between a hadronic or leptonic origin of the TeV emission, based on the H.E.S.S. measurement alone, we opted for a mu… ▽ More The unidentified TeV source HESS J1702-420 has recently been proposed as a new hadronic PeVatron candidate, based on the discovery of a small-scale emission sub-region with extremely hard gamma-ray spectrum up to 100 TeV (named HESS J1702-420A). Given the difficulty to discriminate between a hadronic or leptonic origin of the TeV emission, based on the H.E.S.S. measurement alone, we opted for a multi-wavelength approach. A deep X-ray observation was carried out using the XMM-Newton satellite, with the goal of probing a possible association with a hidden leptonic accelerator. No evidence of a clear counterpart for HESS J1702-420A was found in the X-ray data. After excluding an association with all nearby X-ray point sources, we derived strict upper limits on the diffuse X-ray emission and average magnetic field in the HESS J1702-420A region. We additionally report the serendipitous discovery of a new extended X-ray source, whose association with HESS J1702-420A is not obvious but cannot be ruled out either. A set of scripts dedicated to the multi-wavelength modeling of X-ray and gamma-ray data, based on Gammapy, Naima and Xspec, was developed in the context of this work and is made publicly available along with this paper. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: Accepted for publication in Astronomy & Astrophysics, 13 pages, 11 figures, 3 tables

The Athena X-ray Integral Field Unit: a consolidated design for the system requirement review of the preliminary definition phase

Authors: Didier Barret, Vincent Albouys, Jan-Willem den Herder, Luigi Piro, Massimo Cappi, Juhani Huovelin, Richard Kelley, J. Miguel Mas-Hesse, Stéphane Paltani, Gregor Rauw, Agata Rozanska, Jiri Svoboda, Joern Wilms, Noriko Yamasaki, Marc Audard, Simon Bandler, Marco Barbera, Xavier Barcons, Enrico Bozzo, Maria Teresa Ceballos, Ivan Charles, Elisa Costantini, Thomas Dauser, Anne Decourchelle, Lionel Duband , et al. (274 additional authors not shown)

Abstract: The Athena X-ray Integral Unit (X-IFU) is the high resolution X-ray spectrometer, studied since 2015 for flying in the mid-30s on the Athena space X-ray Observatory, a versatile observatory designed to address the Hot and Energetic Universe science theme, selected in November 2013 by the Survey Science Committee. Based on a large format array of Transition Edge Sensors (TES), it aims to provide sp… ▽ More The Athena X-ray Integral Unit (X-IFU) is the high resolution X-ray spectrometer, studied since 2015 for flying in the mid-30s on the Athena space X-ray Observatory, a versatile observatory designed to address the Hot and Energetic Universe science theme, selected in November 2013 by the Survey Science Committee. Based on a large format array of Transition Edge Sensors (TES), it aims to provide spatially resolved X-ray spectroscopy, with a spectral resolution of 2.5 eV (up to 7 keV) over an hexagonal field of view of 5 arc minutes (equivalent diameter). The X-IFU entered its System Requirement Review (SRR) in June 2022, at about the same time when ESA called for an overall X-IFU redesign (including the X-IFU cryostat and the cooling chain), due to an unanticipated cost overrun of Athena. In this paper, after illustrating the breakthrough capabilities of the X-IFU, we describe the instrument as presented at its SRR, browsing through all the subsystems and associated requirements. We then show the instrument budgets, with a particular emphasis on the anticipated budgets of some of its key performance parameters. Finally we briefly discuss on the ongoing key technology demonstration activities, the calibration and the activities foreseen in the X-IFU Instrument Science Center, and touch on communication and outreach activities, the consortium organisation, and finally on the life cycle assessment of X-IFU aiming at minimising the environmental footprint, associated with the development of the instrument. Thanks to the studies conducted so far on X-IFU, it is expected that along the design-to-cost exercise requested by ESA, the X-IFU will maintain flagship capabilities in spatially resolved high resolution X-ray spectroscopy, enabling most of the original X-IFU related scientific objectives of the Athena mission to be retained. (abridged). △ Less

Submitted 28 November, 2022; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: 48 pages, 29 figures, Accepted for publication in Experimental Astronomy with minor editing

Search for new cosmic-ray acceleration sites within the 4FGL catalog Galactic plane sources

Authors: Fermi-LAT Collaboration, S. Abdollahi, F. Acero, M. Ackermann, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, B. Berenji, A. Berretta, E. Bissaldi, R. D. Blandford, R. Bonino, P. Bruel, S. Buson, R. A. Cameron, R. Caputo, P. A. Caraveo, D. Castro, G. Chiaro, N. Cibrario, S. Ciprini, J. Coronado-Blázquez, M. Crnogorcevic , et al. (95 additional authors not shown)

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

Submitted 6 May, 2022; originally announced May 2022.

Comments: Accepted for publication in The Astrophysical Journal

Incremental Fermi Large Area Telescope Fourth Source Catalog

Authors: Fermi-LAT collaboration, :, Soheila Abdollahi, Fabio Acero, Luca Baldini, Jean Ballet, Denis Bastieri, Ronaldo Bellazzini, Bijan Berenji, Alessandra Berretta, Elisabetta Bissaldi, Roger D. Blandford, Elliott Bloom, Raffaella Bonino, Ari Brill, Richard J. Britto, Philippe Bruel, Toby H. Burnett, Sara Buson, Rob A. Cameron, Regina Caputo, Patrizia A. Caraveo, Daniel Castro, Sylvain Chaty, Teddy C. Cheung , et al. (116 additional authors not shown)

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

Submitted 10 May, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: accepted in ApJS; follow-up paper to 1902.10045

Journal ref: ApJS 260, 53 (2022)

Characterization of the GeV emission from the Kepler supernova remnant

Authors: Fabio Acero, Marianne Lemoine-Goumard, Jean Ballet

Abstract: The Kepler supernova remnant (SNR) is the only historic supernova remnant lacking a detection at GeV and TeV energies which probe particle acceleration. A recent analysis of Fermi-LAT data reported a likely GeV gamma-ray candidate in the direction of the SNR. Using approximately the same dataset but with an optimized analysis configuration, we confirm the gamma-ray candidate to a solid $>6σ$ detec… ▽ More The Kepler supernova remnant (SNR) is the only historic supernova remnant lacking a detection at GeV and TeV energies which probe particle acceleration. A recent analysis of Fermi-LAT data reported a likely GeV gamma-ray candidate in the direction of the SNR. Using approximately the same dataset but with an optimized analysis configuration, we confirm the gamma-ray candidate to a solid $>6σ$ detection and report a spectral index of $2.14 \pm 0.12_{\rm stat} \pm 0.15_{\rm syst}$ for an energy flux above 100 MeV of $(3.1 \pm 0.6_{\rm stat} \pm 0.3_{\rm syst}) \times 10^{-12}$ erg~cm$^{-2}$~s$^{-1}$. The gamma-ray excess is not significantly extended and is fully compatible with the radio, infrared or X-ray spatial distribution of the SNR. We successfully characterized this multi-wavelength emission with a model in which accelerated particles interact with the dense circumstellar material in the North-West portion of the SNR and radiate GeV gamma-rays through $π^{o}$ decay. The X-ray synchrotron and inverse-Compton (IC) emission mostly stem from the fast shocks in the southern regions with a magnetic field B$\sim$100 $μ$G or higher. Depending on the exact magnetic field amplitude, the TeV emission could arise from either the South region (IC dominated) or the interaction region ($π^{o}$ decay dominated). △ Less

Submitted 14 January, 2022; originally announced January 2022.

Comments: accepted for publication in A&A

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

Possible Detection of X-Ray Emitting Circumstellar Material in the Synchrotron-Dominated Supernova Remnant RX J 1713.7-3946

Authors: Dai Tateishi, Satoru Katsuda, Yukikatsu Terada, Fabio Acero, Takashi Yoshida, Shin-ichiro Fujimoto, Hidetoshi Sano

Abstract: We report on a discovery of an X-ray emitting circumstellar material knot inside the synchrotron dominant supernova remnant (SNR) RX J1713.7-3946. This knot was previously thought to be a Wolf-Rayet star (WR 85), but we realized that it is in fact $\sim$40$^{\prime\prime}$ away from WR 85, indicating no relation to WR 85. We performed high-resolution X-ray spectroscopy with the Reflection Grating… ▽ More We report on a discovery of an X-ray emitting circumstellar material knot inside the synchrotron dominant supernova remnant (SNR) RX J1713.7-3946. This knot was previously thought to be a Wolf-Rayet star (WR 85), but we realized that it is in fact $\sim$40$^{\prime\prime}$ away from WR 85, indicating no relation to WR 85. We performed high-resolution X-ray spectroscopy with the Reflection Grating Spectrometer (RGS) on board XMM-Newton. The RGS spectrum clearly resolves a number of emission lines, such as N Ly$α$, O Ly$α$, Fe XVIII, Ne X, Mg XI, and Si XIII. The spectrum can be well represented by an absorbed thermal emission model with a temperature of $k_{\rm B}T_{\rm e} = 0.65\pm 0.02$ keV. The elemental abundances are obtained to be ${\rm N/H} = 3.5\pm 0.8{\rm \left(N/H\right)_{\odot}}$, ${\rm O/H} = 0.5\pm0.1{\rm \left(O/H\right)_{\odot}}$, ${\rm Ne/H} = 0.9\pm0.1{\rm \left(Ne/H\right)_{\odot}}$, ${\rm Mg/H} = 1.0\pm0.1{\rm \left(Mg/H\right)_{\odot}}$, ${\rm Si/H} = 1.0\pm0.2{\rm \left(Si/H\right)_{\odot}}$, and ${\rm Fe/H} = 1.3\pm0.1{\rm \left(Fe/H\right)_{\odot}}$. The enhanced N abundance with others being about the solar values allows us to infer that this knot is circumstellar material ejected when the progenitor star evolved into a red supergiant. The abundance ratio of N to O is obtained to be $\rm N/O = 6.8_{-2.1}^{+2.5}\left(N/O\right)_{\odot}$. By comparing this to those in outer layers of red supergiant stars expected from stellar evolution simulations, we estimate the initial mass of the progenitor star to be $15\, \rm M_{\odot} \lesssim \rm M \lesssim 20\, \rm M_{\odot}$. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: 13 pages, 12 figures, 4 tables, accepted for publication in ApJ

Multiple accelerated particle populations in the Cygnus Loop with Fermi-LAT

Authors: A. Tutone, J. Ballet, F. Acero, A. D'Aì, G. Cusumano

Abstract: The Cygnus Loop (G74.0-8.5) is a very well-known nearby supernova remnant (SNR) in our Galaxy. Thanks to its large size, brightness, and angular offset from the Galactic plane, it has been studied in detail from radio to $γ$-ray emission. The $γ$ -rays probe the populations of energetic particles and their acceleration mechanisms at low shock speeds. We present an analysis of the $γ$-ray emission… ▽ More The Cygnus Loop (G74.0-8.5) is a very well-known nearby supernova remnant (SNR) in our Galaxy. Thanks to its large size, brightness, and angular offset from the Galactic plane, it has been studied in detail from radio to $γ$-ray emission. The $γ$ -rays probe the populations of energetic particles and their acceleration mechanisms at low shock speeds. We present an analysis of the $γ$-ray emission detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope over 11 years in the region of the Cygnus Loop. We performed detailed morphological and spectral studies of the $γ$-ray emission toward the remnant from 100 MeV to 100 GeV and compared it with X-ray, UV, optical, and radio images. The higher statistics with respect to the previous studies enabled us to decompose the emission from the remnant into two morphological components to model its nonthermal multiwavelength emission. The extended $γ$-ray emission is well correlated with the thermal X-ray and UV emission of the SNR. Our morphological analysis reveals that a model considering two contributions from the X-ray and the UV emission regions is the best description of the $γ$-ray data. Both components show a curved spectrum, but the X-ray component is softer and more curved than the UV component, suggesting a different physical origin. The multiwavelength modeling of emission toward the SNR suggests that the nonthermal radio and $γ$-ray emission associated with the UV component is mostly due to the reacceleration of preexisting cosmic rays by radiative shocks in the adjacent clouds, while the nonthermal emission associated with the X-ray component arises from freshly accelerated cosmic rays. △ Less

Submitted 9 December, 2021; v1 submitted 30 September, 2021; originally announced September 2021.

Comments: 15 pages, 12 figures, accepted to A&A, LaTeX typos corrected

Journal ref: A&A 656, A139 (2021)

Learning Perceptual Locomotion on Uneven Terrains using Sparse Visual Observations

Authors: Fernando Acero, Kai Yuan, Zhibin Li

Abstract: To proactively navigate and traverse various terrains, active use of visual perception becomes indispensable. We aim to investigate the feasibility and performance of using sparse visual observations to achieve perceptual locomotion over a range of common terrains (steps, ramps, gaps, and stairs) in human-centered environments. We formulate a selection of sparse visual inputs suitable for locomoti… ▽ More To proactively navigate and traverse various terrains, active use of visual perception becomes indispensable. We aim to investigate the feasibility and performance of using sparse visual observations to achieve perceptual locomotion over a range of common terrains (steps, ramps, gaps, and stairs) in human-centered environments. We formulate a selection of sparse visual inputs suitable for locomotion over the terrains of interest, and propose a learning framework to integrate exteroceptive and proprioceptive states. We specifically design the state observations and a training curriculum to learn feedback control policies effectively over a range of different terrains. We extensively validate and benchmark the learned policy in various tasks: omnidirectional walking on flat ground, and forward locomotion over various obstacles, showing high success rate of traversability. Furthermore, we study exteroceptive ablations and evaluate policy generalization by adding various levels of noise and testing on new unseen terrains. We demonstrate the capabilities of autonomous perceptual locomotion that can be achieved by only using sparse visual observations from direct depth measurements, which are easily available from a Lidar or RGB-D sensor, showing robust ascent and descent over high stairs of 20 cm height, i.e., 50% leg length, and robustness against noise and unseen terrains. △ Less

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

Comments: Video summary can be found at https://youtu.be/vtp43jYQ5w4

Learning Vision-Guided Dynamic Locomotion Over Challenging Terrains

Authors: Zhaocheng Liu, Fernando Acero, Zhibin Li

Abstract: Legged robots are becoming increasingly powerful and popular in recent years for their potential to bring the mobility of autonomous agents to the next level. This work presents a deep reinforcement learning approach that learns a robust Lidar-based perceptual locomotion policy in a partially observable environment using Proximal Policy Optimisation. Visual perception is critical to actively overc… ▽ More Legged robots are becoming increasingly powerful and popular in recent years for their potential to bring the mobility of autonomous agents to the next level. This work presents a deep reinforcement learning approach that learns a robust Lidar-based perceptual locomotion policy in a partially observable environment using Proximal Policy Optimisation. Visual perception is critical to actively overcome challenging terrains, and to do so, we propose a novel learning strategy: Dynamic Reward Strategy (DRS), which serves as effective heuristics to learn a versatile gait using a neural network architecture without the need to access the history data. Moreover, in a modified version of the OpenAI gym environment, the proposed work is evaluated with scores over 90% success rate in all tested challenging terrains. △ Less

Submitted 9 September, 2021; originally announced September 2021.

Comments: 9 pages, 27 figures, 1 table

Survey of the Galactic Plane with the CherenkovTelescope Array

Authors: Q. Remy, L. Tibaldo, F. Acero, M. Fiori, J. Knödlseder, B. Olmi, P. Sharma

Abstract: Observations with the current generation of very-high-energy gamma-ray telescopes have revealed an astonishing variety of particle accelerators in the Milky Way, such as supernova remnants, pulsar wind nebulae, and binary systems. The upcoming Cherenkov Telescope Array (CTA) will be the first instrument to enable a survey of the entire Galactic plane in the energy range from a few tens of GeV to 3… ▽ More Observations with the current generation of very-high-energy gamma-ray telescopes have revealed an astonishing variety of particle accelerators in the Milky Way, such as supernova remnants, pulsar wind nebulae, and binary systems. The upcoming Cherenkov Telescope Array (CTA) will be the first instrument to enable a survey of the entire Galactic plane in the energy range from a few tens of GeV to 300 TeV with unprecedented sensitivity and improved angular resolution. In this contribution we will revisit the scientific motivations for the survey, proposed as a Key ScienceProject for CTA. We will highlight recent progress, including improved physically-motivated models for Galactic source populations and interstellar emission, advance on the optimization of the survey strategy, and the development of pipelines to derive source catalogues tested on simulated data. Based on this, we will provide a new forecast on the properties of the sources thatCTA will detect and discuss the expected scientific return from the study of gamma-ray source populations. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Evidence of 100 TeV $γ$-ray emission from HESS J1702-420: A new PeVatron candidate

Authors: H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun, P. Brun, M. Bryan , et al. (211 additional authors not shown)

Abstract: The identification of PeVatrons, hadronic particle accelerators reaching the knee of the cosmic ray spectrum (few $10^{15}$ eV), is crucial to understand the origin of cosmic rays in the Galaxy. We provide an update on the unidentified source HESS J1702-420, a promising PeVatron candidate. We present new observations of HESS J1702-420 made with the High Energy Stereoscopic System (H.E.S.S.), and p… ▽ More The identification of PeVatrons, hadronic particle accelerators reaching the knee of the cosmic ray spectrum (few $10^{15}$ eV), is crucial to understand the origin of cosmic rays in the Galaxy. We provide an update on the unidentified source HESS J1702-420, a promising PeVatron candidate. We present new observations of HESS J1702-420 made with the High Energy Stereoscopic System (H.E.S.S.), and processed using improved analysis techniques. The analysis configuration was optimized to enhance the collection area at the highest energies. We applied a three-dimensional (3D) likelihood analysis to model the source region and adjust non thermal radiative spectral models to the $γ$-ray data. We also analyzed archival data from the Fermi Large Area Telescope (LAT) to constrain the source spectrum at $γ$-ray energies >10 GeV. We report the detection of a new source component called HESS J1702-420A, that was separated from the bulk of TeV emission at a $5.4σ$ confidence level. The power law $γ$-ray spectrum of HESS J1702-420A extends with an index of $Γ=1.53\pm0.19_\text{stat}\pm0.20_\text{sys}$ and without curvature up to the energy band 64-113 TeV, in which it was detected by H.E.S.S. at a $4.0σ$ confidence level. This brings evidence for the source emission up to $100\,\text{TeV}$, which makes HESS J1702-420A a compelling candidate site for the presence of extremely high energy cosmic rays. Remarkably, in a hadronic scenario, the cut-off energy of the proton distribution powering HESS J1702-420A is found to be higher than 0.5 PeV at a 95% confidence level. HESS J1702-420A becomes therefore one of the most solid PeVatron candidates detected so far in H.E.S.S. data, altough a leptonic origin of its emission could not be ruled out either. △ Less

Submitted 14 June, 2021; v1 submitted 11 June, 2021; originally announced June 2021.

Comments: Accepted for publication in the 2. Astrophysical processes section of Astronomy & Astrophysics

Origin and role of relativistic cosmic particles

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

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

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

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

Constraining the origin of the puzzling source HESS J1640-465 and the PeVatron candidate HESS J1641-463 using Fermi-LAT observations

Authors: A. Mares, M. Lemoine-Goumard, F. Acero, C. J. Clark, J. Devin, S. Gabici, J. D. Gelfand, D. A. Green, M. -H. Grondin

Abstract: There are only few very-high-energy sources in our Galaxy which might accelerate particles up to the knee of the cosmic-ray spectrum. To understand the mechanisms of particle acceleration in these PeVatron candidates, \textit{Fermi}-LAT and H.E.S.S. observations are essential to characterize their $γ$-ray emission. HESS J1640$-$465 and the PeVatron candidate HESS J1641$-$463 are two neighboring (\… ▽ More There are only few very-high-energy sources in our Galaxy which might accelerate particles up to the knee of the cosmic-ray spectrum. To understand the mechanisms of particle acceleration in these PeVatron candidates, \textit{Fermi}-LAT and H.E.S.S. observations are essential to characterize their $γ$-ray emission. HESS J1640$-$465 and the PeVatron candidate HESS J1641$-$463 are two neighboring (\ang[astroang]{0.25}) $γ$-ray sources, spatially coincident with the radio supernova remnants (SNRs) G338.3$-$0.0 and G338.5+0.1. Detected both by H.E.S.S. and \textit{Fermi}-LAT, we present here a morphological and spectral analysis of these two sources using 8 years of \textit{Fermi}-LAT data between 200 \si{\mega\electronvolt} and 1 \si{\tera\electronvolt} with multi-wavelength observations to assess their nature. The morphology of HESS J1640$-$465 is described by a 2D Gaussian ($σ=$ \ang[astroang]{0.053} $\pm$ \ang[astroang]{0.011}$_{stat}$ $ \pm$ \ang[astroang]{0.03}$_{syst}$) and its spectrum is modeled by a power-law with a spectral index $Γ= 1.8\pm0.1_{\rm stat}\pm0.2_{\rm syst}$. HESS J1641$-$463 is detected as a point-like source and its GeV emission is described by a logarithmic-parabola spectrum with $α= 2.7 \pm 0.1_ {\rm stat} \pm 0.2_ {\rm syst} $ and significant curvature of $β= 0.11 \pm 0.03_ {\rm stat} \pm 0.05_ {\rm syst} $. Radio and X-ray flux upper limits were derived. We investigated scenarios to explain their emission, namely the emission from accelerated particles within the SNRs spatially coincident with each source, molecular clouds illuminated by cosmic rays from the close-by SNRs, and a pulsar/PWN origin. Our new \emph{Fermi}-LAT results and the radio and flux X-ray upper limits pose severe constraints on some of these models. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 18 pages 7 figures To be published in ApJ

Discovery of Double-Ring Structure in the Supernova Remnant N103B: Evidence for Bipolar Winds from a Type Ia Supernova Progenitor

Authors: Hiroya Yamaguchi, Fabio Acero, Chuan-Jui Li, You-Hua Chu

Abstract: The geometric structure of supernova remnants (SNR) provides a clue to unveiling the pre-explosion evolution of their progenitors. Here we present an X-ray study of N103B (0509-68.7), a Type Ia SNR in the Large Magellanic Cloud, that is known to be interacting with dense circumstellar matter (CSM). Applying our novel method for feature extraction to deep Chandra observations, we have successfully… ▽ More The geometric structure of supernova remnants (SNR) provides a clue to unveiling the pre-explosion evolution of their progenitors. Here we present an X-ray study of N103B (0509-68.7), a Type Ia SNR in the Large Magellanic Cloud, that is known to be interacting with dense circumstellar matter (CSM). Applying our novel method for feature extraction to deep Chandra observations, we have successfully resolved the CSM, Fe-rich ejecta, and intermediate-mass element (IME) ejecta components, and revealed each of their spatial distribution. Remarkably, the IME ejecta component exhibits a double-ring structure, implying that the SNR expands into an hourglass-shape cavity and thus forms bipolar bubbles of the ejecta. This interpretation is supported by more quantitative spectroscopy that reveals a clear bimodality in the distribution of the ionization state of the IME ejecta. These observational results can be naturally explained if the progenitor binary system had formed a dense CSM torus on the orbital plane prior to the explosion, providing further evidence that the SNR N103B originates from a single-degenerate progenitor. △ Less

Submitted 12 March, 2021; originally announced March 2021.

Comments: 12 pages, 12 figures, accepted for publication in ApJ Letters (originally submitted to ApJ, but offered a transfer to ApJL)

Three-dimensional morphological asymmetries in the ejecta of Cassiopeia A using a component separation method in X-rays

Authors: Adrien Picquenot, Fabio Acero, Tyler Holland-Ashford, Laura A. Lopez, Jérôme Bobin

Abstract: Recent simulations have shown that asymmetries in the ejecta distribution of supernova remnants (SNR) can still reflect asymmetries from the initial supernova explosion. Thus, their study provides a great means to test and constrain model predictions in relation to the distributions of heavy elements or the neutron star kicks, both of which are key to better understanding the explosion mechanisms… ▽ More Recent simulations have shown that asymmetries in the ejecta distribution of supernova remnants (SNR) can still reflect asymmetries from the initial supernova explosion. Thus, their study provides a great means to test and constrain model predictions in relation to the distributions of heavy elements or the neutron star kicks, both of which are key to better understanding the explosion mechanisms in core-collapse supernovae. The use of a novel blind source separation method applied to the megasecond X-ray observations of the well-known Cassiopeia A SNR has revealed maps of the distribution of the ejecta endowed with an unprecedented level of detail and clearly separated from continuum emission. Our method also provides a three-dimensional view of the ejecta by disentangling the red- and blue-shifted spectral components and associated images of the Si, S, Ar, Ca and Fe, providing insights into the morphology of the ejecta distribution in Cassiopeia A. These mappings allow us to thoroughly investigate the asymmetries in the heavy elements distribution and probe simulation predictions about the neutron star kicks and the relative asymmetries between the different elements. We find in our study that most of the ejecta X-ray flux stems from the red-shifted component, suggesting an asymmetry in the explosion. In addition, the red-shifted ejecta can physically be described as a broad, relatively symmetric plume, whereas the blue-shifted ejecta is more similar to a dense knot. The neutron star also moves directly opposite to the red-shifted parts of the ejecta similar to what is seen with 44Ti. Regarding the morphological asymmetries, it appears that heavier elements have more asymmetrical distributions, which confirms predictions made by simulations. This study is a showcase of the capacities of new analysis methods to revisit archival observations to fully exploit their scientific content. △ Less

Submitted 24 May, 2021; v1 submitted 2 February, 2021; originally announced February 2021.

Comments: revised version

Journal ref: A&A Volume 646, February 2021, Article number A82 (for the published version)

The Crab nebula variability at short timescales with the Cherenkov Telescope Array

Authors: Enrique Mestre, Emma de Ona Wilhelmi, Dmitry Khangulyan, Roberta Zanin, Fabio Acero, Diego F. Torres

Abstract: Since 2009, several rapid and bright flares have been observed at high energies (>100 MeV) from the direction of the Crab Nebula. Several hypotheses have been put forward to explain this phenomenon, but the origin is still unclear. The detection of counterparts at higher energies with the next generation of Cherenkov telescopes will be determinant to constrain the underlying emission mechanisms. W… ▽ More Since 2009, several rapid and bright flares have been observed at high energies (>100 MeV) from the direction of the Crab Nebula. Several hypotheses have been put forward to explain this phenomenon, but the origin is still unclear. The detection of counterparts at higher energies with the next generation of Cherenkov telescopes will be determinant to constrain the underlying emission mechanisms. We aim at studying the capability of the Cherenkov Telescope Array (CTA) to explore the physics behind the flares, by performing simulations of the Crab Nebula spectral energy distribution, both in flaring and steady state, for different parameters related to the physical conditions in the nebula. In particular, we explore the data recorded by Fermi during two particular flares that occurred in 2011 and 2013. The expected GeV and TeV gamma-ray emission is derived using different radiation models. The resulting emission is convoluted with the CTA response and tested for detection, obtaining an exclusion region for the space of parameters that rule the different flare emission models. Our simulations show different scenarios that may be favourable for achieving the detection of the flares in Crab with CTA, in different regimes of energy. In particular, we find that observations with low sub-100 GeV energy threshold telescopes could provide the most model-constraining results. △ Less

Submitted 9 October, 2022; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: 11 pages, 9 figures; typos corrected

Journal ref: Volume 501, Issue 1, February 2021, Pages 337-346

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

ALMA CO Observations of Gamma-Ray Supernova Remnant N132D in the Large Magellanic Cloud: Possible Evidence for Shocked Molecular Clouds Illuminated by Cosmic-Ray Protons

Authors: H. Sano, P. P. Plucinsky, A. Bamba, P. Sharda, M. D. Filipovic, C. J. Law, R. Z. E. Alsaberi, Y. Yamane, K. Tokuda, F. Acero, M. Sasaki, J. Vink, T. Inoue, S. Inutsuka, J. Shimoda, K. Tsuge, K. Fujii, F. Voisin, N. Maxted, G. Rowell, T. Onishi, A. Kawamura, N. Mizuno, H. Yamamoto, K. Tachihara , et al. (1 additional authors not shown)

Abstract: N132D is the brightest gamma-ray supernova remnant (SNR) in the Large Magellanic Cloud (LMC). We carried out $^{12}$CO($J$ = 1-0, 3-2) observations toward the SNR using the Atacama Large Millimeter/submillimeter Array (ALMA) and Atacama Submillimeter Telescope Experiment. We find diffuse CO emission not only at the southern edge of the SNR as previously known, but also inside the X-ray shell. We s… ▽ More N132D is the brightest gamma-ray supernova remnant (SNR) in the Large Magellanic Cloud (LMC). We carried out $^{12}$CO($J$ = 1-0, 3-2) observations toward the SNR using the Atacama Large Millimeter/submillimeter Array (ALMA) and Atacama Submillimeter Telescope Experiment. We find diffuse CO emission not only at the southern edge of the SNR as previously known, but also inside the X-ray shell. We spatially resolved nine molecular clouds using ALMA with an angular resolution of $5''$, corresponding to a spatial resolution of $\sim$1 pc at the distance of the LMC. Typical cloud sizes and masses are $\sim$2.0 pc and $\sim$100 $M_\odot$, respectively. High-intensity ratios of CO $J$ = 3-2 / 1-0 $> 1.5$ are seen toward the molecular clouds, indicating that shock-heating has occurred. Spatially resolved X-ray spectroscopy reveals that thermal X-rays in the center of N132D are produced not only behind a molecular cloud, but also in front of it. Considering the absence of a thermal component associated with the forward shock towards one molecular cloud located along the line of sight to the center of the remnant, this suggests that this particular cloud is engulfed by shock waves and is positioned on the near side of remnant. If the hadronic process is the dominant contributor to the gamma-ray emission, the shock-engulfed clouds play a role as targets for cosmic-rays. We estimate the total energy of cosmic-ray protons accelerated in N132D to be $\sim$0.5-$3.8 \times 10^{49}$ erg as a conservative lower limit, which is similar to that observed in Galactic gamma-ray SNRs. △ Less

Submitted 10 November, 2020; v1 submitted 15 July, 2020; originally announced July 2020.

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

Journal ref: The Astrophysical Journal, 902, 53 (2020)

PHEMTO : Polarimetric High Energy Modular Telescope Observatory

Authors: P. Laurent, F. Acero, V. Beckmann, S. Brandt, F. Cangemi, M. Civitani, M. Clavel, A. Coleiro, R. Curado, P. Ferrando, C. Ferrigno, F. Frontera, F. Gastaldello, D. Götz, C. Gouiffès, V. Grinberg, L. Hanlon, D. Hartmann, P. Maggi, F. Marin, A. Meuris, T. Okajima, G. Pareschi, G. W. Pratt, N. Rea , et al. (5 additional authors not shown)

Abstract: With the opening of the X and gamma--ray windows in the sixties, thanks to to sounding rockets and satellite-borne instruments, extremely energetic and violent phenomena were discovered and subsequently found to be ubiquitous in the Universe. Observations in the high energy domain are fundamental for understanding how matter is organized and behaves around black holes; unravelling how these extrem… ▽ More With the opening of the X and gamma--ray windows in the sixties, thanks to to sounding rockets and satellite-borne instruments, extremely energetic and violent phenomena were discovered and subsequently found to be ubiquitous in the Universe. Observations in the high energy domain are fundamental for understanding how matter is organized and behaves around black holes; unravelling how these extreme objects influence their environments on a very large scale; and finding the still elusive obscured massive objects in the centre of galaxies. Other major problems in contemporary astrophysics, such as the understanding of acceleration processes at shocks of all sizes (those of pulsar wind nebulae, supernova remnants, but also at larger scales those of Active Galactic Nuclei radio lobes) in relation to the origin of cosmic-rays, or the definitive characterization of the debated non-thermal X-ray energy content of clusters of galaxies, also requires observations at very high energies. An observatory type medium mission operating from around 1 keV to about 600 keV can provide direct insights into these major questions. The essential characteristics will be coverage of the full energy range by telescopes featuring a large throughput and arc-second resolution optics, coupled to a compact focal plane assembly, with excellent imaging resolution and spectroscopy. In addition, the mission will provide unique polarimetry measurements in the hard X-ray domain, an important new diagnostic tool at energies for which the non-thermal processes dominate. The Polarimetric High-Energy Modular Telescope Observatory (PHEMTO) is designed to have performance several orders of magnitude better than the present hard X-ray instruments. This gives to PHEMTO the improvements in scientific performance needed for a mission in the 2050 era. △ Less

Submitted 22 August, 2019; originally announced August 2019.

Comments: 20 pages, White Paper submitted for High Energy Astrophysics to the Voyage 2050 ESA Call

A novel method for component separation of extended sources in X-ray astronomy

Authors: Adrien Picquenot, Fabio Acero, Jérôme Bobin, Pierre Maggi, Jean Ballet, Gabriel W. Pratt

Abstract: In high-energy astronomy, spectro-imaging instruments such as X-ray detectors allow investigation of the spatial and spectral properties of extended sources including galaxy clusters, galaxies, diffuse interstellar medium, supernova remnants and pulsar wind nebulae. In these sources, each physical component possesses a different spatial and spectral signature, but the components are entangled. Ext… ▽ More In high-energy astronomy, spectro-imaging instruments such as X-ray detectors allow investigation of the spatial and spectral properties of extended sources including galaxy clusters, galaxies, diffuse interstellar medium, supernova remnants and pulsar wind nebulae. In these sources, each physical component possesses a different spatial and spectral signature, but the components are entangled. Extracting the intrinsic spatial and spectral information of the individual components from this data is a challenging task. Current analysis methods do not fully exploit the 2D-1D (x,y,E) nature of the data, as the spatial and spectral information are considered separately. Here we investigate the application of a Blind Source Separation algorithm that jointly exploits the spectral and spatial signatures of each component in order to disentangle them. We explore the capabilities of a new BSS method (General Morphological Component Analysis; GMCA), initially developed to extract an image of the Cosmic Microwave Background from Planck data, in an X-ray context. The performance of GMCA on X-ray data is tested using Monte-Carlo simulations of supernova remnant toy models, designed to represent typical science cases. We find that GMCA is able to separate highly entangled components in X-ray data even in high contrast scenarios, and can extract with high accuracy the spectrum and map of each physical component. A modification is proposed to improve the spectral fidelity in the case of strongly overlapping spatial components, and we investigate a resampling method to derive realistic uncertainties associated to the results of the algorithm. Applying the modified algorithm to the deep Chandra observations of Cassiopeia A, we are able to produce detailed maps of the synchrotron emission at low energies (0.6-2.2 keV), and of the red/blue shifted distributions of a number of elements including Si and Fe K. △ Less

Submitted 17 July, 2019; v1 submitted 24 May, 2019; originally announced May 2019.

Journal ref: A&A 627, A139 (2019)

An improved method for the estimation of the Gumbel distribution parameters

Authors: Rubén Gómez González, M. Isabel Parra, Francisco Javier Acero, Jacinto Martín

Abstract: Usual estimation methods for the parameters of extreme values distribution employ only a few values, wasting a lot of information. More precisely, in the case of the Gumbel distribution, only the block maxima values are used. In this work, we propose a method to seize all the available information in order to increase the accuracy of the estimations. This intent can be achieved by taking advantage… ▽ More Usual estimation methods for the parameters of extreme values distribution employ only a few values, wasting a lot of information. More precisely, in the case of the Gumbel distribution, only the block maxima values are used. In this work, we propose a method to seize all the available information in order to increase the accuracy of the estimations. This intent can be achieved by taking advantage of the existing relationship between the parameters of the baseline distribution, which generates data from the full sample space, and the ones for the limit Gumbel distribution. In this way, an informative prior distribution can be obtained. Different statistical tests are used to compare the behaviour of our method with the standard one, showing that the proposed method performs well when dealing with very shortened available data. The empirical effectiveness of the approach is demonstrated through a simulation study and a case study. Reduction in the credible interval width and enhancement in parameter location show that the results with improved prior adapt to very shortened data better than standard method does. △ Less

Submitted 21 February, 2019; originally announced February 2019.

Comments: 21 pages, 10 figures

Spatially resolved broad-band synchrotron emission from the non-thermal limbs of SN1006

Authors: Jiang-Tao Li, Jean Ballet, Marco Miceli, Ping Zhou, Jacco Vink, Yang Chen, Fabio Acero, Anne Decourchelle, Joel N. Bregman

Abstract: We present ~400ks NuSTAR observations of the northeast (NE) and southwest (SW) non-thermal limbs of the Galactic SNR SN1006. We discovered three sources with X-ray emission detected at >50keV. Two of them are identified as background AGN. We extract the NuSTAR spectra from a few regions along the non-thermal limbs and jointly analyze them with the XMM-Newton spectra and the radio data. The broad-b… ▽ More We present ~400ks NuSTAR observations of the northeast (NE) and southwest (SW) non-thermal limbs of the Galactic SNR SN1006. We discovered three sources with X-ray emission detected at >50keV. Two of them are identified as background AGN. We extract the NuSTAR spectra from a few regions along the non-thermal limbs and jointly analyze them with the XMM-Newton spectra and the radio data. The broad-band radio/X-ray spectra can be well described with a synchrotron emission model from a single population of CR electrons with a power law energy distribution and an exponential cutoff. The power law index of the electron particle distribution function (PDF) is ~1.88-1.95 for both the NE and SW limbs, and we do not find significant evidence for a variation of this index at different energy (curvature). There are significant spatial variations of the synchrotron emission parameters. The highest energy electrons are accelerated in regions with the lowest expansion velocity, which is opposite to what has been found in the Tycho's SNR. In addition to a gradual steepening of synchrotron emission from the center of the non-thermal limbs to larger azimuthal angles, we also find that both the emission spectrum and the PDF are significantly flatter in three regions in the SW limb where the shock encounters higher density ambient medium. The NE limb also shows significantly higher cutoff energy in the PDF than the SW limb. By comparing with the roughly symmetric TeV emission and largely asymmetric GeV emission from the two non-thermal limbs, we conclude that the asymmetry in the ambient medium and magnetic fields may have largely modified the acceleration and emission of CR leptons. △ Less

Submitted 25 July, 2018; originally announced July 2018.

Comments: 21 pages, 14 figures, 7 tables, ApJ in press

Disentangling hadronic from leptonic emission in the composite SNR G326.3$-$1.8

Authors: J. Devin, F. Acero, J. Ballet, J. Schmid

Abstract: G326.3$-$1.8 (also known as MSH 15$-$56) has been detected in radio as a middle-aged composite supernova remnant (SNR) consisting of an SNR shell and a pulsar wind nebula (PWN), which has been crushed by the SNR's reverse shock. Previous $γ$-ray studies of SNR G326.3$-$1.8 revealed bright and extended emission with uncertain origin. Understanding the nature of the $γ$-ray emission allows probing t… ▽ More G326.3$-$1.8 (also known as MSH 15$-$56) has been detected in radio as a middle-aged composite supernova remnant (SNR) consisting of an SNR shell and a pulsar wind nebula (PWN), which has been crushed by the SNR's reverse shock. Previous $γ$-ray studies of SNR G326.3$-$1.8 revealed bright and extended emission with uncertain origin. Understanding the nature of the $γ$-ray emission allows probing the population of high-energy particles (leptons or hadrons) but can be challenging for sources of small angular extent. With the recent $\textit{Fermi}$ Large Area Telescope data release Pass 8, we investigate the morphology of this SNR to disentangle the PWN from the SNR contribution. We perform a morphological and spectral analysis from 300 MeV to 300 GeV. We use the reconstructed events with the best angular resolution to separately investigate the PWN and the SNR emissions, which is crucial to accurately determine the spectral properties of G326.3$-$1.8 and understand its nature. The centroid of the $γ$-ray emission evolves with energy and is spatially coincident with the radio PWN at high energies (E $>$ 3 GeV). The morphological analysis reveals that a model considering two contributions from the SNR and the PWN reproduces the $γ$-ray data better than a single-component model. The associated spectral analysis using power laws shows two distinct spectral features, a softer spectrum for the remnant ($Γ$ = 2.17 $\pm$ 0.06) and a harder spectrum for the PWN ($Γ$ = 1.79 $\pm$ 0.12), consistent with hadronic and leptonic origin for the SNR and the PWN respectively. Focusing on the SNR spectrum, we use one-zone models to derive some physical properties and, in particular, we find that the emission is best explained with a hadronic scenario in which the large target density is provided by radiative shocks in HI clouds struck by the SNR. △ Less

Submitted 28 May, 2018; originally announced May 2018.

Journal ref: A&A 617, A5 (2018)

Discovery of a radio nebula around PSR J0855-4644

Authors: C. Maitra, S. Roy, F. Acero, Y. Gupta

Abstract: We report the discovery of a diffuse radio emission around PSR J0855--4644 using an upgraded GMRT (uGMRT) observation at 1.35 GHz. The radio emission is spatially coincident with the diffuse X-ray pulsar wind nebula (PWN) seen with XMM but is much larger in extent compared to the compact axisymmetric PWN seen with Chandra. The morphology of the emission, with a bright partial ring-like structure a… ▽ More We report the discovery of a diffuse radio emission around PSR J0855--4644 using an upgraded GMRT (uGMRT) observation at 1.35 GHz. The radio emission is spatially coincident with the diffuse X-ray pulsar wind nebula (PWN) seen with XMM but is much larger in extent compared to the compact axisymmetric PWN seen with Chandra. The morphology of the emission, with a bright partial ring-like structure and two faint tail-like features strongly resembles a bow shock nebula, and indicates a velocity of 100 km/s through the ambient medium. We conclude that the emission is most likely to be associated with the radio PWN of PSR J0855-4644. From the integrated flux density, we estimate the energetics of the PWN. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: 4 pages, 2 figures, accepted for publication in MNRAS Letters