Search

The Large Magellanic Cloud (LMC) presents a unique environment for pulsar population studies due to its distinct star formation characteristics and proximity to the Milky Way. As part of the TRAPUM (TRAnsients and PUlsars with MeerKAT) Large Survey Project, we are using the core array of the MeerKAT radio telescope (MeerKAT) to conduct a targeted search of the LMC for radio pulsars at L-band frequencies, 856-1712$\,$MHz. The excellent sensitivity of MeerKAT, coupled with a 2-hour integration time, makes the survey 3 times more sensitive than previous LMC radio pulsar surveys. We report the results from the initial four survey pointings which has resulted in the discovery of seven new radio pulsars, increasing the LMC radio pulsar population by 30 per cent. The pulse periods of these new pulsars range from 278 to 1690$\,$ms, and the highest dispersion measure is 254.20$\,$pc$\,$cm$^{-3}$. We searched for, but did not find any significant pulsed radio emission in a beam centred on the SN$\,$1987A remnant, establishing an upper limit of 6.3$\,{\mu}$Jy on its minimum flux density at 1400$\,$MHz., Comment: 12 pages, 4 figures, 4 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We report new radio timing solutions from a three-year observing campaign conducted with the MeerKAT and Murriyang telescopes for nine Small Magellanic Cloud pulsars, increasing the number of characterised rotation-powered extragalactic pulsars by 40 per cent. We can infer from our determined parameters that the pulsars are seemingly all isolated, that six are ordinary pulsars, and that three of the recent MeerKAT discoveries have a young characteristic age of under 100 kyr and have undergone a spin-up glitch. Two of the sources, PSRs J0040$-$7337 and J0048$-$7317, are energetic young pulsars with spin-down luminosities of the order of 10$^{36}$ erg s$^{-1}$. They both experienced a large glitch, with a change in frequency of about 30 $\mu$Hz, and a frequency derivative change of order $-10^{-14}$ Hz s$^{-1}$. These glitches, the inferred glitch rate, and the properties of these pulsars (including potentially high inter-glitch braking indices) suggest these neutron stars might be Vela-like repeating glitchers and should be closely monitored in the future. The position and energetics of PSR J0048$-$7317 confirm it is powering a new Pulsar Wind Nebula (PWN) detected as a radio continuum source; and similarly the association of PSR J0040$-$7337 with the PWN of Supernova Remnant (SNR) DEM S5 (for which we present a new Chandra image) is strengthened. Finally, PSR J0040$-$7335 is also contained within the same SNR but is a chance superposition. It has also been seen to glitch with a change of frequency of $10^{-2}$ $\mu$Hz. This work more than doubles the characterised population of SMC radio pulsars., Comment: 20 pages, 13 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

The Sculptor Galaxy (NGC 253), located in the Southern Hemisphere, far off the Galactic Plane, has a relatively high star-formation rate of about 7 M$_{\odot}$ yr$^{-1}$ and hosts a young and bright stellar population, including several super star clusters and supernova remnants. It is also the first galaxy, apart from the Milky Way Galaxy to be associated with two giant magnetar flares. As such, it is a potential host of pulsars and/or fast radio bursts in the nearby Universe. The instantaneous sensitivity and multibeam sky coverage offered by MeerKAT therefore make it a favourable target. We searched for pulsars, radio-emitting magnetars, and fast radio bursts in NGC 253 as part of the TRAPUM large survey project with MeerKAT. We did not find any pulsars during a four-hour observation, and derive a flux density limit of 4.4 $\mu$Jy at 1400 MHz, limiting the pseudo-luminosity of the brightest putative pulsar in this galaxy to 54 Jy kpc$^2$. Assuming universality of pulsar populations between galaxies, we estimate that detecting a pulsar as bright as this limit requires NGC 253 to contain a pulsar population of $\gtrsim$20 000. We also did not detect any single pulses and our single pulse search flux density limit is 62 mJy at 1284 MHz. Our search is sensitive enough to have detected any fast radio bursts and radio emission similar to the brighter pulses seen from the magnetar SGR J1935+2154 if they had occurred during our observation., Comment: 6 pages, 3 figures, added citation in section 3, corrected values in section 4, published in MNRAS

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realised the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform $E$- and $B$-fields and a static dipole $B$-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand-Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytic Aristotelian Electrodynamics expectations., Comment: 22 pages, 21 figures, accepted by journal

The Tarantula Nebula in the Large Magellanic Cloud is known for its high star formation activity. At its center lies the young massive star cluster R136, providing a significant amount of the energy that makes the nebula shine so brightly at many wavelengths. Recently, young massive star clusters have been suggested to also efficiently produce high-energy cosmic rays, potentially beyond PeV energies. Here, we report the detection of very-high-energy $\gamma$-ray emission from the direction of R136 with the High Energy Stereoscopic System, achieved through a multicomponent, likelihood-based modeling of the data. This supports the hypothesis that R136 is indeed a very powerful cosmic-ray accelerator. Moreover, from the same analysis, we provide an updated measurement of the $\gamma$-ray emission from 30 Dor C, the only superbubble detected at TeV energies presently. The $\gamma$-ray luminosity above $0.5\,\mathrm{TeV}$ of both sources is $(2-3)\times 10^{35}\,\mathrm{erg}\,\mathrm{s}^{-1}$. This exceeds by more than a factor of 2 the luminosity of HESS J1646$-$458, which is associated with the most massive young star cluster in the Milky Way, Westerlund 1. Furthermore, the $\gamma$-ray emission from each source is extended with a significance of $>3\sigma$ and a Gaussian width of about $30\,\mathrm{pc}$. For 30 Dor C, a connection between the $\gamma$-ray emission and the nonthermal X-ray emission appears likely. Different interpretations of the $\gamma$-ray signal from R136 are discussed., Comment: 10+11 pages, 4+6 figures. Corresponding authors: L. Mohrmann, N. Komin

PSR B1259-63 is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 year period, around an O9.5Ve star. At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of ~100 hours of data taken from $t_p-24$ days to $t_p+127$ days around the system's 2021 periastron passage. We also present the timing and spectral analyses of the source. The VHE light curve in 2021 is consistent with the stacked light curve of all previous observations. Within the light curve, we report a VHE maximum at times coincident with the third X-ray peak first detected in the 2021 X-ray light curve. In the light curve -- although sparsely sampled in this time period -- we see no VHE enhancement during the second disc crossing. In addition, we see no correspondence to the 2021 GeV flare in the VHE light curve. The VHE spectrum obtained from the analysis of the 2021 dataset is best described by a power law of spectral index $\Gamma = 2.65 \pm 0.04_{\text{stat}}$ $\pm 0.04_{\text{sys}}$, a value consistent with the previous H.E.S.S. observations of the source. We report spectral variability with a difference of $\Delta \Gamma = 0.56 ~\pm~ 0.18_{\text{stat}}$ $~\pm~0.10_{\text{sys}}$ at 95% c.l., between sub-periods of the 2021 dataset. We also find a linear correlation between contemporaneous flux values of X-ray and TeV datasets, detected mainly after $t_p+25$ days, suggesting a change in the available energy for non-thermal radiation processes. We detect no significant correlation between GeV and TeV flux points, within the uncertainties of the measurements, from $\sim t_p-23$ days to $\sim t_p+126$ days. This suggests that the GeV and TeV emission originate from different electron populations., Comment: accepted to A&A