Register to reply

Another question on Lagrangians (sorry)

by vertices
Tags: lagrangians
Share this thread:
vertices
#1
Dec31-09, 10:08 AM
P: 230
Sorry about the endless stream of questions about Lagrangians. I am actually beginning to detest them a bit;p

Anyway, if we have a Lagrangian in three dimensional space:

[tex]L=\frac{1}{2}m\dot{\vec{x}}^{2}+e\vec{A}.\dot{\vec{x}}[/tex]

where [tex]A_{i}=\epsilon_{ijk}B_{j}x_{k}[/tex] and B is just a constant (magnetic field).

The question is to find the equations of motion. To do this we can just use the Euler-Lagrange equations thus:

[tex]\frac{{\partial}L}{{\partial}q^{i}}-\frac{d}{dt}\frac{{\partial}L}{{\partial}{\dot{q}}^{i}}=0[/tex]

Now,

[tex]\frac{{\partial}L}{{\partial}x^{i}}=0[/tex]

for all i.

So the EOM are given by

[tex]\frac{d}{dt}\frac{{\partial}L}{{\partial}{\dot{x}}^{i}}=0[/tex]

Now
[tex]\frac{{\partial}L}{{\partial}{\dot{x}}^{i}}=m\dot{x}^{i}+eA_{i}[/tex]

so d/dt of this should just be:

[tex]\frac{d}{dt}\frac{{\partial}L}{{\partial}{\dot{x}}^{i}}=m\ddot{x}^{i}+e \dot{A}_{i}[/tex]

Since [tex]A_{i}=\epsilon_{ijk}B_{j}x_{k}\Rightarrow\dot{A}_{i}=\epsilon_{ijk}B_{j }\dot{x}_{k}[/tex]

so we should get this:

[tex]\frac{d}{dt}\frac{{\partial}L}{{\partial}{\dot{x}}^{i}}=m\ddot{x}^{i}+e \epsilon_{ijk}B_{j}\dot{x}_{k}[/tex]

but according to the solution to the problem it should be:

[tex]\frac{d}{dt}\frac{{\partial}L}{{\partial}{\dot{x}}^{i}}=m\ddot{x}^{i}+2 e\epsilon_{kji}B_{j}\dot{x}_{k}[/tex]

Note the factor of 2.

I posted a similar problem and I was lead to believe that I need to differentiate wrt another superscript - I am not sure if this is at all relavent here since we are differentiating wrt time.

Would be really grateful if someone could offer me some advice...

Thanks..
Phys.Org News Partner Physics news on Phys.org
New approach to form non-equilibrium structures
Nike krypton laser achieves spot in Guinness World Records
Unleashing the power of quantum dot triplets
diazona
#2
Dec31-09, 02:14 PM
HW Helper
P: 2,155
Seems like
[tex]\frac{\partial L}{\partial x^i} \neq 0[/tex]
because the Lagrangian depends on x through the vector potential A.
vertices
#3
Dec31-09, 03:19 PM
P: 230
Quote Quote by diazona View Post
Seems like
[tex]\frac{\partial L}{\partial x^i} \neq 0[/tex]
because the Lagrangian depends on x through the vector potential A.

thanks diazona. You are right:

[I will switch to subscripts, I hope that is okay]

[tex]
L_{i}=\frac{1}{2}m\dot{x}_{i}}^{2}+e\vec{A}.\dot{\vec {x}}
[/tex]

where

[tex]
e\vec{A}.\dot{\vec {x}}=[\epsilon_{123}B_{2}x_{3}+\epsilon_{132}B_{3}x_{2}]\dot{x}_{1} + [\epsilon_{213}B_{1}x_{3}+\epsilon_{231}B_{3}x_{1}]\dot{x}_{2} + [\epsilon_{312}B_{1}x_{2}+\epsilon_{321}B_{2}x_{1}]\dot{x}_{3} [/tex]

So take

[tex]\frac{\partial L}{\partial x_1} = [\epsilon_{231}B_{3}\dot{x}_2+\epsilon_{321}B_{2}\dot{x}_3][/tex]

Now,

[tex]
\frac{d}{dt}\frac{{\partial}L}{{\partial}{\dot{x}_1}}
[/tex]

