Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Corrections to Narlikar paper

  1. Aug 23, 2014 #1


    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    Corrections to Narlikar paper

    The paper by Jayant N Narlikar ’Spectral shifts
    in General Relativity’
    in the American Journal of Physics October 1994 is extremely useful in explaining the 'meaning' of cosmological and gravitational redshift, as well as making sense of the difference between peculiar motion and 'motion' arising from the expansion of space.

    However the calculations in the paper have a few technical errors that can be confusing for beginners like me. These calculations correct the invalidity of Narlikar’s proof arising from his use of bases that do not exist, as well as an error in Equation [12] that generates errors in equations 13-17. 20 and 21.

    Equation [11] for the case i=0 should be as follows:
    0 &= \frac{d^2t}{du^2} - \frac{1}{2}g_{kl,0}\frac{dx^k}{du}\frac{dx^l}{du} \\
    &= \frac{d^2t}{du^2} - \frac{1}{2}g_{00,0}(\frac{dt}{du})^2 - \frac{1}{2}g_{11,0}(\frac{dr}{du})^2\text{ [since }\frac{d\theta}{du}\text{ and }\frac{d\phi}{du}\text{ must be zero]}\\
    &= \frac{d^2t}{du^2} - \frac{1}{2}\frac{\partial(\frac{-a(t)^2}{1-kr^2})}{\partial t}(\frac{dr}{du})^2 \text{ [since }g_{00}\text{ is constant at 1] }\\
    &= \frac{d^2t}{du^2} + (\frac{dr}{du})^2 \frac{a(t)a'(t)}{1-kr^2}\ \ \ \ \textbf{[12]}\\
    &= \frac{d^2t}{du^2} + (\frac{dr}{dt})^2 (\frac{dt}{du})^2 \frac{a(t)a'(t)}{1-kr^2}\\
    &= \frac{d^2t}{du^2} + \frac{1-kr^2}{a(t)^2}(\frac{dt}{du})^2 \frac{a(t)a'(t)}{1-kr^2}\text{ [by 10] }\\
    &= \frac{d^2t}{du^2} + (\frac{dt}{du})^2 \frac{a'(t)}{a(t)}\\
    &= \frac{1}{a(t)} \frac{d(a(t)\frac{dt}{du})}{du}\\

    Hence, as [itex]a(t)>0[/itex] we have [itex]0 = \frac{d(a(t)\frac{dt}{du})}{du}[/itex], whence:
    [itex]a(t)\frac{dt}{du} = A[/itex] for some constant [itex]A[/itex].[13]
    \begin{align*}\frac{dr}{du} &= \frac{dr}{dt}\frac{dt}{du}\\
    &= \frac{A}{a(t)} \frac{\sqrt{1-kr^2}}{a(t)}\ \ \text{[by 13 and 10]}\\
    &= \frac{A\sqrt{1-kr^2}}{a(t)^2} \ \ \ \textbf{[14]}\\
    This enables u to be determined as a function of r (as the light moves from S to O) and hence of t. It is convenient and permissible to set u=0 at S and u=1 at O.
    Integrating [13] we obtain:
    [itex]A = \int_0^1A\ du = \int_{t(O)}^{t(S)}a(t)dt[/itex] [15]
    Let the tangent vector to the geodesic [itex]\Lambda[/itex] at [itex]\lambda(u)\text{ be }U(u)[/itex]. In the comoving coordinates this has components [itex][\frac{dt}{du}, \frac{dr}{du}, \frac{d\theta}{du}, \frac{d\phi}{du}][/itex]. The last two are zero and the first two are given by [13] and [14], hence the components are:
    [itex]U(u) = [\frac{A}{a(t)}, \frac{A\sqrt{1-kr^2}}{a(t)^2}, 0, 0][/itex] [16a]
    and at O we have:
    [itex]U(1) = [\frac{A}{a_O}, \frac{A\sqrt{1-kr_O^2}}{a_O^2}, 0, 0][/itex] [16b]
    where [itex]r_O[/itex] is the radial comoving coordinate of [itex]O[/itex]. This is the first equation that assumes a specific centre ([itex]S[/itex]) for the FLRW system.

    Let the vector [itex]\vec{V}_S[/itex] parallel transported from [itex]S[/itex] to [itex]\lambda(u)[/itex] be [itex]\vec{V}^*(u)[/itex], and let us denote [itex]\vec{V}^*(1)[/itex] by [itex]\vec{\bar{V}}_S[/itex]. Then, as parallel transport preserves magnitude and direction, both [itex]\lvert\lvert\vec{V}^*(u)\rvert\rvert = g(\vec{V}^*(u),\vec{V}^*(u))\text{ and }g(\vec{V}^*(u),\vec{U}(u))[/itex] must be constant over u. [17]

    Hence [itex]\lvert\lvert\vec{\bar{V}}_S\rvert\rvert = \lvert\lvert\vec{V}^*(1)\rvert\rvert = \lvert\lvert\vec{V}^*(0)\rvert\rvert = \lvert\lvert\vec{V}_S\rvert\rvert = 1[/itex] [by 9a]. [18a]

    And [itex]g(\vec{\bar{V}}_S,\vec{U}(1)) = g(\vec{V}^*(1),\vec{U}(1)) = g(V^*(0),U(0)) = g(\vec{V}_S,\vec{U}(0))[/itex] [18b]

    Now we must have [itex]V^{*2}(u) = V^{*3}(u) = 0[/itex] because otherwise the parallel transportation establishes a preferred direction in space, circumferential to [itex]S[/itex], which contradicts the isotropy assumption. [This is a bit hand-wavy. Seek to make it more rigorous]

    From [18a] we get:
    [itex]1 = \lVert\vec{\bar{V}}_S\rVert = g(O)(\vec{\bar{V}}_S, \vec{\bar{V}}_S) = g_{ik}(O)\bar{V}^i_S\bar{V}^k_S= g_{00}(O)(\bar{V}^0_S)^2 + g_{11}(O)(\bar{V}^0_S)^2\\[/itex]
    [as off-diagonal elements of [itex]g[/itex] are zero everywhere in spacetime [from 7] and
    [itex]\bar{V}^2_S = \bar{V}^3_S = V^{*2}(1) = V^{*3}(1) = 0]

    = (\bar{V}^0_S)^2 - \frac{a_O^2}{1-kr_O^2} (\bar{V}^1_S)^2[/itex]

    Hence [itex](\bar{V}^0_S)^2 - \frac{a_O^2}{1-kr_O^2} (\bar{V}^1_S)^2 = 1[/itex] [19]
    From [18b] we get [itex]g_O(\vec{\bar{V}}_S, \vec{U}(1)) = g_S(\vec{V}_S,\vec{U}(0))[/itex]
    hence [itex]g_{ik}(O) \bar{V}^i_S U^k(1) = g_{ik}(S) V^i_S U^k(0)[/itex]

    The left-hand side is:
    [itex] = g_{00}(O)\bar{V}^0_S\ U^0(1) + g_{11}(O)\bar{V}^1_S\ U^1(1)[/itex]
    [since [itex]\bar{V}^2_S = \bar{V}^3_S = 0[/itex]]
    [itex] = \bar{V}^0_S\ U^0(1) - \frac{a_O^2}{1-kr_O^2} \bar{V}^1_S\ U^1(1)
    = \bar{V}^0_S \frac{A}{a_O} - \frac{a_O^2}{1-kr_O^2}\frac{A\sqrt{1-kr_O^2}}{a_O^2} \bar{V}^1_S[/itex] [by 16b]
    [itex] = A(\frac{\bar{V}^0_S}{a_O} - \frac{\bar{V}^1_S}{\sqrt{1-kr_O^2}}) [/itex]

    The right-hand side is (note that, since we are operating in [itex]T_SM[/itex] here, we have to switch to [itex]FO[/itex], the FLRW basis centred at [itex]O[/itex], in order for the components to be well-defined):

    [itex] = {}_{FO}g_{00}(S) {}_{FO}V^0_S {}_{FO}U(0)^0 + {}_{FO}g_{11}(S) {}_{FO}V^1_S {}_{FO}U(0)^1[/itex]

    [since [itex]U(0)^2=U(0)^3=0[/itex] because the light ray is radial. \hl{This is also a bit hand-wavy}]
    [itex]\ = {}_{FO}g_{00}(S) {}_{FO}U(0)^0 [/itex]

    [since the FLRW spatial coordinates of [itex]S[/itex] are constant, so [itex]\vec{V}_S = [1,0,0,0][/itex] in the [itex]FO[/itex] frame]
    [itex]\ = {}_{FO}U(0)^0 \text{ [by 7a]}
    = \frac{A}{a_S}[/itex] [by 16a]

    Hence, equating the right and left sides we get: [itex]\frac{\bar{V}^0_S}{a_O} - \frac{\bar{V}^1_S}{\sqrt{1-kr_O^2}} = \frac{1}{a_S}[/itex] [20]

    Next we note that the four-velocity of the observer has components [itex][1,0,0,0][/itex] in any FLRW frame (because the observer has zero spatial coordinate velocity in that frame) and also in the observer's inertial frame at O. Hence the time basis vectors of the [itex]MO[/itex] and [itex]FS[/itex] frames must be identical: [itex]\vec{e}_0(O) = {}_{MO}\vec{e}_0[/itex].
    Now [itex]\vec{\bar{V}}_S = {}_{MO}\bar{V}^i_S \ {}_{MO}\vec{e}_i = \bar{V}^i_S\ \vec{e}_i(O) = \bar{V}^0_S \ \vec{e}_0(O) + \frac{a_O^2}{\sqrt{1-kr_O^2}}^1_S\ \vec{e}_1(O)\\ = \bar{V}^0_S\ {}_{MO}\vec{e}_0 + \bar{V}^1_S\ \vec{e}_1(O)[/itex]

    Hence, since both bases are orthogonal, we have [itex]{}_{MO}\bar{V}^0_S = \bar{V}^0_S[/itex] and, by rotating the Lorentz frame appropriately around [itex]O[/itex], we can without loss of generality choose our basis vector [itex]{}_{MO}\vec{e}_1[/itex] so that it aligns with [itex]\vec{e}_1(O)[/itex]. Then, since [itex]\vec{\bar{V}}_S = [\bar{V}_S^0, \bar{V}_S^1, 0, 0][/itex] in the [itex]FS[/itex] basis, we can write [itex]\vec{\bar{V}}_S = [\gamma, \gamma\ V, 0, 0][/itex] in the [itex]MO[/itex] basis, where [itex]{}_{MO}\bar{V}^0_S = \gamma = \frac{1}{\sqrt{1-V^2}} = \bar{V}^0_S[/itex] and [itex]V[/itex] has the same sign as [itex]\bar{V}^1_S[/itex] (because [itex]\gamma[/itex] is positive and the corresponding basis vectors, [itex]{}_{MO}\vec{e}_1[/itex] and [itex]\vec{e}_1(O)[/itex], point in the same direction).

    Hence [itex]\lvert\lvert\vec{\bar{V}}_S\rvert\rvert = g(\vec{\bar{V}}_S,\vec{\bar{V}}_S) = \gamma^2 - (\gamma V)^2 = \gamma^2 - \frac{a_O^2}{1-kr_O^2} (\bar{V}^1_S)^2[/itex]
    where we calculate the magnitude in the [itex]MO[/itex] and [itex]FS[/itex] bases and equate the results.\\
    Hence [itex](\gamma V)^2 = \frac{a_O^2}{1-kr_O^2} (\bar{V}^1_S)^2[/itex] and so
    [itex]\gamma V = \frac{a_O}{\sqrt{1-kr_O^2}} \bar{V}^1_S[/itex] [21]
    where the sign is positive because [itex]\gamma[/itex] is positive and [itex]V[/itex] and [itex]\bar{V}^1_S[/itex] have the same sign.

    Now the red shift is given by:
    [itex]1+z = \frac{a_O}{a_S}[/itex] [8]
    [itex] = a_O(\frac{\bar{V}^0_S}{a_O} - \frac{\bar{V}^1_S}{\sqrt{1-kr_O^2}})[/itex] [by 20]
    [itex] = \gamma - \frac{a_O}{\sqrt{1-kr_O^2}}\ \bar{V}^1_S = \gamma-\gamma V[/itex] [by 21]
    [itex] = \gamma(1-V) = \frac{1-V}{\sqrt{1-V^2}} =
    \sqrt{\frac{(1-V)^2}{(1-V)(1+V)} }\\
    = \sqrt{\frac{1 - V}{1+V}}\ \ \ \ \ [/itex] [22]
    This is the reciprocal of the formula for a Doppler shift within a Lorentz frame.

  2. jcsd
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Can you offer guidance or do you also need help?
Draft saved Draft deleted