🔭 astrophysics

Investigation of gravitational stability of protoplanetary disks based on statistical analysis of their masses

By statistically analyzing a sample of 1,155 protoplanetary disks, the study finds that only a small fraction appear gravitationally unstable, suggesting that observational limitations likely lead to a systematic underestimation of disk masses and that the true prevalence of unstable systems is significantly higher.

Original authors: Sophia A. Drobchik, Sergey A. Khaibrakhmanov

Published 2026-04-14

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Sophia A. Drobchik, Sergey A. Khaibrakhmanov

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant cosmic kitchen. In this kitchen, stars are being born, and swirling around them are giant, spinning pancakes of gas and dust. These are protoplanetary disks, the nurseries where planets like Earth, Jupiter, and Mars are eventually baked.

For a long time, astronomers have wondered: How do these pancakes turn into planets?

There are two main recipes:

The Slow Build: Tiny dust grains stick together like LEGO bricks, slowly building up into pebbles, then rocks, then planets. This is called core accretion.
The Big Collapse: Sometimes, the pancake gets so heavy and unstable that it rips apart on its own, forming giant clumps that instantly become massive planets. This is called gravitational instability.

The paper you shared is a massive statistical investigation into Recipe #2. The authors, Sophia and Sergey, wanted to know: How often does the "Big Collapse" actually happen in real life?

The Detective Work: Counting the Pancakes

To answer this, the authors didn't just look at one or two disks. They went on a cosmic scavenger hunt, gathering data on 1,155 different protoplanetary disks from ten different star-forming regions (like the Orion Nebula and the Taurus cloud). It's like if a chef wanted to know how often a soufflé collapses, they wouldn't just check one kitchen; they'd check 1,000 kitchens across the country.

They used a mathematical tool called the Toomre Parameter (Q). Think of this as a "Stability Scorecard."

High Score (Q > 2): The pancake is happy, stable, and spinning smoothly. It's not going to collapse.
Medium Score (1 ≤ Q ≤ 2): The pancake is wobbly. It's on the edge of falling apart.
Low Score (Q < 1): The pancake is doomed. It's about to rip into pieces and form giant planets.

The Surprise: A Very Quiet Kitchen

When the authors calculated the scores for all 1,155 disks, they found something strange.

Only about 1.2% of the disks had a "Low Score" (Q < 1).
In other words, out of 1,000 cosmic pancakes, only about 12 were actually unstable enough to collapse on their own.

This was a shock because computer simulations (the theoretical recipes) predicted that gravitational instability should be much more common. The models suggested that massive disks should be ripping apart all the time. But in reality, the kitchen looked surprisingly calm.

The "Underestimation" Mystery

Why is there such a big gap between the theory (which says "lots of collapsing!") and the observation (which says "almost nothing collapsing!")?

The authors propose a clever explanation: We are likely blind to the true weight of these pancakes.

Imagine you are trying to weigh a cloud. If you look at it from the outside, it looks fluffy and light. But if you could see inside, you might realize it's actually a dense, heavy storm cloud.

In astronomy, measuring the mass of a disk is incredibly hard because:

It's Opaque: The dust is so thick (like a foggy window) that we can't see the whole thing. We only see the surface, so we think the disk is lighter than it really is.
The Gas is Hiding: Most of the disk is gas (mostly Hydrogen and Carbon Monoxide). But the gas is often "frozen out" onto dust grains or turned into other chemicals, making it invisible to our telescopes.
The Dust is Growing: The dust grains might be bigger than we think, which changes how they glow.

The authors suggest that because of these "foggy windows," we are underestimating the mass of these disks by a factor of 10 or more.

The Analogy:
If you think a pancake weighs 10 grams, you'd say it's too light to collapse. But if you realized it actually weighs 100 grams, you'd realize, "Oh! That's heavy enough to collapse!"

The authors conclude that the "stable" disks we see might actually be heavily unstable if we could just see their true mass. The "Big Collapse" recipe might be happening much more often than we think; we just can't see the ingredients clearly enough to count them.

What About the Shapes We See?

The paper also looked at pictures of these disks taken by the powerful ALMA telescope.

Stable Disks: Look like smooth, perfect rings.
Unstable/Marginally Stable Disks: Look like they have spiral arms (like a galaxy) or giant gaps.

They found that the disks with the lowest "Stability Scores" (Q < 2) were the ones showing these dramatic spiral shapes. This confirms that when the math says a disk should be unstable, the pictures show it is behaving unstably.

The Bottom Line

This paper is a massive reality check for astronomers. It says:

Observation: We only see a tiny fraction of disks collapsing (1.2%).
Theory: We expect many more to collapse.
Conclusion: Our telescopes are probably underestimating how heavy these disks are. If we correct for the "fog," the number of collapsing disks likely jumps up, bringing reality closer to the theory.

It's a reminder that in the universe, just because something looks light and stable from the outside, doesn't mean it isn't heavy and chaotic on the inside. We just need better glasses to see the truth.

1. Problem Statement

The paper addresses a significant discrepancy between theoretical predictions and observational data regarding the gravitational stability of protoplanetary disks (PPDs).

Theoretical Expectation: Numerical simulations suggest that massive protoplanetary disks should frequently undergo gravitational instability (GI), leading to fragmentation, spiral arm formation, and episodic accretion. This phase is predicted to last for a significant portion of the disk's lifetime.
Observational Reality: Observations of substructures (rings, gaps, spirals) are common, but direct statistical evidence of disks meeting the criteria for gravitational instability (specifically the Toomre parameter $Q < 1$ ) is scarce.
Core Issue: Current observational methods for estimating disk mass (primarily via millimeter continuum and CO isotopologues) suffer from systematic uncertainties, including optical depth effects, CO depletion, and uncertain gas-to-dust ratios. These factors likely lead to a severe underestimation of actual disk masses, potentially masking the true prevalence of gravitational instability.

