Simulation of a protoplanetary disk accretion activity due to a collision with a gas stream
Pith reviewed 2026-05-08 09:54 UTC · model grok-4.3
The pith
A gas stream colliding with a protoplanetary disk produces accretion rates and time evolution that match observations of FU Orionis stars.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Through three-dimensional numerical gas-dynamic simulation of the collision, the paper establishes that the accretion rate onto the star is governed primarily by the orbital inclination of the infalling stream and the initial mass of the stream material; the resulting maximum accretion rates agree with observational estimates for FU Ori stars, and the time-dependent behavior of the accretion closely tracks available long-term light curves of these objects.
What carries the argument
Three-dimensional numerical gas-dynamic simulation that tracks the density and velocity evolution of the disk and stream after their collision to compute the mass accretion rate onto the central star.
If this is right
- Higher initial mass of the infalling stream produces higher peak accretion rates onto the star.
- The orbital inclination of the stream is the dominant factor shaping how the accretion rate changes with time.
- The simulated accretion rates reach maximum values that are consistent with observational estimates for FU Ori type stars.
- The time evolution of the accretion rate in the simulations reproduces the shape of available long-term light curves for these stars.
Where Pith is reading between the lines
- This collision mechanism may operate across a range of young stellar objects that display episodic accretion, not only FU Ori stars.
- Post-collision disk structures predicted by the model could be targeted in future imaging surveys to identify signs of recent stream infall.
- Adding magnetic fields or radiative transfer would likely modify the detailed flow but leave the strong dependence on inclination and mass intact.
- Light-curve fitting with this model could be used to infer the mass and trajectory of past gas streams that triggered observed outbursts.
Load-bearing premise
The simulation assumes that simplified initial density, velocity, and temperature profiles for the disk and stream, together with pure hydrodynamics, capture the main physical processes that control the collision and subsequent accretion.
What would settle it
High-resolution observations of a FU Ori star during an outburst that show either no infalling gas stream with the modeled velocities and densities or an accretion-rate time series that deviates substantially from the simulated curves for the inferred orbital inclination and stream mass.
Figures
read the original abstract
The consequences of a protoplanetary disk collision with a gas stream are being studied using three-dimensional numerical gas-dynamic simulation. The influence of orbital parameters and the stream mass on the accretion activity of the star is examined. It is shown that the orbital inclination and the initial mass of the infalling material are the most influential parameters in determining the accretion rate. The obtained accretion rate dependencies are compared with actual observational data for two FU~Ori type stars. It turns out that not only is the maximum accretion rate consistent with observational estimates, but the behavior of the accretion rate over time is very similar to available long-term light curves.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents three-dimensional numerical gas-dynamic simulations of a protoplanetary disk colliding with an infalling gas stream. It examines the influence of orbital parameters (particularly inclination) and stream mass on the resulting stellar accretion rate, concluding that inclination and initial stream mass are the dominant factors. The simulated accretion-rate time series is reported to match both the peak value and the temporal behavior of long-term light curves for two FU Ori-type stars.
Significance. If the central results hold under more complete physics, the work identifies a plausible external trigger for episodic accretion outbursts and isolates the key controlling parameters. The forward parameter survey combined with direct comparison to observed light-curve shapes provides a concrete, falsifiable link between simulation and data that could guide future multi-wavelength monitoring and more sophisticated modeling.
major comments (3)
- [§2] §2 (Numerical model): The simulations employ an inviscid, ideal-gas, purely hydrodynamic treatment without magnetic fields or radiative transfer. This omits the magnetorotational instability that supplies the dominant angular-momentum transport in real disks; consequently the post-collision relaxation timescale and the duration of elevated accretion cannot be reliably compared to observed light-curve durations.
- [§3] §3 (Results and comparison): The claim that the time evolution of Ṁ(t) is “very similar” to FU Ori light curves equates the raw mass-accretion-rate history directly to luminosity variations. Because luminosity depends on local dissipation, optical depth, and cooling timescales—none of which are computed—the reported similarity rests on an untested assumption and is therefore not yet load-bearing for the observational consistency argument.
- [§2–§3] §2–§3 (Parameter study): The conclusion that orbital inclination and stream mass are the most influential parameters is drawn from a forward survey of idealized initial conditions. Without reported resolution or convergence tests, or quantitative sensitivity measures, it remains unclear whether the reported ranking of parameters is robust against numerical artifacts or changes in the (unspecified) grid resolution.
minor comments (2)
- [Abstract and §1] The abstract and §1 should explicitly state the numerical code, grid type, and boundary conditions used, as these details are required for reproducibility of the 3-D hydro results.
- [Figures] Figure captions in the results section would benefit from explicit labeling of the time axis units and a brief note on how Ṁ is extracted from the simulation volume.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments. We address each major point below, indicating revisions where the manuscript will be updated to clarify limitations and strengthen the presentation.
read point-by-point responses
-
Referee: [§2] §2 (Numerical model): The simulations employ an inviscid, ideal-gas, purely hydrodynamic treatment without magnetic fields or radiative transfer. This omits the magnetorotational instability that supplies the dominant angular-momentum transport in real disks; consequently the post-collision relaxation timescale and the duration of elevated accretion cannot be reliably compared to observed light-curve durations.
Authors: We agree that the purely hydrodynamic, inviscid model omits MRI-driven transport and radiative effects, limiting the reliability of precise post-collision relaxation timescales. The focus of the work is on the dynamical trigger of the accretion burst by the stream collision and the resulting peak rates. In revision we will expand the methods and discussion sections to explicitly state these limitations and qualify the comparison of burst durations as qualitative rather than quantitative. revision: yes
-
Referee: [§3] §3 (Results and comparison): The claim that the time evolution of Ṁ(t) is “very similar” to FU Ori light curves equates the raw mass-accretion-rate history directly to luminosity variations. Because luminosity depends on local dissipation, optical depth, and cooling timescales—none of which are computed—the reported similarity rests on an untested assumption and is therefore not yet load-bearing for the observational consistency argument.
Authors: The comparison assumes that accretion-rate variations are the primary driver of the observed luminosity changes, which is a standard proxy used in the FU Ori literature. We acknowledge that without radiative transfer this remains indirect. We will revise the results and conclusions to describe the similarity strictly in terms of the accretion-rate time series, note the underlying assumption, and present the match as supportive rather than conclusive evidence. revision: yes
-
Referee: [§2–§3] §2–§3 (Parameter study): The conclusion that orbital inclination and stream mass are the most influential parameters is drawn from a forward survey of idealized initial conditions. Without reported resolution or convergence tests, or quantitative sensitivity measures, it remains unclear whether the reported ranking of parameters is robust against numerical artifacts or changes in the (unspecified) grid resolution.
Authors: The survey was performed at a fixed resolution sufficient to capture the large-scale collision dynamics. While explicit convergence tests and quantitative sensitivity metrics were not included, the dominant trends with inclination and stream mass were reproducible across the explored parameter space. We will add a dedicated paragraph on the numerical resolution and grid setup, together with a statement on the robustness of the parameter ranking based on the existing runs. revision: partial
Circularity Check
No circularity: forward hydrodynamical simulation with parameter variation and post-hoc consistency check
full rationale
The paper conducts 3-D gas-dynamic simulations of a protoplanetary disk colliding with an infalling gas stream, varying orbital inclination and stream mass as inputs to produce accretion-rate time series. These outputs are then compared to FU Ori light curves for consistency, without any fitting of simulation parameters to the observational data or any derivation that reduces the reported accretion behavior to the inputs by construction. No self-citations, ansatzes, or uniqueness theorems are invoked as load-bearing steps in the provided text; the central result is a numerical experiment whose predictions are independent of the target observations.
Axiom & Free-Parameter Ledger
free parameters (2)
- orbital inclination
- initial stream mass
axioms (1)
- domain assumption Ideal-gas hydrodynamics without magnetic fields or radiative transfer suffices to capture the collision and accretion dynamics.
Reference graph
Works this paper leans on
-
[1]
299. 462–478. Hartmann Lee, Kenyon Scott J.The FU Orionis Phenomenon // Ann. Rev. Astron. Astro- phys.. I 1996. 34. 207–240. Heigl S., Hoemann E., Burkert A.Protostellar disk accretion in turbulent filaments // Astron. Astrophys.. VI 2024. 686. A246. Hsieh T.-H., Segura-Cox D. M., Pineda J. E., Caselli P., Bouscasse L., Neri R., Lopez- Sepulcre A., Valdiv...
work page 1996
-
[2]
231. Kolotilov E. A., Petrov P. P.Studies of the Fu-Orionis Stars - Part One - is V1515-CYGNI Starting to Fade // Soviet Astronomy Letters. II 1983. 9. 92–94. Kuffmeier Michael, Frimann Søren, Jensen Sigurd S., Haugbølle Troels. Episodic accretion: the interplay of infall and disc instabilities // Mon. Not. Roy. Astron. Soc.. IV 2018. 475,
work page 1983
-
[3]
Lachaume R.The vertical structure of T Tauri accretion discs
2642–2658. Lachaume R.The vertical structure of T Tauri accretion discs. IV. Self-irradiation of the disc in the FU Orionis outburst phase // Astron. Astrophys.. VII 2004. 422. 171–176. Laznevoi S. I., Akimkin V. V., Pavlyuchenkov Ya. N., Il’in V. B., K´ osp´ al ´A., ´Abrah´ am P. Time-dependent response of protoplanetary disk temperature to an FU Ori-typ...
work page 2004
-
[4]
889, 1. 59. Popham Robert, Kenyon Scott, Hartmann Lee, Narayan Ramesh. Spectra and Line Profiles of FU Orionis Objects: Comparisons between Boundary Layer Models and Observations // Astrophys. J.. XII 1996. 473. 422. Shakura N. I., Sunyaev R. A.Black holes in binary systems. Observational appearance. // Astron. Astrophys. I 1973. 24. 337–355. Stamatellos ...
work page 1996
-
[5]
667. A12. Wachmann A.Das bisherige Verhalten von FU Orionis. Mit 7 Textabbildungen // ZAp. I
-
[6]
35. 74. Williams Jonathan P., Cieza Lucas A.Protoplanetary Disks and Their Evolution // Ann. Rev. Astron. Astrophys.. IX 2011. 49, 1. 67–117
work page 2011
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.