systemic

Worlds worlds worlds5
greg posted in exoplanet detection on May 25th, 2008

On Friday, I flew back from the Boston IAU meeting, still buzzing with excitement. On Saturday, I woke up with what might best be described as a transit-induced hangover (an entirely distinct condition from transit fever). I’d basically allowed all my professorial responsibilities to slide for a week. On my desk is a mountain of work, a preliminary exam to assemble, and a horrifying backlog of e-mail.

Ahh, but like an exotic sports car bought on credit, it was worth it. The meeting was amazing, certainly the most exciting conference that I’ve ever attended. Big ups to the organizers! Planetary transits are no longer the big deal of the future. They’re the big deal of the right here right now. Spitzer, Epoxi, MOST, HST and CoRoT are firing on all cylinders. The ground-based surveys are delivering bizarre worlds by the dozen. And we’re clearly in the midst of very rapid improvement of our understanding of the atmospheres and interiors of the planets that are being discovered.

From a long-term perspective, the conference’s biggest news was probably provided by the Geneva group, in the form of Christophe Lovis’ presentation on Tuesday afternoon. In his 15-minute talk to a packed auditorium, Lovis covered a lot of ground. I scrambled to take notes. My reconstructed summary (hopefully without major errors) runs like this:

The HARPS planet survey of solar-type stars contains ~400 non-active, slowly rotating FGK dwarfs. Observations with the 3.6-meter telescope have been ongoing since 2004, and over time, their emphasis has been progressively narrowed to focus on stars that harbor low-amplitude radial velocity variations with RMS residuals in the 0.5-2.0 m/s range. The current observing strategy is to obtain a nightly multiple-shot composite velocity of an in-play candidate during block campaigns that run for 7-10 nights.

During the first few minutes, Lovis reviewed the current status of the published results. The Mu Arae planets (including the hot Neptune on the 9.6-day orbit, see here and here) are all present and accounted for. The HD 69830 triple-Neptune data set (see here, here and here) now contains twice as many velocities, with virtually no changes to the masses and orbits of the three known planets. Long-term scatter in the HD 69830 data set is at the ~90 cm/sec level, indicating either the effect of residual stellar jitter, or perhaps the presence of additional as-yet uncharacterized bodies.

He then announced that there are currently forty-five additional candidate planets with Msin(i)<30 Earth masses, P<50 days and acceptable orbital solutions. And that’s not counting candidates orbiting red dwarfs.

He then began to highlight specific systems. To say that planets were flying thick and fast is an understatement. Here’s the verbatim text that I managed to type out while simultaneously attempting to focus on the talk:

Rumor has it that some of these systems will be officially unveiled at the upcoming Nantes meeting on Super Earths. Odds-on, with 45 candidates in play, we’ll soon be hearing about a transiting planet with a mass of order ten times Earth’s. I won’t be at the Nantes meeting, but the stands will be harboring agents of the Oklo Corporation.

The talk finished with an overview of the statistics of the warm Neptune population. Most strikingly, a full 80% of the candidates appear to belong to multiple planet systems, but cases of low-order mean motion resonance seem to be rare [as predicted –Ed.] . There is a concentration of these planets near the 10-day orbital period, and the mass function is growing toward lower masses. Significant eccentricities seem to be the rule. And finally, I think it was mentioned that the planet-metallicity correlation is weaker for the warm Neptunes than for the population of higher-mass planets.

Seems like core accretion is standing the test of time.

Note on the images: Gaspar Bakos (of HAT fame) had the cool idea of machining metal models for the planets of known radius which are correct in terms of relative size, and which have the actual density of their namesakes. HAT-P-6, for example, is constructed from a hollow aluminum shell, and with a density of ~0.6 gm/cc it would float like a boat. HAT-P-2b, on the other hand, which packs 8.6 Jupiter masses into less than a Jovian radius, has the density of lead and (not coincidently) is made out of lead. It’s startling to pick it up. CoRoT-Exo-3b, which was announced at the meeting, has a mass of twenty Jovian masses, and a radius just less than Jupiter. I guess that one will have to be made from Osmium.

Earth, at ~5.5 gm/cc, on the other hand, can be readily manufactured from a variety of different alloys.

A Field Guide to the Spitzer Observations10
greg posted in exoplanet detection on May 18th, 2008

Jonathan Fortney has the office next to mine at UCSC, and so we’re always talking about the Spitzer observations of extrasolar planets. The Spitzer Space Telescope has proved to be an extraordinary platform for observing planets in the near infrared, and during the past year, the number of published and planned observations has really been growing rapidly.

Increasingly, with the flood of data, I’ve been finding that I have trouble keeping mental track of all the photometric observations of all the planets that Spitzer has produced. Let’s see, was Tres-1 observed in primary eclipse? Did someone get a 24-micron time series for HD 149026? And so on.

So Jonathan and I decided to put together a poster that aggregates the observations (that we know of) that have either been completed, or which have been scheduled. The relevant information for each campaign includes the star-planet system, the bandpass, and the duration and phase of the observation. We wanted the information for each system to be presented in a consistent manner, in which the orbits, the stars, and the planets are all shown to scale (and at a uniform scale from system to system). As an example, here’s the diagram for HD 189733:

In putting the poster together, we were struck by the variety of different observational programs that have been carried out. Some of the diagrams, furthermore, with text removed, have a delicate insect-like quality.

(The figure just above shows Bryce Croll’s planned 8-micron observations of Transitsearch.org fave HD 17156b. Croll’s campaign will attempt to measure the pseudo-synchronous rotation period of the planet.)

I’m going to Boston next week to attend the IAU transit meeting, and so I printed out a copy of the poster to put up at the meeting:

Here’s a link to the Illustrator file and the .pdf version. Full size, it’s two feet wide and three feet tall. Going forward, I’ll update the files as new observations come in.

Just like in 18465
greg posted in exoplanet detection on April 28th, 2008

Uranus and Neptune have returned to nearly the configuration that they were in at the time of Neptune’s discovery in 1846. Using Solar System Live, it’s easy to see where the planets were located when Galle and d’ Arrest turned the Berlin Observatory’s 9-inch Fraunhofer refractor to the star fields of the ecliptic near right ascension 22 hours:

In 2011, Neptune, with its 165-year period period, will have made one full orbit since its discovery. Uranus, with an 84-year period, will have gone around the Sun almost two times.

Because the planets are fairly close to conjunction, Neptune has recently gone through the phase of its orbit where it exerts its largest perturbation on the motion of Uranus. This was similarly true in the years running up to 1846, and was responsible for LeVerrier’s sky predictions bearing such a stunning proximity to the spot where Neptune was actually discovered by Galle.

LeVerrier (and Adams) were quite fortunate. Without a computer, multi-parameter minimization is hard, and both astronomers cut down on their computational burden by assuming an incorrect distance for Neptune (based on Bode’s “law”). Their solutions were able to compensate for this incorrect assumption by invoking masses for Neptune that were much too large. They carried out remarkable calculations, but nevertheless, luck (in form of the fact that Uranus and Neptune had recently been near conjunction) played a considerable role.

Predictably, as soon as the real orbit of Neptune was determined, the playa haters tried to rush the stage. Benjamin Peirce of Harvard, in the Proceedings of the American Academy of Arts and Sciences 1, 65 (1847) described LeVerrier’s accomplishment as a mere “happy accident”:

I personally think that’s going a bit far. In any case, it’s interesting to compare the two independent predictions with the actual orbit of Neptune. I pulled the LeVerrier and Adams data in the following table from Baum and Sheehan’s book “In Search of Planet Vulcan” :

There’s been no shortage of hard work, and there’s been no shortage of predictions and false alarms, but nevertheless, nobody has managed to discover another solar system planet via analysis of gravitational perturbations. With the extrasolar planets, however, the prospects look a lot better. In particular, the Systemic Backend collaboration can team up with amateur observers to do the trick.

On the Systemic Backend, there are many candidate planets that have had their orbits characterized. As is usually the case with planet predictions, most of the candidates will wind up being spurious, but it’s definitely true that real planets orbiting real stars have been detected by the Backend user base. For example, Gliese 581 c was accurately characterized by the Systemic users several months before it’s announcement by the Swiss (see this post) and the same holds true for 55 Cancri f (see this post).

In the happy circumstance that a candidate planet is part of a system with a known transiting planet, then there’s an increased probability that if the candidate planet exists then it can also be observed in transit. This provides a channel for detection that completely circumvents the need for professional astronomers to carry out confirming radial velocity observations. Amateur observers are currently pushing the envelope down to milli-mag precision. Here’s an out-of-transit observation of the parent star of XO-1b by Bruce Gary:

This photometry is potentially good enough to confirm a Neptune-sized planet in transit across a Solar-type star, which is absolutely amazing.

