Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Solution of Euler-Lagrange equation

  1. May 20, 2009 #1
    I have the following Lagrangian:

    [tex] \mathcal{L} = 1/2 \partial_{\mu} \varphi \partial^{\mu} \varphi - 1/2 b ( \varphi^{2} - a^{2} )^{2} [/tex], where [tex] a,b \in \mathbb{R}_{>0} [/tex] and [tex] \varphi [/tex] is a real (scalar) field and x are spacetime-coordinates.

    I calculated the Euler-Lagrange eq. and get: [tex] \square \varphi + 2 b ( \varphi^{2} - a^{2} ) \varphi = 0 [/tex]

    My problem is now to find constant solutions and static ones like [tex] \varphi(x) = f(x-x_{0}) [/tex] where [tex] x_{0} [/tex] is constant. But I don't know how to solve the differential equation above. Does anyone have an idea?
     
  2. jcsd
  3. May 20, 2009 #2
    Constant solutions are easy to calculate (there are 3 of them), as [itex]\square \varphi[/itex] is zero if [itex]\varphi[/itex] is constant.

    Now, for the other question, I'm assuming that for static solution you mean independent of time. If so, [itex]\square \varphi = \nabla^2 f[/itex], wich leads to the stationary NLS ecuation with constant potential, and there is plenty of literature on how to solve it.
     
  4. May 20, 2009 #3
    I recognise the solution if [itex]\varphi[/itex] is constant.

    [itex]\varphi = a[/itex] , [itex]\varphi = -a[/itex] and [itex]\varphi = 0[/itex]. :biggrin:
     
  5. May 22, 2009 #4
    Thank you for your replies!

    Remark: I think it should be: [tex] \square \varphi = - \nabla^{2} f [/tex]

    But isn't there any other (easier) approach to solve this differential equation? I thought of something like finding a first integral or something like that, but I can't. Or maybe there is a (simple) physical explanation which could lead to a solution???
     
    Last edited: May 22, 2009
  6. May 23, 2009 #5
    There should be a soliton-like solution arising from from energy conservation. For example, let [itex]x \in \mathbb{R}[/itex], then

    [tex]-f'' + 2b(f^2-a^2)f=0.[/tex]

    Multiplying by f' and integrating,

    [tex]-\frac{f'^2}{2} + b(\frac{f^2}{2}-a^2)f^2 = C,[/tex]

    where C is a constant of integration determined by the boundary conditions (i.e. [itex]f,f'\rightarrow 0[/itex] as [itex]x \rightarrow \infty[/itex]). The remainding equation is separable and solvable (in terms of elliptic functions maybe?), and if you study the face plane of the second order equation, you can prove that there is a soliton-like solution (depending of the existence of homoclinic or heteroclinic orbits).

    In more variables, you should be able to build a Hamiltonian and see if that approach gives you some additional info.
     
  7. May 29, 2009 #6
    Yeap, it is really difficult to expressed f in term of x. Soliton?
    wave that maintains its shape while it travels at constant speed. http://en.wikipedia.org/wiki/Soliton" [Broken] without equation!

    But what are face plane? homoclinic or heteroclinic orbits ? :cry:
     
    Last edited by a moderator: May 4, 2017
  8. May 29, 2009 #7
    If C=0 is not that hard to solve, the equation is separable and i think that the change of variable f = sech² x will do the trick. Even if the constant is different than zero, you can solve the integral in terms of elliptic functions (I believe).
    Regarding the phase plane or portrait (my native language betrayed me there), you can read all about it in a basic DE book like the one of Braun or Boyce, homoclinic and heteroclinic orbits can be seen in more advanced books, like the second volume of Hubbard's dynamical systems book.
     
    Last edited: May 29, 2009
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook