We will report recent progress on the on-going NuSTAR observation campaign of 8 TeV-detected pulsar wind nebulae (PWNe). This campaign constitutes a major part of our NuSTAR study of some of the most energetic TeV sources in our galaxy detected by VERITAS and HAWC. NuSTAR is the first focusing X-ray telescope operating above 10 keV in space with sub-arcminute angular resolution. Broad-band X-ray morphology and spectroscopy data, obtained by NuSTAR, allow us to probe sub-PeV electron populations through detecting synchrotron X-ray radiation. Our targets include PeVatron candidates detected by HAWC, the Boomerang nebula, PWNe crushed by supernova remnant shocks (or else relic PWNe) and G0.9+0.1 in the Galactic Center. Combined with Fermi-LAT data and available TeV data, we aim to provide a complete, multi-wavelength view of a diverse class of middle-aged (∼10-100 kyrs old) PWNe. Our NuSTAR analysis detected hard X-ray emission from the Eel and Boomerang PWNe and characterized their broad-band X-ray spectra most accurately. We plan to apply both time-evolution and multi-zone PWNe models to multi-wavelength spectral energy distribution (SED) data over the radio, X-ray, GeV and TeV bands. In this proceeding, we will review our observation campaign and discuss the preliminary results for several PWNe in more detail.

We report on broadband X-ray properties of the Rabbit pulsar wind nebula (PWN) associated with the pulsar PSR J1418−6058 using archival Chandra and XMM-Newton data, as well as a new NuSTAR observation. NuSTAR data above 10 keV allowed us to detect the 110 ms spin period of the pulsar, characterize its hard X-ray pulse profile, and resolve hard X-ray emission from the PWN after removing contamination from the pulsar and other overlapping point sources. The extended PWN was detected up to ∼20 keV and is described well by a power-law model with a photon index Γ ≈ 2. The PWN shape does not vary significantly with energy, and its X-ray spectrum shows no clear evidence of softening away from the pulsar. We modeled the spatial profile of X-ray spectra and broadband spectral energy distribution in the radio to TeV band to infer the physical properties of the PWN. We found that a model with low magnetic field strength (B ∼ 10 μG) and efficient diffusion (D ∼ 1027 cm2 s−1) fits the PWN data well. The extended hard X-ray and TeV emission, associated respectively with synchrotron radiation and inverse Compton scattering by relativistic electrons, suggest that particles are accelerated to very high energies (≳500 TeV), indicating that the Rabbit PWN is a Galactic PeVatron candidate.

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 E max = 2 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.

We present a detailed analysis of broadband X-ray observations of the pulsar PSR J1420−6048 and its wind nebula (PWN) in the Kookaburra region with Chandra, XMM-Newton, and NuSTAR. Using the archival XMM-Newton and new NuSTAR data, we detected 68 ms pulsations of the pulsar and characterized its X-ray pulse profile, which exhibits a sharp spike and a broad bump separated by ∼0.5 in phase. A high-resolution Chandra image revealed a complex morphology of the PWN: a torus-jet structure, a few knots around the torus, one long (∼7′) and two short tails extending in the northwest direction, and a bright diffuse emission region to the south. Spatially integrated Chandra and NuSTAR spectra of the PWN out to 2.′5 are well-described by a power-law model with a photon index Γ ≈ 2. A spatially resolved spectroscopic study, as well as NuSTAR radial profiles of the 3–7 keV and 7–20 keV brightness, showed a hint of spectral softening with increasing distance from the pulsar. A multiwavelength spectral energy distribution (SED) of the source was then obtained by supplementing our X-ray measurements with published radio, Fermi-LAT, and H.E.S.S. data. The SED and radial variations of the X-ray spectrum were fit with a leptonic multizone emission model. Our detailed study of the PWN may be suggestive of (1) particle transport dominated by advection, (2) a low magnetic-field strength (B ∼ 5 μG), and (3) electron acceleration to ∼PeV energies.

We report recent progress on the on-going NuSTAR observational campaign of 8 TeV-detected pulsar wind nebulae (PWNe). This campaign constitutes a major part of our NuSTAR study of some of the most energetic TeV sources in our Galaxy detected by VERITAS and HAWC. NuSTAR is the first focusing X-ray telescope operating above 10 keV in space with sub-arcminute angular resolution. Broad-band X-ray imaging and spectroscopy data, obtained by NuSTAR, allow us to probe sub-PeV electron populations through detecting synchrotron X-ray radiation. Our targets include PeVatron candidates detected by HAWC, the Boomerang nebula, PWNe crushed by supernova remnant shocks (or relic PWNe) and G0.9+0.1 in the Galactic Center. Using Fermi-LAT data and available TeV data, we aim to provide a complete, multi-wavelength view of a diverse class of middle-aged (∼10-100 kyrs old) PWNe. Our NuSTAR analysis detected hard X-ray emission from the Eel and Boomerang PWNe and characterized their broad-band X-ray spectra most accurately. We plan to apply both time-evolution and multi-zone PWN models to multi-wavelength spectral energy distribution (SED) data over the radio, X-ray, GeV and TeV bands. In this proceeding, we will review our observational campaign and discuss the preliminary results for some PWNe.