Thomas Rotation
   HOME

TheInfoList



OR:

In
theoretical physics Theoretical physics is a branch of physics that employs mathematical models and abstractions of physical objects and systems to rationalize, explain and predict natural phenomena. This is in contrast to experimental physics, which uses experim ...
, the composition of two non- collinear
Lorentz boost In physics, the Lorentz transformations are a six-parameter family of linear transformations from a coordinate frame in spacetime to another frame that moves at a constant velocity relative to the former. The respective inverse transformation i ...
s results in a Lorentz transformation that is not a pure boost but is the composition of a boost and a rotation. This rotation is called Thomas rotation, Thomas–Wigner rotation or Wigner rotation. The rotation was discovered and proved by
Ludwik Silberstein Ludwik Silberstein (1872 – 1948) was a Polish-American physicist who helped make special relativity and general relativity staples of university coursework. His textbook '' The Theory of Relativity'' was published by Macmillan in 1914 with a sec ...
in his 1914 book 'Relativity', rediscovered by Llewellyn Thomas in 1926, and rederived by Wigner in 1939. Wigner acknowledged Silberstein. If a sequence of non-collinear boosts returns an object to its initial velocity, then the sequence of Wigner rotations can combine to produce a net rotation called the
Thomas precession In physics, the Thomas precession, named after Llewellyn Thomas, is a relativistic correction that applies to the spin of an elementary particle or the rotation of a macroscopic gyroscope and relates the angular velocity of the spin of a pa ...
. There are still ongoing discussions about the correct form of equations for the Thomas rotation in different reference systems with contradicting results. Goldstein: :''The spatial rotation resulting from the successive application of two non-collinear Lorentz transformations have been declared every bit as paradoxical as the more frequently discussed apparent violations of common sense, such as the
twin paradox In physics, the twin paradox is a thought experiment in special relativity involving identical twins, one of whom makes a journey into space in a high-speed rocket and returns home to find that the twin who remained on Earth has aged more. Thi ...
.'' Einstein's principle of velocity reciprocity (EPVR) reads :''We postulate that the relation between the coordinates of the two systems is linear. Then the inverse transformation is also linear and the complete non-preference of the one or the other system demands that the transformation shall be identical with the original one, except for a change of to '' With less careful interpretation, the EPVR is seemingly violated in some models. There is, of course, no true paradox present. Let it be u the velocity in which the lab reference frame moves respect an object called A and let it be v the velocity in which another object called B is moving, measured from the lab reference frame. If u and v are not aligned the relative velocities of these two bodies will not be opposite, that is \mathbf v_ \neq - \mathbf v_ since there is a rotation between them The velocity that A will measure on B will be: :\mathbf v_ = \frac \left left(1 + \frac\frac \mathbf u \cdot \mathbf v \right) \mathbf u + \frac \mathbf v \right The Lorentz factor for the velocities that either A sees on B or B sees on A: :\gamma = \gamma_ = \gamma_ = \gamma_\mathbf \gamma_\mathbf\left(1+ \frac\right) \,, The angle of rotation can be calculated in two ways: :\cos\epsilon = \frac - 1~, Or: :\cos\epsilon = - \frac And the axis of rotation is: :\mathbf = -\frac~.


Setup of frames and relative velocities between them


Two general boosts

