Equilibrium solutions in double well potential duffing oscillator

Ratpigeon
Messages
52
Reaction score
0

Homework Statement



I am trying to show that for a duffing oscillator described by
x''+2g x'+ax+bx^3=0
with a<0, b>0
the equilibria at x=+- \sqrt{-a/b} are stable

Homework Equations


I used coupled equations, and the characteristic equation of a linear system

The Attempt at a Solution



Coupled equations, x and x' related by
d/dt(x,x&#039;)=[[0,1],[-a-bx^2, -2g]](x,x&#039;)

Setting x=+-\sqrt{-a/b} gives the characteristic is

\lambda^2+2g \lambda+a+b(-a/b)=0
But this can't be right since, i know from literature that the eigens are
\lambda 1=-g + \sqrt{g^2+2a}
\lambda 2=-g- \sqrt{-g^2+2a}
and this gives eigens 0 and -2g
 
Last edited:
Physics news on Phys.org
After you have reduced the equation to the system and found the equilibrium point, you have to linearize the system about the equilibrium point, and then find the eigenvalue of the linear system.
 
Since no one has chimed in yet, I'll throw in my 2 cents (and it might not be worth even that). I believe that one way to examine the behavior near an equilibrium point is to "linearize" the equation in the neighborhood of the equilibrium point. So, you could let z = x-xo where xo is an equilibrium point. Rewrite the differential equation in terms of z but keep only terms linear in z. Then find the characteristic equation for this linearized differential equation.

[Edit: Wow, I did not see voko's response. Sorry.]
 
I'd done all the linearising stuff - I found my problem; I needed to use a Jacobean matrix, instead of just a system
 
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top