M for all and all for M

superEarths and subNeptunes

I’m always impressed by the efficiency with which red dwarfs pack hydrogen, the stuff of flammable zeppelins, into such a small space: Gliese 1214 is more than twice as dense as led. The density of the Sun, on the other hand, is bubblegum by comparison.

Gliese 1214b’s orbital period is a mere 1.58 days. Its 0.014 AU separation from the system barycenter is the smallest yet measured for any planet. Yet because of the high red dwarf density, the star-planet configuration is actually rather spacious. Here’s the system to scale:

Gl1214toScale

It’s interesting to compare this diagram with that of a genuinely close-in planet such as HAT-P-7b, which actually has a somewhat longer 2.2 day orbital period:

At a given period, a red dwarf fills much less of a planetary orbit than does a Sun-like star. If the occurrence rate of planets at a specified period is the same for stars of different masses, then one needs to look at $\sim(M_{\odot}/M_{\rm RD})^{2/3}$ times more red dwarfs than Sun-like stars to find a given number of transits with a particular period.

Gliese 1214b lies at enough stellar radii from Gliese 1214 that its a-priori transit probability was only about 7%. The Mearth survey currently covers only ~2000 stars, and so the fact that the discovery was made so quickly was probably not luck, but rather points to the existence of a very large number of low-mass planets orbiting small stars.

Let’s face it. The big dough goes to chase potentially habitable transiting planets. With this metric, the red dwarfs come out way ahead. If red dwarfs and Sun-like stars have equal occurrence fractions for planets with Earth’s mass and insolation, then a low-mass red dwarf has roughly four times the probability of a Sun-like star of harboring a transiting potentially habitable planet. Twice the temperature means one-sixteenth the area and the square root of sixteen is four. The red dwarfs also present a number of other advantages, see e.g. here, here, and here.

Ryan Montgomery and I have a recent paper out which foreshadows what I think is the inevitability of transit surveys that use the Mearth strategy to target true-Earth analogs the habitable zones of the lowest-mass red dwarf stars. Mearth  is itself very well-positioned to expand in this direction. I also think that a lot of effort will continue to shift toward improved Doppler-velocity capability in the near-infrared (see, e.g. this recent paper by Jacob Bean and collaborators which describes the use of ammonia gas in a glass cell to imprint a forest of fixed reference lines on a K-band stellar spectrum).

A last note: Twelve-Fourteen-b is likely to become a favorite target for small-telescope observers, so I made sure to add it to the Transitsearch.org candidates table. Now that classes are done for the quarter, I’ve been going through the literature and adding or updating one or two planets a day. It’s tedious work, but I’ve noticed some interesting upcoming opportunities, which I’ll be writing about soon. For transit-themed ephemera and the latest celebrity gossip, look no further than the transitsearch twitter stream: http://twitter.com/Transitsearch.

And a postscript: In the comments, reader cwmagee points out that the implication of the post is that the HAT-P-7 and Gl1214 diagrams are to scale which eachother, but that’s not the case. He attached a version which shows a to-scale comparison of both systems:

Red dwarfs are small!

Arcada Fog

From Garth von Ahnen's Arcada Fog
On the Shoulders of Giants: from Garth von Ahnen's Arcada Fog

Adriane Steinacker teaches one of the few undergraduate-level History of Astronomy courses in the country as part of our UCSC atronomy curriculum. She sent me this Youtube link to the work of one of her students — Garth von Ahnen — who is majoring in art and minoring in astronomy.

Garth has made a remarkable video confection that weaves together an interpretation (with artistic license) of the historical trajectory by which the planetary orbits came to be understood. You’ve simply gotta watch it! In Garth’s words:

All Characters, Events, Places and Various Concepts of the structure of the Solar System are entirely non-fictitious. Any similarity to real or once real Characters, Events, Places, Concepts or Mooses are not coincidence, but both purposeful and slanderous, based on historically accurate hearsay, innuendo and exaggeration, except for Newton using a hoolah hoop, which never actually happened according to anyone.

The piece is far richer than the worn-out versions presented in the Astronomy 101 textbooks. I’ll admit I had to consult the wikipedia for the back story on Jost Bürgi. The original source of Tycho Brahe’s Pet Moose, who comes from left-field to play a starring role in von Ahnen’s version of events, is Pierre Gassendi’s 1654 biography Tycho Brahe, the man and his work (original in Latin).

