Astronomy & Celestial Mechanics
the sky, computed
Depth

Series — page 6

A field says what an essay is about. A series follows one idea essay by essay — from the question that introduces it to the one that assumes all the others.
Five rules, and how fast each one leaves the seasons behind. The accumulated difference between a calendar and the seasons, for five leap-year rules, over 4000 years. Each is a straight line whose slope is the rule's fraction minus the tropical year's 0.2421897, and nothing else about a calendar matters to this plot — not which months are long, not where the year starts, not what anything is called. The Julian quarter is out by 11.25 minutes a year, which is one day every 128 years and is why ten days had to be removed in 1582. The Gregorian rule takes 3223 years to lose a day and the Persian 4264, in the other direction — the older rule is the better one, at a denominator of 33 against 400. The 128-year rule is level on this scale: 454545 years to a day, which is longer than any calendar has been kept — and it is not a separate invention but the Julian rule with its own error subtracted, since 1/4 − 1/128 = 31/128 exactly and the Julian drift is exactly one day per 128 years. What the figure cannot show is that the tropical year is itself shortening, by about half a second a century, so the slopes drawn here are the ones at the present epoch and every line is very slightly curved.

Calendars

  1. 1 A year that is not a whole number of days
  2. 2 A month that has to be tabulated
2 essays · sky A separation of 10⁻¹⁰ becomes 0.0285. The distance between two copies of Burrau's problem, started 10⁻¹⁰ apart in one coordinate of one body and then stepped in lockstep — one loop, one step size, both states advanced by it — for 60 time units. The vertical axis is logarithmic, so the straight stretch is exponential growth, and its slope is the Lyapunov exponent: 0.3413 per time unit, fitted over the 604 samples lying between ten times the starting offset and a tenth of the system's own size. That is an e-folding every 2.93 time units, so the separation multiplies by ten every 6.75 and by 7.8·10⁸ across the whole run. It has not saturated within the drawn interval: the growth rate over the last tenth of the run is still 0.6002 per time unit, 176% of the fitted exponent, and the separation has reached only 0.90% of the system's own size. The exponent is what sets a prediction horizon, and this run is drawn short of it on purpose: the straight line is the measurement, and the flat part that follows is arithmetic about how far apart two bounded systems can get.

Chaos

  1. 1 A prediction with an expiry date
  2. 2 Where the chaos comes from
2 essays · gravitation 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.

Collisional cascade

  1. 1 Where the mass is and where the light is
  2. 2 The drag that sorts a disc by size
2 essays · orbits Every orbit one force allows. Circle, ellipse, parabola and hyperbola, all sharing a focus and a closest approach. The eccentricity alone decides which one a body is on, and whether it returns.

Conic sections

  1. 1 Every orbit one force allows, and the number that picks between them
  2. 2 The third thing that is conserved
2 essays · orbits The mass-to-light ratio of NGC 3198, radius by radius. The dynamical mass inside each radius divided by the light inside it, in solar units. Inside two disc scale lengths it is nearly flat at about 2.2 — a galaxy made of stars, weighing what stars weigh. Outside the disc the light stops and the mass does not, so the ratio climbs to 10 by 30 kpc with no sign of turning over. The curve is not fitted: it is the rotation curve's enclosed mass divided by the light profile's enclosed light, both drawn elsewhere in this collection.

Dark matter

  1. 1 The mass that is not the light
  2. 2 Three mass models that fit the same curve
2 essays · galaxies Contrast against separation, which is where direct imaging lives. Planet-to-star brightness ratio against apparent separation, both logarithmic, for a system 10 parsecs away. The reflected-light curves are A_g(R_p/a)² and fall as the inverse square of the orbit; the thermal curve is the ratio of two Planck functions at 10 µm and does not, which is why every imaged planet so far is young and hot rather than merely large. The vertical lines are diffraction limits λ/D — nothing inside a telescope's own line is reachable by it at any contrast at all. An Earth at ten parsecs sits at 6.3e-10, which is 4 orders of magnitude below the faintest planet yet imaged. The four imaged planets are plotted at their measured near-infrared contrasts rather than at 10 µm, because the near infrared is the band they were found in — which is itself part of the argument, since a young planet is hot enough to be bright where its star is not.

