Derivation of the Lorentz Force Law from the Lagrangian: A Comprehensive Tutorial

Derivation of the Lorentz Force Law from the Lagrangian

Introduction

This tutorial gives a step-by-step derivation of the Lorentz force law from the Lagrangian of a charged particle in electromagnetic potentials. Along the way, we compute the Euler-Lagrange terms carefully, explain the total derivative of the vector potential, and verify the vector identity used in the derivation.

Step 1: Define the Lagrangian

For a charged particle of mass $m$ and charge $q$ moving in electromagnetic fields described by a vector potential $\vec{A}(\vec{x},t)$ and a scalar potential $\phi(\vec{x},t)$, the Lagrangian is

\[L = \frac{1}{2}m\dot{\vec{x}}\cdot\dot{\vec{x}} + q\dot{\vec{x}}\cdot\vec{A} – q\phi .\]

Here $\dot{\vec{x}}$ is the particle velocity, and the dot denotes a time derivative.

Step 2: Compute $\partial L / \partial \dot{\vec{x}}$

The derivative of the Lagrangian with respect to the velocity is

\[\frac{\partial L}{\partial \dot{\vec{x}}} = m\dot{\vec{x}} + q\vec{A}.\]

This quantity is the canonical momentum. It differs from the mechanical momentum $m\dot{\vec{x}}$ by the electromagnetic term $q\vec{A}$.

Step 3: Derive the Vector Identity

We will use the identity

\[\nabla(\vec{U}\cdot\vec{V}) = (\vec{U}\cdot\nabla)\vec{V} + (\vec{V}\cdot\nabla)\vec{U} + \vec{U}\times(\nabla\times\vec{V}) + \vec{V}\times(\nabla\times\vec{U}).\]

In this derivation, take $\vec{U}=\dot{\vec{x}}$ and $\vec{V}=\vec{A}$. When taking $\nabla$ with respect to $\vec{x}$, the velocity $\dot{\vec{x}}$ is treated as independent of position, so

\[(\vec{A}\cdot\nabla)\dot{\vec{x}}=0, \qquad \nabla\times\dot{\vec{x}}=0.\]

Therefore the identity reduces to

\[\nabla(\dot{\vec{x}}\cdot\vec{A}) = (\dot{\vec{x}}\cdot\nabla)\vec{A} + \dot{\vec{x}}\times(\nabla\times\vec{A}).\]

Step 4: Compute $\partial L / \partial \vec{x}$

Now compute the partial derivative of $L$ with respect to $\vec{x}$. The kinetic term has no explicit dependence on $\vec{x}$, so only the potential terms contribute:

\[\frac{\partial L}{\partial \vec{x}} = \nabla\left(q\dot{\vec{x}}\cdot\vec{A} – q\phi\right).\]

The scalar potential term gives

\[\nabla(-q\phi) = -q\nabla\phi.\]

Using the reduced vector identity above, the vector-potential term gives

\[\nabla(q\dot{\vec{x}}\cdot\vec{A}) = q(\dot{\vec{x}}\cdot\nabla)\vec{A} + q\dot{\vec{x}}\times(\nabla\times\vec{A}).\]

Combining the two terms,

\[\frac{\partial L}{\partial \vec{x}} = -q\nabla\phi + q(\dot{\vec{x}}\cdot\nabla)\vec{A} + q\dot{\vec{x}}\times(\nabla\times\vec{A}).\]

Step 5: Use the Euler-Lagrange Equation

The Euler-Lagrange equation is

\[\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{\vec{x}}}\right) – \frac{\partial L}{\partial \vec{x}} = 0.\]

From Step 2,

\[\frac{d}{dt}\left(m\dot{\vec{x}}+q\vec{A}\right) = m\ddot{\vec{x}} + q\frac{d\vec{A}}{dt}.\]

Because $\vec{A}$ depends on both $t$ and the particle position $\vec{x}(t)$, its total time derivative is

\[\frac{d\vec{A}}{dt} = \frac{\partial\vec{A}}{\partial t} + (\dot{\vec{x}}\cdot\nabla)\vec{A}.\]

Therefore,

\[\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{\vec{x}}}\right) = m\ddot{\vec{x}} + q\left(\frac{\partial\vec{A}}{\partial t} + (\dot{\vec{x}}\cdot\nabla)\vec{A}\right).\]

Step 6: Substitute and Simplify

Substitute the two computed terms into the Euler-Lagrange equation:

\[m\ddot{\vec{x}} + q\left(\frac{\partial\vec{A}}{\partial t} + (\dot{\vec{x}}\cdot\nabla)\vec{A}\right) – \left[-q\nabla\phi + q(\dot{\vec{x}}\cdot\nabla)\vec{A} + q\dot{\vec{x}}\times(\nabla\times\vec{A})\right] =0.\]

Expanding and cancelling the repeated $(\dot{\vec{x}}\cdot\nabla)\vec{A}$ terms gives

\[m\ddot{\vec{x}} + q\frac{\partial\vec{A}}{\partial t} + q\nabla\phi – q\dot{\vec{x}}\times(\nabla\times\vec{A}) = 0.\]

Rearrange:

\[m\ddot{\vec{x}} = q\left(-\nabla\phi – \frac{\partial\vec{A}}{\partial t}\right) + q\dot{\vec{x}}\times(\nabla\times\vec{A}).\]

Using the standard definitions

\[\vec{E} = -\nabla\phi – \frac{\partial\vec{A}}{\partial t}, \qquad \vec{B} = \nabla\times\vec{A},\]

we obtain the Lorentz force law:

\[m\ddot{\vec{x}} = q\left(\vec{E}+\dot{\vec{x}}\times\vec{B}\right).\]

Equivalently, since $\vec{F}=m\ddot{\vec{x}}$ and $\vec{v}=\dot{\vec{x}}$,

\[\vec{F}=q(\vec{E}+\vec{v}\times\vec{B}).\]

Supplement: Deriving the Total Derivative

The total derivative of a vector field $\vec{A}(t,\vec{x}(t))$ is defined by

\[\frac{d\vec{A}}{dt} = \lim_{\Delta t\to 0} \frac{\vec{A}(t+\Delta t,\vec{x}(t+\Delta t))-\vec{A}(t,\vec{x}(t))}{\Delta t}.\]

Use the first-order Taylor expansion around $(t,\vec{x}(t))$:

\[\vec{A}(t+\Delta t,\vec{x}(t+\Delta t)) \approx \vec{A}(t,\vec{x}(t)) + \Delta t\frac{\partial\vec{A}}{\partial t} + \Delta\vec{x}\cdot\nabla\vec{A},\]

where

\[\Delta\vec{x}=\vec{x}(t+\Delta t)-\vec{x}(t).\]

Substituting into the limit,

\[\frac{d\vec{A}}{dt} = \lim_{\Delta t\to 0} \left( \frac{\partial\vec{A}}{\partial t} + \frac{\Delta\vec{x}}{\Delta t}\cdot\nabla\vec{A} \right).\]

Since

\[\lim_{\Delta t\to 0}\frac{\Delta\vec{x}}{\Delta t}=\dot{\vec{x}},\]

we get

\[\frac{d\vec{A}}{dt} = \frac{\partial\vec{A}}{\partial t} + (\dot{\vec{x}}\cdot\nabla)\vec{A}.\]

Supplement: Verifying the Vector Identity in 3D

Let

\[\vec{U}=U_x\hat{i}+U_y\hat{j}+U_z\hat{k}, \qquad \vec{V}=V_x\hat{i}+V_y\hat{j}+V_z\hat{k},\]

where each component may depend on $(x,y,z)$. The dot product is

\[\vec{U}\cdot\vec{V}=U_xV_x+U_yV_y+U_zV_z.\]

The $x$-component of $\nabla(\vec{U}\cdot\vec{V})$ is

\[\partial_x(U_xV_x+U_yV_y+U_zV_z) = V_x\partial_xU_x+V_y\partial_xU_y+V_z\partial_xU_z + U_x\partial_xV_x+U_y\partial_xV_y+U_z\partial_xV_z.\]

Now compute the $x$-component of the right-hand side:

\[[(\vec{U}\cdot\nabla)\vec{V}]_x = U_x\partial_xV_x+U_y\partial_yV_x+U_z\partial_zV_x,\]
\[[(\vec{V}\cdot\nabla)\vec{U}]_x = V_x\partial_xU_x+V_y\partial_yU_x+V_z\partial_zU_x,\]
\[[\vec{U}\times(\nabla\times\vec{V})]_x = U_y(\partial_xV_y-\partial_yV_x)-U_z(\partial_zV_x-\partial_xV_z),\]

and

\[[\vec{V}\times(\nabla\times\vec{U})]_x = V_y(\partial_xU_y-\partial_yU_x)-V_z(\partial_zU_x-\partial_xU_z).\]

Adding these four expressions cancels the terms containing $\partial_yV_x$, $\partial_zV_x$, $\partial_yU_x$, and $\partial_zU_x$, leaving exactly

\[\partial_x(U_xV_x+U_yV_y+U_zV_z).\]

The $y$- and $z$-components follow in the same way, so the full vector identity is verified.

Conclusion

Starting from the electromagnetic Lagrangian

\[L = \frac{1}{2}m\dot{\vec{x}}^2 + q\dot{\vec{x}}\cdot\vec{A} – q\phi,\]

the Euler-Lagrange equation gives

\[m\ddot{\vec{x}} = q(\vec{E}+\dot{\vec{x}}\times\vec{B}).\]

The electric force comes from the scalar potential and the explicit time dependence of the vector potential, while the magnetic force comes from the curl of the vector potential.

Leave a Reply