Astronomy & Celestial Mechanics
the sky, computed
The collection

Every essay — page 31

One idea per essay, ordered so that the earlier ones set up the later ones — but nothing here depends on being read in sequence. Essays 601–620 of 909.

Stars

The diagram that sorted them, and the one quantity that decides a star's whole life.

Field strength against radius, from the Sun to a magnetar. Surface magnetic field against radius, both on logarithmic scales, radius in solar radii and field in gauss, for five kinds of star: the Sun's dipole, Ap and Bp stars, magnetic white dwarfs, radio pulsars, magnetars. The diagonal lines are constant magnetic flux, πR²B = 10²³, 10²⁵, 10²⁷ gauss square centimetres, along which a star shrinking with its field frozen in would slide. The Ap and Bp stars, the magnetic white dwarfs and the magnetars all fall close to the upper lines: a field of a few kilogauss on a star of three solar radii, squeezed to the size of a white dwarf, becomes a hundred million gauss, and squeezed to twelve kilometres becomes a hundred trillion. The Sun's dipole and the ordinary pulsars carry a hundred to a thousand times less flux.

The same flux in three sizes

A few per cent of hot main-sequence stars carry kilogauss fields that do not change over decades. So do a similar share of white dwarfs, at a million times the strength, and of young neutron stars, at a hundred billion times it. Multiply each field by its star's area and the three classes carry nearly the same magnetic flux — what a field frozen into a star would keep as the star shrank by a factor of a hundred thousand. A collapse makes a field enormous, but it cannot make it important.

6 figures · Flux freezing A white dwarf cooling track in J−K and absolute V. The track a 0.6 solar-mass hydrogen-atmosphere white dwarf follows as it cools, absolute V magnitude against J−K colour, fainter downward. Dashed, a blackbody of the same radius, which slides steadily redward and fainter. Solid, with collision-induced absorption: the track reddens to J−K = 0.40 at 4250 K and absolute magnitude 16.1, then hooks back to the blue, reaching −0.86 at 3000 K and magnitude 18.1. Temperatures are marked along it. The oldest white dwarfs therefore sit in the same colour range as much younger, hotter ones and are told apart only by being fainter — which is why they were missed by surveys that selected by colour.

A colour that runs backwards at the end of cooling

A white dwarf reddens as it cools, like every other hot body — until its surface falls below about four thousand kelvin. Then pairs of hydrogen molecules, colliding in an atmosphere dense enough to crush them together, begin to absorb its infrared light, and the star turns blue again while it goes on getting fainter. The oldest white dwarfs in the Galaxy have the colours of young ones, and for years they were looked for in the wrong place.

6 figures · White dwarf cooling

Starlight

Distance, brightness and colour: what can be measured when nothing can be visited.

The extinction curve of a mixture, and of single sizes. Extinction relative to the V band against inverse wavelength, from the infrared at 0.3 inverse microns to the near ultraviolet at 3.3. The thick grey curve is the observed Galactic average, the Cardelli–Clayton–Mathis law at Rᵥ = 3.1. The coloured curves are computed: MRN, 0.005–0.25 µm, Rᵥ = 2.4; all grains 0.05 µm, Rᵥ = 1.8; all grains 0.3 µm, Rᵥ = 18.4. Grains of any single size make a curve with a knee: steep where the wavelength is longer than the grain, flat where it is shorter. Only a spread of sizes, with many more small grains than large ones, makes the nearly straight rise with inverse wavelength that is observed, and the first model follows the observed curve to an r.m.s. of 0.19 in these units. The ultraviolet bump beyond this range needs a separate population of small carbonaceous grains and is not modelled.

A size read from the slope of a colour

Interstellar dust dims blue light more than red, and the way the dimming rises toward the blue is nearly a straight line in inverse wavelength. No grain of any single size makes that line. A grain only removes light selectively when it is comparable in size to the wavelength, so the extinction curve read across wavelength is the population of grains read across size — and it says there are vastly more small grains than large ones, in the same proportion that grinding rocks together produces.

6 figures · Extinction Hayashi lines: a luminosity that moves and a temperature that does not. The effective temperature a fully convective star must have, against its luminosity, for 0.5, 1, 2 solar masses, drawn as a Hertzsprung–Russell diagram with temperature falling to the right. The temperature follows from matching the star's convective adiabat to a photosphere whose opacity, from the negative hydrogen ion, rises as the ninth power of temperature. Over a factor of thirty thousand in luminosity each line moves by only 363 K, 398 K, 436 K respectively: the lines are nearly vertical, at 3650 K, 4000 K, 4383 K for one solar luminosity. A contracting young star, or a swelling giant, moves up or down its line almost without changing colour. To the right of the lines no star in hydrostatic equilibrium can sit.

A wall built by a steep opacity

A star that is convective all the way to its surface cannot choose its surface temperature. The interior fixes a relation between pressure and temperature, the photosphere fixes another through its opacity, and the negative hydrogen ion's opacity rises so steeply with temperature that the two meet at nearly the same temperature whatever the star's size. Young stars contracting and old stars swelling slide up and down a line in the diagram almost without changing colour — and nothing in equilibrium can sit to its right.

6 figures · Opacity

Cosmology

One object, seen once, from inside — and every number in it the output of a model.

The Hubble constant the microwave background implies, with an early component. The Hubble constant that keeps the microwave background's measured angular sound horizon fixed, against the peak fraction of early dark energy at z = 3500, with the baryon density held at its standard value and the matter fraction at the 0.315 that late-time distances measure, so that the dark-matter density rises with the Hubble constant. With none, it is 67.4 km/s/Mpc and the horizon 144.3 Mpc. At 5 per cent it is 68.7; at 10 per cent 70.2, with the horizon at 138.5 Mpc; the local value of 73 is reached at 19 per cent. The shrinkage of the horizon is what raises the constant: a shorter ruler at the same angle must be nearer, and nearer at a fixed matter fraction means faster expansion. The extra dark matter this requires is the cost the published fits also find: at the value that reaches 73, the physical matter density is 18 per cent higher than in the standard model.

A ruler shortened by something already gone

The microwave background gives the Hubble constant only through a computed length, the distance sound travelled before the universe became transparent. Add a component of energy that exists for a few thousand years before recombination and then vanishes faster than radiation, and the universe expands a little faster while the sound is travelling, the ruler comes out a few per cent shorter, and the inferred expansion rate rises toward the local value. The price is paid in dark matter and in the finest details of the oldest light.

6 figures · Hubble constant How far a void field turns the leptons that make each cascade photon. The deflection angle of a cascade lepton over its cooling length, in radians on a logarithmic scale, against the energy of the cascade photon it produces, for fields of 10⁻¹⁸ G, 10⁻¹⁷ G, 10⁻¹⁶ G, 10⁻¹⁵ G, 10⁻¹⁴ G coherent over at least that length. The angle falls as the square of the lepton energy, so as the inverse of the photon energy: the gigaelectronvolt halo is made by the slowest leptons, which are turned most. At 1 GeV the deflections are 1e-3, 0.01, 0.11, 1.06, 10.62 rad. The dashed line is the jet's half-opening angle, 0.1 rad; a lepton turned further than that sends its photons away from the observer, and the halo at that energy disappears.

A lower limit made of three lengths and a missing echo

The only measurement of a magnetic field in the empty space between galaxies is a lower limit, and it comes from light that did not arrive. Teraelectronvolt gamma rays from distant blazars make electron pairs in intergalactic space, the pairs should re-emit the energy at gigaelectronvolts in a halo around each source, and the halos are missing. Whether that absence weighs a field of a ten-thousandth of a billionth of a microgauss depends on three lengths, the field's own tangle, and how long the blazar has been shining.

6 figures · Primordial fields

Exoplanets

Planets nobody has seen, weighed and measured from a dip, a wobble and a delay.

