9.6 Center of Mass

We have been avoiding an important issue up to now: When we say that an object moves (more correctly, accelerates) in a way that obeys Newton’s second law, we have been ignoring the fact that all objects are actually made of many constituent particles. A car has an engine, steering wheel, seats, passengers; a football is leather and rubber surrounding air; a brick is made of atoms. There are many different types of particles, and they are generally not distributed uniformly in the object. How do we include these facts into our calculations?

Then too, an extended object might change shape as it moves, such as a water balloon or a cat falling (Figure 9.27). This implies that the constituent particles are applying internal forces on each other, in addition to the external force that is acting on the object as a whole. We want to be able to handle this, as well.

A multiple exposure photograph of a cat falling. In the first image, the cat is held by its feet, upside down. It is released from this position and falls, but rotates as it turns so that in the last few images, it is right side up. — **Figure 9.27** As the cat falls, its body performs complicated motions so it can land on its feet, but one point in the system moves with the simple uniform acceleration of gravity.

The problem before us, then, is to determine what part of an extended object is obeying Newton’s second law when an external force is applied and to determine how the motion of the object as a whole is affected by both the internal and external forces.

Be warned: To treat this new situation correctly, we must be rigorous and completely general. We won’t make any assumptions about the nature of the object, or of its constituent particles, or either the internal or external forces. Thus, the arguments will be complex.

Internal and External Forces

Suppose we have an extended object of mass M, made of N interacting particles. Let’s label their masses as $m_{j}$ , where $j equals 1, 2, 3, …, N$ . Note that

$M = \sum_{j = 1}^{N} m_{j} .$

(9.19)

If we apply some net external force ${\vec{F}}_{ext}$ on the object, every particle experiences some “share” or some fraction of that external force. Let:

${\vec{f}}_{j}^{ext} = the fraction of the external force that the j th particle experiences.$

Notice that these fractions of the total force are not necessarily equal; indeed, they virtually never are. (They can be, but they usually aren’t.) In general, therefore,

$f → 1 ext not equal to f → 2 ext not equal to ⋯ not equal to f → N ext.$

Next, we assume that each of the particles making up our object can interact with (apply forces on) every other particle of the object. We won’t try to guess what kind of forces they are; but since these forces are the result of particles of the object acting on other particles of the same object, we refer to them as internal forces ${\vec{f}}_{j}^{int}$ ; thus:

${\vec{f}}_{j}^{int} =$ the net internal force that the jth particle experiences from all the other particles that make up the object.

Now, the net force, internal plus external, on the jth particle is the vector sum of these:

${\vec{f}}_{j} = {\vec{f}}_{j}^{int} + {\vec{f}}_{j}^{ext} .$

where again, this is for all N particles; $j equals 1, 2, 3, and so on, N$ .

As a result of this fractional force, the momentum of each particle gets changed:

$\begin{matrix} {\vec{f}}_{j} & = & \frac{d {\vec{p}}_{j}}{d t} \\ {\vec{f}}_{j}^{int} + {\vec{f}}_{j}^{ext} & = & \frac{d {\vec{p}}_{j}}{d t} . \end{matrix}$

The net force $\vec{F}$ on the object is the vector sum of these forces:

$\begin{matrix} {\vec{F}}_{net} & = \sum_{j = 1}^{N} ({\vec{f}}_{j}^{i n t} + {\vec{f}}_{j}^{e x t}) \\ = \sum_{j = 1}^{N} {\vec{f}}_{j}^{i n t} + \sum_{j = 1}^{N} {\vec{f}}_{j}^{e x t} . \end{matrix}$

(9.22)

This net force changes the momentum of the object as a whole, and the net change of momentum of the object must be the vector sum of all the individual changes of momentum of all of the particles:

${\vec{F}}_{net} = \sum_{j = 1}^{N} \frac{d {\vec{p}}_{j}}{d t} .$

(9.23)

Combining Equation 9.22 and Equation 9.23 gives

$\sum_{j = 1}^{N} {\vec{f}}_{j}^{int} + \sum_{j = 1}^{N} {\vec{f}}_{j}^{ext} = \sum_{j = 1}^{N} \frac{d {\vec{p}}_{j}}{d t} .$

Let’s now think about these summations. First consider the internal forces term; remember that each ${\vec{f}}_{j}^{int}$ is the force on the jth particle from the other particles in the object. But by Newton’s third law, for every one of these forces, there must be another force that has the same magnitude, but the opposite sign (points in the opposite direction). These forces do not cancel; however, that’s not what we’re doing in the summation. Rather, we’re simply mathematically adding up all the internal force vectors. That is, in general, the internal forces for any individual part of the object won’t cancel, but when all the internal forces are added up, the internal forces must cancel in pairs. It follows, therefore, that the sum of all the internal forces must be zero:

$\sum_{j = 1}^{N} {\vec{f}}_{j}^{int} = 0 .$

(This argument is subtle, but crucial; take plenty of time to completely understand it.)

For the external forces, this summation is simply the total external force that was applied to the whole object:

$\sum_{j = 1}^{N} {\vec{f}}_{j}^{ext} = {\vec{F}}_{ext} .$

As a result,

${\vec{F}}_{ext} = \sum_{j = 1}^{N} \frac{d {\vec{p}}_{j}}{d t} .$

(9.25)

This is an important result. Equation 9.25 tells us that the total change of momentum of the entire object (all N particles) is due only to the external forces; the internal forces do not change the momentum of the object as a whole. This is why you can’t lift yourself in the air by standing in a basket and pulling up on the handles: For the system of you + basket, your upward pulling force is an internal force.

Force and Momentum

Remember that our actual goal is to determine the equation of motion for the entire object (the entire system of particles). To that end, let’s define:

${\vec{p}}_{CM} =$ the total momentum of the system of N particles (the reason for the subscript will become clear shortly)

Then we have

$p → sub CM ≡ ∑ j equals 1 N p → sub j,$

and therefore Equation 9.25 can be written simply as

$\vec{F} = \frac{d {\vec{p}}_{CM}}{d t} .$

(9.26)

Since this change of momentum is caused by only the net external force, we have dropped the “ext” subscript.

This is Newton’s second law, but now for the entire extended object. If this feels a bit anticlimactic, remember what is hiding inside it: ${\vec{p}}_{CM}$ is the vector sum of the momentum of (in principle) hundreds of thousands of billions of billions of particles $(6.02 \times 10^{23})$ , all caused by one simple net external force—a force that you can calculate.

Center of Mass

Our next task is to determine what part of the extended object, if any, is obeying Equation 9.26.

It’s tempting to take the next step; does the following equation mean anything?

$\vec{F} = M \vec{a}$

(9.27)

If it does mean something (acceleration of what, exactly?), then we could write

$M \vec{a} = \frac{d {\vec{p}}_{CM}}{d t}$

and thus

$M \vec{a} = \sum_{j = 1}^{N} \frac{d {\vec{p}}_{j}}{d t} = \frac{d}{d t} \sum_{j = 1}^{N} {\vec{p}}_{j} .$

which follows because the derivative of a sum is equal to the sum of the derivatives.

Now, ${\vec{p}}_{j}$ is the momentum of the jth particle. Defining the positions of the constituent particles (relative to some coordinate system) as ${\vec{r}}_{j} = (x_{j}, y_{j}, z_{j})$ , we thus have

${\vec{p}}_{j} = m_{j} {\vec{v}}_{j} = m_{j} \frac{d {\vec{r}}_{j}}{d t} .$

Substituting back, we obtain

$\begin{matrix} M \vec{a} & = \frac{d}{d t} \sum_{j = 1}^{N} m_{j} \frac{d {\vec{r}}_{j}}{d t} \\ = \frac{d^{2}}{d t^{2}} \sum_{j = 1}^{N} m_{j} {\vec{r}}_{j} . \end{matrix}$

Dividing both sides by M (the total mass of the extended object) gives us

$\vec{a} = \frac{d^{2}}{d t^{2}} (\frac{1}{M} \sum_{j = 1}^{N} m_{j} {\vec{r}}_{j}) .$

(9.28)

Thus, the point in the object that traces out the trajectory dictated by the applied force in Equation 9.27 is inside the parentheses in Equation 9.28.

Looking at this calculation, notice that (inside the parentheses) we are calculating the product of each particle’s mass with its position, adding all N of these up, and dividing this sum by the total mass of particles we summed. This is reminiscent of an average; inspired by this, we’ll (loosely) interpret it to be the weighted average position of the mass of the extended object. It’s actually called the center of mass of the object. Notice that the position of the center of mass has units of meters; that suggests a definition:

${\vec{r}}_{CM} \equiv \frac{1}{M} \sum_{j = 1}^{N} m_{j} {\vec{r}}_{j} .$

(9.29)

So, the point that obeys Equation 9.26 (and therefore Equation 9.27 as well) is the center of mass of the object, which is located at the position vector ${\vec{r}}_{CM}$ .

It may surprise you to learn that there does not have to be any actual mass at the center of mass of an object. For example, a hollow steel sphere with a vacuum inside it is spherically symmetrical (meaning its mass is uniformly distributed about the center of the sphere); all of the sphere’s mass is out on its surface, with no mass inside. But it can be shown that the center of mass of the sphere is at its geometric center, which seems reasonable. Thus, there is no mass at the position of the center of mass of the sphere. (Another example is a doughnut.) The procedure to find the center of mass is illustrated in Figure 9.28.

An illustration of finding the center of mass of three particles. Figure a shows the locations of the three particles in he x y plane. m 1 is in the second quadrant. Vector r 1 starts at the origin and extends to the location of m 1. m 2 is in the first quadrant. Vector r 2 starts at the origin and extends to the location of m 2. m 1 is in the fourth quadrant. Vector r 3 starts at the origin and extends to the location of m 3. Vector r 1 is the shortest of the vectors in the diagram, and r 2 is the longest. Figure b shows the vectors m 1 r 1, m 2 r 2 and m 3 r 3. Vector m 1 r 1 points in the same direction as vector r 1 in figure a, but is longer than r 1. Vector m 2 r 2 points in the same direction as vector r 1 in figure a, but is shorter than r 2. . Vector m 3 r 3 points in the same direction as vector r 3 in figure a, but is shorter than r 3. Vector m 1 r 1 is the longest vector in the diagram. Vectors m 2 r 2 and m 3 r 3 appear to be of equal length. Figure c shows the vector sum of m 1 r 1, m2 r 2 and m 3 r 3, which have been drawn in blue and placed head to tail. The red vector m 1 r 1 plus m 2 r 2 plus m 3 r 3 is the vector from the tail of m 1 r 1 to the head of m 3 r 3. Figure d shows the red vector m 1 r 1 plus m 2 r 2 plus m 3 r 3 all divded by the sum m 1 plus m 2 plus m 3. This vector is in the same direction as the vecor m 1 r 1 plus m 2 r 2 plus m 3 r 3 in figure c, but shorter. — **Figure 9.28** Finding the center of mass of a system of three different particles. (a) Position vectors are created for each object. (b) The position vectors are multiplied by the mass of the corresponding object. (c) The scaled vectors from part (b) are added together. (d) The final vector is divided by the total mass. This vector points to the center of mass of the system. Note that no mass is actually present at the center of mass of this system.

Since ${\vec{r}}_{j} = x_{j} \hat{i} + y_{j} \hat{j} + z_{j} \hat{k}$ , it follows that:

$r_{CM, x} = \frac{1}{M} \sum_{j = 1}^{N} m_{j} x_{j}$

(9.30)

$r_{CM, y} = \frac{1}{M} \sum_{j = 1}^{N} m_{j} y_{j}$

(9.31)

$r_{CM, z} = \frac{1}{M} \sum_{j = 1}^{N} m_{j} z_{j}$

(9.32)

and thus

$\begin{matrix} {\vec{r}}_{CM} = r_{CM, x} \hat{i} + r_{CM, y} \hat{j} + r_{CM, z} \hat{k} \\ r_{CM} = | {\vec{r}}_{CM} | = {(r_{CM, x}^{2} + r_{CM, y}^{2} + r_{CM, z}^{2})}^{1 / 2} . \end{matrix}$

Therefore, you can calculate the components of the center of mass vector individually.

Finally, to complete the kinematics, the instantaneous velocity of the center of mass is calculated exactly as you might suspect:

${\vec{v}}_{CM} = \frac{d}{d t} (\frac{1}{M} \sum_{j = 1}^{N} m_{j} {\vec{r}}_{j}) = \frac{1}{M} \sum_{j = 1}^{N} m_{j} {\vec{v}}_{j}$

and this, like the position, has x-, y-, and z-components.

To calculate the center of mass in actual situations, we recommend the following procedure:

Here are two examples that will give you a feel for what the center of mass is.

Example 1 Center of Mass of the Earth-Moon System

Using data from text appendix, determine how far the center of mass of the Earth-moon system is from the center of Earth. Compare this distance to the radius of Earth, and comment on the result. Ignore the other objects in the solar system.

Strategy We get the masses and separation distance of the Earth and moon, impose a coordinate system, and use Equation 9.29 with just $N equals 2$ objects. We use a subscript “e” to refer to Earth, and subscript “m” to refer to the moon.

Solution Define the origin of the coordinate system as the center of Earth. Then, with just two objects, Equation 9.29 becomes

$R = \frac{m_{e} r_{e} + m_{m} r_{m}}{m_{e} + m_{m}} .$

From Appendix D,

$m_{e} = 5.97 \times 10^{24} kg$

$m_{m} = 7.36 \times 10^{22} kg$

$r_{m} = 3.82 \times 10^{8} m .$

We defined the center of Earth as the origin, so $r sub e equals 0 m$ . Inserting these into the equation for R gives

$\begin{matrix} R & = \frac{(5.97 \times 10^{24} kg) (0 m) + (7.36 \times 10^{22} kg) (3.82 \times 10^{8} m)}{5.97 \times 10^{24} kg + 7.36 \times 10^{22} kg} \\ = 4.64 \times 10^{6} m. \end{matrix}$

Significance The radius of Earth is $6.37 \times 10^{6} m$ , so the center of mass of the Earth-moon system is (6.37 − 4.64) $\times 10^{6} m = 1.73 \times 10^{6} m = 1730 km$ (roughly 1080 miles) below the surface of Earth. The location of the center of mass is shown (not to scale).

Example 2 Center of Mass of a Salt Crystal

Figure 9.29 shows a single crystal of sodium chloride—ordinary table salt. The sodium and chloride ions form a single unit, NaCl. When multiple NaCl units group together, they form a cubic lattice. The smallest possible cube (called the unit cell) consists of four sodium ions and four chloride ions, alternating. The length of one edge of this cube (i.e., the bond length) is $2.36 \times 10^{−10} m$ . Find the location of the center of mass of the unit cell. Specify it either by its coordinates $(r_{CM, x}, r_{CM, y}, r_{CM, z})$ , or by $r_{CM}$ and two angles.

Figure 9.29 A drawing of a sodium chloride (NaCl) crystal.

Strategy We can look up all the ion masses. If we impose a coordinate system on the unit cell, this will give us the positions of the ions. We can then apply Equation 9.30, Equation 9.31, and Equation 9.32 (along with the Pythagorean theorem).

Solution Define the origin to be at the location of the chloride ion at the bottom left of the unit cell. Figure 9.30 shows the coordinate system.

Figure 9.30 A single unit cell of a NaCl crystal.

There are eight ions in this crystal, so N = 8:

${\vec{r}}_{CM} = \frac{1}{M} \sum_{j = 1}^{8} m_{j} {\vec{r}}_{j} .$

The mass of each of the chloride ions is

$35.453 u \times \frac{1.660 \times 10^{−27} kg}{u} = 5.885 \times 10^{−26} kg$

so we have

$m_{1} = m_{3} = m_{6} = m_{8} = 5.885 \times 10^{−26} kg .$

For the sodium ions,

$m_{2} = m_{4} = m_{5} = m_{7} = 3.816 \times 10^{−26} kg .$

The total mass of the unit cell is therefore

$M = (4) (5.885 \times 10^{−26} kg) + (4) (3.816 \times 10^{−26} kg) = 3.880 \times 10^{−25} kg .$

From the geometry, the locations are

$\begin{matrix} {\vec{r}}_{1} = 0 \\ {\vec{r}}_{2} = (2.36 \times 10^{−10} m) \hat{i} \\ {\vec{r}}_{3} = r_{3 x} \hat{i} + r_{3 y} \hat{j} = (2.36 \times 10^{−10} m) \hat{i} + (2.36 \times 10^{−10} m) \hat{j} \\ {\vec{r}}_{4} = (2.36 \times 10^{−10} m) \hat{j} \\ {\vec{r}}_{5} = (2.36 \times 10^{−10} m) \vec{k} \\ {\vec{r}}_{6} = r_{6 x} \hat{i} + r_{6 z} \hat{k} = (2.36 \times 10^{−10} m) \hat{i} + (2.36 \times 10^{−10} m) \hat{k} \\ {\vec{r}}_{7} = r_{7 x} \hat{i} + r_{7 y} \hat{j} + r_{7 z} \hat{k} = (2.36 \times 10^{−10} m) \hat{i} + (2.36 \times 10^{−10} m) \hat{j} + (2.36 \times 10^{−10} m) \hat{k} \\ {\vec{r}}_{8} = r_{8 y} \hat{j} + r_{8 z} \hat{k} = (2.36 \times 10^{−10} m) \hat{j} + (2.36 \times 10^{−10} m) \hat{k} . \end{matrix}$

Substituting:

$\begin{matrix} | {\vec{r}}_{CM, x} | & = \frac{1}{M} \sum_{j = 1}^{8} m_{j} {(r_{x})}_{j} \\ = \frac{1}{M} (m_{1} r_{1 x} + m_{2} r_{2 x} + m_{3} r_{3 x} + m_{4} r_{4 x} + m_{5} r_{5 x} + m_{6} r_{6 x} + m_{7} r_{7 x} + m_{8} r_{8 x}) \\ = \frac{1}{3.8804 \times 10^{−25} kg} [(5.885 \times 10^{−26} kg) (0 m) + (3.816 \times 10^{−26} kg) (2.36 \times 10^{−10} m) \\ + (5.885 \times 10^{−26} kg) (2.36 \times 10^{−10} m) \\ + (3.816 \times 10^{−26} kg) (2.36 \times 10^{−10} m) + 0 + 0 \\ + (3.816 \times 10^{−26} kg) (2.36 \times 10^{−10} m) + 0] \\ = 1.18 \times 10^{−10} m. \end{matrix}$

Similar calculations give $r_{CM, y} = r_{CM, z} = 1.18 \times 10^{−10} m$ (you could argue that this must be true, by symmetry, but it’s a good idea to check).

Significance Although this is a great exercise to determine the center of mass given a Chloride ion at the origin, in fact the origin could be chosen at any location. Therefore, there is no meaningful application of the center of mass of a unit cell beyond as an exercise.

Check Your Understanding

Suppose you have a macroscopic salt crystal (that is, a crystal that is large enough to be visible with your unaided eye). It is made up of a huge number of unit cells. Is the center of mass of this crystal necessarily at the geometric center of the crystal?

On a macroscopic scale, the size of a unit cell is negligible and the crystal mass may be considered to be distributed homogeneously throughout the crystal. Thus,
${\vec{r}}_{CM} = \frac{1}{M} \sum_{j = 1}^{N} m_{j} {\vec{r}}_{j} = \frac{1}{M} \sum_{j = 1}^{N} m {\vec{r}}_{j} = \frac{m}{M} \sum_{j = 1}^{N} {\vec{r}}_{j} = \frac{N m}{M} \frac{\sum_{j = 1}^{N} {\vec{r}}_{j}}{N}$
where we sum over the number N of unit cells in the crystal and m is the mass of a unit cell. Because Nm = M, we can write
${\vec{r}}_{CM} = \frac{m}{M} \sum_{j = 1}^{N} {\vec{r}}_{j} = \frac{N m}{M} \frac{\sum_{j = 1}^{N} {\vec{r}}_{j}}{N} = \frac{1}{N} \sum_{j = 1}^{N} {\vec{r}}_{j} .$
This is the definition of the geometric center of the crystal, so the center of mass is at the same point as the geometric center.

Two crucial concepts come out of these examples:

As with all problems, you must define your coordinate system and origin. For center-of-mass calculations, it often makes sense to choose your origin to be located at one of the masses of your system. That choice automatically defines its distance in Equation 9.29 to be zero. However, you must still include the mass of the object at your origin in your calculation of M, the total mass Equation 9.19. In the Earth-moon system example, this means including the mass of Earth. If you hadn’t, you’d have ended up with the center of mass of the system being at the center of the moon, which is clearly wrong.
In the second example (the salt crystal), notice that there is no mass at all at the location of the center of mass. This is an example of what we stated above, that there does not have to be any actual mass at the center of mass of an object.

Center of Mass of Continuous Objects

If the object in question has its mass distributed uniformly in space, rather than as a collection of discrete particles, then $m_{j} \to d m$ , and the summation becomes an integral:

${\vec{r}}_{CM} = \frac{1}{M} \int \vec{r} d m .$

(9.34)

In this context, r is a characteristic dimension of the object (the radius of a sphere, the length of a long rod). To generate an integrand that can actually be calculated, you need to express the differential mass element dm as a function of the mass density of the continuous object, and the dimension r. An example will clarify this.

Example 3 CM of a Uniform Thin Hoop

Find the center of mass of a uniform thin hoop (or ring) of mass M and radius r.

Strategy First, the hoop’s symmetry suggests the center of mass should be at its geometric center. If we define our coordinate system such that the origin is located at the center of the hoop, the integral should evaluate to zero.

We replace dm with an expression involving the density of the hoop and the radius of the hoop. We then have an expression we can actually integrate. Since the hoop is described as “thin,” we treat it as a one-dimensional object, neglecting the thickness of the hoop. Therefore, its density is expressed as the number of kilograms of material per meter. Such a density is called a linear mass density, and is given the symbol $λ$ ; this is the Greek letter “lambda,” which is the equivalent of the English letter “l” (for “linear”).

Since the hoop is described as uniform, this means that the linear mass density $λ$ is constant. Thus, to get our expression for the differential mass element dm, we multiply $λ$ by a differential length of the hoop, substitute, and integrate (with appropriate limits for the definite integral).

Solution First, define our coordinate system and the relevant variables (Figure 9.31).

Figure 9.31 Finding the center of mass of a uniform hoop. We express the coordinates of a differential piece of the hoop, and then integrate around the hoop.

The center of mass is calculated with Equation 9.34:

${\vec{r}}_{CM} = \frac{1}{M} \int_{a}^{b} \vec{r} d m .$

We have to determine the limits of integration a and b. Expressing $\vec{r}$ in component form gives us

${\vec{r}}_{CM} = \frac{1}{M} \int_{a}^{b} [(r cos θ) \hat{i} + (r sin θ) \hat{j}] d m .$

In the diagram, we highlighted a piece of the hoop that is of differential length ds; it therefore has a differential mass $d m = λ d s$ . Substituting:

${\vec{r}}_{CM} = \frac{1}{M} \int_{a}^{b} [(r cos θ) \hat{i} + (r sin θ) \hat{j}] λ d s .$

However, the arc length ds subtends a differential angle $d θ$ , so we have

$d s = r d θ$

and thus

${\vec{r}}_{CM} = \frac{1}{M} \int_{a}^{b} [(r cos θ) \hat{i} + (r sin θ) \hat{j}] λ r d θ .$

One more step: Since $λ$ is the linear mass density, it is computed by dividing the total mass by the length of the hoop:

$λ = \frac{M}{2 π r}$

giving us

$\begin{matrix} {\vec{r}}_{CM} & = \frac{1}{M} \int_{a}^{b} [(r cos θ) \hat{i} + (r sin θ) \hat{j}] (\frac{M}{2 π r}) r d θ \\ = \frac{1}{2 π} \int_{a}^{b} [(r cos θ) \hat{i} + (r sin θ) \hat{j}] d θ . \end{matrix}$

Notice that the variable of integration is now the angle $θ$ . This tells us that the limits of integration (around the circular hoop) are $θ = 0 to θ = 2 π$ , so $a equals 0$ and $b equals 2 π$ . Also, for convenience, we separate the integral into the x- and y-components of ${\vec{r}}_{CM}$ . The final integral expression is

$\begin{matrix} {\vec{r}}_{CM} & = r_{CM, x} \hat{i} + r_{CM, y} \hat{j} \\ = [\frac{1}{2 π} \int_{0}^{2 π} (r cos θ) d θ] \hat{i} + [\frac{1}{2 π} \int_{0}^{2 π} (r sin θ) d θ] \hat{j} \\ = 0 \hat{i} + 0 \hat{j} = \vec{0} \end{matrix}$

as expected.

Center of Mass and Conservation of Momentum

How does all this connect to conservation of momentum?

Suppose you have N objects with masses $m sub 1, m sub 2, m sub 3,... m sub N$ and initial velocities $v → sub 1, v → sub 2, v → sub 3,..., v → sub N$ . The center of mass of the objects is

${\vec{r}}_{CM} = \frac{1}{M} \sum_{j = 1}^{N} m_{j} {\vec{r}}_{j} .$

Its velocity is

${\vec{v}}_{CM} = \frac{d {\vec{r}}_{CM}}{d t} = \frac{1}{M} \sum_{j = 1}^{N} m_{j} \frac{d {\vec{r}}_{j}}{d t}$

and thus the initial momentum of the center of mass is

$\begin{matrix} {[M \frac{d {\vec{r}}_{CM}}{d t}]}_{i} & = & \sum_{j = 1}^{N} m_{j} \frac{d {\vec{r}}_{j, i}}{d t} \\ M {\vec{v}}_{CM,i} & = & \sum_{j = 1}^{N} m_{j} {\vec{v}}_{j, i} . \end{matrix}$

After these masses move and interact with each other, the momentum of the center of mass is

$M {\vec{v}}_{CM,f} = \sum_{j = 1}^{N} m_{j} {\vec{v}}_{j, f} .$

But conservation of momentum tells us that the right-hand side of both equations must be equal, which says

$M {\vec{v}}_{CM,f} = M {\vec{v}}_{CM,i} .$

(9.36)

This result implies that conservation of momentum is expressed in terms of the center of mass of the system. Notice that as an object moves through space with no net external force acting on it, an individual particle of the object may accelerate in various directions, with various magnitudes, depending on the net internal force acting on that object at any time. (Remember, it is only the vector sum of all the internal forces that vanishes, not the internal force on a single particle.) Thus, such a particle’s momentum will not be constant—but the momentum of the entire extended object will be, in accord with Equation 9.36.

Equation 9.36 implies another important result: Since M represents the mass of the entire system of particles, it is necessarily constant. (If it isn’t, we don’t have a closed system, so we can’t expect the system’s momentum to be conserved.) As a result, Equation 9.36 implies that, for a closed system,

${\vec{v}}_{CM,f} = {\vec{v}}_{CM,i} .$

(9.37)

That is to say, in the absence of an external force, the velocity of the center of mass never changes.

You might be tempted to shrug and say, “Well yes, that’s just Newton’s first law,” but remember that Newton’s first law discusses the constant velocity of a particle, whereas Equation 9.37 applies to the center of mass of a (possibly vast) collection of interacting particles, and that there may not be any particle at the center of mass at all! So, this really is a remarkable result.

Example 4 Fireworks Display

When a fireworks rocket explodes, thousands of glowing fragments fly outward in all directions, and fall to Earth in an elegant and beautiful display (Figure 9.32). Describe what happens, in terms of conservation of momentum and center of mass.

Figure 9.32 These exploding fireworks are a vivid example of conservation of momentum and the motion of the center of mass.

The picture shows radial symmetry about the central points of the explosions; this suggests the idea of center of mass. We can also see the parabolic motion of the glowing particles; this brings to mind projectile motion ideas.

Solution Initially, the fireworks rocket is launched and flies more or less straight upward; this is the cause of the more-or-less-straight, white trail going high into the sky below the explosion in the upper-right of the picture (the yellow explosion). This trail is not parabolic because the explosive shell, during its launch phase, is actually a rocket; the impulse applied to it by the ejection of the burning fuel applies a force on the shell during the rise-time interval. (This is a phenomenon we will study in the next section.) The shell has multiple forces on it; thus, it is not in free-fall prior to the explosion.

At the instant of the explosion, the thousands of glowing fragments fly outward in a radially symmetrical pattern. The symmetry of the explosion is the result of all the internal forces summing to zero $open parenthesis ∑ j f → j int equals 0 close parenthesis;$ for every internal force, there is another that is equal in magnitude and opposite in direction.

However, as we learned above, these internal forces cannot change the momentum of the center of mass of the (now exploded) shell. Since the rocket force has now vanished, the center of mass of the shell is now a projectile (the only force on it is gravity), so its trajectory does become parabolic. The two red explosions on the left show the path of their centers of mass at a slightly longer time after explosion compared to the yellow explosion on the upper right.

In fact, if you look carefully at all three explosions, you can see that the glowing trails are not truly radially symmetric; rather, they are somewhat denser on one side than the other. Specifically, the yellow explosion and the lower middle explosion are slightly denser on their right sides, and the upper-left explosion is denser on its left side. This is because of the momentum of their centers of mass; the differing trail densities are due to the momentum each piece of the shell had at the moment of its explosion. The fragment for the explosion on the upper left of the picture had a momentum that pointed upward and to the left; the middle fragment’s momentum pointed upward and slightly to the right; and the right-side explosion clearly upward and to the right (as evidenced by the white rocket exhaust trail visible below the yellow explosion).

Finally, each fragment is a projectile on its own, thus tracing out thousands of glowing parabolas.

Significance In the discussion above, we said, “…the center of mass of the shell is now a projectile (the only force on it is gravity)….” This is not quite accurate, for there may not be any mass at all at the center of mass; in which case, there could not be a force acting on it. This is actually just verbal shorthand for describing the fact that the gravitational forces on all the particles act so that the center of mass changes position exactly as if all the mass of the shell were always located at the position of the center of mass.

You may sometimes hear someone describe an explosion by saying something like, “the fragments of the exploded object always move in a way that makes sure that the center of mass continues to move on its original trajectory.” This makes it sound as if the process is somewhat magical: how can it be that, in every explosion, it always works out that the fragments move in just the right way so that the center of mass’ motion is unchanged? Phrased this way, it would be hard to believe no explosion ever does anything differently.

The explanation of this apparently astonishing coincidence is: We defined the center of mass precisely so this is exactly what we would get. Recall that first we defined the momentum of the system:

${\vec{p}}_{CM} = \sum_{j = 1}^{N} \frac{d {\vec{p}}_{j}}{d t} .$

We then concluded that the net external force on the system (if any) changed this momentum:

$\vec{F} = \frac{d {\vec{p}}_{CM}}{d t}$

and then—and here’s the point—we defined an acceleration that would obey Newton’s second law. That is, we demanded that we should be able to write

$\vec{a} = \frac{\vec{F}}{M}$

which requires that

$\vec{a} = \frac{d^{2}}{d t^{2}} (\frac{1}{M} \sum_{j = 1}^{N} m_{j} {\vec{r}}_{j}) .$

where the quantity inside the parentheses is the center of mass of our system. So, it’s not astonishing that the center of mass obeys Newton’s second law; we defined it so that it would.

Summary

An extended object (made up of many objects) has a defined position vector called the center of mass.
The center of mass can be thought of, loosely, as the average location of the total mass of the object.
The center of mass of an object traces out the trajectory dictated by Newton’s second law, due to the net external force.
The internal forces within an extended object cannot alter the momentum of the extended object as a whole.

Conceptual Questions

Suppose a fireworks shell explodes, breaking into three large pieces for which air resistance is negligible. How does the explosion affect the motion of the center of mass? How would it be affected if the pieces experienced significantly more air resistance than the intact shell?

Problems

Three point masses are placed at the corners of a triangle as shown in the figure below.

A right triangle with sides length 3 c m and 4 c m has masses of 100 g at the vertex between the hypotenuse and the 4 c m side, 75 g at the vertex between the hypotenuse and the 3 c m side, and 150 g at the vertex between the the 3 c m side and the 4 c m side.

Find the center of mass of the three-mass system.

With the origin defined to be at the position of the 150-g mass, $x sub CM equals −1.23 cm$ and $y sub CM equals 0.69 cm$

Two particles of masses $m sub 1$ and $m sub 2$ separated by a horizontal distance D are released from the same height h at the same time. Find the vertical position of the center of mass of these two particles at a time before the two particles strike the ground. Assume no air resistance.

Two particles of masses $m sub 1$ and $m sub 2$ separated by a horizontal distance D are let go from the same height h at different times. Particle 1 starts at $t equals 0$ , and particle 2 is let go at $t = T$ . Find the vertical position of the center of mass at a time before the first particle strikes the ground. Assume no air resistance.

$y_{CM} = {\begin{matrix} \frac{h}{2} - \frac{1}{4} g t^{2}, t < T \\ h - \frac{1}{2} g t^{2} - \frac{1}{4} g T^{2} + \frac{1}{2} g t T, t \geq T \end{matrix}$

Two particles of masses $m sub 1$ and $m sub 2$ move uniformly in different circles of radii $R sub 1$ and $R sub 2$ about origin in the x, y-plane. The x- and y-coordinates of the center of mass and that of particle 1 are given as follows (where length is in meters and t in seconds):
$x_{1} (t) = 4 cos (2 t), y_{1} (t) = 4 sin (2 t)$

and:
$x_{CM} (t) = 3 cos (2 t), y_{CM} (t) = 3 sin (2 t) .$

Find the radius of the circle in which particle 1 moves.
Find the x- and y-coordinates of particle 2 and the radius of the circle this particle moves.

Two particles of masses $m sub 1$ and $m sub 2$ move uniformly in different circles of radii $R sub 1$ and $R sub 2$ about the origin in the x, y-plane. The coordinates of the two particles in meters are given as follows ( $z equals 0$ for both). Here t is in seconds:
$\begin{matrix} x_{1} (t) & = & 4 cos (2 t) \\ y_{1} (t) & = & 4 sin (2 t) \\ x_{2} (t) & = & 2 cos (3 t - \frac{π}{2}) \\ y_{2} (t) & = & 2 sin (3 t - \frac{π}{2}) \end{matrix}$

Find the radii of the circles of motion of both particles.
Find the x- and y-coordinates of the center of mass.
Decide if the center of mass moves in a circle by plotting its trajectory.

a. $R sub 1 equals 4 m$ , $R sub 2 equals 2 m$ ; b. $X_{CM} = \frac{m_{1} x_{1} + m_{2} x_{2}}{m_{1} + m_{2}}, Y_{CM} = \frac{m_{1} y_{1} + m_{2} y_{2}}{m_{1} + m_{2}}$ ; c. yes, with $R = \frac{1}{m_{1} + m_{2}} \sqrt{16 m_{1}^{2} + 4 m_{2}^{2}}$

Find the center of mass of a one-meter long rod, made of 50 cm of iron (density $8 \frac{g}{{cm}^{3}}$ ) and 50 cm of aluminum (density $2.7 \frac{g}{{cm}^{3}}$ ).

Find the center of mass of a rod of length L whose mass density changes from one end to the other quadratically. That is, if the rod is laid out along the x-axis with one end at the origin and the other end at $x = L$ , the density is given by $ρ (x) = ρ_{0} + (ρ_{1} - ρ_{0}) {(\frac{x}{L})}^{2}$ , where $ρ sub 0$ and $ρ sub 1$ are constant values.

$x_{c m} = \frac{3}{4} L (\frac{ρ_{1} + ρ_{0}}{ρ_{1} + 2 ρ_{0}})$

Find the center of mass of a rectangular block of length a and width b that has a nonuniform density such that when the rectangle is placed in the x,y-plane with one corner at the origin and the block placed in the first quadrant with the two edges along the x- and y-axes, the density is given by $ρ open parenthesis x, y close parenthesis equals ρ sub 0 x$ , where $ρ sub 0$ is a constant.

Find the center of mass of a rectangular material of length a and width b made up of a material of nonuniform density. The density is such that when the rectangle is placed in the xy-plane, the density is given by $ρ open parenthesis x, y close parenthesis equals ρ sub 0 x y$ .

$(\frac{2 a}{3}, \frac{2 b}{3})$

A cube of side a is cut out of another cube of side b as shown in the figure below.

A large cube of side b has a cube of side a cut out of its bottom left front corner.

Find the location of the center of mass of the structure. (Hint: Think of the missing part as a negative mass overlapping a positive mass.)

Find the center of mass of a cone of uniform density that has a radius R at the base, height h, and mass M. Let the origin be at the center of the base of the cone and have +z going through the cone vertex.

$open parenthesis x sub CM, y sub CM, z sub CM close parenthesis equals open parenthesis 0,0, h / 4 close parenthesis$

Find the center of mass of a thin wire of mass m and length L bent in a semicircular shape. Let the origin be at the center of the semicircle and have the wire arc from the +x axis, cross the +y axis, and terminate at the −x axis.

Find the center of mass of a uniform thin semicircular plate of radius R. Let the origin be at the center of the semicircle, the plate arc from the +x axis to the −x axis, and the z axis be perpendicular to the plate.

$open parenthesis x sub CM, y sub CM, z sub CM close parenthesis equals open parenthesis 0, 4 R / open parenthesis 3 π close parenthesis, 0 close parenthesis$

Find the center of mass of a sphere of mass M and radius R and a cylinder of mass m, radius r, and height h arranged as shown below.

Figure a has a sphere on top of a vertical cylinder. Figure b has a sphere centered on top of a horizontal cylinder.

Express your answers in a coordinate system that has the origin at the center of the cylinder.