# Problema di Keplero nella relatività generale

Il problema di Keplero nella relatività generale comporta la risoluzione del moto di due corpi sferici che interagiscono tra di loro per mezzo della gravitazione, come descritto dalla teoria della relatività generale.
In genere, e in questo articolo, un corpo viene ipotizzato avente una massa $m$ trascurabile rispetto alla massa $M$ di un altro corpo; questa è una buona approssimazione nel caso di un pianeta che ruota attorno al Sole, o di un fotone che passa da una stella. In tali casi, si può supporre che solo il corpo più pesante contribuisce alla curvatura dello spazio-tempo e che sia fisso nello spazio. Questo spazio-tempo curvo è descritto dalla soluzione di Schwarzschild per le equazioni di Einstein nel vuoto della relatività generale. Il moto del corpo più leggero (successivamente sarà chiamato "particella") è definito per mezzo delle geodetiche dello spazio-tempo della soluzione di Schwarzschild. Viene qui assunto che il corpo più leggero sia simile a un punto in modo che le forze di marea possano essere ignorate.
Queste soluzioni geodetiche giustificano la precessione anomala del pianeta Mercurio, che è la prova principale a sostegno della teoria della relatività generale. Esse descrivono anche la deflessione della luce in un campo gravitazionale, un'altra previsione notoriamente utilizzata come prova per la relatività generale.
Il decadimento orbitale dovuto all'emissione di radiazione gravitazionale non viene definito dalla soluzione di Schwarzschild.

## Contesto storico e comprensione intuitiva

Nel 1859, Urbain Le Verrier scoprì che la precessione orbitale del pianeta Mercurio non era abbastanza rispetto a quella che avrebbe dovuto essere; l'ellisse della sua orbita ruotava (in precessione) leggermente più veloce di quanto previsto dalla teoria tradizionale della gravità newtoniana, anche dopo aver preso in considerazione tutti gli effetti degli altri pianeti. L'effetto è piccolo (grosso modo 43 arcosecondi di rotazione per secolo), ma ben al di sopra dell'errore di misurazione (pressappoco 0,1 arcsecondi per secolo). Le Verrier comprese l'importanza della sua scoperta immediatamente, sfidando gli astronomi e i fisici a metterla in dubbio. Sono state proposte diverse spiegazioni classiche, come la polvere interplanetaria, lo schiacciamento polare inosservato del Sole, una luna di Mercurio ancora non scoperta, o un nuovo pianeta chiamato Vulcano. Dopo che queste spiegazioni vennero ritenute di poco valore, alcuni fisici furono tentati dall'ipotesi più radicale, considerando sbagliata la legge del quadrato inverso della gravitazione di Newton. Per esempio, alcuni fisici proposero una legge di potenza con un esponente che fosse leggermente diverso da 2.
Altri ancora hanno sostenuto che la legge di Newton dovrebbe essere integrata con un potenziale dipendente dalla velocità. Tuttavia, questo implicava un conflitto con la dinamica celeste newtoniana. Nel suo trattato sulla meccanica celeste, Laplace aveva dimostrato che, se l'influenza gravitazionale non agisse istantaneamente, allora il moto dei pianeti stessi non conserverebbe esattamente la quantità di moto (parte del momento deve quindi essere attribuita al mediatore dell'Interazione gravitazionale, analogamente al momento che si attribuisce al mediatore della interazione elettromagnetica). Come si è osservato dal punto di vista newtoniano, se l'influenza gravitazionale si propaga a una velocità finita, allora in tutti i punti del tempo un pianeta viene attratto in un punto in cui il Sole era qualche tempo prima, e non verso la posizione istantanea del sole. In base all'ipotesi dei fondamenti classici, Laplace aveva dimostrato che se la gravità si propagasse a una velocità dell'ordine della velocità della luce allora il sistema solare sarebbe instabile, e non esistere per molto tempo. L'osservazione sul fatto che il sistema solare è vecchio permette di porre un limite inferiore alla velocità gravitazionale che è di molti ordini di grandezza più veloce della velocità della luce. La valutazione di Laplace per la velocità gravitazionale non è esatta, perché in una teoria di campo che rispetti il principio di relatività, l'attrazione di una carica puntiforme che si muove a una velocità costante è verso la posizione istantanea estrapolata, non verso la posizione apparente che sembra occupi quando si osserva.
Per evitare questi problemi, tra il 1870 e il 1900 molti scienziati utilizzavano le leggi dell'elettrodinamica di Wilhelm Eduard Weber, Carl Friedrich Gauß, Bernhard Riemann per produrre orbite stabili e spiegare lo spostamento del perielio dell'orbita di Mercurio. Nel 1890 Lévy riuscì a farlo combinando le leggi di Weber e di Riemann, per cui la velocità gravitazionale è nella sua teoria uguale a quella della luce. E in un altro tentativo Paul Gerber (1898) riuscì ugualmente a ottenere la formula corretta per lo spostamento del perielio (la quale era identica a quella usata più tardi da Einstein). Tuttavia, dato che le leggi fondamentali di Weber e di altri erano errate (per esempio, la legge di Weber venne rimpiazzata dalla teoria di Maxwell), quelle ipotesi vennero rifiutate. Un altro tentativo fatto da Hendrik Lorentz (1900), che aveva già utilizzato la teoria di Maxwell, produsse un spostamento del perielio troppo basso.
Intorno al 1904–1905, i lavori di Hendrik Lorentz, Henri Poincaré e infine la teoria ristretta della relatività di Albert Einstein, escludevano la possibilità di propagazione di qualsiasi effetto più veloce della velocità della luce. Ne consegue che la legge gravitazionale di Newton avrebbe dovuto essere sostituita con un'altra legge, compatibile con il principio di relatività, nonostante il limite newtoniano restasse valido nelle circostanze in cui gli effetti relativistici fossero trascurabili. Tali tentativi sono stati fatti da Henri Poincaré (1905), Hermann Minkowski (1907) e Arnold Sommerfeld (1910). Nel 1907 Einstein arrivò a concludere che per raggiungere questo obiettivo fosse necessaria un erede della relatività ristretta. Dal 1907 al 1915, Einstein tenta gradualmente di spingersi verso una nuova teoria, usando il suo principio di equivalenza come concetto chiave per raggiungere lo scopo. Secondo questo principio, un campo gravitazionale uniforme agisce allo stesso modo su tutto quanto vi è all'interno di esso e, di conseguenza, non può essere rilevato da un osservatore in caduta libera. Al contrario, tutti gli effetti gravitazionali locali dovrebbero essere riproducibili in un sistema di riferimento che acceleri in modo lineare, e viceversa. Così, la gravità agisce come una forza fittizia come la forza centrifuga o la forza di Coriolis, risultante da un sistema di riferimento accelerato; tutte le forze fittizie sono proporzionali alla massa inerziale, proprio come la gravità. Per conciliare la gravità con la relatività speciale incorporando il principio di equivalenza, qualcosa deve essere sacrificato; quel qualcosa fu il presupposto classico di lunga data che considerava il nostro spazio soggetto alle leggi della geometria euclidea, per es., il fatto che il teorema di Pitagora sia vero sperimentalmente. Einstein utilizzò una geometria più generale, la pseudo-geometria di Riemann, per consentire la curvatura di spazio e tempo, condizione necessaria per la riconciliazione; dopo otto anni di lavoro (1907–1915), riuscì a scoprire il modo preciso in cui spazio-tempo sarebbero curvi al fine di riprodurre le leggi fisiche osservate in natura, in particolare la gravitazione. La gravità è distinta dalle forze fittizie, centrifughe e di Coriolis, nel senso che la curvatura dello spazio-tempo è considerato come reale fisicamente, mentre le forze fittizie non sono considerate come forze. Le prime vere soluzioni delle sue equazioni di campo spiegavano la precessione anomala di Mercurio prevedendo un'insolita curvatura della luce, la quale venne confermata dopo che la sua teoria venne pubblicata. Queste soluzioni sono spiegate qui di seguito.

