Arxiv:1702.03077V1 [Quant-Ph] 10 Feb 2017 Deformation and the Equivalence Proposed

Home , Quantum group

Mapping of the 2+1 q-deformed Dirac oscillator onto the q-deformed Jaynes-Cummings model: It’s non-relativistic limit and Zitterbewegung eﬀect

Parisa Majari,1, ∗ Alfredo Luis,2, † and Mohammad Reza Setare1, ‡ 1Department of Science, University of Kurdistan, Sanandaj, Iran 2Department of Optics, University Complutense, Madrid, Spain (Dated: February 13, 2017) We develop the equivalence between the two-dimensional Dirac oscillator and the anti-Jaynes- Cummings model within a q-deformed scenario. We solve the Hamiltonian spectrum and the time evolution for number and coherent initial states, focusing on the appearance of the Zitterbewegung eﬀect. We show the lack of preservation of the q-deformed versions of the total angular momentum that reproduces a collapse-revival structure. We provide suitable relations for the non relativistic limit.

PACS numbers: 42.50.Pq Cavity quantum electrodynamics; micromasers 42.50.Vk Mechanical eﬀects of light on atoms, molecules, electrons, and ions, 03.65.Pm Relativistic wave equations, 03.65.Ca Formalism

I. INTRODUCTION On the other hand, q-deformation provides a quite simple but rich and powerful enough way to describe complex The introduction and proper combination of suitable situations such as electronic conductance in disordered metals and doped semiconductors, phonon spectrum in physical models provide a significant advance in many 4 different areas of physics. In this work we combine the He, oscillatoryrotational spectra of diatomic and multi- relativistic Dirac oscillator, the Jaynes-Cummings model, atomic molecules [11–13]. Deformation is accomplished and the idea of q-deformation. by entering a dimensionless parameter q in normal Weyl algebra. This grants a more general case of the ancient The Dirac oscillator provides a generalization to the theory at the same time it replicates the primary the- relativistic domain of one of the most applied models in ory for q = 1. For example, we may say that relativistic classical and quantum physics [1, 2]. This is a notably quantum mechanics is a q-deformed Newtonian quantum simple solvable equation, linear both in momentum and mechanics, where the deformation parameter is the veloc- position, formally looking as a simple minimal coupling ity ratio q = v/c. Whenever q tends to zero, relativistic of a charged particle with a magnetic field. Accordingly, quantum mechanics convert to the Newtonian quantum it has found applications in nuclear physics, quantum mechanics. chromodynamics, many-body theories, supersymmetric Thus, in this work we investigate the q-deformed per- relativistic quantum mechanics and quantum optics [3– spective of the interplay between Dirac oscillator and 6]. Jaynes-Cummings model following Ref. [14], both in the This relativistic generalization is accompanied of new relativistic domain and in the non-relativistic limit. In effects such as the unavoidable emergences of spin, nega- particular we focus on the preservation in the q-deformed tive energies and rather unexpected behaviors such as the scenario of peculiar relativistic features such as the Zit- Klein tunneling and the Zitterbewegung as the result of terbewegung effect. interference between positive- and negative-energy components. A further key feature, is that the Dirac oscillator contains in itself another germane model in quantum II. THE q-DEFORMED MODEL physics: the Jaynes-Cummings model (JC), or, more rig- orously, the anti-Jaynes-Cummings version (AJC) [7, 8]. Let us present the main results regarding the q- This describes the coupling of a discrete spin-like variable

arXiv:1702.03077v1 [quant-ph] 10 Feb 2017 deformation and the equivalence proposed. with a bosonic degree of freedom [9]. The powerfulness of the model relies on its simplicity and its capacity to embrace very diﬀerent physical situations such as light- A. q-Deformed observables and states matter interaction, dynamics of trapped ions, or Bose- Einstein condensates [10]. In particular this model allows to use quantum-optics based technology or trapped ions A general case of modiﬁed Weyl algebra holds as fol- technology for simulating relativistic quantum mechanics lows [12, 15]: [3, 4]. † † aqaq − qaqaq = 1, for q ≤ 1, (2.1)

† where aq and aq are annihilation and creation operators, ∗Electronic address: [email protected] that operate such that: †Electronic address: alluis@fis.ucm.es ‡ p † p Electronic address: [email protected] aq|ni = [n]|n − 1i, aq|ni = [n + 1]|n + 1i, (2.2) 2 where the number basis |ni is actually defined by the 1.5 eigenvalue equation Q 1.0 † [ˆn]|ni = [n]|ni, [ˆn] ≡ a aq, (2.3) q 0.5 and [n] satisfies this relation [13]: 0.0 1 − qn [n] = . (2.4) -0.5 1 − q Qq 2 For q → 1 naturally [n] → n while for q < 1 we have 0.5 1.0 1.5 2.0 2.5 3.0 | α | [n] < n, with [n] → 1 as q → 0. Alternatively to the q-number operator [ˆn] in Eq. (2.3) we can define also an operatorn ˆ satisfying the standard 2 FIG. 1: Mandel parameters Qq, Q as functions of |α| for commutation relation with aq, i. e., q = 0.25 (solid line) and q = 0.75 (dashed line).

nˆ|ni = n|ni, [aq, nˆ] = aq. (2.5) being We can present two alternative realizations of the q- ∞ deformed algebra. On the one hand we have the following 1 X αn 1 |αi = |ni, |α|2 < , (2.12) differential form for aq p 2 [n]! 1 − q eq(|α| ) n=0 −2iαz iα d −iαz e − e dz e where [n]! = [n][n − 1] ... [1] with [0]! = 1 and the follow- aq = √ , (2.6) −i 1 − e−2α2 ing relation has been used ∞ where α = p− ln q/2 and z is a Cartesian coordinate. X zn eq(z) = . (2.13) On the other hand, after [15] we may express aq also as [n]! n=0 r[n + 1] These relations hold provided that |α|2 < 1/(1 − q) as a a = F a†a a, F (n) = , (2.7) q n + 1 necessary condition such that the number coefficients in Eq. (2.12) tend to zero as n tend to infinity. where naturally a = aq=1. In particular this means Regarding the number statistics let us show that these that q-deformation implies a highly nonlinear interaction. state are sub-Poissonian for the q-deformed number op- The degree of nonlinearity may be reduced to more sim- erator [ˆn] since ple levels if q → 1. Thus, considering q = exp(−), a power series expansion as → 0 leads to hα|[ˆn]|αi = |α|2, hα|[ˆn]2|αi = |α|2 + q|α|4. (2.14) We can introduce the Mandel Q parameter in the q- F a†a ' 1 − a†a. (2.8) 4 deformed scenario as ∆2[ˆn] Therefore, from a quantum-optical perspective q- Q = − 1 = (q − 1)|α|2, (2.15) deformation corresponds to optics in nonlinear media q h[ˆn]i [20]. Also, the limit (2.8) recalls the dynamics of trapped so that Q < 1 ↔ q < 1. On the other hand, regarding ions where a describes the motion of the center of mass, q the other number operatorn ˆ in Eq. (2.5) we have ob- as shown for example in Ref. [21]. tained numerically that the q-deformed coherent states On the other hand, in the opposite limit of large non- (2.12) are always super-Poissonian, that is linearity q → 0 we get ∆2nˆ † 1 1 Q = − 1 > 0. (2.16) F a a → √ , aq → √ a, (2.9) hnî 1 + a†a 1 + a†a This is illustrated in Fig. 1 where we plot both Mandel so that 2 parameters Qq, Q as functions of |α| for q = 0.25 (solid line) and q = 0.75 (dashed line). aq|ni → |n − 1i, aq|0i = 0. (2.10)

This is to say that in this limit aq becomes the Susskind- B. The q-deformed Dirac oscillator Glogower phase operator [22]. That is that aq tends to represent the complex exponential of the oscillator phase. Besides the number states (2.3) we consider some The Hamiltonian of the Dirac oscillator in two spatial analogs of the coherent states deﬁned via the eigenvalue coordinates x, y is given by [14] equation [17]: √ mc2 i 4mc2 ω a† HD = √ ~ ` , (2.17) 2 2 aq|αi = α|αi, (2.11) −i 4mc ~ω a` −mc 3

JC 0 0 + iφ0 − † −iφ0 where m is the rest mass of the electron, ω is the Dirac- H = ~δ σz + ~η (σ ae + σ a e ), (2.25) oscillator frequency, c is the speed of light, and a is the ` 0 0 left-handed combination of x, y dynamics [3] where η, η are coupling constant, δ, δ represent detuning, and φ, φ0 are arbitrary phases. ax + iay Now we q-deforme these JC and AJC Hamiltonians a` = √ . (2.18) † 2 simply by substituting a, a by their the q-deformed versions of creation and annihilation operators:

Now we can consider the q-deformed version of the AJC + † iφ − −iφ above Hamiltonian [14]: Hq = ~δσz + ~η(σ aqe + σ aqe ), (2.26) " √ # mc2 i 4mc2 ω a† JC 0 0 + iφ0 − † −iφ0 HD = √ ~ `,q , (2.19) Hq = ~δ σz + ~η (σ aqe + σ aqe ). (2.27) q −i 4mc2 ω a −mc2 ~ `,q In general terms, q-deformed JC Hamiltonians have been or equivalently, studied in [19]. √ D 2 2 † Hq = mc σz + i 4mc ~ω σ+a`,q − σ−a`,q , (2.20) D. Equivalence where σ refers to the corresponding Pauli matrices. Al- A readily inspection of Eqs. (2.20) and (2.26) reveals ternatively, an exact mapping between them if we properly identify √ † † the complex-amplitude operators D 2 2 0 ax,q − iay,q Hq = mc σz + i 4mc ~ω . −ax,q − iay,q 0 D AJC Hq ≡ Hq aq ≡ a`,q, (2.28) (2.21) Within the q-deformed scenario we can deﬁne the z with the following correspondence of parameters components of the orbital L, spin S, and total J angular √ 2 2 momentum operators: ~η = 4mc ~ω, φ = π/2, ~δ = mc . (2.29) Alternatively, if we consider explicitly the two- † † ~ L = a a − a a ,S = σ ,J = L +S , dimensional nature of our Dirac oscillator the equivalence z ~ r,q r,q `,q `,q z 2 z z z z (2.22) reads where the right-handed mode ar,q is deﬁned as D AJC AJC Hq ≡ Hq,x + Hq,y , (2.30) ax,q − iay,q with ar,q = √ , (2.23) 2 √ 2 ~ηx = ~ηy = 2mc ~ω, and we will always assume that this mode is in the vac- 2 φx = π/2, φy = 0, ~δx = ~δy = mc . (2.31) uum state. D Next we compare the q-deformed Dirac oscillator Hq with the q-deformed Jaynes-Cummings and anti-Jaynes- III. ENERGY SPECTRUM OF THE Cummings models. q-DEFORMED DIRAC OSCILLATOR

Next we derive the energy spectrum of the q-deformed C. The q-deformed JC and AJC Hamiltonians D Dirac oscillator Hq . The spectrum must be calcu- lated from scratch since the commutation relations have The JC model is an example of Rabi model in the changed in comparison with the q = 1 case. The q- framework of quantum electrodynamics that presents the deformed time-independent Dirac equation is: coupling of two energy levels, described by the Pauli ma- D trices σ, with an harmonic oscillator, described by the Hq |ψi = E|ψi, (3.1) complex-amplitude or annihilation operator a. As two D typical realizations the oscillator we can mention a quan- with Hq in Eq. (2.20). Looking for spinor solutions, tum electrodynamic ﬁeld mode, i. e., photons, or the |ψ i motion of the center of mass of a trapped ion, that is |ψi = 1 , (3.2) |ψ i phonons. This last case is specially interesting since it 2 allows to simulate both JC and AJC models simply by we get that Eq. (3.1) is equivalent to the pair of equations properly selecting the frequency tuning of a laser ﬁeld √ 2 2 † illuminating the ion, leading to the AJC and JC Hamil- mc |ψ1i + i 4mc ~ωa`,q|ψ2i = E|ψ1i, tonians [10] (3.3) √ AJC + † iφ − −iφ 2 2 H = ~δσz + ~η(σ a e + σ ae ), (2.24) −mc |ψ2i − i 4mc ~ωa`,q|ψ1i = E|ψ2i. 4

2 We have a readily solution with E = E0 = mc , |ψ2i = 0, This allows us to interpret the whole process as a series and |ψ1i = |0i, that is of Mach-Zehnder interferometers in parallel, one for each subspace Hn. In this picture the basis vectors |ni|χ↑i, |E0i = |0i|χ↑i. (3.4) |n − 1i|χ↓i play the role of input/output beams, the operators exp(±iθσ ) play the role of input/output beam Whenever E 6= mc2 combining Eqs. (3.3) we get : x splitters, while exp(−iEntσz/~) represents the phase difference acquired within the two arms of the interferome- (E2 − m2c4)|ψ i = 4mc2 ωa† a |ψ i, (3.5) 1 ~ `,q `,q 1 ter. In other words, the beam splitters provide the eigenstates of the Hamiltonian, while the phase difference is given by the eigenvalues. This also allows suitable simple (E2 − m2c4)|ψ i = 4mc2 ωa a† |ψ i, (3.6) 2 ~ `,q `,q 2 expressions regarding the approximations involved in the and we recall that in terms of q-deformed number opera- non-relativistic limit to be considered below. † † tors [ˆn] = a`,qa`,q and [ˆn + 1] = a`,qa`,q. From the above equations the energies can be easily obtained as follows: IV. ZITTERBEWEGUNG EFFECT ω E = ±E = ±mc2p1 + 4ξ[n], ξ = ~ , (3.7) n mc2 The Zitterbewegung effect is a trembling motion of relativistic particles that occurs due to the interference be- for natural numbers n 6= 0. The positive and negative of tween positive and negative energies [16]. Here we want the energy eigenstates | ± Eni are: to investigate this effect for the following q-deformed initial state:  q 2  En±mc |ni | ± E i = 2En . (3.8) n  q 2  |ψ0i = |n − 1i|χ↓i = ic−| + Eni − ic+| − Eni, (4.1) ∓i En∓mc |n − 1i 2En for n ≥ 1. The evolution of this state can be obtained The above result can be rewritten in terms of the eigen- by: † † states of σz, |χ↑i = (1, 0) and |χ↓i = (0, 1) as:

s En 1 p 2 2 En ± mc ωn = = mc 1 + 4ξ[n], (4.3) c± = . (3.10) ~ ~ 2En is the frequency of Zitterbewegung oscillation. The os- Let us present an alternative picture of this diagonal- cillation can be well appreciated for example in the ization and time evolution as an interference eﬀect [4]. mean value of the spin and angular momentum opera- To this end we express the total Hilbert space as direct tors (2.22): sum of one- and two-dimensional subspaces as 4ξ~[n] n−1 2 ∞ hLzit = −~[n − 1] − q sin (ωnt), (4.4) H = ⊕n=0Hn (3.11) 1 + 4ξ[n] where the subspaces Hn are deﬁned as: ~ 4ξ~[n] 2 hSzit = − + sin (ωnt), (4.5) Hn = span{|ni|χ↑i, |n − 1i|χ↓i}, H0 = span{|0i|χ↑i}. 2 1 + 4ξ[n] (3.12) and Leaving aside the trivial subspace H0, within each of the two-dimensional spaces Hn the Dirac oscillator (2.20) can 4ξ [n] hJ i = −~(1 + 2[n − 1])+ ~ sin2(ω t)(1−qn−1), be written simply as, in the basis in Eq. (3.12), z t 2 1 + 4ξ[n] n D 2 p (4.6) Hn,q = mc (σz − 4ξ[n]σy). (3.13) or, equivalently,

This is equivalent to ~ hJzit = − 2 (1 + 2[n − 1]) n n−1 HD = E e−iθnσx σ eiθnσx (3.14) 2 (1−q )(1−q ) n,q n z +4ξ~sin (ωnt) (1−q+4ξ−4ξqn) . (4.7) p where tan(2θn) = 4ξ[n]. Thus the unitary time evolu- Several interesting features can be noticed in these re- tion operator within each subspace Hn is given by: sults: i) As an interference eﬀect we can note that the visi- D −iH t/ −iθnσx −iEntσz / iθnσx Un,q = e n,q ~ = e e ~e . (3.15) bility of the time evolution on hLzi and hSzi depends on 5

-1.38 V. NON-RELATIVISTIC LIMIT -1.40 Since typically ξ[n] mc2 it is worth checking the non -1.42 relativistic limit of the Dirac oscillator, and for that we < J z >t -1.44 may use the quasi degenerate perturbation theory [18]. In this spirit and given that we know the exact solutions -1.46 for energies (3.7) and states (3.8), the non relativistic limit can be achieved from a series expansion in powers -1.48 of ξ. For the ﬁrst order on ξ of the eigenvalues we get:

10 20 30 40 50 60 70 2 τ En ' mc (1 + 2ξ[n]), (5.1) while the order ξ0 of the eigenvectors is

FIG. 2: Plot of hJzit in Eq. (4.10) as a function of the dimen- 2 sionless rescaled time τ = mc t/~ for α = 1 and q = 0.75. | + Eni = |ni|χ↑i, | − Eni = |n − 1i|χ↓i. (5.2)

D Thus, within this level of approximation we have Hq ' D q explicitly as well as implicitly thorough the factor [n]. Heff,q with Since q ≤ 1 and [n] ≤ n we get that the visibility for q 6= 1 is lesser than in the undeformed case q = 1. So " 2 † # D mc + 2~ωa`,qa`,q 0 q-deformation diminishes the effect. Heff,q = 2 † . 0 −mc − 2~ωa`,qa ii) At difference with the undeformed case q = 1, in `,q (5.3) the q-deformed scenario J is not constant of the motion z This result is quite natural since from an optical per- except for n = 1, as clearly shown in Eqs. (4.6) and (4.7). spective the condition ξ[n] mc2 actually means a very This is because after the commutation relation (2.1) we large detuning. In such a case there is no real exchange have [HD,J ] 6= 0. q z of photons or phonons and there is just a shift of levels iii) We can recover the constancy of hJzi provided that caused by virtual exchanges. Alternatively, from the in- we consider an alternative expression for Lz, this is that terferometric perspective in this limit the beam splitters Lz ∝ −~nˆ, wheren ˆ is defined in Eq. (2.5). Now the D exp(±iθσx) are replaced by the identity and the phases commutation relation (2.5) grants that [Hq ,Jz] = 0 and are approximate linearly in ξ. Jz is a constant of the motion as in the undeformed case. As a further example of Zitterbewegung we continue by considering a coherent initial state (2.12) so that |ψ i = 0 A. Zitterbewegung in the non-relativistic limit |αi|χ↓i leading to

∞ Lets peruse the existence of Zitterbewegung in the non- X |α|2n hL i = − |α|2 − ~ qnS (t), (4.8) relativistic limit. We can examine it in several steps. z t ~ e (|α|2) [n]! n+1 q n=0 On the one hand, after the form (5.3) for the eﬀective Hamiltonian we get that in this non-relativistic limit

∞ 2n D D X |α| [H ,Sz] = 0, [H ,Lz] = 0, (5.4) hS i = −~ + ~ S (t), (4.9) eﬀ,q eﬀ,q z t 2 e (|α|2) [n]! n+1 q n=0 so that Sz and Lz are constants of the motion and display no Zitterbewegung in this strict limit.

∞ 2n As a further alternative we may examine the evolu- ~ 2 1 − q X |α| tion of more sophisticated observables such as say M = hJzit = − 1 + 2|α| + ~ 2 [n]Sn+1(t), 2 eq(|α| ) [n]! † n=0 σ−a`,q + σ+a`,q. Since the key point of this eﬀect is the (4.10) interference between negative and positive energies, after where Eq. (5.2) we have that the initial state must be diﬀerent from (4.1). Thus let us consider instead for example the 4ξ[n] S (t) = sin2(ω t). (4.11) following initial state: n 1 + 4ξ[n] n |ψ0i = c↑|ni|χ↑i + c↓|n − 1i|χ↓i, (5.5) All above quantities satisfy the proper q = 1 limit, and in particular we have that hJzit is constant for q = 1. where c↑,↓ are real constants. It‘s time evolution results On the other hand, for q 6= 1 these expressions have the in: typical structure leading to an scenario of collapse and −iωnt iωnt revivals [20]. This is clearly shown in Fig. 2. |ψti = c↑e |ni|χ↑i + c↓e |n − 1i|χ↓i, (5.6) 6

