systemic

Vorticity can be thought of as the tendency of a paddlewheel to spin if placed in the flow. High vorticity is a large counter-clockwise spin, zero vorticity is no spin, and a large negative vorticity is a tendency to spin clockwise. Jonathan Langton’s climate models of short-period extrasolar planets show a remarkable variety of vorticity patterns on their surfaces, in keeping with the incredibly stormy and complex nature of their atmospheres. Here’s a gallery of Mercator-projection vorticity maps for the known strongly irradiated Jovian planets that have significant eccentricities. The red arrows indicate the wind speeds and directions across the planetary surfaces. These figures are all from a paper that’s currently under review at the Astrophysical Journal (see here for an overview of the numerical method that we’re using). Also, a shout-out is due to Edward Tufte for advocating the strong graphic-design effect of small spots of saturated color on a gray-scaled backdrop.

HAT-P-2b:

Here are 1.1 MB North Pole, South Pole and Mercator Projection animations of the HAT-P2b vorticity evolution.

HD 80606 b:

1.1 MB Mercator animation here.

HD 185269 b:

1.1 MB Mercator animation here.

HD 108147 b

1.1 MB Mercator animation here.

HD 118203 b:

1.1 MB Mercator animation here. The animations above are hosted on the Oklo Corporation’s servers.

It’s interesting to compare the vorticity maps with the temperature distributions on the planetary surfaces (shown in the same order as above):

Image Source.

The TrES survey announced the discovery of a new transiting planet today, raising the number of known transits to twenty (including Mercury and Venus). The new planet, “TrES-4″, has a mass of order 84% that of Jupiter, and with a radius of 1.67 Rjup, it’s pumped to nearly five times Jupiter’s volume:

The false color image of Jupiter was produced from near-infrared data obtained with the Gemini telescope. The even more luridly false-color representation of TrES-4 is based on a vorticity map from one of Jonathan Langton’s recent simulations.

In order for TrES-4 to be swollen to its current size, it needs to be experiencing heating of order 6×10^27 ergs per second. One way to do this is to have a significant perturbing companion which drives large time-averaged variations in TrES-4′s orbital eccentricity. So far, there are only four published radial velocities for TrES-4, so the orbit could easily be non-circular. More provocatively, if strong orbital forcing is indeed occurring, then there’s a reasonable chance that the perturber might also be observable in transit. I recommend that Transitsearch.org observers keep this bad boy under constant supervision.

Image Source.

HAT-P-2b. The name doesn’t exactly ring of grandeur, but this planet — a product of Gáspár Bakos’ HAT Net transit survey — is poised to give the Spitzer Space Telescope its most dramatic glimpse to date of a hot Jupiter.

HAT-P-2b’s orbit is remarkably eccentric for a planet with an orbital period of only 5.6 days, and by a stroke of luck, periastron is located almost exactly midway between the primary and the secondary transits (as viewed from Earth). The strength of the stellar insolation at periastron is nine times as strong as at apastron, which more than guarantees that the planet will have disaster-movie-ready weather.

On June 6th, Josh Winn and his collaborators used the Keck telescope to obtain 97 radial velocities for HAT-P-2. The observations were timed to occur before, during, and after primary transit, and the Rossiter-McLaughlin effect is clearly visible in their data (preprint here):

The symmetry of the Rossitered points indicates that the angular momentum vector of the planetary orbit is aligned with the spin pole of the star:

This state of affairs also holds true for the other transiting planets — HD 209458b, HD 149026b, HD 189733b — for which the effect has been measured. The observed alignments are evidence in favor of disk migration as the mechanism for producing hot Jupiters.

With its apparent magnitude of V=8.7, the HAT-P-2b parent star is roughly ten times brighter than the average planet-bearing star discovered in a wide-field transit survey. The star is bright enough, in fact, to have earned an entry in both the Henry Draper Catalog (HD 147506) and the Hipparcos Database (HIP 80076), but with its surface temperature of 6300K (F8 spectral type) it was too hot to have been a sure-fire “add” to the ongoing radial velocity surveys. Prior to this May, it had been entirely ignored in the astronomical literature (save a brief mention in this paper from 1969).