Direct imaging

  1. 1 Nine orders of magnitude, half an arcsecond apart
  2. 2 A star subtracted using the star
2 essays · exoplanets Which gradient is steeper, and where. The temperature gradient radiation would need in order to carry the whole luminosity, against the gradient a rising blob of gas follows, through an n = 3 polytrope with a Kramers opacity and an energy generation going as ρT^4. Wherever the first exceeds the second the layer is unstable and convects. This model puts the outer boundary at 89 per cent of the radius; helioseismology measures the base of the Sun's convection zone at 71.3 per cent, and the gap between the two numbers is the price of a polytrope with one opacity law rather than a solar model. The curve also rises above the adiabatic value inside 4 per cent of the radius, which is a convective core — it appears when the energy generation is concentrated enough, and it is the structure a star burning by the CNO cycle has.

Energy transport

  1. 1 The part of a star that boils
  2. 2 A length nobody derived, fitted to one star
2 essays · stars Mars to five metres and Neptune to five thousand kilometres, in the same file. Present-day heliocentric position uncertainty for each planet, in kilometres, with the range component marked separately below it. The two differ because a transponder measures a distance along the line of sight and says nothing about the two directions across it, so a planet with an orbiter is known radially some 17 times better than it is known altogether. Neptune is 10⁶ times less well determined than Mars and only 20 times further away, which is the whole point: the accuracy is a property of the observations, not of the geometry. Mars has carried a transponder almost continuously since 1976; Neptune has been visited once, in 1989, and everything else known about it is meridian-circle astrometry covering 1.07 of one orbit. The two ice giants are the only entries here whose ephemerides are still limited by nineteenth-century technology, and the only cure is a spacecraft.

Ephemerides

  1. 1 The table that is a fit
  2. 2 Where a planet is and where it is seen
2 essays · orbits Energy generation against core temperature. The proton–proton chain and the CNO cycle, in solar units, against core temperature on logarithmic axes. The CNO curve is far steeper, so the two cross at 18.8 million kelvin — above that temperature a star runs mostly on CNO, and below it mostly on pp.

Fusion

  1. 1 The furnace that runs cooler than a compost heap
  2. 2 The resonance that had to exist
2 essays · stars S2's orbit, and the mass it implies. The orbit of S2 about the centre of the Milky Way, drawn from its measured elements: a semi-major axis of 0.1251 arcseconds, which at 8.28 kpc is 1036 AU, an eccentricity of 0.8843, and a period of 16.05 years watched round twice. Kepler's third law in solar units — a³/P² — gives 4.31 million solar masses, by exactly the calculation that weighs a planetary system. At periapsis the star is 120 AU from the focus, 1407 Schwarzschild radii, and the mean density inside that radius is 5.3e+15 M☉ per cubic parsec — a million times the densest star cluster known. That density, not the mass, is the argument: there is no configuration of stars that would fit.

Galactic nuclei

  1. 1 The mass at the centre that is not stars
  2. 2 The last parsec, and the stars that are not there
2 essays · galaxies A bridge and a tail, integrated. A disc of 180 massless particles on circular orbits, and a companion of 1 times the primary's mass on a parabolic orbit with a pericentre of 1.6 disc radii, integrated from before the encounter to well after it. Times are in units of the disc's own outer orbital period, measured from pericentre. The outermost ring, drawn separately, is the one that produces both the bridge and the tail: at 3 its furthest particle is 31.4 disc radii from the centre, having started at 0.9; the panels are scaled to hold ninety per cent of the particles, so the very end of the tail is outside them. Nothing has been ejected and no material is new — every particle is on the orbit its own initial conditions and the two masses give it.

Galaxy interactions

  1. 1 A bridge and a tail drawn by one force
  2. 2 Whether it merges is a ratio of two times
