Uncertainty Principle and Bifurcations in the SU ( 2 ) Nonlinear Semiquantum Dynamics

In this paper, a nonlinear semiquantum Hamiltonian associated to the special unitary group SU(2) Lie algebra is studied so as to analyze its dynamics. The treatment here applied allows for a reduction in: 1) the system’s dimension, as well as 2) the number of system’s parameters (to only three). We can now discern clear patterns in: 1) the complete characterization of the system’s fixed points and 2) their stability. It is shown that the parameter associated to the uncertainty principle, which constitutes a very strong constraint, is the key one in determining the presence of fixed points and bifurcation curves in the parameter’s space.


Introduction
Semiquantum Dynamics (SD) may be used to describe systems in which quantum and classical degrees of freedom coexist.One finds in [1] an exhaustive compilation of physical phenomena and technological applications successfully modeled by SD.It is also possible to encounter situations in which SD is used to describe physical phenomena [2] [3] [4] [5].A nonlinear semiquantum Hamiltonian associated to the SU(2) Lie algebra is very useful to model the problem of quantum confinement, which is of interest for nanotechnology and solid state physics.In particular, if the quantum subsystem is associated to the SU(2) Lie algebra generators { } ˆˆ, , simple form and turns out to be a motion invariant [6].
In [7], the authors consider the following semiquantum Hamiltonian where ˆx σ and ˆz σ are quantum operators, the x and z components of a 1/2 spin particle, while q and p are canonical conjugated classical variables (position and momentum) that obey the Hamilton equations [5] [8].B, C, m, D, and F are positive and constant parameters.The Hamiltonian given by Equation (1) represents a quantum 1/2 spin particle interacting with an external magnetic field [2] [7] (due to the term ˆz Bσ ).The particle is confined by the double well potential ( ) sical and quantum variables are couple in non-linear fashion via the term, ˆx Cqσ , C being the coupling constant.In [7], the authors concentrate on the likely presence of chaotic motion (semiquantum chaos) for special values of the coupling strength.The authors represent the trajectories for different initial conditions by stroboscopic plots, displaying regular and irregular dynamics.This Hamiltonian also may be reduced to the one in [2] (taking 0 F = ).It also can be used to model the semiquantum differential equations of the spin-boson Hamiltonian of [9] (taking 0 D = ).In [10], the authors considered the simplest case of a spin-boson Hamiltonian, i.e., a two level system coupled to one oscillator degree of freedom, and made a posterior semiclassical approximation, to obtain a semiquantum Hamiltonian similar to that given by Equation (1).
We consider that, since the Hamiltonian of Equation ( 1) is able to model the quantum confinement phenomenon, its dynamics deserves an exhaustive analysis.In the present work, we purport to give a full description of its phase space taking into account that, in conservative systems like this one, the motion is fully determined by the amount and disposition of its fixed points.We make a complete characterization of them and determine their stability according to the system's parameters values.In addition, we obtain the bifurcation curves which divide the phase space into the three different zones in which the fixed points are located, according to their stability.We present a dimensionless formulation, and because the uncertainty principle (UP) is an invariant of the motion for the nonlinear semiquantum Hamiltonians associated to the SU(2) Lie algebra [6] [11], we make a change of variables that removes it as external strong constrain to the system's motion equations.The UP provides then just an additional parameter in the concomitant new motion equations.We mention that, in [10], a similar method is used for a coupled quasiparticle-oscillator system, which enables the authors to study the existence of fixed points and bifurcation curves, allowing for a formulation in canonically conjugate variables of the excitonic subsystem.In our case, the change of variables offers some advantages which are highlighted in describing our treatment and summarized in the conclusions.