## Ambito geometrico

Nella normale geometria euclidea, i triangoli obbediscono al teorema di Pitagora, il quale stabilisce che la distanza al quadrato $ds^{2}$ tra due punti nello spazio è la somma dei quadrati delle sue componenti perpendicolari
    $ds^{2}=dx^{2}+dy^{2}+dz^{2}\,\!$
dove $dx$, $dy$ e $dz$ rappresentano le differenze infinitesimali tra i due punti lungo gli assi $x$, $y$ e $z$ di un sistema di coordinate cartesiane (aggiungere qui la figura). Ora immaginate un mondo in cui questo non è del tutto vero, un mondo dove la distanza è invece data da
    $ds^{2}=F(x,y,z)dx^{2}+G(x,y,z)dy^{2}+H(x,y,z)dz^{2}\,\!$
dove $F$, $G$ e $H$ sono funzioni arbitrarie della posizione. Non è difficile immaginare un siffatto mondo: noi ci viviamo. La superficie del mondo è curvo, e ciò rende vana la possibilità di farne una mappa perfettamente precisa e piana. I sistemi di coordinate non cartesiane descrivono bene questo, per esempio, nelle coordinate sferiche ($r$, $\theta$, $\phi$), la distanza euclidea può essere scritta
    $ds^{2}=dr^{2}+r^{2}d\theta ^{2}+r^{2}\sin ^{2}\theta d\varphi ^{2}\,\!$
Un altro esempio sarebbe un mondo in cui i regoli utilizzati per misurare la lunghezza sarebbero inaffidabili, mutando la loro lunghezza con la loro posizione e anche la loro orientazione. Nel caso più generale, si devono consentire i termini trasversali (cross-terms) quando si calcola la distanza $ds$:
    $ds^{2}=g_{xx}dx^{2}+g_{xy}dxdy+g_{xz}dxdz+\cdots +g_{zy}dzdy+g_{zz}dz^{2}\,\!$
dove nove funzioni $g_{xx}$, $g_{xy}$ costituiscono il tensore metrico, il quale definisce la geometria dello spazio nella geometria di Riemann. Nell'esempio delle coordinate sferiche di cui sopra, non vi sono termini trasversali; le sole componenti del tensore metrico diverse da zero sono $g_{rr}=1,g_{\theta \theta }=r^{2}$, e $g_{\phi \phi }=r^{2}sin^{2}{\theta }$.
Nella sua teoria ristretta della relatività, Albert Einstein dimostrava che la distanza $ds$ tra due punti spaziali non è costante, ma dipende dal moto dell'osservatore. Tuttavia, vi è una misura di separazione tra due punti nello spazio-tempo — chiamata "tempo proprio" e indicata con il simbolo $d\tau$ — che è invariante; in altre parole, non dipende dal moto dell'osservatore.
    $c^{2}d\tau ^{2}=c^{2}dt^{2}-dx^{2}-dy^{2}-dz^{2}\,\!$
la quale può essere scritta con le coordinate sferiche nel seguente modo:
    $c^{2}d\tau ^{2}=c^{2}dt^{2}-dr^{2}-r^{2}d\theta ^{2}-r^{2}\sin ^{2}\theta d\varphi ^{2}\,\!$
Questa formula è l'estensione naturale del teorema di Pitagora e analogamente regge solo quando non c'è curvatura dello spazio-tempo. Nella relatività generale, comunque, spazio e tempo possono avere curvatura, quindi questa formula della distanza deve essere modificata in una forma più generale:
    $c^{2}d\tau ^{2}=g_{\mu \nu }dx^{\mu }dx^{\nu }\,\!$
così come abbiamo generalizzato la formula per misurare la distanza sulla superficie della Terra. La forma esatta della metrica $g_{\mu \nu }$ dipende dalla massa gravitante, momento ed energia, come descritta dalle equazioni di campo di Einstein. Einstein sviluppò tali equazioni di campo per soddisfare le leggi della natura allora conosciute ma, ciò nonostante, esse hanno previsto fenomeni mai visti prima (come la curvatura della luce a causa della gravità), confermati in seguito.

## Metrica di Schwarzschild

Una soluzione per le equazioni di campo di Einstein è la metrica di Schwarzschild, corrispondente al campo gravitazionale di un corpo sfericamente simmetrico, non rotante, senza carica e di massa $M$. La soluzione di Schwarzschild può essere scritta come:
    $c^{2}{d\tau }^{2}=\left(1-{\frac {r_{s}}{r}}\right)c^{2}dt^{2}-{\frac {dr^{2}}{1-{\frac {r_{s}}{r}}}}-r^{2}d\theta ^{2}-r^{2}\sin ^{2}\theta \,d\varphi ^{2}$
dove:
    $\tau$ è il tempo proprio (il tempo misurato da un orologio che si muove insieme alla particella) in secondi,
    $c$ è la velocità della luce in metri al secondo,
    $t$ è la coordinata del tempo (misurata da un orologio stazionario all'infinito) in secondi,
    $r$ è la coordinata radiale (circonferenza di un cerchio centrato sulla stella diviso per $2\pi$) in metri,
    $\theta$ è la colatitudine (angolo da Nord) in radianti,
    $\varphi$ è la longitudine in radianti, e
    $r_{s}$ è il raggio di Schwarzschild (in metri) del corpo massivo, il quale è correlato alla sua massa $M$ da
       $r_{s}={\frac {2GM}{c^{2}}}$
    dove $G$ è la costante gravitazionale.
La teoria newtoniana classica della gravità è recuperata nel limite come il rapporto ${\frac {r_{s}}{r}}$ va a zero. In questo limite, la metrica ritorna nella forma data sopra per la relatività ristretta. In pratica, questo rapporto è quasi sempre estremamente piccolo. Per esempio, il raggio di Schwarzschild $r_{s}$ della Terra è grosso modo 9 mm, laddove un satellite in un'orbita geosincrona ha un raggio $r$ che è di circa quattro bilioni di volte più grande, a 42.164 km. Anche per la superficie della Terra, le correzioni per la gravità newtoniana sono soltanto una parte su un bilione. Il rapporto diventa grande solo vicino ai buchi neri e agli oggetti ultra-densi come le stelle di neutroni.

## Equazione geodetica

Secondo la teoria di Einstein della relatività generale, le particelle di massa trascurabile viaggiano lungo le geodetiche nello spazio-tempo. In uno spazio-tempo non curvo, lontano dalla sorgente di gravità, queste geodetiche corrispondono a linee rette; tuttavia, possono deviare dalla linee rette quando lo spazio-tempo è curvo. L'equazione per le linee geodetiche è
    ${\frac {d^{2}x^{\mu }}{dq^{2}}}+\Gamma _{\nu \lambda }^{\mu }{\frac {dx^{\nu }}{dq}}{\frac {dx^{\lambda }}{dq}}=0$
dove $\Gamma$ rappresenta il simbolo di Christoffel e la variabile $q$ parametrizza il percorso della particella attraverso lo spazio-tempo, la sua cosiddetta linea di universo. Il simbolo di Christoffel dipende soltanto dal tensore metrico $g_{\mu \nu }$, o piuttosto da come esso muta con la posizione. La variabile $q$ è una costante multipla del tempo proprio $\tau$ per orbite di tipo tempo (che sono percorse da particelle massive), ed è di solito considerata uguale a essa. Per le orbite di tipo luce (o nulle) (le quali sono percorse da particelle senza massa come i fotoni), il tempo proprio è zero e, a rigore, non può essere usato come la variabile $q$. Ciò nonostante, le orbite di tipo luce possono essere derivate come il limite ultrarelativistico delle orbite di tipo tempo, vale a dire, il limite come la massa $m$ della particella va a zero mentre mantiene fissa la sua energia complessiva.
Possiamo semplificare il problema utilizzando la simmetria per eliminare una variabile presa in considerazione. Poiché la metrica di Schwarzschild è simmetrica di circa $\theta ={\frac {\pi }{2}}$, ogni geodetica che inizia a muoversi in quel piano rimarrà in quel piano indefinitamente (il piano è totalmente geodetico). Quindi, orientiamo il sistema di coordinate in modo che l'orbita della particella venga a trovarsi in quel piano, e fissiamo la coordinata $\theta$ come ${\frac {\pi }{2}}$ in modo che la metrica (di questo piano) si semplifichi in:
    $c^{2}d\tau ^{2}=\left(1-{\frac {r_{s}}{r}}\right)c^{2}dt^{2}-{\frac {dr^{2}}{1-{\frac {r_{s}}{r}}}}-r^{2}d\varphi ^{2}.$
Due costanti di moto possono essere dunque identificate (cfr. la derivazione data nell'approccio lagrangiano)
    $r^{2}{\frac {d\varphi }{d\tau }}={\frac {L}{m}},$
    $\left(1-{\frac {r_{s}}{r}}\right){\frac {dt}{d\tau }}={\frac {E}{mc^{2}}}.$
Sostituendo queste costanti dentro la definizione della metrica di Schwarzschild
    $c^{2}=\left(1-{\frac {r_{s}}{r}}\right)c^{2}\left({\frac {dt}{d\tau }}\right)^{2}-{\frac {1}{1-{\frac {r_{s}}{r}}}}\left({\frac {dr}{d\tau }}\right)^{2}-r^{2}\left({\frac {d\varphi }{d\tau }}\right)^{2},$
si ottiene l'equazione di moto per la particella:
    $\left({\frac {dr}{d\tau }}\right)^{2}={\frac {E^{2}}{m^{2}c^{2}}}-\left(1-{\frac {r_{s}}{r}}\right)\left(c^{2}+{\frac {L^{2}}{m^{2}r^{2}}}\right).$
La dipendenza dal tempo proprio può essere eliminata usando la definizione di $L$
    $\left({\frac {dr}{d\varphi }}\right)^{2}=\left({\frac {dr}{d\tau }}\right)^{2}\left({\frac {d\tau }{d\varphi }}\right)^{2}=\left({\frac {dr}{d\tau }}\right)^{2}\left({\frac {mr^{2}}{L}}\right)^{2},$
che fornisce l'equazione per l'orbita:
    $\left({\frac {dr}{d\varphi }}\right)^{2}={\frac {r^{4}}{b^{2}}}-\left(1-{\frac {r_{s}}{r}}\right)\left({\frac {r^{4}}{a^{2}}}+r^{2}\right)$
dove, per brevità, le due scale di lunghezza, $a$ e $b$, vengono definite da:
    $a={\frac {L}{mc}},$
    $b={\frac {cL}{E}}.$
La stessa equazione può anche essere derivata usando un approccio lagrangiano o l'equazione di Hamilton – Jacobi (vedi sotto). La soluzione dell'equazione orbitale è:
    $\varphi =\int {\frac {dr}{r^{2}{\sqrt {{\frac {1}{b^{2}}}-\left(1-{\frac {r_{s}}{r}}\right)\left({\frac {1}{a^{2}}}+{\frac {1}{r^{2}}}\right)}}}}.$

## Formula approssimata per la curvatura della luce

Nel limite come la massa $m$ della particella va a zero (o, in modo equivalente, come la scala di lunghezza $a$ va all'infinito), l'equazione per l'orbita diventa:
    $\varphi =\int {\frac {dr}{r^{2}{\sqrt {{\frac {1}{b^{2}}}-\left(1-{\frac {r_{s}}{r}}\right){\frac {1}{r^{2}}}}}}}$
Espandendo in potenze di ${\frac {r_{s}}{r}}$, il termine di ordine dominante in questa formula dà la deflessione angolare approssimata $d\varphi$ per una particella senza massa che viene dall'infinito e ritorna all'infinito:
    $\delta \varphi \approx {\frac {2r_{s}}{b}}={\frac {4GM}{c^{2}b}}.$
Qui, $b$ può essere interpretata come la distanza di un approccio più stretto. Sebbene questa formula sia approssimata, essa è precisa per la maggior parte delle misurazioni del fenomeno delle lenti gravitazionali, dovuto alla piccolezza del rapporto ${\frac {r_{s}}{r}}$. Per la luce radente la superficie del sole, la deflessione angolare approssimativa è di circa 1,75 arcosecondi.

## Relazione della meccanica classica e la precessione delle orbite ellittiche

L'equazione di moto per la particella derivata sopra
    $\left({\frac {dr}{d\tau }}\right)^{2}={\frac {E^{2}}{m^{2}c^{2}}}-c^{2}+{\frac {r_{s}c^{2}}{r}}-{\frac {L^{2}}{m^{2}r^{2}}}+{\frac {r_{s}L^{2}}{m^{2}r^{3}}}$
può essere riscritta usando la definizione del raggio di Schwarzschild $r_{s}$ come:
    ${\frac {1}{2}}m\left({\frac {dr}{d\tau }}\right)^{2}=\left[{\frac {E^{2}}{2mc^{2}}}-{\frac {1}{2}}mc^{2}\right]+{\frac {GMm}{r}}-{\frac {L^{2}}{2mr^{2}}}+{\frac {GML^{2}}{c^{2}mr^{3}}}$
che è equivalente a una particella che si muove in un potenziale effettivo unidimensionale
    $V(r)=-{\frac {GMm}{r}}+{\frac {L^{2}}{2mr^{2}}}-{\frac {GML^{2}}{c^{2}mr^{3}}}.$
I primi termini sono energie classiche ben note, essendo la prima l'energia potenziale gravitazionale attrattiva newtoniana e la seconda corrisponde all'energia potenziale "centrifuga" repulsiva; Tuttavia, il terzo termine è un è un'energia attrattiva unica per la relatività generale. Come mostrato sotto e altrove, questa energia cubica-inversa causa gradualmente la precessione delle orbite ellittiche di un angolo $\delta \varphi$ per la rivoluzione:
    $\delta \varphi \approx {\frac {6\pi GM}{c^{2}A\left(1-e^{2}\right)}}$
dove $A$ è il semi-asse maggiore ed $e$ è l'eccentricità.
Il terzo termine è attrattivo e domina in valori di $r$ piccoli, dando un raggio critico interno $r_{interno}$ dove una particella è attratta inesorabilmente verso l'interno per $r=0$; questo raggio interno è una funzione del momento angolare della particella per unità di massa o, in modo equivalente, la scala di lunghezza $a$ definita sopra.

## Orbite circolari e loro stabilità

Il potenziale effettivo $V$ può essere riscritto nei termini delle lunghezze $a$ e $b$
    $V(r)={\frac {mc^{2}}{2}}\left[-{\frac {r_{s}}{r}}+{\frac {a^{2}}{r^{2}}}-{\frac {r_{s}a^{2}}{r^{3}}}\right]$
Le orbite circolari sono possibili quando la forza effettiva è zero:
    $F=-{\frac {dV}{dr}}=-{\frac {mc^{2}}{2r^{4}}}\left[r_{s}r^{2}-2a^{2}r+3r_{s}a^{2}\right]=0$
i.e., quando due forze attrattive — la gravità newtoniana (primo termine) e l'attrazione unica per la relatività generale (terzo termine) — sono esattamente bilanciati dalla forza centrifuga repulsiva (secondo termine). Ci sono due raggi dove questo bilanciamento può accadere, definito qui come $r_{inner}$ e $r_{outer}$:
    $r_{\mathrm {outer} }={\frac {a^{2}}{r_{s}}}\left(1+{\sqrt {1-{\frac {3r_{s}^{2}}{a^{2}}}}}\right)$
    $r_{\mathrm {inner} }={\frac {a^{2}}{r_{s}}}\left(1-{\sqrt {1-{\frac {3r_{s}^{2}}{a^{2}}}}}\right)={\frac {3a^{2}}{r_{\mathrm {outer} }}}$
che sono ottenute usando la formula quadratica. Il raggio interno $r_{inner}$ è instabile, poiché la terza forza attrattiva aumenta molto più velocemente delle altre due forze quando $r$ diventa piccolo; se la particella scivola leggermente verso l'interno da $r_{inner}$ (dove tutte e tre le forze sono in equilibrio), la terza forza domina le altre due, trascinando la particella inesorabilmente verso l'interno per $r=0$. Nel raggio esterno, comunque, le orbite circolari sono stabili; il terzo termine è meno importante e il sistema si comporta in modo più simile al problema di Keplero non-relativistico.
Quando $a$ è molto più grande di $r_{s}$ (caso classico), queste formule diventano approssimativamente:
    $r_{\mathrm {outer} }\approx {\frac {2a^{2}}{r_{s}}}$
    $r_{\mathrm {inner} }\approx {\frac {3}{2}}r_{s}$
Sostituendo le definizioni di $a$ e $r_{s}$ in $r_{outer}$ si ottiene la formula classica per una particella orbitante una massa $M$ in un cerchio:
    $r_{\mathrm {outer} }^{3}\approx {\frac {GM}{\omega _{\varphi }^{2}}}$
dove $\omega _{\varphi }$ è la velocità angolare orbitale della particella. Questa formula è ottenuta nella meccanica non-relativistica impostando la forza centrifuga uguale alla forza gravitazionale newtoniana:
    ${\frac {GMm}{r^{2}}}=\mu \omega r$
dove $\mu$ è la massa ridotta. Nella nostra notazione, la velocità orbitale classica è uguale a:
    $\omega _{\varphi }^{2}\approx {\frac {GM}{r_{\mathrm {outer} }^{3}}}=\left({\frac {r_{s}c^{2}}{2r_{\mathrm {outer} }^{3}}}\right)=\left({\frac {r_{s}c^{2}}{2}}\right)\left({\frac {r_{s}^{3}}{8a^{6}}}\right)={\frac {c^{2}r_{s}^{4}}{16a^{6}}}$
Nell'altro estremo, quando $a^{2}$ avvicina soprattutto a $3r_{s}^{2}$, i due raggi convergono in un unico valore
    $r_{\mathrm {outer} }\approx r_{\mathrm {inner} }\approx 3r_{s}$
Le soluzioni quadratiche sopra assicurano che $r_{outer}$ è sempre maggiore di $3r_{s}$, laddove $r_{inner}$ si trova fra ${\frac {3}{2}}r_{s}$ e $3r_{s}$. Le orbite circolari più piccole di ${\frac {3}{2}}r_{s}$ non sono possibili. Per le particelle senza massa, $a$ va all'infinito, il che implica che vi sia un'orbita circolare per i fotoni per $r_{inner}={\frac {3}{2}}r_{s}$. La sfera di questo raggio è talvolta nota come sfera fotonica.

## Precessione delle orbite ellittiche

La velocità (rate) di precessione orbitale può essere ricavata utilizzando questo potenziale effettivo radiale $V$. Una piccola deviazione radiale da un'orbita circolare di raggio $r_{outer}$ oscillerà stabilmente con una frequenza angolare
    $\omega _{r}^{2}={\frac {1}{m}}\left[{\frac {d^{2}V}{dr^{2}}}\right]_{r=r_{\mathrm {outer} }}$
che è pari a:
    $\omega _{r}^{2}=\left({\frac {c^{2}r_{s}}{2r_{\mathrm {outer} }^{4}}}\right)\left(r_{\mathrm {outer} }-r_{\mathrm {inner} }\right)=\omega _{\varphi }^{2}{\sqrt {1-{\frac {3r_{s}^{2}}{a^{2}}}}}$
Prendendo la radice quadrata di entrambi i lati e espandendo usando il teorema binomiale si ottiene la formula:
    $\omega _{r}=\omega _{\varphi }\left(1-{\frac {3r_{s}^{2}}{4a^{2}}}+\cdots \right)$
Moltiplicando per il periodo $T$ di una rivoluzione si ottiene la precessione dell'orbita per rivoluzione:
    $\delta \varphi =T\left(\omega _{\varphi }-\omega _{r}\right)\approx 2\pi \left({\frac {3r_{s}^{2}}{4a^{2}}}\right)={\frac {3\pi m^{2}c^{2}}{2L^{2}}}r_{s}^{2}$
dove abbiamo usato $\omega _{\varphi }T=2\pi$ e la definizione della scala di lunghezza $a$. Sostituendo la definizione del raggio di Schwarzschild $r_{s}$ si ha:
    $\delta \varphi \approx {\frac {3\pi m^{2}c^{2}}{2L^{2}}}\left({\frac {4G^{2}M^{2}}{c^{4}}}\right)={\frac {6\pi G^{2}M^{2}m^{2}}{c^{2}L^{2}}}$
Questa può essere semplificata usando il semiasse dell'orbita ellittica $A$ e l'eccentricità $e$ correlata tramite la formula:
    ${\frac {L^{2}}{GMm^{2}}}=A\left(1-e^{2}\right)$
per dare la forma più comune dell'angolo di precessione
    $\delta \varphi \approx {\frac {6\pi GM}{c^{2}A\left(1-e^{2}\right)}}$

## Soluzione orbitale usando funzioni ellittiche

L'equazione per l'orbita
    $\left({\frac {dr}{d\varphi }}\right)^{2}={\frac {r^{4}}{b^{2}}}-\left(1-{\frac {r_{s}}{r}}\right)\left({\frac {r^{4}}{a^{2}}}+r^{2}\right)$
può essere semplificata introducendo una variabile senza dimensione:
    $\zeta ={\frac {r_{s}}{4r}}-{\frac {1}{12}}$
in modo che si riduca a:
    $\left({\frac {d\zeta }{d\varphi }}\right)^{2}=4\zeta ^{3}-g_{2}\zeta -g_{3},$
dove la costante, i coefficienti senza dimensione $g_{2}$ e $g_{3}$ sono definiti tramite
    ${\begin{aligned}g_{2}&={\frac {1}{12}}-{\frac {r_{s}^{2}}{4a^{2}}},\\g_{3}&={\frac {1}{216}}+{\frac {r_{s}^{2}}{24a^{2}}}-{\frac {r_{s}^{2}}{16b^{2}}}.\end{aligned}}$
La soluzione di questa equazione orbitale è data da:
    $\varphi -\varphi _{0}=\int {\frac {d\zeta }{\sqrt {4\zeta ^{3}-g_{2}\zeta -g_{3}}}}.$
Ne consegue che, fino ad un cambiamento di fase, $\zeta =\wp (\varphi -\varphi _{0})$, dove $\wp$ è la funzione ellittica di Weierstrass con i parametri $g_{2}$ e $g_{3}$, e $\varphi _{0}$ è una costante di integrazione (possibilmente complessa).

## Carattere qualitativo di possibili orbite

Le orbite sono definite dall'equazione di moto:
    $\left({\frac {d\zeta }{d\varphi }}\right)^{2}=4\zeta ^{3}-g_{2}\zeta -g_{3},$
Se il discriminante $\Delta =g_{2}^{3}-27g_{3}^{2}$ è maggiore di zero, l'equazione cubica
    $G(\zeta )=4\zeta ^{3}-g_{2}\zeta -g_{3}=0\,$
avrà tre radici reali e distinte, $e_{1},\;e_{2}$, ed $e_{3}$, che possono essere elencate in ordine decrescente:
    $e_{1}>e_{2}>e_{3}.$
In tali casi, la soluzione $\zeta =\wp (\varphi -\varphi _{0})$ è una funzione ellittica con due semi-periodi, uno completamente reale
    $\omega _{1}=\int _{e_{1}}^{\infty }{\frac {dz}{\sqrt {4z^{3}-g_{2}z-g_{3}}}}$
e un altro completamente immaginario
    $\omega _{3}=i\int _{-e_{3}}^{\infty }{\frac {dz}{\sqrt {4z^{3}-g_{2}z-g_{3}}}}$.
La radice restante definisce un semi-periodo complesso $\omega _{2}=-\omega _{1}-\omega _{3}$. Questi tre semi-periodi sono correlati alle radici corrispondenti per mezzo dell'equazione $\wp (\omega _{i})=e_{i}$, dove $i$ può essere uguale a $1,\;2$ o $3$. Quindi, $\varphi _{0}$ è posto uguale a ognuno di questi semi-periodi, la derivata di $\zeta$ è zero, indicando un periapside o un apoapside, i punti di approccio rispettivamente più vicini o più lontani.
    ${\frac {d\zeta }{d\phi }}=0\ \mathrm {quando} \ \zeta =\wp (-\omega _{i})=e_{i}$,
dato che
    $\left({\frac {d\zeta }{d\varphi }}\right)^{2}=G(\zeta )=4\zeta ^{3}-g_{2}\zeta -g_{3}=4\left(\zeta -e_{1}\right)\left(\zeta -e_{2}\right)\left(\zeta -e_{3}\right)$.
Le radici $e_{i}$ sono i valori critici di $\zeta$ per quanto riguarda $\varphi$, i.e., i valori dove la derivata è uguale a zero.
Il carattere qualitativo delle orbite dipende dalla scelta di $\varphi _{0}$. Le soluzioni con $\varphi _{0}$ pari a $\omega _{2}$ corrispondono a orbite oscillatorie che variano tra $\zeta =e_{2}$ e $\zeta =e_{3}$, oppure divergono all'infinito $\left(\zeta =-{\frac {1}{12}}\right)$. Diversamente, le soluzioni con $\varphi _{0}$ pari a $\omega _{1}$ (o a ogni altro numero reale) corrispondono a orbite che decadono a raggio zero, dato che (per essere un numero reale) $\zeta$ non può essere inferiore a $e_{1}$ e perciò aumenta inesorabilmente all'infinito.

### Orbite quasi ellittiche

Le soluzioni $\zeta =\wp (\phi -\phi _{0})$ in cui $\varphi _{0}$ pari a $\omega _{2}$ danno un valore reale di $\zeta$ a condizione che l'energia $E$ soddisfi la disuguaglianza $E^{2}<m^{2}c^{4}$. Per tali soluzioni, la variabile $\zeta$ è confinata tra $e_{3}\leq \zeta \leq e_{2}$. Se entrambe le radici sono maggiori di $-{\frac {1}{12}}$, allora $\zeta$ non diventerà mai $-{\frac {1}{12}}$, il punto in cui il raggio $r$ va all'infinito. Perciò, tali soluzioni corrispondono a un'orbita ellittica in graduale precessione. Dato che la particella (o pianeta) ruota intorno alla origine, il suo raggio oscilla tra un raggio minimo
    $r_{min}={\frac {3r_{s}}{1+12e_{2}}}$
e un raggio massimo
    $r_{max}={\frac {3r_{s}}{1+12e_{1}}}$.
che corrisponde al valore di $\zeta$ ai due estremi (extrema) del raggio. La periodicità reale della funzione ellittica di Weierstrass è $2\omega _{1}$; così la particella ritorna allo stesso raggio dopo aver girato con un angolo di $2\omega _{1}$, che può non essere uguale a $2\pi$. Perciò, l'orbita è in precessione. In generale, la quantità di precessione per orbita ($2\omega _{1}-2\pi$) è abbastanza piccola.

#### Orbite circolari stabili

Un caso speciale si verifica quando $2e_{2}=2e_{3}=-e_{3}$, i.e., due delle radici di $G(\zeta )$ sono uguali e negative, mentre la terza è positiva. In questo caso c'è una soluzione con la costante $\zeta$, pari alla radice ripetuta, $e=e_{2}=e_{3}$. Ciò corrisponde a orbite circolari, in particolare, alla soluzione classica $r_{outer}$ derivata sopra; come qui dimostrato, i raggi di queste orbite devono essere maggiori di $3r_{s}$. Tali orbite circolari sono stabili, perché una piccola perturbazione dei parametri separerebbe le radici ripetute, risultanti in un'orbita quasi-ellittica. Ad esempio, dando un piccolo "colpetto" radiale a una particella in una classica orbita circolare la si spinge in un'orbita ellittica che gradualmente va in precessione.

#### Orbite illimitate

Un'orbita illimitata si verifica quando va all'infinito, corrispondente a $\zeta =-{\frac {1}{12}}$. Le orbite illimitate corrispondono a un'orbita oscillatoria in cui $-{\frac {1}{12}}$ cade tra le due radici limitanti, i.e., quando $e_{3}\leq -{\frac {1}{12}}\leq \zeta \leq e_{2}$.

#### Orbite circolari asintotiche

Un altro caso particolare succede quando $-e_{3}=2e_{2}=2e_{1}$, i.e., due delle radici di $G(\zeta )$ sono uguali e positive, mentre la terza radice $e_{3}$ è negativa. Indicando la radice ripetuta con $e={\frac {n^{2}}{3}}$, le orbite sono asintoticamente circolari all'infinito $\varphi$ positivo e negativo:
    $\zeta ={\frac {r_{s}}{4r}}-{\frac {1}{12}}=e-{\frac {n^{2}}{\cosh ^{2}n\varphi }}.$
come può essere verificato per mezzo della sostituzione. Come $\varphi$ va all'infinito positivo o negativo, l'orbita si avvicina asintoticamente al cerchio:
    ${\frac {r_{s}}{4r}}-{\frac {1}{12}}=e.$
In tali casi, il raggio dell'orbita deve rimanere tra $2r_{s}$ e $3r_{s}$.
La formula asintotica può anche essere derivata dall'espressione per la funzione ellittica di Weierstrass in termini di funzioni ellittiche di Jacobi
    $\zeta =\wp (\phi -\phi _{0})=e_{1}+\left(e_{1}-e_{3}\right){\frac {\mathrm {cn} ^{2}w}{\mathrm {sn} ^{2}w}}$
dove $w=(\phi -\phi _{0}){\sqrt {e_{1}-e_{3}}}$ e il modulo pari a:
    $k={\sqrt {\frac {e_{2}-e_{3}}{e_{1}-e_{3}}}}$
Nel limite, come $e_{2}$ si avvicina a $e_{1}$, il modulo tende a $1$ e $w$ va a $n(\varphi -\varphi _{0})$. Infine scegliendo $\varphi _{0}$ come numero immaginario $iK^{\prime }$ (un quarto periodo) dà la formula asintotica sopra.

### Orbite decadenti

Le soluzioni reali per $\zeta =\wp (\phi -\phi _{0})$ in cui $\varphi _{0}$ è pari a $\omega _{1}$ (o a qualche altro numero reale) hanno la proprietà per cui $\zeta$ non è mai minore di $e_{1}$. Dato che l'equazione di moto
    $\left({\frac {d\zeta }{d\varphi }}\right)^{2}=4\zeta ^{3}-g_{2}\zeta -g_{3}=4\left(\zeta -e_{1}\right)\left(\zeta -e_{2}\right)\left(\zeta -e_{3}\right)$
è positiva per tutti i valori di $\zeta >e_{1}$, $\zeta$ aumenta illimitatamente, corrispondente alla particella che cade inesorabilmente per l'origine $r=0$.

## Correzioni alle soluzioni geodetiche

Secondo la relatività generale, due corpi girando uno sull'altro emetteranno radiazione gravitazionale, facendo sì che le orbite siano leggermente diverse rispetto alle geodetiche calcolate sopra. Questo è stato osservato indirettamente in un sistema di stella binaria conosciuto come PSR B1913+16, per cui Russell Alan Hulse e Joseph Hooton Taylor Jr. vennero insigniti nel 1993 del Premio Nobel per la fisica. Le due stelle di neutroni di questo sistema sono estremamente vicine e ruotano quasi l'un l'altra molto velocemente, completando una rivoluzione in più o meno 465 minuti. La loro orbita è fortemente ellittica, con un'eccentricità di 0,62 (62%). Secondo la relatività generale, il breve periodo orbitale e l'alta eccentricità farebbero del sistema un eccellente emettitore di radiazione gravitazionale, perdendo così energia e riducendo ulteriormente il periodo orbitale. La diminuzione osservata nel periodo orbitale per più di trent'anni combacia con le previsioni della relatività generale anche nell'ambito delle più precise misurazioni. La relatività generale prevede che, nel giro di altri 300 milioni di anni, questi due astri si metteranno a spirale l'uno nell'altro.
Sono state calcolate le formule che definiscono la perdita di energia e il momento angolare dovuto alla radiazione gravitazionale dei due corpi del problema di Kepler. Il tasso di energia persa (in media su un'orbita completa) è dato da
    $-\langle {\frac {dE}{dt}}\rangle ={\frac {32G^{4}m_{1}^{2}m_{2}^{2}\left(m_{1}+m_{2}\right)}{5c^{5}a^{5}\left(1-e^{2}\right)^{7/2}}}\left(1+{\frac {73}{24}}e^{2}+{\frac {37}{96}}e^{4}\right)$
dove $e$ è l'eccentricità orbitale e $a$ è il semiasse maggiore dell'orbita ellittica. Le parentesi angolari sul lato sinistro dell'equazione rappresentano la media su una singola orbita. Allo stesso modo, il tasso medio di perdita del momento angolare è uguale a:
    $-\langle {\frac {dL_{z}}{dt}}\rangle ={\frac {32G^{7/2}m_{1}^{2}m_{2}^{2}{\sqrt {m_{1}+m_{2}}}}{5c^{5}a^{7/2}\left(1-e^{2}\right)^{2}}}\left(1+{\frac {7}{8}}e^{2}\right)$
Le perdite in energia e momento angolare aumentano significativamente come l'eccentricità si avvicina a $1$, i.e., come l'ellisse dell'orbita diventa sempre più allungata. Anche le perdite di radiazione aumentano significativamente con il diminuire della dimensione $a$ dell'orbita.

## Derivazioni matematiche dell'equazione orbitale
### Approccio di Hamilton–Jacobi

L'equazione orbitale può essere derivata dall'equazione di Hamilton-Jacobi. Il vantaggio di questo approccio è che si mette sullo stesso piano il moto della particella con la propagazione di un'onda, e porta ordinatamente dentro la derivazione della deflessione della luce per mezzo della gravità nella relatività generale, attraverso il principio di Fermat. L'idea di base è che, a causa del rallentamento gravitazionale del tempo, le parti di un fronte d'onda più vicino a una massa gravitante si muovono più lentamente rispetto a quelle più lontane, curvando così la direzione di propagazione del fronte dell'onda (aggiungi figura).
Usando la covarianza generica, l'equazione di Hamilton-Jacobi per una singola particella di coordinate arbitrarie può essere espressa come:
    $g^{\mu \nu }{\frac {\partial {\mathcal {S}}}{\partial x^{\mu }}}{\frac {\partial {\mathcal {S}}}{\partial x^{\nu }}}=m^{2}c^{2}.$
Usando la metrica di Schwarzschild $g^{\mu \nu }$, questa equazione diventa:
    ${\frac {1}{c^{2}\left(1-{\frac {r_{s}}{r}}\right)}}\left({\frac {\partial {\mathcal {S}}}{\partial t}}\right)^{2}-\left(1-{\frac {r_{s}}{r}}\right)\left({\frac {\partial {\mathcal {S}}}{\partial r}}\right)^{2}-{\frac {1}{r^{2}}}\left({\frac {\partial {\mathcal {S}}}{\partial \varphi }}\right)^{2}=m^{2}c^{2}$
dove abbiamo di nuovo orientato il sistema di coordinate sferiche con il piano dell'orbita. Il tempo $t$ e la longitudine $\varphi$ sono coordinate cicliche, in modo che la soluzione per l'azione possa essere scritta come:
    ${\mathcal {S}}=-Et+L\varphi +{\mathcal {S}}_{r}(r)\,$
dove $E$ e $L$ rappresentano di nuovo rispettivamente l'energia della particella e il momento angolare. L'equazione di Hamilton-Jacobi fornisce una soluzione integrale per la parte radiale dell'azione:
    ${\mathcal {\mathcal {S}}}_{r}(r)=\int {\frac {Ldr}{1-{\frac {r_{s}}{r}}}}{\sqrt {{\frac {1}{b^{2}}}-\left(1-{\frac {r_{s}}{r}}\right)\left({\frac {1}{a^{2}}}+{\frac {1}{r^{2}}}\right)}}.$
Prendendo la derivata dell'azione nel solito modo
    ${\frac {\partial {\mathcal {S}}}{\partial L}}=\varphi +{\frac {\partial {\mathcal {S}}_{r}}{\partial L}}=\mathrm {costante}$
si ottiene l'equazione orbitale derivata precedentemente:
    $\left({\frac {dr}{d\varphi }}\right)^{2}={\frac {r^{4}}{b^{2}}}-\left(1-{\frac {r_{s}}{r}}\right)\left({\frac {r^{4}}{a^{2}}}+r^{2}\right).$
Questo approccio può essere anche utilizzato per derivare la velocità (rate) di precessione orbitale in modo elegante.
Nel limite di massa zero $m$ (o, in modo equivalente, $a$ infinito), la parte radiale della funzione principale di Hamilton $S$ diventa:
    $S_{r}(r)={\frac {E}{c}}\int dr{\sqrt {{\frac {r^{2}}{\left(r-r_{s}\right)^{2}}}-{\frac {b^{2}}{r\left(r-r_{s}\right)}}}}$
da cui può essere derivata l'equazione per la deflessione della luce.

### Approccio lagrangiano

Nella relatività generale, le particelle libere di massa trascurabile $m$ seguono le geodetiche nello spazio-tempo, a causa del principio di equivalenza. Le geodetiche nello spazio-tempo sono definite come curve per cui le piccole variazioni locali nella loro coordinate (tenendo i loro eventi estremi fissi) non apportano modifiche significative nella loro lunghezza complessiva $s$. Ciò può essere espresso matematicamente usando il calcolo delle variazioni:
    $0=\delta s=\delta \int ds=\delta \int {\sqrt {g_{\mu \nu }{\frac {dx^{\mu }}{d\tau }}{\frac {dx^{\nu }}{d\tau }}}}d\tau =\delta \int {\sqrt {2T}}d\tau$
dove $\tau$ è il tempo proprio, $s=c\tau$ è la lunghezza dell'arco nello spazio-tempo e $T$ viene definito come:
    $2T=c^{2}=\left({\frac {ds}{d\tau }}\right)^{2}=g_{\mu \nu }{\frac {dx^{\mu }}{d\tau }}{\frac {dx^{\nu }}{d\tau }}=\left(1-{\frac {r_{s}}{r}}\right)c^{2}\left({\frac {dt}{d\tau }}\right)^{2}-{\frac {1}{1-{\frac {r_{s}}{r}}}}\left({\frac {dr}{d\tau }}\right)^{2}-r^{2}\left({\frac {d\varphi }{d\tau }}\right)^{2}$
in analogia con l'energia cinetica. Se la derivata rispetto al tempo proprio è rappresentata per brevità da un punto
    ${\dot {x}}^{\mu }={\frac {dx^{\mu }}{d\tau }}$
$T$ può essere scritta come:
    $2T=c^{2}=\left(1-{\frac {r_{s}}{r}}\right)c^{2}\left({\dot {t}}\right)^{2}-{\frac {1}{1-{\frac {r_{s}}{r}}}}\left({\dot {r}}\right)^{2}-r^{2}\left({\dot {\varphi }}\right)^{2}$
I fattori della costante (come $c$ o la radice quadrata di $2$) non incidono sulla risposta al problema variazionale, quindi, tenendo la variazione all'interno dell'integrale si ottiene il principio di Hamilton:
    $0=\delta \int {\sqrt {2T}}d\tau =\int {\frac {\delta T}{\sqrt {2T}}}d\tau ={\frac {1}{c}}\delta \int Td\tau .$
La soluzione del problema variazionale è data dalle equazioni di Lagranges:
    ${\frac {d}{d\tau }}\left({\frac {\partial T}{\partial {\dot {x}}^{\sigma }}}\right)={\frac {\partial T}{\partial x^{\sigma }}}.$
Quando applicate a $t$ e $\varphi$, queste equazioni rivelano due costanti di moto
    ${\frac {d}{d\tau }}\left[r^{2}{\frac {d\varphi }{d\tau }}\right]=0,$
    ${\frac {d}{d\tau }}\left[\left(1-{\frac {r_{s}}{r}}\right){\frac {dt}{d\tau }}\right]=0,$
le quali possomo essere scritte come le equazioni per $L$ e $E$
    $r^{2}{\frac {d\varphi }{d\tau }}={\frac {L}{m}},$
    $\left(1-{\frac {r_{s}}{r}}\right){\frac {dt}{d\tau }}={\frac {E}{mc^{2}}}.$
Come mostrato sopra, con la sostituzione di queste equazioni dentro la definizione della metrica di Schwarzschild si ottiene l'equazione per l'orbita.

### Principio di Hamilton

L'integrale dell'azione per una particella influenzata soltanto dalla gravità è:
    $S=\int {-mc^{2}d\tau }=-mc\int {c{\frac {d\tau }{dq}}dq}=-mc\int {{\sqrt {g_{\mu \nu }{\frac {dx^{\mu }}{dq}}{\frac {dx^{\nu }}{dq}}}}dq}$
dove $\tau$ è il tempo proprio e $q$ è ogni parametrizzazione regolare della linea di universo della particella. Se a questa si applica il calcolo delle variazioni, si ottengono di nuovo le equazioni per una geodetica. I calcoli sono semplificati, se prima prendiamo la variazione del quadrato dell'integrando. Per la metrica e le coordinate di questo caso, quella quadrata è:
    $\left(c{\frac {d\tau }{dq}}\right)^{2}=g_{\mu \nu }{\frac {dx^{\mu }}{dq}}{\frac {dx^{\nu }}{dq}}=\left(1-{\frac {r_{s}}{r}}\right)c^{2}\left({\frac {dt}{dq}}\right)^{2}-{\frac {1}{1-{\frac {r_{s}}{r}}}}\left({\frac {dr}{dq}}\right)^{2}-r^{2}\left({\frac {d\varphi }{dq}}\right)^{2}$
Prendendo la variazioni di questa, otteniamo:
    $\delta \left(c{\frac {d\tau }{dq}}\right)^{2}=2c^{2}{\frac {d\tau }{dq}}\delta {\frac {d\tau }{dq}}=\delta \left[\left(1-{\frac {r_{s}}{r}}\right)c^{2}\left({\frac {dt}{dq}}\right)^{2}-{\frac {1}{1-{\frac {r_{s}}{r}}}}\left({\frac {dr}{dq}}\right)^{2}-r^{2}\left({\frac {d\varphi }{dq}}\right)^{2}\right]$
Se variamo soltanto rispetto alla longitudine $\varphi$, otteniamo:
    $2c^{2}{\frac {d\tau }{dq}}\delta {\frac {d\tau }{dq}}=-2r^{2}{\frac {d\varphi }{dq}}\delta {\frac {d\varphi }{dq}}$
dividiamo per $2c{\frac {d\tau }{dq}}$ per ottenere la variazione dell'integrando stesso
    $c\,\delta {\frac {d\tau }{dq}}=-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}\delta {\frac {d\varphi }{dq}}=-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}{\frac {d\delta \varphi }{dq}}$
Perciò abbiamo:
    $0=\delta \int {c{\frac {d\tau }{dq}}dq}=\int {c\delta {\frac {d\tau }{dq}}dq}=\int {-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}{\frac {d\delta \varphi }{dq}}dq}$
L'integrazione delle parti dà
    $0=-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}\delta \varphi -\int {{\frac {d}{dq}}\left[-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}\right]\delta \varphi dq}$
