Abstract

In this contribution we study a superposition of two finite dimensional trio coherent states (FTCS). The state is regarded as a correlated three-mode state in finite dimensional bases. The framework of Pegg and Barnett formalism, and the phase distribution in addition to the Poissonian distribution are examined. It is shown that the eigenvalue of the difference of the photon number (the q-parameter) is responsible for the non-classical phenomenon. Furthermore, the quasi-probability distribution functions (the Wigner and Q-functions) are also discussed. In this case and for the Wigner function the non-classical behavior is only reported for the odd values of the q-parameter.

Keywords

Entropy Squeezing, Atomic Wehrl Entropy, Marginal Atomic Q-Function

Share and Cite:

1. Introduction

The usual CSs introduced by Glauber [1] are eigenstates of the annihilation operator $\hat{a}$ of the harmonic oscillator. Based on Glauber’s work, the even and odd CSs were introduced (Dodonov et al., 1974) [2] . The even (odd) CSs are the symmetric (antisymmetric) combination of the CSs. They are two orthonormalized eigenstates of $\left({\hat{a}}^2\right)$ the square of the annihilation operator $\hat{a}$ and essentially have two kinds of nonclassical effects: the even CS has a squeezing but no antibunching effect, while the odd CS has an antibunching but no squeezing effect (Hillery, 1987; Xia and Guo, 1989) [3] . The simplest archetype examples of nonclassical states are of course number states, whose experimental realization is however difficult to achieve in the optical regime. Moreover, they share very few of the coherence properties that would be desirable both in practical implementations and in fundamental experiments. On the other hand, the coherent state has become a useful and necessary tool for treating ideal boson fields subjected to external pumping sources [1] [4] [5] [6] . Although any quantum state can be described in terms of the coherent states due to their resolution of unity, the coherent states themselves do not show non-classical effects such as antibunching, squeezing, etc. Another way of generalizing the coherent states is the pair coherent states (PCS). They are eigen states of the product of the two annihilation operators of the two modes and the photon number difference of these two modes. They show non-classical effects. Agarwal [7] suggested that the optical pair coherent states (PCS) can be generated via the competition of 4-wave mixing and two-photon absorption in a nonlinear medium. Another scheme has been suggested for generating vibrational pair coherent states via the motion of trapped ions in a two-dimensional trap [8] .

Another generalized is the trio-coherent state TCS [9] , whose various non-classical effects [10] as well as generation scheme [11] have been investigated, is an eigenstate of the operators $\hat{a}\hat{b}\hat{c}$ , ${\hat{n}}_a-{\hat{n}}_c=\hat{P}$ and ${\hat{n}}_b-{\hat{n}}_c=\hat{Q}$ such that

$\hat{a}\hat{b}\hat{c}|\xi ,p,q\rangle f=\xi |\xi ,p,q\rangle f$ , $\hat{Q}|\xi ,p,q\rangle f=q|\xi ,p,q\rangle f$

and $\hat{P}|\xi ,p,q\rangle f=p|\xi ,p,q\rangle f$ ,

where p and q are positive integers with $\xi$ is complex number respectively.

Finite dimensional states have also been introduced. The finite dimensional pair coherent state $\left(FPCS\right)|\xi ,q\rangle$ has been studied recently as the eigenstate of the pair operators $\left({\hat{a}}^{†}\hat{b}+\frac{{\zeta }^{q+1}{\left(\hat{a}{\hat{b}}^{†}\right)}^q}{{\left(q!\right)}^2}\right)$ and the sum of the photon number operators for the two modes $\left(\hat{Q}={\hat{n}}_a+{\hat{n}}_b\right)$ [12] . Also the superposition of FPCS and some of its properties have been studied recently [13] .

In the present communication we develop this idea and introduce the superposition of the finite dimensional trio coherent state (SFTCS). The correlated three-mode for the FTCS is defined as the eigenstate of the following three operators

$\hat{X}=\left({\hat{a}}^{†}\hat{b}\hat{c}+\frac{{\zeta }^{q+1}p!{\left(\hat{a}{\hat{b}}^{†}{\hat{c}}^{†}\right)}^q}{{\left(q!\right)}^2\left(p+q\right)!}\right)$

for the three modes, and the sum of the photon number operators for the three modes ${\hat{a}}^{†}\hat{a}+{\hat{c}}^{†}\hat{c}=\hat{Q}$ and ${\hat{b}}^{†}\hat{b}-{\hat{c}}^{†}\hat{c}=\hat{P}$ , namely [14] :