2 essays · galaxies Period against size for the planets, around the Sun. Orbital period against semi-major axis on logarithmic axes. The line has slope exactly three-halves — the harmonic law — and the measured bodies sit on it.

Harmonic law

  1. 1 The law that links period to size, and weighs everything
  2. 2 The third law is wrong by the mass of the planet
2 essays · orbits The Hill radius against the mass ratio, and where four planets keep their outermost moon. The Hill radius R_H = a(μ/3)^⅓ and the exact L₁ distance from the same quintic, both in units of the planet's own semi-major axis, against the mass ratio μ = m/(M+m) from 10⁻⁸ to 0.12. The two are the same line to the width of the stroke across the whole planetary range — the cube-root formula runs 0.33% short of the exact L₁ distance at the Earth's mass ratio and 10.0% short at the right-hand edge — so the departure between them is drawn on its own percentage scale at the right, which is the only place it can be seen. Below the curves are the numerically determined stability limits, 0.5 R_H for a prograde satellite and 0.7 R_H for a retrograde one, and the outermost known satellite of each of four planets, each at a Hill fraction computed from that planet's own mass and orbit: Moon 0.257, Sinope 0.451, Phoebe 0.198, Neso 0.424. Of the four, three are retrograde — Sinope, Phoebe, Neso — and they reach 0.451 of their planet's Hill radius against 0.257 for the one prograde satellite: 1.8 times as far out, in a diagram whose two stability limits stand in the ratio 0.7 to 0.5. No equilibrium argument predicts that, and it is what the drawing makes plain.

Hill sphere

  1. 1 The region a planet may keep a moon in
  2. 2 A feeding zone, and the spacing it forces
2 essays · gravitation Inside a star that is holding itself up (polytrope n = 3). Temperature, density and pressure through a star, as fractions of their central values, against fractional radius. All three come from one numerical integration of the Lane–Emden equation at index 3, which is the hydrostatic balance written for a gas whose pressure is a power of its density. The inner half of the radius holds 90% of the mass, and the outer half is nearly weightless — which is why the load, and so the temperature, is concentrated where the burning is.

Hydrostatic equilibrium

  1. 1 A star is held up by its own weight
  2. 2 A floor under the centre that assumes nothing
2 essays · stars 1I/ʻOumuamua: an orbit at e = 1.201, and the angle it turned through. The open branch of a conic at eccentricity 1.201 and periapsis 0.2559 AU, with the Sun at the occupied focus. Two numbers fix everything else: |a| = q/(e−1) = 1.273 AU, an impact parameter b = |a|√(e²−1) = 0.847 AU, an asymptote at ν∞ = arccos(−1/e) = 146.37° from periapsis, and a deflection of δ = 2ν∞ − 180° = 112.74° — which is the same statement as sin(δ/2) = 1/e. The asymptotes cross at |a|e = 1.529 AU from the Sun on the periapsis side, where a bound orbit's centre would be on the other. The speed left over at infinity is √(μ/|a|) = 26.40 km/s. Schematic in one respect: the drawing runs at about 119 px to the AU, so the Sun's own disc would be far smaller than its marker.

Hyperbolic orbits

  1. 1 The orbit that has no period
  2. 2 Aiming at a plane instead of at a planet
2 essays · orbits 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.

Impulse approximation

  1. 1 An answer obtained along a path that was not taken
  2. 2 The part of a kick that heats nothing
2 essays · gravitation One population, three integrals, two ends. The same mass function weighted three ways, each normalised so that the area under it is one, against mass on a logarithmic axis — so a share of the page is a share of the total. The light is a zero-age population's, every star still on the main sequence, which is the only age at which the top of the range is present at all. The steep curve on the left is the number of stars, which peaks at 0.08 M☉ and falls away because the function is steeper than m⁻¹; the middle one is the mass they carry; the one on the right is the light they emit, computed from this file's own main-sequence relation and peaking at 55.0 M☉ — 682 times further up the axis. The population that is counted and the population that is seen are two different populations. Above 8 M☉ there are 0.63% of the stars, 21% of the mass and 99% of the light. The slopes are α = 1.3 above 0.08 M☉, α = 2.3 above 0.5 M☉, and the low-mass end is where the honesty runs out: it has never been measured in another galaxy, and every mass inferred from a luminosity assumes it.

