PhyslibAlpha.ClassicalMechanics.NortonDome.Solution
The motions of the Norton dome and the failure of uniqueness
i. Overview
A particle at rest on the apex of the Norton dome satisfies `r̈ = √r` by staying there forever, and also by staying until an arbitrary instant `T ≥ 0` and then sliding off along `r(t) = (t - T)⁴ / 144`: that curve is `C²`, vanishes with its derivative at `T`, and its second derivative `(t - T)² / 12` is `√r`. All these motions have the same initial position and velocity, so the initial data do not determine the motion. This is Norton's observation.
For contrast, the pendulum's phase-space vector field is Lipschitz and Mathlib's `ODE_solution_unique_univ` gives `SimplePendulum.equationOfMotion_unique`; the force of the dome is not Lipschitz at the apex, `NortonDome.not_lipschitzOnWith_force`. Energy conservation does not help either, since every motion in the family has zero energy (`energy_solution`). Read backwards in time, in section E, a particle sliding up the dome may arrive at the apex at rest in finite time and stay there.
ii. Key results
- `NortonDome.delayedQuartic` is the real function `τ ↦ max (τ - T) 0 ⁴ / 144`, with its derivatives `hasDerivAt_delayedQuartic` and `hasDerivAt_deriv_delayedQuartic`, its `C²` regularity `contDiff_delayedQuartic`, and the identity `sqrt_delayedQuartic`. - `NortonDome.solution T` is the motion leaving the apex at the instant `T`. It is a solution, `solution_isSolution`, starts from rest at the apex when `0 ≤ T`, `solution_zero` and `deriv_solution_zero`, and is injective in `T`, `solution_injective`. - `NortonDome.rest_isSolution` is the motion staying at the apex. - `NortonDome.exists_isSolution_ne` is the failure of uniqueness: two distinct solutions with the same initial position and velocity; `NortonDome.not_isSolution_unique` restates it as the negation of the uniqueness property of the pendulum. All these motions have zero energy, `NortonDome.energy_rest` and `NortonDome.energy_solution`. - `NortonDome.IsSolution.comp_neg` is the time-reversal symmetry of the dome, and `NortonDome.arrival T` the time-reversed motion, off the apex before the instant `-T`, `arrival_apply_pos`, and at rest on it from then on, `arrival_of_le`.
iii. Table of contents
- A. The delayed quartic - A.1. Definition and sign - A.2. Derivatives and regularity - B. The motions leaving the apex - B.1. The definition and its values - B.2. Derivatives and regularity - B.3. The equation of motion - B.4. Distinct instants give distinct motions - C. The motion staying at the apex - D. The failure of uniqueness - E. Time reversal and arrival at the apex
iv. References
- Norton, J. D., *The dome: an unexpectedly simple failure of determinism*, Philosophy of Science 75 (2008), 786–798. - `Physlib.ClassicalMechanics.Pendulum.SimplePendulum.Solution` (uniqueness for a Lipschitz force, for contrast).
A. The delayed quartic
As a function of the real time `τ`, the motion leaving the apex at the instant `T` is `max (τ - T) 0 ⁴ / 144`; its calculus is that of the positive-part powers of `NortonDome.PosPartPow`.
A.1. Definition and sign
A.2. Derivatives and regularity
B. The motions leaving the apex
The motion leaving the apex at the instant `T` is the delayed quartic read on `Time`, times the unit vector of the arc-length coordinate.
B.1. The definition and its values
B.2. Derivatives and regularity
B.3. The equation of motion
B.4. Distinct instants give distinct motions
C. The motion staying at the apex
The force vanishes at the apex, so the constant curve there satisfies the equation of motion.
D. The failure of uniqueness
Rest at the apex and departure at any instant `T ≥ 0` are distinct solutions with the same initial data. So the uniqueness property `SimplePendulum.IsSolution.eq_of_initial` of the pendulum fails for the dome, whose force violates the Picard–Lindelöf hypothesis (`not_lipschitzOnWith_force`).
E. Time reversal and arrival at the apex
The equation of motion has no velocity term, so time reversal maps solutions to solutions. Reversing the departure at the instant `T` gives a motion sliding up the dome, arriving at the apex at the instant `-T` with zero velocity, and staying there. For a locally Lipschitz force this is impossible; for the dome it is the non-uniqueness read backwards.
38 declarations
Delayed quartic function
For a given departure time and an instant , the delayed quartic function is defined as This function is zero for all and follows the power law for . It represents the position of a particle on the Norton dome that remains at rest at the apex until time .
For any departure time and any instant , the delayed quartic function is given by the expression
The Delayed Quartic Function is Non-negative
For any real numbers and , the delayed quartic function is non-negative, i.e., where represents the instant the particle leaves the apex and is the current time.
for
For any departure time and any instant such that , the delayed quartic function , defined as , is equal to .
The delayed quartic is positive for
For any departure time and any instant , if , then the delayed quartic function is strictly positive.
The square root of the delayed quartic function is
For any departure time and any instant , the square root of the delayed quartic function is given by This identity shows that the square root of the position function matches its second derivative , thereby satisfying the equation of motion for a particle on the Norton dome.
The derivative of the delayed quartic is
For any departure time and any instant , the delayed quartic function is differentiable at with derivative
The derivative of the delayed quartic equals
For a given departure time , let the delayed quartic function be defined as . The derivative of this function with respect to time is given by: This represents the velocity of a particle on the Norton dome that begins to slide off the apex at time .
The second derivative of the delayed quartic function is
For any departure time and any instant , the second derivative of the delayed quartic function exists and is given by:
The Delayed Quartic Function is
For any , the delayed quartic function is of class (twice continuously differentiable) on .
Motion on the Norton dome leaving the apex at time
For a given departure time , the function maps an instant to the position of a particle on the Norton dome in a one-dimensional Euclidean space . The position is given by the vector: where is the unit basis vector of . This trajectory describes a particle that remains at rest at the apex (arc length ) until time , and subsequently slides away such that its arc length follows the quartic path .
The arc length of the Norton dome solution equals the delayed quartic
For any departure time and any instant , the arc length of the motion leaving the apex of the Norton dome at time (represented as the -th component of the position vector ) is equal to the delayed quartic function evaluated at :
The Norton dome solution is at the apex for
Let be the trajectory of a particle on the Norton dome that leaves the apex at time . For any time such that , the particle remains at the apex, that is, .
For , the motion starts at the apex
Let denote the trajectory of a particle on the Norton dome that is scheduled to leave the apex at time . For any departure time , the position of the particle at the initial time is at the apex, i.e., .
The particle on the Norton dome is off the apex for
For any departure time and any time , the position of the particle on the Norton dome has a strictly positive arc-length coordinate, i.e., .
The Norton dome motion is
For any departure time , the motion leaving the apex of the Norton dome at time is twice continuously differentiable (i.e., it is a function).
The velocity of the Norton dome motion is
For any departure time , the velocity of the motion leaving the apex of the Norton dome at time is given by the time derivative: where is the standard unit basis vector in the one-dimensional Euclidean space .
For any departure time , the second time derivative (acceleration) of the motion on the Norton dome is given by where is the time and is the unit basis vector in the one-dimensional Euclidean space .
Initial velocity of the Norton dome motion is zero for
For any departure time , the motion on the Norton dome, which describes a particle leaving the apex at instant , has a velocity of zero at the initial time . That is, .
The motion satisfies the Norton dome equation of motion
For any departure time , the motion defined as the trajectory leaving the apex of the Norton dome at instant satisfies the equation of motion . Specifically, its acceleration is equal to the square root of its arc length .
The motion leaving the apex at time is a solution for the Norton dome system.
For any real number , the trajectory defined by (where is the unit basis vector in ) is a solution to the equation of motion for the Norton dome. This represents a particle that remains at rest at the apex until time and then begins to slide down the dome.
for Norton dome motions
Let denote the trajectory of a particle on the Norton dome that leaves the apex at time , defined as where is the unit basis vector. For any two departure times such that , the corresponding trajectories and are distinct.
Injectivity of the mapping from departure time to trajectory on the Norton dome
Let denote the trajectory of a particle on the Norton dome that remains at the apex until time and then slides off, defined by the function: where is the unit basis vector in . The mapping from the departure time to the trajectory is injective; that is, distinct departure times result in distinct physical motions .
for the Norton dome solutions
For any departure time , the motion representing a particle that leaves the apex of the Norton dome at time is not identical to the zero motion (the constant trajectory where the particle remains at rest at the apex).
Rest at the Apex is a Solution for the Norton Dome
The constant trajectory , representing a particle remaining at rest at the apex of the Norton dome, satisfies the equation of motion for the system .
The energy of staying at the apex is zero
For the Norton dome system, the motion that consists of the particle remaining at rest at the apex (represented by the constant zero function) has total energy equal to for all times .
Failure of uniqueness for the Norton dome
There exist two trajectories that are solutions to the equations of motion for the Norton dome system such that they have the same initial position and the same initial velocity , but are distinct trajectories ().
Non-uniqueness of Solutions for the Norton Dome
In the Norton dome system , the uniqueness of solutions given initial conditions does not hold. Specifically, it is not the case that for any two solutions , the conditions and imply that .
Existence of solutions for Norton's dome with arbitrary departure time
For any real number , there exists a trajectory that is a solution to the Norton dome's equation of motion such that the particle starts from rest at the apex at time (i.e., and its time derivative ). This trajectory has the property that the particle remains at the apex for all times (satisfying ) and moves away from the apex for all times (satisfying ).
The Total Energy of the Norton Dome Motion is Zero
For any departure time and any instant , the total energy of the motion (the trajectory where a particle leaves the apex of the Norton dome at time ) is equal to .
The second derivative of is
Let be a normed real vector space and be a function representing a curve. If is differentiable and its time derivative is also differentiable (i.e., is twice differentiable), then for any , the second time derivative of the time-reversed curve evaluated at is equal to the second time derivative of evaluated at . In mathematical notation: where the left-hand side represents the second derivative of the composition of with the negation operation.
Time reversal satisfies the Norton dome equation of motion
Let be a Norton dome system and be a trajectory. If is twice differentiable (specifically, and its time derivative are differentiable) and satisfies the equation of motion for , then the time-reversed trajectory also satisfies the equation of motion. This symmetry arises because the equation of motion for the Norton dome has no velocity-dependent terms.
Time reversal of a Norton dome solution is a solution
Let be a Norton dome system. If a motion is a solution to the equations of motion for , then the time-reversed motion defined by is also a solution for .
Motion on the Norton dome arriving at the apex at time
For a given parameter , this function describes the trajectory of a particle on the Norton dome that arrives at the apex at time . It is defined as the time reversal of the motion leaving the apex at time . The position at time in the one-dimensional Euclidean space is given by: where is the unit basis vector. This motion represents a particle sliding up the dome and coming to a permanent rest at the apex at the instant .
The motion arriving at the apex at time is a solution of the Norton dome equations.
For any real number , the trajectory that arrives at the apex of the Norton dome at time , defined as is a solution to the equations of motion for the Norton dome system .
The motion arriving at the apex at is off the apex for
For any , let denote the trajectory of a particle on the Norton dome that arrives at the apex at time . For any time , the particle is off the apex, meaning its radial position coordinate satisfies .
for
For any real parameter and any time , the position of the particle in the motion is . This means that a particle arriving at the apex of the Norton dome at the instant remains at rest at the apex for all subsequent times.
The velocity of is zero for
For any real parameter , let denote the trajectory of a particle on the Norton dome that arrives at the apex at time . For all times , the velocity of the particle (the time derivative of its position) is zero: This implies that a particle arriving at the apex at time does so with zero velocity and remains at rest for all subsequent times.
