PSR J2229+6114: Young Pulsar & PWN Dynamics

• PSR J2229+6114 is a young, rotation-powered pulsar associated with a boomerang-shaped pulsar wind nebula and composite SNR.
• The pulsar's well-measured spin parameters and glitch events, coupled with detailed X-ray and γ-ray pulse profiles, advance our understanding of high-energy pulsar physics.
• Multi-wavelength studies reveal conflicting environmental models and support its potential role as a Galactic PeVatron candidate due to >100 TeV emissions.

PSR J2229+6114 is a young, energetic, rotation-powered pulsar associated with the Boomerang pulsar wind nebula (PWN) and with the composite supernova remnant G106.3+2.7. Across the recent literature it appears simultaneously as a $\gamma$-ray pulsar, a bright non-thermal X-ray source, a driver of extended radio/X-ray structures, a system discussed as a Galactic PeVatron candidate because of $>100$ TeV emission from the surrounding complex, and a glitching neutron star with no measurable radiative response to several spin-up events (Pope et al., 2023, Pétri et al., 2024, Fujita et al., 2021, Xia et al., 13 May 2026).

1. Identification, spin parameters, and reported global properties

PSR J2229+6114 entered the high-energy pulsar literature as the source providing a reliable identification for the previously unidentified EGRET source 3EG 2227+6122. AGILE reported pulsed $\gamma$-ray emission with significance $>4\,\sigma$ by epoch-folding, confirmed with $Z^2_n$ and H-test statistics, and found that the $\gamma$-ray period matches the radio period to better than $10^{-11}\,\mathrm{s}$ (Pilia et al., 2011). Fermi-LAT work subsequently listed J2229+6114 among the new individual $\gamma$-ray pulsars discovered using LAT observations together with radio and X-ray ephemerides, explicitly describing it as a Vela-like pulsar (Parent, 2010).

The spin period is reported as $P=0.0516\,\mathrm{s}$ in the AGILE timing analysis and as a $51.67\,\mathrm{ms}$ spin period in the 2020 NuSTAR observation. Published period derivatives are consistent at the level of $>100$0 or $>100$1. Derived quantities are correspondingly energetic: the spin-down power is given as $>100$2 in the pulse-profile localization study and as $>100$3 in the Boomerang transport model; the surface dipole field is reported as $>100$4, equivalently $>100$5; and the light-cylinder radius is given as $>100$6 (Pétri et al., 2024, Chen et al., 2024, Pilia et al., 2011).

Reported characteristic ages are not uniform across the literature. The pulse-profile localization Letter gives $>100$7, the Boomerang transport model gives $>100$8, and the time-dependent leptonic PWN modeling gives $>100$9 (Pétri et al., 2024, Chen et al., 2024, Sarkar et al., 2022). A comparable divergence appears in the assumed distance: the 2024 spatially dependent Boomerang model adopts $\gamma$0, whereas the 2023 multi-wavelength study states a preference for a much larger distance, $\gamma$1, together with a much lower nebular field (Chen et al., 2024, Pope et al., 2023). This suggests that the pulsar’s basic rotational phenomenology is well constrained, whereas system-scale environmental parameters remain model dependent.

2. Pulse profiles and the localization of non-thermal X-ray emission

In AGILE data, the folded $\gamma$2-ray profile above $\gamma$3 shows a single, relatively narrow peak offset by $\gamma$4 with respect to the main $\gamma$5 radio pulse, with full-width at half-maximum $\gamma$6 in phase and persistence up to $\gamma$7. AGILE further reported a pulsed flux above $\gamma$8 essentially identical to the EGRET value for 3EG 2227+6122, $\gamma$9–$>4\,\sigma$0, and a single power-law spectrum $>4\,\sigma$1 with $>4\,\sigma$2, without a detected cutoff or break below $>4\,\sigma$3 (Pilia et al., 2011).

A more geometrical interpretation was developed in the pulse-profile localization study, which used a three-dimensional force-free dipole magnetosphere. There, the radio emission is assumed to arise in a hollow cone at height $>4\,\sigma$4, while $>4\,\sigma$5-ray emission is placed in the striped-wind current sheet outside the light cylinder. Using the Fermi-LAT double-peaked $>4\,\sigma$6-ray light curve with peak separation $>4\,\sigma$7 and radio/$>4\,\sigma$8 phase lag $>4\,\sigma$9, the work quotes the relation

