Astronomy & Celestial Mechanics
the sky, computed
Concept

Globular cluster — where it appears

A gravitationally bound sphere of a hundred thousand to a million old stars, orbiting a galaxy's halo. Its relaxation time is under a Hubble time, which makes it the standard laboratory for mass segregation, evaporation and core collapse.

Named by 18 essays across 3 fields — each of them below, with the objects they name alongside it.

7 decades of relaxation time, and a Hubble time between them. Crossing time and two-body relaxation time for 5 self-gravitating systems, against the number of bodies in each. The lower marks are the crossing time R/σ, which spans 3 decades; the upper ones are t_relax ≈ 0.1N/ln N times it, which spans 8. The horizontal rule is a Hubble time, and the same expression puts these systems on opposite sides of it. Two of them relax — an open cluster and a globular cluster — so their stars have exchanged energy, their heavy stars have sunk, and their present structure is not the one they were born with. The rest do not: a dwarf spheroidal, an elliptical galaxy, a cluster of galaxies, the elliptical missing it by a factor of 10⁶, which is why a galaxy is modelled as a smooth potential with no stars in it at all. The dependence is on N and hardly at all on anything else, because the Coulomb logarithm makes every decade of impact parameter contribute equally. A cluster that relaxes also evaporates: about a stellar mass in a hundred and forty leaves per relaxation time, so the globular cluster empties in 6.1·10¹¹ years.

How long a system takes to forget

A star crossing a galaxy is deflected by every other star, and the deflections add as a random walk. The time for that walk to change a star's energy by its own amount is a hundred million billion years for a galaxy and a billion for a globular cluster, and everything about how the two are modelled follows from which side of a Hubble time that falls on.

gravitation · Two-body relaxation Three clusters, three ages, one diagram. Isochrones for populations of 100 Myr, 1 Gyr, 12 Gyr, each drawn as the main sequence up to its own turnoff and then the post-main-sequence track of the turnoff mass. The turnoffs are at 5.52 M☉, 2.20 M☉, 0.94 M☉, from t = 10¹⁰ M/L with this file's own mass–luminosity relation. Nothing here is a track along which a star moves. Every point is a different star of a different mass, all the same age, and the bend is simply where the population runs out of stars that have had time to leave. That is why a cluster has an age and a field star does not: the bend needs a population, and one star is not one. The oldest globular clusters sit near the 12 Gyr line, and in the 1990s the same construction gave them ages of 16 to 18 Gyr against a universe measured at 10 — a two-standard-deviation contradiction that was resolved from the distance side, by Hipparcos, and not from this one.

An age read off a bend

A cluster is a population of one age and many masses, so the place where its stars leave the main sequence is a clock. The bend needs a population — which is why a cluster has an age and a single star does not.

stars · The HR diagram A stream that is not the orbit it came from. 456 stars released in pairs from the two saddles of a 10⁵ solar-mass cluster over 4.0 billion years, integrated in a halo whose circular speed is 220 kilometres a second, drawn with the progenitor's own orbit. The cluster runs between 10 and 25 kiloparsecs and the orbit is the thin closed-looking curve; the stars are everything else. The point of the figure is the discrepancy. Stars leaving through the inner saddle are on slightly smaller orbits, turning round at a median of 24.77 kiloparsecs rather than the progenitor's 25.00, and therefore running ahead; stars leaving through the outer saddle reach 25.22 and fall behind. The whole spread is 1.8 per cent of the apocentre, which is the number worth carrying: an offset far too small to see in this drawing builds the entire stream, because it acts for four billion years. The two arms are therefore not merely displaced along the orbit, they are on different orbits, and the track a survey measures is a family of them rather than any single one. Fitting a Galactic potential by demanding that a stream lie along an orbit is wrong by exactly this much, and the size of the error grows with the mass of the progenitor, because the mass is what sets the distance between the two doors.

A stream is not the orbit it came from

A cluster torn apart by a galaxy leaves a thin trail of stars across the sky, and the obvious thing to do with it is fit an orbit. That is wrong by a knowable amount, because stars leave through two doors with a small energy offset and end up on a family of orbits rather than on one.

