Exoplanets

A torque that nearly cancels

A planet embedded in a gas disc pulls on the material inside its orbit and outside it, and the two torques are almost equal and opposite. What survives the subtraction is a per cent of either, and it is still enough to carry a planet from where it formed to its star in less time than the disc lasts.

Assumes Planet migration, Resonance and Planetary rings.

A planet exists where one cannot form, and the first rung of this anchor established the problem: hot Jupiters sit at a hundredth of an astronomical unit, where the disc is too hot for ice and too sparse for enough solids, so they must have arrived from further out. This rung is about the mechanism, and about why it is embarrassingly efficient.

The mechanism is a torque between the planet and the gas it is embedded in. What makes it interesting is that the torque is a difference between two much larger numbers, so its sign and size are decided by quantities that would otherwise be third-order details of a disc model. A calculation whose answer is a one per cent residual is a calculation in which everything matters.

An edge where two torques balance, 21 kilometres from the shepherd. Two torques on the edge of a ring, against distance from a shepherding moon, both axes logarithmic and both scaled by the same combination of surface density, radius and orbital rate so that only their shapes are being compared. The flat line is the viscous torque, which comes from collisions between ring particles and does not care how far away anything is; it is drawn for a kinematic viscosity of 12 square centimetres a second, within the range ring seismology gives. The falling line is the moon's, summed over the first-order resonances that crowd together as the gap narrows, which makes it an inverse cube. A flat curve and an inverse cube cross once, and the crossing is where an edge can sit: closer in the moon wins and pushes the material back, further out viscosity wins and the ring spreads. For Daphnis and the Keeler gap the balance lands 20.8 kilometres out against a measured half-width of 21, which is agreement to well inside the uncertainty on the viscosity — and it is the only handle anybody has on that viscosity, since the quantity being inferred is the collision rate among particles a metre across, a billion kilometres away.
Fig. 1 The solar system’s own worked example of the balance this essay is about, drawn at the numbers of Daphnis and Saturn’s Keeler gap. A satellite embedded in a disc exerts a torque that falls as the inverse cube of the distance; the disc’s own viscosity resists being pushed and does not care about distance at all. Where the two cross, an edge sits. A planet in a protoplanetary disc is the same problem with a gas disc rather than a ring, and the same crossing decides whether an edge exists at all.

Two torques, and the subtraction

A planet on a circular orbit in a disc raises a spiral wake in the gas, and the wake exerts a torque back on the planet. The calculation is done resonance by resonance: at each radius where the gas’s orbital frequency and the planet’s are in a simple ratio, the perturbation is coherent and angular momentum is exchanged.

Material inside the planet’s orbit goes round faster, so the planet pulls it back and it pushes the planet forward: an outward torque. Material outside goes round more slowly, so the planet drags it forward and it pulls the planet back: an inward torque.

The two are almost equal. Each is enormous — enough to move a planet across the disc in a few hundred orbits — and their sum is a per cent or two of either, with a sign that depends on the disc’s density and temperature gradients. Because the residual is a small difference of large quantities, its sign and magnitude are sensitive to everything: the surface density profile, the temperature profile, whether the disc can radiate, and whether the gas immediately around the planet is co-rotating with it. That sensitivity is the whole difficulty of the subject.

A straight line whose slope is a mass per unit area. The measurement that the previous picture is the raw material for. Each point is one crest of the density wave: its number outwards from the resonance along the horizontal axis, and the square of its distance from the resonance up the vertical one. The points lie on a straight line through the origin to within 6.2e-14 per cent, which is what the phase accumulating as the square of the distance means, and the slope of that line is 8π²Gσr_L/3(m−1)Ω_L² — one expression containing one unknown. Solving it gives 15.0 kilograms per square metre. Everything else in the expression is known to several figures from the moon's orbit and the planet's mass, so the error on the ring's mass is essentially the error on reading crest positions off a light curve. A ring a few tens of metres thick, spread over an area larger than the Earth, is weighed by fitting a line to twenty-six numbers.
Fig. 2 The other thing a disc does with an imposed torque, and the reason the ring analogy is not exact. A self-gravitating sheet supports density waves, so a resonance in a massive ring radiates its angular momentum away rather than piling it up at one radius — and the dispersion relation at fifteen kilograms per square metre says how fast. A protoplanetary disc is pressure-supported rather than self-gravitating over most of its life, so it does the same thing through sound waves instead, and the planet’s wake is that wave. The torque this essay is about is the recoil from launching it.