$Z^2_n$0

and reports a formal $Z^2_n$1 minimization with best-fit $Z^2_n$2 and $Z^2_n$3. More broadly, the preferred ranges are $Z^2_n$4–$Z^2_n$5 and $Z^2_n$6–$Z^2_n$7 (Pétri et al., 2024).

Within that framework, the non-thermal X-ray emission is localized by fitting NICER, NuSTAR, and RXTE pulse profiles with thin shells placed along the separatrix between radii $Z^2_n$8 and $Z^2_n$9. The reported result is that the X-ray signal arises from altitudes between $\gamma$0 and $\gamma$1. Lower altitudes would align the X-rays with the radio peak, contrary to the observed $\gamma$2 rotation lag, while higher altitudes would shift curvature-photon energies below the X-ray band and reduce particle density. The paper therefore supports a stratified picture in which radio emission is associated with $\gamma$3, non-thermal X-rays with $\gamma$4–$\gamma$5, and $\gamma$6-rays with the wind zone outside $\gamma$7 (Pétri et al., 2024).

3. The Boomerang PWN and the surrounding radio/X-ray complex

G106.3$\gamma$82.7 is described as a composite SNR characterized by a boomerang-shaped PWN and by two distinct “head” and “tail” regions in the radio band. In the 2023 multi-wavelength investigation, the Boomerang PWN is defined as the $\gamma$9 region around PSR J2229+6114, examined using archival radio and Chandra data, a new NuSTAR observation from 2020, and $10^{-11}\,\mathrm{s}$0-ray upper limits from Fermi and VERITAS. That NuSTAR observation detected the pulsar’s $10^{-11}\,\mathrm{s}$1 spin period and PWN emission described by a power-law model with $10^{-11}\,\mathrm{s}$2 up to $10^{-11}\,\mathrm{s}$3. Contrary to the previous radio study by Kothes et al. 2006, that work preferred a much lower PWN magnetic field, $10^{-11}\,\mathrm{s}$4, and a larger distance, $10^{-11}\,\mathrm{s}$5, arguing from the non-varying X-ray flux over the last two decades, the energy-dependent X-ray PWN size resulting from synchrotron burn-off, and the multi-wavelength spectral energy distribution. Its SED model suggests that the PWN is currently re-expanding after being compressed by the SNR reverse shock $10^{-11}\,\mathrm{s}$6 years ago, and in that interpretation the head region should be formed by GeV--TeV electrons injected earlier by the pulsar propagating into a low-density environment (Pope et al., 2023).

A different but complementary view comes from the Suzaku study of diffuse emission around the SNR and the adjacent pulsar. Using three overlapping source fields and one reference field, that analysis found diffuse X-ray emission represented either by thermal or non-thermal models, but judged the thermal interpretation implausible because the metal abundance is $10^{-11}\,\mathrm{s}$7. In the combined non-thermal fit, the local absorption is $10^{-11}\,\mathrm{s}$8, the photon index is $10^{-11}\,\mathrm{s}$9, and the summed $\gamma$0–$\gamma$1 diffuse flux is $\gamma$2. The diffuse surface brightness declines monotonically from the East field to the West field, increasing toward PSR J2229+6114 in the same way as the radio emission. The X-ray morphology follows the “head” of the $\gamma$3 radio continuum tail, whereas TeV $\gamma$4-ray emission is concentrated farther down the “tail,” coincident with a molecular cloud (Fujita et al., 2021).

The same Suzaku work emphasizes that the X-ray photon index does not change with distance from the pulsar, indicating that radiative cooling is ineffective and particle diffusion is not extremely slow. It also distinguishes the compact Boomerang PWN, which has a much flatter radio index $\gamma$5 and is confined to $\gamma$6, from the extended radio/X-ray tail beyond $\gamma$7–$\gamma$8 with $\gamma$9–$P=0.0516\,\mathrm{s}$0. In that picture, the extended structure arises from electrons diffusing away from the pulsar rather than from fresh acceleration at the SNR shock (Fujita et al., 2021).

4. TeV and ultrahigh-energy emission, PeVatron status, and competing nebular models

The broad interest in the PSR J2229+6114 system is driven by the fact that very-high-energy emission with $P=0.0516\,\mathrm{s}$1 and ultrahigh-energy emission with $P=0.0516\,\mathrm{s}$2 have been reported from the tail region of G106.3$P=0.0516\,\mathrm{s}$32.7, motivating its discussion as a PeVatron candidate (Pope et al., 2023). In the Suzaku analysis, the spatial anti-correlation between X-rays and TeV emission was interpreted as evidence for different particle populations: X-ray and radio emission appear leptonic, with parent electrons likely originating from the pulsar or its wind nebula, whereas the offset $P=0.0516\,\mathrm{s}$4-ray component appears hadronic because it traces a molecular cloud. There, the parent proton population must reach $P=0.0516\,\mathrm{s}$5, while electrons producing X-rays at $P=0.0516\,\mathrm{s}$6 in $P=0.0516\,\mathrm{s}$7 require energies of order $P=0.0516\,\mathrm{s}$8 and remain below a pure-electron PeVatron interpretation (Fujita et al., 2021).

Time-dependent leptonic modeling of the PWN associated with PSR J2229+6114 has nevertheless shown that a one-zone nebular interpretation can reproduce the multi-frequency SED up to the LHAASO band. In that framework, the best-fit solution adopts braking index $P=0.0516\,\mathrm{s}$9, true age $51.67\,\mathrm{ms}$0, present spin-down power $51.67\,\mathrm{ms}$1, present PWN radius $51.67\,\mathrm{ms}$2, and nebular field $51.67\,\mathrm{ms}$3, with $51.67\,\mathrm{ms}$4. In that model, the radio and X-ray bands arise from synchrotron emission, and the GeV–sub-TeV–PeV emission arises from inverse Compton scattering of CMB photons, reproducing the spectrum up to $51.67\,\mathrm{ms}$5. The same study, however, explicitly identified two caveats: the required PWN radius is much larger than the X-ray “Boomerang” nebula, and the inferred field is close to the large-scale Galactic field (Sarkar et al., 2022).

A more spatially resolved 2024 treatment modeled the Boomerang X-ray emission with convective and diffusive transport. It adopts $51.67\,\mathrm{ms}$6 and the transport equation

$51.67\,\mathrm{ms}$7

Across three transport scenarios, the best-fit core field is weak, with $51.67\,\mathrm{ms}$8 in the range $51.67\,\mathrm{ms}$9–$>100$00 and a radial decline to $>100$01–$>100$02 at $>100$03. Those fits reproduce the rapid decline in X-ray surface brightness from $>100$04 to $>100$05 and the softening of photon index from $>100$06 to $>100$07. The same model predicts that inverse Compton emission from the Boomerang PWN may contribute $>100$08–$>100$09 of the flux of LHAASO J2226+6057 at $>100$10 and up to $>100$11 at $>100$12, while leaving open a substantial hadronic contribution because $>100$13 of the pulsar spin-down power may go into relativistic protons or a thermal bath (Chen et al., 2024).

These nebular interpretations should be distinguished from searches for pulsed TeV emission from the neutron-star magnetosphere. VERITAS accumulated $>100$14 of quality-selected data on PSR J2229+6114 and performed both hard-spectrum and soft-spectrum searches. No pulsed signal was found. The hard-spectrum analysis excludes at $>100$15 confidence emission at the mean flux level of the Vela pulsar, while the soft-spectrum analysis rules out Crab-like emission, and even Crab emission scaled by the ratio of the $>100$16–$>100$17 Fermi-LAT fluxes, at $>100$18 near $>100$19. The resulting differential $>100$20 C.L. upper limits extend from $>100$21 to $>100$22 (Wong et al., 23 Sep 2025). A plausible implication is that the system can remain a strong nebular/UHE source candidate even if it lacks a detectable TeV pulsed component.

5. Glitches, braking index, and the absence of radiative outbursts

