{"id":"184550d7-89f6-4d61-91ae-bb418e814b1d","arxiv_id":"2505.11408","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A self-consistent vacuum-gap model with a coordinate-dependent gap height, solved with physics-informed neural networks, gives the first quantitative non-axisymmetric accelerating potential and secondary plasma density profiles for orthogonal pulsars, and shows inverse Compton scattering is…","lead":"The authors present a new numerical method for computing the electric potential and the density of secondary plasma above the polar caps of radio pulsars, including orthogonal pulsars whose magnetic axis is nearly perpendicular to the spin axis. The method could improve models of radio pulse shapes and interpulse statistics.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central quantitative claim depends on the stationary vacuum gap model, which the authors acknowledge is contradicted by PIC simulations; without evidence that the stationary solution approximates the time-averaged gap, the computed potentials and density profiles are not validated for real…","rationale":"The single most load-bearing premise is the validity of the stationary vacuum gap model for producing quantitative predictions about orthogonal pulsars. The authors themselves acknowledge the non-stationarity of pair generation in the same sections where they formulate the model and state the central claim. The manuscript even mentions that plasma periodically leaves the magnetosphere, which suggests the gap is not in a steady state at any time. If the real gap is time-dependent, the computed H_gap and psi are not the physical values, and the density profiles derived from them lose quantitative meaning. This concern directly undermines the strongest claim of the paper and the abstract's implication that the method yields accurate structures for real pulsars. The reader identified the same weakest assumption, and I agree. A concrete PIC comparison would settle whether the stationary solution is a usable approximation; until then, the conditional verdict is appropriate.","tokens_in":13890,"tokens_out":12611,"duration_ms":130915,"concrete_test":"Run a 3D PIC simulation of an orthogonal pulsar (e.g., P=0.3 s, B=1.5e12 G, chi=89.3 deg) with the same parameters as Figure 4, and time-average the accelerating potential and plasma density in the polar cap region. Compare the time-averaged gap height and potential profile with the stationary solution presented here. If the time-averaged potential differs by more than ~30% in the pair-formation region, or if the gap height varies significantly over the discharge cycle, the stationary assumption fails and the claimed quantitative accuracy is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section II the authors state: 'Although the results of numerical simulations indicate that the particle generation is essentially a time-dependent process [13,16], we chose the stationary vacuum gap model [7] as a starting point.' The Conclusion repeats that 'many numerical [10,13,16] and analytical [44] studies indicate a significant non-stationarity of secondary plasma generation.' The central claim—'quantitatively determine the non-axisymmetric structure of the accelerating potential for orthogonal pulsars for the first time'—requires the stationary gap solution to be physically representative. But PIC simulations (Timokhin & Harding 2015; Philippov et al. 2015; Tolman et al. 2022) show that pair creation occurs in episodic discharges with plasma periodically leaving the magnetosphere; the gap is not in a stationary state, and the screening electric field is time-dependent. The paper provides no argument that the stationary solution equals the time average or that non-stationarity is a small correction. The gap height H_gap(rm,phi_m) and potential psi(rm,phi_m,z) are then solutions of a model that may not describe real orthogonal pulsars, so the word 'quantitative' is not justified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a method for computing the accelerating electric potential above the polar caps of radio pulsars with arbitrary inclination angle, using a stationary vacuum-gap model in which the gap height is coordinate-dependent and determined self-consistently with the potential. The Poisson equation is solved with physics-informed neural networks (PINNs) inside an iterative scheme, and the approach is applied to non-orthogonal and orthogonal pulsars, yielding non-axisymmetric potential structures and gap heights for the latter. From these potentials, the authors compute transverse secondary plasma density profiles using a synchrotron-cascade model, argue that inverse Compton scattering is negligible in the regimes considered, and discuss implications for the statistics of orthogonal interpulse pulsars.","tokens_in":14111,"tokens_out":6956,"duration_ms":61455,"significance":"If the results were validated, the method would be a useful extension of the classical Ruderman–Sutherland vacuum-gap analysis to non-axisymmetric geometries, potentially providing input for radio emission profile and interpulse visibility models. The iterative self-consistent treatment of a coordinate-dependent gap height is a sensible formulation, and the analysis of inverse Compton scattering is a useful addition. The paper is also honest about the model's limitations. However, the central claim of having 'quantitatively determined' the non-axisymmetric accelerating potential for orthogonal pulsars rests on two unverified pillars: the representativeness of the stationary vacuum-gap model for a physically non-stationary pair discharge, and the numerical accuracy of the PINN solutions. The manuscript currently lacks the independent checks needed to support this claim.","major_comments":[{"comment":"The central quantitative claim is undermined by the acknowledged non-stationarity of the physical process. In Section II the authors state that 'the particle generation is essentially a time-dependent process [13,16]' yet they choose the stationary vacuum gap model 'as a starting point', and in the Conclusion they repeat that 'many numerical [10,13,16] and analytical [44] studies indicate a significant non-stationarity of secondary plasma generation.' No argument or test is provided that the stationary gap height H_gap(r_m, phi_m) and potential psi(r_m, phi_m, z) correspond to the time-averaged quantities of the real, episodic discharges seen in PIC simulations (e.g., Timokhin & Harding 2015; Philippov et al. 2015; Tolman et al. 2022). The paper should include a quantitative comparison of the computed gap height and potential with time-averaged PIC results for similar parameters, or at least a timescale argument showing that non-stationarity is a small correction. Without such evidence, the word 'quantitative' in the claim 'quantitatively determine the non-axisymmetric structure of the accelerating potential for orthogonal pulsars for the first time' is not justified.","section":"Section II"},{"comment":"The PINN solution is not independently verified. The authors state that an analytical series solution exists but is not used due to a large condition number, and that grid-based methods are poorly suited to the variable-shaped domain. However, no comparison of the PINN result against a known solution is presented. The reported 'relative errors' in Fig. 3(b) and Fig. 4 are norms of the Poisson-equation residuals, which do not guarantee the accuracy of the potential itself; a small PDE residual can coexist with a poorly satisfied boundary condition (8) or with a spurious solution. The manuscript should provide a comparison of the PINN potential with the analytical series [22] (or a conventional finite-difference solver) for at least one axisymmetric case, and should report the residual of the Neumann condition (8) on the converged gap surface. This is necessary to support the claimed 1–2% accuracy.","section":"Section IIIC"},{"comment":"The convergence of the iterative scheme (10) is not established. The scheme introduces a weighting factor w with no analysis of its value or sensitivity, and the stopping criterion is a 1% relative change in the gap height between iterations. No proof or numerical demonstration is given that the fixed-point iteration converges to a unique self-consistent solution, nor that the converged solution is independent of w and of the initial guess. This is load-bearing because the gap height H_gap(r_m, phi_m) is one of the principal outputs and is used in the density calculations. The authors should show results for several values of w and confirm that the converged H_gap and psi are insensitive to these choices, or provide a contraction argument for the mapping H_gap[psi].","section":"Section III"}],"minor_comments":[{"comment":"The phrase 'border conditions' is used repeatedly; it should be replaced with 'boundary conditions' throughout.","section":"Global"},{"comment":"The caption refers to the 'monoenergetic approximation of the synchrotron spectrum (14)', but Eq. (14) defines the characteristic energy of curvature radiation; 'curvature' should replace 'synchrotron'.","section":"Fig. 3 caption"},{"comment":"The quantity Λ is used in Eq. (12) before it is defined in Eq. (13); consider reordering the definitions or adding a parenthetical note.","section":"Eq. (12) and Eq. (13)"},{"comment":"The label 'vacuum gap height normalized to polar gap radius' likely should read 'polar cap radius' for consistency with the text.","section":"Fig. 4"},{"comment":"References [12] and [40] are the same paper (Tchekhovskoy et al. 2016); one duplicate entry should be removed.","section":"References"},{"comment":"In Eq. (35), the quantity cosθ_b is used without restating its definition from Eq. (5); a brief reminder would improve readability.","section":"Section IV, Eq. (35)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a translation of a Russian-language paper (Astronomicheskii Zhurnal) and is reasonably clear, though language editing would help. The key concern is the gap between the acknowledged non-stationarity of the physical system and the 'quantitative' claim made from a stationary model; the requested comparisons with PIC time-averaged results or a clear downgrade of the claim would make the paper publishable as a methodological contribution. The numeric validation points (PINN accuracy, iteration convergence) are standard referee requests and should be addressable in revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The main thing to know: this is a serious modeling paper with a genuinely new computational approach—self-consistent gap height and Poisson equation solved via a PINN—that delivers the first non-axisymmetric potential maps for orthogonal pulsars. But the physical validity of the central result hinges on an assumption the authors themselves admit is questionable: the stationary vacuum gap. If you hope to use these density profiles to interpret MeerKAT/FAST data, treat them as illustrative of a particular model, not as firm predictions.\n\nWhat I liked: the iteration between gap height and potential is sensible, and they actually implement it on a 3D domain where the boundary is free and shape-changing. PINN is a reasonable choice there. They report 1–2% residual errors and test the monoenergetic approximation’s sensitivity (15–20% changes for a factor-two energy shift). They also show inverse Compton scattering is negligible across their parameter range, which is useful. The density profiles carry interesting qualitative predictions—pair multiplicity suppressed by orders of magnitude near χ=90°, and g(r_m) ∼ exp(−a²/r_m²) near the axis—worth testing against observations.\n\nThe soft spots are real but not fatal to the paper as a modeling exercise. First, the stationary gap assumption. They cite Timokhin & Harding and Philippov et al.’s episodic, time-dependent pair cascades and offer no argument that their steady solution is the time average or a small correction. That is the key weakness, and it directly weakens the word “quantitative” in their central claim. Second, they never validate the PINN solver against the analytical series solution they cite in [22], even in the axisymmetric limit where comparison would be trivial. Third, the iterative scheme’s convergence is empirical; the weighting factor w is introduced without discussion of its tuning or robustness. Fourth, no code or data is deposited, so reproducibility rests on the figures.\n\nThat said, I don’t see a fatal internal error. The equations are standard, the derivation is clean, and the authors are appropriately cautious about the model’s scope in the conclusion. The paper deserves a serious referee, but it should be sent back with a request for solver validation against the analytical solution, a sensitivity analysis of the stationary assumption (even a rough comparison with PIC time-averaged profiles), and a clearer account of convergence.\n\nIt is a solid contribution for the pulsar community, especially for people working on interpulse statistics and profile modeling. I wouldn’t cite it as ground truth, but I’d want it in the literature with the caveats highlighted.","headline":"A well-executed new computational method for non-axisymmetric pulsar vacuum gaps, but the physical claim of quantitative accuracy is undermined by the unvalidated stationary gap assumption and missing solver validation; worth refereeing as a modeling contribution.","tokens_in":14626,"tokens_out":2917,"would_cite":false,"duration_ms":30422,"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":"The paper proposes a self-consistent vacuum-gap method that yields the first quantitative non-axisymmetric accelerating potential and secondary plasma density profiles for orthogonal pulsars.","keywords":["radio pulsars","accelerating potential","secondary plasma density","vacuum gap","orthogonal pulsars","physics-informed neural networks","pair cascade","interpulse pulsars"],"falsifier":"A particle-in-cell simulation of an orthogonal pulsar with $B_{12}=1.5$, $P=0.3$ s, and $\\chi=88^\\circ$ that resolves the time-averaged secondary plasma density above the polar cap would settle the point: if the gap height and density profile oscillate substantially, the stationary self-consistent solution cannot describe the real magnetosphere. Alternatively, a statistical study of orthogonal interpulse pulsars with independently estimated surface fields could test the predicted drop in pair multiplicity for $89^\\circ\\lesssim\\chi\\lesssim91^\\circ$.","tokens_in":13690,"feed_emoji":"⚡","tokens_out":8373,"duration_ms":74415,"temperature":0.7,"pith_summary":"The paper argues that the accelerating potential above a pulsar's polar cap can be computed without the usual assumptions that the vacuum gap is thin and axisymmetric. The authors build a self-consistent scheme in which the height and shape of the vacuum gap and the electric potential inside it are solved together, and they use this scheme to obtain, for the first time, the non-axisymmetric accelerating potential of orthogonal pulsars (magnetic axis nearly perpendicular to the rotation axis). From that potential they compute the transverse profiles of the secondary electron-positron plasma flowing along open field lines. These density profiles are what a quantitative theory of observed radio profiles and interpulses needs, since propagation effects such as refraction and absorption depend on the actual plasma distribution. The paper also concludes that inverse Compton scattering is negligible relative to curvature radiation for the polar cap temperatures and parameters considered.","feed_headline":"First map of the accelerating voltage above orthogonal pulsar poles","feed_subtitle":"Self-consistent vacuum-gap model gives the density profiles needed to interpret radio interpulses.","key_machinery":"The central object is the vacuum gap: the region above the polar cap in which the charge density is much smaller than the Goldreich-Julian value, so a longitudinal electric field accelerates primary particles. The argument is carried by a self-consistent iteration in which the potential $\\psi$ is obtained from the Poisson equation in a domain whose upper boundary $H_{\\rm gap}(r_m,\\varphi_m)$ is itself unknown, and the boundary is recalculated from the condition that curvature-radiation photons produced by the accelerated primaries become able to convert into pairs. The Poisson equation is solved with a physics-informed neural network (a fully connected network trained to minimize the PDE residual) because the domain has a variable shape and the alternative series solution is numerically ill-conditioned. Once the potential is known, the secondary plasma density follows from a synchrotron cascade calculation of the number of pairs produced per primary particle, with the cascade spectrum obtained by solving an integral equation for the photon number distribution.","core_discovery":"On the paper's own terms, the central discovery is that a stationary vacuum-gap model, with the gap height $H_{\\rm gap}(r_m,\\varphi_m)$ treated as a free function determined together with the potential $\\psi(r_m,\\varphi_m,z)$, quantitatively reproduces the non-axisymmetric structure of the accelerating potential for orthogonal pulsars. Solving the Poisson equation with zero potential on the stellar surface and the separatrix and zero parallel electric field at the upper boundary, the authors find that for inclination angles near $90^\\circ$ the Goldreich-Julian charge density can change sign inside the polar cap, and with it the sign of the accelerating potential, over a narrow range roughly $88.5^\\circ\\lesssim\\chi\\lesssim91.5^\\circ$. From these potentials they compute the multiplicity and transverse density profile of secondary plasma, and find that for $\\chi\\simeq89^\\circ$–$91^\\circ$ and moderate fields $B_{12}\\sim2$ the pair multiplicity is orders of magnitude lower than for ordinary pulsars, while for orthogonal interpulse pulsars with $B_{12}=7$ plasma generation is stronger. Inverse Compton scattering, both resonant and non-resonant, is shown to be unimportant for the gap height and multiplicity in the parameter ranges studied.","pith_inferences":["The authors note the stationary-gap assumption as a limitation; I infer that the most direct test of the method is a time-resolved particle-in-cell simulation of the same orthogonal-pulsar parameters ($B_{12}=1.5$, $P=0.3$ s, $\\chi=88^\\circ$ and $89.3^\\circ$) that time-averages the pair density profile and compares it with the stationary prediction.","The predicted suppression of multiplicity near $\\chi=90^\\circ$ could be tested statistically: if the method is right, the fraction of interpulse pulsars should decline toward periods and period-derivatives where the inferred surface field drops below a few $10^{12}$ G, unless an evolution model that predicts enhanced fields applies.","Nothing in the method itself requires a dipolar field; the authors assume dipole geometry, but the same self-consistent scheme could be applied to non-dipolar or multipolar near-surface fields. I infer that including such fields could resolve the orthogonal-pulsar contradiction the authors identify, since local curvature and field strength would change pair-production rates."],"forward_implications":["For non-orthogonal pulsars ($\\chi\\lesssim85^\\circ$) the problem reduces to an axisymmetric one, and the computed potential reproduces the expected slot-gap geometry: the gap height diverges near the magnetic axis and near the polar cap edge.","For orthogonal pulsars the accelerating potential is genuinely non-axisymmetric, and within the narrow band $88.5^\\circ\\lesssim\\chi\\lesssim91.5^\\circ$ the Goldreich-Julian density changes sign in the polar cap, so the accelerating electric field reverses direction across the cap.","The pair multiplicity for orthogonal pulsars with $B_{12}\\sim2$ is several orders of magnitude lower than for ordinary pulsars, which is hard to reconcile with the observed abundance of orthogonal interpulse pulsars unless their magnetic fields are higher or the inclination estimates are biased.","The transverse density profile of secondary plasma near the magnetic axis behaves as $g(r_m)\\propto\\exp(-a^2/r_m^2)$ for $r_m<0.03$, while on the larger scale $0.03