Generalized Coordinates and Configuration Space

Ask Newton where a double pendulum is, and he answers with six Cartesian numbers — two masses, three coordinates each — tied together by two constraint equations and policed by rod tensions he must compute along the way. Ask Lagrange, and he answers with two angles. Nothing else. The rods, the tensions, the bookkeeping — all of it vanishes into the choice of description, and what remains is the honest content of the problem: two numbers that can each run around a circle, so that the set of every shape the pendulum can take is a torus. The motion of the pendulum — all of its swings, flips, and chaos — is a single curve winding over that doughnut's surface. That is the move this article teaches, and it is the founding move of analytical mechanics: stop describing systems by where their atoms sit in space, and start describing them by any set of independent numbers that pins down the configuration — angles, arc lengths, ratios, whatever fits. These are generalized coordinates, and the space they chart, one point per possible shape of the system, is configuration space. The payoff is structural, not cosmetic. Constraints stop being equations you enforce and become geometry you inhabit: a bead on a wire is a one-dimensional world, a rigid body's orientations are the rotation group SO(3)SO(3). Mechanics becomes the study of curves on manifolds — a viewpoint so durable it survived intact into quantum mechanics (path integrals sum over curves in configuration space) and now runs the motion planners of every robot arm shipping in 2026. This is the first brick of the Lagrangian edifice; lay it carefully.

Key takeaways
  • Generalized coordinates q1,…,qnq_1,\dots,q_n are any independent parameters that uniquely fix a system's configuration — angles, lengths, ratios; nothing Cartesian is required.
  • For NN particles under kk independent holonomic constraints, the number needed is n=3N−kn = 3N - k, the number of degrees of freedom.
  • The set of all attainable configurations is an nn-dimensional manifold Q\mathcal{Q}, the configuration space: a circle S1S^1 for a pendulum, a torus T2T^2 for a double pendulum, SO(3)SO(3) for a rigid body's orientation.
  • The embedding ra=ra(q1,…,qn,t)\mathbf{r}_a = \mathbf{r}_a(q_1,\dots,q_n,t) carries the abstract description back to physical space; the chain rule gives r˙a\dot{\mathbf{r}}_a linear in the q˙i\dot q_i, and the "cancellation of dots" ∂r˙a/∂q˙i=∂ra/∂qi\partial\dot{\mathbf{r}}_a/\partial\dot q_i = \partial\mathbf{r}_a/\partial q_i becomes the workhorse identity of everything downstream.
  • Kinetic energy in generalized coordinates is T=12Mij(q) q˙iq˙jT = \tfrac12 M_{ij}(q)\,\dot q_i\dot q_j for time-independent constraints — the mass matrix MijM_{ij} is a metric on Q\mathcal{Q}, and free motions are its geodesics.
  • A point of Q\mathcal{Q} is a photograph, not a movie: dynamics needs (q,q˙)(q, \dot q), which live in the 2n2n-dimensional tangent bundle TQT\mathcal{Q}, the home of the Lagrangian.
  • Key idea: pick n=3N−kn = 3N - k independent parameters that build the constraints into the description from the start; the system then lives on an nn-dimensional manifold, and dynamics reduces to free motion of these parameters.

Where this lives today

Configuration space stopped being an abstraction the day machines had to plan their own motions. In robotics the term of art is literally C-space: a seven-joint manipulator arm is a point moving on a seven-dimensional torus-like manifold, workspace obstacles map to forbidden C-space regions, and the planning algorithms that ship in 2026-era industrial arms and humanoids — RRTs, probabilistic roadmaps, trajectory optimizers — all answer one question: find a curve through the allowed region of configuration space. Grasping-and-manipulation research spends most of its effort on the geometry and topology of these spaces before any dynamics enters; gimbal-lock avoidance in the wrist is exactly the chart-singularity problem of SO(3)SO(3) from this article, solved in production code by quaternion charts. The same picture organizes physics at the molecular scale. A protein of NN atoms lives in a 3N3N-dimensional configuration space, and folding is the search for low-energy basins in that space; the energy landscapes community and the AlphaFold generation of structure predictors both work, in different languages, on the geography of Q\mathcal{Q}. And when mechanics went quantum, configuration space kept its job: Feynman's path integral assigns amplitudes to curves in configuration space, and the wavefunction of a many-body system is a function on Q\mathcal{Q}, not on physical space — a fact responsible for most of quantum theory's strangeness and all of quantum chemistry's computational cost.

The physical picture

Thread a bead onto a bent wire and ask: how much information does it take to say where the bead is? Three Cartesian coordinates is the wrong answer — they are shackled by two constraint equations (the bead cannot leave the wire) and carry redundant information. The right answer is one number: the arc length ss measured along the wire. Give me ss and I know everything; change ss and the bead moves in the only way it can. The wire's bends and kinks, which would torture the Cartesian description with position-dependent constraint forces, are silently absorbed into the meaning of ss. That one number is a generalized coordinate, and the wire itself — regarded abstractly as a one-dimensional line the bead lives on — is the configuration space. Now scale the idea up. A planar pendulum needs one angle θ\theta, and since θ\theta and θ+2π\theta + 2\pi are the same configuration, its configuration space is a circle S1S^1. A double pendulum needs two angles, each periodic, so its configuration space is the product of two circles — a torus T2T^2, an honest doughnut. Every possible contortion of the double pendulum is one point on that doughnut; every actual motion is a curve winding across it.

