#### Abstract

A phase transition from integrability to nonintegrability in two-dimensional Hamiltonian mappings is described and characterized in terms of scaling arguments. The mappings considered produce a mixed structure in the phase space in the sense that, depending on the combination of the control parameters and initial conditions, KAM islands which are surrounded by chaotic seas that are limited by invariant tori are observed. Some dynamical properties for the largest component of the chaotic sea are obtained and described in terms of the control parameters. The average value and the deviation of the average value for chaotic components of a dynamical variable are described in terms of scaling laws, therefore critical exponents characterizing a scaling function that describes a phase transition are obtained and then classes of universality are characterized. The three models considered are: The Fermi-Ulam accelerator model, a periodically corrugate waveguide, and variant of the standard nontwist map.

#### 1. Introduction

Dynamical systems described by two-dimensional nonlinear area preserving mappings have been considered for many years as prototype for the study of low-dimensional physical systems. Applications of the description of dynamical systems by mappings are observed in the study of magnetic field lines in toroidal plasma devices with reversed shear (like tokamaks), channel flows, waveguide, Fermi acceleration and many other [1–10].

In this paper we consider a special class of two-dimensional area preserving mappings exhibiting a phase transition from integrable to nonintegrable. For certain values of the control parameter, the phase space of the system has only periodic and quasiperiodic orbits, thus the phase space is filled by straight lines. As this control parameter varies, the phase space exhibits a mixed form containing both Kolmogorov-Arnold-Moser (KAM) islands surrounded by a chaotic sea (conservative chaos) and a set of invariant spanning curves limiting the size of the chaotic sea. This form of the phase space is generic for nondegenerate Hamiltonian systems [11] and can be observed in many different problems including time dependent potential [12–16], tokamaks [17–20] and also in billiards [21–23].

The main approach considered in this paper is the description of the behavior of some average dynamical variables as a function of the control parameter, like the first momenta and the deviation around the average variable. This characterization is made along the chaotic sea near a transition from integrable to nonintegrable. This formalism furnishes critical exponents that can be used to define classes of universality.

This paper is organized as follows. In Section 2, we present and discuss some details needed to construct a generic two-dimensional area preserving mapping. Then application for a specific mapping is made including a discussion of the stability of fixed points and properties of the phase transition from integrable to nonintegrable. In Section 3 the discussion is devoted to a periodically corrugate waveguide. Some results for the phase space are also discussed in connection with the Standard Mapping. The scaling results for the phase transition are obtained and discussed. In Section 4, we extend the applicability of the theory for the one-dimensional Fermi accelerator model confirming that the rather two distinct models discussed in Sections 3 and 4 belong to the same class of universality and experience the same transition . Finally, the conclusion and final remarks are drawn in Section 5.

#### 2. A Generic Two-Dimensional Area-Preserving Map

In this section we present and discuss some dynamical properties for a two-dimension mapping. We assume that there is a two-dimensional integrable system that is slightly perturbed. Then, the Hamiltonian function that describes the system is written as

where the variables and with correspond, respectively, to the action and angle. Considering the Poincaré section given by the plane assuming that is constant (mod ), one can write a modified twist mapping as

where both , , and are, in principle, nonlinear functions of their variables and the index correspond to the th iteration of the mapping. The variables and corresponds indeed to and .

Since the mapping (2.2) must be area preserving (it was obtained from an Hamiltonian) then the functions and must obey some intrinsic relations. The relations are obtained considering that the determinant of the Jacobian matrix is the unity. Then we have that the coefficients are written as

Given the above coefficients, the determinant of the Jacobian matrix is given by

It is easy then to conclude that area preservation will be observed only if the condition

is matched. For many mappings, the function . Therefore, if we keep the function as , just to illustrate applicability of the formalism, the following mappings have already been studied:

(i), then the Taylor-Chirikov's map is recovered;(ii), then the Fermi-Ulam accelerator model is obtained;(iii), with constant, then the bouncer model is found;(iv)for the case of where is a constant, then we recovered the so-called Hybrid-Fermi-Ulam-bouncer model;(v)considering , the logistic twist mapping is obtained.In this section, we shall consider the following expression for the two-dimensional area-preserving mapping:

where the variables and were replaced, by sake of simplicity, by and . The control parameters of the model are and . The modulus function is requested since the argument of the denominator of the second equation above is not defined for nonpositive values. It is easy to see that when the control parameter , then the system is integrable. For , the phase space of the system exhibits chaotic components and periodic orbits, as we can see in Figure 1. The determinant of the Jacobian matrix of the mapping (2.7) is the unity, so area preservation is observed.