An initial proof-of-concept observation has recently been carried out. On the systemic backend, the users have been investigating the HD 17156 system, which contains a known transiting planet. User “japf ” (José Fernandes) found that a lower chi-square fit to the published radial velocity data can be obtained if there’s a 6.2 Earth-mass companion on a 1.23 day orbit.

The best-fit eccentricity of the planet would bring it to a hair-raising 2 stellar radii of HD 17156, and if the planet is made of rock or water, it’ll be too small to detect, but nevertheless, it’s at least worth having a look. Jose sent the ephemeris to Bruce Gary, who observed on the opportunity falling on April 20, 04.5 UT.

No transit detected. This in itself was not at all surprising, given the long-shot nature of this particular candidate planet. What’s exciting, though, is that the full pipeline is now in place. There will definitely be strong candidates emerging over the coming months, and I think it’s quite probable that we’ll see a prediction-confirmation that is at least as good a match as was obtained for Neptune in 1846…

first quarter numbers3
greg posted in exoplanet detection on April 2nd, 2008

Image Source.

Back in 2002, Keith Horne gave a talk at the Frontiers in Research on Extrasolar Planets meeting at the Carnegie Institute in Washington and showed an interesting table:

At that time, there were more than two dozen active searches for transiting extrasolar planets, but only a single transiting planet — HD 209458 b — had been detected. Transits were generating a lot of excitement, but paradoxically, the community was well into its third straight year with no transit detections. The photometric surveys seemed to be just on the verge of really opening the floodgates, with a total theoretical capacity to discover ~200 planets per month.

It’s been six years, and the total transiting planet count is nowhere near 14,000. Most of the surveys on the table have had a tougher-than-expected time with detections because of the large number of false positives, and because of the need to obtain high-precision radial velocities on large telescopes to confirm candidate transiting planets. Indeed, the surveys that were sensitive to dimmer stars have largely faded out. It’s just too expensive to get high-precision velocities for V>15 stars. With the exception of the OGLE survey (which had been set up to look for microlensing during the 1990s, and which had established a robust pipeline early on) none of the surveys that employed telescopes with apertures larger than 12 cm have been successful. The currently productive photometric projects: TrES, XO, HATnet, and SuperWASP all rely on telescopes of 10 to 11 cm aperture to monitor tens of thousands to hundreds of thousands of stars, and all are sensitive to planets transiting stars in the V~10 to V~12 magnitude range. This magnitude range is the sweet spot: there are plenty of stars (and hence plenty of transits) and the stars are bright enough for reasonably efficient radial velocity confirmation.

Yesterday, SuperWASP rolled out 10 new transits at once, dramatic evidence of the trend toward planetary commoditization and of the fact that it’s getting tougher to make a living out on the discovery side. The detection of new planets is growing routine enough that in order to generate a news splash, you need multiple planets, and the more the better. This inflationary situation for new transit news is highly reminiscent of where the Doppler surveys were at seven years ago. For example, on April 4, 2001, the Geneva team put out a press release announcing the discovery of eleven new planets (including current oklo fave HD 80606b).

I’d like to register some annoyance with this latest SuperWASP announcement. There are no coordinates for the new planets, making it impossible to confirm the transits. There is no refereed paper. The data on the website are inconsistent, making it hard to know what’s actually getting announced. I was astonished, for example, that WASP-6 is reported on the website to have a radius 50% that of Jupiter, and a mass of 1.3 Jovian masses:

That’s nuts! If the planet is so small, why is the transit so deep? And a 2200 K surface temperature for a 3.36d planet orbiting a G8 dwarf? Strange. Perhaps the radius and mass have been reversed? In addition, there are weird inconsistencies between the numbers quoted in the media diagram and in the tables. For example, the diagram pegs WASP-7 at 0.67 Jovian masses, whereas the table lists it at 0.86 Jovian masses. WASP-10 has a period of 5.44 days in the table and 3.093 days in the summary diagram. Putting out a press release without the support a refereed paper is never a very good idea, even when there’s a danger that another team will steal your thunder with an even larger batch of planets.

Despite the difficulty in getting accurate quotes from the exchange, it’s interesting to see how the ten new planets stack up in the transit pricing formula. Using the data from the new WASP diagram (except for the 0.66 day period listed for WASP-9) and retaining the assumption that USD 25M has been spent in aggregate on ground-based transit searches, the 46 reported transits come out with the following valuations:

The ten new WASP planets (assuming that the correct parameters have been used) contribute about 1/10th of the total catalog value. There will likely be interesting follow-up opportunities on these worlds from ground and from space, but its unlikely that they’ll rewrite the book on our overall understanding of the field.

It’s interesting to plot the detection rate via transits in comparison to the overall detection rate of extrasolar planets. (The data for the next plot was obtained using the histogram generators at the Extrasolar Planets Encyclopaedia, which are very useful and are always up-to-date.)

It’s a reasonable guess that 2008 will be the first year in which the majority of discoveries arrive via the transit channel, especially if CoRoT comes through with a big crop. Radial velocity, however holds an edge in that it’s surveying the brightest stars, and (so far) has been responsible for progress toward the terrestrial-mass regime. I think that we might be seeing planets of only a few Earth masses coming out of the RV surveys during the coming year. Certainly, everything else being equal, a planet orbiting an 8th magnitude star is far more useful for follow-up characterization than a planet orbiting a 13th magnitude star.

1:1 eccentric2
greg posted in exoplanet detection on March 30th, 2008

Image Source.

The range of planetary orbits that are observed in the wild is quite a bit more varied than the staid e < 0.20 near-ellipses in our own solar system. For regular oklo readers, the mere mention of Gl 876, 55 Cancri, or HD 80606, is enough to bring to mind exotic worlds on exotic orbits.

Non-conventional configurations involving trojan planets have been getting some attention recently from the cognescenti. Even hipper, however, is a configuration that I’ll call the 1:1 eccentric resonance. Two planets initially have orbits with the same semi-major axis, but with very different eccentricities. Conjunctions initially occur close to the moment of apoastron and periastron for the eccentric member of the pair.

Here’s a movie (624 kB Mpeg) of two Jupiter-mass planets participating in this dynamical configuration.

At first glance, the system doesn’t look like it’ll last very long. Remarkably, however, it’s completely stable. Over the course of a 400-year cycle, the two planets trade their angular momentum deficit back and forth like a hot potato and manage to orbit endlessly without anyone getting hurt.

Here’s an animation (1 MB Mpeg) which shows a full secular cycle. The red and the blue dots show the planet positions during the two orbit crossings per orbit made by one of the planets. It’s utterly bizarre.

These animations were made several years ago by UCSC grad student Greg Novak (who’ll be getting his PhD this coming summer with a thesis on numerical simulations of galaxy formation and evolution). As soon as we can get the time, Greg and I are planning to finish up a long-dormant paper that explores the 1:1 eccentric resonance in detail. In short, these configurations might be more than just a curiosity. When planetary systems having three or more planets go unstable, two of the survivors can sometimes find themselves caught in the 1:1 eccentric resonance. The radial velocity signature of the resulting configuration is eminently detectable if the planets can be observed over a significant number of orbital periods.

one seven one five six redux3
greg posted in exoplanet detection on March 22nd, 2008

Image Source.

Stefano, and Eugenio and I have been completely immersed in several time-critical projects during the past few months, and as a result, the frequency of posts here on oklo.org has not been as high as I would like. We’re starting to see our way clear, however, and very shortly, there’ll be a number of significant developments to report. Also in the cards is a major new release of the console, and a refocus on the research being carried out on the systemic backend. In any case, sincere thanks to all the backend participants for their patience.

Oklo regulars will recall all the excitement last fall surrounding the discovery of transits by HD 17156b. The transit was first observed on September 10th by a cadre of small telescope observers, and was then confirmed 21.21 days later on October 1.

Jonathan Irwin at Harvard CfA has led the effort to analyze and publish the October 1 observations of the transit. The work recently cleared the peer-review process, and was posted on the web a few days ago. (Here’s a link to the paper on astro-ph.)

The night of October 1 was plagued by atrociously aphotometric conditions across the North American continent, and most of the observers who tried to catch the transit were clouded out. Southern California, however, had reasonably clear skies, and three confirming time series came from the Golden State. The Mount Laguna observations were taken from SDSU’s Observatory in the mountains east of San Diego, the Las Cumbres observations were made from the parking lot of the LCOGT headquarters in Santa Barbara, and Transitsearch.org participant Don Davis got his photometry from his backyard in suburban Los Angeles.

The aggregate of data from the October 1 transit allowed us to refine the orbital properties of the planet, and additional confirming observations in a paper by Gillon (of ‘436 fame) et al have given a much better characterization of the orbit.

Because of the high orbital eccentricity, the planet should have very interesting weather dynamics on its surface. Jonathan Langton’s model predicts that the planet’s 8-micron flux should peak strongly during the day or so following periastron passage as the heated hemisphere of the planet turns toward Earth.

By measuring the rise and subsequent decay of the planet’s infrared emission, it’ll be possible to get both a measure of the effective radiative time constant in the atmosphere as well as direct information regarding the planet’s rotation rate. Bryce Croll is leading a team that successfully obtained time on the Spitzer telescope to make the observations.

In another interesting development, a paper by Short et al. appeared on astro-ph last week which proposes the existence of a second planet in the HD 17156 system. The Short et al. planet has an Msin(i) of 0.06 Jupiter masses and an orbital period of 111.3 days. It’s quite similar to the slightly more eccentric (and hence dynamically unstable) version of the HD 17156 system proposed by Andy on the Systemic Backend last December, which was based on the radial velocities and transit timing then available:

The existence of a second planet in the HD 17156 system would be extremely interesting! The immediate question, however, is, how likely is it that the second planet is actually there?

To make an independent investigation, it’s straightforward to use the downloadable systemic console to fit to the available published data on HD 17156. I encourage you to fire up a console and follow along. Now that the Irwin et al. paper is on the web, we have the following transit ephemerides:

These can be added to the HD17156.tds transit timing file in the datafiles directory. The file should be edited to look like this:

When the HD17156v2TD system is opened on the console, it shows both the radial velocity and the transit timing data.

It’s quick work to dial in a one planet fit to the RV and transit timing data. I get a system with the following fit statistics:

The required jitter of 2.12 m/s indicates that a one planet fit to the data should still be perfectly adequate, since the star (which is fairly hot and massive) has an expected stellar jitter of order 3 m/s. Nevertheless, the residuals periodogram does show a distinct peak at ~110 days:

Using the 110 day frequency as a starting point, one finds that ~0.1 Mjup planets do indeed lower the chi-square. I’ve uploaded an example two planet fit to the systemic backend that harbors a second planet in a 113 day orbit and a mass of 0.13 Jupiter masses. Its periastron is aligned with that of planet b, and the RMS has dropped down to 3.08 m/s (for a self-consistent, integrated fit). The implied stellar jitter is a bargain-basement 0.59 m/s, which is almost certainly too good to be true.

When I do an F-test between my one and two planet fits, the false alarm probability for planet ccomes in at 38%. It’s thus fairly likely that the second planet is spurious, but nevertheless, it certainly could be there, and it’ll be very interesting to keep tabs on both the transit timing data and the future radial velocity observations of this very interesting system…

Toward Alpha Cen B b24
greg posted in exoplanet detection, worlds on February 24th, 2008

Image Source.

Yesterday, I gave a talk at the JPL Exoplanet Science and Technology Fair, a one-day meeting that showcased the remarkably broad variety of extrasolar planet-related research being carried out at JPL. In keeping with the wide array of projects, the agenda was fast-paced and completely diverse, with talks on theory, observation, instrumentation, and mission planning.

The moment I walked into the auditorium, I was struck by the out-there title on one of the posters: The Ultimate Project: 500 Years Until Phase E, from Sven Grenander and Steve Kilston. Their poster (pdf version here) gives a thumbnail sketch of how a bona-fide journey to a nearby habitable planet might be accomplished. The audacious basic stats include: 1 million travelers, 100 million ton vessel, USD 50 trillion, and a launch date of 2500 CE.

Fifty trillion dollars, which is roughly equivalent to one year of the World GDP, seems surprisingly, perhaps even alarmingly cheap. The Ultimate Project has a website, and for always-current perspective on interstellar travel, it pays to read Paul Gilster’s Centauri Dreams weblog.

Interest in interstellar travel would ramp up if a truly Earth-like world were discovered around one of the Sun’s nearest stellar neighbors. Alpha Centauri, 4.36 light years distant, has the unique allure. Last year, I wrote a series of posts [1, 2, 3, 4] that explored the possibility that a habitable world might be orbiting Alpha Centauri B. In short, the current best-guess theory for planet formation predicts that there should be terrestrial planets orbiting both stars in the Alpha Cen binary. In the absence of non-gaussian stellar radial velocity noise sources, these planets would be straightforward to detect with a dedicated telescope capable of 3 m/s velocity precision.

Over the past year, we’ve done a detailed study that fleshes out the ideas in those original oklo posts. The work was led by UCSC graduate student Javiera Guedes and includes Eugenio, Erica Davis, myself, Elisa Quintana and Debra Fischer as co-authors. We’ve just had a paper accepted by the Astrophysical Journal that describes the research. Javiera will be posting the article to astro-ph in the next day or so, but in the meantime, here is a .pdf version.

Here’s a diagram that shows the sorts of planetary systems one should expect around Alpha Cen B. The higher metallicity of the star in comparison to the Sun leads to terrestrial planets that are somewhat more massive.

We’re envisioning an all-out Doppler RV campaign on the Alpha Cen System. If the stars present gaussian noise, then with 3 m/s, one can expect a very strong detection after collecting data for five years:

Here’s a link to an animation on Javiera’s project website which shows how a habitable planet can literally jump out of the periodogram.

I think the planets are there. The main question in my opinion is whether the stellar noise spectrum is sufficiently Gaussian. It’s worth a try to have a look…

436 again5
greg posted in exoplanet detection on January 22nd, 2008

Image Source.

There’s a provocative paper up on the astro-ph today. Ignasi Ribas and two collaborators are reporting the “possible discovery” of a 4.8 Earth mass planet in an exterior 2:1 mean motion resonance with the transiting hot Neptune Gliese 436b. Planet four three six b is the well-known subject of great consternation, great scientific value, and many an oklo.org post. (For the chronological storyline, see: 1 (for background), 2, 3, 4, 5, 6, 7, 8, 9, and 10.)

Here’s the basic idea. Ribas et al. note that a single-planet fit to the Maness et al. (2007) radial velocity data set (which is listed as gj_436_M07K on the systemic console) has a peak in the residuals periodogram at P~5.1866 days:

Using this periodogram peak as a starting point, they get a keplerian 2-planet fit that lowers the reduced chi-square from ~4.7 to ~3.7. They then point out that this detection can potentially be confirmed by measuring variations in transit timing. In their picture, the presently-grazing transit has come into visibility only within the last 2.5 years or so, as a result of orbital precession. The transit light curve should thus be showing significant variations in duration as well as deviations from a strictly periodic sequence of central transit times.

This will be a huge big deal if the claim holds up. For starters, it’ll provide a natural explanation for Gl 436b’s outsize eccentricity. And everyone’s been on the lookout for a strongly resonant transiting system with a short orbital period. For the time being, though, I’m withholding judgment. As a first point of concern, Ribas et al. are presenting a keplerian fit to the radial velocities. Yet for the orbital configuration they are proposing, it’s absolutely vital to take planet-planet interactions into account. One can see this by entering their fit into the console. (Use a mean anomaly at the first RV epoch 2451552.077 for planet b=40.441 deg, corresponding to their reported time of periastron of Tp_b=HJD 2451551.78, and a mean anomaly for planet c=268.14 deg, corresponding to their reported value of Tp_c=HJD 2451553.4.) One can also dial in a long-term trend if one wants, but this isn’t necessary. Once the fit is entered, the reduced chi-square is 3.7. Activate integration. (Hermite 4th-order is the faster method.) When the planets are integrated, their mutual interactions utterly devastate the fit, driving the reduced chi-square up to 85.018. Using the zoomer and the scroller, you’ll see that the integrated radial velocity curve and the keplerian curve start off as a good match, but then rapidly get completely out of phase.

In order to examine the plausibility of a two-planet fit in 2:1 mean motion resonance, one needs to fit the radial velocity data with integration turned on. It is also important to include the existing transit timing data in the fit (and to do this, it’s best to use the most recent, so-called unstable version of the console). Over at Bruce Gary’s amateur exoplanet archive (AXA), there are now three transit timing measurements listed, with the latest obtained by Bruce himself this past New Years Eve. The HJD measurements of central transit should be added to the gj436.tds file, along with the HJD 2454280.78149 +/- 0.00016 central transit time measured by Spitzer.

Ideally, the Spitzer secondary transit timing data should also be included, but at the moment, the distribution version of the console does not have the capability to incorporate secondary transit measurements. One approach would be to get a self-consistent fit, and then see whether the epoch of secondary transit matches that observed by Spitzer.

Have fun…

transit valuations19
greg posted in exoplanet detection on December 24th, 2007

Image Source

Discoveries relating to transiting extrasolar planets often make the news. This is in keeping both with the wide public interest in extrasolar planets, as well as the effectiveness of the media-relations arms of the agencies, organizations, and universities that facilitate research on planets. I therefore think that funding support for research into extrasolar planets in general, and transiting planets in particular, is likely to be maintained, even in the face of budget cuts in other areas of astronomy and physics. There’s an article in Saturday’s New York Times which talks about impending layoffs at Fermilab, where the yearly budget has just been cut from $342 million to $320 million. It’s often not easy to evaluate how much a particular scientific result is “worth” in terms of a dollar price tag paid by the public, and Sean Carroll over at Cosmic Variance has a good post on this topic.

For the past two years, the comments sections for my oklo.org posts have presented a rather staid, low-traffic forum of discussion. That suddenly changed with Thursday’s post. The discussion suddenly heated up, with some of the readers suggesting that the CoRoT press releases are hyped up in relation to the importance of their underlying scientific announcements.

How much, actually, do transit discoveries cost? Overall, of order a billion dollars has been committed to transit detection, with most of this money going to CoRoT and Kepler. If we ignore the two spacecraft and look at the planets found to date, then this sum drops to something like 25 million dollars. (Feel free to weigh in with your own estimate and your pricing logic if you think this is off base.)

The relative value of a transit depends on a number of factors. After some revisions and typos (see comment section for this post) I’m suggesting the following valuation formula for the cost, C, of a transit:

The terms here are slightly subjective, but I think that the overall multiplicative effect comes pretty close to the truth.

The normalization factor of 580 million out front allows the total value of transits discovered to date to sum to 25 million dollars. The exponential term gives weight to early discoveries. It’s a simple fact that were HD 209458 b discovered today, nobody would party like its 1999 — I’ve accounted for this with an e-folding time of 5 years in the valuation.

Bright transits are better. Each magnitude in V means a factor of 2.5x more photons. My initial inclination was to make transit value proportional to stellar flux (and I still think this is a reasonable metric). The effect on the dimmer stars, though was simply overwhelming. Of order 6 million dollars worth of HST time was spent to find the SWEEPS transits, and with transit value proportional to stellar flux, this assigned a value of two dollars to SWEEPS-11. That seems a little harsh. Also, noise goes as root N.

Longer period transits are much harder to detect, and hence more valuable. Pushing into the habitable zone also seems like the direction that people are interested in going, and so I’ve assigned value in proportion to the square root of the orbital period. (One could alternately drop the square root.)

Eccentricity is a good thing. Planets on eccentric orbits can’t be stuck in synchronous rotation, and so their atmospheric dynamics, and the opportunities they present for interesting follow-up studies make them worth more when they transit.

Less massive planets are certainly better. I’ve assigned value in inverse proportion to mass.

Finally, small stars are better. A small star means a larger transit depth for a planet of given size, which is undeniably valuable. I’ve assigned value in proportion to transit depth, and I’ve also added a term, Np^2, that accounts for the fact that a transiting planet in a multiple-planet system is much sought-after. Np is the number of known planets in the system. Here are the results:

HD 209458 b is the big winner, as well it should be. The discovery papers for this planet are scoring hundreds of citations per year. It essentially launched the whole field. The STIS lightcurve is an absolute classic. Also highly valued are Gliese 436b, and HD 189733b. No arguing with those calls.

Only two planets seem obviously mispriced. Surely, it can’t be true that HAT-P-1 b is 10 times more valuable than HAT-P-2b? I’d gladly pay $85,507 for HAT-P-2b, and I’d happily sell HAT-P-1b for $969,483 and invest the proceeds in the John Deere and Apple Computer corporations.

Jocularity aside, a possible conclusion is that you should detect your transits from the ground and do your follow up from space — at least until you get down to R<2 Earth radii. At that point, I think a different formula applies.

CoRoT-exo-2 c?14
greg posted in systemic faq, exoplanet detection on December 21st, 2007

Image Source.

The CoRoT mission announced their second transiting planet today, and it’s a weird one. The new planet has a mass of 3.53 Jupiter masses, a fleeting 1.7429964 day orbit, and a colossal radius. It’s fully 1.43 times larger than Jupiter.

The surface temperature on this planet is likely well above 1500K. Our baseline theoretical models predict that the radius of the planet should be ~1.13 Jupiter radii, which is much smaller than observed. Interestingly, however, if one assumes that a bit more than 1% of the stellar flux is deposited deep in the atmosphere, then the models suggest that the planet could easily be swollen to its observed size.

The surest way to heat up a planet is via forcing from tidal interactions with other, as-yet unknown planets in the system. If that’s what’s going on with CoRoT-exo-2 b, then it’s possible that the perturber can be detected via transit timing. The downloadable systemic console is capable of fitting to transit timing variations in conjunction with the radial velocity data. All that’s needed is a long string of accurate central transit times.

The parent star for CoRoT-exo-2-b is relatively small (0.94 solar radii) which means that the transit is very deep, of order 2.3%. That means good signal to noise. At V=12.6, the star should be optimally suited for differential photometry by observers with small telescopes. With a fresh transit occurring every 41 and a half hours, data will build up quickly. As soon as the coordinates are announced, observers should start bagging transits of this star and submitting their results to Bruce Gary’s Amateur Exoplanet Archive. (See here for a tutorial on using the console to do transit timing analyses.)

6 Gigabytes. Two Stars. One Planet.6
greg posted in exoplanet detection, worlds on December 14th, 2007

Image Source.

Another long gap between posts. I’m starting to dig out from under my stack, however, and there’ll soon be some very interesting items to report.

As mentioned briefly in the previous post, our Spitzer observations of HD 80606 did indeed occur as scheduled. Approximately 7,800 8-micron 256×256 px IRAC images of the field containing HD 80606 and its binary companion HD 80607 were obtained during the 30-hour interval surrounding the periastron passage. On Nov. 22nd, the data (totaling a staggering 6 GB) was down-linked to the waiting Earth-based radio telescopes of NASA’s Deep Space Network. By Dec 4th, the data had cleared the Spitzer Science Center’s internal pipeline.

We’re living in a remarkable age. When I was in high school, I specifically remember standing out the backyard in the winter, scrutinizing the relatively sparse fields of stars in Ursa Major with my new 20×80 binoculars, and wondering whether any of them had planets. Now, a quarter century on, it’s possible to write and electronically submit a planetary observation proposal on a laptop computer, and then, less than a year later, 6 GB of data from a planet orbiting one of the stars visible in my binoculars literally rains down from the sky.

It will likely take a month or so before we’re finished with the analysis and the interpretation of the data. The IRAC instrument produces a gradually increasing sensitivity with time (known to the cognescenti as “the ramp”). This leads to a raw photometric light curve that slopes upward during the first hours of observation. For example, here’s the raw photometry from our Gliese 436 observations that Spitzer made last Summer. The ramp dominates the time series (although the secondary eclipse can also be seen):

The ramp differs in height, shape, and duration from case to case, but it is a well understood instrumental effect, and so its presence can be modeled out. Drake Deming is a world expert on this procedure, and so the data is in very capable hands. Once the ramp is gone, we’ll have a 2800-point 30 hour time series for both HD 80606 and HD 80607. We’ll be able to immediately see whether a secondary transit occurred (1 in 6.66 chance), and with more work, we’ll be able to measure how fast the atmosphere heats up during the periastron passage. Jonathan Langton is running a set of hydrodynamical simulations with different optical and infrared opacities, and we’ll be able to use these to get a full interpretation of the light curve.

In another exciting development, Joe Lazio, Paul Shankland, David Blank and collaborators were able to successfully observe HD 80606 using the VLA during the Nov. 19-20 periastron encounter! It’s not hard to imagine that there might be very interesting aurora-like effects that occur during the planet’s harrowing periastron passage. If so, the planet might have broadcasted significant power on the decameter band. Rest assured that when that when their analysis is ready, we’ll have all the details here at oklo.org.

planet per week7
greg posted in exoplanet detection on November 30th, 2007

Image Source.

As the academic quarter draws to a close, it gets harder to keep up a regular posting schedule. This year, certainly, the difficulty has nothing to do with a lack of exciting developments associated with extrasolar planets.

A few unrelated items:

It appears that the HD 80606b Spitzer observations went smoothly, and that the data has been safely transmitted to Earth via NASA’s Deep Space Network. It is currently in the processing pipeline at the Spitzer Science Center. When it clears the pipeline, the analysis can start.

Back in September, I wrote a post about Bruce Gary’s Amateur Exoplanet Archive. This is a web-based repository for photometric transit observations by amateurs. With the number of known transits growing by the month, there’s a planet in transit nearly all of the time. Over 90 light curves have been submitted to the archive thus far. For transiting planets such as HD 189733b or HD 209458b, which have significant numbers of published radial velocity data, it’s very interesting to take the transit center measurements from Bruce’s archive and use them as additional orbital constraints within the console. The September post gives a tutorial on how to do this.

It really is turning out to be a banner year for extrasolar planets. As we head into December, this year is averaging more than one planet per week. The detection rate is more than double that of the previous four years.

