Table of Contents

CM.0 Newtonian Mechanics

CM.0.1 Newton's Laws for a Single Particle

The goal of physics is to predict the future from present conditions. In classical mechanics, we usually want to predict the future position of particles given their current positions and velocities. Newton gave three laws of motion for doing so, which I am here quoting from [1]: (1) An object at rest remains at rest, and an object in motion remains in motion at constant speed and in a straight line unless acted on by an unbalanced force. (2) The acceleration of an object depends on the mass of the object and the amount of force applied. (3) Whenever one object exerts a force on another object, the second object exerts an equal and opposite force on the first, acting in the direction pointing from one object to the other.

We now want to state the mathematical implications of these laws and introduce the notion of conserved quantities, which are fundamental to physics. This Chapter will largely follow the first few pages of Goldstein, [2], since it is the simplest way to present this material.

CM.0.1.1 Conservation of Momentum

Newton's second law gives us our equation of motion for a particle:

𝐅=d𝐩dt,(CM.0.1)

where 𝐅 is the total force on the particle, which can be a function of time or many other variables, 𝐩=m𝐯 is the particle's momentum, m is the particle's mass, 𝐯=d𝐫/dt is the particle's velocity, and 𝐫 is the particle's position (bold quantities are vectors). We will assume that the mass of our particle is constant, but the velocity can be a function of time, 𝐯(t), whose derivative is the acceleration 𝐚=d𝐯/dt. Moving the constant mass outside of the time derivative in Equation CM.0.1 gives the familiar form of Newton's Second Law,

𝐅=md𝐯dt=m𝐚.(CM.0.2)

A conservation law is an equation that states that a certain quantity, call it Q is not a function of time, i.e. dQ/dt=0. From Equation CM.0.1, we can see that if the force 𝐅=0, then the momentum of our particle is conserved, that is it is constant in time, d𝐩/dt=0.

CM.0.1.2 Conservation of Angular Momentum

The angular momentum 𝐋 of a particle, with respect to the origin of coordinates being used to define the position 𝐫, is defined as

𝐋=𝐫×𝐩,(CM.0.3)

where × is the vector cross product. We define the torque τ as

τ=𝐫×𝐅.(CM.0.4)

How does the angular momentum of a particle with constant mass change with time? We need to use the product rule to calculate it:

d𝐋dt=d𝐫dt×𝐩+𝐫×d𝐩dt=m𝐯×𝐯+𝐫×𝐅.

The cross product of any vector with itself is zero, so we conclude that

τ=d𝐋dt,(CM.0.5)

where we have used the definition of the torque. Therefore, we can conclude that the angular momentum of our particle is conserved when there is no applied torque, τ=0.

CM.0.1.3 Conservation of Energy

The next familiar quantity we want to define is the work W. It is a measure of how much the force is pushing in the direction of the particle's movement. The work is the sum of the dot product of the force with all the small displacements of the particle along its path, 𝐅·d𝐫. In the limit that the step size goes to zero, the sum turns into a line integral between two points 1 and 2 along the path of the particle:

W12=12𝐅·d𝐫.(CM.0.6)

In general, the work will depend on the path taken between the two points, not just on the endpoints themselves. When the mass is constant, we can use Newton's Law, Equation CM.0.2, to evaluate the integral defining the work:

𝐅·d𝐫=md𝐯dt·d𝐫dtdt=m2d(𝐯·𝐯)dtdt,

where we have also multiplied and divided by dt inside the integral sign and identified the derivative d(𝐯·𝐯)/dt=(1/2)(d𝐯/dt)·(d𝐫/dt). The quantity 𝐯·𝐯 is just the usual magnitude squared of the velocity vector, v2. The integral of the derivative of a function is just the function itself, so we have now shown that the work done on a particle by the force 𝐅 moving it between points 1 and 2 is just

W12=12mv2212mv12=T2T1,(CM.0.7)

where Ti=mvi2/2 is the particle's kinetic energy evaluated at point i.

We now want to identify the special condition under which we have a conserved quantity, which we will call the energy E, involving T. That special condition is when the work is path-independent. That is, the line integral in Equation CM.0.6 is only a function of the endpoints 1 and 2. We now want a mathematical statement that is equivalent to the work only being a function of the endpoints. Take any two paths A and B between the points 1 and 2. If the particle takes path A, then the work is W12. However, if it takes path B, then the work must be the same, W12, because the work only depends on the endpoints. However, consider now the reversed path B¯ taken by going from point 2 to point 1. The work along B¯ must be W21=W12, because integrals are odd functions under exchanging their bounds. The combined path AB¯ is a closed path from point 1 to point 2 and back to point 1; the work done on it must be W11=W12+W21=W12W12=0. However, the work only depends on the endpoints, so we could have taken any closed path that went from point 1 to point 1 and gotten the same answer, W11=0 (If we wanted to be slick, we could have said that the path where we start at point 1 and do not move immediately gives us W11=0).

We have just shown that if the work only depends on the endpoints and not on the path taken, then the work on a closed path must always be zero:

𝐅·d𝐫=0.(CM.0.8)

The converse is also true. If Equation CM.0.8 holds, then the work must be independent of path. The proof is the reverse of the argument we just gave. If W11=0 for any path, then we can choose any other point 2, and any paths A from 1 to 2 and B¯ from 2 to 1, and find W11=W12,A+W21,B¯=0. Therefore, W12,A=W21,B¯=W12,B by reversing the path B¯ to B using the minus sign to flip the integration bounds. This argument holds for any points 1 and 2 and any paths A and B between these points, so the work only depends on the endpoints 1 and 2, not the paths.

When Equation CM.0.8 holds, the force 𝐅 is called a conservative force. If we pick a convenient point 𝐎, we can define the potential energy function V(𝐩) as the negative of the work from 𝐎 to 𝐩:

V(𝐩)=𝐎𝐩𝐅·d𝐫.(CM.0.9)

Finally, we can write the work W12 in terms of the potential energy:

W12=12𝐅·d𝐫=1𝐎𝐅·d𝐫+𝐎2𝐅·d𝐫=V1V2.(CM.0.10)

Combining with Equation CM.0.7, we find

V1V2=W12=T2T1V1+T1=V2+T2.(CM.0.11)

We can define the energy E=T+V, which is conserved along the particle's path when the force is conservative.

We could have immediately found the above results using multivariable calculus, as you probably know. Stokes' Theorem says that conservative forces satisfy

×𝐅=0,(CM.0.12)

their curl is zero, so that they can be written as the gradient of a scalar function V, the exact potential energy that we defined above:

𝐅=V.(CM.0.13)

CM.0.2 Many Particles

The real world has more than one particle. Newton's second law holds for each particle individually. It is convenient to separate the total force 𝐅i on each particle into those that result from external sources, leading to a total force 𝐅i(e) on particle i, and those that result from interactions between pairs of particles in our system, leading to total forces 𝐅ji on particle i from particle j. Then, each particle in our system, with position 𝐫i and corresponding velocity, mass, and momentum, satisfies

d𝐩idt=𝐅i=𝐅i(e)+j𝐅ji.(CM.0.14)

For simplicity, we will assume that all the masses remain constant. We will also assume Newton's third law holds for all forces, so that 𝐅ji=𝐅ij, and that the forces act along the direction pointing from one object to the other, even though the third law is not true in general (e.g. for magnetic forces). We can also note that Fii=0, since particles do not exert forces on themselves.

CM.0.2.1 Conservation of Momentum

We now define the center of mass 𝐑 of our collection of particles:

𝐑=imi𝐫iM,(CM.0.15)

