PSR J2021 Research Articles

First results of the glitching pulsar monitoring programme at the Argentine Institute of Radioastronomy

ABSTRACT We report here on the first results of a systematic monitoring of southern glitching pulsars at the Argentine Institute of Radioastronomy that started in the year 2019. We detected a major glitch in the Vela pulsar (PSR J0835 − 4510) and two small glitches in PSR J1048 − 5832. For each glitch, we present the measurement of glitch parameters by fitting timing residuals. We then make an individual pulse study of Vela in observations before and after the glitch. We selected 6 days of observations around the major glitch on 2021 July 22 and study their statistical properties with machine learning techniques. We use variational autoencoder (VAE) reconstruction of the pulses to separate them clearly from the noise. We perform a study with self-organizing map (SOM) clustering techniques to search for unusual behaviour of the clusters during the days around the glitch not finding notable qualitative changes. We have also detected and confirmed recent glitches in PSR J0742 − 2822 and PSR J1740 − 3015.

Relativistic force-free models of the thermal X-ray emission in millisecond pulsars observed by NICER

ABSTRACTJPSeveral important properties of rotation-powered millisecond pulsars (MSPs), such as their mass-radius ratio, equation of state and magnetic field topology, can be inferred from precise observations and modelling of their X-ray light curves. In the present study, we model the thermal X-ray signals originated in MSPs, all the way from numerically solving the surrounding magnetospheres up to the ray tracing of the emitted photons and the final computation of their light curves and spectra. The magnetosphere is solved by performing general relativistic force-free simulations of a rotating neutron star (NS) endowed with a simple centred dipolar magnetic field, for many different stellar compactness and pulsar misalignments. From these solutions, we derive an emissivity map over the surface of the star, based on the electric currents in the magnetosphere. In particular, the emission regions (ERs) are determined in this model by spacelike four-currents that reach the NS. We show that this assumption, together with the inclusion of the gravitational curvature on the force-free simulations, lead to non-standard ERs facing the closed-zone of the pulsar, in addition to other ERs within the polar caps. The combined X-ray signals from these two kinds of ERs (both antipodal) allow to approximate the non-trivial interpulses found in several MSPs light curves. Our modelled X-ray signals are compared against very accurate NICER observations of four target pulsars: PSR J043-4715, PSR J1231-1411, PSR J2124-3358, and PSR J0030 + 0451; achieving very good simultaneous fits for their light curves and spectral distributions.

High-resolution Radio Study of the Dragonfly Pulsar Wind Nebula Powered by PSR J2021+3651

The Dragonfly Nebula (G75.2+0.1) powered by the young pulsar J2021+3651 is a rare pulsar wind nebula (PWN) that shows double tori and polar jets enclosed by a bow-shock structure in X-rays. We present new radio observations of this source taken with the Very Large Array at 6 GHz. The radio PWN has an overall size about two times as large as the X-ray counterpart, consisting of a bright main-body region in the southwest, a narrow and fainter bridge region in the northeast, and a dark gap in between. The nebula shows a radio spectrum much softer than that of a typical PWN. This could be resulting from compression by the ram pressure as the system travels mildly supersonically in the interstellar medium. Our polarization maps reveal a highly ordered and complex B-field structure. This can be explained by a toroidal field distorted by the pulsar motion.

Discovery of GeV Gamma-Ray Emission from Pulsar Wind Nebula Kes 75 and PSR J1846–0258

We report the detection of gamma-ray emission from pulsar wind nebula (PWN) Kes 75 and PSR J1846−0258. Through modeling the spectral energy distribution incorporating the new Fermi-LAT data, we find that the observed gamma-ray emission is likely a combination of both the PWN and pulsar magnetosphere. The spectral shape of this magnetospheric emission is similar to the γ-ray spectrum of rotation-powered pulsars detected by Fermi-LAT, and the results from our best-fit model suggest that the pulsar’s magnetospheric emission accounts for 1% of the current spin-down luminosity. Prior works attempted to characterize the properties of this system and found a low supernova (SN) explosion energy and low SN ejecta mass. We reanalyze the broadband emission incorporating the new Fermi emission and compare the implications of our results to prior reports. The best-fit gamma-ray emission model suggests a second very hot photon field possibly generated by the stellar wind of a Wolf–Rayet star embedded within the nebula, which supports the low ejecta mass found for the progenitor in prior reports and here in the scenario of binary mass transfer.