Brightness temperature across the infrared, for three profiles. The brightness temperature of the planet's emission against wavelength from 1 to 12 microns, for the three temperature profiles, with water and carbon monoxide bands in the opacity. In a band the opacity is high, so the light comes from higher, lower-pressure layers; in the windows between bands it comes from deeper. Over a profile falling outward the bands are cooler than the windows: at 6.3 microns 1457 K against 1570 K at 3.6. Over the inverted profile they are hotter: 1851 against 1538. Over the isothermal profile every wavelength has the same brightness temperature and the spectrum is a blackbody, whatever molecules are present.

Molecules seen only through a change of temperature

The heat a planet gives off, measured colour by colour as it disappears behind its star, carries the fingerprints of the molecules in its atmosphere — but only because the atmosphere's temperature changes with height. Water absorbs where the air cools upward and emits where it warms, and in an atmosphere at one temperature every molecule vanishes from the spectrum, however much of it there is. An emission spectrum is a sounding of a temperature gradient, and the molecules are what it is taken with.

6 figures · Exoplanet atmospheres What shares a resolution element with an Earth twin. Light collected in one resolution element around an Earth-like planet in the habitable zone of a Sun-like star, against the system's distance, both on logarithmic scales, for a 6 metre telescope at visible wavelengths, in units of the planet's own light at 10 parsecs. The planet and the leaked starlight both dim as the square of distance. The exozodiacal dust, at 3 times the Solar System's level, and the Solar System's own zodiacal light do not dim at all: each is a surface brightness, and a resolution element at a greater distance simply contains more dust. The exozodi equals the planet at 9.1 parsecs and is 11 times it at 30.

Dust that does not dim with distance

An Earth-like planet around a nearby star is seen against the glow of the dust in its own system, and the two behave differently with distance. The planet is a point, and it fades as the square of how far away it is. The dust is a surface, and a surface does not fade at all — a more distant system packs more dust into each resolution element in exactly the proportion its light is diluted. So the dust around a star, not the telescope's contrast, sets how far away an Earth can be found, and how many stars a search can reach.

6 figures · Direct imaging

Orbits

Kepler's three laws, and the family of curves a single force allows.

Debris discs of very different masses, fading to one line. The fractional luminosity of a debris belt at 50 AU — the share of its star's light the dust absorbs and re-radiates — against age, from a million to ten billion years, for belts that start with 1, 10, 100, 1000 Earth masses of planetesimals up to 60 km across. Each stays near its starting brightness until its largest bodies have had time to collide, 16.5 billion years, 1.7 billion years, 165 million years, 17 million years respectively, and then fades as one over the age. A heavier belt starts brighter and grinds itself down faster, and the two effects cancel: after their collisional times all belts lie on the same line, the dashed maximum, whatever their mass. At ten billion years the two heaviest differ by 1 per cent.

A ceiling that forgets how much there was

A belt of planetesimals around a star shows itself only through its dust, and the dust is ground from the belt's largest bodies by collisions. Give the belt more mass and it makes more dust — but it also collides more often, and grinds itself away faster. At late times the two effects cancel exactly, and every belt at a given distance and age lies on the same line of maximum brightness whatever it started with. A disc above that line cannot be a steady cascade, and it is telling of a recent catastrophe.

6 figures · Collisional cascade An inner orbit's inclination, with and without the octupole. The inclination of an inner orbit to the plane of a distant perturber against time, in units of the quadrupole timescale, starting at 65 degrees with a nearly circular orbit, for octupole strengths ε = 0 and 0.03. At ε = 0 — a circular outer orbit — the inclination oscillates in the regular Kozai–Lidov cycle and never crosses ninety degrees, because the angular momentum along the outer orbit's axis is conserved and cannot change sign. With the outer orbit eccentric, the inclination drifts from cycle to cycle and crosses ninety degrees, the orbit turning from prograde to retrograde, first at t ≈ 120, 2 times in all within the time shown. The integration conserves the secular energy to 7·10⁻⁶.

An orbit that turns over

