pith. machine review for the scientific record. sign in

arxiv: 2512.05194 · v2 · submitted 2025-12-04 · 🌌 astro-ph.HE · astro-ph.GA· astro-ph.IM· physics.plasm-ph

Recognition: no theorem link

Lagrangian versus Eulerian Methods for Toroidally-Magnetized Isothermal Disks

Yashvardhan Tomar , Philip F. Hopkins
Authors on Pith no claims yet

Pith reviewed 2026-05-17 00:29 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.GAastro-ph.IMphysics.plasm-ph
keywords toroidal magnetic fieldsisothermal disksLagrangian MHDEulerian MHDaccretion disk simulationsmidplane flux lossresolution convergencemeshless methods
0
0 comments X

The pith

Lagrangian methods drive toroidally-magnetized disks toward midplane flux loss at low resolution where Eulerian methods show no change.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper reruns an idealized test problem of toroidally-magnetized isothermal disks in ideal MHD using two Lagrangian methods instead of the original Eulerian static-mesh approach. Lagrangian simulations match the high-resolution Eulerian outcome of significant toroidal flux loss from the midplane. At low resolution where the thermal scale height is unresolved, Lagrangian methods still evolve toward the converged flux-loss state while Eulerian methods remain frozen with no evolution. The difference is traced to Lagrangian codes following gas flows into arbitrarily thin midplane layers. The authors conclude that sustained midplane toroidal fields in recent multi-scale multi-physics Lagrangian simulations cannot be explained by resolution effects alone.

Core claim

Rerunning the G25 test problem with meshless finite-mass and meshless finite-volume methods reproduces the high-resolution Eulerian results of toroidal flux loss. At low resolution Lagrangian methods lose flux and evolve as close as possible to the converged solution, while Eulerian methods exhibit no evolution. This behavior difference stems from the ability of Lagrangian codes to follow flows to an arbitrarily thin midplane layer, analogous to well-studied differences in Jeans fragmentation problems.

What carries the argument

The ability of Lagrangian methods to follow gas flows to an arbitrarily thin midplane layer

If this is right

  • Lagrangian and Eulerian methods agree on the loss of midplane toroidal flux once the thermal scale height is resolved.
  • Lagrangian methods at low resolution still drive the system toward the same converged flux-loss state.
  • Sustained midplane toroidal fields in multi-scale multi-physics simulations must arise from physical differences with the G25 setup rather than insufficient resolution.
  • The convergence behavior difference is analogous to known issues in Jeans fragmentation where Lagrangian methods allow collapse to thinner layers.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Adding non-ideal MHD effects or radiative cooling to the G25 problem could test whether sustained fields appear only when extra physics is present.
  • Adaptive mesh refinement Eulerian codes might display intermediate behavior between fixed-grid Eulerian and Lagrangian results.
  • Varying the initial field strength or disk thickness in the test problem would show how general the resolution-dependent difference is.

Load-bearing premise

The idealized G25 test problem is representative enough of more complex multi-physics simulations to isolate resolution effects as the sole cause of differing toroidal field behavior.

What would settle it

A simulation of the G25 problem at substantially higher resolution or with an independent code that maintains sustained midplane toroidal fields without adding extra physics would falsify the claim.

Figures

Figures reproduced from arXiv: 2512.05194 by Philip F. Hopkins, Yashvardhan Tomar.

Figure 1
Figure 1. Figure 1: — G25 collapsing disk test problem (§ 2). Top: Time evo￾lution in an example with unresolved thermal scale-height owing to ex￾tremely low-𝛽 (ℎ = 10−4 ), with Lagrangian MFM. Around the circular￾ization radius, we plot the Alfvén scale-height 𝑣𝐴, 𝜙/𝑣K and midplane den￾sity scale-height 𝐻/𝑅 in cylinders, with the Cartesian-equivalent resolu￾tion Δ𝑥 C = (Volume/𝑁cells) 1/3 and median Lagrangian resolution nea… view at source ↗
read the original abstract