$\begin{array}{l}\left({\hat{a}}^{†}\hat{b}\hat{c}+\frac{{\zeta }^{q+1}p!{\left(\hat{a}{\hat{b}}^{†}{\hat{c}}^{†}\right)}^q}{{\left(q!\right)}^2\left(p+q\right)!}\right)|\zeta ,p,q\rangle =\zeta |\zeta ,p,q\rangle \\ \left({\hat{a}}^{†}\hat{a}+{\hat{c}}^{†}\hat{c}\right)|\zeta ,p,q\rangle =q|\zeta ,p,q\rangle \\ \left({\hat{b}}^{†}\hat{b}-{\hat{c}}^{†}\hat{c}\right)|\zeta ,p,q\rangle =p|\zeta ,p,q\rangle \end{array}\}$ (1)

where the parameter $\zeta$ is a complex variable while the parameters $q\ge 0$ and p are integers and the operators $\hat{a}$ , $\hat{b}$ and $\hat{c}$ commute pairwise. The state takes the form,

$|\zeta ,q,p\rangle =N_{q,p}{\displaystyle \underset{n=0}{\overset{q}{\sum }}{\zeta }^n}\sqrt{\frac{\left(q-n\right)!p!}{q!\left(n+p\right)!n!}}|q-n,n+p,n\rangle ,$ (2)

in the three mode states $|n_a,n_b,n_c\rangle =|n_a\rangle \otimes |n_b\rangle \otimes |n_c\rangle$ , where $|n_s\rangle$ is the Fock state for the mode s(s = a, b and c) and the normalization constant $N_{q,p}$ is given by

$N_{q,p}={\left[{\displaystyle \underset{n=0}{\overset{q}{\sum }}{\left|\zeta \right|}^{2n}\frac{\left(q-n\right)!p!}{q!\left(n+p\right)!n!}}\right]}^{\frac{-1}{2}},$ (3)

the parameter p may be a negative integer, in such case the summation starts from $n=\left|p\right|$ , because of the appearance of the operators ${\hat{a}}^{†}\hat{b}\hat{c}$ or $\hat{a}{\hat{b}}^{†}{\hat{c}}^{†}$ in this form it may be legitimate to call it a finite dimensional trio coherent state.

In this paper we address the problem of constructing and discussing some non-classical properties of the superposition of two finite dimensional trio coherent states. We shall examine some values of the parameters of the state obtained.

The correlated three-mode Schrodinger-cat states $|\zeta ,q,p,\varphi \rangle$ are defined as superpositions of three (FTCS) separated in phase by $\varphi$ (Gerry and Grobe, 1995) [15] :

$|\zeta ,q,p,\varphi \rangle =N_{q,p}^{\varphi }\left(|\zeta ,q,p\rangle +e^{i\varphi }|-\zeta ,q,p\rangle \right),$ (4)

where the normalization constant $N_{q,p}^{\varphi }$ is given by

$N_{q,p}^{\varphi }=\frac{1}{\sqrt{2}}{\left[1+\cos \varphi N_{q,p}^2\underset{n=0}{\overset{q}{{\displaystyle \sum }}}{\left(-1\right)}^n{\left|\zeta \right|}^{2n}\frac{\left(q-n\right)!p!}{q!\left(n+p\right)!n!}\right]}^{\frac{-1}{2}}$ (5)

Now we discuss some statistical properties of these correlated two mode states of Equation (4). The results that we are going to present stem from a new approach to the superposing of the finite dimensional state. Subsequently we shall examine the phase distribution in the framework of Pegg and Barnett formalism, the behavior of the sub-Poissonian distribution, the Wigner function and the Q-function of the state (4) are discussed.

2. Phase Distribution

In this section we shall look at the phase distribution of the above mentioned state. In the Pegg-Barnett formalism for the phase, a Hermitian phase operator is defined in a finite dimensional state space [16] . They used the fact that, in this finite dimensional Hilbert space, one can define phase states rigorously. The phase operator is then defined as the projection operator on the particular phase state multiplied by the corresponding value of the phase. The main idea of the formalism consists in evaluation of all expectation values of physical variables in a finite dimensional Hilbert space. These give real numbers, which depend parametrically on the dimension of the space. Because a complete description of the harmonic oscillator involves an infinite number of states to be taken, a limit is taken only after the physical results are evaluated. This leads to proper limits which correspond to the results obtainable in ordinary quantum mechanics. Here, we study phase properties for the state (4) using this formalism.

The Pegg-Barnett phase distribution $P\left(\eta \right)$ is defined [16] as

$N_{q,p}^{\varphi }=\frac{1}{\text{2π}}{\displaystyle \underset{l,m}{\sum }{\rho }_{lm}\exp \left[i\left(l-m\right)\left(\eta -{\eta }_0\right)\right]}$ (6)

the angle ${\eta }_0$ is the phase reference angle and we take it to be zero. In the case of (4) the phase probability distribution is generalized

$\begin{array}{l}{P}_{q,p}\left({\theta }_1,{\theta }_2,{\theta }_3,\zeta \right)\\ =\frac{{\left(N_{q,p}\right)}^2}{{\left(\text{2π}\right)}^3}\{|\underset{n,m=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}{\zeta }^{2n+j}{\left({\zeta }^{*}\right)}^{2m+j}\sqrt{\frac{\left(q-2n-j\right)!p!\left(q-2m-j\right)!p!}{q!\left(2n+j+p\right)!\left(2n+j\right)!q!\left(2m+j+p\right)!\left(2m+j\right)!}}\\ \cdot \exp i[\left(q-2n-j\right){\theta }_1-\left(q-2m-j\right){\theta }_1+\left(2n+p-j\right){\theta }_2\\ \begin{array}{c}\\ \\ \end{array}-\left(2m+p-j\right){\theta }_2+\left(2n-j\right){\theta }_3-\left(2m-j\right){\theta }_3]|\}\end{array}$

Therefore the phases distribution function can be written as

$P_{q,p}\left(\theta ,\zeta \right)=\frac{{\left(N_{q,p}\right)}^2}{{\left(\text{2π}\right)}^3}{\left\{\underset{n,m=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}{\zeta }^{2n+j}\sqrt{\frac{\left(q-2n-j\right)!p!}{q!\left(2n+j+p\right)!\left(2n+j\right)!}}\exp \left[i2n\left(\theta \right)\right]\right\}}^2$ (7)

which is normalized according to ${\displaystyle \int {\displaystyle \int {\displaystyle {\int }_{-\text{π}}^{\text{π}}{P}_{\zeta ,q,p}\left({\theta }_1,{\theta }_2,{\theta }_3,\zeta \right)\text{d}{\theta }_1\text{d}{\theta }_2\text{d}{\theta }_3}}}=1$ . Due to the connection between the three modes, the phase distribution depends on the difference between the phases of the modes. In the figures we plot $P_{\zeta ,q,p}\left(\theta \right)$ against the angle $\theta ={\theta }_2+{\theta }_3-{\theta }_1$ , $-\text{π}\le \theta \le \text{π}$ for different values of the parameter $q,p$ and $\left|\zeta \right|$ .

To demonstrate the behavior of the distribution function $P_{q,p}\left(\theta ,\zeta \right)$ we have plotted Figure 1 for $q=p=2$ and different values of the parameter $\zeta$ adjust $\zeta =3,5$ . For the j = 0 case in the earlier study [13] the distribution function shows Gaussian behavior with one peak centered symmetrically around the center of the base, but in the present contribution the function $P_{q,p}\left(\theta ,\zeta \right)$ decreases its maximum and the peak splits into two separated identical peaks, both are symmetric with respect to the origin, see Figure 1(a). The maxima for the distribution at $\theta =0$ decrease by increasing of the parameter $\zeta$ . In Figure 1(b) we take larger values for the parameter $q\left(q=0\right)$ and the same values of $\zeta$ . We see that the function $P_{\zeta ,q}\left(\theta ,\zeta \right)$ starts at $P_{\zeta ,q}\left(\theta ,\zeta \right)=0.056$ , (0.125) when $\left|\zeta \right|=3$ , (5) respectively. The maxima of the phase distribution are increased by increasing of the parameter $\zeta$ (see Figure 1(b)). However this increase turns to a decrease for larger values of $\zeta$ . The maximum value for $P_{\zeta ,q}\left(\theta ,\zeta \right)$ shifts

to higher values of $\zeta$ as q increases. In general the function decreases its height as the q-parameter increases. This means that an increase in the value of the q-parameter changes the Pegg-Barnett phase distribution function.

3. Sub-Poissonian Distribution

In this section we discuss an example of non-classical effects of physical states. One of the non-classical phenomena of the quantized electromagnetic radiation field is the sub-Poissonian statistics. A state of a single mode, for convenience, displaying sub-Poissonian statistics is characterized by the fact that the variance of the photon number ${\langle \Delta n_i\rangle }^2$ is less than the average photon number $\langle a_i^{†}{a}_i\left(t\right)\rangle =\langle n_i\rangle$ . Sub-Poissonian statistics does not occur in any classical description for the light field, but only occurs when the field is described quaintly [17] . This can be expressed by means of the normalized second-order correlation function as [17] :

$g_z^{\left(2\right)}\left(\zeta \right)=\frac{{}_j\langle \zeta ,p,q|{\hat{n}}_z\left({\hat{n}}_z-1\right){|\zeta ,p,q\rangle }_j}{{}_j\langle \zeta ,p,q|{\hat{n}}_z{|\zeta ,p,q\rangle }_j^2},\forall z=a,b,c$ (8)

where

$\begin{array}{l}{}_j\langle \zeta ,p,q|{\hat{n}}_a\left({\hat{n}}_a-1\right){|\zeta ,p,q\rangle }_j\\ =N_{q,p}^2\underset{n=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}{\left|\zeta \right|}^{2\left(2n+j\right)}\frac{\left(q-2n-j\right)!p!}{q!\left(2n+j+p\right)!\left(2n+j\right)!}\left(q-2n-j\right)\left(q-2n-j-1\right),\\ {}_j\langle \zeta ,p,q|{\hat{n}}_b\left({\hat{n}}_b-1\right){|\zeta ,p,q\rangle }_j\\ =N_{q,p}^2\underset{n=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}{\left|\zeta \right|}^{2\left(2n+j\right)}\frac{\left(q-2n-j\right)!p!}{q!\left(2n+j+p\right)!\left(2n+j\right)!}\left(p+2n-j\right)\left(p+2n-j-1\right),\\ {}_j\langle \zeta ,p,q|{\hat{n}}_c\left({\hat{n}}_c-1\right){|\zeta ,p,q\rangle }_j\\ =N_{q,p}^2\underset{n=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}{\left|\zeta \right|}^{2\left(2n+j\right)}\frac{\left(q-2n-j\right)!p!}{q!\left(2n+j+p\right)!\left(2n+j\right)!}\left(2n-j\right)\left(2n-j-1\right).\end{array}\}$ (9)

and

$\begin{array}{l}{}_j\langle \zeta ,p,q|{\hat{n}}_c{|\zeta ,p,q\rangle }_j^2=N_{q,p}^2\underset{n=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}{\left|\zeta \right|}^{2\left(2n+j\right)}\frac{\left(q-2n-j\right)!p!}{q!\left(2n+j+p\right)!\left(2n+j\right)!}\left(q-2n-j\right),\\ {}_j\langle \zeta ,p,q|{\hat{n}}_c{|\zeta ,p,q\rangle }_j^2=N_{q,p}^2\underset{n=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}{\left|\zeta \right|}^{2\left(2n+j\right)}\frac{\left(q-2n-j\right)!p!}{q!\left(2n+j+p\right)!\left(2n+j\right)!}\left(p+2n-j\right),\\ {}_j\langle \zeta ,p,q|{\hat{n}}_c{|\zeta ,p,q\rangle }_j^2=N_{q,p}^2\underset{n=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}{\left|\zeta \right|}^{2\left(2n+j\right)}\frac{\left(q-2n-j\right)!p!}{q!\left(2n+j+p\right)!\left(2n+j\right)!}\left(2n-j\right).\end{array}\}$ (10)

where the subscript z relates to the zth mode.

A light field has a sub-Poissonian distribution if $g_z^{\left(2\right)}<1$ , which is a nonclassical effect and means that the probability of detecting an incident pair of photons is less than it would be for a coherent field described by the Poissonian distribution; super-Poissonian distribution if $g_z^{\left(2\right)}>1$ , which is a classical effect ( $g_z^{\left(2\right)}=2$ for thermal light) and Poissonian distribution of photons (standard for the coherent state) if $g_z^{\left(2\right)}=1$ .

For fixed the parameter $p=2,j=0$ we examine the correlation function for the first mode $g_z^{\left(2\right)}\left(\zeta \right)$ we have plotted Figure 2. For the first mode the function $g_z^{\left(2\right)}\left(\zeta \right)$ shows sub-Poissonian behavior for all even values of the parameter q, however for all range of the parameter $\zeta$ see Figure 2(a). For the second mode and for same values of the parameter q the distribution sub-Poissonian for a short range of $\zeta$ . When the parameter $\zeta$ is increased further, the trio-coherent state $|\zeta ,q,p\rangle$ exhibits sup-Poissonian behavior. More increases of the variable $\zeta$ the distribution exhibits thermal distribution as well as super-thermal distribution as observed in Figure 2(b). The nonclassical nature of the state is apparent, which $\zeta$ takes on small values shown in Figure 2(b). In the meantime, it is observed that for $q=4$ sub-Poissonian behavior for almost the same range for the cases in which $q=8$ . On the other hand the correlation function $g_z^{\left(2\right)}\left(\zeta \right)$ for $q=8$ exhibits thermal distribution faster than the case in which $q=4$ . Furthermore, we can also observe that the function starts sub-Poissonian distribution

in both modes, but for large values of the parameter $\zeta$ the distribution becomes sub-Poissonian and super thermal states for first and second modes respectively. We note that the behavior of the case $j=1$ for different values of q same as the case $j=0$ with exchanges of the modes.

4. Quasiprobability Distribution Sub-Poissonian Distribution

It is well known that there are three quasiprobability distribution functions: P-representation, W-Wigner, and Q-function. These functions are regarded as important tools to provide insight into the nonclassical features of the radiation fields. In the meantime, they have advantages and disadvantages connected with their use. As a part of disadvantage the P-function (which describes a quantum state in terms of the probability that the system is in a given coherent state) is highly singular or negative for quantum states with no classical analogues. While the Wigner function may become negative for some quantum states, but it has the considerable advantage for squeezed states that its contour map out the variances in the field quadrature. The Q-function is a positive-definite quasiprobability distribution, but its simple relation to anti-normal operator products makes it difficult to interpret in terms of conventional photon counting or squeezing measurements. They are defined by taking the Fourier transforms of their respective characteristic functions:

$\begin{array}{c}I\left(\alpha ,\beta ,\gamma ,s\right)=\frac{1}{{\text{π}}^6}\underset{-\infty }{\overset{\infty }{{\displaystyle \int }}}\exp \left({\alpha }^{*}\eta -\alpha {\eta }^{*}\right)\exp \left({\beta }^{*}\chi -\beta {\chi }^{*}\right)\\ \times \exp \left({\gamma }^{*}\xi -\gamma {\xi }^{*}\right)C\left(\chi ,\eta ,\xi ,s\right){\text{d}}^2\chi {\text{d}}^2\eta {\text{d}}^2\xi \end{array}$ (11)

where $C\left(\chi ,\eta ,\xi ,s\right)$ is the s-parameterized characteristic function with the parameter $s=1,0,-1$ corresponding to P-representation, Wigner, and Q-function, respectively. The characteristic function can be evaluated through the relation

$\begin{array}{c}C\left(\chi ,\eta ,\xi ,s\right)=Tr\{\hat{\rho }\exp \left({\hat{b}}_1^{†}\eta -{\hat{b}}_1{\eta }^{\ast }\right)\exp \left({\hat{b}}_2^{†}\chi -{\hat{b}}_2{\chi }^{\ast }\right)\begin{array}{c}\\ \end{array}\\ \times \exp \left({\hat{b}}_3^{†}\xi -{\hat{b}}_3{\xi }^{*}\right)\left[\frac{s}{2}\left({\left|\eta \right|}^2+{\left|\chi \right|}^2+{\left|\xi \right|}^2\right)\right]\},\end{array}$ (12)

where $\hat{\rho }$ is the density matrix given by

$\hat{\rho }=|q,\zeta ,\xi \rangle \langle q,\zeta ,\xi |$ (13)

In what follows we consider the Wigner, and the Q-function and for this reason we have to evaluate the integral in Equation (11) for $s=0$ and $s=-1$ , respectively. This can be achieved if one manages to calculate the characteristic function. From Equation (12) and after minor algebra we have

$\begin{array}{l}C\left(\chi ,\eta ,\xi ,s\right)\\ =N_q^2\exp \left[-\frac{1}{2}\left(1-s\right)\left({\left|\eta \right|}^2+{\left|\chi \right|}^2+{\left|\xi \right|}^2\right)\right]\underset{n=0}{\overset{q}{{\displaystyle \sum }}}\underset{r=0}{\overset{q}{{\displaystyle \sum }}}{B}_n{B}_r{L}_n\left({\left|\chi \right|}^2\right)L_{p+n}\left({\left|\eta \right|}^2\right)L_{q-n}\left({\left|\xi \right|}^2\right)\end{array}$ (14)

$\begin{array}{l}=\exp \left[-\frac{\left(1-s\right)}{2}\left({\left|\eta \right|}^2+{\left|\chi \right|}^2+{\left|\xi \right|}^2\right)\right]\underset{n=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}\underset{m=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}{B}_n\left(\zeta ,q\right)B_m^{*}\left(\zeta ,q\right)\\ \times \sqrt{\frac{\left(q-2n-j\right)!\left(p+2n+j\right)!\left(2n+j\right)!}{\left(q-2m-j\right)!\left(p+2m+j\right)!\left(2m+j\right)!}}{L}_{q-2n-j}^{2n-2m}\left({\left|\chi \right|}^2\right)L_{2n+j}^{2m-2n}\left({\left|\eta \right|}^2\right)L_{p+2n+j}^{2m-2n}\end{array}$ (15)

where $L_m^n\left(x\right)$ are associated Laguerre polynomials and the coefficient $B_n\left(\zeta ,q\right)$ given by

$\begin{array}{l}{L}_m^n\left(x\right)=\underset{r=0}{\overset{m}{{\displaystyle \sum }}}\left(\begin{array}{c}m+n\\ m-r\end{array}\right)\frac{{\left(-1\right)}^r}{r!}{x}^r\\ B_n\left(\zeta ,q\right)={\zeta }^{2n+j}\sqrt{\frac{\left(q-2n-j\right)!p!}{q!\left(2n+j+p\right)!\left(2n+j\right)!}}\end{array}$ (16)

Having obtained the parameterized characteristic function, we are therefore in a position to find the Wigner, and Q-function. This will be evaluated in the next subsections.

4.1. The Wigner Function

To obtain the Wigner function $W\left(\alpha ,\beta \right)$ we insert Equation (14) into Equation (11) and perform the integral. For $s=0$ thus, we obtain

$\begin{array}{l}W\left(\alpha ,\beta \right)=\frac{4}{{\text{π}}^3}{N}_q^2\exp \left[-2\left({\left|\eta \right|}^2+{\left|\chi \right|}^2+{\left|\xi \right|}^2\right)\right]\\ \times \underset{n=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}\sqrt{\left(2n+j\right)!\left(p+2n+j\right)!\left(q-2n-j\right)!\left(2r+j\right)!\left(p+2r+j\right)!\left(q-2r-j\right)!}\\ \times {\left(-1\right)}^{q+p-2n-j}{\zeta }^{2\left(n+r+j+1\right)}{\left(\frac{\beta \gamma }{\alpha }\right)}^{2\left(r-n\right)}{L}_{q-2n-j}^{2n-2r}\left(4{\left|\alpha \right|}^2\right)L_{2n+j}^{2r-2n}\left(4{\left|\beta \right|}^2\right)L_{p+2n+j}^{2r-2n}\left({\left|\gamma \right|}^2\right)\end{array}$ (17)

In Figure 3 we have plotted the Wigner function against $Re\left(\alpha \right)$ and $Im\left(\alpha \right)$ for fixed values of $\zeta =5$ and $p=2$ , assuming that $\alpha ,\beta$ and $\gamma$ are coincidence i.e. in the subspace $\alpha =\beta =\gamma$ [18] . In the meantime we examined the function for the cases in which $j=0$ and $q=0,4,6$ and 8. When we consider the case in which $q=0$ , the function displays Gaussian shape with a symmetrical behavior around zero.

In this case one can see a sharp peak centered at the middle of the base, see Figure 3(a). As soon as we consider the value of the q-parameter and take, $j=1$ and $q=1$ , the peak of the function gets upsidedown and the nonclassical effect becomes pronounced. This is clear from Figure 3(b) where the negative values of the function are apparent. The spreading of Wigner over the plane is shown as q increases, this is seen for the case in which $j=0$ and $q=6$ . In this case the function increases its value and starts with a doubly folded peak in the center towards the positive side. This indicates that the function gets more sensitive to the variation in the q-parameter and this of course reflects the change from Gaussian to Fock states; see Figure 3(c). For $j=1$ and $q=5$ the function displays the same shape, however, it changes its direction upsidedown and exhibits

