andrewkirk said:
That sounds like a promising approach that ought to work. I tried doing this but, because of the complexity of L, the algebra soon became very ugly and I ran out of paper ( or patience - one or the other).
Maybe I took a wrong turn somewhere. Have you worked it through?
It's easy to get into a blind alley, but I don't think this is that bad.
With [itex]L = \sqrt{g_{\mu \nu} U^\mu U^\nu}[/itex],
- [itex]\dfrac{\partial}{\partial U^\mu} L = \dfrac{1}{L} g_{\mu \nu} U^\nu[/itex]
- [itex]\dfrac{\partial}{\partial x^\mu} L = \dfrac{1}{2L} \dfrac{\partial g_{\mu' \nu}}{\partial x^\mu} U^{\mu'} U^\nu[/itex]
So the Euler-Lagrange equations give:
[itex]\dfrac{1}{L} \dfrac{d}{ds} (g_{\mu \nu} U^\nu) + g_{\mu \nu} U^\nu \dfrac{d}{ds} \dfrac{1}{L} = \dfrac{1}{2L} \dfrac{\partial g_{\mu' \nu}}{\partial x^\mu} U^{\mu'} U^\nu[/itex]
Now, the great simplification comes from assuming [itex]\dfrac{d}{ds} \dfrac{1}{L} = 0[/itex], so [itex]L =[/itex] a constant along the path. With this assumption, the equation simplifies to (multiplying both sides by [itex]L[/itex])
[itex]\dfrac{d}{ds} (g_{\mu \nu} U^\nu) = \dfrac{1}{2} \dfrac{\partial g_{\mu' \nu}}{\partial x^\mu} U^{\mu'} U^\nu[/itex]
We expand the left-hand side to get:
[itex](\dfrac{d}{ds} g_{\mu \nu}) U^\nu + g_{\mu \nu} \dfrac{d}{ds} U^\nu = \dfrac{1}{2} \dfrac{\partial g_{\mu' \nu}}{\partial x^\mu} U^{\mu'} U^\nu[/itex]
Then you use: [itex]\dfrac{d}{ds} g_{\mu \nu} = (\dfrac{\partial}{\partial x^{\mu'}} g_{\mu \nu}) \dfrac{d}{ds} x^{\mu'} = (\dfrac{\partial}{\partial x^{\mu'}} g_{\mu \nu}) U^{\mu'}[/itex] (the chain rule for derivatives)
So we now have:
[itex](\dfrac{\partial}{\partial x^{\mu'}} g_{\mu \nu}) U^{\mu'} U^\nu + g_{\mu \nu} \dfrac{d}{ds} U^\nu = \dfrac{1}{2} \dfrac{\partial g_{\mu' \nu}}{\partial x^\mu} U^{\mu'} U^\nu[/itex]
Rearranging gives:
[itex]g_{\mu \nu} \dfrac{d}{ds} U^\nu = - ( (\dfrac{\partial}{\partial x^{\mu'}} g_{\mu \nu}) U^{\mu'} U^\nu - \dfrac{1}{2} \dfrac{\partial g_{\mu' \nu}}{\partial x^\mu} U^{\mu'} U^\nu)[/itex]
Finally, get rid of the [itex]g_{\mu \nu}[/itex] from the left-hand side by multiplying by the inverse matrix, [itex]g^{\mu \nu'}[/itex] and summing over [itex]\mu[/itex]. this gives:
[itex]\dfrac{d U^{\nu'}}{ds} = - \frac{1}{2} g^{\mu \nu'} (2 \dfrac{\partial g_{\mu \nu}}{\partial x^{\mu'}} - \dfrac{\partial g_{\mu' \nu}}{\partial x^\mu}) U^{\mu'} U^\nu[/itex]
So if we define [itex]Q^{\nu'}_{\nu \mu'} = \frac{1}{2} g^{\mu \nu'} (2 \dfrac{\partial g_{\mu \nu}}{\partial x^{\mu'}} - \dfrac{\partial g_{\mu' \nu}}{\partial x^\mu})[/itex], then this becomes:
[itex]\dfrac{d U^{\nu'}}{ds} = -Q^{\nu'}_{\nu \mu'} U^{\mu'} U^\nu[/itex]
Almost there! The usual connection coefficient is [itex]\Gamma^{\nu'}_{\nu \mu'} = \frac{1}{2} g^{\mu \nu'} (\dfrac{\partial g_{\mu \nu}}{\partial x^{\mu'}} + \dfrac{\partial g_{\mu \mu'}}{\partial x^{\nu}} - \dfrac{\partial g_{\mu' \nu}}{\partial x^\mu})[/itex], So we can write:
[itex]Q^{\nu'}_{\nu \mu'} = \Gamma^{\nu'}_{\nu \mu'} + \frac{1}{2} g^{\mu \nu'} (\dfrac{\partial g_{\mu \nu}}{\partial x^{\mu'}} - \dfrac{\partial g_{\mu \mu'}}{\partial x^{\nu}})[/itex]
Notice that the last term on the right is antisymmetric under the exchange [itex]\nu \Rightarrow \mu'[/itex]. On the other hand, [itex]U^{\mu'} U^\nu[/itex] is symmetric under that exchange. The product of an antisymmetric tensor with a symmetric tensor is zero. So:
[itex]Q^{\nu'}_{\nu \mu'} U^{\mu'} U^\nu = \Gamma^{\nu'}_{\nu \mu'} U^{\mu'} U^\nu[/itex]
So the Euler-Lagrange equations boil down to:
[itex]\dfrac{d U^{\nu'}}{ds} = -\Gamma^{\nu'}_{\nu \mu'} U^{\mu'} U^\nu[/itex]
which is the geodesic equation (whew!). Okay, I guess it was pretty bad...