skip to main content


Title: Global Kinetic Modeling of the Intrabinary Shock in Spider Pulsars
Abstract

Spider pulsars are compact binary systems composed of a millisecond pulsar and a low-mass companion. The relativistic magnetically dominated pulsar wind impacts onto the companion, ablating it and slowly consuming its atmosphere. The interaction forms an intrabinary shock, a proposed site of particle acceleration. We perform global fully kinetic particle-in-cell simulations of the intrabinary shock, assuming that the pulsar wind consists of plane-parallel stripes of alternating polarity and that the shock wraps around the companion. We find that particles are efficiently accelerated via shock-driven reconnection. We extract first-principles synchrotron spectra and light curves, which are in good agreement with X-ray observations: (1) the synchrotron spectrum is nearly flat,Fν∝ const; (2) when the pulsar spin axis is nearly aligned with the orbital angular momentum, the light curve displays two peaks, just before and after the pulsar eclipse (pulsar superior conjunction), separated in phase by ∼0.8 rad; (3) the peak flux exceeds the one at the inferior conjunction by a factor of 10.

 
more » « less
Award ID(s):
1903412
NSF-PAR ID:
10368899
Author(s) / Creator(s):
;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
933
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 140
Size(s):
["Article No. 140"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    In so-called spider pulsars, the X-ray band is dominated by intrabinary shock (IBS) synchrotron emission. While the double-peaked X-ray light curves from these shocks have been well characterized in several spider systems (both black widows and redbacks), the polarization of this emission is yet to be studied. Motivated by the new polarization capability of the Imaging X-ray Polarization Explorer and the confirmation of highly ordered magnetic fields in pulsar wind nebulae, we model the IBS polarization by employing two potential magnetic field configurations: toroidal magnetic fields imposed by the preshock pulsar wind, and tangential shock-generated fields that follow the postshock flow. We find that if IBSs host ordered magnetic fields, the synchrotron X-rays from spider binaries can display a high degree of polarization (≳50%), while the variation in polarization angle provides a good probe of the binary geometry and the magnetic field structure. Our results encourage observational polarization studies of spider pulsars that can distinguish the proposed magnetic models and that constrain the unique properties of these systems better.

     
    more » « less
  2. Abstract HESS J0632+057 belongs to a rare subclass of binary systems that emit gamma rays above 100 GeV. It stands out for its distinctive high-energy light curve, which features a sharp “primary” peak and broader “secondary” peak. We present the results of contemporaneous observations by NuSTAR and VERITAS during the secondary peak between 2019 December and 2020 February, when the orbital phase ( ϕ ) is between 0.55 and 0.75. NuSTAR detected X-ray spectral evolution, while VERITAS detected TeV emission. We fit a leptonic wind-collision model to the multiwavelength spectra data obtained over the four NuSTAR and VERITAS observations, constraining the pulsar spin-down luminosity and the magnetization parameter at the shock. Despite long-term monitoring of the source from 2019 October to 2020 March, the MDM observatory did not detect significant variation in H α and H β line equivalent widths, an expected signature of Be-disk interaction with the pulsar. Furthermore, fitting folded Swift-XRT light-curve data with an intrabinary shock model constrained the orbital parameters, suggesting two orbital phases (at ϕ D = 0.13 and 0.37), where the pulsar crosses the Be-disk, as well as phases for the periastron ( ϕ 0 = 0.30) and inferior conjunction ( ϕ IFC = 0.75). The broadband X-ray spectra with Swift-XRT and NuSTAR allowed us to measure a higher neutral hydrogen column density at one of the predicted disk-passing phases. 
    more » « less
  3. Abstract We present the discovery of a new optical/X-ray source likely associated with the Fermi γ -ray source 4FGL J1408.6–2917. Its high-amplitude periodic optical variability, large spectroscopic radial-velocity semiamplitude, evidence for optical emission lines and flaring, and X-ray properties together imply the source is probably a new black widow millisecond pulsar binary. We compile the properties of the 41 confirmed and suspected field black widows, finding a median secondary mass of 0.027 ± 0.003 M ⊙ . Considered jointly with the more massive redback millisecond pulsar binaries, we find that the “spider” companion mass distribution remains strongly bimodal, with essentially zero systems having companion masses of between ∼0.07 and 0.1 M ⊙ . X-ray emission from black widows is typically softer and less luminous than in redbacks, consistent with less efficient particle acceleration in the intrabinary shock in black widows, excepting a few systems that appear to have more efficient “redback-like” shocks. Together black widows and redbacks dominate the census of the fastest spinning field millisecond pulsars in binaries with known companion types, making up ≳80% of systems with P spin < 2 ms. Similar to redbacks, the neutron star masses in black widows appear on average significantly larger than the canonical 1.4 M ⊙ , and many of the highest-mass neutron stars claimed to date are black widows with M NS ≳ 2.1 M ⊙ . Both of these observations are consistent with an evolutionary picture where spider millisecond pulsars emerge from short orbital period progenitors that had a lengthy period of mass transfer initiated while the companion was on the main sequence, leading to fast spins and high masses. 
    more » « less
  4. Abstract

    SN 2014C was originally classified as a Type Ib supernova, but at phaseϕ= 127 days, post-explosion strong Hαemission was observed. SN 2014C has since been observed in radio, infrared, optical and X-ray bands. Here we present new optical spectroscopic and photometric data spanningϕ= 947–2494 days post-explosion. We address the evolution of the broadened Hαemission line, as well as broad [Oiii] emission and other lines. We also conduct a parallel analysis of all publicly available multiwavelength data. From our spectra, we find a nearly constant HαFWHM velocity width of ∼2000 km s−1that is significantly lower than that of other broadened atomic transitions (∼3000–7000 km s−1) present in our spectra ([Oi]λ6300; [Oiii]λλ4959, 5007; Heiλ7065; [Caii]λλ7291, 7324). The late radio data demand a fast forward shock (∼10,000 km s−1atϕ= 1700 days) in rarified matter that contrasts with the modest velocity of the Hα. We propose that the infrared flux originates from a toroidal-like structure of hydrogen surrounding the progenitor system, while later emission at other wavelengths (radio, X-ray) likely originates predominantly from the reverse shock in the ejecta and the forward shock in the quasi-spherical progenitor He-wind. We propose that the Hαemission arises in the boundary layer between the ejecta and torus. We also consider the possible roles of a pulsar and a binary companion.

     
    more » « less
  5. Abstract

    We present extensive multifrequency Karl G. Jansky Very Large Array (VLA) and Very Long Baseline Array (VLBA) observations of the radio-bright supernova (SN) IIb SN 2004C that span ∼40–2793 days post-explosion. We interpret the temporal evolution of the radio spectral energy distribution in the context of synchrotron self-absorbed emission from the explosion’s forward shock as it expands in the circumstellar medium (CSM) previously sculpted by the mass-loss history of the stellar progenitor. VLBA observations and modeling of the VLA data point to a blastwave with average velocity ∼0.06cthat carries an energy of ≈1049erg. Our modeling further reveals a flat CSM density profileρCSMR−0.03±0.22up to a break radiusRbr≈ (1.96 ± 0.10) × 1016cm, with a steep density gradient followingρCSMR−2.3±0.5at larger radii. We infer that the flat part of the density profile corresponds to a CSM shell with mass ∼0.021M, and that the progenitor’s effective mass-loss rate varied with time over the range (50–500) × 10−5Myr−1for an adopted wind velocityvw= 1000 km s−1and shock microphysical parametersϵe= 0.1,ϵB= 0.01. These results add to the mounting observational evidence for departures from the traditional single-wind mass-loss scenarios in evolved, massive stars in the centuries leading up to core collapse. Potentially viable scenarios include mass loss powered by gravity waves and/or interaction with a binary companion.

     
    more » « less