Relativistic Lagrangian mechanics explained
In theoretical physics, relativistic Lagrangian mechanics is Lagrangian mechanics applied in the context of special relativity and general relativity.
Introduction
The relativistic Lagrangian can be derived in relativistic mechanics to be of the form:
Although, unlike non-relativistic mechanics, the relativistic Lagrangian is not expressed as difference of kinetic energy with potential energy, the relativistic Hamiltonian corresponds to total energy in a similar manner but without including rest energy. The form of the Lagrangian also makes the relativistic action functional proportional to the proper time of the path in spacetime.
In covariant form, the Lagrangian is taken to be:
where
σ is an
affine parameter which parametrizes the spacetime curve.
Lagrangian formulation in special relativity
Lagrangian mechanics can be formulated in special relativity as follows. Consider one particle (N particles are considered later).
Coordinate formulation
If a system is described by a Lagrangian L, the Euler–Lagrange equations
retain their form in
special relativity, provided the Lagrangian generates equations of motion consistent with special relativity. Here /is the
position vector of the particle as measured in some lab frame where
Cartesian coordinates are used for simplicity, and
is the coordinate velocity, the
derivative of position
r with respect to
coordinate time t. (Throughout this article, overdots are with respect to coordinate time, not proper time). It is possible to transform the position coordinates to
generalized coordinates exactly as in non-relativistic mechanics, . Taking the total differential of
r obtains the transformation of velocity
v to the generalized coordinates, generalized velocities, and coordinate time
remains the same. However, the
energy of a moving particle is different from non-relativistic mechanics. It is instructive to look at the total
relativistic energy of a free test particle. An observer in the lab frame defines events by coordinates
r and
coordinate time t, and measures the particle to have coordinate velocity . By contrast, an observer moving with the particle will record a different time, this is the
proper time,
τ. Expanding in a
power series, the first term is the particle's
rest energy, plus its non-relativistic
kinetic energy, followed by higher order relativistic corrections;
} = m_0 c^2 + m_0 \dot^2 (t) + m_0 \frac + \cdots \,,where
c is the
speed of light in vacuum. The
differentials in
t and
τ are related by the
Lorentz factor γ,
[1]
} \,,\quad \dot = \frac \,, \quad \dot^2 (t) = \dot(t) \cdot \dot(t)\,.where · is the
dot product. The relativistic kinetic energy for an uncharged particle of
rest mass m0 is
and we may naïvely guess the relativistic Lagrangian for a particle to be this relativistic kinetic energy minus the potential energy. However, even for a free particle for which
V = 0, this is wrong. Following the non-relativistic approach, we expect the derivative of this seemingly correct Lagrangian with respect to the velocity to be the relativistic momentum, which it is not.
The definition of a generalized momentum can be retained, and the advantageous connection between cyclic coordinates and conserved quantities will continue to apply. The momenta can be used to "reverse-engineer" the Lagrangian. For the case of the free massive particle, in Cartesian coordinates, the x component of relativistic momentum is
and similarly for the
y and
z components. Integrating this equation with respect to
dx/
dt gives
where
X is an arbitrary function of
dy/
dt and
dz/
dt from the integration. Integrating
py and
pz obtains similarly
where
Y and
Z are arbitrary functions of their indicated variables. Since the functions
X,
Y,
Z are arbitrary, without loss of generality we can conclude the common solution to these integrals, a possible Lagrangian that will correctly generate all the components of relativistic momentum, is
where .
Alternatively, since we wish to build a Lagrangian out of relativistically invariant quantities, take the action as proportional to the integral of the Lorentz invariant line element in spacetime, the length of the particle's world line between proper times τ1 and τ2,[1]
S=\varepsilon
d\tau=\varepsilon
, L=
=\varepsilon\sqrt{1-
}\,,where
ε is a constant to be found, and after converting the proper time of the particle to the coordinate time as measured in the lab frame, the integrand is the Lagrangian by definition. The momentum must be the relativistic momentum,
p=
=\left(
| -\varepsilon | \right)\gamma( |
c2 |
=m0\gamma(
,
which requires
ε = −
m0c2, in agreement with the previously obtained Lagrangian.
Either way, the position vector r is absent from the Lagrangian and therefore cyclic, so the Euler–Lagrange equations are consistent with the constancy of relativistic momentum,
which must be the case for a free particle. Also, expanding the relativistic free particle Lagrangian in a power series to first order in,
L=-m0c2\left[1+
\left(-
\right)+ … \right] ≈ -m0c2+
2,
in the non-relativistic limit when
v is small, the higher order terms not shown are negligible, and the Lagrangian is the non-relativistic kinetic energy as it should be. The remaining term is the negative of the particle's rest energy, a constant term which can be ignored in the Lagrangian.
For the case of an interacting particle subject to a potential V, which may be non-conservative, it is possible for a number of interesting cases to simply subtract this potential from the free particle Lagrangian,
and the Euler–Lagrange equations lead to the relativistic version of Newton's second law. The derivative of relativistic momentum with respect to the time coordinate is equal to the force acting on the particle:
assuming the potential
V can generate the corresponding force
F in this way. If the potential cannot obtain the force as shown, then the Lagrangian would need modification to obtain the correct equations of motion.
Although this has been shown by taking Cartesian coordinates, it follows due to invariance of Euler Lagrange equations, that it is also satisfied in any arbitrary co-ordinate system as it physically corresponds to action minimization being independent of the co-ordinate system used to describe it. In a similar manner, several properties in Lagrangian mechanics are preserved whenever they are also independent of the specific form of the Lagrangian or the laws of motion governing the particles. For example, it is also true that if the Lagrangian is explicitly independent of time and the potential V(r) independent of velocities, then the total relativistic energy
is conserved, although the identification is less obvious since the first term is the relativistic energy of the particle which includes the rest mass of the particle, not merely the relativistic kinetic energy. Also, the argument for homogeneous functions does not apply to relativistic Lagrangians.
The extension to N particles is straightforward, the relativistic Lagrangian is just a sum of the "free particle" terms, minus the potential energy of their interaction;
L=-c2
-V(r1,r2,\ldots,
2,\ldots,t),
where all the positions and velocities are measured in the same lab frame, including the time.
The advantage of this coordinate formulation is that it can be applied to a variety of systems, including multiparticle systems. The disadvantage is that some lab frame has been singled out as a preferred frame, and none of the equations are manifestly covariant (in other words, they do not take the same form in all frames of reference). For an observer moving relative to the lab frame, everything must be recalculated; the position r, the momentum p, total energy E, potential energy, etc. In particular, if this other observer moves with constant relative velocity then Lorentz transformations must be used. However, the action will remain the same since it is Lorentz invariant by construction.
A seemingly different but completely equivalent form of the Lagrangian for a free massive particle, which will readily extend to general relativity as shown below, can be obtained by inserting[1]
d\tau=
\sqrt{η\alpha\beta
} dt \,, into the Lorentz invariant action so that
S=\varepsilon
\sqrt{η\alpha\beta
} dt \quad\Rightarrow\quad L = \frac\sqrt where is retained for simplicity. Although the line element and action are Lorentz invariant, the Lagrangian is
not, because it has explicit dependence on the lab coordinate time. Still, the equations of motion follow from
Hamilton's principle
Since the action is proportional to the length of the particle's worldline (in other words its trajectory in spacetime), this route illustrates that finding the stationary action is asking to find the trajectory of shortest or largest length in spacetime. Correspondingly, the equations of motion of the particle are akin to the equations describing the trajectories of shortest or largest length in spacetime, geodesics.
For the case of an interacting particle in a potential V, the Lagrangian is still
} - V,which can also extend to many particles as shown above, each particle has its own set of position coordinates to define its position.
Covariant formulation
In the covariant formulation, time is placed on equal footing with space, so the coordinate time as measured in some frame is part of the configuration space alongside the spatial coordinates (and other generalized coordinates). For a particle, either massless or massive, the Lorentz invariant action is (abusing notation)
S=
Λ(x\nu(\sigma),u\nu(\sigma),\sigma)d\sigma,
where lower and upper indices are used according to
covariance and contravariance of vectors,
σ is an
affine parameter, and is the
four-velocity of the particle.
For massive particles, σ can be the arc length s, or proper time τ, along the particle's world line,
ds2=c2d\tau2=g\alpha\betadx\alphadx\beta.
For massless particles, it cannot because the proper time of a massless particle is always zero;
g\alpha\betadx\alphadx\beta=0.
For a free particle, the Lagrangian has the form
where the irrelevant factor of 1/2 is allowed to be scaled away by the scaling property of Lagrangians. No inclusion of mass is necessary since this also applies to massless particles. The Euler–Lagrange equations in the spacetime coordinates are
which is the geodesic equation for affinely parameterized geodesics in spacetime. In other words, the free particle follows geodesics. Geodesics for massless particles are called "null geodesics", since they lie in a "
light cone" or "null cone" of spacetime (the null comes about because their inner product via the metric is equal to 0), massive particles follow "timelike geodesics", and hypothetical particles that travel faster than light known as
tachyons follow "spacelike geodesics".
This manifestly covariant formulation does not extend to an N-particle system, since then the affine parameter of any one particle cannot be defined as a common parameter for all the other particles.
Examples in special relativity
Special relativistic 1d free particle
For a 1d relativistic free particle, the Lagrangian is
}\,.
This results in the following equation of motion:
Special relativistic 1d harmonic oscillator
For a 1d relativistic simple harmonic oscillator, the Lagrangian is
} - \fracx^2 \,.where
k is the spring constant.
Special relativistic constant force
For a particle under a constant force, the Lagrangian is
} - mgx \,,where
g is the force per unit mass.
This results in the following equation of motion:
Which, given initial conditions of
\begin{align}
x(t=0)&=x0\\
(t=0)&=v0
\end{align}
results in the position of the particle as a function of time being
} - \sqrt\right] \,.
Special relativistic test particle in an electromagnetic field
See main article: Covariant formulation of classical electromagnetism.
In special relativity, the Lagrangian of a massive charged test particle in an electromagnetic field modifies to
} - q \phi + q \dot \cdot \mathbf \,.
The Lagrangian equations in r lead to the Lorentz force law, in terms of the relativistic momentum
}\right) = q \mathbf + q \dot \times \mathbf \,.
In the language of four-vectors and tensor index notation, the Lagrangian takes the form
where
uμ =
dxμ/
dτ is the
four-velocity of the test particle, and
Aμ the
electromagnetic four-potential.
The Euler–Lagrange equations are (notice the total derivative with respect to proper time instead of coordinate time)
obtains
| \mu | \partialA\mu | \partialx\nu |
|
qu | |
=
(mu\nu+qA\nu).
Under the total derivative with respect to proper time, the first term is the relativistic momentum, the second term is
then rearranging, and using the definition of the antisymmetric
electromagnetic tensor, gives the covariant form of the Lorentz force law in the more familiar form,
(mu\nu)=qu\muF\nu\mu, F\nu\mu=
-
.
Lagrangian formulation in general relativity
The Lagrangian is that of a single particle plus an interaction term LI
Varying this with respect to the position of the particle xα as a function of time t gives
\begin{align}\deltaL&=m
\delta\left(g\mu\nu
\right)+\deltaLI\\
&=m
\left(g\mu\nu,\alpha\deltax\alpha
+2g\alpha\nu
\right)+
| \partialLI |
\partialx\alpha |
\deltax\alpha+
\\
&=
g\mu\nu,\alpha\deltax\alpha
-
\left(mg\alpha\nu
\right)\deltax\alpha+
| \partialLI |
\partialx\alpha |
\deltax\alpha-
\left(
\right)\deltax\alpha+
.
\end{align}
This gives the equation of motion
0=
g\mu\nu,\alpha
-
\left(mg\alpha\nu
\right)+f\alpha
where
f\alpha=
| \partialLI |
\partialx\alpha |
-
\left(
\right)
is the non-gravitational force on the particle. (For
m to be independent of time, we must have .)
Rearranging gets the force equation
\left(m
\right)=-m
+g\nu\alphaf\alpha
where Γ is the
Christoffel symbol, which describes the gravitational field.
If we let
be the (kinetic) linear momentum for a particle with mass, then
=-
p\mu
+g\nu\alphaf\alpha
and
hold even for a massless particle.
Examples in general relativity
General relativistic test particle in an electromagnetic field
In general relativity, the first term generalizes (includes) both the classical kinetic energy and the interaction with the gravitational field. For a charged particle in an electromagnetic field, the Lagrangian is given by
)=-mc2\sqrt{-c-2g\mu\nu(x(\tau))
} + q \frac A_(x(\tau))\,.
If the four spacetime coordinates xμ are given in arbitrary units (i.e. unitless), then gμν is the rank 2 symmetric metric tensor, which is also the gravitational potential. Also, Aμ is the electromagnetic 4-vector potential.
There exists an equivalent formulation of the relativistic Lagrangian, which has two advantages:
- it allows for a generalization to massless particles and tachyons;
- it is based on an energy functional instead of a length functional, such that it does not contain a square root.
In this alternative formulation, the Lagrangian is given by
,e)=
g\mu\nu(x(λ))
-
+qA\mu(x(λ))
,where
is an arbitrary affine parameter and
is an auxiliary parameter that can be viewed as an
einbein field along the worldline. In the original Lagrangian with the square root the
energy-momentum relation appears as a
primary constraint that is also a
first class constraint. In this reformulation this is no longer the case. Instead, the
energy-momentum relation appears as the equation of motion for the auxiliary field
. Therefore, the constraint is now a
secondary constraint that is still a
first class constraint, reflecting the invariance of the action under reparameterization of the affine parameter
. After the equation of motion has been derived, one must gauge fix the auxiliary field
. The standard gauge choice is as follows:
, one fixes
. This choice automatically fixes
, i.e. the affine parameter is fixed to be the proper time.
, one fixes
. This choice automatically fixes
, i.e. the affine parameter is fixed to be the proper length.
, there is no choice that fixes the affine parameter
to a physical parameter. Consequently, there is some freedom in fixing the auxiliary field. The two common choices are:
. In this case,
does not carry a dependence on the affine parameter
, but the affine parameter is measured in units of time per unit of mass, i.e.
.
, where
is the energy of the particle. In this case, the affine parameter is measured in units of time, i.e.
, but
retains a dependence on the affine parameter
.
See also
References
- Book: Penrose, Roger. The Road to Reality. Vintage books. 2007. 978-0-679-77631-4.
- Book: Mechanics. limited. 3rd. 134. Landau. L. D.. Lifshitz. E. M.. 15 January 1976. L. D. Landau. E. M. Lifshitz. Butterworth Heinemann. 978-0-7506-2896-9.
- Book: Lev Davidovich Landau. Evgeny Mikhailovich Lifshitz. Landau. Lifshitz. Lev . Evgeny . The Classical Theory of Fields. 1975. Elsevier Ltd.. 978-0-7506-2768-9.
- Book: Analytical Mechanics. limited. 2nd. 23. Hand. L. N.. Finch. J. D.. 13 November 1998. Cambridge University Press. 978-0-521-57572-0.
- Book: Analytical mechanics . Louis N. Hand . Janet D. Finch . 140–141 . 0-521-57572-9 . 1998 . Cambridge University Press.
- Book: Goldstein . Herbert . Herbert Goldstein . Classical Mechanics . limited . 2nd . 1980. 0-201-02918-9. Addison Wesley . San Francisco, CA . 352–353.
- Book: Goldstein . Herbert . Herbert Goldstein . Poole . Charles P. Jr. . Safko . John L. . Classical Mechanics . 3rd . 2002 . 0-201-65702-3 . Addison Wesley . San Francisco, CA . 347–349.
- Book: Lanczos, Cornelius . The variational principles of mechanics . 43 . II §5 Auxiliary conditions: the Lagrangian λ-method . 0-486-65067-7 . Courier Dover . 1986 . Reprint of University of Toronto 1970 4th . https://books.google.com/books?id=ZWoYYr8wk2IC&pg=PA43 .
- Book: Feynman. R. P.. Richard Feynman. Leighton. R. B.. Robert B. Leighton. Sands. M. . Matthew Sands. The Feynman Lectures on Physics. 2. Addison Wesley. 0-201-02117-X. 1964.
- Book: Foster. J. Nightingale. J.D.. 2nd. A Short Course in General Relativity. 1995. Springer. 0-03-063366-4.
- Book: M. P. . Hobson . G. P. . Efstathiou . A. N. . Lasenby . General Relativity: An Introduction for Physicists . 79–80 . 978-0-521-82951-9 . 2006 . Cambridge University Press .
Notes and References
- The line element squared is the Lorentz invariant
c2d\tau2=η\alpha\betadx\alphadx\beta=c2dt2-dr2,
which takes the same values in all inertial frames of reference. Here ηαβ are the components of the Minkowski metric tensor, dxα = (cdt, dr) = (cdt, dx, dy, dz) are the components of the differential position four-vector, the summation convention over the covariant and contravariant spacetime indices α and β is used, each index takes the value 0 for timelike components, and 1, 2, 3 for spacelike components, anddr2\equivdr ⋅ dr\equivdx2+dy2+dz2
is a shorthand for the square differential of the particle's position coordinates. Dividing by c2dt2 allows the conversion to the lab coordinate time as follows,
so thatd\tau=
\sqrt{η\alpha\beta
} dt = \frac \,.