computo ergo sum

2026-10-01 · chapter Yang–Mills and the mass gapPart 2 Physical preparationthe action principle and Noether

The whole book Contents and reading orders
  1. Part 0 Entrance — how to read this book
  2. 0-01 How to read this book
  3. 0-02 One picture
  4. 0-03 Reading the problem statement word by word
  5. 0-04 Four dimensions in mathematics and four dimensions in physics
  6. Part 1 Mathematical preparation
  7. 1-01 Vector spaces and linear maps
  8. 1-02 What a group is
  9. 1-03 Lie groups and Lie algebras
  10. 1-04 SU(2) and SU(3)
  11. 1-05 Manifolds and tangent spaces
  12. 1-06 Differential forms and the exterior derivative
  13. 1-07 Vector bundles and connections
  14. 1-08 Probability and measure
  15. 1-09 Probability in infinite dimensions
  16. 1-10 Hilbert spaces and self-adjoint operators
  17. 1-11 Fourier analysis and distributions
  18. 1-12 Detour: what happens only in four dimensions
  19. Part 2 Physical preparation
  20. 2-01 Lagrangians and the action
  21. 2-02 The idea of a field
  22. 2-03 Electromagnetism is a gauge theory
  23. 2-04 Special relativity and spacetime
  24. 2-05 The bones of quantum mechanics
  25. 2-06 The harmonic oscillator and creation and annihilation
  26. 2-07 The idea of the path integral
  27. 2-08 Statistical mechanics and phase transitions
  28. 2-09 Detour: the gap and the correlation length in the Ising model
  29. 2-10 The relaxation time and the dynamical exponent z
  30. Part 3 Yang–Mills theory (classical)
  31. 3-01 The gauge principle
  32. 3-02 Non-abelian gauge fields
  33. 3-03 The action and the equations
  34. 3-04 Gauge theory as geometry
  35. 3-05 Instantons and topology
  36. 3-06 Detour: Donaldson and exotic R⁴
  37. 3-07 Yang–Mills inside the Standard Model
  38. Part 4 Quantisation
  39. 4-01 Canonical quantisation and the Hamiltonian
  40. 4-02 The path integral and Euclideanisation
  41. 4-03 Perturbation theory and Feynman diagrams
  42. 4-04 Divergences and renormalisation
  43. 4-05 The algebra of divergences
  44. 4-06 Asymptotic freedom
  45. 4-07 Dimensional transmutation and Λ
  46. 4-08 Gauge fixing and Faddeev–Popov
  47. 4-09 The axioms of quantum field theory
  48. 4-10 The Osterwalder–Schrader axioms and reconstruction
  49. 4-11 The definition of the mass gap
  50. 4-12 Detour: the perturbation series diverges
  51. Part 5 Lattice gauge theory
  52. 5-01 Wilson's lattice
  53. 5-02 The strong-coupling expansion
  54. 5-03 Reflection positivity and the transfer matrix
  55. 5-04 The infinite-volume limit and the cluster expansion
  56. 5-05 Weak coupling and the continuum limit
  57. 5-06 How U(1) differs from the non-abelian case
  58. 5-07 The Monte Carlo method
  59. 5-08 Measuring glueballs and the string tension
  60. 5-09 What to hold fixed when taking the limit
  61. 5-10 Lattice gauge theory with finite groups
  62. 5-11 Detour: Bakry–Émery and one-link integration
  63. Part 6 Constructive field theory
  64. 6-01 What constructive field theory is
  65. 6-02 Solvability in two dimensions and the Yang–Mills measure
  66. 6-03 Constructing scalar fields
  67. 6-04 The cluster expansion
  68. 6-05 The sequence of renormalisation group steps
  69. 6-06 Ultraviolet stability of three-dimensional Yang–Mills
  70. 6-07 Four dimensions
  71. 6-08 Triviality of φ⁴ in four dimensions
  72. 6-09 Stochastic quantisation and regularity structures
  73. 6-10 Where four dimensions stops
  74. 6-11 Obstacles other than divergence
  75. 6-12 Detour: "if 1/0 is given a symbol"
  76. Part 7 From the side of physics
  77. 7-01 How physics sees it
  78. 7-02 Mechanisms of confinement
  79. 7-03 The string picture
  80. 7-04 Large N
  81. 7-05 String theory and extra dimensions
  82. 7-06 How extra dimensions become invisible
  83. 7-07 Where gauge fields come from
  84. 7-08 Holography
  85. 7-09 The mass gap becomes geometry
  86. 7-10 Back to four dimensions
  87. 7-11 Supersymmetry and Seiberg–Witten
  88. 7-12 The sentence that plays the role of an equivalence principle
  89. 7-13 Where the way physics digs parts from mathematics
  90. Part 8 A dictionary of two languages — the physical and the mathematical view
  91. 8-01 How to read the dictionary
  92. 8-02 Dictionary A
  93. 8-03 Dictionary B
  94. 8-04 Dictionary C
  95. 8-05 The physical counterparts of our assumptions
  96. Part 9 Where things stand, and open tasks
  97. 9-01 How far the world has come
  98. 9-02 Where the two ways of digging part
  99. 9-03 Candidates for new concepts
  100. 9-04 The list of open tasks
  101. 9-05 Common misconceptions
  102. Part 10 Towards a proof of the mass gap — building the first step of the staircase
  103. 10-00 Entrance — the staircase, and the random walk, the ruler and the clock
  104. 10-01 The line of construction and Assumption H
  105. 10-02 Roads tried on the four-dimensional lattice, and walls
  106. 10-03 The four steps of the staircase — why we start from 0+1 dimensions
  107. 10-04 The path measure of the double well, and kinks
  108. 10-05 The statement of the first step — four conditions and Proposition 4.5
  109. 10-06 (H1) Local Poincaré within a phase
  110. 10-07 The map of (H2) — seven routes and the seam of the ring
  111. 10-08 What computer-assisted means — the conditions on the chain inside the window
  112. 10-09 The sector framework and the three strands of route 5
  113. 10-10 The list of detour plans — types of breakage
  114. 10-11 (H3)(H4) and the constant of the core
  115. 10-12 The integrated table and the dependency graph
  116. 10-13 Lean statements and the limits of machine checking
  117. 10-14 Reconciliation — the numbers of the record and the known numbers
  118. 10-15 What remains and its cost — the distance to a theorem
  119. 10-16 What is not claimed
  120. 10-17 Insights and small questions
  121. Part 11 Further considerations — beyond Assumption H, and other views
  122. 11-00 Further considerations — entrance
  123. 11-01 One form that obstructs H — a nearly conserved global quantity
  124. 11-02 Choosing the ruler and the coordinates again — the length at which conservation breaks
  125. 11-03 Two ways to repair H — the sector H and the local H
  126. 11-04 Which H the one-sided bridge needs
  127. 11-05 Making time continuous — the Hamiltonian lattice and the vacuum measure
  128. 11-06 Paying with rarity instead of convexity — renormalisation on groups and log-Sobolev
  129. 11-07 Because it is four dimensions — the dimensionless coupling and topology
  130. 11-08 Questions, and what is not claimed
  131. 11-A1 Table of notation for Part 11
  132. Appendices
  133. A-01 Table of notation
  134. A-02 Glossary
  135. A-03 Guide to the literature
  136. A-04 Lean and machine checking
  137. A-05 How this book was made
  138. A-06 Index of assumptions

Lagrangians and the action — motion can be written as a stationary point of "one number attached to the whole path"

What this chapter shows — That the equations of motion can be written as the stationarity condition of the action. That the Euler–Lagrange equations come out of it. That continuous symmetries give conserved quantities (Noether's theorem), and local symmetries give constraints.

Prerequisite chapters — 1-01 Vector spaces and linear maps (rotation matrices, eigenvalues of real symmetric matrices). Otherwise, multivariable calculus and Newton's equation of motion.

Said first — We restrict to finitely many degrees of freedom (fields are in 2-02, the Yang–Mills action in 3-03). The proofs stay within smooth paths and variations with fixed ends. We do not answer why the action principle holds.

Order of this chapter
  1. Why start from the action — three advantages
  2. The action, a single number — the integral of kinetic energy − potential energy
  3. Shifting the path a little — discretised, a function of many variables (Figure 1)
  4. The Euler–Lagrange equations — one integration by parts; no choice of coordinates
  5. Not necessarily a minimum — eigenvalues and conjugate points (with computations)
  6. Noether's theorem — conserved quantities from symmetries (Figure 2)
  7. Detour: the discretised action has Noether too — why angular momentum is conserved exactly
  8. Local symmetries give constraints — a foretaste of gauge theory

01

Why start from the action

Newton's equation of motion mq¨=Fm\ddot q = F is a local rule that fixes the acceleration at each moment. In this chapter we write the same motion in another form. Fix the start and the arrival, and give the whole path between them a single number SS (the action). The actual motion is chosen as the path for which SS does not change at first order when the path is shifted slightly. This is the action principle (also called the principle of least action).

There are three advantages. It does not depend on coordinates (§04). Conservation laws come out automatically from symmetries (§06). It is the language of field theory: Yang–Mills theory is defined by an action, the integral of tr⁡F2\operatorname{tr}F^2 (3-03), and the path integral too sums over paths with the action in the exponent (2-07).

As a metaphor, the equation of motion is "step-by-step directions", and the action principle is the view of "comparing whole routes on a map". To be exact, for systems with conservative forces both give the same equation (§04), and the only difference is what is easy to see.

02

The action, a single number

For a particle on a line, let q(t)q(t) be its position at time tt, q˙=dq/dt\dot q = dq/dt its velocity, mm its mass and V(q)V(q) the potential energy. The Lagrangian LL is the kinetic energy minus the potential energy, a different quantity from the energy, which is their sum. The number obtained by integrating it along the path is the action.

L(q,q˙)=12mq˙2−V(q)L(q,\dot q) = \tfrac12 m\dot q^2 - V(q) S[q]=∫t0t1L(q(t),q˙(t)) dtS[q] = \int_{t_0}^{t_1} L\bigl(q(t),\dot q(t)\bigr)\,dt

The square brackets show that it is a map that takes a whole path and returns a number (a functional).

Example. Take V=0V = 0, m=1m = 1, and ends q(0)=0q(0) = 0, q(1)=1q(1) = 1. For the uniform path q=tq = t, S=0.5S = 0.5. For the detour q=t+0.3sin⁡πtq = t + 0.3\sin\pi t joining the same ends, S=12∫01(1+0.3πcos⁡πt)2 dt=12(1+0.09π2/2)≈0.722S = \tfrac12\int_0^1 (1 + 0.3\pi\cos\pi t)^2\,dt = \tfrac12(1 + 0.09\pi^2/2) \approx 0.722. Taking a detour increases the action.

03

Shifting the path a little

An extremum of a function was a point where the first-order change under a small change of the variable is 0. For the action we move the path. Taking a function η(t)\eta(t) that vanishes at the ends, S(ε)=S[q+εη]S(\varepsilon) = S[q + \varepsilon\eta] is an ordinary function of the real number ε\varepsilon. The action principle is the requirement "S′(0)=0S'(0) = 0 for every η\eta with fixed ends", and a path satisfying it is called a stationary point.

Discretised, it becomes a function of finitely many variables

Dividing time into NN equal parts (step hh), making the path a broken line and approximating the integral by the trapezoidal rule, the action becomes a function of many variables, the intermediate positions q1,…,qN−1q_1, \dots, q_{N-1} (the ends q0q_0, qNq_N fixed).

SN=∑k=0N−1[m(qk+1−qk)22h−h2(V(qk)+V(qk+1))]S_N = \sum_{k=0}^{N-1}\Bigl[\frac{m(q_{k+1}-q_k)^2}{2h} - \frac h2\bigl(V(q_k)+V(q_{k+1})\bigr)\Bigr] ∂SN∂qk=−h[m qk+1−2qk+qk−1h2+V′(qk)]\frac{\partial S_N}{\partial q_k} = -h\Bigl[m\,\frac{q_{k+1}-2q_k+q_{k-1}}{h^2} + V'(q_k)\Bigr]

The action principle becomes "the gradient is 0", and this is a difference version of mq¨=−V′(q)m\ddot q = -V'(q) (used again in §07).

Example: a particle on a spring

Take V=12mω2q2V = \tfrac12 m\omega^2 q^2, m=1m = 1, time length T=1T = 1, and ends q(0)=1q(0) = 1, q(1)=−1q(1) = -1. The solution of the equation of motion is qcl(t)=sin⁡(ω(12−t))/sin⁡(ω/2)q_{\mathrm{cl}}(t) = \sin\bigl(\omega(\tfrac12 - t)\bigr)/\sin(\omega/2) (cl stands for "classical path"). For the shift ηn=sin⁡nπt\eta_n = \sin n\pi t (n=1,2,3n = 1, 2, 3), the action is a quadratic in ε\varepsilon.

S(ε)=ωcot⁡ω2+ε24((nπ)2−ω2)S(\varepsilon) = \omega\cot\frac{\omega}{2} + \frac{\varepsilon^2}{4}\bigl((n\pi)^2 - \omega^2\bigr)

The absence of a first-order term is what "stationary" means. The second-order term comes from 12∫(η˙2−ω2η2) dt\tfrac12\int(\dot\eta^2 - \omega^2\eta^2)\,dt. Moving ε\varepsilon in Figure 1 moves the action along a parabola, with slope 0 at ε=0\varepsilon = 0.

Figure 1 Moving ε, the shape n and ωT changes the shifted path (left) and the action S(ε) (right)
0.60

ComputedThe curves draw the formula above in the browser, and the formula was checked against numerical integration on this machine (Table 1). The dotted line is the classical path, and the dots are the variables when discretised with N=10N = 10.

Table 1 Values of the action (m=T=1m = T = 1, q(0)=1q(0) = 1, q(1)=−1q(1) = -1)
ωT/π\omega T/\pinnS(0)S(0)S(1)S(1)Coefficient of the second-order term
0.511.57083.4213+1.851
1.010.00000.00000
1.51−4.7124−7.7966−3.084
1.52−4.7124−0.3944+4.318

ComputedThe difference between numerical integration by Simpson's rule (20000 intervals) and the formula is at most 4×10−154\times10^{-15}. The first-order coefficient is also at most 5×10−135\times10^{-13}.

04

The Euler–Lagrange equations

We compute the condition of §03 without discretising. Differentiate with respect to ε\varepsilon and integrate the second term by parts (the boundary term vanishes because η=0\eta = 0 at the ends).

S′(0)=∫t0t1[∂L∂q η+∂L∂q˙ η˙]dtS'(0) = \int_{t_0}^{t_1}\Bigl[\frac{\partial L}{\partial q}\,\eta + \frac{\partial L}{\partial\dot q}\,\dot\eta\Bigr]dt =∫t0t1[∂L∂q−ddt∂L∂q˙]η dt= \int_{t_0}^{t_1}\Bigl[\frac{\partial L}{\partial q} - \frac{d}{dt}\frac{\partial L}{\partial\dot q}\Bigr]\eta\,dt

If this is 0 for every η\eta, the bracket is 0 at each time (the bracket is continuous, so if it is positive somewhere it is positive nearby, and taking a η\eta positive only nearby makes the integral positive: the fundamental lemma of the calculus of variations). This gives the Euler–Lagrange equations. With several degrees of freedom there is one for each component.

ddt∂L∂q˙−∂L∂q=0(EL)\frac{d}{dt}\frac{\partial L}{\partial\dot q} - \frac{\partial L}{\partial q} = 0 \tag{EL}

Example. If L=12mq˙2−V(q)L = \tfrac12 m\dot q^2 - V(q), then mq¨=−V′(q)m\ddot q = -V'(q): Newton's equation for a conservative force itself.

No choice of coordinates

This derivation does not use what qq represents. It may be an angle (generalised coordinates), and equations of the same form hold Known. This is because the same path gets the same action whatever the coordinates (1-05). Writing a pendulum with the angle θ\theta, L=12ml2θ˙2+mglcos⁡θL = \tfrac12 ml^2\dot\theta^2 + mgl\cos\theta gives ml2θ¨=−mglsin⁡θml^2\ddot\theta = -mgl\sin\theta. The tension does not appear.

Adding dG(q,t)/dtdG(q,t)/dt to LL does not change the equations (the action changes only by a quantity fixed by the end values) Known. The θ term of Yang–Mills theory is the four-dimensional version of this, with the form of a four-dimensional divergence (a total divergence) instead of a total time derivative. To be exact, the θ term does not change the classical equations, but it matters in quantum theory through the topology of the field (3-05).

05

Not necessarily a minimum

The action principle asks only that the first-order change be 0; whether it is a minimum is decided by the second-order term. In the example of §03 it is the sign of the coefficient ((nπ)2−ω2)/4((n\pi)^2 - \omega^2)/4. If ωT<π\omega T \lt \pi it is positive for every nn, and a general η\eta, expanded in a sine series, splits into a sum over each nn, so the classical path really is a minimum (because the action is quadratic).

When ωT\omega T exceeds π\pi, the coefficient of n=1n = 1 becomes negative, and in that direction the action decreases (a saddle point). Moving ωT/π\omega T/\pi either side of 1 in Figure 1 turns the parabola of n=1n = 1 over. At exactly ωT=π\omega T = \pi, moving in the direction of n=1n = 1 does not change the action (it increases in the direction of n≥2n \ge 2), and the solutions joining the same ends form a family cos⁡πt+csin⁡πt\cos\pi t + c\sin\pi t. In general, a time at which a family of nearby solutions leaving the same point at time 0 gathers again at one point is called a conjugate point of 0 Known. The conjugate points in this example are t=π/ω,2π/ω,…t = \pi/\omega, 2\pi/\omega, \dots, and ωT=π\omega T = \pi is the case where the final time falls exactly on the first conjugate point.

The second-order term is an eigenvalue problem for a matrix

Discretised, the second-order term is the Hessian matrix Hjk=∂2SN/∂qj∂qkH_{jk} = \partial^2 S_N/\partial q_j\partial q_k. It is real symmetric, so it can be diagonalised (1-01), and a local minimum means all eigenvalues are positive. The number of negative eigenvalues is the number of directions in which the action can be reduced, called the Morse index (here a name for the discretised matrix). For the harmonic oscillator HH is tridiagonal (diagonal 2m/h−mω2h2m/h - m\omega^2 h, neighbours −m/h-m/h); its eigenvalues divided by hh approach m((nπ/T)2−ω2)m((n\pi/T)^2 - \omega^2), and ⌊ωT/π⌋\lfloor\omega T/\pi\rfloor of them are negative.

Table 2 Hessian matrix of the discretised action (m=T=1m = T = 1, N=200N = 200)
ωT/π\omega T/\piNumber of negative eigenvaluesSmallest eigenvalue / h/\,hπ2−(ωT)2\pi^2 - (\omega T)^2
0.5007.40207.4022
1.051−1.0118−1.0116
2.502−51.8156−51.8154
3.203−91.1953−91.1951

ComputedAll eigenvalues found with numpy's eigvalsh. With N=20N = 20 the numbers are the same, and the smallest eigenvalues differ by about 0.02.

Without discretising too, the directions that reduce the action in this example are those among sin⁡nπt\sin n\pi t with nπ<ωTn\pi \lt \omega T, and their number equals the number of conjugate points kπ/ωk\pi/\omega in (0,T)(0, T). This is an instance of the Morse index theorem (the number of negative directions of the second-order term equals the number of conjugate points in (0,T)(0, T), counted with multiplicity) Known.

This view is extended to infinite dimensions in 1-10. The smallest eigenvalue here (of the second-order term of the action) and the mass gap of 4-11 (the gap in the spectrum of the Hamiltonian) are different things.

06

Noether's theorem

If LL does not contain the coordinate qq, (EL) says that p=∂L/∂q˙p = \partial L/\partial\dot q (the conjugate momentum) is conserved. L=12m(r˙2+r2φ˙2)−V(r)L = \tfrac12 m(\dot r^2 + r^2\dot\varphi^2) - V(r), a central force written in polar coordinates, does not contain φ\varphi, and the angular momentum mr2φ˙mr^2\dot\varphi is conserved. Noether's theorem is this detached from coordinates.

For transformations q↦Φs(q)q \mapsto \Phi_s(q) labelled continuously by a real number ss and equal to the identity at s=0s = 0 (translations, rotations by the angle ss), let K(q)=ddsΦs(q)∣s=0K(q) = \frac{d}{ds}\Phi_s(q)\big|_{s=0} be the direction of motion. Noether's theorem: if L(Φs(q),ddtΦs(q))=L(q,q˙)L\bigl(\Phi_s(q), \tfrac{d}{dt}\Phi_s(q)\bigr) = L(q, \dot q) for every path, then the following QQ is constant along the motion Known.

Q=∑i∂L∂q˙i Ki(q)(N)Q = \sum_i \frac{\partial L}{\partial\dot q_i}\,K_i(q) \tag{N}

Proof. Differentiating the assumption with respect to s=0s = 0 gives ∑i(∂L/∂qi)Ki+∑i(∂L/∂q˙i) K˙i=0\sum_i (\partial L/\partial q_i) K_i + \sum_i (\partial L/\partial\dot q_i)\,\dot K_i = 0. Replacing ∂L/∂qi\partial L/\partial q_i by ddt(∂L/∂q˙i)\frac{d}{dt}(\partial L/\partial\dot q_i) using (EL), the left-hand side is the derivative of a product dQ/dtdQ/dt. ∎

Transformation that does not change the actionKKConserved quantity QQ
Translation x↦x+se\boldsymbol x \mapsto \boldsymbol x + s\boldsymbol ee\boldsymbol eComponent of momentum mx˙⋅em\dot{\boldsymbol x}\cdot\boldsymbol e
Rotation in the plane x↦R(s)x\boldsymbol x \mapsto R(s)\boldsymbol x(−y,x)(-y, x)Angular momentum Lz=m(xy˙−yx˙)L_z = m(x\dot y - y\dot x)
Translation in time (∂L/∂t=0\partial L/\partial t = 0)—Energy E=∑iq˙i ∂L/∂q˙i−LE = \sum_i \dot q_i\,\partial L/\partial\dot q_i - L

The third row follows by differentiating directly, from dE/dt=−∂L/∂tdE/dt = -\partial L/\partial t (E=12mq˙2+VE = \tfrac12 m\dot q^2 + V).

Checking with a figure: breaking the symmetry removes the conserved quantity

Give a particle in the plane (m=1m = 1) the potential V=12r2+14r4(1+δcos⁡4θ)V = \tfrac12 r^2 + \tfrac14 r^4(1 + \delta\cos 4\theta) in polar coordinates. If δ=0\delta = 0 it is rotationally symmetric, and LzL_z is conserved. With δ>0\delta \gt 0 the rotations that leave VV unchanged are only multiples of 90°, and the symmetry drops to a discrete one. There are no infinitesimally small transformations for the theorem to use, so the theorem says nothing, and in fact LzL_z is not conserved. VV does not depend on time, so EE is conserved.

Figure 2 Moving δ makes the angular momentum Lz start to fluctuate, while the energy E hardly fluctuates (about 10⁻⁵)δ = 0.00

ComputedInitial position (1,0)(1, 0), velocity (0,0.9)(0, 0.9). The orbit is computed by the browser with the velocity Verlet method (h=0.01h = 0.01), and Table 3 gives the values from running the same procedure on this machine (1.4e−14 means 1.4×10−141.4\times10^{-14}). The two bands on the right have the same vertical scale.

Table 3 Changes of angular momentum and energy (0≤t≤600 \le t \le 60)
δ\deltahhMinimum of LzL_zRange of LzL_zmax⁡∣E−E0∣\max\lvert E - E_0\rvert
00.010.90001.4e−142.0e−5
00.0050.90001.6e−145.0e−6
0.10.010.85810.0802.8e−5
0.30.01−1.01522.034.3e−5
0.50.01−1.09382.195.8e−5

For δ=0.3\delta = 0.3 and 0.50.5 even the sign of LzL_z changes, and the particle reverses its direction of rotation partway. The deviation of EE is about 10−510^{-5}, and looks flat in the figure.

07

Detour: the discretised action has Noether too

Detour

A section you may skip.

Figure 2 is computed step by step with step size hh. Here position qkq_k and velocity vkv_k are vectors in the plane (qk∈R2q_k \in \mathbb R^2; the qkq_k of §03 were numbers), and they are advanced by the velocity Verlet method using the force F=−∇VF = -\nabla V.

vk+1/2=vk+h2mF(qk),qk+1=qk+h vk+1/2,vk+1=vk+1/2+h2mF(qk+1)v_{k+1/2} = v_k + \frac{h}{2m}F(q_k),\qquad q_{k+1} = q_k + h\,v_{k+1/2},\qquad v_{k+1} = v_{k+1/2} + \frac{h}{2m}F(q_{k+1})

Even so, in Table 3 at δ=0\delta = 0, LzL_z is constant to the level of rounding error, while the energy deviates by 10−510^{-5}. There is a reason.

Eliminating the velocity gives m(qk+1−2qk+qk−1)=h2F(qk)=−h2∇V(qk)m(q_{k+1} - 2q_k + q_{k-1}) = h^2F(q_k) = -h^2\nabla V(q_k). This is exactly the stationarity condition of the discretised action of §03 carried over to the plane Known (Hairer–Lubich–Wanner 2003, §1.6; this form is called the Störmer–Verlet method, and the velocity Verlet method is one way of writing it). If VV is rotationally symmetric, the discretised action too is unchanged by a simultaneous rotation, the computation of §06 goes through with differences, and the following JJ is constant regardless of kk.

J=mh (qk×qk+1)zJ = \frac{m}{h}\,(q_k \times q_{k+1})_z

Here (a×b)z=axby−aybx(a \times b)_z = a_xb_y - a_yb_x. It can also be checked directly. Multiplying the stationarity condition by qk ×q_k\,\times on the left, for a central force F(qk)F(q_k) is parallel to qkq_k, so the right-hand side vanishes and qk×qk+1=qk−1×qkq_k \times q_{k+1} = q_{k-1} \times q_k remains. JJ is 2m/h2m/h times the area of the triangle made by the origin and two neighbouring points, so the area swept in each step is constant.

This JJ is exactly the Lz=m(qk×vk)zL_z = m(q_k \times v_k)_z of Table 3. The first two update formulas give qk+1=qk+h vk+h22mF(qk)q_{k+1} = q_k + h\,v_k + \frac{h^2}{2m}F(q_k). For a central force F(qk)F(q_k) is parallel to qkq_k, so qk×qk+1=h qk×vkq_k \times q_{k+1} = h\,q_k \times v_k, hence J=LzJ = L_z. Theorem 3.5 of the same paper shows in general that quadratic conserved quantities such as angular momentum are preserved exactly by this method Known. And §1.3, citing Feynman's lectures, notes that the argument by which the "Principia" showed Kepler's second law has the same form as one way of writing this method (half a step at constant velocity, a kick towards the centre, half a step at constant velocity).

On the other hand, time has been discretised, so it can be shifted only by integer multiples of hh. The theorem does not work for discrete symmetries, and energy is not conserved exactly. However, if the formula for the energy (a function of position and velocity) is smooth and the numerical solution stays within some bounded region, the following has been shown for any positive integer NN Known (ibid. Theorem 5.1; CC and CNC_N do not depend on tt and hh).

∣Ek−E0∣≤Ch2+CNhNt(0≤t=kh≤h−N)\lvert E_k - E_0\rvert \le Ch^2 + C_N h^N t \qquad (0 \le t = kh \le h^{-N})

The second term is at most CNh2C_N h^2 if t≤h2−Nt \le h^{2-N}. Taking NN large, the error stays of order h2h^2 up to very long times, high powers of h−1h^{-1}. It is not a theorem about infinite time. In Table 3, halving hh makes the deviation about 1/4.

"Keep the symmetry you want to keep exactly, as a symmetry of the discretised action" is also the starting point of Wilson's method in lattice gauge theory (5-01). What is kept exactly there is gauge symmetry, and what is broken are the continuous rotations and translations of spacetime (four dimensions turned into imaginary time; 4-02). The 90° rotations of the lattice directions and translations by the lattice spacing remain.

08

Local symmetries give constraints

The ss of §06 was a constant (a global symmetry). A symmetry under which the action is unchanged even when transforming by a s(t)s(t) that differs from time to time is a local symmetry, and it is at the centre of gauge theory (3-01). We look at a toy with two variables, x(t)x(t) and λ(t)\lambda(t).

L=12m(x˙−λ)2,x↦x+s(t),λ↦λ+s˙(t)L = \tfrac12 m(\dot x - \lambda)^2, \qquad x \mapsto x + s(t),\quad \lambda \mapsto \lambda + \dot s(t)

For any s(t)s(t), x˙−λ\dot x - \lambda does not change. λ\lambda appears without time derivatives, so its equation is a constraint.

(λ)  m(x˙−λ)=0,(x)  ddt[m(x˙−λ)]=0(\lambda)\ \ m(\dot x - \lambda) = 0, \qquad (x)\ \ \frac{d}{dt}\bigl[m(\dot x - \lambda)\bigr] = 0

The second is the derivative of the first, and only one is independent. xx remains as an arbitrary function that the equations do not fix. Noether's second theorem (1918) states that when the action is unchanged under transformations containing arbitrary functions, some of the Lagrange equations follow from the rest Known. To be exact, as many identities as there are arbitrary functions hold among the left-hand sides of the Lagrange equations and their derivatives. In this toy the arbitrary function is the one s(t)s(t), and the identity is also one: "the left-hand side for (x)(x) is the time derivative of the left-hand side for (λ)(\lambda)". The Q=m(x˙−λ)Q = m(\dot x - \lambda) made by (N) from a constant ss is, far from conserved, constrained to be 0.

In electromagnetism the scalar potential ϕ\phi plays the role of λ\lambda, the vector potential A\mathbf A the role of xx, and Gauss's law plays the role of the equation for λ\lambda (2-03, 4-01). To be exact, in electromagnetism s(t)s(t) becomes a function at each point χ(t,x)\chi(t,\mathbf x), A\mathbf A shifts by ∇χ\nabla\chi and ϕ\phi by −∂χ/∂t-\partial\chi/\partial t ((G) of 2-03; the sign of the shift of ϕ\phi is opposite to the λ\lambda of the toy). The degrees of freedom the equations do not fix are a redundancy of description, and handling them is one of the big problems of quantising Yang–Mills theory (4-08).

The toy itself has no physical content (the only solution is λ=x˙\lambda = \dot x, with xx arbitrary).


What this chapter can say

MatterWhat is saidWhat is not said
The Euler–Lagrange equationsKnownDerived for smooth paths with fixed endsVariations with moving ends, the general theory of constraints
Stationary and minimumComputedFor the harmonic oscillator, a minimum if ωT<π\omega T \lt \pi, a saddle point beyond (Tables 1 and 2)The general theory of conjugate points (Jacobi's condition, the proof of the Morse index theorem)
Noether's first theoremKnownProved for finitely many degrees of freedom when LL is exactly invariantThe case where LL changes by a total derivative, the form in field theory (2-02)
Numerical observationComputedIn one example (t≤60t \le 60), the fluctuation and conservation of LzL_z and the error ∝h2\propto h^2 of EE were observed (Table 3)Not a general claim. Long-time estimates rely on theorems in the literature
Noether's second theoremThe redundancy of equations and constraints seen in one toyThe general proof and the form in field theory (3-01)
The origin of the action principle—Physics2-07

Sources and reproduction

ItemKindSource / tools
(EL), generalised coordinates, (N)KnownTextbookL. D. Landau, E. M. Lifshitz, Mechanics, 3rd ed., Pergamon Press, 1976. Bibliographic data only
Noether's two theoremsKnownPrimary sourceE. Noether, Invariante Variationsprobleme, Nachr. Ges. Wiss. Göttingen, Math.-phys. Kl. (1918) 235–257. English translation M. A. Tavel, Invariant variation problems, Transport Theory and Statistical Physics 1(3) (1971) 186–207, doi:10.1080/00411457108231446, arXiv:physics/0503066 (Theorems I and II of §1 checked)
Second variation, conjugate pointsKnownTextbookD. Liberzon, Calculus of Variations and Optimal Control Theory, Princeton University Press, 2012, doi:10.1515/9781400842643 (Jacobi's condition and conjugate points checked in §2.6.2 of the pre-publication version). I. M. Gelfand, S. V. Fomin, Calculus of Variations, 1963. Bibliographic data only
The Morse index theoremKnownTextbookJ. Milnor, Morse Theory, Annals of Mathematics Studies 51, Princeton University Press, 1963 (§15 Theorem 15.1 checked; what is stated is the case of geodesics)
On the Störmer–Verlet method of §07KnownReviewE. Hairer, C. Lubich, G. Wanner, Acta Numerica 12 (2003) 399–450, doi:10.1017/S0962492902000144 (§1.3, §1.6, Theorems 3.5 and 5.1 checked)
Variational integratorsKnownReviewJ. E. Marsden, M. West, Acta Numerica 10 (2001) 357–514. Bibliographic data only
Tables 1, 2, 3Computedpython3 + numpy (the record kept on this machine). action_scan.py (Simpson's rule), hessian.py (eigvalsh), noether_verlet.py (the velocity Verlet method; Figure 2 runs the same procedure in the browser)

Next chapter: 2-02 The idea of a field (replacing qiq_i by a field ϕ(x)\phi(x), the degrees of freedom become infinitely many).

← 1-12 Detour: what happens only in four dimensionsContents2-02 The idea of a field →

Revised 2026-10-01: first version.