Resonance and Fractal Geometry

The phenomenon of resonance will be dealt with from the viewpoint of dynamical systems depending on parameters and their bifurcations. Resonance phenomena are associated to open subsets in the parameter space, while their complement corresponds to quasi-periodicity and chaos. The latter phenomena occur for parameter values in fractal sets of positive measure. We describe a universal phenomenon that plays an important role in modelling. This paper gives a summary of the background theory, veined by examples.


What Is Resonance?
A heuristic definition of resonance considers a dynamical system, usually depending on parameters, with several oscillatory subsystems having a rational ratio of frequencies and a resulting combined and compatible motion that may be amplified as well. Often the latter motion is also periodic, but it can be more complicated as will be shown below. We shall take a rather eclectic point of view, discussing several examples first. Later we shall turn to a number of universal cases, these are context-free models that occur generically in any system of sufficiently high-dimensional state and parameter space.
Among the examples are the famous problem of Huygens's synchronizing clocks and that of the Botafumeiro in the Cathedral of Santiago de Compostela, but also we briefly touch on tidal resonances in the planetary system. As universal models we shall deal with the Hopf-Neȋmark-Sacker bifurcation and the Hopf saddle-node bifurcation for mappings. The latter two examples form 'next cases' in the development of generic bifurcation theory. The term universal refers to the context independence of their occurrence: in any with certain Left: the periodically driven harmonic oscillator (1) and right: the driven Duffing oscillator (2) generic specifications these bifurcations occur in a persistent way. We witness an increase in complexity in the sense that in the parameter space the resonant phenomena correspond to an open & dense subset union of tongues, while in the complement of this a nowhere dense set of positive measure exists, corresponding to multi-or quasi-periodic dynamics. This nowhere dense set has a fractal geometry in a sense that will be explained later. From the above it follows that this global array of resonance tongues and fractal geometry has a universal character. As we shall see, both locally and globally Singularity Theory can give organizing principles. It should be noted at once that next to periodic and quasi-periodic dynamics also forms of chaotic dynamics will show up.
Remark In many cases the resonant bifurcations are repeated at ever smaller scales inside the tongues, leading to an infinite regress. Then we have to extend the notion of open & dense to residual and that of nowhere dense to meagre. Here a residual set contains a countable intersection of open & dense sets, while a meagre set is a countable union of nowhere dense sets. One sometimes also speaks in terms of G δ -or F σ -sets, respectively [73].

Periodically Driven Oscillators
Many of the resonant phenomena of interest to us are modelled by periodically driven or coupled oscillators. To fix thoughts we now present two examples, namely the harmonic and the Duffing oscillator subjected to periodic forcing.

The Driven Harmonic Oscillator
One of the simplest class-room examples of resonance occurs in the harmonic oscillator with periodic forcingẍ = −α 2 x − cẋ + ε sin(ωt) (1) with x, t ∈ R, where c > 0 is the damping and where ε ≥ 0 controls the size of the forcing. The forcing has frequency ω > 0 and therefore has period T = 2π/ω. The 'response' solution is periodic with this same frequency. Its amplitude and phase are given by R = ε (ω 2 − α 2 ) 2 + c 2 ω 2 and tan φ = cω ω 2 − α 2 .
Fixing α, c and ε we consider R as a function of ω, see the amplitude-response diagram in Fig. 1 (left) where the response amplitude R is plotted as a function of ω.

The Driven Duffing Oscillator
As a non-linear variation on the above we consider the Duffing equation with periodic forcingẍ = −α 2 x − cẋ − δx 3 + ε sin(ωt), (2) where α, ω, c, ε and are positive and where δ with 0 < δ 1 is considered as a perturbation parameter. Fixing α, c and ε as before, by successive approximation of a solution we obtain the following relation between ε and R: as δ ↓ 0, with a similar approximation for the phase φ, compare with Stoker [84], also see [44]. In Fig. 1 (right) we depict the corresponding curve in the (ω, R)-plane, which now no longer is a graph.

Remarks
− One of the exciting things about resonance concerns the peaks of the amplitude R that can be quite high, even where ε is still moderate. Systems like (1) and (2) form models or metaphors for various resonance phenomena in daily life. In many cases high resonance peaks one needs to 'detune' away from the resonance value corresponding to the peak, think of a marching platoon of soldiers that have to go out of pace when crossing a bridge. − In other cases, like when 'tuning' the radio receiver to a certain channel, one takes advantage of the peak. − It should be noted that the nonlinear (2) dynamically is far richer than the linear case (1), e.g., see [58] and references therein.

Geometrical Considerations
In both cases of the driven oscillators witnessed above the state space is with a(x, y, δ) = −α 2 x − cy − δx 3 , where δ = 0 in the harmonic example. The second factor is the circle T 1 = R/(2πZ) which takes into account the periodicity of the systems in z. The response motion of the form x(t) = R sin(ωt + φ) then corresponds to a closed curve ⎛ This closed curve, when projected onto the (x, y)-plane, exactly forms an ellipse.

