Bali Hai

No word yet on whether that newly discovered 11 Earth-mass (and possibly rocky) planet orbiting GJ 674 is transiting or not.

The next opportunity is coming up on April 11th 13:17 UT. Given GJ 674’s location in the sky at RA 17:29, Dec -46 54, the South Pacific has by far the best view of the next event. Anyone willing to jump on the next plane to Tahiti with a Meade LX200 and an SBIG ST-7 in their checked baggage?

The current Tahitian weather forecast for the transit window calls for scattered clouds with a 20% chance of rain:

Not exactly the best conditions for obtaining 0.5% photometry, but not completely hopeless, either. I’m interpreting the current forecast as indicating there’s a 1/3rd chance that the weather will be cooperative. This means that if you fly to French Polynesia and set up your telescope in the hotel parking lot, you’ve got a 1 in 60 shot at walking away with the biggest exoplanet discovery of the year.

Even at 60:1 odds, there’s a case to be made that the trip is a good investment. According to the CoRoT website, the CoRoT satellite will detect “a few tens” of large rocky planets for a price tag of roughly 100 Million USD. That’s ~3 million per large rocky transiting planet.

A trip to Tahiti tomorrow, on the other hand, costs out at under 4K, and involves a more clement destination than Baikonur. In fact, when I dialed up a spur-of-the-moment expedition on expedia, I was informed that the price had just gone down:

The expectation value for the Tahiti mission, therefore, is a comparative bargain at $240,000 per transiting planet.

Assuming that you can show up at LAX by ~10pm this evening (Monday) a direct flight on Air Tahiti Nui gets you in to Papeete at 5:10 Tuesday morning. There’s plenty of time to grab a taxi to the luxe Le Meridien Tahiti, where you can take a refreshing nap in your “over water bungalow” set on one of Tahiti’s few sand beaches. Follow your late afternoon dip in the pool with dinner at Restaurant Le Carre, with its trendy atmosphere and refined A la carte dishes. After dinner, there’s still plenty of time for drinks at the L’Astrolabe Bar, where they’ll likely pick up your tab while you regale the hip-yet-distinguished clientele with astronomical bon mots. Indeed, you’ll likely have an admiring circle of new-found friends as you set up your scope in the parking lot and expertly obtain darks, flats, and baseline photometry, prior to observing well into astronomical twilight.

It’ll then be time to retire to your bungalow for some well-deserved rest. You’ll have the rest of the week to analyze your data and hopefully send that discovery e-mail to the IAU. It’ll be impossible for anyone on Earth to scoop your discovery until the next transit window on April 16th, at which point you’ll be flying home (having upgraded to first class for the long-haul flight back to LA).

What’s that you say? No money for your trip? No Problem. As soon as the market opens this morning, just write a few at-the-money April calls on a precariously high-flying tech stock to raise the necessary cash.

GJ 6-7-4

Image Source.

Word up! Chalk a jet-fresh Neptune on the boards — the Swiss’ve done it again.

The red dwarf GJ 674 lies 14.67 light years away. Minus 49 Dec. Only 53 known stars are closer to the Sun, and at V=9.382, GJ 674 is slightly more than twice as bright in the optical as its far more famous cousin GJ 876. With ~35% of the Sun’s mass, it’s packing more heat as well.

According to the Bonfils et al. discovery preprint posted to astro-ph yesterday, GJ 674 is accompanied by a sub-Neptune mass planet on a 4.6938 day orbit. Bucking the recent trends, the paper doesn’t contain a tabulation of the radial velocities. Eugenio, however, made dextrous use of the Dexter to scrape them off the figures, and they’re now safely packaged into the downloadable Systemic Console. The star has also been added to the “Real Stars” catalog on the Systemic Backend. The internal errors on the velocities are mostly below 1 m/s, which is impressive, given that each data point is based on a 15-minute integration of a rather dim star.

This discovery is a exciting for several reasons. Most immediate, is the fact that the planet does not yet seem to have been fully followed up photometrically to check for transits. At first glance, such an effort might appear to be hampered by the fact that the star is young enough to show significant photometric variability in synch with its 35-day rotation period. A central transit, however, would have a duration of only ~80 minutes — much shorter than starspot-induced variations — and would generate a clearly detectable dip of at least ~0.5% photometric depth.

Transitsearch.org has observers in Australia, South Africa, and South America, and so I’m hoping that they can quickly take advantage of this opportunity. The next transit window is centered about 15 hours from now, on April 06, 2007 at 20:38 UT. Here’s looking at you, Perth. The ephemeris table showing all the upcoming opportunities is at transitsearch.org. Based on a radius estimate for the star of 0.35 solar radii, the geometric transit probability is ~5.0%. Roll that twenty-sided die.

It’s fair to say that the next major discovery in the exoplanet game will likely be the detection of transits of a short-period Neptune-mass planet. Quite a few players are scrambling to be the first in the door. If it isn’t done from the ground during the next 6-months, then it’s likely that CoRoT will take the prize.

There’s a large difference in radius between sub-Neptune-mass planets made from rock and iron and sub-Neptunes composed mostly of water:

A Neptune transiting one of the brightest M-dwarfs in the sky would be a huge big deal. Hundreds of citations, Dude. Even if there’s no transit, this planet will likely be an excellent candidate for observation in the long-wavelength Spitzer bands, and fortunately there’s one more GO cycle before that cryogen runs out.

backlog

Image Source.

It was the end of the Winter quarter here at UCSC last week, and then I went on a trip, and then bam! More than a week with no posts… In the interim, there have been a number of interesting developments related to extrasolar planets. Here’s a brief run-down of some topics that I want to look at in more depth in the very near future:

Thanks to continuing efforts from the back-end user base, we’re accumulating a highly useful database of stable, low chi-square fits to the synthetic radial velocity data sets that comprise the Systemic Jr. catalog. Stefano has run a preliminary analysis and interpretation of the data. There are interesting implications for the overall eccentricity distribution of extrasolar planets, and there also appears to be a robust criterion for determining with confidence when you’ve extracted a real, previously unannounced planet from a given data set. We’re putting together a full report, which will appear quite soon. In the meantime, please keep submitting fits for systems that haven’t yet been adequately characterized.

The detection of another Neptune-mass planet orbiting a nearby red-dwarf was announced today. Yet more evidence for the core-accretion theory of planet formation! The discovery paper stops short of tabulating the radial velocities, but as I write this, Eugenio is busy dextering them onto the systemic back-end and onto the downloadable systemic console.

The theoretical case for the existence of Alpha Centauri B b is getting stronger by the day.

This year’s first ‘606 day is coming up next week, with a transit opportunity following on April 17th. I didn’t do enough to get the word out last December, but I’m hoping for good photometric coverage of the star during the upcoming window. A central transit for HD 80606b would last roughly 18 hours, so participation from observers worldwide will be required to definitively rule out transits.

Systemic Jr. Fit Drive

Image Source.

A big thank-you to everyone who’s been participating in the drive to characterize and study the catalog of synthetic “Systemic Jr.” planetary systems on the Systemic Backend. There’s now enough data to indicate that the analysis is going to be very informative. We’re looking forward to revealing the properties of the underlying planetary systems that were used to generate the data. In the meantime, we need your help to adequately characterize all 520 systems. Data in need of better characterization are marked by flags:

Our backend server is now swarming with various hard-working software robots that Stefano has assembled. The 100-year stability bot is rousted out of bed and set to work whenever a new fit is submitted. It reports a quick initial assessment of orbital stability. Planetary systems that pass through the 100-year stability screen are then put in a queue to wait for the attentions of the 1000-year stability bots. Systems that make it through 1000 years with less than a 1% change in semi-major axis of their planets are awarded a snazzy green flag:

Occasionally, systems that are in mean-motion resonance can show periodic semi-major axis variations of more than 1% while still remaining indefinitely stable. A resonance bot that will go through the fits and check for these special cases is currently being readied.

Systems that pass the minimum stability requirement are handed to a bootstrap bot which uses the bootstrap method to estimate uncertainties on the planetary orbital parameters for each stable fit. We’re currently running the bootstrap bot under the assumption that the orbits are pure Keplerian ellipses, and so the calculations are usually quite rapid. Very shortly, the error estimates for the parameters in submitted fits to the real systems and the Systemic Jr. systems will be showing up on the back-end data pages.

Finally, an “F-bot” has been activated which performs successive F-tests on submitted multiple-planet systems. Using its results, we’ll have a better idea of when the addition of a planet to a system is warranted.

HD 118206…

Image Source.

Most hot Jupiters have orbital eccentricities near zero because the tidal forces exerted on them by their parent stars are strong enough to rapidly circularize their orbits. Any planet whose orbit has been circularized should also be spin-synchronous, and so like our Moon with respect to the Earth, it should turn on its axis once every trip around the star. Synchronicity lends each hot Jupiter a permanent day and night side. This likely imparts a profound effect on both the planetary weather, and on the brightness of the planet when viewed in the infrared at different orbital phases.

All of the planets observed so far with the Spitzer Space telescope have nearly circular orbits, and hence are in (or are very near) the spin-synchronous state. We’re waiting to hear the results of our Spitzer GO-4 application to observe the highly eccentric planet HD 80606b, during an upcoming ‘606 day. If our observing proposal gets a thumbs-up, it’ll dramatically broaden the range of conditions under which planets have been observed. Very shortly, I’ll be posting the results of calculations that Jonathan Langton and I have been doing which predict what the light curve of HD 80606 should look like during the periastron passage in the various Spitzer bands. Here’s a sneak preview of how the temperature distribution on the planet might evolve over a 36-hour period as seen from a direction consistent with our line of sight from the Earth:

In looking over the latest officially published additions to the catalog of extrasolar planets, I noticed that there’s a very interesting object — HD 118203b — that straddles the extremes of the circular hot Jupiters and the ultra-eccentric HD 80606b. This planet was discovered in 2005 by the Swiss Team, has an orbital period of 6.13 days, a mass at least twice that of Jupiter, and a well-determined eccentricity, e=0.3. HD 118203b therefore won’t be spin-synchronous. Rather, as is also the case with HD 80606b (see the diagram here), it’ll have been forced into a state of pseudo-synchronous rotation, in which it does its best to keep one face toward the star during the periastron passage. Its day should be 64.8% as long as its year:

Higher resolution .eps version here.

Which raises a rather interesting question: Why is HD 118203’s eccentricity so high?

Assuming that the planet has a similar structure to Jupiter, the equations of tidal dissipation (see here for a discussion) indicate that the planet’s orbit should circularize in a mere 10-20 million years. This time scale is surprisingly short because the parent star is a subgiant with a radius ~1.5 times larger than the Sun. Something must be exerting a very strong perturbation to keep this planet’s e up.

In their discovery paper, Da Silva et al remark that the residuals around the best 1-planet Keplerian fit to the data are very large. It’s quite straightforward to verify this with the downloadable Systemic console (try it!) Da Silva and company were able to improve their fit by including a linear drift of 49.7 meters per second with their one-planet model. This corresponds to adding the effect of an outer planet that has been observed for only a small part of a single orbit. (The 43 published velocities span a period of 1.1 years.) They speculate that an outer as-yet-uncharacterized planet provides the gravitational perturbation that maintains the high eccentricity for the inner planet.

Last year, Fred Adams and I wrote a computer code (see these papers 1, 2) that includes the effect of general relativistic corrections on long-term planet-planet gravitational interactions. It’s easy to use this program to calculate what the long-term influences of various companion planets would have on HD 118203 b’s eccentricity. I ran a few trial cases, and quickly found that the interactions produced by companions that also provide the observed linear drift in the radial velocities don’t seem to be strong enough to explain HD 118203b’s high eccentricity. Could there be another explanation?