PSR J2229+6114 is also a glitching pulsar. Fermi-LAT timing over MJD $>100$23–$>100$24 identified seven glitches, of which three are classified as large with $>100$25, and measured a post-relaxation second frequency derivative $>100$26. Using

$>100$27

the intrinsic braking index is reported as $>100$28. Among the larger events, G5 has $>100$29, G7 has $>100$30, and G8 has $>100$31. Vortex-creep fits to the large glitches infer participating fractional moments of inertia at the level of a few $>100$32 to $>100$33 and characteristic recoupling times of order $>100$34–$>100$35 days, supporting an inner-crust angular-momentum reservoir and post-glitch dynamics consistent with the vortex creep model (Gügercinoğlu et al., 2020).

A subsequent monitoring campaign with daily CHIME/Pulsar observations examined four glitches in PSR J2229+6114 and connected them to nearly contemporaneous X-ray coverage with NICER and NuSTAR. The four events occurred at MJD $>100$36, $>100$37, $>100$38, and $>100$39, with fractional sizes ranging from $>100$40 to $>100$41. The CHIME analysis used daily drift-scan observations in the $>100$42–$>100$43 band and real-time glitch alerts from the CHAMPSS timing pipeline; glitch 4 triggered a NuSTAR target-of-opportunity observation $>100$44 days after the event (Xia et al., 13 May 2026).

The central empirical result of that campaign is negative: there were no measurable changes in radio or X-ray emission. Radio pulse-profile residuals were consistent with white noise, no genuine radio bursts were seen, NICER pulsed flux in $>100$45–$>100$46 showed no significant deviation from the weighted mean $>100$47, and NuSTAR found pre- and post-glitch unabsorbed $>100$48–$>100$49 fluxes of $>100$50 and $>100$51, implying a $>100$52 upper limit of $>100$53 on any flux increase. No convincing X-ray bursts and no significant X-ray pulse-shape changes were found (Xia et al., 13 May 2026).

That null result is explicitly contrasted with the magnetar-like post-glitch X-ray outbursts seen in the high-$>100$54 rotation-powered pulsars PSRs J1846$>100$550258 and J1119$>100$566127, both with $>100$57–$>100$58. For PSR J2229+6114, with $>100$59, the absence of radiative activity despite sizeable glitches is interpreted as support for a magnetic-field-threshold picture in which magnetar-like post-glitch outbursts become common only above $>100$60 (Xia et al., 13 May 2026).

6. Use as a laboratory for high-field spectral tests

Beyond standard pulsar and PWN phenomenology, PSR J2229+6114 has been used as a target in searches for axion-like particles (ALPs) in strong magnetic fields. A NICER analysis used ten observations, for a total cleaned exposure of approximately $>100$61, over the energy range $>100$62–$>100$63. The continuum was modeled with an absorbed power law,

$>100$64

and the best-fit continuum parameters for PSR J2229+6114 were $>100$65 and $>100$66 at $>100$67 (Liu et al., 5 Nov 2025).

The search then applied a sliding-window local fitting method to residuals. No bin showed $>100$68; the maximum observed $>100$69 was $>100$70, judged consistent with statistical fluctuations given the number of bins. Under a small-mixing, uniform-field approximation with effective surface field $>100$71 and path length $>100$72, the resulting $>100$73 C.L. upper limit is

$>100$74

The paper emphasizes that this limit is order-of-magnitude because a real pulsar magnetosphere is at least dipolar and possibly multipolar, but it places PSR J2229+6114 among the neutron stars used to probe relativistic ALPs in a high-magnetic-field regime distinct from helioscope and haloscope searches (Liu et al., 5 Nov 2025).

Taken together, these studies place PSR J2229+6114 at the intersection of several active research programs: high-energy pulse-profile modeling, PWN transport and reverberation, Galactic PeVatron identification, glitch-driven neutron-star interior inference, and strong-field spectral tests. The recurring disagreements over the nebular field strength, system distance, and the relative leptonic and hadronic contributions to the UHE emission are not peripheral details; they are central to how the source is interpreted in current multi-wavelength astrophysics (Pope et al., 2023, Sarkar et al., 2022, Chen et al., 2024).