Astronomy & Celestial Mechanics
the sky, computed
Series

Impulse approximation — the series

5 essays on one idea, from the one that introduces it to the one that assumes the rest.
  1. The whole kick, delivered in about two encounter times. The transverse force a body feels while a mass sweeps past it on a straight line, and the velocity that force has delivered so far, both against time in units of the impact parameter divided by the relative speed. The force is the component of the inverse-square attraction perpendicular to the path, which is the impact parameter over the cube of the distance, and it is drawn at its peak value of one at closest approach. The rising curve is its running integral, scaled by twice the gravitational constant times the mass over the impact parameter and the speed. Two things are visible and both are the point. The integral of the force over all time is exactly two in these units, so the kick is exactly 2GM/bv with no free constant anywhere — an answer to a three-body-shaped question obtained without solving anything. And it arrives quickly: 71 per cent of it within a single encounter time of closest approach and 98.6 per cent within the 6 drawn, which is what licenses calling the whole thing an impulse. On the scale of anything slower, the velocity changes discontinuously.

    An answer obtained along a path that was not taken

    Integrate the force of a passing mass along the straight line the body would have followed if the encounter had not happened, and out comes an exact deflection with no free constant in it. The approximation is circular, it is wrong in a known direction by a known amount, and it is the reason stellar dynamics has closed forms at all.

    part 1 · gravitation
  2. 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.

    part 2 · gravitation
  3. 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.

    part 3 · gravitation
  4. 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.

    part 4 · gravitation
  5. How long a wide pair survives a halo of heavy objects. The time for a binary of one solar mass to be dissolved by passing compact objects, against its separation, both logarithmic, if the whole local dark-matter halo — 0.01 M☉ per cubic parsec, moving past at 300 km/s — is made of objects of 10, 100, 1000 M☉. A pair 0.1 pc apart (about 20,000 AU) lasts 477, 48 and 4.8 billion years respectively; at 1 pc, 47.7, 4.8 and 0.48. The heavier the objects at fixed total density, the fewer and harder the kicks, and the sooner the widest pairs go. The dashed line is 10 billion years, the age of the halo's oldest binaries.

    Pairs of stars too loosely bound to survive heavy dark matter

    Two stars a tenth of a parsec apart are held together by a force so weak that a passing mass of a few hundred suns, at a distance of a parsec, can loosen them. If the dark matter round the Milky Way were made of heavy compact objects, such pairs would have been unbound long ago. Some have survived ten billion years in the halo. The widest of them sets a ceiling on how heavy the dark matter's pieces can be — and the ceiling has moved twice, both times because of one pair.

    part 5 · gravitation

All series