In this subsection I demonstrate the influence of the kick amplitude for the case ; similar results can be obtained for and . The following discussion holds for small values of . The idea is to derive a separatrix mapping that describes the dynamics of the kicked harmonic oscillator in the vicinity of the separatrices and allows to estimate the width of the channels of diffusive dynamics.
The first step is to identify
the
of equation (1.39c) as a
high frequency perturbation to the time averaged Hamiltonian
(1.44).
At kick times
,
, one has
(1.25) |
From the last formula it becomes clear that can indeed be treated as a high frequency perturbation of the HARPER Hamiltonian : the time dependence of is given by cosine terms the largest period of which is , whereas the smallest period of is -- as is shown below on page -- and can thus be made as large as desired in the limit .
In [LL73] it is argued
that in a first approximation
the higher frequency terms of
a perturbation expanded as in equation (1.50)
can be neglected, although
they contribute with
roughly
the same
weight
() as the small frequency terms.
Therefore I can drop all cosine terms with and get the
approximant
Now consider a typical orbit of the web map in the stochastic region of the web. For some the corresponding orbit point will be close to and just below of the midpoint of the HARPER separatrix with , which is displayed in the upper right quadrant of figure 1.13.1.7
For notational convenience I define with some . At , takes on a certain value since the separatrices are characterized by . Let be the initial value (at time ) for the successive application of the fourth iterate of the iteration of the web map (1.29), . The resulting orbit in the -plane is quasiperiodic if is not too near to the separatrix (cf. figure 1.7f). First the orbit follows the separatrix anticlockwise, until it comes close enough to the hyperbolic point , where the orbit turns left and follows the perpendicular separatrix with ; this is shown schematically in figure 1.14. At some time (where is the according number of iterations of ) the orbit reaches the point , which is defined as that iterate of under that comes closest to the centre of the second separatrix. This point again gives rise to a certain value of , . The desired separatrix mapping is now given by the change of the value of during one such quarter revolution and by the corresponding (approximate) quarter period:
In order to determine an approximate expression for
, I first derive an explicit solution of the
HARPER dynamics
on the separatrix for . Such a solution is given by
The rate of change of the value of
during the
-perturbed dynamics can then approximately be
calculated
as
(1.27) |
(1.28) |
(1.30) |
It remains to determine the time interval in equation
(1.52b).
Within the cells, i.e. away from the separatrices, the equations of
motion of the unperturbed HARPER system
can be integrated using
Jacobian
elliptic functions.
With
(1.32) |
(1.33) |
(1.34) |
(1.35) |
For real arguments,
as in the present case,
is periodic with period , being a
complete elliptic integral of the first kind:
Approximate expressions for the integral (1.65) can be found
in [AS72] and yield here, in the vicinity of a
separatrix:
With equations (1.59) and (1.68),
the separatrix mapping (1.52), approximating the dynamics of
near the separatrices, is completely specified.
Note that
the
separatrix mapping can also be obtained in a different way:
up to terms of order , the fourth iterate of the web map
(1.29) is given by
(1.41) |
Using the separatrix mapping, I now can proceed
to the estimation of the width of the channels of diffusive dynamics.
As a criterion for the
border between regular dynamics within the meshes
and stochastic dynamics in the channels,
(1.44) |
Figure 1.15 illustrates this behaviour: for several
values of stochastic webs are obtained numerically by iterating the
web map for a large number of times;
then the channel widths of these webs are measured
-- by determining the largest distance of any two points of the web
which approximately lie on the line and near the point
-- and compared with the numbers given by the approximating formula
(1.73).
The agreement between the analytical formula
and the numerical data is reasonably good for
and improves for . Figure 1.15 also indicates
that for
the precision of the computer
algorithm
does not suffice any more to produce accurate numerical results,
because in this parameter range the web gets very thin
-- the width of the channels
shrinks to less than
here.
Similar results hold for the other values of
.
In this subsection I have discussed how the kick strength of the kicked harmonic oscillator determines the shape of its phase portrait. In particular it has been shown that the periodic kicking acts in a way which is typical for perturbed systems: starting from an integrable system at (here: the harmonic oscillator) the kicking (with ) renders the system nonintegrable, and the area of the phase space region of irregular, chaotic dynamics grows with increasing perturbation parameter .
Having investigated the shape and the size of the channels of irregular motion in this subsection and the previous one, in the next subsection I turn to the discussion of the most characteristic dynamical aspect of the irregular motion.