Physics:Derivations of the Lorentz transformations

There are many ways to derive the Lorentz transformations using a variety of physical principles, ranging from Maxwell's equations to Einstein's postulates of special relativity, and mathematical tools, spanning from elementary algebra and hyperbolic functions, to linear algebra and group theory.

This article provides a few of the easier ones to follow in the context of special relativity, for the simplest case of a Lorentz boost in standard configuration, i.e. two inertial frames moving relative to each other at constant (uniform) relative velocity less than the speed of light, and using Cartesian coordinates so that the x and x′ axes are collinear.

Lorentz transformation

In the fundamental branches of modern physics, namely general relativity and its widely applicable subset special relativity, as well as relativistic quantum mechanics and relativistic quantum field theory, the Lorentz transformation is the transformation rule under which all four-vectors and tensors containing physical quantities transform from one frame of reference to another.

The prime examples of such four-vectors are the four-position and four-momentum of a particle, and for fields the electromagnetic tensor and stress–energy tensor. The fact that these objects transform according to the Lorentz transformation is what mathematically defines them as vectors and tensors; see tensor for a definition.

Given the components of the four-vectors or tensors in some frame, the "transformation rule" allows one to determine the altered components of the same four-vectors or tensors in another frame, which could be boosted or accelerated, relative to the original frame. A "boost" should not be conflated with spatial translation, rather it's characterized by the relative velocity between frames. The transformation rule itself depends on the relative motion of the frames. In the simplest case of two inertial frames the relative velocity between enters the transformation rule. For rotating reference frames or general non-inertial reference frames, more parameters are needed, including the relative velocity (magnitude and direction), the rotation axis and angle turned through.

Historical background

Further information: History of Lorentz transformations

