Period of anharmonic oscillator

  • Context: Undergrad 
  • Thread starter Thread starter ChrisVer
  • Start date Start date
  • Tags Tags
    Oscillator Period
Click For Summary
SUMMARY

The discussion centers on the period of an anharmonic oscillator with a quartic potential, defined as V(x) = x^4. The period T is calculated using the integral T = √(8m) ∫(0 to a) (dx/√(V(a) - V(x))). As the amplitude a approaches zero, the period diverges, leading to T approaching infinity due to the flatness of the potential at small x values. The conversation highlights the limitations of using harmonic oscillator approximations for potentials like V(x) = x^4, emphasizing that such approximations fail when the second derivative at the minimum is zero.

PREREQUISITES
  • Understanding of classical mechanics and oscillatory motion
  • Familiarity with potential energy functions and their properties
  • Knowledge of Taylor series expansions and their applications in physics
  • Experience with numerical integration techniques for solving differential equations
NEXT STEPS
  • Explore the numerical methods for solving integrals in classical mechanics
  • Study the implications of Taylor series expansions for various potential functions
  • Investigate the behavior of oscillators under different potential energy profiles
  • Learn about the mathematical treatment of singularities in physical systems
USEFUL FOR

Physicists, students of classical mechanics, and researchers interested in the dynamics of anharmonic oscillators and potential energy landscapes.

ChrisVer
Science Advisor
Messages
3,372
Reaction score
465
Well I numerically solved for the potential V(x)=x^4, the period:
\begin{equation}
T = \sqrt{8m} \int_0^a \frac{dx}{\sqrt{V(a) - V(x)}}
\end{equation}

where a was the amplitude of the oscillation and m the mass of the particle.
Nevertheless, what I was asked to plot was the above period T(a) for a\in [0.,2.].
My problem however is that the period of small values of a, diverges. I was able to see how this is the case mathematically, by expanding the square root and obtaining an expression for T that goes for V(x)=x^p as:
T \sim \frac{1}{a^{p/2-1}}
which reproduces also the V=x^2 result of constant period.

However I cannot picture what is happening physically. Does it have to do with the flatness of the potential at so small a's and so x's?
Any idea?
 
Physics news on Phys.org
Imagine the potential curve as a literal valley, with the particle a small object sliding along the surface.
You start the object stationary at x=a, let it go, and measure how long it takes to return... and call that number T.
What happens to T as the start position gets closer to zero?
 
  • Like
Likes   Reactions: ChrisVer
Simon Bridge said:
Imagine the potential curve as a literal valley, with the particle a small object sliding along the surface.
You start the object stationary at x=a, let it go, and measure how long it takes to return... and call that number T.
What happens to T as the start position gets closer to zero?

yup, the potential is pretty flat as you move to smaller x...so the body will "hardly" slide...
in explicitly a\rightarrow0 \Rightarrow T \rightarrow \infty

But then again, I remember some quote saying that any potential around a (local) minimum can be described to a good approximation by an Harmonic Oscillator potential with the minimum at that point... (Taylor expanding the potential around the local minimum)...What about its application to this case?
 
The other interpretation to starting at the minimum point is that T=0 because it take no time at all to return to it's start position.
So there is a singularity there :)

An arbitrary potential can be expanded in a Taylor series, and it is common to take low order's of that series as an approximation.
So have you taken the taylor expansion for your potential to see what the first few terms are?
 
Simon Bridge said:
So have you taken the taylor expansion for your potential to see what the first few terms are?
Oh I see... I tried to but I was carried away and didn't think that the 2nd order of the expansion that I [only] kept was also vanishing.
so it's not a rule for all types of potentials, rather than those who move slower than V''(x_0) \ne 0...
so from power potentials V(x)= x^p no p \ge 3 can be described by an harmonic oscillator...
 
V(x) = x^4

expanding as a taylor series about x=a:

##V(x) = 4a^3(x-a) + 6a^2(x-a)^2 + 4a(x-a)^3 + 1(x-a)^4 + 0 + 0 + \cdots##... put a=0.
##V(x) = 0 + 0 + 0 + x^4 + 0 + 0 + \cdots## ... surprise surprise.

This is telling you that the HO approximation to the quartic potential is V=0.
In this approximation: T=0 or infinity everywhere.
The approximation is valid for x "sufficiently close to" 0 ... how close is "sufficient" depends on the sensitivity of your instruments.
It's like how the approximation that the Earth is flat is really good inside a typical room but not if you are looking from orbit.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 4 ·
Replies
4
Views
628
  • · Replies 20 ·
Replies
20
Views
4K
  • · Replies 17 ·
Replies
17
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 41 ·
2
Replies
41
Views
4K