Cycle 18


It’s been rather arduous past few days as the HST Cycle-18 proposal deadline — 5 PM PST, Friday Feb. 26th, to be exact — bore down like a freight train.

During the past year, I’ve become quite intrigued by the remarkable (and well known) HST observation by Vidal Madjar et al. (2003), who discovered that the Lyman-alpha transit depth of HD 209458b is a whopping 15% (as opposed to the mere 1.5% of the star’s light that gets blocked during the optical transit). The implication of this result is that the planet is surrounded by an outflowing, escaping wind of hydrogen, and the discovery has sparked a lot of theoretical work.

A good test for planetary outflow hypotheses is to see what they predict for eccentric planets that undergo drastic changes in stellar heating during the course of the orbit. Fred Adams and I have been working on hydrodynamical models for these situations, and it soon became clear that oklo.org fave HD 17156b, the P=21.2d, e=0.67 transiting planet provides an intriguing observational opportunity for HST/STIS. HD 17156 is currently the fourth-brightest known parent star of a transiting extrasolar planet (after HD 209458, HD 189733 and HD 149026) and it lies in HST’s so-called continuous viewing zone for part of the year. This means that a full transit can be observed without having to take those leisurely once-per-96-minute pit stops every time Earth blocks the view.

The geometry of the transit, furthermore, is such that the planet is getting its maximum sunburn just a few hours after transit egress. Our calculations indicate that it should take the upper atmosphere only a matter of hours to react to the increased heating, so we’re optimistic about the possibility that not only will HST detect a deep transit in the UV, but that it might even be able to detect the Lyman-alpha transit depth increasing during the course of the transit. Here’s the basic idea:

As of a few minutes ago, the proposal was received safe and sound at STScI, so now it’s time to kick back, wait, and see if it passes muster with the TAC…

upcoming event

I’ve got an upcoming event planned in New York City that should be pretty interesting. From the UCSC Newsletter:

UCSC astronomer joins composer Philip Glass to explore music of the universe

UC Santa Cruz Astronomer Gregory Laughlin joins acclaimed composer Philip Glass February 21 in a “Brainwave” discussion at the Rubin Art Museum in New York.

For its third year and in conjunction with the exhibition Visions of the Cosmos, Brainwave is a series of 20 sessions this winter and spring that bring together eminent thinkers from multiple disciplines with neuroscientists and astrophysicists to ponder big thoughts about “things that matter.”

Laughlin and Glass appear in the third Brainwave event titled “How Do We Listen to the Music of the Spheres?”

Laughlin is a professor of astronomy and astrophysics whose research delves into orbital dynamics and the evolution of planetary systems. Glass is one of the most influential composers of the past half-century. Though sometimes called a “minimalist,” Glass describes his compositions as “music with repetitive structures.”

Laughlin said he and Glass will explore commonalities between music and orbital dynamics. The museum’s initiative to pair the two was sparked in part by Laughlin’s articles on his blog oklo.org that delve into ways to “sonify” planetary movements.

He developed software to map planetary systems as audible waveforms. He said he became intrigued by the realization that planetary systems can be used as a type of nonlinear digital synthesizer and can provide an enormous palette of sound — sounds never before heard.

The Laughlin/Glass Brainwave session begins at 6 p.m. Sunday, February 21 at the Rubin Museum of Art at 150 West 17 St., New York City. Admission is $25.

Over the next week, as I’m preparing for the event, I’ll be working extensively with the sonification capability of the systemic console. Just below, is a reprinted post that touches on this very cool, and still relatively unexplored feature. If you’ve worked with the Console’s N-body sonification, and if you’ve found interesting results, feel free to send me .fit files — an extraordinarily effective form of compression(!) — and I may be able to use them in the discussion.

Potentially the most interesting feature on the downloadable systemic console is the “sonify button”, which integrates the model planetary system specified by the state of the console sliders and produces a .wav format CD-quality audio file of the resulting radial velocity waveform. Not interested in planets? The console is a stand-alone non-linear digital synthesizer. It’s capable of producing strange, remarkable, musically useful sounds. They merely need to be located within the uncountable infinity of solutions to the gravitational N-body problem.

First, use the console to build an interesting multi-planet system (for this purpose, there’s no need to try to fit whatever data is in the window.) Then click the sonify button. This brings up a dialogue window which enables the user to make several specifications for the sound file that is produced.

console sonify dialogue box

The most important user-specified parameter is the frequency onto which the orbital period of the shortest-period planet on the console is mapped. If, for example, the innermost planet has a period of 365.25 days, then a 440 Hz map will play 440 years worth of evolution in one second. (440 Hz corresponds to the A below middle C.) Mapping the radial velocity curve onto a high-frequency note extends the total number of orbits that go into the sample, and thus increases the integration time required to produce the sample. You can also specify the length of the sample, and you can exert simple control over the attack and decay rate of the envelope for the overall waveform.