Fig. 2
Huygens's synchronizing clocks [61] We now can describe this motion in terms of a 2-dimensional torus T 2 = T 1 × T 1 , parametrized by two variables ϕ 1 and ϕ 2 in a system of differential equations of the formatφ where ϕ 1 = z and ω 1 = ω and where for ϕ 2 we take the phase of the motion on the ellipse. Thus ϕ 2 exactly is the time parametrization of this motion, scaled to the period 2π : in this case ω 2 = ω 1 = ω. Therefore in this way, the curve (4) can be seen as a 1 : 1 torus knot.
To view resonant motion in terms of torus dynamics turns out to be extremely useful and this can also be applied to coupled oscillators. Here a classical example is given by Christiaan Huygens [61], who in 1665 observed the following phenomenon, see Fig. 2. Two nearly identical pendulum clocks mounted on a not completely rigid horizontal beam tend to synchronize. Moreover, when the pendula both move in the vertical plane through the beam, they have a tendency to synchronize in anti-phase motion. A simple model describes this system in the format (5), where the angles ϕ 1 and ϕ 2 are the phases of the two oscillators and where again ω 1 = ω 2 . Later on we will come back to this and other examples where we will also see other frequency ratios ω 1 : ω 2 .

Torus Flows and Circle Mappings
In this section we turn to the dynamics on the 2-dimensional torus T 2 from Sect. 1.1.3 for its own sake, introducing the weakly coupled systeṁ 3 Poincaré mapping of a torus flow [19] Here f 1 and f 2 are 2π -periodic functions in both variables. Also we use a parameter ε to control the strength of the coupling. For ε = 0 we retrieve the format (5). If T 1 and T 2 are the respective periods of oscillation, then ω 1 = 2π/T 1 and ω 2 = 2π/T 2 . We first define the Poincaré mapping from the circle T 1 to itself and then introduce the rotation number.

The Poincaré Mapping
If in (6) the size |ε| of the coupling is not too large, a first-return Poincaré mapping is defined, as we shall explain now, also see Fig. 3. Without restricting generality we can take the generating circle T 1 × {0} for simplicity baptising ϕ = ϕ 1 . In that case we take the integral curve with initial state (ϕ 1 , ϕ 2 ) = (ϕ, 0) and follow this till (ϕ 1 , ϕ 2 ) = (P (ϕ), 0), counting mod 2πZ. It is easy to see that P − Id should be a periodic function in ϕ, which gives P the general format where α = ω 2 /ω 1 and where f is a 2π -periodic function. Consideration of the T 1 -dynamics generated by iteration of P gives a lot of information about the original T 2 -flow, in particular its asymptotic properties as t → ∞. For instance, a fixed point attractor of P corresponds to an attracting periodic orbit of the flow which form a 1 : 1 torus knot as we saw at the end of Sect. 1.1. Similarly a periodic attractor of P of period q corresponds to an attracting periodic orbit of the flow. In general, periodicity will be related to resonance, but to explain this further we need the notion of rotation number.

Rotation Number
For orientation-preserving homeomorphisms P : T 1 −→ T 1 Poincaré has left us the extremely useful concept of rotation number (P ), which describes the average amount of rotation as follows: HereP : R −→ R is a (non-unique) lift of P which makes the diagram commute, where pr : R −→ T 1 is the natural projection ϕ → e iϕ . This means that in the formula (9) we do not count modulo 2π , but keep counting in R. From [53,72] we quote a number of properties of (P ): 3. If P : ϕ → ϕ + 2πα is a rigid rotation then (P ) = α mod Z. 4. (P ) ∈ Q precisely when P has a periodic point. Moreover, (P ) = p/q with p and q relatively prime corresponds to a p : q torus knot. 5. If P is of class C 2 and (P ) = α for α ∈ R \ Q, then, by a result of Denjoy, the mapping P is topologically conjugated to the rigid rotation ϕ → ϕ + 2πα.
Recall that in that case any orbit {P n (ϕ)} n∈Z forms a dense subset of T 1 . The corresponding dynamics is called quasi-periodic. 6. If P depends continuously on a parameter, then so does (P ).

The Arnold Family of Circle Mappings
A famous example is formed by the Arnold family Periodicity It is instructive to consider its fixed points, given by the equation A brief graphical inspection reveals that mod 2πZ this equation has exactly two solutions for |ε| > 2π|α|.
In Fig. 4 this region, bounded by the two straight lines ε = ±2πα, is depicted for ε > 0. It is not hard to see that one of the fixed points is attracting and the other repelling. At the boundary |ε| = 2π|α| these annihilate one another in a saddle-node bifurcation. For the entire region |ε| ≥ 2π|α| one has (A α,ε ) = 0 mod Z. This region is called the Arnold tongue of rotation number 0. From the properties of Sect. 1.2.2 it follows that, for (α, ε) = (p/q, 0) with p and q relatively prime, one has (A α,ε ) = p/q. One can show that from each (α, ε) = (p/q, 0) an Arnold tongue emanates, in which for all the parameter points (α, ε) one has (A α,ε ) = p/q, see Fig. 4. The 'sharpness', i.e., the order of contact of the boundaries of the p/q-tongue at (α, ε) = (p/q, 0) exactly is of order q, see [1,3,24].
Fixing ε = ε 0 > 0 small, we consider the graph of α → (A α,ε 0 ). By another general property of Sect. 1.2.2, this function is continuous. Moreover, for every rational value p/q it is constant on some plateau, corresponding to the p/q-tongue, see Fig. 4. The total result is a devil's staircase as depicted in Fig. 5.

Quasi-periodicity
In between the tongues the rotation number (A α,ε ) is irrational and by the properties of Sect. 1.2.2 we know that the corresponding iteration dynamics of A is quasi-periodic and that each individual orbits densely fills T 1 .

Open & Dense Versus Nowhere Dense
In general the (α, ε)-plane of parameters contains a catalogue of the circle dynamics. Again fixing ε = ε 0 > 0 small, consider the corresponding horizontal line in the (α, ε)-plane of parameters. We witness the following, also see Fig. 5 and compare with [24] and references therein. The periodic case corresponds to an open & dense subset of the line, and the quasi-periodic case to a nowhere dense subset, which in the 1-dimensional situation is a Cantor set.

Diophantic Rotation Numbers
Quasi-periodicity corresponds to = (P α,ε 0 ) / ∈ Q. If we restrict even further to Diophantine by requiring that for constants τ > 2 and γ > 0, for all rationals p/q the conjugations of P α,ε 0 with the rigid rotation ϕ → ϕ + 2π can be taken smooth [1,19]. The rotation numbers satisfying (11) form a Cantor subset of the former, which has positive Lebesgue measure, which, by choosing γ = γ (ε 0 ) = O(ε 0 ), can be shown to tend to full measure as ε 0 → 0. A fortiori this holds for the original Cantor set given by Fractal Geometry The Cantor sets under consideration, since they have positive Lebesgue measure, have Hausdorff dimension equal to 1. Moreover Cantor sets have topological dimension 0, since they are totally disconnected: every point has arbitrarily small neighbourhoods with empty boundary. The fact that the Hausdorff dimension strictly exceeds the topological dimension is a characterisation of fractals, see page 15 of [65]. So our Cantor sets are fractals. They also show a lot of self-similarity, a property shared with many other fractals.
Beyond the Arnold Family (10) . . . The organisation of the parameter space in an open & dense set on the one hand, versus a nowhere dense, fractal set of positive measure turns out to be universal, also see [24]. To begin with, any arbitrary smooth (Poincaré) circle mapping of the more general format (8) turns out to have an array of resonance tongues similar to the Arnold family (10), forming an open & dense set that corresponds to periodicity, with a fractal complement which is nowhere dense and of positive measure that corresponds to quasi-periodicity. The only point of difference with (10) is formed by the exact 'sharpness' of the tongues, which depends on the Fourier coefficients of the function f . In particular a tongue at the tip (α, ε) = (p/q, 0) has transverse boundaries if and only if the qth Fourier coefficient does not vanish.

Link with Resonance
Returning to the driven oscillator or the two coupled oscillators we now link periodicity of the Poincaré mapping (8) with resonance. For simplicity we keep |ε| sufficiently small to ensure this mapping to be a diffeomorphism.
As observed in Sect. 1.2.2 the fact that (α, ε) belongs to the p/q-tongue, i.e., that (P α,ε ) = p/q, means that the motion takes place on a p : q torus knot. Generically these periodic orbits come in attracting and repelling pairs and the visible motion takes place on such a periodic attractor. In view of our general 'definition' of resonance in that case we say that the oscillators are in p : q resonance, one also speaks of phase-locking or synchronisation. In the case of 1 : 1 sometimes the term entrainment is being used.
If (α, ε) is outside the tongues, by the Denjoy theory mentioned before, the torus motion takes place on a dense orbit. We also call this torus motion quasi-periodic. In this case from KAM Theory [5,19,44] we derive the following. For Diophantine = (P α,ε ), up to a smooth transformation the Poincaré mapping reads ϕ → ϕ + 2π and the two oscillators correspondingly have the familiar formaṫ such that (P α,ε ) = ω 2 /ω 1 for all ε with |ε| 1. In the parameter space this corresponds to a piece of curve through (α, ε) = ( , 0), parametrized by ε. We emphasize that this uncountable union of these curves has positive measure.

Conclusions and Examples
The literature on resonance phenomena is immense, apart from the references already given, for instance see [3,19,58,85,90] and their bibliographies. For even more references see below. The present point of view models resonant systems in terms of dynamical systems depending on parameters, where resonance takes place in a persistent way.
In the parameter space the resonant set is part of the bifurcation set, which forms a catalogue for transitions to various types of dynamics. What we add to the general discussion on this subject is the overall fractal geometry that usually manifests itself in the complement of all the resonances. We now present a couple of examples.
Huygens's Clocks Returning to Huygens's synchronizing clocks we first consider the problem from the torus flow point of view. We have two almost identical oscillators that are weakly coupled. This means that in (5) for the frequencies we have ω 1 ≈ ω 2 and that |ε| is small. For the Poincaré mapping (8) we only make the assumption that the first Fourier coefficient of f = f (ϕ) does not vanish, which ensures that the 1/1-tongue boundaries meet transverse at (α, ε) = (0, 0), where the first 0 has to be taken mod Z. Compare with Fig. 4.
This implies that (α, ε) belongs to the 1/1-tongue, i.e., that the pendulum clocks are in 1 : 1 resonance, a situation described before as entrainment which is a form of synchronization. This gives a partial explanation of the phenomena discovered by Huygens [61].

Remarks
− Note that the 1 : 1 resonance of the two clocks could be obtained under quite weak assumptions. If one also wants to understand the phase and anti-phase motions, the coupling between the clocks has to be included into the dynamics, compare with [9,76] and references therein. − For another application of these ideas in terms of circadian rhythms and the response to stimuli see [8].
Here it turns out that next to the 1 : 1 'entrainment' resonance also certain other resonances have biological significance. − The above ideas can be largely extended to the case of more than two oscillators. For examples in models for the visual neurocortex see [31,32].

Resonances in the Solar System
From ancient times on resonances have been known to occur in solar system, which are more or less in the spirit of the present section. A wellknown example is the orbital 1 : 2 : 4 resonance of Jupiter's moons Ganymede, Europa, and Io which was studied by De Sitter [51,52] using the 'méthodes nouvelles' of Poincaré [77]. The 2 : 5 orbital resonance between Jupiter and Saturn is described by Moser et al. [67,68,81]. These and other resonances by certain authors are being held responsable for gaps in the rings of Saturn and in the asteroid belt.
Another type of resonance is the spin-orbit resonance. As an example thereof, the Moon is 'captured by' the Earth in a 1 : 1 resonance: the lunar day with respect to the Earth is (approximately) equal to one month. Similarly Pluto and Charon have caught each other in such a 1 : 1 resonance: as an approximately rigid body the two orbit around the Sun. Interestingly, the planet Mercury is captured in a 3 : 2 spin-orbit resonance around the Sun [49].

Remarks
− The spin orbit resonances are explained by tidal forces, for instance, the rotation of the Moon has been slowed down to a standstill by tidal friction brought about largely by the reciprocal tidal forces exerted of Earth and Moon. Similarly the rotation of the Earth in the very long run will be put to a stand still by the tidal forces of mainly the Moon. But probably by that time the Sun has already turned into a red giant . . . − This brings us to the subject of adiabatically changing systems as described and summarized by Arnold [2,3] and which may be used to model such slow changes. 2 One may perhaps expect that the 3 : 2 spin orbit resonance of Mercury in the very long run, and after quite a number of transitions, will evolve towards another 1 : 1 resonance. This part of nonlinear dynamical systems still is largely unexplored.

Periodically Driven Oscillators Revisited
We now return to periodically driven oscillators, showing that under certain circumstances exactly the set-up of Sect. 1.2 applies.
As a motivating example we discuss the Botafumeiro in the cathedral of Santiago de Compostela, see Fig. 6. Here a large incense container is suspended by a pully in the dome where it can swing in the longitudinal direction of the church. A few men pull up the container when it approaches the ground and let go after, thereby creating a periodic forcing and in this way creating a stable motion of exactly twice the period of the forcing.

Parametric Resonance
As a another model consider the parametrically driven oscillator with q(t + 2π) ≡ q(t), see [21]. Here a and ε are considered as parameters. 3 For the periodic function f we have studied several examples, namely f (t) = cos t and = cos t + 3 2 cos(2t) and = signum (cos t), corresponding to the Mathieu case, the Mathieu case modified by a higher harmonic, and the square case. In this setting the issue is whether the trivial 2π -periodic solution is elliptic, hence stable, or not. For any of such systems (12) from the points tongues emanate in the (a, ε)-plane, where inside the trivial 2π -periodic solution is hyperbolic, hence unstable. On the tongue boundaries this solution is parabolic. See Fig. 7 and compare with [14,17]. We note that this gives rise to a discrete union of tongues, where again the sharpness is governed by the Fourier coefficients of the periodic function f = f (t). Fig. 7 Stability diagrams of Mathieu, modified Mathieu and square case [17] Subharmonics and Covering Spaces On the tongue boundaries subharmonic bifurcations occur, see [21,22,44] where each bifurcation can be understood in terms of a pitchfork bifurcation on a suitable covering space. Generally, for the p : q resonance this is constructed as follows. Writing where we identify C ∼ = R 2 , the covering map is given by the Van der Pol transformation On the covering space the group D of decktransformations is generated by which means that • T = for any T ∈ D. 4 Note that D is cyclic of order q, i.e., isomorphic to Z q , for details see [21,22]. Instead of the original system in (z, t) we pull the Fig. 8 Orbits of the Poincaré mapping of the swing (12) system back along , so obtaining a Z q -equivariant system on the (ζ, t)-covering space.
In fact this defines a 1 : 1 correspondence of systems and usually it is most convenient to work on the covering space. Generally on the covering space equivariant Singularity Theory can be practised, see [26,56,57] and references therein, as well as equivariant KAM Theory [23,27,33,37,44].
In the present case this construction is only needed for p/q = k/2. The strongest of these resonances occurs inside the tongue labeled by k = 1, where a cylinder with Z 2 -symmetry is the double cover of a Möbius strip, for a remark in this direction see [58]. This is the setting for the simplest subharmonic bifurcation, namely a period doubling bifurcation: inside the tongue the trivial 2π -periodic solution is unstable but a stable periodic orbit occurs of period 4π . At the boundary a period-doubling bifurcation takes place. The corresponding Z 2 -equivariant bifurcation on the covering space exactly is the pitchfork. The ensuing period-doubled periodic motion is exactly the one that occurs in the Botafumeiro example. See Fig. 8 for 'phase portaits' of the Poincaré mapping, left for ε = 0.25 and right for ε = 0.40

Remarks
− The geometric complexity of the individual tongues in Fig. 7 can be described by Singularity Theory; in fact it turns out that we are dealing with type A 2k−1 , see [17,21]. − The parametric 1 : 2 resonance sometimes also is called the parametric roll. By this mechanism ships have been known to capsize . . . − In Fig. 8 also invariant circles can be witnessed. KAM Theory, as discussed before, in particular an application of Moser's Twist Theorem [66], shows that the union of such invariant circles carrying quasi-periodic dynamics has positive measure. − In both cases the cloud of points 5 is formed by just one or two orbits under the iteration of the Poincaré mapping. These clouds are associated to homoclinic orbits related to the upside down unstable periodic solution, which gives rise to horseshoes. Therefore such an orbit is chaotic since it has positive topological entropy, see [19] and references therein. A classical conjecture is that the cloud densely fills a subset of the plane of positive Lebesgue measure on which the Poincaré mapping is ergodic [4]. 6 Fig. 9 Devil's staircase in the Schrödinger equation with quasi-periodic potential: ω 1 = 1 and ω 2 = 1 2 ( √ 5 − 1), see [28] 2

.2 The Hill-Schrödinger Equation
Another famous equation is a linearized version of (12) where the forcing term is quasiperiodic in t :ẍ where now f (t) = F (ω 1 t, ω 2 t, . . . , ω n t) for a function F : T n → R, see [28,54,69]. As in the case of Sect. 1.2 the countable union of tongues again becomes open & dense and separated by a nowhere dense set of positive measure, determined by Diophantine conditions. The geometry of the individual tongues for small |ε| is exactly as in the periodic case. For larger values of |ε| the situation is more complicated also involving non-reducible quasi-periodic tori, compare with [16]. Equation (14) happens to be the eigenvalue equation of the 1-dimensional Schrödinger operator with quasi-periodic potential. We here sketch how our geometric approach fits within the corresponding operator theory. This operator reads with potential εf ; it acts on wave functions x = x(t) ∈ L 2 (R).
We like to note that in the corresponding literature usually the value of ε = ε 0 = 0 is fixed and the intersection of the horizontal line ε = ε 0 with a tongue is referred to as gap: it is a gap in the spectrum of the Schrödinger operator (15). The approach with tongues and the results of [28] regarding the A 2k−1 -singularity therefore leads to a generic gap closing theory.

Remarks
− In the context of Schrödinger operators the letters are chosen somewhat differently. In particular, instead of x(t) one often considers u(x), which gives this theory a spatial interpretation. Also instead of εf (t) one uses V (x), compare with [69]. − For a fixed value ε = ε 0 the Diophantine Cantor set leads to Cantor spectrum. The total picture is illustrated in the devil's staircase of Fig. 9, where we took n = 2, ω 1 = 1 and ω 2 = 1 2 ( . The rotation number is defined almost as before [28] as a function of a. Fig. 10 Phase portrait of the free Van der Pol oscillator [19] − The nonlinear equationẍ with q quasi-periodic is dealt with in [27]. In comparison with the case of periodic f the averaged, approximating situation, is identical. However, the infinite number of resonances and the Cantorization we saw before leads to an infinite regress of the bifurcation scenarios. For this use was made of equivariant Hamiltonian KAM Theory on a suitable covering space [23,44]. As a consequence the resonant set becomes residual and the quasi-periodic set meagre. Compare this with [7,[45][46][47] in the dissipative case.

Driven and Coupled Van der Pol-Like Oscillators
The examples of the driven oscillator in Sect. 1.1 were based on approximations of the damped pendulum, the free oscillation of which always tends to the lower equilibrium x = 0,ẋ = 0. Our present interest is formed by Van der Pol-like oscillators that for |x| and |ẋ| sufficiently small have negative damping, for this approach compare with [19,44]. Such oscillators are known to occur in electronics [80,88,89]. Therefore our starting point is the periodically driven Van der Pol oscillator in a slightly more general formẍ where the function q is 2π -periodic in the time t . Van der Pol originally considered a(x,ẋ) = bx 2ẋ and f (x,ẋ, t; ε) = sin t . We here assume that a and q are sufficiently smooth, say of class C ∞ or real-analytic. The important restriction on a is that the corresponding free oscillator has a hyperbolic attractor. For a phase portrait in the (x,ẋ)-plane, see Fig. 10. Passing to the system forṁ as before, we get a 3-dimensional state space R 2 × T 1 = {(x, y), z}, so with angular variable z. Let us denote the corresponding vector field by X α,ε .
This brings us back to the general setting of a 2-torus flow, with two phase angles ϕ 1 , ϕ 2 , e.g., with ϕ 1 the phase of the free oscillator, i.e., its time parametrization scaled to period 2π , and ϕ 2 = z. For ε = 0 we so obtainφ which is of the familiar format (5). From here the theory of Sect. 1.2 applies in al its complexity, with in the parameter space an open & dense, countable union of resonances and a fractal set of positive measure regarding quasi-periodicity.
Similar results hold for n coupled Van der Pol type oscillators, now with state space T n , the Cartesian product of n copies of T 1 . Next to periodic and quasi-periodic motion, now also chaotic motions occur, see [19] and references therein.

Universal Studies
Instead of studying classes of driven or coupled oscillators we now turn to a few universal cases of 'generic' bifurcations. The first of these is the Hopf-Neȋmark-Sacker bifurcation for diffeomorphisms, which has occurrence codimension 1. This means that the bifurcation occurs persistently in generic 1-parameter families. However, the open & dense occurrence of countably many resonances and the complementary fractal geometry of positive measure in the bifurcation set are only persistent in generic 2-parameter families. A second bifurcation we study is the Hopf saddle-node bifurcation for diffeomorphisms where we use 3 parameters for describing the persistent complexity of the bifurcation set.

The Hopf-Neȋmark-Sacker Bifurcation
We start with the Hopf-Neȋmark-Sacker bifurcation for diffeomorphisms, but also discuss certain consequences for systems of differential equations. As an example to illustrate our ideas consider the following Duffing-Van der Pol-Liénard type driven oscillator the coefficients of which can be considered as parameters. Note that for the free oscillator at ν 1 = 0 the eigenvalues of the linear part at (x,ẋ) = (0, 0) cross the imaginary axis at ±i √ ν 2 . Excluding the strong resonances where √ ν 2 = p/q for p and q relatively prime with q ≤ 5 and assuming that the constants ν 3 and ν 4 are generically chosen, the organization of the (ν 1 , ν 2 )-plane, for ν 1 consists of an open & dense union of countably many resonance tongues separated by a nowhere dense set of positive measure. This situation is comparable to the Arnold family (10) and Fig. 4, at least for large values of q.

The Nondegerate Case
The general set-up just considers a mapping around a fixed point, say P (0) = 0, where the eigenvalues of the linear part read e ν 1 ±2π i √ ν 2 with |ν 1 | small. We consider a fixed resonance (ν 1 , ν 2 ) = (0, p 2 /q 2 ) with p and q relatively Fig. 11 Tongue in nondegenerate resonance q = 7 of the Hopf-Neȋmark-Sacker bifurcation [26,34,38,43] prime, where we will need 2 parameters to versally unfold the linear part [3,44]. Our main interest is with the periodic points of period q, so in solving the equation The zeroes of P q − Id are studied by a Lyapunov-Schmidt reduction, which leads to a Z qequivariant family of functions where z is an appropriate complex variable and where B μ is a polynomial in u = |z| 2 of degree less than (q − 1)/2. We study the corresponding discriminant set given by G μ (z) = 0 and det D z G μ (z) = 0.
Here μ is an unfolding-multiparameter detuning the resonance at hand. The way to study this discriminant set is by Z q -equivariant contact equivalence [26,34,38,40,43]. In the present non-degenerate case (22) can be reduced to the polynomial normal form for a complex parameter σ . In general this set turns out to be a 'tongue' ending in a cusp of sharpness (q − 2)/2, which is part of a familiar bifurcation diagram with two periodic orbits of period q inside that annihilate one another at the tongue boundaries in a saddle-node or fold bifurcation [3]. See Fig. 11 which is embedded in the context of (18), of which P is a Poincaré mapping. Here the dynamics of P also has been described in terms of a Poincaré-Takens interpolating normal form approximations, e.g., see [13,21,22,86].
Globally a countable union of such cusps is separated by a nowhere dense set of positive measure, corresponding to invariant circles with Diophantine rotation number. As before, see Fig. 4, the latter set contains the fractal geometry.

