- #71
Rasalhague
- 1,387
- 2
kev said:That should work now.
It does. It gives the same answer as x*Sqrt[1 - bx^2], and is unaffected by changes in by.
kev said:That should work now.
How does your equation compare numerically with mine?Rasalhague said:For the effect of a general boost on a general space vector, I get:
[tex]\textbf{r}'=\textbf{r}-\frac{\textbf{r}\cdot \pmb{\beta}}{\beta^2}\pmb{\beta}\left ( 1-\frac{1}{\gamma(\beta)} \right )[/tex]
And for the angle,
[tex]\cos(\textbf{r},\textbf{r}')=\frac{r^2 \beta^2-(\textbf{r}\cdot \pmb{\beta})^2}{r\beta^2\sqrt{r^2-(\textbf{r}\cdot \pmb{\beta})^2}} \left ( 1-\frac{1}{\gamma(\beta)} \right )[/tex]
[tex]=\frac{1-\cos^2(\textbf{r},\pmb{\beta})}{1-\beta^2 \cos(\textbf{r},\pmb{\beta})}[/tex]
where, for example, [itex]\cos(\textbf{r},\textbf{r}')[/itex] is the angle between the original space vector [itex]\textbf{r}[/itex] and its boosted counterpart [itex]\textbf{r}'[/itex].
I haven't been following this thread and I'm reading this post out-of-context, but if you're using standard notation, I believe you should have [itex]\gamma[/itex] instead of [itex]1 / \gamma[/itex]. See Lorentz transformation - Matrix form, final equation in section.Rasalhague said:For the effect of a general boost on a general space vector, I get:
[tex]\textbf{r}'=\textbf{r}-\frac{\textbf{r}\cdot \pmb{\beta}}{\beta^2}\pmb{\beta}\left ( 1-\frac{1}{\gamma(\beta)} \right )[/tex]
kev said:How does your equation compare numerically with mine?
In[1]:= th + ArcTan[Cot[th] Sqrt[1 - (b^2)]] - Pi/2
Out[1]:= -0.489976
In[2]:= ArcCos[(1 - Cos[th]^2)/(1 - b^2 Cos[th])]
Out[2]:= 0.817476
Rasalhague said:Code:In[1]:= th + ArcTan[Cot[th] Sqrt[1 - (b^2)]] - Pi/2 Out[1]:= -0.489976 In[2]:= ArcCos[(1 - Cos[th]^2)/(1 - b^2 Cos[th])] Out[2]:= 0.817476
This is with th = ArcTan[.3/.9] = 0.321751. Did I type it right? The angle should be getting bigger, shouldn't it?
kev said:The rotation [itex](\phi)[/itex] of the rod relative to the x' axis in the S' frame is given by:
[tex]\phi = \theta + tan^{-1}\left( \frac{\cos\theta}{\sin\theta } \, \sqrt{(1-\beta_x^2-\beta_y^2)} \right)-\frac{\pi}{2} [/tex]
This rotation is generally away from the line of motion and the orientation of the rod tends towards being orthogonal to the motion at velocities approaching c.
Now that we have the length of the rod [itex]\| L' \|[/itex] in S' and its angle with respect to the x' axis, it is easy enough to work out the projection (L'x) of the rod's length onto the x' axis as:
[tex]L'_x = \| L' \| \cos\phi[/tex]
The angle of the rod with respect to to the line of motion is [itex]-\theta[/itex] in frame S and [itex]-\theta+\phi[/itex] in frame S'. [itex]\phi[/itex] has the opposite sign to [itex]\ \theta[/itex] and so the angle wrt the line of motion is greater in S' than in S.
kev said:Your transcription of my equation looks OK. That equation is for phi which I have defined as the rotation of the rod away from the x' axis. If you look at the bit I have highlighted in the quote from #63 below you will see that I given the total rotation of the rod away from the line of motion as -theta +phi so that you should get something like -0.321751-0.489976 which is a clockwise rotation and I think if done to enough decimal places should be close to your result give or take the sign. Which way the rod rotates depends on which side of the line of motion the rod is on, but the rotation is always away from the line of motion towards being perpendicular to the line of motion, (assuming one end of the rod is on the line of motion).
In the limit as theta goes to zero, the atan function in my expression goes to Pi/2 meaning that there is zero rotation when the by=0. The mathematical software strugglesat this extreme.Rasalhague said:Ah, I see. Curiously, for this example, they're close, but not exactly the same:
ArcCos[(1 - Cos[th]^2)/(1 - b^2 Cos[th])] = 0.817476
ArcTan[Cot[th] Sqrt[1 - (b^2)]] - Pi/2 = ph - th = -0.811726
But in general, they're not the same. When I set by = 0, and all the velocity in the x direction, mine gives pi/2, yours indeterminate, but they should give 0 in that case, shouldn't they? I sometimes get imaginary results with mine, and it doesn't agree with the other equation I had either. Hmm. I'm going to have to leave it for now, and try again tomorrow...
Hi. I think you might have missed a late edit to #76 where I decided that your rotation function should be the same as phi rather than the phi -theta I thought earlier.Rasalhague said:[tex]\cos(\textbf{r},\textbf{r}')=\frac{1-\left ( 1-\gamma^{-1} \right ) \cos^2(\textbf{r},\pmb{\beta})}{\sqrt{1-\beta^2 \cos^2(\textbf{r},\pmb{\beta})}}[/tex]
Well, I picked up a couple of mistakes I made when I tried to test this yesterday, but it's still not agreeing with ArcTan[Cot[th] Sqrt[1 - (b^2)]] - Pi/2.
kev said:Hi. I think you might have missed a late edit to #76 where I decided that your rotation function should be the same as phi rather than the phi -theta I thought earlier.
Try comparing your function to
phi = th + ArcTan[Cot[th] Sqrt[1 - (b^2)]] - Pi/2
Sorry, but I do not have time to check it myself at the moment.
Cheers.
kev said:Your welcome
Oops, I had it working yesrday. It looks like another error crept in when translating my scribbled calculations from paper to latex. I should have squared the gamma factor factor as well as the trigometric functions. It should have been:
In frame S' with relative motion ([itex]\beta_x,\beta_y[/itex]), the apparent length of the rod becomes
[tex]\| L' \| \ = L_x\sqrt{ sin^2\theta + \cos^2\theta(1-\beta_x^2-\beta_y^2)}[/tex]
where [tex]\theta = \tan^{-1}(\beta_y/\beta_x)[/tex]
which is the angle the velocity vector makes with the x' axis.
Rasalhague said:Ah, I see. Curiously, for this example, they're close, but not exactly the same:
ArcCos[(1 - Cos[th]^2)/(1 - b^2 Cos[th])] = 0.817476
ArcTan[Cot[th] Sqrt[1 - (b^2)]] - Pi/2 = ph - th = -0.811726
But in general, they're not the same. When I set by = 0, and all the velocity in the x direction, mine gives pi/2, yours indeterminate, but they should give 0 in that case, shouldn't they? I sometimes get imaginary results with mine, and it doesn't agree with the other equation I had either. Hmm. I'm going to have to leave it for now, and try again tomorrow...
kev said:In frame S' with relative motion ([itex]\beta_x,\beta_y[/itex]), the apparent length of the rod becomes
[tex]\| L' \| \ = L_x\sqrt{ sin^2\theta + \cos^2\theta(1-\beta_x^2-\beta_y^2)}[/tex]
where [tex]\theta = \tan^{-1}(\beta_y/\beta_x)[/tex]
which is the angle the velocity vector makes with the x' axis.
starthaus said:This is incorrect.
kev said:That is not very helpful.
Do you agree that [tex]\| L'\| \ = L_x \sqrt{1-\beta_x^2) [/tex] is correct.
I have checked that is what the the longer equation you said is incorrect simplifies to symbolically. If you agree the short equation is correct, then my original longer equation is correct too.
It is not a guess, it is derived from base principles using simple geometry.starthaus said:No, your "long" equation is just an incorrect guess.
kev said:It is not a guess, it is derived from base principles using simple geometry.
Do you agree my "Long" equation is symbolically and numerically the same as your short equation? If so then you have no reason to say the "long" equation is incorrect. I did not use your Starthaus blueprint method, or Rasalhague's method, but we all obtained the same result. Many ways to skin a cat!
Seriously, what are the chances of guessing:
[tex]\| L '\| \ = L_x\sqrt{ sin^2\theta + \cos^2\theta(1-\beta_x^2-\beta_y^2)}[/tex]
and the related:
[tex]\phi = \theta + tan^{-1}\left( \frac{\cos\theta}{\sin\theta } \, \sqrt{(1-\beta_x^2-\beta_y^2)} \right)-\frac{\pi}{2} [/tex]
and getting both right by pure luck?
SO you are suggesting that Rasalhague, using a completely different method and presumably making a different mistake, arrived at the same (but allegedly incorrect) solution as myself?starthaus said:Problem is, you guessed both wrong.
BTW, your formula for [tex]\theta[/tex] is also wrong.
kev said:SO you are suggesting that Rasalhague, using a completely different method and presumably making a different mistake, arrived at the same (but allegedly incorrect) solution as myself?
starthaus said:You try (and fail) to calculate length contraction. Two very different issues. I suggest that you take a good hard look at your derivation, you will find your mistakes. There are three mistakes.
starthaus said:There are three mistakes.
My derivation shows this formula to be true only if the rod is aligned with the x-axis and only if the two frames have aligned axes. It doesn't say anything about the case of arbitrary alignment. where your attempt at guessing the result fails.kev said:You are missing the very simple and self evident truth that if:
[tex]\| L'\| \ = L_x \sqrt{1-\beta_x^2} [/tex]
is true
and if:
[tex] L_x \sqrt{1-\beta_x^2} = L_x\sqrt{ sin^2\theta + \cos^2\theta(1-\beta_x^2-\beta_y^2)}[/tex]
is also true,
Rasalhague said:When you said "you both guessed wrong"
starthaus said:No, I was referring to kev's two guesses. Your formula is correct but it doesn't bring you any closer to finding the answer to the question. If you complete the calculations in my hint, you'll find the correct formula.
kev said:You are missing the very simple and self evident truth that if:
[tex]\| L'\| \ = L_x \sqrt{1-\beta_x^2} [/tex]
is true...starthaus said:My derivation shows this formula to be true only if the rod is aligned with the x-axis and only if the two frames have aligned axes. It doesn't say anything about the case of arbitrary alignment. where your attempt at guessing the result fails.
which makes it clear I am considering the limited case of the rod being aligned withe x-axis and arbitrary motion with respect to the x and y axes. This is exactly the same limited case that you considered in your https://www.physicsforums.com/blog.php?b=1959 .kev said:...
Let us say we have a rod with rest length Lx aligned with the x-axis in frame S with one end at the origin of the S frame. In frame S' with relative motion ([itex]\beta_x,\beta_y[/itex]), the apparent length of the rod becomes
[tex]\| L' \| \ = L_x\sqrt{ sin^2\theta + \cos^2\theta\ (1-\beta_x^2-\beta_y^2)}[/tex]
where [tex]\theta = \tan^{-1}(\beta_y/\beta_x)[/tex]
This is misleading, because it implies that for any arbitrary orientation of the rod, the length contraction is a function of [tex]\beta_x[/tex] only, which is obviously untrue if the rod is aligned with y or z axis. The correct way to write the expression is the way I did as:starthaus said:... the proof shows clearly:
[tex]L'=L\sqrt{1-\beta_x^2}[/tex]
No dependency whatsoever of [tex]\beta_y[/tex] or [tex]\beta_z[/tex]
I never claimed it was for arbitrary orientation of the rod and it made it clear it wasn't. It isn't, just as your solution in your blog isn't. The motion is for arbitrary motion relative to the x and y axes, and you are always free in the 3 dimensional case to rotate the axes so the motion lies in the x,y plane or for even greater simplification to align the motion with the x axis.starthaus said:But [tex] L_x\sqrt{ sin^2\theta + \cos^2\theta(1-\beta_x^2-\beta_y^2)}[/tex]
isn't the correct formula for arbitrary orientation of rod and arbitrary motion between the frames of reference. You are simply copying the formula for the simpler case, when the rod is aligned with the x axis.
Rasalhague said:Assuming "hint" is the PDF attachment to #84,
and assuming "your formula" is
[tex]\textbf{r}'=\textbf{r}-\frac{\textbf{r}\cdot \pmb{\beta}}{\beta^2}\left ( 1-\frac{1}{\gamma} \right )[/tex]
and assuming "the question" is how to find the angle between r and r' according to the boosted-to coordimate system,
you seem to be saying that my formula for this angle is wrong:
[tex]\frac{1-\cos^2(\textbf{r},\pmb{\beta})(1-1/\gamma)}{\sqrt{1-\beta^2\cos^2(\textbf{r},\pmb{\beta})}}[/tex]
Alternatively, if by "the question" you mean my more general question: "what is the nature of Thomas rotation, how can its effects be calculated", and that my formula does give the correct angle between these 3-vectors according to the boosted-to coordinate system, perhaps you're saying that this is not relevant to Thomas rotation. Then perhaps the "correct formula" you invite me to derive from your hint is a formula for some other quantity. If that's the case, could you elaborate a bit on why that is, and explain what quantity it is that I should be trying to find a formula for?
starthaus said:...
No, the question is finding the length contraction formula for the general case (see the hint file)
...
No (see above)> What I am saying is that , though the angle formula is correct, it doesn't help in finding the length contraction formula.
...
Sure, with pleasure,we were discussing length contraction when kev butted in.
kev said:The heading of the thread is "Stone's derivation of Thomas rotation", not "Stone's derivation of length contraction".
Lets look at the "testimony":starthaus said:Sure but we were discussing length contraction. Testimony is that your only contribution has been guesses about length contraction.
Rasalhague said:Brilliant! Thanks for that (kev). I was having no end of trouble getting through all those various powers of beta and gamma.
Rasalhague said:Thanks kev and starthaus for all your help.
Rasalhague said:Were they the computations Kev showed in #66? These fill in the gap in your PDF that I was having trouble with. I've worked through them, and they seem okay to me. They reach the same conclusion you did.
starthaus said:I have already answered that, it isn't a "pure" rotation. It is simply called a rotation , by abuse of language.
starthaus said:Sure but we were discussing length contraction. Testimony is that your only contribution has been guesses about length contraction.
starthaus said:No, the question is finding the length contraction formula for the general case (see the hint file).
Rasalhague said:That'd just be the norm of r', wouldn't it? When all the cancelling's done: r' = (r2-(r.b)2)1/2. I used it to get my formula for the angle.
kev said:Your above statement suggests that the rotation is somehow not physical
Rasalhague said:For the effect of a general boost on a general space vector, I get:
[tex]\textbf{r}'=\textbf{r}-\frac{\textbf{r}\cdot \pmb{\beta}}{\beta^2}\pmb{\beta}\left ( 1-\frac{1}{\gamma(\beta)} \right )[/tex]