Part III · General Relativity — Chapter 3.1

The Equivalence Principle

Gravity is the one force you can switch off by letting go. What refuses to switch off is the whole of the subject.

Where we are

Chapter 2.6 closed Part II by naming the one thing that had not moved. Everything else had. Space and time are a single fabric with a speed limit built into its geometry. Electricity and magnetism are one field seen from two angles. Energy, momentum and stress are one object. And every law we have written keeps its shape when you change who is watching.

Gravity is still a force reaching across empty space, arriving the moment it is sent. Newton's F=GMm/r2F=GMm/r^{2} contains no cc and no delay. Move the Sun and the Earth's orbit responds in the same instant.

Nothing in Part II permits that, and no small repair is available. Chapter 1.1 already showed that patching an instantaneous force law by hand breaks momentum conservation. Chapter 2.3 showed that "the same instant" is not a frame-independent phrase.

So gravity has to be rebuilt. This chapter does not rebuild it. It does something that comes logically first, and that is more persuasive when it is done properly. It argues that gravity can be geometry. That is, the rebuilding is allowed to take the particular form Part III is going to give it, and no other force in nature is eligible for the same treatment.

Here is the route, announced in advance so you always know which step you are on.

  • An experimental fact about two numbers that had no reason to be equal (§1).
  • The immediate consequence: a falling body's path forgets what is falling (§2). A uniform gravitational field can therefore be deleted by relabelling coordinates.
  • Einstein's enlargement of that claim from mechanics to all of physics (§3).
  • The part of gravity that cannot be deleted, which is the relative acceleration of two nearby falling bodies (§4). This is the physical content of the chapter. Everything else is scaffolding.
  • Three consequences computed with no general relativity at all: how big a falling laboratory may be (§5), why light climbing out of a well is reddened (§6), and by how much a ray of light bends past the Sun (§7).
  • The argument that gravity alone is geometrisable (§8).

One warning about that third consequence. The bending calculation comes out wrong by a factor of two. We will say so plainly when we reach it, and pay the debt back in Chapter 3.8.

Tools you'll need  — Chapter 0.1 §3: the derivative as the coefficient of the best linear approximation. Every appearance of a "tidal" effect below is that sentence applied to a force field. Chapter 0.6 §1 and §6: partial derivatives, the matrix of second partial derivatives (the Hessian), and Taylor's theorem to second order in several variables. Chapter 0.7 §7: the Laplacian 2\nabla^{2}, and the fact that Newtonian gravity obeys 2Φ=4πGρ\nabla^{2}\Phi=4\pi G\rho, which is Poisson's equation, so that 2Φ=0\nabla^{2}\Phi=0 in empty space. Chapter 1.1 for what is unsatisfactory about forces. Chapter 2.3 §5: proper time as the length of a worldline. Chapter 2.5 §7.2: the longitudinal Doppler shift, of which we need only the first-order form. Chapter 1.2's worked example on geodesics is the picture waiting at the end of this one. We will point at it and not spend it.

1 · Two masses that did not have to be equal

Let's begin with a piece of bookkeeping that Newton himself noticed and could not explain. Write down the two places the word "mass" appears in Newtonian physics, and notice that they are two different words wearing the same name.

The first is in the second law. Push on something and it resists:

F  =  mIa. \vv F \;=\; m_{I}\,\vv a . (3.1.1)

The number mIm_{I} is inertial mass. It is defined by an experiment involving no gravity whatsoever. A spring, a puck on ice, a collision: any of those will do. What it measures is one thing only, which is how reluctant the object is to have its velocity changed.

The second is in the law of gravitation. Put an object near a mass MM and it is pulled:

F  =  GMmGr2r^. \vv F \;=\; -\,\frac{GMm_{G}}{r^{2}}\,\hat{\vv r}. (3.1.2)

The number mGm_{G} is gravitational mass. It is defined by an experiment involving no accelerations at all. Hang the object from a balance and read the deflection. What it measures is a completely different thing, which is how strongly the object couples to the gravitational field. It is a charge, the exact analogue of the qq in F=qE\vv F=q\vv E.

There is no reason within Newtonian physics for these two numbers to be related. Electric charge supplies the control case. A proton and an electron have charges equal in magnitude, and inertial masses differing by a factor of 18361836. So charge and inertia are demonstrably independent properties of matter. Nothing in (3.1.1) or (3.1.2) forbids gravitational charge from being independent too.

Now let's combine them. We want the acceleration of a falling body, so substitute (3.1.2) into (3.1.1) and divide both sides by mIm_{I}:

a  =  mGmIGMr2r^. \vv a \;=\; -\,\frac{m_{G}}{m_{I}}\,\frac{GM}{r^{2}}\,\hat{\vv r}. (3.1.3)

Let's read what that says. The acceleration of a falling body is proportional to its ratio of gravitational to inertial mass. Suppose that ratio varied from substance to substance. Then a lead ball and a wooden ball released together would separate visibly as they fell, and the whole of what follows would be impossible.

⚑ Quoted, not derived — the universality of free fall

The equality mG=mIm_{G}=m_{I} is an experimental result and nothing else. It is not derivable from anything in this book, and in Newtonian physics it is an unexplained numerical coincidence. The standard way to report a test is the Eötvös parameter for two materials AA and BB dropped in the same field,

η    2aAaBaA+aB, \eta \;\equiv\; 2\,\frac{a_{A}-a_{B}}{a_{A}+a_{B}},

which is zero exactly when the two ratios mG/mIm_{G}/m_{I} agree. ⚑ The record, in order:

TestBound on η\abs{\eta}
Eötvös torsion balance, 1890–19225×109\sim 5\times10^{-9}
Modern laboratory torsion balances2×1013\sim 2\times10^{-13}
Lunar laser ranging (Earth and Moon falling toward the Sun)1×1013\sim 1\times10^{-13}
MICROSCOPE satellite, final result 2022 (titanium vs platinum)3×1015\sim 3\times10^{-15}

Three parts in a thousand million million. That is among the best-tested statements in physics. It is worth being blunt about the status of that number, because everything Part III builds stands on it. If a material were ever found with η0\eta\neq0, the geometrical picture would be wrong in detail. That is why the experiment keeps being repeated with better apparatus.

Given that result, we may as well choose our unit of gravitational charge so that the ratio is exactly one. From here on there is one mass, written mm, and it cancels out of (3.1.3) altogether:

a  =  g(r),g(r)    GMr2r^. \vv a \;=\; \vv g(\vv r), \qquad \vv g(\vv r)\;\equiv\;-\,\frac{GM}{r^{2}}\,\hat{\vv r}. (3.1.4)

Let's pause and look at what happened there. The left-hand side is kinematics, the acceleration of this object. The right-hand side is a field, a property of the point in space and of nothing else. The object has vanished from its own equation of motion.

In plain terms 3.1.1

Rebuilding gravity is not optional, for the reason the previous part left standing, and the rebuilding starts not with a new equation but with an old and very strange piece of arithmetic. It is an equality between two numbers that had no reason to be equal, noticed by the author of the theory it sits in and never explained by it.

Mass enters the older physics twice, in roles unrelated to each other. In the first it is stubbornness, the reluctance of a thing to have its motion changed, measured with springs and collisions where gravity plays no part. In the second it is a kind of charge, the strength with which a thing answers gravitational attraction, measured by hanging it from a balance and reading a deflection. Electric charge shows that such a pairing need not hold, since a proton and an electron carry equal charge while differing in stubbornness by a factor of nearly two thousand. For gravity the two numbers agree to a few parts in a thousand million million.

The consequence is a cancellation, and it is the whole foundation of what follows. Because one number governs both the pull and the resistance, it divides out, and what remains describes the place rather than the object standing in it. Everything built here rests on that measured coincidence, which is why the experiment keeps being repeated with better apparatus.

2 · The path forgets what is falling

That cancellation has a consequence we should draw out at once. Equation (3.1.4) is a second-order differential equation for r(t)\vv r(t) with no free parameters describing the body. Chapter 0.8 §1 established what such an equation delivers: given a starting position and a starting velocity, the solution exists and is unique. Therefore:

The weak equivalence principle (WEP)

Two bodies released from the same event with the same velocity follow the same trajectory, whatever they are made of. Composition, internal structure, temperature, chemical binding and total mass are all absent from the equation that determines the path.

Nothing about that is true of any other force, and it is worth seeing the contrast immediately. Put a charge in an electromagnetic field and Chapter 2.6 §8 gives

a  =  qm(E+v×B), \vv a \;=\; \frac{q}{m}\Big(\vv E + \vv v\times\vv B\Big), (3.1.5)

Here the factor q/mq/m stubbornly refuses to cancel. It is 1.759×1011 Ckg1-1.759\times10^{11}\ \mathrm{C\,kg^{-1}} for an electron, 9.579×1079.579\times10^{7} for a proton, 4.822×1074.822\times10^{7} for a helium nucleus and exactly zero for a neutron. Release those four particles side by side from rest in the same electric field and they follow four different curves, one of which is a straight line. So there is no such thing as "the trajectory through this field". There is only "the trajectory of this particle through this field."

2.1 · Deleting a uniform field

Here is the first payment universality makes, and it is a large one. Suppose the field is uniform: g\vv g is one and the same constant vector everywhere. Define new coordinates by

r    r    12gt2. \vv r^{\,\prime} \;\equiv\; \vv r \;-\; \tfrac{1}{2}\,\vv g\,t^{2}. (3.1.6)

Notice what that is and what it is not. It is a change of coordinates and nothing more, an instruction for relabelling each event (t,r)(t,\vv r) with a new spatial address, using a rule that depends on tt. We want to know how a free body moves in the new labels, so differentiate the definition twice with respect to tt. The first differentiation gives r˙=r˙gt\dot{\vv r}^{\,\prime}=\dot{\vv r}-\vv g t. The second gives

r¨  =  r¨    g  =  gg  =  0, \ddot{\vv r}^{\,\prime} \;=\; \ddot{\vv r} \;-\; \vv g \;=\; \vv g - \vv g \;=\; 0, (3.1.7)

The last step used (3.1.4). In the primed coordinates every freely falling body moves in a straight line at constant speed. The gravitational field is not weak in these coordinates, and it is not approximately absent. It is gone.

Now notice the step that made this work, because it is easy to read past. The substitution (3.1.6) contains no reference to the body at all. No mm, no qq, no composition. So that single relabelling removes the field from the equation of motion of every body at once.

Try the same trick on (3.1.5) with B=0\vv B=0 and constant E\vv E. The shift you would need is 12(q/m)Et2\tfrac12 (q/m)\vv E t^{2}, which is a different relabelling for each species of particle. A coordinate system is not permitted to depend on what you are looking at, so the trick fails. It fails for exactly the reason §1 identified.

Run the argument backwards and it says something equally strong. Start in an inertial frame with no gravity, and change to coordinates r=r+12a0t2\vv r^{\,\prime\prime}=\vv r+\tfrac12\vv a_{0}t^{2}. Those are the coordinates an observer accelerating uniformly at a0-\vv a_{0} would naturally use. Differentiating twice, free particles now obey r¨=a0\ddot{\vv r}^{\,\prime\prime}=\vv a_{0}. A uniform gravitational field has appeared out of nothing.

So a uniform field and a uniformly accelerated frame are not merely similar. As far as (3.1.4) is concerned, they are the same statement written in two coordinate systems.

⚠ Why this isn't obvious

It is tempting to conclude that gravity is therefore "not real", an artefact of a badly chosen frame, like the centrifugal term Chapter 1.2 §7 dissolved. That conclusion is false, and §4 is entirely devoted to why.

Here is the loophole in the argument above. It assumed a uniform field, and no real gravitational field is uniform. Equation (3.1.4) points toward the centre of the attracting body and weakens with distance, so g\vv g genuinely varies from place to place. Equation (3.1.6) cannot delete a field that is different at different points, because the relabelling would then have to be different at different points. Once you do that, you have changed the geometry rather than merely the labels.

What survives is the local statement. At any single event you may delete the field, and you may keep it deleted over a region whose size §5 computes. The whole of general relativity lives in the gap between "at a point" and "over a region".

In plain terms 3.1.2

The same symbol stands on both sides of the equation for a falling body, so it cancels, and it takes the body's identity with it. What is left determines a path from a starting place and a starting velocity alone, which means that a feather and a cannonball released together in the absence of air trace out one curve rather than two. Every other force keeps a residue of the object in its equation of motion, most visibly electricity, where the surviving ratio of charge to stubbornness varies over three orders of magnitude across ordinary particles and reaches zero for a neutron.

That cancellation buys something startling. If the pull is the same everywhere, then relabelling positions with a rule that accelerates along with the fall removes the pull from the equations entirely, and because the rule mentions nothing about the body, one relabelling does this for everything in the room at the same moment. A field that can be abolished by an act of bookkeeping, for all objects simultaneously, is not behaving like a force.

The reverse reading is worth holding onto. Begin with no gravity at all and describe events using coordinates that accelerate, and a uniform field appears out of nothing. Acceleration and uniform gravity are therefore not two phenomena resembling each other; they are one situation described twice, and nothing measurable distinguishes them.

3 · Einstein's enlargement

Everything in §2 concerned trajectories of test bodies. Einstein's move was to assert that the same holds for every experiment, not merely mechanical ones. That is a large step, and it is worth displaying with its clauses separated. Each clause is testable on its own, and the third is where most of the content sits.

The Einstein equivalence principle (EEP)

In a laboratory falling freely under gravity, confined to a sufficiently small region of space and to a sufficiently short interval of time, the outcome of any local experiment not itself involving gravity is

(i) independent of the composition and structure of the apparatus (the weak principle, §2),
(ii) independent of the velocity of the falling laboratory (local Lorentz invariance),
(iii) independent of where and when the laboratory happens to be (local position invariance).

Equivalently: in such a laboratory, the physics of Part II holds exactly.

Three things must be said about this immediately.

It is a hypothesis, not a theorem. Clause (i) is the experimental result of §1. Clauses (ii) and (iii) are extrapolations from mechanics to optics, electromagnetism, nuclear physics and everything else. They were not derived from anything. They are, however, sharply testable, and §6's redshift is a test of clause (iii). A frequency standard that depended on where it sat in a gravitational field would violate that clause. ⚑ Experiments comparing different atomic clocks at different gravitational potentials constrain such a dependence at the level of a few parts in 10610^{6} of the predicted effect. Local position invariance therefore holds to that accuracy.

The word "local" carries the entire burden. Both qualifications in the statement are doing real work, the small region and the short time alike, and neither can be dropped. Section 5 turns them into an inequality with numbers in it.

There is a stronger version, and the difference matters. The EEP as stated exempts experiments that themselves involve gravity. Drop that exemption, and insist that even the gravitational binding energy of the apparatus falls at the same rate as everything else. What you then have is the strong equivalence principle.

That is not a formality. A body's gravitational self-energy is a genuine contribution to its mass, so a theory in which self-energy gravitates differently from ordinary mass would show the Earth and the Moon falling toward the Sun at slightly different rates. ⚑ Lunar laser ranging bounds that difference at about one part in 10410^{4} of the self-energy contribution. General relativity satisfies the strong principle. Most of its competitors do not, which is why the measurement is made.

In plain terms 3.1.3

Widening a claim about falling stones into a claim about every experiment anybody could perform inside a sealed box is a genuine leap, and it deserves to be recognised as one. What the extended claim says is that a laboratory in free fall, kept small enough and watched briefly enough, is indistinguishable from a laboratory drifting in empty space far from anything, so that chemistry, optics, radioactive decay and the behaviour of clocks all come out as the previous part said they would.

Three separable assertions hide inside that sentence, and separating them is what makes the principle testable rather than rhetorical. The first is the measured fact already in hand, that what the apparatus is built from makes no difference. The second is that how fast the falling laboratory happens to be moving makes no difference. The third, which carries most of the weight, is that where and when the laboratory sits makes no difference, so that no experiment can reveal the local strength of gravity from inside.

None of this was derived, and saying so is not a weakness in the argument but a statement about where its risk lies. Each clause can fail, each has been looked at hard, and the third is checked by the most direct experiment imaginable, which is holding two identical clocks at different heights and asking whether they keep the same time.

4 · What free fall cannot remove

This section is the physical content of the chapter. Everything before it says what gravity is not. This section says what is left over when you have taken away everything a change of coordinates can take away. Here is the destination, announced first:

Where we are going

Two nearby particles, both in free fall, both with the field deleted at the location of the first, do not stay at rest relative to one another. Their relative acceleration is proportional to their separation and to the second derivative of the gravitational potential. It cannot be removed by any relabelling of coordinates, because it is a comparison between two freely falling bodies and refers to no frame at all. That relative acceleration is what gravity is.

4.1 · The tidal equation, in five lines

The plan is to write down the equation of motion for each of two neighbouring particles and then subtract, so that what is common to both drops out. Let r0(t)\vv r_{0}(t) be the position of a reference particle in free fall, and let a second free particle sit at r0(t)+d(t)\vv r_{0}(t)+\vv d(t), where d\vv d is a small separation. Both obey (3.1.4).

Line 1. Start with the reference particle, whose only motion is its own free fall:

r¨0  =  g(r0). \ddot{\vv r}_{0} \;=\; \vv g(\vv r_{0}). (3.1.8)

Line 2. Now the companion. Its acceleration is the field evaluated at its own location, which is displaced from the reference particle by d\vv d:

r¨0+d¨  =  g(r0+d). \ddot{\vv r}_{0} + \ddot{\vv d} \;=\; \vv g(\vv r_{0}+\vv d). (3.1.9)

Line 3. We want the motion of one particle relative to the other, so subtract (3.1.8) from (3.1.9) and keep only the separation. Nothing has been approximated yet:

d¨  =  g(r0+d)    g(r0). \ddot{\vv d} \;=\; \vv g(\vv r_{0}+\vv d) \;-\; \vv g(\vv r_{0}). (3.1.10)

Line 4. The right-hand side is a difference of two nearby values of the same function, so the thing to do is linearise. That is Chapter 0.1 §3 applied in several variables. Writing the ii-th component and using the multivariable linear approximation of Chapter 0.6 §2,

gi(r0+d)  =  gi(r0)  +  gixjr0 ⁣ ⁣dj  +  (smaller than d), g^{i}(\vv r_{0}+\vv d) \;=\; g^{i}(\vv r_{0}) \;+\; \left.\pdv{g^{i}}{x^{j}}\right|_{\vv r_{0}}\!\! d^{j} \;+\; (\text{smaller than } \abs{\vv d}), (3.1.11)

summed over j=1,2,3j=1,2,3. The step is legitimate precisely because d\vv d is small. Everything below is a statement about nearby particles, and nothing more than that is being claimed.

Line 5. Substitute that expansion into (3.1.10). The zeroth-order terms cancel, and that cancellation is the deletion of the field performed in §2.1. What survives is

  d¨i  =  gixjr0dj.   \boxed{\;\ddot{d}^{\,i} \;=\; \left.\pdv{g^{i}}{x^{j}}\right|_{\vv r_{0}} d^{\,j}.\;} (3.1.12)

Let's read (3.1.12) carefully, because it is the sentence Part III is built on. The field g\vv g has disappeared from it. What governs the relative motion of two neighbouring freely falling particles is not the field but the derivative of the field.

Choose coordinates as cleverly as you like and that stays true. The field at a point can always be made zero, since that is three numbers and the shift (3.1.6) has three parameters to spend on them. Its derivative is a different object altogether, and no shift touches it.

4.2 · In terms of the potential

Potentials are easier to differentiate twice than fields are, so let's rewrite the tidal equation in terms of one. Chapter 0.7 §2 established that the Newtonian field is a gradient, g=Φ\vv g=-\nabla\Phi, so gi=iΦg^{i}=-\partial_{i}\Phi and

gixj  =  jiΦ    Φ,ij. \pdv{g^{i}}{x^{j}} \;=\; -\,\partial_{j}\partial_{i}\Phi \;\equiv\; -\,\Phi_{,ij}. (3.1.13)

That array is the Hessian of Chapter 0.6 §6, the matrix of second partial derivatives, and Clairaut's theorem (0.6 §6.1) makes it symmetric. Substituting it, (3.1.12) reads

d¨i  =  Φ,ijdj, \ddot{d}^{\,i} \;=\; -\,\Phi_{,ij}\,d^{\,j}, (3.1.14)

This is a linear equation whose coefficient matrix is symmetric. By the spectral theorem of Chapter 0.5 §6, it therefore has three real eigenvalues and three mutually perpendicular eigendirections. Along each eigendirection the separation obeys d¨=λd\ddot d = -\lambda d, giving exponential growth where λ<0\lambda<0 and oscillation where λ>0\lambda>0.

That is worth saying on its own line. Gravity, seen from inside a falling laboratory, is a symmetric 3×33\times3 matrix.

4.3 · Computing it for a point mass

An abstract matrix is easier to trust once you have seen a concrete one. So take the potential of a point mass, Φ=GM/r\Phi=-GM/r with r=rr=\abs{\vv r}, and differentiate it twice, one step at a time.

First derivative of rr. We will need this repeatedly, so get it out of the way. Since r2=xkxkr^{2}=x^{k}x^{k}, differentiating both sides gives 2rir=2xi2r\,\partial_{i}r = 2x^{i}, so

ir  =  xir    ni, \partial_{i} r \;=\; \frac{x^{i}}{r} \;\equiv\; n^{i}, (3.1.15)

Those are the components of the unit vector pointing outward from the mass.

First derivative of Φ\Phi. With ir\partial_{i}r in hand, the chain rule applied to GMr1-GM\,r^{-1} gives

iΦ  =  +GMr2ir  =  GMxir3. \partial_{i}\Phi \;=\; +\,GM\,r^{-2}\,\partial_{i}r \;=\; \frac{GM\,x^{i}}{r^{3}}. (3.1.16)

