Physics Derivations

1. Coordinates and Galilean Reference Frames

To describe physical phenomena, we construct a coordinate system. A point particle's position is represented by a vector $\mathbf{r}(t) = (x(t), y(t), z(t))$ in Cartesian space relative to an origin. Galileo observed that there is no absolute resting reference frame; uniform linear motion is fundamentally relative. If an observer in reference frame $S'$ moves at a constant velocity $\mathbf{v}$ along the x-axis relative to a "stationary" frame $S$, the coordinates transform via the Galilean Transformations:

$x' = x - vt, \quad y' = y, \quad z' = z, \quad t' = t$

Differentiating these positions with respect to time yields the velocity addition rule: $u'_x = u_x - v$. Differentiating once more, because $\mathbf{v}$ is strictly constant, its derivative vanishes ($\frac{d\mathbf{v}}{dt} = 0$), showing that physical accelerations are identical in all Galilean inertial reference frames: $\mathbf{a}' = \mathbf{a}$. This mathematical equivalence is the cornerstone of the Galilean relativity principle: the laws of mechanics are identical for all observers in uniform motion.

2. Defining the Concept of Force: Galileo's Inertia to Newton's Second Law

In Aristotelian physics, objects were believed to naturally come to rest, requiring a continuous push to maintain a velocity ($\mathbf{F} \propto \mathbf{v}$). Galileo dismantled this through rigorous inclined-plane and double-ramp rolling tests: he observed that a polished bronze ball rolled down one ramp would rise to the exact same height on an opposite ramp. When he flattened the second ramp (reducing its angle to 0°), the ball rolled indefinitely, limited only by the roughness of the wood. This meant objects possess a physical property called **inertia**—they naturally preserve their state of motion.

Newton codified this into three laws, translating qualitative observations into rigorous algebraic dynamics:

• Newton's First Law (Inertia): A body remains at rest or in uniform linear motion unless compelled to change its state by a net external force. This defines the existence of inertial reference frames.
• Newton's Second Law (Dynamics): How do we formulate an equation to predict the exact acceleration caused by a force?
  Let us perform an experiment with a rolling cart. If we pull the cart using a spring stretched to a constant length, it experiences a constant force $\mathbf{F}$, and we observe that it accelerates at a constant rate $\mathbf{a}$.
  If we double the mass of the cart ($m \to 2m$) by stacking an identical block on it, we observe that the same spring pull produces exactly **half** the acceleration ($\mathbf{a} \to \mathbf{a}/2$). This means acceleration is inversely proportional to mass: $\mathbf{a} \propto 1/m$.
  Now, if we double the force ($F \to 2\mathbf{F}$) by pulling with two identical springs side-by-side, we observe that the acceleration of the cart **doubles** ($\mathbf{a} \to 2\mathbf{a}$). This means acceleration is directly proportional to force: $\mathbf{a} \propto \mathbf{F}$.
  Combining these experimental observations yields: $$ \mathbf{a} \propto \frac{\mathbf{F}}{m} \implies \mathbf{F} \propto m\mathbf{a} $$ By defining the unit of force (the "Newton") as the exact amount of force required to accelerate a $1\text{ kg}$ mass at $1\text{ m/s}^2$, the constant of proportionality becomes exactly $1$. Generalizing force as the rate of change of linear momentum $\mathbf{p} = m\mathbf{v}$, we write: $$ \mathbf{F} = \frac{d\mathbf{p}}{dt} = \frac{d(m\mathbf{v})}{dt} = m\mathbf{a} \quad (\text{for constant mass } m) $$
• Newton's Third Law (Reciprocity): When body $A$ exerts a force $\mathbf{F}_{AB}$ on body $B$, body $B$ simultaneously exerts an equal and opposite force $\mathbf{F}_{BA}$ on body $A$: $$ \mathbf{F}_{AB} = -\mathbf{F}_{BA} $$ This is experimentally demonstrated by two interacting carts with a compressed spring between them; when released, the change in their momenta is equal in magnitude and opposite in direction ($\Delta \mathbf{p}_A = -\Delta \mathbf{p}_B$), indicating that the forces acting during the collision were symmetric.

3. Direct Derivations of Kinematics: Linking Galileo's $s \propto t^2$ to Calculus

Through his inclined plane experiments, Galileo observed a precise quantitative rule: the distance $s$ rolled by a ball starting from rest was strictly proportional to the square of the elapsed time ($s \propto t^2$). If we assume a constant force causing a constant physical acceleration $\mathbf{a}(t) = \mathbf{a}_0$, how do we mathematically formulate this observation?