Mearth!

M8arth

Of course, there are still 7 hours and 13 days left until the close of 2009, but I’ve got every confidence that the discovery of the decade has landed on the ground. The Mearth project has found a transiting 6.55 Earth-mass planet in orbit around the nearby red dwarf star GJ 1214. The parent star is bright enough, and the planet-star area ratio is large enough so that direct atmospheric characterization will be possible not just with JWST, but with HST. Incredible. I’m inspired, invigorated, envious. This discovery is a game changer.

The GJ1214 discovery is all over the news today. The coverage is deservedly laudatory, but interestingly, the most dramatic aspect of the detection received rather short schrift. This is easily the most valuable planet yet found by any technique, and the discovery, start to finish, required an investment of ~500K (along with the equivalent of 1-2 nights of HARPS time to do the follow-up confirmation and to measure the planet’s mass). By contrast, well over a billion dollars has been spent on the search for planets.

I’m milking that contrast for drama, of course. It’s true that GJ1214b is low-hanging fruit. The team with the foresight to arrive on the scene first gets to pick it. And the last thing I’m suggesting is a cut in the resources devoted to exoplanet research — it’s my whole world, so to speak. I do think, though, that Mearth epitomizes the approach that will ultimately yield the planets that will give us the answers we want. You search for transits among the brightest stars at given spectral type, and you design your strategy from the outset to avoid the impedance mismatches that produce bottlenecks at the RV-confirmation stage.

There’s a factor-of-fourteen mass gap in our solar system between the terrestrial planets and the ice giants, and so with the discovery of Gl 1214b (and the bizzare CoRoT-7b) we’re getting the “last first look” at a fundamentally new type of planet. CoRoT-7b is clearly a dense iron-silicate dominated object, but it likely didn’t form that way. Gliese 1214b’s radius indicates that it probably contains a lot of water. I think this is going to turn out to be the rule as more transiting objects in the Earth-to-Neptune mass range are detected.

So what next? With a modest increase in capability, Mearth is capable of going after truly habitable planets orbiting the very nearest stars. I think it’s time to put some money down…

Earthbound Planet Search

61 Vir b (simulation by J. Langton, Principia College)
61 Vir b (simulation by J. Langton, Principia College)

