PSR B1259-63 Research Articles

Precise Measurements of TeV Halos around Geminga and Monogem Pulsars with HAWC

We present the most precise measurements to date for the spatial extension and energy spectrum of the γ-ray region between a pulsar’s wind nebula and the interstellar medium, better known as the halo, present around Geminga and PSR B0656+14 (Monogem) using ∼2398 days of >1 TeV data collected by the HAWC observatory. We interpret the data using a physically motivated model for the diffuse γ-ray emission generated from positrons and electrons (e±) injected by the pulsar wind nebula and inverse Compton scattering with interstellar radiation fields. We find the morphologies of the regions inside these halos are characterized by an inhibited diffusion that are approximately three orders of magnitudes smaller than the Galactic average. We also obtain the e± emission efficiencies of 6.6% and 5.1%, respectively, for Geminga and Monogem. These results have remarkable consequences for the study of the particle diffusion in the region between the pulsar wind nebulae and the interstellar medium, and for the interpretation of the flux of positrons measured by the AMS-02 experiment above 10 GeV.

Magnetic Field Geometry of the Gamma-Ray Binary PSR B1259–63 Revealed via X-Ray Polarization

Some X-ray binaries containing an energetic pulsar in orbit around a normal star accelerate particles to high energies in the shock cone formed where the pulsar and stellar winds collide. The magnetic field geometry in the acceleration region in such binaries is unknown. We performed the first measurement of the polarization of the X-ray synchrotron emission from a gamma-ray-emitting binary system. We observed PSR B1259–63 with the Imaging X-ray Polarimetry Explorer during an X-ray bright phase following the periastron passage in 2024 June. X-ray polarization is detected with a polarization degree of 8.3% ± 1.5% at a significance of 5.3σ. The X-ray polarization angle is aligned with the axis of the shock cone at the time of the observation. This indicates that the predominant component of the magnetic field in the acceleration region is oriented perpendicular to the shock-cone axis.

A Growing Braking Index and Spin-down Swings for the Pulsar PSR B0540–69

The way pulsars spin down is not understood in detail, but a number of possible physical mechanisms produce a spin-down rate that scales as a power of the rotation rate ( ν̇∝−νn ), with the power-law index n called the braking index. PSR B0540−69 is a pulsar that in 2011, after 16 yr of spinning down with a constant braking index of 2.1, experienced a giant spin-down change and a reduction of its braking index to nearly zero. Here, we show that, following this episode, the braking index monotonically increased during a period of at least 4 yr and stabilized at ∼1.1. We also present an alternative interpretation of a more modest rotational irregularity that occurred in 2023, which was modeled as an anomalous negative step of the rotation rate. Our analysis shows that the 2023 observations can be equally well described as a transient swing of the spin-down rate (lasting ∼65 days), and the Bayesian evidence indicates that this model is strongly preferred.

Discovery of the First Antiglitch Event in the Rotation-powered Pulsar PSR B0540-69

Using data from the Neutron star Interior Composition ExploreR observatory, we identified a permanent spin frequency decrease of Δν = −(1.04 ± 0.07) × 10−7 Hz around MJD 60132 in the rotation-powered pulsar PSR B0540-69, which exhibits a periodic signal at a frequency of ν ∼ 19.6 Hz. This points to an antiglitch event, a sudden decrease of the pulsar’s rotational frequency without any major alteration in the pulse profile or any significant increase of the pulsed flux. Additionally, no burst activity was observed in association with the antiglitch. To date, observations of the few known antiglitches have been made in magnetars or accreting pulsars. This is the first antiglitch detected in a rotation-powered pulsar. Given its radiatively quiet nature, this antiglitch is possibly of internal origin. Therefore, we tentatively frame this event within a proposed mechanism for antiglitches where the partial “evaporation” of the superfluid component leads to an increase in the normal component’s moment of inertia and a decrease in the superfluid one.

Searching for X-ray counterparts of degree-wide TeV halos around middle-aged pulsars with SRG/eROSITA

