Advanced Classical Mechanics/Many-Particle Systems

Here we will look a group of particles with forces that exert forces among themselves and have external forces exerted upon them as well. If we focus on particle $i$, we have

${\dot{\overrightarrow{p}}}_i=\sum _j{\overrightarrow{F}}_{ji}+F_i^{(e)}$

where the first term is the force that particle $j$ exerts on particle $i$ summed over the other particles. The second term denotes the external force on particle $i$.

Weak law of action and reaction -- total momentum

Let's assume the weak law of action and reaction which states that ${\overrightarrow{F}}_{ij}=-{\overrightarrow{F}}_{ji}$ and sum the change in the momentum of all the particles

$\sum _i{\dot{\overrightarrow{p}}}_i=\frac{d^2}{d{t}^2}\sum _i{m}_i{\overrightarrow{r}}_i=\sum _i{F}_i^{(e)}+\sum _{i,j}{F}_{ji}$

The second term vanishes because each force between each pair of particles cancels out. If we define the position of the center of mass of the system as

$\overrightarrow{R}=\frac{\sum _i{m}_i{\overrightarrow{r}}_i}{\sum _i{m}_i}=\frac{1}{M}\sum _i{m}_i{\overrightarrow{r}}_i$

we have

$M\frac{d^2\overrightarrow{R}}{d{t}^2}=\sum _i{\overrightarrow{F}}_i^{(e)}={\overrightarrow{F}}^{(e)}.$

The total momentum is

$\overrightarrow{P}=\sum _i{m}_i\frac{d{\overrightarrow{r}}_i}{dt}=M\frac{d\overrightarrow{R}}{dt},$

so if ${\overrightarrow{F}}^{(e)}=0$ then $\overrightarrow{P}$ is conserved.

Strong law of action and reaction -- total angular momentum

The total angular momentum of the system is given by

$\overrightarrow{L}=\sum _i{\overrightarrow{r}}_i\times {\overrightarrow{p}}_i$

and its rate of change is

$\dot{\overrightarrow{L}}=\sum _i({\dot{\overrightarrow{r}}}_i\times {\overrightarrow{p}}_i+{\overrightarrow{r}}_i\times {\dot{\overrightarrow{p}}}_i)=\sum _i{\overrightarrow{r}}_i\times {\dot{\overrightarrow{p}}}_i$

where the first term vanishes because the linear momentum and velocity point in the same direction.

Substituting in the forces yields

$\dot{\overrightarrow{L}}=\sum _i{\overrightarrow{r}}_i\times {\overrightarrow{F}}_i^{(e)}+\sum _{i,j}{\overrightarrow{r}}_i\times {\overrightarrow{F}}_{ji}.$

If we look at the final term more closely we have

$\sum _{i,j}{\overrightarrow{r}}_i\times {\overrightarrow{F}}_{ji}=\sum _{i<j}{\overrightarrow{r}}_i\times {\overrightarrow{F}}_{ji}+{\overrightarrow{r}}_j\times {\overrightarrow{F}}_{ij}=\sum _{i<j}({\overrightarrow{r}}_i-{\overrightarrow{r}}_j)\times {\overrightarrow{F}}_{ji},$

so if we assume that the force between the particles points in the same direction as the displacement between the particles (strong law of action and reaction), the forces between the particles do not affect the total angular momentum of the system. This type of force is called a central force. In this case we have $\sum _{i,j}{\overrightarrow{r}}_i\times {\overrightarrow{F}}_{ji}=0$ so

$\dot{\overrightarrow{L}}=\sum _i{\overrightarrow{r}}_i\times {\overrightarrow{F}}_i^{(e)}={\overrightarrow{N}}^{(e)}.$

A central force satisfies the strong law of action and reaction (not necessarily the weak law). Gravity and electric forces are central and also satisfy the weak law of action and reaction but magnetic forces generally satisfy neither.

Let's take a closer look at the total angular momentum of the system of particles by defining the displacement of a particle relative to the center of mass.

${\overrightarrow{r}}_i={\overrightarrow{s}}_i+\overrightarrow{R}\mathrm{a}\mathrm{n}\mathrm{d}{\overrightarrow{\mathrm{v}}}_{\mathrm{i}}\mathrm{=}{\overrightarrow{\mathrm{u}}}_{\mathrm{i}}\mathrm{+}\overrightarrow{\mathrm{V}}\mathrm{w}\mathrm{h}\mathrm{e}\mathrm{r}\mathrm{e}\overrightarrow{\mathrm{V}}\mathrm{=}\frac{d\overrightarrow{R}}{dt}\mathrm{.}$

The total angular momentum is given by

$\overrightarrow{L}=\sum _i{m}_i\overrightarrow{s_i}\times {\overrightarrow{u}}_i+\sum _i{m}_i{\overrightarrow{s}}_i\times \overrightarrow{V}+\sum _i{m}_i\overrightarrow{R}\times {\overrightarrow{v}}_i+\sum _i{m}_i\overrightarrow{R}\times \overrightarrow{V}.$

The two middle terms vanish because

