The orbital parameters of gamma-ray binary PSR~J2032+4127

PSR~J2032+4127 is the only one of gamma-ray binary, that exhibits pulsations in gamma-ray. Previous research has indicated that the pulsar and the Be star MT91 213 orbit each other in a highly eccentric orbit with an extremely long period, with the pulsar reaching its periastron on November 13, 2017. Since its launch, the \fermi{} satellite has been monitoring this pulsar for 16 years, covering the 8 years before and the 8 years after the pulsar passed its periastron. Using these data, we present an analysis of pulse arrival times, and precisely determine the orbital parameters for the first time: the orbital period of $P_{\rm orb} \sim 52.3$ yr, the eccentricity of $e \sim 0.98$, the semimajor axis of $a$sin$i \sim 25.3$ AU, and the orbital inclination of $\sim$ 47.1$^\circ$ – 55.1$^\circ$. We also reveal another small glitch occurred in 2021, MJD $\sim$ 59500.

💡 Research Summary

This paper presents a comprehensive timing analysis of the gamma‑ray binary PSR J2032+4127 using sixteen years of Fermi‑LAT observations (MJD 54682–60949). The authors first reduced the Pass 8 Front+Back events within a 0.8° region of interest and 100 MeV–20 GeV energy range, applying standard zenith‑angle and data‑quality cuts. Photon arrival times were barycentered to the Solar System barycenter using the known sky position of the source.

The data were divided into consecutive 50‑day segments. For each segment a pulsation search was performed using the Pearson χ² statistic, scanning a two‑dimensional grid of spin frequency (ν) and its first derivative (˙ν) to capture the rapid variations induced by orbital motion. The measured ν values displayed a clear sinusoidal modulation, confirming that the dominant effect is the Doppler shift caused by the binary orbit rather than intrinsic spin‑down.

Using the Keplerian formalism, the authors fitted five orbital parameters (orbital period P_orb, eccentricity e, projected semi‑major axis a sin i, longitude of periastron ω, and epoch of periastron T_ω) together with higher‑order spin derivatives (¨ν, ν³, ν⁴, ν⁵). The fitting was carried out with a custom Python implementation of the Kepler equation, which also identified two glitches: the previously reported one at MJD 55810.76 and a new, smaller glitch around MJD ≈ 59500.

To remove the orbital Doppler effect, the derived Keplerian solution was used to apply a binary‑barycentric correction to the photon times. After this correction, a second round of period searches yielded the intrinsic spin frequency evolution, which was used to generate pulse‑profile templates and calculate Times‑of‑Arrival (TOAs) for each 50‑day segment. TOA uncertainties were estimated via 10 000 Monte‑Carlo realizations.

The TOAs were then processed with TEMPO2 using the DD orbital model. An iterative refinement (partially phase‑coherent timing, PPCT, followed by fully phase‑coherent timing, FPCT) allowed the authors to converge on a highly precise orbital solution:

• Orbital period P_orb = 19110.5 ± 2.1 days (≈ 52.3 yr)
• Eccentricity e = 0.979889 ± 0.000003
• Projected semi‑major axis a sin i = 12613.2 ± 1.0 light‑seconds (≈ 25.3 AU)
• Longitude of periastron ω = 29.9680°
• Epoch of periastron T_ω = MJD 58070.7351

From the mass function f_m = ( M_Be sin i )³/(M_NS + M_Be)² = 4π² G⁻¹ (a sin i)³ P_orb⁻², and adopting a canonical neutron‑star mass of 1.4 M⊙ together with the measured Be‑star mass range of 13.1–17.5 M⊙, the orbital inclination i is constrained to 47°–55°. This range is consistent with previous estimates but narrows the allowed geometry substantially.

Glitch analysis revealed:

• Glitch 1 (MJD 55810.76): Δν = 1.9089 × 10⁻⁶ Hz, Δ˙ν = −6.3 × 10⁻¹⁶ Hz s⁻¹
• Glitch 2 (MJD ≈ 59500): Δν = 2.52 × 10⁻⁷ Hz, Δ˙ν = −1.8 × 10⁻¹⁶ Hz s⁻¹

The fractional frequency jumps are 2.74 × 10⁻⁷ and 3.61 × 10⁻⁸ respectively, placing both events among the smallest glitches observed in the pulsar population. The inferred glitch rate for PSR J2032+4127 is ≈ 0.12 yr⁻¹, significantly lower than that of Vela‑like or Crab‑like pulsars.

The authors discuss the astrophysical implications of the refined orbit. The extreme eccentricity and long period mean that the pulsar never enters a sustained accretion or propeller regime; it remains a rotation‑powered pulsar throughout its orbit, similar to PSR B1259‑63. The stable GeV emission observed across periastron, contrasted with variable X‑ray and TeV fluxes, can now be modeled with greater confidence using the precise orbital geometry. The inclination range also informs the viewing angle for multi‑wavelength light‑curve modeling, affecting interpretations of shock geometry between the pulsar wind and the Be‑star’s circumstellar disc.

In summary, by exploiting the continuous, all‑sky monitoring capability of Fermi‑LAT, the authors have achieved the first high‑precision determination of the orbital parameters of PSR J2032+4127 and identified a new small glitch. These results provide essential constraints for theoretical models of particle acceleration and high‑energy emission in long‑period, highly eccentric gamma‑ray binaries, and demonstrate the power of long‑baseline gamma‑ray timing for probing binary dynamics.