Why the timescale is a problem

Put numbers in. For an Earth-mass planet at five astronomical units in a disc of the mass the solar system is thought to have formed from, the residual torque moves it inward on a timescale of a few hundred thousand years. The disc lasts a few million.

So a planet that reaches an Earth mass early in the disc’s life should reach the star before the disc disperses. The migration rate scales with the planet’s mass, so a ten-Earth-mass core migrates ten times faster still.

That is not what is observed. Planetary systems exist, with planets at every separation, and the cores that would become giant planets have to survive long enough to accrete gas. The theory as stated predicts that they do not.

What stops it

Several mechanisms have been proposed and the most robust is the one the hero figure draws: the planet stops migrating smoothly when it becomes massive enough to clear a gap around itself.

The criterion is the same balance as a shepherd moon’s. The planet’s torque tries to push gas away; the disc’s viscosity and pressure try to fill the gap in. Above a critical mass the torque wins, a gap opens — the same way a resonance clears a gap in the asteroid belt and locks a moon elsewhere — and the planet is no longer embedded in a smooth disc — it sits in a hole, with material on both sides that can only reach it slowly.

Once a gap is open the migration changes character. The planet is locked to the disc’s own viscous evolution rather than sliding through it, which is much slower, and the direction depends on which way the disc is spreading.

An edge where two torques balance, 14 kilometres from the shepherd. Two torques on the edge of a ring, against distance from a shepherding moon, both axes logarithmic and both scaled by the same combination of surface density, radius and orbital rate so that only their shapes are being compared. The flat line is the viscous torque, which comes from collisions between ring particles and does not care how far away anything is; it is drawn for a kinematic viscosity of 40 square centimetres a second, within the range ring seismology gives. The falling line is the moon's, summed over the first-order resonances that crowd together as the gap narrows, which makes it an inverse cube. A flat curve and an inverse cube cross once, and the crossing is where an edge can sit: closer in the moon wins and pushes the material back, further out viscosity wins and the ring spreads. For Daphnis and the Keeler gap the balance lands 13.9 kilometres out against a measured half-width of 21, which is agreement to well inside the uncertainty on the viscosity — and it is the only handle anybody has on that viscosity, since the quantity being inferred is the collision rate among particles a metre across, a billion kilometres away.
Fig. 3 And the same balance with the viscosity three times larger. The equilibrium edge moves inward by the cube root of three — the width goes as the inverse cube root of the viscous torque — so a threefold error in the viscosity is a forty-per-cent error in the inferred shepherd mass. That sensitivity is the whole reason the density-wave measurement of viscosity matters here: without it the shepherd mass is a free parameter, and with it the two torques can be compared.

The critical mass depends on the disc’s viscosity and thickness, and for a typical protoplanetary disc it comes out around a Saturn mass. Below that a planet migrates fast; above it, slowly. That is a mechanism for stopping giant planets and not for stopping the cores that have to grow into them, which is the part still unresolved.

