Creo que su confusión queda recogida de forma más concisa en esta parte de su pregunta:
Si insisto en utilizar estas coordenadas para describir mi espacio, no veo cómo puede quedar ninguna libertad galvánica.
Sin entrar en la discusión que ha surgido en los comentarios a tu pregunta sobre si tu elección particular de coordenadas tiene sentido, es generalmente cierto que una vez que eliges un conjunto de coordenadas, haces no que le quede alguna libertad de calibre. La libertad gauge en la relatividad general es exactamente la libertad de elegir sus coordenadas. Ni más ni menos.
Creo que lo más natural es verlo a través de una descomposición ADM a 3+1 dimensiones. (La parte de tu pregunta sobre las ondas gravitacionales también se computa en este formalismo, o en una variante del mismo). Asumiendo por simplicidad en este punto una solución de vacío, si $g_{\mu\nu}$ es la métrica de su espaciotiempo, Arnowitt, Deser y Misner demostraron en los años 50 que si se toma una foliación de hipersuperficies parecidas al espacio, se puede descomponer en una métrica 3 inducida en las superficies $\gamma_{ij}$ más y la función "lapso" $\alpha$ y un vector de desplazamiento $\beta^i$ . Estos están relacionados por $$ \left( \begin{array}{cc} g_{00} & g_{0j} \\ g_{i0} & g_{ij} \end{array} \right) = \left( \begin{array}{cc} \beta_k \beta^k - \alpha^2 & \beta_j \\ \beta_i & \gamma_{ij} \end{array} \right) $$ donde el índice del turno se rebaja con la métrica 3, $\beta_k = \gamma_{ik} \beta^i$ . También hay un campo conjugado $\pi^{ij}$ en este formalismo que permiten las derivadas temporales (de coordenadas) $\partial_t \gamma_{ij}$ para ser expresadas como EDP de primer orden en el tiempo.
No voy a escribir aquí todas las ecuaciones porque están fácilmente disponibles, pero hacer esta división divide las ecuaciones de Einstein en dos tipos de ecuaciones: Hay dos ecuaciones de "evolución" para $\gamma_{ij}$ y para $\pi^{ij}$ y hay dos ecuaciones de restricción. Las restricciones son una restricción escalar llamada restricción hamiltoniana y una restricción de 3 vectores, conocida como restricción de momento. (Esto es análogo, respectivamente, en EM a las dos ecuaciones de Maxwell que tienen derivadas temporales y a las dos que no las tienen).
Tanto en la formulación 3+1 de la RG como en las ecuaciones de Maxwell, se puede demostrar que las "restricciones se propagan". Esto significa que si las restricciones se satisfacen en un momento dado, se satisfacen en todo momento. Véase, por ejemplo, el apéndice de Wald, que trata esto con gran detalle.
Ahora bien, tenga en cuenta que el lapso y el desplazamiento son sus campos "gauge" en este formalismo . No tienen ningún significado físico. En particular, no aparecen en la restricción hamiltoniana ni en la restricción de momento, lo que es análogo al hecho de que los campos gauge EM no aparecen en las ecuaciones de restricción $\nabla \cdot E = 0$ o $\nabla \cdot B = 0$ (de nuevo asumiendo el vacío). El hecho de que no aparezca en las ecuaciones de restricción significa que si tiene alguna ( $\gamma_{ij}$ , $\pi^{ij}$ ) que satisfagan las ecuaciones de restricción en cualquier momento, es decir, que sean "físicas" en cualquier momento, se puede construir una solución a las ecuaciones de Einstein completas para cualquier elección de lapso y turno resolviendo las ecuaciones de evolución como un problema de valor inicial, debido a que las restricciones se propagan, como se mencionó anteriormente.
En segundo lugar, Obsérvese que la elección del lapso y el desplazamiento son esencialmente equivalentes a una elección de coordenadas . Se pueden encontrar figuras de esto en muchas referencias, incluyendo MTW y Wald, pero básicamente se puede demostrar que entre dos superficies "cercanas" $\Sigma_0$ y $\Sigma_1$ un punto $x$ en $\Sigma_0$ tendrá la misma coordenada espacial que un punto $y$ en $\Sigma_1$ si $y \approx x + \alpha n + \beta$ , donde $n$ es normal que $\Sigma_0$ en $x$ . También tienes que el lapso es proporcional al tiempo propio (en lugar del tiempo de coordenadas) entre rodajas.
Así que, en cierto sentido, tienes razón en tu queja. Una vez que fijas las coordenadas, no te queda libertad de calibre. Cuando construyas tu "LIGO personal", según el final de tu pregunta, habrás hecho efectivamente tu elección de coordenadas, ya que las coordenadas "naturales" para tu configuración estarán determinadas por la forma en que dispongas los brazos de tu detector, dándote un sistema de coordenadas rectangulares en tus "rebanadas espaciales". En la región de campo débil de la Tierra, eso será completamente compatible con el lapso unitario y el desplazamiento cero, es decir, el espacio de Minkowski que habrías linealizado para calcular tu solución de onda en el detector.
Incluso con esta comprensión, todavía se pueden encontrar los efectos del calibre. Cuando finalmente detectes una onda en tu LIGO personal, tendrá una polarización, que probablemente descompondrás en polarizaciones "plus" y "cross". Si giras tu detector, lo cual es efectivamente una elección de calibre diferente, mezclarás esas polarizaciones de manera diferente para una señal dada. Este es probablemente el más sencillo de los múltiples factores que tienes que "ajustar" entre tu cálculo teórico de la onda -donde tienes que fijar el calibre- y tu configuración real del detector sobre el terreno. Este es el otro extremo de tu error/confusión. Cuando haces la teoría, si quieres que coincida con tu detector, no sólo tienes que elegir un calibre, sino que tienes que elegir un calibre que sea compatible con la forma en que vas a leer los resultados de tu detector. No puedes elegir cualquier calibre que te guste e intentar utilizar los resultados del cálculo directamente con tu detector.