galaxies · Stellar streams Hard below 8.3 astronomical units, soft above, and nothing settles at the line. The binding energy of a binary of two 0.7 solar-mass stars against its separation, both axes logarithmic, with the mean kinetic energy of a single cluster star at a velocity dispersion of 5 kilometres a second drawn as a level. Where the curve is above the level the binary is bound more tightly than a passing star's motion, and encounters on average take energy out of the field and put it into the pair; where it is below, they do the reverse. The crossing at 8.3 astronomical units is the hard–soft boundary, and it is the entire content of Heggie's law: hard binaries harden and soft binaries soften. The arrows are the direction each side moves, and they point away from the crossing in both directions rather than toward it. That is not a coincidence but a negative heat capacity, the same property that makes a star contract when it radiates: taking energy out of a bound pair moves it closer together and speeds it up, so a hard binary that gives energy to the cluster becomes harder still and gives more. A cluster with binaries in it therefore has a heat source that turns itself up, and the boundary drawn here is a watershed rather than an equilibrium.

A pair that heats what is trying to cool

A star cluster has a negative heat capacity, so it cannot reach equilibrium — its core contracts and gets hotter without limit. What stops it is a binary, and which way a binary exchanges energy with the stars around it is settled by one comparison of two energies.

gravitation · Binary heating Dynamical friction is strongest at 0.97 times the dispersion, and weaker either side. The dynamical friction on a body moving through a Maxwellian sea of stars, against its own speed in units of the square root of two times the velocity dispersion. The dashed curve is the fraction of the field moving slower than the body, which is the only part of the field that contributes: a star overtaking from behind pulls the body forward exactly as often as one being overtaken pulls it back, so the fast tail cancels out entirely. That fraction rises from nothing to one. The drag is that fraction divided by the square of the speed, and the quotient of a saturating numerator and a growing denominator has a single maximum, here at 0.97. The consequences run in both directions. A body moving much faster than the stars around it is barely slowed at all, which is why a galaxy passing through a cluster at a thousand kilometres a second does not sink; and a body already at rest with respect to the field feels nothing either, because there is no wake to be behind it. Sinking is therefore fastest in the middle of the process rather than at the start of it.

A drag that is strongest in the middle

Dynamical friction is not like air resistance. Only the stars moving slower than the body contribute at all, so the drag rises from nothing, peaks when the body has slowed to about the speed of everything around it, and falls away again — which is why a fast satellite is never captured and a slow one sinks in a hurry.

gravitation · Dynamical friction Mestel cooling for white dwarfs of 0.4 to 1 solar masses. Luminosity against age for white dwarfs of 0.4, 0.6, 0.8, 1 solar masses, both axes logarithmic. Each line is the whole of the star's later life: there is no nuclear source, so what is radiated is the thermal energy the ions had when the star was made, leaking out through a thin envelope whose opacity is Kramers'. That one sentence gives L ∝ t^(−7/5), and the lines are straight here because a power law is straight on these axes — the slope is the exponent and can be measured off the picture. Two features are worth more than the numbers. The first is how flat the sequence becomes: a factor of ten in luminosity costs a factor of 5.2 in age, so a white dwarf spends most of its existence faint, and any population of them piles up at the bottom. The second is that the heavier tracks lie above the lighter ones. A heavier white dwarf is smaller, so it has more ions to cool and less surface to lose them through, and at a given age it is the brighter object. The picture does not show what happens at the faint end, where the ions crystallise and the heat capacity stops being constant; that is a separate figure, and it is where this law stops being enough.

A clock with no fuel in it

A white dwarf has nothing left to burn, so its brightness is a record of how long it has been cooling. The faintest ones in the Galaxy are not faint because they are small — they are as faint as anything has had time to become, and the place where they stop is a date.

stars · White dwarf cooling Every pair arrives circular. Eccentricity against gravitational-wave frequency for four binaries of 30 and 30 solar masses, each starting at 0.2 astronomical units with an eccentricity of 0.3, 0.7, 0.9, 0.99. Both axes are logarithmic and the tracks run left to right as the orbit shrinks. The curves are Peters' closed solution a(e), and they are checked against Peters' differential equation at three eccentricities on each track rather than against the integral they came from. The ordering is preserved — a pair that starts rounder stays rounder — but by the time the orbit is radiating at 10 hertz, where a ground-based detector begins to hear it, the four eccentricities are 8.8·10⁻⁸, 5.6·10⁻⁷, 3.6·10⁻⁶, 1.5·10⁻⁴. All four are far below anything a detector could measure. That is the figure's whole content and it is a strong statement: radiation reaction removes angular momentum faster, relative to energy, than a circular orbit would need, so eccentricity is destroyed on the way in. A binary observed to be eccentric in band therefore cannot have spent long shrinking quietly, and must have been put on that orbit recently — by a third body, or in the crowded centre of a cluster. The figure assumes the two bodies are points and nothing else acts on them, which is exactly the assumption an eccentric detection would refute.

