{"id":"f6c62875-c505-4402-8df1-9fa5b3e8ccef","arxiv_id":"2607.29126","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Asymmetric polar-cap plasma loading in the striped-wind current sheet can reproduce most Fermi pulsar gamma-ray light-curve morphologies using low-order spherical-harmonic modes.","lead":"This paper builds a 16.5-million-light-curve atlas for pulsar gamma-ray pulses using a split-monopole current-sheet model with asymmetric plasma loading. It fits 130 Fermi pulsars and argues that low-order asymmetric emissivity can reproduce most observed profile shapes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Two unvalidated premises — the split-monopole sheet locus (Eq. 1) and the polar-cap-to-sheet emissivity mapping (Eq. 4) — carry the physical interpretation; without a check against a dipolar force-free sheet, the fitted morphology atlas does not constrain pair loading.","rationale":"The reader's weakest_assumption lists exactly these two premises: the split-monopole locus standing in for the force-free dipole current sheet and the ad hoc Eq. (4) mapping. My pass finds the same. The paper is self-aware about the phenomenological nature of the emissivity prescription (§4.6) and about the degeneracy of the fitted parameters (§4.3, conclusion), which is why the appropriate verdict is CONDITIONAL rather than ACCEPT or REJECT. However, the paper's strongest-claim phrasing — that morphology encodes pair-loading pattern — depends on both premises, and neither is independently checked. The absence of a citation for the 'several previous works' (§2.1) is a concrete, fixable gap. Because the atlas itself is large and the qualitative diversity result is plausible, this does not overturn the paper; it makes the central physical inference conditional on a test that is currently missing. A targeted force-free comparison is the cheapest decisive check, since it isolates the geometric assumption without requiring a full PIC treatment of pair cascades. If the sheet locus check passes, Eq. (4) still needs its own validation, but that is a second, separate step. I therefore keep the reader's CONDITIONAL verdict (UNCHANGED).","tokens_in":14832,"tokens_out":5756,"duration_ms":65098,"concrete_test":"Take a published force-free oblique-dipole solution (e.g., Kalapotharakos et al. 2012, or equivalent) for α = 30°, 60°, 90°. Locate the B^r = 0 surface in the wind zone for r ∈ [R_L, 10 R_L] and compare its local phase/tilt and radius with Eq. (1) using V = c. If the mean angular deviation exceeds ~5° or the radial offset exceeds ~0.1 R_L, recompute a representative slice of the atlas (say mN = mS = 1, κS/κN = 2, all α, ζ) with the simulated sheet surface and re-run the top-3 fits of §4.3. If the inferred mode/angle distributions change materially (≥10 pulsars switching preferred mode or ≥10° shift in α/ζ), the split-monopole locus is not a safe stand-in and the pair-loading interpretation is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that observed γ-ray pulse morphology encodes the polar-cap pair-loading pattern as transported into the striped-wind current sheet. Every atlas light curve and every fitted (α, ζ, mN, mS, κN/S, φ) is produced from two premises that the paper does not validate. First, §2.1 identifies the emitting surface with the split-monopole locus cosθ cosα + sinθ sinα cosψ = 0, asserting only that 'several previous works showed' the force-free dipolar current sheet lies nearly there; no citation or quantitative comparison is given. A detailed dipolar treatment is explicitly deferred (§2). If the true B^r=0 sheet differs in the near zone, in phase lag, or in tilt for high obliquity, then the entire atlas and all fits are displaced. Second, Eq. (4) assumes the polar-cap spherical-harmonic pattern κ_K |Y^R_{mK,mK}|² r^{-q} propagates unchanged into the sheet, with φ replaced by ψ. The paper itself labels this phenomenological (§2.2) and later concedes it is 'not derived from a self-consistent treatment of pair creation and injection' (§4.6). Because the model family is only this one, a failure of either step changes the inferred parameter distributions — e.g., the 'more than 90%' low-mode statement in §4.5 and the young/MSP mode separation in Fig. 19 — not merely the fine details. The conclusion is honestly hedged as 'within the adopted split-monopole geometry and emissivity prescription,' but the physical interpretation in the strongest claim goes beyond that hedge.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":15309,"tokens_out":4198,"duration_ms":45622,"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":[{"comment":"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.","section":"§2.1, Eq. (1)"},{"comment":"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.","section":"§2.2, Eq. (4)"},{"comment":"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.","section":"§4.2–§4.3, Table 2, Eq. (6)"},{"comment":"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.","section":"§3.1, §4.3"}],"minor_comments":[{"comment":"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.","section":"§2.2, Eq. (2)"},{"comment":"There are several typographical errors: 'reproduces by' should be 'reproduced by', 'wit the mode' should be 'with the mode', and 'depict' should be 'depicted'.","section":"§4.2"},{"comment":"'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.","section":"Fig. 10"},{"comment":"The sentence 'firm conclusions are difficult to drawn' should read 'difficult to draw'.","section":"§4.4"},{"comment":"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.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about its limitations and the atlas is a substantial computational product, but the physical conclusion currently rests on two unvalidated assumptions and on fits whose statistical significance is not established. I do not see this as an irreparable flaw; the claims can be re-scoped to the split-monopole parametric model, or the missing validation can be supplied. I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe thing to know: this is a big engineering extension of Pétri's split-monopole current-sheet light-curve program, and the atlas it produces is genuinely useful. But the physical interpretation—that pulse morphology encodes polar-cap pair loading—rests on two assumptions the paper never validates, and the statistical support for its population claims is softer than the abstract implies.\n\nWhat is new: previous split-monopole papers used symmetric, radial emissivity. Here the emissivity is modulated by squared spherical harmonics with separate north/south modes, amplitudes, phases, and a radial index. The paper computes 16.5M light curves and fits 130 Fermi pulsars. The figures show clearly that asymmetric emissivity produces asymmetric, multi-peaked, and bridge-emission profiles that symmetric models cannot make. That is a real step forward. The paper is also honest in places: it labels the emissivity prescription phenomenological, admits degeneracy between geometry and loading parameters, and concedes some pulsars (Vela, J1231-1411) cannot be reproduced.\n\nThe soft spots are concentrated in the interpretive layer. First, Eq. (1) identifies the emitting surface with the split-monopole current sheet, citing \"several previous works\" with no actual reference or quantitative check. If the force-free dipole sheet differs by tens of degrees for high obliquity, the entire atlas and all fitted angles are offset. Second, Eq. (4) maps the polar-cap spherical-harmonic pattern into the sheet by replacing φ with ψ. That is ad hoc, not derived from cascade or transport physics. The paper admits it, but the abstract's \"possibly related to pair-plasma loading\" leans on it.\n\nStatistically, the fits have no reported uncertainties. Figure 18 shows the top-3 models often have near-identical χ², so the mode distribution in Fig. 13 is fragile. The \"more than 90%\" claim in §4.5 overstates Table 2, which lists top-3 occurrences, not unique pulsars. No code or atlas tables are released, so the population fits cannot be checked in detail.\n\nBottom line: a useful atlas and a clear demonstration that symmetry breaking matters. The physical interpretation should be framed as speculative. A referee should insist on a calibration of the sheet locus against a dipolar force-free solution and a release of the fitting tables. I would not cite it as evidence, but it deserves a serious referee and would make a good reading-group discussion on the dangers of fitting an unvalidated model to a large sample.","headline":"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.","tokens_in":15752,"tokens_out":4261,"would_cite":false,"duration_ms":43745,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["gamma-ray pulsars","light-curve morphology","striped wind","current sheet","pair plasma loading","spherical harmonics","split monopole","pulse profile"],"falsifier":"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.","tokens_in":14684,"feed_emoji":"⚡","tokens_out":10038,"duration_ms":91980,"temperature":0.7,"pith_summary":"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.","feed_headline":"Asymmetric polar-cap loading explains gamma-ray pulse shapes","feed_subtitle":"An atlas of 16.5 million light curves shows the striped wind's pair-loading pattern encodes most observed pulse profiles.","key_machinery":"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","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Current-sheet geometry reproduces pulsar light-curve diversity","Asymmetric pair loading encodes pulsar pulse shapes","Striped wind model fits broad range of gamma-ray pulses","One current sheet explains many pulsar light-curve forms","Emissivity asymmetry in current sheet shapes pulsar pulses"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Current-sheet geometry reproduces pulsar light-curve diversity","Asymmetric pair loading encodes pulsar pulse shapes","Striped wind model fits broad range of gamma-ray pulses","One current sheet explains many pulsar light-curve forms","Emissivity asymmetry in current sheet shapes pulsar pulses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000209,"raw_usage":{"total_tokens":1277,"prompt_tokens":813,"completion_tokens":464,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":384}},"tokens_in":557,"tokens_out":464,"duration_ms":5324,"temperature":1.0,"reasoning_tokens":384,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T13:12:29.414155+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}