Astronomy & Celestial Mechanics
the sky, computed
Concept

Adiabatic invariant — where it appears

A quantity preserved when a system's parameters change slowly compared with its own oscillation. It is what makes a resonance protective, and its failure when two resonances overlap is the mechanism that produces chaos.

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

The last curve across the cylinder, before and after it breaks. Surfaces of section for the standard map at K = 0.5, 0.9716, 1.4, each 22 trajectories iterated 41 times from a column of starting values. Nothing here is placed: every dot is an iterate. At K = 0.5 the islands are separate and the space between them is filled with curves that run all the way round in θ — a trajectory cannot get from one island to the next, and one started beside the hyperbolic fixed point wanders 2.91 in p — 0.46 of a cylinder — and stops. At K = 0.9716, which is Greene's threshold to four figures, the last of those curves is on the point of going and the same trajectory still only reaches 4.84, which is 0.77 of a cylinder. At K = 1.4 it covers 6.1 cylinders: the barrier is gone and there is nothing left to stop it. That transition is what a chaotic zone in the asteroid belt is, drawn without any asteroids.

Where the chaos comes from

A Lyapunov time says how long a prediction lasts. It does not say what destroyed it. The mechanism is two resonances whose libration widths overlap, and the transition can be watched happening on a surface of section as one number is turned up.

gravitation · Chaos Libration and circulation of a resonant angle. The critical angle of a resonance over time, integrated from the pendulum equation it obeys. Below the separatrix the angle oscillates about a fixed value — the body is locked; above it, the angle runs without bound and the body is not.

A chain that could not have been assembled in place

Seven planets whose successive period ratios are all close to small whole numbers. Capture into a resonance requires the orbits to converge slowly, and converging slowly is something planets can only do in a disc.

exoplanets · Resonance 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.

gravitation · Impulse approximation A resonance keeps what convergence brings it and releases what divergence takes away. The resonant angle of a body inside a resonance whose strength is changing, for the two signs of that change. The angle obeys a pendulum, and the strength of the pendulum is set by how close the two orbits are; migration changes it slowly compared with the swing, which is the condition under which the area a trajectory encloses is conserved. Convergent migration strengthens the resonance, so the separatrix grows around a trajectory of fixed area and the swing narrows — from 1.05 radians to 0.77 across the figure, the body ending more deeply locked than it began. Divergent migration weakens it, the separatrix shrinks, and it eventually passes inside the trajectory: the angle stops oscillating and begins to run, 37.5π in the last quarter of the drawing alone, and the lock is gone for good. Nothing here is dissipative and nothing is random. The whole asymmetry is the sign of one derivative, which is why a chain of planets in resonance is direct evidence that they migrated toward one another, and why a chain cannot survive a phase in which they moved apart.

Capture is a direction, not a strength

A resonance holds a body that drifts into it from one side and lets go of one that drifts out the other way, and the asymmetry is not about how strong the resonance is. It is the sign of a derivative — whether the trapped region is growing or shrinking — which is why a chain of planets in resonance is direct evidence that they migrated toward each other.

exoplanets · Resonance Relativity switches the cycle off, halving its reach at a ratio of 0.80. The greatest eccentricity a Kozai–Lidov cycle reaches, against the strength of the orbit's own relativistic pericentre precession, measured in units of the cycle's own precession rate at zero eccentricity. The horizontal axis is logarithmic and spans three decades. At the left the relativistic term is negligible and the cycle reaches 0.838, which is the closed-form value for a start at 65 degrees and is what the integration is checked against. At the right it is gone. The mechanism depends on the pericentre staying put while the outer body pulls on the same side of the orbit for a whole cycle, and the relativistic precession is a competing rotation of that same pericentre; when it is faster, the pull averages away. The threshold sits near one by construction and the transition is sharp rather than gradual, with the reach halved at 0.80. What makes it matter is where the relativistic term is largest: it grows as the pericentre falls, so it strengthens exactly as the cycle drives the orbit inward, and it therefore sets a floor on the pericentre distance that this mechanism can deliver a body to.

The precession that switches the cycle off

A distant companion can trade an orbit's inclination for its eccentricity, over and over, and drive a pericentre almost onto the central body. General relativity's own precession competes with the mechanism for the same pericentre, and when it wins the cycle stops — sharply, at a ratio of one, which puts a floor on how close anything can be delivered.

