Pith. sign in

REVIEW 4 major objections 5 minor 15 references

More than 300 gamma-ray pulsars show a wide variety of pulse shapes; this paper argues the diversity is explained by asymmetric pair-plasma loading in the striped-wind current sheet, with a few physically motivated parameters.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-03 13:12 UTC pith:OIFI5LGR

load-bearing objection A large and useful atlas of split-monopole light curves with asymmetric emissivity, but the polar-cap pair-loading interpretation rests on an unvalidated sheet locus and an ad hoc emissivity mapping. the 4 major comments →

arxiv 2607.29126 v1 pith:OIFI5LGR submitted 2026-07-31 astro-ph.HE

Polar cap plasma loading and the morphology of pulsar γ-ray light curves

J\'er\^ome P\'etri This is my paper
classification astro-ph.HE
keywords gamma-ray pulsarslight-curve morphologystriped windcurrent sheetpair plasma loadingspherical harmonicssplit monopolepulse profile
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

More than 300 gamma-ray pulsars are known, and their pulse profiles come in a wide range of shapes that a single emission geometry has been hard to reproduce. This paper argues that the diversity can be explained if the high-energy emission comes from the equatorial current sheet of the striped pulsar wind, provided the emissivity along that sheet is not uniform but follows the pair-plasma production pattern of the two polar caps. By breaking north-south and azimuthal symmetry with a few parameters — a pair-multiplicity ratio, two spherical-harmonic mode numbers, and phase shifts — the model reproduces symmetric, asymmetric, multi-peaked, and bridge-emission profiles from an atlas of 16.5 million light curves, and fits 130 observed pulsars. The dominant pattern is a dipolar (m=1) loading mode, with quadrupolar modes accounting for more than 90% of the sample; millisecond pulsars prefer somewhat higher modes, hinting at non-dipolar surface fields. If correct, the pulse profile is not arbitrary: it encodes the plasma-loading structure of the two polar caps as transported into the current sheet, making the current sheet a unified framework for high-energy pulsar emission.

Core claim

The paper's central claim is that a spatially dependent emissivity in the split-monopole current sheet — not a collection of separate emission zones — is enough to account for most observed gamma-ray pulse shapes. The current sheet is located where the radial field reverses, at the locus cosθ cosα + sinθ sinα cosψ = 0 (Eq. 1), and each polar cap contributes an emissivity f_K = κ_K |Y^R_{mK,mK}(θ,φ,ϕ_K)|^2 r^{-q}, meaning the pair-multiplicity pattern above each cap is squared and carried into the sheet with a radial falloff. Fitting this model to 130 observed pulsars, the paper finds the symmetric mode (m_N = m_S = 0) almost never works (2% of top-3 fits), the simplest asymmetric dipolar mod

What carries the argument

The machinery is a two-part construction. First, the emitting surface: the split-monopole current sheet, defined analytically by cosθ cosα + sinθ sinα cosψ = 0 with ψ = φ − Ω(t − r/V), the locus of radial magnetic-field reversal, assumed to coincide with the current sheet of a real oblique force-free dipole. Second, the emissivity prescription: f_K(r,θ,φ) = κ_K |Y^R_{mK,mK}(θ,φ,ϕ_K)|^2 r^{-q}, where κ_K sets the pair multiplicity in polar cap K, m_K is the azimuthal spherical-harmonic mode in that cap, ϕ_K is an azimuthal rotation, and q controls the radial falloff. For m=0 with equal amplitudes this reduces to the old uniform-emissivity striped-wind model; for m≥1 the squared spherical harm

Load-bearing premise

The load-bearing premise is that gamma rays are emitted from the split-monopole current-sheet surface given by cosθ cosα + sinθ sinα cosψ = 0, and that the polar-cap pair pattern is transported onto it unchanged; if the true emitting sheet of a dipolar magnetosphere lies elsewhere or the pattern is distorted in transit, every atlas light curve and every fitted obliquity, viewing angle, and loading mode is shifted.

What would settle it