When studying the Thomas rotation at the fundamental level, one typically uses a setup with three coordinate frames, . Frame has velocity relative to frame , and frame has velocity relative to frame . The axes are, by construction, oriented as follows. Viewed from , the axes of and are parallel (the same holds true for the pair of frames when viewed from .) Also viewed from , the spatial axes of and are parallel (and the same holds true for the pair of frames when viewed from .) This is an application of EVPR: If is the velocity of relative to , then is the velocity of relative to . The velocity makes the ''same'' angles with respect to coordinate axes in both the primed and unprimed systems. This does ''not'' represent a snapshot taken in any of the two frames of the combined system at any particular time, as should be clear from the detailed description below. This is possible, since a boost in, say, the positive , preserves orthogonality of the coordinate axes. A general boost can be expressed as , where is a rotation taking the into the direction of and is a boost in the new . Each rotation retains the property that the spatial coordinate axes are orthogonal. The boost will stretch the (intermediate) by a factor , while leaving the and in place. The fact that coordinate axes are non-parallel in this construction after ''two'' consecutive non-collinear boosts is a precise expression of the phenomenon of Thomas rotation.This preservation of orthogonality of ''coordinate axes'' should not be confused with preservation of angles between spacelike vectors taken at one and the same time in one system, which, of course, does not hold. The coordinate axes transform under the ''passive'' transformation presented, while the vectors transform under the corresponding ''active'' transformation. The velocity of as seen in is denoted , where ⊕ refers to the relativistic addition of velocity (and not ordinary vector addition), given by and : \gamma_\mathbf = \frac is the Lorentz factor of the velocity (the vertical bars indicate the magnitude of the vector). The velocity can be thought of the velocity of a frame relative to a frame , and is the velocity of an object, say a particle or ''another frame'' relative to . In the present context, all velocities are best thought of as relative velocities of frames unless otherwise specified. The result is then the relative velocity of frame relative to a frame . Although velocity addition is
nonlinear In mathematics and science, a nonlinear system is a system in which the change of the output is not proportional to the change of the input. Nonlinear problems are of interest to engineers, biologists, physicists, mathematicians, and many other ...
, non-
associative In mathematics, the associative property is a property of some binary operations, which means that rearranging the parentheses in an expression will not change the result. In propositional logic, associativity is a valid rule of replacement f ...
, and non- commutative, the result of the operation correctly obtains a velocity with a magnitude less than . If ordinary vector addition was used, it would be possible to obtain a velocity with a magnitude larger than . The Lorentz factor of both composite velocities are equal, :\gamma = \gamma_ = \gamma_ = \gamma_\mathbf \gamma_\mathbf\left(1+ \frac\right) \,, and the
norm Naturally occurring radioactive materials (NORM) and technologically enhanced naturally occurring radioactive materials (TENORM) consist of materials, usually industrial wastes or by-products enriched with radioactive elements found in the envir ...
s are equal under interchange of velocity vectors :, \mathbf u \oplus \mathbf v, = , \mathbf v \oplus \mathbf u, = \frac\sqrt \,. Since the two possible composite velocities have equal magnitude, but different directions, one must be a rotated copy of the other. More detail and other properties of no direct concern here can be found in the main article.


Reversed configuration

Consider the reversed configuration, namely, frame moves with velocity relative to frame , and frame , in turn, moves with velocity relative to frame . In short, and by EPVR. Then the velocity of relative to is . By EPVR again, the velocity of relative to is then . One finds . While they are equal in magnitude, there is an angle between them. For a single boost between two inertial frames, there is only one unambiguous relative velocity (or its negative). For two boosts, the peculiar result of ''two'' inequivalent relative velocities instead of one seems to contradict the symmetry of relative motion between any two frames. Which is the correct velocity of relative to ? Since this inequality may be somewhat unexpected and potentially breaking EPVR, this question is warranted.This is sometimes called the "Mocanu paradox". Mocanu himself didn't name it a paradox, but rather a "difficulty" within the framework of relativistic electrodynamics in a 1986 paper. He was also quick to acknowledge that the problem is explained by Thomas precession , but the name lingers on.


Formulation in terms of Lorentz transformations


Two boosts equals a boost and rotation

