Astronomy & Celestial Mechanics
the sky, computed
Generator

The family-cluster generator

Six families and four gaps, in coordinates no telescope reports
Six families and four gaps, in coordinates no telescope reports. The main belt in proper elements: semi-major axis against proper sin i, with 620 background bodies and the members of six families. Nothing here is what an orbit catalogue lists. Every coordinate is a free-oscillation amplitude recovered by averaging away the short-period terms and subtracting the forced ones, and the same objects plotted in the elements they are catalogued with fill this frame as a single smooth cloud. What appears once the averaging is done is structure of two opposite kinds. The clumps are collisions: each is the debris of one parent body, and each drawn family scatters between 5.8 and 12.3 times more tightly than the background at its own distance from the Sun. The gaps are resonances — the vertical lanes at 3:1, 5:2, 7:3, 2:1 with Jupiter — and they are empty because proper elements stop existing there: a body in a mean-motion resonance has no quasi-integral to compute, its eccentricity wanders, and it leaves. The two features are the same argument seen twice: the belt keeps a memory of individual events for hundreds of millions of years everywhere except where the memory is erased, and both the memory and the erasure are invisible until the coordinates are changed.

The main belt in proper elements: semi-major axis against proper sin i, with 620 background bodies and the members of six families. Nothing here is what an orbit catalogue lists. Every coordinate is a free-oscillation amplitude recovered by averaging away the short-period terms and subtracting the forced ones, and the same objects plotted in the elements they are catalogued with fill this frame as a single smooth cloud. What appears once the averaging is done is structure of two opposite kinds. The clumps are collisions: each is the debris of one parent body, and each drawn family scatters between 5.8 and 12.3 times more tightly than the background at its own distance from the Sun. The gaps are resonances — the vertical lanes at 3:1, 5:2, 7:3, 2:1 with Jupiter — and they are empty because proper elements stop existing there: a body in a mean-motion resonance has no quasi-integral to compute, its eccentricity wanders, and it leaves. The two features are the same argument seen twice: the belt keeps a memory of individual events for hundreds of millions of years everywhere except where the memory is erased, and both the memory and the erasure are invisible until the coordinates are changed.

5 essays call family-cluster. The drawing above is what it returns with no arguments at all; every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

Where it is called

Every figure listed here is the same construction drawn at different numbers, so a correction to one is a correction to all of them.

A 128.8-million-year-old collision, dated from the shape of a scatter plot. The Erigone family: 165 members drawn at their diameters and their proper semi-major axes, with inverse diameter up the page. The cloud is a V, and the V is a clock. Each member has been drifting in semi-major axis ever since the collision at a rate that goes as one over its diameter, with a sign set by which way it spins — prograde outward, retrograde inward — so after 130 million years the small members have moved far and the large ones have barely moved at all. Plotted against 1/D that envelope is a straight line through the family's centre, and its slope is the drift rate for a one-kilometre body multiplied by the elapsed time. Fitting the two edges of the points actually drawn here returns 128.85 million years against the 130 the members were generated from. The rounding at the bottom is not an artefact: it is the ejection velocity, some 15 metres per second, which every member got at the moment of the collision and which is the same for all sizes. The picture cannot show the interlopers — background asteroids that happen to lie inside the V and have nothing to do with the family — and it cannot show the members that have drifted into a resonance and left the belt entirely, which is the reason the oldest families have the softest edges. Orbits

A collision dated by a scatter plot

Nothing in the solar system carries a date. A collisional family does — because a force that depends on a body's size has been pushing its fragments apart ever since, so the cloud is a V whose slope is an elapsed time, and one of those dates is confirmed by fossil meteorites in Swedish limestone.

 Mass at the top, area at the bottom, 10 decades apart. Two moments of a collisional cascade's size distribution, per logarithmic interval of diameter, over 10 decades from a ten-micron grain to a hundred-kilometre parent body. Both axes are logarithmic and the vertical scale is arbitrary; only the slopes carry the argument. A population in which every collision makes fragments that go on to collide reaches a steady state where the same mass flows through every size per unit time, and that fixes the differential number distribution at an index of 3.5. The two consequences pull opposite ways. Mass per decade goes as the diameter to the power 0.5, so it climbs and almost all the mass is in the largest few bodies. Cross-sectional area per decade goes as the diameter to the power -0.5, so it falls, and almost all the area — which is what scatters light, what is detected, and what anything passing through gets hit by — is in the smallest. Over the range drawn the small end carries 10⁵ times the area of the large end and 10⁻⁵ times its mass. A disc's brightness therefore measures a population whose mass it says nothing whatever about, and the two numbers are connected only through the index of this line. Orbits