Every pair arrives circular

Gravitational radiation drains a binary's energy and its angular momentum at rates that do not keep step, and the mismatch destroys eccentricity far faster than it shrinks the orbit. A pair that starts at 0.99 and spirals in from a fifth of an astronomical unit is round to better than a part in a hundred million by the time a detector can hear it.

gravitation · Gravitational waves Four defensible choices, and a factor of 2.9 between them. The Coulomb logarithm against the ratio of the two impact parameters it is cut off at. The curve is a logarithm, so it is flat — a factor of ten in the ratio buys 2.3 — and that is usually offered as the reason not to worry. The four marked conventions are all in current use and all defensible, and they give ln Λ from 3.4 to 9.9. Since the drag force is proportional to ln Λ and not to its logarithm, that is a factor of 2.9 in every sinking time computed from it. The flatness protects the answer from a wrong guess about the ratio; it does not protect it from there being no correct guess, which is the actual situation.

A drag computed with a logarithm nobody can pin down

Chandrasekhar's drag formula is exact, derived from first principles, and contains a logarithm of a ratio of two lengths that the derivation does not supply. Every sinking time in astronomy is proportional to that logarithm, and the four conventions in current use differ by a factor of three.

gravitation · Dynamical friction The kick falls as b⁻² and the heating as b⁻⁴. Two quantities delivered by the same distant encounter, against impact parameter in units of the target's half-mass radius, both logarithmic and both scaled to cross near one. The upper line is the impulse itself, 2Gm/vb, which every star in a bound system receives almost equally — so it moves the system and changes nothing inside it. The lower line is what is left after the common part is subtracted: the difference of the impulse across the system, which is its gradient multiplied by the system's own size, one power of b smaller and squared in the energy. The slopes are −2 and −4 exactly, and they are measured off the drawn curves rather than quoted. What follows is the point. Encounters at impact parameter b arrive at a rate proportional to b db, so summing the impulse over all of them gives ∫b⁻¹ db, which diverges logarithmically and is the origin of the Coulomb logarithm that appears in every treatment of relaxation. Summing the heating gives ∫b⁻³ db, which converges: extending the population from ten half-mass radii out to 300 multiplies the summed impulse by 2.48 and the summed heating by 1.010. Distant encounters diffuse velocities and heat nothing, and a tidal-heating calculation therefore needs no cutoff at large impact parameter, where a relaxation calculation cannot proceed without one.

The part of a kick that heats nothing

Most of the impulse a passing mass delivers to a star cluster is delivered equally to every star in it, so the cluster moves and nothing inside it changes. What heats it is the difference across it — one power of the impact parameter smaller, squared in the energy — and that one distinction decides which encounters matter and which cannot.

gravitation · Impulse approximation Supply and demand in a cluster's core. The rate at which hard binaries heat a star cluster's core, against the core's central density, for binary fractions of 1%, 5%, 20%, both on logarithmic axes and in units that put the balance for five per cent at one. The heating rate per binary does not depend on how hard the binary is — the energy released per encounter is proportional to its binding, and the chance of an encounter to its size, so the two cancel — and summed over a core whose radius shrinks as its density rises, the total grows as the square root of the density. The horizontal line is the demand: the energy the rest of the cluster carries away by relaxation through its half-mass radius, which the core has no say in. Each supply curve meets it at one density: 25 for 1%, 1.0 for 5%, 0.063 for 20%. The core contracts or expands until it produces exactly what is demanded, and a core rich in binaries does so at a lower density.

A core that makes what the cluster asks for

After core collapse, a star cluster's binaries heat it — but how much they heat it is not theirs to decide. The rest of the cluster carries energy outward at a rate fixed by its own size and mass, and the core contracts or swells until its binaries supply exactly that. So the core's structure follows from the halo's demand, and the cluster expands at a rate that forgets how it was born.

gravitation · Binary heating The same energy given to a cluster two ways. The distribution of the specific energies of 20,000 stars in a Plummer cluster (units G = M = a = 1), before and after an impulsive compressive shock of strength k = 0.06, which changes each star's vertical velocity by k times its height. Stars whose energy ends above the potential at the tidal radius, −0.196, marked by the vertical line, are lost. The thin curve is the cluster before. The dashed curve adds only the half-square term, the heating the textbook integral computes: every star gains a little, and 0.8 per cent of the stars cross the line. The solid curve adds the whole change, including the cross term between each star's own velocity and its kick, which averages to zero — the mean energy input is 1.59e-3 against 1.52e-3 — but broadens the distribution: 1.6 per cent of the stars cross. The same energy, 2.2 times as many escapers.