The answer to the question lies in the Thomas rotation, and that one must be careful in specifying which coordinate system is involved at each step. When viewed from , the coordinate axes of and are ''not'' parallel. While this can be hard to imagine since both pairs and have parallel coordinate axes, it is easy to explain mathematically. Velocity addition does not provide a complete description of the relation between the frames. One must formulate the complete description in terms of Lorentz transformations corresponding to the velocities. A Lorentz boost with any velocity (magnitude less than ) is given symbolically by :X' = B(\mathbf)X where the coordinates and transformation matrix are compactly expressed in
block matrix In mathematics, a block matrix or a partitioned matrix is a matrix that is '' interpreted'' as having been broken into sections called blocks or submatrices. Intuitively, a matrix interpreted as a block matrix can be visualized as the original mat ...
form :X' = \begin c t' \\ \mathbf' \end \quad B(\mathbf) = \begin \gamma_\mathbf & - \dfrac\mathbf^\mathrm \\ -\dfrac\mathbf & \mathbf + \dfrac \dfrac \\ \end \quad X = \begin c t \\ \mathbf \end and, in turn, are column vectors (the matrix transpose of these are row vectors), and is the Lorentz factor of velocity . The boost matrix is a symmetric matrix. The inverse transformation is given by :B(\mathbf )^=B(-\mathbf )\quad \Rightarrow \quad X=B(-\mathbf )X' . It is clear that to each admissible velocity there corresponds a ''pure'' Lorentz boost, :\mathbf \leftrightarrow B(\mathbf). Velocity addition corresponds to the composition of boosts in that order. The acts on first, then acts on . Notice succeeding operators act on the ''left'' in any composition of operators, so should be interpreted as a boost with velocities then , not then . Performing the Lorentz transformations by block matrix multiplication, :X'' = B(\mathbf) X' \,,\quad X' = B(\mathbf) X \quad \Rightarrow \quad X'' = \Lambda X the composite transformation matrix is : \Lambda = B(\mathbf) B(\mathbf) = \begin\gamma & -\mathbf^\mathrm\\ -\mathbf & \mathbf \end and, in turn, :\begin \gamma &= \gamma_\mathbf\gamma_\mathbf\left(1+\frac\right) \\ \mathbf &= \frac \mathbf\oplus\mathbf \,, \quad \mathbf = \frac \mathbf\oplus\mathbf \\ \mathbf &= \gamma_\mathbf\gamma_\mathbf\frac + \left(\mathbf + \frac\frac\right)\left(\mathbf + \frac\frac\right) \end Here is the composite Lorentz factor, and and are 3×1 column vectors proportional to the composite velocities. The 3×3 matrix will turn out to have geometric significance. The inverse transformations are : X = B(-\mathbf) X' \,,\quad X' = B(-\mathbf) X'' \quad \Rightarrow \quad X = \Lambda^ X'' and the composition amounts to a negation ''and'' exchange of velocities, :\Lambda^ = B(-\mathbf) B(-\mathbf) = \begin\gamma & \mathbf^\mathrm\\ \mathbf & \mathbf^\mathrm \end If the relative velocities are exchanged, looking at the blocks of , one observes the composite transformation to be the matrix transpose of . This is not the same as the original matrix, so the composite Lorentz transformation matrix is not symmetric, and thus not a single boost. This, in turn, translates to the incompleteness of velocity composition from the result of two boosts; symbolically, :B(\mathbf u \oplus \mathbf v) \ne B(\mathbf v)B(\mathbf u) \,. To make the description complete, it is necessary to introduce a rotation, before or after the boost. This rotation is the Thomas rotation. A rotation is given by :X' = R(\boldsymbol)X where the 4×4 rotation matrix is :R() = \begin 1 & 0 \\ 0 & \mathbf(\boldsymbol) \end and is a 3×3
rotation matrix In linear algebra, a rotation matrix is a transformation matrix that is used to perform a rotation in Euclidean space. For example, using the convention below, the matrix :R = \begin \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end ...
. In the literature, the 3d rotation matrix may be denoted by other letters, others use a name and the relative velocity vectors involved; e.g., for "Thomas rotation" or for "gyration" (see
gyrovector space A gyrovector space is a mathematical concept proposed by Abraham A. Ungar for studying hyperbolic geometry in analogy to the way vector spaces are used in Euclidean geometry.Abraham A. Ungar (2005), "Analytic Hyperbolic Geometry: Mathematical Fo ...
). Correspondingly the 4d rotation matrix (non-bold italic) in this article may be denoted : \mathrm = \begin1 & 0 \\ 0 & \mathrm mathbf, \mathbf\end \quad \text \quad \mathrm = \begin1 & 0 \\ 0 & \mathrm mathbf, \mathbf\end
In this article the axis-angle representation is used, and is the "axis-angle vector", the angle multiplied by a unit vector parallel to the axis. Also, the ''right-handed'' convention for the spatial coordinates is used (see orientation (vector space)), so that rotations are positive in the anticlockwise sense according to the
right-hand rule In mathematics and physics, the right-hand rule is a common mnemonic for understanding orientation of axes in three-dimensional space. It is also a convenient method for quickly finding the direction of a cross-product of 2 vectors. Most of th ...
, and negative in the clockwise sense. With these conventions; the rotation matrix rotates any 3d vector about the axis through angle anticlockwise (an
active transformation Active may refer to: Music * ''Active'' (album), a 1992 album by Casiopea * Active Records, a record label Ships * ''Active'' (ship), several commercial ships by that name * HMS ''Active'', the name of various ships of the British Royal ...
), which has the equivalent effect of rotating the coordinate frame clockwise about the same axis through the same angle (a passive transformation). The rotation matrix is an orthogonal matrix, its transpose equals its inverse, and negating either the angle or axis in the rotation matrix corresponds to a rotation in the opposite sense, so the inverse transformation is readily obtained by : R(\boldsymbol)^ = R(\boldsymbol)^\mathrm = R(-\boldsymbol) \quad \Rightarrow \quad X = R(-\boldsymbol)X'\,. A boost followed or preceded by a rotation is also a Lorentz transformation, since these operations leave the spacetime interval invariant. The same Lorentz transformation has two decompositions for appropriately chosen rapidity and axis-angle vectors; :\Lambda(\boldsymbol,\mathbf) = R(\boldsymbol)B(\mathbf) :\Lambda(\mathbf,\boldsymbol) = B(\mathbf)R(\boldsymbol) and if these are two decompositions are equal, the two boosts are related by :B(\mathbf) = R(-\boldsymbol)B(\mathbf)R(\boldsymbol) so the boosts are related by a matrix similarity transformation. It turns out the equality between two boosts and a rotation followed or preceded by a single boost is correct: the rotation of frames matches the angular separation of the composite velocities, and explains how one composite velocity applies to one frame, while the other applies to the rotated frame. The rotation also breaks the symmetry in the overall Lorentz transformation making it nonsymmetric. For this specific rotation, let the angle be and the axis be defined by the unit vector , so the axis-angle vector is . Altogether, two different orderings of two boosts means there are two inequivalent transformations. Each of these can be split into a boost then rotation, or a rotation then boost, doubling the number of inequivalent transformations to four. The inverse transformations are equally important; they provide information about what the other observer perceives. In all, there are eight transformations to consider, just for the problem of two Lorentz boosts. In summary, with subsequent operations acting on the left, they are Matching up the boosts followed by rotations, in the original setup, an observer in notices to move with velocity then rotate clockwise (first diagram), and because of the rotation an observer in Σ′′ notices to move with velocity then rotate anticlockwise (second diagram). If the velocities are exchanged an observer in notices to move with velocity then rotate anticlockwise (third diagram), and because of the rotation an observer in notices to move with velocity then rotate clockwise (fourth diagram). The cases of rotations then boosts are similar (no diagrams are shown). Matching up the rotations followed by boosts, in the original setup, an observer in notices to rotate clockwise then move with velocity , and because of the rotation an observer in notices to rotate anticlockwise then move with velocity . If the velocities are exchanged an observer in notices to rotate anticlockwise then move with velocity , and because of the rotation an observer in notices to rotate clockwise then move with velocity .