GeV Gamma-Ray Emission from Pulsar Wind Nebula HESS J1356-645 with Fermi-Large Area Telescope

HESS J1356-645 is considered to be a pulsar wind nebula (PWN) associated with the pulsar PSR J1357-6429. We reanalyze the GeV gamma-ray emission in the direction of HESS J1356-645 with more than 13 yr of Fermi-Large Area Telescope data. The extended gamma-ray emission above 5 GeV is found to be spatially coincident with HESS J1356-645. The spectrum in the energy range of 1 GeV–1 TeV can be described by a power law with an index of Γ = 1.51 ± 0.10. The broadband spectrum of HESS J1356-645 can be reproduced by a leptonic model with a broken power-law electronic spectrum. In addition, we found evidence that the morphology of the GeV emission from HESS J1356-645 varies with energy; a behavior that is similar to that of the PWN Vela-X. More broadband observations will be helpful to study the energy-dependent characteristics of HESS J1356-645.

Investigating the radiative properties of LHAASO J1908 + 0621

ABSTRACT LHAASO J1908+0621 has recently been detected as a source emitting γ-rays with energies above 100 TeV, and multiband observations show that a break around 1 TeV appears in the γ-ray spectrum. We have reanalysed the GeV γ-ray properties of the 100-TeV source using 14 years of data recorded by the Fermi Large Area Telescope (Fermi-LAT). The spectrum in the energy range range 30–500 GeV has an index of 1.50 ± 0.26, which is much smaller than that detected in TeV γ-rays. Additionally, the radiation properties of this source are investigated based on a one-zone time-dependent model. In the model, LHAASO J1908+0621 is associated with a pulsar wind nebula (PWN) powered by the pulsar PSR J1907 + 0602. High-energy particles composed of electrons and positrons are injected into the nebula. Multiband non-thermal emission is produced via synchrotron radiation and inverse Compton scattering (ICS). Taking the effect of radiative energy losses and adiabatic cooling into account, the spectral energy distribution from a model with a broken power law for the distribution of the injected particles can explain the fluxes detected in the γ-ray bands. The results support the idea that LHAASO J1908 + 0621 originates from the PWN powered by PSR J1907 + 0602, and γ-rays with energy above 100 TeV are produced by electrons/positrons in the nebula via ICS.

Bursts from High-magnetic-field Pulsars Swift J1818.0-1607 and PSR J1846.4-0258

The detection of magnetar-like bursts from highly magnetic (B > 1013 G) rotation-powered pulsars (RPPs) opened the magnetar population to yet another group of neutron stars. At the same time the question arose as to whether magnetar-like bursts from high-B RPPs have similar characteristics to bursts from known magnetar sources. We present here our analyses of the Fermi Gamma-ray Burst Monitor (GBM) data from two magnetar candidates, Swift J1818.0−1607 (a radio-loud magnetar) and PSR J1846.4−0258. Both sources entered active bursting episodes in 2020 triggering Fermi-GBM in 2020 and in early 2021. We searched for untriggered bursts from both sources and performed temporal and spectral analyses on all events. Here, we present the results of our comprehensive burst search and analyses. We identified 37 and 58 bursts that likely originated from Swift J1818.0−1607 and PSR J1846.4−0258, respectively. We find that the bursts from these sources are shorter on average than typical magnetar bursts. In addition, their spectra are best described with a single blackbody function with kT ∼ 10–11 keV; several relatively bright events, however, show higher energy emission that could be modeled with a cutoff power-law model. We find that the correlation between the blackbody emitting area and the spectral temperature for the burst ensemble of each pulsar deviates from the ideal Stefan–Boltzmann law, as it does for some burst-active magnetars. We interpret this characteristic as being due to the significant radiation anisotropy expected from optically thick plasmas in very strong magnetic fields.

