Equations of motion

In mathematical physics, equations of motion are equations that describe the behaviour of a physical system in terms of its motion as a function of time. More specifically, the equations of motion describe the behaviour of a physical system as a set of mathematical functions in terms of dynamic variables: normally spatial coordinates and time are used, but others are also possible, such as momentum components and time. The most general choice are generalized coordinates which can be any convenient variables characteristic of the physical system. The functions are defined in a Euclidean space in classical mechanics, but are replaced by curved spaces in relativity. If the dynamics of a system is known, the equations are the solutions to the differential equations describing the motion of the dynamics.

There are two main descriptions of motion: dynamics and kinematics. Dynamics is general, since momenta, forces and energy of the particles are taken into account. In this instance, sometimes the term refers to the differential equations that the system satisfies (e.g., Newton's second law or Euler–Lagrange equations), and sometimes to the solutions to those equations.

However, kinematics is simpler as it concerns only spatial and time-related variables. In circumstances of constant acceleration, these simpler equations of motion are usually referred to as the "SUVAT" equations, arising from the definitions of kinematic quantities: displacement (S), initial velocity (U), final velocity (V), acceleration (A), and time (T). (see below).

Equations of motion can therefore be grouped under these main classifiers of motion. In all cases, the main types of motion are translations, rotations, oscillations, or any combinations of these.

Historically, equations of motion initiated in classical mechanics and the extension to celestial mechanics, to describe the motion of massive objects. Later they appeared in electrodynamics, when describing the motion of charged particles in electric and magnetic fields. With the advent of general relativity, the classical equations of motion became modified. In all these cases the differential equations were in terms of a function describing the particle's trajectory in terms of space and time coordinates, as influenced by forces or energy transformations. However, the equations of quantum mechanics can also be considered equations of motion, since they are differential equations of the wavefunction, which describes how a quantum state behaves analogously using the space and time coordinates of the particles. There are analogs of equations of motion in other areas of physics, notably waves. These equations are explained below.

Introduction

Qualitative

Equations of motion typically involve:

• a differential equation of motion, usually identified as some physical law and applying definitions of physical quantities, is used to set up an equation for the problem,
• setting the boundary and initial value conditions,
• a function of the position (or momentum) and time variables, describing the dynamics of the system,
• solving the resulting differential equation subject to the boundary and initial conditions.

The differential equation is a general description of the application and may be adjusted appropriately for a specific situation, the solution describes exactly how the system will behave for all times after the initial conditions, and according to the boundary conditions.

Quantitative

In Newtonian mechanics, an equation of motion M takes the general form of a second order ordinary differential equation (ODE) in the position r (see below for details) of the object:

$M\left[\mathbf {r} (t),\mathbf {\dot {r}} (t),\mathbf {\ddot {r}} (t),t\right]=0$

where t is time, and each overdot denotes a time derivative.

The initial conditions are given by the constant values at t = 0:

$\mathbf {r} (0),\quad \mathbf {\dot {r}} (0).$

Another dynamical variable is the momentum p of the object, which can be used instead of r (though less commonly), i.e. a second order ODE in p:

${\tilde {M}}\left[\mathbf {p} (t),\mathbf {\dot {p}} (t),\mathbf {\ddot {p}} (t),t\right]=0$

with initial conditions (again constant values)

$\mathbf {p} (0),\quad \mathbf {\dot {p}} (0).$

The solution r (or p) to the equation of motion, combined with the initial values, describes the system for all times after t = 0. For more than one particle, there are separate equations for each (this is contrary to a statistical ensemble of many particles in statistical mechanics, and a many-particle system in quantum mechanics - where all particles are described by a single probability distribution). Sometimes, the equation will be linear and can be solved exactly. However in general, the equation is non-linear, and may lead to chaotic behaviour depending on how sensitive the system is to the initial conditions.

In the generalized Lagrangian mechanics, the generalized coordinates q (or generalized momenta p) replace the ordinary position (or momentum). Hamiltonian mechanics is slightly different, there are two first order equations in the generalized coordinates and momenta:

$M\left[\mathbf {q} (t),\mathbf {\dot {q}} (t),t\right]=0,\quad {\tilde {M}}\left[\mathbf {p} (t),\mathbf {\dot {p}} (t),t\right]=0$

where q is a tuple of generalized coordinates and similarly p is the tuple of generalized momenta. The initial conditions are similarly defined.

Kinematic equations for one particle

Kinematic quantities

Kinematic quantities of a classical particle of mass m: position r, velocity v, acceleration a.

From the instantaneous position r = r (t), instantaneous meaning at an instant value of time t, the instantaneous velocity v = v (t) and acceleration a = a (t) have the general, coordinate-independent definitions;

$\mathbf {v} ={\frac {d\mathbf {r} }{dt}},\quad \mathbf {a} ={\frac {d\mathbf {v} }{dt}}={\frac {d^{2}\mathbf {r} }{dt^{2}}}\,\!$

Notice that velocity always points in the direction of motion, in other words for a curved path it is the tangent vector. Loosely speaking, first order derivatives are related to tangents of curves. Still for curved paths, the acceleration is directed towards the center of curvature of the path. Again, loosely speaking, second order derivatives are related to curvature.

The rotational analogues are the angular position (angle the particle rotates about some axis) θ = θ(t), angular velocity ω = ω(t), and angular acceleration a = a(t):

${\boldsymbol {\omega }}=\mathbf {\hat {n}} {\frac {d\theta }{dt}},\quad {\boldsymbol {\alpha }}={\frac {d\mathbf {\omega } }{dt}}=\mathbf {\hat {n}} {\frac {d^{2}\theta }{dt^{2}}}\,\!$

where

$\mathbf {\hat {n}} =\mathbf {\hat {e}} _{r}\times \mathbf {\hat {e}} _{\theta }\,\!$

is a unit axial vector, pointing parallel to the axis of rotation, $\scriptstyle \mathbf {\hat {e}} _{r}\,\!$ is the unit vector in direction of r, and $\scriptstyle \mathbf {\hat {e}} _{\theta }\,\!$ is the unit vector tangential to the angle. In these rotational definitions, the angle can be any angle about the specified axis of rotation. It is customary to use θ, but this does not have to be the polar angle used in polar coordinate systems.

The following relations hold for a point-like particle, orbiting about some axis with angular velocity ω:

$\mathbf {v} ={\boldsymbol {\omega }}\times \mathbf {r} \,\!$

$\mathbf {a} ={\boldsymbol {\alpha }}\times \mathbf {r} +{\boldsymbol {\omega }}\times \mathbf {v} \,\!$

where r is a radial position, v the tangential velocity of the particle, and a the particle's acceleration. More generally, these relations hold for each point in a rotating continuum rigid body.

Uniform acceleration

Constant linear acceleration

These equations apply to a particle moving linearly, in three dimensions in a straight line, with constant acceleration. Since the vectors are collinear (parallel, and lie on the same line) - only the magnitudes of the vectors are necessary, hence non-bold letters are used for magnitudes, and because the motion is along a straight line, the problem effectively reduces from three dimensions to one.

Derivation

Two arise from integrating the definitions of velocity and acceleration:

${\begin{aligned}\mathbf {v} &=\int \mathbf {a} dt=\mathbf {a} t+\mathbf {v} _{0}\quad [1]\\\mathbf {r} &=\int (\mathbf {a} t+\mathbf {v} _{0})dt={\frac {\mathbf {a} t^{2}}{2}}+\mathbf {v} _{0}t+\mathbf {r} _{0}\quad [2]\\\end{aligned}}$

in magnitudes:

${\begin{aligned}v&=at+v_{0}\quad [1]\\r&={\frac {{a}t^{2}}{2}}+v_{0}t+r_{0}\quad [2]\\\end{aligned}}$

One is the average velocity - since the velocity increases linearly, the average velocity multiplied by time is the distance travelled while increasing the velocity from v₀ to v (this can be illustrated graphically by plotting velocity against time as a straight line graph):

$\mathbf {r} =\left({\frac {\mathbf {v} +\mathbf {v} _{0}}{2}}\right)t\quad [3]\,\!$

in magnitudes

$r=r_{0}+\left({\frac {v+v_{0}}{2}}\right)t\quad [3]\,\!$

From [3]

$t=\left(r-r_{0}\right)\left({\frac {2}{v+v_{0}}}\right)\,\!$

substituting for t in [1]:

${\begin{aligned}v&=a\left(r-r_{0}\right)\left({\frac {2}{v+v_{0}}}\right)+v_{0}\\v\left(v+v_{0}\right)&=2a\left(r-r_{0}\right)+v_{0}\left(v+v_{0}\right)\\v^{2}+vv_{0}&=2a\left(r-r_{0}\right)+v_{0}v+v_{0}^{2}\\v^{2}&=v_{0}^{2}+2a\left(r-r_{0}\right)\quad [4]\\\end{aligned}}$

From [3]:

$2\left(r-r_{0}\right)-vt=v_{0}t\,\!$

substituting into [2]:

${\begin{aligned}r&={\frac {{a}t^{2}}{2}}+2r-2r_{0}-vt+r_{0}\\0&={\frac {{a}t^{2}}{2}}+r-r_{0}-vt\\r&=r_{0}+vt-{\frac {{a}t^{2}}{2}}\quad [5]\end{aligned}}\,\!$

Usually only the first 4 are needed, the fifth is optional.

${\begin{aligned}v&=at+v_{0}\quad [1]\\r&=r_{0}+v_{0}t+{\frac {{a}t^{2}}{2}}\quad [2]\\r&=r_{0}+\left({\frac {v+v_{0}}{2}}\right)t\quad [3]\\v^{2}&=v_{0}^{2}+2a\left(r-r_{0}\right)\quad [4]\\r&=r_{0}+vt-{\frac {{a}t^{2}}{2}}\quad [5]\\\end{aligned}}$

where r₀ and v₀ are the particle's initial position and velocity, r, v, a are the final position (displacement), velocity and acceleration of the particle after the time interval.

Here a is constant acceleration, or in the case of bodies moving under the influence of gravity, the standard gravity g is used. Note that each of the equations contains four of the five variables, so in this situation it is sufficient to know three out of the five variables to calculate the remaining two.

For non collinear vectors

In the case of the acceleration and velocity vectors not being collinear, it is possible to rewrite the equations of motion in terms of vectors. The derivations are essentially the same as in the scalar case, however it is worth considering the derivation of Torricelli's equation, as this contains vector multiplication.

${\begin{aligned}\mathbf {v} &=\mathbf {a} t+\mathbf {v} _{0}\quad [1]\\\mathbf {r} &=\mathbf {r} _{0}+\mathbf {v} _{0}t+{\frac {{\mathbf {a} }t^{2}}{2}}\quad [2]\\\mathbf {r} &=\mathbf {r} _{0}+\left({\frac {\mathbf {v} +\mathbf {v} _{0}}{2}}\right)t\quad [3]\\v^{2}&=v_{0}^{2}+2\mathbf {a} \cdot \left(\mathbf {r} -\mathbf {r} _{0}\right)\quad [4]\\\mathbf {r} &=\mathbf {r} _{0}+\mathbf {v} t-{\frac {{\mathbf {a} }t^{2}}{2}}\quad [5]\\\end{aligned}}$

Where the Torricelli equation can be derived as follows:

$v^{2}=\mathbf {v} \cdot \mathbf {v} =(\mathbf {v} _{0}+\mathbf {a} t)\cdot (\mathbf {v} _{0}+\mathbf {a} t)=v_{0}^{2}+2t(\mathbf {a} \cdot \mathbf {v} _{0})+a^{2}t^{2}$

$(2\mathbf {a} )\cdot (\mathbf {r} -\mathbf {r} _{0})=(2\mathbf {a} )\cdot (\mathbf {v} _{0}t+{\frac {1}{2}}\mathbf {a} t^{2})=2t(\mathbf {a} \cdot \mathbf {v} _{0})+a^{2}t^{2}=v^{2}-v_{0}^{2}$

$\therefore v^{2}=v_{0}^{2}+2(\mathbf {a} \cdot (\mathbf {r} -\mathbf {r} _{0}))$

SUVAT equations

In elementary physics the above formulae are frequently written as:

${\begin{aligned}v&=u+at\quad [1]\\s&=ut+{\frac {1}{2}}at^{2}\quad [2]\\s&={\frac {1}{2}}(u+v)t\quad [3]\\v^{2}&=u^{2}+2as\quad [4]\\s&=vt-{\frac {1}{2}}at^{2}\quad [5]\\\end{aligned}}$

where u has replaced v₀, s replaces r, and s₀ = 0. They are often referred to as the "SUVAT" equations, where "SUVAT" is an acronym from the variables: s = displacement (s₀ = initial displacement), u = initial velocity, v = final velocity, a = acceleration, t = time.

Applications

Elementary and frequent examples in kinematics involve projectiles, for example a ball thrown upwards into the air. Given initial speed u, one can calculate how high the ball will travel before it begins to fall. The acceleration is local acceleration of gravity g. At this point one must remember that while these quantities appear to be scalars, the direction of displacement, speed and acceleration is important. They could in fact be considered as uni-directional vectors. Choosing s to measure up from the ground, the acceleration a must be in fact −g, since the force of gravity acts downwards and therefore also the acceleration on the ball due to it.

At the highest point, the ball will be at rest: therefore v = 0. Using equation [4] in the set above, we have:

$s={\frac {v^{2}-u^{2}}{-2g}}.$

Substituting and cancelling minus signs gives:

$s={\frac {u^{2}}{2g}}.$

Constant circular acceleration

The analogues of the above equations can be written for rotation. Again these axial vectors must all be parallel (to the axis of rotation), so only the magnitudes of the vectors are necessary:

${\begin{aligned}\omega &=\omega _{0}+\alpha t\\\theta &=\theta _{0}+\omega _{0}t+{\tfrac {1}{2}}\alpha t^{2}\\\theta &=\theta _{0}+{\tfrac {1}{2}}(\omega _{0}+\omega )t\\\omega ^{2}&=\omega _{0}^{2}+2\alpha (\theta -\theta _{0})\\\theta &=\theta _{0}+\omega t-{\tfrac {1}{2}}\alpha t^{2}\\\end{aligned}}\,\!$

where α is the constant angular acceleration, ω is the angular velocity, ω₀ is the initial angular velocity, θ is the angle turned through (angular displacement), θ₀ is the initial angle, and t is the time taken to rotate from the initial state to the final state.

General planar motion

Main article: General planar motion

These are the kinematic equations for a particle traversing a path in a plane, described by position r = r(t). They are actually no more than the time derivatives of the position vector in plane polar coordinates using the definitions of physical quantities (like angular velocity ω).

The position, velocity and acceleration of the particle are respectively:

${\begin{aligned}\mathbf {r} &=\mathbf {r} \left(r(t),\theta (t)\right)=r\mathbf {\hat {e}} _{r}\\\mathbf {v} &=\mathbf {\hat {e}} _{r}{\frac {dr}{dt}}+r\omega \mathbf {\hat {e}} _{\theta }\\\mathbf {a} &=\left({\frac {d^{2}r}{dt^{2}}}-r\omega ^{2}\right)\mathbf {\hat {e}} _{r}+\left(r\alpha +2\omega {\frac {dr}{dt}}\right)\mathbf {\hat {e}} _{\theta }\end{aligned}}\,\!$

where $\scriptstyle \mathbf {\hat {e}} _{r},\mathbf {\hat {e}} _{\theta },\,\!$ are the polar unit vectors. Notice for a the components (–rω²) and 2ωdr/dt are the centripetal and Coriolis accelerations respectively.

Special cases of motion described be these equations are summarized qualitatively in the table below. Two have already been discussed above, in the cases that either the radial components or the angular components are zero, and the non-zero component of motion describes uniform acceleration.

State of motion	Constant r	Linear r	Quadratic r	Non-linear r
Constant θ	Stationary	Uniform translation (constant translational velocity)	Uniform translational acceleration	Non-uniform translation
Linear θ	Uniform angular motion in a circle (constant angular velocity)	Uniform angular motion in a spiral, constant radial velocity	Angular motion in a spiral, constant radial acceleration	Angular motion in a spiral, varying radial acceleration
Quadratic θ	Uniform angular acceleration in a circle	Uniform angular acceleration in a spiral, constant radial velocity	Uniform angular acceleration in a spiral, constant radial acceleration	Uniform angular acceleration in a spiral, varying radial acceleration
Non-linear θ	Non-uniform angular acceleration in a circle	Non-uniform angular acceleration in a spiral, constant radial velocity	Non-uniform angular acceleration in a spiral, constant radial acceleration	Non-uniform angular acceleration in a spiral, varying radial acceleration

General 3d motion

Main article: Spherical coordinate system

In 3d space, the equations become more complicated and unwieldy in spherical coordinates (r, θ, ϕ) with corresponding unit vectors $\scriptstyle \mathbf {\hat {e}} _{r},\mathbf {\hat {e}} _{\theta },\mathbf {\hat {e}} _{\phi }\,\!$, the position, velocity, and acceleration are respectively:

${\begin{aligned}\mathbf {r} &=\mathbf {r} \left(t\right)=r\mathbf {\hat {e}} _{r}\\\mathbf {v} &=v\mathbf {\hat {e}} _{r}+r\,{\frac {d\theta }{dt}}\mathbf {\hat {e}} _{\theta }+r\,{\frac {d\phi }{dt}}\,\sin \theta \mathbf {\hat {e}} _{\phi }\\\mathbf {a} &=\left(a-r\left({\frac {d\theta }{dt}}\right)^{2}-r\left({\frac {d\phi }{dt}}\right)^{2}\sin ^{2}\theta \right)\mathbf {\hat {e}} _{r}\\&+\left(r{\frac {d^{2}\theta }{dt^{2}}}+2v{\frac {d\theta }{dt}}-r\left({\frac {d\phi }{dt}}\right)^{2}\sin \theta \cos \theta \right)\mathbf {\hat {e}} _{\theta }\\&+\left(r{\frac {d^{2}\phi }{dt^{2}}}\,\sin \theta +2v\,{\frac {d\phi }{dt}}\,\sin \theta +2r\,{\frac {d\theta }{dt}}\,{\frac {d\phi }{dt}}\,\cos \theta \right)\mathbf {\hat {e}} _{\phi }\end{aligned}}\,\!$

In the case of a constant ϕ this reduces to the planar equations above.

Harmonic motion of one particle

Translation

The kinematic equation of motion for a simple harmonic oscillator (SHO), oscillating in one dimension (the ±x direction) in a straight line is:

${\frac {d^{2}x}{dt^{2}}}=-\omega ^{2}x$

where ω is the angular frequency of the oscillatory motion, related to the general frequency f and the time period T (time taken for one cycle of oscillation):

$\omega =2\pi f=2\pi /T$

Many systems approximately execute simple harmonic motion (SHM). The complex harmonic oscillator is a superposition of simple harmonic oscillators:

${\frac {d^{2}x}{dt^{2}}}=-\sum _{n}\omega _{n}^{2}x$

It is possible for simple harmonic motions to occur in any direction:

${\frac {d^{2}\mathbf {r} }{dt^{2}}}=-\sum _{n}\omega _{n}^{2}\mathbf {r} _{n}$

known as a multidimensional harmonic oscillator. In cartesian coordinates, each component of the position will be a superposition of sinusiodal SHM.

Rotation

The rotational analogue of SHM in a straight line is angular oscillation about an axle or fulcrum:

${\frac {d^{2}\theta }{dt^{2}}}=-\omega ^{2}\theta$

where ω is still the angular frequency of the oscillatory motion - though not the angular velocity which is the rate of change of θ.

This form can be identified (at least approximately) as libration. The complex analogue is again a superposition of simple harmonic oscillators:

${\frac {d^{2}\theta }{dt^{2}}}=-\sum _{n}\omega _{n}^{2}\theta$

Dynamic equations of motion

Newtonian mechanics

Main article: Newtonian mechanics

It may be simple to write down the equations of motion in vector form using Newton's laws of motion, but the components may vary in complicated ways with spatial coordinates and time, and solving them is not easy. Often there is an excess of variables to solve for the problem completely, so Newton's laws are not the most efficient method for generally finding and solving for the motion of a particle. In simple cases of rectangular geometry, the use of Cartesian coordinates works fine, but other coordinate systems can become dramatically complex.

Newton's second law for translation

The first developed and most famous is Newton's second law of motion, there are several ways to write and use it, the most general is:

$\mathbf {F} ={\frac {d\mathbf {p} }{dt}}$

where p = p(t) is the momentum of the particle and F = F(t) is the resultant external force acting on the particle (not any force the particle exerts) - in each case at time t. The law is also written more famously as:

$\mathbf {F} =m\mathbf {a}$

since m is a constant in Newtonian mechanics. However the momentum form is preferable since this is readily generalized to more complex systems, generalizes to special and general relativity (see four-momentum), and since momentum is a conserved quantity; with deeper fundamental significance than the position vector or its time derivatives.

For a number of particles (see many body problem), the equation of motion for one particle i influenced by other particles is:

${\frac {d\mathbf {p} _{i}}{dt}}=\mathbf {F} _{E}+\sum _{i\neq j}\mathbf {F} _{ij}\,\!$

where p_i = momentum of particle i, F_ij = force on particle i by particle j, and F_E = resultant external force (due to any agent not part of system). Particle i does not exert a force on itself.

Newton's(Euler's) second law for rotation

For rigid bodies, Newton's second law for rotation takes the same form as for translation:

$\mathbf {\tau } ={\frac {d\mathbf {L} }{dt}}\,\!$

where L is the angular momentum. Analogous to force and acceleration:

$\mathbf {\tau } =\mathbf {I} \cdot {\boldsymbol {\alpha }}$

where I is the moment of inertia tensor. Likewise, for a number of particles, the equation of motion for one particle i is:

${\frac {d\mathbf {L} _{i}}{dt}}=\mathbf {\tau } _{E}+\sum _{i\neq j}\mathbf {\tau } _{ij}\,\!$

where L_i = angular momentum of particle i, τ_ij = torque on particle i by particle j, and τ_E = resultant external torque (due to any agent not part of system). Particle i does not exert a torque on itself.

Applications

Some examples of Newton's law include describing the motion of a pendulum:

$-mg\sin \theta =m{\frac {d^{2}(\ell \theta )}{dt^{2}}}\quad \Rightarrow \quad {\frac {d^{2}\theta }{dt^{2}}}=-{\frac {g}{\ell }}\sin \theta \,\!$

a damped, driven harmonic oscillator:

$F_{0}\sin(\omega t)=m\left({\frac {d^{2}x}{dt^{2}}}+2\zeta \omega _{0}{\frac {dx}{dt}}+\omega _{0}^{2}x\right),$

or a ball thrown in the air, in air currents (such as wind) described by a vector field of resistive forces R = R(r, t):

$-{\frac {GmM}{|\mathbf {r} |^{2}}}\mathbf {\hat {e}} _{r}+\mathbf {R} =m{\frac {d^{2}\mathbf {r} }{dt^{2}}}+0\quad \Rightarrow \quad {\frac {d^{2}\mathbf {r} }{dt^{2}}}=-{\frac {GM}{|\mathbf {r} |^{2}}}\mathbf {\hat {e}} _{r}+\mathbf {A} \,\!$

where G = gravitational constant, M = mass of the Earth and A = R/m is the acceleration of the projectile due to the air currents at position r and time t. Newton's law of gravity has been used. The mass m of the ball cancels.

Eulerian mechanics

Euler developed Euler's laws of motion, analogous to Newton's laws, for the motion of rigid bodies.

Newton–Euler equations

The Newton–Euler equations combine Euler's equations into one.

Analytical mechanics

Main articles: Analytical mechanics, Lagrangian mechanics and Hamiltonian mechanics

Constraints and motion

Using all three coordinates of 3d space is unnecessary if there are constraints on the system. Generalized coordinates q(t) = [q₁(t), q₂(t) ... q_N(t)], where N is the total number of degrees of freedom the system has, are any set of coordinates used to define the configuration of the system, in the form of arc lengths or angles. They are a considerable simplification to describe motion since they take advantage of the intrinsic constraints that limit the system's motion - i.e. the number of coordinates is reduced to a minimum, rather than demanding rote algebra to describe the constraints and the motion using all three coordinates.

Corresponding to generalized coordinates are:

• their time derivatives, the generalized velocities: $\mathbf {\dot {q}} =d\mathbf {q} /dt$,
• conjugate "generalized" momenta: $\mathbf {p} =\partial L/\partial \mathbf {\dot {q}} =\partial S/\partial \mathbf {q}$,

(see matrix calculus for the denominator notation) where

• the Lagrangian is a function of the configuration q, the rate of change of configuration dq/dt, and time t: $L=L\left[\mathbf {q} (t),\mathbf {\dot {q}} (t),t\right]$,
• the Hamiltonian is a function of the configuration q, motion p, and time t: $H=H\left[\mathbf {q} (t),\mathbf {p} (t),t\right]$, and
• Hamilton's principal function, also called the classical action is a functional of L: $S[\mathbf {q} ,t]=\int _{t_{1}}^{t_{2}}L(\mathbf {q} ,\mathbf {\dot {q}} ,t)dt$.

The Lagrangian or Hamiltonian function is set up for the system using the q and p variables, then these are inserted into the Euler–Lagrange or Hamilton's equations to obtain differential equations of the system. These are solved for the coordinates and momenta.

Generalized classical equations of motion

As the system evolves, q traces a path through configuration space (only some are shown). The path taken by the system (red) has a stationary action (δS = 0) under small changes in the configuration of the system (δq).

Principle of least action

All classical equations of motion can be derived from this variational principle:

$\delta S=0$

stating the path the system takes through the configuration space is the one with the least action.

Euler–Lagrange equations

The Euler–Lagrange equations are:

${\frac {d}{dt}}\left({\frac {\partial L}{\partial \mathbf {\dot {q}} }}\right)={\frac {\partial L}{\partial \mathbf {q} }}$

After substituting for the Lagrangian, evaluating the partial derivatives, and simplifying, a second order ODE in each q_i is obtained.

Hamilton's equations

Hamilton's equations are:

$\mathbf {\dot {p}} =-{\frac {\partial H}{\partial \mathbf {q} }}\quad \mathbf {\dot {q}} =+{\frac {\partial H}{\partial \mathbf {p} }}$

Notice the equations are symmetric (remain in the same form) by making these interchanges simultaneously:

$\mathbf {p} \rightleftharpoons \mathbf {q} ,\quad H\rightarrow -H.$

After substituting the Hamiltonian, evaluating the partial derivatives, and simplifying, two first order ODEs in q_i and p_i are obtained.

Hamilton–Jacobi equation

Hamilton's formalism can be rewritten as:

$-{\frac {\partial S(\mathbf {q} ,t)}{\partial t}}=H\left(\mathbf {q} ,\mathbf {p} ,t\right)$

