{"id":"47451148-12d9-4e72-99de-f9a8267a1f03","arxiv_id":"2412.00867","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A new spherical shallow-water code simulates neutron-star spreading layers and produces a rigidly rotating two-armed pattern and hemispheric asymmetry.","lead":"This paper introduces SPLASH, a new hydrodynamics code for compressible flows on a sphere, and uses it to simulate how matter spreads over a neutron star's surface. It finds a stable two-armed tennis ball pattern and lopsided cyclones that could explain quasi-periodic brightness oscillations seen in X-ray binaries.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed observable flux variability is computed from a pressure-projection proxy (Eq. 63) that the authors concede 'do[es] not necessarily affect the radiative energy loss term'; no check shows that real radiative flux follows this proxy.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing issue: Eq. (63) is a pressure-weighted projected area, while the paper's QPO and flux-variability claims require a relation to actual emitted radiation. This is not a manufactured concern; the manuscript itself states that the proxy 'do[es] not necessarily affect the radiative energy loss term', yet the conclusions and abstract treat Lc as though it were the observed flux. The proposed check is concrete and decisive: recompute the periodograms with the physical radiative sink (Eq. 44) and see whether the claimed peaks survive. Because the reader already assigned CONDITIONAL and my concern reinforces that same condition rather than moving it, the verdict should remain unchanged. I agree with the reader's identification and do not see a separate concern that would justify a different adjustment.","tokens_in":17857,"tokens_out":4916,"duration_ms":51293,"concrete_test":"Recompute the four observer time series from the stored snapshots using the actual radiative loss per cell, L_rad(t) = sum_{i: R_obs·R_i>0} [(c g_eff/κ)(1-β)]_i K_i (R_obs·R_i), optionally with a limb-darkening factor, and repeat the same detrending and Lomb-Scargle analysis as in §5.3. If the peaks near 2ν_cycl and ν_cycl/2 survive with comparable power and with false-alarm probabilities below, say, 1%, the proxy concern is resolved; if they shift or vanish, the QPO interpretation must be restricted to the proxy and the conclusion about observable flux variations withdrawn.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The conclusion that the two-armed pattern 'causes flux variations for an observer at a large inclination' depends entirely on Eq. (63), where Lc = sum_i Π_i K_i (R_obs·R_i). This is a pressure-weighted projected area, not a radiative flux. The physical cooling sink used in the simulation is Eq. (44), dE_rad/dt = -(c g_eff/κ)(1-β), whose angular dependence is controlled by g_eff and β, not directly by Π. The authors flag this at the start of §5.3: such a diagnostic 'outline[s] the most dynamically violent processes in the layer, which do not necessarily affect the radiative energy loss term'. Nevertheless, all periodograms in Fig. 13, the QPO discussion, and the v1 abstract's high-Q claim are built on Lc. The two-armed pattern may well exist in the hydrodynamic fields, but the paper provides no evidence that an actual observer's flux tracks Lc. Therefore the central astrophysical claim—observable quasi-periodic variability at ~2ν_cycl and ν_cycl/2—is unsupported unless Lc is validated against a flux-based diagnostic.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops SPLASH, a new two-dimensional hydrodynamics code that solves compressible shallow-water equations on an arbitrary irregular spherical mesh using a multislope second-order MUSCL scheme with an HLLC+ Riemann solver. After verifying conservation and stability on stationary rotating atmospheres, a shock test, and a split-sphere shear test, the code is applied to an accreting neutron-star spreading layer with sources and sinks for mass, angular momentum, and energy. The simulation develops a convective instability and subsequently a global two-armed 'tennis ball' pattern rotating nearly rigidly, accompanied by long-lived cyclones. Light curves are computed for four observer inclinations via a pressure-weighted projected area diagnostic, and time-resolved Lomb-Scargle periodograms show peaks near harmonics of the mean rotation and cyclone pattern frequencies. The authors conclude that the pattern causes observable flux variations at high inclination and that the code can simulate spreading layers in a realistic parameter range.","tokens_in":17980,"tokens_out":2455,"duration_ms":25115,"significance":"If the numerical method holds up, this is a useful contribution: SPLASH is publicly available, conserves mass and energy to roughly machine precision in the tested stationary case, and demonstrates stable behavior at Mach numbers of order 5-10 on spherical unstructured meshes without polar singularities. The spreading-layer simulation also produces a striking, long-lived two-armed pattern that is physically interesting and worthy of further study. However, the central astrophysical claim—that the pattern produces detectable flux variability for a real observer—rests entirely on a diagnostic, Eq. (63), that the authors themselves state need not track the radiative energy loss. The single simulation also lacks a resolution study, and the paper's own Section 5.1 concedes that the run is far from a realistic steady state and omits surface friction. The observational conclusion is therefore not yet established, even though the code-oriented part of the paper is largely sound.","major_comments":[{"comment":"The 'light curve' Lc is a pressure-weighted projected area, not a radiative flux. The paper explicitly notes that this diagnostic 'outline[s] the most dynamically violent processes in the layer, which do not necessarily affect the radiative energy loss term (equation 44).' Nevertheless, all periodograms in Fig. 13, the QPO discussion, and the Section 7 conclusion that the pattern 'causes flux variations for an observer at a large inclination' are derived from Lc. As written, the paper does not demonstrate that a real observer's flux would follow Lc. I request either a validation of Lc against a flux-based diagnostic (e.g., a light curve constructed from the actual cooling term of Eq. 44, or a simple emission model using the local effective temperature) or a substantial toning-down of the observational claims so that they apply only to the pressure diagnostic.","section":"§5.3, Eq. (63)"},{"comment":"The spreading-layer result is based on a single simulation with one cubic mesh of 24,576 faces. Section 4 tests the scheme on stationary and shock problems, but there is no resolution study for the full accretion run, so it is unknown whether the two-armed pattern, its rotation frequency, and the cyclone lifetimes are converged or mesh-dependent. A coarser and a finer run (or at least a second mesh type at similar resolution) are needed to support the claim that the pattern is a robust physical outcome rather than a low-resolution artifact.","section":"§5.1, §5.2, Fig. 10"},{"comment":"The conclusion that this is 'for the first time to simulate the dynamics of an accretion spreading layer in the realistic parameter range appropriate for a real weakly magnetized accreting NS' contradicts the limitations stated in Section 5.1: the simulation covers only about one second, does not reach mass equilibrium, and does not include friction with the NS surface, which the authors identify as 'a crucial constituent for angular momentum conservation and energy release.' The word 'realistic' should either be removed or the claim restricted to the early dynamical phase that is actually simulated.","section":"§5.1 and §7"},{"comment":"The identification of the low-frequency peaks near νcycl/2 and νcycl as 'probably being a result of some inertial mode' is not tested. Section 6 lists several possible modes (Rossby, gravity, pressure) but does not compare the observed peak frequencies to the dispersion relations, except by qualitative statements about 'likely' origins. A mode identification would require, for example, checking the azimuthal wavenumber and latitudinal structure of the oscillations, or comparing with the Rossby-wave relation of Eq. (64) at specific latitudes. As it stands, the mode assignments in the abstract and conclusion are speculative.","section":"§5.3, §6"}],"minor_comments":[{"comment":"The author list contains a typo: 'Pa velAbolmasov' should presumably read 'Pavel Abolmasov.'","section":"Title page"},{"comment":"The GitHub URL is set in quotes as `https://github.com/TURBOLOSE/SPLASH' with a stray quotation mark; this should be cleaned up.","section":"Section 4, footnote 1"},{"comment":"The physical meaning of the relation between H, Π, Σ, and geff would be clearer if the authors stated explicitly that Π is the vertically integrated pressure before introducing Eq. (21).","section":"Section 2.3, Eq. (21)"},{"comment":"The description of the light-curve normalization is incomplete: the units of Lc in Eq. (63) are not stated, and the figure caption for Fig. 13 does not say whether the vertical axis is normalized frequency or physical frequency. Please make the units and normalization explicit.","section":"Section 5.3"},{"comment":"The Gaussian accretion source uses σα as the width but the text later refers to σ = 6°; for clarity, use the same symbol in both places.","section":"Section 2.5, Eq. (35)"},{"comment":"The text says the increment is 'always imaginary, meaning convective stability,' but the sign of N2 and the condition for instability in the Boussinesq analysis are not fully discussed. Please clarify whether the sign convention is chosen so that N2 < 0 corresponds to instability.","section":"Appendix B, Eq. (B18)"}],"recommendation":"major_revision","confidential_remarks":"The paper's numerical-method contribution is defensible: the code is public, the conservation tests are strong, and the mesh-adaptation details are original. My main concern is that the headline astrophysical result—observable flux variability from the two-armed pattern—rests on a pressure-projection proxy that the authors themselves caution may be unrelated to the radiative loss. The authors can likely fix this within the scope of revision by adding a flux-based light curve and clearly separating the robust hydrodynamic findings from the speculative observational interpretation. The absence of a resolution study is also fixable. I would not reject; the method and the pattern dynamics are worth publishing once the claims are matched to the evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid numerical methods paper with an astrophysical application that is more suggestive than conclusive. The code SPLASH is real, public, and tested: conservation to ~1e-14, stable rigid rotation for ten periods, and a Mach-6 shock test on irregular spherical meshes. That part holds up. The multislope MUSCL plus HLLC+ combination on spherical tessellations appears genuinely new for this problem, and the previous spectral approach could not handle shocks. The tennis-ball pattern and the hemispheric asymmetry are new results from the simulation, not fitting artifacts.\n\nThe soft spots are concentrated in the astrophysical section. The main run is a single simulation, about one second long, with no resolution study and no control without the 6-degree tilt. The authors are upfront about the lack of steady state and the missing friction in Section 5.1, which earns credit. The weakest link is the light-curve claim. Lc in Eq. (63) is a pressure-weighted projected area, and Section 5.3 concedes it does not necessarily affect the radiative energy loss term. All periodograms and the QPO discussion use Lc, so the paper shows a pattern in the hydrodynamic fields but does not show that a real observer's flux tracks it. The periodogram peaks are interpreted without significance thresholds or error bars, and the v1 abstract's 'high-quality-factor QPO' overclaims; the v2 abstract is more careful. The 6-degree tilt is purposeful, but since the asymmetry is described as spontaneous, a zero-tilt control would clarify how much of the lopsidedness is source geometry.\n\nThe citation pattern is reasonable. Prior SL work is cited, and the self-citation to Abolmasov et al. 2020 is appropriate because this directly extends that spectral method. No invented entities or hidden circular fits here.\n\nWho gets value from this: anyone working on numerical shallow-water on spheres, and anyone studying spreading layers or boundary layers in LMXBs. It deserves a serious referee. I would send it to review with the expectation that the astrophysical claims get tightened or the paper is reframed as a methods paper with a preliminary application. The code and tests are worth publishing even if the QPO interpretation fades.","headline":"Solid open-source code and honest methods tests; the spreading-layer QPO claim outruns the evidence.","tokens_in":18651,"tokens_out":1449,"would_cite":true,"duration_ms":14349,"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":"A new spherical hydrodynamics code, SPLASH, solves supersonic neutron-star spreading-layer flows at Mach 5-10 and finds a rigidly rotating two-armed pattern that appears as a quasi-periodic oscillation to inclined observers.","keywords":["Computational astronomy","Hydrodynamics","Astrophysical fluid dynamics","Stellar accretion","Neutron stars","Low-mass X-ray binary stars","spreading layer","quasi-periodic oscillations"],"falsifier":"Replace the light-curve proxy with the actual radiative-loss term $\\dot{E}_{\\rm rad} = -c g_{\\rm eff}(1-\\beta)/\\kappa$ from equation (44), recompute the Lomb-Scargle periodograms from the same snapshots, and check whether the peak at twice the pattern rotation frequency survives; if it vanishes, the claimed quasi-periodic oscillation is an artifact of the diagnostic. A second check is to rerun the accretion case on a finer mesh for several mass-renewal times (≳100 s) and see whether the two-armed pattern and its harmonic persist.","tokens_in":17531,"feed_emoji":"🌀","tokens_out":14717,"duration_ms":120857,"temperature":0.7,"pith_summary":"This paper develops SPLASH, a new numerical code that solves the compressible shallow-water equations on a sphere using a multislope second-order MUSCL scheme and an HLLC+ all-speed Riemann solver on an unstructured spherical mesh. The authors' goal is to simulate the accretion spreading layer of a weakly magnetized neutron star in the physically realistic regime, where the flow is supersonic with Mach numbers of order 5-10, a regime that earlier spectral methods could not reach because of Gibbs phenomena at shocks. When they feed a realistic accretion rate onto a 1.5 solar-mass, 10 km neutron star, the equatorial heating triggers a Rayleigh-Taylor instability that mixes the layer in latitude and eventually organizes it into a rigidly rotating two-armed 'tennis-ball' pattern. The paper claims this pattern produces a high-quality quasi-periodic oscillation in the light curve of a high-inclination observer, with a frequency well below the rotation frequency of the matter in the layer, and that the other variability is likely Rossby-wave-like. If correct, the work provides the first numerically grounded, time-dependent model for the variability of accreting neutron stars and a tool for studying any two-dimensional supersonic flow on a sphere.","feed_headline":"Accretion layer on neutron star forms two-armed 'tennis-ball' pattern","feed_subtitle":"A rigidly rotating two-armed pattern in the simulated layer may drive X-ray quasi-periodic oscillations.","key_machinery":"The carrying machinery is the SPLASH finite-volume scheme: a multislope second-order MUSCL reconstruction with a hybrid limiter, paired with an HLLC+ all-speed approximate Riemann solver, implemented on a spherical tessellation whose edges are great-circle arcs, so the grid has no polar singularity and can capture shocks at Mach 5-10. The conserved variables are the surface density $\\Sigma$, the three Cartesian components of angular-momentum surface density $\\mathbf{l}$, and the surface energy density $E$, with an effective two-dimensional adiabatic index $\\Gamma = 2 - 1/\\gamma$. On the physics side, the argument is carried by the latitudinal entropy gradient produced by accretion heating near the equator: Appendix B derives the Brunt-Väisälä-type criterion for rigid-body rotation, showing that a constant-entropy state is neutral while a constant-density state is convectively unstable, and the paper relates the nonlinear outcome to Rayleigh-Taylor growth with effective gravity $g_{\\theta,\\rm eff} = \\Omega^2 R \\sin\\theta \\cos\\theta$. The observable signal is defined by the light-curve proxy $L_c = \\sum_i \\Pi_i K_i (\\mathbf{R}_{\\rm obs}\\cdot \\mathbf{R}_i)$, a pressure-weighted projected area computed for four observer inclinations.","core_discovery":"The paper makes a two-part claim. First, SPLASH, built on a multislope second-order MUSCL reconstruction and an HLLC+ all-speed Riemann solver on an unstructured spherical tessellation, solves the compressible shallow-water equations with second-order accuracy for Mach numbers up to at least 5-10, as demonstrated by stationary-atmosphere, shock, and split-sphere tests. Second, when applied to constant accretion onto a spherical neutron star, the code shows that equatorial heating triggers a Rayleigh-Taylor instability that mixes the layer in latitude and saturates into a rigidly rotating two-armed 'tennis-ball' pattern. The pattern appears in the computed light curves, defined as the pressure-weighted projected area $L_c = \\sum_i \\Pi_i K_i (\\mathbf{R}_{\\rm obs}\\cdot \\mathbf{R}_i)$, as a high-quality quasi-periodic oscillation for high-inclination observers, at a frequency well below the local matter rotation frequency, with additional variability attributed to Rossby modes. The authors conclude that this is the first simulation of a spreading layer in the realistic parameter range appropriate for a weakly magnetized accreting neutron star.","pith_inferences":["The QPO claim should be tested by replacing the pressure-weighted projected-area diagnostic with a proper radiative-transfer calculation of escaping flux; if the $2\\nu_{\\rm cycl}$ peak disappears, the mode is a property of the diagnostic, not of the star.","The topological selection of azimuthal number $m=2$ suggests a testable family of predictions: narrower or tilted accretion bands, or different initial rotation profiles, should excite $m=1$ or $m=3$ patterns with correspondingly different harmonic structure.","If the pattern survives into a steady state, the ratio of the QPO frequency to the neutron star spin frequency becomes a probe of the layer's rotation profile and of the observer's inclination, which could be compared directly with LMXB timing observations."],"forward_implications":["SPLASH can now be used for any two-dimensional hydrodynamic problem on a sphere with Mach numbers up to at least 5-10, including shock-dominated flows, without the polar singularities or Gibbs oscillations that limited earlier spectral approaches.","If the 'tennis-ball' pattern is generic, accreting neutron stars should show a high-quality quasi-periodic oscillation near twice the pattern rotation frequency to high-inclination observers, at a frequency well below the local Keplerian rate.","The predicted hemisphere asymmetry, about 20% difference in mean rotation rates by the end of the run, implies that variability and light-curve properties may depend on which hemisphere faces the observer.","The same simulation framework, with friction and mass-sink terms balanced, should be able to reach a steady state and test whether the pattern and its oscillation modes survive beyond the initial transient accretion phase."],"supporting_citations":[{"why":"Introduces the spreading-layer concept, the vertically integrated equations, and the long angular-momentum-loss timescale that sets the surface-density scale used in the initial conditions.","marker":"Inogamov & Sunyaev 1999"},{"why":"Gives the preceding spectral-method treatment that fails at supersonic Mach numbers, supplies the physical parameter estimates and radiation-pressure equation of state, and identifies the inertial-oscillation mode discussed as the kHz-QPO candidate.","marker":"Abolmasov et al. 2020"},{"why":"Provides the multislope second-order MUSCL scheme with hybrid limiter for general unstructured meshes that SPLASH adapts to the spherical tessellation.","marker":"Touze et al. 2015"},{"why":"Supplies the HLLC+ all-speed Riemann solver whose low- and high-Mach-number fixes make the scheme accurate for Mach 5-10 flows and near-zero polar velocities.","marker":"Chen et al. 2020"},{"why":"Develops the steady-state analytic spreading-layer models that the numerical simulation is designed to connect to and eventually replace.","marker":"Suleimanov & Poutanen 2006"},{"why":"Provides the analytical stability and oscillation-mode analysis of spreading layers that frames the identification of the simulated variability modes.","marker":"Watts et al. 2003"},{"why":"Contributes the Lomb-Scargle periodogram method used to analyze the unevenly sampled simulated light curves.","marker":"Lomb 1976"},{"why":"Completes the Lomb-Scargle spectral-estimation method used for the dynamic power spectra.","marker":"Scargle 1982"}],"fun_headline_variants":["Two-armed tennis-ball pattern forms in neutron-star spreading layer","Neutron-star spreading layer spins a tennis-ball pattern","Code reveals tennis-ball rotation in neutron-star accretion layer","Tennis-ball pattern spins in neutron-star accretion simulation","Simulation shows neutron-star layer forms a rotating tennis-ball pattern"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the pressure-weighted projected area of the simulated layer, equation (63), behaves like the radiation a real telescope would detect; the paper's timing analysis and QPO interpretation rest entirely on this proxy rather than on the actual radiative loss rate.","fun_headline_variants_meta":{"raw":{"variants":["Two-armed tennis-ball pattern forms in neutron-star spreading layer","Neutron-star spreading layer spins a tennis-ball pattern","Code reveals tennis-ball rotation in neutron-star accretion layer","Tennis-ball pattern spins in neutron-star accretion simulation","Simulation shows neutron-star layer forms a rotating tennis-ball pattern"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001013,"raw_usage":{"total_tokens":4324,"prompt_tokens":1037,"completion_tokens":3287,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":653,"completion_tokens_details":{"reasoning_tokens":3206}},"tokens_in":653,"tokens_out":3287,"duration_ms":22159,"temperature":1.0,"reasoning_tokens":3206,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:54:48.169560+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Replace the light-curve proxy with the actual radiative-loss term $\\dot{E}_{\\rm rad} = -c g_{\\rm eff}(1-\\beta)/\\kappa$ from equation (44), recompute the Lomb-Scargle periodograms from the same snapshots, and check whether the peak at twice the pattern rotation frequency survives; if it vanishes, the claimed quasi-periodic oscillation is an artifact of the diagnostic. A second check is to rerun the accretion case on a finer mesh for several mass-renewal times (≳100 s) and see whether the two-armed pattern and its harmonic persist.","supporting_citations":[{"cited_title":"2020, A&A, 638, A142, doi: 10.1051/0004-6361/201936958","cited_arxiv_id":null,"evidence_quote":"Gives the preceding spectral-method treatment that fails at supersonic Mach numbers, supplies the physical parameter estimates and radiation-pressure equation of state, and identifies the inertial-oscillation mode discussed as the kHz-QPO candidate."},{"cited_title":"L., Murrone, A., & Guillard, H","cited_arxiv_id":null,"evidence_quote":"Provides the multislope second-order MUSCL scheme with hybrid limiter for general unstructured meshes that SPLASH adapts to the spherical tessellation."},{"cited_title":"2020, SIAM Journal on Scientific Computing, 42, B921, doi: 10.1137/18M119032X","cited_arxiv_id":null,"evidence_quote":"Supplies the HLLC+ all-speed Riemann solver whose low- and high-Mach-number fixes make the scheme accurate for Mach 5-10 flows and near-zero polar velocities."}],"review_version":1}