A distant third body makes an inner orbit trade inclination for eccentricity in slow, regular cycles, and one quantity — the inner orbit's angular momentum along the outer orbit's axis — keeps the trade bounded and forbids the orbit from ever tipping past ninety degrees. That conservation holds only if the outer orbit is a circle. Make it eccentric and the quantity drifts, the cycles stop repeating, and at a moment of extreme eccentricity the inner orbit turns over, from prograde to retrograde, with its pericentre pressed a thousand times closer than the regular cycle allowed.

6 figures · Kozai–Lidov

The observed sky

Coordinates, seasons, phases and shadows — geometry seen from inside it.

Mars's tilt over fifty million years, and a steadier twin. The obliquity of a Mars-like planet against time over 50 million years, integrated from 25 degrees with its orbital plane wobbling in five secular modes. With Mars's precession constant, 8.26 arcseconds a year, the axis wanders between 1 and 70 degrees, lingering near some values and lurching between them. The same planet precessing 6 times faster — about the Earth's rate, which the Moon's pull makes possible — stays between 23.2 and 30.8 degrees, because its precession is too fast to resonate with any of the orbital-plane frequencies. The orbital perturbations are illustrative, of the size Mars's own are; the chaos is computed from them.

A past that can only be given as odds

Mars's axis is tilted by twenty-five degrees today, and nothing holds it there. Its slow precession falls among the frequencies at which its orbital plane wobbles, the resonances overlap, and the tilt wanders chaotically between nearly zero and sixty degrees. Two calculations that start a hair apart part company within a few tens of millions of years — so Mars's obliquity a billion years ago is not a number anybody will compute. What can be computed is how much of its history the planet spent at each tilt, and that distribution is a prediction the ice at its poles can test.

6 figures · Obliquity The nutation lines and the core's resonance. The resonance factor |σ₀ / (σ − σ₀)|, σ₀ being the free core nutation's frequency, on a logarithmic scale, against frequency in the celestial frame in cycles per sidereal day, negative for retrograde motion, with the free core nutation at a period of 430.2 days retrograde. The vertical lines are the largest forced nutations the Sun and Moon drive: the 18.6-year, annual, semiannual and fortnightly terms, each prograde and retrograde, and precession at zero. Any part of the Earth's response that involves the core is multiplied by this factor. The retrograde annual term, 365 days, is 5.6 on this scale; the prograde annual term 0.46; the retrograde 18.6-year term 1.07. The annual retrograde nutation is the one the core's resonance distorts most, because a year is the forced period nearest the core's own.

A core that wobbles at its own period

The Earth's liquid core can turn about an axis slightly tilted from the mantle's, and the flattened boundary between them pushes it back. That gives the planet a second free mode of rotation: in space a slow retrograde nutation of about 430 days, on the ground a wobble just faster than once a day. Its period is a measurement of the shape of a surface under 2,900 kilometres of rock, and it is read twice — in the nodding of the Earth's axis against the quasars, and in the rise and fall of the solid ground under the daily tides.

6 figures · Precession

Spaceflight

Celestial mechanics used forwards — where to burn, and what it costs.

How much future debris each removal prevents. The share of the derelict population's expected future fragment production removed, against the number of objects removed, out of 1858, for three orders of removal: by the ranking weight — collision rate times fragments times fragment lifetime — by mass alone, and at random. Removing the 200 top-ranked objects removes 70 per cent; the top fifty alone, 45 per cent; the top ten, 19. Removing by mass does nearly as well at first and falls behind, because mass is most of the weight and altitude the rest. Removing at random, fifty objects remove 3.5 per cent. The population is synthetic, built to the published inventory's broad shape; the concentration of risk in a few tens of heavy upper stages is the finding the real analyses also make.

Fifty objects that hold most of the future

Low Earth orbit holds thousands of dead satellites and spent rocket stages, and the debris of the future will come mostly from their collisions. But not in equal shares. A derelict's contribution is the chance it is hit, times the fragments it would make, times how long they would stay up — and all three favour the heaviest objects in the most crowded, highest orbits. Ranked that way, a few tens of rocket bodies near 850 kilometres carry nearly half the risk, and the arithmetic of cleaning space becomes the arithmetic of choosing them.

