Every method prefers a circle, and not for the same reason
Assumes Detection bias and The ellipse.
The eccentricity distribution of planets is one of the most argued-about numbers in the subject, because it is the clearest surviving trace of how a system was assembled. A planet that arrived by high-eccentricity migration must have spent time on a long ellipse; one that spiralled in through a disc need never have been anything but round.
So it matters a great deal whether the measured distribution is the real one. It is not, and the corrections are three separate effects that do not point the same way.
The probability rises
A planet transits if its orbit carries it across the stellar disc from here, and for a circular orbit the chance of that is — the fraction of viewing directions from which the geometry works.
An eccentric orbit does not have one distance. It spends part of its period inside its own semi-major axis and part outside, and the transit probability depends on the distance at conjunction, which is where the planet happens to be when it crosses the line between here and the star. Averaged over the orientation of the ellipse in its own plane, the probability comes to
which is larger than the circular value, and by a factor that runs away as the eccentricity approaches one.
The reason the average goes that way rather than the other is worth pausing on, because the intuition that an eccentric planet spends most of its time far out is correct and points the wrong direction. The planet does spend most of its time near apoapsis — that is Kepler’s second law — but the transit probability does not average over time. It averages over the orientation of the orbit, which is uniform, and the periapsis passage contributes a large probability from a narrow range of orientations.
The duration falls
A transit’s depth is a ratio of areas and does not depend on the eccentricity at all. Its duration does, and it depends on it twice.
An eccentric planet at conjunction is at a distance that differs from , so the chord it draws across the star differs in length. And it is moving at a speed that differs from the circular speed at that distance, by the amount vis-viva requires. Combining the two, the duration relative to a circular orbit of the same semi-major axis is
and averaged over the argument of periastron the denominator drops out, leaving .
A shorter transit carries less signal. The detection statistic scales as the square root of the time in transit, so a transit half as long has seventy per cent of the signal-to-noise, and near a threshold that is the difference between a detection and nothing.
Which very nearly cancel
Put the two together — the probability of a transit happening, times the signal-to-noise it delivers when it does — and the product at is 1.24 relative to a circle. At it is 1.35.
That is a remarkably small bias for a quantity that runs over a factor of three in each of its two parts. A transit survey, averaged over orientations, has almost no preference between round orbits and eccentric ones, and the cancellation is not approximate for any deep reason — it happens because one effect goes as and the other as after the square root of the duration is taken, and the residual is mild over the range that matters.
Why the cancellation is not a coincidence worth trusting
It is tempting to read the near-cancellation as a piece of luck that lets transit surveys off, and the reading is half right.
The two effects have a common origin: both come from the fact that an eccentric orbit’s conjunction distance differs from its semi-major axis. The probability goes as the inverse of that distance and the transit’s chord-crossing speed goes as roughly its inverse square root, so the exponents are and before the signal-to-noise takes another root. Written that way the partial cancellation is structural rather than accidental.
What is accidental is how complete it is, and that part depends on the noise. The exponent connecting duration to detectability is for white noise and for noise correlated on the transit timescale, and real light curves sit between them. So the residual bias is a small number whose size is set by an instrumental property nobody quotes.
A small residual whose magnitude is unknown is a worse situation than a large one that is known. A factor of three can be corrected; a factor of cannot be, because the uncertainty on the correction is comparable to the correction, and it propagates into every moment of the recovered eccentricity distribution.
The one that does not cancel
Radial velocity is a different situation entirely, and the difference is that the signal is a shape rather than an event.
A circular orbit produces a pure sinusoid in the star’s line-of-sight velocity. An eccentric one does not: the velocity curve is skewed, with a rapid excursion through periastron and a long slow return, and Fourier analysis puts a growing fraction of its variance into the second harmonic and beyond.
A search that fits a sinusoid — which is what a periodogram is — sees only the fundamental. So the amplitude it recovers is smaller than the true semi-amplitude by the square root of the fraction of power remaining in the first harmonic, and that fraction falls steadily with eccentricity.
The fraction is 68 per cent at and 47 per cent at , which is a loss of 18 and 31 per cent in recovered amplitude. Near a detection threshold that is decisive, and unlike the transit case there is nothing pushing the other way.
A full Keplerian fit does better than a periodogram, because it is fitting the right shape. It pays for that in two extra free parameters — the eccentricity and the argument of periastron — and two extra parameters cost significance in any model comparison. So even a search that fits the correct model requires a higher signal-to-noise to claim an eccentric planet than a circular one of the same amplitude.
The bias that is a geometry and the bias that is a search
It is worth setting the three effects against each other by kind, because only two of them are correctable in principle.
The transit probability is a geometric bias. It depends on nothing but the orbit and the observer’s position, it is exact, and correcting for it needs no knowledge of the instrument at all. Every survey that has ever existed shares it, and it will not change.
The transit duration is an instrumental bias. It becomes a bias only because a shorter event carries less signal, and how much less depends on the noise. A survey with infinite sensitivity would have no duration bias whatever.
The harmonic loss is a methodological bias. It exists because a particular search algorithm looks for a particular shape, and a different algorithm has a different bias or none. A search that fitted full Keplerians on a grid of eccentricity would not lose the harmonics; it would lose significance to the extra parameters instead, which is a different penalty of a similar size.
Only the first is a property of the universe. The other two are properties of a decision somebody made, and they change between surveys, between pipelines, and between reanalyses of the same data — which is why two eccentricity distributions derived from the same catalogue by different groups can differ, without either being wrong about the catalogue.
What the corrected distribution says
Correcting for all of this is what turns the catalogue into a statement about formation, and the statement is worth having.
Giant planets on periods beyond about ten days have a broad eccentricity distribution, running from zero to above 0.9, with a median near 0.25. That is nothing like a disc-driven population, which would be nearly circular, and it is the strongest evidence that gravitational interaction between planets — scattering, or secular exchange with a distant companion — is a common part of a system’s history.
Inside ten days the distribution collapses to zero, and that wall is a tidal clock rather than a formation signature.
And the small planets, the ones the transit surveys find in abundance, are a third case: eccentricities that are small but not zero, typically below 0.1, and systematically smaller in systems with more planets than in systems with one. That last correlation is the interesting one, and it is close enough to the detection biases to be worth suspicion — a survey finds more planets per system when the system is flat and calm, which is a selection effect on multiplicity, and flat calm systems are also the ones with low eccentricities.
An asymmetry nobody put there
There is a third consequence of the orientation dependence, and it is the one that produces an effect in a quantity that ought to have none.
The transit probability depends on through , and so does the duration. Both favour periastron pointing towards the observer. So a transit survey does not merely over-represent eccentric planets slightly — it over-represents eccentric planets in a particular orientation, and the orientation is defined relative to the line of sight.
That is a selection effect on a quantity that has no physical meaning. The argument of periastron measured from the sky plane is an accident of where the Earth happens to be, and any real distribution of it must be uniform. So a measured distribution of that is not uniform is a pure measurement of the selection function, with no astrophysics in it whatever, and it can be used as a check on everything above.
The transiting planets do show it. The distribution of the argument of periastron in the transiting sample is skewed towards periastron-on-the-near-side, by about the amount the geometry predicts, and the agreement is one of the few places where a detection bias has been confirmed rather than modelled.
There is one confounder, and it is the orbital decay of the very hottest planets, where a tidal argument also produces an orientation preference by circularising the systems whose periastra come closest. The two effects are separable by period, and the check is made outside the range where tides act.
Exact geometry, an honest quadrature, and the thing neither supplies
The transit probability and duration expressions are geometry and Kepler’s laws, exact within the two-body problem, with no free parameters. They are not approximations in any sense that matters here.
The harmonic decomposition is computed from the velocity curve itself — Kepler’s equation solved at four thousand mean anomalies, the true anomaly taken from the eccentric one, the resulting curve transformed, and the power in the fundamental compared against the total by Parseval’s relation rather than against a series expansion. At it returns exactly one, which is the check that the solver and the transform agree.
What is not measured is the thing the corrections are for, and this is the honest limit of the whole exercise. The detection efficiency of a real pipeline against an eccentric planet is not a formula; it is the same injection experiment that measures a pipeline’s efficiency, run with eccentric signals, and it has been done far less often than the circular version.
Fixed semi-major axis, and what that choice hides
Every curve is drawn at fixed semi-major axis, which is a choice and not the only one. Compared at fixed period the same statements hold, since period and semi-major axis are related by a law with no eccentricity in it. Compared at fixed periastron distance — which is what a tidal argument would want — they do not.
The transit case ignores impact parameter. An eccentric orbit changes the geometry of the crossing as well as its duration, and a grazing transit of an eccentric planet has a different shape from a grazing transit of a circular one. That affects what can be measured from the transit rather than whether it is found.
And nothing here shows the eccentricity distribution itself. The figures give the bias against each eccentricity; converting that into a corrected distribution requires the true distribution’s shape, which is what the exercise is trying to find, so the correction has to be done inside the inference rather than applied to its output.
What a bias in a shape does to a mass
One more consequence deserves separating, because it changes a number rather than a count.
A velocity semi-amplitude gives a minimum mass, and the relation between the two contains the eccentricity: at fixed mass and period, rises as . So an eccentric planet produces a larger velocity amplitude than a circular one of the same mass — which is a bias towards detecting it, running opposite to the harmonic loss.
The two do not cancel, because they act at different stages. The amplitude gain is in the signal; the harmonic loss is in the search. A periodogram searching for sinusoids sees the fundamental of a larger curve, and the product of the two factors is , which stays near one up to and then falls.
So radial velocity has its own near-cancellation, over a narrower range, and the interesting regime is above it. Beyond the harmonic loss wins decisively, which is precisely where the most dynamically informative planets are.
And the mass that comes out of a fit in which the eccentricity is poorly constrained inherits that. Fitting a low-signal eccentric orbit with the eccentricity free returns a posterior that is broad in and correlated with , so the mass and the eccentricity are measured together or not at all — which is the same coupling a transit timing amplitude has, arriving from a completely different observable.
The general shape
The recurring structure is that a selection effect is only correctable when it is known, and the three effects here are known to very different degrees.
Geometry is known exactly. A transit probability is a solid angle and a duration is a chord, and neither depends on anything a survey does. Those corrections can be applied with confidence.
Pipeline response is known only by experiment, and the experiment has usually been run for the circular case alone.
And the correlation between eccentricity and everything else is not known at all. Eccentric planets are found around different stars, at different periods, in different multiplicities — so a correction applied to eccentricity alone is applied to a variable that is not independent of the others. A duration measured against a transit’s expected length gives the eccentricity directly for an individual system, which is the way out, and it works only where the stellar density is independently known.
Why this is not a small correction to a large number
A reader could reasonably conclude from the near-cancellations above that the eccentricity biases are a second-order concern. For the mean of the distribution they nearly are. For its tail they are not, and the tail is the part that carries the physics.
The distinguishing prediction of planet–planet scattering is not a shifted mean; it is the existence of planets at , which disc migration cannot produce at all. The number of such planets in a catalogue is small, the correction to it is large, and a factor of two in the velocity method’s efficiency at moves the inferred fraction of scattered systems by a factor of two.
So the quantity most sensitive to the bias is the one the measurement exists to determine, which is the usual arrangement and is worth stating each time it recurs. The safe part of a measurement is the part nobody needs.
There is a partial escape, and it is the reason the two methods are worth running together. The transit sample’s bias at high eccentricity is mild and the velocity sample’s is severe, so a planet population measured both ways gives a consistency check on the correction: if the corrected distributions agree, the corrections were probably right, and if they disagree, the disagreement localises which method’s correction failed. That check has been run, and the two agree within their uncertainties for giant planets and have too few objects in common to compare for small ones.
Still open: the multiplicity that travels with it
The correlation between low eccentricity and high multiplicity is the loose end. It is a real pattern in the catalogue and it has at least two readings, one astrophysical and one instrumental, and they predict the same thing.
The astrophysical reading is that a system with many surviving planets is one that never scattered, so its planets are both numerous and round. The instrumental reading is that a survey finds many planets in a system only when the system is flat — and a flat system is a dynamically cold one, which is also a circular one.
Separating them requires knowing how a survey’s multiplicity is biased, and that turns out to be the hardest quantity a catalogue reports.
About the same objects
Not linked from either essay — found by the objects both name.
- One timing curve and five planets that could draw it radial velocity · transit probability
- The only stars whose masses are known radial velocity · selection effect
- The planets that were not seen survey completeness · transit probability
What links here
Essays that link to this one from their own argument.
The objects this essay names
Each one links to every other essay that touches it.
Argument of periastronDetection thresholdKepler's equationOrbital eccentricityRadial velocitySelection effectSignal-to-noiseSurvey completenessTransit durationTransit probability