The plot above gives a hint that Saturn-mass planets might wind up being fairly rare, as one might expect from the zeroth-order version of the core accretion theory. (For more information, this series: 1, 2, 3, 4, 5, 6, and 7 of oklo posts compares and contrasts the gravitational instability and core accretion theories for giant planet formation.)
Also, if you give talks, here’s a larger version of the above figure.

Another interesting diagram is obtained by plotting orbital period vs. year of discovery:

It’s possible that this diagram might be hinting that true Jupiter analogs are relatively rare. Could be that the disks around metal-rich stars are able to form Jovian mass planets and then migrate them in, while stars with subsolar metallicity form ice giants beyond the ice line. In this scenario, our solar system lies right on the boundary between the two outcomes.

It could also be the case that there are a whole slew of true-Jupiter analogs just on the verge of being announced. Time will tell.

And as always, it’s interesting to spend time with the correlation diagram tool over at exoplanet.eu.

160 basis points6
greg posted in exoplanet detection on November 25th, 2007

Image Source.

It’s sometimes a little weird to realize that my daily schedule is dictated by the orbits of alien planets. HD 80606b went through periastron passage at 07:00 UT last Tuesday, with the Spitzer Space Telescope’s rattlesnake’s eye vision trained intently upon it. Over the past few days, it’s been hurtling away from the star, gradually reducing its velocity as it climbs up the gravitational potential well of the star.

At 07:45 UT on Monday morning, HD 80606b is scheduled to go through inferior conjunction. In the 1.6% a-priori geometric chance that the orbital plane of the planet is in near-perfect alignment with the line of sight to the solar system, then it will be possible to observe the planet in transit. The 1.6% transit probability is fairly high for a planet with a period of 111 days, but much lower than the 15% probability that a secondary eclipse can be observed. If the planet is undergoing secondary eclipse, then we’ll know as soon as the Spitzer data comes in.

Back in early 2005, Transitsearch.org coordinated a campaign to check for transits of HD 80606b. At that time, there were fewer radial velocities available, and so the transit window was less well constrained. A number of observers got data, and there was no sign of transit, but the coverage was not good enough to rule out a transit. I’m thus encouraging observers to monitor HD 80606 during the next 48 hours on the off chance that it can be observed in transit. Given the small chance involved, it seems appropriate to refer to the transit probability in terms of basis points. As in, “In ‘05, we got about 40 basis points. That means there’s still 120 basis points out there to collect.”

HD 80606 is a visual binary. The companion, HD 80607, provides a good comparison star in telescopes with a large enough aperture under good seeing conditions. For most observers, however, the light from the two stars is combined. A transit by HD 80606b is expected to have a depth of order 1.4%, and (if its a central transit) will last about 14 hours. It’s a long-shot for sure, but worthwhile and fun nonetheless.

Got the ‘606 kickin’ & the 436 written2
greg posted in exoplanet detection on November 20th, 2007

Image Source.

As I write this, it’s JD 2454425.219 (17:16 UT, Nov. 20 2007). HD 80606 b whipped through periastron a little more than 10 hours ago, and the Spitzer Space telescope is literally just finishing its 31-hour observation of the event. Next comes the downlink of the data to Earth on the Deep Space Network, and then the analysis. Definitely exciting!

The Spitzer Space Telescope is scheduled to run out of cryogen in early 2009. When the telescope heats up, we’ll lose our best platform for mid-infrared observations of hot extrasolar planets, and so there was a palpable urgency last week as everyone prepared their proposals to meet the submission deadline for Spitzer’s last general observing cycle. During the next few years, there is going to be intense development of detailed 3D radiation-hydrodynamical models for simulating the time-dependent surface flows on extrasolar planets. These models will need contact points with hard data. It’s thus vital to bank as wide a variety of observations of as wide a variety of actual planets under as wide variety of different conditions as possible. A number of fascinating exoplanet observing proposals were submitted last week by a variety of highly competent teams. I’m urging that they all be accepted!

Most of the exoplanet observations that have been done with Spitzer have focused on tidally locked transiting planets on circular orbits. HD 189733b, HD 209458b, TrES-1 and HD 149026b are the flagship examples of this class. In the past year, however, eccentric transiting planets have started turning up. Gliese 436b (e=0.15) was the first, followed by HAT-P-2b (e=0.5), and HD 17156b (e=0.67).

Drake Deming, Jonathan Langton and I decided that the most interesting proposal that we could make would be for Gliese 436 b. This is the Neptune-mass, Neptune-sized planet transiting a nearby red dwarf star. Here’s the to-scale diagram of the 2.644-day orbit:

After Gliese 436b was discovered to transit last spring, it triggered a Joe Harrington’s standing Target of Opportunity program. Both a primary and a secondary transit were observed (see this post) which confirmed the startlingly high eccentricity, and which allowed an estimate of the planet’s temperature (or, more precisely, the 8-micron brightness temperature). This turned out to be 712±36 K, which is significantly higher than the ~650 K baseline prediction.

The hotter-than-expected temperature measurement could arise from a number of different effects (or combinations of effects). By measuring the secondary eclipse, you strobe one hemisphere of the planet. If there are significant temperature variations across the surface of the planet, then a high reading might arise from chancing on the hotter side of the planet. Alternately, the effective temperature implied by measuring the energy coming out at 8-microns could be seriously skewed if the spectrum of the planet has deep absorption or emission bands at the 8-micron wavelength. Another possibility is that we’re observing tidal heating in action. Gliese 436b is being worked pretty hard in its eccentric orbit, and it should be generating quite a bit of interior luminosity as a result. If its structure is similar to Neptune, then a 712K temperature is completely understandable.

Io, of course, is subject to a similar situation. Here’s a K-band infrared photo of Io in transit in front of Jupiter:

Image Source.

Gliese 436b is in pseudo-synchronous rotation, and spins on its axis every ~2.3 days. The eccentricity of the orbit leads to an 83% variation in the amount of light received from the star over a 1.3 day timescale. This leads to a complicated flow pattern on the surface.

Here’s what Jonathan Langton’s model predicts for the appearance of the hemisphere facing Earth at five successive secondary eclipses:

Globally, the hydrodynamical model produces a statistically steady-state flow pattern that is dominated by a persistent eastward equatorial jet with a zonally averaged speed of ~150 meters per second. This eastward flow in the planet’s frame produces a light curve in the lab frame that has a ~3 day periodicity. This period is significantly longer than both the planet’s orbital period and the planet’s spin period. Our Spitzer proposal is to observe a sequence of 8 secondary transits in hopes of confirming both the amplitude and the periodicity of this light curve.

It’s certainly the case that our current hydrodynamical model is not the definitive explanation of what these planets are doing. I won’t be at all surprised if the flux variation from eclipse to eclipse is more complicated than what we predict. I’m highly convinced, however, that the model is good enough to indicate that the situation on Gliese 436b will be interesting, dynamic, and complex. The actual variation in the real observations will provide an interesting and non-trivial constraint that a definitive model of the planet will need to satisfy. The observations, if approved, will thus be of great use to everyone in the business of constructing GCMs for short period planets.

Stay tuned…

55 Cancri - A tough nut to crack.5
greg posted in exoplanet detection, worlds on November 13th, 2007

Image Source.

As soon as the new data sets for 55 Cancri from the Keck and Lick Observatories were made public last week, they were added to the downloadable systemic console and to the systemic backend. The newly released radial velocities can be combined with existing published data from both ELODIE and HET.

Just as we’d hoped, the systemic backend users got right down to brass tacks. As anyone who has gone up against 55 Cnc knows, it is the Gangkhar Puensum of radial velocity data sets. There are four telescopes, hundreds of velocities, a nearly twenty year baseline, and a 2.8 day inner periodicity. Keplerian models, furthermore, can’t provide fully definitive fits to the data. Planet-planet gravitational perturbations need to be taken into account to fully resolve the system.

Eugenio has specified a number of different incarnations of the data set. It’s generally thought that fits to partial data sets will be useful for building up to a final definitive fit. Here’s a snapshot of the current situation on the backend:

The “55cancriup_4datasets” aggregate contains all of the published data for all four telescopes. This is therefore the dataset that is most in need of being fully understood. The best fit so far has been provided by Mike Hall, who submitted on Nov. 9th. After I wrote to congratulate him, he replied,

Thanks Greg, […] It actually slipped into place very easily. About 13-30 minutes of adding planets and polishing with simple Keplerian, then 25 iterations overnight with Hermite 4th Order.

The problem is that it seemed like I was getting sucked into a very deep chi^2 minimum, so getting alternative fits may be tricky!

Here’s a detail from his fit which illustrates the degree of difference between the Keplerian and the full dynamical model:

and here’s a thumbnail of the inner configuration of the system. It’s basically a self-consistent version of the best 5-Keplerian fit.