(Check: g=Φ\vv g=-\nabla\Phi then points inward with magnitude GM/r2GM/r^{2}, as it must.)

Second derivative. This is the one we actually want. Differentiate (3.1.16) with respect to xjx^{j}, using the product rule on xir3x^{i}\cdot r^{-3}. The first factor gives jxi=δij\partial_{j}x^{i}=\delta_{ij}. The second gives 3r4jr=3xj/r5-3r^{-4}\partial_{j}r=-3x^{j}/r^{5} by (3.1.15). Hence

Φ,ij  =  GM(δijr33xixjr5)  =  GMr3(δij3ninj). \Phi_{,ij} \;=\; GM\left(\frac{\delta_{ij}}{r^{3}} - \frac{3x^{i}x^{j}}{r^{5}}\right) \;=\; \frac{GM}{r^{3}}\Big(\delta_{ij}-3n^{i}n^{j}\Big). (3.1.17)

Before reading anything off, let's give a name to the overall factor, since it is the single number that will appear everywhere below. Note that it has the dimensions of an inverse time squared:

K    GMr3,[K]=s2, K \;\equiv\; \frac{GM}{r^{3}}, \qquad [K]=\mathrm{s^{-2}}, (3.1.18)

Putting (3.1.17) and that definition into the tidal equation (3.1.14), and expanding the bracket so the two pieces stand apart, we get

d¨i  =  K(δij3ninj)dj  =  Kdi  +  3Kni(njdj). \ddot{d}^{\,i} \;=\; -\,K\big(\delta_{ij}-3n^{i}n^{j}\big)\,d^{\,j} \;=\; -K\,d^{\,i} \;+\; 3K\,n^{i}\big(n^{j}d^{\,j}\big). (3.1.19)

Now read off the two cases. Everything hangs on the scalar njdj=n^dn^{j}d^{\,j}=\hat{\vv n}\cdot\vv d, the component of the separation along the line to the mass.

Radial separation. Put d=dn^\vv d = d\,\hat{\vv n}, one particle directly above the other. Then n^d=d\hat{\vv n}\cdot\vv d=d and ni(njdj)=dni=din^{i}(n^{j}d^{\,j}) = d\,n^{i} = d^{\,i}, so the two terms of (3.1.19) combine as Kdi+3Kdi-Kd^{i}+3Kd^{i}:

d¨  =  +2GMr3  d. \ddot{\vv d} \;=\; +\,\frac{2GM}{r^{3}}\;\vv d. (3.1.20)

The sign is positive, so the acceleration points along the separation, away from the reference particle. Radially, the two particles pull apart. The physical reason is plain enough in (3.1.4), since the lower particle is closer to the mass and falls harder. But notice that we did not have to say so. It fell out of the second derivative on its own.

Transverse separation. Put dn^\vv d\perp\hat{\vv n}, the two particles side by side at the same height. Then n^d=0\hat{\vv n}\cdot\vv d=0, the second term of (3.1.19) vanishes outright, and

d¨  =  GMr3  d. \ddot{\vv d} \;=\; -\,\frac{GM}{r^{3}}\;\vv d. (3.1.21)

The sign is negative, so the acceleration points back along the separation, toward the reference particle. Transversely, the two particles come together. Again this is as it must be, since both are falling toward the same centre along converging lines.

Grind box — the tidal matrix in general position, and the size of the next term

1 · The two cases of §4.3 were not special. We evaluated (3.1.19) for separations exactly along and exactly across n^\hat{\vv n}. Show that nothing else happens by diagonalising the matrix Mijδij3ninjM_{ij}\equiv\delta_{ij}-3n^{i}n^{j} directly, for an arbitrary unit vector n^\hat{\vv n}.

Along n^\hat{\vv n}. Apply MM to n^\hat{\vv n} itself, using njnj=1n^{j}n^{j}=1:

Mijnj  =  ni3ni(njnj)  =  ni3ni  =  2ni. M_{ij}n^{j} \;=\; n^{i} - 3n^{i}\big(n^{j}n^{j}\big) \;=\; n^{i}-3n^{i} \;=\; -2\,n^{i}.

So n^\hat{\vv n} is an eigenvector with eigenvalue 2-2.

Across n^\hat{\vv n}. Take any w\vv w with njwj=0n^{j}w^{j}=0. Then the second term dies outright:

Mijwj  =  wi3ni(njwj)  =  wi. M_{ij}w^{j} \;=\; w^{i} - 3n^{i}\big(n^{j}w^{j}\big) \;=\; w^{i}.

So every vector in the plane perpendicular to n^\hat{\vv n} is an eigenvector with eigenvalue +1+1. That plane is two-dimensional, and together with n^\hat{\vv n} it spans all of R3\R^{3}. So we have a complete eigenbasis, and the spectrum is (2,+1,+1)(-2,+1,+1) whatever direction the mass lies in. Multiplying by the overall KK and putting in the minus sign of (3.1.14), the relative acceleration has coefficients (+2K,K,K)(+2K,-K,-K). That is one stretch and two squeezes, always.

Two checks fall out for free. The trace is 2+1+1=0-2+1+1=0, which is (3.1.22). The determinant is (2)(1)(1)=2(-2)(1)(1)=-2, so detΦ,ij=2K3\det\Phi_{,ij}=-2K^{3}, which is not zero. The tidal matrix is therefore invertible, and in particular there is no direction along which two free particles stay put. (Eigenvalues, trace and determinant confirmed numerically for randomly chosen n^\hat{\vv n}.)

2 · What (3.1.24) threw away. The closed-form ellipse treated KK as constant, and it is not. The reference particle falls, rr shrinks and KK grows. Our goal here is the next term. Starting from rest, the reference particle's own free fall gives r0(t)=r12(GM/r2)t2+O(t4)r_{0}(t)=r-\tfrac12(GM/r^{2})t^{2}+O(t^{4}), so

K(t)  =  GMr03  =  K(1Kt22)3  =  K(1+32Kt2+O(K2t4)), K(t) \;=\; \frac{GM}{r_{0}^{3}} \;=\; K\left(1-\frac{K t^{2}}{2}\right)^{-3} \;=\; K\left(1+\tfrac32 K t^{2}+O(K^{2}t^{4})\right),

by the binomial series of Chapter 0.3 §2. We now want the quartic coefficient, so put ξz/ρ=1+Kt2+a4K2t4\xi_{z}/\rho = 1 + Kt^{2} + a_{4}K^{2}t^{4} into ξ¨z=2K(t)ξz\ddot\xi_{z}=2K(t)\xi_{z} and match the coefficient of t2t^{2} on both sides. The left side gives 12a4K212a_{4}K^{2}. The right side gives 2K(32K+K)=5K22K\big(\tfrac32 K + K\big)=5K^{2}, so a4=512a_{4}=\tfrac{5}{12}. The same match on ξ¨x=K(t)ξx\ddot\xi_{x}=-K(t)\xi_{x} gives 112-\tfrac1{12}. Hence

ξzρ=1+Kt2+512(Kt2)2+,ξxρ=112Kt2112(Kt2)2+ \frac{\xi_{z}}{\rho} = 1 + Kt^{2} + \tfrac{5}{12}\big(Kt^{2}\big)^{2}+\ldots, \qquad \frac{\xi_{x}}{\rho} = 1 - \tfrac12 Kt^{2} - \tfrac{1}{12}\big(Kt^{2}\big)^{2}+\ldots

This is exactly the quantity the figure below calls error 2, and it explains the numbers there. At Kt2=0.05Kt^{2}=0.05 the estimate 512(Kt2)2=0.104%\tfrac{5}{12}(Kt^{2})^{2}=0.104\% stands against a measured 0.107%0.107\%. At Kt2=0.25Kt^{2}=0.25 it is 2.60%2.60\% against 3.06%3.06\%. At Kt2=0.55Kt^{2}=0.55 it is 12.6%12.6\% against 19.0%19.0\%, where the series has plainly stopped being a good guide and the numerical integration has to be believed instead.

Note what the correction does not depend on, which is the ring's size. It is an error in the expansion in time, not in the expansion in separation. The figure separates the two so that neither can be mistaken for the other.

Recap — what went in, what came out

In: Newton's law of gravitation, the equality of the two masses, and one linearisation in the separation.

Out: two nearby freely falling particles accelerate relative to each other by 2K2K times their separation when it is radial, stretching them apart, and by KK times their separation when it is transverse, squeezing them together, with K=GM/r3K=GM/r^{3}.

Cost: the result holds only to first order in the separation. Everything here is a statement about neighbouring particles, and Chapter 3.4 will keep exactly that restriction when it does the same calculation covariantly.

4.4 · The trace, and what it is secretly telling you

One number extracted from the tidal matrix carries more than its share of meaning, and it is the simplest one available. Add up the diagonal entries of (3.1.17). Contract ii with jj, using δii=3\delta_{ii}=3 in three dimensions and nini=1n^{i}n^{i}=1 because n^\hat{\vv n} is a unit vector:

Φ,ii  =  GMr3(33)  =  0. \Phi_{,ii} \;=\; \frac{GM}{r^{3}}\big(3-3\big) \;=\; 0. (3.1.22)

Zero. Now look at what that zero is. The quantity Φ,ii\Phi_{,ii} is 2Φ\nabla^{2}\Phi, the Laplacian of Chapter 0.7 §7, and Poisson's equation for Newtonian gravity says 2Φ=4πGρ\nabla^{2}\Phi=4\pi G\rho. So (3.1.22) is not an accident of the point-mass potential. It is the statement that we are computing in empty space, where ρ=0\rho=0.

The tidal matrix of any vacuum region is traceless. Stretch along one direction must be paid for by squeeze along the others, in exactly compensating amounts. The eigenvalues (2K,+K,+K)(-2K,+K,+K) of (3.1.17) sum to zero for that reason and no other.

Hold onto this. Chapter 3.6 will produce field equations whose vacuum form, Rμν=0R_{\mu\nu}=0, is the covariant descendant of (3.1.22), and whose sourced form is the descendant of Poisson's equation.

4.5 · The ring, and the honest arithmetic about its area

Make the shape visible. Release a ring of NN particles, all at rest relative to a central reference particle, arranged in the plane containing the radial direction. Put z^\hat{\vv z} along n^\hat{\vv n} (radially outward) and x^\hat{\vv x} transverse, so the ring is

d(0)  =  ρ(cosθ  x^  +  sinθ  z^),d˙(0)=0. \vv d(0) \;=\; \rho\,\big(\cos\theta\;\hat{\vv x} \;+\; \sin\theta\;\hat{\vv z}\big), \qquad \dot{\vv d}(0)=0. (3.1.23)

The two components decouple, because x^\hat{\vv x} and z^\hat{\vv z} are eigendirections of the tidal matrix. So we can integrate (3.1.20) and (3.1.21) twice from rest, one line each, provided the coefficient KK may be treated as constant over the interval. In general it may not, because the reference particle is itself falling and rr is shrinking. Keeping only the leading behaviour, which means working to first order in the dimensionless combination Kt2Kt^{2},