Once you’ve produced the sound file, it appears in the “soundClips” subdirectory within the systemic parent directory. Both of these directories are automatically created when you download and expand the console — see the instruction set for the downloadable console for more details. With a Macintosh, you get the best results if you play the sample right from the folder. i-Tunes seems to want to convert the samples to .mp3 format in a manner that introduces audible noise, and we’re not yet sure how to resolve this issue.

To the extent that planets orbit independently of one another, the console behaves like a simple additive synthesizer, in which the individual Kepler waveforms add to form a composite sound. Much more interesting, is the situation when planets experience significant gravitational interaction, leading, for example, to resonance and to nonlinear instability (here are examples, 1, 2, from the resources page of both types of waveforms). Close encounters provide discontinuities between individual blocks of sound that resemble the results of granular synthesis.

The strongest 2-planet mean-motion resonances occur when the pair of planets share a common period and engage in a one-to-one resonant motion. There are a variety of different one-to-one resonances, including binary planet orbits (e.g. Earth and Moon), trojan configurations, and generalizations of retrograde satellite orbits. In this last category, one can have two planets with the same semi-major axis, but with different eccentricities. If one starts the planets in the following configuration, then the motion is dynamically stable, and evolves in a complicated way over time.

evolution of eccentric 1:1 resonance

The motion leads to an interesting audio wave-form, in which you can hear the system cycling between configurations in which both planets are modestly eccentric and configurations in which one orbit is nearly circular while the other one is highly eccentric. As a specific example, set the console to the following configuration: P1=P2=10 days, M1=M2=0.3 Mjup, MA1=180., MA2=190., e1=0.9, e2=0.1, long1=0.0, long2=0.0. If you increase MA2 to about 225 degrees while keeping the other parameters fixed, you’ll hear the system go unstable.

Evolving, high-eccentricity orbits tend to have an insect-like quality, which brings to mind the 1986 album, The Insect Musicians, by Greame Revell (formerly of SPK). From the album jacket:

For the two years 1984-85, Graeme Revell traveled from Australia to Europe, to Africa, Indonesia and North America recording and negotiating copyrights of insect sound recordings. It took another full year sampling and metamorphosing some forty sounds thus gathered using the Fairlight Computer Musical Instrument, to produce this record. The only sounds used are those of insects, altered digitally and combined into a unique orchestra of instruments, an orchestra of strange and delicate timbres, music of natural rhythm and texture.

limit cycle

The minimum threshold level for amazement will rise quickly once Kepler’s discoveries start to accumulate, and already, it’s getting very hard to remember which transiting planet is unusual for which reason. Let’s see, was it TrES-4 or WASP-17 that had that styrofoam-like density? Or was it both of them?

Even in a crowded field, though, HAT-P-13 is likely to endure as a touchstone. In the next five years, it’s likely that there will emerge only a select handful of systems in which a well-characterized transiting planet orbiting a relatively bright star is being substantially perturbed by a companion on a well-characterized orbit:

After the HAT-13 system was announced, we showed that the planets “b” and “c” should have evolved to an eccentricity fixed point configuration, in which the orbits’ apsidal lines co-rotate, and in which the orbital eccentricity of planet “b” has a very sensitive dependence on its internal structure. Further modeling, using reasonable assumptions, gives strong limits on the tidal Q of planet “b”. In essence, one can potentially accomplish with an exoplanet a big chunk of what the Juno Mission expects to accomplish at Jupiter at of order a thousandth of the cost.

Our analysis assumed that HAT-P-13 b and c are on co-planar orbits. There’s an interesting new paper by Rosemary Mardling that explores the significantly more complex situation that arises if the orbital planes of the planets are significantly misaligned. In this case, tidal dissipation in the inner planet causes the system to settle into a limit cycle, where the eccentricity and the angle between the apsides circulate on a secular timescale, and the easy insight into the structure of planet b is no longer possible.

Interestingly, however, Mardling’s analysis suggests that if the orbits are misaligned, then the mutual inclination is likely to be in the neighborhood of 45 or 50 degrees. A mutual inclination of, say, 30 degrees is inconsistent with the currently observed parameters of planet b. The following two diagrams (figure 8 a and b) from her paper show how this works:

Within the next few months, we should get improved values for the eccentricity and radius of planet b, which will significantly shrink the size of the peach-colored boxes in the two figures above. HAT-P-13c is also currently headed in for periastron, with the next transit opportunity scheduled for April 12, 2010. A transit by planet c would provide strong evidence that the system is reasonably close to co-planar (and would be quite remarkable in its own right!) Furthermore, during the periastron passage of c, there should be readily detectable transit timing variations for b, which should give us a shot at distinguishing between the co-planar case and the case with a mutual eccentricity of 45-50 degrees. In the next post, I’ll look in detail at the numbers…