The ranks of the super Earths and the sub Neptunes continue to grow! In two papers that have been accepted by the Astrophysical Journal, and which will be coming out in tomorrow’s astro-ph mailing, the Earthbound Planet Search team is announcing the detection of very low mass planets orbiting the nearby solar twins 61 Vir and HD 1461. (Link to paper #1, link to paper #2).

The 61 Vir system is particularly compelling. The radial velocity data for this star indicate that at least three planets are present, with an architecture that’s quite a bit more crowded than the Sun’s terrestrial planet zone:

61Virorbits

The innermost planet, 61 Vir b, with Msin(i)~5.5 Earth masses, has a radial velocity half-amplitude K=2.15 m/s, which puts it in league with Gl581e (with K=1.9 m/s) and HD40307b (with K=2.0 m.s) as the lowest-amplitude Doppler detections to date.

We’ve adopted the Systemic Console software to analyze the Doppler velocities that are produced by the Earthbound Planet Search. I’ve written a tutorial (link here) that explores the 61 Vir dataset in detail, and shows how the planets are extracted.

parallel observing

noisydata

As the decade draws to a close, it’s hard not to be amazed at the progress that’s been made on every research front related to extrasolar planets.

An area that I think is now ripe for progress comprises coordinated multi-observer checks for transits by super-Earth/sub-Neptune planets. There are now over thirty known extrasolar planets with Msin(i)’s less than that of Gliese 436b (which tips the scales at 23 Earth masses). Of these, only CoRoT-7b has so far been observed to transit, and it’s very probable that the current catalog of low-mass RV-detected planets contains one or more transiting members. Needless to say, it’d be very interesting to locate them.

To my knowledge, the lowest-amplitude transits that have been observed by amateur astronomers have been those by HD 149026b. This anomalously dense Saturn-mass planet induces a photometric transit depth of roughly 0.4%.  State-of-the-art amateur detections show the transit very clearly. Here’s an example (the observer was Luboš Brát of the Czech Republic) taken from the TRESCA database:

149026sampletransit

The identification of transits by small planets certainly won’t be a picnic. Super-Earths and  sub-Neptunes orbiting G and K stars present targets that are intrinsically much tougher than HD 149026. Unless the parent star is a red dwarf, the expected transit depths will generally be less than 0.1%, and it’ll be extremely difficult for a single small-telescope observer to obtain a definitive result.

On the other hand, if a platoon of experienced observers mount a coordinated campaign on a single star, then there’s a possibility that a startlingly good composite light curve might be obtained. In theory, if one were to combine the results from sixteen independent observers, one could obtain a light curve of the equal signal-to-noise as the HD 149026b curve shown above, but for a planet with a transit depth of only 0.1%.

I spent time this weekend making sure that the transitsearch.org transit predictions for the known RV-detected low-mass planets are as up-to-date and accurate as possible. I found that HD 7924 is a good candidate star with which to test a coordinated observing strategy. The star harbors a low-mass RV-detected planet was announced earlier this year (discovery paper here):

hd7492

HD 7924b has Msin(i)~10 Earth Masses, a P=5.3978d orbital period, and a 6.7% a-priori chance of being observable in transit. The (folded) photometry in the discovery paper is of quite high quality, and shows that the star is not photometrically variable. The photometry also indicates that transits with depth greater than 0.05% are probably not occurring. The parent star, HD 7924 is a K-dwarf, with a radius of something like 78% that of the Sun, which means that if the planet is a sub-Neptune it’ll have a central transit depth of order 0.075%, whereas if it is a rocky object, the depth will likely be less than 0.05%. The 1-sigma uncertainty on the time of the transit midpoint is about 0.35 days. The parent star has V=7.2, and with Dec=+76 deg, it’s circumpolar for high-latitude observers (RA=01h 21m).

Here are the next predicted transit midpoints (dates and times are UT):

HJD        Y    M  D  H  M
2455182.04 2009 12 16 12 51
2455187.01 2009 12 21 12 14
2455192.41 2009 12 26 21 48
2455197.81 2010  1  1  7 21
2455203.20 2010  1  6 16 54
2455208.60 2010  1 12  2 28

Because HD 7924b’s period is known to an accuracy of 0.0013 days (2 minutes), participating Northern-hemisphere observers can obtain data during any of the upcoming opportunities. Their light curves, once standardized, can in theory be stacked to obtain increased precision. It would be very interesting to get a sense of the practical limits inherent in such an approach. I think the best way to test the limits is to give the observations a try!

that golden age

planetsareeverywhere

I’m nostalgic for ’97, when the discovery of a new extrasolar planet was literally front-page news. What’s now cliche was then fully viable poetic sweep. Epicurus and his multitude of worlds. Bruno burning at the stake. In that frame of mind, it’s fascinating to go back and read John Noble Wilford’s extended New York Times piece, written at the moment when the number of known extrasolar planets equaled the number of planets in our own solar system.

Some of the hyperbole still seems fresh, especially with regard to the frequency and diversity of planetary systems:

And the discoveries may be only beginning. One recent study suggested that planets might be lurking around half the Milky Way’s stars. Astronomers have already seen enough to suspect that their definition of planets may have to be broadened considerably to encompass the new reality. As soon as they can detect several planets around a single star, they are almost resigned to finding that the Sun’s family, previously their only example, is anything but typical among planetary systems.

At the recent Porto conference, the Geneva team not only reiterated their claims regarding the frequency of low-mass planets, but actually upped their yield predictions. According to a contact who heard Stephane Udry’s talk, the latest indication from HARPS is that between 38% (at the low end) and 58% (at the high end) of nearby solar-type stars harbor at least one readily detectable M<50 Earth-mass planet. This is quite extraordinary, especially given the fact that were the HARPS GTO survey located 10 parsecs away and observing the Sun, our own solar system (largely in the guise of Jupiter’s decade-long 12-m/s wobble)  would not yet be eliciting any particular cause for remark.

It also looks like planets beyond the snowline are quite common. In yesterday’s astro-ph listing, there’s a nice microlensing detection of a cold Neptune-like planet orbiting a ~0.65 solar mass star with a semi-major axis of at least 3 AU. The microlensing detections to date indicate that Neptune-mass objects are at least three times as common as Jupiter mass objects when orbital periods are greater than five years or so. Microlensing detections are an extremely cost-effective way to build up the statistics of the galactic planetary census during belt-tightening times. Much of the work is done for free by small telescope observers.

microlens20091208

Yet another dispatch pointing toward a profusion of planets comes from an article posted last week on astro-ph by Brendan Bowler of the IfA in Hawaii. Work that he’s done with John Johnson and collaborators indicates that the frequency of true gas giant planets orbiting intermediate-mass stars (former A-type stars like Sirius that are now in the process of crossing the Hertzsprung gap) is a hefty 26% within ~3 AU.

An embarrassment of riches? Certainly, the outsize planetary frequency means that the cutting-edge of the planet-detection effort will be shifting toward the Sun’s nearest stellar neighbors, as these are the stars that offer by far the best opportunities for follow-up with space-based assets such as HST, Spitzer, JWST et al.

As competition for ground-based large-telescope RV confirmation of run-of-the-mill planet transit candidates orbiting dim stars heats up, the threshold magnitude (at a given bandpass) at which stars become largely too faint to bother with will grow increasingly bright. We’re talking twelve. Maybe nine. Pont et al., in their discovery paper for OGLE-TR-182b refer to this threshold as the “Twilight Zone” of transit surveys:

The confirmation follow-up process for OGLE-TR-182 necessitated more than ten hours of FLAMES/VLT time for the radial velocity orbit, plus a comparable amount of FORS/VLT time for the transit lightcurve. In addition, several unsuccessful attempts were made to recover the transit timing in 2007 with the OGLE telescope, and 7 hours of UVES/VLT were devoted to measuring the spectroscopic parameters of the primary. This represents a very large amount of observational resources, and can be considered near the upper limit of what can reasonably be invested to identify a transiting planet.

Transitsearch back on the air

transitsearchsplashimage

A quick addendum to the previous post. After a rather lengthy and undeserved “vacation”, Transitsearch.org is back on the air. The old website is running as a placeholder, and updated content will follow on soon.

I’ve moved the front-end of the transitsearch site to the hosting service that runs oklo.org, so the real URL is www.oklo.org/transitsearch/ By Dec. 10th, the domain name transfer will be complete, and the old www.transitsearch.org address should properly redirect.

Further updates can be had by subscribing to Transitsearch.org’s twitter stream: http://twitter.com/Transitsearch. We’re planning events to surround the next ‘606 day, and we’re also planning to organize a campaign for the HAT-P-13c transit opportunity that’s centered on April 12, 2010.

Arrived: ETD

Transits come in all shapes and sizes

A recent e-mail from Bruce Gary prompted me to pay a return visit the Exoplanet Transit Database (ETD) which is maintained by the variable star and exoplanet section of the Czech Astronomical Society. I came away both impressed and inspired. The ETD is really leveraging the large, fully global community of skilled small-telescope photometric observers.

There are hundreds of citizen scientists worldwide who have demonstrated the ability to obtain high-quality light curves of transiting extrasolar planets. I’ve developed many contacts with this cohort over the past decade through the Transitsearch.org project, and small-telescope observers played a large role in the discovery of the two longest-period transits, HD 17156b, and HD 80606b.

Once a particular planet has been found to transit, there is considerable scientific value in continued monitoring of the transits. Additional perturbing planets can cause the transit times to deviate slightly from strict periodicity, and a bona-fide case of such transit timing variations (TTVs)  has become something of a holy grail in the exoplanet community. A perturbing body will also produce changes in the depth and duration of transits as a consequence of changes in the orbital inclination relative to the line of sight. Moreover, for favorable cases, a large moon orbiting a transiting planet can produce TTVs detectable with a small telescope from the ground.

New transiting planets are being announced at a rate of roughly one per month. The flow of fresh transits continuously improves the odds that systems with detectable TTVs are in the catalog, but also makes it harder for any single observing group (e.g. the TLC project) to stay on top of all the opportunities.

The Exoplanet Transit Database maintains a catalog of all publicly available transit light curves. At present, there are 1113 data sets distributed over 58 transiting planets. The ETD site provides a facility for photometric observers to upload their data, and also provides online tools for observation scheduling and automated model fitting. Simply put, this is a groundbreaking resource for the community.

The ETD also provides concise summaries of the state of the data sets. Light curves are divided into five quality bins, depending on the noise level, the cadence, and the coverage of the photometry:

Picture 4

It’s interesting to go through the summary reports for each of the transiting planets. Here’s the current plot of predicted and observed transit times for Gliese 436b, the famously transiting hot Neptune:

ETDgl436b

The data show no hint of transit timing variations. (So what’s up with that e?)

In other cases, however, there are hints that either the best-fit orbital period needs adjustment, or that, more provocatively, the TTVs are already being observed. TrES-2 provides an intriguing example:

ETDTres2

In sifting through the database, it looks like XO-1, CoRoT-1, Hat-P-2, OGLE-TR-10, OGLE-TR-132, OGLE-TR-182, TrES-1, TrES-3, and WASP-1 are all worthy of further scrutiny.

Over the past year, as a result of Stefano Meschiari’s efforts, the Systemic Console (latest version downloadable here) has been evolving quite quickly behind the scenes. Stefano and I are working on a paper which illustrates how the console can be used to solve the TTV inverse problem through the joint analysis of radial velocity and transit timing data. In the meantime, it’s worth pointing out that the ETD database lists transit midpoints in HJD for all of the cataloged light curves. These midpoints can easily be added to the .tds files that come packaged with the console.

lithium-induced speculations

Lithium Depletion

Israelian et al’s Nature paper on the planet-stellar lithium correlation (featured in last week’s post) caused quite a stir in the community. The depletion of lithium in the atmosphere of a solar-type star seems to be a prerequisite for the presence of a detectable planetary system. Here’s the paper’s plot again, this time, with Alpha Cen A added for effect. lithiumwalphacen

Had Israelian et al.’s paper come out a decade ago, much of the ensuing hubub would have focused on the fact that low lithium abundance is an effective signpost to planetary systems. Nowadays, though, mere detection of new planets is passé. Everyone knows there are tons of planets out there. Focus is shifting to finding the lowest-mass (and preferably transiting) planets around the brightest M, K, and G main sequence stars in the Sun’s neighborhood. There is a short, highly select, list of worlds that have been, and will eventually be, followed up to great advantage with HST, Warm Spitzer, and JWST.  All of the Sun’s most alluring stellar neighbors are under heavy and ongoing scrutiny, and in fact,  it’s these particular stars (in the form of the HARPS GTO list) that enabled discovery of the planet-lithium correlation.

So planet-finding utility aside, the intense interest in the planet-lithium effect stems from the fact that it’s guaranteed to be imparting an important clue to the planet-formation process.

With over 400 planets known, clear populations are starting to emerge. It’s remarkable that the strength of the lithium-planet correlation seems to be largely independent of the masses and periods of the planets themselves. The mass-period diagram for planets, on the other hand, shows that there are at least three distinct concentrations of planet formation outcomes:

currentpop2009

It’s important to keep in mind that Israelian et al.’s correlation holds over only a very narrow range of stellar temperature. The M-dwarfs (Gliese 581, Gliese 876), the K-dwarfs (HD 69830, Alpha Cen B), and the F-dwarfs (Upsilon Andromedae) all fall outside the band of utility. This dovetails nicely with standard models of stellar evolution that suggest the amount of Lithium depletion in stars with masses very close the the Sun (that is, stars falling in the narrow effective temperature range of the above plot) depends sensitively on both the efficiency of convection and also on rotational mixing. That is, the stars that show the lithium-planet effect, are exactly the stars where subtle differences in properties seem to generate a big effect on lithium abundance.

After writing last week’s post, I got an e-mail from Jonathan Irwin (of MEarth fame) who makes several interesting points:

The low lithium could be more of a coincidence resulting from the long-lived circumstellar disks that are presumably needed to form planets.

Mediation of the stellar rotation rates by long-lived disks is thought to be responsible for generating the wide dispersion in rotation rates observed in open clusters around 100Myr age, and there have been suggestions (e.g. Denissenkov et al.’s paper that appeared on astro-ph 2 weeks ago) that the slowly-rotating stars evolve developing some degree of decoupling of the rotation rates of their radiative core and convective envelope, whereas the rapidly-rotating stars evolve more like solid bodies.

Bouvier (2008) has suggested that the shear at the radiative convective boundary resulting from this could mix lithium into the interior more efficiently, and thus could result in lower lithium for stars that were slow rotators, preserving evidence of their rotational history even though the final rotation rates all converge by the solar age.  Some evidence for this last part exists in the form of a correlation between rotation and lithium in young open clusters such as the Pleiades.

A hypothesis along these lines seems quite appealing to me. As long as a protoplanetary disk is present, and as long as its inner regions are sufficiently ionized, then there’ll be a connection between the stellar magnetic field and the magnetic field of the disk. To a (zeroth) degree of approximation, the equations of ideal MHD allow us to envision the situation as consisting of a rapidly rotating star connected to a slower-rotating disk by lot of weak rubber bands. The net effect will be to slow down the stellar rotation to bring it into synch with the rotation at the inner edge of the disk.

Trying to sound like a tough-guy, I stressed the importance of predictions in last weeks post. If Irwin’s hypothesis is correct, then the formation of the Mayor et al. 2008 planet population is associated with disks that contain lots of gas, even in regions interior to R~0.1 AU. I’d thus expect that the “super Earths” are actually “sub Neptunes”, and that we can expect considerable H-He envelopes for the majority of these planets.

Another speculative prediction concerns the stars that aren’t depleted in lithium. In Irwin’s picture, these stars had short-lived disks and lost their gas relatively rapidly. This shouldn’t hinder the formation of terrestrial planets, but one would expect that the final configurations of the rocky planets would sport higher eccentricities, as there was little or no gas to damp the orbits down during the final stages of terrestrial planet accretion (see this paper for more on this).

Back on line

the bugs have been swept out...

We’re back on line after a skin-crawling attack that exploited the WordPress installation to rebrand the oklo.org name as synonomous with the latest in spamware. I noticed the problem yesterday morning, and took the site offline. Buried in the WordPress .php scripts, I found a piece of code that looked like this:

the bugs have been swept out...

Luckily, the MySQL database seems to have been unaffected, so I did an rm -rf * and started from scratch with the latest WordPress.

It’s been a rather apocalyptic-themed week: Russian hackers attack oklo.org, the University of California is disentigrating under the weight of repeated budget cuts, and on Tuesday, I went to Los Angeles to film a segment for a History Channel episode describing how the Earth would fare in the sudden absence of human presence. My particular interview focused on what would happen to the geostationary satellites over a timescale of weeks to months to years. The filming was done at an abandoned hospital, which was one of the creepiest places I’ve ever seen.

lithium

Mysterious

Diamond prospecting proceeds through the identification of indicator minerals such as specific forms of garnet. The garnets can be traced upstream to the Kimberlite pipes. The Kimberlite pipes contain the sparkling gemstones.

Planet prospecting can be done in similar fashion. If you want to jump-start a new planet search, it’s wise to observe metal-rich stars. Stars with more than twice the Sun’s metal abundance are roughly five times more likely than average to harbor one or more planets in the readily detectable hot Jupiter and Eccentric Giant categories. Histogrammed data from Exoplanet.eu shows the metallicity correlation quite nicely:

Planet Metallicity Correlation

The metallicity correlation can be readily interpreted in the context of the core-accretion paradigm for giant planet formation. In this picture, nascent planets reach the stage of rapid gas accretion when their rocky-icy cores grow to somewhere in the neighborhood of ten Earth masses. The speed with which a core can be assembled in a protoplanetary disk is a very sensitive function of the density of solid material (e.g. ices and dust) in the disk. The density of solids, in turn, scales with metallicity.

If one explains the planet-metallicity correlation with the core-accretion theory, several predictions follow almost immediately. One expects that low-mass stars will show a paucity of readily detectable giant planets, and that high-mass stars will have a larger fraction of giant planets. Observationally, both of these trends have been shown to hold.

A less-well-known prediction is that one also expects that stars with high oxygen (and by proxy, silicon) abundances relative to iron will also show increased planet fractions at given metallicity. Sarah Dodson-Robinson showed this was true as part of her Ph.D. Thesis. Here’s the the key diagram from her paper on the topic:

Silicon-Planet Correlation

A very interesting paper came out in Nature this week which shows an equally compelling, but significantly harder-to-understand abundance correlation. Garik Israelian, and colleagues that include members of the Geneva Team, write (italics are mine):

Here we report Li abundances for an unbiased sample of solar-analogue stars with and without detected planets. We find that the planet-bearing stars have less than one per cent of the primordial Li abundance, while about 50 per cent of the solar analogues without detected planets have on average ten times more Li.

Here’s the graphic from their paper. The filled red circles are planet-bearing stars. Downward arrows indicate that the measurement is an upper limit, and in all likelihood lies at a lower value. Note also, that the y-axis has a logarithmic scale, which de-emphasizes the strength of the effect. To the eye, it’s clear that the lithium abundances of the planet-bearing stars are quite low:

Lithium-Planet Anticorrelation

The effect is dramatic, and yet its origin is mysterious and seems to have gone unpredicted. It’s the best sort of scientific puzzle. Lithium is a rather fragile element, and undergoes nuclear fusion in a star when the temperature reaches ~2.5 million degrees. Lithium depletion in the atmosphere of a star can thus be taken as evidence that the gas that’s currently at the surface has, at one point, been mixed far down enough into the star for the lithium to have burned. This implies that the base of the star’s convective envelope has dipped further into the star than the 2.5 million degree isotherm. (The hot F-type stars on the far right of the diagram have very thin convective envelopes nearly right from the start, and so have been unable to burn their lithium.)

So it seems that somehow, the presence of a planetary system (and even one as wimpy as our own solar system) is enough to alter the evolution of the stellar convective envelope. This, in turn, likely has something to do with angular momentum transfer mediated by planets, but quite frankly the story isn’t very clear. Certainly, there will be papers that explain the effect, and certainly, they are being cranked out even as I write, but unless they make specific, testable, and preferably startling predictions, I’d advise taking them with a grain of lithium chloride.

It’s 5 pm somewhere

Grant-proposal season puts a crimp on one’s style. Despite many interesting developments in the field over the past few weeks, I haven’t had time to write. I’m glad that’ll change shortly.

We’re also very close to getting upgraded versions of the systemic backend and a new Transitsearch-related project on line. In the interim, here’s a link to the old transitsearch.org candidates page. I have it running on our server here at UCO/Lick, and it’s updated every 10 minutes. This information should also soon be available at JPL’s NStED site.

Q: What’s Jupiter’s Q?


With the flood of detail from extrasolar planets, one can forget that our knowledge of the worlds in our own solar system is literally centuries ahead of what we know about planets orbiting other stars. For example, careful naked-eye observations can be used to derive better orbital models for Venus et al. than we currently possess for any exoplanet (assuming, of course, that one owns a good watch and eyesight sufficient to resolve the disk of Venus when it transits the Sun). One of the best ways to learn about what’s out there is to learn as much as we can about what’s right here.

In this vein, an important paper came out in Nature last summer, in which Lainey et al describe a direct and unprecedentedly accurate measurement of the present value of Jupiter’s tidal quality factor, Q. The tidal quality factor encapsulates the ability of an object to dissipate disturbances raised by tidal gravity. The lower the Q, the more capable is the body at damping out the perturbations generated by tidal forcing. Q can depend quite sensitively on the frequency at which perturbations occur, and with a few notable exceptions (for example, the Earth and the Moon), it is notoriously tricky to determine. Previous estimates for Jupiter’s Q ranged from Q~60,000 to over a million. By extension, Q values for Jupiter-mass extrasolar planets are often assumed to lie in this range.

In order to directly measure the Jovian Q, Lainey et al. adopted a procedure that’s conceptually very similar to what goes on inside the systemic console. They first collected measurements of the positions of the galilean satellites that were obtained from 1891 all the way through 2007. They then constructed an orbital N-body model that includes the full gravitational forces acting on Jupiter and the galilean satellites, and which incorporates the non-axisymmetric gravitational pulls exerted by the tidal bulges of Jupiter and Io. The fitted parameters — that is, the initial conditions and undetermined constants — for their model are the osculating orbital elements of the moons, and the values of Q/k2 for Jupiter and Io. (The Love number, k2, is a measure of the degree of central concentration of a body, and has a value of k2~0.37 for Jupiter. For more, see these posts, one, two, from last summer).
Lainey et al. varied the parameters and repeatedly carried out new integrations until the the agreement between where the integrated orbital model said the moons should be located and where they were actually observed was optimized. For this type of direct integrations, goodness-of-fit is highly sensitive to the amount of tidal dissipation in Io and in Jupiter — the larger the dissipation, the larger the effect on the orbit. As a consequence, when a best-fit orbital model is attained, one has direct estimates for the Q‘s of both Jupiter and Io.

And the result? The integrations suggest that the current value of Jupiter’s Q is of order 30,000. This suggests that Jupiter is much more dissipative than has been assumed, and is indeed quite comparable to Neptune or Uranus in terms of its ability to damp out tidal disturbances. The measured Q is low enough, in fact, to suggest that Jupiter currently lies in a state where the tidal forcing by Io is leading to a historically large rate of dissipation. Over the past several billion years, as the orbital frequencies of Io, Europa and Ganymede evolved through a range of values, Jupiter’s Q was on average likely quite a bit higher than it is now.

Jupiter’s low Q hints that the transiting Neptune-mass planet Gliese 436b is even more mysterious than previously though. Gliese 436b has a significantly eccentric orbit whose non-circular figure can only be understood if (1) there’s a suitably influential perturber in the system, or (2) there was a relatively recent disaster, or (3) if the planetary Q has somehow stayed anomalously high through billions of years of orbital evolution. No matter which one of these possibilities turns out to be correct, it’ll be a very interesting story.

the last first look

As is usually the case, there’s been little or no shortage of interesting developments in the field of extrasolar planets. The biggest recent news has been the announcement at the Barcelona conference of a definitive mass for the ultra-short period transiting planet CoRoT-7b. It weighs in at a mere 4.8 Earth Masses (copy of the Queloz et al. preprint here).

Recall that CoRoT-7b caused quite a stir earlier this year with its weird properties. The planet’s year is a fleeting twenty hours and twenty nine minutes, and it induces a tiny transit depth of 0.03%. Unfortunately, the parent star presents a less-than-ideal target for high-precision radial velocity work. It has spots that come and go, and its stellar activity produces frustratingly noisy Doppler measurements. As a result, at the time of CoRoT-7b’s initial announcement, there was no definitive measurement of the planet’s mass.

That’s changed, however, with an unprecedentedly all-out deployment of the HARPS spectrograph. From the Queloz et al. preprint:

A total of 106 measurements between 30 and 60 minute exposure time each were obtained over 4 months, and with sometimes 3 measurements being taken on the same night.

Now in my notoriously biased opinion, such observational enthusiasm is perhaps best reserved for stars such as Alpha Cen B, but a fair argument can be made that the massive investment of time did pay off. Remarkably, the radial velocity data set shows that there are two short-period planets in the CoRoT-7 system. The outer companion, which doesn’t transit, has a period of 3.7 days and at least eight Earth masses. Most dramatically, by combining the mass and radius measurements of CoRoT-7b, one arrives at a density of 5.5 grams per cubic centimeter, essentially identical to that of the Earth, suggesting that the planet is largely composed of refractory materials. (I hesitate to apply the term “rocky” to the CoRoT-7c landscape for the same reason that I’d refrain from describing the Amazon Delta as “icy”.)

In a very real sense, the HARPS campaign on CoRoT-7b has given us our last first look at a fundamentally new category of planet — that is, a world lying in the factor-of-fourteen mass gap spanned by Earth and Uranus. And, from exo-political point of view, the stakes surrounding this discovery were very high. The first density measurement of a planet in this category could just as easily have been made by teams combining high-precision Doppler measurements with either (1) Warm Spitzer, (2) ground-based photometry, (3) Kepler, (4) MOST, (5) HST, or (6) CoRoT. So I can imagine that there was a certain impetus underlying the scheduling of that huge block of HARPS time.

The discovery could, however, still be waiting to be made. Despite all the effort with HARPs, there remains a hefty 70% error on the density determination. This means that there’s a ~16% chance that CoRoT-7b is actually less dense than Neptune.

I’ll go out on a limb: CoRoT-7b’s density will turn out to be anomalously high. More than 90% of “super Earths” will turn out to be “sub-Neptunes” as far as their density is concerned.

campaign mode

Full-resolution Poster-sized .pdf of the above.

The next HD 80606 transit is coming up this week. While the sky position of the star will be much more favorable during the coming January event, observers across the US have an opportunity to get photometric measurements of the ingress early Thursday morning.

The transit begins just after 11 AM UT on Sept. 24, and will unfold over the next 12 hours, meaning that observers in Japan and East Asia will be able to catch the egress.

Josh Winn of MIT is organizing a repeat of the successful June campaign (detailed in this post). If you’re a capable photometric observer, and if you’re interested in participating in the campaign, definitely get in touch with him.