PSR J1640-4631: Young Energetic Pulsar

• PSR J1640-4631 is a young, energetic X-ray pulsar located within SNR G338.3-0.0, noted for its anomalous braking index of 3.15.
• Multiwavelength observations from NuSTAR, XMM-Newton, and Chandra reveal a pulsar wind nebula associated with the TeV source HESS J1640-465, fueling debates over leptonic and hadronic emission models.
• Detailed phase-coherent timing and analysis of braking mechanisms, including magnetospheric torques, gravitational-wave contributions, and magnetic-field decay, offer insights into its unusual spin-down behavior.

Searching arXiv for PSR J1640-4631 and closely related work to ground the article in the literature. PSR J1640-4631 is a young, energetic, 206 ms X-ray pulsar discovered with NuSTAR in spatial association with the shell-type supernova remnant G338.3-0.0, the X-ray point source and putative pulsar wind nebula previously identified in XMM-Newton and Chandra images, and the luminous TeV gamma-ray source HESS J1640-465 (Gotthelf et al., 2014). It is chiefly distinguished by a phase-coherent timing measurement of the braking index, $n = 3.15 \pm 0.03$, the first robust case in which a pulsar braking index is measured to be greater than 3, thereby placing it at the center of debates over magnetospheric torques, inclination-angle evolution, magnetic-field decay, higher multipoles, and gravitational-wave contributions to neutron-star spin-down (Archibald et al., 2016).

1. Discovery, identification, and basic observables

Gotthelf et al. reported PSR J1640-4631 as a pulsar associated with HESS J1640-465 using NuSTAR X-ray observations. The pulsar lies within the shell of G338.3-0.0, whose diameter is $8'$, and coincides with an X-ray point source and surrounding extended emission interpreted as a pulsar wind nebula. The discovery was based on coherent pulsations found in the $3\text{--}25$ keV band, with follow-up observations confirming rapid spin-down (Gotthelf et al., 2014).

The basic measured parameters reported at discovery established PSR J1640-4631 as a high-$\dot{E}$, young system with an unusually luminous nebular and gamma-ray environment.

Quantity	Reported value
Spin period $P$	$0.206443048(33)$ s
Period derivative $\dot P$	$9.758 \times 10^{-13}$ s s$^{-1}$
Spin-down luminosity $\dot E$	$8'$0 erg s$8'$1
Characteristic age $8'$2	3350 yr
Surface dipole magnetic field $8'$3	$8'$4 G

At a distance of 12 kpc to G338.3-0.0, the $8'$5 TeV luminosity of HESS J1640-465 is reported as $8'$6 of the pulsar’s present $8'$7. A search using five years of Fermi LAT data found no gamma-ray pulsations. These facts immediately connected PSR J1640-4631 to the broader problem of whether the GeV–TeV luminosity of the system is powered primarily by a pulsar wind nebula or by hadronic processes in the supernova remnant shell (Gotthelf et al., 2014).

2. Timing solution and the anomalous braking index

Archibald et al. presented a phase-coherent timing solution for PSR J1640-4631 from NuSTAR X-ray timing observations spanning MJD 56463.0–57298.8. At epoch MJD 56741.00000, they reported

$8'$8

$8'$9

$3\text{--}25$0

with a $3\text{--}25$1 upper limit

$3\text{--}25$2

Using the canonical definition

$3\text{--}25$3

they obtained

$3\text{--}25$4

with timing residuals consistent with a unique phase-connected solution and no evidence of pulse-count ambiguities or significant unmodeled noise (Archibald et al., 2016).

The immediate significance of this measurement is that a pure magnetic dipole rotator in vacuum yields $3\text{--}25$5, whereas most other pulsars with measured braking indices have $3\text{--}25$6. PSR J1640-4631 is therefore exceptional in two senses recorded repeatedly in the literature: it is the highest reliably measured braking index among young pulsars, and it is the only one above 3 in the measured sample summarized by these papers (Archibald et al., 2016).