HAT-P-2′s intrisic brightness and its planet’s orbital geometry mean that in a relatively compact 34-hour observation, Spitzer can collect on the most interesting features of the orbit with high signal-to-noise. In particular, there is an excellent opportunity to measure the rate at which the day-side atmosphere heats up during the close approach to the star. The planet, in fact, presents such a remarkable situation that a block of Director’s Discretionary time was awarded so that the observations can be made during the current GO-4 cycle. They’ll be occurring soon.

Both HAT-P-2b and HD 80606 b will provide a crucial ground truth for extrasolar planetary climate simulations. Jonathan Langton’s current model, for example, predicts that that the temperatures on HAT-P-2b will range over more than 1000K. At the four times shown in the above orbital diagram, the hemisphere facing Earth is predicted to show the following appearances:

Spitzer, of course, can’t resolve the planetary disk. It measures the total amount of light coming from the planet in chosen passband. At 8-microns, the planet’s light curve should look like this:

The temperature maps only hint at the complex dynamics of the surface flow. A better indication is given by the distribution of vorticity,

which we’ll pick up in the next post…

One of the systemic project’s most important goals is to reach the point where we can have genuinely realistic and aesthetically satisfying simulated images and animations of extrasolar planets. This will serve the scientific purpose of allowing us to get better comparisons with infrared observational data, and will ultimately allow us to embark on vicarious missions to the new-found worlds of our Galaxy.

In the interim, there’s a lot of coding and computing to do.

As described in this post from earlier this month, I’m advising UCSC Physics graduate student Jonathan Langton on a Ph.D. thesis geared to simulating the atmospheres of irradiated extrasolar planets. Jonathan finally graded his way through a horrific stack of lab reports and final exams, and has now been able to put full focus on the research. Progress is evident in a heavy stream of e-mails containing increasingly detailed animations.

Jonathan is using a numerical technique known as the pseudo-spectral method to do his simulations. The key idea is that the flow pattern on the surface of the simulated planetary sphere can be broken down into a superposition of Fourier modes. For example, as one moves around the planet, the longitudinal variations in the flow can be described in terms of a superposition of sinusoidal patterns. Sinusoids have analytically computable derivatives, which allow one to make a highly accurate representation of changes in the flow without resorting to a cripplingly large amount of computation.

Spectral methods have their drawbacks, however, in the form of high-frequency numerical noise. This noise was evident in the earlier simulations in the form of transient ribbed patterns within the flow. Over the past few days, Jonathan has designed an elegant filtering scheme which seems to be working very well in supressing these spurious features without killing the actual structures in the flow.

The snapshot above is from a test-calculation that implements the new filtering scheme. It’s part of an animation that simulates the development of an initially random vortical flow on the surface of a planet with the radius, mass, and rotation period of HD 209458 b. (Potential vorticity is the quantity plotted, resolution is 256×128, and the simulation runs for 5 rotational periods).

Here is a link to Jonathan’s latest (7MB) animation. It’s hot off the computer.

Regular oklo readers all know that HD 209458 b is a lot bigger than it’s supposed to be.

In a previous post, we saw that the theoretical models that provide reasonable matches for the other 9 transiting planets predict that HD 209458 b’s radius should be slightly larger than Jupiter’s radius. The observations, on the other hand, make it clear that the planet is actually has a diameter about 1.35 times larger than Jupiter. HD 209458 b is by far the best-studied exoplanet, so it’s of more than passing interest to understand why it’s so large.

There’s general agreement that HD 209458 b must be privately tapping an unusual source of internal heat. Somehow, a lot of extra energy is being generated in the planetary interior. The surplus heat allows the planet to maintain an expanded outer envelope, and hence endows the planet with a larger overall size. The big question is: what is the anomalous extra heat source? A few years ago, I was enthusiastic about the idea that there might be a second companion planet that is gravitationally perturbing HD 209458 b, forcing it to maintain a slightly eccentric orbit. This eccentricity would be continually damped as a result of tidal interactions with the parent star, which would generate a sufficient amount of interior heating. Such a state of affairs is analogous to the heating of the inner Jovian satellites. The heat generated by tidal friction lends Io its off-the-hook volcanism, and maintains Europa as the Astrobiology poster world.