where M=mi is the total mass of our system of particles. The center of mass is the (mass-weighted) average position of particles in our collection. Let's try to treat our center of mass as the position of a particle with mass M. What would its equation of motion be? Since all of the masses of our particles are constant, taking two time derivatives of 𝐑 only hits the individual particles' position vectors in the sum. Therefore, we find that Newton's law for the center of mass is

Md2𝐑dt2=imid2𝐫idt2=i𝐅i(e)+i,j𝐅ji.

There are two force sums. The first is the sum over all of the external forces to give the total external force, 𝐅(e)=𝐅i(e). The second sum is the sum over all of the pairwise internal forces. By Newton's third law, this sum should be zero, since each force has an equal and opposite force canceling it. Mathematically, 𝐅ji+𝐅ij=0. So, the right-hand side just becomes the total external force.

We consider now the total momentum 𝐏 of our system, defined as the sum of the individual momenta of each particle:

𝐏=imid𝐫idt=Mddt(imi𝐫iM)=Md𝐑dt,(CM.0.16)

where we multiplied and divided by the total mass M and moved masses in and out of derivatives as necessary. Combining the previous two results, we have shown that

d𝐏dt=𝐅(e),(CM.0.17)

the total momentum of the system and the total external force on all the particles give the usual form of Newton's second law. Therefore, the total momentum of our collection of particles is conserved if the total external force is zero.

CM.0.2.2 Conservation of Angular Momentum

We have defined the total momentum and found a conservation law for it. Now we will define the total angular momentum 𝐋 and find its conservation law. The total angular momentum is

𝐋=i𝐋i=i𝐫i×𝐩i.(CM.0.18)

We should now take the time derivative of the total angular momentum. Using the product rule gives

d𝐋dt=id𝐫idt×𝐩i+i𝐫i×d𝐩idt=i𝐫i×d𝐩idt,

since the individual momenta and velocities point in the same direction, meaning their cross-product is zero. We can expand the derivatives of the momenta using Newton's law, Equation CM.0.14, so that

i𝐫i×d𝐩idt=i𝐫i×𝐅i(e)+i,j𝐫i×𝐅ji.

The first term on the right-hand side looks like the total external torque. We need to do some work on the second term. As before, the trick is to pair term ij with term ji, since 𝐅ji=𝐅ij. Then, term ij is 𝐫i×𝐅ji, while term ji can be rewritten as 𝐫j×𝐅ji. Adding these two terms and factoring out the common force gives the cross product (𝐫i𝐫j)×𝐅ji. While this form may not look any better, it is actually the key. Newton's third law not only guaranteed that the forces were equal and opposite, but it also stipulated that they acted on the line connecting the two particles. That is, 𝐅ji is parallel to the vector (𝐫i𝐫j), since the difference in position vectors is the vector pointing between the two particles. Therefore, by pairing up terms and noticing that they are zero, we have found

i,j𝐫i×𝐅ji=0.

Putting everything together and defining the total external torque τ(e)=τi(e)=𝐫i×𝐅i(e), we find

d𝐋dt=τ(e).(CM.0.19)

Therefore, if the total external torque is zero, then the total angular momentum of our collection of particles is conserved.

Optional: Total angular momentum in terms of the center of mass. We have found the conservation law for angular momentum, but it would be nice to also know an expression for the total angular momentum in terms of the center of mass 𝐑. For the total linear momentum, the expression was as simple as Equation CM.0.16. To find the corresponding equation for the total angular momentum, we will have to directly calculate from Equation CM.0.18. First, we rewrite the position vector of each particle with respect to the center of mass:

𝐫i=ρi+𝐑,(CM.0.20)

where we have defined the vectors ρi=𝐫i𝐑, the position of particle i measured in a coordinate system where the center of mass is at the origin. By taking the time derivative, we get a similar equation for the velocities, defining the new quantity 𝐰i=dρi/dt and the center of mass velocity 𝐕=d𝐑/dt:

𝐯i=𝐰i+𝐕.(CM.0.21)

Using Equations CM.0.20 and CM.0.21 in Equation CM.0.18 and distributing the cross product gives

𝐋=imi𝐫i×𝐯i=i(miρi×𝐰i+miρi×𝐕+mi𝐑×𝐰i+mi𝐑×𝐕);.

The summations can now commute with some of the cross products. In the final term, the sum just collects all of the masses into the total mass M. In the second term, we have (miρi)×𝐕. In the third term, we have 𝐑×d(miρi)/dt. Both of these latter terms rely on the quantity miρi, which we can calculate directly from Equation CM.0.20:

imiρi=imi𝐫iimi𝐑=M𝐑M𝐑=0,

where we have used the definition of the center of mass position. We could have anticipated this result, since miρi is (proportional to) the center of mass's position, but calculated in the coordinate system where the center of mass is at the origin.

Finally, we can state the result:

𝐋=M𝐑×𝐕+imiρi×𝐰i.(CM.0.22)

In words, the total angular momentum is the angular momentum of the center of mass around the origin plus the total angular momentum of the particles moving around the center of mass. The value of the first term, the angular momentum of the center of mass, depends on the origin of our coordinate system. The second term, the total angular momentum of the particles about the center of mass, does not depend on the choice of coordinate system, because it is always calculated in the coordinate frame where the center of mass is at the origin. Therefore, in all coordinate systems where the first term is zero---for example, when the origin of coordinates is moving with the center of mass's velocity 𝐕---the total angular momentum 𝐋 has the same value, equal to the total angular momentum measured in the frame where the center of mass is at the origin.

CM.0.2.3 Conservation of Energy

Finally, we want to determine when energy is conserved in our many particle system. We start again with the work done on all of the particles, which is just the sum of the work done on the individual particles:

W12=i12𝐅i·d𝐫i=i12𝐅i(e)·d𝐫i+i,j12𝐅ji·d𝐫i.(CM.0.23)

We can use the same trick of multiplying and dividing by dt, then substituting out the forces 𝐅i using Newton's law, as we did in Equation CM.0.6, to find that the right-hand side of Equation CM.0.23 turns into the change in total kinetic energy, where the total kinetic energy is a sum over the individual kinetic energies of mivi2/2. So, we can write the work as

W12=T2T1,(CM.0.24)

where

T=imivi22.(CM.0.25)

To have conservation of energy, we need conservative forces. First, we require that the external forces on each particle are conservative. That is, each Fi(e) is the negative gradient of a potential function Vi, so that

12𝐅i(e)·d𝐫i=12iVi·d𝐫i=Vi(1)Vi(2),(CM.0.26)

where the index on the gradient means we are only taking derivatives with respect to the coordinates 𝐫i.

The pairwise forces are slightly more difficult. The external forces came from a potential that only depended on the position of particle i. We expect the pairwise forces to come from a potential that can depend on two particles i and j. We must also require that the forces satisfy Newton's third law of being equal and opposite, lying on the line connecting the two particles:

𝐅ji=𝐅ij=(𝐫i𝐫j)f,(CM.0.27)

where f is a function, so that the directon of the force vector is set by the difference in position vectors. We need a potential function whose negative gradient with respect to either set of coordinates i or j produces this type of force. If the potential is a function of only the magnitude of the separation vector 𝐫ij=𝐫i𝐫j between the particles, then it gives the right form for the forces. Start by taking the gradient using the chain rule:

iVij(|𝐫i𝐫j|)=Vij(|𝐫i𝐫j|)i|𝐫i𝐫j|.

We can calculate the gradient of the length of the separation vector using its definition. In three dimensions, we know

|𝐫i𝐫j|=(xixj)2+(yiyj)2+(zizj)2,(CM.0.28)

so that, for example, the derivative with respect to xi gives (xixj)/(). Therefore, we have shown that

i|𝐫i𝐫j|=𝐫i𝐫j|𝐫i𝐫j|,(CM.0.29)

so that we have proven what we wanted to show, namely that the force resulting from the pairwise potential has the form of Equation CM.0.27

iVij(|𝐫i𝐫j|)=Vij(|𝐫i𝐫j|)𝐫i𝐫j|𝐫i𝐫j|=(𝐫i𝐫j)f.(CM.0.30)

Incidentally, Equations CM.0.29 and CM.0.30 still hold if we replace i with j or ij (the gradient with respect to the coordinates 𝐫ij), which can be proved using the same argument. Note that this result implies that the potentials Vij and Vji are not completely independent, they must be related so that the forces resulting from their gradients are the same.

Now that we have the correct form of the potential function for the pairwise forces, we can try to integrate the pairwise forces' contribution to the work, from Equation CM.0.23. We cannot directly integrate each term separately, because

𝐅ji·d𝐫i=iVij(|𝐫i𝐫j|)·d𝐫i

is not the integral of the derivative of a function of the single position 𝐫i. Instead, we should group pairs of terms. The 𝐅ij term will give

𝐅ij=jVij(|𝐫i𝐫j|)=+iVij(|𝐫i𝐫j|)=𝐅ji,

just the negative of the force 𝐅ji. However, the force 𝐅ji is dotted with the infinitesimal line element d𝐫i, while the force 𝐅ij is dotted with the line element d𝐫j. So, the combined i and j terms of the sum in the work integral is

12𝐅ji·d𝐫i+𝐅ij·d𝐫j=12iVij(rij)·(d𝐫id𝐫j).(CM.0.31)

The differential d is linear, so (d𝐫id𝐫j)=d𝐫ij, and we already proved that iVij=ijVij. We now have what we wanted, the integral of the derivative of one vector variable, 𝐫ij. Evaluating the integral for each pair and summing gives

i,j12𝐅ji·d𝐫i=i<jVij(1)Vij(2),(CM.0.32)

where the sum bounds on the right-hand side make sure that we only count the pairs once. If we wanted to extend the sum to all values of i and j, we can use the symmetry of Vij to see that ij and ji terms would be equal. Therefore, extending the sum overcounts by a factor of 2, which we can cancel by a factor of 1/2:

i<jVij(1)Vij(2)=12i,jVij(1)Vij(2).

Finally, we can see that defining the total potential energy as

V=iVi+12i,jVij(CM.0.33)

allows us to conclude that the total energy E=T+V is conserved, just as before.

Optional: Total kinetic energy in terms of the center of mass. As with the angular momentum, we can rewrite the total kinetic energy in terms of the center of mass's motion. To do so, we use the previously defined relationship for velocity relative to the center of mass, 𝐯i=𝐰i+𝐕, and just calculate vi2:

vi2=(𝐰i+𝐕)·(𝐰i+𝐕)=wi2+V2+2𝐰i·𝐕.

To find the total kinetic energy, we multiply by mi/2 and sum over i:

T=imivi22=12imiwi2+12imiV2+imi𝐰i·𝐕.

By commuting the mass-weighted sum through the time derivative and the dot product, we can see that the final term is zero by the same arguments we made earlier about the position of the center of mass in the ρ coordinates being zero, so we find that the total kinetic energy is

T=12MV2+12imiwi2.(CM.0.34)

The total kinetic energy is the kinetic energy of the center of mass plus the sum of the kinetic energy of each particle about the center of mass. Again, if V=0, then the first term disappears. So, in each coordinate frame where the center of mass has no velocity, the kinetic energy is the same, just the kinetic energy measured in the frame where the center of mass is at the origin.

CM.0.3 References

  1. Newton's laws of motion. https://www1.grc.nasa.gov/beginners-guide-to-aeronautics/newtons-laws-of-motion/.
  2. H. Goldstein. Classical Mechanics. Addison-Wesley series in advanced physics. Addison-Wesley Press, 1950.

Table of Contents

This is a teXsite.