Mike’s fit has a reduced chi-square of 7.72. This would require a Gaussian stellar jitter of 6.53 m/s in order to drop the reduced chi-square to unity. Yet 55 Cancri is an old, inherently quiet star, and so I think it’s possible, even likely, that there is still a considerable improvement to be had. It’s just not clear how to make the breakthrough happen.

This situation is thus what we’ve been hoping for all along with the systemic collaboration: A world-famous star, a high-quality highly complex published data set, a tough unsolved computational problem, and the promise of a fascinating dynamical insight if the problem can be solved.

I’ll end with two comments posted by the frontline crew (Eric Diaz, Mike Hall, Petej, and Chris Thiessen) that I found quite striking. These are part of a very interesting discussion that’s going on right now inside the backend.

When something is this difficult to solve using the ordinary approaches, I start to look to improbable and difficult solutions. In the case of 55C, my hunch is that it’s a system where the integration is necessary, but not sufficient to build a correct solution. I think that the parameter space of solutions is so chaotic that the L-M minimization doesn’t explore it well, or that the inclination of the system is significant enough to skew the planet-to-planet interactions in the console, or both. Trojans or horseshoe orbits would fit these conditions. Perhaps other resonant or eccentric orbits would as well.

I think the high chi square results and flat periodograms after fitting the known planets also point to a 1:1 resonant solution or significant inclination. I just don’t think there’s enough K left to fit another significant planet unless it’s highly interactive with the others.

I’m going to keep working on this system in the hopes that we can find a solution (and because it’s really, really fun), but I suspect that a satisfactory answer won’t be found without a systematic search of the parameter space including inclination.

– Chris

“Nature is not stranger than we imagine but stranger than we can imagine.” Or words to that effect, I can’t remember who said that but in all probability this system shall have more questions answered about it (or not as is often the case!) by direct imaging e.g. such as by the Terrestrial Planet Finder (TPF) mission to show what is really happening (if it is ever launched). The 55 Cancri system is listed as 63 on the top TPF 100 target stars.

In the meantime, we struggle on… I don’t think I can add anything else to what Eric and everyone else has said…

– Petej

The latest on 55 Cancri8
greg posted in systemic faq, exoplanet detection on November 6th, 2007

Image Source.

Here’s a development that systemic regulars will find interesting! In a press release today, came announcement of the detection of a fifth planet in the 55 Cancri system (paper here). The new planet has an Msin(i) of 0.144 Jupiter masses, a 260-day orbital period and a low eccentricity. The detection is based on a really amazing set of additions to the Lick and Keck radial velocities:

For background on the 55 Cancri system, check out this oklo.org post from December 2005.

The outer four planets in the 55 Cancri system all have fairly low eccentricities in the new five-planet model. This leads to a diminished importance for planet-planet interactions, but nevertheless, the system does require a fully integrated fit. Deviations between the Keplerian and integrated models arise primarily from the orbital precessions of planets b, c, and e that occur during the long time frame spanned by the radial velocity observations.

Eugenio has added the velocities onto a fully updated version of the downloadable systemic console. The new version of the console adds a wide variety of new features (including dynamical transit timing) that were formerly available only on the unstable distribution. Check it out, and see the latest news on the console change log and the backend discussion forum. Over the next month, we’ll be talking in detail about the new features on the updated console.

Very shortly, a new entries corresponding to the updated 55 Cancri data sets will be added to the “Real Stars” catalog on the systemic backend. I’ll then upload my baseline integrated 5-planet fit to the joint Keck-Lick data set. I’m almost certain that with some computational work, this baseline model can be improved. Such a task is not for the squeamish, however. Obtaining self-consistent 6-body models to the 55 Cancri data set is a formidable computational task for the console. There are 29 parameters to vary (if the Lick, Keck, ELODIE and HET radial velocity data sets are all included). The inner planet orbits every 2.79 days, and the data spans nearly two decades. Fortunately, Hermite integration is now available on the console. Hermite integration speeds things up by roughly a factor of ten in comparison to Runge Kutta integration.

There have been hints of the 260-day planet for a number of years now because it presents a clear peak in the residuals periodogram. After the 2004 announcement of planet “e” in its short-period 2.8 day orbit, Jack Wisdom of MIT circulated a paper that argued against the existence of planet “e”, and simultaneously argued that there was evidence for a 260-day planet in the data available at that time. More recently, a number of very nice fully self consistent fits to the available data have been submitted to the backend (by, e.g., users thiessen, EricFDiaz, and flanker). Their fits all contain both the 2.8 day and the 260-day planets, and happily, are fully consistent with the new system configuration based on the updated velocities. Congratulations, guys!

Interestingly, the best available self-consistent fits to the system indicate that planets b and c do not have any of the 3:1 resonant arguments in libration. It will be interesting to see whether this continues to be the case as the new fits roll into the systemic backend.

“seventy six seven hundred”6
greg posted in exoplanet detection on October 20th, 2007

The flurry of activity surrounding the detection of the HD 17156b transits, combined with the start of the academic quarter here at UCSC, caused me to fall way behind on my stack. All of a sudden, over two weeks have passed with no oklo.org posts. I think that’s a record, unfortunately.

And it’s not as if there’s been no new exoplanet news. The past two weeks have seen the announcement of two more new transiting planets from Gaspar Bakos’ HATnet project. “Yo, what up TrES?” HAT-P-5b has a period of 2.79 days, and looks good in a 400-pixel wide everything’s to scale diagram:

There are enough transiting extrasolar planets now, so that it’s interesting to look for trends in the planetary properties. At Jean Schneider’s exoplanet.eu site, there’s a nifty set of php routines that make it very easy to dial up all the different correlation diagrams. With the inclusion of HD 17156b, HAT-P-5b, and HAT-P-6b, a plot of planetary radii vs. stellar metallicity is pointing to an interesting trend. It’s quite apparent that the metal-rich stars tend to harbor smaller planets. This seems to be indicating that metal-rich disks are yielding planets that have highly enriched heavy element fractions, which in turn is giving us an important clue into the planet formation process.

If we ignore Gliese 436b, which is far smaller in both mass and radius than all the other known transiting exoplanets, then there’s a fairly obvious hint of two separate sequences in the diagram — a large-radius sequence, and a small-radius sequence. Feel free to voice your opinion in the comments section…

It’d certainly be nice to get more transits by planets orbiting bright parent stars. To that end, it’s important to stress that literally every single planet transiting a V<13 parent star is located north of the the celestial equator. It’s pretty clear that the southern hemisphere Doppler-wobble planets have not been fully followed up with photometric campaigns. I’m thus keen to get the Southern-Hemisphere Transitsearch.org corps out on the sky during the coming austral summer. First on the list is HD 76700b. This 0.20 Jupiter-mass planet orbits with a period of 3.970985 days, and has an extensive and fairly recent set of published radial velocities. I just updated the orbital fit, and found that the transit windows are still quite narrow. A simple bootstrap analysis shows that the uncertainty in the time of the transit midpoint is not much wider than the expected duration of the transit itself. The star is just coming visible in the early morning, and so it should be straightforward to either confirm or rule out a transit for this particular planet.

Confirmed!6
greg posted in exoplanet detection on October 4th, 2007

Image Source.

I’m happy to report that HD 17156b is observable in transit, and that Transitsearch.org observers played the key role in the discovery.

Regular oklo.org readers are familiar with HD 17156 b. This planet has an orbital period of 21.2 days, which is nearly four times longer than any other known transiting planet, and an eccentricity of e=0.67 (even higher than HAT-P-2b’s eccentricity of e=0.5). The geometry of the eccentric orbit has a periastron angle of 121 degrees, which means that the planet is quite close to the star as it perforates the plane containing the line of sight to Earth. Here’s a scale model showing the planet (at equally spaced intervals), the star, and the orbit:

Early last week, I got e-mail from Mauro Barbieri, an Italian post-doctoral researcher who is working on the CoRoT satellite, and who’s based at LAM in Marseille, France. In his spare time, Mauro has been writing articles for popular astronomy magazines in Italy, and has worked to coordinate Italian amateur astronomers for participation in Transitsearch campaigns. On the night of Sept. 9/10, he recruited four Italian amateur astronomers to monitor the transit window.

Two observers in Northern Italy were clouded out prior to the start of the transit, but two others, D. Gasparri and C. Lopresti, were able to observe through much of the night in Central Italy. Their data looked promising, showing clear ingresses at the same time as the ingress was observed by Jose Almenara in the Canary Islands. Ron Bissinger, observing from Pleasanton, California, and seven hours farther west, was able to start observing just after the transit ended.