Figure
The double pendulum and its configuration space — every shape is one point on a torus, every motion a curve on it
physical space: two angles θ₁ θ₂ six Cartesian numbers, two constraints, unknown rod tensions — or just (θ₁, θ₂) configuration space: a torus θ₁ around θ₂ around one motion = one curve on the torus every possible contortion = one point

Three habits of mind make the machinery feel natural. First, a configuration is a photograph, not a movie: a point of configuration space records where everything is, and says nothing about where anything is going. Two identical photographs of a pendulum at its lowest point can hide wildly different fates — hanging at rest, or blasting through at ten radians per second. The dynamical state needs the pair (q,q˙)(q, \dot q), and that pair lives in a bigger space (the tangent bundle) that we will meet below. Second, coordinates are labels, not physics. You may chart the same configuration space by angles, by arc lengths, by any invertible relabeling, and no experiment can tell the difference; the art — most of the craft of problem-solving in this subject — is picking labels that make the equations pretty. Nature does not care that you measured θ\theta from the downward vertical. Third, the space itself carries physics even before any equation of motion is written. That the double pendulum's space is a compact torus rather than an infinite plane already shapes what motions are possible; that a rigid body's orientation space SO(3)SO(3) cannot be covered by one smooth chart is precisely why every flight controller must manage gimbal lock or switch to quaternions. Topology is destiny.

Figure
A field guide to configuration spaces — the same counting, radically different geographies
pendulum Q = S¹ one angle, periodic n = 1 bead on wire Q = the wire arc length s n = 1 rigid body Q = R³ × SO(3) place + orientation n = 6, any size of body the count n says how many numbers; the topology says how they wrap — both matter

The core equations

Everything begins with a declaration of independence. A system's configuration shall be specified by a tuple of independent parameters, whose number is fixed by the counting rule made precise in the degrees-of-freedom article:

q=(q1,q2,…,qn),n=3N−k(1)q = (q_1, q_2, \dots, q_n), \qquad n = 3N - k \qquad (1)

where: qiq_i — the ii-th generalized coordinate (an angle, a length, a ratio — anything); NN — the number of point particles; kk — the number of independent holonomic constraints, each one scalar equation fj(r1,…,rN,t)=0f_j(\mathbf{r}_1,\dots,\mathbf{r}_N,t) = 0 of the kind classified in constraints and their classification. Independent means no qiq_i is a function of the others; complete means the tuple pins the configuration uniquely. The count nn is intrinsic to the system; the choice of which nn parameters is gloriously free. The abstract labels must still connect to physical space — forces act on particles, not on angles. The bridge is the embedding map, which every problem's setup phase constructs explicitly:

ra=ra(q1,…,qn,t),a=1,…,N(2)\mathbf{r}_a = \mathbf{r}_a(q_1, \dots, q_n, t), \qquad a = 1, \dots, N \qquad (2)

where: ra\mathbf{r}_a — the Cartesian position of particle aa, written as a function of the generalized coordinates. The explicit tt appears only when the constraints themselves move — a bead on a rotating wire, a pendulum whose pivot is being shaken (the rheonomic case). For the double pendulum the map reads r1=(ℓ1sin⁡θ1,−ℓ1cos⁡θ1)\mathbf{r}_1 = (\ell_1\sin\theta_1, -\ell_1\cos\theta_1), r2=r1+(ℓ2sin⁡θ2,−ℓ2cos⁡θ2)\mathbf{r}_2 = \mathbf{r}_1 + (\ell_2\sin\theta_2, -\ell_2\cos\theta_2): six Cartesian numbers manufactured from two angles, constraints satisfied automatically, forever, by construction.

Derivation AThe velocity expansion and the cancellation of dots

We are deriving the two identities on which the entire Lagrangian formalism will run: the expansion of a particle's true velocity r˙a\dot{\mathbf{r}}_a in terms of generalized velocities q˙i\dot q_i, and the innocent-looking "cancellation of dots" ∂r˙a/∂q˙i=∂ra/∂qi\partial\dot{\mathbf{r}}_a/\partial\dot q_i = \partial\mathbf{r}_a/\partial q_i. The strategy is nothing but the multivariable chain rule applied to the embedding map (2), followed by careful bookkeeping of what depends on what. The point to watch is exactly that bookkeeping: on the tangent bundle, qiq_i and q˙i\dot q_i are treated as independent variables — a convention that confuses everyone once, and which this derivation makes precise rather than mysterious.

Step 1. Differentiate the embedding along a motion. On an actual trajectory the coordinates are functions of time, ra=ra(q(t),t)\mathbf{r}_a = \mathbf{r}_a(q(t), t), so the chain rule gives

r˙a=∑i=1n∂ra∂qi q˙i+∂ra∂t(3)\dot{\mathbf{r}}_a = \sum_{i=1}^{n} \frac{\partial \mathbf{r}_a}{\partial q_i}\,\dot{q}_i + \frac{\partial \mathbf{r}_a}{\partial t} \qquad (3)

where: q˙i=dqi/dt\dot q_i = \mathrm{d}q_i/\mathrm{d}t — the generalized velocities; ∂ra/∂qi\partial\mathbf{r}_a/\partial q_i — the tangent vectors, pointing along the direction particle aa moves when coordinate qiq_i alone is nudged; ∂ra/∂t\partial\mathbf{r}_a/\partial t — the drift velocity imposed by moving constraints, zero in the scleronomic case. Why this is more than calculus: the nn vectors ∂ra/∂qi\partial\mathbf{r}_a/\partial q_i span, at each configuration, the directions of allowed motion — they are a moving basis for the constraint surface, and every legal velocity is a linear combination of them plus the drift.

Step 2. Read off the structure. The generalized velocities enter (3) linearly, with coefficients ∂ra/∂qi\partial\mathbf{r}_a/\partial q_i that depend only on qq and tt — not on any velocity. The velocity r˙a\dot{\mathbf{r}}_a is thus a fixed, configuration-dependent linear machine fed by the q˙i\dot q_i.

Step 3. Differentiate with respect to a generalized velocity, treating (q,q˙)(q, \dot q) as independent arguments:

∂r˙a∂q˙j=∂∂q˙j(∑i∂ra∂qi q˙i+∂ra∂t)=∂ra∂qj(4)\frac{\partial \dot{\mathbf{r}}_a}{\partial \dot{q}_j} = \frac{\partial}{\partial \dot q_j}\left(\sum_i \frac{\partial\mathbf{r}_a}{\partial q_i}\,\dot q_i + \frac{\partial\mathbf{r}_a}{\partial t}\right) = \frac{\partial \mathbf{r}_a}{\partial q_j} \qquad (4)

where: only the i=ji = j term survives, because the coefficients and the drift are q˙\dot q-independent. Why treating qq and q˙\dot q as independent is legitimate: the derivative asks how r˙a\dot{\mathbf{r}}_a responds when the system passes through the SAME configuration with a slightly different velocity — comparing neighboring states, not neighboring instants of one motion. States with equal qq and unequal q˙\dot q are genuinely different points of the state space, so the partial derivative is well defined.

Step 4. Record the companion identity for later use. Differentiating (3) with respect to qjq_j instead, and comparing with the total time derivative of ∂ra/∂qj\partial\mathbf{r}_a/\partial q_j, shows that ddt∂ra∂qj=∂r˙a∂qj\frac{\mathrm{d}}{\mathrm{d}t}\frac{\partial\mathbf{r}_a}{\partial q_j} = \frac{\partial\dot{\mathbf{r}}_a}{\partial q_j} — time differentiation and qq-differentiation commute on the embedding. Why both identities matter: the pair (4) and this commutation rule are precisely the two moves that will turn Newton's equations into the Euler–Lagrange equations one article downstream — they are the entire technical content of that future derivation, banked here. *Verified: Goldstein §1.4 and Landau–Lifshitz §1–2 state both identities in this form; check on the planar pendulum, r=(ℓsin⁡θ,−ℓcos⁡θ)\mathbf{r} = (\ell\sin\theta, -\ell\cos\theta): eq. (3) gives r˙=ℓθ˙(cos⁡θ,sin⁡θ)\dot{\mathbf{r}} = \ell\dot\theta(\cos\theta, \sin\theta), and indeed ∂r˙/∂θ˙=ℓ(cos⁡θ,sin⁡θ)=∂r/∂θ\partial\dot{\mathbf{r}}/\partial\dot\theta = \ell(\cos\theta,\sin\theta) = \partial\mathbf{r}/\partial\theta.*

One more object completes the kinematic toolkit: the virtual displacement, the difference between neighboring configurations at frozen time,

δra=∑i=1n∂ra∂qi δqi(5)\delta \mathbf{r}_a = \sum_{i=1}^{n} \frac{\partial \mathbf{r}_a}{\partial q_i}\, \delta q_i \qquad (5)

where: δqi\delta q_i — independent, arbitrary variations of the coordinates; the explicit-time drift of (3) is absent because δt=0\delta t = 0 by decree. Freezing time looks like a technicality and is anything but: for moving constraints the real displacement in time dt\mathrm{d}t includes the constraint's own motion, while the virtual one probes only the directions the constraint currently permits. This distinction is the load-bearing wall of the principle of virtual work, where these δqi\delta q_i become the test probes of statics; here we need only the definition.

Which formula answers which question

you wantusenotes
how many coordinates a system needsn=3N−kn = 3N - k, eq. (1)holonomic constraints only; details in the DOF article
to set up a problem in your chosen qiq_ithe embedding ra(q,t)\mathbf{r}_a(q, t), eq. (2)write it once, differentiate forever
particle velocities from generalized oneschain rule, eq. (3)drift term only for moving constraints
the workhorse identity for Lagrangian proofscancellation of dots, eq. (4)plus the commutation rule of Step 4
constraint-respecting test displacementsδra\delta\mathbf{r}_a, eq. (5)frozen time; feeds virtual work
kinetic energy in the qiq_iT=12Mijq˙iq˙jT = \tfrac12 M_{ij}\dot q_i\dot q_j, eq. (7)scleronomic; MijM_{ij} from eq. (8)

The geometry: configuration space and its tangent bundle

The set of every configuration the constraints allow is not just a list — it has a shape. It is a smooth nn-dimensional manifold, the configuration space Q\mathcal{Q}: locally chartable by nn real numbers (that is what a coordinate system is), globally as curved or as wrapped as the constraints dictate. A motion of the system is a curve q(t)q(t) on Q\mathcal{Q}; solving a mechanics problem means finding which curves nature selects. Positions alone, however, do not determine futures. To launch a motion you must supply a configuration and a velocity, and the velocity q˙\dot q at a point qq is a vector in the tangent space at that point. The full state space is therefore the collection of all (point, attached velocity) pairs — the tangent bundle:

dim⁡Q=n,dim⁡TQ=2n(6)\dim \mathcal{Q} = n, \qquad \dim T\mathcal{Q} = 2n \qquad (6)

