## Locally inertial frames (LIFs)

Reference: Moore, Thomas A., A General Relativity Workbook, University Science Books (2013) – Chapter 17; g.

We’ve seen that for a smooth metric, we can define a locally flat coordinate system at a given event in which the metric is the flat spacetime metric ${\eta_{ij}}$. Although we specify the value of the metric tensor at that point, we haven’t imposed any conditions on its derivatives with respect to the various coordinates. If we impose the addition constraint that all the derivatives of the metric are zero, then we have a locally inertial frame or LIF.

It’s not immediately obvious that we can impose this restriction, so let’s see how we can demonstrate this. We start with some arbitrary metric ${g_{ij}}$ which doesn’t satisfy the LIF conditions and consider the transformation to another metric ${g'_{ij}}$ on which we try to impose these conditions. As usual, the transformation is

$\displaystyle g'_{ij}=\partial_{i}^{'}x^{a}\partial_{j}^{'}x^{b}g_{ab} \ \ \ \ \ (1)$

Each of the partial derivatives is a function of the primed coordinates so, for a region close to the event point ${P}$, we can expand these derivatives in Taylor series:

$\displaystyle \partial_{i}^{'}x^{a}\left(x'_{P}+\Delta x'\right)=a_{i}^{a}+b_{ij}^{a}\Delta x'^{j}+c_{ijk}^{a}\Delta x'^{j}\Delta x'^{k}+... \ \ \ \ \ (2)$

where the coefficients are defined in terms of ${\partial'_{i}x^{a}}$ and its derivatives, all evaluated at ${P}$:

 $\displaystyle a_{i}^{a}$ $\displaystyle =$ $\displaystyle \partial_{i}^{'}x^{a}\left(x'_{P}\right)\ \ \ \ \ (3)$ $\displaystyle b_{ij}^{a}$ $\displaystyle =$ $\displaystyle \partial'_{j}\partial_{i}^{'}x^{a}\left(x'_{P}\right)\ \ \ \ \ (4)$ $\displaystyle c_{ijk}^{a}$ $\displaystyle =$ $\displaystyle \frac{1}{2}\partial'_{k}\partial'_{j}\partial_{i}^{'}x^{a}\left(x'_{P}\right) \ \ \ \ \ (5)$

We can also think of the unprimed metric as a function of the primed coordinates, and expand it in a Taylor series as well:

$\displaystyle g_{ab}\left(x'_{P}+\Delta x'\right)=g_{ab}\left(x'_{P}\right)+\Delta x'^{c}\partial'_{c}g_{ab}\left(x'_{P}\right)+\frac{1}{2}\Delta x'^{d}\Delta x'^{c}\partial'_{d}\partial'_{c}g_{ab}\left(x'_{P}\right)+... \ \ \ \ \ (6)$

We can now substitute these two series into 1 and collect terms. Due to the large number of indices floating about, it’s easier to use a condensed notation for this step. We can temporarily drop the indices to get

 $\displaystyle \partial'x$ $\displaystyle =$ $\displaystyle a+b\Delta x'+c\left(\Delta x'\right)^{2}\ \ \ \ \ (7)$ $\displaystyle g$ $\displaystyle =$ $\displaystyle g_{P}+\Delta x'\partial'g_{P}+\frac{1}{2}\left(\Delta x'\right)^{2}\partial'^{2}g_{P} \ \ \ \ \ (8)$

The transformation now becomes, up to second order terms:

 $\displaystyle g'\left(x'_{P}+\Delta x'\right)$ $\displaystyle =$ $\displaystyle \left(\partial'x\right)^{2}g\ \ \ \ \ (9)$ $\displaystyle$ $\displaystyle =$ $\displaystyle \left(a+b\Delta x'+c\left(\Delta x'\right)^{2}\right)\left(a+b\Delta x'+c\left(\Delta x'\right)^{2}\right)\left(g_{P}+\Delta x'\partial'g_{P}+\frac{1}{2}\left(\Delta x'\right)^{2}\partial'^{2}g_{P}\right)\ \ \ \ \ (10)$ $\displaystyle$ $\displaystyle =$ $\displaystyle a^{2}g_{P}+\Delta x'\left[abg_{P}+bag_{P}+a^{2}\partial'g_{P}\right]+\nonumber$ $\displaystyle$ $\displaystyle$ $\displaystyle \left(\Delta x'\right)^{2}\left[\frac{1}{2}a^{2}\partial'^{2}g_{P}+ab\partial'g_{P}+acg_{P}+ba\partial'g_{P}+b^{2}g_{P}+cag_{P}\right] \ \ \ \ \ (11)$

The last equation is the Taylor expansion of ${g'\left(x'_{P}+\Delta x'\right)}$ where every factor is now a function of ${x'}$. That is

$\displaystyle g'_{ij}\left(x'_{P}+\Delta x'\right)=g'_{ij}\left(x'_{P}\right)+\Delta x'^{k}\partial'_{k}g'_{ij}\left(x'_{P}\right)+\Delta x'^{\ell}\Delta x'^{k}\partial'_{\ell}\partial'_{k}g'_{ij}\left(x'_{P}\right)$

We therefore have, taking each term separately and restoring the indices:

 $\displaystyle g'_{ij}\left(x'_{P}\right)$ $\displaystyle =$ $\displaystyle a^{2}g_{P}\ \ \ \ \ (12)$ $\displaystyle$ $\displaystyle =$ $\displaystyle a_{i}^{a}a_{j}^{b}\left(g_{ab}\right)_{P}\ \ \ \ \ (13)$ $\displaystyle \partial'_{k}g'_{ij}\left(x'_{P}\right)$ $\displaystyle =$ $\displaystyle abg_{P}+bag_{P}+a^{2}\partial'g_{P}\ \ \ \ \ (14)$ $\displaystyle$ $\displaystyle =$ $\displaystyle \left[a_{i}^{a}b_{jk}^{b}+b_{ik}^{a}a_{j}^{b}\right]\left(g_{ab}\right)_{P}+a_{i}^{a}a_{j}^{b}\partial'_{k}\left(g_{ab}\right)_{P}\ \ \ \ \ (15)$ $\displaystyle \partial'_{\ell}\partial'_{k}g'_{ij}\left(x'_{P}\right)$ $\displaystyle =$ $\displaystyle \frac{1}{2}a^{2}\partial'^{2}g_{P}+ab\partial'g_{P}+acg_{P}+ba\partial'g_{P}+b^{2}g_{P}+cag_{P}\ \ \ \ \ (16)$ $\displaystyle$ $\displaystyle =$ $\displaystyle \frac{1}{2}a_{i}^{a}a_{j}^{b}\partial'_{k}\partial'_{\ell}\left(g_{ab}\right)_{P}+a_{i}^{a}b_{jk}^{b}\partial'_{\ell}\left(g_{ab}\right)_{P}+a_{i}^{a}c_{jk\ell}^{b}\left(g_{ab}\right)_{P}+\nonumber$ $\displaystyle$ $\displaystyle$ $\displaystyle b_{ik}^{a}a_{j}^{b}\partial'_{\ell}\left(g_{ab}\right)_{P}+b_{ik}^{a}b_{j\ell}^{b}\left(g_{ab}\right)_{P}+c_{ik\ell}^{a}a_{j}^{b}\left(g_{ab}\right)_{P} \ \ \ \ \ (17)$

We won’t actually try to solve these equations, but what we need to do is see if we have enough freedom to impose the conditions above, that is, that ${g'_{ij}\left(x'_{P}\right)=\eta_{ij}}$ and ${\partial'_{k}g'_{ij}\left(x'_{P}\right)=0}$. To do this, we need to see how many independent equations each of these equations provides.

Consider 13 first. Each index ${i}$ and ${j}$ has 4 possible values, but because ${g'_{ij}=g'_{ji}}$ we have to exclude permutations of ${i}$ and ${j}$. We can get these using binomial coefficients:

$\displaystyle N=\binom{4}{1}+\binom{4}{2}=4+6=10$

For 15, the symmetry of ${g'}$ means we have ${4\times10=40}$ (10 choices for ${i}$ and ${j}$, times the 4 different values of ${k}$) independent equations. Finally, for 17, since the order of partial derivatives doesn’t matter and the metric is symmetric, we have ${10\times10=100}$ independent equations.

The original metric ${g_{ij}}$ is assumed to be given, so we can’t change that, so how many variables can we solve for in these equations? The transformation is effectively determined (up to second order) by the coefficients ${a_{i}^{a}}$, ${b_{ij}^{a}}$ and ${c_{ijk}^{a}}$, so we need to figure out how many of these are independent of each other.

For ${a_{i}^{a}}$ every choice of the two indices gives an independent quantity, so there are 16 different variables. Thus we have 6 more degrees of freedom than we need to solve 13 and set ${g'_{ij}=\eta_{ij}}$.

For ${b_{ij}^{a}}$, the order of the partial derivatives doesn’t matter so there are 10 possible choices of ${i}$ and ${j}$ for each value of ${a}$, giving ${10\times4=40}$. Having chosen the ${a_{i}^{a}}$ in solving 13, we can plug these into the 40 equations specified by 15 and solve for the 40 ${b_{ij}^{a}}$ coefficients. Therefore, in principle, we can set ${\partial'_{k}g'_{ij}\left(x'_{P}\right)=0}$.

When it comes to ${c_{ijk}^{a}}$, again the order of the 3 partial derivatives doesn’t matter, so for each value of ${a}$, we can have ${i=j=k}$, or two indices equal with the third different, or all three indices different. The total number of choices for each value of ${a}$ is then

$\displaystyle N_{c}=\binom{4}{1}+2\binom{4}{2}+\binom{4}{3}=20$

To get the middle term, note that ${\binom{4}{2}}$ is the number of ways of choosing 2 numbers out of 4, but it assumes there are only 2 slots into which these numbers can be placed. In our case, one of the numbers is used twice, so this term would imply that, for example, 1-1-2 and 2-2-1 are just permutations of the same combination, whereas we want them to be distinct, so we have to multiply by 2. The total number of independent ${c_{ijk}^{a}}$ is therefore ${4\times20=80}$. The fact that there are 100 independent equations specified by 17 means that we can’t impose values on all of the ${\partial'_{\ell}\partial'_{k}g'_{ij}\left(x'_{P}\right)}$ so in general even if we require 80 of them to be zero, the other 20 may have non-zero values. It turns out that this is a result of the inherent curvature of the spacetime, so it’s not surprising that we can’t get rid of that.

The LIF imposes stricter conditions on the local metric than the locally flat frame we mentioned at the start. The latter just specifies the metric at a given point, while the LIF specifies both the metric and its derivatives.