6 figures · Orbital debris Propellant a moon saves, against how close the flyby comes. The reduction in the Jupiter orbit-insertion burn bought by one flyby of a Galilean moon on the way in, in metres per second, against the flyby altitude, for an arrival at 5.6 km/s, a periapsis of five Jupiter radii and a capture orbit of 200 days, whose insertion alone costs 825 m/s. At 300 km: Io 268 m/s, Europa 122 m/s, Ganymede 191 m/s, Callisto 124 m/s. The saving falls with altitude because the turn angle does, roughly as one over the flyby distance once that is well above the moon's surface. Io saves most: it is deep in Jupiter's well, where it moves fastest, so a given turn of the relative velocity takes more energy from the spacecraft. Ganymede, the most massive, is second.

Braking on a moon before braking with the engine

A gravity assist is usually told as a way to gain speed, but the same exchange runs the other way. A spacecraft falling toward Jupiter can pass behind one of its moons on the way in, lose a few hundred metres per second of its speed to the moon's orbit, and arrive at its capture burn needing a fifth less propellant. Which moon brakes best depends on three things at once — its mass, how fast it is met and how fast it moves — and the best one sits in radiation no spacecraft wants to fly through.

6 figures · Gravity assist

Cosmology

One object, seen once, from inside — and every number in it the output of a model.

How much of a background a survey can resolve, against sources per beam. The fraction of a background resolved into individual sources before confusion stops the count, against X, the number of sources brighter than the break in the counts that fall inside one beam, on a logarithmic axis. The curve is computed for counts falling as S^−2.5 above the break and S^−1.5 below it, with a source counted as recoverable while there is less than one source per 30 beams brighter than it. It is a single curve: telescope, wavelength and integration time enter only through X. Half the background is resolved at X ≈ 0.033, and the fall from nearly all to nearly none takes about four decades. The surveys are placed on it at their own beams and counts: Hubble, optical at X = 8.7·10⁻⁵, 100 per cent; Chandra, soft X-ray at X = 6.4·10⁻⁶, 100 per cent; Spitzer, 24 µm at X = 0.012, 54 per cent; Herschel, 250 µm at X = 0.043, 24 per cent; Herschel, 500 µm at X = 0.18, 11 per cent. The points do not sit exactly on the line because each band has its own slopes.

A background that stops breaking into sources

A deep enough image turns a glowing sky into a list of galaxies. Past a certain depth it stops doing so, and the depth is set by how many sources share one beam — which a longer exposure cannot change.

7 figures · Olbers's paradox The histogram of a confused map, for three faint-end slopes. The distribution of pixel values in a simulated map in which every beam holds 0.2 sources brighter than the break in the counts, for counts identical above the break (falling as S^−3.8) and falling as S^−b below it with b = 1.3, 1.5, 1.7. Pixel values are in units of the break flux per beam, measured from the map's own mean as a real map's are and no instrumental noise has been added. No single source below the break can be picked out of any of these maps, and yet the histograms differ: b = 1.3 gives a peak 0.15 below the mean, a width of 0.566 and a skewness of 2.5; b = 1.5 gives a peak 0.15 below the mean, a width of 0.629 and a skewness of 1.6; b = 1.7 gives a peak 0.13 below the mean, a width of 0.617 and a skewness of 1.5. Shallower faint counts leave a map dominated by the rarer bright sources, with a sharp peak and a long positive tail; steeper ones fill each beam with many faint sources and the histogram moves towards a Gaussian. The shape of the histogram is therefore a measurement of the counts below the level at which anything can be counted.

The histogram of a sky too crowded to count

When every beam holds more than one faint source, a map stops being a list. It remains a measurement — the histogram of its pixels carries the number of sources nobody can pick out, down to the depth where there are so many that their sum turns Gaussian.

7 figures · Olbers's paradox

Exoplanets

Planets nobody has seen, weighed and measured from a dip, a wobble and a delay.

