Early-Phase Circumstellar Disk Dynamics

• Early-phase circumstellar disks are massive, rotationally supported gas and dust structures that emerge around nascent protostars during molecular cloud collapse.
• Simulations show that rapid evolution, influenced by magnetic braking and self-gravity, leads to fragmentation and variable accretion rates impacting stellar assembly.
• High-resolution observations reveal substructures like spiral arms and rings, serving as key diagnostics of disk instability and early chemical evolution.

An early-phase circumstellar disk is the rotationally supported, massive structure of gas and dust formed around a newborn protostar during the collapse of a molecular cloud core, typically within the first 10³–10⁵ years of stellar evolution. These disks emerge immediately following the formation of the first hydrostatic core and evolve rapidly due to interplay between infall, angular momentum transport, magnetic fields, and self-gravity, presenting distinctive thermodynamic, kinematic, and structural properties that profoundly impact pathways of stellar and planetary assembly.

1. Formation Sequence and Initial Disk Properties

During the isothermal collapse of a molecular cloud core ($M\sim1~M_{\odot}$, $T_0=10$ K, $n_0\sim10^6$ cm$^{-3}$), gas density increases until the core becomes optically thick, stalling collapse and generating the first hydrostatic core of radius $\sim$10 AU and mass $M_{\rm disk}\sim10^{-2}$–$10^{-1}\,M_{\odot}$. Conservation of angular momentum ensures rotational support, making this first core the direct seed of the nascent circumstellar disk (Machida et al., 2010, Machida et al., 2010). The protostar (second core, $M_*\sim10^{-3}\,M_{\odot}$, $R_*\sim0.01$ AU) forms as the innermost region heats to $T\sim2000$ K and molecular hydrogen dissociates, but remains embedded in the much more massive disk.

Two disk-forming pathways arise depending on initial core rotation $\beta=E_{\rm rot}/|E_{\rm grav}|$:

• For $\beta\lesssim10^{-4}$, a nested-disk configuration briefly appears, with a small centrifugal structure ($r\sim0.1$–$0.5$ AU) merging with the outer first core in $\sim3$–4 yr.
• For $\beta\sim10^{-3}$–$10^{-2}$ (typical of observed cores), a Keplerian disk of $\sim8$ AU forms contemporaneously with or before protostar birth, rapidly growing in mass and extent (Machida et al., 2010).

Consequently, the disk-to-protostar mass ratio in the initial few $10^3$–$10^4$ yr is generically $\mu\equiv M_{\rm disk}/M_*\sim10$–100 (Machida et al., 2010, Ahmad et al., 2024).

2. Self-Gravity, Disk Instability, and Fragmentation

Disk self-gravity dominates early evolution due to high $M_{\rm disk}/M_*$. The Toomre parameter $Q=c_s\,\kappa/(\pi\,G\,\Sigma)$ (with $c_s$ local sound speed, $\kappa$ epicyclic frequency, and $\Sigma$ surface density) quantifies local gravitational stability:

• $Q\lesssim1$ leads to instability and fragmentation.
• Early-phase disks maintain regions with $Q\sim1$, self-regulated by a balance of heating, cooling, and mass infall (Tsukamoto et al., 2011, Kimura et al., 2020).

Three-dimensional hydrodynamic and radiation-hydrodynamic simulations demonstrate prolific fragmentation under $Q\lesssim1.5$ and efficient cooling, giving rise to numerous bound clumps (Tsukamoto et al., 2013). These clumps, initially of mass $\sim$0.01–0.03 $M_{\odot}$, can be modeled as polytropic spheres with index $n\sim3$ ($T_c\lesssim100$ K) rising to $n\gtrsim4$ ($T_c\gtrsim100$ K) as thermal evolution proceeds, with the maximum clump mass before “second collapse” set by $T_c\sim1000$ K and $M_{\max}\sim0.03~M_\odot$.

Most early-formed clumps undergo rapid inward Type I-like migration, merging with the protostar unless scattered to larger radii. Surviving fragments generally accrete to brown dwarf or low-mass star masses ($\gtrsim0.03\,M_{\odot}$), highlighting the challenge in direct formation of planetary-mass companions by early disk instability (Tsukamoto et al., 2013).

3. Angular Momentum Transport: Magnetic Braking and the Role of Non-Ideal MHD

Ideal magnetohydrodynamics predicts catastrophic magnetic braking, removing disk angular momentum on timescales competitive with main envelope accretion. In magnetized cores with mass-to-flux ratio $\mu\lesssim5$, early disk growth is suppressed, restricting disk radii to $r\lesssim10$ AU as long as the envelope mass exceeds disk mass (Machida et al., 2010, Tsukamoto, 2016, Tomida et al., 2015). The characteristic braking timescale

$$t_{\rm brake} \sim \frac{z_d}{v_{A,{\rm env}}}\left(\frac{\rho_d}{\rho_{\rm env}}\right)$$

matches the accretion time in these regimes, enforcing early-phase disk compactness.

Non-ideal MHD effects—Ohmic dissipation, ambipolar diffusion, and the Hall current—become dynamically significant at $n_{\rm H}\gtrsim10^{10}$–$10^{11}$ cm${}^{-3}$. Ambipolar diffusion and Ohmic resistivity decouple the disk from the field, allowing formation of rotationally supported disks of $r\sim1$–10 AU, even before protostar formation (Tomida et al., 2015, Tsukamoto, 2016). Hall effect introduces bimodality: for anti-parallel alignment of magnetic field and rotation axis (para-disks), the Hall term weakens braking and enables $\gtrsim$20 AU disks, while the parallel case (ortho-disks) can keep disks $<1$ AU (Tsukamoto et al., 2015).

Once envelope mass drops below disk mass, braking torques weaken and outer disk expansion is triggered, rapidly growing the disk to $r\sim100$–200 AU by the end of the Class 0/I phase (Machida et al., 2010, Tomida et al., 2016).

4. Thermodynamic and Kinematic Evolution

The early-phase disk transitions from pressure-supported (first core) to centrifugally supported (Keplerian) as infall continues:

• Initial thickness $H/R\gtrsim0.6$–0.9 decreases to $<0.1$ over $10^3$–$10^5$ yr (Machida et al., 2010, Machida et al., 2010).
• Surface density and temperature profiles generally follow $\Sigma(r)\propto r^{-p}$ ($p\sim1$–1.5) and $T(r)\propto r^{-0.3}$–$-0.75$ (Ahmad et al., 2024, Lee et al., 2017).

Disk self-gravity modifies kinematic structure; rotational velocities deviate from pure Keplerian scaling due to pressure support and mass distribution, e.g., $v_\phi\propto r^{-0.3}$ instead of $r^{-1/2}$ (Ahmad et al., 2024). Toomre-$Q$ remains marginally stable ($Q\gtrsim1$), permitting fragmentation only with continued mass loading or rapid cooling.

Radiative transfer in optically thick disks yields distinctive observational features—e.g., the equatorial dark lane seen in the HH 212 system, reflecting cooler, high-$\tau$ midplane flanked by warmer emitting surfaces (Lee et al., 2017).

5. Morphological and Observational Diagnostics

Early-phase disks display complex substructure:

• Grand-design spiral arms are recurrent manifestations of gravitational instability, efficiently transporting angular momentum and reappearing episodically as fresh infall keeps $Q\sim1$ (Tomida et al., 2016).
• High-resolution ALMA imaging in Class I objects reveals multiple rings and central gaps already at $T_{\rm bol}=235$ K, with features consistent with planet–disk interaction models (e.g., $r_{\rm gap}\sim35$ AU inferring $M_p\sim0.1$–1.8 $M_{\rm Jup}$) and/or dust growth fronts (Shoshi et al., 14 Jan 2026, Ohashi et al., 2020). Magnetic-flux-driven interchange instability further modulates local turbulence and dead zone structure, enabling early planetesimal formation (Shoshi et al., 14 Jan 2026).

Molecular line observations identify chemical youth (e.g., SO detection in AB Aur), as warm, turbulent upper layers prevent rapid freeze-out and grain-surface sequestration of volatiles (Fuente et al., 2010).

6. Environmental Sensitivity, Parameter Dependence, and Evolutionary Modes

The birth and evolution of early-phase disks are sensitive to initial conditions:

• High cosmic-ray ionization rates enhance magnetic coupling and braking, leading to systematically smaller and shorter-lived disks, while the dust fraction (metallicity) exerts a secondary influence (Kobayashi et al., 2023).
• Turbulence in the parent core produces filamentary accretion flows and misaligned disk orientations, driving stochastic variability in disk structure and angular momentum (Tsukamoto et al., 2012).

Multidimensional simulations map disk evolution into a taxonomy spanning massive disk, early/late fragmentation, and protostar-dominant modes depending on the cloud's $\alpha=E_{\rm thermal}/|E_{\rm grav}|$ and $\beta$ parameters. Most models traverse the massive-disk regime ($M_{\rm disk}/M_*\gtrsim1$, no fragmentation), but early/late fragmentation or single-star-dominated pathways are also realized (Tsukamoto et al., 2011). Disks with $\beta\gtrsim3\times10^{-3}$ and $0.4\leq\alpha\leq0.8$ dominate the Class 0/I epoch.

Accretion onto the protostar is highly variable, modulated by both envelope infall and internal disk perturbations (e.g., spiral arms, migrating clumps), resulting in $\dot M_*\sim10^{-6}$–$10^{-4}~M_\odot\,{\rm yr}^{-1}$ with stochastic spikes—potentially observable as FU Ori-type bursts (Machida et al., 2010).

---

The early-phase circumstellar disk sets the physical and chemical initial conditions for planet formation and binary assembly, imprints on observational diagnostics (substructure, chemical tracers, variability), and fundamentally constrains subsequent evolutionary pathways for both the central protostar and any gravitationally bound companions (Machida et al., 2010, Tsukamoto et al., 2013, Tsukamoto et al., 2011, Ahmad et al., 2024).