An update on Fermi-LAT transients in the Galactic plane, including strong activity of Cygnus X-3 in mid-2020

ABSTRACT We present a search for Galactic transient γ-ray sources using 13 yr of the Fermi Large Area Telescope data. The search is based on a recently developed variable-size sliding-time-window (VSSTW) analysis and aimed at studying variable γ-ray emission from binary systems, including novae, γ-ray binaries, and microquasars. Compared to the previous search for transient sources at random positions in the sky with 11.5 yr of data, we included γ-rays with energies down to 500 MeV, increased a number of test positions, and extended the data set by adding data collected between 2020 February and 2021 July. These refinements allowed us to detect additional three novae, V1324 Sco, V5855 Sgr, V357 Mus, and one γ-ray binary, PSR B1259-63, with the VSSTW method. Our search revealed a γ-ray flare from the microquasar, Cygnus X-3, occurred in 2020. When applied to equal quarters of the data, the analysis provided us with detections of repeating signals from PSR B1259-63, LS I +61°303, PSR J2021+4026, and Cygnus X-3. While the Cygnus X-3 was bright in γ-rays in mid-2020, it was in a soft X-ray state and we found that its γ-ray emission was modulated with the orbital period.

Exploring the Hadronic Origin of LHAASO J1908+0621

Recent observations by the Large High Altitude Air Shower Observatory (LHAASO) have paved the way for the observational detection of PeVatrons in the Milky Way, thus revolutionizing the field of γ-ray astrophysics. In this paper, we study one such detected source, LHAASO J1908+0621, and explore the origin of multi-TeV γ-ray emission from this source. A middle-aged radio supernova remnant SNR G40.5–0.5 and a GeV pulsar PSR J1907+0602 are cospatial with LHAASO J1908+0621. Dense molecular clouds are also found to be associated with SNR G40.5–0.5. We explain the multi-TeV γ-ray emission observed from the direction of LHAASO J1908+0621, by the hadronic interaction between accelerated protons that escaped from the SNR shock front and cold protons present inside the dense molecular clouds, and the leptonic emission from the pulsar wind nebula (PWN) associated with the pulsar J1907+0602. Moreover, we explain lower energy γ-ray emission by considering the radiative cooling of the electrons that escaped from SNR G40.5–0.5. Finally, the combined lepto-hadronic scenario was used to explain the multiwavelength spectral energy distribution of LHAASO J1908+0621. Although not yet significant, an IceCube hotspot of neutrino emission is spatially associated with LHAASO J1908+0621, indicating a possible hadronic contribution. In this paper, we show that if a hadronic component is present in LHAASO J1908+0621, then the second-generation IceCube observatory will detect neutrinos from this source.

Single-pulse Variability in Gamma-Ray Pulsars

The Fermi Large Area Telescope receives ≪1 photon per rotation from any γ-ray pulsar. However, out of the billions of monitored rotations of the bright pulsars Vela (PSR J0835−4510) and Geminga (PSR J0633+1746), a few thousand have ≥2 pulsed photons. These rare pairs encode information about the variability of pulse amplitude and shape. We have cataloged such pairs and find the observed number to be in good agreement with simple Poisson statistics, limiting any amplitude variations to <19% (Vela) and <22% (Geminga) at 2σ confidence. Using an array of basis functions to model pulse-shape variability, the observed pulse phase distribution of the pairs limits the scale of pulse-shape variations of Vela to <13%, while for Geminga we find a hint of ∼20% single-pulse-shape variability most associated with the pulse peaks. If variations last longer than a single rotation, more pairs can be collected, and we have calculated upper limits on amplitude and shape variations for assumed coherence times up to 100 rotations, finding limits of ∼1% (amplitude) and ∼3% (shape) for both pulsars. Because a large volume of the pulsar magnetosphere contributes to γ-ray pulse production, we conclude that the magnetospheres of these two energetic pulsars are stable over one rotation and very stable on longer timescales. All other γ-ray pulsars are too faint for similar analyses. These results provide useful constraints on rapidly improving simulations of pulsar magnetospheres, which have revealed a variety of large-scale instabilities in the thin equatorial current sheets where the bulk of GeV γ-ray emission is thought to originate.

Scattering Searches for Dark Matter in Subhalos: Neutron Stars, Cosmic Rays, and Old Rocks.

In many cosmologies dark matter clusters on subkiloparsec scales and forms compact subhalos, in which the majority of Galactic dark matter could reside. Null results in direct detection experiments since their advent four decades ago could then be the result of extremely rare encounters between the Earth and these subhalos. We investigate alternative and promising means to identify subhalo dark matter interacting with standard model particles: (1)subhalo collisions with old neutron stars can transfer kinetic energy and brighten the latter to luminosities within the reach of imminent infrared, optical, and ultraviolet telescopes; we identify new detection strategies involving single-star measurements and Galactic disk surveys, and obtain the first bounds on self-interacting dark matter in subhalos from the coldest known pulsar, PSR J2144-3933; (2)subhalo dark matter scattering with cosmic rays results in detectable effects; (3)historic Earth-subhalo encounters can leave dark matter tracks in Paleolithic minerals deep underground. These searches could discover dark matter subhalos weighing between gigaton and solar masses, with corresponding dark matter cross sections and masses spanning tens of orders of magnitude.

Magnetar XTE J1810–197: Spectro-temporal Evolution of Average Radio Emission

We present the long-term spectro-temporal evolution of the average radio emission properties of the magnetar XTE J1810−197 (PSR J1809−1943), following its most recent outburst in late 2018. We report the results from a 2.5 yr monitoring campaign with the upgraded Giant Metrewave Radio Telescope, carried out over the frequency range of 300–1450 MHz. Our observations show intriguing time variability in the average profile width, flux density, spectral index, and broadband spectral shape. While the average profile width appears to gradually decrease at later epochs, the flux density shows multiple episodes of radio rebrightening over the course of our monitoring. Our systematic monitoring observations reveal that the radio spectrum has steepened over time, resulting in evolution from a magnetar-like spectrum to a more pulsar-like spectrum. A more detailed analysis reveals that the radio spectrum has a turnover, and that this turnover shifts toward lower frequencies with time. We present the details of our analysis leading to these results, and discuss our findings in the context of magnetar radio emission mechanisms, as well as potential manifestations of the intervening medium. We also briefly discuss whether an evolving spectral turnover could be a ubiquitous property of radio magnetars.

The Eel Pulsar Wind Nebula: A PeVatron-candidate Origin for HAWC J1826−128 and HESS J1826−130

HAWC J1826−128 is one of the brightest Galactic TeV γ-ray sources detected by the High Altitude Water Cherenkov (HAWC) observatory, with photon energies extending up to nearly ∼100 TeV. This HAWC source spatially coincides with the H.E.S.S. TeV source HESS J1826−130 and the “Eel” pulsar wind nebula (PWN), which is associated with the GeV pulsar PSR J1826−1256. In the X-ray band, Chandra and XMM-Newton revealed that the Eel PWN is composed of both a compact nebula (∼15″) and diffuse X-ray emission (∼6′ × 2′) extending away from the pulsar. Our NuSTAR observation detected hard X-ray emission from the compact PWN up to ∼20 keV and evidence of the synchrotron burn-off effect. In addition to the spatial coincidence between HESS J1826−130 and the diffuse X-ray PWN, our multiwavelength spectral energy distribution (SED) analysis using X-ray and γ-ray data establishes a leptonic origin of the TeV emission associated with the Eel PWN. Furthermore, our evolutionary PWN SED model suggests (1) a low PWN B-field of ∼1 μG, (2) a significantly younger pulsar age (t ∼ 5.7 kyr) than the characteristic age (τ = 14.4 kyr), and (3) a maximum electron energy of PeV. The low B-field, as well as the putative supersonic motion of the pulsar, may account for the asymmetric morphology of the diffuse X-ray emission. Our results suggest that the Eel PWN may be a leptonic PeVatron particle accelerator powered by the ∼6 kyr old pulsar PSR J1826−1256 with a spin-down power of 3.6 × 1036 erg s−1.