Equations of Motion
If we consider the generators of the SU(2) Lie algebra, { } ˆˆ, , x y z σ σ σ , it can be easily seen that they close a partial Lie algebra under commutation with the Hamiltonian of Equation ( 1), since the commutator of any ˆi σ with Ĥ may be expressed as a linear superposition of these generators.Following the prescription given in [6] [12] [13], we obtain the equations of motion for the system given by Equation (1): and they must obey the uncertainty relation [6] [11] which, for the SU(2) Lie algebra case, adopts the form [6]: We will consider the whole range of values that the polarization vector σ can achieve in the interval ( ) must remain a constant of motion during the whole temporal evolution [6].

The System's Fixed Points
In order to determine the system's fixed points we proceed, as usual, by equating ( 11)- (15) to zero.Note that the existence of them does not depend on the parameter ε , but only on the values of r and s.From Equation (13), there exist two situations: the cases 0 q * = and 0 q * ≠ .

3.2.
Case q * ≠ 0 From Equations ( 11)-( 15) it follows that these kinds of fixed points must fulfill: ( ) 0, 0, and , y p x q s q z s q * * * = = = − = − (18) so they adopt the appearance ( ) , and must obey the strong constraint given by Equation ( 16), the uncertainty principle, which in terms of q * , reads: , q q r s * * = are the roots of Equation ( 19), which should be tackled numerically.The system's phase space, given by Equations ( 11)-( 15), is five-dimensional.However, the Jacobian matrix at the fixed point must have, at least, one null eigenvalue, since the uncertainty condition of Equation ( 16) is an external constraint added to the equations of motion.
This means that, in fact, the solutions lie on a 4D invariant manifold,  .
Accordingly, we represent the quantum degrees of freedom, [ ] , , x y z , in spherical coordinates: , , , , , , Using the relations: we find: The system of Equations ( 11)-( 16) becomes now: where ( ) ( ) The system's energy given Applied Mathematics by Equation ( 17) becomes: We claim that this change of variables (CV) offers some advantages: 1) The quantum variables θ and ϕ obey the relationship: as if they were canonical spherical-conjugates, meanwhile the classical ones, q and p, obey, as usual, the Hamilton's equations [6] [13]: Thus, the fixed points * ξ of the system given by Equations ( 23)-( 26), are the critical points of the energy function h of Equation ( 27), since: 2) The uncertainty relation in Equation ( 16), now of the form ( ) , is incorporated into the new system's equations in a natural way by reducing the system dimension and removing a superfluous null eigenvalue.One may speak of linearization of the original system at each fixed point.
3) Our CV provides a better characterization of the fixed points, since with this change we find a generic expression for them explicitly written in terms of the variable q.This fact facilitates the study of 1) the presence of bifurcations curves in the parameter's space, and 2) the stability analysis of the fixed points.
To obtain the fixed points, From 0 q = in Equation (25), it follows that * 0 p = , and from 0 ϕ = in Equation ( 24), we are led to 0 0 2 π k

Case k
Then, in this case, * ξ adopts the generic form ( ) From Equation (32), ( ) ( ) ( ) ( ) and then: ( ) ( ) Replacing Equation (33) into Equation (35), it follows that * 0 q must be the solution of the following equation: . 1 r q q s q q Since * 0 q ≠ , the right side on Equation (36) is positive.Thus, the possible values of * 0 q are restricted to the range ( ) ( ) In each case, the amount of solutions will be obtained graphically by means of the intersection points of two curves, namely: , 0 and ., 0 1 q s q r f q g q s q q q 1 depicts the intersection between the left branch of ( ) f q and ( ) g q .It is possible to see that for all , 0 r s > , there is always a unique solution, ( ) , in the range of interest, and so a corresponding fixed point, ( ) , is obtained.Note also that this point does not bifurcate in the ( ) , r s -parameter space.
• 0 q > .Here, we are looking for the intersection points of the right branch of ( ) f q and ( ) g q .The situations for different values of r and s are depicted in Figure 2, by considering r as a "fixed'' parameter and s going down from s r > to s r < .q existing for 0 q < .Applied Mathematics Decreasing the s value, the curves become tangent at 0 q = , at the critical value s r = .Note that s is a "rigid'' parameter since only produces f graph shifts, in contrast to the "flexible'' parameter r which bends the graph of g.The scenario splits into two ones when s r < , depending on the convexity of g as compared to the f value at 0 q = , i.e., depending on the absolute values of the second derivatives: Thus, for , see Figure 2(c)), there is no intersection point between f and g.However, for s r , see Figure 2 ξ , is created (the term "imperfect" means that only one point is created/extinguished, instead of two, as would happen in a "perfect" saddle-node.This is due to the loss of the quadratic symmetry of 2 s q − when considering only its right branch ( 0 q > ) [15].This last scenario ( 2 r > ) persists until the two graphs becomes tangent again at the critical value ( ) q and 3 * 0 q meet each other at ( ) , , q f q q g q = the tangency point.For s below this critical value s , the two graphs do not match any longer, and the fixed points From this, we obtain the relation between q and r: ( ) valid for 2 r > .Evaluating f and g at Equation (40), and so equating ( ) g q , we find the relation between r and s: ( ) or, equivalently, ( ) Replacing Equation (41) into Equation (40), we find the relation between q and s , that allows to better fit the possible ranges for 2 * 0 q and 3 * 0 q , within the range ( ) 0, s , which will be useful later on: or, equivalently:

. 3
Thus, as , we have for s r < : and , . 3 3 Note that the bifurcation curve ( ) s r meets the line s r = tangentially at a codimension-2 point ( ) 2, 2 in the ( ) , r s -parameter space (Figure 3).

Case k
Accordingly, in this case * ξ adopts the generic form ( ) 33) is still valid, but from 0 p = in Equation ( 26 ( ) ( ) Replacing ( 33) into (47), it follows that * π q must be the solutions of the equation: ( ) Accordingly, in this case the possible values of * π q are restricted to the range . Therefore, the solutions should now be obtained by the intersection points of the functions ( ) f q − and ( ) g q (defined in (37)) according to the r and s values.It is easy to see that, for symmetry reasons, the analysis is completely analogous to the one developed in the previous Section.
Then, the solutions thus obtained are: ( )  Summing up, the amount of fixed points in the ( ) , r s -plane is (Figure 3):

The Stability of the Fixed Points
The system given by Equation (1) (and then by Equation ( 27)) is conservative.Thus, the local behavior at the fixed point may be studied by considering the energy h of Equation ( 27) as a function of

[ ]
, , , which is a constant of motion, since 0 h = .The h-level sets, ( ) 3D invariant manifolds containing orbits given by the equations of motion.The h-level sets help in understanding the structure of the system's phase space.The local behavior around a fixed point will be determined by the index k of * ξ as a non-degenerate critical point of the Morse function h [16].
Definition 1.The instability index ( ) is the number of negative eigenvalues of the Hessian matrix of h at * ξ , ( ) ξ is a non-degenerate critical point, i.e.
( ) This is to say, k is the number of independent directions along which h decrease from * ξ .Therefore, if ( ) , the h-level set around * ξ , corres- ponds to a positive definite quadratic form, * ξ is a local minimum of h and a system's nonlinear center.If, in our case, ( ) , * ξ it is unstable, and of the saddle type.

Obtaining the Index k of * ξ
From Equation (27) the Hessian matrix of h at * ξ is: Since H h * ξ is a square-block diagonal matrix, its eigenvalues are those from the block submatrices.Taking into account that for a 2 2 × -matrix the eigenvalues may be written in terms of its trace and determinant, we have, then, that the four eigenvalues of H h * ξ are: ( ) ( ) ( )  ( ) ( )( ) for all * ξ , as befits to a conservative system which cannot posses a completely unstable critical point.
Figure 4 and Figure 5 illustrate a case relative to the orange zone in Figure 3, for which the six fixed points coexist ( 11, 9 r s = = ).Table 1.The nondegenerate critical points * ξ and their instability index k:

Note
The other two fixed points of the system, [ ] 0, 0, , 0, 0 r ± , which are not described by the change of variables ( [ ] , , , q p θ ϕ = ξ ), result degenerate critical points of the energy h of Equation ( 17) as a function of [ ] , , , , x y z q p , since they force the vanishing of the Hessian matrix of h .

Conclusions
The dimensionless formulation of Equations ( 2)-( 6), given by Equations ( 11)-( 15), allows for reduction in the number of system's parameters to three: r, s, and ε .The posterior change of variables (CV) has additional advantages as highlighted before: the reduction of the system's dimension through proper consideration of the uncertainty principle constraint.It also affords a convenient way to obtain the fixed points, providing a better characterization of them, according to the system's parameters, and depending only on the classical degree of freedom, q.We can display a relation between the new quantum degrees of freedom, θ and ϕ , which turn out to be canonically "spherical-conjugates", making the system's fixed points to be the critical points of the energy function h .In addition, the CV illustrates the role of each parameter in the system dynamics in a very clear fashion: ε can only play a role in the stability analysis, ensuring that the unstable fixed points are saddles.The parameter s accounts for the range of possible values that q can achieve at each fixed point, dividing the phase space into three regions.Finally, the parameter r forces the uncertainty principle to play a fundamental role in the appearance/disappearance of fixed points, thus governing the presence of bifurcations in the SU(2) nonlinear semiquantum dynamics.
The putative presence of chaotic dynamics in this system, for some ( ) , r sparameter's region, is still under consideration and will be reported elsewhere.
Remark: by virtue of the SU(2) Lie algebra, the uncertainty principle becomes a constant of the motion.Thus, we claim that the methodology used in the present work applies even in the case in which the quantum subsystem of Equation (1) was nonlinear in the spin variables (as in the non-dissipative Hamiltonian case treated in [17], which can be used to model a SQUID).Despite the nonlinearity in the quantum subsystem, the uncertainty principle would remain there a constant of the motion.This and other topics related to semiquantum dynamics will be part of future work.
, the uncertainty principle (UP) adopts a very R. Hansen et al.DOI: 10.4236/am.2018.91001 2 Applied Mathematics generated by a classical particle of mass m and it undergoes elastic reflections between the moving frontier, V ∂ , of the double well potential.The term 2 2 p m represents the classical kinetic energy.The clas-

==
in Equation(11), it follows 0 y = .From d 0 d y τ = in Equation (12), it follows 0 x = .So, from Equation (16), we are led to z r = ± .From Applied Mathematics in Equation (14), it follows 0 p = .Therefore, in this case, one obtains two fixed points, and south pole of the dimensionless Bloch sphere.
the system's fixed points by means of the four-dimensional variables,

Figure 1 .
Figure 1.The graph of the functions g and the left branch of f of Equation (37) for , 0 r s > , showing the unique intersection point 1 * 0

Figure 2 .q
Figure 2. The graph of the functions g and the right branch of f of Equation (37) for , 0 r s > .The figures in the first column ((a)-(c)) depict the case where f is more convex than g at 0 q = ( (d)), the graphs intersect at a new point, to another solution in the range of interest.This fact makes the line s r = to be a codimension-1 bifurcation curve in the ( ) , r s -parameter space: an imperfect saddle-node bifurcation occurs for 2 r < , since the fixed point 2 * 0 ξ is lost, and another imperfect saddle-node bifurcation R. Hansen et al.DOI: 10.4236/am.2018.

ξ
are mutually destroyed.The system undergoes a ("perfect") saddle-node bifurcation at s , and this makes the relation curve in the ( ) , r s -parameter space.Note the fundamental role that plays the uncertainty principle parameter r in the coming into being of these bifurcation curves.To find the relation ( ) s s r = at q , we proceed by equating the first derivatives:

Figure 3 .
Figure 3.The ( ) , r s -parameter space and the codimension-1 bifurcation curves: the and the same bifurcation curves, s r = and , the matrix in Equation (49) for * 0 ξ , is thus expressed in terms of * q :