**(a)**

**(b)**

The fixed points of the mapping (2.7) are obtained by matching the following conditions: and with . After solving these two conditions, the two sets of fixed points are given by: fixed point 1 ; fixed point 2 . The fixed point 1 is always hyperbolic while fixed point 2 is elliptic for

From now on, we shall consider as fixed the control parameter . Results for will be presented elsewhere.

We can see in Figure 1 that the size of the chaotic sea varies as the control parameter varies. As rises, the position of the lowest invariant spanning curve also rises. It is then expected that the average value of , say changes as a function of the control parameter. To obtain the average value of , we proceed as follows. First we obtain the average value along an orbit

and then we evolve the average over an ensemble of different initial conditions along the chaotic sea. The average value is then written as

where is a sample of the ensemble. It is shown in Figure 2 the behavior of for three different control parameters. We can see that starts growing and then, after reaching a characteristic crossover iteration number, it bends towards a regime of convergence. We see also that the three curves start growing at different positions. Therefore, a transformation coalesces all curves at small iterations, as we can see in Figure 2(b).

**(a)**

**(b)**

The behavior shown in Figure 2 allow us to suppose the following.

(i)For low iterations and after a short initial transient, that is, , the average is described as , where is a critical exponent. (ii)For large enough , say , then , where is a critical exponent. (iii)The number of iterations which characterizes the changeover from growth to the saturation is given where is also a critical exponent.Plotting the behavior of and we can obtain the critical exponents and . The corresponding plots of and are shown in Figure 3. After applying a power law fitting we obtain that and . The critical exponent was obtained for several different time series. Its average value is around .

**(a)**

**(b)**

Rescaling properly, we can see that all curves coalesce together, as they are shown in Figure 4. Such collapse is a clear evidence of the criticality for the chaotic sea near a transition from integrability to nonintegrability. It must be emphasized that the values obtained for the critical exponents are useful to define classes of universality. As we will see in the two next sections, for the periodically corrugate waveguide and for the one-dimensional Fermi accelerator model, the exponents obtained here are different from the critical exponents obtained there, thus confirming that, despite the similarity of the two-dimensional mappings, the classes of universality of such models are different.

#### 3. Periodically Corrugate Waveguide

Let us now discuss in this section some properties of a periodically corrugated waveguide. This section is a review of papers [24, 25]. The interest in describing the problem of * guiding a light ray* inside a periodically corrugated boundary has increased along last years, mainly because the topic is applicable in so many different fields of science. Applications involve, for example, ray chaos in underwater acoustics [26–28], quantum transport in ballistic cavities [29, 30], transport through a finite hetero-structure [31], quantized ballistic conductance in a periodically modulated quantum channel [32], comparisons of classical versus quantum behavior in periodic mesoscopic systems [33, 34], scattering of a quantum particle in a rippled waveguide [35] and anomalous wave transmittance in the stop band of a corrugated parallel-plane waveguide [36].

There are many different ways to describe problems involving waveguides. One of them consists in considering the well-known billiard approach (A billiard problem consists of a system in which a point-like particle moves freely inside a bounded region and suffers specular reflections with the boundaries.). Generally, the dynamics is described using the formalism of discrete maps. Therefore, depending on the combinations of the control parameters as well as on the initial conditions, the phase space for such mappings might be included in three different classes namely: (i) regular, (ii) ergodic, and (iii) mixed. Roughly speaking, the integrability of the regular cases is generally related to the angular momentum preservation and the static circular billiard is a typical example. On the other hand, for the completely ergodic billiards, only chaotic and therefore unstable periodic orbits are present in the dynamics. The well-known Bunimovich stadium [37] and the Sinai billiard [38] (sometimes also-called as Lorenz gas) are typical examples of case (ii). On these two kinds of systems, the time evolution of a single initial condition, for the appropriated combinations of control parameters, is enough to fill up ergodically the entire phase space. Finally, there are many different billiards that present mixed phase space structure [39–43], which have control parameters with different physical significance. Depending on the combination of both initial conditions and control parameters, the phase space presents a very rich structure which contains invariant spanning curves (sometimes also-called as invariant tori), KAM islands and chaotic seas.

In this section, scaling arguments are used to describe the behavior of the deviation of the average reflection angle within the chaotic sea, for the problem of a * classical* light ray inside a periodically corrugated waveguide suffering specular reflections. The model consists of a classical light ray which is specularly reflected between a flat plane at and a corrugated surface given by , where is the average distance between the flat plane and the corrugated surface, is the amplitude of the corrugation and is the wave number. The dynamical variables, which are used in the description of the problem are the angle of the trajectory measured from the positive horizontal axis and the corresponding value of the coordinate at the instant of the reflection. It must be emphasized that the mapping is iterated when the ray light * hits* the flat plane, thus multiple reflections with the periodically corrugated surface are neglected. (This is a standard approximation which is very useful to speed up the numerical simulations since no transcendental equations must be solve, as they have to be in the corresponding complete version of the problem).

##### 3.1. The Model and the Mapping

Let us discuss now all the details needed to describe the model and obtain the equations that describe the dynamics of the problem. The model thus consists in obtaining a mapping , given the initial conditions as shown in Figure 5(a). The mapping is obtained from purely geometrical considerations of Figure 5(b). It is easy to see that

**(a)**

**(b)**

and in a similar way

We stress that the term gives the exact location of the collision on the corrugated surface. On the other hand, the angle is written as

where gives the slope of the corrugated surface at the point . The condition of specular reflection (A specular reflection is that in which the angle of incidence is equals to the angle of reflection. See details in Figure 5(b)) was used in the derivation of the (3.3). The Equations (3.2) and (3.3) correspond to the exact form of the mapping. However, I will discuss in this paper only a simplified version. To do this, the following approximations are assumed.

(1)The relative corrugation is assumed to be small, so that , yielding consequently . (2)In the limit of small corrugation, it is also assumed that .Let us now argue on the approximation of small relative corrugation. It is easy to see that, for , the system is integrable, since the two boundaries are parallel plates and the phase space shows only straight lines. However, if the phase space birth as mixed therefore containing both chaos, invariant spanning curves and KAM islands. Thus, there is an abrupt transition from integrability to nonintegrability when the control parameter goes from to . This is the transition we shall address in this section. We stress also that the formalism can be extendable to many other systems and even to other different kinds of transitions. Before write the equations of the mapping, we can see that there are an excessive number of control parameters, namely , , . The dynamics however does not depend on all of them. It is therefore convenient to define the following dimensionless and more convenient variables , and finally . Thus the only relevant control parameter is . With these new variables, the two-dimensional nonlinear mapping is given by

This mapping preserves the phase space measure since that the determinant of the Jacobian matrix is the unity.

It is shown in Figure 6, the phase space generated by the iteration of mapping (3.4). We can see that there is a complex structure which contains KAM islands surrounded by a large chaotic sea that is confined between two invariant spanning curves. The chaotic sea was characterized using the well-known Lyapunov exponent. Such exponent shows great applicability as a practical tool that can quantify the average expansion or contraction rate for a small volume of initial conditions. As discussed in [44], the Lyapunov exponents are defined as

where are the eigenvalues of and is the Jacobian matrix evaluated over the orbit . It is important to stress that a direct implementation of a computational algorithm to evaluate (3.5) has a severe limitation to obtain . Even in the limit of very short , the components of the matrix can assume very different orders of magnitude for chaotic orbits and periodic attractors, yielding impracticably the implementation of the algorithm. In order to avoid such problem, we note that can be written as where is an orthogonal matrix and is a right triangular matrix. Thus we rewrite as , where . A product of defines a new . In a next step, it is easy to show that . The same procedure can be used to obtain and so on. Using this procedure, the problem is reduced to evaluate the diagonal elements of . Finally, the Lyapunov exponents are now given by

If at least one of the is positive then the orbit is classified as chaotic. Therefore we define as the larger value of the . It is shown in Figure 7 the behavior of the positive Lyapunov exponent for an ensemble of 10 different initial conditions on the chaotic sea for the control parameter . The initial conditions were chosen such that and different initial values of were uniformly distributed along the interval . Each initial condition was iterated up to times. The average value obtained was where the error corresponds to the standard deviation of the samples.

Let us now discuss one of the consequences of the variation of the control parameter on the location of the invariant spanning curves. Since the two first invariant spanning curves define the limit for the chaotic sea, the variation of consequently change the edges of the chaotic sea. Before do that, I must first give a brief review in the well-known * Standard Map* (see also [45] for specific details). It is defined as

