Unit 8: Physical Applications: Continuum Mechanics, Electromagnetism & Gravitation
Application of tensor analysis to the fundamental laws of classical physics and general relativity: Cauchy stress and the stress-energy-momentum tensor T^{\mu\nu}, relativistic conservation laws \nabla_\mu T^{\mu\nu} = 0, covariant formulation of Maxwell electrodynamics with electromagnetic field tensor F_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\mu and Lorentz force, the geodesic principle and gravitational time dilation / redshift, derivation and physical interpretation of the Einstein Field Equations G_{\mu\nu} + \Lambda g_{\mu\nu} = (8\pi G / c^4) T_{\mu\nu}, the static spherically symmetric Schwarzschild solution, Mercury's anomalous perihelion precession, and gravitational deflection of light.
§8.1 Stress-Energy-Momentum Tensors & Conservation Laws in Curved Coordinates
1. Cauchy Stress Tensor in Continuum Mechanics
In classical 3D continuum mechanics, internal contact forces distributed across an infinitesimal area element $d\mathbf{S} = \mathbf{n} \, dS$ are characterized by the Cauchy stress tensor $\sigma^{ij}$:
- Conservation of Linear Momentum: In the absence of body forces:
where $\nabla_j$ is the covariant derivative associated with the spatial Riemannian metric $g_{ij}$.
- Conservation of Angular Momentum: By requiring the net torque on an infinitesimal volume element to vanish, the stress tensor must be symmetric:
2. The Relativistic Stress-Energy-Momentum Tensor $T^{\mu\nu}$
In 4-dimensional pseudo-Riemannian spacetime $(M, g_{\mu\nu})$, energy, momentum density, shear stress, and isotropic pressure are unified into a symmetric rank-2 contravariant tensor field $T^{\mu\nu}$:
Perfect Fluid Stress-Energy Tensor:
For an ideal fluid characterized by proper energy density $\rho(x)$, isotropic pressure $p(x)$, and 4-velocity field $u^\mu = \frac{dx^\mu}{d\tau}$ (normalized such that $g_{\mu\nu} u^\mu u^\nu = -c^2$):
3. Local Conservation Laws
In special relativity (flat Minkowski spacetime $\eta_{\mu\nu}$), conservation of energy and momentum is expressed as $\partial_\mu T^{\mu\nu} = 0$. According to the principle of minimal gravitational coupling (equivalence principle), the physical conservation law in curved spacetime is obtained by replacing ordinary derivatives with covariant derivatives:
Fundamental Conservation Law:
Projecting along and orthogonal to the 4-velocity $u^\nu$ yields the relativistic continuity equation and the relativistic Euler equation of fluid dynamics!
§8.2 Covariant Formulation of Maxwell's Electrodynamics ($F_{\mu\nu}$, $\nabla_\mu F^{\mu\nu} = \mu_0 J^\nu$)
1. The 4-Potential and Field Strength Tensor
Classical electrodynamics is governed by electric and magnetic vector fields $\mathbf{E}$ and $\mathbf{B}$, which mix under Lorentz coordinate transformations. Tensor analysis unifies them into a single antisymmetric rank-2 covariant tensor field $F_{\mu\nu}$, called the Faraday (electromagnetic field) tensor.
Let $A_\mu = (-\phi/c, \mathbf{A})$ be the covariant electromagnetic 4-potential.
Definition 8.1 (Electromagnetic Field Tensor):
Because the connection coefficients $\Gamma^\lambda_{\mu\nu} = \Gamma^\lambda_{\nu\mu}$ are symmetric, the Christoffel symbols cancel identically!
In matrix form with metric signature $(-, +, +, +)$:
By raising indices with $g^{\mu\alpha} g^{\nu\beta}$, we obtain the contravariant tensor $F^{\mu\nu}$.
2. Covariant Maxwell Equations
Theorem 8.1 (Maxwell's Equations in Curved Spacetime): The four classical Maxwell equations reduce to two elegant, manifestly covariant tensor equations:
- Inhomogeneous Maxwell Equations (Sources):
where $J^\nu = (\rho c, \mathbf{J})$ is the 4-current density.
- Homogeneous Maxwell Equations (Bianchi Identity):
Automatic Charge Conservation:
Taking the covariant divergence of the inhomogeneous equation:
because $F^{\mu\nu}$ is antisymmetric while the commutator contraction is symmetric. Thus, electric charge conservation $\nabla_\nu J^\nu = 0$ is an exact geometric identity!
3. The Lorentz Force Density
The electromagnetic 4-force density acting on charged matter is:
which unifies Coulomb electric force and magnetic Lorentz force $\mathbf{F} = q(\mathbf{E} + \mathbf{v} \times \mathbf{B})$.
§8.3 The Geodesic Principle & Gravitational Time Dilation / Redshift
1. The Geodesic Hypothesis
In Einstein's theory of General Relativity, gravity is not a physical Newtonian force, but the curvature of spacetime. Free particles experiencing no non-gravitational forces follow geodesics of the spacetime metric:
where $\tau$ is the proper time experienced by a clock moving along the worldline:
2. Gravitational Time Dilation
Consider a static gravitational field with metric:
For a stationary observer situated at fixed spatial coordinates ($dx^i = 0$):
Let two clocks be stationed at positions $A$ and $B$ where the metric components are $g_{00}(A)$ and $g_{00}(B)$:
In a weak Newtonian gravitational potential $\Phi(x)$, $g_{00} \approx -\left( 1 + \frac{2\Phi}{c^2} \right)$:
Clocks situated deeper in a gravitational well ($\Phi < 0$) run measurably slower than clocks higher up!
3. Gravitational Redshift of Light
A light wave emitted at position $A$ with frequency $\nu_A$ and received at position $B$ undergoes a gravitational frequency shift:
For light escaping from the surface of a star of mass $M$ and radius $R$ ($g_{00} = -\left(1 - \frac{2GM}{c^2 R}\right)$) to a distant observer ($g_{00}(\infty) = -1$):
The light is redshifted ($\Delta \lambda > 0$), directly confirming the curvature of the metric tensor!
§8.4 The Einstein Field Equations $G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}$
1. The Search for the Gravitational Field Equations
In Newtonian gravity, the gravitational potential $\Phi$ satisfies Poisson's equation:
Einstein sought a tensorial equation of the form:
where $\mathcal{G}_{\mu\nu}$ is a symmetric rank-2 geometric tensor constructed from the metric $g_{\mu\nu}$ and its first and second derivatives.
Essential Constraints:
- Since energy-momentum is conserved ($\nabla^\mu T_{\mu\nu} = 0$), the geometric tensor must satisfy the divergence-free identity:
- In 1915, David Hilbert and Albert Einstein proved that the unique symmetric tensor containing at most second derivatives of the metric that satisfies this condition (Lovelock's theorem in 4D) is:
2. The Einstein Field Equations
Definition 8.2 (The Einstein Field Equations):
where:
- $G_{\mu\nu} = R_{\mu\nu} - \frac{1}{2} R g_{\mu\nu}$ is the Einstein tensor.
- $\Lambda$ is the cosmological constant.
- $G = 6.674 \times 10^{-11} \, \text{m}^3 \text{kg}^{-1} \text{s}^{-2}$ is Newton's gravitational constant.
- $c$ is the speed of light.
- $\kappa = \frac{8\pi G}{c^4} \approx 2.076 \times 10^{-43} \, \text{N}^{-1}$ is Einstein's gravitational coupling constant.
3. Trace-Reversed Form and Vacuum Equations
Taking the trace with $g^{\mu\nu}$:
Substituting back into $G_{\mu\nu}$:
Vacuum Field Equations ($\Lambda = 0, T_{\mu\nu} = 0$):
In empty spacetime outside mass distributions:
Spacetimes satisfying $R_{\mu\nu} = 0$ are called Ricci-flat. Crucially, $R_{\mu\nu} = 0$ does NOT imply flat spacetime! The full Riemann tensor $R^\rho_{\mu\sigma\nu}$ can still be non-zero through its Weyl tensor component, propagating gravitational waves and mediating gravitational attraction!
§8.5 The Schwarzschild Metric, Planetary Perihelion Precession & Gravitational Deflection
1. The Schwarzschild Solution (1916)
Karl Schwarzschild derived the exact vacuum solution ($R_{\mu\nu} = 0$) for a static, spherically symmetric mass $M$:
Definition 8.3 (Schwarzschild Metric):
where $r_s = \frac{2GM}{c^2}$ is the Schwarzschild radius (event horizon). For the Sun: $r_s \approx 2.95 \text{ km}$; for the Earth: $r_s \approx 8.87 \text{ mm}$.
2. Relativistic Planetary Orbits and Perihelion Precession
Consider motion confined to the equatorial plane $\theta = \frac{\pi}{2}$. From the geodesic equations, energy per unit mass $E$ and angular momentum per unit mass $L$ are conserved:
Substituting into the metric $g_{\mu\nu} \dot{x}^\mu \dot{x}^\nu = -c^2$:
where the relativistic effective potential is:
The third term $-\frac{GML^2}{c^2 r^3}$ is the general relativistic correction to Newtonian gravity.
Using $u = 1/r$, the orbital Binet equation becomes:
By perturbation theory, an elliptical orbit precesses by an angle per revolution:
For Mercury: $a = 5.79 \times 10^{10} \text{ m}, e = 0.2056$. This yields $\Delta\phi = 42.98'' \text{ per century}$, matching the observed anomalous 43 arcseconds per century unexplained by Newtonian planetary perturbations!
3. Gravitational Deflection of Starlight
For null geodesics (photons, $ds^2 = 0$), the orbital equation with impact parameter $b$ gives:
Integrating the perturbation from $\phi = -\pi/2$ to $\pi/2$ for light grazing the Sun at radius $R$:
Deflection Angle Formula:
For the Sun ($M = M_\odot, R = R_\odot$):
Exactly twice the naive Newtonian corpuscular value ($0.875''$)! This famous prediction was confirmed by Arthur Eddington's 1919 solar eclipse expedition, providing historic observational confirmation of general relativity.
Rigorous Tiered Solved Examination Problems
Step-by-step unskipped derivations, complete proofs, and verification across Foundational, Advanced, and Honors tiers.
Consider a pressureless cloud of non-interacting particles ('dust') in spacetime with proper mass density $\rho_0$ and 4-velocity field $u^\mu = \frac{dx^\mu}{d\tau}$. Its stress-energy tensor is given by:
- Verify that $T^{\mu\nu}$ is symmetric.
- In the rest frame of the dust, where $u^\mu = (c, 0, 0, 0)$ in Minkowski spacetime $\eta_{\mu\nu} = \text{diag}(-1, 1, 1, 1)$, write down the matrix components of $T^{\mu\nu}$ and $T^\mu_\nu$.
- Compute the scalar trace $T = g_{\mu\nu} T^{\mu\nu}$.
- Show that the conservation law $\nabla_\mu T^{\mu\nu} = 0$ decomposes into the continuity equation $\nabla_\mu (\rho_0 u^\mu) = 0$ and the geodesic equation $u^\mu \nabla_\mu u^\nu = 0$.
1. Symmetry
Symmetric by definition. $\blacksquare$
2. Matrix Components in Rest Frame
With $u^\mu = (c, 0, 0, 0)$:
Lowering the second index with $\eta_{\nu\lambda} = \text{diag}(-1, 1, 1, 1)$:
3. Trace $T$
4. Decomposition of Conservation Law
Expand $\nabla_\mu T^{\mu\nu} = 0$:
Contract this equation with the 4-velocity covector $u_\nu$:
Since $u^\nu u_\nu = -c^2 = \text{const}$:
Thus, the second term vanishes identically:
This is the relativistic mass-continuity equation (conservation of particle number).
Now substitute $\nabla_\mu (\rho_0 u^\mu) = 0$ back into the original conservation equation:
This is the geodesic equation for the dust worldlines $\frac{D u^\nu}{d\tau} = 0$! Energy-momentum conservation automatically forces pressureless dust particles to follow spacetime geodesics! $\blacksquare$
The electromagnetic stress-energy tensor in 4D spacetime is defined as:
- Prove that $T^{\mu\nu}_{\text{EM}}$ is strictly symmetric.
- Prove that in $n = 4$ spacetime dimensions, the trace of $T^{\mu\nu}_{\text{EM}}$ vanishes identically: $T = g_{\mu\nu} T^{\mu\nu}_{\text{EM}} = 0$.
- Use Maxwell's equations $\nabla_\mu F^{\mu\nu} = \mu_0 J^\nu$ and $\nabla_{[\lambda} F_{\mu\nu]} = 0$ to calculate $\nabla_\mu T^{\mu\nu}_{\text{EM}}$ and show that $\nabla_\mu T^{\mu\nu}_{\text{EM}} = -F^{\nu\alpha} J_\alpha$.
- What physical principle does $\nabla_\mu (T^{\mu\nu}_{\text{matter}} + T^{\mu\nu}_{\text{EM}}) = 0$ express?
1. Proof of Symmetry
Using antisymmetry $F^{\mu\alpha} = -F^{\alpha\mu}$ and $F^{\nu\beta} = -F^{\beta\nu}$:
Thus $T^{\mu\nu}_{\text{EM}} = T^{\nu\mu}_{\text{EM}}$. $\blacksquare$
2. Trace-Free Property in 4 Dimensions
Evaluate each term:
- $g_{\mu\nu} F^{\mu\alpha} F^\nu_{\ \alpha} = F_\nu^{\ \alpha} F^\nu_{\ \alpha} = F^{\alpha\nu} F_{\alpha\nu} = F_{\alpha\beta} F^{\alpha\beta}$.
- In 4D, $g_{\mu\nu} g^{\mu\nu} = \delta^\mu_\mu = 4$.
Substituting:
Because $T \equiv 0$, pure radiation has zero gravitational trace, meaning $R = 0$ in the presence of electromagnetic fields!
3. Covariant Divergence
Using $\nabla_\mu F^{\mu\alpha} = \mu_0 J^\alpha$:
The remaining terms cancel by the homogeneous Maxwell equation $\nabla_{[\mu} F_{\alpha\beta]} = 0$:
Therefore:
4. Physical Conservation Principle
Since $\nabla_\mu T^{\mu\nu}_{\text{matter}} = +F^{\nu\alpha} J_\alpha$ (Lorentz force acting on matter):
This expresses total energy-momentum conservation: whatever energy and momentum the electromagnetic field loses is transferred into kinetic energy and momentum of the charged matter! $\blacksquare$
Consider a general static, spherically symmetric spacetime metric in spherical coordinates $(t, r, \theta, \phi) = (x^0, x^1, x^2, x^3)$:
where $\Phi(r)$ and $\Lambda(r)$ are unknown radial metric functions.
- Compute the non-zero Christoffel symbols $\Gamma^\lambda_{\mu\nu}$ for this metric.
- Compute the Ricci tensor components $R_{00}$, $R_{11}$, and $R_{22}$.
- Impose the vacuum Einstein field equations $R_{\mu\nu} = 0$:
- Combine $R_{00}$ and $R_{11}$ to show that $\Phi'(r) + \Lambda'(r) = 0 \implies \Phi(r) = -\Lambda(r) + \text{const}$.
- Use $R_{22} = 0$ to solve for $e^{-2\Lambda(r)}$.
- Determine the integration constant from the Newtonian weak-field limit and arrive at the exact Schwarzschild metric.
1. Christoffel Symbols
Using $\Gamma^\lambda_{\mu\nu} = \frac{1}{2} g^{\lambda\sigma} (\partial_\mu g_{\nu\sigma} + \partial_\nu g_{\mu\sigma} - \partial_\sigma g_{\mu\nu})$:
- $\Gamma^0_{01} = \Phi'$, $\Gamma^1_{00} = c^2 \Phi' e^{2(\Phi - \Lambda)}$
- $\Gamma^1_{11} = \Lambda'$, $\Gamma^1_{22} = -r e^{-2\Lambda}$, $\Gamma^1_{33} = -r \sin^2\theta e^{-2\Lambda}$
- $\Gamma^2_{12} = \frac{1}{r}$, $\Gamma^2_{33} = -\sin\theta\cos\theta$
- $\Gamma^3_{13} = \frac{1}{r}$, $\Gamma^3_{23} = \cot\theta$
2. Ricci Tensor Components
Evaluating $R_{\mu\nu} = \partial_\lambda \Gamma^\lambda_{\mu\nu} - \partial_\nu \Gamma^\lambda_{\mu\lambda} + \Gamma^\lambda_{\mu\nu} \Gamma^\sigma_{\lambda\sigma} - \Gamma^\lambda_{\mu\sigma} \Gamma^\sigma_{\nu\lambda}$:
3. Vacuum Field Equations $R_{\mu\nu} = 0$
Step A: Sum of $R_{00}$ and $R_{11}$
Consider the linear combination:
Since $r \ne 0$:
As $r \to \infty$, spacetime must be flat Minkowski space: $\Phi(\infty) = 0$ and $\Lambda(\infty) = 0 \implies C_0 = 0$. Thus:
Step B: Solve $R_{22} = 0$
Substituting $\Phi' = -\Lambda'$ into $R_{22}$:
Notice that:
Therefore, $R_{22} = 0$ is equivalent to:
Integrating with respect to $r$:
where $r_s$ is a constant of integration. Since $e^{2\Phi} = e^{-2\Lambda}$:
4. Determination of Integration Constant $r_s$
At large distances $r \gg r_s$, the Newtonian limit requires:
Comparing with $e^{2\Phi} = 1 - \frac{r_s}{r}$:
Substituting back into the line element:
This completes the exact, closed-form derivation of the Schwarzschild Metric! $\blacksquare$