Context. Extended gamma-ray TeV emission (TeV halos) has been detected around middle-aged pulsars. A proposed model to explain these TeV halos is that electrons from a degree-wide pulsar wind nebula (PWN) get up-scattered by cosmic microwave background photons through inverse Compton processes. However, no X-ray degree-wide faint diffuse PWNe have been found around these middle-aged pulsars in previous X-ray observations. Aims. We performed a search for degree-wide PWNe around Geminga, PSR B0656+14, B0540+23, J0633+0632, and J0631+1036 using data from the first four consecutive Spectrum Roentgen Gamma/eROSITA all-sky surveys. In order to better understand the mechanisms underlying the formation of TeV halos, we investigated the magnetic field strength within an area of 1° in radius centred on those pulsars. Methods. To achieve our goals, we selected a list of suitable candidate pulsars in the eROSITA-DE part of the sky and applied data reduction techniques to process the eROSITA data. We then performed a spatial analysis of the regions around selected pulsars. Results. We do not detect any degree-wide diffuse emission driven by rotation-powered pulsars around Geminga, PSR B0656+14, B0540+23, J0633+0632, or J0631+1036. Indeed, close inspection of the data shows that the pulsars of interest are all embedded in diffuse emission from supernova remnants such as the Monogem Ring or the Rosetta Nebula, while PSR B0540+23 is located ~2.5° away from the Crab pulsar, which is exceptionally bright, causing its X-ray emission to spread over a significant area up to the position of PSR B0540+23 and therefore strongly diminishes our capacity to search for degree-wide bright diffuse X-ray emission. Conclusions. Despite the non-detection of any degree-wide PWN surrounding the analysed pulsars, we set flux upper limits to provide useful information on magnetic field strength and its spatial distribution around those pulsars, providing additional constraints on the proposed theory for the formation of TeV halos around pulsars.

Multiwavelength Pulsations and Surface Temperature Distribution in the Middle-aged Pulsar B1055–52

We present a detailed study of the X-ray emission from PSR B1055–52 using XMM-Newton observations from 2019 and 2000. The phase-integrated X-ray emission from this pulsar is poorly described by existing models of neutron star atmospheres. Instead, we confirm that, similar to other middle-aged pulsars, the best-fitting spectral model consists of two blackbody components, with substantially different temperatures and emitting areas, and a nonthermal component characterized by a power law. Our phase-resolved X-ray spectral analysis using this three-component model reveals variations in the thermal emission parameters with the pulsar’s rotational phase. These variations suggest a nonuniform temperature distribution across the neutron star’s surface, including the cold thermal component and probable hot spot(s). Such a temperature distribution can be caused by external and internal heating processes, likely a combination thereof. We observe very high pulse fractions, 60%–80% in the 0.7–1.5 keV range, dominated by the hot blackbody component. This could be related to temperature nonuniformity and potential beaming effects in an atmosphere. We find indication of a second hot spot that appears at lower energies (0.15–0.3 keV) than the first hot spot (0.5–1.5 keV) in the X-ray light curves and is offset by about half a rotation period. This finding aligns with the nearly orthogonal rotator geometry suggested by radio observations of this interpulse pulsar. If the hot spots are associated with polar caps, a possible explanation for their temperature asymmetry could be an offset magnetic dipole and/or an additional toroidal magnetic field component in the neutron star crust.

First Detection of Polarization in X-Rays for PSR B0540-69 and Its Nebula

We report on X-ray polarization measurements of the extragalactic Crab-like PSR B0540-69 and its Pulsar Wind Nebula (PWN) in the Large Magellanic Cloud, using a ∼850 ks Imaging X-ray Polarimetry Explorer (IXPE) exposure. The PWN is unresolved by IXPE. No statistically significant polarization is detected for the image-averaged data, giving a 99% confidence polarization upper limit (MDP99) of 5.3% in the 2–8 keV energy range. However, a phase-resolved analysis detects polarization for both the nebula and pulsar in the 4–6 keV energy range. For the PWN defined as the off-pulse phases, the polarization degree (PD) of (24.5 ± 5.3)% and polarization angle (PA) of (78.1 ± 6.2)° is detected at 4.6σ significance level, consistent with the PA observed in the optical band. In a single on-pulse window, a hint of polarization is measured at 3.8σ with PD of (50.0 ± 13.1)% and PA of (6.2 ± 7.4)°. A “simultaneous” PSR/PWN analysis finds two bins at the edges of the pulse exceeding 3σ PD significance, with PD of (68 ± 20)% and (62 ± 20)%; intervening bins at 2–3σ significance have lower PD, hinting at additional polarization structure.

Constraining the Orbital Inclination and Companion Properties of Three Black Widow Pulsars Detected by FAST

Black widows (BWs) are millisecond pulsars ablating their companion stars. The out-flowing material from the companion can block the radio emission of the pulsar, resulting in eclipses. In this paper, we construct a model for the radio eclipse by calculating the geometry of the bow shock between the winds of the pulsar and companion, where the shock shapes the eclipsing medium but had not been described in detail in previous works. The model is further used to explain the variations of the flux density and dispersion measure of three BW pulsars (i.e., PSR B1957+20, J2055+3829, and J2051−0827) detected by the Five-hundred-meter Aperture Spherical radio Telescope. Consequently, we constrained the parameters of the three BW systems such as the inclination angles and true anomalies of the observer as well as the mass-loss rates and wind velocity of the companion stars. With the help of these constraints, it is expected that magnetic fields of companion stars and even masses of pulsars could further be determined as some extra observation can be achieved in the future.

The Polarized Cosmic Hand: IXPE Observations of PSR B1509−58/MSH 15−52

We describe IXPE polarization observations of the pulsar wind nebula (PWN) MSH 15−52, the “Cosmic Hand.” We find X-ray polarization across the PWN, with B-field vectors generally aligned with filamentary X-ray structures. High-significance polarization is seen in arcs surrounding the pulsar and toward the end of the “jet,” with polarization degree PD > 70%, thus approaching the maximum allowed synchrotron value. In contrast, the base of the jet has lower polarization, indicating a complex magnetic field at significant angle to the jet axis. We also detect significant polarization from PSR B1509−58 itself. Although only the central pulse phase bin of the pulse has high individual significance, flanking bins provide lower-significance detections and, in conjunction with the X-ray image and radio polarization, can be used to constrain rotating vector model solutions for the pulsar geometry.

X-Ray and Near-Infrared Observations of the Middle-aged Pulsar B1055–52, Its Multiwavelength Spectrum, and Proper Motion* * Based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. Based also on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by

Previous observations of the middle-aged γ-ray, X-ray, and radio pulsar B1055–52 indicated some peculiarities, such as a suspected changing of the X-ray flux and spectral parameters, a large excess of the alleged thermal component of the UV spectrum over the Rayleigh–Jeans extension of the X-ray thermal spectrum, and a possible double break in the nonthermal spectral component between the optical and X-ray bands. We observed PSR B1055–52 with the XMM-Newton observatory in X-rays and the Hubble Space Telescope in near-infrared (NIR). The analysis of the XMM-Newton observations does not support the notion of long-term changes in the X-ray flux and broadband X-ray spectrum of the pulsar. Using an observing mode less affected by background noise than the previous XMM-Newton observations, we constrain the power-law (PL) spectral index as (F ν ∝ ν α ) in the energy band 3–10 keV. From the NIR–optical data we obtain a PL slope α O = −0.24 ± 0.10 for the color index E(B − V) = 0.03 mag. The slopes and fluxes of the NIR–optical and X-ray nonthermal spectra suggest that the NIR through X-ray emission can be described by the same PL and is generated by the same mechanism, unlike the pulsar’s γ-ray emission. The excess of the UV thermal component over the extension of the X-ray thermal component became smaller but did not disappear, indicating a nonuniformity of the bulk surface temperature. The NIR data also enable us to measure the proper motion accurately, with values μ α = 47.5 ± 0.7 mas yr−1 and μ δ = −8.7 ± 0.7 mas yr−1.

From the External to the Internal Dynamics of the Neutron Star: The Exotic Braking Indices of PSR B0540−69

The braking index is of great importance for interpreting the spin-down mechanism of pulsar rotation. The sudden spin-down rate transition of PSR B0540−69, the lowest braking index n = 0.031, and its variations up to 1.2 in its later phases, without glitches or changes in X-ray pulsed flux or shape, are the most enigmatic problem that challenges our understanding of the correlation between pulsar spin-down and magnetospheric emission. Here we discuss the above issue on the external and internal grounds of pulsar dynamics. It is found that the increase of the open field line region of the pulsar magnetosphere would give a plausible explanation for the state transition and the braking index of 0.031 of PSR B0540−69, and changes in the moment of inertia of PSR B0540−69 after the state transition may account for the variable braking indices in its subsequent phases. Results indicate that, on the one hand, a change in the magnetosphere size of a pulsar would influence the external braking torque and have a substantial impact on the observed braking index; and, on the other hand, a sudden change in external torque may trigger the mechanism that could slowly increase the moment of inertia of the pulsar and cause an observable effect on the spin frequency second derivatives. This is hardly explained under the regular glitch hypothesis. In this respect, PSR B0540−69 would be the ideal candidate to study the inside and outside dynamics of a pulsar.