This is the sort of situation where the collaborative systemic back-end is extremely useful. I had a look at the stable fits that have been submitted so far for HD 118203. The best stable, self-consistent fit was uploaded back in October by the user Flanker, and has a reduced chi-square statistic of 1.96:

null.

This fit might point in an interesting direction for some further inquiry. Instead of using a linear trend to soak up the residuals to the one-planet fit, Flanker added two additional planets. One of the planets has a mass of 0.3 Jupiter masses and is orbiting with a period of 15 days. Its periapse is nearly aligned with the periapse of the inner planet. The resulting short-period secular interaction may well be strong enough to keep the eccentricity of the innermost planet high in the face of tidal dissipation. Flanker’s model also contains an outer planet with an orbital period of 244 days and a minimum mass 0.6 times that of Jupiter.

I think it’s worthwhile to explore additional models for this system that contain planets with short enough periods to intereact strongly with the 6.13 day planet. If the perturbing body has a relatively short-period orbit, then its presence will not be hard to verify with additional radial velocity observations of the star. And also, if Spitzer’s cryogen holds out, HD 118203b might be a very interesting target for a full-phase campaign.

Bootstrap

Image Source.

Stefano and Eugenio have been making quite a bit of development progress on the downloadable systemic console. A new version of the console (available for beta testing on the systemic backend) is now capable of providing an estimate of the uncertainties on the orbital elements associated with a fit to a particular data set.

Radial velocity data don’t provide an exact determination of planetary orbits. The most obvious shortcoming is that Keplerian orbital fits can’t determine the inclination of the planetary orbits, and so for a given system, we’re only able to measure then mass of the planet multiplied by the sine of the unknown inclination angle. Furthermore, the stellar radial velocity signal created by a planetary system is corrupted by astrophysical noise introduced by the parent star, as well as by noise introduced during the measurement process here on Earth.

Determination of the true uncertainties in a planetary orbital model is a subtle problem (for more detail, see Eric Ford‘s recent work in this area). As a first straightforward step, we’ve implemented the so-called “bootstrap” method of error estimation into the console. The bootstrap works by taking the original data set, and then successively redrawing time + velocity + uncertainties triples from the data with replacement. This procedure creates alternate realizations of the original data set in which some of the original measurements appear more than once, and in which some don’t appear at all. The best-fit parameters obtained by the console are then used as a starting guess to fit the bootstraped data sets. The standard deviations measured from the distributions of orbital elements thus obtained give error estimates for the parameters of the original fit.

The bootstrap routine is menu-accessed, and is simple to use. First, create a fit to a dataset. In the example just below, I’ve fitted to the data for HD 80606:

Once the fit has been polished, the bootstrap can be run. In the default configuration it uses Keplerian fitting and does 100 trials.

HD 80606 has been observed for nearly 20 orbital periods, and velocities have been obtained at a wide variety of orbital phases. As a result, the orbit is very well constrained. The bootstrap indicates that the uncertainty on the e=0.932 eccentricity is only 0.003. For other systems, such as hd 20782, which also seems to have a high eccentricity:

the uncertainty on the parameters is much larger:

Give the routine a try! In upcoming posts, we’ll talk more about how uncertainty estimates will be incorporated into the planetary catalogs on the backend.

Impact!

Image Source.

The news out of the Planetary Defense Conference is that NASA has the ability to locate all the potentially dangerous asteroids in our solar system by the year 2020, but that the cash to carry out the project is not currently forthcoming.

CNN and many other news outlets carried the story, along with a dramatic artist’s rendition of an asteroid striking the Earth.

That’s a rather large asteroid.