where: TQT\mathcal{Q} — the tangent bundle, coordinatized by (q1,…,qn,q˙1,…,q˙n)(q_1,\dots,q_n,\dot q_1,\dots,\dot q_n). This is the natural home of the Lagrangian function L(q,q˙,t)L(q, \dot q, t), and the doubling from nn to 2n2n is why the equations of motion will be second-order in time: a point of Q\mathcal{Q} underdetermines the future exactly by the missing nn velocities. (One more doubling awaits much later: replacing velocities by momenta trades TQT\mathcal{Q} for the phase space T∗QT^*\mathcal{Q} of Hamiltonian mechanics — same dimension, profoundly different geometry.)

Figure
Configuration space versus state space — a photograph underdetermines the movie
configuration space Q (dim n) the point q — a photograph q̇ fast q̇ slow same point, different attached velocities → different futures. state = (q, q̇) lives in the tangent bundle TQ, dimension 2n — the Lagrangian's home ground

Distance on Q\mathcal{Q} is not measured with a ruler — it is measured with mass. The kinetic energy, rewritten in generalized coordinates, hands the manifold a natural metric, and that rewriting is short enough to do in full:

Go deeper · Derivation B · 2 minKinetic energy as a quadratic form and the mass metric
Derivation BKinetic energy as a quadratic form and the mass metric

The goal is the master formula for kinetic energy in generalized coordinates: for time-independent constraints, TT is a homogeneous quadratic form 12Mij(q) q˙iq˙j\tfrac12 M_{ij}(q)\,\dot q_i \dot q_j, whose coefficient matrix — the mass matrix — is symmetric, positive-definite, and configuration-dependent. The strategy is direct substitution of the velocity expansion (3) into the Cartesian kinetic energy, then reading the resulting structure with a geometer's eye. Watch two things: where the scleronomic assumption enters (it kills the drift term, and with it any pieces linear in or independent of q˙\dot q), and the fact that MijM_{ij} depends on qq but never on q˙\dot q — that is what "quadratic form" means and what makes the geometry clean.

Step 1. Start from the Cartesian kinetic energy of NN particles, T=12∑ama r˙a⋅r˙aT = \tfrac12\sum_a m_a\,\dot{\mathbf{r}}_a\cdot\dot{\mathbf{r}}_a. Why this is the honest starting point: whatever exotic coordinates we adopt, the physical kinetic energy is defined in the inertial Cartesian frame; the generalized-coordinate expression must be derived, not decreed.

Step 2. Substitute the velocity expansion. With scleronomic constraints, ∂ra/∂t=0\partial\mathbf{r}_a/\partial t = 0, so (3) collapses to r˙a=∑i(∂ra/∂qi) q˙i\dot{\mathbf{r}}_a = \sum_i (\partial\mathbf{r}_a/\partial q_i)\,\dot q_i — homogeneous linear in the generalized velocities. Why the drift's absence matters: keep it, and squaring would generate terms linear in q˙\dot q and a term with no velocities at all — the rheonomic complications the second worked example exhibits concretely.

Step 3. Expand the dot product and swap summation order:

T=12∑i,j=1n(∑ama ∂ra∂qi⋅∂ra∂qj)q˙i q˙jT = \tfrac{1}{2}\sum_{i,j=1}^{n}\left(\sum_a m_a\,\frac{\partial\mathbf{r}_a}{\partial q_i}\cdot\frac{\partial\mathbf{r}_a}{\partial q_j}\right) \dot q_i\,\dot q_j

Step 4. Name the object in parentheses:

T=12∑i,j=1nMij(q) q˙iq˙j,Mij(q)=∑ama ∂ra∂qi⋅∂ra∂qj(7), (8)T = \tfrac{1}{2}\sum_{i,j=1}^{n} M_{ij}(q)\,\dot q_i \dot q_j, \qquad M_{ij}(q) = \sum_{a} m_a\,\frac{\partial\mathbf{r}_a}{\partial q_i}\cdot\frac{\partial\mathbf{r}_a}{\partial q_j} \qquad (7),\ (8)

where: MijM_{ij} — the mass matrix (or mass metric): symmetric by inspection, positive-definite because T>0T > 0 for any nonzero velocity, and a function of configuration alone. Why "metric" is the right word: 2T dt2=Mij dqi dqj2T\,\mathrm{d}t^2 = M_{ij}\,\mathrm{d}q_i\,\mathrm{d}q_j assigns a squared length to every small step in configuration space, weighted by how much mass each direction drags. With no forces acting, the system coasts along geodesics of this metric — straightest-possible lines in the mass geometry — which is the earliest glimpse of mechanics becoming geometry.

Step 5. Sanity-check on one particle in polar coordinates. The embedding r=(rcos⁡θ,rsin⁡θ)\mathbf{r} = (r\cos\theta, r\sin\theta) gives tangent vectors ∂r/∂r=(cos⁡θ,sin⁡θ)\partial\mathbf{r}/\partial r = (\cos\theta, \sin\theta) and ∂r/∂θ=(−rsin⁡θ,rcos⁡θ)\partial\mathbf{r}/\partial\theta = (-r\sin\theta, r\cos\theta); their dot products give Mrr=mM_{rr} = m, Mθθ=mr2M_{\theta\theta} = mr^2, Mrθ=0M_{r\theta} = 0, so T=12m(r˙2+r2θ˙2)T = \tfrac12 m(\dot r^2 + r^2\dot\theta^2) — the formula every textbook quotes, now manufactured rather than memorized. Verified: Landau & Lifshitz §5 and Goldstein §1.4 give (7)–(8) in identical form; the polar-coordinate special case and the positive-definiteness argument match Arnold Ch. 4, where the geodesic reading is developed.

Go deeper · Derivation C · 2 minThe mass matrix of the double pendulum
Derivation CThe mass matrix of the double pendulum

The goal is to see the machinery of the previous derivation produce something with real teeth: the mass matrix of the planar double pendulum (masses m1,m2m_1, m_2, rod lengths ℓ1,ℓ2\ell_1, \ell_2, angles θ1,θ2\theta_1, \theta_2 from the vertical), whose off-diagonal entry — the dynamical coupling between the two arms — is the seed of the system's celebrated chaos. The strategy is the assembly line just built: write the embedding, differentiate for the tangent vectors, take dot products weighted by mass. The step to watch is the cross term: it appears because both angles move the second mass, and its size depends on the relative angle θ1−θ2\theta_1 - \theta_2, which is what makes the matrix genuinely configuration-dependent rather than constant.

Step 1. Write the embedding, measuring yy downward for convenience: r1=ℓ1(sin⁡θ1,cos⁡θ1)\mathbf{r}_1 = \ell_1(\sin\theta_1, \cos\theta_1) and r2=r1+ℓ2(sin⁡θ2,cos⁡θ2)\mathbf{r}_2 = \mathbf{r}_1 + \ell_2(\sin\theta_2, \cos\theta_2). Why mass 2 is the interesting one: it is carried by both coordinates, so both tangent vectors at r2\mathbf{r}_2 are nonzero — the anatomical source of the coupling.

Step 2. Differentiate. For mass 1: ∂r1/∂θ1=ℓ1(cos⁡θ1,−sin⁡θ1)\partial\mathbf{r}_1/\partial\theta_1 = \ell_1(\cos\theta_1, -\sin\theta_1), ∂r1/∂θ2=0\partial\mathbf{r}_1/\partial\theta_2 = 0. For mass 2: ∂r2/∂θ1=ℓ1(cos⁡θ1,−sin⁡θ1)\partial\mathbf{r}_2/\partial\theta_1 = \ell_1(\cos\theta_1, -\sin\theta_1), ∂r2/∂θ2=ℓ2(cos⁡θ2,−sin⁡θ2)\partial\mathbf{r}_2/\partial\theta_2 = \ell_2(\cos\theta_2, -\sin\theta_2).

Step 3. Assemble MijM_{ij} from (8). The diagonal entries collect every mass moved by each angle:

M11=(m1+m2) ℓ12,M22=m2 ℓ22M_{11} = (m_1 + m_2)\,\ell_1^2, \qquad M_{22} = m_2\,\ell_2^2

Step 4. Compute the off-diagonal entry — only mass 2 contributes, and the dot product of the two unit tangent directions is cos⁡θ1cos⁡θ2+sin⁡θ1sin⁡θ2=cos⁡(θ1−θ2)\cos\theta_1\cos\theta_2 + \sin\theta_1\sin\theta_2 = \cos(\theta_1 - \theta_2):

M12=M21=m2 ℓ1ℓ2 cos⁡(θ1−θ2)(9)M_{12} = M_{21} = m_2\,\ell_1 \ell_2\,\cos(\theta_1 - \theta_2) \qquad (9)

where: the coupling is strongest when the arms are parallel (θ1=θ2\theta_1 = \theta_2, straight or folded) and vanishes when they are perpendicular — momentarily, the two angles then move mass 2 in orthogonal directions and cannot trade kinetic energy. Why this matters downstream: a configuration-dependent, non-diagonal mass matrix is precisely what no rescaling of coordinates can globally flatten, and the resulting velocity-dependent couplings are the mechanical origin of the double pendulum's chaotic dynamics.

Step 5. Sanity checks. Set m2→0m_2 \to 0: the matrix collapses to M11=m1ℓ12M_{11} = m_1\ell_1^2, a lone pendulum, as it must. Freeze θ1\theta_1: the remaining T=12m2ℓ22θ˙22T = \tfrac12 m_2\ell_2^2\dot\theta_2^2 is a simple pendulum hung from a fixed point. Positive-definiteness: det⁡M=m2ℓ12ℓ22 [m1+m2sin⁡2(θ1−θ2)]>0\det M = m_2\ell_1^2\ell_2^2\,[m_1 + m_2\sin^2(\theta_1-\theta_2)] > 0 always. Verified: matches the double-pendulum kinetic energy in Goldstein Ch. 1 problems and Morin Ch. 6; the determinant is positive for all angles, confirming positive-definiteness, and both limiting checks reduce to the correct single-pendulum results.

Limiting cases

CaseResultMeaning
no constraintsn=3Nn = 3N, Q=R3N\mathcal{Q} = \mathbb{R}^{3N}Cartesian coordinates are already generalized coordinates; nothing to gain
one particle on a surfacen=2n = 2the surface itself is the configuration space; its curvature enters MijM_{ij}
rigid body (NN large, distances all fixed)n=6n = 6, Q=R3×SO(3)\mathcal{Q} = \mathbb{R}^3 \times SO(3)∼1023\sim 10^{23} particles, six numbers — the greatest data compression in physics
scleronomic constraintsT=12Mijq˙iq˙jT = \tfrac12 M_{ij}\dot q_i\dot q_j, purely quadraticenergy bookkeeping stays simple; E=T+VE = T + V downstream
rheonomic constraintsTT gains linear and q˙\dot q-free termsmoving constraints pump energy; E≠T+VE \ne T + V in general
nonholonomic constraintscoordinates cannot be eliminatedvelocity restrictions do not shrink Q\mathcal{Q}; see the constraints article

Worked examples

1 · Charting the spherical pendulum (confidence builder)

A mass m=0.20 kgm = 0.20\ \mathrm{kg} hangs from a rigid massless rod of length ℓ=0.50 m\ell = 0.50\ \mathrm{m}, free to swing in any direction. Count its degrees of freedom, name its configuration space, and compute its kinetic energy at the instant when the rod makes θ=60°\theta = 60° with the vertical, with θ˙=1.2 rad/s\dot\theta = 1.2\ \mathrm{rad/s} and azimuthal rate ϕ˙=2.0 rad/s\dot\phi = 2.0\ \mathrm{rad/s}.

Step 1. Count. One particle, 3N=33N = 3; the rigid rod fixes one distance, k=1k = 1; so n=3−1=2n = 3 - 1 = 2. Natural coordinates: polar angle θ\theta from the downward vertical and azimuth ϕ\phi. The configuration space is the sphere S2S^2 of radius ℓ\ell — every point the bob can occupy.

Step 2. Build the mass matrix from the embedding r=ℓ(sin⁡θcos⁡ϕ, sin⁡θsin⁡ϕ, −cos⁡θ)\mathbf{r} = \ell(\sin\theta\cos\phi,\ \sin\theta\sin\phi,\ -\cos\theta). The tangent vectors have squared lengths ∣∂r/∂θ∣2=ℓ2|\partial\mathbf{r}/\partial\theta|^2 = \ell^2 and ∣∂r/∂ϕ∣2=ℓ2sin⁡2θ|\partial\mathbf{r}/\partial\phi|^2 = \ell^2\sin^2\theta, and they are orthogonal, so Mθθ=mℓ2M_{\theta\theta} = m\ell^2, Mϕϕ=mℓ2sin⁡2θM_{\phi\phi} = m\ell^2\sin^2\theta, Mθϕ=0M_{\theta\phi} = 0 — the standard metric of a sphere, weighted by mass.

Step 3. Kinetic energy from equation (7): T=12mℓ2(θ˙2+sin⁡2θ ϕ˙2)T = \tfrac12 m\ell^2\big(\dot\theta^2 + \sin^2\theta\,\dot\phi^2\big). Insert numbers: 12mℓ2=0.5×0.20×0.25=0.025 kg m2\tfrac12 m\ell^2 = 0.5 \times 0.20 \times 0.25 = 0.025\ \mathrm{kg\,m^2}, and θ˙2+sin⁡260° ϕ˙2=1.44+0.75×4.0=4.44 s−2\dot\theta^2 + \sin^2 60°\,\dot\phi^2 = 1.44 + 0.75\times 4.0 = 4.44\ \mathrm{s^{-2}}. So T=0.025×4.44=0.111 JT = 0.025 \times 4.44 = \boxed{0.111\ \mathrm{J}}.

Step 4. Read the structure: ϕ\phi never appears in MijM_{ij} (only ϕ˙\dot\phi does) — the azimuth is on its way to becoming a cyclic coordinate with a conserved momentum, the payoff banked for a later article. The coordinate choice did that, not the physics.

2 · Bead on a spinning hoop — a rheonomic embedding (exam level)

A bead of mass m=0.10 kgm = 0.10\ \mathrm{kg} slides on a circular hoop of radius R=0.30 mR = 0.30\ \mathrm{m} that is forced to spin about its vertical diameter at constant ω=8.0 rad/s\omega = 8.0\ \mathrm{rad/s}. (a) Count the degrees of freedom and write the embedding. (b) Find the bead's kinetic energy at polar angle θ=45°\theta = 45° (from the downward vertical) when θ˙=3.0 rad/s\dot\theta = 3.0\ \mathrm{rad/s}. (c) Identify which term violates the "purely quadratic" rule and why.

Step 1. Count. One particle, 3N=33N = 3; two constraints (on the hoop's sphere of radius RR, and in the hoop's instantaneous plane), k=2k = 2; so n=1n = 1. One generalized coordinate: the polar angle θ\theta. But the hoop's plane is prescribed to rotate, so the embedding carries explicit time: r=R (sin⁡θcos⁡ωt, sin⁡θsin⁡ωt, −cos⁡θ)\mathbf{r} = R\,(\sin\theta\cos\omega t,\ \sin\theta\sin\omega t,\ -\cos\theta) — the rheonomic case of equation (2).

Step 2. Differentiate via equation (3). The θ˙\dot\theta part (along the hoop) has speed Rθ˙R\dot\theta; the drift part ∂r/∂t\partial\mathbf{r}/\partial t (azimuthal, from the forced spin) has speed ωRsin⁡θ\omega R\sin\theta. The two directions are orthogonal, so v2=R2θ˙2+ω2R2sin⁡2θv^2 = R^2\dot\theta^2 + \omega^2 R^2\sin^2\theta.

Step 3. Put in the numbers. Along-hoop: Rθ˙=0.30×3.0=0.90 m/sR\dot\theta = 0.30 \times 3.0 = 0.90\ \mathrm{m/s}. Drift: ωRsin⁡45°=8.0×0.30×0.7071=1.70 m/s\omega R \sin 45° = 8.0 \times 0.30 \times 0.7071 = 1.70\ \mathrm{m/s}. Then v2=0.81+2.88=3.69 m2/s2v^2 = 0.81 + 2.88 = 3.69\ \mathrm{m^2/s^2} and T=12×0.10×3.69=0.18 JT = \tfrac12 \times 0.10 \times 3.69 = \boxed{0.18\ \mathrm{J}}.

Step 4. Diagnose the structure: T=12mR2θ˙2+12mω2R2sin⁡2θT = \tfrac12 mR^2\dot\theta^2 + \tfrac12 m\omega^2R^2\sin^2\theta. The first term is the honest quadratic form in the single generalized velocity. The second contains no θ˙\dot\theta at all — it is the drift energy the motor forces on the bead, the fingerprint of a rheonomic constraint predicted in the kinetic-energy derivation. (No term linear in θ˙\dot\theta appears here only because drift and hoop directions happen to be orthogonal; a skewed drive would produce one.)

Step 5. The physical consequence, stated for later: because the motor continuously does work through the constraint, the bead's T+VT + V is not conserved — energy bookkeeping for rheonomic systems must wait for the energy-function machinery of a downstream article. The counting n=1n = 1, however, was never in danger: rheonomic or not, one number locates the bead.

Origins · 1 minJoseph-Louis Lagrange and how this was discovered
JL
Joseph-Louis Lagrange
1736–1813 · Turin → Berlin → Paris

Lagrange taught himself mathematics as a teenager in Turin, and by nineteen he had mailed Euler a solution of the isoperimetric problem so elegant that Euler delayed publishing his own to let the young man claim the calculus of variations. It was the perfect apprenticeship for the book that defines this article's subject: the Mécanique analytique of 1788, which he boasted contained not a single diagram — mechanics rebuilt as pure analysis, with configurations described by arbitrary parameters and the equations of motion generated by a universal recipe rather than drawn force-by-force. He succeeded Euler in Berlin at Frederick the Great's invitation ("the greatest king in Europe" required "the greatest mathematician"), survived the French Revolution as the Republic's favorite savant while colleagues lost their heads, chaired the commission that created the metric system, and taught at the new École Polytechnique. Napoleon made him a count. The coordinates qiq_i of this article are called "Lagrangian coordinates" in older books with justice: the freedom to describe a system by any independent parameters — and the guarantee that mechanics survives the relabeling — is his single most consequential idea, the one on which the next several articles are built.

History
1687
— Newton's Principia states mechanics in geometric, particle-by-particle language; constraints must be handled by explicit forces.
1743
— d'Alembert's Traité de dynamique recasts dynamics through constraint-respecting displacements, the seed of the virtual-displacement idea in equation (5).
1788
— Lagrange's Mécanique analytique rebuilds mechanics on arbitrary independent parameters — generalized coordinates — famously without a single figure.
1834
— Hamilton reformulates the theory yet again; states in configuration space acquire their velocity (and later momentum) partners, prefiguring the tangent-bundle picture.
1974
— Arnold's Mathematical Methods of Classical Mechanics completes the geometrization: configuration spaces as manifolds, kinetic energy as a Riemannian metric, dynamics as geodesic flow.

Going further

References. Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., §1.3–1.4 · Landau & Lifshitz, Mechanics, 3rd ed., §1 · Morin, Introduction to Classical Mechanics, Ch. 6 · Taylor, Classical Mechanics, Ch. 7 · Arnold, Mathematical Methods of Classical Mechanics, Ch. 4

Provenance. DOF formula n=3N−kn=3N-k and configuration-space discussion: Goldstein §1.3–1.4. Kinetic-energy quadratic form and mass metric: Landau & Lifshitz §5, Goldstein §1.4. Tangent bundle / manifold viewpoint and SO(3)SO(3) example: Arnold Ch. 4. A verifier should recompute the chain-rule velocity expression, the 'cancellation of dots' identity ∂r˙a/∂q˙i=∂ra/∂qi\partial\dot{\mathbf r}_a/\partial\dot q_i=\partial\mathbf r_a/\partial q_i, the DOF counts in the worked examples (pendulum n=1n=1 on S1S^1, double pendulum n=2n=2 on T2T^2, rigid body n=6n=6), and confirm T=12m(r˙2+r2θ˙2)T=\tfrac12 m(\dot r^2+r^2\dot\theta^2) in polar coordinates as a special case of the quadratic-form derivation. Worked-example numbers on this page (T=0.111T = 0.111 J for the spherical pendulum; v2=3.69 m2/s2v^2 = 3.69\ \mathrm{m^2/s^2} and T=0.18T = 0.18 J for the bead on the spinning hoop; M12=m2ℓ1ℓ2cos⁡(θ1−θ2)M_{12} = m_2\ell_1\ell_2\cos(\theta_1-\theta_2) for the double pendulum) are illustrative arithmetic recomputed from the stated formulas.

Builds on: Constraints and Their Classification · Degrees of Freedom. Leads to: The Lagrangian Function · The Euler-Lagrange Equations · Generalized Forces and Generalized Momenta · Cyclic Coordinates and Conservation of Conjugate Momenta · Virtual Displacements and the Principle of Virtual Work.

Is configuration space a real place or just bookkeeping?

It is exactly as real as any coordinate description of physics — which is to say, the structure is real even though the labels are conventions. No experiment can distinguish two chartings of the same Q\mathcal{Q}, but the manifold's intrinsic properties — its dimension, its topology, its mass metric — have measurable consequences. The double pendulum's torus is compact, which bounds its excursions; SO(3)SO(3)'s failure to admit one global chart forces every satellite attitude controller to handle gimbal lock or adopt quaternions; the metric MijM_{ij} decides which motions are "straight." A good test for what is real: anything that survives every change of coordinates. Dimension, topology, and the metric survive; the particular angles you chose do not.

Why exactly is the double pendulum's configuration space a torus and not, say, a square?