Archibald et al. also examined whether the result could be a timing artifact. Their simulations of red timing noise found that only $3\text{--}25$7 of simulated cases produced a false $3\text{--}25$8 with the observed significance and noise properties, and that the probability dropped to $3\text{--}25$9 when the constraint on $\dot{E}$0 was included. They could not rule out contamination from an unseen glitch recovery, but the recovery timescale would have to be longer than most yet observed. In the radio band, Parkes observations yielded a $\dot{E}$1 upper limit on the pulsed flux at 1.4 GHz of $\dot{E}$2 mJy, reinforcing the source’s radio-quiet characterization in the observational literature (Archibald et al., 2016).

3. Inclination-angle evolution and magnetospheric interpretations

One class of explanations attributes $\dot{E}$3 to evolution of the magnetic inclination angle. In the plasma-filled magnetosphere model, the spin-down and alignment torques are written as

$\dot{E}$4

$\dot{E}$5

leading to

$\dot{E}$6

Within this framework, the allowed braking-index range is $\dot{E}$7, so the observed $\dot{E}$8 is directly admissible. Solving the $\dot{E}$9–$P$0 relation for PSR J1640-4631 gives two possible inclination angles,

$P$1

with the smaller angle preferred because the X-ray pulse profile is single-peaked and the radio output is weak. The inferred inclination-angle evolution is

$P$2

that is, secular alignment of the magnetic and rotation axes (Ekşi et al., 2016).

This interpretation was sharpened by comparison with the vacuum model. In vacuum,

$P$3

which can reproduce $P$4 only for an almost orthogonal rotator. The single-peaked pulse profile makes that geometry disfavored in the plasma-filled interpretation, and the vacuum explanation is described as incompatible with the observed pulse morphology (Ekşi et al., 2016).

A more explicit internal-mechanism model was proposed by applying the two-dipole model of Hamil et al. In that construction, $P$5 is a weak, centered, spin-axis-aligned dipole associated with the rotation effect of a charged sphere, while $P$6 is a strong, localized, off-centered, mobile dipole associated with ferromagnetically ordered material. The magnetic inclination angle is

$P$7

and the equilibrium position of the mobile dipole is

$P$8

Because the potential energy has a minimum, $P$9 rotates toward equilibrium in a pendulum-like manner. The model calculation indicates that $0.206443048(33)$0 would evolve toward alignment with the spin axis for PSR J1640-4631, thereby decreasing $0.206443048(33)$1 and producing $0.206443048(33)$2. The paper further states that the single peak in the pulse profile favors a relatively low change rate of the magnetic inclination angle, and that model rates in the range $0.206443048(33)$3 to $0.206443048(33)$4 can reproduce the observed braking index (Shi et al., 2019).

4. Alternative torque models and implications for neutron-star interior physics

A second class of explanations invokes additional torques beyond standard dipole braking. One possibility is a combined magnetic-dipole and gravitational-wave spin-down model. In that formulation, pure magnetic dipole braking gives $0.206443048(33)$5, pure gravitational-wave braking gives $0.206443048(33)$6, and a combined model yields an intermediate value. For PSR J1640-4631, the gravitational-wave fraction is

$0.206443048(33)$7

so $0.206443048(33)$8 of the power loss is attributed to gravitational waves in that model. The inferred ellipticity is

$0.206443048(33)$9

and the predicted gravitational-wave strain is

$\dot P$0

approximately a factor of four lower than the pure spin-down-limit estimate. This scenario was judged inaccessible to aLIGO at the relevant frequency, but potentially detectable by the Einstein Telescope with one year of integration (Araujo et al., 2016).

A more restrictive gravitational-wave interpretation was advanced in the low-mass-neutron-star scenario. There, conventional neutron star and low-mass quark star candidates require ellipticities much larger than the theoretical estimated maximum value, while a low-mass neutron star with $\dot P$1 yields a required ellipticity $\dot P$2, approximately equal to the theoretical estimated maximum value. To explain the radio-quiet nature within the vacuum gap model, the inclination angle is then constrained to

$\dot P$3

This proposal therefore links the high braking index and radio quietness to a low-mass neutron star with a large inclination angle, while explicitly stating that future gravitational-wave detection and long-term timing are needed to confirm or confute the scenario (Chen, 2016).

A third family of models attributes $\dot P$4 to magnetic-field decay. In the exponential dipole-decay picture,

$\dot P$5

Applying this to PSR J1640-4631, Gao et al. estimated an initial spin period $\dot P$6 ms, a moment of inertia $\dot P$7, a true age $\dot P$8 yr, a decay timescale $\dot P$9 yr, an initial surface dipole field $9.758 \times 10^{-13}$0 G, and a present field-decay rate $9.758 \times 10^{-13}$1 (Gao et al., 2017).

A further extension combined gravitational-wave emission from magnetic deformation with dipole-field decay and tilt-angle evolution to constrain the number of precession cycles, $9.758 \times 10^{-13}$2. In this framework, a future measurement of the tilt angle $9.758 \times 10^{-13}$3 becomes diagnostically important. The reported conclusion is that $9.758 \times 10^{-13}$4 would be larger than previous estimates unless a tiny angle $9.758 \times 10^{-13}$5 is observed, and that $9.758 \times 10^{-13}$6 would indicate $9.758 \times 10^{-13}$7, at least ten times larger than earlier suggestions. This suggests that PSR J1640-4631 may constrain internal dissipation efficiency and internal magnetic-field geometry, provided $9.758 \times 10^{-13}$8 can eventually be measured (Cheng et al., 2019).

Taken together, these models show that the braking index alone does not uniquely identify the dominant torque. The same observed $9.758 \times 10^{-13}$9 has been used to motivate secular alignment in a plasma-filled magnetosphere, a pendulum-like evolution of an internal mobile dipole, long-term dipole-field decay, mixed dipole–gravitational-wave braking, and more specific low-mass-star hypotheses (Ekşi et al., 2016, Shi et al., 2019, Gao et al., 2017, Araujo et al., 2016, Chen, 2016).

5. Supernova-remnant environment, pulsar wind nebula, and the GeV–TeV emission debate

PSR J1640-4631 sits inside a complex multiwavelength environment in which source association and emission mechanism have been debated extensively. The pulsar is located within the shell of G338.3-0.0 and is spatially coincident with HESS J1640-465, described in one study as the most luminous gamma-ray source associated with a supernova remnant in the Galaxy. HI absorption studies place the SNR and pulsar in the distance range $^{-1}$0 kpc, though several models adopt 12 kpc as a working value (Supan et al., 2016, Gotthelf et al., 2014).

The presence of PSR J1640-4631 naturally favored leptonic interpretations in which a pulsar wind nebula powers the GeV–TeV emission through inverse Compton scattering. Gotthelf et al. used the pulsar energetics to revise evolutionary one-zone PWN models and found that the broadband spectrum of HESS J1640-465 could be accommodated provided the pulsar’s braking index is approximately 2 and its initial spin period was $^{-1}$1 ms. In that model, the true age is $^{-1}$2 yr, the initial spin-down power is $^{-1}$3 erg s$^{-1}$4, the explosion energy is $^{-1}$5 erg, and the ambient ISM density is $^{-1}$6 (Gotthelf et al., 2014).

At the same time, a hadronic interpretation has remained viable. Based on HI and $^{-1}$7 data, one analysis measured the ambient proton density in the G338.3-0.0/HESS J1640-465 region to be $^{-1}$8. Using that density and an updated proton–proton gamma-ray cross-section, the authors found that the total energy in accelerated protons required to fit the gamma-ray data is $^{-1}$9 erg for 8.5 kpc and $\dot E$0 erg for 13 kpc. They concluded that a pure hadronic scenario is feasible, while also stating explicitly that a contribution from the PWN associated with PSR J1640-4631 cannot be discarded by the present data (Supan et al., 2016).