dz(t)=ρsinθ(1+Kt2),dx(t)=ρcosθ(112Kt2). d^{z}(t) = \rho\sin\theta\,\Big(1+K t^{2}\Big), \qquad d^{x}(t) = \rho\cos\theta\,\Big(1-\tfrac12 K t^{2}\Big). (3.1.24)

Those are the parametric equations of an ellipse, stretched along the field line and squeezed across it. That shape is what a tidal field is, and Chapter 3.4 will recover this very ellipse from the curvature tensor.

Now let's say plainly what was assumed, because the figure below measures it. Equation (3.1.24) carries two approximations, and they are logically independent of each other.

  • One in the separation, since the tidal equation itself is only first order in d\vv d.
  • One in the time, since only the first two terms of the expansion in Kt2Kt^{2} have been kept.

The figure reports them separately so that neither can hide behind the other, and the grind box above computes the leading size of the second.

It is often said that the ring deforms at constant area. Do the arithmetic before believing it. The semi-axes are a=ρ(1+Kt2)a=\rho(1+Kt^{2}) and b=ρ(112Kt2)b=\rho(1-\tfrac12 Kt^{2}), so the enclosed area is

πabπρ2  =  (1+Kt2)(112Kt2)  =  1  +  12Kt2    12K2t4, \frac{\pi a b}{\pi\rho^{2}} \;=\; \big(1+Kt^{2}\big)\big(1-\tfrac12 Kt^{2}\big) \;=\; 1 \;+\; \tfrac12 K t^{2} \;-\; \tfrac12 K^{2}t^{4}, (3.1.25)

That grows, at the same order in tt as the deformation itself. So the slogan is wrong in two dimensions. In three it survives, and the reason is worth seeing. Release a small ball of particles rather than a ring. One direction stretches and both transverse directions squeeze, so the volume ratio is

(1+Kt2)(112Kt2)2  =  1    34K2t4  +  14K3t6, \big(1+Kt^{2}\big)\big(1-\tfrac12 Kt^{2}\big)^{2} \;=\; 1 \;-\; \tfrac34 K^{2}t^{4} \;+\; \tfrac14 K^{3} t^{6}, (3.1.26)

The t2t^{2} term has cancelled identically. A small ball of freely falling particles keeps its volume to second order in time while changing its shape.

And the reason it does is (3.1.22). The rate of fractional volume change is the sum of the three fractional stretch rates, which is the trace, which vanishes in vacuum. The volume statement and Laplace's equation are the same fact.

K t² = 0.250
ρ/r = 1.0e-3
closed-form semi-axes: radial 1.2500 ρ transverse 0.8750 ρ
area 1.093750 × initial volume of the ball 0.957031 × initial
error 1 · exact − tidal equation: 0.0476% of ρ
error 2 · tidal equation − closed form: 3.06% of ρ
that is t = 402.7 s at the Earth’s surface, or 0.434 s at 10 000 km from a 10 M☉ black hole
The shape of gravity, in the frame that has deleted it. The view is from inside the falling laboratory: the central dot is the reference particle, held fixed by construction, and the mass lies far below the bottom of the frame. The faint dashed circle is where the ring started. Three descriptions of the same thing are drawn on top of one another, and the figure exists to separate them. (i) The purple dots come from integrating the exact Newtonian equation r¨=GMr/r3\ddot{\vv r}=-GM\vv r/r^{3} for every particle with a fourth-order Runge–Kutta step and then subtracting the reference particle's own motion. No formula from the text is used to place them. (ii) The thin green ellipse is the tidal equation (3.1.19) integrated numerically, with the coefficient K=GM/r3K=GM/r^{3} recomputed at each step as the reference particle falls. It is exact in time and first order in the separation. (iii) The thick blue ellipse is the closed form (3.1.24), which additionally keeps only the first two terms of the expansion in tt. The two error readouts measure the two approximations separately, each as the largest gap between the corresponding pair of curves. Error 1 is the price of linearising in the separation: shrink the ring and it falls by a factor of ten for every factor of ten in ρ/r\rho/r, running 5.7%0.48%0.048%0.0048%5.7\%\to0.48\%\to0.048\%\to0.0048\% as you drag that slider left. Error 2 is the price of truncating the expansion in time, and it does not move at all while you do so. Push the elapsed slider instead and the second grows steeply, 0.11%3.1%19%0.11\%\to3.1\%\to19\%, while the first barely stirs. The arrows are the relative acceleration Φ,ijdj-\Phi_{,ij}d^{j}, outward at top and bottom, inward at the sides. Finally, compare the area and volume readouts: at Kt2=0.05Kt^{2}=0.05 the enclosed area has grown by 2.4%2.4\% while the volume of the ball obtained by spinning the ellipse about the vertical axis has moved by only 0.18%0.18\%, because the area's change begins at order t2t^{2} and the volume's only at order t4t^{4}. That cancellation is (3.1.22), which is Laplace's equation, which in Chapter 3.6 becomes the vacuum field equations.

4.6 · The thesis of Part III, stated early

Collect what §2 and §4 have between them established.

  • The gravitational field at a point can always be set to zero by a change of coordinates, and can be set to zero for all bodies at once. It is therefore not an invariant feature of the situation. It is a fact about the coordinates, in exactly the sense Chapter 2.4 taught us to be suspicious of.
  • The derivative of the field cannot be set to zero. It shows up as the relative acceleration of two freely falling bodies, which is measured by comparing two objects with each other and never by referring to a frame.

So if gravity is going to be described by something coordinate-independent, that something must be built out of second derivatives of the potential, and it must vanish when and only when the tidal effects vanish.

Chapter 3.4 constructs exactly such an object, a tensor with four indices built from second derivatives of the metric. Chapter 3.4 §4 will then derive an equation of the form ξ¨Rξ\ddot\xi \sim -R\,\xi, which is (3.1.14) with Φ,ij-\Phi_{,ij} replaced by that tensor. Set the two side by side and the identification is unavoidable: tidal acceleration is curvature. That is the thesis of Part III, and it was decided here.

In plain terms 3.1.4

So far the story has been about what falling removes; here is what it leaves behind, and it is the only part of gravity really there. Take two neighbouring specks of dust, released and left alone, and ask how they move relative to each other rather than to a frame. Subtract their equations of motion and the pull cancels, being nearly the same at both places; what survives is the rate at which the pull changes from place to place, times the gap.

Work that out for one attracting body and the picture is worth having. Two specks one above the other drift apart, the lower being nearer and falling harder; two side by side drift together, both falling toward one centre along converging lines. A ring of dust stretches along the pull, narrows across it, becoming an ellipse. This is the tide, and the ocean does it twice a day for the same reason.

The stretching and squeezing balance, and that is no coincidence. A small ball of dust changes shape while holding its volume fixed, and that sentence is the law of gravity in empty space in plain clothes. It is not the phase-space volume theorem in new clothes: that one concerned possible states, this one concerns actual dust, and each is a flow with nothing made or lost. No coordinates remove this relative drift, so whatever gravity ultimately is, this is it.

5 · How big may a falling laboratory be?

The equivalence principle was stated with two qualifications, "sufficiently small" and "sufficiently short". Section 4 supplies everything we need to replace those two words with an inequality. Here is the destination: we want a bound on the size and duration of a region in which the tidal drift stays below the precision of whatever instrument is in the room.

Take the worst case, which is radial separation. From (3.1.20), two particles a distance \ell apart, released at relative rest, have relative acceleration 2K2K\ell. Because \ell barely changes over the interval, the acceleration is very nearly constant, and integrating twice from rest gives a relative displacement after time TT of

δ  =  12(2K)T2  =  KT2. \delta \;=\; \tfrac12\,(2K\ell)\,T^{2} \;=\; K\,\ell\,T^{2}. (3.1.27)

The frame is good enough as long as that drift stays below what the instrument in the room can resolve. Call that resolution ε\varepsilon, which is a length, and demand that the drift stay under it:

  T2  <  εK  =  εr3GM.   \boxed{\;\ell\,T^{2} \;\lt\; \frac{\varepsilon}{K} \;=\; \frac{\varepsilon\,r^{3}}{GM}.\;} (3.1.28)

Three things are worth extracting from (3.1.28).

The constraint is on a product, not on a length. A local inertial frame is a region of spacetime. You may have a large laboratory for a short time or a small one for a long time, and the trade is T2\ell\propto T^{-2}. This is why general relativity is a differential theory. Special relativity is exact in the infinitesimal and only there.

If your precision is fractional rather than absolute, the size drops out. Dividing (3.1.27) by \ell, the fractional drift is δ/=KT2\delta/\ell=KT^{2}, independent of how big the laboratory is. Requiring δ/<εfrac\delta/\ell<\varepsilon_{\text{frac}} then gives T<εfrac/KT<\sqrt{\varepsilon_{\text{frac}}/K}, which is a pure statement about duration.

KK is the local answer to "how curved is it here". Far from everything, K0K\to0 and the frame may be enormous. Close to a compact mass, KK is large and the frame is tiny. And KK is a second derivative of the potential, which is §4's object all over again.

5.1 · Numbers

At the Earth's surface, with M=5.972×1024 kgM=5.972\times10^{24}\ \mathrm{kg} and r=R=6.371×106 mr=R_{\oplus}=6.371\times10^{6}\ \mathrm{m},

K  =  GMR3  =  gR  =  1.541×106 s2, K \;=\; \frac{GM}{R_{\oplus}^{3}} \;=\; \frac{g}{R_{\oplus}} \;=\; 1.541\times10^{-6}\ \mathrm{s^{-2}}, (3.1.29)

where the middle expression follows because g=GM/R2g=GM/R_{\oplus}^{2}. Feeding this into (3.1.28):

LaboratoryResolution ε\varepsilonSize \ellTime available
Optical bench1 μm1\ \mu\mathrm{m}1 m1\ \mathrm{m}0.81 s0.81\ \mathrm{s}
Optical bench1 μm1\ \mu\mathrm{m}1 cm1\ \mathrm{cm}8.1 s8.1\ \mathrm{s}
Atom interferometer1 nm1\ \mathrm{nm}1 cm1\ \mathrm{cm}0.25 s0.25\ \mathrm{s}
Gravitational-wave detector1018 m10^{-18}\ \mathrm{m}4 km4\ \mathrm{km}1.3×108 s1.3\times10^{-8}\ \mathrm{s}

The last row is the interesting one. A gravitational wave is a travelling tidal field, so a detector built specifically to measure tidal effects has a local inertial frame lasting about thirteen nanoseconds. That is not a defect of the instrument. It is the statement that the instrument works.

5.2 · A debt from Chapter 1.1, collected

Chapter 1.1 §1 identified a circularity at the foundation of Newtonian mechanics that no amount of care could remove. An inertial frame was defined as one in which force-free bodies move uniformly. A force-free body was one that moves uniformly in an inertial frame. Newton's escape was to postulate absolute space, an entity he conceded was unobservable. That chapter promised the circle would be broken here, by changing the theory rather than by being more careful.

