In this study we quantitatively examine the manner pulsar wind, supernova ejecta, and defunct stellar wind materials distribute and melt together into plerions. We performed 2.5D magnetohydrodynamical simulations of the entire evolution of their stellar surroundings and different scenarios are explored, whether the star dies as a red supergiant and Wolf-Rayet supernova progenitors, and whether it moved with velocity 20 or 40kms(-1) through the interstellar medium. Within the post-explosion, early 10kyr, the H-burning-products rich red supergiant wind only mixes by <= 20 per cent, due to its dense circumstellar medium filling the progenitor's bow shock trail, still unaffected by the supernova blast wave. Wolf-Rayet materials, enhanced in C, N, O elements, distribute circularly for the 35M(circle dot) star moving at 20kms(-1) and oblongly at higher velocities, mixing efficiently up to 80 per cent. Supernova ejecta, filled with Mg, Si, Ca, Ti, and Fe, remain spherical for longer times at 20kms(-1) but form complex patterns at higher progenitor speeds due to earlier interaction with the bow shock, in which they mix more efficiently. The pulsar wind mixing is more efficient for Wolf-Rayet (25 per cent) than red supergiant progenitors (20 per cent). This work reveals that the past evolution of massive stars and their circumstellar environments critically shapes the internal distribution of chemical elements on plerionic supernova remnants, and, therefore, governs the origin of the various emission mechanisms at work therein. This is essential for interpreting multifrequency observations of atomic and molecular spectral lines, such as in optical, infrared, and soft X-rays.

We introduce and apply a methodology based on dynamic time warping to compare the whole set of gamma-ray light curves reported in the Third Fermi-Large Area Telescope Pulsar Catalogue. Our method allows us to quantitatively measure the degree of global similarity between two light curves beyond comparing indicators such as how many peaks there are, what their separation is, and their width and height. Once the morphology of the light curve is showcased via background subtraction, min-max scaler normalization and rotations are considered to take into account that phase 0 is arbitrary, and the level of detail with which light curves of different pulsars appear is revealed. In many cases, their similarity is striking and occurs disregarding any other timing, physical, or spectral property. In particular, some millisecond pulsars and young pulsars share detailed light-curve morphology.

Following a recent detection of TeV radiation by the Large High Altitude Air Shower Observatory and the High-Altitude Water Cherenkov Observatory, coincident with the direction of the microquasar V4641 Sgr, we search for possible GeV emission from this source. We explored the morphology and temporal features of the source as well as two nearby unassociated point sources, which could be a part of the extended structure of V4641 Sgr, and compared results with corresponding X-ray and TeV emissions. No significant gamma-ray signal linked to V4641 Sgr was detected in the energy range of 1–300 GeV. The 95% confidence level upper limits for the flux from the source, assuming both point and extended source models, were 5.38 × 10 ^−13 erg cm ^−2 s ^−1 and 1.12 × 10 ^−12 erg cm ^−2 s ^−1 , respectively. Additionally, no correlation between the gamma-ray light curve and X-ray outbursts was observed.

An increasing number of pulsar wind nebulae (PWNe) are being identified in the TeV band by ground-based Imaging Air Cherenkov Telescopes such that they constitute the dominant source class of Galactic TeV emitters. However, MeV-GeV PWN counterparts are still largely lacking. To date, only a dozen PWNe are identified by the Fermi-Large Area Telescope (LAT) in the MeV-GeV band. Most PWNe are located along the Galactic plane embedded within the prominent, diffuse Galactic gamma-ray emission, which makes these sources difficult to disentangle from the bright diffuse background. We present a systematic search for gamma-ray counterparts to known PWNe in the 300MeV-2TeV energy band using the Fermi-LAT. We target locations of previously identified PWNe that lack detected Fermi-LAT pulsars to minimize associated pulsar contamination. The sample includes 6 previously identified Fermi-LAT PWNe and 8 Fermi-LAT sources associated with PWNe. We report the analysis of 58 regions of interest and classify detected sources as either a likely PWN or a candidate PWN counterpart based on their morphological and spectral characteristics across the broadband spectrum. There are 9 unidentified Fermi-LAT sources that we consider as likely PWN counterparts, which, if confirmed to be PWNe, would greatly increase the PWN population detected by the Fermi-LAT from 12 to 21. The remaining Fermi-LAT detected sources are considered weaker PWN candidates. A second approach in the systematic search for gamma-ray emitting PWNe will involve studying the off-pulse phases of Fermi-LAT pulsars for the presence of an obscured PWN and will be reported in a subsequent paper.

On October 9th, 2022, the brightest gamma-ray burst (GRB) since the first GRB observation in the late sixties was detected by the Fermi-GBM and Swift-BAT telescopes (GRB 221009A). The outstanding characteristics of this GRB triggered extensive follow-up observations of the source across all wavebands, including at very-high-energy (VHE) gamma rays with the Large-Sized Telescope prototype (LST-1) of the upcoming Cherenkov Telescope Array Observatory (CTAO). In this contribution, we present the analysis and results of the LST-1 observation campaign in October 2022, focusing on the data taken under nominal observing conditions and above 200 GeV.

Fast radio bursts (FRBs) represent one of the most intriguing phenomena in modern astrophysics. However, their classification into repeaters and nonrepeaters is challenging. Here, we present the application of the graph theory minimum spanning tree (MST) methodology as an unsupervised classifier of repeater and nonrepeater FRBs. By constructing MSTs based on various combinations of variables, we identify those that lead to MSTs that exhibit a localized high density of repeaters at each side of the node with the largest betweenness centrality. Comparing the separation power of this methodology against known machine learning methods, and with the random expectation results, we assess the efficiency of the MST-based approach to unravel the physical implications behind the graph pattern. We finally propose a list of potential repeater candidates derived from the analysis using the MST.

The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high energy (VHE) gamma-rays, as well as details of the individual observations and light curves. We also conduct phenomenological modelling to investigate the basic source properties. We present the first VHE gamma-ray flare from M87 detected since 2010. The flux above 350 GeV has more than doubled within a period of about 36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE gamma-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and emphasises the need for combined image and spectral modelling.

Recent observations of a few young pulsar wind nebulae (PWNe) have revealed their morphologies in some detail. Given the availability of spatio-spectral-temporal data, we use our multi-zone (1D) leptonic emission code to model the PWNe associated with G29.7-0.3 (Kes 75) and G21.5-0.9 (G21.5) and obtain (by-eye) constraints on additional model parameters compared to spectral-only modelling. Kes 75 is a Galactic composite supernova remnant (SNR) with an embedded pulsar, PSR J1846-0258. X-ray studies reveal rapid expansion of Kes 75 over the past two decades. PWN G21.5 is also a composite SNR, powered by PSR J1833-1034. For Kes 75, we study a sudden plasma bulk speed increase that may be due to the magnetar-like outbursts of the central pulsar. An increase of a few percent in this speed does not result in any significant change in the model outputs. For G21.5, we investigate different diffusion coefficients and pulsar spin-down braking indices. We can reproduce the broadband spectra and X-ray surface brightness profiles for both PWNe, and the expansion rate, flux over different epochs, and X-ray photon index vs epoch and central radius for Kes 75 quite well. The latter three features are also investigated for G21.5. Despite obtaining reasonable fits overall, some discrepancies remain, pointing to further model revision. We find similar values to overlapping parameters between our 1D code and those of an independent 0D dynamical code (TIDE). Future work will incorporate spatial data from various energy wavebands to improve model constraints.

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

We present the first X-ray observation of the energetic millisecond pulsar binary PSR J1431-4715, performed with XMM-Newton and complemented with fast optical multi-band photometry acquired with the ULTRACAM instrument at ESO-NTT. It is found as a faint X-ray source without a significant orbital modulation. This contrasts with the majority of systems that instead display substantial X-ray orbital variability. The X-ray spectrum is dominated by non-thermal emission and, due to the lack of orbital modulation, does not favour an origin in an intrabinary shock between the pulsar and companion star wind. While thermal emission from the neutron star polar cap cannot be excluded in the soft X-rays, the dominance of synchrotron emission favours an origin in the pulsar magnetosphere that we describe at both X-ray and gamma-ray energies with a synchro-curvature model. The optical multi-colour light curve folded at the 10.8 h orbital period is double-humped and dominated by ellipsoidal effects, but also affected by irradiation. The ULTRACAM light curves are fit with several models encompassing direct heating and a cold spot, or heat redistribution after irradiation either through convection or convection plus diffusion. Despite the inability to constrain the best irradiation models, the fits provide consistent system parameters, giving an orbital inclination of 59 +/- 6 degrees and a distance of 3.1 +/- 0.3 kpc. The companion is found to be an F-type star, underfilling its Roche lobe (f(RL) = 73 +/- 4%) with a mass of 0.20 +/- 0.04 M-circle dot, confirming the redback status, but hotter than the majority of redbacks. The stellar dayside and nightside temperatures of 7500 K and 7400 K, respectively, indicate a weak irradiation effect on the companion, likely due to its high intrinsic luminosity. Although the pulsar mass cannot be precisely derived, a heavy (1.8-2.2 M-circle dot) neutron star is favoured.