A later reanalysis of eight years of Fermi-LAT Pass 8 data separated the GeV field into two overlapping sources and identified the hard component, “Source B,” as the one spatially matching HESS J1640-465. Its spectrum in the $\dot E$1 GeV band is a power law with photon index

$\dot E$2

and an extended $\dot E$3D Gaussian morphology with

$\dot E$4

The GeV spectrum connects smoothly to the TeV spectrum, and the broadband SED can be fit by a one-zone leptonic model with a broken power-law electron spectrum, an exponential cutoff at $\dot E$5 TeV, and a magnetic field $\dot E$6. That analysis argued that morphology, spectral continuity, and population trends favor a PWN origin powered by PSR J1640-4631 rather than an SNR-shell origin (Xin et al., 2018).

The multiwavelength picture changed again with the 2025 report of diffuse radio emission inside G338.3-0.0 using MeerKAT at 816 MHz and 1.4 GHz. The detected radio emission is centrally peaked, elliptical, and well contained within the SNR shell, spatially overlapping the X-ray PWN and HESS J1640-465. At 816 MHz, the diffuse PWN-region flux density is reported as

$\dot E$7

with semi-major axis $\dot E$8 and semi-minor axis $\dot E$9. The absence of mid- and far-infrared counterparts and of catalogued H II regions argues against a thermal origin, while the morphology and radial profile are described as suggestive of a PWN origin powered by PSR J1640-4631. Time-dependent one-zone modeling of the radio, X-ray, and gamma-ray emission further suggests that the PWN is currently interacting with the reverse shock of its host SNR and may be releasing high-energy $8'$00 into the ISM (Abdelmaguid et al., 26 Aug 2025).

6. Observational constraints, controversies, and future tests

Several recurrent misconceptions are explicitly countered in the literature. First, $8'$01 is not treated as proof of a single mechanism. The observational papers state that mass or magnetic quadrupoles could be important if the braking index is stable, while the theoretical papers show that decreasing inclination angle, field decay, mixed dipole–gravitational-wave torque, or combinations of these can all be constructed to match the measured value (Archibald et al., 2016, Ekşi et al., 2016, Gao et al., 2017, Araujo et al., 2016).

Second, the high braking index is not considered secure against every systematic, even though it is robust against ordinary timing noise in the simulations. The possibility of contamination by an unseen glitch recovery is retained in the discovery timing paper, although the required recovery timescale would have to be longer than most yet observed. This keeps long-baseline timing central to the source’s interpretation (Archibald et al., 2016).

Third, the geometry of the star remains underconstrained. The source is radio quiet to the Parkes limit of $8'$02 mJy at 1.4 GHz, and no gamma-ray pulsations were found in five years of Fermi LAT data. Because of this, the actual inclination angle $8'$03, its secular derivative, and the tilt angle $8'$04 are not directly measured. Multiple theoretical papers therefore identify future gamma-ray, X-ray, polarization, and timing observations as the critical tests that could discriminate among alignment-torque, field-decay, two-dipole, and gravitational-wave scenarios (Gotthelf et al., 2014, Ekşi et al., 2016, Shi et al., 2019, Cheng et al., 2019).

Finally, the broader system remains a live case study in coupled pulsar–PWN–SNR evolution. The new radio detection suggests a reverse-shock interaction stage for the PWN, whereas the older gamma-ray literature leaves room for both hadronic and leptonic contributions. A plausible implication is that PSR J1640-4631 is important not only because of its braking index, but also because it connects neutron-star torque physics to reverse-shock-modified PWN evolution and to the origin of the exceptionally luminous GeV–TeV emission from the G338.3-0.0/HESS J1640-465 complex (Supan et al., 2016, Xin et al., 2018, Abdelmaguid et al., 26 Aug 2025).