pith. machine review for the scientific record. sign in

arxiv: 2602.20691 · v2 · submitted 2026-02-24 · 🌌 astro-ph.SR · astro-ph.GA

Recognition: 1 theorem link

· Lean Theorem

Multiple protostellar outflows from a single protostar with a misaligned disk

Shingo Hirano , Yuri Aikawa , Masahiro N. Machida
Authors on Pith no claims yet

Pith reviewed 2026-05-15 20:07 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.GA
keywords protostellar outflowsmagnetic misalignmentdisk windsspiral flowsMHD simulationsstar formationcore collapse
0
0 comments X

The pith

Misaligned protostellar disks launch both a disk wind and a spiral flow outflow.

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

This paper examines how misalignment between a protostar's angular momentum and the magnetic field leads to multiple outflow components from a single system. Simulations reveal that a disk wind forms along the disk normal in all cases, but large misalignments also produce a spiral flow parallel to the disk plane. For angles around 60 degrees, the spiral flow can dominate in mass and extent, coexisting with the disk wind. This offers an explanation for observed misaligned secondary outflows without needing multiple protostars. Readers care because it simplifies the interpretation of complex outflow structures in star-forming regions.

Core claim

Three-dimensional nonideal MHD simulations of rotating magnetized cores with misalignment angles from 0 to 90 degrees show that all models produce a magnetocentrifugal disk wind along the disk normal. For misalignments of 60 degrees or larger, a spiral flow component emerges parallel to the disk plane and can become more massive and extended than the disk wind, with intermittent coexistence. The paper concludes that secondary misaligned outflows observed in protostellar systems correspond to this spiral flow, while the primary bipolar outflow is the disk wind from the same misaligned configuration.

What carries the argument

The spiral flow component that develops parallel to the disk plane due to misalignment between core angular momentum and the large-scale magnetic field.

If this is right

  • Mass and size ratios of spiral flow to disk wind increase for misalignment angles of 60 degrees or more.
  • The system can transition from disk wind-dominated to spiral flow-dominated outflow phases.
  • Observed secondary misaligned outflows can be explained as the spiral flow from a single misaligned protostar.
  • Relative lifetimes of the two outflow components vary with the misalignment angle.

Where Pith is reading between the lines

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

  • This could explain outflow multiplicity in regions without clear binary signatures.
  • High-resolution observations might detect the parallel velocity component to confirm the mechanism.
  • It suggests that magnetic field and rotation misalignments are common enough to influence early stellar evolution broadly.

Load-bearing premise

That the spiral flow component in the simulations corresponds directly to observable secondary outflows, and that the initial conditions adequately represent real collapsing cores.

What would settle it

High-resolution observations of a protostellar system showing a secondary outflow component with velocity structure parallel to the disk plane matching the simulated spiral flow characteristics.

Figures

Figures reproduced from arXiv: 2602.20691 by Masahiro N. Machida, Shingo Hirano, Yuri Aikawa.

