{"id":"8931659d-ba7c-456a-8bc6-7f4990a75463","arxiv_id":"2506.12423","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"FAST and MeerKAT pulsar surveys show O-mode mean profiles are wider than X-mode profiles and that orthogonal interpulse pulsar statistics favor the classical BGI vacuum-gap evolutionary model over MHD models.","lead":"This paper compares new radio pulsar polarization data from FAST and MeerKAT against theoretical predictions for two wave modes, and argues that the fraction and periods of orthogonal interpulse pulsars support a specific evolutionary model. The authors report that O-mode profiles are wider than X-mode profiles and that the observed statistics favor the classical Ruderman-Sutherland (BGI) vacuum-gap model, where inclination angles evolve toward 90 degrees.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Interpulse-fraction test is circular and underpowered: model predictions come from in-house fits (death-line constant k fitted in Novoselov+2020), the FAST fraction (6%, 14/233) sits at/above the BGI band, and total counts are small.","rationale":"The reader's stated weakest assumption is the mode-classification sign relation, which affects only the first (width) conclusion. I agree that is a real concern—the 5% cross-check compares two applications of the same theoretical rule, not an independent validation—but the second, more consequential conclusion about pulsar evolution rests on a shakier foundation: a circular, in-house, post-fit model comparison with small number counts. The paper explicitly acknowledges it cannot directly test theoretical beam widths, so the width claim is modest; the evolutionary claim is the headline and the one that 'challenges the currently accepted belief.' Setting the verdict to UNCHANGED reflects that the reader's conditional verdict already demands the external, model-independent calculation that would settle this concern. I do not see an internally inconsistent calculation, nor a reason to reject outright; the paper is honest about several limitations. Credit is due for using homogeneous FAST and MeerKAT samples and for reporting both automated and manual mode classifications, which strengthens the width result within its stated scope. The missing piece is an independent population synthesis for interpulse fractions, which is exactly the condition the reader set.","tokens_in":11877,"tokens_out":8823,"duration_ms":109791,"concrete_test":"Have an independent group (not affiliated with the BGI authors) run a population synthesis for orthogonal interpulse pulsars using a standard MHD torque/current model from 3D PIC simulations (e.g., Philippov et al. 2014) and a plausible beaming model, without fitting any death-line constant to the FAST/MeerKAT interpulse data. Compare the resulting predicted fraction and period distribution in 0.033–0.5 s to Tables 3 and 4. If the MHD prediction rises above ~2% or the BGI prediction falls below ~3%, the claimed discrimination collapses.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's most consequential claim—that FAST/MeerKAT interpulse statistics confirm the BGI/Ruderman–Sutherland model and χ→90°—depends on comparing Table 3 and Table 4 to predictions from Novoselov et al. (2020). That comparison is not an independent test: Novoselov et al. is by the same group and calibrates the death-line constant k in equation (7) to similar data (kBGI=0.5±0.1, kMHD=1.0±0.1). The stated BGI range 2.5–5.5% and MHD upper bound <1% are therefore post-fit, not true predictions. The new FAST fraction in the key period range is 6.0% (14/233), which exceeds the upper end of the BGI band, and MeerKAT gives 4.2% (25/590); the combined counts are tens of objects, so Poisson and survey-selection uncertainties are large. The paper itself labels the FAST count a lower limit, so the true fraction could move further from the BGI band. Meanwhile the MHD 'prediction' is derived from the authors' own interpretation of MHD current/potential structure rather than an independent simulation-based population synthesis. Without an external, uncalibrated calculation of interpulse fractions under both models, the claimed confirmation of the vacuum-gap evolutionary model is not established.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses polarization data for 682 FAST and 1170 MeerKAT pulsars to classify mean pulse profiles as ordinary (O) or extraordinary (X) mode via the theoretical sign relation between the derivative of the position angle and the Stokes V parameter. It reports that O-mode profiles are significantly wider than X-mode profiles (Tables 1 and 2), and compiles orthogonal interpulse pulsar counts from FAST, MeerKAT, and BSA (Tables 3 and 4) to argue that the observed fraction and period distribution favor the BGI/Ruderman-Sutherland evolutionary model, in which the inclination angle tends toward 90 degrees, over the MHD model, in which it decreases. The paper also discusses the death-line relation (7) and the predicted period peak for orthogonal pulsars around 0.1–0.3 s.","tokens_in":12173,"tokens_out":7119,"duration_ms":82892,"significance":"If the O-mode width difference is robust, it provides a valuable large-sample confirmation of propagation effects in pulsar magnetospheres; the statistical tests are appropriate, and the use of both manual and automated classification is a strength. The interpulse counts compiled from homogeneous surveys are a useful observational resource. However, the evolutionary-model comparison is not an independent test: the predicted fractions are taken from Novoselov et al. (2020), which calibrates the death-line constant to pulsar data, and the paper explicitly leaves unmodeled the viewing-angle and beam-intensity distributions needed for a quantitative comparison with Eqs. (1)–(3). The paper is therefore best read as presenting strong new observational constraints and a qualitative consistency check, not a decisive test between evolution models.","major_comments":[{"comment":"The central claim that FAST and MeerKAT data confirm the BGI model over the MHD model is not established as stated, because the model predictions are not independent of the data used to calibrate them. The predicted fractions (BGI 2.5–5.5%, MHD <1%) are taken from Novoselov et al. (2020), and the paper itself notes in the next section that Novoselov et al. fitted the death-line constant in Eq. (7) to the pulsar distribution, obtaining kBGI = 0.5 ± 0.1 and kMHD = 1.0 ± 0.1. A post-fit comparison of a new sample to a model whose parameter was tuned to earlier pulsar data cannot be presented as confirmation. Moreover, the new FAST value 14/233 = 6.0% lies above the stated BGI upper bound of 5.5%, and the paper labels the FAST count as a lower limit, so even the direction of the discrepancy is not controlled. The authors should either obtain predictions from an independently calibrated model or explicitly recast the result as a consistency check with large systematic uncertainty.","section":"Statistics of orthogonal interpulse pulsars: Total number, Tables 3–4"},{"comment":"The paper stops short of a quantitative test of the theoretical widths. It states that direct comparison of the measured W10 P^{1/2} distributions with Eqs. (1)–(3) requires modeling the unknown viewing and inclination angle distributions and non-uniform beam patterns, and that such a model is beyond the scope of the paper. In that case, the abstract's and Discussion's statements that the data 'fully correspond to the predictions' and that the theory has received 'confident confirmation' go beyond what the statistical tests establish: the tests in Table 2 show only that the X- and O-mode width distributions differ in the predicted direction, not that their magnitudes match Eqs. (1)–(3).","section":"Widths of the radio pulsar mean profiles, after Eq. (3) and Discussion"},{"comment":"The claimed <5% mode-misclassification rate is estimated from the agreement between the manual and automated classifications, both of which implement the same theoretical sign rule (sign V = ± sign dp.a./dφ). This does not validate the sign rule itself. Since the entire width comparison rests on the physical identification of O and X modes, the paper would be strengthened by an external check, such as comparison with independent mode diagnostics or a test on simulated polarization profiles with known modes; without such a check, a systematic failure of the sign rule for a subset of pulsars could bias the comparison.","section":"Widths of the radio pulsar mean profiles, Eq. (5) and threshold eta_cr"},{"comment":"The statistical evidence distinguishing BGI from MHD is weaker than the text implies. The counts are small (14/233 for FAST and 25/590 for MeerKAT in the key period range), and the FAST sample is explicitly a lower limit because faint interpulses may be missing from Fig. A6 of Wang et al. (2023). With Clopper-Pearson 95% confidence intervals, the FAST fraction is consistent with values up to about 10%, which would be difficult to reconcile with the BGI band, while the MeerKAT fraction alone is consistent with the entire BGI band. The authors should report confidence intervals for the observed fractions and discuss survey selection effects before concluding that the models are distinguished at high confidence.","section":"Statistics of orthogonal interpulse pulsars: Total number, Table 3"}],"minor_comments":[{"comment":"The model is called 'BIG' in the Introduction and 'BGI' elsewhere; use one abbreviation consistently.","section":"Introduction"},{"comment":"The header 'P (с)' uses a Cyrillic 'с' instead of the Latin 's'.","section":"Table 4"},{"comment":"The sentence contains a duplicated citation: 'Novoselov et al. Novoselov et al. (2020)' should be 'Novoselov et al. (2020)'.","section":"Section 4"},{"comment":"Table 1 would be easier to read if the four subcategories (Xs, Xd, Os, Od) were arranged as separate columns with clear headers; the current layout is difficult to parse.","section":"Table 1"},{"comment":"The notation is inconsistent: the threshold is introduced as eta_cr but the text then says 'the value of eta was chosen to be eta = 0.4'; use eta_cr consistently.","section":"Widths of the radio pulsar mean profiles"},{"comment":"The summation condition should be typeset as I >= 0.1 I_max for clarity.","section":"Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The second half of the paper is essentially a re-affirmation of the authors' own prior results (Novoselov et al. 2020, Istomin et al. 2024) using new catalogs. The new observational numbers are worth publishing, but the interpretive claims need to be separated from the data presentation and substantially toned down or backed by an independent model calculation. This is a scope and calibration concern, not a question of author intent."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe part of this paper worth your time is the first half: the O-mode/X-mode width comparison on the FAST and MeerKAT samples. The mode classification is transparent—manual plus an automated eta criterion with a stated 5% misclassification target, applied to public polarization profiles—and the Anderson–Darling and permutation tests genuinely reject equal distributions. The paper also admits that the quantitative width relations (1)–(3) cannot be tested directly without modeling the unknown viewing-angle and beam-intensity distributions. That honesty is to its credit. The result that O-mode profiles are wider than X-mode, especially for double-humped profiles, is a solid observational statement, modulo the assumed sign relation between d.p.a./dphi and V that defines O vs X.\n\nThe second half, on orthogonal interpulse pulsars, is much weaker than the abstract's 'confirm the validity' language. The comparison fractions in Table 3 are tested against predictions from Novoselov et al. (2020), which is the same group's model with the death-line constant k fitted to similar data. That is not an independent test. The MHD 'prediction' of <1% is also the authors' own interpretation of MHD potential drops, not an independent population synthesis. And the numbers themselves are not clean: the FAST fraction in the key period range is 6.0%, above the upper end of the BGI band (2.5–5.5%), and the paper concedes it is a lower limit, so the true value could move further out. The combined sample is a few tens of objects, so the statistical constraints are loose. The period distribution in Table 4 is suggestive but again compared to in-house BGI predictions. The death line relation (7) is treated as established based on the same group's earlier fit, which is circular.\n\nSo: the paper ships a reproducible classification and a new catalog-based measurement—that is real. But the interpulse-evolution conclusion is not supported by the evidence presented. A referee should ask for (a) an external, uncalibrated calculation of interpulse fractions for both evolutionary models, (b) a forward model of the width distributions rather than medians alone, and (c) language that matches the actual strength of the tests.\n\nThis paper deserves a serious referee: the width result is worth publishing, and the interpulse question is important even if this attempt doesn't settle it. I would send it to review with a major-revision recommendation, not desk-reject it.","headline":"The O-mode/X-mode width difference on the FAST and MeerKAT samples is a solid, reproducible result; the interpulse-pulsar argument for the BGI model and chi -> 90 deg is overinterpreted and rests on in-house predictions that are not an independent test.","tokens_in":12723,"tokens_out":3288,"would_cite":true,"duration_ms":39415,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Using FAST and MeerKAT polarization samples, the paper confirms that ordinary-mode pulsar profiles are wider than extraordinary-mode ones and that orthogonal interpulse statistics favor the vacuum-gap model of magnetic-axis evolution…","keywords":["radio pulsars","mean pulse profiles","polarization modes","ordinary and extraordinary modes","orthogonal interpulse pulsars","pulsar evolution","death line"],"falsifier":"A complete, sensitivity-matched survey that turned up several orthogonal interpulse pulsars with periods above 0.5 s at a rate comparable to the 0.1–0.3 s peak would contradict the death-line prediction and support the competing alignment model; the absence of such long-period orthogonal pulsars would strengthen the paper's conclusion.","tokens_in":11660,"feed_emoji":"📡","tokens_out":20541,"duration_ms":210634,"temperature":0.7,"pith_summary":"Using the large polarized samples from the FAST and MeerKAT surveys, the paper re-examines two basic predictions of pulsar theory. It finds that mean profiles emitted in the ordinary (O) mode are statistically wider than those emitted in the extraordinary (X) mode, with median scaled widths $W_{10}P^{1/2}$ about 15 versus 12 degrees in FAST and 15 versus 11 degrees in MeerKAT, matching the expectation that O-mode rays are refracted and broadened in the magnetosphere. It also finds that the fraction of orthogonal interpulse pulsars in homogeneous samples, about 4–6 percent for periods $0.033<P<0.5$ s, is far above the MHD model's $\\lesssim1$ percent prediction and consistent with the classical vacuum-gap (BGI) model, in which the inclination angle $\\chi$ evolves toward $90^\\circ$. Their period distribution, peaking near 0.1–0.3 s and almost vanishing above 0.5 s, matches a death-line restriction on orthogonal rotators in the BGI model. If both conclusions hold, the common interpretation of pulsar spin-axis alignment with age would need to be revised.","feed_headline":"Pulsar data confirm wider O-mode beams and 90-degree axis evolution","feed_subtitle":"Large polarization samples favor the classical vacuum-gap model over spin-alignment models.","key_machinery":"Two pieces of theory carry the argument. The first is the mode-classification sign rule: in pulsar mean profiles, the sign of the derivative of the position-angle swing $\\mathrm{d\\,p.a.}/\\mathrm{d}\\varphi$ is the same as the sign of circular polarization $V$ for the extraordinary mode and opposite for the ordinary mode; the paper converts this into a normalized score $\\eta$, the mean of sign($V$) times sign($\\mathrm{d\\,p.a.}/\\mathrm{d}\\varphi$) over points with intensity above 10 percent of the maximum, and classifies pulsars as X-mode for $\\eta>0.4$, O-mode for $\\eta<-0.4$, with misclassification estimated below 5 percent. The second is the death line of the classical vacuum-gap (BGI) model, $\\cos\\chi> k\\,P^{15/7}B^{-8/7}$ with $k\\approx1$, which for orthogonal rotators (small $\\cos\\chi$) implies a period cap $P<0.2\\,B_{12}^{16/37}$ s. Together these tools let the paper separate the two mode populations by profile width and use the abundance and period distribution of orthogonal interpulse pulsars to discriminate between the BGI and MHD evolutionary models.","core_discovery":"At the core of the paper is a sign rule for pulsar polarization: for a mean profile, the sign of the derivative of the position angle $\\mathrm{d\\,p.a.}/\\mathrm{d}\\varphi$ relative to the sign of the Stokes $V$ parameter identifies the mode—same sign is the extraordinary mode, opposite sign the ordinary mode. Applying this rule to the FAST and MeerKAT samples gives median widths $W_{10}P^{1/2}$ of about $12.5^\\circ$ (X) versus $15.4^\\circ$ (O) for FAST and $10.6^\\circ$ (X) versus $15.3^\\circ$ (O) for MeerKAT, and Anderson–Darling and permutation tests reject equal distributions in every comparison. For orthogonal interpulse pulsars, the homogeneous samples give relative fractions of $6.0\\%$ and $4.2\\%$ in the period range $0.033<P<0.5$ s, far above the MHD model's $\\lesssim1\\%$ prediction; the period distribution, peaking near $0.1$–$0.3$ s and almost absent above $0.5$ s, matches the death-line restriction that in the vacuum-gap (BGI) model limits orthogonal rotators to short periods. The authors conclude that the classical vacuum-gap evolutionary model, in which the inclination angle $\\chi$ increases toward $90^\\circ$, is favored over the MHD model, and that the observed decline of average inclination angle with period is a selection effect of the death line rather than evidence that individual pulsars align with age.","pith_inferences":["A direct test the paper leaves implicit: compare the $\\eta=0.4$ mode classification against independent mode identifications from orthogonal-mode jumps to measure the true misclassification rate on a pulsar-by-pulsar basis.","The authors do not split their sample by magnetic-field strength, but the death-line expression implies that high-field orthogonal pulsars should be able to radiate at longer periods; checking that trend would sharpen the model comparison.","The FAST sample contains an 11 percent orthogonal-interpulse fraction below $P<0.033$ s, which the paper sets aside because millisecond pulsars evolve differently; whether recycled pulsars obey the same death line is an open extension of the analysis."],"forward_implications":["The statistically robust width difference between O- and X-mode profiles directly confirms the refraction picture in which ordinary-mode rays are broadened as they traverse the pulsar magnetosphere.","The observed fraction of orthogonal interpulse pulsars, 4–6 percent for periods between 0.033 and 0.5 s, is too high for the MHD evolutionary model's at-most-1 percent prediction, so that model's monotonic decrease of inclination angle is disfavored.","The scarcity of orthogonal interpulse pulsars with periods above 0.5 s supports the death-line restriction specific to the vacuum-gap model and contradicts the MHD expectation of comparable numbers at $P\\sim1$ s.","The decline of average inclination angle with pulsar period can be produced by the death-line selection of short-period objects with small $\\cos\\chi$, so this observed trend does not by itself prove that individual pulsars align with age."],"supporting_citations":[{"why":"Establishes the sign relation between the position-angle derivative and circular polarization that the paper uses to sort pulsars into X- and O-modes.","marker":"Andrianov and Beskin (2010)"},{"why":"Develops the propagation theory predicting O-mode broadening and links the sign rule to mean-profile morphology.","marker":"Beskin and Philippov (2012)"},{"why":"Extends the sign-rule method to complex mean profiles, supporting the paper's manual and automated classification.","marker":"Hakobyan et al. (2017)"},{"why":"Provides the theoretical width scalings for X-mode and wider O-mode profiles that the FAST and MeerKAT data are compared against.","marker":"Beskin et al. (1988)"},{"why":"Supplies the FAST polarization database of 682 pulsars used for both the width comparison and the interpulse sample.","marker":"Wang et al. (2023)"},{"why":"Supplies the MeerKAT thousand-pulsar-array polarization data used for widths and interpulse statistics.","marker":"Johnston et al. (2023)"},{"why":"MeerKAT pulse-width catalogue from which the paper draws part of the homogeneous orthogonal-interpulse sample.","marker":"Posselt et al. (2021)"},{"why":"Gives the earlier heterogeneous interpulse statistics and the BGI-model fraction predictions (2.5–5.5 percent) and MHD prediction (about 1 percent) that the new data are tested against.","marker":"Novoselov et al. (2020)"},{"why":"Defines the vacuum-gap accelerating potential and the death-line condition that underlies the BGI model's period and fraction predictions.","marker":"Ruderman and Sutherland (1975)"},{"why":"MHD magnetosphere simulations predicting that inclination angle decreases with time, the competing evolutionary model disfavored here.","marker":"Philippov et al. (2014)"}],"fun_headline_variants":["Pulsar polarization sign rule backs classical vacuum-gap model","FAST and MeerKAT data: O-mode beams wider, axes evolve to 90°","Survey data favor vacuum-gap over MHD for pulsar evolution","Pulsar mode widths and axis tilt confirm vacuum-gap theory","Orthogonal pulsars rare, supporting vacuum-gap evolution"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the sign of the swing of the polarization angle relative to the sign of the circular polarization correctly identifies the emission mode for almost every pulsar in the sample; if that sign relation fails for a substantial fraction, or if visual inspection biases which pulsars are kept, both the width comparison and the orthogonal-pulsar statistics lose their foundation.","fun_headline_variants_meta":{"raw":{"variants":["Pulsar polarization sign rule backs classical vacuum-gap model","FAST and MeerKAT data: O-mode beams wider, axes evolve to 90°","Survey data favor vacuum-gap over MHD for pulsar evolution","Pulsar mode widths and axis tilt confirm vacuum-gap theory","Orthogonal pulsars rare, supporting vacuum-gap evolution"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000622,"raw_usage":{"total_tokens":2908,"prompt_tokens":1000,"completion_tokens":1908,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":616,"completion_tokens_details":{"reasoning_tokens":1814}},"tokens_in":616,"tokens_out":1908,"duration_ms":13197,"temperature":1.0,"reasoning_tokens":1814,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:49:18.243301+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A complete, sensitivity-matched survey that turned up several orthogonal interpulse pulsars with periods above 0.5 s at a rate comparable to the 0.1–0.3 s peak would contradict the death-line prediction and support the competing alignment model; the absence of such long-period orthogonal pulsars would strengthen the paper's conclusion.","supporting_citations":[],"review_version":1}