Fig. 12
Two-dimensional tomography in a mildly degenerate resonance q = 7 of the Hopf-Neȋmark-Sacker bifurcation [26,34,38,43] − The strong resonances with q = 1, 2, 3 and 4 form a completely different story where the Singularity Theory is far more involved [3,86]. Still, since the higher order resonances accumulate at the boundaries, there is fractal geometry around, always of positive measure. − Regarding structural stability of unfoldings of P as in (19) under topological conjugation, all hopes had already disappeared since [71].

A Mildly Degenerate Case
In the mildly degenerate case the 'Hopf' coefficient in the previous example vanishes and is introduced as an extra parameter. This leads to another normal form which is now parametrized over C 2 ∼ = R 4 , hence of codimension 4. As before [26,40] the normal form (24) is structurally stable under Z q -equivariant contact equivalence. Here the Singularity Theory is more complex, involving folds, cusps, a swallowtail and a Whitney umbrella. The complete bifurcation diagram is more involved, even at the level of Poincaré-Takens normal form vector field approximations. In Fig. 12, again for the case q = 7, a 2dimensional tomography is shown of the bifurcation set, in which still a tongue-like structure can be discerned, for more details see [38,43] where also 3-dimensional tomographies are presented.
To illustrate a mildly degenerate case of the Hopf-Neȋmark-Sacker bifurcation one may well consider the preceding Duffing-Van der Pol-Liénard type driven oscillator (18) where we need all four parameters.  Table 1 for the color coding

Concluding Remarks
For both cases of the Hopf-Neȋmark-Sacker bifurcation we have a good grip on the part of the bifurcation set that governs the number of periodic points. The full bifurcation set is far more involved and the corresponding dynamics is described only at the level of Poincaré-Takens normal-form vector fields [13,38,43]. We note that homo-and heteroclinic phenomena occur at a flat distance in terms of the bifurcation parameters [15,18,75].

The Hopf Saddle-Node Bifurcation for Diffeomorphisms
As a continuation of the above programme, we now consider the Hopf saddle-node (or fold Hopf) bifurcation for diffeomorphisms [35,36,41], in which the central singularity is a fixed point of a 3-dimensional diffeomorphism, such that the eigenvalues of the linear part at bifurcation are 1 and e 2π iα , where e 2πniα = 1 for n = 1, 2, 3 and 4, so excluding strong resonances as in the Hopf-Neȋmark-Sacker case of Sect. 3.1. The Hopf saddle-node bifurcation for flows is well-known [20,58], especially because of the subordinate Hopf-Neȋmark-Sacker and Šilnikov homoclinic bifurcation. Our main interest is how the Hopf-Neȋmark-Sacker bifurcation is being changed into one of the simplest quasiperiodic bifurcations near a 2 : 5 resonance.  Fig. 13, see [35,41]. The attractors are classified by means of the

From Vector Fields to Mappings
The linear part of the vector field at bifurcation has eigenvalues 0 and ±iα. This linear part generates an axial symmetry that in a normal form procedure can be pushed stepwise over the entire Taylor series, see [13] and references. This makes it possible to first consider axially symmetric systems, that turn out to be topologically determined by their 3rd-order truncation given byẇ where w ∈ C and z ∈ R and where β 1 and β 2 are unfolding parameters [58,62]. A scaling β 1 = γ 2 , β 2 = γ 2 μ leads to a vector field From this an axially symmetric map is obtained in a kind of Euler step. To study a 2 : 5 resonance we take α 0 = 2/5, writing α = α 0 + γ δ, and perturb to by adding axially non-symmetric order r resonant terms w 4 ∂ ∂w and Re w 5 ∂ ∂z .
The scaling parameter γ and the other constants are fixed suitably.