Stars lost to a kick that adds no energy

A tidal shock changes each star's energy by two terms. One is the heating every textbook integral computes, always positive and second order in the shock. The other is the product of the star's own velocity and its kick, first order in the shock and exactly zero on average — and for the weak shocks a cluster meets most often it decides who leaves, doubling the losses without adding any energy at all.

gravitation · Impulse approximation A slow shock still heats the cluster. The energy a Plummer star cluster gains from a compressive tidal pulse, as a fraction of what the same impulse delivered instantaneously would give, against the pulse's duration times the orbital frequency at the half-mass radius. The points come from integrating 1200 stars through the pulse one by one. The lines are prescriptions applied to the same stars: an oscillator's exponential suppression and a power law, each applied star by star at the star's own orbital frequency, and the shortcut of one power-law factor for the whole cluster, evaluated at the half-mass radius. At the longest pulse drawn the cluster keeps 0.28 of the impulsive energy; the star-by-star rules give 0.24 and 0.36, bracketing it, and the half-mass shortcut gives 0.017, 17 times too little. The total is set by the outermost stars, which carry most of the impulsive energy and are always slower than the shock.

A slow shock heats the halo and robs the core

The textbook treatment of a slow tidal shock multiplies the instantaneous answer by a suppression factor that falls steeply as the shock lengthens. Integrated star by star through shocks of intermediate length, a cluster behaves differently. Its total heating falls slowly, because the outer stars are always slower than any shock and carry most of the energy, and its core does not merely escape the heating: on balance its stars lose energy, which no factor that multiplies a positive number can describe.

gravitation · Impulse approximation Brighter by mass, fainter by temperature. Main sequences from 0.5 to 1.3 solar masses in luminosity against effective temperature, hotter to the left, for solar composition, [Fe/H] = −1.5, from homology applied to the solar-composition relation. The dots join stars of the same mass: lowering the metals moves each one up and to the left, so a 0.8 solar-mass star goes from 0.41 to 0.66 solar luminosities and from 5049 to 5879 K. The main sequence is steep, so a move of that size to the left carries a star past solar stars much brighter than itself: at a fixed temperature of 5800 K the metal-poor sequences lie 0.59 magnitudes below the solar one. A metal-poor star is brighter than a solar one of the same mass and fainter than a solar one of the same temperature — which is why such stars were called subdwarfs.

Brighter for its mass and fainter for its colour

Two stars of the same mass are almost the same star only if they are made of the same material. Take most of the metals out and the opacity falls, the radiation leaks out faster and the star settles brighter and much hotter — so far hotter that, at a given colour, it sits below the ordinary main sequence. The relation between mass and light has composition inside it, and every mass, age and distance read off it inherits that.

stars · The mass–luminosity relation The binaries in a fed core are spread evenly in energy. The steady distribution of binding energies among the hard binaries of a core where new binaries enter at x = 5 and leave by ejection from a cluster whose escape speed is 4 times its dispersion: the number per unit binding energy, against binding energy in units of mσ², from 20,000 simulated binaries weighted by the time each spends at each energy. Between twice the entry energy and a third of the typical ejection energy (about 120) it is flat to within 5 per cent, and it falls away smoothly near ejection because a single encounter can throw a binary out from anywhere above about 24. A flat distribution is what a population with a constant flux through energy and a constant drift must have: binaries pass every energy at the same rate, and spend the same time at each.

The binaries a core keeps are spread evenly in energy

A hard binary in a star cluster gains binding energy at the same average rate however hard it already is. So in a core that is fed new binaries and loses hardened ones, every binding energy is passed through at the same speed, and the steady population is flat in energy — which, counted per factor of energy, means most binaries are within a factor of ten of being thrown out. Where they were born hardly matters; where they leave decides everything.

gravitation · Binary heating Ejected before touching, or touching before ejected. The separation at which a hardening pair of 0.7-solar-mass main-sequence stars is typically ejected — where one encounter's recoil exceeds an escape speed 4 times the dispersion, at binding energy about 120 mσ² — against the cluster's velocity dispersion, both axes logarithmic, beside the separation at which the two stars fill their Roche lobes and begin to merge, 0.0082 au (1.8 solar radii). The ejection separation falls as the inverse square of the dispersion, because a hotter cluster needs a harder binary to throw it out. The lines cross at a dispersion of about 18 km/s. In clusters slower than that — open clusters and nearly all globulars — a binary is ejected while its stars are still apart. In faster ones — the most massive globulars and the nuclear star clusters at galaxies' centres — it is driven into contact first, and its two stars become one.

