{"id":"7e014d4f-de40-4855-a38a-56a7bdfae4f3","arxiv_id":"2609.12998","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A fluid-frame pertabative formulation of resistive GRMHD avoids matrix inversions and reduces cost compared to IMEX methods while capturing leading-order resistive effects.","lead":"The paper presents a new formulation for including finite resistivity in general relativistic magnetohydrodynamic simulations of binary neutron star mergers. It reduces computational cost compared to standard implicit methods by using a peturbative expansion around the ideal limit and a fluid-frame approach that avoids matrix inversions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek/deepseek-v4-flash","headline":"Neglect of time derivatives in fluid-frame REGIME correction is untested in the dynamic binary merger scenario, undermining claim of accuracy.","rationale":"The reader's weakest assumption is the same as ours: the quasi-static approximation for neglecting time derivatives. Our analysis reinforces that this assumption is load-bearing and untested in the dynamical merger case. The reader's verdict was CONDITIONAL, which we maintain, but the condition should be made more explicit: validation against a full resistive simulation. The paper provides no evidence that the correction term is accurate under merger dynamics, only that it is efficient and reproduces ideal and simple resistive tests. Therefore, the verdict remains CONDITIONAL, as the central claim requires further validation.","tokens_in":26898,"tokens_out":2550,"duration_ms":23975,"concrete_test":"Perform a binary neutron star merger simulation with the same initial data and resistivity (η=10^-6) using a full IMEX resistive GRMHD code (e.g., those from [45,64,67]) at the same resolution (∆x=156.8 m on finest grid). Compare the post-merger gravitational-wave phase difference, maximum density evolution, and magnetic energy evolution between REGIME and IMEX up to 20 ms post-merger. If the differences between REGIME and IMEX are of the same order as the differences between REGIME and ideal MHD, or if they exceed typical numerical uncertainty, the neglected time derivatives are significant and the central claim of accuracy is unsupported.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that REGIME reproduces resistive GRMHD accuracy in BNS mergers rests on the approximation Eq. (33) which neglects time-derivative terms V^f in the fluid-frame electric field correction (Eqs. 30-33, Section II D 2). The authors justify this by a quasi-static assumption and near-coplanarity of acceleration, magnetic field, and velocity, but provide no quantification of when these conditions hold or fail. The validation tests do not probe this: the current sheet has zero fluid velocity (so the neglected terms vanish identically), and the TOV star is static. In the binary merger, where velocities are large and the merger is highly dynamical, these assumptions are likely violated at least during merger and early post-merger. The authors present a qualitative comparison between ideal and resistive REGIME simulations, but crucially do not compare against a full resistive GRMHD (IMEX) simulation or an alternative that includes the time derivatives. Without such a comparison, it is impossible to determine whether the observed differences are physical or artifacts of the approximation. The paper itself acknowledges this caveat only in passing and does not estimate the magnitude of neglected terms.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek/deepseek-v4-flash","summary":"This paper presents a fluid-frame formulation of the Resistive Extension upGrade for Ideal MagnEtohydrodynamics (REGIME) to approximate resistive GRMHD in the high-conductivity regime. The method adds a leading-order correction term to the vector-potential evolution equation, computed via spatial derivatives of the ideal-MHD Faraday tensor, avoiding implicit time integration. The implementation in IllinoisGRMHD is tested on a 1D current sheet, a magnetised TOV star, and a binary neutron star merger. A comparison between ideal and resistive binary simulations shows differences in the post-merger gravitational-wave phase, remnant compactness, and magnetic-field structure.","tokens_in":27340,"tokens_out":1933,"duration_ms":19881,"significance":"If the approximation captures the dominant resistive effects at the claimed accuracy, REGIME would offer a computationally cheaper alternative to IMEX-based resistive GRMHD for binary neutron star mergers, enabling routine inclusion of resistivity in parameter studies. The method is derived from first principles and avoids phenomenological tuning. However, the validation is incomplete: the central approximation (neglect of time derivatives in the fluid-frame correction, Eq. 33) is untested in the dynamical merger environment where it is most needed, and no comparison against a full resistive (IMEX) simulation is provided for a problem with non-zero fluid velocity.","major_comments":[{"comment":"The central approximation of the paper—neglecting the time-derivative terms V^f in Eq. (31) to obtain Eq. (33)—is justified by a quasi-static assumption and near-coplanarity, but no quantitative estimate of the neglected terms is given for the binary merger scenario. The validation tests do not probe this: the current sheet (Fig. 1) has v^i = 0 so V^f vanishes identically; the TOV star (Fig. 6) is static. The binary simulation (Sect. V) compares REGIME only against ideal MHD, not against a full resistIVe IMEX run or an alternative that retains the time-derivative terms. Without such a comparison, the claim that 'the new formulation reproduces the expected resisTIve behaviour with acuracy comparable to IMEX-based methods' (abstract) is unsubstantiated for the very system the method is designed for. This should be addressed either by a dedicated test with non-zero velocity (e.g., a moving current sheet or a rotating magnetised star) or by directly computing the magnitude of the neglected terms in the binary simulation.","section":"§II.D.2 (Eqs. 30–33)"},{"comment":"The L1-error scaling for the current sheet shows first-order convergence (p ≈ 1.17), which is lower than the second-order reported for IMEX-based implementations (e.g., Ref. [47,63,65]). The authors attribute this to simple interpolation and differencing, but it is not shown whether the model error (the neglect of higher-order terms in the perturbative expansion) contributes to this order reduction. Since the method is intended to reproduce resistive behaviour 'with accuracy comparable to IMEX,' the origin of the first-order rate should be clarified and a path to second-order accuracy should be discussed.","section":"§IV.A, Fig. 3"}],"minor_comments":[{"comment":"The abstract and introduction state that REGIME 'reproduces the expected resistive behaviour with accuracy comparable to IMEX-based methods.' This claim is too strong given that the binary validation compares only to ideal MHD. I suggest rewording to 'demonstrates consistency with expected resistive behaviour in test problems and produces measurable differences in binary merger simulations.'","section":"Abstract and §I"},{"comment":"The resistivity profile in the TOV simulations is density-dependent, but in the binary simulation (Sect.V) a uniform η=10⁻⁶ is used. The motivation for this choice and its relation to the physical conductivity expected in the merger should be discussed. Also, the binary simulation uses a fixed resistivity, while the TOV test shows the profile implementation works; why not also use a density-dependent profile for the merger?","section":"§III.B.3, Eq. (46)"},{"comment":"The visual comparison of magnetic field lines is qualitative. Consider adding a quantitative measure of ordering (e.g., magnetic energy spectra or a correlation length) to support the claim that resistivity promotes large-scale structures.","section":"§V.C.3, Fig. 14"},{"comment":"The paper acknowledges that 'the neglect of the time-derivative terms… may be inaccurate' in the binary case, but this caveat is buried in the discussion of results rather than prominently featured in the abstract or conclusions. Given the centrality of this approximation, the limitations should be stated more plainly in the summary.","section":"§VI"},{"comment":"The decomposition of ∇_b into D_b and a time derivative along n^a uses n^e ∇_e. The sign convention in the last line of Eq. (31) should be double-checked for consistency with the definition of the acceleration a_b.","section":"§II.D.2, Eq. (31)"}],"recommendation":"major_revision"},"author_rebuttal":null,"desk_editor":{"model":"deepseek/deepseek-v4-flash","letter":"This paper does something genuinely useful: it reformulates the REGIME resistive-GRMHD correction in the fluid frame, eliminating the matrix inversions that made the earlier Eulerian version (Wright 2020) clunky. The algebra is clean, the implementation is straightforward, and the cost overhead (~10%) is small enough that adding resistivity to an existing ideal-GRMHD code becomes practical. The current-sheet and TOV tests behave sensibly. The binary merger simulation shows interesting differences in gravitational-wave phase and remnant structure—plausible signs that resistivity matters even in the high-conductivity regime.\n\nThe soft spot is that the central dynamical test remains a proof-of-principle, not a validation. The derivation neglects time-derivative terms in the fluid-frame correction (Eq. 33), justified by a quasi-static argument and near-coplanarity of acceleration, magnetic field, and velocity. But the current-sheet has zero fluid velocity (so those terms vanish identically), and the TOV star is static. The binary merger simulation, where the approximation is supposed to matter, is never compared against a full resistive-GRMHD or an IMEX run. Without that comparison, we cannot tell whether the observed differences are physical or artifacts of the neglected terms. The paper acknowledges the approximation but does not quantify its error. The first-order convergence is also lower than what IMEX methods typically achieve, which weakens confidence in the accuracy of the correction.\n\nThis is a solid methods paper with an honest limitation statement. The fluid-frame approach is genuinely new and makes resistive simulations cheaper to implement. The missing validation in the dynamical regime is the main issue, but it is a fixable one. I would cite this after seeing a stronger dynamical test, but not yet.\n\nWho this is for: groups already running BNS mergers with ideal GRMHD who want a cheap way to estimate resistive effects. They should be aware that the method's accuracy in the merger phase is unquantified.\n\nRecommendation: send to peer review. The paper is worth refereeing, but the referee should ask the authors to either estimate the size of the neglected time-derivative terms in a representative dynamical configuration, or run a comparison with a full resistive (IMEX) simulation for a simpler but still dynamical problem (e.g., a collapsing magnetized star) where the approximation can be tested directly.","headline":"New fluid-frame REGIME formulation is clever and computationally efficient, but the accuracy of its key approximations in dynamical binary mergers remains unverified.","tokens_in":27687,"tokens_out":1835,"duration_ms":23280},"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek/deepseek-v4-flash","headline":"The REGIME formulation approximates resistive GRMHD in binary neutron star mergers with fluid-frame variables, avoiding the computational cost of implicit-explicit methods.","keywords":["resistive GRMHD","binary neutron star mergers","perturbative method","fluid frame","computational efficiency","gravitational waves","magnetic diffusion","REGIME"],"falsifier":"Compare REGIME predictions with full IMEX resistive GRMHD in a system with rapid magnetic evolution, such as a forced reconnection layer or a turbulent Kelvin-Helmholtz unstable region, where time derivatives of the leading-order fields are large. If the phase differences between the two solutions exceed model error, the quasi-static assumption is violated. Alternatively, run a binary merger simulation with a much higher resistivity (e.g., $\\eta=10^{-4}$) where the perturbation expansion should break down, and verify that the discrepancy with IMEX grows as expected.","tokens_in":28501,"feed_emoji":"💥","tokens_out":5750,"duration_ms":50825,"temperature":0.7,"pith_summary":"Binary neutron star mergers require simulations with finite electrical resistivity, but standard implicit-explicit (IMEX) methods are computationally expensive and complex. This paper presents REGIME, a perturbative method that captures leading-order resistive corrections around the ideal-MHD limit by working in the fluid frame, avoiding matrix inversions. The method adds only a modest overhead to ideal-MHD codes and reproduces expected resistive behaviour in current-sheet tests, magnetised neutron stars, and merger simulations. In binary merger simulations with resistivity, the gravitational-wave signal and remnant compactness differ measurably from the ideal-MHD case, showing that even small resistivity matters for post-merger dynamics.","feed_headline":"Fluid-frame technique models resistivity in neutron-star mergers cheaply","feed_subtitle":"New method avoids costly implicit solvers, finds measurable gravitational-wave phase shifts from finite conductivity.","key_machinery":"The central mechanism is a Chapman-Enskog-like perturbative expansion in the fluid frame. The fluid-frame electric field, which vanishes in ideal MHD, is expanded to first order in the resistivity $\\eta$. Using the Maxwell equations and assuming that time derivatives of the leading-order magnetic field can be neglected (a quasi-static approximation justified for binary neutron star remnants), the first-order correction is expressed as a spatial derivative of the ideal Faraday tensor, projected orthogonally to the fluid four-velocity. This correction is then converted to the Eulerian frame through simple algebraic relations involving the fluid Lorentz factor and the fluid-frame electric field. The resulting evolution equation for the vector potential includes an explicit $\\eta$-dependent correction term that captures leading-order resistive effects without stiff source terms.","core_discovery":"The paper asserts that the fluid-frame REGIME formulation provides a computationally efficient approximation to resistive GRMHD in the high-conductivity regime by expanding the electric field in the fluid frame around its ideal-MHD limit, deriving a leading-order correction that requires only spatial derivatives and simple algebraic transformations back to the Eulerian frame. This avoids the matrix inversions necessary in earlier Eulerian formulations and eliminates the need for implicit time integration. The method is validated against analytical current-sheet solutions, magnetised TOV stars, and binary neutron star mergers, demonstrating accuracy comparable to IMEX methods at substantially reduced complexity and cost. A key finding is that finite resistivity in binary mergers produces cumulative phase differences in the gravitational-wave signal of up to 10 radians and alters the compactness and magnetic-field structure of the remnant, indicating that resistive effects can be observable even in the high-conductivity regime.","pith_inferences":["If the quasi-static approximation used in the derivation (neglect of time derivatives of the leading-order fields) is tested in more turbulent regimes, the method's accuracy may degrade, limiting its applicability to very dynamical phases of the merger.","The uniform resistivity $\\eta=10^{-6}$ used in the merger simulation is likely much larger than actual neutron-star conductivities; systematic studies with smaller $\\eta$ are needed to assess whether such large phase differences persist in physically realistic cases.","The method's first-order spatial accuracy in current-sheet tests may limit its use in problems requiring high precision, but higher-order discretisations could mitigate this without changing the core algorithm."],"forward_implications":["Resistive effects can be incorporated into existing ideal-GRMHD codes with minimal modifications, enabling large-scale parameter studies of magnetised mergers at modest computational cost.","The gravitational-wave phase difference of up to 10 radians between resistive and ideal simulations suggests that future gravitational-wave detectors may constrain the resistivity of neutron-star matter from post-merger signals.","The method naturally handles non-uniform and data-dependent resistivity profiles, as the resistivity appears only algebraically in the update, allowing more realistic microphysical models to be tested.","The suppression of small-scale magnetic structure and the more ordered field in the resistive case implies that resistivity influences angular-momentum transport and disc formation, with implications for kilonova and short gamma-ray burst modelling."],"supporting_citations":[{"why":"provides the foundational Chapman-Enskog-type reduction used to derive the perturbative expansion.","marker":"[72]"},{"why":"introduced the REGIME approach in special-relativistic MHD, which the present work extends to general relativity.","marker":"[75]"},{"why":"developed an earlier Eulerian-frame REGIME formulation that required matrix inversions, motivating the simpler fluid-frame version.","marker":"[77]"},{"why":"one of the first resistive GRMHD implementations using IMEX schemes, used as a reference for comparison.","marker":"[63]"},{"why":"another key IMEX-based resistive GRMHD code, used for validation in TOV star tests.","marker":"[64]"},{"why":"the IllinoisGRMHD code into which REGIME is implemented.","marker":"[24]"}],"fun_headline_variants":["Fluid-frame trick sidesteps costly solvers for neutron-star mergers","Resistive merger signals captured without implicit time integration","New REGIME cuts cost, reveals conductivity’s gravity-wave fingerprint","No matrix inversions: faster resistive GRMHD for binary mergers"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The method assumes that time derivatives of the leading-order magnetic field are negligible and that the lapse gradient (normal-observer acceleration) is nearly coplanar with the magnetic field and velocity; if these do not hold, the resistive correction becomes inaccurate.","fun_headline_variants_meta":{"raw":{"variants":["Fluid-frame trick sidesteps costly solvers for neutron-star mergers","Resistive merger signals captured without implicit time integration","New REGIME cuts cost, reveals conductivity’s gravity-wave fingerprint","No matrix inversions: faster resistive GRMHD for binary mergers"]},"model":"deepseek/deepseek-v4-flash","effort":"low","cost_usd":0.000196,"raw_usage":{"cost":0.000196168,"is_byok":false,"cost_details":{"upstream_inference_cost":0.000196168,"upstream_inference_prompt_cost":0.000101808,"upstream_inference_completions_cost":9.436e-05},"total_tokens":1269,"prompt_tokens":932,"completion_tokens":337,"prompt_tokens_details":{"audio_tokens":0,"video_tokens":0,"cached_tokens":256,"cache_write_tokens":0},"completion_tokens_details":{"audio_tokens":0,"image_tokens":0,"reasoning_tokens":263}},"tokens_in":932,"tokens_out":337,"cost_source":"openrouter_usage","duration_ms":4234,"image_pages":[],"temperature":1.0,"reasoning_tokens":263,"transport_provider":"openrouter","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-09-14T02:43:35.837817+00:00","model_set":{"reader":"deepseek/deepseek-v4-flash"},"pipeline_version":"pith-pipeline@v1.4.1-daily-deepseek-or-medium","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":"Compare REGIME predictions with full IMEX resistive GRMHD in a system with rapid magnetic evolution, such as a forced reconnection layer or a turbulent Kelvin-Helmholtz unstable region, where time derivatives of the leading-order fields are large. If the phase differences between the two solutions exceed model error, the quasi-static assumption is violated. Alternatively, run a binary merger simulation with a much higher resistivity (e.g., $\\eta=10^{-4}$) where the perturbation expansion should break down, and verify that the discrepancy with IMEX grows as expected.","supporting_citations":[{"cited_title":"A fully covariant mean-field dynamo closure for numerical 3+1 resistive GRMHD","cited_arxiv_id":"1205.2951","evidence_quote":"provides the foundational Chapman-Enskog-type reduction used to derive the perturbative expansion."},{"cited_title":"Classical and Quantum Gravity , keywords =","cited_arxiv_id":null,"evidence_quote":"introduced the REGIME approach in special-relativistic MHD, which the present work extends to general relativity."},{"cited_title":"Resistive relativistic MHD simulations of astrophysical jets","cited_arxiv_id":"2308.09477","evidence_quote":"one of the first resistive GRMHD implementations using IMEX schemes, used as a reference for comparison."},{"cited_title":"Magnetic dissipation in short gamma-ray burst jets. I. Resistive relativistic MHD evolution in a model environment","cited_arxiv_id":"2407.11581","evidence_quote":"another key IMEX-based resistive GRMHD code, used for validation in TOV star tests."},{"cited_title":"3D magnetised jet break-out from neutron-star binary merger ejecta: afterglow emission from the jet and the ejecta","cited_arxiv_id":"2009.09714","evidence_quote":"the IllinoisGRMHD code into which REGIME is implemented."}],"lane_accuracy":{"status":"not yet measured","lane":{"model":"deepseek/deepseek-v4-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}