Finding the axis and angle of the Thomas rotation

The above formulae constitute the relativistic velocity addition and the Thomas rotation explicitly in the general Lorentz transformations. Throughout, in every composition of boosts and decomposition into a boost and rotation, the important formula :\mathbf = \mathbf + \frac1\mathbf^\mathrm holds, allowing the rotation matrix to be defined completely in terms of the relative velocities and . The angle of a rotation matrix in the axis–angle representation can be found from the trace of the rotation matrix, the general result for ''any'' axis is . Taking the trace of the equation gives :\cos\epsilon = \frac - 1 The angle between and is ''not'' the same as the angle between and . In both frames Σ and Σ′′, for every composition and decomposition, another important formula :\mathbf = \mathbf holds. The vectors and are indeed related by a rotation, in fact by the same rotation matrix which rotates the coordinate frames. Starting from , the matrix rotates this into anticlockwise, it follows their
cross product In mathematics, the cross product or vector product (occasionally directed area product, to emphasize its geometric significance) is a binary operation on two vectors in a three-dimensional oriented Euclidean vector space (named here E), and is ...
(in the right-hand convention) :\mathbf\times\mathbf = \frac \mathbf\times\mathbf defines the axis correctly, therefore the axis is also parallel to . The magnitude of this pseudovector is neither interesting nor important, only the direction is, so it can be normalized into the
unit vector In mathematics, a unit vector in a normed vector space is a vector (often a spatial vector) of length 1. A unit vector is often denoted by a lowercase letter with a circumflex, or "hat", as in \hat (pronounced "v-hat"). The term ''direction vecto ...
:\mathbf = \frac which still completely defines the direction of the axis without loss of information. The rotation is simply a "static" rotation and there is no relative rotational motion between the frames, there is relative translational motion in the boost. However, if the frames accelerate, then the rotated frame rotates with an angular velocity. This effect is known as the
Thomas precession In physics, the Thomas precession, named after Llewellyn Thomas, is a relativistic correction that applies to the spin of an elementary particle or the rotation of a macroscopic gyroscope and relates the angular velocity of the spin of a pa ...
, and arises purely from the kinematics of successive Lorentz boosts.