Measuring the circumscribed circles, it appears that the impacting body is 1/10th the Earth’s diameter, or approximately 600 km in radius. A bolide of this magnitude would currently rank as the ~25th largest object in the solar system, larger than Uranus’s moon Umbriel (584 km radius), smaller than Saturn’s moon Iapetus (736 km in radius), and nearly exactly the same size as Pluto’s moon Charon (606 km in radius). Ceres, the largest object in the main asteroid belt, by contrast, has a radius of ~475 km.

Based on the location of the late-afternoon catastrophe relative to the day-night terminator, the impactor seems to have had an orbit that was highly inclined relative to the solar system’s angular momentum plane. Perhaps it was undergoing Jupiter-driven Kozai oscillations prior to striking the Earth.

The last impact of the magnitude shown in the illustration was probably the Moon-forming impact ~4.4 Billion years ago, in which a Mars-sized body struck the Earth. Kevin Zahnle of the NASA Ames Research Center has estimated the distribution of giant impacts after the Moon-forming event. It’s likely that the largest strikes were by bodies with roughly half the radius of the object shown in the above picture.

The consequences of even a 300 km object hitting the Earth are severe. Such an impact is energetic enough to entirely vaporize the Earth’s oceans and create a temporary rock-vapor atmosphere with a surface pressure of ~100 bars. For a period of several months to a few years, the Earth would radiate with a temperature of ~2000 Kelvin — hot enough to glow brightly in the visible region of the spectrum. And depressingly, our planet would be fully sterilized by the event.

This week’s crop

Image Source.

The year 2007 is off to a reasonably good start. Three more planets were announced by the Geneva Planet Search team at a conference in Chile, bringing the total planet crop for ’07 up to seven.

The rate of planet discovery, however, has definitely leveled off. For the past four years, the detection rate has remained fixed at 26 new planets per year. The low-hanging fruit — the 51 Pegs, the 47 Ursae Majorii, the Upsilon Andromedaes — have all been harvested from the bright nearby stars, and increasingly extractive methods are being brought to bear. Transits are starting to contribute significantly to the overall detection rate. Radial velocity is pushing to planets with progressively lower masses. Surveys such as N2K are rapidly screening metal-rich stars that have high a-priori probabilities for harboring readily detectable planets. The neccessity of finding more planets is driving up the average metallicity of the known planet-bearing stars:

The three new planets, HD 100777b, HD 190647b, and HD 221287b are quite ordinary as far as extrasolar planets go. They all have masses somewhat greater than Jupiter, and they all take more than a year to orbit their parent stars. Their discovery seems not to have registered with the news media:

HD 100777 b, however, is actually deserving of some attention. Its orbital period of 383.7 days places it squarely in the habitable zone of its parent star. The eccentricity, e=0.36, is fairly high, and likely leads to interesting seasonal effects in the atmosphere of the planet.

HD 100777 b lies a regime where we expect to see white water clouds forming in the visible atmosphere. The planet is probably very reflective in the optical region of the spectrum (quite unlike the hot Jupiters, which are likely cloud-free, and which are known to absorb almost all of the starlight that strikes them). Convection of interior heat to the surface of HD 100777b is almost certainly driving collossal thunderstorms, and the atmospheric disturbances created by the thunderstorms likely feed giant vortical storms similar to Jupiter’s great red spot.

It’s also possible that the atmosphere is much clearer in regions where air wrung dry by rainfall is downwelling. This phenomenon occurs on Jupiter, where highly transparent patches occur over several percent of the Jovian surface:

Image Source.

The Galileo entry probe went right into one of these regions, and sampled very dry air. On HD 100777, the regions of high atmospheric transparency will probably preferentially absorb red and green light (as a result of Rayleigh scattering of incoming photons). The surface, then, in the vicinity of a downwelling region may look something like this:

flyby

Image: NASA New Horizons Spacecraft (false color by oklo).

One day, one hour, and nine minutes ago, the New Horizons spacecraft sailed flawlessly through its closest approach to Jupiter. A day later, Jupiter still looms large in New Horizon’s field of view, with an angular size more than five times greater than the size of the full moon in our sky.