so feeding these in to the EL equation we get:

[tex]
m\ddot{x}_{1}+e(\epsilon_{123}B_{2}\dot{x}_3+\epsilon_{132}B_{3}\dot{x} _3) - [ \epsilon_{231}B_{3}\dot{x}_2+\epsilon_{321}B_{2}\dot{x}_3 ] = 0
[/tex]

I can't see why we end up up

[tex]
m\ddot{x}_{i}+2e\epsilon_{kji}B_{j}\dot{x}_{ k}
[/tex]

diazona
#4
Dec31-09, 04:54 PM
HW Helper
P: 2,155
Another question on Lagrangians (sorry)

No problem with the subscripts By the way, you're missing a factor of [itex]e[/itex] (charge) in the last term there.

Anyway, big hint: what happens when you put in numbers for the [itex]\epsilon_{ijk}[/itex]'s?
vertices
#5
Dec31-09, 05:18 PM
P: 230
Quote Quote by diazona View Post
No problem with the subscripts By the way, you're missing a factor of [itex]e[/itex] (charge) in the last term there.

Anyway, big hint: what happens when you put in numbers for the [itex]\epsilon_{ijk}[/itex]'s?
Yes, it works! (I think I made a mistake, which was confusing me somewhat!) Many thanks Diazona:)

Is decomposing into components the only way to see the levi-cevita (if that's what you call it) connection?
diazona
#6
Dec31-09, 06:32 PM
HW Helper
P: 2,155
I'm not really sure. I guess you could try it in general tensor notation:
[tex]L = \frac{1}{2}m \delta_{jk} \dot x^j\dot x^k + e \epsilon_{ijk}B^i x^j \dot x^k[/tex]
then
[tex]\frac{\partial L}{\partial x^j} = e\epsilon_{ijk} B^i \dot x^k = -e\epsilon_{ikj} B^i \dot x^k[/tex]
and
[tex]\frac{\mathrm{d}}{\mathrm{d}t}\frac{\partial L}{\partial \dot x^k} = m \delta_{jk} \ddot x^j + e\epsilon_{ijk} B^i \dot x^j[/tex]
so when you set those equal to each other (or subtract one from the other), you wind up with
[tex]m \delta_{jk} \ddot x^j + 2e\epsilon_{ijk} B^i \dot x^j = 0[/tex]

You could also make an attempt at it in vector notation,
[tex]L = \frac{1}{2}m (\dot{\vec{x}})^2 + e(\vec{B}\times\vec{x})\cdot\dot{\vec{x}}[/tex]
but you'd have to use a vector identity,
[tex](\vec{A}\times\vec{B})\cdot\vec{C} = (\vec{B}\times\vec{C})\cdot\vec{A} = (\vec{C}\times\vec{A})\cdot\vec{B}[/tex]
so you could write
[tex]L = \frac{1}{2}m (\dot{\vec{x}})^2 + e(\dot{\vec{x}}\times\vec{B})\cdot\vec{x}[/tex]
Then you get
[tex]\frac{\partial L}{\partial \vec{x}} = e \dot{\vec{x}}\times\vec{B} = -e \vec{B}\times\dot{\vec{x}}[/tex]
and
[tex]\frac{\mathrm{d}}{\mathrm{d}t}\frac{\partial L}{\partial \dot{\vec{x}}} = \frac{\mathrm{d}}{\mathrm{d}t}\left[m \dot{\vec{x}} + e \vec{B}\times\vec{x}\right] = m \ddot{\vec{x}} + e \vec{B}\times\dot{\vec{x}}[/tex]
and again, when you put them together, you get
[tex]m \ddot{\vec{x}} + 2e \vec{B}\times\dot{\vec{x}} = 0[/tex]

It's the same math going on underneath, it's just that there are several different ways to write it out. If I didn't already have a lot of experience with using vector or tensor notation, I'd probably do it the way you did it, manually, one component at a time.
vertices
#7
Dec31-09, 07:42 PM
P: 230
Hi Diazona - this is really, really helpful:)

Can I ask a really stupid question though:

How do you get from this (I've omitted the term in B):

[tex]
L = \frac{1}{2}m \delta_{jk} \dot x^j\dot x^k[/tex]

to this (this what the third line implies, ie. before d/dt'ing it):

[tex]
\frac{\partial L}{\partial \dot x^k} = m \delta_{jk} \dot x^j[/tex]

Where has the factor of 1/2 gone?

Thanks again:)
diazona
#8
Jan1-10, 04:55 PM
HW Helper
P: 2,155
Ah, well in that one you have to remember that j and k are dummy indices, so that [itex]\dot{x}^j[/itex] and [itex]\dot{x}^k[/itex] are really the same thing. So the derivative applies to both of them equally.

If you want to show it explicitly, you can apply the product rule. I'm going to relabel the index of the [itex]\dot{x}^k[/itex] in the derivative so there's no confusion (of course, the label makes no difference because it is a dummy index):
[tex]\frac{\partial}{\partial \dot{x}^a}\biggl[\frac{1}{2}m\delta_{jk}\dot{x}^j\dot{x}^k\biggr] = \frac{1}{2}m\delta_{jk}\frac{\partial \dot{x}^j}{\partial \dot{x}^a}\dot{x}^k + \frac{1}{2}m\delta_{jk}\dot{x}^j\frac{\partial \dot{x}^k}{\partial \dot{x}^a}[/tex]
Now, do you accept that [itex]\partial \dot{x}^j/\partial \dot{x}^a = \delta^{j}_{a}[/itex]? Then
[tex]\begin{align*}
\frac{\partial}{\partial \dot{x}^a}\biggl[\frac{1}{2}m\delta_{jk}\dot{x}^j\dot{x}^k\biggr] &= \frac{1}{2}m\delta_{jk}\delta^{j}_{a}\dot{x}^k + \frac{1}{2}m\delta_{jk}\dot{x}^j\delta^{k}_{a} \\
&= \frac{1}{2}m\delta_{ak}\dot{x}^k + \frac{1}{2}m\delta_{ja}\dot{x}^j \\
&= \frac{1}{2}m\delta_{ja}\dot{x}^j + \frac{1}{2}m\delta_{ja}\dot{x}^j \\
&= m\delta_{ja}\dot{x}^j
\end{align*}[/tex]
I went from [itex]\delta_{ja}\dot{x}^j[/itex] to [itex]\delta_{ak}\dot{x}^k[/itex] just by relabeling that index from j to k, and then switching the order of indices in the Kronecker delta because it's a symmetric tensor.
vertices
#9
Jan1-10, 05:45 PM
P: 230
Aahh that is clear now - it's so much better knowing how things work, than just accepting them as facts... thanks ever so much diazona:)

I have one final question (I promise I won't bother you anymore after this!) - Infact, it's in this thread, but I'll reproduce it here:

If we have the following Lagrangian:

[tex]
L={\frac{m}{2}}{g_{ij}(x)}.{\dot{x^{i}}{\dot{x^{j} }
[/tex]

I want to show the momentum is this:

[tex]
(\frac{{\partial}L}{{\partial}\dot{x^{i}}})=mg_ {ij}x^{j}

[/tex]

So I differentiate wrt another superscript (as per George's advice):

[tex]
\frac{{\delta}L}{d \dot{x}^k}=\frac{m}{2}g_{ij}[\frac{d \dot{x}^i}{d \dot{x}^k}{\dot{x}^j+\frac{d \dot{x}^j}{d \dot{x}^k}{\dot{x}^i]
[/tex]

This is equal to:

[tex]
\frac{m}{2}g_{ij}[{\delta}_{ik}\dot{x}^{j}+{\delta}_{jk}\dot{x}^{i}]
[/tex]

Now replacing k with i we find

[tex]
\frac{{\delta}L}{d \dot{x}^i}=\frac{m}{2}g_{ij}\dot{x}^{j}
[/tex]

ie. still have a factor of 1/2. My concern is that we can't relabel indices (as you did above), as there are no repeated indices.

Can you shed any light on where I have gone wrong?

Thanks..


Register to reply

Related Discussions
Rigid Body Lagrangians Advanced Physics Homework 0
Lagrangians and P-Conservation Advanced Physics Homework 2
Quantization without Lagrangians. General Physics 19