# 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 {\vec {p}}}_{i}=\sum _{j}{\vec {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 ${\vec {F}}_{ij}=-{\vec {F}}_{ji}$ and sum the change in the momentum of all the particles
$\sum _{i}{\dot {\vec {p}}}_{i}={\frac {d^{2}}{dt^{2}}}\sum _{i}m_{i}{\vec {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
${\vec {R}}={\frac {\sum _{i}m_{i}{\vec {r}}_{i}}{\sum _{i}m_{i}}}={\frac {1}{M}}\sum _{i}m_{i}{\vec {r}}_{i}$
we have
$M{\frac {d^{2}{\vec {R}}}{dt^{2}}}=\sum _{i}{\vec {F}}_{i}^{(e)}={\vec {F}}^{(e)}.$
The total momentum is
${\vec {P}}=\sum _{i}m_{i}{\frac {d{\vec {r}}_{i}}{dt}}=M{\frac {d{\vec {R}}}{dt}},$
so if ${\vec {F}}^{(e)}=0$ then ${\vec {P}}$ is conserved.

## Strong law of action and reaction -- total angular momentum

The total angular momentum of the system is given by
${\vec {L}}=\sum _{i}{\vec {r}}_{i}\times {\vec {p}}_{i}$
and its rate of change is
${\dot {\vec {L}}}=\sum _{i}\left({\dot {\vec {r}}}_{i}\times {\vec {p}}_{i}+{\vec {r}}_{i}\times {\dot {\vec {p}}}_{i}\right)=\sum _{i}{\vec {r}}_{i}\times {\dot {\vec {p}}}_{i}$
where the first term vanishes because the linear momentum and velocity point in the same direction.
Substituting in the forces yields
${\dot {\vec {L}}}=\sum _{i}{\vec {r}}_{i}\times {\vec {F}}_{i}^{(e)}+\sum _{i,j}{\vec {r}}_{i}\times {\vec {F}}_{ji}.$
If we look at the final term more closely we have
$\sum _{i,j}{\vec {r}}_{i}\times {\vec {F}}_{ji}=\sum _{i<j}{\vec {r}}_{i}\times {\vec {F}}_{ji}+{\vec {r}}_{j}\times {\vec {F}}_{ij}=\sum _{i<j}\left({\vec {r}}_{i}-{\vec {r}}_{j}\right)\times {\vec {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}{\vec {r}}_{i}\times {\vec {F}}_{ji}=0$ so
${\dot {\vec {L}}}=\sum _{i}{\vec {r}}_{i}\times {\vec {F}}_{i}^{(e)}={\vec {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.
${\vec {r}}_{i}={\vec {s}}_{i}+{\vec {R}}~{\rm {{and}~{\vec {v}}_{i}={\vec {u}}_{i}+{\vec {V}}~{\rm {{where}~{\vec {V}}={\frac {d{\vec {R}}}{dt}}.}}}}$
The total angular momentum is given by
${\vec {L}}=\sum _{i}m_{i}{\vec {s_{i}}}\times {\vec {u}}_{i}+\sum _{i}m_{i}{\vec {s}}_{i}\times {\vec {V}}+\sum _{i}m_{i}{\vec {R}}\times {\vec {v}}_{i}+\sum _{i}m_{i}{\vec {R}}\times {\vec {V}}.$
The two middle terms vanish because
$\sum _{i}m_{i}{\vec {s}}_{i}=\sum _{i}m_{i}\left({\vec {r}}_{i}-{\vec {R}}\right)=M\left({\vec {R}}-{\vec {R}}\right)=0$
from the definition of the center of mass. We are left with
${\vec {L}}=\sum _{i}m_{i}{\vec {s_{i}}}\times {\vec {u}}_{i}+M{\vec {R}}\times {\vec {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
${\vec {N}}^{(e)}=\sum _{i}m_{i}{\vec {r}}_{i}\times {\vec {g}}=M{\vec {R}}\times {\vec {g}}.$
Let's calculate the change in the velocity of the center of mass
${\dot {\vec {V}}}={\frac {1}{M}}\sum _{i}m_{i}{\dot {\vec {v}}}_{i}={\frac {1}{M}}\sum _{i}m_{i}{\vec {g}}={\vec {g}}$
and the change in the angular momentum
${\dot {\vec {L}}}=M{\vec {R}}\times {\vec {g}}={\frac {d}{dt}}\left(\sum _{i}m_{i}{\vec {s_{i}}}\times {\vec {u}}_{i}\right)+M\left({\vec {V}}\times {\vec {V}}+{\vec {R}}\times {\dot {\vec {V}}}\right)={\frac {d}{dt}}\left(\sum _{i}m_{i}{\vec {s_{i}}}\times {\vec {u}}_{i}\right)+M{\vec {R}}\times {\vec {g}}$
so
${\frac {d}{dt}}\left(\sum _{i}m_{i}{\vec {s_{i}}}\times {\vec {u}}_{i}\right)=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}\int _{1}^{2}{\vec {F}}_{i}\cdot d{\vec {s}}_{i}=\sum _{i}\int _{1}^{2}{\vec {F}}_{i}^{(e)}\cdot d{\vec {s}}_{i}+\sum _{i,j}\int _{1}^{2}{\vec {F}}_{ji}\cdot d{\vec {s}}_{i}.$
We also have
$W_{12}=\sum _{i}\int _{1}^{2}{\vec {F}}_{i}\cdot d{\vec {s}}_{i}=\sum _{i}\int _{1}^{2}m_{i}{\dot {\vec {v}}}_{i}\cdot {\vec {v}}_{i}dt=\sum _{i}\int _{1}^{2}d\left({\frac {1}{2}}m_{i}v_{i}^{2}\right)$
so
$W_{12}=T_{2}-T_{1}~{\rm {where}}~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}\left({\vec {u}}_{i}+{\vec {V}}\right)^{2}={\frac {1}{2}}\sum m_{i}{\vec {u}}_{i}^{2}+\sum m_{i}u_{i}\cdot {\vec {V}}+{\frac {1}{2}}MV^{2}.$
The middle term vanishes for the reasons outlined earlier.

## Potential Energy

First let's look at the external forces
$\sum _{i}\int _{1}^{2}{\vec {F}}_{i}^{(e)}\cdot d{\vec {s}}_{i}=-\sum _{i}\int _{1}^{2}{\vec {\nabla }}_{i}V_{i}\cdot d{\vec {s}}_{i}=-\sum _{i}\left.V_{i}\right|_{1}^{2}$
where $\nabla _{i}\equiv {\frac {d}{d{\vec {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(\left|{\vec {r}}_{i}-{\vec {r}}_{j}\right|\right).$
Such a potential has two nice properties. First,
${\vec {F}}_{ji}=-{\vec {\nabla }}_{i}V_{ij}={\vec {\nabla }}_{j}V_{ij}=-{\vec {F}}_{ij}$
so it satisfies the weak law of action-reaction. Second we have
${\vec {\nabla }}_{i}V_{ij}\left(\left|{\vec {r}}_{i}-{\vec {r}}_{j}\right|\right)=\left({\vec {r}}_{i}-{\vec {r}}_{j}\right)f\left(\left|{\vec {r}}_{i}-{\vec {r}}_{j}\right|\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}\int _{1}^{2}{\vec {F}}_{ji}\cdot d{\vec {r}}_{i}=-\sum _{i<j}\int _{1}^{2}\left({\vec {\nabla }}_{i}V_{ij}\cdot d{\vec {r}}_{i}+\nabla _{j}V_{ij}\cdot d{\vec {r}}_{j}\right)$
Let's define ${\vec {r}}_{ij}={\vec {r}}_{i}-{\vec {r}}_{j}$. We have
${\frac {dV}{d{\vec {r}}_{ij}}}={\frac {dV}{d{\vec {r}}_{i}}}=\nabla _{i}V=-\nabla _{j}V$
and
$d{\vec {r}}_{ij}=d\left({\vec {r}}_{i}-{\vec {r}}_{j}\right)=d{\vec {r}}_{i}-d{\vec {r}}_{j}.$
If we substitute these results into the expression for the work we have
$-\sum _{i<j}\int _{1}^{2}\left({\vec {\nabla }}_{i}V_{ij}\cdot d{\vec {r}}_{i}+\nabla _{j}V_{ij}\cdot d{\vec {r}}_{j}\right)=-\sum _{i<j}\int _{1}^{2}\left({\vec {\nabla }}_{ij}V_{ij}\cdot dr_{ij}\right)=-\left.\sum _{i<j}V_{ij}\right|_{1}^{2}=-{\frac {1}{2}}\left.\sum _{i\neq j}V_{ij}\right|_{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\neq 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}({\vec {r}}_{i},{\vec {r}}_{j})=k\left|{\vec {r}}_{ij}\right|^{n}.$
and that there is no external potential. This means than
$\nabla _{i}V_{ij}=nV_{ij}{\frac {{\vec {r}}_{ij}}{\left|{\vec {r}}_{ij}\right|^{2}}}.$
We can define
$K=\sum _{i}{\frac {1}{2}}m_{i}{\vec {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}{dt^{2}}}=\sum _{i}\left(m_{i}{\dot {\vec {r}}}_{i}^{2}+m_{i}{\vec {r}}_{i}\cdot {\ddot {\vec {r}}}_{i}\right)=2T-\sum _{i}{\vec {r}}_{i}\cdot \sum _{j}\nabla _{i}V_{ij}.$
Again let's group the particles in the final summation in pairs
${\frac {d^{2}K}{dt^{2}}}=2T-\sum _{i<j}\left({\vec {r}}_{i}-{\vec {r}}_{j}\right)\nabla _{i}V_{ij}=2T-\sum _{i<j}\left({\vec {r}}_{i}-{\vec {r}}_{j}\right)nV_{ij}{\frac {{\vec {r}}_{i}-{\vec {r}}_{j}}{\left|{\vec {r}}_{i}-{\vec {r}}_{j}\right|^{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=\left(1+{\frac {2}{n}}\right)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.