{"id":"eca56ffe-2b6a-47c6-ac4b-9a44b8ecc376","arxiv_id":"2507.08509","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"Accretion mounds on neutron stars spread equator-ward as their mass grows, measurably burying the dipole field for low-field stars, and a new current-free boundary condition removes dependence on the outer domain size.","lead":"This paper simulates how matter piling onto a neutron star's magnetic pole bends and buries the star's magnetic field, using a new current-free outer boundary condition and high-resolution numerical solutions. For low-field neutron stars, the dipole moment can drop by about a third as the accretion mound spreads toward the equator, which is relevant to millisecond pulsar field decay and gravitational wave estimates.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Headline 37% dipole reduction is quoted at a resolution-dependent numerical maximum mass; convergence of this value is not demonstrated","rationale":"The paper's central claim as summarized in the strongest_claim has two components: (1) ring mounds spread equator-ward with increasing mass, and (2) for B_d=1e9 G, ζ=0.6 the dipole moment ratio reaches 0.627, demonstrating onset of field burial. Component (1) is supported by the resolution studies in Sec. 4.2.1 and Appendix C. Component (2), however, is quoted at the NMM, and Appendix C explicitly demonstrates that NMM is resolution-dependent (Table C1). Since the dipole moment ratio monotonically decreases with mass (Fig. 7 bottom, Fig. 9), the value 0.627 is the burial at the largest mass accessible at 5000^2 resolution, not at a physical or converged mass limit. The authors acknowledge this ('this NMM is limited by resolution'), but they do not provide a resolution-convergence test for the headline ratio. This is the most load-bearing weakness because if the ratio continues to drop at higher resolution, the paper's quantitative claim of '37% reduction' is an underestimate of the effect; if it rises, the claim is an overestimate. In either case, the current number is not a prediction. This is distinct from the reader's chosen weakest assumption about r0(ψ): that affects which physical system is being solved; the resolution issue affects whether the reported numerical result is even converged. Because the paper is candid about both, the CONDITIONAL verdict stands, but the condition should include a resolution-convergence check of μ1(Rout)/μ1(R*).","tokens_in":36801,"tokens_out":6385,"duration_ms":70805,"concrete_test":"Recompute the B_d=1e9 G, ζ=0.6 ring-shaped mound at NMM for three resolutions (5000, 10000, 20000) with otherwise identical setup; plot μ1(Rout)/μ1(R*) (Eq. 20) and the mound mass against resolution. If the ratio changes by more than 5% between 5000 and 10000, the 0.627 value is not converged and the quantitative field-burial claim is unresolved. Also check whether NMM converges to a finite mass; if it keeps increasing, the spreading solutions are cut off by the grid, not by physics.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative result, μ1(Rout)/μ1(R*) = 0.627 for B_d=1e9 G, ζ=0.6 (Fig. 9, Table 5), is obtained at the Numerical Maximum Mass (NMM) as defined in Appendix C. The NMM is not a physical mass limit: Appendix C and Table C1 show it increases monotonically with resolution, from 3.34e-13 M_sun at 700^2 to 5.65e-13 M_sun at 5000^2, and the authors state that NMM is 'limited by resolution.' Because the normalized dipole moment decreases monotonically with accreted mass (Figs. 7 and 9), a higher resolution that permits a larger NMM could produce a dipole moment ratio substantially below 0.627; conversely, if the trend reverses, the quoted burial would be overestimated. The paper does not report how μ1(Rout)/μ1(R*) changes with resolution for this case, so the specific '37% reduction' is not a converged numerical prediction. This is independent of the profile-function caveat: even keeping r0(ψ) fixed, the headline number is bound to a numerical cutoff. The qualitative spreading may survive, but the quantitative field-burial claim — the paper's primary evidence for low-field fast burial — is not yet established as a physical result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper solves the axisymmetric Grad-Shafranov (GS) equation for magnetically confined accretion mounds on neutron stars, using stretched spherical coordinates and a new multipolar current-free outer boundary condition. With a ring-shaped mound profile r0(psi) whose truncation angle is tied to the Alfven radius, the authors find that as mound mass increases, the mound spreads latitudinally toward the equator and the normalized dipole moment at the outer radius decreases. For B_d = 10^9 G and zeta = 0.6, the reported ratio is mu_1(Rout)/mu_1(R*) = 0.627, which they interpret as the onset of field burial. They also model mounds on a pre-existing ocean/envelope and quadrudipolar inner boundary configurations, which produce asymmetric polar mounds. The solver is validated against a Green's-function boundary condition (Appendix B3) and a perturbation-based uniqueness test defines a resolution-dependent Numerical Maximum Mass (NMM).","tokens_in":37105,"tokens_out":5048,"duration_ms":62220,"significance":"If the central result holds, the paper would overturn the common interpretation that closed magnetic loops set an intrinsic mass ceiling in GS mound calculations, and would support rapid field burial for low-field neutron stars in the early accretion phase. The new current-free boundary condition, the large parameter suite, the validation against a Green's-function volume boundary condition (relative difference below 5e-5 in Appendix B3), and the clearly described perturbation tests for numerical uniqueness are all valuable and reproducible contributions. However, because the headline dipole reduction is quoted at a resolution-limited numerical cutoff and the mass-loading profile is prescribed, the quantitative field-burial claim is not yet established as a converged physical prediction.","major_comments":[{"comment":"The central quantitative result, mu_1(Rout)/mu_1(R*) = 0.627 for B_d = 10^9 G and zeta = 0.6, is the value at the Numerical Maximum Mass, and NMM is not a physical mass limit. Table C1 shows that NMM increases monotonically with resolution, from 3.34e-13 M_sun at 700^2 to 5.65e-13 M_sun at 5000^2, and the text states that NMM is 'limited by resolution.' Because the dipole ratio decreases monotonically with accreted mass (Figs. 7 and 9), a higher-resolution run that supports a larger NMM could shift the ratio substantially below 0.627. The paper does not report how mu_1(Rout)/mu_1(R*) itself converges with resolution for this case, so the specific '37% reduction' is not a converged numerical prediction.","section":"Sec. 4.2.3, Table 5; Appendix C, Table C1"},{"comment":"The profile function r0(psi) in Eq. (16) is prescribed, and the limitations acknowledge that no fixed mass-loading per flux tube dM/dpsi is imposed between different solutions. Consequently, the sequence of solutions in Fig. 6 and the mass-dipole relations in Figs. 8-10 are not a single accretion sequence: each snapshot is an independent equilibrium with a different dM/dpsi. The equator-ward spreading and the 0.627 dipole ratio therefore measure properties of the chosen profile family, not of a physically self-consistent loading history. A model that fixes dM/dpsi during mound growth, or a detailed justification that the adopted r0(psi) follows from disk-magnetosphere interaction, is needed to support the field-burial interpretation.","section":"Sec. 2.5, Eq. (16); Limitations (ii)"},{"comment":"The text notes that matter beyond theta_t is supported by vacuum magnetic fields (or by a low-density envelope in the ocean model) and states that the dynamical stability of these spreading mounds should be verified. Since the central claim is that a mound can persist while spreading and burying the field, a static GS equilibrium alone does not demonstrate field burial unless the configuration is stable or the instability growth time exceeds the accretion time. Given the known interchange and pressure-driven instabilities of such mounds, the burial result should be presented conditional on a stability check.","section":"Sec. 4.2.1; Limitations (v)"},{"comment":"The authors report that general-relativistic GS solutions show approximately three times less screening than the Newtonian solutions (Rossetto et al. 2023). This directly affects the quantitative claim: the 37% reduction is a Newtonian value and may overestimate burial for the same mound parameters. The limitation is acknowledged, but the abstract and Section 4 present the Newtonian number without this caveat; the discussion in Section 6 should carry the GR caveat whenever the 0.627 figure is quoted.","section":"Limitation (iv)"}],"minor_comments":[{"comment":"The column header 'R_b (km)' lists values such as 10, 9.98, and 9.93, which could be confused with R*; clarify that the dipole moment ratio in the last column is measured with respect to the inner boundary R_b, not the stellar surface.","section":"Table 3"},{"comment":"The symbols alpha, epsilon, and h/r are used before their adopted values are introduced; the text should define these disk parameters explicitly when they first appear.","section":"Eq. (29)"},{"comment":"These fits are empirical descriptions of the authors' own computed points; the text should state explicitly that they are not derived relations and should not be extrapolated outside the fitted mass and height ranges.","section":"Eqs. (23), (25), (26)"},{"comment":"The phrase 'Masses with an order of 10^-8 M_sun can be simulated' should read 'masses of order 10^-8 M_sun' to avoid grammatical ambiguity.","section":"Sec. 4.1"}],"recommendation":"major_revision","confidential_remarks":"The main barrier to acceptance is the resolution dependence of the NMM. The authors should either provide a convergence study of mu_1(Rout)/mu_1(R*) at fixed r0(psi) or present the dipole-reduction result as a qualitative trend rather than a headline number. I do not see circularity: the dipole moment is an output, not an input, and the limitations are stated openly. The paper is within the journal's scope, but the quantitative burial claim needs revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this is a useful methods paper. The multipolar current-free outer boundary condition (Sec. 2.4.2, Eq. 15) is a real step beyond the fixed, outflow, and reduced-dipole options, and the validation is solid — a Green's function volume boundary check with relative difference below 5e-5, plus a domain-independence test that the older BCs fail. The solver is standard SOR but carefully documented, the grid stretching is sensible, and the resolution study (up to 10000^2) is honest. The perturbation-based uniqueness test for defining the Numerical Maximum Mass is a clear idea, and the result that earlier mass limits were partly numerical artifacts is worth having. The new ocean-mound and quadrudipolar profiles broaden the scope.\n\nThe soft spots are real but mostly acknowledged by the authors. The headline 37% dipole reduction for B_d = 1e9 G, zeta=0.6 is quoted at the NMM, and Table C1 shows that mass increases monotonically with resolution, from 3.34e-13 to 5.65e-13 M_sun. Since the dipole moment falls with mass, the 0.627 ratio is not a converged number; higher resolution could make it lower. The paper never reports how the dipole moment itself behaves with resolution, so that one number should be cited as evidence of onset of burial, not as a quantitative prediction. The hand-specified profile r0(psi) also gives different dM/dpsi for each solution, so the equator-ward spreading is a property of the assumed mass-loading family, not a self-consistent accretion history. The authors state this limitation. Newtonian gravity and unverified dynamical stability are also acknowledged.\n\nNo code or data is public, which limits reproducibility, but the appendices give enough detail to reimplement the solver.\n\nNone of this sinks the paper. The qualitative picture — ring-shaped mounds spreading equator-ward and burying a low dipole — is plausible and now at least demonstrated in a well-documented solver. The contribution is the BC and the numerical framework, not the specific 37%.\n\nWho is it for: neutron-star astro, accretion-induced field decay, mound equilibria, CGW ellipticity. It deserves a serious referee. Recommendation: send to peer review, with the request that the authors either provide a resolution-convergence test of the dipole moment at fixed r0(psi) or soften the quantitative claim to a demonstration of onset.","headline":"Useful methods paper with a genuinely new current-free BC; the 37% burial headline is real but resolution-bound, so treat it as evidence of onset, not a converged number.","tokens_in":37707,"tokens_out":3346,"would_cite":true,"duration_ms":38719,"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":"Mound spreading buries neutron star dipole fields by up to 37 percent","keywords":["neutron stars","accretion mounds","magnetic field burial","Grad-Shafranov equation","magnetostatic equilibria","current-free boundary condition","multipolar magnetic fields","mass ellipticity"],"falsifier":"Run the same GS solver while conserving $dM/d\\psi$ between solutions, or perform a time-dependent axisymmetric MHD simulation of accretion onto a $10^9$ G neutron star with mass loaded up to $r_t = 0.6 R_A$, tracking the mound for hours: if the mound does not spread to the equator and the normalized dipole moment does not decline to about $0.627$, the burial result is an artifact of the assumed profile function.","tokens_in":36527,"feed_emoji":"🧲","tokens_out":10397,"duration_ms":101332,"temperature":0.7,"pith_summary":"The paper claims that magnetically confined accretion mounds on neutron stars do not stay where they are placed: as they grow, ring-shaped mounds spread latitudinally toward the equator, and this spreading is what buries the star's dipole field. Solving the Grad-Shafranov magnetostatic equations at high resolution with a new current-free outer boundary condition, the authors find that a $10^9$ G neutron star with a mound loaded near the Alfvén truncation radius reduces its normalized dipole moment at the outer radius to $0.627$, a 37 percent reduction, with matter reaching the equator. Previous solutions stopped at lower masses because of closed magnetic loops; here the loops are shown to be numerical artifacts of low resolution or of the filled mound profile. If correct, the result means low-field neutron stars can bury their fields within hours of accretion, supporting the accretion-reprocessing scenario for millisecond pulsars.","feed_headline":"Mound spreading buries neutron star dipole fields by up to 37 percent","feed_subtitle":"Mounds slide equatorward as they grow, shrinking the dipole and supporting rapid field burial in low-field pulsars.","key_machinery":"The machinery is the axisymmetric, zero-toroidal-field Grad-Shafranov equation $\\Delta_2 \\psi = K(\\psi,r,\\theta)$ with $K = -4\\pi r^2 \\sin^2\\theta\\, \\rho g\\, dr_0/d\\psi$, where $\\psi$ is the poloidal flux function and $r_0(\\psi)$ is a prescribed profile giving the height of each flux surface; for the ring-shaped mound, $r_0(\\psi)$ is the inverted parabola of Eq. 16, which loads mass only up to the truncation angle $\\theta_t$. The new ingredient is a multipolar current-free boundary condition (Eq. 15) that updates the outer radial boundary at every SOR iteration by decomposing the interior solution into associated Legendre multipoles ($\\ell \\le 33$) and matching to a source-free vacuum exterior, so the dipole and higher moments are free to evolve rather than being pinned to their surface values. This boundary condition is validated by showing that solutions become independent of the outer radius choice. The stretched radial coordinate ($y = \\log((r-aR_*)/(R_*(1-a)))$ with $a=0.999$) and grid resolutions up to $5000^2$–$12000^2$ are what let the solver follow the large gradients at the mound edge and the subsequent equatorward spread.","core_discovery":"The paper's central claim is that magnetically confined accretion mounds, when computed at sufficient resolution with a boundary condition that lets multipoles evolve freely, do not hit a hard mass ceiling: instead the mound spreads equatorward as its mass grows, and this spreading is the physical mechanism of field burial. For a $10^9$ G neutron star with a mound loaded out to $r_t = 0.6 R_A$, the solution spreads all the way to the equator and the normalized dipole moment at the outer radius falls to $0.627$, a 37 percent reduction, with an octupole contribution rising to $0.12$. The paper argues this is the onset of field burial that earlier semi-analytic work anticipated, and that the closed magnetic loops which stopped previous solutions near $10^{-7}$–$10^{-8}$ M$\\odot$ were numerical artifacts: they disappear at high resolution or with the ring-shaped profile, and the true limit is a resolution-dependent Numerical Maximum Mass beyond which solutions become non-unique. It also shows that the same physics applies when the mound sits on a pre-existing ocean, where sinking reduces ellipticity and dipole moment, and that quadru-dipolar surface fields make the two polar mounds asymmetric, with the asymmetry growing until matter funnels through the quadrupolar field.","pith_inferences":["The 37 percent burial figure is tied to the assumed $r_0(\\psi)$; a self-consistent model that conserves $dM/d\\psi$ between solutions could reduce or amplify the spreading, so the number should be read as a proof of mechanism rather than a precise prediction.","If the current-free boundary result transfers to full time-dependent MHD, the same equator-ward spreading might operate in accreting millisecond pulsars on timescales of hours, making magnetic field evolution observable in a single outburst rather than over Gyr.","The octupole enhancement (0.12 of the dipole) in the $10^9$ G burial case suggests that the surface field complexity inferred from X-ray pulse-profile modeling could be a direct signature of ongoing burial, not just an intrinsic multipole.","The quadru-dipolar funneling limit predicts a qualitative change in accretion geometry, hot spots near the equator rather than the pole, which could be searched for in X-ray light curves of accreting pulsars with known field geometries."],"forward_implications":["For a $10^9$ G neutron star with $\\zeta=0.6$, a mound accreted in about 9 hours spreads to the equator and cuts the normalized dipole moment to $0.627$, implying short-term field burial on outburst timescales.","The resolution dependence of the Numerical Maximum Mass means published caps such as $10^{-7}$ M$\\odot$ for a relativistic degenerate EOS are not physical limits; denser or better-resolved grids should push them higher.","Mass ellipticity grows with accreted mass, reaching $\\sim 3.75\\times10^{-10}$ at $10^{12}$ G, and turns over once a mound spreads to the equator, a concrete signature for continuous gravitational wave searches.","Mounds that sink into a pre-existing ocean retain most of their field-burial effect: with a helium composition and the crystallization-depth inner boundary, dipole reduction is slightly weaker and ellipticity is lower than for a pure hard-crust mound.","Quadru-dipolar surface fields create asymmetric polar mounds with different heights and masses at the two poles, and beyond a field-dependent quadrupole fraction ($f_{qd} \\gtrsim 1$ at $10^{10}$ G, $\\gtrsim 5$ at $10^{12}$ G) accretion funnels directly through the quadrupolar field, outside the model's validity."],"supporting_citations":[{"why":"Supplies the Grad-Shafranov formalism, the Green's function solution, and the isothermal EOS baseline that this work extends.","marker":"Payne & Melatos 2004"},{"why":"Introduced the ring-shaped (hollow) mound profile and the fixed dipole boundary; the mass limits this paper revisits.","marker":"Mukherjee & Bhattacharya 2012"},{"why":"Provides the ring-shaped profile, Paczynski EOS, and review of limitations (closed loops, boundary conditions) that motivate the new current-free boundary condition.","marker":"Mukherjee 2017"},{"why":"Previous adiabatic-EOS mound solutions with outflow boundary; the mass limits (1e-7...3e-8 M_sun) shown here to be resolution-dependent.","marker":"Priymak et al. 2011"},{"why":"Spherical-coordinate GS solutions with multipole moments and multipolar initial fields; used for comparison of multipole behavior and the 1/cosh profile.","marker":"Fujisawa et al. 2022"},{"why":"Reduced dipole boundary condition and general-relativistic GS solutions; provides context for why the current-free boundary is needed and the GR limitation.","marker":"Rossetto et al. 2023"},{"why":"Time-dependent MHD simulations of mound sinking onto a fluid base; used for comparison of ellipticity reduction due to sinking.","marker":"Wette et al. 2010"},{"why":"MHD simulations giving truncation radius r_t ~ 0.5-1 R_A; motivates the zeta = 0.6-1 parameter range.","marker":"Parfrey et al. 2017"}],"fun_headline_variants":["Neutron star mounds slide equatorward, burying their dipole fields","Spreading accretion mounds shrink neutron star dipole moment by 37%","Field burial on neutron stars explained by mound spreading, not mass caps","High-res simulations reveal how accretion mounds bury neutron star fields","Mound spreading, not mass limit, drives neutron star field burial"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The equatorward spreading and the 37 percent dipole reduction depend on the hand-specified ring-shaped profile $r_0(\\psi)$, which fixes the mass loading per flux surface; the paper explicitly does not impose a fixed $dM/d\\psi$ between different solutions, so a physically self-consistent accretion model could yield a different profile and materially different burial efficiency.","fun_headline_variants_meta":{"raw":{"variants":["Neutron star mounds slide equatorward, burying their dipole fields","Spreading accretion mounds shrink neutron star dipole moment by 37%","Field burial on neutron stars explained by mound spreading, not mass caps","High-res simulations reveal how accretion mounds bury neutron star fields","Mound spreading, not mass limit, drives neutron star field burial"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001226,"raw_usage":{"total_tokens":5080,"prompt_tokens":1027,"completion_tokens":4053,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":643,"completion_tokens_details":{"reasoning_tokens":3963}},"tokens_in":643,"tokens_out":4053,"duration_ms":33230,"temperature":1.0,"reasoning_tokens":3963,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:17:41.048841+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same GS solver while conserving $dM/d\\psi$ between solutions, or perform a time-dependent axisymmetric MHD simulation of accretion onto a $10^9$ G neutron star with mass loaded up to $r_t = 0.6 R_A$, tracking the mound for hours: if the mound does not spread to the equator and the normalized dipole moment does not decline to about $0.627$, the burial result is an artifact of the assumed profile function.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the ring-shaped profile, Paczynski EOS, and review of limitations (closed loops, boundary conditions) that motivate the new current-free boundary condition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Spherical-coordinate GS solutions with multipole moments and multipolar initial fields; used for comparison of multipole behavior and the 1/cosh profile."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reduced dipole boundary condition and general-relativistic GS solutions; provides context for why the current-free boundary is needed and the GR limitation."}],"review_version":1}