Compute the null surface where the radial magnetic field reverses in a force-free oblique-dipole magnetosphere and compare it with the split-monopole locus cosθ cosα + sinθ sinα cosψ = 0 across the full obliquity range; if the two surfaces differ by more than the current-sheet thickness, the atlas light curves and fitted (α, ζ, κ, m, ϕ) parameters are systematically displaced. A secondary check is to measure whether the pair density in the current sheet of a kinetic simulation actually tracks the polar-cap loading pattern, or is scrambled by reconnection.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If the current-sheet scenario with asymmetric loading is correct, a pulsar's gamma-ray pulse profile directly encodes the pair-loading pattern of its two polar caps as transported into the wind.
  • The symmetric uniform current-sheet model is essentially ruled out as a general explanation: it appears in the top-3 fits for only about 2% of the 130-pulsar sample.
  • More than 90% of the sample is fitted with dipolar (m=1) or quadrupolar (m=2) loading patterns, implying large-scale asymmetries dominate the emitting sheet over small-scale structure.
  • Millisecond pulsars systematically require higher-order loading modes than young pulsars, a statistical hint of non-dipolar magnetic structures near their surfaces — though the paper stresses this is not yet direct evidence.
  • Because the model uses a single emission surface, the magnetic obliquity α and viewing angle ζ can be extracted from peak separation and shape, with the loading parameters controlling asymmetry and substructure.

Where Pith is reading between the lines

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

  • Editorial inference: If the loading pattern truly travels from the polar caps to the current sheet, pulse-profile fitting becomes an indirect probe of the pair-creation cascade near the neutron-star surface, complementing radio drifting-subpulse observations.
  • Editorial inference: The model's preference for phase alignments near ϕ_N=0 and 0.5 is a testable prediction; force-free or particle-in-cell simulations of polar-cap current patterns could confirm or refute that preferred orientation, since the paper itself leaves the causal link unresolved.
  • Editorial inference: The model's flexibility implies a strong degeneracy between geometry (α, ζ) and loading parameters (κ, m, ϕ); population-level conclusions about spin-axis geometry should marginalize over loading parameters before interpreting fitted angles.
  • Editorial inference: A direct numerical comparison of the split-monopole null surface with the true current-sheet null surface of a force-free oblique dipole — especially near the light cylinder and at high obliquity — would quantify the largest systematic error in this atlas.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The manuscript constructs an analytic model of pulsed gamma-ray emission in which the emitting region is the split-monopole striped-wind current sheet (Eq. 1) and the local emissivity is modulated by squared spherical-harmonic patterns associated with the north and south polar caps (Eq. 4). The author computes a large atlas of light curves, fits the model to 130 Fermi 3PC pulsars, and reports that low-order asymmetric plasma-loading modes, especially (m_N, m_S) = (1, 1), reproduce a large fraction of observed profiles, with millisecond pulsars preferring somewhat higher-order patterns than young pulsars. The central claim is that pulse morphology encodes the polar-cap pair-loading pattern as transported into the current sheet.

Significance. If the two key premises of the model were validated—namely, that the gamma-ray emitting sheet coincides with the split-monopole locus and that the polar-cap pattern is advected unchanged into the sheet—this would be a useful, economical unified framework for Fermi pulsar light curves and would connect macroscopic pulse morphology to polar-cap pair physics. The paper is computationally thorough and transparent about several limitations, and the atlas itself may be a useful reference. However, the physical interpretation goes beyond what is demonstrated: the fit quality is assessed with an ad hoc score, no parameter uncertainties are reported, and the two load-bearing geometric/emissivity assumptions are asserted rather than tested. The conclusions are hedged in places, but the abstract and parts of §4.5 make stronger causal claims than the model currently supports.

