AP Physics C: Mechanics — Simple Harmonic Motion & Calculus-Driven Oscillations
Target Audience: Caltech-Bound Candidates Aiming for a Score of 5
Institutional Target: Caltech Physics 1a (Ph 1a) Advanced Placement / Exemption Exam
1. Introduction & AP Exam Weight
Simple Harmonic Motion (SHM) and calculus-driven oscillations account for 6%–14% of the AP Physics C: Mechanics exam. While introductory physics treats SHM as a plug-and-chug exercise with sine and cosine equations, the AP Physics C exam—and institutional placement exams at elite universities like Caltech—treats SHM as a gateway to differential equations, Lagrangian mechanics, and Taylor-series approximations of continuous potential energy functions $U(x)$.
To earn a 5 on the AP Physics C exam and position yourself to pass Caltech’s Ph 1a Advanced Placement Examination, you must move beyond memorizing $T = 2\pi\sqrt{\frac{m}{k}}$ and $T = 2\pi\sqrt{\frac{\ell}{g}}$. You must be prepared to: 1. Formulate second-order differential equations from linear forces, torsional forces, and variable mass distributions. 2. Linearize non-linear systems near stable equilibrium points using Taylor expansions ($U''(x_0) > 0$). 3. Integrate continuous density distributions to determine moments of inertia ($I$) and center-of-mass locations ($r_{cm}$) before setting up physical pendulum dynamics. 4. Apply conservation of energy to differential forms $\frac{dE}{dt} = 0$ to solve for oscillatory parameters.
At Caltech, the simple harmonic oscillator is not merely a topic in classical mechanics; it serves as the foundational model for quantum state bound systems (Ph 12a), molecular vibrational modes, and quantum field theory.
2. Deep Concept Breakdown
2.1 The Differential Equation of SHM
The canonical hallmark of Simple Harmonic Motion is a restoring force directly proportional to displacement. Applying Newton's Second Law:
$$F_{net} = -kx \implies m\frac{d^2x}{dt^2} + kx = 0$$
Dividing by $m$ yields the homogeneous second-order linear differential equation:
$$\frac{d^2x}{dt^2} + \omega^2 x = 0 \quad \text{where} \quad \omega = \sqrt{\frac{k}{m}}$$
General Solution Derivation via Complex Exponentials
Assume a trial solution of the form $x(t) = C e^{rt}$. Substituting into the differential equation:
$$r^2 C e^{rt} + \omega^2 C e^{rt} = 0 \implies r^2 + \omega^2 = 0 \implies r = \pm i\omega$$
Thus, the general complex solution is:
$$x(t) = C_1 e^{i\omega t} + C_2 e^{-i\omega t}$$
Using Euler's formula ($e^{i\theta} = \cos\theta + i\sin\theta$), this transforms into real trigonometric form:
$$x(t) = A \cos(\omega t + \phi)$$
Where $A$ is the amplitude and $\phi$ is the phase constant determined by boundary conditions at $t = 0$.
2.2 Oscillations from Arbitrary Potential Energy Functions $U(x)$
For a conservative system, force is related to potential energy by $F(x) = -\frac{dU}{dx}$. An equilibrium point $x_0$ satisfies:
$$\left. \frac{dU}{dx} \right|_{x=x_0} = 0$$
To analyze small displacements $x = x_0 + \eta$ (where $\eta \ll 1$), expand $U(x)$ as a Taylor series around $x_0$:
$$U(x) \approx U(x_0) + \left.\frac{dU}{dx}\right|{x_0}(x - x_0) + \frac{1}{2}\left.\frac{d^2U}{dx^2}\right|{x_0}(x - x_0)^2 + \mathcal{O}((x - x_0)^3)$$
Since $\left.\frac{dU}{dx}\right|_{x_0} = 0$, and defining the potential energy baseline $U(x_0) = 0$:
$$U(\eta) \approx \frac{1}{2} \left( \left.\frac{d^2U}{dx^2}\right|_{x_0} \right) \eta^2$$
Differentiating to find the effective restoring force:
$$F(\eta) = -\frac{dU}{d\eta} = -\left( \left.\frac{d^2U}{dx^2}\right|_{x_0} \right) \eta$$
This matches Hooke's Law $F = -k_{eff} \eta$, where:
$$k_{eff} = \left. \frac{d^2U}{dx^2} \right|_{x = x_0}$$
The angular frequency of small oscillations about any stable equilibrium ($U''(x_0) > 0$) is rigorously given by:
$$\omega = \sqrt{\frac{\left.\frac{d^2U}{dx^2}\right|_{x_0}}{m}}$$
2.3 Physical Pendulums with Variable Density
For a rigid body pivoted about a fixed horizontal axis, torque analysis gives:
$$\Sigma \tau = I \alpha \implies -m g d \sin\theta = I \frac{d^2\theta}{dt^2}$$
where: * $I$ is the moment of inertia about the pivot point (calculated using $I = \int r^2 dm$ or $I = I_{cm} + Mh^2$), * $d$ is the distance from the pivot to the center of mass ($r_{cm}$).
Applying the small-angle approximation $\sin\theta \approx \theta - \frac{\theta^3}{6} \approx \theta$:
$$\frac{d^2\theta}{dt^2} + \left( \frac{m g d}{I} \right) \theta = 0$$
This gives the period of a physical pendulum:
$$\omega = \sqrt{\frac{m g d}{I}} \implies T = 2\pi \sqrt{\frac{I}{m g d}}$$
2.4 Numerical Modeling of Non-Linear Oscillations (Python Implementation)
When displacements $\theta$ are large, the small-angle approximation fails, leading to non-linear behavior:
$$\frac{d^2\theta}{dt^2} + \omega_0^2 \sin\theta = 0$$
The script below uses the numerical Fourth-Order Runge-Kutta (RK4) integration scheme to model exact pendulum dynamics versus the simple harmonic approximation.
import numpy as np
import matplotlib.pyplot as plt
def pendulum_derivs(t, state, omega0_sq):
"""
State vector: state = [theta, omega]
d(theta)/dt = omega
d(omega)/dt = -omega0_sq * sin(theta)
"""
theta, omega = state
dtheta_dt = omega
domega_dt = -omega0_sq * np.sin(theta)
return np.array([dtheta_dt, domega_dt])
def rk4_step(dt, state, omega0_sq):
k1 = pendulum_derivs(0, state, omega0_sq)
k2 = pendulum_derivs(0, state + 0.5 * dt * k1, omega0_sq)
k3 = pendulum_derivs(0, state + 0.5 * dt * k2, omega0_sq)
k4 = pendulum_derivs(0, state + dt * k3, omega0_sq)
return state + (dt / 6.0) * (k1 + 2*k2 + 2*k3 + k4)
# Parameters
g = 9.81
L = 1.0
omega0_sq = g / L
dt = 0.001
t_max = 10.0
time_steps = int(t_max / dt)
# Initial conditions: Large angle displacement (70 degrees)
theta0 = np.radians(70.0)
omega0 = 0.0
state_rk4 = np.array([theta0, omega0])
t_vec = np.linspace(0, t_max, time_steps)
theta_exact = np.zeros(time_steps)
theta_shm = theta0 * np.cos(np.sqrt(omega0_sq) * t_vec) # Small-angle approximation
for i in range(time_steps):
theta_exact[i] = state_rk4[0]
state_rk4 = rk4_step(dt, state_rk4, omega0_sq)
# Plotting phase trajectory comparison
plt.figure(figsize=(10, 5))
plt.plot(t_vec, np.degrees(theta_exact), label="Exact Numerical (RK4: $\\sin\\theta$)", color="crimson")
plt.plot(t_vec, np.degrees(theta_shm), label="Linear SHM Approximated ($\\theta$)", color="navy", linestyle="--")
plt.title("Caltech Baseline: Non-Linear Oscillation vs. Linear SHM Approximation")
plt.xlabel("Time (s)")
plt.ylabel("Angle $\\theta$ (degrees)")
plt.grid(True)
plt.legend()
plt.savefig("pendulum_comparison.png")
plt.show()
3. Common AP Exam Pitfalls & Score 5 Scoring Rubric Nuances
Scoring Breakdown: Score 4 vs. Score 5 Performance
| Analytical Concept | Score 4 Response Pattern | Score 5 Response Pattern |
|---|---|---|
| Derivation of Period | Recalls pre-derived formula $T = 2\pi\sqrt{I/mgd}$ and plugs in numbers directly. | Derives the torque balance equation $\Sigma \tau = I\alpha$, proves harmonic form $\ddot{\theta} + \omega^2\theta = 0$, explicitly identifies $\omega^2$, and states $T = \frac{2\pi}{\omega}$. |
| Arbitrary Potentials | Guesses harmonic motion frequency by setting $U(x) = \frac{1}{2}kx^2$ at arbitrary boundaries. | Computes $x_0$ via $U'(x_0) = 0$, evaluates $k_{eff} = U''(x_0)$, proves stability ($U''(x_0) > 0$), and constructs $\omega = \sqrt{U''(x_0)/m}$. |
| Moment of Inertia Integration | Uses tabular values $I = \frac{1}{3}ML^2$ even when density $\lambda(x)$ is non-uniform. | Sets up spatial integrals $M = \int \lambda(x) dx$, $x_{cm} = \frac{1}{M}\int x \lambda(x) dx$, and $I = \int x^2 \lambda(x) dx$ cleanly from limits. |
| Phase Constant ($\phi$) Formulation | Assumes $\phi = 0$ or $\phi = -\pi/2$ without checking both initial position $x(0)$ and initial velocity $v(0)$. | Uses $\tan\phi = -\frac{v(0)}{\omega x(0)}$ and checks quadrant constraints explicitly to fix $\phi$. |
AP Rubric Points Breakdown (Free-Response Scoring Insights)
- The Differential Equation Point: You must explicitly show the second derivative term ($\frac{d^2x}{dt^2}$ or $\frac{d^2\theta}{dt^2}$) set equal to restoring forces/torques. Simply writing $F = -kx$ is insufficient; you must write $m\frac{d^2x}{dt^2} = -kx$.
- Small Angle Approximation Point: You must explicitly state $\sin\theta \approx \theta$ for $\theta \ll 1$ rad. Omitting this step causes an automatic 1-point deduction on pendulum derivations.
- Parallel Axis Application: When analyzing physical pendulums rotated away from the center of mass, $I = I_{cm} + M d^2$ must be shown algebraically before substituting numerical parameters.
4. Caltech Placement Pathway
Exemption Criteria: Waiving Physics 1a (Ph 1a)
Caltech requires all incoming undergraduates to take or place out of Ph 1a (Classical Mechanics). To earn a waiver through the placement exam administered during Orientation Week: * A score of 5 on AP Physics C: Mechanics is a prerequisite to sit for the advanced standing exam. * The exemption exam tests classical mechanics at the level of Analytical Mechanics (Fowles & Cassiday or Morin).
Transition to Ph 1b / Ph 12a
By passing the Ph 1a Exemption Exam, high-performing students bypass introductory kinematics and accelerate directly into: * Ph 1b (Electromagnetism): Utilizes dynamic field equations and driven harmonic oscillators under Lorentz forces. * Ph 12a (Quantum Mechanics): Covers the Quantum Harmonic Oscillator via ladder operators ($a, a^\dagger$):
$$\hat{H} = \hbar \omega \left( a^\dagger a + \frac{1}{2} \right)$$
Mastering Taylor expansions of classical potentials $U(x)$ directly translates to constructing quantum potential wells and computing zero-point energies.
5. High-Yield Practice Problem
Problem Statement
A non-uniform, thin rigid rod of length $L$ and total mass $M$ has a variable linear mass density given by:
$$\lambda(x) = \lambda_0 \left(1 + \frac{x}{L}\right)$$
where $x$ is measured from an axle located at the end $x = 0$. The rod is free to rotate without friction in a vertical plane about a horizontal axis passing through $x = 0$.
Pivot (Axle)
O=====================================================> x
| Rod Length L
(x=0) (x=L)
|
v Gravity (g)
- Express the constant $\lambda_0$ in terms of $M$ and $L$.
- Determine the distance $d$ from the pivot at $x = 0$ to the center of mass of the rod.
- Calculate the moment of inertia $I$ of the rod about the pivot at $x = 0$.
- Derive the second-order differential equation governing small angular displacements $\theta$ from vertical equilibrium.
- Derive an expression for the period of small oscillations $T$.
- At $t = 0$, the rod is displaced by an angle $\theta_0$ ($\theta_0 \ll 1$ rad) and released with an initial angular velocity $\Omega_0$ directed toward the equilibrium position. Determine the complete equation for angular position $\theta(t)$.
Complete Step-by-Step Solution Checklist
Step 1: Calculate Mass Density Constant $\lambda_0$
The total mass $M$ is the integral of the mass element $dm = \lambda(x) dx$:
$$M = \int_0^L \lambda(x) \, dx = \int_0^L \lambda_0 \left(1 + \frac{x}{L}\right) dx$$
$$M = \lambda_0 \left[ x + \frac{x^2}{2L} \right]_0^L = \lambda_0 \left( L + \frac{L}{2} \right) = \frac{3}{2} \lambda_0 L$$
$$\implies \lambda_0 = \frac{2M}{3L}$$
Step 2: Determine Distance to Center of Mass $d$
The center of mass position $x_{cm} = d$ is:
$$d = \frac{1}{M} \int_0^L x \, dm = \frac{1}{M} \int_0^L x \cdot \lambda_0 \left(1 + \frac{x}{L}\right) dx$$
Substitute $\lambda_0 = \frac{2M}{3L}$:
$$d = \frac{1}{M} \left( \frac{2M}{3L} \right) \int_0^L \left( x + \frac{x^2}{L} \right) dx$$
$$d = \frac{2}{3L} \left[ \frac{x^2}{2} + \frac{x^3}{3L} \right]_0^L = \frac{2}{3L} \left( \frac{L^2}{2} + \frac{L^2}{3} \right) = \frac{2}{3L} \left( \frac{5L^2}{6} \right) = \frac{5}{9}L$$
Step 3: Compute Moment of Inertia $I$ about the Pivot
$$I = \int_0^L x^2 \, dm = \int_0^L x^2 \lambda_0 \left(1 + \frac{x}{L}\right) dx$$
$$I = \left(\frac{2M}{3L}\right) \int_0^L \left(x^2 + \frac{x^3}{L}\right) dx$$
$$I = \left(\frac{2M}{3L}\right) \left[ \frac{x^3}{3} + \frac{x^4}{4L} \right]_0^L = \left(\frac{2M}{3L}\right) \left( \frac{L^3}{3} + \frac{L^3}{4} \right)$$
$$I = \left(\frac{2M}{3L}\right) \left( \frac{7L^3}{12} \right) = \frac{7}{18} M L^2$$
Step 4: Derive Differential Equation for Oscillations
Apply Newton's Second Law for Rotation:
$$\Sigma \tau_{pivot} = I \alpha = I \frac{d^2\theta}{dt^2}$$
The restoring torque exerted by gravity acts through the center of mass at distance $d$:
$$\tau_{restoring} = - M g d \sin\theta$$
Equating torque to angular acceleration:
$$- M g d \sin\theta = I \frac{d^2\theta}{dt^2} \implies \frac{d^2\theta}{dt^2} + \frac{M g d}{I} \sin\theta = 0$$
Apply the small-angle approximation ($\sin\theta \approx \theta$ for $\theta \ll 1$):
$$\frac{d^2\theta}{dt^2} + \left(\frac{M g d}{I}\right) \theta = 0$$
Step 5: Calculate Period of Oscillation $T$
From the differential equation, identify $\omega^2$:
$$\omega^2 = \frac{M g d}{I}$$
Substitute $d = \frac{5}{9}L$ and $I = \frac{7}{18} M L^2$:
$$\omega^2 = \frac{M g \left(\frac{5}{9}L\right)}{\frac{7}{18} M L^2} = \frac{\frac{5}{9}}{\frac{7}{18}} \cdot \frac{g}{L} = \left(\frac{5}{9} \cdot \frac{18}{7}\right) \frac{g}{L} = \frac{10}{7} \frac{g}{L}$$
$$\omega = \sqrt{\frac{10 g}{7 L}}$$
The period $T$ is:
$$T = \frac{2\pi}{\omega} = 2\pi \sqrt{\frac{7 L}{10 g}}$$
Step 6: Construct Complete Solution $\theta(t)$ with Phase Boundary Conditions
The general harmonic solution is:
$$\theta(t) = A \cos(\omega t + \phi)$$
Taking the time derivative to get angular velocity $\omega_{ang}(t)$:
$$\omega_{ang}(t) = \frac{d\theta}{dt} = -A \omega \sin(\omega t + \phi)$$
Apply initial conditions at $t = 0$: 1. Position: $\theta(0) = \theta_0 \implies A \cos\phi = \theta_0$ 2. Velocity: $\omega_{ang}(0) = -\Omega_0$ (negative because initial motion is toward equilibrium $x=0$):
$$-A \omega \sin\phi = -\Omega_0 \implies A \sin\phi = \frac{\Omega_0}{\omega}$$
Squaring and adding both initial condition equations ($A^2\cos^2\phi + A^2\sin^2\phi = A^2$):
$$A = \sqrt{\theta_0^2 + \left(\frac{\Omega_0}{\omega}\right)^2}$$
Dividing the sine equation by the cosine equation:
$$\tan\phi = \frac{\Omega_0}{\omega \theta_0} \implies \phi = \arctan\left(\frac{\Omega_0}{\omega \theta_0}\right)$$
Substituting $\omega = \sqrt{\frac{10g}{7L}}$ into the general expression:
$$\theta(t) = \sqrt{\theta_0^2 + \frac{7L\Omega_0^2}{10g}} \, \cos \left( \sqrt{\frac{10g}{7L}} \, t + \arctan\left( \Omega_0 \sqrt{\frac{7L}{10g \, \theta_0^2}} \right) \right)$$
Master Review Verification Criteria
- Units Check: $\left[\frac{g}{L}\right] = \frac{\text{m/s}^2}{\text{m}} = \text{s}^{-2}$, giving $\omega$ in rad/s. Correct.
- Mass Invariance: $\omega$ is independent of total mass $M$, as expected for linear gravitational restoring torques.
- Limit Check: As density becomes uniform ($\lambda(x) = \text{constant}$), $d \to \frac{1}{2}L$ and $I \to \frac{1}{3}ML^2$, yielding $\omega^2 \to \frac{M g (L/2)}{(1/3)M L^2} = \frac{3g}{2L}$. In our result, setting spatial weighting density uniformly reproduces classic rigid rod dynamics.