Vertical Shear Instability (VSI): How Disk Temperature Gradients Drive Turbulence

What the VSI Is: Vertical Differential Rotation as Fuel

A protoplanetary disk is not a rigid body: its gas orbits at close to the local Keplerian angular velocity, Ω_K = √(GM/R³). But real disks are pressure-supported and radially stratified in temperature, so the balance of gravity, centrifugal force, and pressure gradient forces Ω to depend weakly on height above the midplane as well as radius. This vertical shear, ∂Ω/∂z ≠ 0, is a store of free energy.

For a locally isothermal disk with a radial temperature power law T ∝ R^q, the vertical gradient scales as

• R·(∂Ω/∂z) ≈ q · (z/R) · Ω_K · (h/R), where h is the pressure scale height,

so the shear is a small fraction — a few percent — of the orbital velocity. A rotating fluid is normally protected against overturning by the Rayleigh criterion (angular momentum rising outward) and by vertical buoyancy. The VSI is the instability that defeats buoyancy and exploits that tiny vertical shear to drive motion.

The Mechanism: Cooling Kills Buoyancy, Shear Does the Rest

The key is double diffusion — inherited from the stellar Goldreich–Schubert–Fricke (GSF) instability. Displace a fluid ring slightly vertically. In a stably stratified disk, buoyancy (set by the Brunt–Väisälä frequency N) pushes it back, stabilizing the flow. But if the disk radiates away temperature perturbations faster than the ring can oscillate, the displaced gas thermally equilibrates with its surroundings and loses its buoyant restoring force.

With buoyancy neutralized, the ring conserves angular momentum instead of entropy. The vertical shear then means that a ring moved along a suitably tilted, radially elongated path finds itself surrounded by gas of different angular momentum — and it keeps going. Growth requires:

• t_cool · Ω ≲ |q| · (h/R) — cooling must be fast, typically t_cool ≲ 0.1 Ω⁻¹;
• modes with short radial wavelength and long vertical extent, which minimize the stabilizing epicyclic response.

The result is columnar, corrugating motions spanning many scale heights that break down into turbulence.

Characteristic Numbers and a Worked Estimate

Consider a solar-mass star (M = 1 M_sun) with a disk at R = 50 au. The Keplerian period is P = 2π/Ω_K ≈ 354 years, so Ω_K ≈ 5.6 × 10⁻¹⁰ s⁻¹. The VSI's fast-cooling condition, t_cool ≲ 0.1/Ω, demands cooling in under about 6 years at that radius — easily met in the optically thin, dust-rich outer disk, but violated near the optically thick midplane inside ~5–10 au.

• Growth timescale: a few to tens of orbits, giving growth rates σ ~ 0.1–0.3 (h/R) Ω.
• Aspect ratio: h/R ≈ 0.05 at 50 au, so vertical velocities reach v_z ~ 0.01–0.1 of the sound speed c_s (tens of m/s).
• Angular-momentum transport: α ≈ 10⁻⁴ to 10⁻³ in simulations — modest, but enough to influence dust.

The turbulence is strongly anisotropic: vertical stirring dominates over radial, which is exactly what makes it observable through how it lofts dust.

How It Is Detected: Dust Settling and CO Kinematics

The VSI leaves two observable fingerprints. First, its vigorous vertical motions loft small dust grains off the midplane, opposing gravitational settling. ALMA images of edge-on disks such as HL Tau and Oph 163131 show razor-thin dust rings, implying vertical diffusion coefficients α_z of only a few × 10⁻⁴ down to ≲ 10⁻⁵ — lower than unimpeded VSI predicts, hinting the instability is damped or dust loading suppresses it.

Second, the corrugation mode produces large-scale meridional circulation — coherent up-and-down gas flows that imprint characteristic wiggles on molecular-line velocity maps. Surveys of CO isotopologue kinematics (e.g., the MAPS program and dedicated VSI-observability studies) search for these ~10–100 m/s corrugation signatures.

• A 2025 A&A analysis of 33 disks found outer regions vertically settled, with α_z typically ≤ radial turbulence.

Together these place real, if still loose, upper limits on VSI strength in nature.

VSI Versus Its Cousins: A Family of Hydro Instabilities

The VSI belongs to a family of hydrodynamic instabilities that became prominent once simulations showed the magnetorotational instability (MRI) is quenched across large "dead zones" where the gas is too poorly ionized for magnetic fields to couple. Each family member exploits a different gradient and thrives in a different cooling regime:

• VSI needs fast cooling (t_cool ≲ 0.1 Ω⁻¹) and feeds on vertical shear — favored in outer, irradiated disks.
• Convective Overstability (COS) and the Subcritical Baroclinic Instability (SBI) use a radial entropy gradient at intermediate cooling and tend to spawn long-lived vortices.
• Zombie Vortex Instability (ZVI) prefers slow, near-adiabatic cooling and vertical buoyancy.

Because cooling time varies with radius, a single disk can host different instabilities at different distances — VSI in the cold outer disk, COS/ZVI further in — making the disk a patchwork of turbulence regimes rather than one uniform state.

Why It Matters and What Remains Open

The VSI matters because turbulence controls planet formation. It sets how efficiently millimeter grains settle into a thin midplane layer, whether the streaming instability can concentrate solids into planetesimals, and how fast gas accretes onto the star. Too much VSI turbulence stirs dust and delays coagulation; too little lets solids settle and clump.

Key open questions:

• Real cooling times. VSI strength hinges on radiative cooling, which depends on dust opacity and grain growth — quantities that VSI turbulence itself alters, creating a feedback loop still being modeled.
• Dust back-reaction. Loading the gas with settled dust can stabilize or self-limit the VSI, possibly explaining the very thin observed rings.
• Interplay with magnetic winds and the MRI at disk surfaces, where recent MHD studies show the picture is more nuanced than pure-hydro models.

Whether the VSI, its cousins, or magnetized winds dominate angular-momentum transport in real disks is one of the central unresolved debates in planet formation.

VSI compared with the other leading purely hydrodynamic disk instabilities. Each exploits a different gradient and a different thermal-relaxation regime.

Instability	Free-energy source	Cooling requirement	Typical α transport
Vertical Shear Instability (VSI)	Vertical shear ∂Ω/∂z from radial temperature gradient	Fast cooling, t_cool ≲ 0.1 Ω⁻¹	10⁻⁴–10⁻³
Convective Overstability (COS)	Radial entropy gradient + epicyclic oscillation	Intermediate, t_cool ~ Ω⁻¹	10⁻⁵–10⁻³ (drives vortices)
Zombie Vortex Instability (ZVI)	Vertical buoyancy (Brunt–Väisälä) + noise	Slow cooling / near-adiabatic	~10⁻⁴
Subcritical Baroclinic Instability (SBI)	Radial entropy gradient (finite amplitude)	Intermediate cooling	Sustains vortices, α ~ 10⁻⁴
Magnetorotational (MRI, for contrast)	Magnetic tension in weak field + shear	Independent of cooling; needs ionization	10⁻³–10⁻² where active