Ejected from a globular, merged in a nucleus

A hard binary in a star cluster hardens until something ends it: an encounter that throws it out, or its own stars touching. Which comes first is decided by one comparison — the separation at which the cluster's escape speed is exceeded against the separation at which the stars meet — and the first shrinks as the square of the cluster's velocity dispersion. In a globular cluster binaries leave while their stars are still apart; in the nucleus of a galaxy they merge first, and the heating that holds the cluster up becomes a factory of stars more massive than they should be.

gravitation · Binary heating A cluster sinking by Chandrasekhar's formula alone. The orbital radius of a 2.0 × 10⁵ solar-mass globular cluster released on a circular orbit at 1 kpc in the two Fornax halo models, decaying under Chandrasekhar's dynamical friction with a velocity dispersion of 11.7 km/s and a Coulomb logarithm of 3. In the cusp it reaches the centre in 6.7 billion years. In the core the formula gives no arrival within 12 billion years, crossing into the core at 2.8. The formula already slows the cluster inside the core, because the circular speed there drops below the stars' velocity dispersion and few stars move slower than the cluster; but it still lets the cluster creep inward without limit. N-body experiments of clusters sinking into cored haloes find something stronger: the decay stops altogether near the core, which the local formula cannot produce at any setting of its logarithm.

A drag that stops where every orbit keeps time

Chandrasekhar's formula reads the density and the velocities where a sinking body is, and nothing about the orbits those stars are on. Inside a constant-density core every orbit has the same frequency, the differential motion that drags a body inward disappears, and the sinking stops — which is why the five globular clusters of the Fornax dwarf galaxy are still where they are.

gravitation · Dynamical friction Which populations can share a cluster's temperature. The share of a cluster's mass in a heavy population against the mass of its members relative to an average star of 0.5 solar masses, both on logarithmic scales, with Spitzer's limit: a heavy component can reach energy equipartition with the stars only below the line where its share of the mass times its members' mass ratio to the power three halves equals 0.16. Above it the heavy bodies sink, cannot slow down enough to match the stars' kinetic temperature, and decouple into a subsystem of their own. Stellar black holes, 10 solar masses and about 2% of the mass, sit at 1.8 — 11 times beyond the limit. Neutron stars sit at 0.05 and massive white dwarfs at 0.14, below it; binaries, at a tenth of the mass, sit at 0.28, just beyond it, which is part of why they gather toward the centre. Only the black holes are beyond it by an order of magnitude.

A population that cannot share the cluster's temperature

Encounters between stars push a cluster toward equal kinetic energies, so heavy bodies slow down and sink. Spitzer showed in 1969 that past a certain point they cannot slow down enough: the heavy population falls away from the rest and becomes a cluster of its own. Stellar black holes are past that point by a factor of ten, in every globular cluster that kept them — and the subsystem they form, not the stars' binaries, is what holds most clusters' cores open until the black holes have been spent.

gravitation · Binary heating The clusters scale with the halo and the stars do not. The total mass in a galaxy's globular clusters (lower, with its ±0.3 dex scatter shaded) and its stellar mass (upper), against the mass of its dark-matter halo, all logarithmic. The clusters follow a straight line of slope one, 3.9×10⁻⁵ of the halo, across five and a half decades. The stars follow a curve that bends at both ends round a peak near 5.8×10¹¹ M☉, because star formation is suppressed in small halos and in large ones. Stars outweigh clusters by 12 times in a halo of 10¹⁰ M☉, 879 times at 10¹², 31 times at 3.2×10¹⁴.

A halo weighed by counting its oldest clusters

The globular clusters round a galaxy add up to about four parts in a hundred thousand of the mass of its dark-matter halo — in a dwarf, in the Milky Way, in a giant elliptical and in a cluster of galaxies, across five decades of halo mass. The galaxy's own stars follow no such rule. So a count of a few dozen old star clusters weighs a halo that its stars cannot, and the strangest thing about the rule may be that it is less strange than it looks.

galaxies · Dark matter

Named alongside it

The objects these essays reach for when they reach for this one.

Binary heatingCore collapseCoulomb logarithmImpulse approximationRelaxation timeTwo-body relaxationDark matter haloDynamical frictionEscape velocityHard binaryMass segregationPlummer sphere

All concepts