HAWC Study of the Ultra-high-energy Spectrum of MGRO J1908+06

Abstract We report TeV gamma-ray observations of the ultra-high-energy source MGRO J1908+06 using data from the High Altitude Water Cherenkov Observatory. This source is one of the highest-energy known gamma-ray sources, with emission extending past 200 TeV. Modeling suggests that the bulk of the TeV gamma-ray emission is leptonic in nature, driven by the energetic radio-faint pulsar PSR J1907+0602. Depending on what assumptions are included in the model, a hadronic component may also be allowed. Using the results of the modeling, we discuss implications for detection prospects by multi-messenger campaigns.

Modelling the correlated keV/TeV light curves of Be/gamma-ray binaries

Be/γ-ray binaries comprise a confirmed or presumptive pulsar orbiting a Be star and emit luminous γ-rays. Non-thermal emissions are thought to arise from synchrotron radiation and inverse-Compton (IC) scattering in the shock where the pulsar wind is terminated by the stellar outflow. We study wind interactions and shock radiations from such systems and show that the bimodal structures observed in keV/TeV light curves are caused by enhanced synchrotron radiation and IC scattering during disc passages. We use a simple radiation model to reproduce orbital modulations of keV X-ray and TeV γ-ray flux and compare with two confirmed pulsar/Be star binaries (i.e. PSR B1259-63/LS 2883 and PSR J2032+4127/MT91 213), and two candidates (i.e. HESS J0632+057 and LS I +61°303). We find that the keV/TeV light curves of the former two binaries can be well explained by the inclined disc model, while modelling the modulated emissions of the latter two sources remains challenging with current orbital solutions. Therefore, we propose alternative orbital geometries for HESS J0632+057 and LS I +61°303. We estimate the positions and inclination angles of Be discs by fitting correlated keV/TeV light curves. Our results could be beneficial for future measurements of orbital parameters and searches for radio pulsations from presumed pulsars.

Search for Very High-energy Emission from the Millisecond Pulsar PSR J0218+4232

Abstract PSR J0218+4232 is one of the most energetic millisecond pulsars known and has long been considered as one of the best candidates for very high-energy (VHE; &gt;100 GeV) γ-ray emission. Using 11.5 yr of Fermi Large Area Telescope (LAT) data between 100 MeV and 870 GeV, and ∼90 hr of Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in the 20 GeV to 20 TeV range, we searched for the highest energy γ-ray emission from PSR J0218+4232. Based on the analysis of the LAT data, we find evidence for pulsed emission above 25 GeV, but see no evidence for emission above 100 GeV (VHE) with MAGIC. We present the results of searches for γ-ray emission, along with theoretical modeling, to interpret the lack of VHE emission. We conclude that, based on the experimental observations and theoretical modeling, it will remain extremely challenging to detect VHE emission from PSR J0218+4232 with the current generation of Imaging Atmospheric Cherenkov Telescopes, and maybe even with future ones, such as the Cherenkov Telescope Array.

Is the birth of PSR J0538+2817 accompanied by a gamma-ray burst?

ABSTRACT Recently, the Five-hundred-meter Aperture Spherical radio Telescope measured the three-dimensional velocity of PSR J0538+2817 with respect to its associated supernova remnant S147 and found a possible spin–velocity alignment for this pulsar. Here, we show that the high velocity and the spin–velocity alignment of this pulsar can be explained by the so-called electromagnetic rocket mechanism. In this framework, the pulsar is kicked in the direction of the spin axis, which naturally explains the spin–velocity alignment. We scrutinize the evolution of the pulsar and show that the kick process can create a highly relativistic jet at the opposite direction of the kick velocity. The lifetime and energetics of the jet is estimated. It is argued that the jet can generate a gamma-ray burst (GRB). The long-term dynamical evolution of the jet is calculated. It is found that the shock radius of the jet should expand to about 32 pc at present, which is well consistent with the observed radius of the supernova remnant S147 (32.1 ± 4.8 pc). Additionally, our calculations indicate that the current velocity of the GRB remnant should be about 440 km s−1, which is also roughly consistent with the observed blast-wave velocity of the remnant of S147 (500 km s−1).

