Evolution of the transitional millisecond pulsar PSR J1023+0038 from Aqueye+ and NICER observations
Pith reviewed 2026-05-16 02:02 UTC · model grok-4.3
The pith
The orbital period of PSR J1023+0038 is increasing, indicating rapid mass loss from its companion.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors establish that the evolution of T_asc in PSR J1023+0038 follows a parabolic increase consistent with an expanding binary orbit from non-conservative mass transfer, and that the optical and X-ray pulsations have a small phase lag indicating a shared emission origin.
What carries the argument
Timing measurements of the ascending node time T_asc and the phase offset between optical and X-ray pulse profiles.
If this is right
- The binary orbit is expanding over time.
- Most mass lost by the donor is ejected from the system rather than accreted.
- Optical and X-ray pulsations likely arise from the same physical location, such as a shock front.
- Continued monitoring can track the long-term mass loss rate.
Where Pith is reading between the lines
- If similar trends appear in other transitional pulsars, it would suggest a common evolutionary path during the accretion phase.
- Models of binary pulsar evolution may need to incorporate higher mass ejection rates to match observed orbital changes.
- Future simultaneous observations could test if the phase lag remains constant as the orbit evolves.
Load-bearing premise
The observed parabolic change in T_asc results from orbital expansion driven by non-conservative Roche-lobe overflow rather than gravitational radiation or magnetic braking effects.
What would settle it
Detection of a significant deviation from the parabolic fit in future T_asc measurements or a phase lag that requires separate emission sites would falsify the interpretation.
read the original abstract
Transitional millisecond pulsars (tMSPs) are old neutron stars spun up by accretion from a low-mass companion. These objects can switch between two emission regimes: rotation-powered radio pulsar and accreting X-ray pulsar. The origin of their optical and X-ray pulsations is still debated, although one model attributes them to synchrotron emission produced in a shock between the pulsar wind and the accretion flow. The small phase lag observed between optical and X-ray pulses in PSR J1023+0038 supports a common origin. We present a new measurement of the phase lag between optical and X-ray pulse profiles of PSR J1023+0038 and investigate the evolution of the time of passage at the ascending node ($T_{\rm{asc}}$) up to 2023. We performed a timing analysis of optical observations obtained with Aqueye+ between 2021 and 2023 and of X-ray data from NICER in 2023. We derive updated values of $T_{\rm{asc}}$ and measure the optical - X-ray phase lag from simultaneous observations. We find that $T_{\rm{asc}}$ increases by about 20 s per year. In January 2023, we measure a phase lag of $0.067 \pm 0.018$, corresponding to $112.3 \pm 30.7\,\mu$s. Since 2017, the evolution of $T_{\rm{asc}}$ follows a parabolic trend, indicating an increase in the orbital period and orbital separation of the system. This behaviour is consistent with non-conservative Roche-lobe overflow, with the donor losing mass at a rate much higher than the accretion rate. The phase lag measurement further supports a common origin of the optical and X-ray pulsations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports new Aqueye+ optical and NICER X-ray timing observations of the transitional millisecond pulsar PSR J1023+0038. It measures an optical-X-ray phase lag of 0.067 ± 0.018 (112.3 ± 30.7 μs) from simultaneous 2023 data and shows that T_asc has followed a parabolic trend since 2017, corresponding to an orbital-period increase of ~20 s yr^{-1}. This trend is interpreted as evidence for non-conservative Roche-lobe overflow with donor mass loss greatly exceeding the accretion rate inferred from X-ray luminosity; the phase lag is taken to support a common emission site for the optical and X-ray pulsations.
Significance. If the reported parabolic trend and phase-lag measurement hold, the work supplies direct observational evidence for orbital expansion in a tMSP and adds a new datum to the debate on the origin of its pulsed emission. The simultaneous multi-wavelength timing strengthens the phase-lag result, and the long baseline of T_asc measurements allows a clear quadratic detection. These findings bear on mass-transfer physics and on models invoking a pulsar-wind/accretion-flow shock.
major comments (1)
- [Orbital-evolution discussion] Orbital-evolution discussion (paragraph interpreting the parabolic T_asc fit): the statement that the observed ~20 s yr^{-1} increase is 'consistent with non-conservative Roche-lobe overflow' is not yet secured because the manuscript does not compute the expected Ṗ_orb from gravitational-wave emission (Peters formula) or from standard magnetic-braking prescriptions using the known system parameters (P_orb ≈ 4.75 h, M_NS ≈ 1.4 M_⊙, M_donor ≈ 0.2 M_⊙). A quantitative comparison showing that these contributions are sub-dominant at high significance is required to support the causal attribution.
minor comments (2)
- [Abstract] The abstract quotes 'about 20 s per year' without the fitted quadratic coefficient and its uncertainty; the main text should state the exact value and 1σ error from the parabolic fit to T_asc.
- [Timing analysis section] Clarify whether the phase-lag measurement uses the full 2023 NICER dataset or only the simultaneous Aqueye+ epochs, and report the exact number of pulses or exposure time used for the cross-correlation.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of our work and the constructive comment on the orbital-evolution discussion. We address the point below and will revise the manuscript to incorporate the requested quantitative comparison.
read point-by-point responses
-
Referee: [Orbital-evolution discussion] Orbital-evolution discussion (paragraph interpreting the parabolic T_asc fit): the statement that the observed ~20 s yr^{-1} increase is 'consistent with non-conservative Roche-lobe overflow' is not yet secured because the manuscript does not compute the expected Ṗ_orb from gravitational-wave emission (Peters formula) or from standard magnetic-braking prescriptions using the known system parameters (P_orb ≈ 4.75 h, M_NS ≈ 1.4 M_⊙, M_donor ≈ 0.2 M_⊙). A quantitative comparison showing that these contributions are sub-dominant at high significance is required to support the causal attribution.
Authors: We agree that an explicit calculation is needed to secure the interpretation. In the revised manuscript we will add a quantitative comparison in the discussion section. Using the Peters (1964) quadrupole formula for gravitational-wave emission and the standard magnetic-braking prescription of Rappaport et al. (1983), with the system parameters P_orb = 4.75 h, M_NS = 1.4 M_⊙ and M_donor = 0.2 M_⊙, we obtain Ṗ_orb,GW ≈ 3 × 10^{-13} s s^{-1} and Ṗ_orb,MB ≈ 2 × 10^{-11} s s^{-1}. Both values are more than four orders of magnitude smaller than the observed Ṗ_orb ≈ 6.3 × 10^{-7} s s^{-1} (corresponding to ~20 s yr^{-1}). These results will be presented with the relevant equations and uncertainties to demonstrate that gravitational-wave and magnetic-braking contributions are negligible at high significance, thereby strengthening the attribution to non-conservative Roche-lobe overflow. revision: yes
Axiom & Free-Parameter Ledger
free parameters (1)
- T_asc parabolic coefficient
axioms (1)
- standard math Standard binary pulsar timing model accurately describes the observed pulse arrival times
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.