Euler-Lagrange Equations for Two Body Problem

[_Kenny_]

Homework Statement

I'm trying to do a little review of Lagrangian Mechanics through studying the two-body problem for a radial force. I have the Lagrangian of the system L=\frac{1}{2}m_1\dot{\vec{r_1}}^{2}+\frac{1}{2}m_2\dot{\vec{r_2}}^{2}-V(|{\vec{r_1}-\vec{r_2}}|). Now I'm trying to find the Euler-Lagrange Equations for r_1 and r_2 but I'm confused about taking the derivative of the potential portion with respect to either r_1 or r_2. Please call me stupid and then tell me why I'm being stupid here.

Homework Equations

L=\frac{1}{2}m_1\dot{\vec{r_1}}^{2}+\frac{1}{2}m_2\dot{\vec{r_2}}^{2}-V(|{\vec{r_1}-\vec{r_2}}|)

\frac{dL}{dq}=\frac{d}{dt}\frac{dL}{d\dot{q}}

[/B]

The Attempt at a Solution

\frac{\partial L}{\partial r_1}=-\frac{\partial V(|{\vec{r_1}-\vec{r_2}}|)}{\partial r_1}=...?[/B]

 

[ShayanJ]

You're a stupid because you're doing it in terms of \vec r_1 and \vec r_2. Change coordinates to \vec R=\frac{m_1\vec r_1+m_2\vec r_2}{m_1+m_2} and \vec r=\vec r_1-\vec r_2(centre of mass coordinates) and you'll see that there is only one variable here!

 

[jk86]

In central force problems, it is useful to instead use the coordinates:
<br /> \begin{align}<br /> \vec{r} &amp;= \vec{r}_2-\vec{r}_1\\<br /> \vec{R} &amp;= \frac{m_1\vec{r}_1 + m_2\vec{r}_2}{M}<br /> \end{align}<br />
The reason is that with these coordinates the Lagrangian reduces to something that looks like a single body with a "reduced mass" \mu:
<br /> L=\frac{1}{2}\underbrace{\frac{m_1m_2}{m_1+m_2}}_{\mu}\dot{\vec{r}}^2 + \underbrace{\frac{1}{2}M\dot{\vec{R}}^2}_{0} - U(\vec{r}) \ ,<br />

 

[_Kenny_]

I'm not switching to the easier coordinates just yet, although I do know they simplify things greatly. What would the Euler-Lagrange Equations be if I remained with the non-reduced coordinates for right now?

 

[ShayanJ]

I'm not switching to the easier coordinates just yet, although I do know they simplify things greatly. What would the Euler-Lagrange Equations be if I remained with the non-reduced coordinates for right now?

In contrast to what you may think, its not only a bit harder, its much harder. Because |\vec r_1-\vec r_2| depends on all six components of the position vectors and so you should consider all the components in both the kinetic and potential terms. Just expand the position vectors w.r.t. their Cartesian components. So you'll have 6 EL equations.

 

[voko]

\frac{\partial L}{\partial r_1}=-\frac{\partial V(|{\vec{r_1}-\vec{r_2}}|)}{\partial r_1}=...?

This is not entirely correct. $\vec r_1$ was a the position of a body. What is $r_1$? Usually, that denotes the magnitude of the corresponding vector. You could use this as one of the generalized coordinates, but then you would need two other for the body, and what would those be?

I interpret your messages as if you really want to use Cartesian coordinates, in which case it is not $r_1$ that you should be using. It is $r_{11}, \ r_{12}, \ r_{13}$. Then you need to compute, for example, $$ \partial V(|\vec r_1 - \vec r_2|) \over \partial r_{11} . $$ The question is how?

Observe that $V(s)$ is a scalar function of one scalar argument, so it has some derivative $V'(s)$. Then you could use the chain rule and obtain $$ {\partial V(|\vec r_1 - \vec r_2|) \over \partial r_{11} } = V'(|\vec r_1 - \vec r_2|) {\partial |\vec r_1 - \vec r_2| \over \partial r_{11} }. $$

Can you continue from here?