PhyslibAlpha.ClassicalMechanics.NortonDome.Basic
The Norton dome
i. Overview
The Norton dome is a point mass `m` sliding without friction, under gravity `g`, on a rotationally symmetric dome whose surface lies a depth `h(r) = (2/(3g)) r^{3/2}` below its apex at arc length `r`. Its force is continuous, yet a mass at rest on the apex may stay there forever or start sliding down at any later instant: the equation of motion has more than one solution with given initial data. It is the standard example of the failure of determinism in Newtonian mechanics and, mathematically, of the failure of uniqueness for an ODE whose right-hand side is continuous but not Lipschitz.
Newton's laws, as stated, do not determine the motion. A Newtonian system that is to be deterministic must add something the laws do not state, a regularity condition on the force or a restriction on the admissible motions; which addition belongs to Newtonian mechanics is debated, and the folder does not decide it. It proves what each candidate condition buys and that the dome fails each. Reading order:
- `NortonDome.Basic` (this file): the dome, its energies, force, equation of motion `r̈ = √r`, and the failure of the Lipschitz condition at the apex. - `NortonDome.Solution`: rest at the apex and departure at any instant `T` are all solutions from the same initial data, `exists_isSolution_ne`. - `NortonDome.PhysicalSpace`: the arc-length chart as a point mass constrained to the surface. - `NortonDome.NewtonianSystem`: the minimal `NewtonianSystem`; determinism, the first law and the regularity conditions as predicates; a locally Lipschitz force gives determinism. - `NortonDome.Determinism`: the dome is not deterministic, violates the first law, and fails each regularity condition. - `NortonDome.PeanoExistence`, `NortonDome.PosPartPow`, `NortonDome.Sqrt`: supporting analysis. Peano's theorem is pending in Mathlib and marked `@[sorryful]`; all else is proved.
The configuration is the arc length `r`, carried as for the pendulum on the Euclidean lift `Time → EuclideanSpace ℝ (Fin 1)`. The physical dome is `r ≥ 0`; the chart is all of `ℝ`, and for `r < 0` the truncated square root makes the force vanish. With `T = ½ m ṙ²` and `V = -m g h(r) = -(2m/3) r^{3/2}` the equation of motion is `m r̈ = -dV/dr = m √r`, so `r̈ = √r`: `m` and `g` cancel. Norton's profile presumes units in which the constant fixing the size of the dome is `1`.
The force `√r` is continuous but not Lipschitz at the apex, and the potential is `C¹` but not `C²` there; see `NortonDome.Determinism`. So the Picard–Lindelöf hypothesis, which the pendulum satisfies, fails for the dome. For the same reason the Euler–Lagrange theorem `euler_lagrange_varGradient`, which needs a smooth Lagrangian, does not apply, and the equation of motion is identified only with the vanishing of the pointwise Euler–Lagrange operator.
ii. Key results
- `NortonDome` holds the input data, the mass `m` and the gravitational acceleration `g`. - `NortonDome.height` is the depth `(2/(3g)) r^{3/2}` of the dome below its apex, with its derivative `hasDerivAt_height` and `C¹` regularity `height_contDiff_one`. - `NortonDome.kineticEnergy`, `NortonDome.potentialEnergy` and `NortonDome.energy` are the energies, with the gradient `gradient_potentialEnergy` of the potential and the time derivatives `kineticEnergy_deriv`, `potentialEnergy_deriv` and `energy_deriv` along twice differentiable curves. - `NortonDome.lagrangian` is the Lagrangian `T - V`, with its partial gradients `gradient_lagrangian_position_eq` and `gradient_lagrangian_velocity_eq`. - `NortonDome.force` is the generalized force `m √r` conjugate to the arc length. It is continuous, `force_continuous`, but not Lipschitz on any closed ball about the apex, `not_lipschitzOnWith_force`. - `NortonDome.EquationOfMotion` is the equation of motion `m r̈ = F(r)`, with its scalar form `equationOfMotion_iff_scalar`; `NortonDome.IsSolution` is a twice differentiable solution. - `NortonDome.equationOfMotion_iff_eulerLagrangeOp_zero` identifies the equation of motion, for curves with differentiable velocity, with the vanishing of the Euler–Lagrange operator. - `NortonDome.IsSolution.energy_eq` is the conservation of energy along a solution.
iii. Table of contents
- A. The input data - B. The profile of the dome - B.1. The height below the apex - B.2. The derivative of the profile - C. The energies - C.1. The definitions of the energies - C.2. Differentiability and the gradient of the potential - C.3. Time derivatives of the energies - D. The Lagrangian - E. The force and the equation of motion - E.1. The force - E.2. Regularity of the force - E.3. The equation of motion - E.4. Solutions - F. The Euler–Lagrange operator - G. Energy conservation
iv. References
- Norton, J. D., *The dome: an unexpectedly simple failure of determinism*, Philosophy of Science 75 (2008), 786–798. - Malament, D. B., *Norton's slippery slope*, Philosophy of Science 75 (2008), 799–816.
A. The input data
The dome is specified by the mass of the particle and the gravitational acceleration; its shape is fixed by Norton's formula.
B. The profile of the dome
Norton's profile is the depth `h(r) = (2/(3g)) r^{3/2}` of the surface below the apex at arc length `r`, written here as `(2/(3g)) √r ^ 3` so that it is defined, and vanishes, for `r ≤ 0`.
B.1. The height below the apex
B.2. The derivative of the profile
The slope `√r / g` of the profile is continuous and vanishes at the apex, so the profile is `C¹`. The slope itself is not differentiable at the apex, so the profile is not `C²`.
C. The energies
The kinetic energy is `½ m ṙ²`, `r` being the arc length so that the speed is `|ṙ|`, and the potential energy is `-m g h(r)`, normalized to vanish at the apex.
C.1. The definitions of the energies
C.2. Differentiability and the gradient of the potential
The potential energy is differentiable, with gradient `-m √r` times the unit vector of the arc-length coordinate. It is `C¹` but not `C²`, since `√r` is not differentiable at `0`.
C.3. Time derivatives of the energies
Along a twice differentiable curve, solution or not, the time derivatives of the energies are inner products against the velocity; the equation of motion makes the two contributions cancel.
D. The Lagrangian
The Lagrangian `L = ½ m ṙ² + m g h(r)` is a function on phase space. Unlike those of the harmonic oscillator and the pendulum it is only `C¹`; its partial gradients still exist, and are the force of section E and the momentum `m ṙ`.
E. The force and the equation of motion
Gravity exerts the generalized force `m √r` conjugate to the arc length, balanced against `m r̈`. As for the pendulum this pointwise relation is the definition of the equation of motion; the variational derivative of the action is not available, the Lagrangian not being smooth.
E.1. The force
E.2. Regularity of the force
The force is continuous. It is not Lipschitz on any closed ball about the apex, since the square root is not Lipschitz on any `[0, ε]`: the Picard–Lindelöf hypothesis fails, which is the source of the non-uniqueness.
E.3. The equation of motion
E.4. Solutions
A solution is a twice differentiable curve, position and velocity both differentiable, satisfying the equation of motion: the regularity of `ClassicalMechanics.ReferenceFrame.Particle` and the least under which the acceleration exists. Regularity is part of the definition because the bare pointwise equation, whose derivatives are zero wherever they do not exist, admits rough accidental solutions. Smoothness is not demanded: the motions leaving the apex are `C³` but not `C⁴`.
F. The Euler–Lagrange operator
With the gradients of section D, the pointwise Euler–Lagrange operator along any curve with differentiable velocity is the force minus `m r̈`, so its vanishing is the equation of motion. The identification of this operator with the variational derivative of the action, `euler_lagrange_varGradient`, needs a smooth Lagrangian and so does not apply to the dome.
G. Energy conservation
Along any twice differentiable curve satisfying the equation of motion the energy is constant. This follows from `energy_deriv`. Conservation of energy does not restore uniqueness: the motion at rest at the apex and every motion leaving it all have zero energy, which is proved in `NortonDome.Solution` as `energy_rest` and `energy_solution`.
50 declarations
In the Norton dome system, the mass of the particle is non-zero, .
For a Norton dome system, the gravitational acceleration satisfies .
Height of the Norton dome
The function represents the depth of the Norton dome's surface below its apex at an arc length from the apex. It is defined by the formula: where is the gravitational acceleration. The height vanishes at the apex () and for all negative arguments .
Explicit formula for the Norton dome height
For any real number , the depth of the Norton dome's surface below its apex at arc length is given by the formula where is the gravitational acceleration.
Non-negativity of the Norton dome depth
For any real number , the depth of the Norton dome below its apex is non-negative, i.e., . Here, is defined as , where is the gravitational acceleration.
For a Norton dome, the depth of the dome's surface below its apex at arc length is equal to zero if and only if .
The derivative of the dome height is
For a Norton dome with gravitational acceleration , let be the function representing the depth of the surface below the apex at arc length . For any , the function is differentiable at with derivative where is the square root of (defined as for ).
The height function of the Norton dome is differentiable
The height function of the Norton dome, defined by where is the gravitational acceleration and is the arc length from the apex, is differentiable for all .
For the Norton dome, the derivative of the height function with respect to the arc length is given by where is the gravitational acceleration.
The depth of the Norton dome is
Let be the function representing the depth of the Norton dome's surface below its apex at an arc length , defined as where is the gravitational acceleration. Then is a continuously differentiable function (of class ) on .
Kinetic energy
For a particle of mass on the Norton dome, given a trajectory representing the arc length, the kinetic energy is the function of time defined as: where denotes the time derivative (velocity) of the trajectory at time , and denotes the standard inner product on the 1-dimensional Euclidean space.
Potential energy of the Norton dome
The potential energy of a particle of mass on the Norton dome at a position is defined as where is the gravitational acceleration, is the arc length from the apex, and is the depth of the dome's surface below the apex. The potential is normalized to vanish at the apex ().
Energy of a particle on the Norton dome
For a particle on the Norton dome, given a trajectory where represents the position (arc length from the apex), the energy at time is the sum of its kinetic energy and its potential energy :
Kinetic Energy for the Norton Dome
For a particle of mass on a Norton dome, given a trajectory representing the arc length from the apex, the kinetic energy at time is given by: where denotes the time derivative of the trajectory at time , and denotes the standard inner product on the 1-dimensional Euclidean space.
Potential Energy for the Norton Dome
For a particle of mass on a Norton dome, given a position where represents the arc length from the apex, the potential energy is given by:
Total Energy for the Norton Dome
For a particle on the Norton dome with a trajectory representing the arc length from the apex, the total energy at time is the sum of the kinetic energy and the potential energy :
for the Norton dome potential energy
For any position in the 1-dimensional Euclidean space (representing the arc-length coordinate on the dome), the potential energy of the particle is non-positive, i.e., . This reflects the fact that the apex (where ) is the highest point of the dome.
for the Norton dome potential energy
For a particle on the Norton dome at position in the 1-dimensional Euclidean space , the potential energy is equal to zero if and only if the arc-length coordinate satisfies . This condition corresponds to the particle being at or "behind" the apex of the dome.
The function is differentiable everywhere
For every point in the 1-dimensional Euclidean space (representing the arc-length coordinate), the function is differentiable at , where is the coordinate of .
The Potential Energy is Differentiable
The potential energy of a particle on the Norton dome is a differentiable function of the arc-length coordinate .
The gradient of the potential energy of the Norton dome is
For a particle of mass on the Norton dome, let denote its position, where is the arc-length coordinate measured from the apex. The gradient of the potential energy at position is given by where is the unit basis vector of the 1-dimensional arc-length coordinate.
Kinetic Energy is Differentiable along a Curve with Differentiable Velocity
Let be a trajectory representing the arc length of a particle on the Norton dome. If the velocity is differentiable with respect to time, then the kinetic energy is a differentiable function of time.
Potential Energy is Differentiable along a Differentiable Curve
Let be a differentiable curve representing the position (arc length) of a particle on the Norton dome. Then the potential energy is a differentiable function of time .
Total Energy is Differentiable along a Twice Differentiable Curve
Let be a trajectory representing the position (arc length) of a particle on the Norton dome. If is differentiable and its velocity is also differentiable (meaning the curve is twice differentiable), then the total energy is a differentiable function of time.
Time Derivative of Kinetic Energy equals
Let be a curve representing the arc length of a particle of mass on a Norton dome. If the velocity is differentiable, then the rate of change of the kinetic energy is given by the inner product of the velocity and the mass times the acceleration: where denotes the time derivative and is the standard inner product on the one-dimensional Euclidean space.
Time Derivative of Potential Energy equals
Let be a differentiable curve representing the position of a particle on the Norton dome (parameterized by arc length). The rate of change of the potential energy along this curve is given by the inner product of the velocity and the gradient of the potential : where denotes the time derivative and is the standard inner product on the one-dimensional Euclidean space.
Time Derivative of Energy equals
Let be a curve representing the trajectory (arc length from the apex) of a particle of mass on the Norton dome. If is twice differentiable (i.e., both the position and the velocity are differentiable), then the time derivative of the total energy is the inner product of the velocity and the sum of the mass-acceleration term and the gradient of the potential energy : where is the velocity, is the acceleration, and denotes the standard inner product on the one-dimensional Euclidean space.
Lagrangian of the Norton dome
For a particle of mass on the Norton dome, the Lagrangian is a function of time , position (represented as the arc length from the apex), and velocity , defined as the difference between the kinetic energy and the potential energy: where is the square of the norm of the velocity and is the potential energy at position . The Lagrangian is independent of the time .
for the Norton Dome
For any trajectory representing the arc length from the apex of the Norton dome and for any time , the Lagrangian evaluated at position and velocity is equal to the kinetic energy at time minus the potential energy at position : where denotes the time derivative of at .
for the Norton Dome
For a particle of mass on a Norton dome, let be the Lagrangian of the system, where is the time, is the position vector (with representing the arc length from the apex), and is the velocity vector. The gradient of the Lagrangian with respect to the position is given by where is the unit basis vector in . This gradient is equal to the generalized force , where is the potential energy.
for the Norton dome
For a particle of mass on the Norton dome, let be the Lagrangian at time , position , and velocity . The gradient of the Lagrangian with respect to the velocity is equal to the momentum , which is expressed as:
Generalized force of the Norton dome
For a position in the one-dimensional configuration space (representing the arc length from the apex), the generalized force acting on the particle on the Norton dome is defined as the negative gradient of the potential energy , namely .
The generalized force on the Norton dome is
For a Norton dome with mass , the generalized force acting on a particle at a position in the one-dimensional configuration space (where the coordinate represents the arc length from the apex) is given by where is the unit basis vector in . The force points away from the apex and vanishes at .
The component of the generalized force is
For a Norton dome with mass , the -th component of the generalized force at a position is given by where is the coordinate of representing the arc length from the apex.
The Generalized Force Vanishes at the Apex of the Norton Dome
For a Norton dome , the generalized force acting on the particle at the apex () is zero, i.e., .
The Generalized Force is Continuous
The generalized force acting on a particle on the Norton dome is continuous with respect to the arc length .
for the arc-length unit vector
In a one-dimensional Euclidean space, let be the unit vector (represented by `EuclideanSpace.single 0 1`). For any two real numbers and , the distance between the vectors and is equal to the absolute difference between the coefficients, .
The generalized force of the Norton dome is not Lipschitz continuous on any neighborhood of the apex
For any non-negative Lipschitz constant and any radius , the generalized force of the Norton dome is not Lipschitz continuous with constant on the closed ball centered at the apex.
Equation of motion for the Norton dome
A trajectory satisfies the equation of motion for a particle of mass on the Norton dome if, at every instant , the product of the mass and the acceleration equals the generalized force . That is, where represents the arc length coordinate from the apex and denotes the second derivative with respect to time.
The Norton dome equation of motion is equivalent to
For a trajectory representing the arc length of a particle of mass on a Norton dome, the equation of motion is satisfied if and only if at every instant , the sum of the mass multiplied by the acceleration and the gradient of the potential energy vanishes: Here, denotes the second derivative with respect to time, and is the potential energy of the system.
Scalar Equation of Motion for the Norton Dome
For a trajectory representing the arc length from the apex of a Norton dome, the equation of motion is satisfied if and only if for every instant , the second time derivative of the position satisfies the scalar relation: In this form, the mass and the gravitational acceleration have cancelled out.
is a solution for the Norton dome
A trajectory , representing the arc length coordinate from the apex, is a solution for the Norton dome if it satisfies the following conditions: 1. The position function is differentiable with respect to time. 2. The velocity function is differentiable with respect to time (making twice differentiable). 3. The trajectory satisfies the equation of motion for the Norton dome, , where is the mass and is the generalized force.
Solutions for the Norton Dome are Differentiable
For a Norton dome , let be a trajectory representing the arc length coordinate from the apex. If is a solution for the system , then the position function is differentiable with respect to time.
The velocity of a Norton dome solution is differentiable
Let be a Norton dome and be a trajectory representing the arc length from the apex. If is a solution to the equations of motion for the dome, then its velocity is differentiable with respect to time.
Solutions to the Norton dome satisfy the equation of motion
Let be a Norton dome with mass and generalized force . If a trajectory , representing the arc length coordinate from the apex, is a solution for the Norton dome, then it satisfies the equation of motion at every instant , where is the second derivative of with respect to time.
Euler-Lagrange Operator for the Norton Dome equals
Consider a particle of mass on a Norton dome. Let be a curve representing the arc length from the apex, and suppose that the velocity is differentiable. Then the Euler-Lagrange operator applied to the Lagrangian along the trajectory is equal to the generalized force minus the product of the mass and the acceleration : where the Euler-Lagrange operator is defined as .
Equation of Motion Euler-Lagrange Operator Vanishes for the Norton Dome
Let be a Norton dome with mass . For a trajectory representing the arc length from the apex, suppose that the velocity is differentiable. Then the trajectory satisfies the equation of motion if and only if the Euler-Lagrange operator of the Lagrangian vanishes along the curve, that is: where and are the first and second time derivatives of , is the generalized force, and the Lagrangian is defined as .
Vanishing of the Time Derivative of Energy for the Norton Dome Equation of Motion
Let be a trajectory representing the arc length from the apex of a Norton dome . Suppose that is twice differentiable (specifically, that and its velocity are differentiable). If satisfies the equation of motion , then the time derivative of the total energy of the system is zero for all : where is the sum of the kinetic and potential energies at time .
Conservation of Energy for the Norton Dome Equation of Motion:
Consider a particle on a Norton dome . Let be a trajectory representing the arc length from the apex. Suppose that is twice differentiable and satisfies the equation of motion . Then for any time , the total energy of the particle is equal to its initial value .
The energy of a Norton dome solution at any time equals its initial energy
For a Norton dome , let be a trajectory representing the arc length from the apex. If is a solution to the equation of motion (as defined by `IsSolution`), then for any time , the total energy of the particle is equal to its initial energy at time :