Finding the Thomas rotation

In principle, it is pretty easy. Since every Lorentz transformation is a product of a boost and a rotation, the consecutive application of two pure boosts is a pure boost, either followed by or preceded by a pure rotation. Thus, suppose :\Lambda = B(\mathbf w)R. The task is to glean from this equation the boost velocity and the rotation from the matrix entries of . The coordinates of events are related by :x'^\mu = _\nu x^\nu. Inverting this relation yields :_\mu_\rho x^\rho = _\mu x'^\mu, or :x^\nu = ^\nu x'^\mu. Set Then will record the spacetime position of the origin of the primed system, :x^\nu = ^\nu x'^0, or : x = \beginct \\ x_1 \\ x_2 \\ x_3\end = \begin^0 ct' \\ ^1 ct' \\ ^2 ct' \\ ^3 ct'\end. . But :\Lambda^ = (B(\mathbf w)R)^ = R^B(-\mathbf). Multiplying this matrix with a pure rotation will not affect the zeroth columns and rows, and : x = \beginct \\ x_1 \\ x_2 \\ x_3\end = \begin\gamma ct' \\ \gamma\beta_x ct' \\ \gamma\beta_y ct' \\ \gamma\beta_z ct'\end = \begin\gamma ct' \\ \gamma w_x t' \\ \gamma w_y t' \\ \gamma w_z t'\end = \gamma \begin ct' \\ w_x t' \\ w_y t' \\ w_z t'\end, which could have been anticipated from the formula for a simple boost in the -direction, and for the relative velocity vector : \frac\mathbf x = \frac = \boldsymbol \beta = \begin \frac \\ \frac \\ \frac \end = \begin \beta_x \\ \beta_y \\ \beta_z \end = \begin^1/^0 \\ ^2/^0 \\ ^3/^0 \end. Thus given with , one obtains and by little more than inspection of . (Of course, can also be found using velocity addition per above.) From , construct . The solution for is then :R = B(-\mathbf w)\Lambda. With the ansatz :\Lambda = RB(\mathbf w), one finds by the same means :R = \Lambda B(-\mathbf w). Finding a formal solution in terms of velocity parameters and involves first ''formally'' multiplying , formally inverting, then reading off form the result, ''formally'' building from the result, and, finally, formally multiplying . It should be clear that this is a daunting task, and it is difficult to interpret/identify the result as a rotation, though it is clear a priori that it is. It is these difficulties that the Goldstein quote at the top refers to. The problem has been thoroughly studied under simplifying assumptions over the years.


Group theoretical origin

Another way to explain the origin of the rotation is by looking at the generators of the Lorentz group.


Boosts from velocities

The passage from a velocity to a boost is obtained as follows. An arbitrary boost is given by : e^, where is a triple of real numbers serving as coordinates on the boost subspace of the Lie algebra spanned by the matrices :(K_1, K_2, K_3) = \left( \left \begin0&1&0&0\\1&0&0&0\\0&0&0&0\\0&0&0&0 \end \right \left \begin0&0&1&0\\0&0&0&0\\1&0&0&0\\0&0&0&0 \end \right \left \begin0&0&0&1\\0&0&0&0\\0&0&0&0\\1&0&0&0 \end \right\right). The vector :\boldsymbol \zeta = \frac \tanh^\beta is called the ''boost parameter'' or ''boost vector'', while its norm is the
rapidity In relativity, rapidity is commonly used as a measure for relativistic velocity. Mathematically, rapidity can be defined as the hyperbolic angle that differentiates two frames of reference in relative motion, each frame being associated with di ...
. Here is the ''velocity parameter'', the magnitude of the vector . While for one has , the parameter is confined within , and hence . Thus : e^ = e^ \equiv B(\mathbf u)~. The set of velocities satisfying is an open ball in and is called the space of admissible velocities in the literature. It is endowed with a hyperbolic geometry described in the linked article.