orbits · Kozai–Lidov Prograde and retrograde, the same encounter. The same encounter twice: a companion of 1 times the primary's mass passing at 1.5 disc radii, with the disc spinning the same way the companion orbits and then the other way. Nothing else differs. The prograde disc grows a bridge towards the companion and a tail away from it; the retrograde one is barely disturbed. The reason is a resonance rather than a force: in the prograde case the outer particles keep pace with the companion for a substantial part of the encounter and are pulled the whole time, and in the retrograde case they sweep past it and the impulses cancel.

Whether it merges is a ratio of two times

Two galaxies passing each other either merge or do not, and what decides it is not how close they come. It is whether the encounter lasts long enough for the internal motions to respond — a slow, prograde passage transfers orbital energy into stellar orbits and the pair is bound; a fast one leaves both galaxies heated and still moving.

galaxies · Galaxy interactions Three motions, three decades apart, and that is the point. The three periods of a trapped 1-MeV electron's motion against the shell it is trapped on, on a logarithmic time axis. It gyrates about a field line in 0.22 milliseconds, bounces between its mirror points in 0.33 seconds, and drifts right round the planet in 16 minutes — ratios of 1498 and 3015 at L = 4. Each motion carries a conserved quantity: the magnetic moment, the longitudinal invariant, and the magnetic flux the drift shell encloses. The separation is what makes them conserved. An invariant survives anything that changes slowly compared with its own period, so a disturbance lasting minutes destroys the third and leaves the first two untouched — and a particle that keeps its magnetic moment while being moved inward to a stronger field must gain energy. That is not a loophole; it is how the belts are filled.

Three clocks and nothing to fall onto

A charged particle in a dipole field gyrates, bounces and drifts, on timescales a millisecond, a second and a quarter of an hour. The three periods are three decades apart, and that separation is not a curiosity — it is the reason each motion has a conserved quantity, and the reason a magnetic storm can accelerate particles rather than merely stir them.

spaceflight · Radiation belts 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 The same well after the primary has lost 45 per cent of its mass. Two effective potentials for one body: the solid curve before the primary loses mass and the faint one after, both at the same angular momentum, because a central force of any strength exerts no torque. The well shallows and its floor moves out from r = 1.00 to 1.82. The body's own level moves with it — from E = -0.420 to -0.127 — and the two horizontal lines are drawn where the radial action is conserved, which puts the turning points at 0.71–1.67 before and 1.30–3.03 after. The ratio between them is 2.3333 in both, so the orbit is the same shape at a larger size: everything about the body's path has scaled and its eccentricity of 0.4 has not moved. That is what a slow change leaves behind, and it is not what a sudden one leaves.

The well moves, and the body keeps its share of it

When the Sun becomes a white dwarf it will throw away half its mass, and every planet's orbit will swell by the same factor. Their eccentricities will not change at all — provided the loss is slow, and the only meaning "slow" has here is slow compared with one orbital period.

orbits · Effective potential A cusp pulled deeper by the disc inside it. The dark-matter density of a halo of 10¹² solar masses and concentration 10 before and after an exponential disc of 5·10¹⁰ solar masses and scale length 3 kpc settles slowly into its centre, computed shell by shell from the adiabatic invariant r·M(r). Far outside the disc the two profiles coincide — 4.8% apart at 83 kpc — because contraction moves mass inward without changing how much there is. Inside, the dark matter is pulled in after the baryons: at half a disc scale length the density is 2.9 times higher, and at 1 kpc the logarithmic slope goes from −1.09 to −1.25. The cusp that dark-matter-only simulations predict is therefore the shallowest the inner halo of a disc galaxy can be before feedback acts, not the profile feedback starts from. The calculation assumes circular orbits and slow settling, and simulations that include gas find the true response weaker than this.

The halo a galaxy drags inward

The cusp that dark-matter simulations predict is a halo with no galaxy in it. When gas cools and settles into a disc at the centre, the dark matter is pulled in after it, and a quantity conserved on slowly changing orbits says by how much. In a large disc galaxy the pull triples the dark matter at the centre and multiplies tenfold the energy needed to dig a core; in a dwarf it barely registers.

