Magnetorotational Instability Growth Rate Calculator

JJ Ben-Joseph headshot JJ Ben-Joseph

Introduction to Magnetorotational Instability Growth Rates

The magnetorotational instability, usually shortened to MRI, is a central mechanism in accretion-disk physics because it explains how a weak magnetic field can unlock turbulence in a differentially rotating flow. In an idealized disk, the field couples neighboring rings of gas so that the inner ring gives up angular momentum while the outer ring gains it, and that exchange lets matter drift inward. In that sense, the MRI is the bridge between smooth rotation and the turbulent transport that powers accretion.

This calculator focuses on the textbook MRI setup: an axisymmetric perturbation in an incompressible disk threaded by a uniform vertical magnetic field. With the local angular velocity Ω, field strength B, and mass density ρ, it estimates the Alfvén speed, the most unstable wavelength, the maximum linear growth rate, and the matching e-folding time. Those four outputs are enough to judge whether the MRI is likely to fit inside the local disk and how rapidly it would amplify once it starts.

The scope is intentionally narrow. It does not model full nonlinear turbulence, vertical stratification, or non-ideal magnetohydrodynamic effects such as resistivity, Hall drift, or ambipolar diffusion. Instead, it gives a clean first-pass answer from the standard vertical-field MRI result that is commonly used for order-of-magnitude estimates.

How to Use the MRI Growth Rate Calculator

Use this MRI growth-rate calculator by entering three positive quantities in SI units. The first input is the angular velocity Ω in radians per second, which sets the local orbital timescale and captures the differential rotation that drives the instability. In a Keplerian disk, Ω decreases outward, and that shear is the background the MRI feeds on.

The second input is the vertical magnetic field B in tesla. The field can be weak; the important point is that it provides magnetic tension between neighboring fluid elements. The third input is the mass density ρ in kilograms per cubic meter. Density matters because it controls inertia and therefore the Alfvén speed, which is the characteristic propagation speed for magnetic disturbances in the gas.

After you press the compute button, the calculator reports vA, λmax, γmax, and the e-folding time tgrowth. The e-folding time is often the quickest way to interpret the result: if a perturbation grows like eγt, then one e-folding time is the interval needed for the amplitude to rise by a factor of about e, or 2.718. A shorter time means a faster-growing MRI mode.

For the formulas used here to remain valid, keep the units consistent and enter strictly positive values. Zero or negative inputs do not describe the ideal MRI case this page is built around, so the calculator will reject them.

Formula for the Magnetorotational Instability Growth Rate

The MRI growth rate used on this page comes from the classic linear dispersion relation. For the axisymmetric, incompressible, vertical-field case, perturbations of the form exp(γt + ikz) satisfy the relation below. It shows how rotation, magnetic tension, and the epicyclic response of the disk combine to determine whether the disturbance grows.

γ 4 + ( κ 2 + 2 k 2 v A 2 ) · γ 2 + k 2 v A 2 ( k 2 v A 2 + d Ω 2 d ln r ) = 0

Here, κ is the epicyclic frequency, k is the vertical wavenumber, and vA is the Alfvén speed. For a Keplerian disk, the epicyclic frequency equals the orbital angular velocity, so κ = Ω. In that common MRI limit, the fastest-growing mode has a maximum growth rate of

γmax = (3/4)Ω

and occurs at a wavenumber

kmax = Ω / (√2 vA).

The corresponding wavelength is λmax = 2π / kmax, and the e-folding time is tgrowth = 1 / γmax. The Alfvén speed is given by the standard magnetohydrodynamic expression

v A = B μ 0 ρ

where μ0 is the permeability of free space. These are the exact relations used by the calculator script, so the displayed values come directly from the ideal vertical-field MRI model rather than from a fitted approximation.

MRI Growth Rate Interpretation

The MRI outputs are easiest to read when each one is tied to a physical question about the disk. The Alfvén speed tells you how efficiently the magnetic field can transmit a force through the gas. A stronger field or a lower density gives a larger vA. The most unstable wavelength tells you the vertical size of the fastest-growing disturbance. If that wavelength is much larger than the local disk thickness, the ideal mode may not fit in the actual system. The growth rate tells you how quickly the perturbation amplifies, while the e-folding time translates that rate into a more intuitive timescale.

A useful MRI rule of thumb is that the maximum growth rate in a Keplerian disk is tied directly to the orbital frequency. Since γmax = 0.75Ω, the instability grows on a dynamical timescale rather than a slow diffusive one. That is why the MRI is so effective as an accretion driver: once the conditions are right, only a handful of orbits are needed for a small perturbation to become dynamically important.

MRI Example: a Protoplanetary Disk at 1 AU

Consider a simple MRI example in a protoplanetary disk. Suppose the local angular velocity is Ω = 2 × 10−7 rad/s, the vertical magnetic field is B = 1 × 10−7 T, and the density is ρ = 1 × 10−9 kg/m³. These are rough order-of-magnitude values often used for illustrative work near 1 AU in a young stellar disk.

First compute the Alfvén speed. With μ0 = 4π × 10−7 H/m, the result is about 2.82 m/s. Next compute the most unstable wavenumber using kmax = Ω/(√2 vA). That gives a wavelength λmax of roughly 1.25 × 108 m. The maximum growth rate is γmax = 0.75Ω = 1.5 × 10−7 s−1, so the e-folding time is about 6.67 × 106 s, or roughly 77 days.

What does that mean physically? It means that, under these ideal MRI assumptions, a tiny perturbation would grow by a factor of e in a little over two months. After several e-foldings, the linear description would break down and the disturbance could move into the nonlinear turbulent regime. On the timescales of a disk that persists for thousands to millions of years, that is extremely rapid.

Reference MRI Scenarios

The table below compares a few illustrative MRI environments. The numbers are not intended to be exact models of any one object. Instead, they show how the same growth-rate formula behaves across very different astrophysical settings, from cool protoplanetary disks to compact-object accretion flows.

Scenario Ω (rad/s) B (T) ρ (kg/m³) λₘₐₓ (m) t_growth (s)
Protoplanetary disk at 1 AU 2×10⁻⁷ 1×10⁻⁷ 1×10⁻⁹ 3×10⁸ 6×10⁶
AGN disk at 10 r_s 1×10⁻³ 1×10⁻³ 1×10⁻⁴ 2×10⁵ 1×10³
X-ray binary disk 1 1 1 9 1.3

Limitations and Assumptions for MRI Growth Rates

This calculator assumes the ideal, axisymmetric, vertical-field MRI in a local Keplerian setting. That is a standard and very useful approximation, but it is still an approximation. Real disks can be vertically stratified, compressible, partially ionized, and threaded by magnetic fields that are not purely vertical. In some environments, especially protoplanetary disks, non-ideal effects such as Ohmic resistivity, Hall drift, and ambipolar diffusion can strongly alter or even suppress MRI growth. If the gas is poorly ionized, the magnetic field may not couple well enough to the fluid for the ideal formulas to remain accurate.

The calculator also reports the fastest-growing wavelength from local linear theory. In a real disk, that wavelength has to fit inside the available geometry. If λmax is larger than the local scale height or larger than the region over which the background conditions are roughly uniform, then the ideal fastest-growing mode may not be physically realized. Likewise, the growth rate shown here applies only during the linear phase. Once the perturbation becomes large, nonlinear saturation determines the eventual turbulence level, and that requires simulations or more advanced modeling.

Another practical limitation is that the result assumes a Keplerian-like rotation law with outwardly decreasing angular velocity. If the rotation profile differs substantially, the exact growth rate and unstable range of wavelengths can change. The page therefore works best as a quick MRI estimate for standard accretion-disk conditions rather than as a complete stability solver for arbitrary rotating plasmas.

Why This MRI Calculator Is Useful

Despite those limitations, a compact MRI calculator is valuable because it turns abstract linear theory into immediate scale estimates. Researchers, students, and science communicators often want to know whether a chosen field strength is dynamically relevant, whether the unstable wavelength is microscopic or macroscopic, and whether the instability grows in seconds, days, or years. Those questions can be answered quickly with just three inputs. The result is especially helpful when comparing MRI to other timescales such as the orbital period, sound-crossing time, or viscous evolution time.

In short, this page gives a practical entry point into one of the central instabilities of astrophysical fluid dynamics. If the output shows a short e-folding time and a wavelength that fits within the system, then the MRI is likely to be dynamically important. If not, you may need to consider whether the field is too weak, the density too high, the geometry too constrained, or whether non-ideal effects dominate the local plasma behavior.

Enter positive SI values to estimate the fastest-growing MRI mode in the ideal vertical-field case.

Use the local orbital angular velocity in radians per second. Enter the vertical magnetic field strength in tesla. Enter the local mass density in kilograms per cubic meter.
Enter disk parameters to estimate the fastest-growing MRI mode.