An edge where two torques balance, 43 kilometres from the shepherd. Two torques on the edge of a ring, against distance from a shepherding moon, both axes logarithmic and both scaled by the same combination of surface density, radius and orbital rate so that only their shapes are being compared. The flat line is the viscous torque, which comes from collisions between ring particles and does not care how far away anything is; it is drawn for a kinematic viscosity of 12 square centimetres a second, within the range ring seismology gives. The falling line is the moon's, summed over the first-order resonances that crowd together as the gap narrows, which makes it an inverse cube. A flat curve and an inverse cube cross once, and the crossing is where an edge can sit: closer in the moon wins and pushes the material back, further out viscosity wins and the ring spreads. For Daphnis and the Keeler gap the balance lands 42.9 kilometres out against a measured half-width of 40, which is agreement to well inside the uncertainty on the viscosity — and it is the only handle anybody has on that viscosity, since the quantity being inferred is the collision rate among particles a metre across, a billion kilometres away.
Fig. 4 The same balance for a shepherd three times heavier, which is the gap-opening criterion read in the direction the planet problem needs. The equilibrium edge moves out to forty-three kilometres — the two-thirds power again — so the gap a body clears grows more slowly than its mass. That is the awkward scaling for planet formation: a core has to grow by a large factor to buy a modest increase in how far it holds the gas off, and until the gap is clean the planet is still embedded and still migrating at the embedded rate. The transition is gradual in exactly the mass range the cores have to pass through.

The other ways out

Three further mechanisms are worth naming because each solves a different part of the problem.

The first is a trap. If the disc has a place where the density or temperature gradient reverses — at the inner edge, at the boundary between magnetically active and inactive regions, or at the ice line — the residual torque changes sign there, and a planet arriving from either side stops. The second is the disc’s dispersal. A disc does not fade uniformly; it is cleared from the inside out by photoevaporation, opening a hole that expands. A planet outside the hole stops migrating when the gas inside it is gone.

The third is that the planets stop each other. Two migrating planets converging in period ratio are convergently migrating, and convergent migration captures into resonance. A captured pair migrates as a unit, more slowly than either would alone, and a chain of three or four is slower still.

The regime that runs away

There is a third mode of migration between the two described above, and it is the one that behaves worst.

A planet massive enough to perturb the disc substantially but not massive enough to clear a gap sits in an awkward intermediate state: it has carved a partial gap, so the surface density around it is no longer smooth, and the gas on horseshoe orbits in its co-orbital region is being processed at a rate comparable to the planet’s own migration.

In that situation the torque acquires a feedback. As the planet migrates, material passes across its orbit from one side to the other, and the mass flux through the co-orbital region exerts a torque proportional to how fast the planet is already moving. If the coefficient exceeds one, the migration accelerates itself.

The result is runaway: a planet can cross a substantial fraction of the disc in a few tens of orbits, in whichever direction it happened to be drifting when the feedback engaged. It is called Type III migration, it is fast enough to be an instability rather than a drift, and its direction is set by an initial condition rather than by the disc’s gradients.

An edge where two torques balance, 30 kilometres from the shepherd. Two torques on the edge of a ring, against distance from a shepherding moon, both axes logarithmic and both scaled by the same combination of surface density, radius and orbital rate so that only their shapes are being compared. The flat line is the viscous torque, which comes from collisions between ring particles and does not care how far away anything is; it is drawn for a kinematic viscosity of 4 square centimetres a second, within the range ring seismology gives. The falling line is the moon's, summed over the first-order resonances that crowd together as the gap narrows, which makes it an inverse cube. A flat curve and an inverse cube cross once, and the crossing is where an edge can sit: closer in the moon wins and pushes the material back, further out viscosity wins and the ring spreads. For Daphnis and the Keeler gap the balance lands 29.9 kilometres out against a measured half-width of 30, which is agreement to well inside the uncertainty on the viscosity — and it is the only handle anybody has on that viscosity, since the quantity being inferred is the collision rate among particles a metre across, a billion kilometres away.
Fig. 5 The same shepherd at a third of the viscosity, which is the direction the runaway regime lives in. A less viscous disc resists the torque less, so the gap is wider for the same mass — thirty kilometres rather than twenty-one — and the co-orbital region holds proportionally more of the mass deficit the feedback needs. Low viscosity is therefore where a partial gap is easiest to open and hardest to refill, which is exactly the condition the runaway requires. That the same parameter controls both is why the assessment of Type III migration turns on a quantity nobody has measured in any protoplanetary disc.