Supernovae producing unbound binaries and triples

ABSTRACT The fraction of stars that are in binaries or triples at the time of stellar death and the fraction of these systems that survive the supernova explosion are crucial constraints for evolution models and predictions for gravitational wave source populations. These fractions are also subject to direct observational determination. Here, we search 10 supernova remnants containing compact objects with proper motions for unbound binaries or triples using Gaia EDR3 and new statistical methods and tests for false positives. We confirm the one known example of an unbound binary, HD 37424 in G180.0−01.7, and find no other examples. Combining this with our previous searches for bound and unbound binaries, and assuming no bias in favour of finding interacting binaries, we find that 72.0 per cent (52.2–86.4 per cent, 90 per cent confidence) of supernova producing neutron stars are not binaries at the time of explosion, 13.9 per cent (5.4–27.2 per cent) produce bound binaries, and 12.5 per cent (2.8–31.3 per cent) produce unbound binaries. With a strong bias in favour of finding interacting binaries, the medians shift to 76.0 per cent were not binaries at death, 9.5 per cent leave bound binaries, and 13.2 per cent leave unbound binaries. Of explosions that do not leave binaries, ${\lt}18.9{{\ \rm per\ cent}}$ can be fully unbound triples. These limits are conservatively for $M\gt 5\, \mathrm{M}_\odot$ companions, although the mass limits for some individual systems are significantly stronger. At birth, the progenitor of PSR J0538+2817 was probably a 13–$19\, \mathrm{M}_\odot$ star, and at the time of explosion, it was probably a Roche limited, partially stripped star transferring mass to HD 37424 and then producing a Type IIL or IIb supernova.

Radio study of HESS J1857+026: Gamma-rays from a superbubble?

Aims. We provide new insights into the nature of HESS J1857+026, a very-high-energy γ-ray source whose complex morphology in the TeV band was attributed to the superposition of two distinct sources. Methods. We performed radio continuum observations to look for the pulsar wind nebula and the supernova remnant associated with the pulsar PSR J1856+0245, which might be powering part of the γ-ray emission. We observed HESS J1857+026 with the Karl G. Jansky Very Large Array (VLA) at 1.5 GHz in the C configuration. In addition, using the same array configuration, we observed a region of 0.4° × 0.4° towards PSR J1856+0245 at 6.0 GHz. We obtained complementary data for the neutral hydrogen and molecular gas emission from public surveys in order to investigate the properties of the interstellar medium in the direction of HESS J1857+026. Results. The new observations at 1.5 GHz do not show evidence of emission above the noise level of 0.7 mJy beam−1 that could be associated with either HESS J1857+026 or PSR J1856+0245. Also, in the new image at 6.0 GHz we do not detect radio emission from a pulsar wind nebula powered by PSR J1856+0245. The neutral gas analysis shows the existence of a superbubble in the direction of the γ-ray source. We suggest that this structure is located at ~5.5 kpc, compatible with the distance to the pulsar PSR J1856+0245. Conclusions. We conclude that TeV emission from HESS J1857+026 originates in a superbubble, arguing in favour of a single γ-ray source rather than the superposition of two distinct sources. The pulsar PSR J1856+0245 could also be contributing as a source of γ-rays within the bubble.

A Search for TeV Gamma-Ray Emission from Pulsar Tails by VERITAS

Abstract We report on the search for very-high-energy gamma-ray emission from the regions around three nearby supersonic pulsars (PSR B0355+54, PSR J0357+3205, and PSR J1740+1000) that exhibit long X-ray tails. To date there is no clear detection of TeV emission from any pulsar tail that is prominent in X-ray or radio. We provide upper limits on the TeV flux, and luminosity, and also compare these limits with other pulsar wind nebulae detected in X-rays and the tail emission model predictions. We find that at least one of the three tails is likely to be detected in observations that are a factor of 2–3 more sensitive. The analysis presented here also has implications for deriving the properties of pulsar tails, for those pulsars whose tails could be detected in TeV.