Where the mass is and where the light is

A population that grinds itself up settles into a size distribution with a fixed slope, and that slope puts almost all the mass in the largest bodies and almost all the cross-section in the smallest. So a debris disc's brightness measures a population whose mass it says nothing about, and the two are connected only by the exponent.

 A branch that adds nothing above 1 km and everything below it. Cumulative crater counts on a 3.5-billion-year-old surface, with the population split into the craters made by objects arriving from outside and the craters made by blocks thrown out of larger ones on the same surface. The two are indistinguishable in a photograph and completely different as a statistic. Secondaries stop at about 1 kilometre, because that is the largest crater a block leaving at a few hundred metres a second can excavate, so the upper half of the plot is unaffected. Below it they are steeper — slope -3.2 against the primaries' -2 — and by the smallest diameter drawn they outnumber the primaries 292 to one. A count taken at 100 metres and read through the primary production curve returns an age of 4.34 billion years for ground that is 3.5, and it returns it with a small formal error, because the counting statistics are excellent. The error is not in the counting. Orbits

The craters that were not primary

Counting craters dates a surface, and the method works because impacts from space arrive at a known rate. Some of the holes were not made from space. They were made by rock thrown out of the larger holes on the same surface, and they are far more numerous than anything that arrived.

 Below 0.46 microns a grain is not in orbit at all. The ratio of the radiation force to the gravitational force on a dust grain, against the grain's radius, for three densities. Every line has slope exactly −1 because gravity acts on the mass and radiation on the cross-section, and the ratio of a volume to an area is a length. Two horizontal lines matter and they are different statements. At β = 1 the star does not attract the grain at all. At β = 1/2 a grain released at rest from a circular orbit is already unbound, because it keeps the speed appropriate to the full stellar mass while feeling only half of it — and since dust is made by breaking up larger bodies that were on circular orbits, the lower line is the one that applies. For rock at 2500 kilograms a cubic metre that is 0.46 microns; for ice it is 1.15, and for iron 0.15. A collisional cascade that grinds material finer runs into this floor and stops, and the material that would have been finer leaves the system on a hyperbola. Orbits

The drag that sorts a disc by size

Starlight does three different things to a dust grain depending on how big it is — blows it out of the system, drags it inward over millennia, or ignores it entirely. Which one happens is decided by a single length, and whether it happens at all is decided by how crowded the disc is.

 A family of 276, or of 582, depending on one number. Membership of a synthetic family against the cutoff velocity used to define it, for a family of 260 genuine fragments ejected at 15 metres a second sitting in a background of 417 unrelated bodies. The clustering is single linkage in the standard proper-element metric, started from one object and grown until nothing more is within the cutoff of anything already absorbed. Three curves: the number of objects claimed, scaled to its largest value; the fraction of the real family recovered; and the fraction of the claim that is background. The shape is the whole difficulty. At low cutoff the family is fragmented and only its core is found. There is then a plateau — near 58 metres a second here, giving 276 members of which 6 per cent are background — and that plateau is what every published family list is chosen at. Past it the count runs away, because single linkage absorbs an object and then searches from the object it has just absorbed: one chance interloper bridges the family to the belt and the algorithm returns half the main belt. Completeness and contamination rise together, so there is no cutoff at which both are good and the plateau is a compromise rather than a discovery. Orbits

A family whose size is a choice

Nothing observable distinguishes a collisional fragment from an asteroid that happens to be nearby. Membership is assigned by clustering at a cutoff velocity, and the same family has 276 members or 582 according to which cutoff is used — with completeness and contamination rising together, so no cutoff is good.

The whole library · All essays