# PlanetPhysics/Total Energy of a System of Particles

Let us multiply the equation of motion of the k th particle scalarly with ${\frac {d\mathbf {r} }{dt}}$, and sum over all the particles. Then
   $\sum _{k}m_{k}{\frac {d^{2}\mathbf {r} _{k}}{dt^{2}}}{\frac {d\mathbf {r} _{k}}{dt}}={\frac {d}{dt}}{\frac {1}{2}}\sum _{k}m_{k}\left({\frac {d\mathbf {r} _{k}}{dt}}\right)^{2}=\sum _{k}\mathbf {F} _{k}{\frac {d\mathbf {r} _{k}}{dt}}+\sum _{k}\sum _{j}\epsilon _{jk}\mathbf {F} _{jk}{\frac {d\mathbf {r} _{k}}{dt}}$
Integrating between the times $t_{0}$ and $t$:
   ${\frac {1}{2}}\sum _{k}m_{k}\left({\frac {d\mathbf {r} _{k}}{dt}}\right)_{t}^{2}-{\frac {1}{2}}\sum _{k}m_{k}\left({\frac {d\mathbf {r} _{k}}{dt}}\right)_{t_{0}}^{2}=\int _{r_{k}(t_{0})}^{r_{k}(t)}\sum _{k}\mathbf {F} _{k}d\mathbf {r} _{k}+\int _{r_{k}(t_{0})}^{r_{k}(t)}\sum _{k}\sum _{j}\mathbf {F} _{jk}d\mathbf {r} _{k}$
The left member represents the total change in kinetic energy of the system, the right member gives the work done by the internal and external forces. But it is by no means the case that the work done by the internal forces cancels out in calculating the energy, as one might expect it to do. The kinetic energy may be divided into two parts, each of which has a physical meaning. If we introduce a second coordinate system, whose origin $O^{\prime }$ is at the center of gravity of the system, and if we denote all radius vectors referred to this system by primes, we have
$\mathbf {r} _{k}={\bar {\mathbf {r} }}+\mathbf {r} _{k}^{\prime }$
Then, identically,
   $\sum _{k}{\frac {1}{2}}m_{k}\left({\frac {d\mathbf {r} _{k}}{dt}}\right)^{2}={\frac {1}{2}}\left({\frac {d{\bar {\mathbf {r} }}}{dt}}\right)^{2}\sum _{k}m_{k}+{\frac {d{\bar {\mathbf {r} }}}{dt}}\sum _{k}m_{k}{\frac {d\mathbf {r_{k}^{\prime }} }{dt}}+{\frac {1}{2}}\sum _{k}m_{k}\left({\frac {d\mathbf {r_{k}^{\prime }} }{dt}}\right)^{2}$
The second sum on the right vanishes, however, since $\sum m_{k}\mathbf {r} _{k}/M$ is, by equation (3), the radius vector of the center of gravity, and this, by hypothesis, is zero in the primed coordinates. The first term on the right represents the kinetic energy of the system, considering the entire mass to be concentrated at the center of gravity. The last term gives the kinetic energy of motion of the system referred to the center of gravity, when considered at rest. Thus, we may say:
The total kinetic energy is equal to the \htmladdnormallink{translational kinetic energy {http://planetphysics.us/encyclopedia/KineticEnergy.html} of the entire mass, considered concentrated at the center of gravity, plus the energy of motion of the parts of the system relative to the center of gravity}.
We further assume that the internal forces are such that they are derivable from a potential. The potential of the force operating between the points $j$ and $k$ is a function of the distance between the two points, and therefore of their coordinates:
   $U_{jk}=U_{jk}(\mathbf {r} _{jk})=U_{jk}\left({\sqrt {(x_{j}-x_{k})^{2}+(y_{j}-y_{k})^{2}+(z_{j}-z_{k})^{2}}}\right)$
The force acting on $k$ is obtained by taking $j$ to be fixed, and considering $k$ to move in the potential field given by the point function $U_{jk}$; i.e. we consider the coordinates of $j$ to be fixed, those of $k$ to be variable. Then
   $\mathbf {F} _{jk}=-{\hat {i}}{\frac {\partial U_{jk}}{\partial x_{k}}}-{\hat {j}}{\frac {\partial U_{jk}}{\partial y_{k}}}-{\hat {k}}{\frac {\partial U_{jk}}{\partial z_{k}}}=-\nabla _{k}U_{jk}$
in like manner,
   $\mathbf {F} _{kj}=-{\hat {i}}{\frac {\partial U_{jk}}{\partial x_{j}}}-{\hat {j}}{\frac {\partial U_{jk}}{\partial y_{j}}}-{\hat {k}}{\frac {\partial U_{jk}}{\partial z_{j}}}=-\nabla _{j}U_{jk}=-\mathbf {F} _{jk}$
The work done in causing small displacements of $j$ and $k$ is
   $\mathbf {F} _{jk}d\mathbf {r} _{k}+\mathbf {F} _{kj}d\mathbf {r} _{j}=-\left({\frac {\partial U_{jk}}{\partial x_{k}}}dx_{k}+{\frac {\partial U_{jk}}{\partial y_{k}}}dy_{k}+{\frac {\partial U_{jk}}{\partial z_{k}}}dz_{k}+{\frac {\partial U_{jk}}{\partial x_{j}}}dx_{j}+{\frac {\partial U_{jk}}{\partial y_{j}}}dy_{j}+{\frac {\partial U_{jk}}{\partial z_{j}}}dz_{j}\right)=-dU_{jk}$
The negative of the sum of $\mathbf {F} _{jk}d\mathbf {r} _{k}$ and $\mathbf {F} _{kj}d\mathbf {r} _{j}$ is therefore obtained by forming the total differential of $U_{jk}$, defined as a funtion of the six coordinates of the two points, in (11). If, then, we wish to introduce the internal potential into the right member of equation (9), we must write
   $\sum _{k}\sum _{j}\epsilon _{jk}\mathbf {F} _{jk}d\mathbf {r} _{k}=-{\frac {1}{2}}\epsilon _{jk}\sum _{k}\sum _{j}dU_{jk}$
It is readily seen that the factor $1/2$ enters: If we start with point $1$, and calculate the mutual energy $U_{jk}$ between this and all the other points, $k$ runs from $2$ to $N$; but when we take point 2, we must start counting with $3$, since the mutual effect of points $1$ and $2$ was already taken into account in dealing with point $1$, and so on. Thus, in extending the summation over all combinations $j$ and $k$, we must divide by two.
If the external forces have also a potential, the energy equation (9) becomes
   $T+\sum _{k}U_{k}+{\frac {1}{2}}\sum _{k}\sum _{j}\epsilon _{jk}U_{jk}=T^{(0)}+\sum _{k}U_{k}^{(0)}+{\frac {1}{2}}\sum _{k}\sum _{j}\epsilon _{jk}U_{jk}^{(0)}=const.$
where $T$ denotes the kinetic energy. The sum of the kinetic energy and of the external and internal potential energy of a system is constant, if the external as well as the internal forces are conservative .