Toomre Q Disk Stability Calculator

What this calculator does

The Toomre Q parameter is one of the most widely used local stability checks in disk astrophysics. It asks whether a rotating, self-gravitating disk can resist axisymmetric collapse at a given location. This page lets you enter the local surface density, sound speed, and epicyclic frequency, then computes the dimensionless stability parameter $Q$, a characteristic critical wavelength, and an approximate fragment mass scale. The result is a fast first-pass diagnostic that is useful in protoplanetary disk work, galactic dynamics, and classroom examples where you want to connect physical intuition to a numerical estimate.

The physical picture is simple but powerful. Gravity tries to pull matter together. Pressure or random motions resist compression. Rotation and shear also help prevent collapse by changing how perturbations move through the disk. The Toomre criterion combines those competing effects into one compact quantity. When $Q$ is small, self-gravity is strong compared with pressure support and rotational stabilization, so the disk is more vulnerable to local collapse. When $Q$ is large, the disk is harder to compress and is more likely to remain locally stable.

This calculator is intentionally focused on the classic thin-disk, local, axisymmetric version of the problem. That makes it easy to use and easy to interpret, but it also means the output should be treated as a diagnostic rather than a complete physical model. Real disks can be thick, turbulent, magnetized, cooling, multi-component, or globally unstable in ways that a one-number local criterion cannot fully capture. Even so, the Toomre parameter remains a standard starting point because it quickly tells you whether a region deserves closer attention.

How to use the inputs

The form below asks for three quantities. The first is the surface density $\Sigma$, entered in kilograms per square meter. This represents how much mass is spread across the disk surface locally. Larger values strengthen self-gravity and tend to lower $Q$. The second is the sound speed $c_s$, entered in kilometers per second. In gas disks this may be the thermal sound speed, while in some practical applications it can stand in for an effective velocity dispersion. Larger values increase pressure support and tend to raise $Q$. The third is the epicyclic frequency $\kappa$, entered in inverse seconds. This measures the local restoring effect of orbital motion and differential rotation.

When you press the compute button, the script preserves the original calculator behavior. It reads the three fields, converts the sound speed from kilometers per second to meters per second, evaluates the formulas, and writes the result into the output area. The output includes the numerical value of $Q$, the critical wavelength in kilometers, the fragment mass in kilograms, and a simple status label: Unstable, Marginal, or Stable. The status thresholds are practical labels for quick interpretation, not a replacement for detailed modeling.

It is worth pausing on units because unit mistakes are the most common source of confusion. The sound speed field is labeled in kilometers per second, but the equations are evaluated in SI units. The script therefore multiplies the entered sound speed by 1000 before using it. If you accidentally type meters per second into that field, the calculator will treat the number as kilometers per second and overestimate the stabilizing pressure term by a factor of 1000. Surface density should be entered directly in kilograms per square meter, and epicyclic frequency should be entered in inverse seconds. Scientific notation such as 1e-15 is acceptable in the numeric fields.

In a nearly Keplerian disk, such as an idealized disk orbiting a dominant central mass, the epicyclic frequency is often close to the orbital angular frequency, so one may use $\kappa =\Omega$ as a good approximation. In galactic disks, however, $\kappa$ is usually derived from the rotation curve and need not equal $\Omega$. The calculator does not derive $\kappa$ from radius or velocity data, so the quality of the result depends on supplying a physically appropriate value for the system you are studying.

Core formulas

The main quantity computed here is the classic Toomre stability parameter:

In this expression, $c_s$ and $\kappa$ appear in the numerator because pressure support and orbital restoring forces help resist collapse, while $\Sigma$ appears in the denominator because a denser disk has stronger self-gravity. The gravitational constant $G$ sets the scale for how strongly mass attracts mass. The result is dimensionless, which is one reason the parameter is so convenient for comparing very different astrophysical environments.

The local axisymmetric stability analysis is often summarized by the thin-disk dispersion relation:

Here $\omega$ is the perturbation frequency and $k$ is the radial wavenumber. The first term is stabilizing rotation, the second is destabilizing self-gravity, and the third is stabilizing pressure. Instability becomes possible when ${\omega }^2$ drops below zero for some wavelength. The Toomre threshold condenses that condition into the simpler statement that local axisymmetric instability is expected when in the idealized thin-disk limit.

The calculator also reports a characteristic critical wavelength:

and an approximate fragment mass:

These extra outputs are useful because a stability label alone does not tell you the scale of the structures involved. A disk can be unstable in principle, but the preferred unstable wavelength may correspond to very different physical objects depending on the environment. In one setting it might suggest a clump on a planetary scale; in another it might point toward a giant cloud complex or a much larger galactic feature. The fragment mass estimate should be read as an order-of-magnitude guide based on a circular patch of radius $\frac{{\lambda }_{\text{crit}}}{2}$, not as a precise prediction of what nature must produce.

For context, the epicyclic frequency itself is often related to the angular frequency by

That relation is included here because it explains where the input $\kappa$ often comes from in practice, especially in galactic dynamics. The calculator does not evaluate this derivative expression; it simply uses the value you provide.

Worked example

Suppose you enter a surface density of 1000 kg/m², a sound speed of 1 km/s, and an epicyclic frequency of 1×10⁻¹⁵ s⁻¹. The script first converts the sound speed to 1000 m/s. It then evaluates the Toomre parameter using the formula above. With those numbers, the result is about 0.00 when rounded to two decimal places in the original script output, because the chosen epicyclic frequency is extremely small compared with the self-gravity term for this surface density in SI units. The exact unrounded value is still physically meaningful, but the displayed output follows the preserved JavaScript formatting.

The same inputs also produce a critical wavelength and a fragment mass estimate. Those values can be very large because the characteristic scale depends on ${c_s}^2$ divided by $G\Sigma$. If the surface density is modest and the sound speed is not tiny, the implied unstable scale can span a large physical region. That does not mean a real disk must fragment exactly at that size. Instead, it tells you the order of magnitude of the scale that the simplified local analysis considers dynamically important.

It is also helpful to think about how each input changes the answer. If you double the surface density while keeping the other two inputs fixed, the denominator of the Toomre expression doubles, so $Q$ is cut in half. If you double the sound speed, the numerator doubles and $Q$ doubles. If you increase the epicyclic frequency, the disk becomes harder to collapse locally because orbital restoring forces are stronger. These direct proportionalities make the parameter intuitive and useful for quick sensitivity checks.

How to interpret the result

The simplest rule of thumb is that indicates local axisymmetric instability in the ideal thin-disk model, while $Q>1$ indicates stability. This page preserves a practical middle label of Marginal for values between 1 and 1.5. That middle category is not a universal law of nature; it is simply a useful warning that the system lies near the threshold and may be sensitive to assumptions, measurement uncertainty, or additional physics not included in the one-component model.

The critical wavelength should be interpreted as a characteristic scale associated with the instability criterion, not as a guaranteed fragment diameter. Real disks have finite thickness, radial gradients, turbulence, magnetic fields, and cooling processes that can shift the preferred scale or suppress fragmentation entirely. Likewise, the fragment mass is best treated as an order-of-magnitude estimate. It is valuable for intuition, comparison, and rough planning, but it should not be mistaken for a detailed prediction from a full simulation.

In practical work, the most useful way to read the output is often comparative rather than absolute. You might evaluate several radii in a disk model, several times in an evolving simulation snapshot, or several observationally inferred parameter sets. The absolute number matters, but the trend can matter even more. A region where $Q$ drops steadily toward unity is often more interesting than a region that is stably above the threshold everywhere. This calculator is therefore especially handy for quick scans and sanity checks before moving on to more detailed analysis.

Assumptions and limitations

The classic Toomre analysis assumes a thin disk and local perturbations. It is designed for axisymmetric disturbances, which means it does not directly describe every kind of structure seen in real disks. Spiral arms, bars, swing amplification, and other non-axisymmetric effects can be important even when the local axisymmetric criterion suggests stability. A disk can therefore look dynamically active while still having a local $Q$ value above unity in some regions.

Cooling is another major caveat, especially in protoplanetary disks. A low $Q$ value may indicate that self-gravity is strong enough to matter, but fragmentation can still depend on whether the gas can lose heat quickly enough. If cooling is inefficient, the disk may develop spiral structure and transport angular momentum without breaking into bound clumps. If cooling is rapid, fragmentation becomes easier. This calculator does not include thermal timescales, radiative transfer, or opacity effects.

There is also a modeling choice hidden in the sound speed input. In some contexts, users enter a true thermal sound speed. In others, they use an effective dispersion that includes turbulence or random stellar motions. That can be a sensible approximation, but it changes the interpretation of the result. The calculator does not distinguish among those cases; it simply applies the number you provide in the standard formula. For stellar disks or multi-component gas-plus-star systems, more advanced effective stability criteria are often used.

Finally, this tool is best viewed as a first-order diagnostic. It is excellent for teaching, quick estimates, parameter sweeps, and checking whether a region is plausibly near the instability threshold. It is not a substitute for hydrodynamic simulations, N-body calculations, or detailed observational modeling. If your result is close to the threshold, the safest conclusion is usually that the region deserves deeper study rather than that the answer is settled.

Reference comparison

The short table below is included only as a supporting aid. It shows how the classification changes when surface density and sound speed vary while the epicyclic frequency is held fixed. The exact numbers are illustrative, but the trend is the important part: increasing $\Sigma$ tends to lower $Q$, while increasing $c_s$ tends to raise it.

That simple competition between self-gravity, pressure support, and rotational stabilization is the heart of the Toomre criterion. Even when a full research problem requires more sophisticated tools, this local parameter remains a useful way to organize intuition and communicate what is driving a disk toward or away from instability.

Additional formula notes

Because the original page included MathML-based notation, the formulas are preserved here in MathML rather than being replaced with plain text. That matters for accessibility, machine readability, and consistency with the scientific content. The following compact identities restate the same relationships in slightly different forms that readers often encounter in textbooks and lecture notes:

$Q$ is dimensionless.

$\Sigma$ carries units of mass per area.

$c_s$ is a speed.

$\kappa$ is a frequency.

$G$ is the gravitational constant.

$\omega$ describes perturbation frequency.

$k$ is the radial wavenumber.

$\lambda$ and $k$ are related through wavelength and wavenumber concepts.

${\kappa }^2$ contributes a stabilizing term.

$\pi G\Sigma$ sets the self-gravity scale in the denominator of $Q$.

${c_s}^2$ appears because pressure support depends on the square of the effective sound speed.

$M_{\text{frag}}$ is only an approximate characteristic mass.

${\lambda }_{\text{crit}}$ is a characteristic scale, not a guaranteed observed size.

$\Omega$ is the angular frequency often used when discussing nearly Keplerian disks.

$R$ is the cylindrical radius appearing in the epicyclic relation.

These notes do not change the calculator logic, but they make the page easier to read for students and researchers who want a quick reminder of what each symbol means before entering values.