The usual treatment (e.g., Albert Einstein's original work) is based on the invariance of the speed of light. However, this is not necessarily the starting point: indeed (as is described, for example, in the second volume of the Course of Theoretical Physics by Landau and Lifshitz), what is really at stake is the locality of interactions: one supposes that the influence that one particle, say, exerts on another can not be transmitted instantaneously. Hence, there exists a theoretical maximal speed of information transmission which must be invariant, and it turns out that this speed coincides with the speed of light in vacuum. Newton had himself called the idea of action at a distance philosophically "absurd", and held that gravity had to be transmitted by some agent according to certain laws.

Michelson and Morley in 1887 designed an experiment, employing an interferometer and a half-silvered mirror, that was accurate enough to detect aether flow. The mirror system reflected the light back into the interferometer. If there were an aether drift, it would produce a phase shift and a change in the interference that would be detected. However, no phase shift was ever found. The negative outcome of the Michelson–Morley experiment left the concept of aether (or its drift) undermined. There was consequent perplexity as to why light evidently behaves like a wave, without any detectable medium through which wave activity might propagate.

In a 1964 paper, Erik Christopher Zeeman showed that the causality-preserving property, a condition that is weaker in a mathematical sense than the invariance of the speed of light, is enough to assure that the coordinate transformations are the Lorentz transformations. Norman Goldstein's paper shows a similar result using inertiality (the preservation of time-like lines) rather than causality.

Physical principles

Einstein based his theory of special relativity on two fundamental postulates. First, all physical laws are the same for all inertial frames of reference, regardless of their relative state of motion; and second, the speed of light in free space is the same in all inertial frames of reference, again, regardless of the relative velocity of each reference frame. The Lorentz transformation is fundamentally a direct consequence of this second postulate.

The second postulate

Assume the second postulate of special relativity stating the constancy of the speed of light, independent of reference frame, and consider a collection of reference systems moving with respect to each other with constant velocity, i.e. inertial systems, each endowed with its own set of Cartesian coordinates labeling the points, i.e. events of spacetime. To express the invariance of the speed of light in mathematical form, fix two events in spacetime, to be recorded in each reference frame. Let the first event be the emission of a light signal, and the second event be it being absorbed.

Pick any reference frame in the collection. In its coordinates, the first event will be assigned coordinates [math]\displaystyle{ x_1, y_1, z_1, ct_1 }[/math], and the second [math]\displaystyle{ x_2, y_2, z_2, ct_2 }[/math]. The spatial distance between emission and absorption is [math]\displaystyle{ \sqrt{(x_2-x_1)^2+(y_2-y_1)^2+(z_2-z_1)^2} }[/math], but this is also the distance [math]\displaystyle{ c(t_2-t_1) }[/math] traveled by the signal. One may therefore set up the equation

[math]\displaystyle{ c^2(t_2 - t_1)^2 - (x_2 - x_1)^2 - (y_2 - y_1)^2 - (z_2 - z_1)^2 = 0. }[/math]

Every other coordinate system will record, in its own coordinates, the same equation. This is the immediate mathematical consequence of the invariance of the speed of light. The quantity on the left is called the spacetime interval. The interval is, for events separated by light signals, the same (zero) in all reference frames, and is therefore called invariant.

Invariance of interval

For the Lorentz transformation to have the physical significance realized by nature, it is crucial that the interval is an invariant measure for any two events, not just for those separated by light signals. To establish this, one considers an infinitesimal interval,

[math]\displaystyle{ ds^2 = c^2 dt^2 - dx^2 - dy^2 - dz^2, }[/math]

as recorded in a system [math]\displaystyle{ K }[/math]. Let [math]\displaystyle{ K' }[/math] be another system assigning the interval [math]\displaystyle{ ds'^2 }[/math] to the same two infinitesimally separated events. Since if [math]\displaystyle{ ds^2 = 0 }[/math], then the interval will also be zero in any other system (second postulate), and since [math]\displaystyle{ ds^2 }[/math] and [math]\displaystyle{ ds'^2 }[/math] are infinitesimals of the same order, they must be proportional to each other, [math]\displaystyle{ ds^2 = a ds'^2. }[/math]

On what may [math]\displaystyle{ a }[/math] depend? It may not depend on the positions of the two events in spacetime, because that would violate the postulated homogeneity of spacetime. It might depend on the relative velocity [math]\displaystyle{ V' }[/math] between [math]\displaystyle{ K }[/math] and [math]\displaystyle{ K' }[/math], but only on the speed, not on the direction, because the latter would violate the isotropy of space.

Now bring in systems [math]\displaystyle{ K_1 }[/math] and [math]\displaystyle{ K_2 }[/math], [math]\displaystyle{ ds^2 = a(V_1)ds_1^2, \quad ds^2 = a(V_2)ds_2^2, \quad ds_1^2 = a(V_{12})ds_2^2. }[/math] From these it follows, [math]\displaystyle{ \frac{a(V_2)}{a(V_1)} = a(V_{12}). }[/math]

Now, one observes that on the right-hand side that [math]\displaystyle{ V_{12} }[/math] depend on both [math]\displaystyle{ V_1 }[/math] and [math]\displaystyle{ V_2 }[/math]; as well as on the angle between the vectors [math]\displaystyle{ \textbf{V}_1 }[/math] and [math]\displaystyle{ \textbf{V}_2 }[/math]. However, one also observes that the left-hand side does not depend on this angle. Thus, the only way for the equation to hold true is if the function [math]\displaystyle{ a(V) }[/math] is a constant. Further, by the same equation this constant is unity. Thus, [math]\displaystyle{ ds^2 = ds'^2 }[/math] for all systems [math]\displaystyle{ K' }[/math]. Since this holds for all infinitesimal intervals, it holds for all intervals.

Most, if not all, derivations of the Lorentz transformations take this for granted. In those derivations, they use the constancy of the speed of light (invariance of light-like separated events) only. This result ensures that the Lorentz transformation is the correct transformation.

Rigorous Statement and Proof of Proportionality of ds² and ds′²

Theorem: Let [math]\displaystyle{ n,p\geq 1 }[/math] be integers, [math]\displaystyle{ d:= n+p }[/math] and [math]\displaystyle{ V }[/math] a vector space over [math]\displaystyle{ \Reals }[/math] of dimension [math]\displaystyle{ d }[/math]. Let [math]\displaystyle{ h }[/math] be an indefinite-inner product on [math]\displaystyle{ V }[/math] with signature type [math]\displaystyle{ (n,p) }[/math]. Suppose [math]\displaystyle{ g }[/math] is a symmetric bilinear form on [math]\displaystyle{ V }[/math] such that the null set of the associated quadratic form of [math]\displaystyle{ h }[/math] is contained in that of [math]\displaystyle{ g }[/math] (i.e. suppose that for every [math]\displaystyle{ v \in V }[/math], if [math]\displaystyle{ h(v,v) = 0 }[/math] then [math]\displaystyle{ g(v,v)=0 }[/math]). Then, there exists a constant [math]\displaystyle{ C\in\Reals }[/math] such that [math]\displaystyle{ g = Ch }[/math]. Furthermore, if we assume [math]\displaystyle{ n\neq p }[/math] and that [math]\displaystyle{ g }[/math] also has signature type [math]\displaystyle{ (n,p) }[/math], then we have [math]\displaystyle{ C\gt 0 }[/math].

Remarks.

• In the section above, the term "infinitesimal" in relation to [math]\displaystyle{ ds^2 }[/math] is actually referring (pointwise) to a quadratic form over a four-dimensional real vector space (namely the tangent space at a point of the spacetime manifold). The argument above is copied almost verbatim from Landau and Lifshitz, where the proportionality of [math]\displaystyle{ ds^2 }[/math] and [math]\displaystyle{ ds'^2 }[/math] is merely stated as an 'obvious' fact even though the statement is not formulated in a mathematically precise fashion nor proven. This is a non-obvious mathematical fact which needs to be justified; fortunately the proof is relatively simple and it amounts to basic algebraic observations and manipulations.
• The above assumptions on [math]\displaystyle{ h }[/math] means the following: [math]\displaystyle{ h:V\times V\to\Reals }[/math] is a bilinear form which is symmetric and non-degenerate, such that there exists an ordered basis [math]\displaystyle{ \{v_1,\dots, v_n,v_{n+1},\dots, v_d\} }[/math] of [math]\displaystyle{ V }[/math] for which [math]\displaystyle{ h(v_a,v_b) = \begin{cases} -1 & \text{if } a = b \in \{1,\dots, n\}\\ 1 & \text{if } a = b \in \{n+1,\dots, d\}\\ 0&\text{ otherwise} \end{cases} }[/math] An equivalent way of saying this is that [math]\displaystyle{ h }[/math] has the matrix representation [math]\displaystyle{ \begin{pmatrix} -I_n & 0 \\0 & I_p\end{pmatrix} }[/math] relative to the ordered basis [math]\displaystyle{ \{v_1,\dots, v_d\} }[/math].
• If we consider the special case where [math]\displaystyle{ n = 1, p = 3 }[/math] then we're dealing with the situation of Lorentzian signature in 4-dimensions, which is what relativity is based on (or one could adopt the opposite convention with an overall minus sign; but this clearly doesn't affect the truth of the theorem). Also, in this case, if we assume [math]\displaystyle{ g }[/math] and [math]\displaystyle{ h }[/math] both have quadratics forms with the same null-set (in physics terminology, we say that [math]\displaystyle{ g }[/math] and [math]\displaystyle{ h }[/math] give rise to the same light cone) then the theorem tells us that there is a constant [math]\displaystyle{ C\gt 0 }[/math] such that [math]\displaystyle{ g = Ch }[/math]. Modulo some differences in notation, this is precisely what was used in the section above.

Proof of Theorem.

Fix a basis [math]\displaystyle{ \{v_1,\dots, v_d\} }[/math] of [math]\displaystyle{ V }[/math] relative to which [math]\displaystyle{ h }[/math] has the matrix representation [math]\displaystyle{ [h]= \begin{pmatrix} -I_n&0\\ 0&I_p \end{pmatrix} }[/math]. The point is that the vector space [math]\displaystyle{ V }[/math] can be decomposed into subspaces [math]\displaystyle{ V^- }[/math] (the span of the first [math]\displaystyle{ n }[/math] basis vectors) and [math]\displaystyle{ V^+ }[/math] (then span of the other [math]\displaystyle{ p }[/math] basis vectors) such that each vector in [math]\displaystyle{ V }[/math] can be written uniquely as [math]\displaystyle{ v + w }[/math] for [math]\displaystyle{ v \in V^- }[/math] and [math]\displaystyle{ w \in V^+ }[/math]; moreover [math]\displaystyle{ h(v,v) \leq 0 }[/math], [math]\displaystyle{ h(w,w) \geq 0 }[/math] and [math]\displaystyle{ h(v,w) = 0 }[/math]. So (by bilinearity)

[math]\displaystyle{ h(v+w,v+w) = h(v,v) + h(w,w) }[/math] Since the first summand on the right in non-positive and the second in non-negative, for any [math]\displaystyle{ v \in V^- }[/math] and [math]\displaystyle{ w \in V^+ }[/math], we can find a scalar [math]\displaystyle{ \alpha }[/math] such that [math]\displaystyle{ h(v + \alpha w, v + \alpha w) = 0 }[/math].

From now on, always consider [math]\displaystyle{ v \in V^- }[/math] and [math]\displaystyle{ w \in V^+ }[/math]. By bilinearity [math]\displaystyle{ \begin{align} g(v+w,v+w) &= g(v,v) + g(w,w) + 2g(v,w) \\ g(v-w,v-w) &= g(v,v) + g(w,w) - 2g(v,w) \end{align} }[/math]

If [math]\displaystyle{ h(v+w,v+w) = 0 }[/math], then also [math]\displaystyle{ h(v-w, v-w) = 0 }[/math] and the same is true for [math]\displaystyle{ g }[/math] (since the null-set of [math]\displaystyle{ h }[/math] is contained in that of [math]\displaystyle{ g }[/math]). In that case, subtracting the two expression above (and dividing by 4) yields

[math]\displaystyle{ 0 = g(v, w) }[/math]

As above, for each [math]\displaystyle{ v \in V^- }[/math] and [math]\displaystyle{ w \in V^+ }[/math], there is a scalar [math]\displaystyle{ \alpha }[/math] such that [math]\displaystyle{ h(v + \alpha w, v + \alpha w) = 0 }[/math], so [math]\displaystyle{ g(v, \alpha w) = 0 }[/math], which by bilinearity means [math]\displaystyle{ g(v,w) = 0 }[/math].

Now consider nonzero [math]\displaystyle{ v, v' \in V^- }[/math] such that [math]\displaystyle{ h(v, v) = h(v', v') }[/math]. We can find [math]\displaystyle{ w \in V^+ }[/math] such that [math]\displaystyle{ 0 = h(v + w, v + w) = h(v,v) + h(w,w) = h(v' + w, v' + w) }[/math]. By the expressions above, [math]\displaystyle{ g(v,v) = -g(w,w) = g(v', v') }[/math] Analogically, for [math]\displaystyle{ w, w' \in V^+ }[/math], one can show that if [math]\displaystyle{ h(w,w)=h(w',w') }[/math], then also [math]\displaystyle{ g(w,w)=g(w',w') }[/math]. So it holds for all vectors in [math]\displaystyle{ V }[/math].

For [math]\displaystyle{ u, u' \in V }[/math], if [math]\displaystyle{ g(u, u) = Ch(u, u) \neq 0 }[/math], [math]\displaystyle{ g(u', u') = C'h(u',u') \neq 0 }[/math] for some [math]\displaystyle{ C, C' \in \mathbb{R} }[/math], we can (scaling one of the if necessary) assume [math]\displaystyle{ h(u, u) = h(u',u') }[/math], which by the above means that [math]\displaystyle{ C = C' }[/math]. So [math]\displaystyle{ g = Ch }[/math].

Finally, if we assume that [math]\displaystyle{ g,h }[/math] both have signature types [math]\displaystyle{ (n,p) }[/math] and [math]\displaystyle{ n\neq p }[/math] then [math]\displaystyle{ C \gt 0 }[/math] (we can't have [math]\displaystyle{ C = 0 }[/math] because that would mean [math]\displaystyle{ g = 0 }[/math], which is impossible since having signature type [math]\displaystyle{ (n,p) }[/math] means it is a non-zero bilinear form. Also, if [math]\displaystyle{ C\lt 0 }[/math], then it means [math]\displaystyle{ g }[/math] has [math]\displaystyle{ n }[/math] positive diagonal entries and [math]\displaystyle{ p }[/math] negative diagonal entries; i.e. it is of signature [math]\displaystyle{ (p,n)\neq (n,p) }[/math], since we assumed [math]\displaystyle{ n\neq p }[/math], so this is also not possible. This leaves us with [math]\displaystyle{ C \gt 0 }[/math] as the only option). This completes the proof of the theorem.

Standard configuration

The spacetime coordinates of an event, as measured by each observer in their inertial reference frame (in standard configuration) are shown in the speech bubbles.
Top: frame F′ moves at velocity v along the x-axis of frame F.
Bottom: frame F moves at velocity −v along the x′-axis of frame F′.

The invariant interval can be seen as a non-positive definite distance function on spacetime. The set of transformations sought must leave this distance invariant. Due to the reference frame's coordinate system's cartesian nature, one concludes that, as in the Euclidean case, the possible transformations are made up of translations and rotations, where a slightly broader meaning should be allowed for the term rotation.

The interval is quite trivially invariant under translation. For rotations, there are four coordinates. Hence there are six planes of rotation. Three of those are rotations in spatial planes. The interval is invariant under ordinary rotations too.

It remains to find a "rotation" in the three remaining coordinate planes that leaves the interval invariant. Equivalently, to find a way to assign coordinates so that they coincide with the coordinates corresponding to a moving frame.

The general problem is to find a transformation such that [math]\displaystyle{ \begin{align} & c^2(t_2 - t_1)^2 - (x_2 - x_1)^2 - (y_2 - y_1)^2 - (z_2 - z_1)^2 \\ ={}& c^2(t_2' - t_1')^2 - (x_2' - x_1')^2 - (y_2' - y_1')^2 - (z_2' - z_1')^2. \end{align} }[/math]

To solve the general problem, one may use the knowledge about invariance of the interval of translations and ordinary rotations to assume, without loss of generality, that the frames F and F′ are aligned in such a way that their coordinate axes all meet at t = t′ = 0 and that the x and x′ axes are permanently aligned and system F′ has speed V along the positive x-axis. Call this the standard configuration. It reduces the general problem to finding a transformation such that

[math]\displaystyle{ c^2(t_2 - t_1)^2 - (x_2 - x_1)^2 = c^2(t_2' - t_1')^2 - (x_2' - x_1')^2. }[/math]

The standard configuration is used in most examples below. A linear solution of the simpler problem

[math]\displaystyle{ (ct)^2 - x^2 = (ct')^2 - x'^2 }[/math]

solves the more general problem since coordinate differences then transform the same way. Linearity is often assumed or argued somehow in the literature when this simpler problem is considered. If the solution to the simpler problem is not linear, then it doesn't solve the original problem because of the cross terms appearing when expanding the squares.

The solutions

As mentioned, the general problem is solved by translations in spacetime. These do not appear as a solution to the simpler problem posed, while the boosts do (and sometimes rotations depending on angle of attack). Even more solutions exist if one only insist on invariance of the interval for lightlike separated events. These are nonlinear conformal ("angle preserving") transformations. One has

Lorentz transformations ⊂ Poincaré transformations ⊂ conformal group transformations.

Some equations of physics are conformal invariant, e.g. the Maxwell's equations in source-free space, but not all. The relevance of the conformal transformations in spacetime is not known at present, but the conformal group in two dimensions is highly relevant in conformal field theory and statistical mechanics. It is thus the Poincaré group that is singled out by the postulates of special relativity. It is the presence of Lorentz boosts (for which velocity addition is different from mere vector addition that would allow for speeds greater than the speed of light) as opposed to ordinary boosts that separates it from the Galilean group of Galilean relativity. Spatial rotations, spatial and temporal inversions and translations are present in both groups and have the same consequences in both theories (conservation laws of momentum, energy, and angular momentum). Not all accepted theories respect symmetry under the inversions.

Using the geometry of spacetime

Landau & Lifshitz solution

These three hyperbolic function formulae (H1–H3) are referenced below:

1. [math]\displaystyle{ \cosh^2\Psi - \sinh^2 \Psi = 1, }[/math]
2. [math]\displaystyle{ \sinh \Psi= \frac{\tanh\Psi}{\sqrt{1-\tanh^2\Psi}}, }[/math]
3. [math]\displaystyle{ \cosh\Psi = \frac{1}{\sqrt{1-\tanh^2\Psi}}, }[/math]

The problem posed in standard configuration for a boost in the x-direction, where the primed coordinates refer to the moving system is solved by finding a linear solution to the simpler problem

[math]\displaystyle{ (ct)^2 - x^2 = (ct')^2 - x'^2. }[/math]

The most general solution is, as can be verified by direct substitution using (H1),

[math]\displaystyle{ x = x'\cosh\Psi + ct'\sinh \Psi, \quad ct = x'\sinh \Psi + ct'\cosh \Psi. }[/math]

(1)

To find the role of Ψ in the physical setting, record the progression of the origin of F′, i.e. x′ = 0, x = vt. The equations become (using first x′ = 0),

[math]\displaystyle{ x = ct'\sinh \Psi, \quad ct = ct'\cosh \Psi. }[/math]

Now divide:

[math]\displaystyle{ \frac{x}{ct} = \tanh \Psi = \frac{v}{c} \Rightarrow \quad \sinh \Psi= \frac{\frac{v}{c}}{\sqrt{1-{\frac{v^2}{c^2}}}}, \quad \cosh\Psi = \frac{1}{\sqrt{1-{\frac{v^2}{c^2}}}}, }[/math]

where x = vt was used in the first step, (H2) and (H3) in the second, which, when plugged back in (1), gives

[math]\displaystyle{ x=\frac{x'+vt'}{\sqrt{1-\frac{v^2}{c^2}}}, \quad t= \frac{t'+\frac{v}{c^2}x'}{\sqrt{1-\frac{v^2}{c^2}}}, }[/math]

or, with the usual abbreviations,

[math]\displaystyle{ x=\gamma(x'+vt'),\,\, t=\gamma\left(t'+\frac{vx'}{c^2}\right),\quad x'=\gamma(x-v t),\,\, t'=\gamma\left(t-\frac{vx}{c^2}\right). }[/math]

This calculation is repeated with more detail in section hyperbolic rotation.

Hyperbolic rotation

The Lorentz transformations can also be derived by simple application of the special relativity postulates and using hyperbolic identities.

Relativity postulates

Start from the equations of the spherical wave front of a light pulse, centred at the origin:

[math]\displaystyle{ (ct)^2 - (x^2+y^2+z^2) = (ct')^2 - (x'^2+y'^2+z'^2) =0 }[/math]

which take the same form in both frames because of the special relativity postulates. Next, consider relative motion along the x-axes of each frame, in standard configuration above, so that y = y′, z = z′, which simplifies to

[math]\displaystyle{ (ct)^2 - x^2 = (ct')^2 - x'^2 }[/math]

Linearity

Now assume that the transformations take the linear form:

[math]\displaystyle{ \begin{align} x' & = Ax + Bct \\ ct' & = Cx + Dct \end{align} }[/math]

where A, B, C, D are to be found. If they were non-linear, they would not take the same form for all observers, since fictitious forces (hence accelerations) would occur in one frame even if the velocity was constant in another, which is inconsistent with inertial frame transformations.

Substituting into the previous result:

[math]\displaystyle{ (ct)^2 - x^2 = [(Cx)^2 + (Dct)^2 + 2CDcxt] - [(Ax)^2 + (Bct)^2 + 2ABcxt] }[/math]

and comparing coefficients of x², t², xt:

[math]\displaystyle{ \begin{align} - 1 = C^2 - A^2 & \Rightarrow & A^2 - C^2 = 1 \\ c^2 = (Dc)^2 - (Bc)^2 & \Rightarrow & D^2 - B^2 = 1 \\ 2CDc - 2ABc = 0 & \Rightarrow & AB = CD \end{align} }[/math]

Hyperbolic rotation

The equations suggest the hyperbolic identity [math]\displaystyle{ \cosh^2\phi-\sinh^2\phi=1 . }[/math]

Introducing the rapidity parameter ϕ as a hyperbolic angle allows the consistent identifications

[math]\displaystyle{ A = D = \cosh\phi\,,\quad C = B = -\sinh\phi }[/math]

where the signs after the square roots are chosen so that x and t increase. The hyperbolic transformations have been solved for:

[math]\displaystyle{ \begin{align} x' & = x \cosh\phi - ct \sinh\phi \\ ct' & = -x \sinh\phi + ct \cosh\phi \end{align} }[/math]

If the signs were chosen differently the position and time coordinates would need to be replaced by −x and/or −t so that x and t increase not decrease.

To find how ϕ relates to the relative velocity, from the standard configuration the origin of the primed frame x′ = 0 is measured in the unprimed frame to be x = vt (or the equivalent and opposite way round; the origin of the unprimed frame is x = 0 and in the primed frame it is at x′ = −vt):

[math]\displaystyle{ 0 = vt \cosh\phi - ct \sinh\phi \, \Rightarrow \, \tanh\phi = \frac{v}{c} = \beta }[/math]

and hyperbolic identities [math]\displaystyle{ \sinh \Psi= \frac{\tanh\Psi}{\sqrt{1-\tanh^2\Psi}},\, \cosh\Psi = \frac{1}{\sqrt{1-\tanh^2\Psi}} }[/math] leads to the relations between β, γ, and ϕ,

[math]\displaystyle{ \cosh\phi = \gamma,\,\quad \sinh\phi=\beta\gamma \,. }[/math]

From physical principles

The problem is usually restricted to two dimensions by using a velocity along the x axis such that the y and z coordinates do not intervene, as described in standard configuration above.

Time dilation and length contraction

The transformation equations can be derived from time dilation and length contraction, which in turn can be derived from first principles. With O and O′ representing the spatial origins of the frames F and F′, and some event M, the relation between the position vectors (which here reduce to oriented segments OM, OO′ and O′M) in both frames is given by:

OM = OO′ + O′M.

Using coordinates (x,t) in F and (x′,t′) in F′ for event M, in frame F the segments are OM = x, OO′ = vt and O′M = x′/γ (since x′ is O′M as measured in F′): [math]\displaystyle{ x = vt + x'/\gamma. }[/math] Likewise, in frame F′, the segments are OM = x/γ (since x is OM as measured in F), OO′ = vt′ and O′M = x′: [math]\displaystyle{ x / \gamma = vt' + x'. }[/math] By rearranging the first equation, we get [math]\displaystyle{ x' = \gamma(x - vt), }[/math] which is the space part of the Lorentz transformation. The second relation gives [math]\displaystyle{ x = \gamma(x' + vt'), }[/math] which is the inverse of the space part. Eliminating x′ between the two space part equations gives

[math]\displaystyle{ t' = \gamma t + \frac{ \left( 1 - { \gamma^2} \right)x}{ \gamma v}. }[/math]

that, if [math]\displaystyle{ \gamma^2=\frac{1}{1-v^2/c^2} }[/math], simplifies to:

[math]\displaystyle{ t' = \gamma(t-vx/c^2), }[/math] which is the time part of the transformation, the inverse of which is found by a similar elimination of x: [math]\displaystyle{ t = \gamma(t' + vx'/c^2). }[/math]

Spherical wavefronts of light

The following is similar to that of Einstein. As in the Galilean transformation, the Lorentz transformation is linear since the relative velocity of the reference frames is constant as a vector; otherwise, inertial forces would appear. They are called inertial or Galilean reference frames. According to relativity no Galilean reference frame is privileged. Another condition is that the speed of light must be independent of the reference frame, in practice of the velocity of the light source.

Consider two inertial frames of reference O and O′, assuming O to be at rest while O′ is moving with a velocity v with respect to O in the positive x-direction. The origins of O and O′ initially coincide with each other. A light signal is emitted from the common origin and travels as a spherical wave front. Consider a point P on a spherical wavefront at a distance r and r′ from the origins of O and O′ respectively. According to the second postulate of the special theory of relativity the speed of light is the same in both frames, so for the point P: [math]\displaystyle{ \begin{align} r &= ct \\ r' &= ct'. \end{align} }[/math]

The equation of a sphere in frame O is given by [math]\displaystyle{ x^2 + y^2 + z^2 = r^2. }[/math] For the spherical wavefront that becomes [math]\displaystyle{ x^2 + y^2 + z^2 = (ct)^2. }[/math] Similarly, the equation of a sphere in frame O′ is given by [math]\displaystyle{ x'^2 + y'^2 + z'^2 = r'^2, }[/math] so the spherical wavefront satisfies [math]\displaystyle{ x'^2 + y'^2 + z'^2 = (ct')^2. }[/math]

The origin O′ is moving along x-axis. Therefore, [math]\displaystyle{ \begin{align} y' &= y \\ z' &= z. \end{align} }[/math]

x′ must vary linearly with x and t. Therefore, the transformation has the form [math]\displaystyle{ x' = \gamma x + \sigma t. }[/math] For the origin of O′ x′ and x are given by [math]\displaystyle{ \begin{align} x' &= 0 \\ x &= vt, \end{align} }[/math] so, for all t, [math]\displaystyle{ 0 = \gamma vt + \sigma t }[/math] and thus [math]\displaystyle{ \sigma = -\gamma v. }[/math] This simplifies the transformation to [math]\displaystyle{ x' = \gamma \left( x - v t \right) }[/math] where γ is to be determined. At this point γ is not necessarily a constant, but is required to reduce to 1 for v ≪ c.

The inverse transformation is the same except that the sign of v is reversed: [math]\displaystyle{ x = \gamma \left( x' + v t' \right). }[/math]

The above two equations give the relation between t and t′ as: [math]\displaystyle{ x = \gamma \left[ \gamma \left( x - v t \right) + v t' \right] }[/math] or [math]\displaystyle{ t' = \gamma t + \frac{ \left( 1 - { \gamma^2} \right)x}{ \gamma v}. }[/math]

Replacing x′, y′, z′ and t′ in the spherical wavefront equation in the O′ frame, [math]\displaystyle{ x'^2 + y'^2 + z'^2 = (ct')^2, }[/math] with their expressions in terms of x, y, z and t produces: [math]\displaystyle{ {\gamma^2} \left( x - v t \right)^2 + y^2 + z^2 = c^2 \left[ \gamma t + \frac{ \left( 1 - { \gamma^2} \right)x}{ \gamma v} \right]^2 }[/math] and therefore, [math]\displaystyle{ \gamma^2 x^2 + \gamma^2 v^2 t^2 - 2 \gamma^2 v t x + y^2 + z^2 = c^2 {\gamma^2} t^2 + \frac{ \left( 1 - {\gamma^2} \right)^2 c^2 x^2}{ {\gamma^2} v^2} + 2 \frac{ \left( 1 - {\gamma^2} \right) t x c^2}{ v} }[/math] which implies, [math]\displaystyle{ \left[ {\gamma^2} - \frac{ \left( 1 - {\gamma^2} \right)^2 c^2}{ {\gamma^2} v^2} \right] x^2 - 2 {\gamma^2} v t x + y^2 + z^2 = \left( c^2 {\gamma^2} - v^2 {\gamma^2} \right) t^2 + 2 \frac{ \left[ 1 - {\gamma^2} \right] t x c^2}{ v} }[/math] or [math]\displaystyle{ \left[ {\gamma^2} - \frac{ \left( 1 - {\gamma^2} \right)^2 c^2}{ {\gamma^2} v^2} \right] x^2 - \left[ 2 {\gamma^2} v + 2 \frac{ \left( 1 - {\gamma^2} \right) c^2}{ v}\right] t x + y^2 + z^2 = \left[ c^2 {\gamma^2} - v^2 {\gamma^2} \right] t^2 }[/math]

Comparing the coefficient of t² in the above equation with the coefficient of t² in the spherical wavefront equation for frame O produces: [math]\displaystyle{ c^2 {\gamma^2} - v^2 {\gamma^2} = c^2 }[/math] Equivalent expressions for γ can be obtained by matching the x² coefficients or setting the tx coefficient to zero. Rearranging: [math]\displaystyle{ {\gamma^2} = \frac{1}{1 - \frac{v^2}{c^2}} }[/math] or, choosing the positive root to ensure that the x and x' axes and the time axes point in the same direction, [math]\displaystyle{ {\gamma} = \frac{1}{ \sqrt{1 - \frac{v^2}{c^2}}} }[/math] which is called the Lorentz factor. This produces the Lorentz transformation from the above expression. It is given by [math]\displaystyle{ \begin{align} x' &= \gamma \left( x - v t \right)\\ t' &= \gamma \left( t - \frac{vx}{c^2} \right) \\ y' &= y \\ z' &= z \end{align} }[/math]

The Lorentz transformation is not the only transformation leaving invariant the shape of spherical waves, as there is a wider set of spherical wave transformations in the context of conformal geometry, leaving invariant the expression [math]\displaystyle{ \lambda\left(\delta x^2 + \delta y^2 + \delta z^2 - c^2\delta t^2\right) }[/math]. However, scale changing conformal transformations cannot be used to symmetrically describe all laws of nature including mechanics, whereas the Lorentz transformations (the only one implying [math]\displaystyle{ \lambda = 1 }[/math]) represent a symmetry of all laws of nature and reduce to Galilean transformations at [math]\displaystyle{ v \ll c }[/math].

Galilean and Einstein's relativity

Galilean reference frames

In classical kinematics, the total displacement x in the R frame is the sum of the relative displacement x′ in frame R′ and of the distance between the two origins x − x′. If v is the relative velocity of R′ relative to R, the transformation is: x = x′ + vt, or x′ = x − vt. This relationship is linear for a constant v, that is when R and R′ are Galilean frames of reference.

In Einstein's relativity, the main difference from Galilean relativity is that space and time coordinates are intertwined, and in different inertial frames t ≠ t′.

Since space is assumed to be homogeneous, the transformation must be linear. The most general linear relationship is obtained with four constant coefficients, A, B, γ, and b: [math]\displaystyle{ x'=\gamma x + b t }[/math] [math]\displaystyle{ t'=A x + B t. }[/math] The linear transformation becomes the Galilean transformation when γ = B = 1, b = −v and A = 0.

An object at rest in the R′ frame at position x′ = 0 moves with constant velocity v in the R frame. Hence the transformation must yield x′ = 0 if x = vt. Therefore, b = −γv and the first equation is written as [math]\displaystyle{ x'=\gamma \left(x - v t\right) . }[/math]

Using the principle of relativity

According to the principle of relativity, there is no privileged Galilean frame of reference: therefore the inverse transformation for the position from frame R′ to frame R should have the same form as the original but with the velocity in the opposite direction, i.o.w. replacing v with -v: [math]\displaystyle{ x=\gamma\left(x' - (-v)t'\right), }[/math] and thus [math]\displaystyle{ x=\gamma\left(x' + vt'\right) . }[/math]

Determining the constants of the first equation

Since the speed of light is the same in all frames of reference, for the case of a light signal, the transformation must guarantee that t = x/c when t′ = x′/c.

Substituting for t and t′ in the preceding equations gives: [math]\displaystyle{ x'= \gamma\left(1 - v/c\right) x , }[/math] [math]\displaystyle{ x= \gamma\left(1 + v/c\right) x'. }[/math] Multiplying these two equations together gives, [math]\displaystyle{ xx' = \gamma^2 \left(1 - v^2/c^2\right) xx'. }[/math] At any time after t = t′ = 0, xx′ is not zero, so dividing both sides of the equation by xx′ results in [math]\displaystyle{ \gamma=\frac{1}{\sqrt{1 - \frac{v^2}{c^2}}}, }[/math] which is called the "Lorentz factor".

When the transformation equations are required to satisfy the light signal equations in the form x = ct and x′ = ct′, by substituting the x and x'-values, the same technique produces the same expression for the Lorentz factor.

Determining the constants of the second equation

The transformation equation for time can be easily obtained by considering the special case of a light signal, again satisfying x = ct and x′ = ct′, by substituting term by term into the earlier obtained equation for the spatial coordinate [math]\displaystyle{ x'=\gamma (x - v t), \, }[/math] giving [math]\displaystyle{ ct'=\gamma \left(ct - \frac{v}{c} x\right), }[/math] so that [math]\displaystyle{ t'=\gamma \left(t - \frac{v}{c^2} x\right), }[/math] which, when identified with [math]\displaystyle{ t'=A x + B t, \, }[/math] determines the transformation coefficients A and B as [math]\displaystyle{ A= -\gamma v/c^2,\, }[/math] [math]\displaystyle{ B=\gamma. \, }[/math] So A and B are the unique constant coefficients necessary to preserve the constancy of the speed of light in the primed system of coordinates.

Einstein's popular derivation

In his popular book Einstein derived the Lorentz transformation by arguing that there must be two non-zero coupling constants λ and μ such that

[math]\displaystyle{ \begin{cases} x' - ct' = \lambda \left( x - ct \right) \\ x' + ct' = \mu \left( x + ct \right) \, \end{cases} }[/math]

that correspond to light traveling along the positive and negative x-axis, respectively. For light x = ct if and only if x′ = ct′. Adding and subtracting the two equations and defining

[math]\displaystyle{ \begin{cases} \gamma = \left( \lambda + \mu \right) /2 \\ b = \left( \lambda - \mu \right) /2 ,\, \end{cases} }[/math]

gives

[math]\displaystyle{ \begin{cases} x' = \gamma x - bct \\ ct' = \gamma ct - bx . \, \end{cases} }[/math]

Substituting x′ = 0 corresponding to x = vt and noting that the relative velocity is v = bc/γ, this gives

[math]\displaystyle{ \begin{cases} x' = \gamma \left( x - vt \right) \\ t' = \gamma \left( t - \frac{v}{c^2} x \right) \, \end{cases} }[/math]

The constant γ can be evaluated by demanding c²t² − x² = c²t′² − x′² as per standard configuration.

Using group theory

From group postulates

Following is a classical derivation (see, e.g., [1] and references therein) based on group postulates and isotropy of the space.

Coordinate transformations as a group

The coordinate transformations between inertial frames form a group (called the proper Lorentz group) with the group operation being the composition of transformations (performing one transformation after another). Indeed, the four group axioms are satisfied:

1. Closure: the composition of two transformations is a transformation: consider a composition of transformations from the inertial frame K to inertial frame K′, (denoted as K → K′), and then from K′ to inertial frame K′′, [K′ → K′′], there exists a transformation, [K → K′] [K′ → K′′], directly from an inertial frame K to inertial frame K′′.
2. Associativity: the transformations ( [K → K′] [K′ → K′′] ) [K′′ → K′′′] and [K → K′] ( [K′ → K′′] [K′′ → K′′′] ) are identical.
3. Identity element: there is an identity element, a transformation K → K.
4. Inverse element: for any transformation K → K′ there exists an inverse transformation K′ → K.

Transformation matrices consistent with group axioms

Consider two inertial frames, K and K′, the latter moving with velocity v with respect to the former. By rotations and shifts we can choose the x and x′ axes along the relative velocity vector and also that the events (t, x) = (0,0) and (t′, x′) = (0,0) coincide. Since the velocity boost is along the x (and x′) axes nothing happens to the perpendicular coordinates and we can just omit them for brevity. Now since the transformation we are looking after connects two inertial frames, it has to transform a linear motion in (t, x) into a linear motion in (t′, x′) coordinates. Therefore, it must be a linear transformation. The general form of a linear transformation is [math]\displaystyle{ \begin{bmatrix} t' \\ x' \end{bmatrix} = \begin{bmatrix} \gamma & \delta \\ \beta & \alpha \end{bmatrix} \begin{bmatrix} t \\ x \end{bmatrix}, }[/math] where α, β, γ and δ are some yet unknown functions of the relative velocity v.

Let us now consider the motion of the origin of the frame K′. In the K′ frame it has coordinates (t′, x′ = 0), while in the K frame it has coordinates (t, x = vt). These two points are connected by the transformation [math]\displaystyle{ \begin{bmatrix} t' \\ 0 \end{bmatrix} = \begin{bmatrix} \gamma & \delta \\ \beta & \alpha \end{bmatrix} \begin{bmatrix} t \\ vt \end{bmatrix}, }[/math] from which we get [math]\displaystyle{ \beta = -v\alpha \,. }[/math] Analogously, considering the motion of the origin of the frame K, we get [math]\displaystyle{ \begin{bmatrix} t' \\ -vt' \end{bmatrix} = \begin{bmatrix} \gamma & \delta \\ \beta & \alpha \end{bmatrix} \begin{bmatrix} t \\ 0 \end{bmatrix}, }[/math] from which we get [math]\displaystyle{ \beta=-v\gamma \,. }[/math] Combining these two gives α = γ and the transformation matrix has simplified, [math]\displaystyle{ \begin{bmatrix} t' \\ x' \end{bmatrix} = \begin{bmatrix} \gamma & \delta \\ -v\gamma & \gamma \end{bmatrix} \begin{bmatrix} t \\ x \end{bmatrix}. }[/math]

Now consider the group postulate inverse element. There are two ways we can go from the K′ coordinate system to the K coordinate system. The first is to apply the inverse of the transform matrix to the K′ coordinates:

[math]\displaystyle{ \begin{bmatrix} t \\ x \end{bmatrix} = \frac{1}{\gamma^2+v\delta\gamma} \begin{bmatrix} \gamma & -\delta \\ v\gamma & \gamma \end{bmatrix} \begin{bmatrix} t' \\ x' \end{bmatrix}. }[/math]

The second is, considering that the K′ coordinate system is moving at a velocity v relative to the K coordinate system, the K coordinate system must be moving at a velocity −v relative to the K′ coordinate system. Replacing v with −v in the transformation matrix gives:

[math]\displaystyle{ \begin{bmatrix} t \\ x \end{bmatrix} = \begin{bmatrix} \gamma(-v) & \delta(-v) \\ v\gamma(-v) & \gamma(-v) \end{bmatrix} \begin{bmatrix} t' \\ x' \end{bmatrix}, }[/math]

Now the function γ can not depend upon the direction of v because it is apparently the factor which defines the relativistic contraction and time dilation. These two (in an isotropic world of ours) cannot depend upon the direction of v. Thus, γ(−v) = γ(v) and comparing the two matrices, we get [math]\displaystyle{ \gamma^2 + v\delta\gamma = 1. }[/math]

According to the closure group postulate a composition of two coordinate transformations is also a coordinate transformation, thus the product of two of our matrices should also be a matrix of the same form. Transforming K to K′ and from K′ to K′′ gives the following transformation matrix to go from K to K′′:

[math]\displaystyle{ \begin{align} \begin{bmatrix} t'' \\ x'' \end{bmatrix} & = \begin{bmatrix} \gamma(v') & \delta(v') \\ -v'\gamma(v') & \gamma(v') \end{bmatrix} \begin{bmatrix} \gamma(v) & \delta(v) \\ -v\gamma(v) & \gamma(v) \end{bmatrix} \begin{bmatrix} t \\ x \end{bmatrix}\\ & = \begin{bmatrix} \gamma(v')\gamma(v)-v\delta(v')\gamma(v) & \gamma(v')\delta(v)+\delta(v')\gamma(v) \\ -(v'+v)\gamma(v')\gamma(v) & -v'\gamma(v')\delta(v)+\gamma(v')\gamma(v) \end{bmatrix} \begin{bmatrix} t\\x \end{bmatrix}. \end{align} }[/math]

In the original transform matrix, the main diagonal elements are both equal to γ, hence, for the combined transform matrix above to be of the same form as the original transform matrix, the main diagonal elements must also be equal. Equating these elements and rearranging gives:

[math]\displaystyle{ \begin{align} \gamma(v')\gamma(v)-v\delta(v')\gamma(v) &= -v'\gamma(v')\delta(v)+\gamma(v')\gamma(v) \\ v\delta(v')\gamma(v) &= v'\gamma(v')\delta(v) \\ \frac{\delta(v)}{v\gamma(v)} &= \frac{\delta(v')}{v'\gamma(v')}. \end{align} }[/math]

The denominator will be nonzero for nonzero v, because γ(v) is always nonzero;

[math]\displaystyle{ \gamma^2 + v \delta \gamma = 1. }[/math]

If v = 0 we have the identity matrix which coincides with putting v = 0 in the matrix we get at the end of this derivation for the other values of v, making the final matrix valid for all nonnegative v.

For the nonzero v, this combination of function must be a universal constant, one and the same for all inertial frames. Define this constant as δ(v)/v γ(v) = κ, where κ has the dimension of 1/v². Solving [math]\displaystyle{ 1 = \gamma^2 + v\delta\gamma = \gamma^2 (1 + \kappa v^2) }[/math] we finally get [math]\displaystyle{ \gamma = 1/\sqrt{1 + \kappa v^2} }[/math] and thus the transformation matrix, consistent with the group axioms, is given by

[math]\displaystyle{ \begin{bmatrix} t' \\ x' \end{bmatrix} = \frac{1}{\sqrt{1 + \kappa v^2}} \begin{bmatrix} 1 & \kappa v \\ -v & 1 \end{bmatrix} \begin{bmatrix} t \\ x \end{bmatrix}. }[/math]

If κ > 0, then there would be transformations (with κv² ≫ 1) which transform time into a spatial coordinate and vice versa. We exclude this on physical grounds, because time can only run in the positive direction. Thus two types of transformation matrices are consistent with group postulates:

1. with the universal constant κ = 0, and
2. with κ < 0.

Galilean transformations

If κ = 0 then we get the Galilean-Newtonian kinematics with the Galilean transformation, [math]\displaystyle{ \begin{bmatrix} t' \\ x' \end{bmatrix} = \begin{bmatrix} 1 & 0 \\ -v & 1 \end{bmatrix} \begin{bmatrix} t \\ x \end{bmatrix}\;, }[/math] where time is absolute, t′ = t, and the relative velocity v of two inertial frames is not limited.

Lorentz transformations

If κ < 0, then we set [math]\displaystyle{ c = 1/\sqrt{-\kappa} }[/math] which becomes the invariant speed, the speed of light in vacuum. This yields κ = −1/c² and thus we get special relativity with Lorentz transformation [math]\displaystyle{ \begin{bmatrix} t' \\ x' \end{bmatrix} = \frac{1}{\sqrt{1 - {v^2 \over c^2}}} \begin{bmatrix} 1 & {- v \over c^2} \\ -v & 1 \end{bmatrix} \begin{bmatrix} t \\ x \end{bmatrix}\;, }[/math] where the speed of light is a finite universal constant determining the highest possible relative velocity between inertial frames.

If v ≪ c the Galilean transformation is a good approximation to the Lorentz transformation.

Only experiment can answer the question which of the two possibilities, κ = 0 or κ < 0, is realized in our world. The experiments measuring the speed of light, first performed by a Danish physicist Ole Rømer, show that it is finite, and the Michelson–Morley experiment showed that it is an absolute speed, and thus that κ < 0.

Boost from generators

Using rapidity ϕ to parametrize the Lorentz transformation, the boost in the x direction is

[math]\displaystyle{ \begin{bmatrix} c t' \\ x' \\ y' \\ z' \end{bmatrix} = \begin{bmatrix} \cosh\phi & -\sinh\phi & 0 & 0\\ -\sinh\phi & \cosh\phi & 0 & 0\\ 0&0&1&0\\ 0&0&0&1\\ \end{bmatrix} \begin{bmatrix} c\,t \\ x \\ y \\ z \end{bmatrix} , }[/math]

likewise for a boost in the y-direction

[math]\displaystyle{ \begin{bmatrix} c t' \\ x' \\ y' \\ z' \end{bmatrix} = \begin{bmatrix} \cosh\phi & 0 & -\sinh\phi & 0\\ 0&1&0&0\\ -\sinh\phi &0& \cosh\phi &0\\ 0&0&0&1\\ \end{bmatrix} \begin{bmatrix} c\,t \\ x \\ y \\ z \end{bmatrix} , }[/math]

and the z-direction

[math]\displaystyle{ \begin{bmatrix} c t' \\ x' \\ y' \\ z' \end{bmatrix} = \begin{bmatrix} \cosh\phi &0&0&-\sinh\phi\\ 0&1&0&0\\ 0&0&1&0\\ -\sinh\phi&0&0&\cosh\phi\\ \end{bmatrix} \begin{bmatrix} c\,t \\ x \\ y \\ z \end{bmatrix} \,. }[/math]

where e_x, e_y, e_z are the Cartesian basis vectors, a set of mutually perpendicular unit vectors along their indicated directions. If one frame is boosted with velocity v relative to another, it is convenient to introduce a unit vector n = v/v = β/β in the direction of relative motion. The general boost is

[math]\displaystyle{ \begin{bmatrix} c\,t' \\ x' \\ y' \\ z' \end{bmatrix} = \begin{bmatrix} \cosh\phi & -n_x\sinh\phi & -n_y\sinh\phi & -n_z\sinh\phi\\ -n_x\sinh\phi & 1+(\cosh\phi-1)n_x^2&(\cosh\phi-1)n_x n_y&(\cosh\phi-1)n_x n_z\\ -n_y\sinh\phi & (\cosh\phi-1)n_y n_x&1+(\cosh\phi-1)n_y^2&(\cosh\phi-1)n_y n_z\\ -n_z\sinh\phi & (\cosh\phi-1)n_z n_x&(\cosh\phi-1)n_z n_y&1+(\cosh\phi-1)n_z^2\\ \end{bmatrix} \begin{bmatrix} c\,t \\ x \\ y \\ z \end{bmatrix}\,. }[/math]

Notice the matrix depends on the direction of the relative motion as well as the rapidity, in all three numbers (two for direction, one for rapidity).

We can cast each of the boost matrices in another form as follows. First consider the boost in the x direction. The Taylor expansion of the boost matrix about ϕ = 0 is

[math]\displaystyle{ B(\mathbf{e}_x,\phi)=\sum_{n=0}^{\infty}\frac{\phi^n}{n!}\left.\frac{\partial^nB(\mathbf{e}_x,\phi)}{\partial\phi^n}\right|_{\phi=0} }[/math]

where the derivatives of the matrix with respect to ϕ are given by differentiating each entry of the matrix separately, and the notation |_{ϕ = 0} indicates ϕ is set to zero after the derivatives are evaluated. Expanding to first order gives the infinitesimal transformation

[math]\displaystyle{ B(\mathbf{e}_x,\phi)=I+\phi\left.\frac{\partial B}{\partial\phi}\right|_{\phi=0} = \begin{bmatrix} 1 & 0 & 0 & 0\\ 0 & 1 & 0 & 0\\ 0 & 0 & 1 & 0\\ 0 & 0 & 0 & 1 \end{bmatrix}-\phi \begin{bmatrix} 0 & 1 & 0 & 0\\ 1 & 0 & 0 & 0\\ 0 & 0 & 0 & 0\\ 0 & 0 & 0 & 0 \end{bmatrix} }[/math]

which is valid if ϕ is small (hence ϕ² and higher powers are negligible), and can be interpreted as no boost (the first term I is the 4×4 identity matrix), followed by a small boost. The matrix

[math]\displaystyle{ K_x=\begin{bmatrix}0 & 1 & 0 & 0\\ 1 & 0 & 0 & 0\\ 0 & 0 & 0 & 0\\ 0 & 0 & 0 & 0 \end{bmatrix} }[/math]

is the generator of the boost in the x direction, so the infinitesimal boost is

[math]\displaystyle{ B(\mathbf{e}_x,\phi)=I-\phi K_x }[/math]

Now, ϕ is small, so dividing by a positive integer N gives an even smaller increment of rapidity ϕ/N, and N of these infinitesimal boosts will give the original infinitesimal boost with rapidity ϕ,

[math]\displaystyle{ B(\mathbf{e}_x,\phi)=\left(I-\frac{\phi K_x}{N}\right)^N }[/math]

In the limit of an infinite number of infinitely small steps, we obtain the finite boost transformation

[math]\displaystyle{ B(\mathbf{e}_x,\phi)=\lim_{N\to\infty}\left(I-\frac{\phi K_x}{N}\right)^{N}=e^{-\phi K_x} }[/math]

which is the limit definition of the exponential due to Leonhard Euler, and is now true for any ϕ.

Repeating the process for the boosts in the y and z directions obtains the other generators

[math]\displaystyle{ K_y=\begin{bmatrix}0 & 0 & 1 & 0\\ 0 & 0 & 0 & 0\\ 1 & 0 & 0 & 0\\ 0 & 0 & 0 & 0 \end{bmatrix}\,,\quad K_z=\begin{bmatrix}0 & 0 & 0 & 1\\ 0 & 0 & 0 & 0\\ 0 & 0 & 0 & 0\\ 1 & 0 & 0 & 0 \end{bmatrix} }[/math]

and the boosts are

[math]\displaystyle{ B(\mathbf{e}_y,\phi)=e^{-\phi K_y}\,,\quad B(\mathbf{e}_z,\phi)=e^{-\phi K_z}\,. }[/math]

For any direction, the infinitesimal transformation is (small ϕ and expansion to first order)

[math]\displaystyle{ B(\mathbf{n},\phi)=I+\phi\left.\frac{\partial B}{\partial\phi}\right|_{\phi=0} = \begin{bmatrix} 1 & 0 & 0 & 0\\ 0 & 1 & 0 & 0\\ 0 & 0 & 1 & 0\\ 0 & 0 & 0 & 1 \end{bmatrix}-\phi\begin{bmatrix}0 & n_x & n_y & n_z\\ n_x & 0 & 0 & 0\\ n_y & 0 & 0 & 0\\ n_z & 0 & 0 & 0 \end{bmatrix} }[/math]

where

[math]\displaystyle{ \begin{bmatrix}0 & n_x & n_y & n_z\\ n_x & 0 & 0 & 0\\ n_y & 0 & 0 & 0\\ n_z & 0 & 0 & 0 \end{bmatrix} = n_x K_x + n_y K_y + n_z K_z = \mathbf{n}\cdot\mathbf{K} }[/math]

is the generator of the boost in direction n. It is the full boost generator, a vector of matrices K = (K_x, K_y, K_z), projected into the direction of the boost n. The infinitesimal boost is

[math]\displaystyle{ B(\mathbf{n},\phi)=I-\phi(\mathbf{n}\cdot\mathbf{K}) }[/math]

Then in the limit of an infinite number of infinitely small steps, we obtain the finite boost transformation

[math]\displaystyle{ B(\mathbf{n},\phi)=\lim_{N\to\infty}\left(I-\frac{\phi(\mathbf{n}\cdot\mathbf{K})}{N}\right)^{N}=e^{-\phi(\mathbf{n}\cdot\mathbf{K})} }[/math]

which is now true for any ϕ. Expanding the matrix exponential of −ϕ(n ⋅ K) in its power series

[math]\displaystyle{ e^{-\phi\mathbf{n}\cdot\mathbf{K}} = \sum_{n=0}^{\infty}\frac{1}{n!}(-\phi \mathbf{n}\cdot\mathbf{K})^n }[/math]

we now need the powers of the generator. The square is

[math]\displaystyle{ (\mathbf{n}\cdot\mathbf{K})^2=\begin{bmatrix}1 & 0 & 0 & 0\\ 0 & n_x^2 & n_x n_y & n_x n_z\\ 0 & n_y n_x & n_y^2 & n_y n_z\\ 0 & n_z n_x & n_z n_y & n_z^2 \end{bmatrix} }[/math]

but the cube (n ⋅ K)³ returns to (n ⋅ K), and as always the zeroth power is the 4×4 identity, (n ⋅ K)⁰ = I. In general the odd powers n = 1, 3, 5, ... are

[math]\displaystyle{ (\mathbf{n}\cdot\mathbf{K})^n = (\mathbf{n}\cdot\mathbf{K}) }[/math]

while the even powers n = 2, 4, 6, ... are

[math]\displaystyle{ (\mathbf{n}\cdot\mathbf{K})^n = (\mathbf{n}\cdot\mathbf{K})^2 }[/math]

therefore the explicit form of the boost matrix depends only the generator and its square. Splitting the power series into an odd power series and an even power series, using the odd and even powers of the generator, and the Taylor series of sinh ϕ and cosh ϕ about ϕ = 0 obtains a more compact but detailed form of the boost matrix

[math]\displaystyle{ \begin{align} e^{-\phi\mathbf{n}\cdot\mathbf{K}}& = -\sum_{n=1,3,5\ldots}^{\infty}\frac{1}{n!}\phi^n(\mathbf{n}\cdot\mathbf{K})^n+\sum_{n=0,2,4\ldots}^{\infty} \frac{1}{n!}\phi^n(\mathbf{n}\cdot\mathbf{K})^n\\ & = - \left[\phi+\frac{\phi^3}{3!}+\frac{\phi^5}{5!} +\cdots \right] (\mathbf{n}\cdot\mathbf{K}) + I +\left[ -1 + 1 + \frac{1}{2!} \phi^2+\frac{1}{4!}\phi^4 +\frac{1}{6!}\phi^6 +\cdots\right](\mathbf{n}\cdot\mathbf{K})^2 \\ & = - \sinh\phi(\mathbf{n}\cdot\mathbf{K}) + I +(-1+\cosh\phi)(\mathbf{n}\cdot\mathbf{K})^2 \end{align} }[/math]

where 0 = −1 + 1 is introduced for the even power series to complete the Taylor series for cosh ϕ. The boost is similar to Rodrigues' rotation formula,

[math]\displaystyle{ B(\mathbf{n},\phi) = e^{-\phi\mathbf{n}\cdot\mathbf{K}} = I -\sinh\phi(\mathbf{n}\cdot\mathbf{K})+(\cosh\phi-1)(\mathbf{n}\cdot\mathbf{K})^2 \,. }[/math]

Negating the rapidity in the exponential gives the inverse transformation matrix,

[math]\displaystyle{ B(\mathbf{n},-\phi) = e^{\phi\mathbf{n}\cdot\mathbf{K}} = I +\sinh\phi(\mathbf{n}\cdot\mathbf{K})+(\cosh\phi-1)(\mathbf{n}\cdot\mathbf{K})^2 \,. }[/math]

In quantum mechanics, relativistic quantum mechanics, and quantum field theory, a different convention is used for the boost generators; all of the boost generators are multiplied by a factor of the imaginary unit i = √−1.

From experiments

Howard Percy Robertson and others showed that the Lorentz transformation can also be derived empirically. In order to achieve this, it's necessary to write down coordinate transformations that include experimentally testable parameters. For instance, let there be given a single "preferred" inertial frame [math]\displaystyle{ X, Y, Z, T }[/math] in which the speed of light is constant, isotropic, and independent of the velocity of the source. It is also assumed that Einstein synchronization and synchronization by slow clock transport are equivalent in this frame. Then assume another frame [math]\displaystyle{ x, y, z, t }[/math] in relative motion, in which clocks and rods have the same internal constitution as in the preferred frame. The following relations, however, are left undefined:

• [math]\displaystyle{ a(v) }[/math] differences in time measurements,
• [math]\displaystyle{ b(v) }[/math] differences in measured longitudinal lengths,
• [math]\displaystyle{ d(v) }[/math] differences in measured transverse lengths,
• [math]\displaystyle{ \varepsilon(v) }[/math] depends on the clock synchronization procedure in the moving frame,

then the transformation formulas (assumed to be linear) between those frames are given by:

[math]\displaystyle{ \begin{align} t & =a(v)T+\varepsilon(v) x\\ x & =b(v)(X-vT)\\ y & =d(v)Y\\ z & =d(v)Z \end{align} }[/math]

[math]\displaystyle{ \varepsilon (v) }[/math] depends on the synchronization convention and is not determined experimentally, it obtains the value [math]\displaystyle{ -v/c^{2} }[/math] by using Einstein synchronization in both frames. The ratio between [math]\displaystyle{ b(v) }[/math] and [math]\displaystyle{ d(v) }[/math] is determined by the Michelson–Morley experiment, the ratio between [math]\displaystyle{ a(v) }[/math] and [math]\displaystyle{ b(v) }[/math] is determined by the Kennedy–Thorndike experiment, and [math]\displaystyle{ a(v) }[/math] alone is determined by the Ives–Stilwell experiment. In this way, they have been determined with great precision to [math]\displaystyle{ 1/a(v)=b(v)=\gamma }[/math] and [math]\displaystyle{ d(v)=1 }[/math], which converts the above transformation into the Lorentz transformation.

See also

• Lorentz group
• Noether's theorem
• Poincaré group
• Proper time
• Relativistic metric
• Spinor

References

• Greiner, W.; Bromley, D. A. (2000). Relativistic Quantum Mechanics (3rd ed.). springer. ISBN 9783540674573. https://books.google.com/books?id=a6_rCAAAQBAJ&q=walter+greiner+theoretical+physics&pg=PR5. 
• Landau, L.D.; Lifshitz, E.M. (2002). The Classical Theory of Fields. Course of Theoretical Physics. 2 (4th ed.). Butterworth–Heinemann. ISBN 0-7506-2768-9. 
• Weinberg, S. (2002), The Quantum Theory of Fields, 1, Cambridge University Press, ISBN 0-521-55001-7, https://archive.org/details/quantumtheoryoff00stev