major comments (4)
  1. [§2.1, Eq. (1)] The emitting surface is identified with the split-monopole locus cosθ cosα + sinθ sinα cosψ = 0. The text states that 'several previous works showed' the force-free dipole sheet lies nearly at this location, but no citation or quantitative comparison is given. Every fitted α, ζ, and all mode counts depend on this locus. If the true current sheet in a dipolar magnetosphere differs, e.g., in the near zone or for high obliquity, all fitted parameters and the population-level distributions in Fig. 19 shift. This premise needs to be supported by a citation and a quantitative error estimate, or the paper must explicitly restrict its physical interpretation to the split-monopole toy model rather than claiming that observed morphology encodes the pair-loading pattern.
  2. [§2.2, Eq. (4)] The central emissivity prescription f_K = κ_K |Y^R_{mK,mK}|² r^{-q} assumes that the polar-cap spherical-harmonic pattern propagates unchanged into the current sheet, with φ replaced by ψ. The paper itself calls this 'phenomenological' and later concedes it is 'not derived from a self-consistent treatment of pair creation and injection' (§4.6). This mapping is load-bearing for the conclusion that pulse morphology reflects polar-cap pair loading. Since only this one emissivity family is used, the successful fits do not test the transport assumption. I request either a physical derivation/justification of Eq. (4) or a robustness check with an alternative ansatz (e.g., localized spots, Gaussian columns, or different harmonic orders) to show that the inferred mode distribution is not an artifact of the specific prescription.
  3. [§4.2–§4.3, Table 2, Eq. (6)] The population-level claims (e.g., the (1,1) mode in the top-3 for 60% of pulsars, and the young/MSP separation in Fig. 19) are based on χ² minimization with the modified score S_mod of Eq. (6), but no uncertainties on the fitted parameters are reported and no model-selection or cross-validation is performed. Figure 17 shows that many fits, especially for bright pulsars, have very large χ²/d.o.f. values; the 'top-3' criterion selects among poorly fitting models. The statement in §4.5 that 'more than 90%' of pulsars are reproduced by dipolar/quadrupolar modes needs an explicit goodness-of-fit threshold. Without parameter errors or a validation scheme, the mode-count distributions may be dominated by noise and degeneracy rather than by physically meaningful structure.
  4. [§3.1, §4.3] The atlas is described as containing 16,549,260 light curves, but the fitting section uses only a coarse grid of κ_S/κ_N ∈ {1,2,5,10}, m up to 4, and fixed phase grids. The paper does not discuss the uniqueness or degeneracy of the best-fit parameters, nor the sensitivity of the mode ranking to the chosen grid resolution. Since multiple nearly degenerate solutions are acknowledged (Fig. 18 and the top-3 discussion), the reported histograms of m_N, m_S, φ_N, and Δφ should be accompanied by an assessment of how representative the selected best fits are, for example by showing the spread of parameters among all statistically acceptable fits.
minor comments (5)
  1. [§2.2, Eq. (2)] The notation arctan(A,B) should be defined explicitly as the two-argument arctangent, and the range of the angle used in cos(m arctan(A,B)) should be stated to avoid ambiguity.
  2. [§4.2] There are several typographical errors: 'reproduces by' should be 'reproduced by', 'wit the mode' should be 'with the mode', and 'depict' should be 'depicted'.
  3. [Fig. 10] 'translucence blue square' and 'translucence red cross' should be 'translucent'; the caption could state briefly how the α and ζ values were obtained for the plotted pulsars.
  4. [§4.4] The sentence 'firm conclusions are difficult to drawn' should read 'difficult to draw'.
  5. [References] The claim in §2.1 that previous works located the force-free dipole sheet near the split-monopole locus is made without citation; please add the relevant references or remove the claim.

Circularity Check

0 steps flagged

No significant circularity: the atlas is a forward model and the fitted loading parameters are explicitly labeled phenomenological, not independent predictions.

full rationale

The claimed derivation chain is: adopt the split-monopole current-sheet locus (Eq. 1), prescribe an emissivity (Eq. 4), compute light curves over a grid, and fit those light curves to Fermi pulsar profiles. The paper does not claim that the observed light curves are derived from first principles; it explicitly describes the emissivity as "phenomenological but physically motivated" (§2.2) and later concedes that it "is not derived from a self-consistent treatment of pair creation and injection" (§4.6). The statement that asymmetric emissivity produces asymmetric profiles is a mathematical consequence of the forward model, not a circular inference from the data: the asymmetry parameters were introduced before the light-curve computation, and the atlas is generated independently of the 3PC data before fitting. The fitted mode/phase/asymmetry distributions are parameter estimates under an assumed model, not independent predictions; they should be weighed for model validity and degeneracy, but that is a correctness/robustness concern, not circularity in the derivation. The main load-bearing premise — that the force-free dipolar current sheet lies near the split-monopole locus — is asserted in §2.1 with "Several previous works showed" but without a citation or quantitative comparison, and a dipolar treatment is deferred; this is a missing-support/correctness risk, not a circular equivalence. Self-citations to Pétri (2011, 2024, 2026) and Pétri & Mitra (2021) are used for the base geometry and method, not to forbid alternatives or to import a uniqueness theorem. The paper's conclusions are explicitly hedged "within the adopted split-monopole geometry and emissivity prescription," so it does not overclaim an independent verification of pair loading. No equation or fitted parameter is, on inspection, identical to the claimed result by construction in a way that qualifies as circularity under the rubric.

Axiom & Free-Parameter Ledger

7 free parameters · 5 axioms · 0 invented entities

The central claim rests on a split-monopole geometrical proxy, an assumed current-sheet emission site, and an ad hoc spherical-harmonic mapping from polar-cap pair supply to emissivity. The fitted parameters (α, ζ, κ, φ, m, q) are numerous enough that the atlas's ability to match observed shapes is largely expected from model flexibility rather than from independent physical constraint.

free parameters (7)
  • magnetic obliquity α = best-fit per pulsar from a 35-value grid (5°–175°); not tabulated
    Scanned and inferred per pulsar; controls peak separation and overall morphology.
  • observer line-of-sight angle ζ = best-fit per pulsar from a 91-value grid (0°–180°); not tabulated
    Scanned and inferred per pulsar; degenerate with α in many fits.
  • north-south emissivity ratio κ_S/κ_N = discrete grid {1,2,5,10}; distribution in Fig. 16
    Controls relative peak amplitudes and asymmetry of the profile.
  • spherical-harmonic modes m_N, m_S = 0–4 with m_S ≤ m_N; top-3 ranking in Table 2
    Discrete model complexity per polar cap; central to the population-level claims.
  • north-pole emissivity phase φ_N = 11 grid values; distribution in Fig. 15
    Azimuthal orientation of the northern emissivity pattern.
  • south-pole emissivity phase φ_S (or Δφ) = 11×11 grid; distribution in Fig. 15
    Relative azimuthal offset of the southern pattern; nearly uniform after degeneracy removal.
  • radial emissivity index q = not reported / not fitted
    Power-law decline r^{-q} in Eq. (4); introduced by hand without constraint or fitted value.
axioms (5)
  • domain assumption The force-free split-monopole solution is a sufficient proxy for the global pulsar magnetosphere.
    All light curves are computed from this geometry; dipolar solutions are explicitly deferred (§2, intro and §2.1).
  • domain assumption Gamma-ray emission is produced only in a thin current sheet where B_r reverses sign, and the force-free dipole sheet lies nearly at the split-monopole surface.
    Sect. 2.1; the equivalence to the dipole is asserted with 'several previous works' but no citation or derivation is given.
  • ad hoc to paper Current-sheet emissivity is proportional to |shifted spherical harmonic Y^R_{m,m}|² times r^{-q}, with φ replaced by the phase ψ, as a proxy for pair plasma loading.
    Eqs. (2)–(4); no derivation from pair cascades or particle transport; the paper calls it phenomenological.
  • ad hoc to paper Restricting to ℓ=m and m≤4 preserves the physically relevant angular structure of the loading pattern.
    §2.2 and §3.4; chosen for tractability and coherence with expected pair-cascade scales, not derived.
  • domain assumption A line-of-sight integration of the emissivity with relativistic beaming, ignoring spectral and polarization details, captures pulse morphology.
    §2.1; the radiation mechanism is unspecified, which is acceptable for morphology but limits physical completeness.

pith-pipeline@v1.3.0-daily-deepseek · 14457 in / 22939 out tokens · 237035 ms · 2026-08-03T13:12:29.414155+00:00 · methodology

0 comments
read the original abstract

The discovery of more than 300 $\gamma$-ray pulsars by the Fermi Large Area Telescope (Fermi-LAT) has revealed a rich diversity of light-curve morphologies that remains challenging to reproduce within purely magnetospheric emission models, suggesting that particle acceleration and radiation may occur, at least in part, in the striped-wind current sheet. We compute a comprehensive atlas of pulsar $\gamma$-ray light curves based on the split-monopole current sheet geometry, with the goal of quantifying the role of pair plasma loading in shaping the observed pulse morphology. The current-sheet surface is described analytically and assumed to be the dominant site of high-energy photon emission. Spatially dependent emissivity prescriptions are introduced through a minimal set of parameters that explicitly break the north-south and azimuthal symmetries of the emitting plasma. With only a few physically motivated parameters, the model is able to reproduce a broad range of the observed $\gamma$-ray light-curve morphologies, including asymmetric, multi-peaked, and highly structured profiles. The results show that asymmetric emissivity, possibly related to pair-plasma loading, provides a natural and economical mechanism capable of reproducing a substantial fraction of the observed diversity of $\gamma$-ray pulse profiles. The current-sheet scenario therefore appears to provide a promising and physically motivated framework for high-energy pulsar emission.

