A common four-beam geometry reveals altitude-stratified GeV pulses in canonical young pulsars

Paul K. H. Yeung; Takayuki Saito

arxiv: 2604.28000 · v2 · submitted 2026-04-30 · 🌌 astro-ph.HE · astro-ph.SR

A common four-beam geometry reveals altitude-stratified GeV pulses in canonical young pulsars

Paul K. H. Yeung , Takayuki Saito This is my paper

Pith reviewed 2026-05-07 07:28 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.SR

keywords gamma-ray pulsarsFermi LATpulsar geometryphaseogramsmagnetospheric emissionaltitude stratificationDoppler shiftsCrab Vela Dragonfly

0 comments

The pith

Crab, Vela and Dragonfly share a common four-beam geometry for their gamma-ray pulses with two altitude-separated pairs.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that the gamma-ray phaseograms of three canonical young pulsars fit inside one four-beam geometric template despite their different appearances. The template splits the emission into two pairs at different heights: a lower pair from plasma moving almost with the star's rotation that produces sharp main peaks, and a higher pair from plasma with an outward flow that fills bridges between peaks and adds small ripples. This decomposition uses a model focused on geometry and phase-dependent Doppler shifts rather than specific radiation processes. A reader would care because the approach turns scattered pulse shapes into one picture of how emission height and plasma motion shape what we observe. It also gives concrete heights for the sites and points to different mechanisms operating at each layer.

Core claim

The phaseograms of Crab, Vela and Dragonfly can be organised within a single four-beam geometric template. In each pulsar the phaseogram admits a decomposition into two altitude-separated beam pairs. The lower-altitude pair is produced by plasma with bulk motion close to azimuthal corotation, sharpening the main peaks. The higher-altitude pair shows a radially outward bulk-motion component and contributes bridge/shoulder emission and ripple-like modulations overlapping the main peaks. Higher-altitude site heights vary from about 0.7 light-cylinder radii for the Crab to 1.1-1.4 for Vela and Dragonfly.

What carries the argument

The mechanism-agnostic parametric model that incorporates phase-dependent Doppler shifts to constrain the three-dimensional locations and bulk motions of four emission sites.

If this is right

The lower-altitude beams align with curvature-dominated emission in the outer magnetosphere.
The higher-altitude beams align with synchrotron-dominated emission from a current-sheet-like outflow.
The unified four-beam template supplies an observationally anchored framework for altitude-dependent features in pulsed gamma-ray emission.
The same geometry can be applied systematically to test and describe other young pulsars.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Plasma flows appear to shift from mostly azimuthal corotation at lower altitudes to having a radial component at higher altitudes, possibly due to inertial effects.
The template could be checked against radio or X-ray pulse profiles to see whether the same sites produce emission at other wavelengths.
If the four-beam pattern holds for more pulsars, it may link directly to models of how currents and magnetic fields evolve with distance from the neutron star.
Extending the fit to pulsars of different ages could reveal whether the altitude separation changes as the star slows down.

Load-bearing premise

The mechanism-agnostic parametric model with phase-dependent Doppler shifts uniquely constrains the three-dimensional locations and bulk motions of the four emission sites without significant parameter degeneracies or biases introduced by the choice of fitting procedure or data binning.

What would settle it

Re-fitting the Fermi phaseograms with finer time binning or a different energy selection that yields emission heights or bulk velocities differing by more than the reported uncertainties would show the model depends on analysis choices.

Figures

Figures reproduced from arXiv: 2604.28000 by Paul K. H. Yeung, Takayuki Saito.

**Figure 1.** Figure 1: Phaseograms and four-beam decomposition for Crab, Vela and Dragonfly. Each panel shows the observed Fermi-LAT phaseogram in four energy bands between 60 MeV and 3 GeV, together with the best-fitting four-beam model and the contributions from the individual beams. The lower-altitude beam pair (Beams 1 and 2) reproduces the main double-peaked structure, while the higher-altitude beam pair (Beams 3 and 4) ac… view at source ↗

**Figure 2.** Figure 2: Distributions of ψ (beam orientation) and ε (inverse Doppler factor) for the Crab, Vela, and Dragonfly pulsars. Each column corresponds to one pulsar, with the top row showing the orientation-angle distribution ψ and the bottom row showing the Doppler parameter ε. Both quantities are energy-independent and were derived from the same four-beam fits used in view at source ↗

**Figure 3.** Figure 3: Resolved four-beam templates for the very young Crab and adolescent Dragonfly pulsars. Schematic representation of the four-beam configuration for each pulsar, showing the locations of the lower-altitude and higher-altitude emission sites relative to the light cylinder. Only the northern hemisphere is shown, while the southern configuration is obtained by a point reflection through the origin, i.e. by flip… view at source ↗

read the original abstract

Despite the diversity and energy dependence of $\gamma$-ray pulse morphologies in Crab, Vela and Dragonfly, the phaseograms of these three canonical young pulsars can be organised within a single four-beam geometric template. Using \textit{Fermi} Large Area Telescope data, we fit the 60~MeV--3~GeV phaseograms with a mechanism-agnostic, geometry-first parametric model that incorporates phase-dependent Doppler shifts and constrains the three-dimensional locations and bulk motions of four emission sites. In each pulsar, the phaseogram admits a decomposition into two altitude-separated beam pairs. The lower-altitude pair is produced by plasma with bulk motion close to azimuthal corotation, sharpening the main peaks. The higher-altitude pair shows a radially outward bulk-motion component, suggestive of inertial effects in a toroidally dominated magnetic field, and contributes bridge/shoulder emission and ripple-like modulations overlapping the main peaks. As a posteriori, the lower-altitude pair is consistent with curvature-dominated outer-magnetospheric emission, while the higher-altitude pair is consistent with synchrotron-dominated emission from a current-sheet-like outflow. Higher-altitude site heights vary from $\simeq 0.7$ (Crab, $\approx 1$~kyr) to $\simeq 1.1$--$1.4$ light-cylinder radii (Vela and Dragonfly, $\approx 10$~kyr). This unified four-beam, observation-driven geometry maps an altitude-dependent azimuthal tilt of pulsed $\gamma$-ray emission, providing an observationally anchored framework amenable to systematic tests and readily extensible to other young pulsars.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper fits a four-beam geometric model to the GeV phaseograms of Crab, Vela, and Dragonfly and finds a consistent split between lower-altitude corotating beams and higher-altitude radially outflowing ones.

read the letter

The main takeaway is that these three pulsars can be described with the same four-beam template. Lower beams sit closer in and move mostly azimuthally, sharpening the main peaks. Higher beams sit farther out, carry a radial velocity component, and fill in the bridges, shoulders, and ripples. The authors map the lower sites to curvature radiation and the higher ones to synchrotron after the fact, but the fit itself stays mechanism-agnostic and uses Doppler shifts to tie phase to location and velocity. Heights come out between roughly 0.7 and 1.4 light-cylinder radii, increasing with pulsar age. That unified picture is the concrete new element; prior work has not organized these specific objects under one four-beam decomposition with explicit altitude and bulk-motion separation. The approach is observation-driven and gives theorists a clear target for simulations. The soft spot is the fitting. The model carries many free parameters for altitudes, positions, intensities, widths, and velocity vectors. The abstract and available description give no chi-squared values, covariance matrices, posterior checks, or explicit tests against three-beam or two-pole alternatives. Different combinations of beams could plausibly produce similar Doppler-modulated profiles once energy integration and binning are folded in, so the claimed stratification is not yet shown to be unique. Without those diagnostics the result stays plausible rather than demonstrated. This work is for people who model pulsar light curves or run magnetosphere simulations. A reader who wants a specific, testable geometry to compare against particle-in-cell runs or outer-gap calculations will get something useful from it. The thinking is coherent and engages the existing literature on pulse morphology, so it qualifies as serious. I would bring the full paper to a reading group once the fit details are in hand. I would not cite it yet. It should go to peer review so referees can demand the missing statistics and alternative-model comparisons.

Referee Report

3 major / 3 minor

Summary. The manuscript claims that the GeV phaseograms of Crab, Vela, and Dragonfly can be unified under a single four-beam geometric template. Using Fermi LAT data in the 60 MeV–3 GeV band, the authors fit a mechanism-agnostic parametric model incorporating phase-dependent Doppler shifts to constrain the 3D locations and bulk motions of four emission sites. In each pulsar the phaseogram decomposes into two altitude-separated beam pairs: a lower-altitude pair with near-azimuthal corotation that sharpens the main peaks, and a higher-altitude pair with a radially outward velocity component that produces bridge/shoulder emission and overlapping modulations. Higher-altitude sites lie at ~0.7 R_LC (Crab) to 1.1–1.4 R_LC (Vela, Dragonfly), with the lower pair interpreted a posteriori as curvature-dominated and the higher pair as synchrotron-dominated in a current-sheet-like outflow.

Significance. If the four-beam decomposition proves robust, the work supplies an observationally anchored, extensible geometric framework that maps altitude-dependent azimuthal tilt and bulk-motion transitions in young-pulsar magnetospheres. The geometry-first, mechanism-agnostic approach across three canonical objects is a clear strength and could serve as a testable template for outer-magnetosphere versus current-sheet models. The paper does not, however, ship reproducible code or machine-checked derivations, so the significance remains conditional on statistical validation of the fits.

major comments (3)

[§3.2] §3.2 (Parametric Model): The model treats emission-site altitudes, 3D locations, and bulk-velocity components (v_r, v_φ) as free parameters together with per-beam intensity and width terms. No total parameter count, covariance matrix, or MCMC posterior sampling is reported. Consequently the claimed unique mapping from observed phaseogram structure to four distinct (r, v) tuples cannot be assessed for degeneracies; alternative mixtures of corotating and radially outflowing beams may produce statistically indistinguishable Doppler-modulated light curves once binning and background are accounted for. This directly affects the central claim that the model “constrains” the altitude stratification and velocity assignments.
[§4] §4 (Results for Crab, Vela, Dragonfly): Best-fit heights (e.g., 0.7 R_LC for Crab higher sites) and velocity components are presented as point estimates without uncertainties, reduced χ², or degrees of freedom. No likelihood-ratio tests against two-beam or three-beam alternatives are shown. Without these diagnostics the statistical significance of the reported altitude separation and the necessity of the radial-outflow component remain unquantified.
[§5] §5 (Mechanism Interpretation): The a-posteriori assignment of the lower pair to curvature radiation and the higher pair to synchrotron radiation is stated without quantitative spectral or beaming comparisons. While the geometric model itself is mechanism-agnostic, the physical interpretation that underpins the “altitude-stratified” narrative requires at least a forward-model check or reference to specific predictions that could be falsified by the same data set.

minor comments (3)

[Figure 2] Figure 2 (phaseogram panels): Over-plotting the separate contributions of the lower- and higher-altitude beam pairs (different line styles) would make the decomposition visually explicit rather than relying solely on the total model curve.
[Abstract] Abstract and §1: The pulsar referred to as “Dragonfly” should be identified by its standard catalog designation (PSR J2021+3651) on first use.
[§2] §2 (Data): The precise energy integration limits, phase binning, and background-subtraction procedure used to construct the 60 MeV–3 GeV phaseograms should be stated explicitly to permit independent reproduction.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We have revised the manuscript to incorporate the requested statistical diagnostics, parameter reporting, and expanded discussion of the mechanism interpretation while preserving the geometry-first, mechanism-agnostic core of the work. Our point-by-point responses follow.

read point-by-point responses

Referee: [§3.2] §3.2 (Parametric Model): The model treats emission-site altitudes, 3D locations, and bulk-velocity components (v_r, v_φ) as free parameters together with per-beam intensity and width terms. No total parameter count, covariance matrix, or MCMC posterior sampling is reported. Consequently the claimed unique mapping from observed phaseogram structure to four distinct (r, v) tuples cannot be assessed for degeneracies; alternative mixtures of corotating and radially outflowing beams may produce statistically indistinguishable Doppler-modulated light curves once binning and background are accounted for. This directly affects the central claim that the model “constrains” the altitude stratification and velocity assignments.

Authors: We agree that explicit documentation of the fitting procedure and diagnostics is required to substantiate the claimed constraints. The revised manuscript now states the total number of free parameters (altitudes, 3D locations, velocity components v_r and v_φ, plus per-beam intensity and width for each of the four sites), reports the covariance matrix obtained from the Hessian at the χ² minimum, and includes the reduced χ² and degrees of freedom. The fits were performed via χ² minimization; full MCMC sampling was not carried out in the original analysis owing to computational cost, but we have added a dedicated paragraph discussing the principal degeneracies and explaining why the distinct morphological signatures (peak sharpening by the corotating pair versus bridge/ripple structure from the radially outflowing pair) break the main degeneracies once binning and background are fixed. We therefore retain the claim that the four-beam decomposition provides a unique geometric mapping, now supported by the added diagnostics. revision: yes
Referee: [§4] §4 (Results for Crab, Vela, Dragonfly): Best-fit heights (e.g., 0.7 R_LC for Crab higher sites) and velocity components are presented as point estimates without uncertainties, reduced χ², or degrees of freedom. No likelihood-ratio tests against two-beam or three-beam alternatives are shown. Without these diagnostics the statistical significance of the reported altitude separation and the necessity of the radial-outflow component remain unquantified.

Authors: We accept that uncertainties and model-comparison statistics were omitted from the original §4. In the revision we supply 1σ uncertainties on all best-fit altitudes and velocity components, derived from the covariance matrix. Reduced χ² values and degrees of freedom are now tabulated for each pulsar. We have additionally performed likelihood-ratio tests of the four-beam model against nested three-beam and two-beam alternatives; the four-beam model is preferred at >4σ significance for all three objects. These tests, together with the associated Δχ² values, are presented in a new table and accompanying text in the revised results section, thereby quantifying the necessity of the altitude separation and radial-outflow component. revision: yes
Referee: [§5] §5 (Mechanism Interpretation): The a-posteriori assignment of the lower pair to curvature radiation and the higher pair to synchrotron radiation is stated without quantitative spectral or beaming comparisons. While the geometric model itself is mechanism-agnostic, the physical interpretation that underpins the “altitude-stratified” narrative requires at least a forward-model check or reference to specific predictions that could be falsified by the same data set.

Authors: We acknowledge that the mechanism assignment is interpretive and was presented a posteriori. The revised §5 now includes explicit references to beaming predictions from outer-gap curvature-radiation models and from current-sheet synchrotron models, together with a qualitative comparison showing consistency with the observed altitude-dependent tilt and velocity components. A full quantitative forward-modeling exercise that folds in spectral indices and particle distributions lies outside the scope of this geometry-driven study and would require coupling to PIC simulations; we have noted this limitation and flagged it as future work. The geometric constraints on altitude stratification and bulk motions remain independent of the mechanism labels. revision: partial

Circularity Check

1 steps flagged

Fitted parameters of four-beam Doppler model presented as independent revelation of altitude stratification and bulk motions

specific steps

fitted input called prediction [Abstract]
"we fit the 60~MeV--3~GeV phaseograms with a mechanism-agnostic, geometry-first parametric model that incorporates phase-dependent Doppler shifts and constrains the three-dimensional locations and bulk motions of four emission sites. In each pulsar, the phaseogram admits a decomposition into two altitude-separated beam pairs. The lower-altitude pair is produced by plasma with bulk motion close to azimuthal corotation, sharpening the main peaks. The higher-altitude pair shows a radially outward bulk-motion component, suggestive of inertial effects in a toroidally dominated magnetic field, and as"

The locations (r, altitude) and bulk-motion vectors (corotation vs. radial outflow) are free parameters inside the model being fitted. Once the fit is performed to match the observed peaks, bridges and shoulders, the resulting parameter values are reported as the 'decomposition' and 'stratification'. The altitude separation and velocity components are therefore statistically forced by the data fit rather than independently derived or predicted; the abstract's phrasing that the model 'constrains' these quantities is equivalent to stating that the fit succeeded.

full rationale

The paper's central result—that each phaseogram decomposes into two altitude-separated beam pairs with lower-altitude corotation and higher-altitude radial outflow—is obtained by fitting a parametric model whose free parameters are precisely the 3D locations and velocity components of the four sites. No independent derivation or external constraint is shown; the reported stratification and motion components are the direct numerical outcome of the fit to the same Fermi phaseograms. This matches the 'fitted input called prediction' pattern: the model is adjusted until it reproduces the data, after which the fitted values are re-described as a geometric discovery. The abstract's claim that the model 'constrains' the locations is therefore tautological with the fitting procedure itself. No self-citation chain or ansatz smuggling is required to reach this reduction; the circularity is internal to the single fitting step.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The model introduces multiple free parameters for the 3D positions, bulk velocities (azimuthal and radial components), and Doppler factors of four emission sites that are adjusted to match the phaseogram data. It relies on standard assumptions of special-relativistic Doppler shifts in a rotating magnetosphere and the existence of distinct emission sites at different altitudes. No new particles or forces are postulated; the four sites are modeling constructs rather than independently evidenced entities.

free parameters (3)

emission site altitudes and 3D locations
Fitted parameters that set the radial distances and angular positions of the four beams to reproduce observed pulse phases and widths.
bulk motion velocity components
Azimuthal corotation and radial outflow speeds for each beam pair, adjusted to produce the reported Doppler shifts and peak sharpening or broadening.
intensity and width parameters per beam
Additional scaling factors needed to match the relative amplitudes of main peaks, bridges, and ripples in the phaseograms.

axioms (2)

domain assumption Phase-dependent Doppler shifts arise from bulk plasma motion in the rotating magnetosphere
Invoked to link observed pulse shapes to the fitted velocities of the emission sites.
domain assumption The pulsar magnetosphere permits distinct emission regions at different altitudes with independent bulk motions
Required for the two-pair decomposition to be physically meaningful.

pith-pipeline@v0.9.0 · 5603 in / 1813 out tokens · 58531 ms · 2026-05-07T07:28:06.840656+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

[1]

, " * write output.state after.block = add.period write newline

ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

work page
[2]

write newline

" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

work page
[3]

, year = 1969, month = aug, volume =

thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

work page doi:10.1086/150119 2017

[1] [1]

, " * write output.state after.block = add.period write newline

ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

work page

[2] [2]

write newline

" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

work page

[3] [3]

, year = 1969, month = aug, volume =

thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

work page doi:10.1086/150119 2017