Archive for August, 2013

Central Limit Theorem

August 21st, 2013 Comments off


We’re putting the finishing touches on a new research paper that deals with an old favorite: HD 80606b. The topic is the Spitzer Telescope’s 4.5-micron photometry taken during the interval surrounding the planet’s scorching periastron passage, including the secondary eclipse that occurs several hours prior to the moment of closest approach (see the diagram just below). I’ll write a synopsis of what we’ve found as soon as the paper has been refereed.


In writing the conclusion for the paper, we wanted to try to place our results in perspective — the Warm Mission has been steadily accumulating measurements of secondary eclipses. There are now over 100 eclipse depth measurements for over 30 planets, in bandpasses ranging from the optical to the infrared.

A set of secondary eclipse measurements at different bandpasses amount to a low-resolution dayside emission spectrum of an extrasolar planet. When new measurements of secondary eclipse depths for an exoplanet are reported, a direct comparison is generally made to model spectra from model atmospheres of irradiated planets. Here is an example from a recent paper analyzing Warm Spitzer’s measurements of WASP-5:


Dayside planet/star flux ratio vs. wavelength for three model atmospheres (Burrows et al. 2008) with the band-averaged flux ratios for each model superposed (colored circles). Stellar fluxes were calculated using a 5700 K ATLAS stellar atmosphere model (Kurucz 2005). The observed contrast ratios are overplotted as the black circles, with uncertainties shown. The model parameter kappa is related to the atmosphere’s opacity, while p is related to the heat redistribution between the day and night sides of the planet (Pn = 0.0 indicates no heat redistribution, and Pn = 0.5 indicates complete redistribution).

As is certainly the case in the figure just above, the atmospheric models that are adopted for comparison often have a high degree of sophistication, and are informed by a substantial number of free parameters and physical assumptions. In most studies, some of the atmospheric parameters, such as the presence or absence of a high-altitude inversion-producing absorber, or the global average efficiency of day-to-night side heat redistributions are varied, whereas others, such as the assumption of hydrostatic equilibrium and global energy balance, are assumed to be settled. Invariably, the number of implicit and explicit parameter choices tend to substantially exceed the number of measurements. This makes it very hard to evaluate the degree to which a given, highly detailed, planetary atmospheric model exhibits any actual explanatory power.

The central limit theorem states that any quantity that is formed from a sum of n completely independent random variables will approach a normal (Gaussian) distribution as n becomes large. By extension, any quantity that is the product of a large number of random variables will be distributed approximately log-normally. We’d thus expect that if a large number of independent processes contribute to a measured secondary eclipse depth, then the distribution of eclipse depth measurements should be either normally (or possibly log-normally) distributed. The “independent processes” in question can arise from measurement errors or from systematic observational issues, as well as from the presence of any number of physical phenomena on the planet itself (such as the presence or absence of a temperature inversion layer, or MHD-mediated weather, or a high atmospheric C/O ratio, etc.).

The plot just below consolidates more than 100 existing secondary eclipse measurements onto a single diagram. Kudos to for tracking the secondary eclipse depths and maintaining a parseable database! The observed systems are ordered according to the specific orbit-averaged flux as expressed by the planetary equilibrium temperaturs — the nominal black-body temperature of a zero-albedo planet that uniformly re-radiates its received orbit-averaged stellar energy from its full four-pi worth of surface area. The secondary eclipse depths in the various bands are transformed to flux ratios, F, relative to what would be emitted from a black-body re-radiator. If all of the measurements were perfect, and if all of the planets were blackbodies, all of the plotted points would lie on the horizontal line F=1.


It’s somewhat startling to see that there is little or no systematic degree of similarity among the measurements. One is hard pressed to see any trends at all. Taken together, the measurements are consistent with a normal distribution of flux ratios relative to a mean value F=1.5, and with standard deviation of 0.65:


This impression is amplified by the diagram just below, which is a quantile-quantile plot comparing the distribution of F values to an N(0,1) distribution.


The nearly gaussian distribution of flux ratios suggests that the central limit theorem may indeed find application, and imparts a bit of uneasiness about comparing highly detailed models to secondary eclipse measurements. I think we might know less about what’s going on on the hot Jupiters than is generally assumed…

Categories: worlds Tags:


August 10th, 2013 3 comments


One prediction regarding exoplanets that did hold true was the Moore’s-Law like progression toward the detection of planets of ever-lower mass. More than seven years ago, not long after the discovery of Gliese 876 d, the plot of Msin(i) vs. year of discovery looked like this:


With a logarithmic scale for the y-axis, the lower envelope of masses adhered nicely to a straight line progression, pointing toward the discovery of the first Earth-mass exoplanet sometime shortly after 2010. The honors went, rather fittingly, last year, to Alpha Cen B b. Here’s an update to the above plot. Planets discovered via Doppler velocity only are indicated in gray, transiting planets are shown in red…


The data for the plot were parsed out of the very useful exoplanets.csv file published at

And wait, what’s going on with that point in 1993? See

Categories: worlds Tags:


August 6th, 2013 7 comments


I think it’s worth making an attempt to coin a term for these “ungiant” planets that are, effectively by default, largely being referred to as super-Earths, a term which brings to mind Voltaire’s remark regarding the Holy Roman Empire.

Planets in the category:

1. Have masses between ~1%  and ~10% of Jupiter’s mass.
2. Have unknown composition, even if their density is known.


Ideally, a term for such planets would:

3. Have a satisfying etymology springing from the ancient Greek.
4. Not be pretentious, or, much more critically, not be seen as being pretentious.

Simultaneously satisfying conditions 3 and 4 is certainly not easy, and indeed, may not be possible. (See, e.g.,

I’ve noticed that the esoteric efforts to describe the interiors of these planets — in the absence of any data beyond bulk density — effectively boil down to Robert Fludd’s 1617 macrocosm of the four classical elemental spheres:


This led me to look into Empedocles’ four elements themselves, see, e.g., here. Specifically, can a term of art for the planets of interest be constructed from the original Greek roots?

The following table on p. 23 of Wright, M. R., Empedocles: The Extant Fragments, Yale University Press, 1981, contains various, possibly appropriate, possibilities:


To get going, I had to refer to the rules for romanization of Greek. Initial attempts to coin names (while abundantly satisfying requirement #3 above) have so far failed miserably on requirement #4: chonthalaethian planets, ambroaethic planets, gaiapontic planets. Yikes!

The Tetrasomia, or Doctrine of the Four Elements, alludes to the secure fact that these planets are unknown compounds of metal, rock, ices, and gas. Tetrian planets, maybe? Suggestions welcome…

Categories: worlds Tags: