Velocity-addition formula
In relativistic physics, a velocity-addition formula is an equation that specifies how to combine the velocities of objects in a way that is consistent with the requirement that no object's speed can exceed the speed of light. Such formulas apply to successive Lorentz transformations, so they also relate different frames. Accompanying velocity addition is a kinematic effect known as Thomas precession, whereby successive non-collinear Lorentz boosts become equivalent to the composition of a rotation of the coordinate system and a boost.
Standard applications of velocity-addition formulas include the Doppler shift, Doppler navigation, the aberration of light, and the dragging of light in moving water observed in the 1851 Fizeau experiment.<ref>Template:Harvnb</ref>
The notation employs Template:Math as velocity of a body within a Lorentz frame Template:Math, and Template:Math as velocity of a second frame Template:Math, as measured in Template:Math, and Template:Math as the transformed velocity of the body within the second frame.
HistoryEdit
The speed of light in a fluid is slower than the speed of light in vacuum, and it changes if the fluid is moving along with the light. In 1851, Fizeau measured the speed of light in a fluid moving parallel to the light using an interferometer. Fizeau's results were not in accord with the then-prevalent theories. Fizeau experimentally correctly determined the zeroth term of an expansion of the relativistically correct addition law in terms of Template:Math as is described below. Fizeau's result led physicists to accept the empirical validity of the rather unsatisfactory theory by Fresnel that a fluid moving with respect to the stationary aether partially drags light with it, i.e. the speed is Template:Math instead of Template:Math, where Template:Math is the speed of light in the aether, Template:Math is the refractive index of the fluid, and Template:Math is the speed of the fluid with respect to the aether.
The aberration of light, of which the easiest explanation is the relativistic velocity addition formula, together with Fizeau's result, triggered the development of theories like Lorentz aether theory of electromagnetism in 1892. In 1905 Albert Einstein, with the advent of special relativity, derived the standard configuration formula (Template:Math in the Template:Nowrap) for the addition of relativistic velocities.<ref name=Einstein_1905>Template:Harvnb</ref> The issues involving aether were, gradually over the years, settled in favor of special relativity.
Galilean relativityEdit
It was observed by Galileo that a person on a uniformly moving ship has the impression of being at rest and sees a heavy body falling vertically downward.<ref>Template:Harvnb</ref> This observation is now regarded as the first clear statement of the principle of mechanical relativity. Galileo saw that from the point of view of a person standing on the shore, the motion of falling downwards on the ship would be combined with, or added to, the forward motion of the ship.<ref>Template:Harvnb Galileo used this insight to show that the path of the weight when seen from the shore would be a parabola.</ref> In terms of velocities, it can be said that the velocity of the falling body relative to the shore equals the velocity of that body relative to ship plus the velocity of the ship relative to the shore.
In general for three objects A (e.g. Galileo on the shore), B (e.g. ship), C (e.g. falling body on ship) the velocity vector <math>\mathbf{u}</math> of C relative to A (velocity of falling object as Galileo sees it) is the sum of the velocity <math>\mathbf{u'}</math> of C relative to B (velocity of falling object relative to ship) plus the velocity Template:Math of B relative to A (ship's velocity away from the shore). The addition here is the vector addition of vector algebra and the resulting velocity is usually represented in the form
<math display="block"> \mathbf{u} = \mathbf{v} + \mathbf{u'}.</math>
The cosmos of Galileo consists of absolute space and time and the addition of velocities corresponds to composition of Galilean transformations. The relativity principle is called Galilean relativity. It is obeyed by Newtonian mechanics.
Special relativityEdit
According to the theory of special relativity, the frame of the ship has a different clock rate and distance measure, and the notion of simultaneity in the direction of motion is altered, so the addition law for velocities is changed. This change is not noticeable at low velocities but as the velocity increases towards the speed of light it becomes important. The addition law is also called a composition law for velocities. For collinear motions, the speed of the object, <math>u'</math>, e.g. a cannonball fired horizontally out to sea, as measured from the ship, moving at speed <math>v</math>, would be measured by someone standing on the shore and watching the whole scene through a telescope as<ref>Template:Cite book Extract of page 367</ref> <math display="block"> u = {v+u'\over 1+(vu'/c^2)} . </math> The composition formula can take an algebraically equivalent form, which can be easily derived by using only the principle of constancy of the speed of light,<ref>Template:Harvnb</ref> <math display="block"> {c-u \over c+u} = \left({c-u' \over c+u'}\right)\left({c-v \over c+v}\right).</math> The cosmos of special relativity consists of Minkowski spacetime and the addition of velocities corresponds to composition of Lorentz transformations. In the special theory of relativity Newtonian mechanics is modified into relativistic mechanics.
Standard configurationEdit
The formulas for boosts in the standard configuration follow most straightforwardly from taking differentials of the inverse Lorentz boost in standard configuration.<ref name=LL>Template:Harvnb</ref><ref>Template:Harvnb</ref> If the primed frame is travelling with speed <math>v</math> with Lorentz factor <math display="inline">\gamma_{_v} = 1/\sqrt{1-v^2/c^2}</math> in the positive Template:Nowrap relative to the unprimed frame, then the differentials are
<math display="block">dx = \gamma_{_v}(dx' + vdt'), \quad dy = dy', \quad dz = dz', \quad dt = \gamma_{_v}\left(dt' + \frac{v}{c^2}dx'\right).</math>
Divide the first three equations by the fourth,
<math display="block">\frac{dx}{dt} = \frac{\gamma_{_v}(dx' + vdt')}{\gamma_{_v}(dt' + \frac{v}{c^2}dx')},
\quad \frac{dy}{dt} = \frac{dy'}{\gamma_{_v}(dt' + \frac{v}{c^2}dx')}, \quad \frac{dz}{dt} = \frac{dz'}{\gamma_{_v}(dt' + \frac{v}{c^2}dx')},</math>
or
<math display="block">u_x = \frac{dx}{dt} = \frac{\frac{dx'}{dt'} + v}{(1 + \frac{v}{c^2}\frac{dx'}{dt'})}, \quad u_y = \frac{dy}{dt} = \frac{\frac{dy'}{dt'}}{\gamma_{_v} \ (1 + \frac{v}{c^2}\frac{dx'}{dt'})}, \quad u_z = \frac{dz}{dt} = \frac{\frac{dz'}{dt'}}{\gamma_{_v} \ (1 + \frac{v}{c^2}\frac{dx'}{dt'})},</math>
which is
Template:Equation box 1}{1 + \frac{v}{c^2}u_x'}, \quad u_y' = \frac{u_y\sqrt{1-\frac{v^2}{c^2}}}{1 - \frac{v}{c^2}u_x},</math> <math display="block">u_z = \frac{u_z'\sqrt{1-\frac{v^2}{c^2}}}{1 + \frac{v}{c^2}u_x'}, \quad u_z' = \frac{u_z\sqrt{1-\frac{v^2}{c^2}}}{1 - \frac{v}{c^2}u_x},</math> |cellpadding= 6 |border |border colour = #0073CF |bgcolor=#F9FFF7}}
in which expressions for the primed velocities were obtained using the standard recipe by replacing Template:Math by Template:Math and swapping primed and unprimed coordinates. If coordinates are chosen so that all velocities lie in a (common) Template:Math plane, then velocities may be expressed as <math display="block">u_x = u\cos \theta, u_y = u\sin \theta,\quad u_x' = u'\cos \theta', \quad u_y' = u'\sin \theta',</math> (see polar coordinates) and one finds<ref name=Einstein_1905 /><ref>Template:Harvnb</ref> Template:Equation box 1{1 + \frac{v}{c^2}u'\cos \theta'},</math> <math>\tan \theta = \frac{u_y}{u_x} = \frac{\sqrt{1-\frac{v^2}{c^2}}u_y'}{u_x' + v} = \frac{\sqrt{1-\frac{v^2}{c^2}}u'\sin \theta'}{u'\cos \theta' + v}.</math> |cellpadding= 6 |border |border colour = #0073CF |bgcolor=#F9FFF7}}
Template:Hidden begin <math display="block">\begin{align} u &= \sqrt{u_x^2 + u_y^2} = \frac{\sqrt{(u_x'+v)^2 + (1-\frac{v^2}{c^2})u_y'^2}}{1 + \frac{v}{c^2}u_x'} = \frac{\sqrt{u_x'^2+v^2+2u_x'v + (1-\frac{v^2}{c^2})u_y'^2}}{1 + \frac{v}{c^2}u_x'}\\ &=\frac{\sqrt{u'^2 \cos^2 \theta'+v^2+2vu'\cos \theta' + u'^2\sin^2\theta' - \frac{v^2}{c^2}u'^2\sin^2\theta'}}{1 + \frac{v}{c^2}u_x'}\\ &=\frac{\sqrt{u'^2 +v^2+2vu'\cos \theta' - (\frac{vu'\sin\theta'}{c})^2}}{1 + \frac{v}{c^2}u'\cos \theta'} \end{align}</math>
The proof as given is highly formal. There are other more involved proofs that may be more enlightening, such as the one below. Template:Math proof
General configurationEdit
Starting from the expression in coordinates for Template:Math parallel to the Template:Nowrap, expressions for the perpendicular and parallel components can be cast in vector form as follows, a trick which also works for Lorentz transformations of other 3d physical quantities originally in set up standard configuration. Introduce the velocity vector Template:Math in the unprimed frame and Template:Math in the primed frame, and split them into components parallel (∥) and perpendicular (⊥) to the relative velocity vector Template:Math (see hide box below) thus <math display="block">\mathbf{u} = \mathbf{u}_\parallel + \mathbf{u}_\perp,\quad \mathbf{u}' = \mathbf{u}'_\parallel + \mathbf{u}'_\perp ,</math> then with the usual Cartesian standard basis vectors Template:Math, set the velocity in the unprimed frame to be <math display="block">\mathbf{u}_\parallel = u_x \mathbf{e}_x,\quad \mathbf{u}_\perp = u_y \mathbf{e}_y + u_z \mathbf{e}_z ,\quad \mathbf{v} = v\mathbf{e}_x,</math> which gives, using the results for the standard configuration, <math display="block">\mathbf u_\parallel = \frac{\mathbf u_\parallel' + \mathbf v}{1 + \frac{\mathbf v \cdot \mathbf u_\parallel'}{c^2}}, \quad \mathbf u_\perp = \frac{\sqrt{1-\frac{v^2}{c^2}}\mathbf u_\perp'}{1 + \frac{\mathbf v\cdot \mathbf u_\parallel'}{c^2}}. </math> where Template:Math is the dot product. Since these are vector equations, they still have the same form for Template:Math in any direction. The only difference from the coordinate expressions is that the above expressions refers to vectors, not components.
One obtains <math display="block">\mathbf{u} = \mathbf u_\parallel + \mathbf u_\perp = \frac{1}{1+\frac{\mathbf{v}\cdot\mathbf{u}'}{c^{2}}}\left[\alpha_v\mathbf{u}'+ \mathbf{v} + (1-\alpha_v)\frac{(\mathbf{v}\cdot\mathbf{u}')}{v^{2}}\mathbf{v}\right] \equiv \mathbf v \oplus \mathbf u',</math> where Template:Math is the reciprocal of the Lorentz factor. The ordering of operands in the definition is chosen to coincide with that of the standard configuration from which the formula is derived. Template:Hidden begin <math display="block">\begin{align} \frac{\mathbf u'_\parallel + \mathbf v}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}} + \frac{\alpha_v \mathbf u'_\perp}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}} &= \frac{\mathbf v + \frac{\mathbf v \cdot \mathbf u'}{v^2}\mathbf v}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}} + \frac{\alpha_v \mathbf u' - \alpha_v\frac{\mathbf v \cdot \mathbf u'}{v^2}\mathbf v}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\\
&=\frac{1 + \frac{\mathbf v \cdot \mathbf u'}{v^2}(1 - \alpha_v)}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\mathbf v + \alpha_v\frac{1}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\mathbf u'\\
&=\frac{1}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\mathbf v + \alpha_v\frac{1}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\mathbf u' + \frac{1}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\frac{\mathbf v \cdot \mathbf u'}{v^2}(1 - \alpha_v)\mathbf v\\
&=\frac{1}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\mathbf v + \alpha_v\frac{1}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\mathbf u' + \frac{1}{c^2}\frac{1}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\frac{\mathbf v \cdot \mathbf u'}{v^2/c^2}(1 - \alpha_v)\mathbf v\\
&=\frac{1}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\mathbf v + \alpha_v\frac{1}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\mathbf u' + \frac{1}{c^2}\frac{1}{1 + \frac{\mathbf v \cdot \mathbf u'}{c^2}}\frac{\mathbf v \cdot \mathbf u'}{(1-\alpha_v)(1+\alpha_v)}(1 - \alpha_v)\mathbf v\\
&= \frac{1}{1+\frac{\mathbf{v}\cdot\mathbf{u}'}{c^{2}}}\left[\alpha_v\mathbf{u}'+ \mathbf{v} + (1-\alpha_v)\frac{(\mathbf{v}\cdot\mathbf{u}')}{v^{2}}\mathbf{v}\right].
\end{align} </math>
Template:Hidden end Template:Hidden begin Either the parallel or the perpendicular component for each vector needs to be found, since the other component will be eliminated by substitution of the full vectors.
The parallel component of Template:Math can be found by projecting the full vector into the direction of the relative motion <math display="block">\mathbf{u}'_\parallel = \frac{\mathbf{v} \cdot \mathbf{u}'}{v^2}\mathbf v,</math> and the perpendicular component of Template:Math can be found by the geometric properties of the cross product (see figure above right), <math display="block">\mathbf{u}'_\perp = - \frac{\mathbf{v} \times (\mathbf{v} \times \mathbf{u}')}{v^2}.</math>
In each case, Template:Math is a unit vector in the direction of relative motion.
The expressions for Template:Math and Template:Math can be found in the same way. Substituting the parallel component into <math display="block">\mathbf u = \frac{\mathbf u_\parallel' + \mathbf v}{1 + \frac{\mathbf v \cdot \mathbf u_\parallel'}{c^2}} +\frac{\sqrt{1-\frac{v^2}{c^2}}(\mathbf u' - \mathbf u_\parallel')}{1 + \frac{\mathbf v\cdot \mathbf u_\parallel'}{c^2}}, </math>
results in the above equation.<ref>Template:Harvnb</ref>
Using an identity in <math>\alpha_v</math> and <math>\gamma_v</math>,<ref>Template:Harvnb</ref><ref group=nb>These formulae follow from inverting Template:Math for Template:Math and applying the difference of two squares to obtain Template:Block indent so that Template:Block indent</ref>
Template:Equation box 1\left[\mathbf v + \frac{\mathbf u'}{\gamma_v} + \frac{1}{c^2}\frac{\gamma_v}{1+\gamma_v}(\mathbf u' \cdot \mathbf v)\mathbf v\right]\\ &= \frac{1}{1 + \frac{\mathbf u' \cdot \mathbf v}{c^2}}\left[\mathbf v + \mathbf u' + \frac{1}{c^2}\frac{\gamma_v}{1+\gamma_v} \mathbf v \times(\mathbf v \times \mathbf u')\right], \end{align}</math> and in the forwards (v positive, S → S') direction <math display="block">\begin{align} \mathbf v \oplus \mathbf u \equiv \mathbf u' &=\frac{1}{1 - \frac{\mathbf u \cdot \mathbf v}{c^2}}\left[\frac{\mathbf u}{\gamma_v} - \mathbf v + \frac{1}{c^2}\frac{\gamma_v}{1+\gamma_v}(\mathbf u \cdot \mathbf v)\mathbf v\right]\\ &= \frac{1}{1 - \frac{\mathbf u \cdot \mathbf v}{c^2}}\left[ \mathbf u - \mathbf v + \frac{1}{c^2}\frac{\gamma_v}{1+\gamma_v} \mathbf v \times(\mathbf v \times \mathbf u)\right] \end{align}</math> |cellpadding= 6 |border |border colour = #0073CF |bgcolor=#F9FFF7}}
where the last expression is by the standard vector analysis formula Template:Math. The first expression extends to any number of spatial dimensions, but the cross product is defined in three dimensions only. The objects Template:Math with Template:Math having velocity Template:Math relative to Template:Math and Template:Math having velocity Template:Math relative to Template:Math can be anything. In particular, they can be three frames, or they could be the laboratory, a decaying particle and one of the decay products of the decaying particle.
PropertiesEdit
The relativistic addition of 3-velocities is non-linear, so in general <math display="block">(\lambda \mathbf{v}) \oplus (\lambda \mathbf{u}) \neq \lambda (\mathbf{v} \oplus \mathbf{u}) , </math> for real number Template:Math, although it is true that <math display="block">(-\mathbf{v}) \oplus (-\mathbf{u}) = - (\mathbf{v} \oplus \mathbf{u}) , </math>
Also, due to the last terms, is in general neither commutative <math display="block">\mathbf v \oplus \mathbf u \ne \mathbf u \oplus \mathbf v, </math> nor associative <math display="block">\mathbf v \oplus (\mathbf u \oplus \mathbf w) \ne (\mathbf v \oplus \mathbf u) \oplus \mathbf w. </math>
It deserves special mention that if Template:Math and Template:Math refer to velocities of pairwise parallel frames (primed parallel to unprimed and doubly primed parallel to primed), then, according to Einstein's velocity reciprocity principle, the unprimed frame moves with velocity Template:Math relative to the primed frame, and the primed frame moves with velocity Template:Math relative to the doubly primed frame hence Template:Math is the velocity of the unprimed frame relative to the doubly primed frame, and one might expect to have Template:Math by naive application of the reciprocity principle. This does not hold, though the magnitudes are equal. The unprimed and doubly primed frames are not parallel, but related through a rotation. This is related to the phenomenon of Thomas precession, and is not dealt with further here.
The norms are given by<ref>Template:Harvnb Equation (12.6) This is derived quite differently by consideration of invariant cross sections.</ref> <math display="block">| \mathbf u |^2 \equiv |\mathbf v \oplus \mathbf u'|^2 = \frac{1}{\left(1+\frac{\mathbf v \cdot \mathbf u'}{c^2}\right)^2}\left[\left(\mathbf v + \mathbf u' \right)^2 - \frac{1}{c^2}\left(\mathbf v \times \mathbf u'\right)^2 \right] = |\mathbf u' \oplus \mathbf v|^2.</math> and <math display="block">| \mathbf u' |^2 \equiv |\mathbf v \oplus \mathbf u|^2 = \frac{1}{\left(1-\frac{\mathbf v \cdot \mathbf u}{c^2}\right)^2}\left[\left(\mathbf u - \mathbf v \right)^2 - \frac{1}{c^2}\left(\mathbf v \times \mathbf u\right)^2 \right] = |\mathbf u \oplus \mathbf v|^2.</math>
Template:Math proof It is clear that the non-commutativity manifests itself as an additional rotation of the coordinate frame when two boosts are involved, since the norm squared is the same for both orders of boosts.
The gamma factors for the combined velocities are computed as <math display="block">\gamma_u = \gamma_{\mathbf v \oplus \mathbf u'} =\left[ 1 - \frac{1}{c^2}\frac{1}{(1+\frac{\mathbf v \cdot \mathbf u'}{c^2})^2} \left( (\mathbf v + \mathbf u')^2 - \frac{1}{c^2}(v^2u'^2 - (\mathbf v \cdot \mathbf u')^2)\right)\right]^{-\frac{1}{2}}=\gamma_v\gamma_u'\left(1+\frac{\mathbf v \cdot \mathbf u'}{c^2}\right), \quad \quad \gamma_u' = \gamma_v\gamma_u\left(1-\frac{\mathbf v \cdot \mathbf u}{c^2}\right)</math>
&=\left[ \frac{c^2(1+\frac{\mathbf v \cdot \mathbf u'}{c^2})^2 - (\mathbf v + \mathbf u')^2 + \frac{1}{c^2}(v^2u'^2 - (\mathbf v \cdot \mathbf u')^2)}{c^2(1+\frac{\mathbf v \cdot \mathbf u'}{c^2})^2} \right]^{-\frac{1}{2}}\\
&=\left[ \frac{c^2(1+2\frac{\mathbf v \cdot \mathbf u'}{c^2} + \frac{(\mathbf v \cdot \mathbf u')^2}{c^4}) - v^2 - u'^2 - 2(\mathbf v \cdot \mathbf u') + \frac{1}{c^2}(v^2u'^2 - (\mathbf v \cdot \mathbf u')^2)}{c^2(1+\frac{\mathbf v \cdot \mathbf u'}{c^2})^2} \right]^{-\frac{1}{2}}\\
&=\left[ \frac{1+2\frac{\mathbf v \cdot \mathbf u'}{c^2} + \frac{(\mathbf v \cdot \mathbf u')^2}{c^4} - \frac{v^2}{c^2} - \frac{u'^2}{c^2} - \frac{2}{c^2}(\mathbf v \cdot \mathbf u') + \frac{1}{c^4}(v^2u'^2 - (\mathbf v \cdot \mathbf u')^2)}{(1+\frac{\mathbf v \cdot \mathbf u'}{c^2})^2} \right]^{-\frac{1}{2}}\\
&=\left[ \frac{1 + \frac{(\mathbf v \cdot \mathbf u')^2}{c^4} - \frac{v^2}{c^2} - \frac{u'^2}{c^2} + \frac{1}{c^4}(v^2u'^2 - (\mathbf v \cdot \mathbf u')^2)}{(1+\frac{\mathbf v \cdot \mathbf u'}{c^2})^2} \right]^{-\frac{1}{2}}\\
&=\left[ \frac{\left(1-\frac{v^2}{c^2}\right)\left(1-\frac{u'^2}{c^2}\right)}{\left(1+\frac{\mathbf v \cdot \mathbf u'}{c^2}\right)^2} \right]^{-\frac{1}{2}} =\left[ \frac{1}{\gamma_v^2\gamma_u'^2\left(1+\frac{\mathbf v \cdot \mathbf u'}{c^2}\right)^2} \right]^{-\frac{1}{2}}\\
&=\gamma_v\gamma_u' \left(1+\frac{\mathbf v \cdot \mathbf u'}{c^2}\right) \end{align}</math>
Reverse formula found by using standard procedure of swapping Template:Math for Template:Math and Template:Math for Template:Math. }}
Notational conventionsEdit
Notations and conventions for the velocity addition vary from author to author. Different symbols may be used for the operation, or for the velocities involved, and the operands may be switched for the same expression, or the symbols may be switched for the same velocity. A completely separate symbol may also be used for the transformed velocity, rather than the prime used here. Since the velocity addition is non-commutative, one cannot switch the operands or symbols without changing the result.
Examples of alternative notation include:
- No specific operand
- Template:Harvtxt (using units where c = 1) <math display="block">| \mathbf{v_{rel}} |^2 =\frac{1}{(1 - \mathbf{v_1}\cdot\mathbf{v_2})^2}\left[(\mathbf{v_1}-\mathbf{v_2})^2 - (\mathbf{v_1} \times \mathbf{v_2})^2\right] </math>
- Left-to-right ordering of operands
- Template:Harvtxt <math display="block">\mathbf{u}\oplus\mathbf{v} = \frac{1}{1+\frac{\mathbf{u}\cdot\mathbf{v}}{c^2}}\left[\mathbf{v}+\mathbf{u}+\frac{1}{c^2}\frac{\gamma_\mathbf{u}}{\gamma_\mathbf{u}+1}\mathbf{u}\times(\mathbf{u}\times\mathbf{v})\right] </math> Template:BrTemplate:Harvtxt <math display="block">\mathbf{u}*\mathbf{v}=\frac{1}{1+\frac{\mathbf{u}\cdot\mathbf{v}}{c^2}}\left[\mathbf{v}+\mathbf{u}+\frac{1}{c^2}\frac{\gamma_\mathbf{u}}{\gamma_\mathbf{u}+1}\mathbf{u}\times(\mathbf{u}\times\mathbf{v})\right] </math>
- Right-to-left ordering of operands
- Template:Harvtxt <math display="block">\mathbf{w}\circ\mathbf{v}=\frac{1}{1+\frac{\mathbf{v}\cdot\mathbf{w}}{c^{2}}}\left[\frac{\mathbf{w}}{\gamma_\mathbf{v}}+\mathbf{v}+\frac{1}{c^{2}}\frac{\gamma_\mathbf{v}}{\gamma_\mathbf{v}+1}(\mathbf{w}\cdot\mathbf{v})\mathbf{v}\right]</math>
ApplicationsEdit
Some classical applications of velocity-addition formulas, to the Doppler shift, to the aberration of light, and to the dragging of light in moving water, yielding relativistically valid expressions for these phenomena are detailed below. It is also possible to use the velocity addition formula, assuming conservation of momentum (by appeal to ordinary rotational invariance), the correct form of the Template:Math-vector part of the momentum four-vector, without resort to electromagnetism, or a priori not known to be valid, relativistic versions of the Lagrangian formalism. This involves experimentalist bouncing off relativistic billiard balls from each other. This is not detailed here, but see for reference Template:Harvtxt Wikisource version (primary source) and Template:Harvtxt.
Fizeau experimentEdit
{{#invoke:Labelled list hatnote|labelledList|Main article|Main articles|Main page|Main pages}} When light propagates in a medium, its speed is reduced, in the rest frame of the medium, to Template:Math, where Template:Math is the index of refraction of the medium Template:Math. The speed of light in a medium uniformly moving with speed Template:Math in the positive Template:Math-direction as measured in the lab frame is given directly by the velocity addition formulas. For the forward direction (standard configuration, drop index Template:Mvar on Template:Math) one gets,<ref>Template:Harvnb</ref> <math display="block">\begin{align} c_m &= \frac{V + c_m'}{1 + \frac{Vc_m'}{c^2}} = \frac{V + \frac{c}{n}}{1 + \frac{Vc}{nc^2}} = \frac{c}{n} \frac{1 + \frac{nV}{c}}{1 + \frac{V}{nc}}\\ & = \frac{c}{n} \left(1 + \frac{nV}{c}\right) \frac{1}{1 + \frac{V}{nc}} = \left(\frac{c}{n} + V\right) \left(1 - \frac{V}{nc} + \left(\frac{V}{nc}\right)^2 - \cdots\right). \end{align}</math>
Collecting the largest contributions explicitly, <math display="block">c_m = \frac{c}{n} + V\left(1 - \frac{1}{n^2} - \frac{V}{nc} + \cdots\right).</math> Fizeau found the first three terms.<ref>Template:Harvnb</ref><ref>Template:Harvnb</ref> The classical result is the first two terms.
Aberration of lightEdit
{{#invoke:Labelled list hatnote|labelledList|Main article|Main articles|Main page|Main pages}}
Another basic application is to consider the deviation of light, i.e. change of its direction, when transforming to a new reference frame with parallel axes, called aberration of light. In this case, Template:Math, and insertion in the formula for Template:Math yields <math display="block">\tan \theta = \frac{\sqrt{1-\frac{V^2}{c^2}}c\sin \theta'}{c\cos \theta' + V} = \frac{\sqrt{1-\frac{V^2}{c^2}}\sin \theta'}{\cos \theta' + \frac{V}{c}}.</math>
For this case one may also compute Template:Math and Template:Math from the standard formulae,<ref>Template:Harvnb</ref> <math display="block">\sin \theta =\frac{\sqrt{1-\frac{V^2}{c^2}}\sin \theta'}{1+\frac{V}{c}\cos \theta'},</math>
Template:Hidden begin <math display="block">\begin{align}\frac{v_y}{v} &= \frac{\frac{\sqrt{1-\frac{V^2}{c^2}}v_y'}{1 + \frac{V}{c^2}v_x'}}{\frac{\sqrt{v'^2 +V^2+2Vv'\cos \theta' - (\frac{Vv'\sin\theta'}{c})^2}}{1 + \frac{V}{c^2}v'\cos \theta'}}\\
&= \frac{c\sqrt{1-\frac{V^2}{c^2}}\sin \theta'}{\sqrt{c^2 +V^2+2Vc\cos \theta' - V^2\sin^2\theta'}}\\
&= \frac{c\sqrt{1-\frac{V^2}{c^2}}\sin \theta'}{\sqrt{c^2 +V^2+2Vc\cos \theta' - V^2(1 - \cos^2\theta')}} = \frac{c\sqrt{1-\frac{V^2}{c^2}}\sin \theta'}{\sqrt{c^2 +2Vc\cos \theta' + V^2\cos^2\theta'}}\\
&= \frac{\sqrt{1-\frac{V^2}{c^2}}\sin \theta'}{1+\frac{V}{c}\cos \theta'}, \end{align}</math>
Template:Hidden end <math display="block">\cos \theta = \frac{\frac{V}{c} + \cos \theta'}{1+\frac{V}{c}\cos \theta'},</math>
the trigonometric manipulations essentially being identical in the Template:Math case to the manipulations in the Template:Math case. Consider the difference,
<math display="block">\begin{align}\sin \theta - \sin \theta' &= \sin \theta'\left(\frac{\sqrt{1 - \frac{V^2}{c^2}}}{1 + \frac{V}{c} \cos \theta'} - 1\right)\\ &\approx \sin \theta'\left(1 -\frac{V}{c} \cos \theta' - 1\right) = -\frac{V}{c}\sin\theta'\cos\theta',\end{align}</math> correct to order Template:Math. Employ in order to make small angle approximations a trigonometric formula, <math display="block">\sin \theta' - \sin \theta = 2\sin \frac{1}{2}(\theta'-\theta)\cos\frac{1}{2}(\theta + \theta') \approx (\theta' - \theta)\cos\theta', </math> where Template:Math were used.
Thus the quantity <math display="block">\Delta \theta \equiv \theta' - \theta = \frac{V}{c}\sin \theta',</math> the classical aberration angle, is obtained in the limit Template:Math.
Relativistic Doppler shiftEdit
{{#invoke:Labelled list hatnote|labelledList|Main article|Main articles|Main page|Main pages}}
Here velocity components will be used as opposed to speed for greater generality, and in order to avoid perhaps seemingly ad hoc introductions of minus signs. Minus signs occurring here will instead serve to illuminate features when speeds less than that of light are considered.
For light waves in vacuum, time dilation together with a simple geometrical observation alone suffices to calculate the Doppler shift in standard configuration (collinear relative velocity of emitter and observer as well of observed light wave).
All velocities in what follows are parallel to the common positive Template:Nowrap, so subscripts on velocity components are dropped. In the observers frame, introduce the geometrical observation <math display="block">\lambda = -sT + VT = (-s + V)T</math> as the spatial distance, or wavelength, between two pulses (wave crests), where Template:Math is the time elapsed between the emission of two pulses. The time elapsed between the passage of two pulses at the same point in space is the time period Template:Mvar, and its inverse Template:Math is the observed (temporal) frequency. The corresponding quantities in the emitters frame are endowed with primes.<ref>Template:Harvnb In the reference, the speed of an approaching emitter is taken as positive. Hence the sign difference.</ref>
For light waves <math display="block">s = s' = -c,</math> and the observed frequency is<ref name=Einstein_1905 /><ref>Template:Harvnb</ref><ref>Template:Harvnb</ref> <math display="block">\nu = {-s \over \lambda} = {-s \over (V-s)T} = {c \over (V+c)\gamma_{_V} T'} = \nu'\frac{c\sqrt{1 - {V^2 \over c^2}}}{c+V} = \nu'\sqrt{\frac{1-\beta}{1+\beta}}\,.</math> where Template:Math is standard time dilation formula.
Suppose instead that the wave is not composed of light waves with speed Template:Mvar, but instead, for easy visualization, bullets fired from a relativistic machine gun, with velocity Template:Math in the frame of the emitter. Then, in general, the geometrical observation is precisely the same. But now, Template:Math, and Template:Math is given by velocity addition, <math display="block">s = \frac{s' + V}{1+{s'V\over c^2}}.</math>
The calculation is then essentially the same, except that here it is easier carried out upside down with Template:Math instead of Template:Mvar. One finds
Template:Equation box 1\right), \quad \nu = \gamma_{_V}\nu'\left(1+{V\over s'}\right)</math> |cellpadding=6|border|border colour = #0073CF|bgcolor=#F9FFF7}}
Template:Hidden begin <math display="block">\begin{align}{L\over -s} &= \frac{\left(\frac{-s'-V}{1+{s'V\over c^2}} + V\right) T} {\frac{-s'-V}{1+{s'V\over c^2}}}\\
&={\gamma_{_V} \over \nu'}\frac{-s'-V + V(1+{s'V\over c^2})}{-s'-V}\\ &={\gamma_{_V} \over \nu'}\left(\frac{s'\left(1-{V^2\over c^2}\right)}{s'+V}\right)\\ &={\gamma_{_V} \over \nu'}\left(\frac{s'\gamma^{-2}}{s'+V}\right)\\ &={1 \over \gamma_{_V}\nu'}\left(\frac{1}{1+{V\over s'}}\right).\\ \end{align}</math>
Observe that in the typical case, the Template:Math that enters is negative. The formula has general validity though.<ref group=nb>Note that Template:Math is negative in the sense for which that the problem is set up, i.e. emitter with positive velocity fires fast bullets towards observer in unprimed system. The convention is that Template:Math should yield positive frequency in accordance with the result for the ultimate velocity, Template:Math. Hence the minus sign is a convention, but a very natural convention, to the point of being canonical.Template:PbThe formula may also result in negative frequencies. The interpretation then is that the bullets are approaching from the negative Template:Math-axis. This may have two causes. The emitter can have large positive velocity and be firing slow bullets. It can also be the case that the emitter has small negative velocity and is firing fast bullets. But if the emitter has a large negative velocity and is firing slow bullets, the frequency is again positive.Template:PbFor some of these combination to make sense, it must be required that the emitter has been firing bullets for sufficiently long time, in the limit that the Template:Math-axis at any instant has equally spaced bullets everywhere.</ref> When Template:Math, the formula reduces to the formula calculated directly for light waves above,
<math display="block">\nu = \nu'\gamma_{_V}(1-\beta) = \nu'\frac{1-\beta}{\sqrt{1-\beta}\sqrt{1+\beta}}=\nu'\sqrt{\frac{1-\beta}{1+\beta}}\,.</math>
If the emitter is not firing bullets in empty space, but emitting waves in a medium, then the formula still applies, but now, it may be necessary to first calculate Template:Math from the velocity of the emitter relative to the medium.
Returning to the case of a light emitter, in the case the observer and emitter are not collinear, the result has little modification,<ref name=Einstein_1905 /><ref>Template:Harvnb</ref><ref>Template:Harvnb</ref> <math display="block">\nu = \gamma_{_V}\nu' \left(1+\frac{V}{s'}\cos\theta\right),</math> where Template:Mvar is the angle between the light emitter and the observer. This reduces to the previous result for collinear motion when Template:Math, but for transverse motion corresponding to Template:Math, the frequency is shifted by the Lorentz factor. This does not happen in the classical optical Doppler effect.
Hyperbolic geometryEdit
Associated to the relativistic velocity <math>\boldsymbol \beta</math> of an object is a quantity <math>\boldsymbol{\zeta}</math> whose norm is called rapidity. These are related through <math display="block">\mathfrak{so}(3,1) \supset \mathrm{span}\{K_1, K_2, K_3\} \approx \mathbb{R}^3 \ni \boldsymbol{\zeta} = \boldsymbol{\hat{\beta}} \tanh^{-1}\beta, \quad \boldsymbol{\beta} \in \mathbb{B}^3,</math> where the vector <math>\boldsymbol \zeta</math> is thought of as being Cartesian coordinates on a 3-dimensional subspace of the Lie algebra <math>\mathfrak{so}(3, 1)</math> of the Lorentz group spanned by the boost generators <math>K_1, K_2, K_3</math>. This space, call it rapidity space, is isomorphic to Template:Math as a vector space, and is mapped to the open unit ball, <math> \mathbb B^3</math>, velocity space, via the above relation.<ref>Template:Harvnb</ref> The addition law on collinear form coincides with the law of addition of hyperbolic tangents <math display="block">\tanh(\zeta_v + \zeta_{u'}) = {\tanh \zeta_v + \tanh \zeta_{u'} \over 1+ \tanh \zeta_v \tanh \zeta_{u'}}</math> with <math display="block">\frac{v}{c} = \tanh \zeta_v \ , \quad \frac{u'}{c} = \tanh \zeta_{u'} \ , \quad\, \frac{u}{c} = \tanh(\zeta_v + \zeta_{u'}).</math>
The line element in velocity space <math>\mathbb B^3</math> follows from the expression for relativistic relative velocity in any frame,<ref>Template:Harvnb</ref> <math display="block">v_{r} = \frac{\sqrt{(\mathbf{v_1}-\mathbf{v_2})^2 - (\mathbf{v_1} \times \mathbf{v_2})^2}}{1 - \mathbf{v_1}\cdot\mathbf{v_2}},</math> where the speed of light is set to unity so that <math>v_i</math> and <math>\beta_i</math> agree. It this expression, <math>\mathbf{v}_1</math> and <math>\mathbf{v}_2</math> are velocities of two objects in any one given frame. The quantity <math>v_r</math> is the speed of one or the other object relative to the other object as seen in the given frame. The expression is Lorentz invariant, i.e. independent of which frame is the given frame, but the quantity it calculates is not. For instance, if the given frame is the rest frame of object one, then <math>v_r = v_2</math>.
The line element is found by putting <math>\mathbf v_2 = \mathbf v_1 + d\mathbf v</math> or equivalently <math>\boldsymbol \beta_2 = \boldsymbol \beta_1 + d\boldsymbol \beta</math>,<ref>Template:Harvnb</ref> <math display="block">dl_\boldsymbol{\beta}^2 = \frac{d\boldsymbol \beta^2 - (\boldsymbol \beta \times d\boldsymbol \beta)^2}{(1-\beta^2)^2} = \frac{d\beta^2}{(1-\beta^2)^2} + \frac{\beta^2}{1-\beta^2}(d\theta^2 + \sin^2\theta d\varphi^2),</math> with Template:Math and Template:Mvar the usual spherical angle coordinates for <math>\boldsymbol \beta</math> taken in the Template:Math-direction. Now introduce Template:Mvar through <math display="block">\zeta = |\boldsymbol \zeta| = \tanh^{-1}\beta,</math> and the line element on rapidity space <math>\mathbb R^3</math> becomes <math display="block">dl_{\boldsymbol \zeta}^2 = d\zeta^2 + \sinh^2\zeta(d\theta^2 + \sin^2\theta d\varphi^2).</math>
Relativistic particle collisionsEdit
In scattering experiments the primary objective is to measure the invariant scattering cross section. This enters the formula for scattering of two particle types into a final state <math>f</math> assumed to have two or more particles,<ref>Template:Harvnb</ref> <math display="block">dN_f = R_f \, dV \, dt = \sigma F \, dV \, dt</math> or, in most textbooks, <math display="block">dN_f = \sigma n_1 n_2 v_r \, dV \, dt</math> where
- <math>dVdt</math> is spacetime volume. It is an invariant under Lorentz transformations.
- <math>dN_f</math> is the total number of reactions resulting in final state <math>f</math> in spacetime volume <math>dVdt</math>. Being a number, it is invariant when the same spacetime volume is considered.
- <math>R_f = F\sigma</math> is the number of reactions resulting in final state <math>f</math> per unit spacetime, or reaction rate. This is invariant.
- <math>F = n_1n_2v_{r}</math> is called the incident flux. This is required to be invariant, but isn't in the most general setting.
- <math>\sigma</math> is the scattering cross section. It is required to be invariant.
- <math>n_1, n_2 </math> are the particle densities in the incident beams. These are not invariant as is clear due to length contraction.
- <math>v_{r} = |\mathbf v_2 - \mathbf v_1|</math> is the relative speed of the two incident beams. This cannot be invariant since <math>F = n_1n_2v_{r}</math> is required to be so.
The objective is to find a correct expression for relativistic relative speed <math>v_\text{rel}</math> and an invariant expression for the incident flux.
Non-relativistically, one has for relative speed <math>v_r = |\mathbf v_2 - \mathbf v_1|</math>. If the system in which velocities are measured is the rest frame of particle type <math>1</math>, it is required that <math>v_\text{rel} = v_r = |\mathbf v_2|.</math> Setting the speed of light <math>c = 1</math>, the expression for <math>v_\text{rel}</math> follows immediately from the formula for the norm (second formula) in the general configuration as<ref name="Cannoni 2017 4">Template:Harvnb</ref><ref>Template:Harvnb</ref> <math display="block">v_\text{rel} =\frac{\sqrt{(\mathbf{v_1}-\mathbf{v_2})^2 - (\mathbf{v_1} \times \mathbf{v_2})^2}}{1 - \mathbf{v_1}\cdot\mathbf{v_2}}.</math>
The formula reduces in the classical limit to <math>v_r = |\mathbf v_1 - \mathbf v_2|</math> as it should, and gives the correct result in the rest frames of the particles. The relative velocity is incorrectly given in most, perhaps all books on particle physics and quantum field theory.<ref name="Cannoni 2017 4" /> This is mostly harmless, since if either one particle type is stationary or the relative motion is collinear, then the right result is obtained from the incorrect formulas. The formula is invariant, but not manifestly so. It can be rewritten in terms of four-velocities as <math display="block">v_\text{rel} = \frac{\sqrt{(u_1 \cdot u_2)^2 - 1}}{u_1 \cdot u_2}.</math>
The correct expression for the flux, published by Christian Møller<ref>Template:Harvnb</ref> in 1945, is given by<ref>Template:Harvnb</ref> <math display="block">F = n_1n_2\sqrt{(\mathbf v_1 - \mathbf v_2)^2 - (\mathbf v_1 \times \mathbf v_2)^2} \equiv n_1n_2\bar v.</math>
One notes that for collinear velocities, <math>F = n_1n_2|\mathbf v_2 - \mathbf v_1| = n_1n_2v_r</math>. In order to get a manifestly Lorentz invariant expression one writes <math>J_i = (n_i, n_i\mathbf v_i)</math> with <math>n_i = \gamma_i n_i^0</math>, where <math>n_i^0</math> is the density in the rest frame, for the individual particle fluxes and arrives at<ref name="Cannoni 2017 13">Template:Harvnb</ref> <math display="block">F = (J_1 \cdot J_2) v_\text{rel}.</math>
In the literature the quantity <math>\bar v</math> as well as <math>v_r</math> are both referred to as the relative velocity. In some cases (statistical physics and dark matter literature), <math>\bar v</math> is referred to as the Møller velocity, in which case <math>v_r</math> means relative velocity. The true relative velocity is at any rate <math>v_\text{rel}</math>.<ref name="Cannoni 2017 13"/> The discrepancy between <math>v_\text{rel}</math> and <math>v_r</math> is relevant though in most cases velocities are collinear. At LHC the crossing angle is small, around Template:Val, but at the old Intersecting Storage Ring at CERN, it was about 18°.<ref>Template:Harvnb</ref>
With rapidityEdit
{{#invoke:Labelled list hatnote|labelledList|Main article|Main articles|Main page|Main pages}}
For collinear velocities interpreted as rapidity, the sum formula is simple addition. According to Minkowski, the time-like vector (zeitartiger Vektor) for a given duration lies on a hyperbola. Since the hyperbola is traced by a hyperbolic angle, and velocity determines a point on the hyperbola, there is a hyperbolic angle called rapidity associated with a particular velocity. For a unit of duration, the unit hyperbola provides a reference, where hyperbolic angle forms a one-parameter group isomorphic to the real number line under addition. The form of the formula for velocity addition can then be accounted in terms of the hyperbolic tangent function tanh which takes hyperbolic angle (rapidity) as an argument. In fact, the hyperbolic tangent of rapidity is the ratio of velocity to the speed of light in vacuum.<ref>Template:Cite book</ref>
See alsoEdit
RemarksEdit
NotesEdit
ReferencesEdit
- Template:Cite journal
- Template:Cite journal
- Template:Cite book
- Template:Cite book
- Template:Cite book
- Template:Cite book (graduate level)
- Template:Cite book (introductory level)
- Template:Cite book (graduate level)
- Template:Cite book
- Template:Cite book
- Template:Cite journal
- Template:Cite journal
- Template:Cite book
- Template:Cite book
- Template:Cite book
- Template:Cite book
- Template:Cite journal
- Historical
- Template:Cite journal
- Template:Citation
- Template:Cite journal
- Template:Cite journal
- Template:Cite journal
- Template:Cite journal
- Template:Cite book
- Template:Cite book
- Template:Cite journal Wikisource version