Unfortunately, however, the perturbing companion model no longer seems to be a viable explanation of HD 209458 b’s large size. That is, if you use the systemic console to fit to the HD 209458 b radial velocities, you’ll find that there is very little latitude for inserting significant extra planets. Try it and see, and upload your fits to the systemic back-end.

Josh Winn (MIT) and Matthew Holman (Harvard-Smithsonian CfA) have written a paper that presents an interesting hypothesis for resolving the HD 209458 b radius dilemma. Winn and Holman propose that the planet is caught in a so-called Cassini state, which is a resonance between spin precession and orbital precession. In a future post, I’ll give a heuristic discussion of the dynamics of how this situation can arise, and how the Cassini states work, but in short, if HD 209458b is trapped in the “Cassini state 2″, then its spin axis will lie almost in the orbital plane. Like all hot Jupiters, the planet will spin once per orbit, but it will literally be lying on its side as it orbits the parent star. A synchronous planet in Cassini state 2 will experience a large amount of tidal heating, even in the complete absence of any other planets in the system.

I like the Winn-Holman hypothesis because it’s potentially testable. If the planet is in Cassini state 2, then the pattern of illumination on the surface, and hence the time-dependant global infrared signature, will be very different than if it is locked into the standard upright configuration. In the standard scenario, a hot Jupiter has a fixed substellar point on its equator that does not wander significantly as the planet executes its orbit. One hemisphere of the planet is in perpetual day, while the other hemisphere experiences an endless night. Hydrodynamic calculations by James Cho and his collaborators (link), and by Adam Showman and his students (link), suggest that hot Jupiters should have a single strong equatorial jet that advects heat from the hot dayside to the cool night side. The oklo splash image has been adapted from Cho’s calculations, and shows this jet in action (see this post for more discussion).

I’ve been advising UCSC Physics graduate student Jonathan Langton, who has recently begun a study of what the flow pattern on a hot Jupiter should look like if the planet is caught in Cassini state 2. If the planet’s rotation axis lies in the orbital plane, and if the planet spins on its axis once per orbit, then the play of light and shadow across the planetary orb has a pattern that is totally unlike our seasons here on Earth. At the north and south poles, of a spin-synchronous Cassini-state-2 planet, the parent star rises, passes directly overhead, and then sets once per orbit. At one special spot on the equator, on the other hand, the star is always visible, and additionally passes directly overhead once per orbit. At the opposite spot on the equator (which we’ll call the anti-stellar point), the star never fully rises, but rather peeks half of its diameter above opposite horizons once per orbit.

Jonathan has made two short .avi format animations that help to illustrate the situation. In the first animation, we hover above the point on the equator that receives maximum illumination. In the second animation, we hover above the point on the equator that receives the least illumination. The mythology on such a world would likely be pretty interesting.

When a spin-synchronous planet is illuminated in this bizarre manner, the flow pattern on its surface should be very different than the flow pattern that would occur if the planet is in the standard upright configuration. Jonathan has finished a preliminary set of simulations using the so-called shallow water approximation which indicate that this is indeed the case. (The shallow water approximation is a 2-dimensional method for simulating atmospheric dynamics on the surface of the planet under the assumption that the depth of the fluid is much smaller than the horizontal scales of interest. Use of this approximation doesn’t require us to assume that the red-hot Jupiter is actually covered with water!)

Here are two of Jonathan’s .avi format animations that show the (still very) preliminary results. The first animation [11 MB, modem users watch out!] shows the evolution of the temperature distribution (on the anti-stellar hemisphere) for a planet in Cassini state 2. Here’s a snapshot at a particular moment in time:

The second animation [12 MB] shows the distribution of vorticity across the planet surface. The vorticity at a particular spot in a fluid flow can be thought of as the ability of the flow to cause a tiny imaginary paddle-wheel to spin. Here’s a snapshot from the animation. The high-vorticity orange structure is a giant fiery hurricane-like storm on the surface of the planet: