arxiv: 2602.11818 · v2 · submitted 2026-02-12 · 🌌 astro-ph.EP

Recognition: 2 theorem links

· Lean Theorem

Global magnetohydrodynamic simulations of the inner regions of protoplanetary discs. II. Vertical-net-flux regime

Matthew J. O. Roberts , Henrik N. Latter , Geoffroy Lesur

Authors on Pith no claims yet

Pith reviewed 2026-05-16 02:55 UTC · model grok-4.3

classification 🌌 astro-ph.EP

keywords protoplanetary discsdead zoneactive zonemagnetic fluxMHD simulationsdisc windsvorticesinterface

0 comments

The pith

The active-dead zone interface in protoplanetary discs acts as a one-way barrier that depletes large-scale magnetic flux from the inner active region.

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

Global magnetohydrodynamic simulations model the transition between the well-ionized turbulent inner zone and the poorly ionized outer dead zone in protoplanetary discs threaded by vertical net flux. The interface prevents inward transport of magnetic flux from the dead zone, so flux either drains to the inner boundary or piles up just inside the interface. This setup generates strong variability at the boundary because maintaining a steady vertically integrated current is difficult across evolving magnetic states. A pressure maximum develops at the interface and triggers Rossby-wave vortices, while multiple outflow zones form at all radii instead of a single flux-free transition region. The hot corona keeps the active disc from expanding vertically.

Core claim

In five three-dimensional global simulations with physically motivated resistivity profiles and a cool-disc hot-corona thermodynamic model, the active-dead zone interface functions as a one-way barrier to inward advection of large-scale poloidal magnetic flux. Flux therefore depletes throughout most of the active zone, either advecting to the inner numerical boundary or accumulating immediately inside the interface. Two sources of strong variability appear at the interface because a constant vertically integrated electrical current cannot be sustained across distinct and evolving magnetic-field states. Despite a weak magnetothermal wind in the dead zone, a pressure maximum forms at the dead–

What carries the argument

The active-dead zone interface acting as a one-way magnetic-flux barrier, set by radially varying ionization and resistivity that prevent steady current continuity.

Load-bearing premise

The use of specific Ohmic and ambipolar resistivity profiles together with a simplified cool-disc hot-corona thermodynamic model.

What would settle it

A radial map of magnetic flux density or accretion stress in a real disc showing whether flux is systematically depleted interior to the expected dead-zone edge.

read the original abstract

The inner regions of protoplanetary discs, which encompass the putative habitable zone, are dynamically complex, featuring a relatively well-ionised, turbulent active zone located interior to a poorly ionised 'dead' zone. In this second paper, we investigate a model of the magnetohydrodynamic processes around the interface between these two regions, using five three-dimensional global magnetohydrodynamic simulations of discs threaded by a large-scale poloidal-net-flux magnetic field. We employ physically motivated profiles for Ohmic resistivity and ambipolar diffusion, alongside a simplified thermodynamic model comprising a cool disc and hot corona. Our results show that, first, the interface acts as a one-way barrier to inward transport of large-scale magnetic flux from the dead zone. This leads to magnetic flux depletion throughout most of the active zone, whereby it either advects inwards to the inner numerical boundary or accumulates just inside the interface. Second, two sources of strong variability emerge from the interface due to the difficulty of maintaining a constant, vertically integrated electrical current across distinct and evolving magnetic-field states. Third, despite the weak magnetothermal wind in the dead zone, a pressure maximum forms at the interface, leading to Rossby-wave-induced vortices. Fourth, unlike the model of Iwasaki et al. (2024), there is no 'transition zone' devoid of magnetic flux and magnetic winds. Instead, multiple outflow zones span all disc radii reflecting the radially varying launch conditions, with an inner turbulent wind impinging upon an outer, more laminar one. Fifth, a heated corona prevents the 'puffing up' of poloidal-net-flux, active disc regions.

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.

Desk Editor's Note private letter to a colleague

The simulations find the active-dead zone interface acting as a one-way barrier to large-scale flux, depleting the active zone, but this rests on fixed resistivity profiles that may not generalize.

read the letter

Hi, the central result here is that in these five global 3D MHD runs the interface between the active and dead zones functions as a one-way barrier to inward transport of large-scale poloidal flux. Flux either piles up just inside the interface or gets advected to the inner boundary, leaving most of the active zone depleted. Two variability sources appear at the interface because a constant vertically integrated current is hard to sustain across evolving field states, a pressure maximum forms and spawns Rossby-wave vortices, multiple outflow zones appear at all radii rather than a clean transition zone, and the hot corona keeps the active regions from puffing up. These outcomes differ from Iwasaki et al. (2024) in the absence of a flux-free transition region and in the radially varying launch conditions for the winds. The work uses physically motivated resistivity and ambipolar profiles plus a simple cool-disc/hot-corona thermodynamics, which is a reasonable step forward for global models of the inner disc. The limitation is that everything is shown for one specific set of profiles and only five runs; a different dead-zone width or ambipolar scaling could allow net inward flux motion or remove the depletion. Without reported resolution or convergence tests it is also hard to judge how sensitive the barrier is to numerics. This is the sort of targeted simulation study that people building planet-formation or disc-evolution models would want to read. It is worth sending to referees so they can check the implementation details and ask for profile variations.

Referee Report

2 major / 2 minor

