Delayed Feedback

Last weekend, there was an engrossing article in the New York Times.

Titled, ‘A Frankenstein’ that Never Lived, the piece’s top-line summary runs, “On Jan. 4, 1981, the effects-heavy production opened and closed on the same night. Forty years later, the creators revisit a very expensive Broadway flop.

According to the article, prior to the open, hopes for success of the big-budget enactment of Mary Shelley’s 1818 classic had run high, but the show’s prospects were dashed in large and immediate measure by Frank Rich’s dreadful review in the Times. Rich’s write-up brims with both arch sarcasm and gut-punch lines such as, “we feel nothing except the disappointment that comes from witnessing an evening of misspent energy.” Reading the article, I can feel a queasy, visceral sympathy for everyone who worked on that production.

At the end of the last century, Fred Adams and I were riding high on our forecasts for the denouement of the entire cosmos. Our trade book — The Five Ages of the Universe — had, to our great thrill, just been published by an imprint of Simon and Schuster, and we were enjoying the modest acclaim that had proceeded from dividing the entire past and the entire future of the Universe into five thumbnail-friendly eras of time.

The sales, the buzz, and the idle dreams of pop-culture stardom all came to a crashing halt from one day to the next when the New York Times ran its review of our book. I can recall the sinking feeling at the moment when our agent called Fred to let us know what had happened. “Brace yourself.” I think that was what she said.

Dick Teresi’s review was entitled, “The First Squillion Years”, and the title telegraphed the intent. It starts out bad, it gets progressively worse, and finally ends rather crushingly, with “Imagine an astronomer looking back at us ages hence. Perhaps he can read the bound galleys in my hand. Maybe he will be kinder than I have been here.” Yikes.

It’s amazing, however, how quickly one recalibrates. The book tour that had been in the cards was scaled way back to a handful of local appearances, with Borders Books in Ann Arbor standing in as the high point. Sales dried up. Just like that, I was back to trying to find bugs in Fortran codes. Over time, it became clear to me that it was almost certainly just as well.

Rather amusingly, the distant future is currently having something of a moment. There have been several new trade books on the topic, including at least one that got a review in the Times that was entirely different in tone that the one we drew. And then, a fantastic coda for the whole episode arrived unexpectedly a week or two ago as I was leafing through a recent issue of the New Yorker.

I’ve always enjoyed Roz Chast’s cartoons. She has a great sense of humor. It was oddly gratifying to see that she’d somehow been seized to lift out our eras and our thumbnails, elevating them to a full-page cartoon.

The Black Clouds

Time has a way of sliding by. More than five years ago, I wrote a blog post with a suggestion for a new cgs unit, the oklo, which describes the rate per unit mass of computation done by a given system:

1 oklo = 1 bit operation per gram per second

With limited time (and limited expertise) it can be tricky to size up the exact maximum performance in oklos of that new iPhone 12 that I’ve been eyeing. Benchmarking sites suggest that Apple’s A14 “Bionic” SoC runs at 824 32-bit Gflops, so (with Tim Cook in charge of the rounding) an iPhone runs at roughly a trillion oklos. Damn.

Computational energy efficiency also keeps improving. Proof-of-work based cryptocurrencies such as Bitcoin have laid the equivalency of power, money and bit operations into stark relief. The new Antminer S19 Pro computes SHA-256 hashes at an energy cost of about one erg per five million bit operations, a performance that’s only six million times worse than the Landauer limit. There is still tremendous opportunity for efficiency improvements, but the hard floor is starting to come into view.

For over a century, it’s been straightforward to write isn’t-it-amazing-what-they-can-do-these-days posts about the astonishing rate of technological progress.

Where a calculator like ENIAC today is equipped with 18,000 vacuum tubes and weighs 30 tons, computers in the future may have only 1000 vacuum tubes and perhaps weigh only 1½ tons.

And indeed, there’s a certain intellectual laziness associated with taking the present-day state-of-the-art and the projecting it forward into the future. Exponential growth in linear time has a way of making one’s misses seem unremarkable. “But I was only off by a few years!”

Fred Adams was visiting, and we were sitting around the kitchen table talking about the extremely distant future.

We started batting around the topic of computational growth. “Given the current rate of increase in of artificial irreversible computation, and given the steady increase in computational energy efficiency, how long will it be before the 1017 Watts that Earth gets from the Sun is all used in service of bit operations?”

As with many things exponential, the time scale is startlingly sooner than one might expect: about 90 years. In other words, the long-term trend of terrestrial economic growth will end within a lifetime. A computational energy crisis looms.

In any conversation along such lines, it’s a small step from Dyson Spheres to Kardashev Type N++. If you want to do a lot of device-based irreversible computation, what is the most effective strategy?

A remarkably attractive solution is to catalyze the outflow from a post-main-sequence star to produce a dynamically evolving wind-like structure that carries out computation. Extreme AGB (pre-planetary nebula phase) stars have lifetimes of order ten thousand years, generate thousands of solar luminosities, produce prodigious quantities of device-ready graphene, and have photospheres near room temperature.

We (Fred and I, along with Khaya Klanot and Darryl Seligman) wrote a working paper that goes into detail. In particular, the hydrodynamical solution that characterizes the structures is set forth. We’re posting it here for anyone interested in reading.

As we remark in the abstract, Possible evidence for past or extant structures may arise in presolar grains within primitive meteorites, or in the diffuse interstellar absorption bands, both of which could display anomalous entropy signatures.

Proximate

For well over a year now, oklo.org has languished on Bluehost’s servers. No new posts. No maintenance whatsoever. The PHP installation was so dated that for months, visits to the site received nothing but a cryptic SQL fail. An hour with technical support was sufficient to bring the site back up, but the result was disheartening. The URL was landing hard, with an insecure SSL warning, a creaky WordPress theme that shouted “mid-aughts”, and an alarmingly wanton disrespect for the aspect ratios of the once-carefully embedded images.

New decade. New resolve.

The name on the site header now matches the site URL. Systemic referred to the long-lapsed Systemic Console project which Aaron Wolf and I started in 2004. Later, Stefano Meschiari took over the code base, and developed the software to pearlescent degrees of refinement. Yet the intangibles that generated the thrill in the Doppler velocities have all but faded, now that an essentially Earth-mass planet orbits at an essentially temperate distance from Proxima Centauri, and rumors circulate of a signal of order a GHz rather than a meter per second.

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.

A Final Exam…

Latest Deep Space Climate Observatory Image

Over the past several semesters at Yale, I’ve been working out a new take on the standard “Astronomy 101” class for non-science majors. Broadly, the goal is to stage a wide-angle view of the Anthropocene, thereby forging an understanding of how Earth fits into its broader cosmic context. Economics, Political Science, and History constituted the largest groups of majors in the class.

I’m working on getting the class notes, problem sets and readings into a form that’s distributable. In the interim, I’ll cut right to the chase with the final exam. Per Yale’s official instructions:

Final examinations normally last either two or three hours but, in either case, students are permitted to take an additional half hour before being required to turn in their answers. This additional time is given for improving what has already been written, rather than for breaking new ground.

Link to .pdf version of the AY 105 Final Exam

This one was set as the three-hour variety. (Didn’t seem like any utility would be gained by imposing time pressure.)

Intercept

‘Oumuamua breezed in unexpectedly and it left in a rush. Faded now, to twentynine, soaring up and out over Jupiter’s orbit. No sum, it seems, sufficient to compel it to pick up the phone, to give us a call.

Maybe it was a one time fluke — a color out of space, but it’s also possible that it was unexceptional, a mundane representative from a vast distribution. If so, what can we do to be ready for the next one?

Darryl Seligman has a new paper up on arXiv that outlines a plan. Had ‘Oumuamua been spotted on its way in, and if a probe had been loitering in anticipation, fueled and ready to go at L1, it would have been an easy thing (energetically at least) to rocket over and intercept it, Deep Impact style, in a blaze of glory.

With LSST set to start monitoring the skies, there should be an opportunity every decade or so to “get interstellar” by barely leaving home.

Until Time Without End

‘Oumuamua’s encounter with the inner solar system is dying down on Twitter, yet still it bristles with consequence and the uneasiness of unanswered questions. Why no coma?

Occam’s razor is a dull instrument that points almost unerringly to the mundane (as opposed to pointing to interstellar probes). One thus draws several conclusions. (1) ‘Oumuamua’s aspect ratio is substantially less than 10:1. (2) Billions of years in the interstellar environment lead to the buildup of a tarry crust that resists temporary heating, and this process is enhanced for comet-like planetesimals that form in systems with supersolar C/O ratios. (3) Most stars have true-Neptune analogs.

The resulting prediction is that slightly tweaked ongoing surveys, and soon LSST, should start turning up interstellar asteroids and perhaps interstellar comets with some frequency. If another one is found in the near-term, it would be interesting to look at the optimal mission designs that could accomplish an opportunistic sample-return.

From ‘Oumuamua’s perspective, the close encounter with the Sun was a near-indescribable stroke of luck. To scale, the stars of the galactic disk are like grains of sand separated by miles and crawling through space at a few feet per year. The Galaxy is the archetypal collisionless fluid. Vaulting from ‘Oumuamua’s current encounter to its next connects the all too human interval of waking-up-at-3AM anxieties — the scale of days and months — to the frigid waste of a quadrillion years.

Why cold? When fusion has ended, dark matter annihilation and proton decay take over, and both (while uncertain) are certainly slow processes. Grand Unified Theories predict that proton decay should occur, but so far, there is no experimental evidence. The lower bounds on the proton half-life are ~10^34 years via the sluggishly competing processes of positron and muon decay.