where is a control parameter. This mapping can also shows mixed phase space for . Moreover, one of the most important properties is that it shows a transition from local to global chaos (see [45, 46] for specific details) when the control parameter matches the condition . After this transition, no invariant spanning curves are present and consequently the chaotic sea spreads over the phase space. The connection of this result with the phase space of the mapping that describes the behavior of the light ray in a periodically corrugate waveguide suppose that, near the invariant spanning curves that limit the chaotic sea, the dynamics could be locally described by using the standard mapping. Thus we can suppose that near the invariant spanning curve, the reflection angle can be written as

where is a typical value of the reflection angle along the invariant spanning curve, and is a small perturbation of the angle. After defining , the first equation of the mapping (3.4) is written as

Using (3.8), we can rewrite (3.9) as

Expanding (3.10) in Taylor series, we obtain that

Considering the first representative term in the series, we obtain that

On the other hand, the second equation of mapping (3.4) may be rewritten as

Multiplying both sides of (3.13) by and adding again in both sides the term , defining the term

and calling , we obtain the corresponding mapping

An immediate comparison of (3.7) and (3.15) allow us to conclude that there is an effective control parameter given by

Since the transition from local to global chaos occurs at , we obtain that the two invariant spanning curves that limit the chaotic sea are given by

Therefore we can conclude that the *size* of the chaotic sea is proportional to [25].

##### 3.2. Scaling Approach

This section is devoted to discuss a scaling property which is present in the chaotic sea. As it was discussed in Section 2, an initial condition in the region of the chaotic sea wanders in a chaotic way along the accessible region, however it is always confined between two invariant spanning curves. The location of the first positive and negative invariant spanning curves depend on the value of the control parameter . As a consequence, the “amplitude” of a chaotic time series is dependent on the control parameter too. We will now explore the behavior of the deviation of the average value for the angle , an observable defined as

where it is used as

The (3.19) gives the correspondent average reflection angle averaged over the orbit, while the in (3.18) corresponds to an ensemble of different initial conditions. Therefore, (3.18) was iterated using an ensemble of different initial conditions. The variable was kept fixed as while the values of were uniformly distributed along . It is shown in Figure 8(a) the behavior of three different values of for different control parameters. It is easy to see that grows for short iterations and then suddenly it bends towards a regime of saturation for large enough iterations. The changeover from growth to the saturation regime is marked by a characteristic crossover iteration number denoted as . It is worth stressing that different values for the control parameters generate different curves for for short iteration . This might be an indicative that is not a “good” scaling variable. Thus, a trick transformation coalesces all the curves for short iteration number, as can be seen in Figure 8(b). Based on the behavior of Figure 8, the following three scaling hypotheses are proposed.

**(a)**

**(b)**

(i) For a short iteration number, say , behaves according to the power law where is a critical exponent.

(ii) For large enough iteration number, say , approaches a regime of saturation marked by a constant “plateau’’ given by where the exponent is a critical exponent;

(iii) The number of iterations that characterizes the crossover, that is, the iterations that marks the change from growth to the saturation is written as
The exponent is defined as the * dynamical* exponent.

This scaling formalism is commonly used in the description of phase transitions in critical phenomena. It is also very useful in studies of problems of surface sciences (see, e.g., [47]). Considering these three initial suppositions, can now be formally described in terms of a scaling function of the type

where is a scaling factor, and are the so-called scaling exponents. Moreover, the exponents and must be necessarily related to the critical exponents , and . Since is a scaling factor, it could be chosen that and then (3.23) is rewritten as

where the function is assumed to be constant on the limit of . A comparison of (3.24) and (3.20) allows us to conclude that . Choosing now that , (3.23) is given by

where the function is assumed to be constant for . Comparing now (3.25) and (3.21), it is easy to see that . Given the two different expressions of the scaling factor , we can obtain the relation between the critical exponents as

Note however that the scaling exponents are all determined if the critical exponents and were numerically obtained. Thus, it is shown in Figure 9 the behavior of (a) and (b) . The point which is important to stress is that the saturation values were obtained via extrapolation since that even after almost , the saturation value for has not yet been reached. A power law on Figures 9(a) and 9(b) gives that , and . Since the dynamical exponent can also be obtained of (3.26), from evaluation of the numerical values of and , it was found that . Such result is indeed in good agreement with the numerical result obtained in Figure 9(b).

**(a)**

**(b)**

Given that the values of the critical exponents are now obtained, the scaling hypotheses can be verified. In this sense, it is shown in Figure 10 the collapse for three different curves generated from different values of the control parameters onto a single and universal plot. The control parameters used in Figure 10(a) were , and . The collapse obtained in Figure 10(b) reinforces that the scaling suppositions are indeed correct and therefore confirms a criticality near .