In the Product of State Space and Parameter Space
In Fig. 13 a Lyapunov diagram is depicted in the parameter plane of the mapping family G. Table 1 contains the corresponding color code. In Fig. 14 we show the dynamics corresponding to two values of (μ, δ/2π). Let us discuss these numerical data.
The Parameter Space On the right-hand-side of the figure this method detects an attracting invariant circle of focus type (blue). In the gaps larger resonances are visible, compare with Fig. 4 for a fixed value of ε. Moving to the left, in the neighbourhood of the line indicated by H a quasi-periodic Hopf bifurcation occurs from a circle attractor to a 2-torus attractor (green). Also here the parameter space is interspersed with a resonance web of which the larger lines are visible. The remaining features, among other things, indicate invariant tori and strange attractors of various types and also more invariant circles.
The State Space The upper two figures of Fig. 14 show an invariant circle, once seen from the z-direction and once from some w-direction. The lower two figures indicate how this circle has become a strange attractor, from the same two points of view. Part of these results can be justified mathematically, as seen from the Perturbation Theory point of view. The invariant circles all have one Lyapunov exponent equal to 0 and these are quasi-periodic, perturbations of closed integral curves of a vector field (averaging) approximation, whence their existence can be proven by KAM Theory [5,23]. A similar statement can be made about the 2-tori with two Lyapunov exponents equal to 0. In fact the transition is a quasi-periodic Hopf bifurcation as discussed by Broer et al. [10,11,19,44].
By the same references, this also holds for the quasi-periodic invariant circle in the upper half of Fig. 14. The lower half of this figure is conjectured to show a quasi-periodic Hénon-like attractor, which is the closure of the unstable manifold of an unstable quasiperiodic invariant circle. This is the previous quasi-periodic circle that has become unstable through a quasi-periodic saddle-node bifurcation [11]. For this kind of strange attractor the mathematical background theory largely fails, so the results must remain experimental; for indications in this direction however see [41] and references therein.
For a detailed, computer-assisted bifurcation analysis of the 2 : 5 resonance 'bubble' we refer to [35]. Compare with earlier work of Chenciner [45][46][47]. We like to note that the family of mappings G forms a concrete model for the Ruelle-Takens scenario regarding the onset of turbulence. In fact it also illustrates how the earlier scenario of Hopf-Landau-Lifschitz is also included: the present multi-parameter set-up unifies both approaches. For details and background see [19,44,59,60,63,64,70,78,79].
Resonance and Fractal Geometry Interestingly, the blue colors right and left correspond to quasi-periodic circle attractors. The fact that the corresponding regions of the plane look like open sets is misleading. In reality these are meagre sets, dense veined by the residual sets associated to periodicity. These details are just too fine to be detected by the computational precision used.
Particularly in the latter case, in the left half of the diagram, we are dealing with the Arnold resonance web, for a detailed analysis see [36].

Conclusions
We discuss a number of consequences of the present paper in terms of modelling of increasing complexity.

'Next Cases'
The Hopf saddle-node bifurcation for maps, see Sect. 3.2, can be viewed as a 'next case' in the systematic study of bifurcations as compared to, e.g., [58,62] and many others. The nowhere dense part of parameter space, since it lacks interior points is somewhat problematic to penetrate by numerical continuation methods. Nevertheless, from the 'physical' point of view, this part surely is visible when its measure is positive or, as in the present examples, even close to full measure. Needless to say that this observation already holds for the Hopf-Neȋmark-Sacker bifurcation as described in Sect. 3.1.
Other 'next cases' are formed by the quasi-periodic bifurcations which is a joint application of Kolmogorov-Arnold-Moser Theory [5,11,19,30,44] and Singularity Theory [55][56][57]87]. For overviews see [12,48,91]. The quasi-periodic bifurcations are inspired by the classical ones in which equilibria or periodic orbits are replaced by quasi-periodic tori. As an example, in the Hopf saddle node of Sect. 3.2 we met quasi-periodic Hopf bifurcation for mappings from circles to a 2-tori in a subordinate way. Here we witness a global geometry inspired by the classical Hopf bifurcation, which concerns the quasi-periodic dynamics associated to the fractal geometry in the parameter space, compare with Fig. 13. The gaps or tongues in between concern the resonances inside, within which we notice a further 'fractalization' or 'Cantorization'.
A similar 'next case' in complexity is given by the parametrically forced Lagrange top [29,30], in which a quasi-periodic Hamiltonian Hopf bifurcation occurs. Indeed, we recall from [50] that in the Lagrange top a Hamiltonian Hopf bifurcation occurs, the geometry of which involves a swallowtail catastrophy. By the periodic forcing this geometry is 'Cantorized' yielding countably many tongues with fractal geometry in between.

Remarks
− As said before, in cases with infinite regress the fractal complement is a meagre set which has positive measure. Simon [82] describes a similar situation for 1-dimensional Schrödinger operators. Also see [6]. − It is an interesting property of the real numbers to allow for this kind of dichotomy in measure and topology, compare with Oxtoby [73]. Interestingly, although these properties in the first half of the 20th century were investigated for theoretical reasons, they here naturally show up in the context of resonances and spectra.

Modelling
We like to note that our investigations on the Hopf saddle-node bifurcation for mappings were inspired by climate models [25,42,83], where in about 80-dimensional Galerkin projections of PDE models such bifurcations were detected in 3-dimensional center manifolds. Generally speaking there exists a large-scale programme of modelling in terms of dynamical systems depending on parameters, with applications varying from climate research to mathematical physics and biological cell systems. These models often are high-dimensional and their complexity is partly explained by mechanisms of the present paper, also see [19,78,79,91]. In general such models exhibit the coexistence of periodicity (including resonance), quasi-periodicity and chaos, best observed in the product of state-and parameter space.