The absence of magnetic monopoles is expressed by . In a regular simply connected region, a magnetic vector potential may be introduced so that
Its units are . It is not unique: adding a gradient leaves the magnetic field unchanged because . In time-dependent electromagnetism,
is a gauge transformation preserving both physical fields. Stokes theorem gives for a contour bounding an appropriate surface.
For localized steady currents, choose the Coulomb gauge . Substitution into Ampère's law and the curl-of-curl identity gives
Boundary conditions are needed in addition to the gauge choice. For a uniform magnetic field one convenient choice is . For a magnetic dipole, outside its source. The magnetic vector potential simplifies field calculations and appears directly in the Lagrangian of a charged particle, . In multiply connected regions, a field-free region can retain nontrivial circulation of , relevant to the Aharonov-Bohm effect; a local gauge argument cannot remove every global circulation.
In Lagrangian mechanics, generalized coordinates describe the independent configurations allowed by the constraints. For particles with ordinary conservative forces in an inertial frame, use , with kinetic energy and potential energy . Holonomic constraints can often be built into the coordinates, avoiding explicit constraint forces.
The principle of stationary action makes the action stationary for variations vanishing at the endpoints. Integrating the velocity variation by parts gives
Since each independent variation is arbitrary, the Euler-Lagrange equations are
If nonconservative generalized forces are present, their components appear on the right. The formalism is invariant under suitable coordinate changes and adding a total time derivative to leaves the fixed-endpoint equations unchanged.
The generalized momentum is . If is cyclic, and is conserved. If has no explicit time dependence, the conserved energy function is
For a kinetic energy quadratic in velocities and a velocity-independent potential energy, this equals . This connects symmetry to conservation laws.
For example, a simple pendulum of mass and length has
Its tension does not appear because the fixed-length constraint is already imposed. Expanding for small gives simple harmonic motion with . More generally, expansion of and to quadratic order about a stable equilibrium leads directly to the normal-mode eigenvalue problem.
Definition and equations. A normal mode is a motion in which every component of a linear system has the same time dependence with fixed amplitude ratios and phases. For a conservative system near equilibrium, the most useful form is
The vector describes the mode shape, while is its angular frequency. Linear combinations of modes describe arbitrary sufficiently small motions when the modes form a complete basis.
Let be displacements of independent generalized coordinates from equilibrium. Expand the Lagrangian to quadratic order:
The mass matrix defines a positive-definite quadratic form on the independent coordinates, while the stiffness matrix is the Hessian of the potential energy. Substituting a mode gives the generalized eigenvalue problem for small oscillations
A practical calculation identifies the equilibrium, chooses independent coordinates, computes the quadratic kinetic and potential energies, solves for the eigenvalues and then obtains their eigenvectors. One can equivalently diagonalize the real symmetric matrix . This guarantees real and a complete orthogonal basis in mass-weighted coordinates.
For distinct eigenvalues, the fact that and are real symmetric matrices implies
Thus distinct modes are orthogonal in the inner product weighted by the mass matrix. Choosing decouples the motion into , with total energy . A positive-definite matrix gives positive eigenvalues and oscillatory stability. Zero modes can describe rigid translations or rotations without restoring force; a negative eigenvalue gives exponential instability rather than an oscillatory mode.
Symmetry. A symmetry preserving both quadratic energies maps modes into modes of the same frequency. Coordinates adapted to symmetry divide the eigenvalue problem into smaller independent blocks. For two identical subsystems, exchanging them gives symmetric and antisymmetric modes. The degeneracy of normal modes means that independent modes share one frequency. Higher-dimensional group representations of a symmetry can enforce such degeneracies; symmetry does not require every pair of modes to be degenerate, and accidental degeneracies can also occur. Within a degenerate eigenspace, orthogonal linear combinations are equally valid mode shapes.
For two equal masses , each attached to a fixed support by a spring of stiffness and coupled to the other by a spring of stiffness ,
The symmetric mode leaves the coupling spring unstretched and has . The antisymmetric mode stretches it and has . Superposing the two gives beats when their frequencies are close.
Continuous systems. The same idea replaces matrices by differential operators and mode vectors by functions satisfying boundary conditions. For a taut string of length , tension and mass per unit length , the wave equation and fixed endpoints give
Orthogonality of the sine functions allows an arbitrary initial displacement to be decomposed into modes. Changing the endpoints changes the spectrum: boundary conditions are part of the physical mode problem.
A uniform chain of equal masses joined by springs , with spacing , has travelling normal modes with
Their quantization gives phonons; the long-wavelength limit is a sound wave. Rectangular electromagnetic cavity modes satisfy Maxwell's equations and conductor boundary conditions, with in vacuum. An elastic beam instead has a fourth-order spatial equation and flexural dispersion . These examples show how restoring forces and geometry determine the frequency spectrum.
Real systems. Small-amplitude motion is often approximately linear even when the full equations are nonlinear. Large amplitudes produce frequency shifts, mode coupling and transfer of energy between modes. Friction broadens resonances and makes free oscillations decay. A damping matrix proportional to or can retain the undamped mode shapes; general damping need not, and the appropriate eigenfrequencies and modes can be complex. External periodic forcing predominantly excites modes near resonance, with amplitude set by coupling to the forcing and dissipation. Imperfections split ideal degeneracies, while radiation or coupling to a surrounding fluid can remove energy from the chosen subsystem. Normal modes remain useful when linearization, the boundary conditions and the treatment of damping accurately describe the motion being measured.
The dynamical method. In Lagrangian mechanics, choose independent generalized coordinates consistent with the constraints, form the Lagrangian , and solve the Euler-Lagrange equations,
This form applies to a Newtonian system with ideal constraints, so constraint forces do no allowed virtual work, and with the remaining forces derivable from the stated potential. Smooth holonomic constraints can be built into the coordinates. Additional nonconservative generalized forces would appear on the right. Prescribed shaft rotation is a driven constraint; its constant is a parameter, not a conserved angular momentum.
Energy of the two balls. The diagram labels each ball , so their combined mass is . It also shows between an arm and the downward vertical, rather than the full angle between the arms. Each ball has horizontal radius , vertical coordinate , and squared speed . Thus
and
The massless connecting rods and sliding collar enforce the geometry but contribute no energy.
Effective potential. Define . At fixed ,
This effective potential determines the reduced dynamics. Its conserved effective energy is not the full lab-frame mechanical energy: the shaft drive can exchange work and angular momentum with the balls.
The change of stability. Near ,
For , the downward state is a minimum. At , its quadratic stiffness vanishes but the quartic coefficient is . Above that speed, the downward state becomes a maximum and stable nonzero openings appear. This is the gravity-restored centrifugal governor bifurcation.
Figure 1.
Original calculated effective potentials below, at and above the critical shaft speed; signed angles display the symmetric mathematical extension of the opening coordinate
. Details.
The physical opening is nonnegative; the two minima in the signed extension represent opposite transverse displacements. The stable opening increases continuously from zero toward a right angle as grows. The apparatus is frictionless in the model, so a displaced system oscillates rather than automatically settling at the equilibrium shown by the minima.
Equilibria. Differentiation gives
Hence the downward equilibrium is , and the open equilibria satisfy
The signed extension includes the negative of this angle. Formally is also stationary, but it is an unstable inverted configuration and lies outside the usual opening range of the mechanism.
Small oscillations. Write . The linear equation is , where
At the downward state and at a nonzero open state, respectively,
These are angular frequencies in ; ordinary frequencies are . Negative indicates exponential instability rather than an oscillation. At , . At the critical speed, the linear frequency is zero and the positive quartic potential gives amplitude-dependent nonlinear restoring motion.
Inertia is a tensor. For a rigid body, distances between material points are fixed. Write each position relative to the centre of mass as , so its rotational velocity is . Integrating gives angular momentum
The rotational kinetic energy is
The inertia tensor is real and symmetric, so orthogonal principal axes diagonalize it. Angular momentum is parallel to angular velocity only along a principal axis, or within a subspace of equal principal moments. For translation plus rotation about the centre of mass, . Moving the reference to a parallel axis adds to the tensor, the tensor form of the parallel-axis theorem.
Euler's equations. The external torque about the centre of mass obeys . In body axes, the transport rule gives
The principal moments are constant in that frame, hence
These are Euler's equations. With zero torque, spatial angular momentum and kinetic energy are conserved, although the body components generally vary. For distinct moments, small perturbations about rotation around the largest or smallest principal moment oscillate, whereas those about the intermediate moment grow. This is the intermediate-axis instability.
A lamina. Choose the lamina's plane as . Then , so the normal is a principal axis, and
This perpendicular-axis theorem applies to a planar mass distribution; it is not a general identity for a three-dimensional body. Rotate the in-plane axes to diagonalize the remaining symmetric two-by-two block. For a uniform rectangular lamina with sides and mass , its central principal moments are , and . If , then : torque-free rotation about the long in-plane axis or the normal is stable, while that about the short in-plane axis is unstable. Coincident moments in a square require the degenerate, rather than strict intermediate-axis, case. Under a constraint to rotate only about the normal, the motion reduces to the scalar equation .
A symmetric top without torque. Let and . Euler's equations reduce to
Thus rotates at signed body-frame frequency , while the axial component stays fixed. The symmetry axis has constant angle to fixed spatial angular momentum. Indeed , so
It precesses in space about with angular speed . This spatial precession and the body-frame frequency describe different motions and should not be identified.
A heavy top. If a symmetric top is pivoted a distance from its centre of mass, gravity gives torque . Here the relevant moments are about the fixed pivot. The axial angular momentum and the vertical angular-momentum component are conserved, because gravity supplies neither an axial torque nor a torque about the vertical. Using tilt , precession angle and spin angle , its Lagrangian is
Eliminating the two conserved cyclic momenta leaves one-dimensional motion in an effective tilt potential; bounded oscillation in tilt is nutation. Steady precession at and fixed tilt satisfies
For a rapidly spinning top on the slow branch, the first term is small and
Gravity then changes the direction of angular momentum rather than simply tipping the top straight down. The approximation fails near conditions where nutation or the fast-precession branch matters; the full equations retain those possibilities.
Lagrangian mechanics expresses dynamics through generalized coordinates that respect the constraints. For a system with conservative, velocity-independent forces and standard kinetic energy, the Lagrangian is . Its equations follow from the principle of stationary action:
for variations that vanish at the endpoints. Varying the integral and integrating the term by parts gives
Independent variations therefore yield the Euler-Lagrange equations:
Additional nonconservative generalized forces appear on the right-hand side as . The action is stationary, not necessarily a minimum. For ideal holonomic constraints, using suitable generalized coordinates removes the need to solve explicitly for constraint forces.
The conjugate momentum is . If a coordinate is cyclic, and its conjugate momentum is conserved. If has no explicit time dependence, the energy function
is conserved. It equals for the usual quadratic kinetic energy with time-independent geometry, but that identification is not universal for every Lagrangian or constraint.
For example, a simple pendulum with angle , length and bob mass has
The tension does no virtual work in the allowed angular displacement and never needs to be calculated to obtain the angular equation. The method scales naturally to coupled systems and coordinates with complicated geometry.
The Lagrangian is . Away from a collision, its equations of motion are
Summing forces proves conservation of total linear momentum. The total mechanical energy is also conserved: its derivative is , because .
For the first proposed angular momentum,
which need not vanish. Thus particle one's angular momentum about the fixed origin is not generally conserved. Its force is central relative to particle two, not necessarily relative to that origin.
For small oscillations about equilibrium, put and . The potential has no linear term at equilibrium; to quadratic order the Lagrangian is . Its equations are . A normal mode satisfies the generalized eigenvalue problem . For a positive-definite mass matrix, diagonalization gives the squared frequencies and independent normal coordinates. A strictly positive Hessian gives positive frequencies; zero eigenvalues require separate higher-order stability analysis.
Here the spring extensions are , , and . Thus
The stiffness matrix is and . Its eigenvectors and have stiffness eigenvalues and , giving
In the first mode the masses move together with equal displacements and the middle spring keeps its length. In the second they move oppositely with equal amplitudes, stretching and compressing the middle spring twice as much as either end spring.
Write , and . The Euler-Lagrange equations give . A normal mode is a motion with fixed nonzero vector ; its normal frequency is . Substitution gives . Thus the squared normal frequencies are the assumed positive distinct roots of
For either root, solve the associated two linear equations for the mode ratio.
For the rods, the small angles are and . The first rod has kinetic energy . The second rod's center has horizontal velocity and its angular velocity is . Its center-of-mass moment of inertia is , giving
Adding these kinetic energies gives . The first center rises by , and the second by . Multiplying by the respective masses and gives
Therefore the stated Lagrangian is .
Set . The mode equation is
Hence
The mode ratio is , so the ratio of the rods' angles is . It is positive for and negative for . Thus the lower-frequency normal mode has the rods in phase, and the higher-frequency normal mode has them out of phase.
The lower mass has coordinates . Its squared speed is . Adding the bead's kinetic energy gives the kinetic part of the supplied Lagrangian. Equal spring prestresses give spring energy up to a constant: the two squared extensions are proportional to . The gravitational potential is . Thus has exactly the stated form, with constant energy offsets discarded.
Since has no explicit time dependence, the conserved energy is
The two Euler–Lagrange equations are
At equilibrium, and . There are two geometrically distinct positions, modulo . The downward one has positive potential Hessian and is stable; the inverted one has a negative angular curvature and is unstable.
Keeping only quadratic terms about the downward equilibrium gives
For unit parameters, the mass and stiffness matrices are and . The equation is . Therefore the normal angular frequencies are
By reflection symmetry, is the midpoint of the cylinder axis. With position measured from , define . Symmetry makes the off-diagonal entries zero. Uniform density gives and , by integrating over a disk and the interval . Hence
For parallel axes based at , the tensor form of the parallel-axis theorem is , where . This follows by substituting in the defining integral and using . Here is an axial vector of length , so
The Euler angles are cyclic coordinates. Their conserved conjugate momenta are and . The Lagrangian has no explicit time dependence, so the conserved energy is kinetic plus gravitational potential, namely the expression given with constant .
Eliminate the angular velocities using and . The energy equation becomes
Thus
The elimination is valid away from the coordinate poles; the multiplied equation extends by continuity to any pole reached by a regular motion.
Take a uniform nearest-neighbour harmonic chain with displacement , conjugate momentum , mass , spring stiffness , lattice spacing and periodic length . These specify a standard model for the requested notes, not additional numerical exam data. Its Lagrangian and Hamiltonian are
Variation gives . For a mode , periodicity requires and
The lattice makes periodic modulo , so choose one Brillouin zone. At long wavelength , where is the sound speed. The continuum limit of a harmonic chain expands and replaces by . Defining mass density and one-dimensional elastic coefficient gives
The momentum density is , and the Hamiltonian density is . Translational invariance of makes total momentum conserved. The continuum approximation describes acoustic wavelengths large compared with ; it is not a justification for extending the lattice dispersion and zero-point sums to arbitrarily large wavevector.[ref]
canonical quantization promotes the coordinates and momenta to operators with and the other equal-time commutators zero. In the continuum, . The nonzero normal modes can be represented by independent phonon operators obeying . A useful lattice expansion is
where is the centre-of-mass displacement. The corresponding momentum is
Including the centre-of-mass term restores the full site commutator: the nonzero-mode sum alone gives . Substitution diagonalizes the Hamiltonian,
A phonon is one quantum of a collective normal mode, carrying energy and lattice momentum modulo a reciprocal-lattice vector. Reality of relates its Fourier coordinates by ; it does not remove the independent oppositely travelling phonon modes. Using one complete Brillouin zone avoids double-counting real standing-wave coordinates.[ref]
The free centre of mass has no restoring force and must not be treated as an oscillator of zero frequency. Fixing it, weakly pinning it, or specifying a finite centre-of-mass volume defines the corresponding zero-mode factor. With a fixed centre of mass, the internal thermal partition function and occupations are
The internal energy is ; its high-temperature heat capacity approaches . For a large one-dimensional chain at temperatures below the lattice cutoff but above its finite-size mode spacing,
This one-dimensional acoustic phonon heat capacity shows that quantization turns classical normal waves into bosonic excitations with a linear low-temperature heat capacity in one dimension. The Euclidean field action is Gaussian, with periodic imaginary-time boundary conditions; its Fourier modes reproduce these oscillators. Absolute displacement correlations require the zero-mode convention, and the lattice regulates short-distance fluctuations.
Keep the electromagnetic units of the supplied Lagrangian, including its factor , and regard as the signed coupling used there. Set . The canonical momentum is , so the Legendre transform gives
Single-valued periodic wavefunctions have normalized eigenstates with integer . Defining dimensionless flux gives , and the trace is therefore
The flux-threaded quantum rotor spectrum is periodic under after relabelling , corresponding to flux period in these units. A local gauge change can remove the magnetic vector potential from the differential operator but then imposes a twisted boundary condition with phase ; it cannot remove the flux dependence from both the Hamiltonian and boundary condition simultaneously. In SI conventions the flux coupling has no , so one must not mix the two unit systems.
Differentiate the Lagrangian with respect to the generalized velocities:
The Legendre transform cancels the vector-potential energy term, giving . Eliminating velocity with therefore yields
For the given gauge, set , , . Hamilton's equations give , , , . Hence the conserved quantities are
The two transverse combinations are components of magnetic pseudomomentum. The printed conserved combinations omit the electric charge factor . With the stated vector potential and general electric charge they need the factors displayed above; as printed they hold only if or if electric charge is already absorbed into the field convention.
The initial data give . Thus , throughout, and
Solving with the initial point gives
Consequently : the guiding centre is the origin and the signed cyclotron frequency fixes the rotation direction. The energy is the constant .
Put and let dots now denote an affine parameter derivative. The null geodesic lagrangian in the orbital plane is . The time and azimuthal Euler-Lagrange equations integrate once; the null normalization gives
For completeness, the differential radial geodesic equation is
the remaining differential equations are and . Spherical symmetry permits the equatorial choice.
Eliminating yields . For , divide by :
The impact parameter is unchanged by a common rescaling of the affine parameter. At infinity it is times angular momentum divided by energy. Radial photons have and obey the undivided equation .
The effective potential is zero at , positive outside, and tends to zero as . Its derivative is . Thus
The circular null geodesic lies at a maximum, since . It is unstable: a small radial displacement away from the maximum permits motion away from the orbit. For incoming photons from infinity in the exterior of a black hole, an outer turning point exists when ; for the photon can cross the event horizon. For a material central body, the vacuum circular orbit is relevant only if lies outside its surface.
Figure 1.
Photon potential and horizon-penetrating light cones
. Details.
For the new coordinates, set , so . The geodesic Lagrangian is . Cyclic-coordinate conservation and the null condition give the required analogous equations:
Here denotes the same conserved time-translation quantity as before. The transformation is . In particular, the displayed time coordinate belongs to Ingoing Kerr-Schild coordinates for Schwarzschild spacetime; the usual advanced null coordinate of ingoing Eddington-Finkelstein coordinates is . Using a metric tensor with a vanishing coefficient while retaining the name would therefore solve a different coordinate representation.
An elimination that stays finite at the horizon is especially simple:
The middle equality uses . Hence the radial energy equation and are unchanged, including the horizon limit. The radial metric tensor block has determinant and is regular at . Its radial null directions are and , as shown in the original schematic.
This coordinate invariance of the Schwarzschild photon radial potential follows because the transformation changes the time slicing but preserves the areal radius and the stationary Killing vector field . Consequently the conserved Killing energy, angular momentum and impact parameter are the same. photon capture, the circular orbit and the radial turning points cannot depend on coordinate labels; coordinate velocities can.
Use an affine parameter , with dots denoting its derivatives, and take so the metric is nondegenerate for . The geodesic Lagrangian is . The Euler-Lagrange equations give
or . The conserved cyclic angular momentum is
This is the cyclic angular momentum in a rotational metric, associated with the rotational Killing vector . Explicitly, by the angular equation.
Differentiating the Lagrangian and substituting both equations gives
Thus the squared tangent norm is constant. If is arc length, it equals one; its numerical value otherwise depends on affine normalization. An arbitrary nonaffine parametrization would not preserve this constant-speed statement.
Define the canonical momentum and the energy function . Along an Euler-Lagrange equation solution,
Thus an autonomous Lagrangian has conserved . For a simple harmonic oscillator, and is constant. For ordinary time-independent mechanical coordinates this is the total mechanical energy.
The orbiting planet produces a potential . It has explicit time dependence and exerts an azimuthal torque, so generally and are nonzero. Nevertheless a simultaneous time-translation symmetry and rotation leaves the Lagrangian and unchanged. Noether's theorem gives and hence
Equivalently implies . This Jacobi integral assumes a rigidly rotating potential with constant and negligible back-reaction from the test particle.
The inertial Lagrangian is . Put and define . The transport identity and invariance of the Euclidean norm under a rotation give
Its expansion contains the ordinary kinetic energy, , and . These velocity and position terms generate the apparent forces in the rotating frame.
For a regular Lagrangian, form the canonical momenta , invert for , and take the Legendre transform . On the upper cone let be distance along a generator from the apex and the azimuth. Then the horizontal radius is and height is , so
Consequently