In physics, kinematics describes how the position of a body changes with time, without addressing the forces responsible for that motion. Starting from the position vector $\boldsymbol{x}(t)$, successive time derivatives define velocity, acceleration, and higher-order kinematic quantities. The first derivative,

$$\boldsymbol{v}=\frac{\mathrm{d}\boldsymbol{x}}{\mathrm{d}t},$$

gives the instantaneous rate of change of position and the direction of motion. The second derivative,

$$\boldsymbol{a} =\frac{\mathrm{d}\boldsymbol{v}}{\mathrm{d}t} =\frac{\mathrm{d}^2\boldsymbol{x}}{\mathrm{d}t^2},$$

gives the instantaneous rate of change of velocity, including changes in either speed or direction.

Higher derivatives, however, play a markedly different role. Jerk,

$$\boldsymbol{j} =\frac{\mathrm{d}\boldsymbol{a}}{\mathrm{d}t} =\frac{\mathrm{d}^3\boldsymbol{x}}{\mathrm{d}t^3},$$

is mathematically well defined and important in many applications, yet it rarely appears as an independent state variable or as the highest-order derivative in a fundamental equation of motion. The same is true, to an even greater extent, of still higher kinematic derivatives, such as the fourth time derivative of position, commonly called snap or jounce in engineering. Successive time derivatives of force are also occasionally assigned special names—for example, the first time derivative of force is sometimes called yank—but this terminology is far less standardized than the language of velocity, acceleration, and jerk.

The issue is not whether higher derivatives can be defined, but which variables specify a complete state. In many common models of classical mechanics, positions and velocities—or generalized coordinates and canonical momenta—already form a closed initial-value problem. Acceleration and higher derivatives then follow from the state and the equations of motion. The sections below explain this through initial-value theory, Lagrangian and Hamiltonian mechanics, local field theory, and higher-derivative stability.

1. The Initial-Value Structure of Second-Order Dynamics

One of the central tasks of physics is to predict a system’s subsequent evolution from its present state. For a classical particle of constant mass, Newton’s second law is

$$\boldsymbol{F} =m\boldsymbol{a} =m\frac{\mathrm{d}^2\boldsymbol{x}}{\mathrm{d}t^2}.$$

In many standard models of classical mechanics, the force can be expressed as a function of position, velocity, and time:

$$\boldsymbol{F} =\boldsymbol{F}(\boldsymbol{x},\boldsymbol{v},t).$$

The equation of motion therefore takes the form

$$m\frac{\mathrm{d}^2\boldsymbol{x}}{\mathrm{d}t^2} = \boldsymbol{F}\left( \boldsymbol{x}, \frac{\mathrm{d}\boldsymbol{x}}{\mathrm{d}t}, t \right).$$

For a system with $N$ configurational degrees of freedom, let

$$\boldsymbol{q}\in\mathbb{R}^N, \qquad \boldsymbol{u}=\dot{\boldsymbol{q}}.$$

A second-order system,

$$\ddot{\boldsymbol{q}} = \boldsymbol{f}(t,\boldsymbol{q},\dot{\boldsymbol{q}}),$$

can then be rewritten as the first-order system

$$\frac{\mathrm{d}}{\mathrm{d}t} \begin{pmatrix} \boldsymbol{q}\\ \boldsymbol{u} \end{pmatrix} = \begin{pmatrix} \boldsymbol{u}\\ \boldsymbol{f}(t,\boldsymbol{q},\boldsymbol{u}) \end{pmatrix}.$$

If the vector field on the right-hand side is locally Lipschitz continuous with respect to the state variables $(\boldsymbol{q},\boldsymbol{u})$, together with the usual continuity assumptions, the Picard–Lindelöf theorem implies that specifying an initial time $t_0$ and a complete initial state,

$$\boldsymbol{q}(t_0)=\boldsymbol{q}_0, \qquad \dot{\boldsymbol{q}}(t_0)=\boldsymbol{u}_0,$$