2 where here ωn = mc (1 + 2ξ[n])/~. So we may easily have highlined the algebraic modifications introduced by calculate hMi leading to the q-deformation. In particular we have shown the lack of preservation of the q-deformed versions of the total p hMit = 2c↑c↓ [n] cos (2ωnt) , (5.7) angular momentum. Actually we have shown it leads to a time evolution mimicking the well-known collapse- displaying the desired oscillating fully analogous Eqs. revival structure of the population inversion or spin in (4.4) to (4.6). standard undeformed Jaynes-Cummings model. More- On the other hand we may consider also further correc- over, we have provided suitable relations for the non rel- tions to the non-relativistic limit. For example, regard- ativistic limit. ing Zitterbewegung we may consider the same initial state (4.1) including also the first approximation in powers of ξ in exp(±iθσx), leading to Acknowledgments hJ i = −~ (1 + 2[n − 1]) + 4ξ [n](1 − qn−1) sin2 (ω t) , A. L. acknowledges financial support from Spanish z t 2 ~ n (5.8) Ministerio de Econom´ıa y Competitividad Project No. for the same ω above. FIS2016-75199-P, and from the Comunidad Autónoma n de Madrid research consortium QUITEMAD+ Grant No. S2013/ICE-2801. VI. CONCLUSIONS

We have shown that the equivalence between the Appendix A: Differential forms two-dimensional Dirac oscillator and the anti-Jaynes- Cummings model extends to a q-deformed scenario. Af- Here we provide explicit differential expressions for the ter this fundamental equivalence we have investigated Hamiltonians of the q-deformed Dirac oscillator and the from first principles the Hamiltonian spectrum and time Jaynes-Cummings and anti-Jaynes Cummings models. evolution for initial number and coherent states, focus- For the Dirac oscillator, we have after Eqs. (2.6) and ing on the appearance of the Zitterbewegung effect. We (2.20):

d d HD = mc2σ + G[σ [(2 cos(2xα) + 2 sin(2yα) − cos(xα + α ) − i sin(xα + α ) q z x dx dx d d d d + i cos(yα + α ) − sin(yα + α ) − cos(α − xα) − i sin(α − xα) dy dy dx dx d d d − i cos(α − yα) + sin(α − yα)] + σ [−2 sin(2xα) + 2 cos(2yα) − i cos(xα + α ) dy dy y dx d d d d d − cos(yα + α ) + sin(xα + α ) − i sin(yα + α ) + i cos(α − xα) − sin(α − xα) dy dx dy dx dx d d − cos(α − yα) − i sin(α − yα)]], (A1) dy dy where G satisﬁes the relation √ mω G = c ~ . (A2) p2(1 − e−2α2 ) For the JC Hamiltonian (2.27): η0 d d HJC = √ ~ [σ [−2 sin(φ0 − 2r α) − i cos(φ0 − r α + α ) + sin(φ0 − r α + α ) q,i 2 x i i i 2 1 − e−2α dri dri 0 d 0 d 0 +i cos(φ − α − riα) + sin(φ − α − riα)]+ ([σy[−2 cos(φ − 2αri) dri dri 0 d 0 d 0 d −i sin(φ − α − riα) + i sin(φ + α − riα) + cos(φ − α − riα) dri dri dri 0 d 0 +cos(φ + α − riα )])] + (~δ σz) , (A3) dri 7 where ri=x,y = x, y, while for the A JC Hamiltonian (2.26):

η d d HAJC = √ ~ [σ (2 sin(φ + 2r α) + i cos(φ + r α + α ) − sin(φ + r α + α ) q,i 2 x i i i 2 1 − e−2α dri dri d d d − i cos(φ + α + riα) − sin(φ − α + riα)]+[σy(2 cos(φ + 2αri)−i sin(φ + α + riα) dri dri dri d d d − cos(φ + α + riα) − cos(φ − α + riα) + +i sin(φ − α + riα]] + ~δσz. (A4) dri dri dri

[1] D. Itô,K. Mori and E. Carrieri, An example of dynami- oscillator, coherent states, and the uncertainty relation, cal systems with linear trajectory, Nuovo Cimento 51A, Theor. Math. Phys. 147, 709 (2006). 1119 (1967). [14] N. Hatami and M. R. Setare, The q-deformed Dirac os- [2] M. Moshinsky and A. Szczepaniak, The Dirac oscillator, cillator in 2 + 1 dimensions, Phys. Lett. A 380, 3469 J. Phys. A 22, L817 (1989). (2016). [3] A. Bermudez, M. A. Martin-Delgado, and E. Solano, Ex- [15] V. V. Dodonov, ‘Nonclassical’ states in quantum optics: act mapping of the 2+1 Dirac oscillator onto the Jaynes- a ‘squeezed’ review of the first 75 years, J. Opt. B: Quan- Cummings model: Ion-trap experimental proposal, Phys. tum Semiclass. Opt. 4, R1 (2002). Rev. A 76, 041801(R) (2007). [16] W. Greiner, Relativistic Quantum Mechanics: Wave [4] A. Bermudez, M.A Martin-Delgado, and A. Luis, Non- Equations (Springer, Berlin, 2000). relativistic limit in the 2+1 Dirac oscillator: A Ramsey- [17] D. I. Fivel, Quasi-coherent states and the spectral reso- interferometry effect, Phys. Rev. A 77, 033832 (2008). lution of the q-Bose field operator, J. Phys. A 24, 3575 [5] J. A. Franco-Villafañe, E. Sadurn´ı, S. Barkhofen, U. (1991); V. Spiridonov, Coherent states of the q-Weyl al- Kuhl, F. Mortessagne, and T. H. Seligman, First Exper- gebra, Lett. Math. Phys. 35, 179 (1995). imental Realization of the Dirac Oscillator, Phys. Rev. [18] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, Lett. 111, 170405 (2013). Atom-Photon Interactions. Basic Processes and Applica- [6] S. Longhi, Photonic realization of the relativistic Dirac tions, (Wiley-VCH, Weinheim, 2004). oscillator, Opt. Lett. 35, 1302 (2010). [19] E. A. Kochetov, Exactly solvable non-linear generalisa- [7] Y. Luo, Y. Cui, Z. Long, and J. Jing, 2+1 Dimensional tions of the Jaynes-Cummings model, J. Phys. A: Math. Noncommutative Dirac Oscillator and (Anti)-Jaynes- Gen. 20, 2433 (1987); J. Crnugelj,ˇ M. Martinis, and Cummings Models, Int. J. Theor. Phys. 50, 2992 (2011). V. Mikuta-Martinis, Properties of a deformed Jaynes- [8] A. Jellal, A. El Mouhafid and M. Daoud, Massless Dirac Cummings model, Phys. Rev. A 50, 1785 (1994). fermions in an electromagnetic field, J. Stat. Mech. 2012, [20] S. Cordero and J. Récamier, Selective transition and P01021 (2012). complete revivals of a single two-level atom in the Jaynes- [9] E. T. Jaynes and W. F. Cummings, Comparison of quan- Cummings Hamiltonian with an additional Kerr medium, tum and semiclassical radiation theories with application J. Phys. B: At. Mol. Opt. Phys. 44, 135502 (2011) ; O. de to the beam maser, Proc. IEEE 51, 89 (1963); B. W. los Santos-Sánchez and J. Récamier,The f-deformed Jay- Shore and P. L. Knight , The Jaynes-Cummings Model, nesCummings model and its nonlinear coherent states, J. J. Mod. Opt. 40, 1195 (1993). Phys. B: At. Mol. Opt. Phys. 45, 015502 (2012). [10] D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, [21] W. Vogel and R. L. de Matos Filho, Nonlinear Jaynes- Quantum dynamics of single trapped ions, Rev. Mod. Cummings dynamics of a trapped ion, Phys. Rev. A 52, Phys. 75, 281(2003). 4214 (1995); R. L. de Matos Filho and W. Vogel, Non- [11] A. J. Macfarlane, On q-analogues of the quantum har- linear coherent states, Phys. Rev. A 54, 4560 (1996). monic oscillator and the quantum group SU(2)q, J. Phys. [22] R. Lynch, The quantum phase problem: a critical review, A 22, 4581 (1989); L. S. Biedenharn, The quantum group Phys. Rep. 256, 367 (1995); V. Peˇrinová,A. Lukˇsand SUq(2) and a q-analogue of the boson operators, J. Phys. J. Peˇrina, Phase in Optics (World Scientific, Singapore, A 22, L873 (1989); P. P. Kulish and E. V. Damaskinsky, 1998); P. Carruthers and M. M. Nieto, Phase and angle On the q oscillator and the quantum algebra suq(1, 1), J. variables in quantum mechanics, Rev. Mod. Phys. 40, Phys. A 23, L415 (1990). 411 (1968); D. T. Pegg and S. M. Barnett, Quantum [12] A. Lorek, A. Ruffing and J. Wess, A q-Deformation of optical phase, J. Mod. Opt. 44, 225 (1997). the Harmonic Oscillator, Z. Phys. C 74, 369 (1997). [13] V.V. Eremin, A.A. Meldianov, The q-deformed harmonic