Journal of Chaos
Volume 2013 (2013), Article ID 587548, 7 pages
Difficulties in Evaluating Lyapunov Exponents for Lie Governed Dynamics
1Ciclo Basico Comun, Catedra de Fisica, Universidad de Buenos Aires, Argentina and Laboratorio de Sistemas Complejos, FI-UBA, Argentina
2Instituto de Física La Plata-CCT-CONICET, Universidad Nacional de La Plata, C.C. 727, 1900 La Plata, Argentina
3Physics Department and IFISC-CSIC, University of the Balearic Islands, 07122 Palma de Mallorca, Spain
4Laboratorio de Sistemas Complejos, FI-UBA, Argentina and GIDyCSI-Universidad Nacional de Lanús, Argentina
Received 26 February 2013; Revised 17 June 2013; Accepted 18 June 2013
Academic Editor: Svetoslav Nikolov
Copyright © 2013 C. M. Sarris et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
We consider here an environment in which the fact that a semiquantum Hamiltonian obeys SU(2) symmetries poses serious difficulties if one wants to compute Lyapunov exponents.
An extreme complexity of phase space's trajectories that are very sensitive to small changes in the initial conditions is the signature of classical chaos, accompanied by (i) an ostensibly random allotment of phase points on a Poincare's surface of section and (ii) an exponentially rapid separation of two initially close trajectories . Instead, it is clear that the state vector of a closed quantum system cannot exhibit chaotic motion in Hilbert space. The interaction between a quantum system and a classical one may instead lead to authentic chaotic motion of the quantum component, a phenomenon known as semiquantum chaos [2, 3]. Remark that the vocable semiquantum is reserved to systems for which neither the quantum part nor the classical part would be chaotic by themselves. If chaos ensues, this happens because of the classical-quantum coupling. Semiquantumness implies that one part is treated classically and the other one in quantal fashion .
We consider here an environment in which the semiquantum Hamiltonian obeys symmetries and consider the difficulties that arise if one wants to compute Lyapunov exponents .
Consider a system that possesses both quantum and classic degrees of freedom, with a coupling amongst them, that we call semiquantum [3, 6–10]. The associated Hamiltonian is of the general form  where , , and are the quantum, classical, and interaction parts, respectively. There exist many situations in which a semiquantum description has been attempted [6, 9, 12]. Porter  made an exhaustive compilation of physical systems for which this kind of description is relevant. One may highlight vibrating quantum billiards as a useful abstraction of the ensuing semiquantum dynamics . Indeed, many semiquantum Hamiltonians are found in the literature [3, 14–23].
In this work, via the Maximum Entropy Principle (MEP) vantage point , we show that serious difficulties arise, in the case of these systems, if one wants to compute Lyapunov exponents, because the quantum degrees of freedom of the system must abide by the generalized uncertainty principle (GUP)  which in the case of an underlying Lie algebra becomes a dynamic invariant. In turn, the other dynamic invariants (the energy) involve both quantum and classical degrees of freedom. Accordingly, the initial condition pertaining to at least one of the classical degrees of freedom will depend not only upon the energy value but also on the uncertainty relation as well. This will be illustrated below.
These invariants allow us to (i) adequately define initial conditions and (ii) follow the details of the temporal evolution, from regular regimes to irregular ones. Our main invariant is none other than the generalized uncertainty principle (GUP) . From a MEP viewpoint, in the analysis of a semiquantum dynamics, one deals with a peculiar dynamic “working” space spanned by the variables , . The first variables are the mean values of a set of noncommuting observable closing a partial Lie semialgebra under commutation with the Hamiltonian. We insist on the fact that these quantum variables must obey the uncertainty relation not only at the instant but also for any other . The , are classical variables of the system. As in the full quantum case [25, 26], we will see for semiquantum nonlinear dynamics; the closure condition defines an dynamic matrix which governs the dynamics of the system's quantum degrees of freedom. This matrix is now of a semiquantal nature since it also depends upon classical degrees of freedom.
We follow the prescription given in [8, 9, 11, 12, 18, 19]; that is, the energy of the system is taken to coincide with the expectation value of the Hamiltonian, , traced over the quantum state, . In turn, generates the time evolution of the classical degrees of freedom in the orthodox classical mechanics' fashion. The MEP’s point of view approach to semiquantal systems takes advantage of two facts:(1)it is possible to describe the time evolution of the quantum degrees of freedom in the dual space of Lagrange multipliers associated to the quantum observable [25–27] and(2)the existence of the GUP invariant, say, makes it possible to analyze the dynamics of the system in different regimes (irregular and regular) by varying 's value through the initial conditions of the system.
3. Quantum Plus Classical Hamiltonians and Maximum Entropy Approach
Consider the Hamiltonian (1). The classical degrees of freedom are the canonical conjugate variables . We choose for the Hamiltonian (1) the specific form where the first term includes both the (nonchaotic) and (chaos-generator) ingredients, with the ’s being quantum operators, while the last term is a purely classical and nonchaotic, with functions of the canonically conjugate classical variables . We have in fact a family of Hamiltonians with a classical phase-space substratum.
The MEP formalism [25–27] deals with the quantum degrees of freedom of the systems (2) and provides a density operator [25–27] in terms of Lagrange multipliers (’s), where is the identity operator and is determined through the normalization condition. Alhassid and Levine demonstrated [25, 26] that the sufficient number of constraints Is necessary in order for the state operator (3) to be valid for any time deriveing from the closure condition Equation (5) defines a matrix whose coefficients may depend upon the time if the Hamiltonian of the system is time dependent. With the help of (5) and considering that the set of relevant operators defined through (5) are linearly independent, it is possible to obtain the equations of motion of the Lagrange multipliers [25–27] (see (8)) below. Note that (5) involves the classical variables in the guise of parameters, on account of the coefficients intervening in (2). Thus, (5) adopts the appearance with . The ’s are quantum operators closing a partial Lie algebra under commutation with (2) [25, 26]. The defines a semiquantum matrix .
Since we are going to describe the quantum state of the system (2) by means of the quantum statistical operator given by (3), it is required that this operator be represented by (3) for all so that the entropy [25, 26] is a constant of the motion [25, 26]. The Lagrange multipliers, in turn, fulfill the following equations of motion (see  for more details): The semiquantum matrix defined through (6) enables us to obtain the equations of motion of the mean values for the quantum degrees of freedom via  Equations (8) and (9) display the interplay between the quantum and classical degrees of freedom. In order to describe our classical degrees of freedom, we need to know the equation of motion for the classical variables .
To do so, we consider the mean value of (2), traced over the quantum state (3), as a classical Hamiltonian function so that generates the temporal evolution of the classical degrees of freedom [8, 9, 11, 12, 18, 19] where indicates Poisson brackets.
The generalized uncertainty principle (GUP) in the present framework reads 
The GUP relationship given by (11) is obtained as a sum over the principal minors (of order 2) of the correlation matrix , corresponding to the observable set defined by the closure condition (6). The correlation matrix is a definite positive one, and its elements are defined via 
3.1. The SU(2) Instance
Proposition 1. If a set of operators, which fulfills the commutation relation (13), closes a commutation algebra with a Hamiltonian of the type (14), then the semiquantum matrix of the system, defined by means of the closure condition (6), is an antisymmetric one.
Every Hamiltonian that closes an algebra with the generators is accompanied by the GUP invariant (11) which for the Lie algebra adopts the fashion  consider the following: Because of the uncertainty principle (11), , and Schwarz' inequality entails , that is, the uncertainty principle for the Lie algebra, that can be expressed in the guise defining the celebrated Bloch sphere of the system.
4. First Difficulty: The Normalization Difficulty
In evaluating Lyapunov exponents (LE), we must consider the action of the motion invariants (MI). Any “orbit” calculation must “respect them.” In particular, we must take care of not violating the uncertainty principle as the system evolves (same for the normalization procedure) in the process of Lyapunov exponents' evaluation. Once we have chosen a GUP-value at via (16), this value must remain the same at any other time so as not to violate the uncertainty principle. More specifically, we are speaking of the fiduciary trajectory and of the ones infinitely close to it that one keeps renormalizing when the Euclidean distance between them grows. If the uncertainty principle ((11) or (16)) is violated during the normalization procedure, the concept of semiquantum system (2) will become a mere “mathematical artifact” (in the words of ) but cease to have physical meaning. To illustrate what we mean, we consider the semiquantum Hamiltonian  whose equations of motion for quantum degrees of freedom (mean values and Lagrange multipliers) and classical degrees of freedom are
Our motion invariants (MIs) for the system given by (17) are (i)generalized uncertainty principle (GUP) (see (11)) (ii)energy () (iii)evolution vector (iv)entropy () with invariant in -space whose value is not an arbitrary one. The ’s values are adjusted after setting initial conditions (IC) for the UP , since the invariants' expressions in the two reciprocal spaces expectation values one and one are related by  All these invariants' values are closely linked to the uncertainty principle's value given by (26) and, therefore, cannot be chosen disregarding it. Accordingly, one must mind the evolution equations for , , , , , , , and given by (18)–(25). Let be a vector field belonging to the solution space of (18)–(20) and (24)-(25) and let be the flow induced by , that is, as it reads in : “for any let , where is an integral curve of the vector field such that .” Now let us choose two neighboring initial conditions and not lying on the same trajectory as prescribed by , but minding the conservation of the invariants given by (26) to (31) over the whole evolution process. Thus, we choose In (33) and (34), the initial conditions must be set as follows: (i) choose any value for our GUP invariant (26) that it does not violate the uncertainty principle (16), that is, . Once this value is chosen, it must remain the same during the whole evolution of the system not only for the fiduciary orbit (33) but also for any other infinitely close to it that one keeps renormalizing (in this way, we assure that the quantum state is always evolving on the same Bloch sphere, i.e., the same radius, in order not to violate the uncertainty principle, see (34)); (ii) choose arbitrarily two of the three following values: , , or , and the third must be obtained from the relationship given by (26); (iii) choose a value for the energy which is also another invariant of the motion; (iv) choose an initial value for one of the classical variables, for example, ; (v) it is clear that the value cannot be fixed at will since is of the form and its initial value is constrained to the value and to the uncertainty principle value! This is a curious feature of the semiquantum dynamics: the values of the classical degrees of freedom are influenced by the uncertainty principle followed by the quantum degrees of freedom of the system (see  for more details). The steps (i) to (v) given above give conditions to be followed by not only on the fiducial trajectory (33) but also by all trajectories infinitely close to it. Otherwise, the uncertainty principle will be violated. With such restrictions in mind, we choose ; ; ; ; .
We denote, as in (), by the length of the segment joining to () at , so that where is the Euclidean norm. Let and . Reference  prescribes: “denote by the unique point on the half-line issuing from and containing such that ” in order to be able to iterate this procedure and define and “and find the unique point on the half-line issuing from and containing such that ” in order to get a sequence which enables one to evaluate .
However, here the nonlinear semiquantum dynamics reigns. We must ensure that those required points are to be chosen in such a way that the invariants given by (26) to (31) are always respected. How could one ensure that imposing on the system the normalization condition at every step does not violate the uncertainty principle? Enters now our crucial caveat: this prescription may be inconsistent with that given by . Such constraint does not exist for a classical system since there is no obligation to strictly abide by the uncertainty principle. Such an obligation is not a mere mathematical feature but a strong physical constraint . Unfortunately, all we can expect is to be able to find a point (or for each iteration) connecting to with a distance from as small as possible, such that the function is minimal. This must comply with the uncertainty principle (a dynamic invariant for the Lie algebra) whose value was set at by means of the initial conditions.
Summing up, in order to find the initial condition of the fiduciary's very close trajectory (in each iteration path), we MUST externalize the distance function , subject to the following constraints:(i) should respect the fixed value of GUP invariant (see (26)) so as not to violate the uncertainty principle,(ii) does the same with the invariant (see (27)) so as not to violate energy conservation,(iii) does the same with the invariant (see (29)) so as not to violate entropy conservation,(iv) does the same with the invariant (see (28)) so as not to violate the conservation of the evolution vector,(v) must lie in the half-line starting at and containing . If we do mind all this panoply at each iteration path, we can never be sure that we are evaluating in a correct fashion. In general, we face an insurmountable task.
5. Second Difficulty: Numerical Simulations' Cumulative Error
If we evaluate Lyapunov exponents with the so-called one step method Runge-Kutta 4 (RK4), whose local truncation error is of the order (given that RK4 is a fourth order method), we face the following situation: suppose we choose an integration path . This entails that our approximation to the actual solution starts off (at ) with a local truncation error . However, as pointed out by , it is necessary to take into account the error arising as a consequence of cumulative effects coming from all local truncation errors, namely, the global truncation error . After a few millions of iteration paths, the global truncation error is so important that the temporal series' approximation that we have obtained becomes meaningless. It is the same for the swarm of points in Poincaré surfaces of section. Accordingly, such evaluation of the Lyapunov exponents numerical is also meaningless (remember that even if the normalized trajectory begins at every and does not accumulate much error, the fiduciary trajectory is evolving from and does accrue a large global cumulative truncation error ). To illustrate our claim, let us return to the semiquantum Hamiltonian given by (17): if we take ;;;;;;;, an integration path , and integrate (18)–(20) and (24)-(25), using RK4 from to , we obtain the Poincaré surface of section versus and the Bloch sphere of Figure 1, while if we do the same from to , using the same integration path, we obtain the Poincaré surface of section versus and Bloch sphere depicted in Figure 2.
We have shown that the evaluation of Lyapunov exponents on semiquantum systems like (17) may become an impossible task given that these kinds of Hamiltonian systems contain both a classical component and a quantum component which must obey the uncertainty principle. We appreciate via (18)–(25) the interplay between classical and quantum degrees of freedom. This means that once (i) the value is chosen (without violating the uncertainty principle) together with (ii) the energy value (another dynamic invariant), the classical variables become strongly constrained in the values they can attain by two facts: (i) the constant energy of the system determines the energetically accessible regions (this fact is shared by classical Hamiltonians) and (ii) the uncertainty principle, that imposes strong restrictions on their evolution (see  for more details). Thus, the value cannot be fixed without minding the condition that follows from the uncertainty principle. This fact is encountered only in dealing with semiquantum dynamics. Accordingly, one does not preserve the uncertainty principle's dynamic invariant value GUP and one will get nonsensical results. This is a distinctive feature of semiquantum dynamics that has no classical counterpart. Evaluating Lyapunov exponents for semiquantum systems is an entirely different problem for semiquantum systems than for classical ones.
In other words, as the MEP approach is able to convert the evolution equation of motion of the quantum state of a semiquantum system into a set of nonlinear semiquantum differential equations (see [21, 22] for more details), the quantum nature of the system rises its head in setting the initial conditions for the quantum degrees of freedom. One must not violate the uncertainty principle. At every when we obtain the values infinitely close to the fiduciary trajectory, we must be sure that the system abides by the uncertainty principle from the very beginning. In this sense, calculating the Lyapunov exponents for a semiquantum system is in no way a similar task to that of doing it classically.
This work was partially supported by the Project PIP1177 of CONICET (Argentina) and the Projects FIS2008-00781/FIS (MICINN)-FEDER (EU) (Spain, EU).
- A. J. Lichtenberg and M. A. Lieberman, Regular and Chaotic Dynamics, Springer, New York, NY, USA, 1992.
- R. Bluemel and W. P. Reinhardt, Chaos in Atomic Physics, Cambridge University Press, Cambridge, Mass, USA, 1997.
- L. E. Ballentine, “Is semiquantum chaos real?” Physical Review E, vol. 63, no. 5, Article ID 056204, 7 pages, 2001.
- G. Benettin, L. Galgani, and J.-M. Strelcyn, “Kolmogorov entropy and numerical experiments,” Physical Review A, vol. 14, no. 6, pp. 2338–2345, 1976.
- C. M. Sarris and A. N. Proto, “Generalized metric phase space for quantum systems and the uncertainty principle,” Physica A, vol. 377, no. 1, pp. 33–42, 2007.
- M. A. Porter, “Nonadiabatic dynamics in semiquantal physics,” Reports on Progress in Physics, vol. 64, no. 9, pp. 1165–1189, 2001.
- R. Blumel and B. Esser, “Quantum chaos in the Born-Oppenheimer approximation,” Physical Review Letters, vol. 72, no. 23, pp. 3658–3661, 1994.
- H. Schanz and B. Esser, “Mixed quantum-classical versus full quantum dynamics: coupled quasiparticle-oscillator system,” Physical Review A, vol. 55, no. 5, pp. 3375–3387, 1997.
- A. M. Kowalski, A. Plastino, and A. N. Proto, “Semiclassical model for quantum dissipation,” Physical Review E, vol. 52, no. 1, pp. 165–177, 1995.
- J. Ma and R. K. Yuan, “Semiquantum Chaos,” Journal of the Physical Society of Japan, vol. 66, no. 8, pp. 2302–2307, 1997.
- A. M. Kowalski, M. T. Martin, J. Nuñez, A. Plastino, and A. N. Proto, “Quantitative indicator for semiquantum chaos,” Physical Review A, vol. 58, no. 3, pp. 2596–2599, 1998.
- A. M. Kowalski, A. Plastino, and A. N. Proto, “A semiclassical statistical model for quantum dissipation,” Physica A, vol. 236, no. 3-4, pp. 429–447, 1997.
- M. A. Porter and R. L. Liboff, “Vibrating quantum billiards on Riemannian manifolds,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 11, no. 9, pp. 2305–2315, 2001.
- T. C. Blum and H.-T. Elze, “Semiquantum chaos in the double well,” Physical Review E, vol. 53, no. 4, pp. 3123–3133, 1996.
- L. L. Bonilla and F. Guinea, “Collapse of the wave packet and chaos in a model with classical and quantum degrees of freedom,” Physical Review A, vol. 45, no. 11, pp. 7718–7728, 1992.
- R. I. Cukier and M. Morillo, “Comparison between quantum and approximate semiclassical dynamics of an externally driven spin-harmonic oscillator system,” Physical Review A, vol. 61, no. 2, 4 pages, 2000.
- A. K. Pattanayak and W. C. Schieve, “Semiquantal dynamics of fluctuations: ostensible quantum chaos,” Physical Review Letters, vol. 72, no. 18, pp. 2855–2858, 1994.
- A. M. Kowalski, M. T. Martin, J. Nuñez, A. Plastino, and A. N. Proto, “Semiquantum chaos and the uncertainty principle,” Physica A, vol. 276, no. 1, pp. 95–108, 2000.
- A. M. Kowalski, M. T. Martin, A. Plastino, A. N. Proto, and O. A. Rosso, “Wavelet statistical complexity analysis of the classical limit,” Physics Letters, Section A, vol. 311, no. 2-3, pp. 180–191, 2003.
- F. Cooper, J. F. Dawson, D. Meredith, and H. Shepard, “Semiquantum chaos,” Physical Review Letters, vol. 72, no. 9, pp. 1337–1340, 1994.
- A. M. Kowalski, A. Plastino, and A. N. Proto, “Classical limits,” Physics Letters, Section A, vol. 297, no. 3-4, pp. 162–172, 2002.
- C. M. Sarris and A. N. Proto, “Information entropy and nonlinear semiquantum dynamics,” International Journal of Bifurcation and Chaos, vol. 19, no. 10, pp. 3473–3484, 2009.
- A. M. Kowalski, M. T. Martin, A. Plastino, and A. N. Proto, “Classical limit and chaotic regime in a semi-quantum Hamiltonian,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 13, no. 8, pp. 2315–2325, 2003.
- C. M. Sarris, F. Caram, and A. N. Proto, “The uncertainty principle as invariant of motion for time-dependent Hamiltonians,” Physics Letters, Section A, vol. 324, no. 1, pp. 1–8, 2004.
- Y. Alhassid and R. D. Levine, “Entropy and chemical change. III. The maximal entropy (subject to constraints) procedure as a dynamical theory,” The Journal of Chemical Physics, vol. 67, no. 10, 1977.
- Y. Alhassid and R. D. Levine, “Connection between the maximal entropy and the scattering theoretic analyses of collision processes,” Phys. Rev. A, vol. 18, no. 1, pp. 89–116, 1978.
- D. Otero, A. Plastino, A. N. Proto, and G. Zannoli, “Ehrenfest theorem and information theory,” Physical Review A, vol. 26, no. 3, pp. 1209–1217, 1982.
- E. Merzbaher, Quantum Mechanics, Wiley, 1963.
- W.-K. Tung, Group Theory in Physics, World Scientific Publishing, Singapore, 1985.
- W. Cheney and D. Kincaid, Métodos Numéricos y Computación, CENGAE Learning, Mexico, 2011.