As it happened, Roi Alonso, another CoRoT postdoc, is good friends with Jose Almenara, who made the Canary Island observations on Sep 9/10. Barbieri and Alonso did a careful analysis of the three transit-bearing data sets, and concluded that the transit is present at 3-sigma, 5.3-sigma, and 7.9-sigma, respectively. Almenara’s data, in particular, is excellent, despite of the fact that high winds occluded part of the mid-transit time series. By last Friday, on the strength of the detections, we had begun drafting a paper that discusses the discovery.

The Almenara egress is particularly convincing (the bottom time series shows the result of subtracting out the best-fit transit signal):

Barbieri’s and Alonso’s fit to the data from Sep. 9/10 implies the following properties for the planet and the transit. The model consistently takes into account the eccentric character of the planetary orbit:

The fit to the data suggests that the radius of HD 17156b is just a bit larger than the radius of Jupiter. This is fully in line with our theoretical models of the planet. HD 17156b experiences strong tidal forces during its periastron passages. This tidal heating might be observable in the form of excess infrared radiation, but it is not serving to inflate the planet beyond its expected radius.

Needless to say, we were quite excited by the quality of the fit. Everything seemed to hang together quite well, but confirmation was essential. The full transit would be visible across the United States and Canada. I wrote to the transitsearch.org mailing list, urging observers to monitor the star through the night of Sep. 30/Oct. 1. Dave Charbonneau of Harvard had been following oklo.org, and was impressed by the Sep 9/10 data. He worked very hard to organize and coordinate observers, and he’ll be leading a follow-up paper that uses data from all the transits to improve the planetary and the orbital characterization. Dave was very generous in offering to notify us if a transit was confirmed by the cadre of observers that he’d recruited.

Sunday morning was perfectly clear in Santa Cruz. Ron Bissinger lives in Pleasanton, just forty miles to the Northeast, and was ready to observe. His pipeline is quite automated, and so if he could observe, I knew that we would rapidly rule out or confirm a transit.

Late Sunday afternoon, I went running, and noticed that a gloomy bank of clouds was visible over the Pacific to the west:

At dusk, it was still clear, but the weather forecast did not look good. The bands of clouds that had stayed offshore to the Northwest were beginning to slide across the skies. By midnight, it was evident that the Bay Area would not be producing useful photometry. Furthermore, all of Arizona and New Mexico were rained out. Observer after observer reported in to say that they had not gotten data. The only positive report of clear skies came from Dave, graduate student Philip Nutzman, and postdoc Jonathan Irwin, who had set up a small telescope on the roof of a Harvard/CfA building in Cambridge MA.

By Monday morning, however, it was clear that several observers had indeed managed to obtain data. In addition to the Harvard roof observations, Bill Welsh and Abhijith Rajan had obtained usable photometry from the Mt. Laguna Observatory run by San Diego State University. Dave reported to me that on Sunday, while Dave was visiting the Zoo with his daughter, Bill had called with the news that they had managed to secure a night on the telescope and were at that moment driving up the mountain. In addition, a report came from Tim Brown of Las Cumbres observatory that while their Hawaii site had been weathered out, observations had been successfully made from parking lot of the observatory headquarters in Santa Barbara. And in addition, Don Davies, a Transitsearch.org observer in Torrance California had obtained a 10,000 CCD frames under good sky conditions.

On Tuesday evening, I got a phone call from Dave. The transit was clearly visible in both the Cambridge and the Mt. Laguna data. Thirty minutes later, I got an e-mail from Don Davies, who, in the early stages of analysis was seeing a clear transit-like signal at the expected time. We signed on Davies as a co-author, added Dave’s personal communication to the paper draft, and submitted the discovery. The paper will be showing up on astro-ph today. Here’s a link to a .pdf version:

Barbieri et al. 2007, AA, submitted (157 kb)

As for HD 17156b itself, the transit should present a number of exciting opportunities for follow-up observations. The large orbit leads to a 26-fold orbital variation in the amount of flux received from the parent star. This should drive complex weather on the surface, and indeed, even the night side of the planet should be glowing from its own radiation. Here’s a frame from Jonathan Langton’s most recent simulation of the planet which shows the night side hemisphere:

And here’s a one-orbit 1 MB animation of the surface flow patterns, glowing and roiling with their own emitted heat.

Tonight’s the Night6
greg posted in exoplanet detection on September 30th, 2007

Image Source.

Tonight, Sept. 30/Oct 1, is the night to follow-up to confirm whether HD 17156b can be observed in transit. Earlier this morning, I sent the following e-mail to the Transitsearch Observers list:

Hello Everyone,

I’d like to alert you to an important follow-up opportunity TONIGHT to observe HD 17156 for a possible transit by its companion planet. North American Observers are best situated for the event.

HD 17156 b has been the topic of several blog posts on oklo.org, see: [1], [2], [3], [4], [5].

Photometry taken by Jose Manuel Alemenara Villa on the Sept. 9/10 opportunity was suggestive of a possible transit with duration 169 minutes, a photometric depth of 0.007 magnitudes, and a mid-transit time of HJD ~ 2454353.614. These values are all quite close to what one would expect if HD 17156b is really transiting.

If the event observed by Alemenara Villa is due to a transit, then the next transit will be centered at HJD~2454374.83 (CE 2007 October 01 07:55 UT Monday) with the transit beginning at about 06:30 UT.

Observing should start as soon as possible this evening, and observers are encouraged to take photometry for as long as possible.

My fit to the published radial velocities predicts a transit midpoint centered at HJD 2454374.87 (CE 2007 Oct. 01 08:52 UT Monday), with a +/- 0.3d uncertainty in the time of central transit. The Alemenara Villa event sits nicely inside this window.

Thanks very much!
best regards,
Greg

It looks like much of the Southwest is clouded out, and although the skies outside here in Santa Cruz are currently cobalt blue, it’s predicted that clouds and even rain will materialize after midnight. SoCal, however, and many locations in the midwest and east look good to go. Here’s a selection of California predictions from the clear sky clock. This is a cool graphical tool for use in scheduling observations. Dark blue is good, white is bad.

follow-up still in order…2
greg posted in exoplanet detection on September 18th, 2007

Image Source.

In the last post, I pretty much wrote off HD 17156 b, which was the subject of last week’s transitsearch.org photometric follow-up campaign. Ron Bissinger observed the star during the latter part of the transit window, and saw no evidence of a transit. Tonny Vanmunster wrote with the news that Belgium was clouded out.

Soon after the post went up, however, Jose Manuel Almenara Villa of the Instituto de Astrofisica de Canarias posted a comment:

Hi Greg,

I observed HD17156 in the transit window. Unfortunately the night was windy, affecting the small telescope so the photometry is not so clear as we would wish. Anybody else observe?

It’s possible that I have a central transit. I can show you some plots if you want. I will try to observe again on December 3 (I think that is my next opportunity).

Regards,
Jose

On Saturday, Jose sent me his photometric plots, I should point out that he emphasized once again that the night was windy. In his plots (I’ve rewritten the labels in illustrator so that they show up better on the narrow blog-page format) the black dots are individual observations (R filter, 7 s exposures), the red dots bin 6 observations, and the blue dots bin 12 observations.

On the night before the night of the transit window, he got baseline photometry which shows considerably less scatter, and which does a nice job of showing his excellent photometric technique:

He fit a simple trapezoidal transit template to his data. The resulting fit has a duration of 169 minutes, a depth of 0.007 magnitudes, and a mid-transit time (HJD) ~ 2454353.614. These values are all quite close to what one would expect if HD 17156 b really is transiting. The possible event ends just prior to the start of Ron Bissinger’s time series.

So what to think? It’s most important to reiterate Jose’s point that the weather was not particularly good, and that a block of critical data is missing during the event itself. I myself have contracted transit fever several times in the past, and have built up sufficient immunity to refrain from getting too excited. I think a conservatively realistic assessment would say that there’s still an 11% chance that HD 17156 b transits are occurring, and that the uncertainty in the window has been narrowed down significantly. Over the long run, if transitsearch.org is going to succeed, then its important to stay cautiously optimistic. The good thing about a transit is that it repeats with clockwork regularity (barring the unlikely, but tough-luck situation where dynamically induced precession of the node induces transit seasons.) The next chance to observe HD 17156 during the transit window falls to North America on Oct. 1, where hopefully there’ll be multiple observers on the sky. We’re bad - We’re Nationwide…

To end on a heartfelt note, I think that the global collaborative efforts that go into these transitsearch campaigns have been both fun and inspiring, even when the result is the high-probability flat-line light curve. It would be exciting, though, if Jose ends up leading a discovery paper with the other participating observers as co-authors.

Results4
greg posted in exoplanet detection on September 12th, 2007

Image Source.

It’s not looking good for transits by HD 17156 b. Ron Bissinger of Pleasanton, California obtained a block of photometric data that covered a significant chunk of the transit window. His time series lasts from JD 2454353.68 through 2452353.88, and shows no hint of an event:

His observations were taken just after the peak of the transit midpoint histogram:

No word yet on whether anyone in Europe or the eastern US were able to observe during the first half of the window. If you got data, let me know.

Also, the Gliese 176 window has opened up. If you’ve got a telescope, a CCD, and a free evening, you know what to do!

Discover a planet3
greg posted in exoplanet detection, non-technical posts on September 11th, 2007

Image Source.

My tight 30-minute layover in Denver turned into an eight-hour delay yesterday when a solenoid somewhere in our Boeing 777 malfunctioned just prior to pushback, giving me an unexpected opportunity to attempt to catch up on all the work that’s been piling up.

After 6 hours of tapping on the laptop, I’d exhausted my effectiveness, so I bought glossy magazines from the airport newstand. In the latest issue of Portfolio from Conde Nast, you can read an in-depth Vanity Fair style puff piece on ex-Tyco CFO Mark Swartz’s life in the Big House, and, in one of the advertisements, you’re encouraged to use a Visa “Signature” card to charge up some of the finer experiences in life. Quite to my surprise, #17 on a list that includes “See the Tony Awards live”, and “Test-drive a supercar”, is “Discover a planet”.

Now regular visitors to oklo.org all know that you can get your planet-discovery experience right here on the systemic backend without ever having to reach for your wallet. In fact, just yesterday, we learned from Gregory’s latest preprint on astro-ph that Eric Diaz (and a number of other systemic users) appear to have made the first characterizations of the most statistically probable planetary system fits to the HD 11964 radial velocity data set.

The HD 11964 data set was published by Butler et al. (2006). Two planets are already known to orbit this star. HD11964 b has roughly 1/3rd of a Saturn mass and a ~38-day orbit, whereas HD 11964 c is a sub-Jovian mass planet on a ~2110-day orbit. There’s a wide dynamically stable gap between the two planets, making this system a fertile hunting ground for additional companions.

Gregory does an extensive statistical analysis and argues that there’s strong evidence for a sub-Saturn mass planet on a year-long orbit. Eric Diaz’s version of this planet shows up in the fit that he submitted to systemic back in July 2007:

Eric also suggests the presence of a 12.4-day planet in the system. The Gregory analysis suggests that this planet is not statistically significant, but I’m going to add it to the transitsearch.org unpublished candidates list. There’s certainly no reason not to have a look-see if anyone has unused photometric capability.

HD 17156 at inferior conjunction (right now!)5
greg posted in exoplanet detection on September 10th, 2007

Image Source.

It’s 01:58 UT Sep. 10, and HD 17156 has moved into its transit window. Hopefully photometric transit observers across Europe have clear skies. If you’re collecting data, drop us a note on the comments page!

Sep 09, 2007 Europe Satellite Map Source.

Most of California looks pretty good for catching the latter part of the transit window once it gets dark tonight. I was up on Mt. Hamilton last night, and even though it was clear, there was a strong smell of smoke in the air. Bits of gray ash from the nearby forest fires were floating down like snow, and so they couldn’t open the dome of the 36-inch.

Image Source.

The odds of a HD 17156 transit are 10.9%, so its best not to get hopes up too high. Its always good to have the next candidate ready to go, and as luck has it, there’s another good one in the hopper.

Endl et al. have published a preprint describing the discovery of a Neptune-mass planet orbiting the nearby red dwarf star Gliese 176 (aka HD 285968). This discovery is further evidence in favor of the core-accretion prediction that Neptune planets should be common around low-mass stars whereas Jovian-mass planets should be relatively rare. Endl et al.’s new planet has an orbital period of 10.24 days, an a-priori transit probability of 3%, and an expected transit depth of 0.4%. This is a low-amplitude signal, but it is nevertheless accessible to many experienced amateur astronomers. The discovery paper makes no mention of a photometric transit search, making this planet a very attractive Transitsearch.org candidate. The star is located at RA 04:43, Dec +18:57, and the next transit window is centered on Sep. 15, 2007.

HD 17156 b5
greg posted in exoplanet detection, worlds on September 4th, 2007

Image Source.

Last week, I wrote a post introducing HD 17156 b, a Jovian planet on a highly eccentric 21.2-day orbit around a V=8.17 solar-type star lying 250 light-years away in Cassiopeia (RA=2h 50m, Dec=72 deg).

A photometric check for transits by HD 17156 b was reported in the discovery paper, but due to the nearly three-week orbital period, it was only possible to rule out about 25% of the transit window. Given the highly favorable geometry of the planetary orbit, this means that there’s an impressive ~11% chance (8.25% if you take the discount) that the planet can be observed in transit. The expected transit depth is a very respectable 1%, and given the bright parent star, it’s a straightforward detection for small-telescope observers everywhere in the Northern Hemisphere.

What’s it worth to catch HD 17156 b in transit? From a crass cash-money standpoint, one can estimate a dollar value. Because the planet has a long period and an eccentric orbit, it would be the first transiting example of its kind, and would thus be expected to generate a fairly large number of citations. From a career standpoint, an ADS citation is worth at least $100 (see, e.g. here). Based on the citation count for the TrES-1 discovery paper (144 citations in three years) it’s reasonable to expect that at one decade out, a HD 17156 b transit would garner of order 200 citations, for a conservative total value of 20K. Given the 10% probability of the transit coming through, the resulting expectation value is equivalent to having twenty Benjamins floating down from the black velvet of the night sky.

I used the systemic console’s bootstrap utility to generate a set of orbital fits to the published radial velocities for HD 17156. Each orbital fit describes a unique sequence of central transit times. For a particular transit opportunity, the aggregate of predicted central transit times from the different fits can be plotted as a histogram. Here’s the resulting plot for the transit opportunity that’ll occur next Monday (HJD 2454353.68):

The uncertainty in the time of central transit is ~0.3 days. A window this narrow is rare for a planet that hasn’t yet been thoroughly checked. In fact, as far as Transitsearch.org opportunities are concerned, it doesn’t get much better than this. Extending our opportunity cost analysis, the expected monetary return for observations within the 1-sigma transit window is an impressive $114 per hour. (Only rarely does the expected return per hour exceed minimum wage for existing transit opportunities.)

Scientifically, a transit by HD 17156 b would certainly be very exciting. The planet should be heating up very rapidly during its periastron passage, which should spur the generation of hemispheric-scale vortices and an 8-micron light curve that’s detectable with the Spitzer telescope. Observation of the secondary eclipse (assuming it occurs) would allow for a measurement of the global planetary temperature near the orbital apastron.

The frame above is from a hydrodynamical study of HD 17156 b that Jonathan Langton has just finished computing. If all the talk of dollars, ephemerides, opportunity cost, and expectation value is leaving you stressed out, then just kick back with this fat 1.0 MB .mov of the simulation and get your groove on.

Countdown1
greg posted in exoplanet detection, worlds on August 20th, 2007

Image Source.

August 1st marked the most recent ‘606 day, which came and went without wide remark. Perhaps this was because in late Summer, HD 80606 rises and sets in near-synch the Sun, and is thus lost from the Earth’s night skies.

At the moment, HD 80606b is headed back out toward apastron.

The global storms and shockwaves that were unleashed at the beginning of August are dissipating rapidly, and the flux of heat from the planet is likely fading back down to the sullen baseline glow that arises from tidal heating.

HD 80606’s next periastron passage occurs on November 20th, and the Spitzer Space Telescope is scheduled to observe the whole event (details here). It’s going to be a big deal. Spitzer can only observe HD 80606 during two three-week windows each year, and fortunately, the Nov. 20th Periastron passage occurs during one of these windows. It’s literally the only opportunity to catch HD 80606 b’s big swing before Spitzer’s cryogen runs out in 2009.

The orbital geometry of the periastron passage looks like this:

Each marker of the orbit is separated by one hour. The prediction for the pseudo-synchronous rotation of the planet is also indicated. The planet should be spinning with a period of 36.8 hours. Jonathan Langton’s hydrodynamics code predicts what the temperature distribution on the planet should look like at each moment from Spitzer’s viewpoint in our solar system:

Transitsearch.org observers have covered a number of the HD 80606 b transit opportunities, and it seems pretty certain that the planet doesn’t transit. This isn’t surprising. The geometry of the orbit is such that when the planet crosses the plane containing the line of sight to the Earth, it’s quite a distance away from the star. Not so, however, for the secondary transit. There’s a very respectable 15% chance that Spitzer will detect a secondary transit centered two hours prior to the periastron passage.

Even if the planet doesn’t transit, we should be able to get a good sense of the orbital inclination from the shape of the light curve. If the orbit is nearly in the plane of the sky, then we should see a steady rise followed by a plateau in the 8-micron flux coming from the planet. For more nearly edge-on configurations, the flux peak should be clearly discernable. The observations are scheduled to start 20 hours prior to periastron and end 10 hours after.