Is the EOM for a Lagrangian with a scalar operator always 2D phi = 0?

  • Context: Graduate 
  • Thread starter Thread starter RedX
  • Start date Start date
  • Tags Tags
    Lagrangian
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
5 replies · 2K views
RedX
Messages
963
Reaction score
3
If you have a Lagrangian of the form:

[tex]L=\phi \partial^2 \phi[/tex]

how would you derive its equation of motion? All the books seem to say to treat this Lagrangian as if it were only a function of the field, and not derivatives of the field.

So to calculate this they seem to do a product rule:

[tex]\partial^2 \phi+\phi \partial^2=0[/tex]

The latter term is somehow equal to the first term, so you get:

[tex]2\partial^2 \phi=0[/tex]

Is this generally true, that if you have some scalar operator D sandwiched between two fields:

[tex]L=\phi D \phi[/tex]

then the EOM is:

[tex]2 D \phi=0[/tex] ?
 
Physics news on Phys.org
RedX said:
If you have a Lagrangian of the form:

[tex]L=\phi \partial^2 \phi[/tex]

how would you derive its equation of motion? All the books seem to say to treat this Lagrangian as if it were only a function of the field, and not derivatives of the field.

So to calculate this they seem to do a product rule:

[tex]\partial^2 \phi+\phi \partial^2=0[/tex]

The latter term is somehow equal to the first term, so you get:

[tex]2\partial^2 \phi=0[/tex]

Yes, if you gain some experience in these calculations, you can use such things as a shortcut. However, if this is all new to you I suggest that you try writing everything out explicitly to try and get some insight in what is happening. First of all, what you are calling a Lagrangian is actually a Lagrangian density,
[tex]\mathcal L = \phi \partial_\mu \partial^\mu \phi[/tex]
(summation convention in use). This is however not a physical quantity, the relevant thing is the action,
[tex]S = \int d^4x \mathcal L[/tex].
This means that you are allowed to do partial integration and distribute the derivatives any way you like:
[tex]S = B_1 - \int d^4x (\partial_\mu\phi)(\partial^\mu\phi) = B_2 + \int d^4x (\partial_\mu \partial^\mu \phi) \phi[/tex]
The B1,2 are boundary terms, which are usually neglected (e.g. we assume that the field vanishes at the boundary of the integration; although sometimes we cannot make such assumptions or at least the boundary conditions are not trivial).

Anyway, hardly anyone derives the equations of motion directly from varying the Lagrangian; instead one usually takes the Euler-Lagrange equations
[tex]\frac{\partial L}{\partial \phi} - \partial_\mu \frac{\partial L}{\partial (\partial_\mu\phi)} = 0[/tex]
as a starting point. You should have seen this equation in a course on classical mechanics, if not, look it up somewhere (for example in Fowles & Cassiday's "Analytical Mechanics"). This equation is very important and very useful, and if you are going to do this kind of calculations more often you should derive them yourself at least once. In your case, if you do the partial integration once, you get
[tex]\mathcal L = -(\partial_\mu \phi)(\partial^\mu \phi) = - \eta^{\mu\nu} (\partial_\mu\phi)(\partial_\nu\phi)[/tex].
This does not depend on [itex]\phi[/itex] but it depends on [itex]\partial_\mu\phi[/itex]. Let me rename the dummy summation index:
[tex]\mathcal L = - \eta^{\rho\sigma} (\partial_\rho\phi)(\partial_\sigma\phi)[/tex].
If you do the differentiation, you will have to use the product rule and take into account that
[tex]\frac{\partial_\lambda \phi}{\partial_\mu\phi} = \delta_\lambda^\mu[/tex];
then in the end you will find for the equations of motion
[tex]0 - \partial_\mu( 2 \partial^\mu\phi) = -\partial^2 \phi = 0.[/tex]
(Again, I urge you to do the calculation and write out all the steps; it's actually quite wonderful how it all fits together).

RedX said:
Is this generally true, that if you have some scalar operator D sandwiched between two fields:

[tex]L=\phi D \phi[/tex]

then the EOM is:

[tex]2 D \phi=0[/tex] ?
Yes, this is generally true and if D is an order 2 differential operator, like [itex]\partial^2\phi[/itex], then you can easily prove this by varying [itex]\phi \to \phi + \delta\phi[/itex] (where [itex]\delta\phi^2 = 0[/itex]) and doing some partial integrations. In fact there is a whole theory using so-called Greens functions, which can be applied to such cases. If you take a good course on quantum field theory, you will be taught about them (actually I first learned about them in a course called "topics in mathematical physics"). However, until you acquire some fluency doing these calculations I suggest sticking to the Euler-Lagrange formalism (derive the Euler-Lagrange equations so you don't have to vary the action every time you get a Lagrangian) so you can see in each case how the equations of motion appear.
 
Thanks. So the trick is to integrate by parts first, and then discard boundary terms from the Lagrangian. I was aware of the equations of motion, but couldn't quite figure out how to apply them to [tex]\phi \partial^2 \phi[/tex], since naively taking the derivative with respect to the derivative of the field gives you [tex]\phi \partial[/tex] acting on nothing.
 
Keep remembering that [itex]\partial^2[/itex] is shorthand for [itex]\partial_\mu \partial^\mu[/itex], so you would get [itex]\phi \partial^\mu[/itex]
I suppose you can show that for higher derivatives, you can extend the Euler-Lagrange equations to
[tex]\frac{\partial L}{\partial \phi} - \partial_\mu \frac{\partial L}{\partial (\partial_\mu \phi)} + \partial_\mu \partial_\nu \frac{\partial L}{\partial (\partial_\mu \partial_\nu \phi)} - \partial_\mu \partial_\nu \partial_\lambda \frac{\partial L}{\partial (\partial_\mu \partial_\nu \partial_\lambda \phi)} + \cdots = 0[/tex]
and apply it directly to the Lagrangian written as [itex]\mathcal L = g^{\mu\nu} \phi \partial_\mu \partial_\nu \phi[/itex].
Of course, to show the extended version of the E-L equations, you still have to perform some partial integrations so this does not "circumvent" the boundary conditions.