Since acceleration is defined as the rate of change of velocity ($\mathbf{a}_0 = \frac{d\mathbf{v}}{dt}$), we separate variables and integrate both sides from initial time $t = 0$ to final time $t$: $$ \int_{\mathbf{v}_0}^{\mathbf{v}(t)} d\mathbf{v} = \int_0^t \mathbf{a}_0 dt' \implies \mathbf{v}(t) - \mathbf{v}_0 = \mathbf{a}_0 t \implies \mathbf{v}(t) = \mathbf{v}_0 + \mathbf{a}_0 t $$ This shows that velocity increases linearly with time. Now, since velocity is the rate of change of position ($\mathbf{v}(t) = \frac{d\mathbf{r}}{dt}$), we integrate once more: $$ \int_{\mathbf{r}_0}^{\mathbf{r}(t)} d\mathbf{r} = \int_0^t (\mathbf{v}_0 + \mathbf{a}_0 t') dt' \implies \mathbf{r}(t) - \mathbf{r}_0 = \left[ \mathbf{v}_0 t' + \frac{1}{2}\mathbf{a}_0 t'^2 \right]_0^t \implies \mathbf{r}(t) = \mathbf{r}_0 + \mathbf{v}_0 t + \frac{1}{2}\mathbf{a}_0 t^2 $$ This calculus derivation beautifully reproduces Galileo's experimental discovery: if an object starts from rest ($\mathbf{v}_0 = 0$) at the origin ($\mathbf{r}_0 = 0$), its position is $\mathbf{r}(t) = \frac{1}{2}\mathbf{a}_0 t^2$, mathematically proving that $s \propto t^2$.

To obtain the 1D speed-displacement relationship without explicit time variables, we write acceleration using the chain rule as $a = \frac{dv}{dt} = \frac{dv}{dx}\frac{dx}{dt} = v\frac{dv}{dx}$. Separating variables and integrating: $$ \int_{x_0}^x a \, dx' = \int_{v_0}^v v' \, dv' \implies a(x - x_0) = \left[ \frac{1}{2}v'^2 \right]_{v_0}^v \implies a\Delta x = \frac{1}{2}v^2 - \frac{1}{2}v_0^2 \implies v^2 = v_0^2 + 2a\Delta x $$

4. Energy and Momentum Conservation

We define the mechanical **Work** $W$ done by a force on a particle traveling along a path $C$ as the line integral: $$ W = \int_C \mathbf{F} \cdot d\mathbf{r} = \int_{t_0}^{t_1} \mathbf{F} \cdot \frac{d\mathbf{r}}{dt} dt = \int_{t_0}^{t_1} \mathbf{F} \cdot \mathbf{v} dt $$

Substituting Newton's second law $\mathbf{F} = m\frac{d\mathbf{v}}{dt}$: $$ W = \int_{t_0}^{t_1} m\frac{d\mathbf{v}}{dt} \cdot \mathbf{v} dt = \int_{\mathbf{v}_0}^{\mathbf{v}_1} m \mathbf{v} \cdot d\mathbf{v} = \left[ \frac{1}{2}m(\mathbf{v} \cdot \mathbf{v}) \right]_{\mathbf{v}_0}^{\mathbf{v}_1} = \frac{1}{2}mv_1^2 - \frac{1}{2}mv_0^2 = \Delta K $$ This is the **Work-Energy Theorem**: the work done by a net force equals the change in the particle's **Kinetic Energy** $K$.

If the force is conservative (path-independent), we can define a scalar **Potential Energy** function $U(\mathbf{r})$ such that $\mathbf{F} = -\nabla U$. Thus: $$ W = \int_{\mathbf{r}_0}^{\mathbf{r}_1} \mathbf{F} \cdot d\mathbf{r} = -\int_{\mathbf{r}_0}^{\mathbf{r}_1} \nabla U \cdot d\mathbf{r} = -(U_1 - U_0) = -\Delta U $$ Combining the two definitions: $\Delta K = -\Delta U \implies \Delta(K + U) = 0 \implies E = K + U = \text{Constant}$.

5. Gravitational Fields & Kepler's Laws: Deducing the Inverse-Square Law

Terrestrial gravity causes objects to fall with a constant acceleration of $g \approx 9.8\text{ m/s}^2$. In the celestial realm, Johannes Kepler synthesized decades of meticulous astronomical observations of planetary orbits, identifying a striking empirical harmony (Kepler's Third Law): the square of a planet's orbital period $T$ is directly proportional to the cube of its orbital radius ($T^2 \propto r^3$).

How did Isaac Newton formulate a single gravitational equation to explain both phenomena?
Let us model a planet of mass $m$ in a stable circular orbit of radius $r$ around the sun. The centripetal force holding the planet in its orbit is: $$ F_c = \frac{m v^2}{r} $$ Since the planet's orbital speed is the circumference divided by the period ($v = \frac{2\pi r}{T}$), its square is $v^2 = \frac{4\pi^2 r^2}{T^2}$.
Substituting Kepler's experimental observation that $T^2 = k r^3$ (where $k$ is a constant of nature): $$ v^2 = \frac{4\pi^2 r^2}{k r^3} = \left(\frac{4\pi^2}{k}\right) \frac{1}{r} $$ This indicates that planetary speed squared drops off inversely with distance. Plugging this speed back into our centripetal force equation: $$ F_c = \frac{m v^2}{r} \propto \frac{m \left(\frac{1}{r}\right)}{r} = \frac{m}{r^2} $$ To explain Kepler's astronomical observations, the gravitational force **must** decay with the inverse square of the distance: $F \propto 1/r^2$.
Furthermore, by Newton's Third Law (Reciprocity), gravity must be mutual: if the Sun pulls the planet, the planet pulls the Sun. The force must therefore be proportional to both the Sun's mass $M$ and the planet's mass $m$. Combining these experimental and logical constraints yields: $$ \mathbf{F} = -G \frac{M m}{r^2} \hat{\mathbf{r}} $$ where $G$ is the universal gravitational constant. Applying this central force model to Newton's Second Law in general elliptical coordinates, we mathematically derive Kepler's three empirical laws:

• Kepler's First Law (Law of Ellipses): The orbit of a planet is an ellipse with the Sun at one of the two foci.
  To derive this, we set up the equation of motion in plane polar coordinates. The radial component of acceleration is $a_r = \ddot{r} - r\dot{\theta}^2$. Under gravitational attraction: $$ m(\ddot{r} - r\dot{\theta}^2) = -G\frac{Mm}{r^2} \implies \ddot{r} - r\dot{\theta}^2 = -\frac{GM}{r^2} $$ Since the central force is radial, angular momentum $L = mr^2\dot{\theta}$ is conserved, which means $\dot{\theta} = \frac{L}{mr^2}$. Making the standard orbital substitution $u = 1/r$ and writing $r$ derivatives in terms of angle $\theta$: $$ \frac{dr}{dt} = \frac{d}{dt}\left(\frac{1}{u}\right) = -\frac{1}{u^2}\frac{du}{d\theta}\frac{d\theta}{dt} = -\frac{1}{u^2}\frac{du}{d\theta}\left(\frac{Lu^2}{m}\right) = -\frac{L}{m}\frac{du}{d\theta} $$ Differentiating again with respect to time: $$ \ddot{r} = -\frac{L}{m}\frac{d^2u}{d\theta^2}\dot{\theta} = -\frac{L^2u^2}{m^2}\frac{d^2u}{d\theta^2} $$ Substituting $\ddot{r}$ and $\dot{\theta}$ back into the radial equation of motion yields Binet's orbital equation: $$ -\frac{L^2u^2}{m^2}\frac{d^2u}{d\theta^2} - \frac{1}{u}\left(\frac{Lu^2}{m}\right)^2 = -GMu^2 \implies \frac{d^2u}{d\theta^2} + u = \frac{GMm^2}{L^2} $$ This is a second-order linear differential equation. Its general physical solution is: $$ u(\theta) = \frac{GMm^2}{L^2} \left[ 1 + e \cos(\theta - \theta_0) \right] $$ Inverting back to $r(\theta)$ and setting $\theta_0 = 0$: $$ r(\theta) = \frac{\alpha}{1 + e\cos\theta} \quad \text{where} \quad \alpha = \frac{L^2}{GMm^2} $$ This represents the polar coordinate equation of a conic section with eccentricity $e$. For bounded planetary systems where $0 \le e < 1$, this equation dictates an **ellipse**, proving Kepler's First Law.
• Kepler's Second Law (Law of Equal Areas): A line segment joining a planet and the Sun sweeps out equal areas during equal intervals of time.
  The differential area $dA$ swept out by the position vector $\mathbf{r}$ as it rotates through a small angle $d\theta$ is: $$ dA = \frac{1}{2} r (r d\theta) = \frac{1}{2} r^2 d\theta $$ Dividing by $dt$ to find the rate of area swept over time (areal velocity): $$ \frac{dA}{dt} = \frac{1}{2} r^2 \frac{d\theta}{dt} $$ Since angular momentum is $L = m r^2 \dot{\theta}$, we substitute $\dot{\theta} = \frac{L}{mr^2}$: $$ \frac{dA}{dt} = \frac{1}{2} r^2 \left( \frac{L}{mr^2} \right) = \frac{L}{2m} = \text{Constant} $$ Because angular momentum $L$ is conserved in a central field, $\frac{dA}{dt}$ is strictly constant, demonstrating that planets sweep out equal orbital areas in equal times.
• Kepler's Third Law (Law of Harmonies): The square of the orbital period $T$ of a planet is directly proportional to the cube of the semi-major axis $a$ of its orbit.
  Integrating the constant areal velocity $\frac{dA}{dt} = \frac{L}{2m}$ over one full orbital period $T$ must yield the total area of the ellipse $A = \pi a b$, where $b$ is the semi-minor axis: $$ A = \int_0^T \frac{dA}{dt} dt = \frac{L}{2m} T \implies \pi a b = \frac{L}{2m} T \implies T = \frac{2m\pi ab}{L} $$ Squaring both sides gives: $$ T^2 = \frac{4\pi^2 m^2 a^2 b^2}{L^2} $$ The geometry of an ellipse relates the axes by $b^2 = a^2(1 - e^2)$. Recall from the First Law derivation that the semi-latus rectum is $\alpha = \frac{L^2}{GMm^2} = a(1-e^2)$. Rearranging terms: $$ 1 - e^2 = \frac{L^2}{GMm^2 a} \implies b^2 = a^2 \left( \frac{L^2}{GMm^2 a} \right) = \frac{a L^2}{GMm^2} \implies \frac{m^2 b^2}{L^2} = \frac{a}{GM} $$ Substituting this relationship into our squared period expression: $$ T^2 = 4\pi^2 a^2 \left( \frac{m^2 b^2}{L^2} \right) = 4\pi^2 a^2 \left( \frac{a}{GM} \right) = \left( \frac{4\pi^2}{GM} \right) a^3 $$ Thus, we derive that the square of the orbital period is proportional to the cube of the semi-major axis ($T^2 \propto a^3$), completing the mechanical synthesis.

1. The Speed of Light Paradox

Maxwell's wave equation derived the velocity of light $c = 1/\sqrt{\mu_0\epsilon_0}$ strictly from universal parameters. However, in Galilean reference frames, speeds are relative ($c' = c - v$). Under Galilean relativity, there could only be *one* privileged frame where the speed of light is exactly $c$—the frame of the medium carrying light waves, dubbed the **Luminiferous Aether**.

2. Michelson-Morley Interferometer (1887)

Albert A. Michelson and Edward W. Morley attempted to measure the velocity of Earth through this Aether wind by splitting a light beam into two perpendicular arms, bouncing them off mirrors, and recombining them to form interference fringes.

Let one arm (length $L_1$) lie parallel to Earth's velocity vector $\mathbf{v}$ through the Aether, and the other arm (length $L_2$) lie perpendicular. The classical time of round-trip travel along the parallel arm is calculated via velocity addition: $$ t_{\text{para}} = \frac{L_1}{c - v} + \frac{L_1}{c + v} = \frac{2 L_1 c}{c^2 - v^2} = \frac{2 L_1}{c} \frac{1}{1 - v^2/c^2} $$ While the cross-wind travel along the perpendicular arm requires diagonal aiming, yielding: $$ t_{\text{perp}} = \frac{2 L_2}{\sqrt{c^2 - v^2}} = \frac{2 L_2}{c} \frac{1}{\sqrt{1 - v^2/c^2}} $$

As the Earth rotated, the phase difference between the arms should shift, shifting the interference pattern. However, despite extreme sensitivity, **no fringe shift was detected**. The velocity of Earth relative to the Aether was zero at all seasons, showing that the speed of light is isotropic (constant in all directions).

3. Postulates of Special Relativity

Albert Einstein resolved this crisis by discarding the Aether and reforming the concepts of space and time. He established two postulates:

1. The Principle of Relativity: The laws of physics are identical in all inertial reference frames.
2. The Constancy of the Speed of Light: Light propagates in empty space with a constant velocity $c$ independent of the motion of the emitting body.

4. Derivation of Time Dilation via Light-Clock

Consider a light clock consisting of two parallel mirrors separated by a distance $D$ in a stationary frame $S_0$. A light pulse bounces between them. The proper time interval $\Delta t_0$ for a round-trip bounce is: $$ \Delta t_0 = \frac{2D}{c} $$

Now, let this clock move horizontally at constant velocity $v$ relative to an observer in frame $S$. To this observer, the light pulse must travel a diagonal path of length $2L$ to trace the moving mirrors. In the time interval $\Delta t$ of the round trip, the clock travels a horizontal distance $v \Delta t$.

Applying the Pythagorean theorem to the triangle formed by the light path: $$ L^2 = D^2 + \left( \frac{v \Delta t}{2} \right)^2 $$

Since the speed of light is constant ($c$) in frame $S$, the distance $L$ traveled by the pulse is $L = c \frac{\Delta t}{2}$. Substituting: $$ \left( \frac{c \Delta t}{2} \right)^2 = D^2 + \left( \frac{v \Delta t}{2} \right)^2 $$

Multiplying by 4 and rearranging terms to isolate $\Delta t$: $$ c^2 \Delta t^2 - v^2 \Delta t^2 = 4D^2 \implies \Delta t^2(c^2 - v^2) = 4D^2 \implies \Delta t^2 = \frac{4D^2}{c^2 - v^2} = \frac{4D^2}{c^2 \left(1 - \frac{v^2}{c^2}\right)} $$ Taking the square root and substituting proper time $\Delta t_0 = \frac{2D}{c}$: $$ \Delta t = \frac{2D}{c} \frac{1}{\sqrt{1 - v^2/c^2}} = \gamma \Delta t_0 $$ where $\gamma = 1/\sqrt{1 - v^2/c^2}$ is the relativistic **Lorentz Factor**.

5. The Lorentz Transformations: Deriving Spacetime Transformations from Constant $c$