Initial mass function

  1. 1 Most stars are small, and most of the light is not
  2. 2 A mass function corrected by an age
2 essays · stars The transform plane an array of 9 actually samples. Every point in this plane is a spatial frequency the array has measured, in thousands of wavelengths. The 9 antennas make 36 pairs, each pair measures one point at any instant, and turning the Earth sweeps each of them along an ellipse — so eight hours of tracking turns 36 measurements into the arcs drawn here. Two properties are structural rather than chosen. The plane is Hermitian: a real sky forces V(−u,−v) = V(u,v), so half the points are free and the coverage is symmetric through the origin. And every ellipse has axis ratio exactly sin δ = 0.707* at this declination, measured off the longest track as 0.707 — an array is squashed in one direction by where the source is in the sky, and at the equator the tracks collapse to lines whatever the array. What the figure cannot show is the hole in the middle: no baseline is shorter than an antenna is wide, so the largest structures on the sky are simply not measured, and no processing recovers them.

Interferometry

  1. 1 Resolution without a mirror
  2. 2 A phase that survives what corrupts it
2 essays · starlight The Sun's internal rotation: differential above 0.693 R, rigid below, and the 15 nHz between two modes. Angular velocity against fractional radius, in nanohertz, from the inversion of hundreds of thousands of measured rotational splittings. Two features had to be discovered rather than deduced. The convection zone rotates differentially in latitude — 473 nanohertz at the equator, 330 at sixty degrees — and that latitude dependence persists all the way down through it, in surfaces that are very nearly radial rather than the cylinders a rotating convective fluid was expected to produce. And below about 0.693 of the radius the latitude dependence stops: the radiative interior turns as a single rigid body at 430 nanohertz, which is a remarkable thing for a fluid with no strength to do, and requires something — most likely a weak internal magnetic field — to be enforcing it. The two regimes are joined by a shear layer a few per cent of the radius thick, the tachocline, and it is where the solar magnetic cycle is generally thought to be generated, because it is the only place in the Sun with the shear a dynamo of that strength needs. The horizontal marks are what individual modes would report: the rotation averaged over the cavity each one occupies, weighted by the time the wave spends at each radius. A mode of degree 100 is trapped near the surface and reports 473 nanohertz; one of degree 1 passes through the deep interior and reports 458. The 15-nanohertz difference between them is the measurement, and splittings are measured to about a nanohertz — which is why the profile can be resolved at all. The kernels used here have the right support and the right sense; a real inversion uses computed eigenfunctions, and its resolution below 0.2 R is poor for the reason the previous figure gives.

Internal rotation

  1. 1 A fluid that turns as one piece
  2. 2 A shear layer that should have spread
2 essays · stars Two densities in thermal balance at one pressure, 121 times apart. Thermal equilibrium for interstellar gas, drawn as pressure against density with both axes logarithmic. Every point on the curve is a temperature between 40 and 9000 K at which cooling exactly balances the 2·10⁻²⁶ erg per second per hydrogen nucleus that grain photoelectrons deliver. The curve is not monotonic: it rises to 5007 K cm⁻³, falls to 1597, and rises again, so a horizontal line anywhere between those two crosses it three times. At the 3000 K cm⁻³ of the local medium the three crossings are a warm phase at 0.47 cm⁻³ and 6354 K, a cold phase at 57 cm⁻³ and 52 K, and one in between drawn dashed because it cannot survive — a parcel there that is squeezed cools and keeps contracting, and one that expands heats and keeps expanding. The two survivors differ by a factor of 121 in density and by exactly the same factor in temperature — necessarily the same, because their product is the pressure both are held at — and yet they press on one another equally, which is why they can share the same volume of the disc indefinitely rather than mixing.

Interstellar medium

  1. 1 Two temperatures, and nothing in between
  2. 2 A disc held open by what cannot be photographed
2 essays · galaxies

All essays