Vertical Structure of Protoplanetary Disks in Scattered Light: A large sample analysis

This paper presents a large-scale analysis of 92 protoplanetary disks using a new ellipse-fitting algorithm on VLT/SPHERE scattered-light images to derive vertical height profiles, revealing that while flared geometries are common, a single power-law model fails to describe the entire sample and significant correlations with system properties are only found in extended disks with radii exceeding 150 au.

The Gist

Imagine a protoplanetary disk not as a flat, featureless pancake, but as a giant, cosmic flying saucer made of dust and gas, swirling around a baby star. For a long time, astronomers knew these disks existed, but they were like trying to guess the shape of a cloud by looking at a 2D photograph. We knew they were there, but we didn't know how "tall" they were or if they were flat like a dinner plate or flared out like a trumpet.

This paper is a massive detective story where the authors tried to measure the vertical height of 92 of these cosmic saucers to see if they follow a specific rule.

Here is the breakdown of their adventure, using some everyday analogies:

1. The Mission: Measuring the "Trumpet" Shape

The main goal was to figure out how the height of the dust changes as you move away from the star.

• The Analogy: Imagine a trumpet. Near the mouthpiece (the star), it's narrow. As you move toward the bell (the outer edge), it flares out wide. Astronomers call this "flaring."
• The Problem: We can't see the disk from the side (unless it's perfectly edge-on). We usually see them from above or at an angle. It's like trying to guess the height of a hill by looking at its shadow on the ground.
• The Solution: The team used a clever computer algorithm called SEEF (Structure Extraction and Ellipse Fitting). Think of this like a digital "cookie cutter." They took high-resolution photos of the dust, found the edges of the rings, and fitted an oval (ellipse) around them. By measuring how much the center of that oval was shifted from the star, they could mathematically calculate how "tall" the dust wall was.

2. The Big Discovery: One Rule Doesn't Fit All

The authors hoped that all these disks would follow the same simple mathematical rule (a "power law"), just like how all human heights might roughly follow a bell curve.

• The Result: They found that no single rule works for everyone.
• The Analogy: It's like trying to describe the shape of all vehicles in a parking lot with one sentence. You can't say "all vehicles are tall and thin" because you have sedans, trucks, and motorcycles.
• The Exception: However, they found one specific group that did follow a perfect rule: the Extended Disks.
  • The Analogy: These are the "Super-Trucks" of the disk world. They are huge (over 150 astronomical units wide). When the authors looked only at these giant disks, they lined up perfectly like soldiers, showing a clear, predictable flaring shape. It seems that when a disk gets big enough, it settles into a very orderly, trumpet-like shape.

3. The "Why" Mystery: What Makes Them Flare?

The team tried to figure out why some disks flare more than others. They looked at three main suspects:

1. How heavy is the dust? (Dust Mass)
2. How old is the star? (Stellar Age)
3. How heavy is the star? (Stellar Mass)

• The Result: No clear connection.
• The Analogy: Imagine trying to figure out why some people are tall. You check their diet (dust mass), their age, and their parents' height (stellar mass). But you find that tall people can be young or old, eat anything, and have parents of any height.
• The "Self-Shadowing" Clue: The authors suspect a sneaky culprit called self-shadowing.
  • The Analogy: Imagine a lighthouse. If the tower is too tall and puffy at the bottom, it casts a shadow on the path in front of it. Similarly, if the inner part of the disk is "puffed up," it blocks the star's light from hitting the outer parts. This makes the outer disk colder and flatter.
  • The Theory: Younger disks might be "puffed up" and shadowing themselves, making them look flatter. Older disks might have cleared out their inner puffy parts, letting the light hit the outer edges, making them flare out beautifully (like the Extended Disks).

4. The "Planet Hunters" Application

Finally, the authors used these height measurements to guess the size of hidden planets.

• The Analogy: Imagine a river (the disk) flowing smoothly. If you see a gap in the water, you know a rock (a planet) is sitting there, pushing the water aside. The size of the gap tells you how big the rock is.
• The Calculation: By measuring the width of the gaps in the dust rings and knowing how "tall" the water (dust) is, they calculated the mass of the planets hiding inside.
• The Surprise: The planets they found are likely small to medium-sized (like Jupiter or smaller) and are currently invisible to our telescopes because they are too faint. They are like "ghosts" in the machine, carving out paths we can see but can't yet see the carvers themselves.

The Bottom Line

This paper is a massive census of 92 baby solar systems. It tells us that:

1. There is no "one size fits all" shape for these disks; they are a diverse bunch.
2. The biggest disks are the most orderly and follow a perfect flaring rule.
3. We still don't know exactly what controls the shape (dust, age, or mass don't seem to be the main drivers), but shadowing might be the key.
4. We can now use these measurements to estimate the size of invisible planets hiding in the gaps.

It's a bit like realizing that while all houses have roofs, they aren't all built the same way, and the biggest mansions seem to follow a very specific architectural blueprint that the smaller cottages don't.

Technical Summary

Here is a detailed technical summary of the paper "Vertical Structure of Protoplanetary Disks in Scattered Light: A large sample analysis" by Byrne et al. (2026).

1. Problem Statement

Protoplanetary disks are the birthplaces of planets, and their vertical geometry (specifically the height of the dust scattering surface) is critical for understanding disk evolution, dust properties, and planet formation mechanisms. While high-resolution imaging has revealed complex morphologies (rings, spirals, gaps), previous studies have been limited to small samples or individual systems.

• The Gap: There is a lack of large-scale, systematic analysis regarding the vertical height structure across diverse disk morphologies.
• The Challenge: Determining the vertical height ($h$) from 2D scattered-light images is geometrically complex due to projection effects. Furthermore, it is unclear how disk flaring correlates with system properties (stellar mass, age, dust mass) or if a universal power-law description exists for all disk types.
• The Goal: To develop a robust methodology for measuring vertical height profiles across a large sample (92 disks) to identify trends, test correlations with system properties, and constrain the masses of potential gap-opening planets.

2. Methodology

The authors utilized a dataset of 92 unique near-infrared (NIR) polarimetric images of protoplanetary disks obtained with the VLT/SPHERE instrument (IRDIS subsystem) between 2015 and 2024. The data covers J, H, and K-bands, with a dominance of H-band observations.

The SEEF Algorithm (Structure Extraction and Ellipse Fitting):
The core of the study is a custom algorithm designed to extract vertical height profiles:

1. Pre-processing: Images are rotated, and optional high-pass filtering (unsharp masking) is applied to enhance edges for high Signal-to-Noise (S/N) disks. An interactive masking tool allows users to exclude artifacts or isolate specific rings in multi-ringed systems.
2. Structure Extraction: Two methods are employed to locate the disk structure:
   • Edge Detection: Identifies the outermost signal exceeding a threshold (typically $5\sigma$) along radial cuts. A step-function criterion filters single-pixel noise.
   • Gaussian Fitting: Fits a Gaussian profile to the intensity peak of the ring structure. This targets the centroid rather than the edge, mitigating some effects of the scattering phase function but requiring higher S/N.
3. Ellipse Fitting:
   • Free Fit: Uses the Direct Least Squares (DLS) method (Fitzgibbon et al. 1999) to fit an ellipse to the extracted points.
   • Constrained Fit: Uses the Levenberg-Marquardt algorithm with priors from ALMA observations (inclination and Position Angle) to stabilize fits for noisy or partially visible structures.
4. Height Calculation: The vertical height ($h$) is derived from the geometric offset ($u$) between the fitted ellipse center and the star position along the minor axis:
   $$h = \frac{u}{\sin(i)}$$
   where $i$ is the disk inclination. The aspect ratio is calculated as $h/r$.

3. Key Contributions