La variazione della longitudine è assunta uguale a zero nei punti finali, così il primo termine scompare. L'integrale può essere fatto diverso da zero tramite una scelta ostinata (perverse) di $\delta \varphi$ salvo che l'altro fattore interno sia ovunque zero. Così otteniamo l'equazione di moto:
    ${\frac {d}{dq}}\left[-{\frac {r^{2}}{c}}{\frac {d\varphi }{d\tau }}\right]=0$
Se variamo soltanto rispetto al tempo $t$, otteniamo:
    $2c^{2}{\frac {d\tau }{dq}}\delta {\frac {d\tau }{dq}}=2\left(1-{\frac {r_{s}}{r}}\right)c^{2}{\frac {dt}{dq}}\delta {\frac {dt}{dq}}$
dividiamo per $2c{\frac {d\tau }{dq}}$ per ottenere la variazione dell'integrando stesso
    $c\delta {\frac {d\tau }{dq}}=c\left(1-{\frac {r_{s}}{r}}\right){\frac {dt}{d\tau }}\delta {\frac {dt}{dq}}=c\left(1-{\frac {r_{s}}{r}}\right){\frac {dt}{d\tau }}{\frac {d\delta t}{dq}}$
Così abbiamo:
    $0=\delta \int {c{\frac {d\tau }{dq}}dq}=\int {c\left(1-{\frac {r_{s}}{r}}\right){\frac {dt}{d\tau }}{\frac {d\delta t}{dq}}dq}$
L'integrazione delle parti dà
    $0=c\left(1-{\frac {r_{s}}{r}}\right){\frac {dt}{d\tau }}\delta t-\int {{\frac {d}{dq}}\left[c\left(1-{\frac {r_{s}}{r}}\right){\frac {dt}{d\tau }}\right]\delta tdq}$,
fornendo l'equazione di moto:
    ${\frac {d}{dq}}\left[c\left(1-{\frac {r_{s}}{r}}\right){\frac {dt}{d\tau }}\right]=0$
Se integriamo queste equazioni di moto e determiniamo le costanti di integrazione, otteniamo di nuovo le equazioni
    $r^{2}{\frac {d\varphi }{d\tau }}={\frac {L}{m}},$
    $\left(1-{\frac {r_{s}}{r}}\right){\frac {dt}{d\tau }}={\frac {E}{mc^{2}}}.$
Queste due equazioni per le costanti di moto $L$ ed $E$ possono essere combinate per formare un'equazione che è valida anche per i fotoni e altre particelle senza massa per cui il tempo proprio lungo una geodetica è zero.
    ${\frac {r^{2}}{bc}}{\frac {d\varphi }{dt}}=1-{\frac {r_{s}}{r}}$