Gravitational force and acceleration in General Relativity.

In summary: According to General Relativity, the coordinate acceleration (measured by an observer at infinity) is:a= \frac{GM}{r^2}\left(1-\frac{2GM}{rc^2}\right)This is very wrong. You need to start with the metric:ds^2=\alpha dt^2-\frac{1}{\alpha}dr^2-r^2d\phi^2\alpha=1-\frac{2m}{r}From this you construct the Lagrangian:L=\alpha (\frac{dt}{ds})
  • #106
Your equation for the coordinate acceleration gives a value that is larger than the proper acceleration and that is exactly the opposite to the generally perceived wisdom.
 
Physics news on Phys.org
  • #107
kev said:
Your equation for the coordinate acceleration gives a value that is larger than the proper acceleration .

You don't know what you are doing. It is really simple:

[tex]a_{proper}=\frac{d\Phi}{dr'}=\frac{d\Phi}{dr}*\frac{dr}{dr'}[/tex]
[tex]a_{coordinate}=\frac{d\Phi}{dr}[/tex]

If [tex]\frac{dr}{dr'}=\sqrt{1-\frac{2GM}{rc^2}[/tex] by what factor do the two accelerations differ?

and that is exactly the opposite to the generally perceived wisdom

LOL
 
  • #108
kev said:
Are you saying that the proper acceleration is:

[tex]\frac{d\Phi}{dr'} = \frac{d\Phi}{dr/\sqrt{(1-2M/R)}} = \frac{GM}{r^2}\left(1-\frac{2GM}{R}\right)^{-1/2} [/tex]

(assuming [tex]dr/dr' = \sqrt{(1-2M/R)} [/tex] where dr/dr' is the ratio of coordinate length to local (proper) length in the radial direction)?

You still have a mistake:

[tex]\frac{d\Phi}{dr'} = \frac{d\Phi}{dr/\sqrt{(1-2M/r)}} = \frac{GM}{r^2}\left(1-\frac{2GM}{rc^2}\right)^{-1/2} [/tex]

You can't get anything right.
 
Last edited:
  • #109
Man, tough crowd. LOL!
 
  • #110
starthaus said:
You still have a mistake:

[tex]\frac{d\Phi}{dr'} = \frac{d\Phi}{dr/\sqrt{(1-2M/r)}} = \frac{GM}{r^2}\left(1-\frac{2GM}{rc^2}\right)^{-1/2} [/tex]

You can't get anything right.

LOL. Great coming from someone who made a complete hash of their calculations in this thread https://www.physicsforums.com/showthread.php?t=403978&page=5 and when I clearly pointed out your errors in #71 you continued arguing the mistake was mine until finally admitting your error in #81.

There is no error in my calculation above, unless you are making a petty reference to inconsistent use of explicitly stating G and sometimes using G=1. In that case you have made the a similar error of explicitly stating c^2 in the final term and using implied c^2 =1 in the middle term. We are all big boys here and occasionally lapse between using explicit c or G and sometimes using units of G=c=1. Get over it. We all do it. I suspect this is yet another red herring to distract us from a major blunder in your work.
 
Last edited:
  • #111
kev said:
There is no error in my calculation above.

Of course not, only about several errors still linger around, it's excusable, you just learned how to do the correct calculations. Unfortunately, you still have no clue what you have calculated even after I took you through the derivation step by step.
 
  • #112
kev said:
LOL. Great coming from someone who made a complete hash of their calculations in this thread https://www.physicsforums.com/showthread.php?t=403978&page=5 and when I clearly pointed out your errors in #71 you continued arguing the mistake was mine until finally admitting your error in #81.

There is no error in my calculation above, unless you are making a petty reference to inconsistent use of explicitly stating G and sometimes using G=1. In that case you have made the a similar error of explicitly stating c^2 in the final term and using implied c^2 =1 in the middle term. We are all big boys here and occasionally lapse between using explicit c or G and sometimes using units of G=c=1. Get over it. We all do it. I suspect this is yet another red herring to distract us from a major blunder in your work.

Nah, you are mixing [tex]R[/tex] with [tex]r[/tex]. There is no [tex]R[/tex]. But worse, you have no clue what you have calculated, even after I took you through the derivation step by step.
 
  • #113
starthaus said:
Nah, you are mixing [tex]R[/tex] with [tex]r[/tex]. There is no [tex]R[/tex]. But worse, you have no clue what you have calculated, even after I took you through the derivation step by step.

Yep, I was right, you were being petty. Sorry for inadvertently using capital R instead of lowercase r. FoR a minute theRe, I thought you had something impoRtant to say.

The important issue is why your equation for coordinate acceleration is so wrong.
 
  • #114
kev said:
Assuming Schwarzschild geometry (uncharged non-rotating body) and defining the gravitational gamma factor as:

[tex] \gamma = \frac{1}{\sqrt{1-\frac{2GM}{rc^2}}} [/tex][tex]a_0 = a \gamma^3[/tex]

Nope, this is just as wrong as the rest. It is only true for [tex]\gamma(v)=\frac{1}{\sqrt{1-(v/c)^2}}[/tex] (and only when [tex]\vec{v}[/tex] and [tex]\vec{a}[/tex] are co-linear)There is no reason for it to be true for

[tex] \gamma = \frac{1}{\sqrt{1-\frac{2GM}{rc^2}}} [/tex]

Numerology and physics are two very different fields.
 
  • #115
kev said:
The important issue is why your equation for coordinate acceleration is so wrong.

Because you don't know what you are doing, even after someone leads you through the derivation step by step.
You really need to take a break from numerology and reflect on what I explained to you in posts 107 and 114.
 
  • #116
starthaus said:
kev said:
Your equation for the coordinate acceleration gives a value that is larger than the proper acceleration .
You don't know what you are doing. It is really simple:

[tex]a_{proper}=\frac{d\Phi}{dr'}=\frac{d\Phi}{dr}*\frac{dr}{dr'}[/tex]
[tex]a_{coordinate}=\frac{d\Phi}{dr}[/tex]

If [tex]\frac{dr}{dr'}=\sqrt{1-\frac{2GM}{rc^2}[/tex] by what factor do the two accelerations differ?
So are you claiming coordinate acceleration is greater than proper acceleration or not? Your response here seems to agree with kev's statement.
 
  • #117
Al68 said:
So are you claiming coordinate acceleration is greater than proper acceleration or not? Your response here seems to agree with kev's statement.

That's what the derivation shows.
 
  • #118
starthaus said:
Because you don't know what you are doing, even after someone leads you through the derivation step by step.
You really need to take a break from numerology and reflect on what I explained to you in posts 107 and 114.

I can do a derivation based on the exact Schwarzschild solution that gives equations that agree exactly with equations in #1 and disagree significantly with your conclusion based on a field aproximation.
 
  • #119
Why don' you two guys do the following: You must arrive to Newtons law of Gravitation in the weak field also. Suppose two spherical masses of 1 and 2 kg at a distance of 1m.For example made of iron with density 7.9 g/cm^3.
Write on the left hand side :
m1 = 2 kg
m2 = 1 kg
r1 = 1m
G = 6,67 *10^-11*m^3/(kg*s^2)

F = G* m1*m2/r^2

F = 1.334*10^-10 N

Now write down on the right hand side every single step you need to calculate to get this single number in General relativity including conversions for G = c = 1 etc. So this would settle this discussion (hopefully!)
 
Last edited:
  • #120
kev said:
I can do a derivation based on the exact Schwarzschild solution that gives equations that agree exactly with equations in #1 and disagree significantly with your conclusion based on a field aproximation.
Starthaus, I must say that this is the part that troubles me most about the potential approach. It is indeed more elegant, but this field approximation step makes it very suspicious to me. Particularly in light of the results for the rotating reference frame where we found that obtaining the correct answer with the potential approach depended critically on whether you were using a strong-field approximation or a weak-field approximation. We should have been able to tell, from the beginning, that we needed to use the strong-field approximation.

Can you explain more about that step? When is the weak field approximation safe, under what conditions is the strong field approximation required, and when would even the strong field approximation fail? Without that information it seems that we have to do the brute-force approach anyway, just to check if the field approximations were good.
 
  • #121
DaleSpam said:
Starthaus, I must say that this is the part that troubles me most about the potential approach. It is indeed more elegant, but this field approximation step makes it very suspicious to me. Particularly in light of the results for the rotating reference frame where we found that obtaining the correct answer with the potential approach depended critically on whether you were using a strong-field approximation or a weak-field approximation. We should have been able to tell, from the beginning, that we needed to use the strong-field approximation.

First things first.

1. Do you agree that post 1 of the OP is a mess?
2. That the results are put in by hand, none of them is derived? kev claims that he sent the derivation to you and DrGreg, is this true? (the derivation cannot be valid given the mistakes in the post).
3. That the only correct formula is the one for proper acceleration [itex]a_0=\frac{GM}{r^2}(1-\frac{GM}{rc^2})^{-1/2}[/itex]?.
5. Since kev did not derive the above for proper acceleration until very late under my guidance, it isn't clear if he din't simply lucked out in the OP.

Now, to the justification of the derivation based on the strong field:

6. Rindler gives the justification for this approach in chapter 9.6.
7. Rindler uses it for deriving proper acceleration for rotation in 9.7.
8. Rindler uses it for deriving proper acceleration in a gravitational field in 11.2.
9. What troubled me was something different, the fact that he uses the approach by equating the strong field metric

[tex]ds^2=e^{2 \Phi/c^2} dt^2-...[/tex]

with the Newtonian approximation for the weak field [tex]ds^2=(1-\frac{2m}{r})dt^2-...[/tex]. I wrote to him about this (in conjunction to criticizing his circular derivation of the equations of accelerated motion in 3.7). He (Rindler) got very defensive and said (textually) that he prefers simpler proofs to the more rigorous ones, so I let the issue drop.

Can you explain more about that step? When is the weak field approximation safe,

See point 9 above. Both fields are required.

under what conditions is the strong field approximation required, and when would even the strong field approximation fail?

Both fields are required (see point 9 , above). To my knowledge, Rindler's approach works for all the examples.
Without that information it seems that we have to do the brute-force approach anyway, just to check if the field approximations were good.

Not necessarily. Rindler's approach gives the correct answers. I had to resort to this simpler approach because there was no hope in teaching kev how to approach the problem from the Lagrangian angle. The other approach that always works is the covariant derivative. So, we have three different approaches that work.
For the rotation case we have the method described in Nikolic's (Demystifier) paper. I much prefer that approach as I demonstrated in the attempt to teach kev. So, we really have 4 different ways to solving the rotating motion. There are more, there is an excellent chapter on the equations of rotating frames in Moller's "Theory of RElativity".
 
Last edited:
  • #122
kev said:
I can do a derivation based on the exact Schwarzschild solution that gives equations that agree exactly with equations in #1

Really, does this include correcting the blunder about [tex]a_0=a\gamma^3[/tex]? I just taught you how to do the derivations starting from the Schwarzschild solution and you are all ready to do it? Then, let's see it. For once, do a derivation rather than put expressions in by hand.
and disagree significantly with your conclusion based on a field aproximation.

Let's see it. Drop the bluster and do the derivations. Once you publish your results, I suggest that you contact Rindler.
 
Last edited:
  • #125
starthaus said:
First things first.
1. Do you agree that post 1 of the OP is a mess?
Post 1 is just a collection of statements and assumptions requesting input from other members of this forum. Whether it is a mess or not, depends on whether the statements are true or not and that remains to be determined.
starthaus said:
2. That the results are put in by hand, none of them is derived?
I never claimed that #1 was a formal derivation. (See accusation 1.)
starthaus said:
3. That the only correct formula is the one for proper acceleration [itex]a_0=\frac{GM}{r^2}(1-\frac{GM}{rc^2})^{-1/2}[/itex]?.
We agree on the equation for proper acceleration. I disgree about your claim that the rest of the equations are wrong.
starthaus said:
5. Since kev did not derive the above for proper acceleration until very late under my guidance, it isn't clear if he din't simply lucked out in the OP.
If the statement's in #1 are correct, it does not really matter how I obtained them, but for the record they are substantially based on the mathpages relativity website which is usually fairly reliable and rigorous.
starthaus said:
Now, to the justification of the derivation based on the strong field:

6. Rindler gives the justification for this approach in chapter 9.6.
7. Rindler uses it for deriving proper acceleration for rotation in 9.7.
8. Rindler uses it for deriving proper acceleration in a gravitational field in 11.2.
From the above it is clear that Rindler has not directly derived the coordinate gravitational acceleration and so the incorrect equation you have obtained for coordinate acceleration, is your own hash up rather than Rindler's.

starthaus said:
9. What troubled me was something different, the fact that he uses the approach by equating the strong field metric

[tex]ds^2=e^{2 \Phi/c^2} dt^2-...[/tex]

with the Newtonian approximation for the weak field [tex]ds^2=(1-\frac{2m}{r})dt^2-...[/tex]. I wrote to him about this (in conjunction to criticizing his circular derivation of the equations of accelerated motion in 3.7). He (Rindler) got very defensive and said (textually) that he prefers simpler proofs to the more rigorous ones, so I let the issue drop.

Allow me to clear up some of your misconceptions. The metric [tex]ds^2=(1-\frac{2m}{r})dt^2-...[/tex] is not an approximation. It is the exact Schwarzschild solution (with units of c=1). The weak field approximation is the assumption that the gravitational potential is given by [tex] \Phi = GM/r [/tex] (same as the Newtonian potential) and strong field approximation is the assumption that the potential is given by [tex] \Phi = (1/2)c^2 \log(1-2GM/(rc^2)) [/tex]

The trouble with the field approximations, is that while it fairly clear what the limitations of the weak field are, it is not clear how strong the field has to be before the strong field approximation breaks down (as Dalespam has pointed out).

Here is another misconception of yours:
starthaus said:
That's too bad, if you end up buying Rindler, you will find out that, contrary to your post 1 (and to your incorrect claims above),

[tex]\vec{f} = \frac{GM}{r^2}\left(1-\frac{2GM}{rc^2}\right)^{-1} [/tex]

is indeed the coordinate force per unit mass. You can multiply by [tex]m_0[/tex] all by yourself. Besides, you have now the results derived rather than put in by hand as you did in post 1.
You seem to think that coordinate force is obtained by multiplying the coordinate acceleration by the rest mass. That is simply wrong in both SR and GR. It is very easy to demonstrate that your assumption is false in SR.
starthaus said:
...
Both fields are required (see point 9 , above).
This is just all wrong. See my comments above. It is the strong field approximation of potential substituted into the Schwarzschild metric. It is not a combination of weak and strong field approximations. Does Rindler really give the exact equation you have quoted above for coordinate acceleration or is it just a botched extrapolation of what he said, by you?

P.S. I think it is time to end the feud and if you agree to stop being insulting, demeaning, evasive, inflamatory and confrontational, then I will agree to be civil with you too.
 
Last edited:
  • #126
kev said:
Post 1 is just a collection of statements and assumptions requesting input from other members of this forum. Whether it is a mess or not, depends on whether the statements are true or not and that remains to be determined.

I never claimed that #1 was a formal derivation. (See accusation 1.)

No, it is just a collection of errors put in by hand. That much we have established.
We agree on the equation for proper acceleration. I disgree about your claim that the rest of the equations are wrong.

How about the blunder:

[tex]a_0=a\gamma^3[/tex] ?

If the statement's in #1 are correct,

They aren't, see bove.
it does not really matter how I obtained them,

Sure it does, physics is about deriving results. You are just putting in results by hand. Sometimes you hit but most often you miss.
 
  • #127
kev said:
Allow me to clear up some of your misconceptions. The metric [tex]ds^2=(1-\frac{2m}{r})dt^2-...[/tex] is not an approximation. It is the exact Schwarzschild solution (with units of c=1). The weak field approximation is the assumption that the gravitational potential is given by [tex] \Phi = GM/r [/tex] (same as the Newtonian potential) and strong field approximation is the assumption that the potential is given by [tex] \Phi = (1/2)c^2 \log(1-2GM/(rc^2)) [/tex]

LOL.

[tex]ds^2=e^{2\Phi/c^2}dt^2-...[/tex]

is approximated in the weak field by:

[tex]ds^2=(1+\frac{2\Phi}{c^2})dt^2-...[/tex]

(simple Taylor expansion)

Now, substitute [tex]\Phi=-mc^2/r[/tex] above.
 
  • #128
starthaus said:
[tex]ds^2=e^{2\Phi/c^2}dt^2-...[/tex]

is approximated in the weak field by:

[tex]ds^2=(1+\frac{2\Phi}{c^2})dt^2-...[/tex]

(simple Taylor expansion)

Now, substitute [tex]\Phi=-mc^2/r[/tex] above.

I concede that

[tex]ds^2= (1+\frac{2\Phi}{c^2})dt^2 - (1+\frac{2\Phi}{c^2})^{-1}dr^2 - ...[/tex]

approximates to the Newtonian gravitational field when (1-GM/r) is approximately equal to unity and when [tex]\Phi=-GMc^2/r[/tex] is substituted. It is also independent of the value or units used for c.

I was thrown off course by this contradictory post https://www.physicsforums.com/showpost.php?p=2730381&postcount=81 by you in another thread:
starthaus said:
Start with the Schwarzschild solution in the weak field approximation:

[tex](cd\tau)^2=(1+\frac{2\Phi}{c^2})(cdt)^2+(1+\frac{2\Phi}{c^2})^{-1}(dr)^2+...[/tex]
...
...
At the Earth surface :

[tex]\Phi_1=-\frac{GM}{R}[/tex]

The result of the statements you make in the above quote is that the weak field approximation is:

[tex](cd\tau)^2=(1-\frac{2GM}{Rc^2})(cdt)^2+(1-\frac{2GM}{Rc^2})^{-1}(dr)^2+...[/tex]

which is not correct.

At least the equation for the weak field you give in this thread does approximate to the Newtonian one, as shown here:

[tex]d\tau^2= (1-\frac{2GM}{r})dt^2 - (1+\frac{2GM}{r})^{-1}dr^2 [/tex]

Assume radial motion only and assume dr/dt is small so that
[tex]d\tau \approx dt[/tex]

[tex]dt^2 = (1-\frac{2GM}{r})dt^2 - (1+\frac{2GM}{r})^{-1}dr^2 [/tex]

[tex]\frac{dr}{dt} = \sqrt{ \frac{4G^2M^2}{r^2}-\frac{2GM}{r} }[/tex]

[tex]\frac{d^2r}{dt^2} = \frac{d(dr/dt)}{dt} = \frac{d(dr/dt)}{dr}\frac{dr}{dt}= \frac{GM}{r^2}(1-2GM/r) - \frac{2G^2M^2}{r^3} [/tex]

In the weak field (1-2GM/r) is approximately unity so:

[tex]\frac{d^2r}{dt^2} = \frac{GM}{r^2}(1-2GM/r) - \frac{2G^2M^2}{r^3} \approx \frac{GM}{r^2} - \frac{2G^2M^2}{r^3} \approx \frac{GM}{r^2} (1-2GM/r) \approx \frac{GM}{r^2}[/tex]

starthaus said:
...
How about the blunder:

[tex]a_0=a\gamma^3[/tex] ?

That is not a mistake. I still maintain the equations in #1 are correct. I can admit and correct my mistakes as I have done in this post. I learn and move on. You should try doing the same. It makes for much better progress all round.

Anyway, the weak field approximation is not important to this thread. We are more interested in where the fields are much stronger.

However, this post does establish that your derivation uses the strong field potential approximation substituted in the Schwarzschild metric (See post #95 https://www.physicsforums.com/showpost.php?p=2731106&postcount=95) and not a combination of the strong and weak field approximations as you claim.
 
Last edited:
  • #129
kev said:
I concede that

[tex]ds^2= (1+\frac{2\Phi}{c^2})dt^2 - (1+\frac{2\Phi}{c^2})^{-1}dr^2 - ...[/tex]

approximates to the Newtonian gravitational field when (1-GM/r) is approximately equal to unity and when [tex]\Phi=-GMc^2/r[/tex] is substituted. It is also independent of the value or units used for c.

I was thrown off course by this contradictory post https://www.physicsforums.com/showpost.php?p=2730381&postcount=81 by you in another thread:


The result of the statements you make in the above quote is that the weak field approximation is:

[tex](cd\tau)^2=(1-\frac{2GM}{Rc^2})(cdt)^2+(1-\frac{2GM}{Rc^2})^{-1}(dr)^2+...[/tex]

which is not correct.

THat was another thread. That was a typo. That has no bearing on the final result, a result that to this day you still do not understand since you don't know that [tex]d\tau_1/d\tau_2[/tex] is not the ratio of frquencies. You really need to read on this subject, I can recommend pages 346-7 in Moller's "Theory of Relativity".

For the time being, stick to the subject.


Assume radial motion only and assume dr/dt is small so that
[tex]d\tau \approx dt[/tex]

But this is a horrible hack that is not true. You know very well the value of the raitio [tex]d\tau/dt=\sqrt{1-2m/r}[/tex]





That is not a mistake. I still maintain the equations in #1 are correct. I can admit and correct my mistakes as I have done in this post. I learn and move on. You should try doing the same. It makes for much better progress all round.

LOL. You are plugging in formulas that don't apply because you don't know how they were derived
 
  • #130
kev said:
...
Assume radial motion only and assume dr/dt is small so that
[tex]d\tau \approx dt[/tex]
starthaus said:
But this is a horrible hack that is not true. You know very well the value of the raitio [tex]d\tau/dt=\sqrt{1-2m/r}[/tex]

I know that, but approximation solutions always involve "horrible hacks" and should be used with caution. If you do not like hacks, then you should probably avoid derivations based on the weak or strong field approximations. They start with failure built in.

There is no notion of any distinction between proper time and coordinate time in Newtonian equations and from that point of view all Newtonian equations are a horrible hack as you put it.

The formula you base you derivation on, of

[tex]\vec{F}=-grad\Phi[/tex]

is true in Newtonian physics and happens to also be true for proper acceleration and force in relativity, but it is not true for coordinate acceleration without modification.
 
Last edited:
  • #131
This is a question for Dalespam.

I know you are an expert in 4 vectors and have already obtained a correct solution for the proper acceleration. Would you be able to obtain a solution for the coordinate gravitational acceleration for us, using those methods? It would be great help to this thread.
 
  • #132
pepeherborn said:
Why don' you two guys do the following: You must arrive to Newtons law of Gravitation in the weak field also. Suppose two spherical masses of 1 and 2 kg at a distance of 1m.For example made of iron with density 7.9 g/cm^3.
Write on the left hand side :
m1 = 2 kg
m2 = 1 kg
r1 = 1m
G = 6,67 *10^-11*m^3/(kg*s^2)

F = G* m1*m2/r^2

F = 1.334*10^-10 N

Now write down on the right hand side every single step you need to calculate to get this single number in General relativity including conversions for G = c = 1 etc. So this would settle this discussion (hopefully!)

Hi pepeherborn and welcome to PF :smile:. Sorry, didn't mean to ignore you, but we were caught up in a rather heated argument. (You might have noticed LOL). Unfortunately the solution is not as simple as you suggest and there are many ways of interpreting mass and distance (r) in General Relativity and that is probably why we often appear to be in conflict here. For simplicity it is also probably best to use to use one mass that is significantly larger than the other. If we ever come to an agreement in this thread, we may be able to answer you question in your terms.
 
  • #133
kev said:
kev said:
...
Assume radial motion only and assume dr/dt is small so that
[tex]d\tau \approx dt[/tex]I know that, but approximation solutions always involve "horrible hacks" and should be used with caution. If you do not like hacks, then you should probably avoid derivations based on the weak or strong field approximations. They start with failure built in.

Physics and numerology are different disciplines. What you are doing isn't physics. I suggest that you invest in some good books , I recommended quite a few.(Rindler,Moller)
There is no notion of any distinction between proper time and coordinate time in Newtonian equations and from that point of view all Newtonian equations are a horrible hack as you put it.

While true this is not a valid argument for hacking GR. There are perfectly correct solutions to the problem that don't involve hacks. You have been given at least two in this thread, I don't see the point in your continuing to produce hacks.
The formula you base you derivation on, of

[tex]\vec{F}=-grad\Phi[/tex]

is true in Newtonian physics and happens to also be true for proper acceleration and force in relativity, but it is not true for coordinate acceleration without modification.

Are you reading these things in any book or are you simply making them up?
 
  • #134
kev said:
Assume radial motion only and assume dr/dt is small so that
[tex]d\tau \approx dt[/tex]

[tex]dt^2 = (1-\frac{2GM}{r})dt^2 - (1+\frac{2GM}{r})^{-1}dr^2 [/tex]

[tex]\frac{dr}{dt} = \sqrt{ \frac{4G^2M^2}{r^2}-\frac{2GM}{r} }[/tex]

[tex]\frac{d^2r}{dt^2} = \frac{d(dr/dt)}{dt} = \frac{d(dr/dt)}{dr}\frac{dr}{dt}= \frac{GM}{r^2}(1-2GM/r) - \frac{2G^2M^2}{r^3} [/tex]

In the weak field (1-2GM/r) is approximately unity so:

[tex]\frac{d^2r}{dt^2} = \frac{GM}{r^2}(1-2GM/r) - \frac{2G^2M^2}{r^3} \approx \frac{GM}{r^2} - \frac{2G^2M^2}{r^3} \approx \frac{GM}{r^2} (1-2GM/r) \approx \frac{GM}{r^2}[/tex]

If the above were correct (it isn't) and since we have already established that :

[tex]a_0= \frac{GM}{r^2}(1-2GM/(rc^2))^{-1/2}[/tex]

it becomes clearly apparent that you can't have

[tex]a_0=a\gamma^3[/tex]

So, is the your coordinate acceleration expression false or is the relationship between proper and coordinate acceleration false?
Hint: they are both false.
 
Last edited:
  • #135
starthaus said:
If the above were correct (it isn't) and since we have already established that :

[tex]a_0= \frac{GM}{r^2}(1-2GM/(rc^2))^{-1/2}[/tex]

it becomes clearly apparent that you can't have

[tex]a_0=a\gamma^3[/tex]

So, is the your coordinate acceleration expression false or is the relationship between proper and coordinate acceleration false?
Hint: they are both false.

Wrong again. As I said a few posts before, if you start a derivation from an approximation, you have failure built in from the start, if you try to come to any general conclusions outside the scope of the approximation.


This is the (long awaited) correct way to derive the coordinate acceleration:

(It is the substance of the derivation I sent to Dalespam and DrGreg, but tidied up and debugged. Thanks are due to Dalespam for looking over it and picking up some of the original typos. )

Derivation of the coordinate acceleration of in Schwarzschild spacetime.

Starting with the radial Schwarzschild metric:

[tex]ds^2=\alpha dt^2-\frac{1}{\alpha}dr^2[/tex]

[tex]\alpha=1-\frac{2M}{r}[/tex]

Solve for dr/dt:

[tex] \frac{dr}{dt} \;=\; \sqrt{\alpha^2 - \alpha\left( \frac{ds}{dt}\right )^2}[/tex]

From the Lagragian and the equations of motion we have the conserved quantity k defined as:

[tex]\frac{dt}{ds} \;=\; \frac{k}{\alpha}[/tex]

(See #2 by Starthaus https://www.physicsforums.com/showthread.php?p=2710570#post2710570 )

(Note that the equations of motion are only valid for a particle in freefall.)

This value [tex](k/\alpha)[/tex] is substituted for dt/ds in the equation above it.

[tex] \frac{dr}{dt} \;=\; \sqrt{\alpha^2 - \frac{\alpha^3}{k^2}}[/tex]

Now dr/dt is differentiated again:

[tex]\frac{d^2r}{dt^2} \;=\; \frac{d(dr/dt)}{dr}\frac{dr}{dt}\;=\; -\frac{M}{r^2}\left(\frac{3\alpha^2}{k^2}-2\alpha \right )[/tex]

For a particle falling from infinity, k has the value 1, but for a particle at its apogee, [tex]k = \sqrt{\alpha} .[/tex]

(See https://www.physicsforums.com/showpost.php?p=2714572&postcount=4 for how the value of k is obtained).

Substituting [tex]k=\sqrt{\alpha}[/tex] gives the coordinate acceleration of a particle at its apogee at r as:

[tex]\frac{d^2r}{dt^2} \;=\;-\frac{M}{r^2}\alpha \;=\; -\frac{M}{r^2}\left(1-\frac{2M}{r}\right )[/tex]

Note that the particle is stationary, in the sense that it is at the peak of its trajectory and is in freefall.

If a value of k=1 is used, then the coordinate acceleration of a particle falling from infinity is obtained.

The derivation of the proper acceleration is given in the next post.
 
Last edited:
  • #136
Derivation of the proper acceleration of a particle in Schwarzschild spacetime.

Starting with the radial Schwarzschild metric:

[tex]ds^2=\alpha dt^2-\frac{1}{\alpha}dr^2[/tex]

[tex]\alpha=1-\frac{2M}{r}[/tex]

Solve the metric for dr/ds:

[tex]\frac{dr}{ds} = \sqrt{\alpha^2\frac{dt^2}{ds^2} - \alpha}[/tex]

Substitute [tex]k/\alpha[/tex] for dt/ds: (See previous post).

[tex]\frac{dr}{ds} = \sqrt{k^2 - \alpha}[/tex]

Differentiate again with respect to s:

[tex]\frac{d^2r}{ds^2} = \frac{d(dr/ds)}{dr} \frac{dr}{ds} = -\frac{M}{r^2}[/tex]

For a particle in freefall at r, but stationary or nearly stationary (i.e at the apogee of its trajectory):

[tex]\frac{dr'}{dr} = \frac{1}{\sqrt{1-2M/r}}[/tex]

and the proper acceleration is:

[tex]\frac{d^2r'}{ds^2} = \frac{d(dr)}{ds^2} \frac{dr'}{dr} = \frac{d^2r}{ds} \frac{1}{\sqrt{1-2M/r}} = -\frac{M}{r^2}\frac{1}{ \sqrt{1-2M/r}}[/tex]

(Note that proper acceleration is independent of the k parameter and therefore the proper acceleration of a free falling particle is indepenent of its velocity.)

This is the initial acceleration of a particle released from radial coordinate r as measured by a local observer at r and is equal in magnitude to the proper acceleration of a particle at rest at r, as measured by an accelerometer.
 
Last edited:
  • #137
kev said:
Now dr/dt is differentiated again:

[tex]\frac{d^2r}{dt^2} \;=\; \frac{d(dr/dt))}{dr}\frac{dr}{dt}\;=\; -\frac{m}{r^2}\left(\frac{3\alpha^2}{k^2}-2\alpha \right )[/tex]

For a particle falling from infinity, k has the value 1, but for a particle at rest at a given radial coordinate, [tex]k = \sqrt{\alpha} .[/tex]

(See https://www.physicsforums.com/showpost.php?p=2714572&postcount=4 for how the value of k is obtained).

Substituting [tex]k=\sqrt{\alpha}[/tex] gives the coordinate acceleration of a particle at its apogee at r as:

Ummm, no. You found the calculations "ready made" by a very good guy (prof. Kevin Brown).
But he goofed badly on his page, one of the rare instances that he made a very serious mistake.

Here is the problem : when you (or more exactly he) differentiated wrt [tex]r[/tex] he did it as if [tex]k[/tex] is a constant, i.e. it does not depend of [tex]r[/tex].
Unfortunately he later "forgets" about [tex]k[/tex] being a constant when he makes [tex]k=\sqrt{\alpha}=\sqrt{1-2m/r}[/tex] dependent on [tex]r[/tex]. Even the very best make mistakes and this is a very bad one. So, prof. Brown's proof is invalid.
 
Last edited:
  • #138
kev said:


Substitute [tex]k/\alpha[/tex] for dt/ds: (See previous post).

[tex]\frac{dr}{ds} = \sqrt{k^2 - \alpha}[/tex]
Nope, coupled with [tex]k=\sqrt{\alpha}[/tex] (see your previous post) that would make

[tex]\frac{dr}{ds} = 0[/tex]

This is really bad.

Differentiate again with respect to ds:

[tex]\frac{d^2r}{ds^2} = \frac{d(dr/ds)}{dr} \frac{dr}{ds} = -\frac{m}{r^2}[/tex]

Nope, same mistake as before, you can't forget that [tex]k[/tex] is a function of [tex]r[/tex] and differentiate as if [tex]k[/tex] is a constant. (it isn't).

I think I have Kevin's email somewhere, I used to email him suggestions how to improve his book.
 
  • #139
starthaus said:
Nope, coupled with [tex]k=\sqrt{\alpha}[/tex] (see your previous post) that would make

[tex]\frac{dr}{ds} = 0[/tex]

This is really bad.
Why is that bad? Shouldn't [tex]\frac{dr}{ds} = 0[/tex] be true for a particle at apogee?

Not trying to butt in, just trying to follow along. :smile:
 
  • #140
starthaus said:
Nope, coupled with [tex]k=\sqrt{\alpha}[/tex] (see your previous post) that would make

[tex]\frac{dr}{ds} = 0[/tex]

This is really bad.
Al68 said:
Why is that bad? Shouldn't [tex]\frac{dr}{ds} = 0[/tex] be true for a particle at apogee?

Not trying to butt in, just trying to follow along. :smile:
You are quite right Al. I think the objection that Starthaus is making is that this implies the next step is differentiating zero with respect to s. The clever part of Prof Brown's derivation is that by using the constant k he can get around this and carry out the differentiation. If Starthaus was as good at calculus as he makes out he is, he would realize there is no problem with the step taken by K. Brown. You only have to read Prof Brown's website to realize he very good at calculus and physics.

Here is another way of looking at it.

[tex]x = \sqrt{\frac{1}{2}-\left(1-\frac{2M}{r}\right)}[/tex]

[tex]\frac{dx}{dr} = \frac{M}{r^2\sqrt{2M/r-1/2}}[/tex]

The above is a perfectly valid derivative, even though for a certain value of r, x=0.

The above can rephrased using [tex]k^2=1/2[/tex] and [itex]\alpha = (1-2M/r)[/itex] as:

[tex]x = \sqrt{k^2-\alpha}[/tex]

and as long as it understood that k is a constant and [tex]\alpha[/tex] is a function of the variable r, there is no problem differentiating it.
 
Last edited:

Similar threads

Replies
2
Views
1K
Replies
16
Views
2K
Replies
11
Views
934
Replies
3
Views
690
Back
Top