determines a unique local solution near $t_0$. Since the state space is typically $2N$-dimensional, one must specify $2N$ independent initial data.

For a regular Lagrangian system, the same initial state may instead be represented by generalized coordinates and canonical momenta:

$$\bigl(\boldsymbol{q}(t_0),\dot{\boldsymbol{q}}(t_0)\bigr) \quad\longleftrightarrow\quad \bigl(\boldsymbol{q}(t_0),\boldsymbol{p}(t_0)\bigr),$$

provided that the Legendre transformation between velocities and momenta is non-degenerate. These are two parameterizations of the same initial-value problem, not interchangeable versions of an initial-value problem and a boundary-value problem.

A boundary-value problem has a different mathematical structure. For example, prescribing

$$\boldsymbol{q}(t_1)=\boldsymbol{q}_1, \qquad \boldsymbol{q}(t_2)=\boldsymbol{q}_2,$$

does not in general guarantee either existence or uniqueness. Hamilton’s principle does vary the action over paths with fixed endpoints, but this does not imply that arbitrary endpoints determine a unique classical trajectory. There may be no classical solution, or there may be several. Boundary data should therefore not be treated as simply another configuration equivalent to Cauchy initial data.

For the second-order system above, once $\boldsymbol{q}(t_0)$ and $\dot{\boldsymbol{q}}(t_0)$ are specified, the initial acceleration is fixed by the equation of motion:

$$\ddot{\boldsymbol{q}}(t_0) = \boldsymbol{f}\bigl( t_0, \boldsymbol{q}(t_0), \dot{\boldsymbol{q}}(t_0) \bigr).$$

Acceleration is therefore not an additional initial degree of freedom in this class of second-order models, and jerk follows by differentiating the equation of motion along the solution. Once the model and state are specified, higher time derivatives are normally derived quantities rather than independent state variables.

Third- and higher-order equations are mathematically possible. A genuinely third-order equation generally requires additional initial data, corresponding to extra dynamical degrees of freedom or to redundant variables removed by constraints. Such a theory must be judged by its degrees of freedom, stability, causal structure, and empirical predictions, not merely by the existence of solutions.

2. State Space in Lagrangian and Hamiltonian Mechanics

The initial-value argument states what data a second-order equation requires. Analytical mechanics explains why this structure is common. For a Lagrangian that depends on coordinates and velocities,

$$L=L(\boldsymbol{q},\dot{\boldsymbol{q}},t),$$

the Euler–Lagrange equations are

$$\frac{\mathrm{d}}{\mathrm{d}t} \left( \frac{\partial L}{\partial\dot{\boldsymbol{q}}} \right) - \frac{\partial L}{\partial\boldsymbol{q}} =0.$$

When the Hessian of the Lagrangian with respect to the velocities is non-degenerate, these equations can typically be solved for $\ddot{\boldsymbol{q}}$, yielding a second-order system.

Defining the canonical momenta,

$$\boldsymbol{p} = \frac{\partial L}{\partial\dot{\boldsymbol{q}}},$$

and performing a Legendre transformation gives the Hamiltonian,

$$H(\boldsymbol{q},\boldsymbol{p},t) = \boldsymbol{p}\cdot\dot{\boldsymbol{q}}-L.$$

The evolution is governed by Hamilton’s canonical equations:

$$\dot{\boldsymbol{q}} = \frac{\partial H}{\partial\boldsymbol{p}}, \qquad \dot{\boldsymbol{p}} = -\frac{\partial H}{\partial\boldsymbol{q}}.$$

For a regular system with $N$ configurational degrees of freedom, phase space is locally coordinatized by $(\boldsymbol{q},\boldsymbol{p})$, has dimension $2N$, and carries the canonical symplectic form

$$\omega = \sum_{i=1}^{N}\mathrm{d}q^i\wedge\mathrm{d}p_i.$$

Once an initial point in phase space is specified, the Hamiltonian vector field determines the local evolution. Acceleration and jerk are not additional phase-space coordinates; they are derived by repeatedly differentiating the position variables along the Hamiltonian flow.

If a Lagrangian depends explicitly on higher time derivatives, the Ostrogradsky construction introduces additional coordinates and momenta on an enlarged phase space. For a non-degenerate higher-derivative Lagrangian, the Hamiltonian is linear in some canonical momenta and is generally unbounded below. Degenerate theories may avoid this instability through constraints that reduce the number of degrees of freedom.

Second-order equations are common in classical mechanics because a regular Lagrangian depending on coordinates and velocities naturally yields second-order Euler–Lagrange equations and a closed system on a $2N$-dimensional phase space. Higher-derivative theories enlarge that state space and require a separate stability analysis.

The same general idea—that a theory must identify a complete state and then prescribe its evolution—also appears in quantum mechanics. The Schrödinger equation is first order in time:

$$\mathrm{i}\hbar\frac{\partial}{\partial t} \lvert\psi(t)\rangle = \hat{H}\lvert\psi(t)\rangle.$$

Given a quantum state at one time, together with the Hamiltonian operator and its domain, its subsequent unitary evolution is determined in principle. The first-order temporal structure of quantum mechanics is not simply a continuation of classical second-order mechanics, but both theories illustrate the same organizing principle: one first specifies what constitutes a complete state, and the evolution equation then propagates that state.

3. Local Field Theory, Causal Propagation, and the Time Dependence of Force

The preceding sections concern finite-dimensional mechanical systems. Field theory has infinitely many degrees of freedom but follows the same organizing principle: define the complete local state and its evolution, then derive particle forces and higher kinematic quantities. A charged particle, for example, experiences the Lorentz force

$$\boldsymbol{F} = q\left( \boldsymbol{E} + \boldsymbol{v}\times\boldsymbol{B} \right),$$

which depends on the electromagnetic field at the particle’s spacetime position and on its instantaneous velocity. Here, local does not mean that influences propagate instantaneously. The electromagnetic field itself obeys Maxwell’s equations, and changes in the sources propagate causally at finite speed.

With the field retained, particles and fields form a local coupled system. If the field is integrated out, the effective particle equation may depend on its past through delays or memory kernels. The instantaneous variables $\boldsymbol x(t)$ and $\boldsymbol v(t)$ may then be insufficient. The form $\boldsymbol F=\boldsymbol F(\boldsymbol x,\boldsymbol v,t)$ is common, but it is not universal.

When the force does have the form

$$\boldsymbol{F} = \boldsymbol{F}(\boldsymbol{x},\boldsymbol{v},t),$$

its total derivative along the particle’s trajectory is

$$\frac{\mathrm{d}\boldsymbol{F}}{\mathrm{d}t} = \frac{\partial\boldsymbol{F}}{\partial t} + (\boldsymbol{v}\cdot\nabla_{\boldsymbol{x}})\boldsymbol{F} + (\boldsymbol{a}\cdot\nabla_{\boldsymbol{v}})\boldsymbol{F}.$$

The three terms have distinct origins:

• $\partial\boldsymbol{F}/\partial t$ describes the explicit time dependence of the force field;
• $(\boldsymbol{v}\cdot\nabla_{\boldsymbol{x}})\boldsymbol{F}$ describes the change caused by motion through a spatially non-uniform region;
• $(\boldsymbol{a}\cdot\nabla_{\boldsymbol{v}})\boldsymbol{F}$ accounts for the velocity dependence of the force.

Only when the force is independent of velocity, so that $\boldsymbol{F}=\boldsymbol{F}(\boldsymbol{x},t)$, does the last term vanish:

$$\frac{\mathrm{d}\boldsymbol{F}}{\mathrm{d}t} = \frac{\partial\boldsymbol{F}}{\partial t} + (\boldsymbol{v}\cdot\nabla)\boldsymbol{F}.$$

For a particle of constant mass, $\boldsymbol{F}=m\boldsymbol{a}$ implies

$$m\boldsymbol{j} = \frac{\mathrm{d}\boldsymbol{F}}{\mathrm{d}t}.$$

Jerk can therefore be computed from the rate of change of force. This does not imply, however, that $\dot{\boldsymbol{F}}$ is always an algebraic function of the instantaneous position and velocity. Its form depends on which degrees of freedom are retained in the model and on whether the dynamics includes retardation, self-interaction, dissipation, or memory effects.

4. Practical Roles of Higher-Order Derivatives

Jerk is usually not an independent state variable, but it directly measures the smoothness of motion and appears in some effective equations. Engineering control and radiation reaction illustrate these two roles.

4.1 Ride Comfort and Biomechanics

The human body responds differently to sustained and rapidly changing acceleration. Large jerk can produce impact sensations and transient loads even when acceleration remains within an acceptable range. Tolerance also depends on magnitude, direction, duration, and point of application, so jerk is one metric rather than a complete measure of comfort or injury. Trains, elevators, roller coasters, and vehicle controllers commonly limit both acceleration and jerk.

4.2 Robotics, Machine Tools, and Trajectory Planning

In industrial robots, precision machine tools, and motion platforms, discontinuous velocity or acceleration commands can produce impulsive acceleration or very large jerk, thereby exciting flexible structural modes and causing vibration, tracking error, and mechanical wear. Practical motion planning therefore often uses jerk-limited S-curves, and in more demanding applications may also constrain snap.

Several distinct idealizations should be kept separate. A step in acceleration corresponds to a finite but discontinuous acceleration and gives jerk a Dirac-$\delta$ contribution in the idealized description. A step in force means that the force itself is discontinuous. Real systems have finite bandwidth and cannot realize truly infinite rates of change, but the idealized analysis correctly identifies the source of high-frequency excitation.

4.3 Radiation Reaction in Classical Electrodynamics

For a non-relativistic point charge, the Abraham–Lorentz radiation-reaction force is

$$\boldsymbol{F}_{\mathrm{rad}} = \frac{q^2}{6\pi\varepsilon_0c^3} \frac{\mathrm{d}\boldsymbol{a}}{\mathrm{d}t} = \frac{\mu_0q^2}{6\pi c} \frac{\mathrm{d}\boldsymbol{a}}{\mathrm{d}t}.$$

Including this term raises the equation of motion from second to third order. Its general solution requires additional initial data and includes extra branches, among them exponentially growing runaway solutions and pre-acceleration before an applied force begins. The non-relativistic equation is called Abraham–Lorentz; its relativistic extension is Lorentz–Dirac.

These problems arise from the classical point-particle model and do not make every third-order equation unphysical. Treating radiation reaction as a small correction and reducing the order gives approximations such as the Landau–Lifshitz equation. Under the required scale separation, it avoids explicit runaway solutions while preserving the original predictions to the chosen order. Quantum effects become necessary at sufficiently high energy or short distance.

Conclusion

No physical principle forbids jerk or higher time derivatives. They are rarely independent state variables because positions and velocities—or generalized coordinates and momenta—already specify the state of many mechanical systems, and regular Lagrangians depending on coordinates and velocities naturally give second-order equations. Higher derivatives can then be obtained from the state and the dynamics.

Putting higher derivatives into a fundamental equation usually enlarges the state space and can introduce the Ostrogradsky instability in non-degenerate theories. Degenerate theories, delay systems, and radiation reaction require more specific treatment, while jerk remains useful in trajectory design and transient-load analysis. The derivative order is ultimately fixed by what constitutes a complete state and whether the resulting evolution is stable, causal, and consistent with experiment.

Send Feedback