Although the equation has a simple form, it's actually a non-linear PDE, first order in N + 1 variables, rather than 2N such equations. Due to the action S, it can be used to identify conserved quantities for mechanical systems, even when the mechanical problem itself cannot be solved fully, because any differentiable symmetry of the action of a physical system has a corresponding conservation law, a theorem due to Emmy Noether.

Electrodynamics

Lorentz force f on a charged particle (of charge q) in motion (instantaneous velocity v). The E field and B field vary in space and time.

In electrodynamics, the force on a charged particle of charge q is the Lorentz force:

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

Combining with Newton's second law gives a first order differential equation of motion, in terms of position of the particle:

$m{\frac {d^{2}\mathbf {r} }{dt^{2}}}=q\left(\mathbf {E} +{\frac {d\mathbf {r} }{dt}}\times \mathbf {B} \right)\,\!$

or its momentum:

${\frac {d\mathbf {p} }{dt}}=q\left(\mathbf {E} +{\frac {\mathbf {p} \times \mathbf {B} }{m}}\right)\,\!$

The same equation can be obtained using the Lagrangian (and applying Lagrange's equations above) for a charged particle of mass m and charge q:

$L={\frac {m}{2}}\mathbf {\dot {r}} \cdot \mathbf {\dot {r}} +q\mathbf {A} \cdot \mathbf {\dot {r}} -q\phi$

where A and ϕ are the electromagnetic scalar and vector potential fields. The Lagrangian indicates an additional detail: the canonical momentum in Lagrangian mechanics is given by:

$\mathbf {P} ={\frac {\partial L}{\partial \mathbf {\dot {r}} }}=m\mathbf {\dot {r}} +q\mathbf {A}$

instead of just mv, implying the motion of a charged particle is fundamentally determined by the mass and charge of the particle. The Lagrangian expression was first used to derive the force equation.

Alternatively the Hamiltonian (and substituting into the equations):

$H={\frac {\left(\mathbf {P} -q\mathbf {A} \right)^{2}}{2m}}-q\phi \,\!$

can derive the Lorentz force equation.

General relativity

Geodesic equation of motion

Geodesics on a sphere are arcs of great circles (yellow curve). On a 2d–manifold (such as the sphere shown), the direction of the accelerating geodesic is uniquely fixed if the separation vector ξ is orthogonal to the "fiducial geodesic" (green curve). As the separation vector ξ₀ changes to ξ after a distance s, the geodesics are not parallel (geodesic deviation).

Main articles: Geodesics in general relativity and Geodesic equation

The above equations are valid in flat spacetime. In curved space spacetime, things become mathematically more complicated since there is no straight line; this is generalized and replaced by a geodesic of the curved spacetime (the shortest length of curve between two points). For curved manifolds with a metric tensor g, the metric provides the notion of arc length (see line element for details), the differential arc length is given by:

$ds={\sqrt {g_{\alpha \beta }dx^{\alpha }dx^{\beta }}}$

and the geodesic equation is a second-order differential equation in the coordinates, the general solution is a family of geodesics:

${\frac {d^{2}x^{\mu }}{ds^{2}}}=-\Gamma ^{\mu }{}_{\alpha \beta }{\frac {dx^{\alpha }}{ds}}{\frac {dx^{\beta }}{ds}}$

where Γ^μ_αβ is a Christoffel symbol of the second kind, which contains the metric (with respect to the coordinate system).

Given the mass-energy distribution provided by the stress–energy tensor T^αβ, the Einstein field equations are a set of non-linear second-order partial differential equations in the metric, and imply the curvature of space time is equivalent to a gravitational field (see principle of equivalence). Mass falling in curved spacetime is equivalent to a mass falling in a gravitational field - because gravity is a fictitious force. The relative acceleration of one geodesic to another in curved spacetime is given by the geodesic deviation equation:

${\frac {d^{2}\xi ^{\alpha }}{ds^{2}}}=-R^{\alpha }{}_{\beta \gamma \delta }{\frac {dx^{\alpha }}{ds}}\xi ^{\gamma }{\frac {dx^{\delta }}{ds}}$

where ξ^α = (x₂)^α − (x₁)^α is the separation vector between two geodesics, D/ds (not just d/ds) is the covariant derivative, and R^α_βγδ is the Riemann curvature tensor, containing the Christoffel symbols. In other words, the geodesic deviation equation is the equation of motion for masses in curved spacetime, analogous to the Lorentz force equation for charges in an electromagnetic field.

For flat spacetime, the metric is a constant tensor so the Christoffel symbols vanish, and the geodesic equation has the solutions of straight lines. This is also the limiting case when masses move according to Newton's law of gravity.

Spinning objects

In general relativity, rotational motion is described by the relativistic angular momentum tensor, including the spin tensor, which enter the equations of motion under covariant derivatives with respect to proper time. The Mathisson–Papapetrou–Dixon equations describe the motion of spinning objects moving in a gravitational field.

Analogues for waves and fields

Unlike the equations of motion for describing particle mechanics, which are systems of coupled ordinary differential equations, the analogous equations governing the dynamics of waves and fields are always partial differential equations, since the waves or fields are functions of space and time. Sometimes in the following contexts, the wave or field equations are also called "equations of motion".

Field equations

Equations that describe the spatial dependence and time evolution of fields are called field equations. These include

• Maxwell's equations for the electromagnetic field,
• Poisson's equation for Newtonian gravitational or electrostatic field potentials,
• the Einstein field equation for gravitation (Newton's law of gravity is a special case for weak gravitational fields and low velocities of particles).

This terminology is not universal: for example although the Navier–Stokes equations govern the velocity field of a fluid, they are not usually called "field equations", since in this context they represent the momentum of the fluid and are called the "momentum equations" instead.

Wave equations

Equations of wave motion are called wave equations. The solutions to a wave equation give the time-evolution and spatial dependence of the amplitude. Boundary conditions determine if the solutions describe traveling waves or standing waves.

From classical equations of motion and field equations; mechanical, gravitational wave, and electromagnetic wave equations can be derived. The general linear wave equation in 3d is:

${\frac {1}{v^{2}}}{\frac {\partial ^{2}X}{\partial t^{2}}}=\nabla ^{2}X$

where X = X(r, t) is any mechanical or electromagnetic field amplitude, say:

• the transverse or longitudinal displacement of a vibrating rod, wire, cable, membrane etc.,
• the fluctuating pressure of a medium, sound pressure,
• the electric fields E or D, or the magnetic fields B or H,
• the voltage V or current I in an alternating current circuit,

and v is the phase velocity. Non-linear equations model the dependence of phase velocity on amplitude, replacing v by v(X). There are other linear and non-linear wave equations for very specific applications, see for example the Korteweg–de Vries equation.

Quantum theory

In quantum theory, the wave and field concepts both appear.

In quantum mechanics, in which particles also have wave-like properties according to wave–particle duality, the analogue of the classical equations of motion (Newton's law, Euler–Lagrange equation, Hamilton–Jacobi equation, etc.) is the Schrödinger equation in its most general form:

$i\hbar {\frac {\partial \Psi }{\partial t}}={\hat {H}}\Psi \,,$

where Ψ is the wavefunction of the system, ${\hat {H}}$ is the quantum Hamiltonian operator, rather than a function as in classical mechanics, and ħ is the Planck constant divided by 2π. Setting up the Hamiltonian and inserting it into the equation results in a wave equation, the solution is the wavefunction as a function of space and time. The Schrödinger equation itself reduces to the Hamilton–Jacobi equation in when one considers the correspondence principle, in the limit that ħ becomes zero.

Applying special relativity to quantum mechanics results in their unification as relativistic quantum mechanics; this is achieved by inserting relativistic Hamiltonians into the Schrödinger equation, leading to relativistic wave equations.

In the context of relativistic and non-relativistic quantum field theory, in which particles are interpreted and treated as fields rather than waves, the Schrödinger equation above has solutions Ψ which are interpreted as fields.

Throughout all aspects of quantum theory, relativistic or non-relativistic, there are various formulations alternative to the Schrödinger equation that govern the time evolution and behavior of a quantum system, for instance:

• the Heisenberg equation of motion resembles the time evolution of classical observables as functions of position, momentum, and time, if one replaces dynamical observables by their quantum operators and the classical Poisson bracket by the commutator,
• the phase space formulation closely follows classical Hamiltonian mechanics, placing position and momentum on equal footing,
• the Feynman path integral formulation extends the principle of least action to quantum mechanics and field theory, placing emphasis on the use of a Lagrangians rather than Hamiltonians.

See also

• Scalar (physics)
• Vector
• Distance
• Displacement
• Speed
• Velocity
• Acceleration
• Angular displacement
• Angular speed
• Angular velocity
• Angular acceleration

• Equations for a falling body
• Parabolic trajectory
• Curvilinear coordinates
• Orthogonal coordinates
• Newton's laws of motion
• Torricelli's equation
• Euler–Lagrange equation
• Generalized forces
• Defining equation (physics)
• Newton–Euler laws of motion for a rigid body

External links

• Equations of Motion Applet