galaxies · Halo profiles A band and a strip, and the constant where they cross. The plane of the two quantities a planetary ephemeris cannot separate: the fractional rate of change of the gravitational constant, vertically, and of the Sun's mass, horizontally, both in units of 10⁻¹⁴ per year. The orbits respond only to their sum, the rate of change of GM☉, so a measurement of that sum — drawn here as −9.0 ± 3.0 × 10⁻¹⁴ per year, the precision recent ephemerides reach, not any one published value — is a diagonal band. On its own it allows any Ġ/G at all, provided the Sun's mass is changing to compensate. The vertical strip is the Sun's mass loss computed from its luminosity and its wind, −9.9 to −8.8 × 10⁻¹⁴ per year. Where the two overlap, Ġ/G lies between −3.2 and 3.9 × 10⁻¹⁴ per year. The bound on a change in a fundamental constant is therefore set as much by how well the solar wind is measured as by how well the planets are.

A constant measured through a Sun that is getting lighter

Whether the gravitational constant changes with time is tested by watching the planets, and the planets report only the product of G with the Sun's mass. The Sun is losing mass — four million tonnes a second as light, a million and a half more as wind — at nearly a ten-trillionth of itself a year, and every bound on a changing G is a statement about the solar wind as much as about gravity.

gravitation · Gravitational constant A survival probability with a step in the middle of the spectrum. The probability that an electron neutrino made in the core of the Sun is still an electron neutrino at the Earth, against its energy, from the resonant conversion it undergoes in the Sun's own electrons: mixing angle sin²θ₁₂ = 0.307, splitting Δm² = 7.4·10⁻⁵ eV², and an electron density at the point of production of 90 moles per cubic centimetre. At low energy the matter potential is small beside the vacuum term and the probability is the vacuum average, 0.54 at 0.3 MeV (the flat limit is 0.55). At high energy the matter term dominates, the neutrino is born almost entirely in the heavier mass state, and the probability falls to 0.32 at 10 MeV, close to sin²θ₁₂. The two terms are equal at 2.1 MeV, which falls in the gap between the ⁷Be line and the bulk of the ⁸B spectrum. The points are the measured survival probabilities for each source separately — pp, ⁷Be, pep and ⁸B from Borexino, and ⁸B from SNO — and every one of them lies within about two standard deviations of the curve: low-energy neutrinos survive about half the time, high-energy ones about a third.

A survival that depends on the energy

A low-energy neutrino from the Sun reaches the Earth still an electron neutrino about half the time, and a high-energy one about a third of the time. The difference is made by the Sun's own electrons, and the energy at which it switches over — near two million electronvolts — sits in the gap between the neutrinos that never meet a resonance and the ones that all do.

stars · Solar neutrinos Two clocks for a radiation-belt electron, at Kp 2. Two times against drift shell L, on a logarithmic scale, for electrons of about 1 MeV. The first is how long radial diffusion takes to move an electron by one Earth radius, the inverse of an empirical diffusion coefficient that grows as the tenth power of L — at Kp 2, 550 years at L = 2, 10 years at L = 3, 21 days at L = 5. The second is how long the electron survives before it is scattered into the loss cone. Inside the plasmapause, here at L = 4.7, that is set by plasmaspheric hiss: years at L = 1.3, 14 days at L = 2.5, 5.2 days at L = 3; outside it a single lifetime of 30 days stands for every other loss. Inside L = 4.8 the electron is lost before diffusion can carry it further in, which is the shaded region; outside, diffusion wins and the belt is fed. The two curves cross at a single place, and that crossing is the inner edge of the outer belt.

The gap between the belts is a race between two clocks

Electrons reach the inner magnetosphere by diffusing inward from the outer edge, and they are lost there to waves that scatter them into the atmosphere. The diffusion slows as the tenth power of distance and the loss does not, so the two clocks cross at one sharp place. Inside it lies the slot, and at megaelectronvolt energies the inner electron belt is not a steady population at all but what the last great storm left behind.

spaceflight · Radiation belts Phase-space density during a storm fed from inside. Electron phase-space density at fixed first and second adiabatic invariants, against the drift-shell parameter L, before a two-day storm (grey) and 0.5, 1, 2 days after it begins, from a model of radial diffusion with a diffusion rate rising as L¹⁰ and a loss lifetime of 4 days. The boundary is held fixed and a source centred at L = 4.5, standing for acceleration by whistler-mode chorus waves, acts for the storm's two days; the density grows a peak at L = 4.5, well inside the boundary, and diffusion then spreads it both inward and outward. A peak in phase-space density inside the outer boundary cannot be produced by diffusion from outside while the boundary is rising — diffusion only smooths gradients — so its presence then is the signature of local acceleration, and its absence is consistent with transport.

A peak that says where the electrons were made

