Two-mode photon-added entangled coherent-squeezed states: their entanglement and nonclassical properties (original) (raw)
Two-mode photon-added entangled coherent-squeezed states: their entanglement and nonclassical properties
Amir Karimi
Applied Physics B
Lasers and Optics
ISSN 0946-2171
Volume 123
Number 6
Appl. Phys. B (2017) 123:1-9
DOI 10.1007/s00340-017-6757-0
Applied Physics B
Lasers
and Optics
(2) Springer
Your article is protected by copyright and all rights are held exclusively by SpringerVerlag Berlin Heidelberg. This e-offprint is for personal use only and shall not be selfarchived in electronic repositories. If you wish to self-archive your article, please use the accepted manuscript version for posting on your own website. You may further deposit the accepted manuscript version in any repository, provided it is only made publicly available 12 months after official publication or later and provided acknowledgement is given to the original source of publication and a link is inserted to the published article on Springer’s website. The link must be accompanied by the following text: “The final publication is available at link.springer.com”.
Two-mode photon-added entangled coherent-squeezed states: their entanglement and nonclassical properties
Amir Karimi 1{ }^{1}
Abstract
Recently, we introduced and generated new types of two-mode entangled states named “entangled coherent-squeezed states”. In these states, two common states for quantum information processing, coherent state, and squeezed states have been used. Now, based on the generated entangled states, we introduce “two-mode photon-added entangled coherent-squeezed states”. These states are obtained via the iterated action of two creation operators on the two modes of the “entangled coherentsqueezed states”. Next, we study the amount of entanglement of the introduced states using concurrence criterion. In the continuation, some of the nonclassical features such as photon-statistics, second-order correlation function, and quadrature squeezing are considered. In addition, we study the influence of photon-addition of two modes on the mentioned properties of the introduced states. We will observe that the entanglement of the introduced photon-added entangled states increases more rapidly as photon-addition of two modes increases. Moreover, some of the nonclassical features for the first mode of the introduced states such as sub-Poissonian photon-statistics and squeezing in pp appear and disappear by photon-addition of two modes, respectively.
1 Introduction
Quantum entanglement can be mentioned as the main concept and most important basis in quantum information
[1]science such as cryptography [1], quantum teleportation [2], quantum computation [3], and so on. This notion (entanglement or entangled states) which is related to special superposition of several different states was first used by Schrödinger in 1935. He called the wave function of two particles as an entangled state if it cannot be written as the product of the wave functions of two particles [4]. In another word, it can be said that the concept of entanglement is referred to the inseparability feature of a multi-mode system to multiplication of its single-mode subsystems.
In the last two decades, various types of entangled quantum states have been proposed and studied which, for example, can be mentioned entangled coherent states [5-7], entangled squeezed states [8-11], and entangled cat states [12, 13]. Soon, after the introduction of these entangled states, the corresponding photon-added and photon-subtracted entangled states attracted a great deal of attention and some of their nonclassical properties have been considered [14-24]. Moreover, many studies have been done to generate such entangled states theoretically and experimentally [25-33].
Recently, we introduced and generated new types of entangled states namely “entangled coherent-squeezed states”, and studied some of their nonclassical properties [34]. It is worth noting that to generate the entangled states, we used a theoretical scheme based on the interaction between a Λ\Lambda-type three-level atom and a two-mode quantized field under certain conditions (For more information, see [34]). Now, the increasing importance to achieve new entangled states with different nonclassical features motivated us to introduce “two-mode photon-added entangled coherent-squeezed states” and study some of their nonclassicalities.
- ㅁ Amir Karimi amirkarimi.phy@gmail.com; akarimi@iauabadeh.ac.ir
1 Physics Department, Abadeh Branch, Islamic Azad University, Fars, Iran ↩︎
The remainder of paper organizes as follows: in the next section, we have a brief review on the “entangled coherent-squeezed states”. Next, by making use of these states, we introduce the “two-mode photon-added entangled coherent-squeezed states” in Sect. 3. In Sect. 4, we study the amount of entanglement of the introduced pho-ton-added entangled states by evaluating concurrence. In the continuation, we consider some of the nonclassical features of the introduced photon-added entangled states such as quantum statistics, second-order correlation function, and quadrature squeezing, as well as the influence of two-mode photon-addition on the mentioned nonclassical properties in Sect. 5. Finally, we end our paper by presenting a summary in Sect. 6.
2 Entangled coherent-squeezed states: a brief review
The special form of the “entangled coherent-squeezed states” has recently been introduced and generated in the following forms [34]:
∣ψ(α,ξ)⟩i=Mi(α,ξ)(∣α⟩⊗∣−ξ⟩i+∣−α⟩⊗∣ξ⟩ii=e,o\begin{gathered} |\psi(\alpha, \xi)\rangle_{i}=M_{i}(\alpha, \xi)(|\alpha\rangle \otimes|-\xi\rangle_{i}+|-\alpha\rangle \otimes|\xi\rangle_{i} \\ i=\mathrm{e}, \mathrm{o} \end{gathered}
where Me(α,ξ)M_{\mathrm{e}}(\alpha, \xi) and Mo(α,ξ)M_{\mathrm{o}}(\alpha, \xi) are normalization constants. The states ∣α⟩|\alpha\rangle and ∣ξ⟩i(i=e,o)|\xi\rangle_{i}(i=\mathrm{e}, \mathrm{o}) are, respectively, coherent state and squeezed states that will be explained in what follows.
The coherent state ∣α⟩|\alpha\rangle that known as the minimum uncertainty state can be generated by displacing the vacuum state ∣α⟩=D(α)∣0⟩|\alpha\rangle=D(\alpha)|0\rangle where D(α)=eαa†−α∗aD(\alpha)=e^{\alpha a^{\dagger}-\alpha^{*} a} is displacement operator ( aa and a†a^{\dagger} are the bosonic annihilation and creation operators, respectively, and α=reiθ\alpha=r e^{i \theta} is a complex number). Moreover, the coherent state can be obtained by solving the eigen-value equation a∣α⟩=α∣α⟩a|\alpha\rangle=\alpha|\alpha\rangle as [35]:
∣α⟩=e−∣α∣22∑p=0∞αpp!∣p⟩|\alpha\rangle=e^{-\frac{|\alpha|^{2}}{2}} \sum_{p=0}^{\infty} \frac{\alpha^{p}}{\sqrt{p!}}|p\rangle.
Two types of squeezed states ∣ξ⟩e|\xi\rangle_{\mathrm{e}} and ∣ξ⟩o|\xi\rangle_{\mathrm{o}} that for which one of the quadrature fluctuation becomes less than that of the vacuum state can be, respectively, obtained by squeezing the vacuum state ∣ξ⟩e=S(ξ)∣0⟩|\xi\rangle_{\mathrm{e}}=S(\xi)|0\rangle and one-photon state of the field ∣ξ⟩o=S(ξ)∣1⟩|\xi\rangle_{\mathrm{o}}=S(\xi)|1\rangle in the following forms [36]:
∣ξ⟩e=Ne(ξ)∑q=0∞(−1)q(2q)!2qq!(eiθtanhr)q∣2q⟩|\xi\rangle_{\mathrm{e}}=N_{\mathrm{e}}(\xi) \sum_{q=0}^{\infty}(-1)^{q} \frac{\sqrt{(2 q)!}}{2^{q} q!}\left(e^{i \theta} \tanh r\right)^{q}|2 q\rangle,
∣ξ⟩o=No(ξ)∑q=0∞(−1)q2qq!(2q+1)!(eiθtanhr)q∣2q+1⟩|\xi\rangle_{\mathrm{o}}=N_{\mathrm{o}}(\xi) \sum_{q=0}^{\infty}(-1)^{q} \frac{2^{q} q!}{\sqrt{(2 q+1)!}}\left(e^{i \theta} \tanh r\right)^{q}|2 q+1\rangle,
where S(ξ)=eξ22a2−ξ2a2S(\xi)=e^{\frac{\xi^{2}}{2} a^{2}-\frac{\xi}{2} a^{2}} is squeezing operator. The operators aa and a†a^{\dagger} have the same definitions of coherent state and ξ=reiθ\xi=r e^{i \theta} is a complex parameter, too.
Due to the fact that only the even and odd number states appear, respectively, in Eqs. (3) and (4), these squeezed states are, respectively, called even and odd squeezed states and labeled by e and o.
The normalization coefficients Ne(ξ)N_{\mathrm{e}}(\xi) and No(ξ)N_{\mathrm{o}}(\xi) can be easily obtained as:
Ne−2(ξ)=∑q=0∞(2q)!22qq!2(tanhr)2qN_{\mathrm{e}}^{-2}(\xi)=\sum_{q=0}^{\infty} \frac{(2 q)!}{2^{2 q} q!^{2}}(\tanh r)^{2 q},
No−2(ξ)=∑q=0∞23qq!2(2q+1)!(tanhr)2qN_{\mathrm{o}}^{-2}(\xi)=\sum_{q=0}^{\infty} \frac{2^{3 q} q!^{2}}{(2 q+1)!}(\tanh r)^{2 q}.
It is worth noting that the even (odd) squeezed states ∣ξ⟩e|\xi\rangle_{\mathrm{e}} and ∣−ξ⟩e(∣ξ⟩o|-\xi\rangle_{\mathrm{e}}\left(|\xi\rangle_{\mathrm{o}}\right. and ∣−ξ⟩o|-\xi\rangle_{\mathrm{o}} ) in Eq. (1) can be obtained by setting θ=0\theta=0 and θ=π\theta=\pi in Eqs. (3) and (4), respectively. In the same way, the coherent states ∣±α⟩| \pm \alpha\rangle are obtained from Eq. (2), too.
Now, using the coherent state in Eq. (2) and squeezed states in Eqs. (3) and (4), we can derive the number bases representation of the introduced “entangled coherentsqueezed states” in Eq. (1) in the following forms:
∣ψ(α,ξ)⟩e=Me(α,ξ)Ne(ξ)e−τ22×∑p=0∞∑q=0∞((−1)p+(−1)q)×rpp!(2q)!2qq!(tanhr)q∣p⟩∣2q⟩∣ψ(α,ξ)⟩o=Mo(α,ξ)No(ξ)e−τ22×∑p=0∞∑q=0∞((−1)p+(−1)q)×rpp!2qq!(2q+1)!(tanhr)q∣p⟩∣2q+1⟩\begin{aligned} |\psi(\alpha, \xi)\rangle_{\mathrm{e}}= & M_{\mathrm{e}}(\alpha, \xi) N_{\mathrm{e}}(\xi) e^{-\frac{\tau^{2}}{2}} \times \sum_{p=0}^{\infty} \sum_{q=0}^{\infty}\left((-1)^{p}+(-1)^{q}\right) \\ & \times \frac{r^{p}}{\sqrt{p!}} \frac{\sqrt{(2 q)!}}{2^{q} q!}(\tanh r)^{q}|p\rangle|2 q\rangle \\ |\psi(\alpha, \xi)\rangle_{\mathrm{o}}= & M_{\mathrm{o}}(\alpha, \xi) N_{\mathrm{o}}(\xi) e^{-\frac{\tau^{2}}{2}} \times \sum_{p=0}^{\infty} \sum_{q=0}^{\infty}\left((-1)^{p}+(-1)^{q}\right) \\ & \times \frac{r^{p}}{\sqrt{p!}} \frac{2^{q} q!}{\sqrt{(2 q+1)!}}(\tanh r)^{q}|p\rangle|2 q+1\rangle \end{aligned}
The normalization factors Me(α,ξ)M_{\mathrm{e}}(\alpha, \xi) and Mo(α,ξ)M_{\mathrm{o}}(\alpha, \xi) are given by:
Me−2(α,ξ)=Ne2(ξ)e−r2∑p=0∞∑q=0∞(2+2(−1)p+q)×r2pp!(2q)!22qq!2(tanhr)2q\begin{aligned} M_{\mathrm{e}}^{-2}(\alpha, \xi)= & N_{\mathrm{e}}^{2}(\xi) e^{-r^{2}} \sum_{p=0}^{\infty} \sum_{q=0}^{\infty}\left(2+2(-1)^{p+q}\right) \\ & \times \frac{r^{2 p}}{p!} \frac{(2 q)!}{2^{2 q} q!^{2}}(\tanh r)^{2 q} \end{aligned}
Mo−2(α,ξ)=No2(ξ)e−r2∑p=0∞∑q=0∞(2+2(−1)p+q)×r2pp!22qq!2(2q+1)!(tanhr)2q\begin{aligned} M_{\mathrm{o}}^{-2}(\alpha, \xi)= & N_{\mathrm{o}}^{2}(\xi) e^{-r^{2}} \sum_{p=0}^{\infty} \sum_{q=0}^{\infty}\left(2+2(-1)^{p+q}\right) \\ & \times \frac{r^{2 p}}{p!} \frac{2^{2 q} q!^{2}}{(2 q+1)!}(\tanh r)^{2 q} \end{aligned}
where Ne(ξ)N_{\mathrm{e}}(\xi) and No(ξ)N_{\mathrm{o}}(\xi) have been obtained in Eqs. (5) and (6), respectively.
It is worth mentioning that here, ∣α∣|\alpha| and ∣ξ∣|\xi| have been supposed the same and equal to rr.
3 Two-mode photon-added entangled coherent-squeezed states
Now, we introduce the “two-mode photon-added entangled coherent-squeezed states”. These states can be obtained via the iterated actions of two creation operators a†a^{\dagger} and b†b^{\dagger}, respectively, on the first and second modes of the “entangled coherent-squeezed states” in Eq. (1) in the following forms:
∣ψ(α,ξ,m,n)⟩i=Mi(α,ξ,m,n)(a†)m(b†)n(∣α⟩⊗∣−ξ⟩i+∣−α⟩⊗∣ξ⟩i),i=e,o\begin{aligned} & |\psi(\alpha, \xi, m, n)\rangle_{i}=M_{i}(\alpha, \xi, m, n)\left(a^{\dagger}\right)^{m}\left(b^{\dagger}\right)^{n} \\ & \left(|\alpha\rangle \otimes|-\xi\rangle_{i}+|-\alpha\rangle \otimes|\xi\rangle_{i}\right), \quad i=\mathrm{e}, \mathrm{o} \end{aligned}
Using the even and odd “entangled coherent-squeezed states” in Eqs. (7) and (8), we can get the number basis representations of the introduced photon-added entangled states in Eq. (11) as:
∣ψ(α,ξ,m,n)⟩e=Me(α,ξ,m,n)Ne(ξ)e−r22×∑p=0∞∑q=0∞((−1)p+(−1)q)(p+m)!p!(2q+n)!2qq!×rp(tanhr)q∣p+m⟩∣2q+n⟩∣ψ(α,ξ,m,n)⟩o=Mo(α,ξ,m,n)No(ξ)e−r22×∑p=0∞∑q=0∞((−1)p+(−1)q)(p+m)!p!2qq!(2q+n+1)!(2q+n+1)!(2q+1)!×rp(tanhr)q∣p+m⟩∣2q+n+1⟩\begin{aligned} & |\psi(\alpha, \xi, m, n)\rangle_{\mathrm{e}}=M_{\mathrm{e}}(\alpha, \xi, m, n) N_{\mathrm{e}}(\xi) e^{-\frac{r^{2}}{2}} \\ & \quad \times \sum_{p=0}^{\infty} \sum_{q=0}^{\infty}\left((-1)^{p}+(-1)^{q}\right) \frac{\sqrt{(p+m)!}}{p!} \frac{\sqrt{(2 q+n)!}}{2^{q} q!} \\ & \quad \times r^{p}(\tanh r)^{q}|p+m\rangle|2 q+n\rangle \\ & |\psi(\alpha, \xi, m, n)\rangle_{\mathrm{o}}=M_{\mathrm{o}}(\alpha, \xi, m, n) N_{\mathrm{o}}(\xi) e^{-\frac{r^{2}}{2}} \\ & \quad \times \sum_{p=0}^{\infty} \sum_{q=0}^{\infty}\left((-1)^{p}+(-1)^{q}\right) \frac{\sqrt{(p+m)!}}{p!} \frac{2^{q} q!}{\sqrt{(2 q+n+1)!}} \frac{\sqrt{(2 q+n+1)!}}{(2 q+1)!} \\ & \quad \times r^{p}(\tanh r)^{q}|p+m\rangle|2 q+n+1\rangle \end{aligned}
where the normalization constants Me(α,ξ,m,n)M_{\mathrm{e}}(\alpha, \xi, m, n) and Mo(α,ξ,m,n)M_{\mathrm{o}}(\alpha, \xi, m, n) are easily obtained as:
Me−2(α,ξ,m,n)=Ne2(ξ)e−r2∑p=0∞∑q=0∞(2+2(−1)p+q)×(p+m)!p!2(2q+n)!22qq!2r2p(tanhr)2q\begin{aligned} M_{\mathrm{e}}^{-2}(\alpha, \xi, m, n)= & N_{\mathrm{e}}^{2}(\xi) e^{-r^{2}} \sum_{p=0}^{\infty} \sum_{q=0}^{\infty}\left(2+2(-1)^{p+q}\right) \\ & \times \frac{(p+m)!}{p!^{2}} \frac{(2 q+n)!}{2^{2 q} q!^{2}} r^{2 p}(\tanh r)^{2 q} \end{aligned}
Mo−2(α,ξ,m,n)=No2(ξ)e−r2×∑p=0∞∑q=0∞(2+2(−1)p+q)(p+m)!p!222qq!2(2q+n+1)!(2q+1)!2r2p(tanhr)2q\begin{gathered} M_{\mathrm{o}}^{-2}(\alpha, \xi, m, n)=N_{\mathrm{o}}^{2}(\xi) e^{-r^{2}} \times \sum_{p=0}^{\infty} \sum_{q=0}^{\infty}\left(2+2(-1)^{p+q}\right) \\ \frac{(p+m)!}{p!^{2}} \frac{2^{2 q} q!^{2}(2 q+n+1)!}{(2 q+1)!^{2}} r^{2 p}(\tanh r)^{2 q} \end{gathered}
In these relations, the coefficients Ne(ξ)N_{\mathrm{e}}(\xi) and No(ξ)N_{\mathrm{o}}(\xi) have been, respectively, determined in Eqs. (5) and (6).
4 Entanglement of the introduced states
In this section, we want to study the amount of entanglement of the introduced “two-mode photon-added entangled coherent-squeezed states” by evaluating “concurrence”.
This criterion that known as one of the appropriate measures to study the entanglement degree of the introduced states is defined for discrete-variable entangled states in terms of the Pauli matrix σy\sigma_{y} as [38]:
C=∣⟨Ψ∣σy⊗σy∣Ψ∗⟩∣C=\left|\left\langle\Psi \mid \sigma_{y} \otimes \sigma_{y} \mid \Psi^{*}\right\rangle\right|,
where ∣Ψ∗⟩\left|\Psi^{*}\right\rangle stands for the complex conjugate of ∣Ψ⟩|\Psi\rangle. Notice that the maximum value of the concurrence is 1 for a maximum entangled state.
To calculate the concurrence of the introduced entangled states, we consider a two-mode entangled state in the following general form:
∣Ψ⟩=1N(μ∣η⟩⊗∣γ⟩+v∣ξ⟩⊗∣δ⟩)|\Psi\rangle=\frac{1}{N}(\mu|\eta\rangle \otimes|\gamma\rangle+v|\xi\rangle \otimes|\delta\rangle),
where μ\mu and vv are complex numbers. The first (second) mode of such entangled state has been shown with ∣η⟩|\eta\rangle and ∣ξ⟩⟨∣γ⟩|\xi\rangle\langle | \gamma\rangle and ∣δ⟩)|\delta\rangle).
For such state, using the transformation of the continu-ous-variables-type components to the discrete orthogonal basis and the Schmidt decomposition, the concurrence of the entangled state (17) is determined by [39-42]:
C=2∣μ∣∣v∣N2(1−∣P(1)∣2)(1−∣P(2)∣2)C=\frac{2|\mu||v|}{N^{2}} \sqrt{\left(1-\left|P^{(1)}\right|^{2}\right)\left(1-\left|P^{(2)}\right|^{2}\right)}, P(1)=⟨η∣ξ⟩,P(2)=⟨δ∣γ⟩P^{(1)}=\langle\eta \mid \xi\rangle, P^{(2)}=\langle\delta \mid \gamma\rangle.
Now, to calculate the concurrence of the “two-mode photonadded entangled coherent-squeezed states” in Eq. (11), we first rewrite the introduced “two-mode photon-added entangled coherent-squeezed states” in Eq. (11) in the explicit form. For this purpose, we use the photon-added coherent states introduced first by Agarwal and Tara [37] as:
∣±α,m⟩=d(α,m)(a†)m∣±α⟩=d(α,m)e−∣α∣22×∑p=0∞(p+m)!p!(±α)p∣p+m⟩\begin{aligned} | \pm \alpha, m\rangle= & d(\alpha, m)\left(a^{\dagger}\right)^{m}| \pm \alpha\rangle=d(\alpha, m) e^{-\frac{|\alpha|^{2}}{2}} \\ & \times \sum_{p=0}^{\infty} \frac{\sqrt{(p+m)!}}{p!}( \pm \alpha)^{p}|p+m\rangle \end{aligned}
where the coefficient d(α,m)d(\alpha, m) is given by:
d−2(α,m)=e−∣α∣2∑p=0∞(p+m)!p!2∣α∣2pd^{-2}(\alpha, m)=e^{-|\alpha|^{2}} \sum_{p=0}^{\infty} \frac{(p+m)!}{p!^{2}}|\alpha|^{2 p}.
In the similar way, we can define the even and odd pho-ton-added squeezed states, respectively, in the following forms:
±ξ,n⟩e=de(ξ,n)(b†)n∣±ξ⟩e=de(ξ,n)Ne(ξ)∑q=0∞(∓1)q(2q+n)!2qq!×(tanhr)q∣2q+n⟩,±ξ,n⟩o=do(ξ,n)(b†)n∣±ξ⟩o=do(ξ,n)No(ξ)∑q=0∞(∓1)q2qq!(2q+n+1)!(2q+1)!×(tanhr)q∣2q+n+1⟩.\begin{aligned} & \left. \pm \xi, n\right\rangle_{\mathrm{e}}=d_{\mathrm{e}}(\xi, n)\left(b^{\dagger}\right)^{n}| \pm \xi\rangle_{\mathrm{e}} \\ & =d_{\mathrm{e}}(\xi, n) N_{\mathrm{e}}(\xi) \sum_{q=0}^{\infty}(\mp 1)^{q} \frac{\sqrt{(2 q+n)!}}{2^{q} q!} \\ & \times(\tanh r)^{q}|2 q+n\rangle, \\ & \left. \pm \xi, n\right\rangle_{\mathrm{o}}= d_{\mathrm{o}}(\xi, n)\left(b^{\dagger}\right)^{n}| \pm \xi\rangle_{\mathrm{o}} \\ & =d_{\mathrm{o}}(\xi, n) N_{\mathrm{o}}(\xi) \sum_{q=0}^{\infty}(\mp 1)^{q} \frac{2^{q} q!\sqrt{(2 q+n+1)!}}{(2 q+1)!} \\ & \times(\tanh r)^{q}|2 q+n+1\rangle . \end{aligned}
The coefficients de(ξ,n)d_{\mathrm{e}}(\xi, n) and do(ξ,n)d_{\mathrm{o}}(\xi, n) are determined by the normalization condition as:
de−2(ξ,n)=Ne2(ξ)∑q=0∞(2q+n)!22qq!2(tanhr)2qd_{\mathrm{e}}^{-2}(\xi, n)=N_{\mathrm{e}}^{2}(\xi) \sum_{q=0}^{\infty} \frac{(2 q+n)!}{2^{2 q} q!^{2}}(\tanh r)^{2 q},
do−2(ξ,n)=No2(ξ)∑q=0∞22qq!2(2q+n+1)!(2q+1)!2(tanhr)2qd_{\mathrm{o}}^{-2}(\xi, n)=N_{\mathrm{o}}^{2}(\xi) \sum_{q=0}^{\infty} \frac{2^{2 q} q!^{2}(2 q+n+1)!}{(2 q+1)!^{2}}(\tanh r)^{2 q}.
By making use of the introduced photon-added coherent states in Eq. (19) and photon-added squeezed states in Eqs. (21) and (22), we can write the explicit form of the “two-mode photon-added entangled coherent-squeezed states” in Eq. (11) in the following forms:
∣ψ(α,ξ,m,n)⟩i=Mi(α,ξ,m,n)d−1(α,m)di−1(ξ,n)×(∣α,m⟩⊗∣−ξ,n⟩i+∣−α,m⟩⊗∣ξ,n⟩i),i=e,o\begin{aligned} & |\psi(\alpha, \xi, m, n)\rangle_{i}=M_{i}(\alpha, \xi, m, n) d^{-1}(\alpha, m) d_{i}^{-1}(\xi, n) \\ & \quad \times\left(|\alpha, m\rangle \otimes|-\xi, n\rangle_{i}+|-\alpha, m\rangle \otimes|\xi, n\rangle_{i}\right), \quad i=\mathrm{e}, \mathrm{o} \end{aligned}
Now, we can obtain the concurrence of the introduced “two-mode photon-added entangled coherent-squeezed states” in Eq. (11) by comparing Eq. (25) with Eq. (17) and using Eq. (18) as:
Ci(α,ξ,m,n)=2Mi2(α,ξ,m,n)d2(α,m)di2(ξ,n)×(1−∣P(1)(α,m)∣2)(1−∣Pi(2)(ξ,n)∣2),i=e,o\begin{aligned} C_{i}(\alpha, \xi, m, n)= & \frac{2 M_{i}^{2}(\alpha, \xi, m, n)}{d^{2}(\alpha, m) d_{i}^{2}(\xi, n)} \\ & \times \sqrt{\left(1-\left|P^{(1)}(\alpha, m)\right|^{2}\right)\left(1-\left|P_{i}^{(2)}(\xi, n)\right|^{2}\right)}, \\ & i=\mathrm{e}, \mathrm{o} \end{aligned}
Using the Eqs. (19) and (21), we can, respectively, calculate the quantities P(1)(α,m)P^{(1)}(\alpha, m) and Pe(2)(ξ,n)P_{\mathrm{e}}^{(2)}(\xi, n) of the even “two-mode photon-added entangled coherent-squeezed states” in Eq. (25) in the following forms:
P(1)(α,m)=⟨−α,m∣α,m⟩=d2(α,m)e−∣α∣2∑p=0∞(−1)p(p+m)!p!2∣α∣2pPe(2)(ξ,n)=e⟨−ξ,n∣ξ,n⟩e=de2(ξ,n)Ne2(ξ)∑q=0∞(−1)q(2q+n)!22qq!2(tanhr)2q\begin{aligned} P^{(1)}(\alpha, m) & =\langle-\alpha, m \mid \alpha, m\rangle \\ & =d^{2}(\alpha, m) e^{-|\alpha|^{2}} \sum_{p=0}^{\infty}(-1)^{p} \frac{(p+m)!}{p!^{2}}|\alpha|^{2 p} \\ P_{\mathrm{e}}^{(2)}(\xi, n) & =\mathrm{e}\langle-\xi, n \mid \xi, n\rangle_{\mathrm{e}} \\ & =d_{\mathrm{e}}^{2}(\xi, n) N_{\mathrm{e}}^{2}(\xi) \sum_{q=0}^{\infty}(-1)^{q} \frac{(2 q+n)!}{2^{2 q} q!^{2}}(\tanh r)^{2 q} \end{aligned}
For the odd “two-mode photon-added entangled coher-ent-squeezed states”, the quantity P(1)(α,m)P^{(1)}(\alpha, m) is the same with that of the even photon-added entangled states which has been obtained in Eq. (27). Similarly, using Eq. (22), the quantity Po(2)(ξ,n)P_{\mathrm{o}}^{(2)}(\xi, n) of the odd “two-mode pho-ton-added entangled coherent-squeezed states” can be easily obtained as:
Po(2)(ξ,n)=o⟨−ξ,n∣ξ,n⟩o=do2(ξ,n)No2(ξ)∑q=0∞(−1)q22qq!2(2q+n+1)!(2q+1)!2×(tanhr)2q+1\begin{aligned} P_{\mathrm{o}}^{(2)}(\xi, n)= & { }_{\mathrm{o}}\langle-\xi, n \mid \xi, n\rangle_{\mathrm{o}} \\ = & d_{\mathrm{o}}^{2}(\xi, n) N_{\mathrm{o}}^{2}(\xi) \sum_{q=0}^{\infty}(-1)^{q} \frac{2^{2 q} q!^{2}(2 q+n+1)!}{(2 q+1)!^{2}} \\ & \times(\tanh r)^{2 q+1} \end{aligned}
Now, we plot the concurrence of the introduced “twomode photon-added entangled coherent-squeezed states” in Eq. (11) versus r(∣α∣r(|\alpha| and ∣ξ∣)|\xi|) for the different values of photon-additions of two-mode, mm and nn, in Fig. 1.
It can be observed that the concurrence of the even and odd “two-mode photon-added entangled coherentsqueezed states” in Eq. (11) shows similar behaviours.
The concurrence of the introduced photon-added entangled states starts from zero, increases monotonically and tends to 1 in the region of r<2r<2, and remains fixed. Clearly, it can be found that by photon-addition and also increasing the number of photon-addition of two modes, mm and nn, the concurrence increases more rapidly. Moreover, it is worth noting that the photon-addition of the second mode (squeezed states ∣±ξ⟩e,o| \pm \xi\rangle_{\mathrm{e}, \mathrm{o}} ) of the introduced photon-added entangled states has more impact on the increasing rate of concurrence.
5 Other nonclassical properties of the introduced states
In this section, some of the nonclassical properties of the introduced “two-mode photon-added entangled coherentsqueezed states” in Eq. (11) are studied. For this purpose,
Author’s personal copy
Fig. 1 The concurrence C(r)C(r) of the a\mathbf{a} even and b\mathbf{b} odd “twomode photon-added entangled coherent-squeezed states” ∣ψ(α,ξ,m,n)⟩e,α|\psi(\alpha, \xi, m, n)\rangle_{e, \alpha} versus r⟨∣α∣r\langle | \alpha| and ∣ξ⟩|\xi\rangle for the different values of photon-additions m,n=0m, n=0 (solid), m,n=2m, n=2 (dashed), m=4,n=2m=4, n=2 (dash-dot), and m=2,n=4m=2, n=4 (dot)
we consider a few important nonclassical features such as “quantum statistics”, “second-order correlation function”, and “quadrature squeezing”. It should be noted that all necessary expectation values in this section will be obtained using the expanded forms of the introduced states in Eqs. (12) and (13).
5.1 Quantum statistics
The quantum statistics of any quantum state is usually investigated by studying the two criteria of that state as sub-Poissonian photon-statistics and oscillatory behaviour of photon distribution.
5.1.1 Sub-Poissonian photon-statistics
First, we study the sub-Poissonian photon-statistics of the introduced photon-added entangled states in Eq. (11). To achieve this goal, we calculate the Mandel parameter [43] for both modes which is defined as [43]:
Qi=⟨ni2⟩−⟨ni⟩2⟨ni⟩−1,i=a,bQ_{i}=\frac{\left\langle n_{i}^{2}\right\rangle-\left\langle n_{i}\right\rangle^{2}}{\left\langle n_{i}\right\rangle}-1, \quad i=a, b,
where na=a†an_{a}=a^{\dagger} a and nb=b†bn_{b}=b^{\dagger} b are the number operators corresponding to the first and second modes of the introduced entangled states, respectively.
The positivity and negativity of the Mandel parameter indicate, respectively, to super-Poissonian and subPoissonian photon-statistics. Moreover, this parameter becomes zero for states with Poissonian photon-statistics such as coherent state. It should be noted that one of the nonclassical features of each quantum state is specified by negative values of Mandel parameter [44].
By performing the necessary calculations, the Mandel parameter of the first and second modes of the even (odd) “two-mode photon-added entangled coherent-squeezed states” in Eq. (11) is plotted versus r⟨∣α∣r\langle | \alpha| and ∣ξ⟩|\xi\rangle for different values of photon-additions of two modes, mm and nn, in Fig. 2 (Fig. 3).
Clearly, the Mandel parameter of the first mode Qa(r)Q_{a}(r) of the introduced even and odd “two-mode photon-added entangled coherent-squeezed states” in Figs. 2a and 3a shows the same behaviours. It can be observed that the first mode of the even and odd “entangled coherentsqueezed states” (m,n=0)(m, n=0) has super-Poissonian photonstatistics over all of the regions of rr. While the first mode of the even and odd “two-mode photon-added entangled coherent-squeezed states” becomes sub-Poissonian for all chosen values of photon-additions (m,n≠0)(m, n \neq 0) over all of the regions of rr. Moreover, it can be found from Figs. 2a and 3a that the sub-Poissonian photon-statistics of the first mode of the introduced even and odd photon-added entangled states increases by increasing the photon-addition of the first mode (m)(m) of the “two-mode photon-added entangled coherent-squeezed states”.
In the case of the second mode of the introduced entangled states, the Mandel parameter Qb(r)Q_{b}(r) shows no sub-Poissonian photon-statistics for the even “entangled coherent-squeezed state” (no photon-additions m,n=0m, n=0 ) in Fig. 2b while for the odd “entangled coher-ent-squeezed state” (no photon-additions m,n=0m, n=0 ) subPoissonian photon-statistics is seen over the region of small r⟨∣α∣r\langle | \alpha| and ∣ξ⟩|\xi\rangle in Fig. 3b. By photon-addition of two modes (m,n≠0)(m, n \neq 0), the sub-Poissonian photon-statistics is observed for the second mode of both states (even and odd) “two-mode photon-added entangled coherentsqueezed states” in Figs. 2b and 3b over the region of very small r⟨∣α∣r\langle | \alpha| and ∣ξ⟩|\xi\rangle ).
It is worth noting that the region of sub-Poissonian photon-statistics of the second mode of the even “twomode photon-added entangled coherent-squeezed states” becomes more slightly wider with increasing the pho-ton-addition of the second mode ( nn ) of the introduced entangled states (Fig. 2b). While, for the second mode of the odd “two-mode photon-added entangled coherentsqueezed states”, the region of sub-Poissonian photonstatistics becomes slightly narrower with increasing nn (Fig. 3b).
Fig. 2 Mandel parameter Qp(r)Q_{p}(r) of the a first mode and b\mathbf{b} second mode of the even “two-mode photon-added entangled coherentsqueezed states” ∣ψ(α,ξ,m,n)⟩0|\psi(\alpha, \xi, m, n)\rangle_{0} versus r⟨∣α∣r\langle | \alpha \mid and ∣ξ⟩|\xi\rangle for the different values of photon-additions m,n=0m, n=0 (solid), m,n=2m, n=2 (dashed),
m=4,n=2m=4, n=2 (dash-dot), and m=2,n=4m=2, n=4 (dot) (the m,n=2m, n=2 and m=2,n=4m=2, n=4 curves in a and the m,n=2m, n=2 and m=4,n=2m=4, n=2 curves in b coincide with each other)
Fig. 3 Mandel parameter Qo(r)Q_{o}(r) of the a first mode and b\mathbf{b} second mode of the odd “twomode photon-added entangled coherent-squeezed states” ∣ψ(α,ξ,m,n)⟩0|\psi(\alpha, \xi, m, n)\rangle_{0} versus r⟨∣α∣r\langle | \alpha \mid and ∣ξ⟩|\xi\rangle for the different values of photon-additions mm and nn (all of the parameters and explanations are the same as Fig. 2)
5.1.2 Photon number distribution
As it is pointed out before, the photon number distribution is another criterion to study the quantum statistics which its oscillatory behaviour can show the nonclassical property of each quantum state.
The photon number distribution P(j,k)P(j, k) which is defined as the probability of finding jj and kk photons, respectively, in the first and second modes of the introduced “two-mode photon-added entangled coherentsqueezed states” in Eq. (11) is given by:
P1(j,k)=∣⟨j,k∣ψ(α,ξ,m,n)⟩1∣2,i=e,oP_{1}(j, k)=\left|\langle j, k \mid \psi(\alpha, \xi, m, n)\rangle_{1}\right|^{2}, \quad i=\mathrm{e}, \mathrm{o}.
By calculating Pi(j,k)P_{i}(j, k), we plotted the three-dimensional photon number distribution of the introduced even (odd) photon-added entangled states in Eq. (11) versus the photon numbers of the first mode jj and second mode kk in Fig. 4 (Fig. 5) by fixing r=1r=1.
Clearly, an oscillatory behaviour can be observed in the photon number distribution of the introduced
photon-added entangled states which indicates to the nonclassical feature of these states.
5.2 Two-mode correlation function
In this subsection, the second-order correlation function of the introduced “two-mode photon-added entangled coherent-squeezed states” in Eq. (11) is considered. This feature which has been defined as [45]:
gnb(3)(0)=⟨a†ab†b⟩⟨a†a⟩⟨b†b⟩g_{n b}^{(3)}(0)=\frac{\left\langle a^{\dagger} a b^{\dagger} b\right\rangle}{\left\langle a^{\dagger} a\right\rangle\left\langle b^{\dagger} b\right\rangle},
becomes less than 1 for a nonclassical state [46]. By plotting the second-order correlation function of the introduced photon-added entangled states versus r⟨∣α∣r\langle | \alpha \mid and ∣ξ⟩|\xi\rangle in Fig. 6, this nonclassical feature is observed for none of the values of photon-addition.
5.3 Quadrature squeezing
As the last proposed nonclassical feature, we study the quantum fluctuations of the quadratures of the
introduced photon-added entangled states. For this purpose, the two Hermitian operators x=(a+a†)/2x=\left(a+a^{\dagger}\right) / \sqrt{2} and p=(a−a†)/i2p=\left(a-a^{\dagger}\right) / i \sqrt{2} have been investigated which satisfy the commutation relation [x,p]=i[x, p]=i.
Using the common definitions of the variances of these operators, position (Δx)2(\Delta x)^{2} and momentum (Δp)2(\Delta p)^{2}, the quadrature squeezing is defined as sw=2(Δw)2−1,(w=x,p)s_{w}=2(\Delta w)^{2}-1,(w=x, p) or equivalently:
sx(p)=2⟨a†a⟩±⟨a2⟩±⟨a†2⟩∓⟨a⟩2∓⟨a†⟩2−2⟨a⟩⟨a†⟩s_{x(p)}=2\left\langle a^{\dagger} a\right\rangle \pm\left\langle a^{2}\right\rangle \pm\left\langle a^{\dagger 2}\right\rangle \mp\langle a\rangle^{2} \mp\left\langle a^{\dagger}\right\rangle^{2}-2\langle a\rangle\left\langle a^{\dagger}\right\rangle.
A squeezed state in x(p)x(p) satisfies the inequality −1<sx(p)<0-1<s_{x(p)}<0.
To observe the squeezing in xx or pp corresponding to the two modes of the introduced photon-added entangled states in Eq. (11), we have plotted the xx - and pp - squeezing
Fig. 4 Photon number distribution Pn(j,k)P_{n}(j, k) (the probability of finding jj photons in the first mode and kk photons in the second mode) of the even “two-mode photon-added entangled coherent-squeezed states”
∣ψ(α,ξ,m,n)⟩n|\psi(\alpha, \xi, m, n)\rangle_{n} with r=1r=1 for the different values of photon-additions a m,n=0m, n=0, b m,n=2m, n=2, c m=4,n=2m=4, n=2, and d m=2,n=4m=2, n=4
Fig. 5 Photon number distribution Pn(j,k)P_{n}(j, k) of the odd “two-mode photon-added entangled coherent-squeezed states” ∣ψ(α,ξ,m,n)⟩n|\psi(\alpha, \xi, m, n)\rangle_{n} for the different photon-additions mm and nn with r=1r=1 (all of the parameters are the same as Fig. 4)
Fig. 6 Second-order correlation function gnk(2)(0)g_{n k}^{(2)}(0) of the a even and b\mathbf{b} odd “two-mode photon-added entangled coherent-squeezed states” ∣ψ(α,ξ,m,n)⟩e,n|\psi(\alpha, \xi, m, n)\rangle_{e, n} versus r⟨∣α∣r\langle | \alpha \mid and ∣ξ⟩|\xi\rangle for the different values of photon-additions m,n=0m, n=0 (solid), m,n=2m, n=2 (dashed), m=4,n=2m=4, n=2 (dash-dot), and m=2,n=4m=2, n=4 (dot)
parameters of the even (odd) “two-mode photon-added entangled coherent-squeezed states” versus r(∣α∣r(|\alpha| and ∣ξ∣)|\xi|) in Fig. 7 (Fig. 8). It is worth noting that in Figs. 7b and 8 b , the sps_{p} and sxs_{x} curves coincide with each other.
As it can be observed in Figs. 7a and 8a, squeezing occurs only in the pp component for the first mode of the even and odd “entangled coherent-squeezed states” (no photon-additions m,n=0m, n=0 ) in Eq. (11) in the region of small r(∣α∣r(|\alpha| and ∣ξ∣)|\xi|). Moreover, it can be seen that the squeezing in pp-component is disappeared with photon-addition.
6 Summary and conclusion
In this paper, we started with a brief review on new types of entangled states named “entangled coherent-squeezed states” which have recently been introduced and generated by us. Next, using these entangled states, we introduced “two-mode photon-added entangled coherentsqueezed states”. These states can be obtained via the iterated action of two creation operators a†a^{\dagger} and b†b^{\dagger} corresponding to the first and second modes of the “entangled
coherent-squeezed states”, respectively. In the continuation, we studied the amount of the entanglement by evaluating concurrence and some of the nonclassical features of the introduced photon-added entangled states such as quantum statistics, second-order correlation function, and quadrature squeezing, too. Moreover, we considered the influence of two-mode photon-addition on the mentioned properties.
In conclusion, the outcome results are briefly summarized as follows:
- The concurrence (that showing the amount of entanglement) of the even and odd “two-mode photonadded entangled coherent-squeezed states” starts from zero, increases suddenly in the region of small r(r<2)r(r<2), and finally gets its maximum value 1 for all of the considered values of photon-additions mm and nn. Moreover, the rate of increasing of the entanglement increased with two-mode photon-addition.
- The sub-Poissonian photon-statistics appears for the first mode of the even and odd “two-mode photonadded entangled coherent-squeezed states” for all con-
Fig. 7 Squeezing parameters sp(r)s_{p}(r) (thin) and sx(r)s_{x}(r) (thick) of a the first mode and b\mathbf{b} the second mode of the even “two-mode photonadded entangled coherent-squeezed states” ∣ψ(α,ξ,m,n)⟩n|\psi(\alpha, \xi, m, n)\rangle_{n} versus r(∣α∣r(|\alpha| and ∣ξ∣)|\xi|) for the different values of photon-additions m,n=0m, n=0
(solid), m,n=2m, n=2 (dashed), m=4,n=2m=4, n=2 (dash-dot), and m=2,n=4m=2, n=4 (dot) (in b, the sps_{p} and sxs_{x} curves coincide with each other. Moreover, the m,n=2m, n=2 and m=2,n=4m=2, n=4 curves in a\mathbf{a} and the m,n=2m, n=2 and m=4,n=2m=4, n=2 curves in b\mathbf{b} coincide with each other)
Fig. 8 Squeezing parameters sp(r)s_{p}(r) (thin) and sx(r)s_{x}(r) (thick) of a the first mode and b\mathbf{b} the second mode of the odd “twomode photon-added entangled coherent-squeezed states” ∣ψ(α,ξ,m,n)⟩n|\psi(\alpha, \xi, m, n)\rangle_{n} versus r(∣α∣r(|\alpha| and ∣ξ∣)|\xi|) for the different photonadditions mm and nn (all of the parameters and explanations are the same as Fig. 7)
sidered values of photon-additions mm and nn over all of the regions of rr. Furthermore, the sub-Poissonian characteristics of the first mode of these entangled states increase with increasing the photon-addition of the first mode (m)(m). The second mode of the odd “entangled coherent-squeezed state” (no photon-additions m,n=0m, n=0 ) and also the second mode of the even and odd “two-mode photon-added entangled coherentsqueezed states” show sub-Poissonian photon-statistics for all considered values of photon-additions mm and nn in the region of very small r(∣α∣r(|\alpha| and ∣ξ∣)|\xi|). For the second mode of the even (odd) “two-mode photonadded entangled coherent-squeezed states”, the region of sub-Poissonian photon-statistics becomes more slightly wider (slightly narrower) with increasing the photon-addition of the second mode (n)(n).
- The oscillatory behaviour is clearly observed in the photon number distribution of both even and odd “two-mode photon-added entangled coherentsqueezed states”.
- Studying the second-order correlation function does not show any nonclassical feature for the even and odd “two-mode photon-added entangled coherentsqueezed states” for all considered values of photonadditions mm and nn.
- Since the squeezed states (which are the second mode of the introduced entangled states) have no quadrature squeezing, this nonclassical feature is observed only in the pp-component for the first mode (which is the coherent state) of the even and odd “entangled coher-ent-squeezed states” (no photon-additions m,n=0m, n=0 ) in the region of small r(∣α∣r(|\alpha| and ∣ξ∣)|\xi|).
Acknowledgements The author would like to express his utmost gratitude to Prof. Mohammad Kazem Tavassoly at Yazd university who taught him many things and also his family specially his wife, Masi for their encouragement, patience, and support throughout his study.
References
- C.H. Bennett, Phys. Rev. Lett. 68, 3121 (1992)
- S.L. Braunstein, H.J. Kimble, Phys. Rev. Lett. 80, 869 (1998)
- D.P. DiVincenzo, Science 270, 255 (1995)
- Schröodinger, E. Naturwissenschaften 23, 807-812, 823-828, 844-849 (1935)
- B.C. Sanders, Phys. Rev. A 45, 6811 (1992)
- S.J. Van Enk, O. Hirota, Phys. Rev. A 64, 022313 (2000)
- X. Wang, B.C. Sanders, Phys. Rev. A 65, 012303 (2001)
- X.H. Cai, L.M. Kuang, Chin. Phys. 11, 876 (2002)
- L. Zhou, L.M. Kuang, Phys. Lett. A 302, 273 (2002)
- H. Lu, L. Chen, J. Lin, Chin. Opt. Lett. 2, 618 (2004)
- S. Dey, V. Hussin, Phys. Rev. D 91, 124017 (2015)
- S. Dey, Phys. Rev. D 91, 044024 (2015)
- S. Dey, A. Fring, V. Hussin, Int. J. Mod. Phys. B 31, 1650248 (2017)
- L. Xu, L.M.J. Kuang, Phys. A. Math. Gen. 39, L191 (2006)
- R. Zhen-Zhong, J. Rui, Z. Xian-Zhou, Chin. Phys. Lett. 25, 3562 (2008)
- H.M. Li, H.C. Yuan, H.Y. Fan, Int. J. Theor. Phys. 48, 2849 (2009)
- Z. Dong-Lin, K. Le-Man, Chin. Phys. B 18, 1328 (2009)
- A. Karimi, M.K. Tavassoly, Phys. Scr. 90, 015101 (2014)
- A. Karimi, M.K. Tavassoly, Int J Theor Phys 55, 563 (2016)
- L. Hong, Phys. Lett. A 264, 265 (1999)
- A. Biswas, G.S. Agarwal, Phys. Rev. A 75, 032104 (2007)
- E. Dibakar, A. Karimi, M.K. Tavassoly, Phys. Scr. 90, 085102 (2015)
- A. Karimi, M.K. Tavassoly, Chin. Phys. B 25, 040303 (2016)
- S. Dey, V. Hussin, Phys. Rev. A 93, 053824 (2016)
- C.C. Gerry, Phys. Rev. A 55, 2478 (1997)
- M. Paternostro, M.S. Kim, B.S. Ham, Phys. Rev. A 67, 023811 (2003)
- L.M. Kuang, Z.B. Chen, J.W. Pan, Phys. Rev. A 76, 052324 (2007)
- A. Karimi, M.K. Tavassoly, J. Opt. Soc. Am. B 31, 2345 (2014)
- A. Karimi, M.K. Tavassoly, Commun. Theor. Phys. 64, 341 (2015)
- L.M. Kuang, A.H. Zeng, Z.H. Kuang, Phys. Lett. A 319, 24 (2003)
- Z. Zhi-Ming, Y. Jian, Y. Ya-Fei, Chin. Phys. Lett. 24, 352 (2007)
- S.B. Zheng, Z.B. Yang, Y. Xia, Phys. Rev. A 81, 015804 (2010)
- A. Karimi, M.K. Tavassoly, Laser Phys. 25, 115201 (2015)
- A. Karimi, M.K. Tavassoly, Quant. Inf. Proc. 15, 1513 (2016)
- J. Glauber, Phys. Rev. 131, 2766 (1963)
- C.C. Gerry, P.L. Knight, Introductory Quantum Optics (Cambridge University Press, NewYork, 2005)
- G.S. Agarwal, K. Tara, Phys. Rev. A 43, 492 (1991)
- S. Hill, W.K. Wootters, Phys. Rev. Lett. 78, 5022 (1997)
- A. Mann, B.C. Sanders, W.J. Munro, Phys. Rev. A 50, 989 (1995)
- P. Rungta, V. Buzěk, C.M. Caves, M. Hillery, G.J. Milburn, Phys. Rev. A 64, 042315 (2001)
- X.G. Wang, J. Phys. A. Math. Gen. 35, 165 (2002)
- L.M. Kuang, L. Zhou, Phys. Rev. A 68, 043606 (2003)
- L. Mandel, Opt. Lett. 4, 205 (1979)
- L. Davidovich, Rev. Mod. Phys. 68, 127 (1996)
- H. Paul, Rev. Mod. Phys. 54, 1061 (1982)
- T.Q. Song, H.Y. Fan, J. Phys. A. Math. Gen. 35, 1071 (2002)