Small planets and giant planets per star, against the mass of the star. Planets per star against the host star's mass, both on logarithmic axes, for small planets (one to four Earth radii, orbital periods under 200 days) and giant planets (above a third of Jupiter's mass, within 2.5 astronomical units). The small-planet rate is drawn as 1.1 per star at one solar mass rising as M^−0.6; the giant rate as 7 per cent at one solar mass falling towards small stars as M^1. Both are representative of published fits, not fits of their own. At 0.3 solar masses a star has 2.3 small planets and a 2.1 per cent chance of a giant; at one solar mass 1.1 and 7 per cent. The ratio of small planets to giants changes by a factor of 7 between the two, so a statement about how common planets are depends on which star it is about as much as on which planet.

More small planets round smaller stars

A red dwarf has about twice as many small planets as the Sun's kind of star and a third as many giants. The two populations move in opposite directions along the same axis, and the small planets were built from discs several times poorer.

5 figures · Occurrence rates The chance that a candidate is a planet, against how strong its signal is. The reliability of a transit candidate — the fraction of candidates at a given detection statistic that are planets rather than false alarms — against the multiple-event statistic, the combined signal-to-noise of all its transits, for orbital periods of 10, 100, 350 days. Below the 7.1 threshold nothing is a candidate. At a statistic of 9 the reliability is 100 per cent, 80 per cent, 23 per cent; at 12 it is 100 per cent, 95 per cent, 44 per cent. The true planets and the false alarms are modelled, not fitted: the planets as a falling power law above the survey's completeness curve, the false alarms as an exponential tail above the threshold that grows longer and taller with period, because a long-period signal is two or three transits and two or three transits average away little of the non-Gaussian noise that makes false alarms. At short periods a candidate near threshold is nearly always real; at a year, a candidate near threshold is usually not.

At a year, most candidates are noise

A transit survey's threshold is set so that Gaussian noise yields less than one false planet in the whole catalogue. The noise is not Gaussian, and at year-long periods the candidates just above the threshold are more often artefacts than planets — so the rate of Earth analogues is partly a measurement of the noise.

7 figures · Occurrence rates Giant planets per logarithmic interval of orbital distance, and who can see where. Giant planets (a third of Jupiter's mass to thirteen Jupiters) per star per unit of the natural logarithm of orbital distance, against distance, round Sun-like stars. The curve is a broken power law representative of the long radial-velocity surveys: rising as a^0.8 to a peak at 3 AU and falling beyond it as a^−1, which integrates to 2.7 per cent of stars inside 1 AU and 12 per cent in all. The shaded bands are where each method has the sensitivity to measure it: radial velocities out to about 7 AU, where two decades of monitoring cover one orbit, and direct imaging beyond about 10 AU, where a young giant can be separated from its star. The peak sits near the outer edge of the radial-velocity window, where a planet has completed one or two orbits in the data, and the fall beyond it lies entirely outside that window — so the peak's position is where the two methods have to be joined rather than a place either can measure alone.

A peak that neither survey can see alone

Giant planets become commoner with distance from their star out to a few astronomical units, and then rarer. The rise is measured by radial velocities and the fall by direct imaging, and the peak between them sits where neither method is comfortable.

6 figures · Occurrence rates

Galaxies

Where the unknown stops being a number and becomes a profile — and most of it is not light.

Three initial mass functions, normalised to the same light. The number of stars per logarithmic interval of mass, on a logarithmic axis, for three initial mass functions — Kroupa, Salpeter, bottom-heavy — each normalised to the same number of stars above one solar mass, so that the stars which make an old population's light are the same in each. They differ only below a solar mass. In a 10 billion year old population the stars below half a solar mass hold 62 per cent, 78 per cent, 91 per cent of the living stellar mass, and they give 15.8, 24.2, 41.7 per cent of the light of the main sequence, which is itself a third of the population's. The difference between them is almost entirely mass that emits almost nothing.

Stars that weigh and do not shine

Most of the mass in an old stellar population is in stars too faint to contribute to its light. How many of them there are is set by the mass function stars are born with, and the mass-to-light ratio every dark-matter measurement divides by can be wrong by a factor of two without a single colour changing.

5 figures · Dark matter

Every sequence · Every named object · Search