Two things make it hard to assess. It requires the co-orbital mass deficit — the difference between what is in the partial gap and what would be there in a smooth disc — to exceed the planet’s own mass, which is a condition on the disc’s structure rather than on the planet. And it is a numerical result: the linear theory that gives the Type I torque cannot contain it, because the feedback is intrinsically nonlinear, so its existence and its rate come from simulations whose resolution near the planet is the limiting factor.

A migration mode whose speed depends on the speed it already has is not a correction to the others, and whether real planets pass through it decides whether the intermediate mass range is a place a planet can linger or a place it crosses in a hurry.

The number the whole thing turns on

It is worth writing the residual out, because the smallness is the point and the smallness is not an accident.

The torque from a single Lindblad resonance carries the square of the planet’s mass, the local surface density, and a function of the disc’s aspect ratio. Summing the inner resonances gives a positive number; summing the outer gives a negative one of nearly the same size. The reason they nearly cancel is geometric: the resonances are placed symmetrically about the planet in frequency, and the disc’s properties vary slowly across the region they occupy.

What breaks the symmetry is that the outer resonances sit slightly closer to the planet in radius than the inner ones do, because of the way the epicyclic frequency runs, and that the surface density and temperature are not constant. Both effects are of order the disc’s aspect ratio, which is a few per cent — so the residual is a few per cent of either sum, and its sign depends on gradients rather than on values.

A wave whose crests count out 70 kilograms per square metre of ring. Optical depth across a spiral density wave in Saturn's A ring, drawn outwards from the 5:3 inner Lindblad resonance with Mimas at 131,988 km. The wave is launched at the resonance on the left and damps away to the right. Its wavelength is not constant: it starts at about 7.9 km and has shortened to 2.26 km by the last crest drawn, because the wavenumber grows in proportion to the distance from resonance and the ring's own self-gravity is the only restoring force in the dispersion relation. That makes the pattern a chirp with exactly one unknown in it. Fitting the 12 crest positions actually drawn here — the square of each one's distance from resonance against its number, which is a straight line — returns a surface density of 70.0 kilograms per square metre against the 70 the profile was built from. The rings have been weighed this way rather than by anything touching them: the mass per unit area follows from counting bright bands in a light curve as a star sets behind the ring.
Fig. 6 A density wave at twice the surface density, and the reason the ring analogy is worth pursuing rather than merely invoking. The wavelength of the crests is set by the surface density, so counting them weighs the material — seventy kilograms per square metre here against thirty-five before — and it is the one place in this whole family of problems where the disc’s own mass is measured rather than assumed. The planet-migration calculation needs exactly that number and cannot get it, which is why the residual torque is quoted as a rate per unit surface density and the surface density is a model.

That is why the subject is so sensitive to disc models. A change in the temperature profile that would be invisible in any observation flips the sign of the residual, and with it the direction a planet moves.

An edge where two torques balance, 10 kilometres from the shepherd. Two torques on the edge of a ring, against distance from a shepherding moon, both axes logarithmic and both scaled by the same combination of surface density, radius and orbital rate so that only their shapes are being compared. The flat line is the viscous torque, which comes from collisions between ring particles and does not care how far away anything is; it is drawn for a kinematic viscosity of 12 square centimetres a second, within the range ring seismology gives. The falling line is the moon's, summed over the first-order resonances that crowd together as the gap narrows, which makes it an inverse cube. A flat curve and an inverse cube cross once, and the crossing is where an edge can sit: closer in the moon wins and pushes the material back, further out viscosity wins and the ring spreads. For Daphnis and the Keeler gap the balance lands 10.0 kilometres out against a measured half-width of 10, which is agreement to well inside the uncertainty on the viscosity — and it is the only handle anybody has on that viscosity, since the quantity being inferred is the collision rate among particles a metre across, a billion kilometres away.
Fig. 7 The same balance for a shepherd three times lighter. The equilibrium gap narrows by the two-thirds power — twenty-one kilometres to about ten — while the satellite’s own diameter falls by only the cube root, so a lighter shepherd sits proportionally closer to the edge it maintains. That scaling is how the masses of the unseen shepherds were estimated: the gap is measured, the viscosity is measured from a density wave, and the mass is what is left.