Here is the break. Freely falling bodies define the inertial frames, and free fall is identifiable by an experiment performed inside, with no reference to anything external. Release two objects and see whether they stay put relative to each other.

Gravity leaves the list of forces entirely. A body under gravity alone is not being pushed. It is the standard of not being pushed. So "force-free" acquires an independent meaning, and the definition stops chasing its own tail.

The price is (3.1.28). The test "do the two objects stay put" succeeds only within the region that inequality allows, so inertial frames now exist locally and not globally.

Whether the local frames can be knitted together into one global frame is a separate question with a definite answer. The obstruction to knitting them is exactly what Chapter 3.4 will call curvature, and Newton's absolute space is the assumption that the knitting always succeeds.

In plain terms 3.1.5

The word local has been carrying a great deal of weight without being asked to pay for it, and it can now be handed a number. Free fall abolishes gravity at one place and one moment exactly; move away from that place, or wait, and the drift computed a moment ago begins to show. Requiring the drift to stay under whatever the instruments in the room can resolve produces an inequality, and the inequality constrains the size of the room multiplied by the square of the time spent in it.

That the constraint falls on a product rather than on a length is the useful part. A falling laboratory is not a small box; it is a small patch of space and time together, and one may be traded for the other. A metre-wide bench near the Earth's surface behaves as though gravity had been switched off for about a second if you can measure to a thousandth of a millimetre, and for thirteen nanoseconds if you are working at the precision of a gravitational-wave detector.

Underneath the numbers sits a structural point about the shape the eventual theory must take. The physics of the previous part is exact only in the limit of a vanishing region, so what replaces Newtonian gravity has to hold in the small and be stitched together across large regions rather than written down globally at a stroke.

6 · Light climbing out of a well

Now three consequences, each obtained with nothing but the equivalence principle and special relativity. No curvature, no metric, no field equations. Here is the destination for this one: a light signal sent upward through a height hh in a gravitational field gg arrives with its frequency lowered by the fraction gh/c2gh/c^{2}, and identical clocks at the two ends therefore run at different rates.

6.1 · The cabin

Put a sealed cabin in deep space, far from any mass, and accelerate it with constant proper acceleration gg along its own axis, which we call upward. A source is bolted to the floor and a detector to the ceiling, a height hh above it. This is a problem in special relativity alone. There is no gravity anywhere in it.

Work in the inertial frame SS in which the cabin is momentarily at rest at the instant of emission, and call that instant t=0t=0. The source emits light of frequency νem\nu_{\text{em}} as measured by the source's own clock, which at t=0t=0 is a clock at rest in SS.

Familiar ground — this cabin was drawn in Chapter 2.2

The worldline of a body with constant proper acceleration gg was computed in Chapter 2.2's Problem 4, the rocket at constant proper acceleration, and it came out a hyperbola, (x+c2/g)2c2t2=(c2/g)2\big(x+c^{2}/g\big)^{2}-c^{2}t^{2}=\big(c^{2}/g\big)^{2}. That chapter promised the crew would find themselves in something remarkably like a gravitational field, with a horizon behind them. Both halves of the promise are now collectable.

The field is the content of §2.1. The cabin's own coordinates are precisely r=r+12gt2\vv r^{\,\prime\prime}=\vv r+\tfrac12 \vv g t^{2}, in which free particles fall.

The horizon is a property of the hyperbola. Its asymptotes are the null lines x=±ctc2/gx=\pm ct-c^{2}/g, and a light signal emitted from behind the asymptote never catches the cabin, however long it chases. The cabin's speed approaches cc and the gap never closes. So a permanently accelerated observer has a region of spacetime from which no signal can ever reach them, at a distance c2/gc^{2}/g behind. For g=9.81 ms2g=9.81\ \mathrm{m\,s^{-2}} that is about a light-year.

The equivalence principle then says a static observer in a gravitational field should have the same feature. It does not say where the horizon is, because gg is not uniform and c2/gc^{2}/g is only the answer for the uniform case. Chapter 3.8 §6 locates it properly, at r=2GM/c2r=2GM/c^{2}, and shows that it is a coordinate artefact rather than a place where anything is singular. The Rindler horizon of Chapter 2.3's Problem 3 is an artefact in exactly the same way, produced by insisting on riding the hyperbola forever.

The flight time. In SS the light travels at cc and must cover a distance of hh plus whatever the ceiling has moved. So