negative values. In the meantime, the width of the base increases and the doubly folded peak gets more pronounced, see Figure 3(d). This indicates that the nonclassical behavior appears only for the odd numbers of the q-parameter while it disappears for the even numbers. This means that the q-parameter plays a role of changing the nonclassical behavior.

4.2. Q-Function

In the following we concentrate on one of the quasiprobability distribution functions, that is the Q-function [3] [19] . The Q-function is a convenient tool to calculate expectation values of anti-normally ordered products of operators. It is well known that the Q-function can be defined in terms of diagonal elements of the density operator in the coherent state. Therefore we use the density operator $\rho$ for the state (4) to study its quasi-probability distribution. The Q-function in this case is given by

$Q\left(\alpha ,\beta ,\gamma \right)=\frac{1}{{\text{π}}^3}{\left|\langle \alpha ,\beta ,\gamma |\zeta ,q,p\rangle \right|}^2,$ (18)

where $\alpha ,\beta ,\gamma \in ℂ$ , $ℂ$ is a complex number and $|\alpha ,\beta ,\gamma \rangle =|\alpha \rangle |\beta \rangle |\gamma \rangle$ , with $|\alpha \rangle ,|\beta \rangle$ and $|\gamma \rangle$ are the usual coherent states. Generally there are six variables associated with the real and imaginary parts of $\alpha ,\beta ,\gamma$ . For visualization let us confine ourselves to a subspace determined by $\alpha =\beta =\gamma$ [18] . In that subspace the Q-function for the state (4) is calculated to be

$Q\left(\alpha ,\beta ,\gamma \right)=\frac{\exp \left[-\frac{3}{2}\left(x^2+y^2\right)\right]}{{\text{π}}^3}{\left|N_{q,p}\underset{n=0}{\overset{\left[\frac{q-j}{2}\right]}{{\displaystyle \sum }}}\frac{{\zeta }^{2n+j}{\alpha }^{q+p+2n+j}\sqrt{p!}}{\sqrt{q!}\left(p+2n+j\right)!\left(2n+j\right)!}\right|}^2,$ (19)

where $x=Re\left(\alpha \right)$ and $y=Im\left(\alpha \right)$ .

Since the maximization or minimization of the Q-function depends on the parameter q. Therefore, our main task is to examine the behavior of the Q-function result of the variation in the q-parameter. For this reason we plot Figure 4 for different values of the q-parameter keeping all other parameters unchanged as in the Wigner function case. For instance when we consider $j=0$ and $q=0$ , the function exhibits Gaussian behavior similar to that of the Wigner function, see Figure 4(a). In this case there exists a unique maximum value at the origin and consequently the function $Q\left(x,y\right)$ for the state $|0,\zeta \rangle$ has one peak centered at $\left(x,y\right)=\left(0,0\right)$ . For the case in which $q>0$ , more precisely $q=4$ , there exists maxima and minima at origin. For the case in which $q=4$ , we observe that the function changes its behavior from Gaussian to almost the Fock state. This is quite obvious from the observation of the doubly folded peaks beside an increase in the width of the base, see Figure 4(b). The increase in the value of q-parameter leads to an increase in the width of the base, this in addition to a decrease in the maximum value of the function. Thus we may conclude that, the shape of the Q-function for the state $|q,\zeta \rangle$ where $q\ne 0$ is sensitive to the change in q-parameter where a crater is apparent in the center. However, if we increase q the crater-like at the center spreads out in the phase space and the diameter increases as the q.

5. Conclusion

In the present paper we have studied the superposition of two entangled finite dimensional trio coherent states. These states can be produced by processes in which there is a strong competition between a three mode parametric conversions. In the mean time, we have employed the Glauber second order-correlation function to examine the nonclassical properties of the state. We have shown that the nonclassical as well as the classical behaviors are apparent in both modes for large values of $\zeta$ . However, the nonclassical is more pronounced in the first mode case while the classical behavior is pronounced in the second mode. We have also considered the quasiprobability distribution functions (the Wigner and Q-functions) where observation of nonclassical properties is reported for the odd values of the q-parameter. In the meantime the Q-function displays Gaussian

behavior and tends to Fock state behavior as q increases. Finally, we have examined the properties of the present state in terms of the phase distribution function introduced by Barnett and Pegg. In this case, regardless of the value of q, the function is symmetric about zero. However, as the q-parameter increases the value of the function decreases but without broking the symmetry.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References