In general relativity, force is a non-tensor, the time derivative of (kinetic) linear momentum,
![{\displaystyle F_{\alpha }={\frac {dp_{\alpha }}{dt}}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7c227188001358c76946fb8fb601174e90207cf1)
where t is any 'time' coordinate which parameterizes the trajectory of the particle (it does not have to be the same as any of x0, x1, x2, or x3).
To use Newton's third law of motion, both forces must be defined as the rate of change of momentum with respect to the same time coordinate. Consequently, the differentiation cannot be with respect to proper time for both particles since it is usually different for two particles.
The linear momentum of a particle is a covariant tensor. For particles which have mass, linear momentum is
![{\displaystyle p_{\alpha }=m\,g_{\alpha \beta }\,{\frac {dx^{\beta }}{d\tau }}=m\,g_{\alpha \beta }\,v^{\beta }{\frac {dt}{d\tau }}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e6186f543419c28bbe642d867ee3ebe72fd8a243)
where: m is the mass of the particle; gαβ is the metric tensor which is also the gravitational potential field; xβ is the position vector in 4D space-time; and τ is the proper time measured along the trajectory of the particle
![{\displaystyle {\frac {d\tau }{dt}}={\sqrt {{\frac {-1}{c^{2}}}g_{\mu \nu }{\frac {dx^{\mu }}{dt}}{\frac {dx^{\nu }}{dt}}}}\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/05e3c5954123c367a58a5dc6ed490c6cf0c154e1)
For particles without mass, all we can say about momentum is that it is parallel to velocity
![{\displaystyle p_{\alpha }\parallel g_{\alpha \beta }\,{\frac {dx^{\beta }}{dt}}\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d124271072a7ebd975dc3e2f4adea5e583300b20)
If t = x0, then
![{\displaystyle v^{\alpha }={\frac {dx^{\alpha }}{dt}}={\frac {v^{\alpha }}{v^{0}}}={\frac {g^{\alpha \beta }p_{\beta }}{g^{0\gamma }p_{\gamma }}}\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/007dcdb10cd4f10b958b15a256f4b28693c4ba50)
Force in special relativity
[edit]
To use Newton's third law of motion, both forces must be defined as the rate of change of momentum with respect to the same time coordinate. Unfortunately, the resulting object is not a tensor even with respect to Lorentz transformations.
In a continuous medium, the 3D density of force combines with the density of power to form a covariant 4-vector density. The spatial part is the result of dividing the force on a small cell (in 3-space) by the volume of that cell. The time component is negative of the power transferred to that cell divided by the volume of the cell.
Weight and fictitious force
[edit]
A point particle in free fall obeys the equation
![{\displaystyle {\frac {dp_{\mu }}{dt}}=\Gamma _{\mu \nu }^{\lambda }p_{\lambda }{\frac {dx^{\nu }}{dt}}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c41217e2192b09258b17fe76fcbc1a3ae9c3e5fe)
where pμ is the momentum 4-vector and Γλμν is the Christoffel symbol (which is the gravitational force field). Any separation of the inertial force on the right hand side into a part which is weight and a part which is fictitious force must be based on an arbitrary choice of a preferred reference frame where there are no fictitious forces by definition.
![{\displaystyle \Gamma _{\mu \nu }^{\lambda }p_{\lambda }{\frac {dx^{\nu }}{dt}}={\big (}\Gamma _{\mu \nu }^{\lambda }-{\widehat {\Gamma }}_{\mu \nu }^{\lambda }{\big )}p_{\lambda }{\frac {dx^{\nu }}{dt}}+{\widehat {\Gamma }}_{\mu \nu }^{\lambda }p_{\lambda }{\frac {dx^{\nu }}{dt}}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7b92cd63afa47a145c5f1b6c0cf5e9051573bdc7)
where the first term on the right side is the weight (gravitational force, a tensor (except for the dt)) and the second term on the right side is the fictitious force (a non-tensor which is zero in the preferred reference frame). The hatted Christoffel symbol is defined by using the Minkowski metric in the preferred reference frame rather than the physical metric used by the regular Christoffel symbol.
Classical force equation
[edit]
The non-relativistic version of the force equation is
![{\displaystyle {\frac {d{\vec {p}}}{dt}}=m{\vec {g}}+q{\vec {E}}+{\text{ other forces }}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4990b28ba58d31cc05be17c8c114645df60edbe0)
where
![{\displaystyle {\vec {p}}=m{\vec {v}}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/10dc8ad00c653c2fb15fb7628c0c1ba92a8fd8e4)
is the linear momentum,
![{\displaystyle {\vec {g}}=-{\vec {\nabla }}\zeta \,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/28e55dbf47387b16c69caa2bfd280d9780544c4d)
is the negative of the gradient of the gravitational potential field, and
![{\displaystyle {\vec {E}}=-{\vec {\nabla }}\phi \,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3f9033e392ffef340f49a0e5bf457f89451382c2)
is the negative of the gradient of the electrostatic potential field.
A test particle obeys the equation
![{\displaystyle {\frac {dp_{\mu }}{dt}}=\Gamma _{\mu \nu }^{\lambda }p_{\lambda }{\frac {dx^{\nu }}{dt}}+F_{\mu \nu }q{\frac {dx^{\nu }}{dt}}+{\text{ other forces }}\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/72e2169d05dbcb2060684877e2637bae4a1791ba)
The linear momentum of a particle is a covariant tensor. For particles which have mass, linear momentum is
![{\displaystyle p_{\alpha }=m\,g_{\alpha \beta }\,{\frac {dx^{\beta }}{d\tau }}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8a94e494c5b017f8fdbbbac2212bceebbfe608ee)
where: m is the mass of the particle; gαβ is the metric tensor which is also the gravitational potential field; xβ is the position vector in 4D space-time; and τ is the proper time measured along the trajectory of the particle
![{\displaystyle {\frac {d\tau }{dt}}={\sqrt {{\frac {-1}{c^{2}}}g_{\mu \nu }{\frac {dx^{\mu }}{dt}}{\frac {dx^{\nu }}{dt}}}}\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/05e3c5954123c367a58a5dc6ed490c6cf0c154e1)
Γλμν is the Christoffel symbol (which is the gravitational force field)
![{\displaystyle \Gamma _{\mu \nu }^{\lambda }={\frac {1}{2}}g^{\lambda \sigma }\left(g_{\sigma \mu ,\nu }+g_{\nu \sigma ,\mu }-g_{\mu \nu ,\sigma }\right)\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9fc35b95c34b1a54c3614f307748bec35b3de103)
q is the electric charge of the particle; and Fμν is the electromagnetic field
![{\displaystyle F_{\alpha \beta }\,=\,A_{\beta ,\alpha }\,-\,A_{\alpha ,\beta }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3aea885f94345a1104c7df3628314860aae13f5e)
where Aβ is the electromagnetic potential 4-vector.
Derivation from Lagrangian
[edit]
The Lagrangian of a general relativistic test particle in an electromagnetic field is
![{\displaystyle L[t]=-mc^{2}{\frac {d\tau [t]}{dt}}+q{\frac {dx^{\gamma }[t]}{dt}}A_{\gamma }[x[t]]\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/508dcffc4976739134c2034ebd8cdf9729c340b0)
where
![{\displaystyle {\frac {d\tau [t]}{dt}}={\sqrt {{\frac {-1}{c^{2}}}\,g_{\alpha \beta }[x[t]]{\frac {dx^{\alpha }[t]}{dt}}{\frac {dx^{\beta }[t]}{dt}}}}\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b538f38387ea0d1e154bd3df3c01cc2ca4d9881a)
Taking the variation of the Lagrangian with respect to variations, δxσ, in the trajectory of the particle, we get
![{\displaystyle \delta L=-mc^{2}{\frac {dt}{2d\tau }}\delta \left({\frac {-1}{c^{2}}}\,g_{\alpha \beta }{\frac {dx^{\alpha }}{dt}}{\frac {dx^{\beta }}{dt}}\right)+q{\frac {d\delta x^{\gamma }}{dt}}A_{\gamma }+q{\frac {dx^{\gamma }}{dt}}\delta A_{\gamma }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ce2947c81c80c12aacc4e1bdbcac9aa5953a63a6)
![{\displaystyle \delta L={\frac {mdt}{2d\tau }}\left(\delta g_{\alpha \beta }{\frac {dx^{\alpha }}{dt}}{\frac {dx^{\beta }}{dt}}+g_{\alpha \beta }{\frac {d\delta x^{\alpha }}{dt}}{\frac {dx^{\beta }}{dt}}+g_{\alpha \beta }{\frac {dx^{\alpha }}{dt}}{\frac {d\delta x^{\beta }}{dt}}\right)+q{\frac {d\delta x^{\gamma }}{dt}}A_{\gamma }+q{\frac {dx^{\gamma }}{dt}}\delta A_{\gamma }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/23b57b76d6e47a4acb0466a1624fc9b382770015)
![{\displaystyle \delta L={\frac {mdt}{2d\tau }}\left(g_{\alpha \beta ,\sigma }\delta x^{\sigma }{\frac {dx^{\alpha }}{dt}}{\frac {dx^{\beta }}{dt}}+2g_{\sigma \beta }{\frac {d\delta x^{\sigma }}{dt}}{\frac {dx^{\beta }}{dt}}\right)+q{\frac {d\delta x^{\sigma }}{dt}}A_{\sigma }+q{\frac {dx^{\gamma }}{dt}}A_{\gamma ,\sigma }\delta x^{\sigma }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/43e1eff8d2ce98591bd718af488816e48e478673)
![{\displaystyle \delta L={\frac {m}{2}}g_{\alpha \beta ,\sigma }{\frac {dt}{d\tau }}\delta x^{\sigma }{\frac {dx^{\alpha }}{dt}}{\frac {dx^{\beta }}{dt}}-{\frac {d}{dt}}\left(mg_{\sigma \beta }{\frac {dx^{\beta }}{d\tau }}\right)\delta x^{\sigma }-q{\frac {dA_{\sigma }}{dt}}\delta x^{\sigma }+q{\frac {dx^{\gamma }}{dt}}A_{\gamma ,\sigma }\delta x^{\sigma }+{\frac {d(\ldots )}{dt}}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/19b58acd03a46d586f71c93efde038d8875a3191)
![{\displaystyle \delta L=mg_{\rho \beta }\Gamma _{\alpha \sigma }^{\rho }{\frac {dt}{d\tau }}\delta x^{\sigma }{\frac {dx^{\alpha }}{dt}}{\frac {dx^{\beta }}{dt}}-{\frac {d}{dt}}\left(mg_{\sigma \beta }{\frac {dx^{\beta }}{d\tau }}\right)\delta x^{\sigma }-qA_{\sigma ,\gamma }{\frac {dx^{\gamma }}{dt}}\delta x^{\sigma }+q{\frac {dx^{\gamma }}{dt}}A_{\gamma ,\sigma }\delta x^{\sigma }+{\frac {d(\ldots )}{dt}}\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/38947b7a8879f3bc33297e5fad361527a5341c7b)
Thus the equation of motion is
![{\displaystyle 0=mg_{\rho \beta }\Gamma _{\alpha \sigma }^{\rho }{\frac {dt}{d\tau }}{\frac {dx^{\alpha }}{dt}}{\frac {dx^{\beta }}{dt}}-{\frac {d}{dt}}\left(mg_{\sigma \beta }{\frac {dx^{\beta }}{d\tau }}\right)-qA_{\sigma ,\gamma }{\frac {dx^{\gamma }}{dt}}+q{\frac {dx^{\gamma }}{dt}}A_{\gamma ,\sigma }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6d8d2e28eae8e76d02197d8cf098711f2f0c301c)
![{\displaystyle {\frac {dp_{\sigma }}{dt}}=\Gamma _{\alpha \sigma }^{\rho }p_{\rho }{\frac {dx^{\alpha }}{dt}}+qF_{\sigma \gamma }{\frac {dx^{\gamma }}{dt}}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d3a4fa1fd3c829496e0d11076a63b86464afd1ef)
Changing notation and letting fσ be the total of all non-gravitational forces, we get
![{\displaystyle {\frac {d}{dt}}\left(mg_{\sigma \beta }\gamma v^{\beta }\right)=mg_{\rho \beta }\Gamma _{\alpha \sigma }^{\rho }\gamma v^{\alpha }v^{\beta }+f_{\sigma }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/271504acb97c69c27e807414e77eef35fb5354d4)
![{\displaystyle mg_{\sigma \beta ,\alpha }v^{\alpha }\gamma v^{\beta }+mg_{\sigma \beta }{\frac {d\gamma }{dt}}v^{\beta }+mg_{\sigma \beta }\gamma a^{\beta }=mg_{\rho \beta }\Gamma _{\alpha \sigma }^{\rho }\gamma v^{\alpha }v^{\beta }+f_{\sigma }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/57e675da712a9ff4cd5c3585ad6da9d101504ae2)
![{\displaystyle mg_{\sigma \beta }\gamma a^{\beta }=-mg_{\sigma \mu ,\alpha }v^{\alpha }\gamma v^{\mu }-mg_{\sigma \beta }{\frac {d\gamma }{dt}}v^{\beta }+mg_{\rho \mu }\Gamma _{\alpha \sigma }^{\rho }\gamma v^{\alpha }v^{\mu }+f_{\sigma }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2630d510ba96604aa6146d0890ad5a116bd6332f)
![{\displaystyle a^{\beta }=-g^{\beta \sigma }g_{\sigma \mu ,\alpha }v^{\alpha }v^{\mu }-{\frac {d\gamma }{\gamma dt}}v^{\beta }+g^{\beta \sigma }g_{\rho \mu }\Gamma _{\alpha \sigma }^{\rho }v^{\alpha }v^{\mu }+{\frac {1}{m\gamma }}g^{\beta \sigma }f_{\sigma }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/56bfdb9ced452edf2721065aef3094e5c3a0776c)
using
![{\displaystyle g_{\sigma \mu ,\alpha }=g_{\rho \mu }\Gamma _{\sigma \alpha }^{\rho }+g_{\sigma \rho }\Gamma _{\mu \alpha }^{\rho }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bbb97e37ecade17a575a95b6ac308472ad5cd2a4)
gives
![{\displaystyle a^{\beta }=-\Gamma _{\alpha \mu }^{\beta }v^{\alpha }v^{\mu }-{\frac {d\gamma }{\gamma dt}}v^{\beta }+{\frac {1}{m\gamma }}g^{\beta \sigma }f_{\sigma }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9846f0b11ba49ad12dc3b8421d20a49329975831)
substituting β=0 and v0=1 and a0=0 gives
![{\displaystyle 0=-\Gamma _{\alpha \mu }^{0}v^{\alpha }v^{\mu }-{\frac {d\gamma }{\gamma dt}}+{\frac {1}{m\gamma }}g^{0\sigma }f_{\sigma }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/88ad95a0c61137e83bcd5bf6d3b70b8b823d3450)
subtracting vβ times this from the previous equation gives
![{\displaystyle a^{\beta }=-\Gamma _{\alpha \mu }^{\beta }v^{\alpha }v^{\mu }+\Gamma _{\alpha \mu }^{0}v^{\alpha }v^{\mu }v^{\beta }+{\frac {1}{m\gamma }}g^{\beta \sigma }f_{\sigma }-{\frac {1}{m\gamma }}g^{0\sigma }f_{\sigma }v^{\beta }\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/09dcda8e5692baf7bbab6911da82cd92c2af603c)
which is the formula for the acceleration of a particle in general relativity.
An alternative derivation
[edit]
Let the Lagrangian be:
![{\displaystyle L=p_{\mu }[t]\,v^{\mu }[t]+\Lambda [t]\,(p_{\alpha }[t]\,g^{\alpha \beta }[x^{.}[t]]\,p_{\beta }[t]+m^{2}c^{2})+L_{\text{int}}[x^{.}[t],v^{.}[t]]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7635a8b68f4821fa3dee8a7c60966654d807d629)
where
.
Variation of the integral of L with respect to p gives:
![{\displaystyle 0=v^{\mu }+2\,\Lambda \,g^{\mu \beta }\,p_{\beta }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1946ab330bcaacb701b9aaf13d8bf26d38b482ea)
.
Variation of the integral of L with respect to Λ gives the constraint equation:
![{\displaystyle 0=p_{\alpha }\,g^{\alpha \beta }\,p_{\beta }+m^{2}c^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5883b6bce4f45e9fa4e4fd0256ca09424dd0ed39)
![{\displaystyle -m^{2}c^{2}\,\Lambda ^{2}={\frac {1}{4}}v^{\alpha }\,g_{\alpha \beta }\,v^{\beta }={\frac {-c^{2}}{4}}\left({\frac {d\tau }{dt}}\right)^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f655a14de9cc24a49d76f0cf11df09b06b22b248)
![{\displaystyle m^{2}\,\Lambda ^{2}={\frac {1}{4}}\left({\frac {d\tau }{dt}}\right)^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5ea6386352269d43f7edf74cecaa07af75be8ab3)
![{\displaystyle m\,\Lambda =-{\frac {1}{2}}{\frac {d\tau }{dt}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e80dc2ddc16b1ca8c2ed792e3c1516e2ccbb2614)
.
Variation of the integral with respect to x gives:
![{\displaystyle 0=-{\frac {dp_{\mu }}{dt}}+\Lambda \,p_{\alpha }\,g^{\alpha \beta }{}_{,\mu }\,p_{\beta }+{\frac {\partial L_{\text{int}}}{\partial x^{\mu }}}-{\frac {d}{dt}}{\frac {\partial L_{\text{int}}}{\partial v^{\mu }}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3610689db3a9b2b7483fba2c50617f230b95ad6f)
![{\displaystyle {\frac {dp_{\mu }}{dt}}=\Lambda \,p_{\alpha }\,\left(-\Gamma _{\rho \mu }^{\alpha }g^{\rho \beta }-\Gamma _{\sigma \mu }^{\beta }g^{\alpha \sigma }\right)\,p_{\beta }+f_{\mu }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/53d73c4f5a51df54a9f6f30278e19bf5d66ad235)
where the non-gravitational forces on the particle are:
![{\displaystyle f_{\mu }={\frac {\partial L_{\text{int}}}{\partial x^{\mu }}}-{\frac {d}{dt}}{\frac {\partial L_{\text{int}}}{\partial v^{\mu }}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/26bb1e0610b822f11bd50ee8e132f6be46237458)
![{\displaystyle {\frac {dp_{\mu }}{dt}}=\Gamma _{\rho \mu }^{\alpha }\,p_{\alpha }\,v^{\rho }+f_{\mu }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bc7ae59611307b991fe191de17c16e2e65b0f017)
Is the equation of motion consistent with the constraint?
![{\displaystyle =2f_{\alpha }\,g^{\alpha \beta }\,p_{\beta }+\Gamma _{\rho \alpha }^{\sigma }\,p_{\sigma }\,v^{\rho }\,g^{\alpha \beta }\,p_{\beta }+p_{\alpha }\,\left(-\Gamma _{\rho \gamma }^{\alpha }g^{\rho \beta }-\Gamma _{\sigma \gamma }^{\beta }g^{\alpha \sigma }\right)v^{\gamma }\,p_{\beta }+p_{\alpha }\,g^{\alpha \beta }\,\Gamma _{\rho \beta }^{\sigma }\,p_{\sigma }\,v^{\rho }=2f_{\alpha }\,g^{\alpha \beta }\,p_{\beta }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e3140d47ab89c57b67154105131e7eafc96c76b9)
![{\displaystyle 0=f_{\alpha }v^{\alpha }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/17136c810e8d95b8576e340206fa4425abfc14a0)
Only if the non-gravitational force is orthogonal to the velocity 4-vector! That is, in an inertial frame where the particle is instantaneously at rest, the non-gravitational force has no time component.
The only interaction Lagrangian consistent with this restriction on the resulting non-gravitational forces is a sum of terms of the form
![{\displaystyle qA_{\alpha }[x^{.}[t]]v^{\alpha }[t]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e36f2ecfb8522578eac5bc2c0d0a0094e7675f70)
where q is a constant charge or coupling constant, and A is a 4-vector potential field. In this case, the resulting forces will have the same form as the Lorentz force.
Null geodesics from what?
[edit]
Assume that a future-oriented path (hopefully a geodesic) is smoothly parameterized by strictly increasing λ, t=x0, and θ. Consider the integral
.
Suppose its variations with respect to xρ and θ are zero (while holding the end-points fixed):
![{\displaystyle 0=g_{\alpha \beta ,\rho }{\frac {dx^{\alpha }}{d\lambda }}{\frac {dx^{\beta }}{d\lambda }}{\frac {d\lambda }{d\theta }}-2{\frac {d}{d\lambda }}\left(g_{\rho \beta }{\frac {dx^{\beta }}{d\lambda }}{\frac {d\lambda }{d\theta }}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a681e4cdf270ce5dbafb6d60ee839a96b584af46)
.
The second equation is satisfied, if there is a constant K such that:
![{\displaystyle g_{\alpha \beta }{\frac {dx^{\alpha }}{d\theta }}{\frac {dx^{\beta }}{d\theta }}=K}](https://wikimedia.org/api/rest_v1/media/math/render/svg/42da64586e29a553d3ffe9e55f57c4f0eeadcbfa)
that is,
.
The first equation gives
![{\displaystyle 0={\tfrac {1}{2}}g_{\alpha \beta ,\rho }{\frac {dx^{\alpha }}{d\lambda }}{\frac {dx^{\beta }}{d\lambda }}{\frac {d\lambda }{d\theta }}-g_{\rho \beta ,\alpha }{\frac {dx^{\alpha }}{d\lambda }}{\frac {dx^{\beta }}{d\lambda }}{\frac {d\lambda }{d\theta }}-g_{\rho \beta }{\frac {d^{2}x^{\beta }}{{d\lambda }^{2}}}{\frac {d\lambda }{d\theta }}+g_{\rho \beta }{\frac {dx^{\beta }}{d\lambda }}{\frac {d\lambda }{d\theta }}{\frac {d\lambda }{d\theta }}{\frac {d^{2}\theta }{{d\lambda }^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a132ffa2404dd29c9164bf9a43b3e69d91049fcf)
.
Use the Christoffel symbols:
.
Use the chain rule to put this in terms of t:
![{\displaystyle 0=-g_{\rho \sigma }\Gamma _{\alpha \beta }^{\sigma }{\frac {dx^{\alpha }}{dt}}{\frac {dx^{\beta }}{dt}}-g_{\rho \beta }{\frac {d^{2}x^{\beta }}{{dt}^{2}}}-g_{\rho \beta }{\frac {dx^{\beta }}{dt}}{\frac {d^{2}t}{{d\lambda }^{2}}}{\frac {d\lambda }{dt}}{\frac {d\lambda }{dt}}+g_{\rho \beta }{\frac {dx^{\beta }}{d\theta }}{\frac {d^{2}\theta }{{dt}^{2}}}+g_{\rho \beta }{\frac {dx^{\beta }}{d\theta }}{\frac {d\theta }{dt}}{\frac {d^{2}t}{{d\lambda }^{2}}}{\frac {d\lambda }{dt}}{\frac {d\lambda }{dt}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/48b6c5405cd08923fb7dd41ae6455cf1ad4cb9b8)
.
.
.
Replacing σ by 0, x0 by t, dt/dt by 1, d2t/dt2 by 0 gives:
.
.
Subtract this from the second previous equation to get
![{\displaystyle {\frac {d^{2}x^{\sigma }}{{dt}^{2}}}=-\Gamma _{\alpha \beta }^{\sigma }{\frac {dx^{\alpha }}{dt}}{\frac {dx^{\beta }}{dt}}+\Gamma _{\alpha \beta }^{0}{\frac {dx^{\alpha }}{dt}}{\frac {dx^{\beta }}{dt}}{\frac {dx^{\sigma }}{dt}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b77aee7e9c158744c286e0051f47e0ff09a4b231)
which is the geodesic equation. We also get a formula for θ:
.
.