2. Methodology

The authors conducted a large-scale statistical analysis to quantify the stability of PPDs using the Toomre parameter ( $Q$ ).

Sample Construction:
- Compiled a dataset of 1,706 objects from ten major astronomical surveys covering diverse star-forming regions (Orion, Lupus, Ophiuchus, Taurus, Chamaeleon, etc.).
- Applied strict selection criteria: availability of measured dust mass, known distance, and estimated stellar mass.
- Final Sample: After filtering out objects with extremely low surface densities ( $Q > 200$ ) which could not be corrected to unstable states even with massive mass adjustments, the final sample consisted of 1,155 disks.
Theoretical Framework:
- Stability Criterion: Used the Toomre parameter $Q = \frac{c_s \kappa}{\pi G \Sigma}$ $Q = \frac{c _{s} κ}{π G Σ}$ .
  - $Q < 1$ : Gravitational instability (fragmentation).
  - $1 \le Q \le 2$ : Marginal stability.
  - $Q > 2$ : Gravitational stability.
- Disk Model: Assumed a power-law surface density profile $\Sigma(r) \propto r^{-1}$ and a constant temperature $T = 20$ K (typical for outer disk regions).
- Calculation: Derived $Q$ for each disk using the formula:
  $Q \approx 18 \left(\frac{r_0}{100 \text{ au}}\right)^{1/2} \left(\frac{M_{disk}}{0.01 M_\odot}\right)^{-1} \left(\frac{M_*}{1 M_\odot}\right)^{1/2}$
- Mass Distribution: Fitted the disk mass distribution ( $M_{disk} > 10 M_\oplus$ ) to a power law: $dN/dM \propto M^{-\beta}$ , finding $\beta = 1.36 \pm 0.14$ .
Validation:
- Compared calculated $Q$ values with observed morphological features (spirals, rings, gaps) from the DSHARP survey and other high-resolution ALMA data.
- Cross-referenced results with detailed modeling of specific well-studied systems (e.g., IM Lup, HL Tau).

3. Key Contributions

Largest Statistical Sample: Doubled the size of previous major studies (e.g., Manara et al. 2023) by integrating data from ten distinct catalogs, providing a more robust statistical basis for stability analysis.
Quantification of the "Missing" Instability: Systematically quantified the fraction of formally unstable disks and highlighted the gap between this fraction and theoretical expectations.
Morphological Correlation: Established a link between low $Q$ values and specific morphological features (spiral arms), validating $Q$ as a diagnostic tool despite observational limitations.
Systematic Error Analysis: Provided a comprehensive breakdown of why observational masses are likely underestimated (optical depth, CO depletion, gas-to-dust ratio uncertainties) and how this impacts stability conclusions.

4. Key Results

Stability Distribution:
- Formally Unstable ( $Q < 1$ ): Only 1.2% (14 disks) of the sample.
- Marginally Stable ( $1 \le Q \le 2$ ): 1.7% (20 disks).
- Stable ( $Q > 2$ ): The vast majority (97.1%).
Mass Dependence: Unstable disks are intrinsically massive. The median mass of stable disks is $3.4 \times 10^{-3} M_\odot$ , whereas unstable systems have a median mass of $2.6 \times 10^{-1} M_\odot$ .
Mass Ratio Threshold: The transition to instability ( $Q < 1$ ) occurs when the disk-to-star mass ratio ( $M_{disk}/M_*$ ) exceeds approximately 0.132.
Morphological Consistency:
- Disks with $Q \lesssim 1.7$ (e.g., IM Lup, HL Tau) frequently exhibit spiral structures or complex ring systems consistent with GI models.
- High- $Q$ disks can still show rings (likely driven by planets or other mechanisms), but global spiral arms are predominantly associated with low- $Q$ systems.
Power-Law Fit: The mass distribution of the sample follows $dN/dM \propto M^{-1.36}$ , with a high coefficient of determination ( $R^2 = 0.885$ ).

5. Significance and Discussion

The "Mass Underestimation" Hypothesis: The paper concludes that the extremely low observed fraction of unstable disks (1.2%) contradicts theoretical models predicting longer-lived GI phases. The authors argue this is not because GI is rare, but because current observational methods systematically underestimate disk masses by an order of magnitude or more.
- Contributing factors: Optical depth in millimeter continuum, CO freeze-out/depletion, and uncertainty in the gas-to-dust ratio (assumed to be 100).
Implications for Planet Formation: If actual masses are significantly higher, a much larger fraction of PPDs may be undergoing gravitational instability. This suggests GI could be a more dominant mechanism for giant planet formation (especially super-Jupiters) and disk evolution than currently believed.
Future Directions: The authors recommend:
- Using rare CO isotopologues and HD molecules to bypass optical depth and depletion issues.
- Moving from global $Q$ estimates to radial maps of stability.
- Focusing on Class 0/I disks where GI conditions are most favorable.
- Improving resolution to accurately measure gas disk radii (which are likely larger than dust radii).

Conclusion: The study confirms that while only a tiny fraction of observed disks formally meet the criteria for gravitational instability, this is likely an artifact of observational limitations. The statistical methodology provides a robust framework for identifying candidate systems for GI, suggesting that the true prevalence of gravitational instability in the early universe is significantly higher than current data indicates.