**(a)**

**(b)**

#### 4. One-Dimensional Fermi Accelerator Model

In this section we will use the formalism described in Section 3 to characterize the same transition in a rather different model, the well-known one-dimensional Fermi accelerator model (FAM). This section is a short review of [48]. The physical motivation of this problem is far away different of the corrugated waveguide, however and as we will see, the two different model belong to the same class of universality. The origin of the FAM back to earlier 1949 when Enrico Fermi [49] proposed a model in order to study the cosmic ray acceleration. After that, his original model was modified and studied using many different approaches [50–54]. A version which we address now, as it is shown in Figure 11, consists in a classical particle which is confined between two rigid walls and suffers elastic collisions with them. One of the walls is fixed at and the other one is periodically time varying . The dynamics is described in terms of a two-dimensional nonlinear mapping for the variables and , where the index corresponds to the *n*th collision with the moving wall. Before obtain the equations of the mapping, let us argue on a simplification commonly used. We will suppose that both walls are fixed but that, when the particle suffers a collision with one of the walls, it exchanges momentum as if the wall was moving. This simplification carries the huge advantage of allowing us substantially to speed up the numerical simulations compared with the complete model (In the complete model, the movement of the oscillating wall is taken into account. Therefore, the instant of the collisions is obtained via a solution of a transcendental equation.) It is usefully applicable because the main dynamical properties of the system are preserved under such conditions. Incorporating this simplification in the model and using dimensionless variables , and measuring the time in term of , the mapping is written as

The term specifies the length of time during which the particle travels between collisions, while gives the corresponding fraction of velocity gained or lost in the collision. The modulus function is introduced to avoid the particle leaving the region between the walls. We stress that the approximation of using the simplified FAM is valid in the limit of small . In such a limit, the same transition observed for the periodically corrugated waveguide is observed in the FAM [48].

The phase space generated for iteration of the mapping (4.1) is shown in Figure 12. It is easy to see that there is a large chaotic sea surrounding KAM islands. The chaotic sea is also limited by an invariant spanning curve. This model exhibits also scaling properties near the phase transition from integrable to nonintegrable like those of the periodically corrugate waveguide. However, the observable is now obtained as

where the average velocity is defined as

and is an ensemble of different initial conditions.

It is shown in Figure 13(a), the behavior of for two different control parameters, while Figure 13(b) shows the curves of after a suitable transformation .

**(a)**

**(b)**

The exponent is obtained in the asymptotic limit of large iteration number and it is independent of . Figure 14(a) illustrates an attempt to characterize this exponent using the extrapolated saturation of . Extrapolation is also required because, even after iterations, has still not quite reached saturation. From a power law fit, we obtain . After averaging over different values of the control parameter in the range , we then obtain . It is found that the scaling relation for the exponent is . Considering the previous values of both and , it is found . The exponent can be also obtained numerically. Figure 14(b) shows the behavior of the crossover iteration number as function of the control parameter . The power law fit gives us that , in good accord with the scaling result. The scaling for is demonstrated in Figure 15, where the three different curves for in Figure 15(a) are very well collapsed onto the universal curve seen in Figure 15(b).

**(a)**

**(b)**

**(a)**

**(b)**

Since the critical exponents are the same of those obtained for the corrugated waveguide, we can conclude that the one-dimensional Fermi accelerator model belongs to the same class of universality of the corrugated waveguide and experiences the same transition from integrability to nonintegrability when the control parameter goes from to . However the critical exponents obtained for a variant of a standard nontwist map are different from those obtained in Sections 3 and 4. Despite the similarities of the model presented in Section 2 and those models discussed in Sections 3 and 4, it belongs to a different class of universality from the corrugate waveguide and Fermi accelerator model since the critical exponents are rather different.

#### 5. Conclusion

In summary, we have studied a phase transition from integrability to nonintegrability in three different two-dimensional mappings. Critical exponents for the behavior of the chaotic time series were obtained and used to define classes of universality. Even though the mappings are similar in their forms, the models discussed in Sections 3 and 4 belong to the same class of universality but then, they do not belong to the same class of universality of a variant of a standard nontwist map since the critical exponents are rather different. Indeed the term that seem to define the universality class is the expression of (see (2.2)) and the exponent of the variable in that equation.

#### Acknowledgments

Financial support from CNPq, FAPESP and FUNDUNESP, Brazilian agencies are gratefully acknowledged. The present paper contains totally new results, those shown in Section 2 and also bring a short review of the articles [24], [25], and [48].