$\sum _i{m}_i{\overrightarrow{s}}_i=\sum _i{m}_i({\overrightarrow{r}}_i-\overrightarrow{R})=M(\overrightarrow{R}-\overrightarrow{R})=0$

from the definition of the center of mass. We are left with

$\overrightarrow{L}=\sum _i{m}_i\overrightarrow{s_i}\times {\overrightarrow{u}}_i+M\overrightarrow{R}\times \overrightarrow{V}.$

That is the angular momentum of a system is the angular momentum of the center of mass plus the angular momentum about the center of mass.

Uniform Gravitational Field

A familar possibility for the external forces and torques is a uniform gravitational field. In this case we have

${\overrightarrow{N}}^{(e)}=\sum _i{m}_i{\overrightarrow{r}}_i\times \overrightarrow{g}=M\overrightarrow{R}\times \overrightarrow{g}.$

Let's calculate the change in the velocity of the center of mass

$\dot{\overrightarrow{V}}=\frac{1}{M}\sum _i{m}_i{\dot{\overrightarrow{v}}}_i=\frac{1}{M}\sum _i{m}_i\overrightarrow{g}=\overrightarrow{g}$

and the change in the angular momentum

$\dot{\overrightarrow{L}}=M\overrightarrow{R}\times \overrightarrow{g}=\frac{d}{dt}(\sum _i{m}_i\overrightarrow{s_i}\times {\overrightarrow{u}}_i)+M(\overrightarrow{V}\times \overrightarrow{V}+\overrightarrow{R}\times \dot{\overrightarrow{V}})=\frac{d}{dt}(\sum _i{m}_i\overrightarrow{s_i}\times {\overrightarrow{u}}_i)+M\overrightarrow{R}\times \overrightarrow{g}$

so

$\frac{d}{dt}(\sum _i{m}_i\overrightarrow{s_i}\times {\overrightarrow{u}}_i)=0.$

In an uniform gravitational field, the internal angular momentum of a body is conserved. Think about footballs and frisbees!

Work and Kinetic Energy

Let's examine the total work performed on the system of particles. We have

$W_{12}=\sum _i\underset{1}{\overset{2}{\int }}{\overrightarrow{F}}_i\cdot d{\overrightarrow{s}}_i=\sum _i\underset{1}{\overset{2}{\int }}{\overrightarrow{F}}_i^{(e)}\cdot d{\overrightarrow{s}}_i+\sum _{i,j}\underset{1}{\overset{2}{\int }}{\overrightarrow{F}}_{ji}\cdot d{\overrightarrow{s}}_i.$

We also have

$W_{12}=\sum _i\underset{1}{\overset{2}{\int }}{\overrightarrow{F}}_i\cdot d{\overrightarrow{s}}_i=\sum _i\underset{1}{\overset{2}{\int }}{m}_i{\dot{\overrightarrow{v}}}_i\cdot {\overrightarrow{v}}_{i}dt=\sum _i\underset{1}{\overset{2}{\int }}d\left(\frac{1}{2}{m}_i{v}_i^2\right)$

so

$W_{12}=T_2-T_1\mathrm{w}\mathrm{h}\mathrm{e}\mathrm{r}\mathrm{e}T=\frac{1}{2}\sum _i{m}_i{v}_i^2.$

We can also look at the energy using the center of mass to get

$T=\frac{1}{2}\sum _i{m}_i{({\overrightarrow{u}}_i+\overrightarrow{V})}^2=\frac{1}{2}\sum m_i{\overrightarrow{u}}_i^2+\sum m_i{u}_i\cdot \overrightarrow{V}+\frac{1}{2}M{V}^2.$

The middle term vanishes for the reasons outlined earlier.

Potential Energy

First let's look at the external forces

$\sum _i\underset{1}{\overset{2}{\int }}{\overrightarrow{F}}_i^{(e)}\cdot d{\overrightarrow{s}}_i=-\sum _i\underset{1}{\overset{2}{\int }}{\overrightarrow{\mathrm{\nabla }}}_i{V}_i\cdot d{\overrightarrow{s}}_i=-\sum _i{V_i|}_1^2$

where ${\mathrm{\nabla }}_i\equiv \frac{d}{d{\overrightarrow{r}}_i}$ if the external forces are conservative. N.B.: in this section, $V$ denotes potential energy not the velocity of the center of mass.

To look at the interparticle forces, let's assume that they too are conservative and that the potential only depends on the distance between the particles so

$V_{ij}=V_{ij}\left(|{\overrightarrow{r}}_i-{\overrightarrow{r}}_j|\right).$

Such a potential has two nice properties. First,

${\overrightarrow{F}}_{ji}=-{\overrightarrow{\mathrm{\nabla }}}_i{V}_{ij}={\overrightarrow{\mathrm{\nabla }}}_j{V}_{ij}=-{\overrightarrow{F}}_{ij}$

so it satisfies the weak law of action-reaction. Second we have

${\overrightarrow{\mathrm{\nabla }}}_i{V}_{ij}\left(|{\overrightarrow{r}}_i-{\overrightarrow{r}}_j|\right)=({\overrightarrow{r}}_i-{\overrightarrow{r}}_j)f\left(|{\overrightarrow{r}}_i-{\overrightarrow{r}}_j|\right)$

so it is central and satifies the strong law of action-reaction. If we use this particular type of force we have

$\sum _{i,j}\underset{1}{\overset{2}{\int }}{\overrightarrow{F}}_{ji}\cdot d{\overrightarrow{r}}_i=-\sum _{i<j}\underset{1}{\overset{2}{\int }}({\overrightarrow{\mathrm{\nabla }}}_i{V}_{ij}\cdot d{\overrightarrow{r}}_i+{\mathrm{\nabla }}_j{V}_{ij}\cdot d{\overrightarrow{r}}_j)$

Let's define ${\overrightarrow{r}}_{ij}={\overrightarrow{r}}_i-{\overrightarrow{r}}_j$. We have

$\frac{dV}{d{\overrightarrow{r}}_{ij}}=\frac{dV}{d{\overrightarrow{r}}_i}={\mathrm{\nabla }}_{i}V=-{\mathrm{\nabla }}_{j}V$

and

$d{\overrightarrow{r}}_{ij}=d({\overrightarrow{r}}_i-{\overrightarrow{r}}_j)=d{\overrightarrow{r}}_i-d{\overrightarrow{r}}_j.$

If we substitute these results into the expression for the work we have

$-\sum _{i<j}\underset{1}{\overset{2}{\int }}({\overrightarrow{\mathrm{\nabla }}}_i{V}_{ij}\cdot d{\overrightarrow{r}}_i+{\mathrm{\nabla }}_j{V}_{ij}\cdot d{\overrightarrow{r}}_j)=-\sum _{i<j}\underset{1}{\overset{2}{\int }}({\overrightarrow{\mathrm{\nabla }}}_{ij}{V}_{ij}\cdot d{r}_{ij})=-{\sum _{i<j}{V}_{ij}|}_1^2=-\frac{1}{2}{\sum _{i\ne j}{V}_{ij}|}_1^2$

so we can definte the total potential energy to be the sum of the potential energies of each particle due to the external force plus the sum of the potential energies of each pair of particles:

$V=\sum _i{V}_i+\frac{1}{2}\sum _{i\ne j}{V}_{ij}.$

The Virial Theorem

Let's look at particular and common form of central force. In particular, let's assume that

$V_{ij}({\overrightarrow{r}}_i,{\overrightarrow{r}}_j)=k{\left|{\overrightarrow{r}}_{ij}\right|}^n.$

and that there is no external potential. This means than

${\mathrm{\nabla }}_i{V}_{ij}=n{V}_{ij}\frac{{\overrightarrow{r}}_{ij}}{{\left|{\overrightarrow{r}}_{ij}\right|}^2}.$

We can define

$K=\sum _i\frac{1}{2}{m}_i{\overrightarrow{r}}_i^2$

as a measure of the size of the system. Let's calculate the second derivative of $K$ with respect to time

$\frac{d^{2}K}{d{t}^2}=\sum _i(m_i{\dot{\overrightarrow{r}}}_i^2+m_i{\overrightarrow{r}}_i\cdot {\ddot{\overrightarrow{r}}}_i)=2T-\sum _i{\overrightarrow{r}}_i\cdot \sum _j{\mathrm{\nabla }}_i{V}_{ij}.$

Again let's group the particles in the final summation in pairs

$\frac{d^{2}K}{d{t}^2}=2T-\sum _{i<j}({\overrightarrow{r}}_i-{\overrightarrow{r}}_j){\mathrm{\nabla }}_i{V}_{ij}=2T-\sum _{i<j}({\overrightarrow{r}}_i-{\overrightarrow{r}}_j)n{V}_{ij}\frac{{\overrightarrow{r}}_i-{\overrightarrow{r}}_j}{{|{\overrightarrow{r}}_i-{\overrightarrow{r}}_j|}^2}=2T-nV.$

If we have a bunch of particles in equilibrium so they are not expanded or contracting on average than the left-hand side of the equation above vanishes. The most familar example of this is the force of gravity with $n=-1$, so for a self-gravitating system in equilibrium we have $2T+V=0$. This relationship is used to weigh distant galaxies and clusters of galaxies.

If $n<0$ the potential energy of the system vanishes as the particles get further apart, so the total energy of the system $E=T+V$ is the negative of the amount of energy needed to break up the system and leave the particles stationary. Let's calculate the total energy of a system in equilibrium in terms of its kinetic energy. We have

$E=T+V=(1+\frac{2}{n})T$

We find that only if $-2<n<0$ is the energy of the system negative. For $n=-2$ the total energy of the equilibrium system is zero, so the equilibrium is marginally stable; even a tiny bit of additional energy is enough to break it apart. For $n<-2$ the equilibrium is unstable. The grouping of particles has more energy than if they separate.

Thankfully the force of gravity has $n=-1$. If $n$ were less than -2, there would be no self-gravitating equilibrium systems such as stars, galaxies etc. In fact theories of gravity in higher dimensions find that $n=2-D$ where $D$ is the number of spatial dimensions. Apparently $D=3$ is kind of special. Template:Subpage navbar