Figure 1
Figure 1. Figure 1: Three-dimensional structure of the magnetized protostellar system in the misaligned model T60 with an initial angle of θ0 = 60◦ . We show the density contours of the disk (blue), the disk wind (DW; yellow), and the spiralflow component (SF; red). The white lines show the magnetic field lines. 10-5 10-4 10-3 10-2 10-1 100 101 0 1 2 3 4 5 6 7 8 9 10 (a) Star Disk DW SF Mass [M⊙ ] Time, tps [104 yr] 101 102 1… view at source ↗
Figure 2
Figure 2. Figure 2: Time evolution of the outflow properties in the misaligned model T60. Panel (a): Masses of the protostar, circumstellar disk, disk wind (DW), and spiralflow (SF) components as functions of time after protostar formation, tps. Panel (b): Corresponding characteristic extents of the DW and SF, where zDW is the maximum vertical distance along the disk axis and rSF is the maximum radial distance in the disk pla… view at source ↗
Figure 3
Figure 3. Figure 3: Time evolution of the outflow properties for all misalignment models T00–T90. This figure generalizes the single-model evolution shown in [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Mass fractions of the circumstellar components (disk, DW, and SF) as functions of the protostellar mass for models T00–T90 (from top to bottom). while for the SF, we define the maximum radial extent in the disk plane, rSF. We define an outflow component to be “on” when its instantaneous outflow mass exceeds a fiducial thresh￾old, Mout > 10−3 M⊙. Because our analysis focuses on the inner (near-source) outfl… view at source ↗
Figure 5
Figure 5. Figure 5: Dependence of the relative importance of the SF component on the initial misalignment angle θ0. Panel (a): Maximum mass ratio (MSF/MDW)max and maximum size ratio (rSF/zDW)max attained during the evolution for models T00–T90. The horizontal gray line indicates unity. Panel (b): Fractional durations of the simulation during which the DW mass, the SF mass, and both components simultaneously exceed a fiducial … view at source ↗
Figure 6
Figure 6. Figure 6: Simulation visualization of the misaligned protostellar system and the resulting multi-component out￾flows. A protostar is surrounded by a tilted circumstellar disk (blue). From the inner disk, a disk wind (DW; yel￾low lobes) is launched roughly along the local disk axis, fol￾lowing large-scale magnetic-field lines (gray curves). At the same time, twisted magnetic-field lines in and around the disk plane l… view at source ↗
Figure 7
Figure 7. Figure 7: Disk-frame diagrams showing the time evolution of the DW (left) and SF (right) in the representative model T60. Panels (a) and (b) show the plasma beta, panels (c) and (d) the Maxwell stress, panels (e) and (f) the Reynolds stress, and panels (g) and (h) the stress ratio log10(|TR|/|TM|). The left column is plotted against z ′ /h for the DW, and the right column against R ′ /Rdisk for the SF. White regions… view at source ↗
read the original abstract

We investigate how misalignment between the core angular momentum and the large-scale magnetic field affects protostellar outflows, and whether a single protostellar system can drive multiple outflow components. We perform three-dimensional nonideal magnetohydrodynamic simulations of magnetized rotating cores, focusing on the formation of a protostar, a circumstellar disk, and magnetically driven outflows. The initial angle between the core angular-momentum vector and the magnetic field is systematically varied from $0^\circ$ to $90^\circ$. All models launch a classical magnetocentrifugal disk wind (DW) roughly along the local disk normal. For large misalignment, the system also develops a spiralflow (SF) component that propagates parallel to the disk plane. In a representative model with a $60^\circ$ misalignment, the outflow transitions from a DW-dominated to an SF-dominated phase, with the SF becoming more massive and more extended than the DW, and the two components intermittently coexisting. Across the model suite, the maximum mass and size ratios of SF to DW, as well as the relative lifetimes of the two components, increase for misalignment angles $\gtrsim60^\circ$. We propose that secondary, misaligned outflows (or their fossil remnants) observed in some protostellar systems can be interpreted as the SF component, while the main bipolar outflow traces the DW from the same misaligned system.

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 paper reports 3D nonideal MHD simulations of collapsing magnetized cores in which the initial angle between the core angular-momentum vector and the large-scale magnetic field is varied from 0° to 90°. All runs produce a classical magnetocentrifugal disk wind (DW) aligned with the local disk normal; for misalignments ≳60° an additional spiral-flow (SF) component develops parallel to the disk plane. In the 60° case the outflow transitions from DW- to SF-dominated, with the SF becoming more massive and spatially extended; the maximum SF/DW mass and size ratios and the relative lifetimes of the two components increase with misalignment angle. The authors propose that the SF component (or its fossil) can account for observed secondary misaligned outflows while the main bipolar outflow traces the DW from the same single protostar.

Significance. If the SF–DW distinction survives quantitative observational tests, the work supplies a single-source mechanism for the multiple, misaligned outflow components seen in some Class 0/I systems, reducing reliance on binary or multiple-protostar interpretations. The systematic parameter study across misalignment angle is a clear strength; the absence of resolution studies, error analysis, and synthetic observables, however, keeps the result provisional.

major comments (2)
  1. [§4] §4 (results on 60° model): the transition from DW- to SF-dominated outflow and the reported mass/size ratios are presented without resolution-convergence tests or quantitative error estimates, so it is unclear whether the SF dominance is numerically robust or sensitive to the adopted grid or nonideal MHD coefficients.
  2. [Discussion] Discussion section: the central interpretive claim—that the simulated SF component corresponds to observed secondary misaligned outflows—rests on qualitative morphological similarity alone; no synthetic position-velocity diagrams, moment maps, or proper-motion predictions are shown to demonstrate that the SF produces distinguishable kinematic signatures separable from the DW or from projection effects.
minor comments (2)
  1. Figure captions and text should explicitly state the spatial scale and time at which the SF/DW mass and size ratios are measured, and whether the ratios are time-averaged or instantaneous.
  2. Notation: the abbreviation 'SF' is introduced late; early sections should define 'spiral flow (SF)' on first use for consistency.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We address each major comment below, indicating the revisions we will incorporate.

read point-by-point responses

  1. Referee: [§4] §4 (results on 60° model): the transition from DW- to SF-dominated outflow and the reported mass/size ratios are presented without resolution-convergence tests or quantitative error estimates, so it is unclear whether the SF dominance is numerically robust or sensitive to the adopted grid or nonideal MHD coefficients.

    Authors: We agree that explicit resolution-convergence tests and quantitative error estimates are absent from the current manuscript. The grid resolution and nonideal MHD coefficients were selected following standard practices in the literature for comparable protostellar collapse simulations, and the SF component emerges consistently for all models with misalignment ≳60°. In the revised version we will add a dedicated paragraph in §4 that describes the numerical setup, discusses potential sensitivities to grid scale and resistivity coefficients, and provides a qualitative robustness argument based on the systematic parameter study. We will also flag the lack of formal convergence tests as a limitation. revision: yes

  2. Referee: [Discussion] Discussion section: the central interpretive claim—that the simulated SF component corresponds to observed secondary misaligned outflows—rests on qualitative morphological similarity alone; no synthetic position-velocity diagrams, moment maps, or proper-motion predictions are shown to demonstrate that the SF produces distinguishable kinematic signatures separable from the DW or from projection effects.

    Authors: We concur that the current interpretation relies primarily on morphological similarity and that synthetic observables would strengthen the link to observations. Producing full position-velocity diagrams and moment maps lies outside the scope of this initial parameter study. In the revised Discussion we will add a paragraph that qualitatively contrasts the expected kinematics of the SF (predominantly azimuthal and radial motions parallel to the disk plane) with the DW (collimated along the local disk normal), including how projection and inclination effects could separate the two components. We will also outline the requirements for future synthetic observations. revision: partial

Circularity Check

0 steps flagged

No circularity: results follow directly from explicit simulation parameter sweeps

full rationale

The paper performs a suite of 3D nonideal MHD simulations in which the sole varied parameter is the initial misalignment angle (0°–90°) between core angular momentum and the large-scale B-field. The disk wind (DW) and spiral flow (SF) components, their mass/size ratios, and the transition to SF dominance at ≳60° are outputs of those runs, not quantities fitted to data and then relabeled as predictions. The final interpretive sentence (“We propose that secondary, misaligned outflows … can be interpreted as the SF component”) is an after-the-fact mapping suggestion; it does not close any equation or self-define any quantity used in the simulation itself. No self-citations, uniqueness theorems, or ansatzes are invoked to justify the core results. The derivation chain is therefore self-contained computational experiment rather than circular.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on standard nonideal MHD equations and typical initial conditions for rotating magnetized cores; no new free parameters or invented entities are introduced in the abstract.

axioms (1)
  • standard math Nonideal magnetohydrodynamic equations govern the dynamics of the collapsing core
    Invoked as the basis for all simulations described

pith-pipeline@v0.9.0 · 5557 in / 1107 out tokens · 37100 ms · 2026-05-15T20:07:50.572173+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

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

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.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages · 1 internal anchor

  1. [1]

    P., Tollerud, E

    Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167, doi: 10.3847/1538-4357/ac7c74

    work page doi:10.1051/0004-6361/201322068 2013

  2. [2]

    2016, ARA&A, 54, 491, doi: 10.1146/annurev-astro-081915-023341

    Bally, J. 2016, ARA&A, 54, 491, doi: 10.1146/annurev-astro-081915-023341

    work page doi:10.1146/annurev-astro-081915-023341 2016

  3. [3]

    D., & Payne, D

    Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883, doi: 10.1093/mnras/199.4.883

    work page doi:10.1093/mnras/199.4.883 1982

  4. [4]

    L., Davidson, J

    Chapman, N. L., Davidson, J. A., Goldsmith, P. F., et al. 2013, ApJ, 770, 151, doi: 10.1088/0004-637X/770/2/151

    work page doi:10.1088/0004-637x/770/2/151 2013

  5. [5]

    A., et al

    Collette, A., Kluyver, T., Caswell, T. A., et al. 2023, h5py/h5py: 3.8.0-aarch64-wheels, 3.8.0-aarch64-wheels Zenodo, doi: 10.5281/zenodo.594310

    work page doi:10.5281/zenodo.594310 2023

  6. [6]

    2024, AJ, 167, 72, doi: 10.3847/1538-3881/ad152b

    Dutta, S., Lee, C.-F., Johnstone, D., et al. 2024, AJ, 167, 72, doi: 10.3847/1538-3881/ad152b

    work page doi:10.3847/1538-3881/ad152b 2024

  7. [7]

    Early Planet Formation in Embedded Disks (eDisk). XIX. Structures of molecular outflows

    Feeney-Johansson, A., Aikawa, Y., Takakuwa, S., et al. 2025, arXiv e-prints, arXiv:2512.21454, doi: 10.48550/arXiv.2512.21454

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2512.21454 2025

  8. [8]

    R., Millman, K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357, doi: 10.1038/s41586-020-2649-2

    work page doi:10.1038/s41586-020-2649-2 2020

  9. [9]

    2009, A&A, 506, L29, doi: 10.1051/0004-6361/200913008

    Hennebelle, P., & Ciardi, A. 2009, A&A, 506, L29, doi: 10.1051/0004-6361/200913008

    work page doi:10.1051/0004-6361/200913008 2009

  10. [10]

    Hirano, S., Aikawa, Y., & Machida, M. N. 2025, ApJ, 982, 166, doi: 10.3847/1538-4357/adb96f

    work page doi:10.3847/1538-4357/adb96f 2025

  11. [11]

    Hirano, S., Tsukamoto, Y., Basu, S., & Machida, M. N. 2020, ApJ, 898, 118, doi: 10.3847/1538-4357/ab9f9d

    work page doi:10.3847/1538-4357/ab9f9d 2020

  12. [12]

    Hull, C. L. H., Plambeck, R. L., Bolatto, A. D., et al. 2013, ApJ, 768, 159, doi: 10.1088/0004-637X/768/2/159

    work page doi:10.1088/0004-637x/768/2/159 2013

  13. [13]

    Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90, doi: 10.1109/MCSE.2007.55

    work page doi:10.1109/mcse.2007.55 2007

  14. [14]

    2012, A&A, 543, A128, doi: 10.1051/0004-6361/201118730

    Joos, M., Hennebelle, P., & Ciardi, A. 2012, A&A, 543, A128, doi: 10.1051/0004-6361/201118730

    work page doi:10.1051/0004-6361/201118730 2012

  15. [15]

    N., Hirano, S., & Kitta, H

    Machida, M. N., Hirano, S., & Kitta, H. 2020, MNRAS, 491, 2180, doi: 10.1093/mnras/stz3159

    work page doi:10.1093/mnras/stz3159 2020

  16. [16]

    N., & Hosokawa, T

    Machida, M. N., & Hosokawa, T. 2020, MNRAS, 499, 4490, doi: 10.1093/mnras/staa3139

    work page doi:10.1093/mnras/staa3139 2020

  17. [17]

    N., Inutsuka, S.-i., & Matsumoto, T

    Machida, M. N., Inutsuka, S.-i., & Matsumoto, T. 2007, ApJ, 670, 1198, doi: 10.1086/521779

    work page doi:10.1086/521779 2007

  18. [18]

    N., Matsumoto, T., Hanawa, T., & Tomisaka, K

    Machida, M. N., Matsumoto, T., Hanawa, T., & Tomisaka, K. 2006, ApJ, 645, 1227, doi: 10.1086/504423

    work page doi:10.1086/504423 2006

  19. [19]

    2011, ApJ, 728, 47, doi: 10.1088/0004-637X/728/1/47

    Matsumoto, T., & Hanawa, T. 2011, ApJ, 728, 47, doi: 10.1088/0004-637X/728/1/47

    work page doi:10.1088/0004-637x/728/1/47 2011

  20. [20]

    N., & Inutsuka, S.-i

    Matsumoto, T., Machida, M. N., & Inutsuka, S.-i. 2017, ApJ, 839, 69, doi: 10.3847/1538-4357/aa6a1c

    work page doi:10.3847/1538-4357/aa6a1c 2017

  21. [21]

    Maury, A., Hennebelle, P., & Girart, J. M. 2022, Frontiers in Astronomy and Space Sciences, 9, 949223, doi: 10.3389/fspas.2022.949223

    work page doi:10.3389/fspas.2022.949223 2022

  22. [22]

    2024, ApJ, 963, 106, doi: 10.3847/1538-4357/ad1990

    Misugi, Y., Inutsuka, S.-i., Arzoumanian, D., & Tsukamoto, Y. 2024, ApJ, 963, 106, doi: 10.3847/1538-4357/ad1990

    work page doi:10.3847/1538-4357/ad1990 2024

  23. [23]

    2021, ApJ, 910, 11, doi: 10.3847/1538-4357/abddb1

    Okoda, Y., Oya, Y., Francis, L., et al. 2021, ApJ, 910, 11, doi: 10.3847/1538-4357/abddb1

    work page doi:10.3847/1538-4357/abddb1 2021

  24. [24]

    Parker, E. N. 1958, ApJ, 128, 664, doi: 10.1086/146579

    work page doi:10.1086/146579 1958

  25. [25]

    2021, A&A, 648, A45, doi: 10.1051/0004-6361/202038429

    Podio, L., Tabone, B., Codella, C., et al. 2021, A&A, 648, A45, doi: 10.1051/0004-6361/202038429

    work page doi:10.1051/0004-6361/202038429 2021

  26. [26]

    E., & Norman, C

    Pudritz, R. E., & Norman, C. A. 1986, ApJ, 301, 571, doi: 10.1086/163924

    work page doi:10.1086/163924 1986

  27. [27]

    N., et al

    Sai, J., Yen, H.-W., Machida, M. N., et al. 2024, ApJ, 966, 192, doi: 10.3847/1538-4357/ad34b7

    work page doi:10.3847/1538-4357/ad34b7 2024

  28. [28]

    N., et al

    Sato, A., Tokuda, K., Machida, M. N., et al. 2023, ApJ, 958, 102, doi: 10.3847/1538-4357/ad0132

    work page doi:10.3847/1538-4357/ad0132 2023

  29. [29]

    1994, ApJ, 429, 781, doi: 10.1086/174363

    Shu, F., Najita, J., Ostriker, E., et al. 1994, ApJ, 429, 781, doi: 10.1086/174363

    work page doi:10.1086/174363 1994

  30. [30]

    Spruit, H. C. 1996, in NATO Advanced Study Institute (ASI) Series C, Vol. 477, Evolutionary Processes in Binary Stars, ed. R. A. M. J. Wijers, M. B. Davies, & C. A. Tout, 249–286

    work page 1996

  31. [31]

    Spruit, H. C. 2009, in IAU Symposium, Vol. 259, Cosmic Magnetic Fields: From Planets, to Stars and Galaxies, ed. K. G. Strassmeier, A. G. Kosovichev, & J. E. Beckman, 61–74, doi: 10.1017/S1743921309030075

    work page doi:10.1017/s1743921309030075 2009

  32. [32]

    Tomida, K., Okuzumi, S., & Machida, M. N. 2015, ApJ, 801, 117, doi: 10.1088/0004-637X/801/2/117

    work page doi:10.1088/0004-637x/801/2/117 2015

  33. [33]

    N., & Inutsuka, S

    Tsukamoto, Y., Iwasaki, K., Okuzumi, S., Machida, M. N., & Inutsuka, S. 2015, MNRAS, 452, 278, doi: 10.1093/mnras/stv1290

    work page doi:10.1093/mnras/stv1290 2015

  34. [34]

    N., & Inutsuka, S

    Tsukamoto, Y., Okuzumi, S., Iwasaki, K., Machida, M. N., & Inutsuka, S. 2018, ApJ, 868, 22, doi: 10.3847/1538-4357/aae4dc

    work page doi:10.3847/1538-4357/aae4dc 2018

  35. [35]

    2023, in Astronomical Society of the Pacific Conference Series, Vol

    Tsukamoto, Y., Maury, A., Commercon, B., et al. 2023, in Astronomical Society of the Pacific Conference Series, Vol. 534, Protostars and Planets VII, ed. S. Inutsuka, Y. Aikawa, T. Muto, K. Tomida, & M. Tamura, 317, doi: 10.48550/arXiv.2209.13765

    work page doi:10.48550/arxiv.2209.13765 2023

  36. [36]

    J., & Davis, Jr., L

    Weber, E. J., & Davis, Jr., L. 1967, ApJ, 148, 217, doi: 10.1086/149138

    work page doi:10.1086/149138 1967