Commutators

The generators of boosts, , in different directions do not commute. This has the effect that two consecutive boosts is not a pure boost in general, but a rotation preceding a boost. Consider a succession of boosts in the x direction, then the y direction, expanding each boost to first order :e^e^ = (I - \zeta_y K_y + \cdots)(I - \zeta_x K_x + \cdots) = I - \zeta_x K_x - \zeta_y K_y + \zeta_x\zeta_y K_y K_x + \cdots then :e^e^ = I + \zeta_x K_x + \zeta_y K_y + \zeta_x\zeta_y K_y K_x + \cdots and the group commutator is :\begin e^e^e^e^ = I &+ \zeta_x \zeta_y _y, K_x- (\zeta_x K_x)^2 - (\zeta_y K_y)^2 \\ &+ \zeta_x^2 \zeta_y _x, K_yK_x + \zeta_x\zeta_y^2 K_y _y, K_x\ &+ (\zeta_x\zeta_y)^2 K_y K_x K_y K_x + \cdots \end Three of the
commutation relation In mathematics, the commutator gives an indication of the extent to which a certain binary operation fails to be commutative. There are different definitions used in group theory and ring theory. Group theory The commutator of two elements, a ...
s of the Lorentz generators are :\begin[] [J_x, J_y] &= J_z \\[] _x, K_y&= -J_z \\[] [J_x, K_y] &= K_z \end where the bracket is a binary operation known as the ''
commutator In mathematics, the commutator gives an indication of the extent to which a certain binary operation fails to be commutative. There are different definitions used in group theory and ring theory. Group theory The commutator of two elements, a ...
'', and the other relations can be found by taking cyclic permutations of x, y, z components (i.e. change x to y, y to z, and z to x, repeat). Returning to the group commutator, the commutation relations of the boost generators imply for a boost along the x then y directions, there will be a rotation about the z axis. In terms of the rapidities, the rotation angle is given by : \tan \frac=\tanh \frac\tanh \frac , equivalently expressible as : \tan \theta=\frac ~.


Spacetime diagrams for non-collinear boosts

The familiar notion of vector addition for velocities in the
Euclidean plane In mathematics, the Euclidean plane is a Euclidean space of dimension two. That is, a geometric setting in which two real quantities are required to determine the position of each point ( element of the plane), which includes affine notions of ...
can be done in a triangular formation, or since vector addition is commutative, the vectors in both orderings geometrically form a parallelogram (see " parallelogram law"). This does not hold for relativistic velocity addition; instead a hyperbolic triangle arises whose edges are related to the rapidities of the boosts. Changing the order of the boost velocities, one does not find the resultant boost velocities to coincide.


See also

* Bargmann-Michel-Telegdi equation * Pauli–Lubanski pseudovector * Velocity-addition formula#Hyperbolic geometry


Footnotes


References

* *Sexl Urbantke mention on p. 39 Lobachevsky geometry needs to be introduced into the usual Minkowski
spacetime diagram A spacetime diagram is a graphical illustration of the properties of space and time in the special theory of relativity. Spacetime diagrams allow a qualitative understanding of the corresponding phenomena like time dilation and length contractio ...
s for non-collinear velocities. *. * * * * Ferraro, R., & Thibeault, M. (1999). "Generic composition of boosts: an elementary derivation of the Wigner rotation". ''European journal of physics'' 20(3):143. * * (free access) * * * * * * * * * * * * * * Thomas L.H The kinematics of an electron with an axis,Phil. Mag. 7, 1927 http://www.clifford.org/drbill/csueb/4250/topics/thomas_papers/Thomas1927.pdf * Silberstein L. The Theory of Relativity, MacMillan 1914


Further reading

*''Relativistic velocity space, Wigner rotation, and Thomas precession'' (2004) John A. Rhodes and Mark D. Semon *''The Hyperbolic Theory of Special Relativity'' (2006) by J.F. Barrett {{relativity Special relativity Coordinate systems Theory of relativity Theoretical physics Mathematical physics