• Large-Scale Dataset: The analysis of 92 unique systems represents the largest sample of scattered-light vertical height measurements to date.
• Robust Methodology: The SEEF algorithm provides a standardized, reproducible pipeline for extracting vertical structures from polarimetric images, handling diverse morphologies (compact, extended, multi-ringed, spirals).
• Morphological Categorization: The sample was subdivided into categories (Extended, Compact, Uniform, Spiral, Inclined, Face-on) to test for morphology-dependent trends.
• Wavelength Dependence Analysis: The study explicitly addresses how different filters (H vs. K-band) probe different vertical layers of the disk, using the WISPIT 2 system as a case study.

4. Key Results

A. Vertical Structure and Power-Law Trends

• Full Sample: The vertical height profile of the entire sample cannot be described by a single power-law relation ($h/r \propto r^{\alpha-1}$). The coefficient of determination ($R^2$) for aspect ratio vs. radius is low ($R^2 = 0.158$), indicating significant scatter.
• Extended Disks ($r_{outer} \ge 150$ au): This subgroup is the only category showing a strong correlation with a single power law ($R^2 = 0.732$ for aspect ratio).
  • Fitted relation: $h = 0.0043 (r/1\text{au})^{1.70}$.
  • This suggests extended disks follow a consistent flaring geometry, likely because they have cleared inner cavities (reducing self-shadowing) and are less affected by the factors causing scatter in other morphologies.
• Other Morphologies: Compact, uniform, and spiral disks show no strong correlation with a single power law, suggesting either distinct vertical structures for each morphology or the influence of unidentified factors.

B. System Property Correlations

• Dust Mass: No correlation found between aspect ratio and dust mass (tested in 40–60 au brackets). Uncertainties in dust mass calculations (optical depth assumptions) may obscure trends.
• Stellar Age: No significant correlation found. However, a weak trend suggests older systems have slightly higher aspect ratios. The authors hypothesize this is due to the reduction of self-shadowing as inner cavities clear in older transition disks, allowing outer regions to heat up and flare more.
• Stellar Mass: No correlation found between aspect ratio and stellar mass.

C. Multi-Ringed Systems & WISPIT 2

• Multi-ringed disks generally show flaring consistent with extended disks.
• WISPIT 2 Case Study: This system was observed in both H and K bands.
  • H-band: Shows a strong correlation with the extended disk trend ($R^2 \approx 0.99$).
  • K-band: Shows a weaker correlation and systematically lower aspect ratios.
  • Implication: Longer wavelengths (K-band) probe deeper, optically thicker layers closer to the midplane, while shorter wavelengths (H-band) trace higher, optically thick scattering surfaces.

D. Planet Mass Constraints

• Using the measured aspect ratios and gap widths, the authors estimated the masses of planets opening gaps using the Kanagawa et al. (2016) empirical model.
• Results: Inferred planet masses range from 0.02 to 2 $M_J$ (Jupiter masses), depending on the assumed disk viscosity ($\alpha_{visc}$).
• Detection Limits: These masses are generally below current direct detection limits for embedded planets. The results suggest these planets may either be a population of wide-separation, low-mass planets or are currently undergoing significant inward migration.

5. Significance and Conclusions

• Methodological Advancement: The paper establishes a consistent framework for measuring disk vertical structure, moving beyond case studies to population-level analysis.
• Morphology Matters: The lack of a universal power law for the full sample implies that disk vertical structure is highly dependent on morphology. Extended disks are the only class exhibiting a "clean" flaring geometry, likely due to the clearing of inner cavities.
• Physical Drivers: The scatter in non-extended disks suggests that factors like self-shadowing (where inner rims block light to outer regions) or varying dust settling/turbulence play a dominant role in shaping vertical structure, overriding simple correlations with stellar mass or age.
• Future Implications: The derived height profiles provide crucial constraints for hydrodynamic simulations of planet-disk interactions. The inferred planet masses, while tentative, point to a population of planets that are currently undetectable by direct imaging but are shaping the disk architecture.

In summary, this work demonstrates that while a universal flaring law does not exist for all protoplanetary disks, extended disks follow a predictable power-law trend. The vertical structure of other disk types is likely governed by complex, morphology-specific interactions involving self-shadowing and dust dynamics.