First Ejection from the PSR B1259-63/LS 2883 High Mass Gamma-Ray Binary Detected During the 2021–2024 Binary Cycle

LS 2883/PSR B1259-63 is a high mass, eccentric gamma-ray binary that has previously been observed to eject X-ray emitting material. We report the results of recent Chandra observations near binary apastron in which a new X-ray emitting clump of matter was detected. The clump has a high projected velocity of v ⊥ ≈ 0.07c and hard X-ray spectrum, which fits an absorbed power-law model with Γ = 1.1 ± 0.3. Although clumps with similar velocities and spectra were detected in some of the previous binary cycles, no resolved clumps were seen near apastron in the preceding cycle of 2017–2021.

Radio Study of the Pulsar Wind Nebula Powered by PSR B1706-44

PSR B1706−44 is an energetic gamma-ray pulsar located inside supernova remnant (SNR) G343.1−2.3 and it powers a compact pulsar wind nebula (PWN) that shows torus and jet structure in X-rays. We present a radio study of the PWN using Australia Telescope Compact Array observations at 3, 6, 13, and 21 cm. We found an overall arc-like morphology at 3 and 6 cm, and the “arc” shows two distinct peaks at 6 cm. The radio emission is faint inside the X-ray PWN and only brightens beyond that. We develop a thick torus model with Doppler boosting effect to explain the radio PWN structure. The model suggests a bulk flow speed of ∼0.2c, which could indicate significant deceleration of the flow from the X-ray emitting region. Our polarization result reveals a highly ordered toroidal B field in the PWN. Its origin is unclear given that the supernova reverse shock should have interacted with the PWN. At a larger scale, the 13 and 21 cm radio images detected a semicircular rim and an east-west ridge of G343.1−2.3. We argue that the latter could possibly be a pulsar tail rather than a filament of the SNR, as supported by the flat radio spectrum and the alignment between the magnetic field and its elongation.

Neutron star mass estimates from gamma-ray eclipses in spider millisecond pulsar binaries

Reliable neutron star mass measurements are key to determining the equation of state of cold nuclear matter, but such measurements are rare. Black widows and redbacks are compact binaries consisting of millisecond pulsars and semi-degenerate companion stars. Spectroscopy of the optically bright companions can determine their radial velocities, providing inclination-dependent pulsar mass estimates. Although inclinations can be inferred from subtle features in optical light curves, such estimates may be systematically biased due to incomplete heating models and poorly understood variability. Using data from the Fermi Large Area Telescope, we have searched for gamma-ray eclipses from 49 spider systems, discovering significant eclipses in 7 systems, including the prototypical black widow PSR B1957+20. Gamma-ray eclipses require direct occultation of the pulsar by the companion, and so the detection, or significant exclusion, of a gamma-ray eclipse strictly limits the binary inclination angle, providing new robust, model-independent pulsar mass constraints. For PSR B1957+20, the eclipse implies a much lighter pulsar (1.81 ± 0.07 solar masses) than inferred from optical light curve modelling.

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.

Plasma lensing near the eclipses of the Black Widow pulsar B1957+20

ABSTRACT Recently, several eclipsing millisecond pulsars have been shown to experience strong and apparent weak lensing from the outflow of their ionized companions. Lensing can be a powerful probe of the ionized plasma, with the strongest lenses potentially resolving emission regions of pulsars. Understanding lensing in the ‘laboratory-like’ conditions of an eclipsing pulsar may be analogously applied to fast radio bursts, many of which reside in dense, magnetized environments. We examined variable dispersion measure (DM), absorption, scattering, and flux density in the original Black Widow pulsar PSR B1957+20 through an eclipse at the Arecibo Observatory at $327\, {\rm MHz}$. We discovered clear evidence of the two regimes of lensing, strong, and apparent weak. We show that the flux density variations in the apparently weak lensing regime can be modelled directly from variations of DM, using geometric optics. The mean effective velocities in the ingress, $(954\pm 99)\, {\rm km\, s^{-1}}$, and egress $(604\pm 47)\, {\rm km\, s^{-1}}$ cannot be explained by orbital motions alone, but are consistent with significant outflow velocity of material from the companion. We also show that geometric optics can predict when and where the lensing regime-change between weak and strong occurs, and argue that the apparent weak lensing is due to averaging many images. Our framework can be applied in any source with variable electron columns, measuring their relative velocities and distances. In other eclipsing pulsars, this provides a unique opportunity to measure companion outflow velocity, predict regions of weak and strong lensing, and in principle independently constrain orbital inclinations.