Because each angle is periodic: θ1=0\theta_1 = 0 and θ1=2π\theta_1 = 2\pi are the same physical configuration, so the coordinate space must glue those edges together. Start with the square [0,2π]×[0,2π][0, 2\pi] \times [0, 2\pi], glue the left edge to the right (periodicity of θ1\theta_1) to get a cylinder, then glue the top to the bottom (periodicity of θ2\theta_2): the result is a torus. The gluing is not pedantry — it changes the physics of paths. On the true torus, a motion in which the inner arm swings all the way over the top is a loop that cannot be shrunk to a point, a genuinely different class of motion from small oscillations. On a square, every loop is shrinkable, and that classification would be lost.

What is the difference between configuration space, state space, and phase space?

Three rungs of one ladder. Configuration space Q\mathcal{Q} (nn dimensions): positions only — where the system is. State space (the tangent bundle TQT\mathcal{Q}, 2n2n dimensions): positions plus velocities (q,q˙)(q, \dot q) — enough to determine the future, and the natural home of the Lagrangian L(q,q˙,t)L(q, \dot q, t). Phase space (T∗QT^*\mathcal{Q}, also 2n2n dimensions): positions plus momenta (q,p)(q, p) — the Hamiltonian description, coming several articles downstream. The middle and last rungs have equal dimension but different geometry, and the differences only start to matter once you meet transformations that mix coordinates and momenta. For now, the crucial boundary is the first one: Q\mathcal{Q} alone cannot run a movie.

Do generalized coordinates have to have units of length?

No — and most of the useful ones do not. Angles (dimensionless), arc lengths (m), areas (m²), even voltage-like variables in electromechanical analogies are all legitimate, provided the set is independent and complete. The units migrate elsewhere: in T=12Mijq˙iq˙jT = \tfrac12 M_{ij}\dot q_i\dot q_j, the mass-matrix entries absorb whatever is needed to make energy come out in joules — for an angle, Mθθ=mℓ2M_{\theta\theta} = m\ell^2 carries kg m2\mathrm{kg\,m^2}, a moment of inertia. Corollary worth flagging early: the "momentum conjugate to qiq_i" that appears downstream is ∂L/∂q˙i\partial L/\partial\dot q_i, and for an angle it is an angular momentum. If your generalized momentum fails a naive kg m/s\mathrm{kg\,m/s} unit check, that is the formalism working, not failing.

What actually goes wrong at gimbal lock — does the physics become singular?

Nothing physical goes wrong; a chart dies. Euler angles (ϕ,θ,ψ)(\phi, \theta, \psi) label orientations the way longitude and latitude label the globe, and they inherit the same disease: at the poles (θ=0\theta = 0 or π\pi for Euler angles), distinct values of ϕ\phi and ψ\psi describe the same orientation, the Jacobian of the labeling drops rank, and two of your three knobs momentarily steer the same way. The rigid body itself rotates serenely through the "lock" — ask it, and nothing happened. The lesson generalizes: SO(3)SO(3) provably cannot be covered by a single nonsingular three-parameter chart, so any three-angle attitude description must go singular somewhere. Engineering fixes are chart-management: switch between overlapping Euler charts, or embed in four parameters (quaternions) with one constraint — the standard choice in every modern flight computer.

If I choose terrible coordinates, do I get wrong physics?

No — you get correct physics wearing uncomfortable clothes. Any invertible, smooth relabeling Qi=Qi(q1,…,qn,t)Q_i = Q_i(q_1,\dots,q_n,t) (a point transformation) describes the same configurations, and the formalism being built in these articles is expressly designed so its equations keep their form under every such change — that form-invariance is the deep reason for preferring this machinery to raw Newton. What a bad choice costs is labor: coupled equations that better coordinates would separate, hidden symmetries that adapted coordinates would expose as cyclic. The spherical pendulum's ϕ\phi was invisible-by-design in Cartesian coordinates and conservation-law-adjacent in polar ones. Same physics; radically different Tuesday afternoon.

Why do rolling constraints not reduce the number of coordinates the way a rigid rod does?

Because they constrain velocities, not configurations. A rigid rod says "these two points are always 0.5 m apart" — an equation among coordinates, usable to eliminate one of them. Rolling without slipping says "the contact point's velocity is zero right now" — a condition on q˙\dot q that cannot be integrated into a relation among the qq alone. The startling consequence: a rolling sphere can reach every position with every orientation (park a car anywhere at any heading despite never sliding sideways), so its configuration space stays full-dimensional even though its instantaneous freedoms are fewer. Such constraints are called nonholonomic, and they get a full classification — and their own warning labels — in the constraints article.

Is the mass matrix the same as the moment of inertia?

The moment of inertia is one corner of it. For a single coordinate that happens to be a rotation angle, Mθθ=∑ama ∣∂ra/∂θ∣2M_{\theta\theta} = \sum_a m_a\,|\partial\mathbf{r}_a/\partial\theta|^2 is exactly the moment of inertia about that axis — the spherical pendulum's mℓ2m\ell^2, for instance. But MijM_{ij} is the general-purpose object: it handles mixed coordinate types (an angle and a length in the same system), it varies with configuration (the double pendulum's cos⁡(θ1−θ2)\cos(\theta_1 - \theta_2) coupling), and its off-diagonal entries encode inertia shared between coordinates, for which elementary mechanics has no standard name. Think of the familiar scalars — mass for translation, moment of inertia for rotation — as the diagonal, constant special cases of one configuration-dependent metric tensor.

Before this

Where this leads