Jupiter, during its 4.5 billion year history, has been visited by at least seven other probes. Pioneer 10, Pioneer 11, Voyager 1, Voyager 2, Ulysses, Galileo, and Cassini have all successfully made the journey. This latest encounter was buried beneath the news of a 500-point drop in the Dow Jones Industrial Average. The flyby, in fact, hasn’t even made it onto the Astronomy Picture of the Day!

A decade ago, many of the metal atoms in the New Horizons spacecraft were still buried in the Earth’s crust. A bit more than a year ago, the assembled spacecraft was flown, in a sealed pressurized container, to Cape Canaveral for launch. All through the past several weeks, it’s been taking pictures of the Jovian system. Most of the data will be radioed back to Earth over the coming months. The image above was taken on Monday, and shows a Von Karman vortex sheet trailing away from the Little Red Spot, currently the second-largest storm in the Solar System.

In a sense, the Jupiter encounter was mostly utilitarian. It boosted the spacecraft’s heliocentric velocity (at the expense of Jupiter’s orbital energy) and cut down the travel time to Pluto.

The next scheduled mission to Jupiter is Juno, the Jupiter Polar Orbiter, which is scheduled to arrive at the Jovian system in 2016.

spectra

Image Source.

So, uhh, yeah, the oklo blog went through a dry spell with no posts last week. This was primarily a consequence of the fact that Spitzer GO-4 proposals were due last Friday. I teamed up with Drake Deming of GSFC and UCSC physics grad student Jonathan Langton to propose a 30-hour observation of HD 80606b during the ‘606 day that’ll occur next November 20th. In an upcoming article, I’ll be pushing the reasons why we’re really excited about the possibility of observing HD 80606 b during its big periastron swing.

The Spitzer Space Telescope has turned out to be a regular wellspring of exo-planet results. It’s providing very interesting and often surprising constraints on the weather conditions at the surfaces of the hot Jupiters, and another big new result was announced today. Three different teams released the first-ever observations of emergent infrared spectra from two observational runs on HD 189733 and HD 209458.

The transiting planet HD 189733 b was discovered by the Swiss team in 2005. Of the fourteen known transiting planets, HD 189733 b is the best-suited for detailed follow-up observations. The parent star lies only 19 parsecs away, the orbital period is a skimpy 2.1 days, and the 1.15 Jupiter-mass planet has a radius fully 15% the size of the primary star’s radius. Like the other transiting systems, the planet, the orbit, and the star can all be drawn completely to scale on a “saved for web” diagram that’s only 420 pixels across:


Grillmair et al.’s Spitzer spectrum of HD 189733 was obtained with 12 hours of observation, in which the brightness of the star at infrared wavelengths between 7 and 14 microns is compared in and out of the secondary transit:

The observed flux distribution from the planet is nearly completely flat as a function of wavelength! Models of the atmospheres of hot Jupiters had all predicted that the presence of water vapor in the planetary atmosphere would lead to a prominent absorption feature at ~8 microns. No hint of the predicted dip was seen. The overall amount of infrared light coming from the planet during the secondary transit indicates that heat is probably being efficiently redistributed between the day and night sides of the planet.

A Nature paper by Jeremy Richardson and collaborators reported a very similar result for HD 209458 b. Their spectrum runs between 7.5 and 13.2 microns, is similarly devoid of absorption features, and also suggests a modest day-night temperature difference.

So how to interpret these results? One possibility is that the lack of absorption lines is caused by a high, uniformly emitting cloud layer, perhaps made of silicate grains. A problem with this interpretation, however, is that the cloud decks would have to be extremely dark and unreflective in the optical. Hot Jupiters absorb nearly 95% of the radiation that they receive from their parent stars. Another possibility, put forward by Jonathan Fortney, is that the atmospheres of these planets are isothermal down to large optical depths. Because we can’t actually see to a hotter underlying layer, there’s no mechanism for deep absorption lines to form.

And finally, another, rather startling, interpretation of the results was offered several hours ago by CBS News:

Image Source.