Figures

Figures reproduced from arXiv: 2607.29126 by J\'er\^ome P\'etri.

Figure 1
Figure 1. Figure 1: An example of emissivity for an obliquity α = 60◦ and the modes mN = 1 and mS = 2 in the (θ, φ) plane. The blue colour means zero density whereas red means maximal density normalised to one. The black undulating line depicts the magnetic equator located at θ ′ = π/2. to be confused with the rotational equator θ = π/2) [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Same as [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Representative γ-ray light curves for each morphological class with the parameters (mN, mS) = (0, 0) and κS/κN = 1 i.e. a fully sym￾metric profile. The separation is forced into six classes. A map of each morphological class in the (α, ζ) plane for (mN, mS) = (0, 0) is shown in Fig.4. Two sharp peaks well sepa￾rated by almost half a period requires α ∼ 90◦ or ζ ∼ 90◦ . Two overlapping peaks requires α + ζ … view at source ↗
Figure 5
Figure 5. Figure 5: Representative γ-ray light curves sample for (mN, mS) = (1, 1) and κS/κN = 1. The asymmetric nature of these light curves is readily visible. exactly the same number of phase points as the theoretical pro￾files. This allows us to use straightforwardly the fast Fourier transform for the cross-correlation as explained in Lorange et al. (2026). Our sample shows the wealth of profile variety: one symmetric or … view at source ↗
Figure 6
Figure 6. Figure 6: A representative subset of the variety of γ-ray pulsar light curves extracted from our sample of 130 pulsars. The best fits for the symmetric mode mN = mS = 0 are shown in Fig.7. Although some pulsars can be explained by this sim￾ple model, the vast majority requires at least asymmetric shapes reproduces by the mode mN = mS = 1 as shown in Fig.8. Fi￾nally, Fig.9 shows the best fit obtained by imposing mN ≤… view at source ↗
Figure 8
Figure 8. Figure 8: Same as Fig.7 but wit the mode mN = mS = 1 [PITH_FULL_IMAGE:figures/full_fig_p007_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Same as Fig.7 but wit the mode satisfying mN ≤ 4 and mS ≤ 4. with bridge emission, impossible to reproduce with any values of mN and mS modes. Even more difficult is J1536-4948 with for peaks and significant bridge emission. The best we could do is to use a (mN, mS) = (3, 1) mode, getting four peaks but not all aligned with observations. Next PSR J0205+6449, like the Crab, shows two well separated narrow p… view at source ↗
Figure 12
Figure 12. Figure 12: Same as [PITH_FULL_IMAGE:figures/full_fig_p008_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Distribution of the emissivity modulation complexity required to reproduce the observed γ-ray light curves of the full sample irrespec￾tive of the young or MSP nature, summarised by the north and south polar cap modes mN and mS, imposing also mS ≤ mN due to the north￾south symmetry. Therefore, to better estimate the number of modes required for explaining the light curves, we extract the top-3 models for … view at source ↗
Figure 14
Figure 14. Figure 14: Same as [PITH_FULL_IMAGE:figures/full_fig_p009_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Distribution of the north pole emissivity phase alignment ϕN, on the left histogram, and distribution of the phase shift ∆ϕ between the south and north pole, on the right histogram, for the young pulsars in blue and for the millisecond pulsars in red. Phase alignment values are favoured around ϕN = 0 and ϕN = 0.5 for both classes whereas the ∆ϕ distribution is more compatible with a uniform distribution … view at source ↗
Figure 16
Figure 16. Figure 16: Distribution of emissivity asymmetry κS/κN between the south and north pole for the two classes of pulsars: young in blue and MSP in red. The gaps are artifact due to our chosen discrete values being "only" κS/κN = {1, 2, 5, 10}. curves. Low-order modes with m ≤ 1 reproduce a large sample of the observed light curves, suggesting possibly that the pair loading modulation may be dominated by large-scale asy… view at source ↗
Figure 17
Figure 17. Figure 17: quantifies the decrease in χ 2 /d.o.f. value when shift￾ing from the symmetric mode (mN, mS) = (0, 0), in blue, to the simplest asymmetric mode (mN, mS) = (1, 1), in red, and finally to the best fit solution (mN ≤ 4, mS ≤ 4), in black. The brightest pulsars are generally on the right sight of this plot, and are badly fitted because of the small uncertainties in the intensity. This ex￾plains their high χ 2… view at source ↗
Figure 18
Figure 18. Figure 18: Diagram showing the χ 2 ratios for the top-3 best solutions for each pulsar, on the x-axis χ 2 top2 /χ2 top1 and on the y-axis χ 2 top3 /χ2 top2 . Blue dots correspond to young pulsars and red triangles to MSP. The χ 2 val￾ues are all comparable. The question therefore arises on which model to pick out according to supplementary constraints. One of the main crite￾ria is certainly physical simplicity. Cons… view at source ↗
Figure 19
Figure 19. Figure 19: Distribution of the minimal modes (mN, mS) required to rea￾sonably fit separately the young pulsars, top panel, and the millisecond pulsars, bottom panel. places the dominant emission region outside the light cylinder, consistent with the pulsar striped wind high energy emission model, and reproduces a wide range of observed light curve mor￾phologies with only a few free parameters, making it both eco￾nom… view at source ↗
Figure 20
Figure 20. Figure 20: Distribution of the minimal modes (mN, mS) required to rea￾sonably fit separately the young pulsars, top panel, and the millisecond pulsars, bottom panel. radiation processes, and should therefore be regarded as a geo￾metrical, parametric framework rather than a completely phys￾ical model. In addition, its applicability is restricted by the as￾sumed magnetic field structure, namely a split monopole config… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

15 extracted references · 1 linked inside Pith

  1. [1]

    & Spitkovsky, A

    Bai, X.-N. & Spitkovsky, A. 2010, ApJ, 715, 1282

    2010

  2. [2]

    A., & Spitkovsky, A

    Cerutti, B., Philippov, A. A., & Spitkovsky, A. 2016, MNRAS, 457, 2401

    2016

  3. [3]

    2024, ApJL, 974, L32

    Chernoglazov, A., Philippov, A., & Timokhin, A. 2024, ApJL, 974, L32

    2024

  4. [4]

    Coroniti, F. V . 1990, ApJ, 349, 538

    1990

  5. [5]

    E., Lupsasca, A., & Philippov, A

    Gralla, S. E., Lupsasca, A., & Philippov, A. 2017, ApJ, 851, 137

    2017

  6. [6]

    2012, MNRAS, 420, 2793

    Kalapotharakos, C., Contopoulos, I., & Kazanas, D. 2012, MNRAS, 420, 2793

    2012

  7. [7]

    2026, A&A, 707, A395

    Lorange, A., Pétri, J., Sautron, M., & Vigon, V . 2026, A&A, 707, A395

    2026

  8. [8]

    Michel, F. C. 1994, ApJ, 431, 397 Pétri, J. 2011, MNRAS, 412, 1870 Pétri, J. 2024, A&A, 687, A169 Pétri, J. 2026, A&A, 706, A215 Pétri, J. & Mitra, D. 2021, A&A, 654, A106

    1994

  9. [9]

    Ruderman, M. A. & Sutherland, P. G. 1975, ApJ, 196, 51

    1975

  10. [10]

    A., Abdollahi, S., Ajello, M., et al

    Smith, D. A., Abdollahi, S., Ajello, M., et al. 2023, ApJ, 958, 191

    2023

  11. [11]

    2006, ApJ, 648, L51

    Spitkovsky, A. 2006, ApJ, 648, L51

    2006

  12. [12]

    Timokhin, A. N. & Arons, J. 2013, MNRAS, 429, 20

    2013

  13. [13]

    Uzdensky, D. A. & Spitkovsky, A. 2014, ApJ, 780, 3

    2014

  14. [14]

    2026, arXiv:2607.08596 [astro-ph.HE]

    Zhang, J., Zhu, H., Lin, G., et al. 2026, arXiv:2607.08596 [astro-ph.HE]

    Pith/arXiv arXiv 2026

  15. [15]

    W., Lin, L

    Zhao, J., Ng, C. W., Lin, L. C. C., et al. 2017, ApJ, 842, 53 Article number, page 12 of 15 J. Pétri: Polar cap plasma loading and the morphology of pulsarγ-ray light curves Appendix A: Light curve fitting results: full sample In this appendix, we show the full sample of fitted light curves, concerning the 130 pulsars we selected. Fig. A.1 shows the light...

    2017