A number of simulations have seen the emergence of strongly-toroidally-magnetized accretion disks from interstellar medium inflows. Recently, Guo et al. 2025 (G25) studied an idealized test problem of toroidally-magnetized disks in isothermal ideal MHD with an Eulerian static-mesh method, and argued the midplane behavior changes qualitatively (with a significant loss of toroidal magnetic flux) when the the thermal scale-length is resolved ($\Delta x < H_{\rm thermal}$). We rerun the G25 test problem with two Lagrangian methods: meshless finite-mass, and meshless finite-volume. We show that Lagrangian methods reproduce the high-resolution ($\Delta x \ll H_{\rm thermal}$) Eulerian G25 results. At low resolution ($\Delta x \gg H_{\rm thermal}$), behaviors differ: Lagrangian methods still lose flux and evolve 'as close as possible' to the converged solution, while Eulerian methods show no evolution. We argue this difference in convergence behavior is related to the ability of Lagrangian codes to follow flows to an arbitrarily thin midplane layer, analogous to the well-studied difference in Jeans fragmentation problems. This and results from other higher-resolution simulations and different codes suggest that the sustained midplane toroidal fields seen in recent Lagrangian multi-scale, multi-physics simulations cannot be a numerical resolution effect, and some physical difference between those simulations and the G25 test problem explains their different behaviors.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript compares Lagrangian meshless finite-mass (MFM) and meshless finite-volume (MFV) methods against Eulerian static-mesh methods on the G25 idealized test problem of toroidally-magnetized isothermal disks in ideal MHD. It reports that Lagrangian and Eulerian results agree at high resolution (Δx ≪ H_thermal), but diverge at low resolution (Δx ≫ H_thermal): Lagrangian methods lose toroidal flux and evolve toward the converged high-resolution state, while Eulerian methods show no evolution. This difference is attributed to Lagrangian codes' ability to follow flows into arbitrarily thin midplane layers (analogous to Jeans fragmentation). The authors conclude that sustained midplane toroidal fields in recent multi-scale, multi-physics Lagrangian simulations cannot be explained as a resolution artifact and must arise from physical differences with the G25 setup.

Significance. If the numerical comparisons hold, the work clarifies method-dependent convergence behaviors in ideal MHD disk problems and provides a concrete explanation for why Lagrangian codes can maintain thin midplane structures even at modest resolution. This has direct implications for interpreting toroidal-field results from complex galaxy-formation and accretion-disk simulations, and the side-by-side testing of independent codes on a shared benchmark is a useful contribution to the numerical astrophysics literature.

major comments (2)
  1. The central claim that sustained midplane toroidal fields in multi-physics Lagrangian simulations are not a resolution effect rests on the assumption that the idealized isothermal G25 test is representative. Additional physics (cooling, self-gravity, star formation, non-isothermal thermodynamics) can alter the effective scale height, thin-layer dynamics, and numerical flux diffusion; without explicit tests showing that the low-resolution Lagrangian convergence behavior persists when these processes are included, the extrapolation remains suggestive rather than definitive.
  2. Quantitative support for the reported flux loss and convergence is limited. The manuscript would benefit from explicit plots or tables of toroidal magnetic flux evolution versus time (or versus resolution) with error estimates or convergence metrics for both Lagrangian and Eulerian runs at matched resolutions, rather than qualitative descriptions of behavior.
minor comments (2)
  1. Clarify the precise definitions and units of Δx and H_thermal in the methods section, and state the initial conditions (density profile, magnetic field strength, velocity field) explicitly so that the test can be reproduced independently.
  2. The analogy to Jeans fragmentation is useful but would be strengthened by a brief reference to the relevant literature on Lagrangian vs. Eulerian behavior in self-gravitating collapse problems.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their constructive and insightful comments on our manuscript. We address each major comment below, indicating the revisions we will incorporate where appropriate.

read point-by-point responses

  1. Referee: The central claim that sustained midplane toroidal fields in multi-physics Lagrangian simulations are not a resolution effect rests on the assumption that the idealized isothermal G25 test is representative. Additional physics (cooling, self-gravity, star formation, non-isothermal thermodynamics) can alter the effective scale height, thin-layer dynamics, and numerical flux diffusion; without explicit tests showing that the low-resolution Lagrangian convergence behavior persists when these processes are included, the extrapolation remains suggestive rather than definitive.

    Authors: We agree that the G25 test problem is idealized and that additional physical processes could modify the effective scale height, thin-layer formation, and numerical diffusion of flux. Our primary objective is to isolate the numerical convergence difference between Lagrangian and Eulerian methods in a controlled isothermal MHD setup, demonstrating that Lagrangian methods can follow flows into arbitrarily thin midplane layers even at low resolution (analogous to Jeans fragmentation). This explains the observed flux loss at low resolution in Lagrangian runs, which approach the high-resolution converged state, while Eulerian runs remain static. We have revised the manuscript to include an expanded discussion section that explicitly acknowledges the limitations of the idealized test and cites supporting results from other higher-resolution and multi-physics Lagrangian simulations in the literature. We do not claim the G25 test is fully representative of all multi-physics cases but argue that the demonstrated numerical behavior makes a resolution artifact explanation for sustained toroidal fields in those simulations unlikely. revision: partial

  2. Referee: Quantitative support for the reported flux loss and convergence is limited. The manuscript would benefit from explicit plots or tables of toroidal magnetic flux evolution versus time (or versus resolution) with error estimates or convergence metrics for both Lagrangian and Eulerian runs at matched resolutions, rather than qualitative descriptions of behavior.

    Authors: We thank the referee for this suggestion to strengthen the quantitative presentation. In the revised manuscript we have added two new figures: one showing the time evolution of the total toroidal magnetic flux (normalized to the initial value) for MFM, MFV, and Eulerian runs at both high and low resolutions, and a second showing flux versus resolution at fixed late time. These include shaded regions for the standard deviation across an ensemble of runs and a table summarizing the fractional flux loss at t = 10 orbital periods together with a simple convergence metric (L2 difference relative to the high-resolution reference solution). revision: yes

standing simulated objections not resolved
  • Explicit tests incorporating additional physics (cooling, self-gravity, star formation, non-isothermal thermodynamics) to confirm that the low-resolution Lagrangian convergence behavior persists under those conditions.

Circularity Check

0 steps flagged

No significant circularity; results from independent numerical experiments

full rationale

The paper's central claims rest on direct reruns of the G25 isothermal ideal-MHD test problem using meshless finite-mass and finite-volume methods, with explicit comparisons to published Eulerian static-mesh results at varying resolutions. Convergence behaviors, flux loss at low resolution, and the analogy to Jeans fragmentation are derived from these new simulations rather than from any self-referential fitting, parameter definitions, or load-bearing self-citations. The conclusion that sustained toroidal fields in multi-physics runs are not a resolution artifact follows from the observed differences in method behavior on the shared test problem, without reducing any prediction to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard assumptions of ideal MHD and an isothermal equation of state for the test problem, with no new free parameters, axioms beyond domain standards, or invented entities.

axioms (2)
  • domain assumption Ideal MHD approximation is valid for the disk plasma.
    Standard assumption invoked for modeling magnetic fields in astrophysical accretion disks.
  • domain assumption Isothermal equation of state applies to the G25 test problem.
    Directly adopted from the referenced G25 setup being re-run.

pith-pipeline@v0.9.0 · 5570 in / 1367 out tokens · 129600 ms · 2026-05-17T00:29:19.457138+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Thick Disks, Thin Hopes: Suppressed Capture and Merger Rates in AGN

    astro-ph.HE 2026-01 unverdicted novelty 6.0

    AGN disk capture and merger rates scale as (H/R)^{-8} and drop by 10-20 orders of magnitude in thick magnetically supported disks compared to thin thermal disks.

Reference graph

Works this paper leans on

5 extracted references · 5 canonical work pages · cited by 1 Pith paper

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

    work page

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

    work page

  3. [3]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence := #2 '...

    work page

  4. [4]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

    work page

  5. [5]

    write newline

    " write newline "" before.all 'output.state := FUNCTION fin.entry write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION new.sentence output.state after.block = 'skip output.state before.all = 'skip after.sentence 'output.state := if if FUNCTION not #0 #1 if FUNCTION and 'skip pop #0 if FUNCTION or pop #1...

    work page