Interstellar Scintillation and Polarization of PSR B0656+14 in the Monogem Ring

High-sensitivity interstellar scintillation and polarization observations of PSR B0656+14 made at three epochs over a year using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) show that the scattering is dominated by two different compact regions. We identify the one nearer to the pulsar with the shell of the Monogem Ring, thereby confirming the association. The other is probably associated with the Local Bubble. We find that the observed position angles of the pulsar spin axis and the spatial velocity are significantly different, with a separation of 19.°3 ± 0.°8, inconsistent with a previously published near-perfect alignment of 1° ± 2°. The two independent scattering regions are clearly defined in the secondary spectra, which show two strong forward parabolic arcs. The arc curvatures imply that the scattering screens corresponding to the outer and inner arcs are located approximately 28 pc from PSR B0656+14 and 185 pc from the Earth, respectively. Comparison of the observed Doppler profiles with electromagnetic simulations shows that both scattering regions are mildly anisotropic. For the outer arc, we estimate the anisotropy A R to be approximately 1.3, with the scattering irregularities aligned parallel to the pulsar velocity. For the outer arc, we compare the observed delay profiles with delay profiles computed from a theoretical strong-scattering model. Our results suggest that the spatial spectrum of the scattering irregularities in the Monogem Ring is flatter than Kolmogorov, but further observations are required to confirm this.

A Quarter Century of Guitar Nebula/Filament Evolution

We have collected a new deep Chandra X-Ray Observatory exposure of PSR B2224+65 and the “Guitar Nebula,” mapping the complex X-ray structure. This is accompanied by a new Hubble Space Telescope (HST) Hα image of the head of the Guitar. Comparing the HST and Chandra structures in four epochs over 25 yr, we constrain the evolution of the TeV particles that light up the filament. Cross-field diffusion appears to be enhanced, likely by the injected particles, behind the filament’s sharp leading edge, explaining the filament width and its evolving surface brightness profile.

TRAPUM upper limits on pulsed radio emission for SMC X-ray pulsar J0058−7218

ABSTRACT The TRAPUM collaboration has used the MeerKAT telescope to conduct a search for pulsed radio emission from the young Small Magellanic Cloud pulsar J0058−7218 located in the supernova remnant IKT 16, following its discovery in X-rays with XMM–Newton. We report no significant detection of dispersed, pulsed radio emission from this source in three 2-h L-band observations using the core dishes of MeerKAT, setting an upper limit of 7.0 $\mu$Jy on its mean flux density at 1284 MHz. This is nearly seven times deeper than previous radio searches for this pulsar in Parkes L-band observations. This suggests that the radio emission of PSR J0058−7218 is not beamed towards Earth or that PSR J0058−7218 is similar to a handful of Pulsar Wind Nebulae systems that have a very low radio efficiency, such as PSR B0540−6919, the Large Magellanic Cloud Crab pulsar analogue. We have also searched for bright, dispersed, single radio pulses and found no candidates above a fluence of 93 mJy ms at 1284 MHz.

Combined significance of spatial coincidence of high energy neutrinos from PSR B1509-58 by Super-Kamiokande and MACRO

In their searches for astrophysical point sources of high energy neutrinos, both the Super-Kamiokande and MACRO neutrino detectors saw the largest angular excess from the same source, viz. PSR B1509-58. We estimate the probability for the observed number of events by both Super-Kamiokande and MACRO to be a chance coincidence due to atmospheric neutrino background. We find that this probability is about 0.4%, corresponding to 2.6σ significance. We also propose some additional tests to ascertain if this excess corresponds to an astrophysical signal or is only a background event.