Center of mass - cant find error in derivation

AI Thread Summary
The discussion revolves around the derivation of the center of mass for a binary system and the confusion arising from the results obtained. The original poster correctly derives the center of mass position but encounters contradictions when analyzing the distances from the masses to the center of mass, leading to the conclusion that these distances must be equal, which is incorrect. The introduction of a third mass at the center of mass raises questions about its motion and forces acting on it, with participants clarifying that while the center of mass remains constant, the third mass will experience unequal gravitational forces from the other two masses, causing it to move. The conversation emphasizes the importance of correctly applying the principles of gravity and the dynamics of multiple bodies in a system. Ultimately, the discussion highlights the complexities of gravitational interactions and the need for careful mathematical treatment in deriving physical laws.
Kyouran
Messages
70
Reaction score
10
First of all, I am sorry if this is in the wrong subforum; I'm not sure if astrodynamics should be under astronomy. The post is also fairly long. My problem is that in my derivation of the center of mass of a binary system, I get a conflicting result for the distances of the center of mass to the masses. I must have made an error in my derivation, and I guessed the best way to ask for help is to just post the entire derivation, so here goes:

Suppose I have a system of 2 masses (for example, planets). The position vector of the center of mass of such a system with respect to some arbitrary inertial reference frame with origin O is given by:

\vec{r_{cm}} = \frac{m_1 \vec{r_1} + m_2 \vec{r_2}}{m_1 + m_2} (eq. 1)

Now, I can easily prove that this point is located on the line segment connecting the 2 masses:

\vec{r_{cm}} = \frac{m_1 \vec{r_1}}{m_1 + m_2} + \frac{m_2 \vec{r_2}}{m_1 + m_2}
\vec{r_{cm}} = \frac{m_1}{m_1 + m_2}\vec{r_1} + (1 - \frac{m_1}{m_1 + m_2})\vec{r_2}
\vec{r_{cm}} = \alpha \vec{r_1} + (1 - \alpha)\vec{r_2}

Similarly, setting \beta = 1-\alpha
I get:
\vec{r_{cm}} = (1-\beta) \vec{r_1} + \beta \vec{r_2}

From these relationships, I can find the distance between the 2 masses and the center of mass as:

\vec{r_{cm,1}} = \beta \vec{r_21} \rightarrow \vec{r_{cm,1}} = \frac{m_2}{m_1 + m_2} \vec{r_{21}}
\vec{r_{cm,2}} = \alpha \vec{r_12} \rightarrow \vec{r_{cm,2}} = \frac{m_1}{m_1 + m_2} \vec{r_{12}}

Where r_ij is the vector from i to j. Now taking the norm of the vectors and taking into consideration that ||\vec{r_{12}} || = ||\vec{r_{21}} ||, i get:

|| \vec{r_{cm,1}} || = \frac{m_2}{m_1 + m_2} || \vec{r_{12}} ||
|| \vec{r_{cm,2}} || = \frac{m_1}{m_1 + m_2} || \vec{r_{12}} ||

Which gives me the relation:

\frac{|| \vec{r_{cm,1}} ||}{|| \vec{r_{cm,2}} ||} = \frac{m_2}{m_1} (eq. 2)

This seems to be logical; if the second mass is heavier then the center of mass will be closer to the second mass. If I now go back to the first equation and differentiate it twice with respect to time I get:

\frac{d^2\vec{r_{cm}}}{dt^2} = \frac{m_1 \frac{d^2\vec{r_1}}{dt^2} + m_2 \frac{d^2\vec{r_2}}{dt^2}}{m_1 + m_2}

The only forces on the masses 1 and 2 is the gravitational force, therefore:

m_1 \frac{d^2\vec{r_1}}{dt^2} = G \frac{m_1 m_2}{|| \vec{r_{12}} ||^3} \vec{r_{12}}
m_2 \frac{d^2\vec{r_2}}{dt^2} = - G \frac{m_1 m_2}{|| \vec{r_{12}} ||^3} \vec{r_{12}}

Therefore, I find that:

\frac{d^2\vec{r_{cm}}}{dt^2} = \vec{0} (eq. 3)

Which means that the center of mass does not accelerate. If I now place a mass m3 at the center of mass of this binary system, the center of mass of the entire system doesn't change, and it can be shown that equations (1), (2) and (3) still hold. Also, it can be shown that this mass will stay at the center of mass, and that it experiences no net acceleration as well. When I thus calculate the forces on this third mass I get:

m_3 \frac{d^2 \vec{r_{cm}}}{dt^2} = G m_3 ( \frac{m_1}{|| \vec{r_{cm,1}}||^3} \vec{r_{cm,1}} + \frac{m_2}{|| \vec{r_{cm,2}} ||^3} \vec{r_{cm,2}} ) = \vec{0}

from which I find:

\frac{|| \vec{r_{cm,1}}||^2}{|| \vec{r_{cm,2}}||^2} = \frac{m_1}{m_2} (eq. 4)

From equation 2 and 4 I then find:

||\vec{r_{cm,1}} ||^ 3 = ||\vec{r_{cm,2}}||^3

Now this would imply that the distance of the 2 objects to the center of mass is the same. I know this can't be true, so there must somewhere be an error in my derivation or thinking. I really don't see where this could be, so can anyone point out the mistake I made?
 
Last edited:
Astronomy news on Phys.org
Kyouran said:
If I now place a mass m3 at the center of mass of this binary system, the center of mass of the entire system doesn't change,

why not?

well,, if i now take the center of mass frame in which origin is at the center of mass, then position of center of mass is:

(m1r1 + m2r2 + m3*0)/(m1+m2+m3) which is a different position from that of earlier center of mass

EDIT: the conclusion is wrong since the above equation actually gives answer =0 which implies center of mass remains there itself even after the introduction of 3rd mass.
 
Last edited:
Kyouran said:
Also, it can be shown that this mass will stay at the center of mass, and that it experiences no net acceleration as well
How can m3 experience no force? if m3 is between m1 ans m2, if experiences a mutually opposite gravitational force which are not equal,

EDIT: but i think i hold this argument
 
Last edited:
ash64449 said:
why not?

Intuitively, it shouldn't change. But, mathematically. If you substitute:

\vec{r_3} = \frac{m_1 \vec{r_1} + m_2 \vec{r_2}}{m_1 + m_2}

into the equation for the center of mass of a three body system:

\vec{r_{cm}} = \frac{m_1 \vec{r_1} + m_2 \vec{r_2} + m_3 \vec{r_3}}{m_1 + m_2 + m_3}

and you do some algebra magic (I just checked with maple), you will find:

\vec{r_{cm}} = \frac{m_1 \vec{r_1} + m_2 \vec{r_2} + m_3 \vec{r_3}}{m_1 + m_2 + m_3} = \frac{m_1 \vec{r_1} + m_2 \vec{r_2}}{m_1 + m_2}
 
  • Like
Likes ash64449
sorry.. i understand the mistake before... but just as i tried to post u posted along!
 
ash64449 said:
How can m3 experience no force? if m3 is between m1 ans m2, if experiences a mutually opposite gravitational force which are not equal,

EDIT: but i think i hold this argument

Well, if I differentiate r3 with respect to time twice, it comes out zero according to my derivation, so somehow the net force on it must be zero.
 
Even if you introduce m3 in the center of mass position, center of mass of the new system remains constant.

But you can't say that m3 will not accelerate just because it is kept there. m3 here experience only mutually opposite gravitational force which are not equal.. hence it will accelerate.

however everything will move in such a manner that center of mass remains there again in the new system since sum total of all forces experienced by all masses is zero.
 
Kyouran said:
Well, if I differentiate r3 with respect to time twice, it comes out zero according to my derivation, so somehow the net force on it must be zero.

did you differentiate r3 w.r.t time twice in your force on m3=0 derivation or Rcm w.r.t time twice? when i read your derivation ,i see latter.
 
Never... force on m3 cannot be zero... i think that is why an error has come in your derivation. it makes sense why that error has come.
 
  • #10
ash64449 said:
did you differentiate r3 w.r.t time twice in your force on m3=0 derivation or Rcm w.r.t time twice? when i read your derivation ,i see latter.

Maybe the issue is that I assumed that r3 stays in the center of mass, and maybe it is only valid at t = 0.

EDIT: I find it weird though that a mass placed at the center of mass would start to move away from it, but it could be if m1 and m2 move accordingly.
 
  • #11
Kyouran said:
Maybe the issue is that I assumed that r3 stays in the center of mass, and maybe it is only valid at t = 0.

you are right.. why it moves afterwards? it has a net force acting on it.. this is what i said.
 
  • #12
Kyouran said:
I find it weird though that a mass placed at the center of mass would start to move away from it, but it could be if m1 and m2 move accordingly.

why weird? yes,m1 and m2 will move accordingly in such a manner that center of mass remains where it is?

motion of m1 and m2 different in first case and second case. why?

because in the second case m3 also produces a force on them which produces the different motion. but center of mass remain where it is in both the cases.
 
  • #13
Well, because i thought everything is supposed to orbit the center of mass of the system, from which it would seem that an object which is placed at the center of mass would orbit it with a radius of 0, meaning it would stay there. After all, that's the idea of gravity: everything falls towards the center of mass.
 
Last edited:
  • #14
Kyouran said:
After all, that's the idea of gravity: everything falls towards the center of mass.
That is *not* the idea of gravity. The idea of gravity is that everything falls towards everything else, and except for the Newtonian two body point mass problem, that in general is not toward the center of mass.
 
Back
Top