What was actually measured

Nothing about the torque. Protoplanetary discs are observed, and planets in them are being observed for the first time, but the exchange of angular momentum between them is inferred entirely from structure.

The strongest evidence is the rings and gaps in discs themselves. Submillimetre imaging of nearby young stars shows discs with concentric gaps, at radii of tens of astronomical units, in objects only a million years old. The natural interpretation is that planets have opened them, and the gap widths imply masses of tens of Earth masses.

That interpretation is not forced. Gaps can also be produced by condensation fronts, by magnetic effects, and by instabilities in the dust rather than the gas. The discriminant would be detecting the planet, and two have been detected so far, in one system. The third piece is the resonant chains, and their fragility is informative. Chains are common in young systems and rare in mature ones, which says migration happened and something later broke most of its products.

Reading a mass off a gap

The resolved discs turned the subject from a theoretical one into an observational one, and it is worth setting out how a gap becomes a planet mass, because the chain is short and every link is contested.

A gap’s width is set by the balance the hero figure draws: the planet’s torque opening it against the disc’s viscosity and pressure closing it. Wider gap, more massive planet — with the width depending also on the disc’s aspect ratio and its viscosity, neither of which is measured directly.

A gap’s depth is set by the same balance and by how much material can flow through the gap region. Deeper gap, more massive planet, with the same two unmeasured quantities entering differently.

So width and depth give two equations, and with three unknowns — planet mass, viscosity, aspect ratio — the system does not close. What is done in practice is to assume a viscosity, which is the least well determined of the three, and read off a mass with an uncertainty that spans an order of magnitude.

The masses obtained that way for the rings in the best-studied discs come out between a few Earth masses and a Jupiter, depending on the assumption. That range is wide enough to include and exclude several of the interesting possibilities, which is why the direct detections matter so much: one planet imaged in a gap calibrates every gap in that disc.

A wave whose crests count out 35 kilograms per square metre of ring. Optical depth across a spiral density wave in Saturn's A ring, drawn outwards from the 5:3 inner Lindblad resonance with Mimas at 131,988 km. The wave is launched at the resonance on the left and damps away to the right. Its wavelength is not constant: it starts at about 4.5 km and has shortened to 0.75 km by the last crest drawn, because the wavenumber grows in proportion to the distance from resonance and the ring's own self-gravity is the only restoring force in the dispersion relation. That makes the pattern a chirp with exactly one unknown in it. Fitting the 36 crest positions actually drawn here — the square of each one's distance from resonance against its number, which is a straight line — returns a surface density of 35.0 kilograms per square metre against the 35 the profile was built from. The rings have been weighed this way rather than by anything touching them: the mass per unit area follows from counting bright bands in a light curve as a star sets behind the ring.
Fig. 8 The same wave at a higher azimuthal harmonic, which is the second thing a resolved image could in principle deliver. A wave’s crests are counted radially to weigh the disc, and its pattern — how many arms it has — says which resonance launched it, and therefore where the perturber is. In Saturn’s rings both are measured and the answer is checked against a moon that can be seen. In a protoplanetary disc the arms are marginally resolved, the perturber is usually not detected, and the inversion runs the other way: the pattern is used to place a planet that nothing else has found.

