Lecture 22a: Nonlinear Systems, Period Doubling, Chaos

Michael Fowler

This is a preliminary version: more will be added. It should be read along with the applet.


This lecture is an extension of lecture 22, which followed Landau’s treatment of the driven damped nonlinear oscillator.  In that lecture, after reviewing Landau’s analysis, we illustrated and confirmed his conclusions with an applet (which also gave a guide to the precision of his predictions, they were within 10% to 20%).  Here we shall use essentially the same applet, suitably modified, to look at related systems that have spectacular new features: period doubling, and approach to chaos. This is now a huge subject, we provide only a glimpse here. We follow part of the treatment in Taylor’s Classical Mechanics, with the addition of our applets for illustrative purposes.   We verify the accuracy of our applet by reproducing some of Taylor’s graphs.

Driven Damped Pendulum: The Road to Chaos

After the anharmonic oscillator treated in the previous lecture, our second example of a one-dimensional nonlinear system is a driven damped pendulum: that just means replacing the potential term ω 0 2 x  in the linear oscillator with ω 0 2 sinx,  or rather ω 0 2 sinϕ,  to make clear we have an angular system.  

For small amplitude oscillations, the free pendulum’s frequency is almost independent of amplitude, as has been known for centuries, and the first correction on increasing the amplitude is from a quartic addition to the potential, just as in the nonlinear case discussed earlier, except now with a definite value, and an opposite sign, so the (undamped, undriven) frequency decreases with amplitude.  

Provided the amplitude remains small, the driven damped pendulum behaves like the driven damped linear oscillator we discussed earlier.  But things get much more interesting as the amplitude increases.

For the driven pendulum, the natural measure of the driving force is its ratio to the weight mg,  Taylor calls this the drive strength, so for driving force

F t = F 0 cosωt,

the drive strength

γ= F 0 /mg, ω 0 = g/L .

The equation of motion (with resistive damping force bv  and hence resistive torque b L 2 ϕ ˙  )  is:

m L 2 ϕ ¨ =b L 2 ϕ ˙ mgLsinϕ+LF t ,

Dividing by m L 2  and writing the damping term b/m=2β  (to coincide with Taylor’s notation, his equation 12.12) we get

ϕ ¨ +2β ϕ ˙ + ω 0 2 sinϕ=γ ω 0 2 cosωt.

(Compare with our equation in the previous lecture for the anharmonic oscillator in Landau notation:   x ¨ +2λ x ˙ + ω 0 2 x= f/m cosγtβ x 3 .  )

Behavior on Gradually Increasing the Driving Force: Period Doubling

The driving force γ= F 0 /mg  is the dimensionless ratio of drive strength to weight, so if this is small the pendulum will not be driven to large amplitudes, and indeed we find that after initial transients it settles to motion at the driving frequency, similar to the linearized case. We would expect things to get more complicated when the oscillations have amplitude of order a radian, meaning a driving force comparable to the weight.  And indeed they do.

Here we’ll show that our applet reproduces the sequence found by Taylor as the driving strength is increased.  We therefore choose his values ω 0 =1.5,2β=0.75,ω=1,  and gradually increase γ  from 0.9 to 1.0829, where chaos begins.

For γ=0.9,  the oscillation (after brief initial transients) looks like a sine wave, although it’s slightly flatter, and notice the amplitude (turquoise in box) is greater than π/2,  the positive and negative swings are equal in magnitude to five figure accuracy after five or so oscillations.

For γ=1.06,  there is a larger initial transient, vaguely resembling the three-cycle we’ll meet shortly.  In fact, this can be suppressed by starting at φ=π/2,  but there are still transients: the peaks are not uniform to five-figure accuracy until about forty cycles.


For γ=1.0662,  there are very long transients: not evident on looking at the graph, but revealed by monitoring the amplitude (turquoise figures in the box on the graph) which reads the value at the latest extremum. 

To get past these transients, set the applet speed to 50 (this does not affect accuracy, only prints a lot less dots), run 4500 cycles (won’t take long!) then pause and go to speed = 1.  You will find the peaks now all have the same height to the five-figure accuracy displayed.

Now repeat for γ=1.0663  (use the button to nudge up by one in that last digit).  Go at speed 50 to 3000 cycles, then you’ll find on pausing and going back to speed 1 that the peaks have period 2, with a 0.01 radian difference, approximately.  But this is a false alarm: it’s an incredibly long transient. The system finally settles at period 1 (to five figures) after 12,000 cycles.  We are apparently right at the bifurcation point.

Going to 1.0664 and beyond you’ll see both peaks and dips separating, for a well-defined period 2, which lasts until 1.0792.


For γ=1.0793,  the period doubles again: successive peaks are now of order 0.4 radians apart, but successive “high” peaks are suddenly about 0.01 radians apart.

For γ=1.0821  there is a further doubling to an 8 cycle, then at 1.0827 to a 16 cycle.

Look at the graph for 1.0826: in particular, the faint red horizontal line near the bottom. Every fourth dip goes below the line, but the next lowest dips alternate between touching the line and not quite reaching it.  This is the 8 cycle.


The intervals between successive doublings decrease by the factor δ=4.6692,  found universally in period doubling cascades, and called the Feigenbaum number, after its discoverer. Our five-figure accuracy is too crude to follow this sequence further, but we can establish (or at least make very plausible!) that beyond the geometric series limit at γ c =1.0829  there is no periodicity, the system is chaotic. (Of course, the values of γ  depend on the chosen damping parameter, etc., only the ratio of doubling intervals is universal.)  However, the full picture is incredibly complex: there are further intervals of periodicity, for example a 6 cycle at γ=1.0845,  pictured below.

Different Attractors

The periodic solutions described above are called “attractors”: configurations where the system settles down after initially wandering around. 

Clearly the attractors change with the driving strength, what is less obvious is that they may be different for different initial conditions.  Taylor shows that for γ=1.077  taking φ 0 =0  gives a 3-cycle after transients, but φ 0 =π/2  gives a 2-cycle. (Easily checked with the applet!)

 Looking at the initial wanderings, which can be quite different for very small changes in the driving strength (compare 1.0730 to 1.0729 and 1.0731, use speed 5, it doesn’t affect the accuracy). But you can see these initial wanderings include elements from both attractors.

Exercise:  Use the applet to find at what angle the transition from one attractor to the other takes place.  And, explore what happens over a range of driving strengths.

These are the simplest attractors: there are far more complex entities, called strange attractors, we may discuss later.

Exercises:  Try different values of the damping constant and see how this affects the bifurcation sequence.