On ‘Oumuamua

‘Oumuamua slipped our grasp.

With JPL’s mission design tool, one can scroll wistfully through a fine-grained list of opportunities lost — a spaceport departure board overflowing with missed flights, all to a single exotic destination.

What could have been -- 2017-era economy-class tickets to 'Oumuamua.

By any measure, ‘Oumuamua imparted an influence that exceeded its size. Its approximately 260-meter greatest extent is dwarfed by the recently visited 2014 MU69.

Which itself is not very large.


It’s interesting to look at the Google Trends listing for worldwide interest in ‘Oumuamua.

In mid-November 2017, a two-lobed burst of interest surrounded first the naming and then the release of ESO’s iconic artist’s impression of a menacing starship-like shard. A barely perceptible blip followed last Summer’s announcement that ‘Oumuamua’s trajectory was non-Keplerian, and then, in Fall 2018, boom, a veritable megastructure of attention.

‘Oumuamua’s visit was short enough so that the data obtained (mostly in the five compressed weeks following its initial discovery) is readily summarized. Our first detected interstellar visitor arrived with almost exactly the speed and direction that one would guess a priori for an object having the average local orbit in the disk of the Milky Way. In a sense, the Solar System ran across ‘Oumuamua rather than the other way around.

The spectrum of sunlight reflecting off ‘Oumuamua’s surface is skewed smoothly to the red, with no identifying bumps or dips. The first spectrum, obtained by Joseph Masiero, who was observing at the Palomar 200-inch telescope when word of the interstellar object was released, was typical.

‘Oumuamua’s colors fall readily onto the locus defined by the various small-body constituents of the outer Solar System, which tend to list toward the reddish and the dark. Saturn’s moon Phoebe is the type of object one can bring to mind in this regard.

Solar System comets invariably spew out micron-sized particles, which are entrained in the gas (mostly water vapor) that erupts from their sun-warmed surfaces. It was thus generally expected that comets arriving from other solar systems would behave in similar fashion. Yet even in the deepest exposures — stacked together from multiple tracking images — ‘Oumuamua appeared completely point-like. Its lack of any observable coma placed stringent, spoonful-per-second limits on the amount of powdery dust emanating from its surface.

The diagram just below (with portions of the graphics taken from last October’s Sky and Telescope article) shows the highest-resolution portion ‘Oumuamua’s light curve stitched together from the data obtained with eight different telescopes. The reflected light signal varied periodically, with a dip-to-dip time scale of about 3.6 hours. There are lots of frequencies (and aliases) in the periodogram of the irregularly sampled data, the data gets proportionally noisier when the signal is at its dimmest, and the light curve does not quite repeat from pulse to pulse. The large peak-to-trough variation, and the more-or-less repetitive pattern suggest a highly elongated object, one that was tumbling chaotically end over end. The rotation period, moreover, was short enough so that ‘Oumumua would need to have some mild degree of physical strength to resist flying apart as it spins. It seems to be a solid object rather than a loosely consolidated rubble pile. The famous artist’s impression that populates a Google search on ‘Oumuamua proceeds from these basic features of the light curve.

Each image from each telescope contributed to the determination of ‘Oumuamua’s trajectory. A careful analysis of all the data leads to the remarkable conclusion that ‘Oumuamua departed the Solar System more quickly than an object starting on the same orbit and that was subject only to the Sun’s gravity. In other words, an additional, anomalous acceleration acted on ‘Oumuamua. The effect was small, apparently exactly akin to reducing the mass of the Sun by 0.1%, but its effect was nonetheless very evident in the trajectory.

Comets routinely display non-gravitational accelerations of the type and the general magnitude that was observed for ‘Omuamua. Jets of sublimating gas (mostly water vapor) stream off the surface, and perform like stochastic rocket engines. The close-range photographs of comet 67/P Churyumov-Gerasimenko (the dramatic subject of the Rosetta Mission) show the process in action.

The catch, however, is that if outgassing is responsible for ‘Oumuamua’s acceleration, there was exceedingly little fine dust entrained in the gas. This disconnect has led to two proposals in which the acceleration arose not from a gas jet, but from solar radiation pressure. The first proposal (by Shmuel Bialy and Avi Loeb) contained the much-discussed speculation that ‘Oumuamua is an engineered solar sail. A paper published several weeks ago (to less media attention) by Amaya Moro-Martin suggests that ‘Oumuamua is a naturally formed fractal, a dust bunny-like aggregate formed in a protostellar disk with an extremely low density and an extremely high porosity. Such a structure would be readily accelerated by radiation pressure, and like the sail model, would have no need to spew gas-entrained micron-sized dust.

One can argue that the lack of micron-sized dust in an outgassing scenario isn’t all that unexpected. An icy body from an alien environment may have undergone compositional processing that was absent in the Solar System. Moreover, if the entrained dust particles were larger, say 100 microns in size, they would have gone undetected.

The Spitzer Space Telescope’s non-detection of re-radiated heat constrains ‘Oumuamua’s reflectivity and its physical size. Spitzer’s wavebands of observation rendered it particularly sensitive to emission from carbon dioxide and carbon monoxide gas in ‘Oumuamua’s vicinity, but neither species of molecule was detected. This is odd. Jets from Solar System comets are primarily composed of water vapor, but typically also contain carbon-based constituents. A water vapor jet from ‘Oumuamua is not directly ruled out by Spitzer, but one needs to argue that the water, in addition to being free of small dust, had a very low level of carbonation. A glass of melted ‘Oumuamua would have to be still, not sparkling.The jet interpretation for ‘Oumuamua’s acceleration can thus wriggle out of the Spitzer non-detection by positing a composition of icy material with a very low carbon-to-oxygen ratio, and it can simultaneously wriggle out of the coma non-detection by positing that very little micron-sized dust was embedded in the ice, but that’s an admittedly uncomfortable amount of wriggling. Requiring a pair of caveats is not very satisfactory, but it allows us the hypothesis that we were visited by a relatively mundane object rather than an exotic natural object or an artificial object. Moreover, with dim prospects for the emergence of additional observational data from a long-departed object that was faint to begin with, it’s unlikely that any of the three classes of competing explanations for ‘Oumuamua’s behavior will ever gain full credence or confirmation.

Interestingly, we do have some direct information on the composition of extrasolar planetesimals. This is obtained obtained by looking at accretion events onto white dwarfs, e.g. this paper by Wilson et al. Among the measurements of this type that have been made, there are a number of cases where the C/O ratios are below the naive limit for ‘Oumuamua (plotted below on Figure 3 of Wilson et al.)

C/O ratios of planetesimals that accrete onto white dwarfs. Data and figure from Wilson et al.( 2016). The upper limit on the C/O ratio for ‘Oumuamua is shown.

A paper published by Roman Rafikov last August raised an additional potential problem with the jet model for ‘Oumuamua’s anomalous acceleration. If a jet acts continuously (or on average, acts continuously) to exert the acceleration at a single point on the surface, then the force that provides the acceleration will also exert a torque, unless it is pushing exactly along one of the principal axes of inertia. This torque would have caused ‘Oumuamua to noticeably increase its rotation rate during the days that it was under close observation.

In a new paper lead-authored by Darryl Seligman and just posted to arXiv, we revisit the jet model, and imagine that the jet is created by ice that lies just below or right at the surface of ‘Oumuamua. If this is the true situation, then ‘Oumuamua’s jet will migrate rapidly across the body, and at any given moment, it will be strongest at the spot where the Sun’s rays are directly perpendicular to the surface. If we model ‘Oumuamua as an ellipsoid, the situation looks like the illustration just below, with comet 67/P Churyumov-Gerasimenko and its sub-solar jet shown to help motivate what we’re proposing:

In an idealized situation in which ‘Oumuamua is a perfect triaxial ellipsoid with a long, an intermediate and a short axis, and where two of the axes lie exactly in ‘Oumuamua’s orbital plane, and where the jet is always located at the spot where the Sun is directly overhead, ‘Oumuamua’s motion resembles that of a pendulum. Here’s a movie made with the open-source ray tracing software POV-ray which implements the appropriate rigid-body dynamics:

In the example shown above, we’ve assigned our idealized ‘Oumuamua a surfboard-like 9:4:1 long-intermediate-short axis ratio, and we’ve positioned the short axis so that it points perpendicular to the orbital plane. The action of the sublimation jet in this model produces a pendulum-like rotation of the body and implies a long semi-axis, \(a\sim 5A_{\rm ng}P^2/4\pi^2 \sim 260\,{\rm m}\), where \(A_{\rm ng}=2.5\times10^{-4}\,{\rm cm\,s^{-2}}\) was the observed non-gravitational acceleration during the high-cadence October observations, and \(P\sim8\,{\rm h}\) was the observed peak-to-trough-to-peak-to-trough-to-peak period of a full \(-\pi \rightarrow \pi \rightarrow -\pi\) oscillation.

Interestingly, the size \(a\sim260\,{\rm m}\) implied by 8 hour period and the \(0.001 g_{\odot}\) acceleration agrees perfectly with the completely independent estimates of `Oumuamua’s size that stemmed from its measured brightness, if we assume the 10% reflectivity that is appropriate to ices that have undergone long-duration exposure to the interstellar cosmic ray flux. The implication is that ‘Oumuamua was either like a pendulum that was rocking back and forth, or perhaps one that was barely swinging “over the top”. A more complete analysis allows us to make a snapshot of the range of possible motions for the idealized case. The angular position of the idealized ‘Oumuamua during its swing is on the x-axis, the varying angular speed of the swing is on the y-axis, and the color coding shows the period for each allowed trajectory. A nice feature of the motion is that as long as the aspect ratio, a:b:c has a\(\gg\)b, the period depends only very weakly on b and c.

With our sub-stellar jet model, ‘Oumuamua doesn’t change its light curve period appreciably during the period that it was observed. This is because the torque from the sub-stellar jet spends equal time working with the instantaneous spin and working against it.

Assuming that ‘Oumuamua was a natural object, it likely had an uneven surface, possibly with regions of varying reflectivity, and it most certainly was not spinning with a principle moment of inertia aligned along its angular momentum vector. Moreover, there was likely a time-varying lag in the response of the jet, and the strength of the jet likely varied stochastically. With our ray-tracing model, we can readily gin up such complications and watch how they affect the motion. Here’s a version of the dynamics where ‘Oumuamua starts in a random orientation, and has a rough mottled surface pattern:

A more realistic model for 'Oumuamua. The sub-solar jet location and direction is shown by the wandering blue pole. (The light curve at the upper left charts the brightness of the body only, not the model axes!)

In the simulation shown just above, we’re tracking ‘Oumuamua, and watching it tumble as it flies through its escaping orbit. The point of view is from an ultra-resolving telescope orbiting along with Earth. This means that the Sun’s illumination of the ‘Oumuamua model is consistent with what an observer on Earth would have seen had they had sufficient magnification. Removing the axes and the jet, we can sum up the brightness of all the pixels in each image to get a frame-by-frame light curve. Sampling the light curve at the actual cadence of the observations, and adding the correct amount of noise permits direct comparisons between various versions of the model and the observations. The first figure below shows the real data, and then various versions of what one would get with our ellipsoid models. No attempt at curve-fitting has been made here, just a proof of concept. The second figure below shows the power spectra of the data, our models, and the zeroed-out observations. It’s clear that if one wanted to, it would be quite possible to fit the light curve pretty well with models of this general type.

Real and synthetic observations of Oumuamua from October 25th-28th 2017. The rows show the real observations (top), synthetic observations for 9:4:1 (upper middle) , 10:1:1 (middle), and 10:10:1 (lower middle) models, as well as an unvarying light curve. The solid lines show the underlying light curve for each model, and the transparent points show synthetic observations sampled at the same epochs that 'Oumuamua was observed, and perturbed with magnitude-dependant Gaussian noise inferred from the observations. The figure below shows the power spectrum for the data in each of the above plots.

There’s a strong possibility that more interstellar objects will be detected fairly soon. A necessarily shaky extrapolation from the “statistics” of one object detected over several years by the Pan-STARRS survey implies that interstellar space is teeming with ‘Oumuamuas. The numbers are impressive — of order \(10^{26}\) bodies total in the Milky Way, totaling a hundred billion Earth masses. At any given moment, roughly one such object should be in the process of threading the 1AU sphere that encompasses Earth’s orbit around the Sun.


Sculpting the orbital distributions

Saturn’s pole

Saturn’s polar axis is tilted relative to the plane in which Saturn orbits the Sun, and the plane of Saturn’s orbit is itself tilted with respect to the averaged orbital plane of the Solar System’s planets.

Pieces of popular scientific writing often start with an engaging “hook”, but the foregoing statement doesn’t do a particularly good job. Saturn and its rings do, however, do a good job of showing off their tilts — their obliquities, to use the vernacular. Saturn gradually shifts in appearance as the Sun’s illumination angle changes, and over time the creeping ring shadows even affect Saturn’s climate. Certainly, at the moment when the rings slice edge-on to the solar rays, the system presents a very different appearance than when the ring plane is inclined.

The geometry was first understood by Christiaan Huygens. By the mid-1600s, he had drawn a clear diagram showing how Saturn’s tilted pole points in a fixed direction as the planet traces its three-decade orbit.

The obliquities of Saturn and Neptune (26.7\(^{\circ}\) and 29.6\(^{\circ}\) respectively) seem odd. Uranus, tipped to its side and then some (97.9\(^{\circ}\)) is odder still. Naively, one might have expected a Solar System forming from a flat, orbiting disk of gas and dust to have ended with the equators of the giant planets lying in the average plane of the planetary orbits. Jupiter, with its axis tilted at only 3.1\(^{\circ}\) is indeed fairly conforming, but the others are all badly out of alignment. Why?

Moreover, with literally thousands of worlds now in the catalog, one also wonders if spin misalignments are rampant among the extrasolar planets. Could it be possible to infer obliquities even if we have no method to photographically resolve the planets themselves?

Left to orbit freely around a star, the tilted spin pole of an isolated planet will precess like a gyroscope. The cycling of precession is slow in comparison to the rate of the spin itself, and it stems from the torque exerted on the planet by the parent star. If the planet — like Saturn — has satellites, they orbit quickly enough to act as if they were a contributing part of the planet, and the joint set-up of planet plus moons precesses as a unit. The moons and rings stay locked to the orbital plane, and the net effect of the satellites, as far as precession is concerned, is to speed up the rate at which the cycle occurs.

Wheels within wheels… The situation grows more complicated if the orbital plane of a precessing planet also precesses. An orbit whose own pole (or orbit normal) traces a slow overhead circle presents, in effect, its planet with a moving target. With the complication induced by the precessing orbit, how does the spin axis respond?

Some intuition can be obtained by experimenting with tops. Orbital precession can be mimicked by placing a spinning top on a plate and evenly rotating the tilt of the plate, much as if panning for gold.


After some practice, one finds that the plate-top system can be transiently “locked” into a state where the spin axis precesses in the opposite direction — but at the same rate — as the plate. When this happens it’s oddly satisfying, imparting a tactile clue that resonance between the precession of an orbit and the precession of a planet’s spin might be capable of playing a dynamical role.

Another clue is supplied by the motion of Earth’s Moon. In 1693, Giovanni Domenico Cassini, who logged careful observations from the Paris Observatory, concluded that the plane of the Moon’s orbit, in the course of executing a 18.6 year precession cycle about Earth’s equator, maintains a constant angle with respect to the ecliptic (Earth’s orbital plane). He also found that the small obliquity of the moon, which is only 1.5\(^{\circ}\) (compared to 23.5\(^{\circ}\) for Earth) precesses at the same rate of one full revolution per 18.6 years.

Giovanni Domenico Cassini (1625-1712). Prior to holding the directorship of the Paris Observatory, he was the highest paid astronomer at the University of Bologna, having been appointed to his professorship by the Pope.

In other words, an arrow pointing out of the Moon’s pole always lies in the plane formed by Earth’s orbit normal and the Moon’s orbit normal. This remarkable co-precession of the Moon’s spin axis and its orbit normal received little attention until a landmark 1966 paper in the Astronomical Journal by Giussepe Colombo — who separately, achieved fame for discovering that Mercury exists in 3:2 spin-orbit resonance, turning three times on its axis for every two trips around the Sun. A few years after Colombo’s paper appeared, Stan Peale coined the term Cassini State to describe the dynamical configuration.

That something so fundamental to the motion of the Earth-Moon system was left apparently unexplained from its discovery in 1693 through 1966 seemed puzzling, so I sent a text to Konstantin Batygin.

Indeed they do:

Issues of obscurity, precedence and priority aside, Colombo used Newton’s laws of motion and gravity to demonstrate that the in-sync cycling of the Moon’s spin and orbital axes represents a minimum-energy configuration. In a frame of reference synchronized to its orbital precession, the Moon’s spin pole is analogous to a marble that frictional dissipation has brought to rest at the bottom of a curved bowl. The Moon is said to be locked in a secular spin-orbit resonance.

If a perturbation — a kick — imparts energy to the marble, it will roll around in the bowl. Likewise, if the spin pole of a body in secular spin-orbit resonance is perturbed, the direction that the pole points will wander when viewed in the precessing frame. In his 1966 paper, Colombo worked out how the trajectories traced by a perturbed spin pole will behave.

Colombo also pointed out that a simple geometric construction can be used to illustrate how the spin pole moves. First, imagine the sphere defining all the possible directions that a spin pole can point. The sphere is oriented so that its own coordinate poles are perpendicular to the overall plane of the system under consideration. In the minimum energy configuration, the spin direction, the coordinate pole, and the direction of the orbit normal all lie in a single plane that slices the sphere in half.

If the planet’s spin axis is not damped to the minimum energy configuration, it will slowly trace a path on the sphere when followed in the co-precessing frame. A detailed analysis (which appears in Colombo’s 1966 paper, and which Stan Peale augmented and corrected in 1969) shows that the set of allowed paths (level curves of the Hamiltonian) are determined by the set of possible intersections between a parabolic cylinder and the sphere. Individual paths, corresponding to individual energies of perturbation, are defined by moving the parabolic cylinder back and forth.

Remarkably, the projection of these intersection paths onto the ecliptic displays a characteristic structure that arises repeatedly in problems involving resonance. The state of co-precession keeps the spin pole of the planet fixed in the perfectly damped configuration, in which the vertex of the parabola just touches the sphere. This vertex point lies at the center of a set of banana-shaped trajectories. Then, when the spin pole of the planet is perturbed, the motion follows a trajectory where the spin pole travels along one of the banana-shaped curves. Inside the shaded region, a full traversal of the curve never entails a full 360\(^{\circ}\) accumulation of angle, and the pole is said to be librating in the resonance.

The Solar System provides several examples of secular spin-orbit resonance. Most prominent from our Earth-bound viewpoint, is the motion of the Moon. The action of tides has damped the motion of its spin pole so that it lies at the core of its sequence of banana-shaped level curves.In a pair of articles [1, 2] published in 2004, Bill Ward and Douglas Hamilton raised the remarkable possibility that Saturn’s spin pole might be librating in a frame that co-precesses with the orbital inclination of the planet Neptune. On the surface, a Saturn-Neptune secular spin-orbit resonance seems nearly unbelievable. I recall hearing Bill Ward describe the work at a conference, a year or so before the papers came out. At that time, I have to admit, I didn’t really understand the details of the talk, other than the the take-away that Neptune was somehow responsible for tipping Saturn over. Although Jupiter and Uranus exhibit substantially larger gravitational perturbations on Saturn than does Neptune, the frequency at which Neptune’s line of nodes regresses happens to almost exactly match the precession rate of Saturn’s pole. Neptune’s orbit, in a sense, shifts at a rate that cuts through the noise to provide a controlling influence that adds up for Saturn.

Once the precession rates of an orbit and a spin pole are locked together, the lock will be maintained even when the orbit’s precession rate slowly changes. As a consequence, if the rate at which the orbit precesses slows down, the planet’s spin pole will slowly tip over so that its precession rate can decline in sync. When that happens on a habitable planet, it’s time to set solar sail for the stars. Closer to home, the long-ago dispersal of planetesimals in the Kuiper Belt led to a slowing of Neptune’s orbital precession. Remarkably, this seemingly minor slowdown seems to have forced Saturn from a small initial obliquity to its current 26.7\(^\circ\).

In addition to affecting Earth’s Moon and Saturn, secular spin-orbit resonance also plays a likely role in the tilts of both Jupiter and Mars (and quite possibly Uranus). It’s fully separate from the phenomenon of spin-orbit resonance, which, for example, maintains Mercury’s spin period at an average rate that is exactly 3/2 times its orbital period.

In the Solar System, the planetary obliquities are readily measured, and have been accurately known for centuries. Orbital precession rates can be calculated either from the well-established techniques of celestial mechanics, or from direct numerical N-body integrations. Even so, secular spin-orbit resonance didn’t garner attention until Colombo’s and Peale’s papers in the 1960s, and even then, it received only limited press. In Murray and Dermott’s Solar System Dynamics, which has become a standard text, the authors state at the outset that Cassini states are not covered in their book. The possible enforced match between Saturn’s polar tilt and Neptune’s orbit went unnoticed until 2004. It thus seems like a long-shot that secular spin-orbit resonances among extrasolar planets have much chance of being a “thing”.

For a planet like Saturn, the slight decrease in the Sun’s gravity from the sub-solar point to the anti-solar point on Saturn’s surface leads to a small tidal deformation of the planet. Friction within Saturn causes Saturn’s rotation to pull this tidal deformation slightly out of alignment, with the net result being a slow decrease in Saturn’s spin rate. The rate of decrease, however, is negligible. It would take many times the current age of the Solar System for Saturn’s spin period to be tangibly modified by this effect.

Tidal forces, however, have an extraordinarily steep fall-off with distance. If Saturn were moved a hundred times closer to the Sun, to a distance where the extrasolar planets are routinely found, the Sun’s tidal influence on Saturn’s spin would be ten billion times stronger.

In the presence of strong tidal forces, the spin period of a planet on a circular or near-circular orbit is brought into sync with the planet’s orbital period. That’s the situation that the Moon finds itself in, and it is also thought to be the case for most of the shorter-period transiting planets that have been discovered by the various ground-based surveys as well as by the Kepler Mission.

In addition to synchronizing the spin, tidal forces also act to align a spin pole with the orbit normal. If, however, a planet is in secular spin orbit resonance of the type we’ve been discussing, the resonant torques can potentially balance the dissipative torques and prevent the planet from being righted.

Tidal dissipation is normally quite self-regulating. If the dissipation caused by tides is strong, then synchronization ensues, and the energy that the dissipation generates drops. If, however, a mechanism exists to thwart synchronization then significant evolution can occur. Io (and to a lesser extent Europa) provide examples. As a consequence of having its eccentricity forced by the resonant interaction between the three inner Galilean satellites, Io undergoes strong tidal dissipation, leading to the famous volcanoes that cover its surface, and to the heavy loss over time of its volatile constituents.

The famous Peale, Cassen and Reynolds article that describes Io’s dissipation belongs near the very top of a list of admired papers. It presents clean dynamical arguments that draw on disparate aspects of geophysics and celestial mechanics to make a non-trivial prediction. And indeed, the paper’s two-sentence abstract is the very model of brevity:

The dissipation of tidal energy in Jupiter’s satellite Io is likely to have melted a major fraction of the mass. Consequences of a largely molten interior may be evident in pictures of Io’s surface returned by Voyager I.

Just days after the March 2, 1979 publication of the paper, Voyager 1 flew through the Jovian system, and recorded Io’s hyperactive volcanism. Here’s a recent photo of Io from NASA’s Juno probe. The picture was taken in the infrared, where it’s pretty clear what’s going on.

In short, the Peale et al. 1979 paper is a tremendous inspiration. For years, I’ve been thinking, could something similar be done with the extrasolar planets?

The Kepler data is certainly the best place to look for opportunity. The precise timing of the planets in Kepler’s multi-planet systems gives the possibility for finding subtle effects that go beyond simple Keplerian orbital motion.

It’s well known that Kepler detected lots of multiple-planet multiple-transiting systems. The planets in these systems tend to lie in the super-Earth/sub-Neptune radius range, and typically have masses of order 5 to 10 times the mass of Earth. A zeroth-order question is what these planets are like and how they got to where they are currently observed.

There is an interesting unexplained clue in the data. One can take pairs of adjacent planets from the Kepler catalog, and plot the period ratios. What one sees is that in the vicinity of low-order orbital commensurabilities, there is a statistically structure in the distribution:

There is an overabundance, or a “pile up” of planets with orbital period ratios that are a few percent larger than the perfect 3:2 and 2:1 orbital commensurabilities, and a relative lack of planet pairs that have orbital period ratios just less than the commensurabilities. It’s as if some mechanism is acting to pry the pairs apart. Moreover, if one looks at the individual sizes of the planets in the distribution, one sees that on average, the radii of the planets that lie just wide of the commensurability are larger than the radii of the planets that have slightly smaller period ratios.

Several theorists have written papers that show this structure, termed “resonant repulsion” can be understood if the participating planets are experiencing a very high rate of tidal dissipation. The difficulty, however, has been that if the standard rate of interior energy dissipation is used, then the rate of dissipation would have to be very high. The planets would have to be extremely inelastic. Earth for example, does fall into this inelastic category because the ocean tides efficiently dissipate energy along shorelines. Most bodies in the Solar System, however, and especially the massive planets – Uranus, Neptune, Saturn and Jupiter – are far less dissipative. In the case of the Solar System’s giant planets, this difference with Earth is of order a factor of a thousand or more.

In a new paper appearing in Nature Astronomy and lead-authored by Yale graduate student Sarah Millholland, we propose a solution. If one or both planets in a pair that has a period ratio lying just outside the low-order commensurability is in secular spin-orbit resonance, and if the spin obliquity is high, then the dissipation within the planet will be large, and indeed large enough to account for the observed effect.

In many respects, the regular satellites of the jovian planets in our solar system resemble the multiple-transiting multiple planet population that was found by the Kepler Mission. Orbital inclinations and eccentricities are small in both types of systems. The orbital periods typically range from days to weeks in both cases, and the mass ratios of satellites to primaries typically tend toward one part in ten thousand. It is thus reasonable to ask why the phenomenon of resonant repulsion die to secular spin-orbit resonance is not found among the jovian satellites, all of which have tiny tilts for their spin poles.

The answer lies in the spin rates of the giant planets, all of which spin relatively rapidly, causing them to bulge significantly at their equators. Jupiter does a full turn in only 9 hours 55 minutes and is noticeably squashed when viewed through a telescope. The quadrupole moment is the jargon for the quantified degree of spin-induced structural flattening. The giant planets’ large quadrupole moments force rapid orbital precession of their satellites. The frequency is substantially higher than the spin precession rates of the satellites can keep up with. As a consequence, all of the major regular satellites of the Jovian planets have their spin axes aligned with their orbit normals.

The parent stars of the Kepler multi-transiting, multiple-planet systems spin much more slowly than Jupiter or the Solar System’s other giant planets. The stars have lost the majority of their spin angular momentum through the process of magnetic braking. The quadrupole moments of the G, K, and M stars hosting Kepler-multiple planet systems are quite small. Our own Sun spins on its axis with a 27-day period, which is fairly typical, and red dwarfs tend to spin even more slowly. As a consequence, the precession periods of the Kepler planet’s orbits are driven primarily by planet-planet interactions and not by the stellar equatorial bulges.

In the plot just below, the natural spin precession frequencies, \(\alpha\)‘s, and the orbital precession frequencies, \(\vert g \vert\)‘s, for the planets in Kepler’s multiple-transiting systems are tallied into histograms. The rate, \(\alpha\) of a planet’s spin precession depends on its internal structure, so that a planet that is highly centrally concentrated (a low \(k_2\)) precesses more slowly than one whose mass is more extended (a high \(k_2\)). The histograms for the spin and orbit rates (\(\alpha\)‘s and \(\vert g \vert\)‘s) show substantial overlap, and both reach peaks near a period of about 3,000 years.

In short, it is a suggestive coincidence that the orbital periods, the masses, the radii and the separations of the Kepler planets combine to generate similar rates of orbital precession and spin precession. This means that capture into spin orbit resonances may be quite likely for these planets.

Capture of a planet into secular spin-orbit resonance will naturally occur if the ratio of the planet’s orbital precession frequency to its spin precession frequency is slowly brought down to unity from above, that is, if \(\vert g \vert/\alpha \rightarrow 1\). This can happen if the planets in a system migrate toward orbital commensurability. This schematic diagram from our paper shows how the process works:

Simulations that track the orbits and the spins of the planets show that the spin precession and orbit precession lock into sync remarkably easily and naturally. Our paper charts several example evolutionary trajectories that look like this one:

In this particular simulation, two 5 \(M_{\oplus}\) super-Earths experience mild disk-driven migration which slowly pushes their orbits together, and, after \(\sim\)1.3 million years, binds them into a 3:2 orbital mean-motion resonance. As this mean-motion resonance capture occurs, the inner planet of the pair finds that its orbital precession rate has slowed to match its spin precession rate; it is caught in secular spin-orbit resonance. Thereafter, as the orbital precession slows still further (as a consequence of the protostellar disk dissipating), the inner planet’s axis is compelled to precess more slowly as well. In order for the planet to slow its spin precession, it is forced over on its side, to a final obliquity of more than 50\(\circ\).

The simulation charted above runs for just a few million years, but the planetary systems that Kepler observed are generally a thousand times older. The outer planet in the simulation, whose obliquity is traced with the green line in the upper panel, sees its tilt kicked up when the ratio \(\vert g \vert/\alpha\) passed through unity from below but does not end up in spin-orbit resonance. Its perturbed obliquity will drop back to zero after a few tens of milions of years. For the inner planet, however, the situation is different. Torque from the tidal dissipation in the planet balances torque from the precessing orbit, the obliquity remains constant, and an Io-like situation is produced. Obliquity-juiced tides generate ~3 million Gigawatts within the planet, roughly a thousand times the total power that Io produces, and roughly three times more energy per unit mass. The relentless dissipation draws energy from the orbit, forcing the period ratio, over time, to creep up from the initial 3:2 ratio.

The net result of this process, replicated again and again in the Kepler sample, can explain the lack of worlds near the exact m:n integer period ratios and simultaneously account for the pile-ups seen just wide of the perfect commensurabilities.

A nice feature of the theory is that it makes some predictions.

,Capture into secular spin-orbit resonance is easier if a planet has a larger radius. As a consequence, if dissipating oblique planets are what drive the Kepler pairs apart, then the planets on the right side of the period ratio gaps should be (on average) larger than those to the left sides of the gap. Pleasingly, this is exactly what is seen in the data, and it’s a feature that has gone unnoticed until now:

Moreover, larger planets are more dissipative, and so statistically, the radii of planets in the member pairs should increase as the period ratios increase. This effect, while subtler, is also present in the data.

Given the actual structure of the period ratio diagram, one can work out the amount of dissipation required to explain each pair if the ages of the systems are known. Statistically, this allows us to determine what kind of planets we’re dealing with. The details are explained in our paper, but the take away is that the planets in the Kepler-multiple systems likely tend to resemble Uranus and Neptune as far as their internal structures are concerned.

And finally, one last, as-yet untested prediction. If a planet with an orbital period in the range spanned by the Kepler-multiple planet systems has significant satellites, its precession rate will be too rapid for secular spin-orbit resonance to work. As a consequence, oblique planets driving resonant repulsion won’t have significant moons of the type seen orbiting the giant planets in the Solar System.