Determination of Riemann curvature tensor from tidal forces

jdstokes
Messages
520
Reaction score
1
Hi,

Given a large number of test particles N, it should be possible to determine the Riemann curvature tensor by tracking their motion as they undergo geodesic deviation.

Is there a minimum number N that will achieve this in any situation, or does it vary from problem to problem?

How difficult is it to compute this least N?
 
Physics news on Phys.org
I have a partial answer to your question. I found this reasoning in https://www.amazon.com/dp/0716703440/?tag=pfamazon01-20 so don't hesitate to read the page 72 if you find my explanation too messy, or simply wrong.

The geodesic deviation equation is

\frac{D^{2}\xi_{\alpha}}{d\tau^{2}} = -R_{\alpha\beta\gamma\delta}u^{\beta}\xi^{\gamma}u^{\delta} in any coordinate system.

The idea is to use the method of least squares to determine the components of the Riemann tensor.

We define the 'error' between the observed value and the predicted value for the Kth particle by

\ddot{\xi}_{\alpha}^{\phantom{\alpha}K} + R_{\alpha\beta\gamma\delta}u^{\beta,K}\xi^{\gamma,K}u^{\delta,K} = \delta a_{\alpha}^{\phantom{\alpha}K}

The expression you now have to minimize is the sum of the squares of the \delta a_{\alpha}^{\phantom{\alpha}K}

S = \sum_{K}\eta^{\alpha\beta}\delta a_{\alpha}^{\phantom{\alpha}K}\delta a_{\beta}^{\phantom{\beta}K}

In this expression, everything is regarded as known, except the components R_{\alpha\beta\gamma\delta}. Minimize with respect to these components. It gives a system of 20 equations (the 20 independant components out of the 256).

It should be possible to determine the minimum number of particles by checking for which value of K the determinant of the 20 equations system 'stops to vanish'. However, this is not very elegant and might prove difficult to achieve (I haven't tried).
 
Last edited by a moderator:
Interesting. I will try to do this after trying the Lorentz force version.

By the way shouldn't the expression for the sum of the square involve a g^{\alpha\beta} rather than \eta^{\alpha\beta}.
 
I'm a little bit confused by this box in MTW.

In the first column they seem to make it out as though 2 particles is insufficient to determine all the components of F^{\alpha\beta} (by boosting to one particle's frame and then rotating spatial axes so that one axis is aligned with the motion of the second particle). Does this mean that 3 is sufficient?
 
Yes, the box is a bit confusing. At the end of the box, the author says that with 2 particles the determinant vanishes. I suppose it means that for 3 particles the system becomes solvable (the determinant is different from 0).

For the g^{\alpha\beta} instead of the \eta^{\alpha\beta} I think you're right.
 
I started reading a National Geographic article related to the Big Bang. It starts these statements: Gazing up at the stars at night, it’s easy to imagine that space goes on forever. But cosmologists know that the universe actually has limits. First, their best models indicate that space and time had a beginning, a subatomic point called a singularity. This point of intense heat and density rapidly ballooned outward. My first reaction was that this is a layman's approximation to...
Thread 'Dirac's integral for the energy-momentum of the gravitational field'
See Dirac's brief treatment of the energy-momentum pseudo-tensor in the attached picture. Dirac is presumably integrating eq. (31.2) over the 4D "hypercylinder" defined by ##T_1 \le x^0 \le T_2## and ##\mathbf{|x|} \le R##, where ##R## is sufficiently large to include all the matter-energy fields in the system. Then \begin{align} 0 &= \int_V \left[ ({t_\mu}^\nu + T_\mu^\nu)\sqrt{-g}\, \right]_{,\nu} d^4 x = \int_{\partial V} ({t_\mu}^\nu + T_\mu^\nu)\sqrt{-g} \, dS_\nu \nonumber\\ &= \left(...
In Philippe G. Ciarlet's book 'An introduction to differential geometry', He gives the integrability conditions of the differential equations like this: $$ \partial_{i} F_{lj}=L^p_{ij} F_{lp},\,\,\,F_{ij}(x_0)=F^0_{ij}. $$ The integrability conditions for the existence of a global solution ##F_{lj}## is: $$ R^i_{jkl}\equiv\partial_k L^i_{jl}-\partial_l L^i_{jk}+L^h_{jl} L^i_{hk}-L^h_{jk} L^i_{hl}=0 $$ Then from the equation: $$\nabla_b e_a= \Gamma^c_{ab} e_c$$ Using cartesian basis ## e_I...
Back
Top