The outer radiation belt refills with electrons of millions of electronvolts within hours of a storm. Whether they were carried inward from the edge of the magnetosphere or accelerated where they are found shows in one place only — the density of electrons in the coordinates their motion conserves, which a detector does not measure directly.

spaceflight · Radiation belts One history of mass loss, three orbits, three outcomes. The eccentricity of three planets' orbits against the fraction of the star's mass lost so far, as the star sheds 45 per cent of its mass in 6 equal pulses at irregular intervals. Each planet starts with e = 0.05. The pulses are the same events in the star's own time, and the planets differ only in their periods, so each pulse lasts 3, 0.3, 0.03 orbital periods respectively. For the innermost planet, to which every pulse is slow, the eccentricity ends at 0.050: the orbit swells and keeps its shape. For the outermost, to which every pulse is sudden, each pulse kicks the eccentricity by an amount set by where the planet happened to be, and it ends at 0.298. The middle planet ends at 0.090. The same total loss leaves three different shapes, and only one of them is predicted by either limit alone.

The size forgets how the mass was lost, and the shape does not

A dying star loses half its mass, and a planet round it keeps its orbit's shape if the loss is slow and has it scrambled if the loss is sudden. Real stars lose mass in pulses — each one sudden or slow depending on the planet — and what survives a train of pulses is neither answer: the orbit's size settles to the slow value quickly, while its eccentricity random-walks and keeps a record of the history.

orbits · Effective potential 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 The plane a moon's orbit turns about, between equator and ecliptic. The inclination of the Laplace plane — the plane about which a satellite's orbit precesses — to the Earth's equator, against distance from the Earth in Earth radii, for an obliquity of 23.44°. Close in, the Earth's equatorial bulge dominates and the plane is the equator; far out, the Sun dominates and it is the ecliptic, 23.44° away. The dashed curve uses today's bulge, which makes the transition at 8.4 Earth radii. The solid curve lets the bulge grow with the faster spin the Earth had when the Moon was closer — a 5.3-hour day at 5 Earth radii, from the system's angular momentum alone — which moves the transition out to about 14. The Moon today, at 60.3 Earth radii, has a Laplace plane within 0.0022° of the ecliptic, so its 5.15° inclination is entirely free: nothing forces it.

A tilt the Moon should not have

A moon's orbit precesses about a plane that sits on its planet's equator close in and on the planet's own orbit far out. A Moon born in the Earth's equator and carried outward by the tides should have followed that plane all the way to the ecliptic, because the angle it keeps to it is an adiabatic invariant. It is inclined by five degrees instead, and slow recession cannot have put them there.

orbits · Angular momentum A cusp steepened into a spike by the black hole growing in it. The dark-matter density near the Galaxy's central black hole of 4.3·10⁶ solar masses, in solar masses per cubic parsec on logarithmic axes. The dashed line is the halo's NFW cusp, falling inward as one over the radius. If the black hole grew slowly at the centre of that cusp, the dark matter would have followed it inward on adiabatically shrinking orbits, steepening the cusp inside about a fifth of the black hole's sphere of influence — here 0.4 parsecs, with the sphere of influence at 1.8 — into a spike falling as r to the minus 7/3, measured off the drawn curve. Close in, the dark matter's own annihilation would cap the density at the particle mass over the annihilation rate times the Galaxy's age — for a 100 GeV particle with the thermal cross-section, 2.8·10⁸ solar masses per cubic parsec — the flat plateau. Between the spike's edge and its plateau the density rises by 6 orders of magnitude.

A cusp a black hole would sharpen

A black hole that grows slowly at the centre of a dark-matter cusp drags the dark matter in after it, turning a density that rises as one over the radius into a spike rising as the radius to the minus seven-thirds. Whatever the cusp was, the spike comes out nearly the same — and the annihilation signal it would produce is set less by the particle than by whether the spike survived the stars around it.

galaxies · Halo profiles How far inside the door a day of compression reaches. The innermost drift shell at which phase-space density has halved after a day of compression, against where the compression put the last closed drift shell, for diffusion at Kp = 4, 5, 6, 7. With the door at L = 5 and Kp = 6 the loss reaches L = 3.8; at Kp = 4 it reaches only 4.3. The diagonal is the door itself: everything above the curves and below the diagonal was emptied without ever being opened. A strong storm both brings the door in and speeds the diffusion that drains the shells behind it, and the two together carry the loss a full Earth radius or more inside the shells that were opened.