Context. Pulsar wind nebulae (PWNe) are a source of very high energy radiation that can reach up to tera-electron volts and even peta-electron volts. Our work uses the pulsar tree, a graph theory tool recently presented to analyze the pulsar population and select candidates of interest. Aims. We aim to discover detectable PWNe. We also aim to test to what extent the pulsar tree is able to group detectable PWNe despite only considering the intrinsic properties of pulsars. Methods. We selected four pulsars as tera-electron volt PWNe candidates based on their positions in the pulsar tree. Using observed and assumed ranges of values for relevant parameters, we anticipated the possible spectral energy distributions of the PWNe of four pulsars (PSR J1208-6238, J1341-6220, J1838-0537, and J1844-0346) via a detailed time-dependent leptonic model that was already found to be appropriate for describing almost all other detected nebulae. Results. We estimated the likelihood of detection for the four candidates we studied by comparing the TeV fluxes predicted by the possible models with the sensitivities of different observatories. In doing so, we provide context for analyzing the advantages and caveats of using the pulsar tree position as a marker for properties that go beyond the intrinsic features of pulsars that are considered in producing the pulsar tree.

The recent detection of extended γ-ray emission around middle-aged pulsars is interpreted as inverse-Compton scattering of ambient photons by electron-positron pairs escaping the pulsar wind nebula, which are confined near the system by unclear mechanisms. This emerging population of γ-ray sources was first discovered at TeV energies and remains underexplored in the GeV range. To address this, we conducted a systematic search for extended sources along the Galactic plane using 14 years of Fermi-LAT data above 10 GeV, aiming to identify a number of pulsar halo candidates and extend our view to lower energies. The search covered the inner Galactic plane (| l|≤ 100^∘, | b|≤ 1^∘) and the positions of known TeV sources and bright pulsars, yielding broader astrophysical interest. We found 40 such sources, forming the Second Fermi Galactic Extended Sources Catalog (2FGES), most with 68 1.0^∘ and relatively hard spectra with photon indices below 2.5. We assessed detection robustness using field-specific alternative interstellar emission models and by inspecting significance maps. Noting 13 sources previously known as extended in the 4FGL-DR3 catalog and five dubious sources from complex regions, we report 22 newly detected extended sources above 10 GeV. Of these, 13 coincide with H.E.S.S., HAWC, or LHAASO sources; six coincide with bright pulsars (including four also coincident with TeV sources); six are associated with 4FGL point sources only; and one has no association in the scanned catalogs. Notably, six to eight sources may be related to pulsars as classical pulsar wind nebulae or pulsar halos.

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 $\alpha_{\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 $\alpha_{\rm CRp}$ down to about $\Delta\alpha_{\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 $\tau_{\chi}>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.

Many models of dark matter (DM) are now widely considered and probed intensively with accelerators, underground detectors, and astrophysical experiments. Among the various approaches, high-energy astrophysical observations are extremely useful to complement laboratory searches for some DM candidates. In the near future, the Cherenkov Telescope Array (CTA) should enable us to access much heavier weakly interacting massive particles, as well as a broad range of other DM candidates. In this talk, we describe DM searches with CTA.

During its early evolution, a pulsar wind nebula (PWN) sweeps the inner part of the supernova ejecta and forms a thin massive shell. Later on, when the shell has been reached by the reverse shock of the supernova remnant, the evolution becomes more complex, in most cases reverting the expansion into a compression: this later phase is called "reverberation". Computations done so far to understand this phase have been mostly performed in the thin-shell approximation, where the evolution of the PWN radius is assimilated to that of the swept-up shell under the effect of both the inner pressure from the PWN, and the outer pressure from the supernova remnant. Despite the thin-shell approach seems rather justifiable, its implementations have so far been inaccurate, and its correctness, never tested. The outer pressure was naively assumed to be scaled according to the Sedov solution (or a constant fraction of it) along the entire evolution. The thin-shell assumption itself fails along the process, being the shell no longer thin in comparison with the size of the PWN. Here, through a combination of numerical models, dimensional arguments, and analytic approximations, we present a detailed analysis of the interaction of the PWN with the supernova remnant. We provide a new analytic approximation of the outer pressure, beyond the Sedov solution, and a revised "thin-shell" able to reproduce results from numerical simulations. Finally, we compute the efficiency by which the PWN is compressed during reverberation over a wide population of sources.

A decade has passed since high-energy astrophysical neutrinos have been discovered by IceCube, however the corresponding sources have not been fully identified yet. The reported coincidence of the high-energy IceCube-170922A with the gamma-ray blazar TXS 0506+056 is not enough to claim that blazars are the dominant high-energy neutrino emitters in the Universe. In fact, recently IceCube announced a second correlation with NGC 1068, a nearby Seyfert galaxy, which is significantly different from a gamma-emitting blazars. The hunt for counterparts of the IceCube neutrinos using gamma-ray telescopes started in 2012. Nonetheless, these efforts will continue with the next-generation gamma-ray telescopes, such as the CTA Large Size Telescopes (LSTs) and other telescopes, by means of an improved and revised observation strategy. These new observations will allow us to detect enough sources in order to elucidate the mystery of the neutrino emitters. In this contribution, we introduces the efforts made thus far in the search for gamma-ray counterpart of high-energy IceCube events using the current generation IACTs, focusing on alerts made of multiple neutrinos events, and present an idea to improve in the observational strategies with the CTA LSTs that will become operational in the coming decade. We will discuss how to reduce the bias to gamma-ray emitters in order to search for possible neutrino counterparts.

Despite the advancement in background rejection techniques, observation of the very-high-energy gamma-ray sky by imaging atmospheric Cherenkov telescopes (IACTs) are subject to an irreducible background from gamma-like hadron- or electron-induced air showers. The determination of this residual background is crucial for accurate spectral and spatial measurements. The Cherenkov Telescope Array (CTA) will become the next generation of IACTs. To unveil its full potential, the improved reconstruction performance of CTA needs to be coupled with a reliable background estimate across the entire field of view. This may become especially important in the case of the planned surveys of large areas of the sky. In this contribution we will present pybkgmodel, an open-source python software package developed for CTA. It aims at providing in a consistent way the various background modelling methods, based on the experience from current IACTs such as H.E.S.S, MAGIC, and VERITAS. It is designed as a toolbox allowing a user to easily choose the optimal reconstruction approach for various target regions or a combination of several algorithms. We will introduce the design of the package as well as demonstrate its performance using simulations for the CTA Large-Sized Telescope prototype (LST-1).

The prototype Large-Sized Telescope (LST-1) of the Cherenkov Telescope Array Observatory (CTAO) is in commissioning phase at the Observatorio del Roque de Los Muchachos at 2200 m a.s.l. in La Palma (Canary Islands, Spain). LST-1 is a 23-m diameter telescope and is the first of four that will compose the LST part of the CTAO Northern array. The LST subarray is optimized to provide the best sensitivity for gamma rays in the 20 GeV - 200 GeV energy range. The MAGIC telescopes, which are located approximately 100 m from the LST-1, is operating as a two 17-m telescope stereoscopic system for more than 14 years. LST-1 and MAGIC routinely perform joint observations of gamma-ray sources to exploit the potential of the three-telescope system. This contribution describes the analysis pipeline and evaluates the performance of the system using Monte Carlo simulations and data on the Crab Nebula. The sensitivity achieved during joint observations with MAGIC and LST-1 is about 30% higher than that of MAGIC alone.

The prototype Large-Sized Telescope (LST-1) of the Cherenkov Telescope Array (CTA) was inaugurated on La Palma, Canary Islands, in 2018. Since then, the telescope is in the commissioning phase and takes regular gamma-ray data on astrophysical sources while waiting for other CTA telescopes in La Palma to be constructed. Here we present the status of the commissioning, lessons learned, the telescope performance, and scientific highlights achieved in the last couple of years. The science results include the detection of active galactic nuclei flares, studies of pulsar wind nebulae and pulsars, the detection of a Nova, and searches for gamma-ray emission from gamma-ray bursts. We also present the status of the construction of the three further LST telescopes (LST-2-4) in La Palma and plans for their commissioning.

Clusters of galaxies are the largest gravitationally-bound structures in the Universe. They are composed of galaxies and gas (approximately 15% of the total mass) mostly dark matter (DM, accounts up to 85% of the total mass). If the DM is composed of Weakly Interacting Massive Particles (WIMPs), galaxy clusters represent one of the best targets to search for gamma-ray signals induced by the decay of WIMPs, with masses around the TeV scale. Due to its sensitivity and energy range of operation (from 20 GeV to 300 TeV), the Cherenkov Telescope Array (CTA) Observatory has a unique opportunity to test WIMPs with masses close to the unitarity limit. This will complement the searches for DM from other gamma-ray observatories as well as direct and collider experiments. The CTA Observatory is planning to search for gamma-ray emission, either its origin may be cosmic-ray (CR) or DM related, in the Perseus galaxy cluster during the first years of operation. In this poster, we will present the software created to perform the analysis using the \texttt{ctools} software and the corresponding results.