tarr  =  hc  +  1c12gtarr2ceiling’s rise  =  hc(1+gh2c2+). t_{\text{arr}} \;=\; \frac{h}{c} \;+\; \frac{1}{c}\cdot\underbrace{\tfrac12 g\,t_{\text{arr}}^{2}}_{\text{ceiling's rise}} \;=\; \frac{h}{c}\left(1 + \frac{gh}{2c^{2}} + \ldots\right). (3.1.30)

We are going to work to first order in the small dimensionless quantity gh/c2gh/c^{2}, and the correction in (3.1.30) is itself of that order. Since it multiplies a quantity that already carries one factor of gg, it contributes at second order and may be dropped. Say this out loud, so it is not mistaken for carelessness: we use tarr=h/ct_{\text{arr}}=h/c and the error is O ⁣((gh/c2)2)O\!\big((gh/c^{2})^{2}\big).

The detector's speed on arrival. Starting from rest in SS and accelerating at gg, the ceiling's velocity when the light reaches it is

v  =  gtarr  =  ghc. v \;=\; g\,t_{\text{arr}} \;=\; \frac{gh}{c}. (3.1.31)

(Proper and coordinate acceleration agree at first order in v/cv/c, which is the order we are keeping. Chapter 2.5 §6 recorded the difference.)

Doppler. The detector is receding from the emission event at speed vv, in the same direction the light is travelling. Chapter 2.5 §7.2 derived the longitudinal Doppler formula, so take it from there and expand it to first order in v/cv/c:

νrecνem  =  1v/c1+v/c  =  1vc+O ⁣(v2c2). \frac{\nu_{\text{rec}}}{\nu_{\text{em}}} \;=\; \sqrt{\frac{1-v/c}{1+v/c}} \;=\; 1 - \frac{v}{c} + O\!\left(\frac{v^{2}}{c^{2}}\right). (3.1.32)

We know the speed already, so put (3.1.31) into that expansion and write the answer as a fractional shift in frequency:

  Δνν    νrecνemνem  =  ghc2.   \boxed{\;\frac{\Delta\nu}{\nu} \;\equiv\; \frac{\nu_{\text{rec}}-\nu_{\text{em}}}{\nu_{\text{em}}} \;=\; -\,\frac{g\,h}{c^{2}}.\;} (3.1.33)

Light climbing inside an accelerating cabin arrives reddened. That is a result in special relativity and it required nothing else.

6.2 · Handing it to gravity

Now invoke the equivalence principle in the direction established at the end of §2.1: a cabin accelerating at gg in empty space and a cabin held stationary in a uniform gravitational field gg are the same situation described in two coordinate systems. Every experiment performed inside must give the same answer. Hence, with no further calculation,

Δνν  =  ghc2for light climbing a height h in a static field g. \frac{\Delta\nu}{\nu} \;=\; -\,\frac{g h}{c^{2}} \qquad\text{for light climbing a height } h \text{ in a static field } g. (3.1.34)

6.3 · Why this is a statement about clocks

Equation (3.1.34) is usually read as "light loses energy climbing out". Resist that reading for a moment, because there is a sharper one available and it costs three lines.

Suppose the source emits NN wave crests over an interval that its own clock records as τem\tau_{\text{em}}, so the emitted frequency is νem=N/τem\nu_{\text{em}}=N/\tau_{\text{em}}. Those crests travel up.

Now use the fact that the situation is static. The field does not change with time, the source and the detector do not move, and the region between them is unchanging. Therefore the number of crests in flight at any moment is constant, and no crest is created or destroyed on the way.

So the detector receives exactly NN crests, over an interval its own clock records as τrec\tau_{\text{rec}}, giving νrec=N/τrec\nu_{\text{rec}}=N/\tau_{\text{rec}}. Divide the one frequency by the other and the count NN drops out:

νrecνem  =  τemτrec  =  1ghc2. \frac{\nu_{\text{rec}}}{\nu_{\text{em}}} \;=\; \frac{\tau_{\text{em}}}{\tau_{\text{rec}}} \;=\; 1 - \frac{gh}{c^{2}}. (3.1.35)

Let's look at what that line is actually saying. The same NN crests took longer to arrive than to leave, as measured by the two local clocks. Neither clock is moving. Neither is defective, since by the equivalence principle both are ideal. The only consistent reading left is that the lower clock runs slow relative to the upper one, by the fraction gh/c2gh/c^{2}:

dτlowerdτupper  =  1ghc2. \frac{\dd\tau_{\text{lower}}}{\dd\tau_{\text{upper}}} \;=\; 1 - \frac{gh}{c^{2}}. (3.1.36)

6.4 · Non-uniform fields, by chaining cabins

Equation (3.1.34) assumed uniformity, which is only good over a height satisfying §5's bound. Remove the assumption by applying the local result over an infinitesimal rise dh\dd h and accumulating, which is Chapter 0.2's business. Over dh\dd h the fractional change in frequency is

dνν  =  g(h)dhc2. \frac{\dd\nu}{\nu} \;=\; -\,\frac{g(h)\,\dd h}{c^{2}}. (3.1.37)

Let's fix the sign convention once, out loud, so it does not have to be checked again. Let hh measure height upward and let Φ\Phi be the gravitational potential, which by (3.1.16) increases with height. The downward field strength is then g=+dΦ/dhg=+\dd\Phi/\dd h, so gdh=dΦg\,\dd h=\dd\Phi. The frequency therefore falls as the potential rises, which is the right direction.

With the sign settled, substitute dΦ\dd\Phi into (3.1.37) and integrate all the way from emitter to receiver:

lnνrecνem  =  ΦrecΦemc2Δνν    ΔΦc2 \ln\frac{\nu_{\text{rec}}}{\nu_{\text{em}}} \;=\; -\,\frac{\Phi_{\text{rec}}-\Phi_{\text{em}}}{c^{2}} \qquad\Longrightarrow\qquad \frac{\Delta\nu}{\nu} \;\simeq\; -\,\frac{\Delta\Phi}{c^{2}} (3.1.38)

to first order in ΔΦ/c2\Delta\Phi/c^{2}. The uniform result is recovered when ΔΦ=gh\Delta\Phi=gh. This form is the one used in practice, because real potentials are not uniform.

6.5 · Numbers, and one very large consequence

SituationΔΦ/c2\Delta\Phi/c^{2}In seconds per year
Top vs bottom of a 1 m1\ \mathrm{m} table1.09×10161.09\times10^{-16}3.4 ns3.4\ \mathrm{ns}
33 cm33\ \mathrm{cm}, the Chou et al. optical-clock experiment ⚑3.60×10173.60\times10^{-17}1.1 ns1.1\ \mathrm{ns}
22.5 m22.5\ \mathrm{m} Harvard tower, Pound and Rebka ⚑2.46×10152.46\times10^{-15}77 ns77\ \mathrm{ns}
GPS orbit vs the ground (26560 km26\,560\ \mathrm{km} vs 6371 km6371\ \mathrm{km})5.29×10105.29\times10^{-10}16.7 ms16.7\ \mathrm{ms}

⚑ Pound and Rebka measured the 22.522.5-metre shift in 1960 using the Mössbauer effect and confirmed it. Later refinements of the same experiment agree with the prediction to about 1%1\%. ⚑ In 2010 a comparison of two aluminium-ion optical clocks resolved the shift over a height difference of 33 cm33\ \mathrm{cm}. Both are quoted as experimental results.

Now the consequence, and it is a large one. Two clocks, at rest with respect to each other, never moving, separated only in height, tick at different rates. In special relativity every rate difference came from relative motion. Here there is no relative motion at all.

So whatever describes spacetime in the presence of a mass, it cannot be ημν\eta_{\mu\nu} everywhere, because ημν\eta_{\mu\nu} assigns the same relation between coordinate time and proper time at every point.

We can go one step further and read off a component of the eventual answer in advance. Chapter 2.3 §5 wrote proper time along a worldline at rest as dτ2=ds2/c2\dd\tau^{2}=\dd s^{2}/c^{2}. Chapter 3.3 will replace ημν\eta_{\mu\nu} by a position-dependent gμν(x)g_{\mu\nu}(x) in that expression, and for a static observer only the 0000 component survives, giving dτ=g00  dt\dd\tau=\sqrt{g_{00}}\;\dd t.

Now compare that with (3.1.38), which says the rate ratio is 1+Φ/c21+\Phi/c^{2} up to a common constant. Matching the two expressions for the rate gives

g00  =  (1+Φc2)2  =  1+2Φc2+O ⁣(Φ2c4). g_{00} \;=\; \left(1+\frac{\Phi}{c^{2}}\right)^{2} \;=\; 1 + \frac{2\Phi}{c^{2}} + O\!\left(\frac{\Phi^{2}}{c^{4}}\right). (3.1.39)

Hold two things about (3.1.39) in mind. First, it will be derived from the field equations in Chapter 3.6 rather than guessed, and the agreement will be a check. Second, and much more importantly for the next section: this argument has said absolutely nothing about the spatial components of the geometry. Clocks probe g00g_{00} and nothing else. Keep that grievance in your pocket, because §7 is where it bites.

In plain terms 3.1.6

Two clocks that never move relative to each other, one on a shelf and one on the floor beneath it, do not agree, and the argument reaching that conclusion uses no general relativity at all. Put a windowless cabin in empty space and accelerate it; send a flash from the floor to the ceiling; the ceiling has picked up speed away from the flash during the crossing, so the light arrives reddened by the ordinary Doppler effect of the previous part. Then apply the principle that a steadily accelerating cabin and a cabin standing still in gravity are one situation described twice.

Reading the result as light losing energy on the climb is the weaker option. The stronger reading counts wave crests. Nothing about the arrangement changes with time, so crests are neither created nor destroyed between floor and ceiling, and exactly as many arrive each second as leave each second. If the receiver nonetheless counts fewer of them in each of its own seconds, the only possibility left is that its seconds and the emitter's seconds are of different lengths.

That is a serious problem for the geometry inherited from the previous part, in which the relation between clock readings and coordinate time is fixed once and holds everywhere. One further caution will matter almost immediately: this argument constrains the timekeeping part of the geometry and says nothing whatever about distances.

7 · Light bending, and a debt of exactly one half

The same cabin, turned through a right angle. Here is the destination in advance, including the bad news: we will obtain a deflection of 2GM/bc22GM/bc^{2} for a ray passing a mass MM at distance bb, which for a ray grazing the Sun is 0.870.87 arcseconds, and the measured value is 1.751.75 arcseconds. Our answer will be exactly half. We will say so and identify what is missing.

7.1 · A ray crossing the cabin

Cabin again, accelerating at gg upward, but now the light enters horizontally through one wall and crosses a width LL to the opposite wall. In the momentarily comoving inertial frame the light travels in a perfectly straight horizontal line, because that frame is inertial and the light is free. It takes time t=L/ct=L/c to cross, again to first order.

During that time the cabin rises by 12gt2\tfrac12 g t^{2}. So relative to the cabin the light has descended by

Δy  =  12gt2  =  gL22c2. \Delta y \;=\; -\,\tfrac12 g t^{2} \;=\; -\,\frac{g L^{2}}{2c^{2}}. (3.1.40)

More useful than the drop is the change in direction. In the inertial frame the light has zero vertical velocity, while the cabin has acquired vertical velocity +gt+gt. So relative to the cabin the light has vertical velocity gt-gt while retaining horizontal velocity cc, and the angle by which its path has turned is

Δθ  =  gtc  =  gLc2. \Delta\theta \;=\; \frac{gt}{c} \;=\; \frac{gL}{c^{2}}. (3.1.41)

Divide by the path length LL to get the result in the form we can integrate. Writing gg_{\perp} for the component of the gravitational field perpendicular to the ray, and \ell for arc length along it, the equivalence principle converts (3.1.41) into

dθd  =  gc2. \dv{\theta}{\ell} \;=\; \frac{g_{\perp}}{c^{2}}. (3.1.42)

A ray bends toward a mass at a rate equal to the transverse field divided by c2c^{2}.

7.2 · Integrating past a mass

Send a ray past a point mass MM with impact parameter bb, meaning the closest approach of the undeflected line. Set up coordinates with the ray along the xx axis and the mass at (0,b)(0,-b), so that a point of the ray at abscissa xx is at distance x2+b2\sqrt{x^{2}+b^{2}} from the mass.

The field magnitude there is GM/(x2+b2)GM/(x^{2}+b^{2}), and its transverse (here, downward) component is obtained by multiplying by the direction cosine b/x2+b2b/\sqrt{x^{2}+b^{2}}:

g(x)  =  GMx2+b2bx2+b2  =  GMb(x2+b2)3/2. g_{\perp}(x) \;=\; \frac{GM}{x^{2}+b^{2}}\cdot\frac{b}{\sqrt{x^{2}+b^{2}}} \;=\; \frac{GM\,b}{(x^{2}+b^{2})^{3/2}}. (3.1.43)

Let's name the approximation before making it. The total deflection comes out at order 10610^{-6} radians, so we may integrate (3.1.42) along the undeflected straight line and take d=dx\dd\ell=\dd x.

This is first-order perturbation theory. You evaluate a small correction along the uncorrected path, so that the piece you neglect is the correction to the correction. The error here is second order in GM/bc2GM/bc^{2}, which for the Sun is one part in 10610^{6} of an already tiny angle.

With that settled, the total deflection α\alpha is the bending rate accumulated along the whole ray, from far before the mass to far after it:

α  =  1c2GMb  dx(x2+b2)3/2. \alpha \;=\; \frac{1}{c^{2}}\int_{-\infty}^{\infty} \frac{GM\,b\;\dd x}{(x^{2}+b^{2})^{3/2}}. (3.1.44)

Do the integral by the trigonometric substitution x=btanux=b\tan u, so that dx=bsec2udu\dd x=b\sec^{2}u\,\dd u and x2+b2=b2sec2ux^{2}+b^{2}=b^{2}\sec^{2}u, whence (x2+b2)3/2=b3sec3u(x^{2}+b^{2})^{3/2}=b^{3}\sec^{3}u. The limits x=x=\mp\infty become u=π/2u=\mp\pi/2:

bdx(x2+b2)3/2  =  π/2π/2bbsec2ub3sec3udu  =  1bπ/2π/2cosudu  =  2b. \int_{-\infty}^{\infty}\frac{b\,\dd x}{(x^{2}+b^{2})^{3/2}} \;=\; \int_{-\pi/2}^{\pi/2}\frac{b\cdot b\sec^{2}u}{b^{3}\sec^{3}u}\,\dd u \;=\; \frac{1}{b}\int_{-\pi/2}^{\pi/2}\cos u\,\dd u \;=\; \frac{2}{b}. (3.1.45)

The integral came out as 2/b2/b. Restoring the constants that were carried outside it gives the total deflection:

αEP  =  2GMbc2. \alpha_{\text{EP}} \;=\; \frac{2GM}{b\,c^{2}}. (3.1.46)

7.3 · The number, and the confession

For a ray grazing the Sun, b=R=6.957×108 mb=R_{\odot}=6.957\times10^{8}\ \mathrm{m} and GM=1.327×1020 m3s2GM_{\odot}=1.327\times10^{20}\ \mathrm{m^{3}\,s^{-2}}:

αEP  =  2(1.327×1020)(6.957×108)(8.988×1016)  =  4.245×106 rad  =  0.875. \alpha_{\text{EP}} \;=\; \frac{2(1.327\times10^{20})}{(6.957\times10^{8})(8.988\times10^{16})} \;=\; 4.245\times10^{-6}\ \mathrm{rad} \;=\; 0.875''. (3.1.47)
⚠ This answer is wrong, by a factor of exactly two

⚑ The measured deflection at the solar limb is 1.751.75''. Eddington's expeditions established it during the eclipse of 1919 at a precision of some tens of per cent. Very-long-baseline radio interferometry has since confirmed the general-relativistic value to about one part in 10410^{4}. Our (3.1.47) is short by a factor of 2.0002.000.

This is not a rounding error, an arithmetic slip, or a subtlety about which radius to use. It is a structural failure of the argument. It was also Einstein's own published prediction in 1911, and he corrected it to 1.751.75'' in 1915 when the field equations were in place. We record it here as a debt, payable in Chapter 3.8.

The reason is visible if you look at what the cabin argument used. Every step of §7.1 was about time. How long the light takes to cross. How much speed the cabin picks up in that time. The whole derivation is a consequence of the fact established in §6, that the rate of a clock depends on where it sits. That is a fact about g00g_{00} alone, exactly as (3.1.39) warned. Nothing in the cabin argument constrains how spatial distances are measured near a mass, because a cabin is a rigid box and we assumed its width was LL without asking any further.

In the full theory the spatial part of the geometry is also modified, by an amount of the same order. A ray of light covers equal amounts of space and time, since it moves at cc, so it picks up an equal contribution from that modification. Two equal contributions, hence the factor two.

A slow-moving planet, by contrast, covers a great deal of time and very little space per unit of proper time, so it barely notices the spatial part. That is exactly why the Newtonian limit is recovered from g00g_{00} alone (Chapter 3.6 §5), and why Newtonian gravity works so well for everything except light.

The register of debts

Two things have now been promised to later chapters, and both will be collected explicitly.

(a) The factor of two in (3.1.46) is paid in Chapter 3.8 §4, which derives α=4GM/bc2\alpha=4GM/bc^{2} from the Schwarzschild solution and names the missing half as spatial curvature.

(b) The identification of tidal acceleration (3.1.14) with a coordinate-independent object is paid in Chapter 3.4 §4.

Being caught out later by a factor of two would cost this book more than admitting it now.

In plain terms 3.1.7

Honesty is cheaper now than it would be later, so here is the result together with what is wrong with it. Turn the accelerating cabin on its side and send a flash across it. In the frame where nobody accelerates the flash goes straight, but the cabin rises while the flash is in transit, so anyone inside sees the beam bend downward. Light falls. Adding that small bending up along a ray skimming past the Sun gives a deflection of about nine tenths of an arcsecond.

The measured value is one and three quarter arcseconds, so the calculation is short by a factor of two, and the factor is exact rather than approximate. Nothing has gone wrong arithmetically. What has gone wrong is that every step of the cabin argument concerned durations: how long the crossing takes, how much speed the cabin picks up meanwhile. It therefore uses only the part of the eventual geometry that governs clocks, and no information whatever about how distances are measured near a heavy body.

Light spends its budget evenly between space and time, since it moves at the limiting speed, so the distance part contributes exactly as much as the time part and doubles the answer. A planet crawling along at a millionth of that speed barely samples the distance part, which is why the old theory works so well for everything except light.

8 · Why gravity, and only gravity, can be geometry

Everything is now in place for the argument the chapter exists to make. It has three steps.

Step 1. The paths belong to the arena. Section 2 established that through each event, and for each initial velocity, there passes exactly one free-fall trajectory, and that trajectory is the same for every body. So the assignment

(event, initial four-velocity)    a unique curve \big(\text{event},\ \text{initial four-velocity}\big) \;\longmapsto\; \text{a unique curve} (3.1.48)

is a well-defined rule that mentions nothing about matter. It is therefore permissible to regard it as a property of spacetime rather than of the objects moving through it. Permissible, note, and not yet compulsory.

Step 2. That is exactly what a geometry supplies. Chapter 1.2's second worked example took a rule for the length of a curve, extremised it, and produced from that rule alone a family of "straightest available" curves: straight lines in the plane, great circles on a sphere. The input was the metric and nothing else, with no reference to what was travelling. The output was precisely an assignment of the form (3.1.48), one curve through each point in each direction. So a geometry generates the right kind of object, and it does so for free.

Chapter 2.2 §7 already named which extremum it will have to be, and the answer is the strange one. Because of the minus sign in the interval, the straight worldline between two events carries more proper time than any bent one. So the curve a free body follows is the one of maximal elapsed proper time rather than minimal length.

That chapter's closing remark is the statement we are heading for: that a satellite's orbit is the path maximising the proper time of the clock aboard it. Chapter 3.3 §8 turns it into an equation by feeding dτ\int\dd\tau into the Euler–Lagrange machinery of Chapter 1.2, with dτ\dd\tau now computed from a position-dependent metric. Not a line of that machinery changes. Only dτ\dd\tau does.

Step 3. The geometry must be on spacetime, not on space. The curve in (3.1.48) depends on the initial velocity. A ball thrown gently and a bullet fired hard from the same window follow wildly different arcs through space. So the data selecting a curve is a point plus a direction in four dimensions, not in three, and the arena carrying the geometry has to be spacetime.

This is where §6 rejoins the argument. What the redshift showed is that the rule for the length of a worldline, the proper time of Chapter 2.3 §5, already varies from place to place. A position-dependent rule for the length of a worldline is exactly what Chapter 3.3 will call a metric field.

8.1 · Why the same treatment fails for every other force

The contrast is quantitative, and it is worth doing with numbers rather than in words. Consider four particles released from rest at the same point in the same uniform electric field EE. By (3.1.5) their accelerations are in the ratio of their charge-to-mass ratios:

Particleq/mq/m (C kg⁻¹)Acceleration relative to the proton
Electron1.7588×1011-1.7588\times10^{11}1836-1836
Proton+9.5788×107+9.5788\times10^{7}+1+1
Helium nucleus+4.8225×107+4.8225\times10^{7}+0.5035+0.5035
Neutron0000

Four particles, four different curves through the same point with the same initial velocity, one of them a straight line. No assignment of the form (3.1.48) exists here, because the map is not single-valued. You must know what is moving before you can say where it goes. There is nothing for a geometry to be.

Electromagnetism gets a different treatment. It becomes a connection on an internal space rather than on spacetime itself, which is Part VI, and the reason the treatments differ is the table above.

The same test disqualifies the strong and weak nuclear interactions for the same reason and disqualifies air resistance, which depends on cross-sectional area. Gravity is alone in the class, and it is alone in it because of §1, an experimental fact with error bars.

⚠ What would break this

The argument above is only as strong as the universality it rests on. Suppose a new interaction were discovered that coupled to, say, baryon number rather than to mass–energy. Its effects would masquerade as gravity but would differ slightly between materials with different neutron-to-proton ratios, η\eta would be nonzero, and step 1 would fail: there would no longer be a single family of free-fall curves for the geometry to encode. General relativity would then be an approximation valid to whatever precision η\eta has been bounded, which as of §1 is a few parts in 101510^{15}. This is why the Eötvös experiment is still being improved a century and a quarter later, and it is the honest statement of what Part III is standing on.

In plain terms 3.1.8

The case for treating gravity as the shape of the arena rather than as a force within it rests on a single experimental fact and collapses without it. Because everything falls the same way, exactly one path leaves each event in each direction, and that path can be described without naming what travels along it. A rule of that kind is what a geometry hands you free of charge, as the earlier chapter on stationary action showed when it produced straight lines on a flat sheet and great circles on a globe from nothing but a recipe for measuring length.

Electricity fails the same test, and loudly rather than marginally. Release an electron, a proton, a helium nucleus and a neutron together in one electric field and they trace four different curves, one of them straight. There is no single family of paths belonging to the region, so there is nothing for a geometry to be, and the electric force needs a quite different device that a later part supplies.

Two closing observations point forward. The paths depend on how fast you were going when you set off, so the geometry must live on space and time together rather than on space alone. And the rule for the length of a path was already shown to vary from place to place by the business with the clocks.

9 · Worked examples

Worked example 1 — where a falling body is torn apart

A person of height =1.8 m\ell=1.8\ \mathrm{m} falls feet-first toward a black hole of mass MM. Take "torn apart" to mean that the head-to-toe relative acceleration reaches 10g10g. At what distance does that happen, and how does that distance compare with the horizon radius rs=2GM/c2r_{s}=2GM/c^{2}? Do it for M=10MM=10M_{\odot} and for M=4×106MM=4\times10^{6}M_{\odot}, the mass of the hole at the centre of our galaxy.

The head-to-toe separation is radial, so (3.1.20) applies:

Δa  =  2GMr3  =  10gr  =  (2GM10g)1/3. \Delta a \;=\; \frac{2GM\,\ell}{r^{3}} \;=\; 10g \qquad\Longrightarrow\qquad r \;=\; \left(\frac{2GM\ell}{10g}\right)^{1/3}.

Ten solar masses. GM=10×1.3271×1020=1.3271×1021 m3s2GM = 10\times1.3271\times10^{20}=1.3271\times10^{21}\ \mathrm{m^{3}s^{-2}}, so

r=(2(1.3271×1021)(1.8)98.07)1/3=(4.872×1019)1/3=3.65×106 m3650 km, r = \left(\frac{2(1.3271\times10^{21})(1.8)}{98.07}\right)^{1/3} = \big(4.872\times10^{19}\big)^{1/3} = 3.65\times10^{6}\ \mathrm{m} \approx 3650\ \mathrm{km},

while rs=2GM/c2=2.95×104 m=29.5 kmr_{s}=2GM/c^{2}=2.95\times10^{4}\ \mathrm{m}=29.5\ \mathrm{km}. The tidal limit is 124 times further out than the horizon. You are destroyed long before you arrive.

Four million solar masses. GM=5.3085×1026GM=5.3085\times10^{26}, giving r=2.69×108 mr=2.69\times10^{8}\ \mathrm{m}, while rs=1.18×1010 mr_{s}=1.18\times10^{10}\ \mathrm{m}. Now the horizon is 44 times larger than the tidal limit, so an infalling observer crosses it feeling almost nothing. The head-to-toe stretch at the horizon is 2GM/rs3=1.16×103 ms22GM\ell/r_{s}^{3}=1.16\times10^{-3}\ \mathrm{m\,s^{-2}}, about 1.2×1041.2\times10^{-4} of gg. For comparison, the same calculation for the ten-solar-mass hole gives 1.9×107g1.9\times10^{7}g at its horizon.

Why the scaling is what it is. The tidal field goes as M/r3M/r^{3} and the horizon as MM, so the tidal field at the horizon goes as M/M3=M2M/M^{3}=M^{-2}. Bigger holes are gentler at the edge. That is a genuinely counter-intuitive consequence of a completely elementary scaling, and it is worth checking against §5. A large hole has small KK near its horizon, hence a large local inertial frame, hence nothing locally remarkable happens there. Chapter 3.8 §6 makes that precise by showing the horizon is a coordinate artefact and not a singularity.

Worked example 2 — why GPS would not work

A GPS satellite orbits at r=26560 kmr=26\,560\ \mathrm{km} from the Earth's centre, and the receiver sits on the ground at R=6371 kmR_{\oplus}=6371\ \mathrm{km}. Compute the rate difference between the two clocks, and the positioning error that would accumulate in one day if it were ignored. Use GM=3.9860×1014 m3s2GM_{\oplus}=3.9860\times10^{14}\ \mathrm{m^{3}s^{-2}}.

The gravitational part is (3.1.38) with Φ=GM/r\Phi=-GM/r:

ΔΦc2=GMc2(1R1r)=4.435×103(1.5696×1073.7651×108)=5.291×1010. \frac{\Delta\Phi}{c^{2}} = \frac{GM_{\oplus}}{c^{2}}\left(\frac{1}{R_{\oplus}}-\frac{1}{r}\right) = 4.435\times10^{-3}\big(1.5696\times10^{-7}-3.7651\times10^{-8}\big) = 5.291\times10^{-10}.

Positive, so the satellite clock runs fast. Over a day, 86400×5.291×1010=45.7 μs86400\times5.291\times10^{-10} = 45.7\ \mu\mathrm{s}.

The kinematic part is Chapter 2.3's time dilation, not this chapter's business, but the two must be combined. Circular orbit speed is v=GM/r=3874 ms1v=\sqrt{GM/r}=3874\ \mathrm{m\,s^{-1}}, so to first order the moving clock runs slow by v2/2c2=8.349×1011v^{2}/2c^{2}=8.349\times10^{-11}, which is 7.2 μs-7.2\ \mu\mathrm{s} per day.

Adding is legitimate here only because both effects are of order 101010^{-10} and their product is of order 102010^{-20}, far below anything measurable. That makes this a first-order combination. Chapter 3.8 §5 is the place where it is done properly, with a single metric rather than two separate small corrections.

net =45.77.2=38.5 μs per day. \text{net } = 45.7 - 7.2 = 38.5\ \mu\mathrm{s\ per\ day}.

The positioning error. GPS works by timing signals, so a clock error δt\delta t becomes a range error cδtc\,\delta t:

c×38.5 μs  =  (2.998×108)(3.85×105)  =  1.15×104 m    11.5 km per day. c\times 38.5\ \mu\mathrm{s} \;=\; (2.998\times10^{8})(3.85\times10^{-5}) \;=\; 1.15\times10^{4}\ \mathrm{m} \;\approx\; 11.5\ \mathrm{km\ per\ day}.

The satellites carry oscillators deliberately offset before launch to compensate. Two remarks are worth making. First, the gravitational term is six times the kinematic one, so this is predominantly a test of this chapter rather than of Part II. Second, we have neglected the receiver's motion with the Earth's rotation, worth v2/2c2=1.2×1012v^{2}/2c^{2}=1.2\times10^{-12} or 0.1 μs0.1\ \mu\mathrm{s} per day at the equator. That is a real correction in the operational system, and it is small enough to leave out of a calculation quoting three figures.

10 · Your turn

Problem 1 — the tide, from scratch

The Moon has Mm=7.342×1022 kgM_{\text{m}}=7.342\times10^{22}\ \mathrm{kg} and orbits at dm=3.844×108 md_{\text{m}}=3.844\times10^{8}\ \mathrm{m}. The Sun has GM=1.3271×1020 m3s2GM_{\odot}=1.3271\times10^{20}\ \mathrm{m^{3}s^{-2}} at d=1.496×1011 md_{\odot}=1.496\times10^{11}\ \mathrm{m}. Using (3.1.20), compute the tidal acceleration each raises across the Earth's radius, and their ratio. The Sun is 2.7×1072.7\times10^{7} times more massive than the Moon, so explain in one sentence how it nonetheless loses.

Solution

The relative acceleration between the Earth's centre and a point on its near surface, a radial separation =R\ell=R_{\oplus}, is 2GMR/d32GM\,R_{\oplus}/d^{3}.

Moon: 2(6.674×1011)(7.342×1022)(6.371×106)(3.844×108)3=1.099×106 ms2, \text{Moon: } \frac{2(6.674\times10^{-11})(7.342\times10^{22})(6.371\times10^{6})}{(3.844\times10^{8})^{3}} = 1.099\times10^{-6}\ \mathrm{m\,s^{-2}}, Sun: 2(1.3271×1020)(6.371×106)(1.496×1011)3=5.051×107 ms2. \text{Sun: } \frac{2(1.3271\times10^{20})(6.371\times10^{6})}{(1.496\times10^{11})^{3}} = 5.051\times10^{-7}\ \mathrm{m\,s^{-2}}.

The ratio is 2.182.18: the lunar tide is a little over twice the solar one, which is why spring and neap tides differ by roughly a factor of three in range rather than being equal or wholly lunar.

The one sentence. The force goes as M/d2M/d^{2} and the tide as M/d3M/d^{3}, so moving the Sun 389389 times further away than the Moon costs it an extra factor of 389389 beyond its mass advantage, and 2.7×107/3893=0.462.7\times10^{7}/389^{3}=0.46.

Note what this problem is really showing. The Sun's pull on the Earth is about 180180 times the Moon's, and is completely invisible in the tides, because a uniform pull is exactly what §2.1 deletes. Only the gradient survives, and gradients weigh distance more heavily.

Problem 2 — sizing your own inertial frame

An atom interferometer drops a cloud of atoms down a 10 m10\ \mathrm{m} tower at the Earth's surface and reads out phase differences corresponding to a position resolution of ε=1 nm\varepsilon=1\ \mathrm{nm} over a baseline of =1 cm\ell=1\ \mathrm{cm}. (a) For how long is the apparatus a local inertial frame at that precision? (b) The drop takes 2×10/g=1.43 s\sqrt{2\times10/g}=1.43\ \mathrm{s}. Is the tidal effect detectable? (c) The tidal field is not uniform over the tower either. Estimate the fractional change in KK between the top and bottom of a 10 m10\ \mathrm{m} drop.

Solution

(a) From (3.1.28) with K=1.541×106 s2K=1.541\times10^{-6}\ \mathrm{s^{-2}}, ε=109 m\varepsilon=10^{-9}\ \mathrm{m}, =102 m\ell=10^{-2}\ \mathrm{m}:

T<εK=109(1.541×106)(102)=6.49×102=0.255 s. T \lt \sqrt{\frac{\varepsilon}{K\ell}} = \sqrt{\frac{10^{-9}}{(1.541\times10^{-6})(10^{-2})}} = \sqrt{6.49\times10^{-2}} = 0.255\ \mathrm{s}.

(b) Very. The drop lasts 1.43 s1.43\ \mathrm{s}, which is 5.65.6 times longer, and since the drift goes as T2T^{2} the tidal displacement is (5.6)2=31(5.6)^{2}=31 times the resolution, or about 31 nm31\ \mathrm{nm}. Such instruments measure the tidal gradient deliberately. It is the signal, not the noise.

(c) K=GM/r3K=GM/r^{3}, so dK/K=3dr/r=3(10)/(6.371×106)=4.7×106\dd K/K = -3\,\dd r/r = -3(10)/(6.371\times10^{6}) = -4.7\times10^{-6}. Five parts in a million is negligible here, and it is a reminder that the hierarchy continues: the field has a gradient, the gradient has a gradient, and each is suppressed by another power of the size of the apparatus over the distance to the source.

Problem 3 — a clock on a mountain

(a) An optical clock is carried from sea level to the summit of a 4000 m4000\ \mathrm{m} mountain. By how much does it gain per year? (b) Modern optical clocks reach a fractional frequency stability of about 1×10181\times10^{-18}. What height difference does that correspond to at the Earth's surface? (c) Comment on what that means for the definition of "sea level".

Solution

(a) ΔΦ/c2=gh/c2=(9.807)(4000)/(8.988×1016)=4.365×1013\Delta\Phi/c^{2}=gh/c^{2}=(9.807)(4000)/(8.988\times10^{16})=4.365\times10^{-13}. Over a year of 3.156×107 s3.156\times10^{7}\ \mathrm{s}, the gain is 1.38×105 s1.38\times10^{-5}\ \mathrm{s}, or about 14 μs14\ \mu\mathrm{s} per year.

(b) Setting gh/c2=1018gh/c^{2}=10^{-18} gives h=(1018)(8.988×1016)/9.807=9.2×103 mh = (10^{-18})(8.988\times10^{16})/9.807 = 9.2\times10^{-3}\ \mathrm{m}, about 9 millimetres.

(c) A clock is now a better altimeter than an altimeter. "Sea level" is properly defined as a surface of constant gravitational potential, which is called the geoid, and (3.1.38) says two clocks tick at the same rate precisely when they sit on the same equipotential. So a network of optical clocks measures the geoid directly, to centimetre accuracy, without anyone having to survey anything. This is called relativistic geodesy, and it is a case of a supposedly exotic effect becoming an instrument.

Problem 4 — could you geometrise electromagnetism if you tried?

Suppose someone proposes that the electric field is also "really" a curvature of spacetime, so that charged particles follow paths determined by the geometry alone. (a) Show that the proposal is already refuted by two particles, and state the minimum experiment. (b) A defender replies: "restrict attention to a single species, say protons, and the paths are universal." Is the reply sound? (c) In a uniform field E=106 Vm1E=10^{6}\ \mathrm{V\,m^{-1}}, compute how far apart a proton and a deuteron released together are after 1 ns1\ \mathrm{ns}. Take md=2mpm_{d}=2m_{p} to two figures.

Solution

(a) Geometrisation requires the map (3.1.48) to be single-valued: one curve per event per initial four-velocity. Release a proton and a neutron from the same point with the same velocity in an electric field. The neutron goes straight and the proton curves. Two curves, same initial data, so the map is not a function of the initial data and there is nothing for a geometry to encode. The minimum experiment is exactly that pair.

(b) Partly sound and entirely useless. Within one species q/mq/m is fixed, so the paths are universal and one could absorb q/mq/m into a geometry. But it would be a different geometry for each species, and spacetime does not come in one copy per particle type. The point of geometrisation is that the arena is shared. This also identifies exactly what §1's experiment is testing: not that some ratio exists, but that one ratio serves for everything.

(c) ap=(q/mp)E=(9.579×107)(106)=9.579×1013 ms2a_{p}=(q/m_{p})E = (9.579\times10^{7})(10^{6}) = 9.579\times10^{13}\ \mathrm{m\,s^{-2}}, and ad=ap/2=4.789×1013a_{d}=a_{p}/2=4.789\times10^{13}. After t=109 st=10^{-9}\ \mathrm{s},

Δx=12(apad)t2=12(4.789×1013)(1018)=2.4×105 m, \Delta x = \tfrac12(a_{p}-a_{d})t^{2} = \tfrac12(4.789\times10^{13})(10^{-18}) = 2.4\times10^{-5}\ \mathrm{m},

That is about 24 μm24\ \mu\mathrm{m}, which is vast on the scale of any detector. Compare it with the corresponding gravitational separation, which by §1 is bounded by η3×1015\eta\lesssim3\times10^{-15} times the common displacement. The two cases are not close.

The brick you just laid

You now hold the argument that gravity is eligible to be geometry, and you hold it in the form of derivations rather than slogans. The equality of inertial and gravitational mass (⚑ measured, to three parts in 101510^{15}) makes a falling body's path independent of the body. That in turn makes a uniform field deletable by the coordinate change (3.1.6), and deletable for all bodies at once, which is the part no other force can match.

What survives the deletion is the tidal equation (3.1.14), d¨i=Φ,ijdj\ddot d^{\,i}=-\Phi_{,ij}d^{\,j}: stretch by 2GM/r32GM/r^{3} along the field, squeeze by GM/r3GM/r^{3} across it, traceless in vacuum, and therefore volume-preserving. That equation is the whole physical content of gravity, and no change of coordinates removes it.

Where this gets spent. Chapter 3.2 builds the arena that can carry a geometry without a background to draw it in. Chapter 3.3 supplies the metric field whose position dependence §6 already forced, and its geodesics are the curves of (3.1.48). Chapter 3.4 constructs the Riemann tensor and derives an equation that is (3.1.14) written covariantly, at which point tidal acceleration and curvature become the same thing. Chapter 3.6's field equations have (3.1.22) as their Newtonian shadow, and read off 8πG/c48\pi G/c^{4} by matching to (3.1.39). Chapter 3.8 pays the factor of two.

The shape of the argument, for the second time. Chapter 2.1 ended by predicting this chapter's structure and it was right. A principle held on excellent grounds met a fact held on equally good grounds. The principle was that the laws are the same for everyone, with a finite speed built into the geometry. The fact was that everything falls the same way.

The two are incompatible, since one forbids instantaneous influence and the other is stated in a theory built on it. As in Part II, the resolution is not a new force but the removal of a piece of assumed structure that nobody had noticed was an assumption. Last time it was absolute simultaneity. This time it is the flatness of spacetime.

What you should not yet believe. Nothing above shows that gravity is geometry. It shows only that it may be, and that nothing else may. The demonstration requires field equations, and those need five more chapters of machinery.