If the proton were completely stable, the end states of stars present a curious state of affairs. Black holes of stellar mass, which are much more tightly bound than degenerate stars, will evaporate through the Hawking effect with a lifetime of “only” 10^66 years Although this time scale is aggressively long compared to the current 13.8-billion year age of the universe, it would be odd if black holes are ephemeral while white dwarfs and neutron stars are forever.

While jarring, this possible divergence of lifetimes is not exactly a matter of pressing concern. Two decades, ago, however, Fred Adams and I had priorities that were definitely skewed toward the really long term. Along with Manasse Mbonye and Malcom Perry, we looked into how quantum tunneling into black holes can erode white dwarfs. In Freeman Dyson’s 1979 article, Time Without End, it is pointed out that an otherwise stable white dwarf will spontaneously tunnel into a black hole on a time scale of order 10^10^76 (!) years. In our article, we argued that the whole star need not make the plunge at once, and that a 10^45 year half-life is a plausible value for black-hole induced proton decay. This has the added benefit of enabling a Hertzsprung-Russell diagram that traces stellar evolution to its absolute bitter end.

Visitors

‘Oumuamua. Up close and alongside, in the vastness of interstellar space, its hurtling bulk imparts no sense of motion as it turns imperceptibly on its axis, blotting out the stars.

For a hundred years, the point-like Sun grew steadily brighter against its frigid airless horizons. First came light, then warmth, and finally searing illumination of the tarry reddish expanse, blistering sluggishly beneath a September Noon far more intense than any summer of Earth.

`Oumuamua is departing the solar system as rapidly as it arrived, heading outward at a current rate of 2.5 million miles a day. Our tiny chance of sending a probe to catch it diminishes with each lagging tick of inactivity. Nonetheless, world-wide interest is mounting, in part as a consequence of two new articles reporting detailed observations. The first, by Jewitt et al. was posted to arXiv last week, while the second, by Meech et al. (which independently comes largely to the same overall conclusions), appeared in Nature earlier this week. Nature being Nature, the Meech et al. article was accompanied by a media push, spearheaded by an extraordinary piece of space art.

Maybe it’s press release fatigue from one “habitable” world after another — a monotony of warm suns glinting off imaginary oceans — that makes this image so arresting.

The observational facts remain stark and limited. `Oumuamua’s double-peaked light curve suggests that it has a large aspect ratio, perhaps as high as 10:1. Assuming that it’s a poor reflector, it’s several hundred meters on its long axis. Its overall color is reddish. It has to have physical strength, or its 7-hour rotation period would be enough to overcome its negligible self-gravity and tear it apart. Most alarmingly, it shows no sign of a coma. At most, less than a sugar cube’s worth of cometary dust per second was emanating from it as it tore through the inner solar system. (As a matter of fact, ‘Oumuamua as observed is entirely consistent with Tintin’s rocket.)

For more on ‘Oumuamua, I have a blog post up at Scientific American.

Interstellar Asteroids

This was no fruit of such worlds and suns as shine on the telescopes and photographic plates of our observatories. This was no breath from the skies whose motions and dimensions our astronomers measure or deem to vast to measure. It was just a colour out of space — a frightful messenger from unformed realms of infinity…

Aww, come off it.

Wild-eyed extravagances aside, A/2017 U1 — the asteroid-like visitor from interstellar space — is an extraordinary object. In traversing the gulfs, its next encounter with a star that is as close as last month’s encounter with the Sun likely won’t occur for another quadrillion years, and so the mere fact that it zipped through suggests that quite a few interstellar asteroids are out there. And this, in turn, has some remarkable consequences. A straightforward cross-section based estimate suggests that the galaxy contains of order a hundred billion earth masses of A/2017 U1-like planetesimals. Hot Jupiters, terrestrial planets, and super-Earths are all incapable of using gravity-assist to eject bodies out of their parent systems, leaving the strong hint that as-yet undetected Neptune-like planets must be extremely common.

In general, extrapolations from a sample size of one don’t have a good track record. Exhibit A would be our own Solar System — hot Jupiters were discovered at better than 100-sigma significance because solar-system expectations had been projected throughout the galaxy; proper planetary systems should have terrestrial bodies near 1 AU and gas giants at 10 AU.

The arrival of A/2017 U1 seems nicely timed to revival of the AAS’ new low-maintenance communication channel, the “Research Note“:

The purpose of the Research Notes is to provide a home for short submissions that are not suitable for publication as a journal article, but are likely to be interesting or useful to members of our community. Appropriate submissions would include brief summaries of work in progress, comments and clarifications, null results, and timely reports of observations (such as the spectrum of a supernova), as well as results that would not traditionally merit a full paper (e.g., the discovery of a single unremarkable exoplanet, a spectrum of a meteor, or contributions to the monitoring of variable sources).

I especially like the part about “single unremarkable exoplanets” being equivalenced to the “spectrum of a meteor”. In any event, Prof. K. Batygin and I have just submitted a research note that gives our take on the implications of A/2017 U1. Here’s a link to a draft of the note, which we’ll also post on the arXiv within the next several days.