There is a further complication that has emerged as the images improved. Some discs have gaps at radii where nothing plausible could have formed and where the migration timescale from anywhere else is longer than the disc’s age. Those are read either as evidence for a mechanism other than a planet — a condensation front, a magnetically dead zone — or as evidence that the discs are older than their apparent ages, and both readings have consequences for everything else in this essay.

There is a statistical version of the same question that avoids committing to any one disc. Count the gaps across a large sample and ask whether their radii cluster: planets should be distributed smoothly, while a condensation front sits at a temperature and therefore at a radius that scales with the star’s luminosity. The observed radii do not obviously cluster once the luminosity scaling is applied, which is weak evidence for planets and is a long way from a demonstration.

The comparison that would settle it is between the gap radii and the stellar luminosities in a sample large enough to have a range of both, and the samples are only now becoming large enough to attempt it.

The awkward summary

The state of the subject is worth stating plainly rather than smoothing over, because it is unusual.

The mechanism is not in doubt: a planet in a gas disc exchanges angular momentum with it, the calculation of how much has been done many times by different methods, and the answers agree. The problem is that the answer is too large. A theory that predicts migration timescales shorter than disc lifetimes by an order of magnitude, in a universe full of planets that did not migrate all the way in, is a theory with something missing rather than a theory that is wrong.

What is missing is most likely in the co-orbital region — the horseshoe of gas that shares the planet’s orbit, whose contribution to the torque is not a small residual and whose behaviour depends on how efficiently that gas can radiate. Torques from that region can be positive, and calculating them requires radiation hydrodynamics rather than a resonance sum. There is a second reason to state the difficulty this way rather than as a failure. The migration calculation is one of the few places in planet formation where the physics is clean enough to give a number with no free parameters, and the number disagrees with the world. That is a far more useful situation than a calculation with three tunable constants that can be made to agree with anything.

What the disagreement has produced, over two decades, is a list of specific things that must be true: the disc must have structure, or the co-orbital region must contribute more than the resonances do, or planets must form later than assumed, or all three. Each of those is testable against observations that did not exist when the problem was posed — the resolved discs, the resonant chains, the radius distribution — and each of them is being tested now. A theory that was merely approximately right would have generated none of that, and the list would still be a list of adjustable choices rather than of consequences.

Worth adding, because it is the part most often left out of the summary: the difficulty is not that the residual is small but that it is small and the two large numbers are known well. Each Lindblad torque can be computed to a few per cent from linear theory, checked against numerical hydrodynamics, and it agrees. So the residual is not lost in the uncertainty of the pieces; it is a genuine, computable quantity whose value happens to depend on gradients that a disc model has to supply and observations cannot. That is a different situation from a calculation that is merely imprecise, and it points somewhere specific. The place it points is the corotation region, where the torque is not a sum over resonances at all but an average over material on horseshoe orbits, and where the answer depends on whether that material can radiate away the heat it gains as it turns. A disc that cools efficiently gives one answer and one that does not gives another of the opposite sign, and the transition between the two regimes happens at radii and surface densities that real discs occupy.

Where the ladder goes

The next rung is the co-orbital torque itself, which is where the residual’s sign is actually decided and which requires the disc’s thermodynamics rather than its geometry.

The other direction is the alternative route entirely: high-eccentricity migration, in which a planet is delivered by a distant companion driving its eccentricity up and tides circularising it afterwards. That mechanism needs no gas, works long after the disc is gone, and predicts misaligned orbits — and it is switched off entirely by the inner orbit’s own relativistic precession in part of the parameter space, which is why the observed population contains both aligned and misaligned hot Jupiters.

About the same objects

Not linked from either essay — found by the objects both name.

What links here

Essays that link to this one from their own argument.

The objects this essay names

Each one links to every other essay that touches it.

Angular momentumGap openingHot jupiterLindblad resonancePlanet migrationProtoplanetary discResonanceThe snow lineSurface densityTorqueViscosity