To ensure that the speed of light $c$ is strictly invariant for all inertial observers, the Galilean transformations must be discarded. How can we formulate an equation to transform coordinates while preserving the isotropic speed of light?
Let us assume a linear coordinate transformation between frame $S$ and frame $S'$ moving at relative velocity $v$ along the x-axis. To ensure uniform motion transforms to uniform motion, the transformation must be linear: $$ x' = \gamma (x - vt) $$ By the principle of relativity, the inverse transformation must carry the exact same mathematical form, differing only by the sign of the relative velocity: $$ x = \gamma (x' + vt') $$ Now, let us perform a thought experiment. A light pulse is emitted from the origin of both frames at the instant they coincide ($t = t' = 0$). Since the speed of light is observed as exactly $c$ in **both** frames, the position of the wave front must satisfy: $$ x = ct \quad \text{and} \quad x' = ct' $$ Let us substitute these isotropic wave equations into our transformation equations: $$ ct' = \gamma (c - v) t \quad \text{and} \quad ct = \gamma (c + v) t' $$ To solve for the scaling factor $\gamma$, we multiply these two expressions together: $$ c^2 t t' = \gamma^2 (c^2 - v^2) t t' $$ Dividing both sides by $t t'$ and isolating $\gamma$: $$ c^2 = \gamma^2 (c^2 - v^2) \implies \gamma^2 = \frac{c^2}{c^2 - v^2} = \frac{1}{1 - v^2/c^2} \implies \gamma = \frac{1}{\sqrt{1 - v^2/c^2}} $$ This is the exact same Lorentz factor derived from our geometric light-clock, proving that the constant speed of light mathematically necessitates this spatial stretch.
To find how time transforms ($t'$), we substitute the spatial equation $x' = \gamma(x-vt)$ into our reciprocal equation $x = \gamma(x' + vt')$: $$ x = \gamma \left[ \gamma(x - vt) + vt' \right] = \gamma^2 x - \gamma^2 vt + \gamma vt' $$ Isolating the $t'$ term: $$ \gamma vt' = x - \gamma^2 x + \gamma^2 vt \implies vt' = \frac{x}{\gamma} - \gamma x + \gamma vt $$ Dividing by $v$ and substituting $\frac{1}{\gamma} = \gamma(1 - v^2/c^2)$: $$ t' = \gamma t + \frac{x}{v} \left( \frac{1}{\gamma} - \gamma \right) = \gamma t + \frac{x}{v} \left( \gamma\left(1 - \frac{v^2}{c^2}\right) - \gamma \right) = \gamma t + \frac{x}{v} \left( \gamma - \gamma\frac{v^2}{c^2} - \gamma \right) $$ Simplifying the expression: $$ t' = \gamma\left(t - \frac{vx}{c^2}\right) $$ This yields the complete and elegant **Lorentz Transformations**:

$x' = \gamma(x - vt), \quad y' = y, \quad z' = z, \quad t' = \gamma\left(t - \frac{vx}{c^2}\right)$

These equations demonstrate that space and time are not independent absolutes, but rather orthogonal coordinates in a unified 4D spacetime continuum. A consequence of this is the relativity of simultaneity: events that occur simultaneously at different spatial locations in frame $S$ will not occur at the same time in frame $S'$.

6. 4-Vectors: Unifying Spacetime and Energy-Momentum

In special relativity, the separate three-dimensional vectors of space and scalar quantities of time are unified into four-dimensional vectors, or **4-vectors**, which transform under Lorentz boosts. The position 4-vector $X^\mu$ is defined as: $$ X^\mu = (X^0, X^1, X^2, X^3) = (ct, x, y, z) $$ where $\mu = 0,1,2,3$ denotes the spacetime components ($0$ is temporal, $1,2,3$ are spatial). Under a Lorentz boost at velocity $v$ along the x-axis, this vector transforms as $X'^\mu = \Lambda^\mu_{\ \nu} X^\nu$, which is exactly the matrix representation of the Lorentz transformations: $$ \begin{pmatrix} ct' \\ x' \\ y' \\ z' \end{pmatrix} = \begin{pmatrix} \gamma & -\beta\gamma & 0 & 0 \\ -\beta\gamma & \gamma & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} ct \\ x \\ y \\ z \end{pmatrix} $$ where $\beta = v/c$. To compute physical quantities that are identical for all observers, we define the scalar product of two 4-vectors using the flat **Minkowski metric** $\eta_{\mu\nu} = \text{diag}(-1, 1, 1, 1)$: $$ A \cdot B = \eta_{\mu\nu} A^\mu B^\nu = -A^0 B^0 + A^1 B^1 + A^2 B^2 + A^3 B^3 $$ The scalar product of the position 4-vector with itself yields the **invariant spacetime interval** $ds^2$: $$ ds^2 = X \cdot X = -c^2 dt^2 + dx^2 + dy^2 + dz^2 $$ Because $\eta_{\mu\nu}$ is invariant under Lorentz boosts, $ds^2$ has the exact same value in all inertial reference frames ($ds^2 = ds'^2$), representing the absolute geometry of spacetime.

Similarly, the three-dimensional momentum vector $\mathbf{p}$ and scalar relativistic energy $E$ are unified into the **energy-momentum 4-vector** $P^\mu$: $$ P^\mu = \left(\frac{E}{c}, p_x, p_y, p_z\right) $$ Taking the Lorentz-invariant scalar product of $P^\mu$ with itself: $$ P \cdot P = -\frac{E^2}{c^2} + |\mathbf{p}|^2 = -m^2 c^2 $$ where $m$ is the particle's invariant rest mass. Rearranging this yields Einstein's celebrated energy-momentum relation: $$ E^2 = (pc)^2 + (mc^2)^2 $$ For a particle at rest ($\mathbf{p} = 0$), this simplifies directly to the most famous equation in physics: $E = mc^2$. Thus, 4-vectors provide the ultimate coordinate-free representation of relativistic dynamics.

1. Cavity Radiators and the Classical Failure: Deducing the Density of States

At the close of the 19th century, physicists sought to compute the exact spectrum of electromagnetic radiation emitted by a thermal cavity (an ideal blackbody). Lord Rayleigh and Sir James Jeans applied classical statistical mechanics to the radiation field inside a conducting cavity of volume $V = L^3$.
They treated the radiation field as a collection of electromagnetic standing waves with nodes at the metallic walls. By solving Maxwell's equations under these conducting boundary conditions, they deduced the number of electromagnetic standing wave modes per unit frequency interval $d\nu$ (the **Density of States**): $$ g(\nu)d\nu = \frac{8\pi V \nu^2}{c^3} d\nu $$

According to the classical statistical **Equipartition Theorem**, every degree of freedom in thermal equilibrium at temperature $T$ must carry an average energy of: $$ \langle E \rangle = k_B T $$ where $k_B$ is the Boltzmann constant. This continuous average energy is derived by integrating the Boltzmann probability distribution over all possible continuous energy states from $0$ to $\infty$: $$ \langle E \rangle = \frac{\int_0^\infty E e^{-\frac{E}{k_B T}} dE}{\int_0^\infty e^{-\frac{E}{k_B T}} dE} = k_B T $$

Multiplying the density of states $g(\nu)$ by the continuous average energy $k_B T$ yields the classical **Rayleigh-Jeans Spectral Radiance** (expressed in wavelength $\lambda = c/\nu$): $$ I(\lambda, T) = \frac{8\pi k_B T}{\lambda^4} $$

2. Planck's Quanta Hypothesis and the Resolution: Formulating Discrete Averages

In 1900, Max Planck resolved this crisis by introducing a radical physical assumption. He realized that to suppress the infinite divergence at high frequencies, there must be a mechanism that prevents high-frequency electromagnetic oscillators from absorbing thermal energy at a given temperature $T$.
If energy is continuous, any oscillator—no matter how high its frequency—can absorb an infinitesimally small amount of energy to participate in thermal equilibrium. But what if energy is strictly **discrete**?
Planck postulated that the microscopic oscillators in the cavity walls exchange energy with the radiation field *only* in discrete packets called **quanta**: $$ E = n h \nu = n \frac{h c}{\lambda}, \quad n \in \{0, 1, 2, 3, \dots\} $$ where $h$ is a new constant of nature (Planck's constant).
Under this discrete assumption, the average energy of an oscillator mode must be computed using a discrete sum over Maxwell-Boltzmann probabilities rather than a continuous integration: $$ \langle E \rangle = \frac{\sum_{n=0}^\infty (nh\nu) e^{-\frac{nh\nu}{k_B T}}}{\sum_{n=0}^\infty e^{-\frac{nh\nu}{k_B T}}} $$ Let us define a dimensionless parameter $x = \frac{h\nu}{k_B T}$. The denominator of this expression is a standard infinite geometric series: $$ S = \sum_{n=0}^\infty e^{-nx} = 1 + e^{-x} + e^{-2x} + \dots = \frac{1}{1 - e^{-x}} $$ The numerator can be rewritten as a derivative of this geometric series: $$ \sum_{n=0}^\infty (nh\nu) e^{-nx} = -h\nu \frac{d}{dx} \left( \sum_{n=0}^\infty e^{-nx} \right) = -h\nu \frac{d}{dx} \left( \frac{1}{1 - e^{-x}} \right) = -h\nu \left( \frac{-e^{-x}}{(1 - e^{-x})^2} \right) = h\nu \frac{e^{-x}}{(1 - e^{-x})^2} $$ Dividing our formulated numerator by the denominator yields the average energy of a quantized mode: $$ \langle E \rangle = \frac{h\nu \frac{e^{-x}}{(1 - e^{-x})^2}}{\frac{1}{1 - e^{-x}}} = h\nu \frac{e^{-x}}{1 - e^{-x}} = \frac{h\nu}{e^x - 1} = \frac{h\nu}{e^{\frac{h\nu}{k_B T}} - 1} $$ Let us analyze what this equation says:
* At low frequencies ($h\nu \ll k_B T$), we can Taylor-expand the denominator: $e^{\frac{h\nu}{k_B T}} - 1 \approx 1 + \frac{h\nu}{k_B T} - 1 = \frac{h\nu}{k_B T}$. The average energy becomes $\langle E \rangle \approx \frac{h\nu}{h\nu/k_B T} = k_B T$, which is exactly the continuous classical Equipartition result!
* At high frequencies ($h\nu \gg k_B T$), the denominator grows exponentially: $e^{\frac{h\nu}{k_B T}} \to \infty$. This forces the average energy to plummet exponentially to zero: $\langle E \rangle \to 0$. Because high-frequency quantum states require a large packet of energy ($h\nu$) to be excited, the available thermal energy ($k_B T$) is simply too small to excite them, mathematically freezing them out!
Multiplying our density of states by this discrete average energy yields the celebrated **Planck's Radiation Law**: $$ I(\lambda, T) = \frac{8\pi hc}{\lambda^5} \frac{1}{e^{\frac{hc}{\lambda k_B T}} - 1} $$ This beautiful formula perfectly matched the experimental curves, successfully predicting the spectral radiance peak and resolving the Ultraviolet Catastrophe by dropping to zero at short wavelengths: $\lim_{\lambda \to 0} I(\lambda, T) = 0$.

1. Einstein and the Photoelectric Effect (1905): The Collision of Light Quanta

Classical wave electrodynamics predicted that light, as a continuous wave, deposits energy over time. Therefore, the kinetic energy of ejected photoelectric electrons should scale strictly with the **intensity** of the light (which represents the field amplitude squared, $|\mathbf{E}|^2$). Under classical physics, a very dim light should require a long delay (minutes to hours) for an electron to accumulate enough wave energy to break free, whereas extremely bright light should eject incredibly fast electrons.

However, experiments conducted by Philipp Lenard and Robert Millikan revealed three shocking observations that completely shattered classical wave expectations:
1. **Frequency Dependency:** The kinetic energy of the ejected electrons depended strictly on the **frequency** $\nu$ (color) of the incident light, completely independent of its intensity!
2. **Zero Time Delay:** Ejection was strictly instantaneous—even for extremely dim light, electrons were emitted immediately with zero delay.
3. **Threshold Frequency:** Below a specific cutoff frequency $\nu_0$ (the work function threshold), **no** electrons were emitted, regardless of how bright or intense the light source was.
To explain these point-like, all-or-nothing collisions, Albert Einstein extended Planck's quanta hypothesis. He proposed that light propagates not as continuous waves, but as localized packets of energy called **photons**, each carrying a discrete energy: $$ E = h\nu $$ When a photon strikes an electron in a metal, it transfers its entire energy packet in a single, instantaneous collision. A fixed amount of energy, called the **Work Function** $\Phi$, is absorbed to overcome the electrostatic binding forces of the metal lattice. The remaining energy is converted into the electron's maximum kinetic energy.
By applying the principle of conservation of energy to this single-photon collision, we formulate Einstein's celebrated photoelectric equation: $$ K_{\text{max}} = h\nu - \Phi $$ This simple algebraic balance perfectly explained all experimental observations: the threshold frequency is where $h\nu_0 = \Phi$, the speed scales with frequency because energy scales with $\nu$, and the process is instant because it is a point-like collision.

2. Bohr's Quantized Hydrogen Atom (1913): Spectral Lines and Stable Orbits

Ernest Rutherford's alpha-particle scattering experiments proved that atoms consist of lightweight, negative electrons orbiting a dense, positively charged nucleus. However, classical electrodynamics dictated that accelerating charges must continuously emit electromagnetic radiation. An electron orbiting classically would continuously bleed away its kinetic energy, spiraling into the positive nucleus within $10^{-11}$ seconds. Classical physics predicted that all matter in the universe should instantly collapse!
Furthermore, when an electric discharge is passed through hydrogen gas, it does not emit a continuous thermal spectrum. Instead, it emits light at **discrete, sharp spectral lines** (the Balmer series), described by the empirical Rydberg formula: $$ \frac{1}{\lambda} = R \left( \frac{1}{n_{\text{lower}}^2} - \frac{1}{n_{\text{upper}}^2} \right) $$

Niels Bohr resolved both the collapse paradox and the discrete spectrum by postulating that electrons orbit only in non-radiating, stationary states, and that their **angular momentum** $L$ is restricted to integer multiples of $\hbar = h/2\pi$: $$ L = m_e v r = n \hbar, \quad n \in \{1, 2, 3, \dots\} $$ Let us derive how this quantization condition dictates orbital mechanics. Under circular motion, we balance the electrostatic attractive force with the centripetal force: $$ \frac{1}{4\pi\epsilon_0} \frac{e^2}{r^2} = \frac{m_e v^2}{r} \implies m_e v^2 = \frac{e^2}{4\pi\epsilon_0 r} $$ From the angular momentum quantization condition, we solve for velocity: $v = \frac{n\hbar}{m_e r}$. Substituting this into our force balance: $$ m_e \left( \frac{n\hbar}{m_e r} \right)^2 = \frac{e^2}{4\pi\epsilon_0 r} \implies \frac{n^2 \hbar^2}{m_e r^2} = \frac{e^2}{4\pi\epsilon_0 r} $$ Solving for the quantized orbit radius $r_n$: $$ r_n = \frac{4\pi\epsilon_0 \hbar^2}{m_e e^2} n^2 = a_0 n^2 $$ where $a_0 \approx 0.529 \text{ Å}$ is the **Bohr Radius** ($n=1$). This shows that the electron's orbit cannot decay continuously; it can only occupy discrete, stable shells.
The total mechanical energy is the sum of kinetic and electrostatic potential energy: $$ E = K + U = \frac{1}{2}m_e v^2 - \frac{e^2}{4\pi\epsilon_0 r} $$ Substituting our force balance expression $\frac{1}{2}m_e v^2 = \frac{e^2}{8\pi\epsilon_0 r}$: $$ E = \frac{e^2}{8\pi\epsilon_0 r} - \frac{e^2}{4\pi\epsilon_0 r} = -\frac{e^2}{8\pi\epsilon_0 r} $$ Substituting our quantized radius $r_n$ yields the quantized **Bohr Energy Levels**: $$ E_n = -\left(\frac{m_e e^4}{32 \pi^2 \epsilon_0^2 \hbar^2}\right) \frac{1}{n^2} \approx -\frac{13.606 \text{ eV}}{n^2} $$ When an electron transitions from a higher energy orbit $n_{\text{upper}}$ to a lower orbit $n_{\text{lower}}$, it emits a single photon carrying the exact energy difference: $\Delta E = E_{\text{upper}} - E_{\text{lower}} = h\nu = \frac{hc}{\lambda}$.
Substituting our quantized energy states: $$ \frac{1}{\lambda} = \frac{E_{\text{upper}} - E_{\text{lower}}}{hc} = \left( \frac{m_e e^4}{64\pi^3 \epsilon_0^2 \hbar^3 c} \right) \left( \frac{1}{n_{\text{lower}}^2} - \frac{1}{n_{\text{upper}}^2} \right) $$ This derived coefficient is exactly the Rydberg constant $R \approx 1.097 \times 10^7 \text{ m}^{-1}$, mathematically proving Bohr's quantization condition.

3. Wave-Particle Duality and the de Broglie Waves

In 1924, Louis de Broglie proposed a beautiful symmetry: if light waves behave as particles (photons), then matter particles must also possess a wave character. He defined the **matter wavelength** as: $$ \lambda = \frac{h}{p} $$ This wave-like behavior was subsequently confirmed experimentally by Clinton Davisson and Lester Germer, who fired a beam of electrons at a nickel crystal and observed sharp diffraction rings on a detector—a phenomenon unique to wave interference.

4. Postulating the Schrödinger Wave Equation: Deducing Operators from a Plane Wave

If particles behave as waves, what wave equation governs them? Erwin Schrödinger formulated his monumental wave equation using classical energy conservation and a plane wave ansatz.
For a non-relativistic particle, the total energy is the sum of kinetic and potential energy: $$ E = \frac{p^2}{2m} + U(x) $$ To describe this matter wave, Schrödinger assumed a free particle propagating along the x-axis is represented by a plane wave function: $$ \Psi(x, t) = A e^{i(kx - \omega t)} $$ Applying de Broglie's relation $p = \hbar k$ and Planck's relation $E = \hbar\omega$, we rewrite the wave in terms of physical observables: $$ \Psi(x, t) = A e^{\frac{i}{\hbar}(px - Et)} $$ To convert our classical energy equation into a differential wave equation, we find mathematical operators that can extract $E$ and $p^2$ from this wave function:
1. **Energy Operator:** Taking a first derivative of the wave function with respect to time $t$: $$ \frac{\partial\Psi}{\partial t} = -\frac{i E}{\hbar} \Psi \implies i\hbar \frac{\partial\Psi}{\partial t} = E\Psi $$ This maps the physical energy $E$ to the temporal differential operator $\hat{E} = i\hbar\frac{\partial}{\partial t}$.
2. **Momentum Operator:** Taking a first spatial derivative with respect to $x$: $$ \frac{\partial\Psi}{\partial x} = \frac{i p}{\hbar} \Psi \implies -i\hbar \frac{\partial\Psi}{\partial x} = p\Psi $$ This maps momentum $p$ to the spatial operator $\hat{p} = -i\hbar\frac{\partial}{\partial x}$. Differentiating twice with respect to $x$: $$ \frac{\partial^2\Psi}{\partial x^2} = -\frac{p^2}{\hbar^2} \Psi \implies -\hbar^2 \frac{\partial^2\Psi}{\partial x^2} = p^2 \Psi \implies -\frac{\hbar^2}{2m} \frac{\partial^2\Psi}{\partial x^2} = \frac{p^2}{2m} \Psi $$ This maps the kinetic energy term $\frac{p^2}{2m}$ to the spatial operator $-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}$.
Now, let us multiply our classical energy conservation equation $E = \frac{p^2}{2m} + U(x)$ by the wave function $\Psi$: $$ E\Psi = \frac{p^2}{2m}\Psi + U(x)\Psi $$ Substituting our temporal and spatial operators directly into this relation yields the monumental **Time-Dependent Schrödinger Equation**: $$ i\hbar \frac{\partial}{\partial t}\Psi(x,t) = \left[ -\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2} + U(x) \right]\Psi(x,t) $$ This equation replaced classical deterministic trajectories with a **wave function** $\Psi$, whose magnitude squared $|\Psi(\mathbf{r}, t)|^2$ determines the probability density of finding a particle in space.