Summary. The paper presents five 3D global MHD simulations of protoplanetary discs with vertical net poloidal flux, using physically motivated Ohmic resistivity and ambipolar diffusion profiles together with a cool-disc/hot-corona thermodynamic model. It claims that the active-dead zone interface functions as a one-way barrier to inward large-scale magnetic flux transport, producing flux depletion in the active zone, strong variability from vertically integrated current maintenance, a pressure maximum with Rossby-wave vortices, multiple radially varying outflow zones without a flux-free transition region, and suppression of poloidal-field puffing by the heated corona. These results are contrasted with the Iwasaki et al. (2024) model.

Significance. If the one-way barrier and associated flux depletion hold under the reported conditions, the work provides valuable numerical evidence for how non-ideal MHD effects and zone interfaces regulate magnetic flux evolution and wind launching in the inner disc. The global 3D approach with realistic resistivity profiles and the identification of multiple outflow zones and interface-driven variability represent concrete advances over local or ideal-MHD models, with direct relevance to accretion and planet-formation processes in the habitable zone.

major comments (2)

[Results (first enumerated claim)] Results, first claim: The one-way barrier to inward large-scale flux transport from the dead zone is demonstrated only for the specific Ohmic resistivity and ambipolar diffusion profiles adopted in the five runs; no additional simulations with varied dead-zone width, ambipolar scaling, or ionization structure are reported, leaving open whether bidirectional flux motion could occur under plausible alternative profiles.
[Numerical methods] Numerical methods: Convergence with respect to grid resolution and explicit diffusivity values is not demonstrated, which is load-bearing for the reported turbulent wind structures, current-sheet variability, and precise location of flux accumulation just inside the interface.

minor comments (2)

[Abstract] Abstract: The statement that 'multiple outflow zones span all disc radii' would benefit from a brief quantitative indication of the radial extent or number of distinct launch regions observed.
[Figures] Figure captions: Captions for the magnetic-field and density slices should explicitly note the time averaging interval and the precise definition of the interface radius used for the one-way barrier analysis.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive report and positive assessment of the significance of our results. We address each major comment below, providing clarifications and indicating where revisions will be made.

read point-by-point responses

Referee: Results, first claim: The one-way barrier to inward large-scale flux transport from the dead zone is demonstrated only for the specific Ohmic resistivity and ambipolar diffusion profiles adopted in the five runs; no additional simulations with varied dead-zone width, ambipolar scaling, or ionization structure are reported, leaving open whether bidirectional flux motion could occur under plausible alternative profiles.

Authors: We agree that the reported one-way barrier is shown for the specific, physically motivated resistivity and ambipolar diffusion profiles described in Section 2, which are based on standard ionization models for the inner disc. The five simulations vary initial magnetic field strength and topology but employ the same profiles, and the barrier behavior is robust across them. We did not perform a broader parameter survey in this work, as the focus was on demonstrating the interface dynamics under representative conditions. We will revise the manuscript to explicitly note this scope limitation and discuss how plausible variations in dead-zone width or ionization could be explored in future studies. revision: partial
Referee: Numerical methods: Convergence with respect to grid resolution and explicit diffusivity values is not demonstrated, which is load-bearing for the reported turbulent wind structures, current-sheet variability, and precise location of flux accumulation just inside the interface.

Authors: We acknowledge that explicit convergence tests with respect to grid resolution and explicit diffusivity are not presented in the current manuscript. The adopted resolution and diffusivity values were chosen based on preliminary lower-resolution runs during code development, in which the key qualitative features (interface barrier, multiple outflow zones, and current-sheet variability) remained present. To address the concern, we will add a dedicated paragraph in the revised numerical methods section summarizing these resolution tests and justifying the final diffusivity choices. revision: partial

Circularity Check

0 steps flagged

No circularity: purely numerical MHD simulation results

full rationale

The paper reports outcomes from five global 3D MHD simulations employing fixed, physically motivated radial profiles for Ohmic resistivity and ambipolar diffusion together with a simplified cool-disc/hot-corona thermodynamics. All central claims—the one-way barrier at the interface, flux depletion in the active zone, variability sources, pressure maximum, and multiple outflow zones—are direct diagnostics extracted from the evolved simulation fields. No analytical derivation chain exists that reduces any prediction to its own inputs by construction, no parameters are fitted to data and then re-predicted, and no load-bearing uniqueness theorem or ansatz is imported via self-citation. The comparison to Iwasaki et al. (2024) is external. The work is therefore self-contained against its stated numerical setup.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claims rest on chosen resistivity and diffusion profiles plus a simplified two-layer thermodynamic model; these are motivated by prior physics but constitute modeling choices rather than derived quantities.

free parameters (2)

Ohmic resistivity profile
Physically motivated functional form selected to represent ionization structure in the disc.
Ambipolar diffusion profile
Physically motivated profile chosen to model non-ideal MHD effects in the dead zone.

axioms (1)

domain assumption Simplified thermodynamic model with cool disc and hot corona
Assumed disc structure used to close the MHD equations.

pith-pipeline@v0.9.0 · 5610 in / 1167 out tokens · 130867 ms · 2026-05-16T02:55:41.125538+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We employ physically motivated profiles for Ohmic resistivity and ambipolar diffusion, alongside a simplified thermodynamic model comprising a cool disc and hot corona.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Rm(R,z) and ΛA(R,z) are prescribed with free parameters Rm0, ΛA0 and interface location RDZI.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.