A belt drained through its outer door

When a storm's solar wind squeezes the magnetopause inward, the outermost drift shells of the radiation belt are cut open and their electrons leave within minutes. The shells inside are untouched by the boundary and by the atmosphere alike, and they empty anyway: the opened region is a sink, and the same diffusion that fills the belt from outside carries electrons out to it. A day of compression drains the belt a full Earth radius inside the shells that were actually opened.

spaceflight · Radiation belts A dropout that partly comes back. The flux of 1 MeV electrons at L = 4.5, relative to before the storm, through a storm whose ring current deepens to Dst = −150 nT over eight hours and decays with a one-day e-folding. Dashed: the adiabatic part alone, with phase-space density unchanged — the flux falls to 18 per cent at the storm's peak and returns completely. Solid: the same with a real loss of 60 per cent of the phase-space density during the main phase, as compression and outward diffusion would cause. The flux falls to 7 per cent and recovers only to 39. Seen in flux at fixed energy, the two causes look the same while the storm is on, and only the recovery — days later — tells how much of the dropout was electrons lost and how much was the field.

A dropout that is only the field

During a storm the ring current weakens the field the radiation belt sits in, and every electron, keeping its adiabatic invariants, loses energy. A detector counting electrons at a fixed energy and place sees the flux fall by half or more — and not one electron has gone anywhere. When the ring current decays, the electrons regain their energy and the flux comes back. Only in the coordinates the electrons themselves conserve can a storm's real losses be told from the field's sleight of hand, and the bigger the storm, the larger the sleight.

spaceflight · Radiation belts How much a dying star stirs its companion's orbit, against their separation. The rms change in orbital eccentricity of a 0.8 M☉ companion when a 2 M☉ star becomes a 0.6 M☉ white dwarf, against the binary's separation before the loss, both logarithmic. Blue: the mass loss alone — twelve thermal-pulse episodes 500 years wide and a superwind 10,000 years wide, each adding a kick of random sign whose size falls off exponentially once it lasts longer than about a sixth of an orbit. Below 100 au the orbit is slow to every episode and the eccentricity is untouched. The pulses switch on near 229 au, and by 3,000 au the change is 0.26 and the superwind too has become sudden. Dashed: the change a velocity kick to the white dwarf of 0.05 or 0.2 km/s would make on its own, which grows as the square root of the separation. A 0.05 km/s kick outweighs the pulses out to 251 au; a 0.2 km/s kick outweighs the pulses at every separation drawn.

A kick that writes over the record of the pulses

A star that becomes a white dwarf loses half its mass in pulses, and a companion far enough away feels each pulse as a sudden jolt to its orbit's shape. The wide binaries of white dwarfs should therefore carry the pulses' record in their eccentricities. They would, except that the white dwarf itself is born moving at a few tenths of a kilometre a second, and a kick that small writes more eccentricity than all the pulses together.

orbits · Effective potential Oceans that a changing G would have boiled. The Earth's surface temperature in a crude climate — equilibrium temperature for a fixed albedo, scaled as the received flux to the quarter power, plus a fixed greenhouse warming of 33 K — against time before the present, if the Sun's luminosity goes as G⁷ and the Earth's orbit as 1/G, so the flux goes as G⁹, on top of the Sun's own brightening. With Dirac's G ∝ 1/t and an age of 13.8 Gyr, the oceans boil 1.9 Gyr ago; with G falling by 10⁻¹¹ of itself a year, the surface is 292 K at 3.8 Gyr. Rocks of that age were laid down in water. Constant G gives the faint young Sun's cooling, 270 K in this climate — a separate problem, which the real atmosphere solved.

A constant bounded by the oceans

In 1937 Dirac proposed that the gravitational constant falls as the universe ages. Eleven years later Teller pointed out what that would have done to the Earth: a star's luminosity depends on G to about the seventh power, the Earth's orbit shrinks as G grows, and the sunlight reaching the Earth goes as the ninth power. Run backwards, Dirac's G boils the oceans less than two billion years ago. The rocks were laid down in water, so G cannot have changed that fast — the first bound on a changing constant came from geology.

gravitation · Gravitational constant

Named alongside it

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

Geomagnetic stormRadiation beltsDrift shellImpulse approximationRadial diffusionEccentricityImpulsive approximationMagnetopause shadowingMass lossPhase-space densityAngular momentumAsymptotic giant branch

All concepts