different aspects of nonclassical properties of electromagnetic field have been studied since the advent of quantum optics .
however , the interest on nonclassical states has been considerably escalated with the progress of interdisciplinary field of quantum computation and quantum communication in recent past as a large number of applications of nonclassical states have been reported in context of quantum computation and quantum communication @xcite .
specifically , it is shown that squeezed states can be used for the implementation of continuous variable quantum cryptography @xcite and teleportation of coherent states @xcite , antibunched states can be used to build single photon sources @xcite , entangled states are essential for the implementation of a set of protocols of discrete @xcite and continuous variable quantum cryptography @xcite , quantum teleportation @xcite , dense - coding @xcite , etc . , states violating bell s inequality are useful for the implementation of protocols of device independent quantum key distribution @xcite .
study of the possibility of generation of nonclassical states ( specially , entanglement and nonlocal characters of quantum states ) in different quantum systems have recently become extremely relevant and important for the researchers working in different aspects of quantum information theory and quantum optics .
existence of entanglement and other nonclassical states in a large number of bosonic systems have already been reported ( see @xcite and references therein ) .
however , it is still interesting to study the possibility of generation of nonclassical states in experimentally realizable simple systems .
a specific system of this kind is a nonlinear optical coupler which can be easily realized using optical fibers or photonic crystals .
optical couplers are interesting for several reasons .
firstly , in a coupler the amount of nonclassicality can be controlled by controlling the interaction length and the coupling constants .
further , optical couplers are of specific interest as recently matthews et al .
have experimentally demonstrated manipulation of multiphoton entanglement in quantum circuits constructed using optical couplers @xcite , mandal and midda have shown that nand gate ( thus , in principle a classical computer ) can be build using nonlinear optical couplers @xcite .
motivated by these facts we aim to systematically investigate the possibility of observing nonclassicality in nonlinear optical couplers . as a first effort
we have reported lower order and higher order entanglement and other higher order nonclassical effects in an asymmetric codirectional nonlinear optical coupler which is prepared by combining a linear waveguide and a nonlinear ( quadratic ) waveguide operated by second harmonic generation @xcite .
waveguides interact with each other through evanescent wave . extending the earlier investigation @xcite here we study a similar asymmetric nonlinear optical coupler for contradirectional propagation of fields .
this type of contradirectional asymmetric optical coupler was studied earlier by some of us @xcite . however , in the earlier studies @xcite intermodal entanglement and some of the higher order nonclassical properties studied here were not studied .
further , in those early studies second - harmonic mode ( @xmath7 was assumed to be pumped with a strong coherent beam . in other words , @xmath8 mode was assumed to be classical and thus it was beyond the scope of the previous studies to investigate single - mode and intermodal nonclassicalities involving this mode .
interestingly , completely quantum mechanical treatment adopted in the present work is found to show intermodal squeezing in a compound mode involving @xmath8 mode .
in addition , conventional short - length solutions of heisenberg s equations of motion were used in earlier studies @xcite , but recently it is established by some of us that using improved perturbative solutions obtained by sen - mandal approach @xcite we can observe several nonclassical characters that are not observed using short - length ( or short - time approach ) @xcite . keeping these facts in mind
present paper aims to study nonclassical properties of this contradirectional asymmetric nonlinear optical coupler with specific attention to entanglement using perturbative solution obtained by sen - mandal technique .
nonclassical properties of different optical couplers are extensively studied in past ( see @xcite for a review ) .
specially , signatures of nonclassicality in terms of photon statistics , phase properties and squeezing were investigated in codirectional and contradirectional kerr nonlinear couplers having fixed and varying linear coupling constant @xcite , raman and brillouin coupler @xcite and parametric coupler @xcite , asymmetric @xcite and symmetric @xcite directional nonlinear coupler etc . here
it would be apt to note that the specific coupler system that we wish to study in the present paper has already been investigated @xcite for contradirectional propagation of classical ( coherent ) input modes and it is shown that depth of nonclassicality may be controlled by varying the phase mismatching @xmath9 @xcite .
existing studies are restricted to the investigation of lower order nonclassical effects ( e.g. , squeezing and antibunching ) under the conventional short - length approximation .
only a few discrete efforts have recently been made to study higher order nonclassical effects and entanglement in optical couplers @xcite , but except a recent study on codirectional nonlinear optical coupler reported by us @xcite , all the other efforts were limited to kerr nonlinear optical coupler .
for example , in 2004 , leonski and miranowicz reported entanglement in kerr nonlinear optical coupler @xcite and pumped kerr nonlinear optical coupler @xcite , subsequently entanglement sudden death @xcite and thermally induced entanglement @xcite were reported in the same system .
amplitude squared ( higher order ) squeezing was also reported in kerr nonlinear optical coupler @xcite .
however , no effort has yet been made to rigorously study the higher order nonclassical effects and entanglement in contradirectional nonlinear optical coupler .
keeping these facts in mind in the present letter we aim to study nonclassical effects ( including higher order nonclassicality and entanglement ) in contradirectional nonlinear optical coupler .
remaining part of the paper is organized as follows . in section [ sec : the - model - and ] , the model momentum operator that represents the asymmetric nonlinear optical coupler is described and perturbative solutions of equations of motion corresponding to different field modes present in the momentum operator are reported . in section [ sec : criteria - of - nonclassicalities ] , we briefly describe a set of criteria of nonclassicality . in section [ sec : nonclassicality - in - contradirectional ] the criteria described in the previous section are used to investigate the existence of different nonclassical characters ( e.g. , lower order and higher order squeezing , antibunching , and entanglement ) in various field modes present in the contradirectional asymmetric nonlinear optical coupler . finally , section [ sec : conclusions ] is dedicated for conclusions .
a schematic diagram of a contradirectional asymmetric nonlinear optical coupler is shown in fig . [
fig : schematic diagram ] . from fig .
[ fig : schematic diagram ] one can easily observe that a linear waveguide is combined with a nonlinear @xmath10 waveguide to constitute the asymmetric coupler of our interest . further , from fig . [ fig : schematic diagram ] we can observe that in the linear waveguide field propagates in a direction opposite to the propagation direction of the field in the nonlinear waveguide . electromagnetic field characterized by the bosonic field annihilation ( creation ) operator @xmath11 propagates through the linear waveguide .
similarly , the field operators @xmath12 corresponds to the nonlinear medium .
specifically , @xmath13 and @xmath14 denote annihilation operators ( wave vectors ) for fundamental and second harmonic modes , respectively .
now the momentum operator for contradirectional optical coupler is @xcite@xmath15 where @xmath16 stands for the hermitian conjugate and @xmath17 represents the phase mismatch between the fundamental and second harmonic beams .
the linear ( nonlinear ) coupling constant , proportional to susceptibility @xmath18 @xmath19 , is denoted by the parameter @xmath20 @xmath21 .
it is reasonable to assume @xmath22 as in a real physical system we usually obtain @xmath23 . as a consequence , in absence of a highly strong pump @xmath24 .
earlier this model of contradirectional optical coupler was investigated by some of the present authors ( @xcite and references therein ) . using ( [ eq:1 ] ) and the procedure described in @xcite we can obtain the coupled differential equations for three different modes as follows @xmath25 here it would be apt to mention that momentum operator for contradirectional asymmetric nonlinear coupler ( [ eq:1 ] ) is same as that of codirectional asymmetric nonlinear optical coupler @xcite . however , for the contradirectional couplers the sign of derivative in the heisenberg s equation of motion of the contra - propagating mode is changed ( i.e. , in the present case @xmath26 is replaced by @xmath27 as mode @xmath28 is considered here as the contra - propagating mode ) .
the method used here to obtain ( [ eq : aaa ] ) is described in refs .
@xcite . further , this particular description of contradirectional coupler is valid only for the situation when the forward propagating waves reach @xmath29 and the counter ( backward ) propagating wave reach @xmath30 .
thus the coupled equations described by ( [ eq : aaa ] ) and their solution obtained below are not valid for @xmath31 @xcite .
earlier these coupled equations ( [ eq : aaa ] ) were solved under short - length approximation . here
we aim to obtain perturbative solution for these equations using sen - mandal method which is already shown to be useful in detecting nonclassical characters not identified by short - length solution @xcite . keeping this in mind
, we plan to solve ( [ eq : aaa ] ) using sen - mandal approach .    ) with a nonlinear ( @xmath32 ) waveguide operated by second harmonic generation .
the fields involved are described by the corresponding annihilation operators , as shown ; @xmath33 is the interaction length . ]    using sen - mandal approach we have obtained closed form perturbative analytic solutions of ( [ eq : aaa ] ) as@xmath34 with@xmath35\sinh|k|l\right\ } , \\
h_{4 } & = & \frac{c^{*}|k|k}{k^{*}}f_{1}\left\ { 2|k|^{2}f_{1}g_{-}^{*}+\delta k\sinh|k|l\left(g_{+}^{*}-1\right)\left[2i|k|+\delta k\tanh|k|l\right]\right\ } , \end{array}\label{eq : soln}\ ] ] where @xmath36 and @xmath37 .
the solution obtained above is verified by escr ( equal space commutation relation ) which implies @xmath38=\left[b_{1}\left(l\right),b_{1}^{\dagger}\left(l\right)\right]=\left[b_{2}\left(l\right),b_{2}^{\dagger}\left(l\right)\right]=1 $ ] while all other equal space commutations are zero .
further , we have verified that the solutions reported here satisfies constant of motion . to be precise , in ref .
@xcite , it was shown that the constant of motion for the present system leads to @xmath39 here we have verified that the solution proposed here satisfies ( [ eq : const - of - motion ] ) .
we have also verified that the solution of contradirectional coupler using short - length solution method reported in @xcite can be obtained as a special case of the present solution .
specifically , to obtain the short - length solution we need to expand the trigonometric functions present in the above solution and neglect all the terms beyond quadratic powers of @xmath33 and consider phase mismatch @xmath40 .
after doing so we obtain@xmath41 which coincides with the short - length solution reported earlier by some of the present authors @xcite .
clearly the solution obtained here is valid and more general than the conventional short - length solution as the solution reported here is fully quantum solution and is valid for any length , restricting the coupling constant only .
further , in case of codirectional optical coupler we have already seen that several nonclassical phenomena not identified by short - length solution are identified by the perturbative solution obtained by sen - mandal method @xcite . keeping this in mind , in what follows we investigate nonclassical characters of the fields that have propagated through a contradirectional asymmetric nonlinear optical coupler .
as the criteria for obtaining signatures of different nonclassical phenomena are expressed in terms of expectation values of functions of annihilation and creation operators of various modes , we can safely state that eqs .
( [ 2 ] ) and ( [ eq : soln ] ) provide us sufficient resource for the investigation of the nonclassical phenomena . to illustrate this point
we may note that the criteria for quadrature squeezing in single - mode ( @xmath42 and compound mode ( @xmath43 are @xcite @xmath44 and @xmath45 where @xmath46 and the quadrature operators are defined as @xmath47 and @xmath48 similarly , the existence of single- and multi - mode nonclassical ( sub - poissonian ) photon statistics can be obtained through the following inequalities@xmath49 and @xmath50 where ( [ eq : antib-1 ] ) provides us the condition for single - mode antibunching and ( [ eq : antib2 ] ) provides us the condition for intermodal antibunching .
many of the early investigations on nonclassicality were limited to the study of squeezing and antibunching only , but with the recent development of quantum computing and quantum information it has become very relevant to study entanglement .
interestingly , there exist a large number of inseparability criteria ( @xcite and references therein ) that can be expressed in terms of expectation values of moments of field operators .
for example , hillery - zubairy criterion i and ii ( hz - i and hz - ii ) @xcite are described as @xmath51 and @xmath52 respectively .
another , interesting criterion of inseparability that can be expressed in terms of moments of creation and annihilation operators is duan et al.s criterion which is described as follows @xcite : @xmath53 where @xmath54 clearly our analytic solution ( [ 2])-([eq : soln ] ) enables us to investigate intermodal entanglement using these criteria and a bunch of other criteria of nonclassicality that are described in ref .
interestingly , all the inseparability criteria described above and in the remaining part of the present work can be viewed as special cases of shchukin - vogel entanglement criterion @xcite . in ref .
@xcite , miranowicz et al .
have explicitly established this point .
so far we have described criteria of nonclassicality that are related to the lowest order nonclassicality .
however , nonclassicality may be witnessed via higher order criterion of nonclassicality , too .
investigations on higher order nonclassical properties of various optical systems have been performed since long .
for example , in late seventies some of the present authors showed the existence of higher order nonclassical photon statistics in different optical systems using criterion based on higher order moments of number operators ( cf .
@xcite and chapter 10 of @xcite and references therein ) .
however , in those early works , higher order antibunching ( hoa ) was not specifically discussed , but existence of higher order nonclassical photon statistics was reported for degenerate and nondegenerate parametric processes in single and compound signal - idler modes , respectively and also for raman scattering in compound stokes - anti - stokes mode up to @xmath55 . criterion for hoa was categorically introduced by c. t. lee @xcite in 1990 . since then hoa is reported in several quantum optical systems ( @xcite and references therein ) and atomic systems @xcite . however , except a recent effort by us @xcite no effort has yet been made to study hoa in optical couplers .
the existence of hoa can be witnessed through a set of equivalent but different criteria , all of which can be interpreted as modified lee criterion . in what follows
we will investigate the existence of hoa in the contradirectional optical coupler of our interest using a simple criterion of @xmath56-th order single - mode antibunching introduced by some of us @xcite @xmath57 here @xmath58 and @xmath59 refer to the usual antibunching and the higher order antibunching , respectively .
similarly , we can also investigate the existence of higher order squeezing , but definition of higher order squeezing is not unique . to be precise
, it is usually investigated using two different criteria @xcite .
the first criterion was introduced by hong and mandel in 1985 @xcite and a second criterion was subsequently introduced by hillery in 1987 @xcite . in hong and mandel criterion @xcite , the reduction of higher order moments of usual quadrature operators with respect to their coherent state counterparts
are considered as higher order squeezing . however , in hillery s criterion , reduction of variance of an amplitude powered quadrature variable for a quantum state with respect to its coherent state counterpart is considered as higher order squeezing . in what follows
we have restricted our study on higher order squeezing to hillery s criterion of amplitude powered squeezing .
specifically , hillery introduced amplitude powered quadrature variables as @xmath60 and @xmath61 it is easy to check that @xmath62 and @xmath63 do not commute and consequently we can obtain an uncertainty relation and thus a criterion of @xmath64 order amplitude squeezing as@xmath65\right\rangle \right|}{\normalcolor < 0}}.}\ ] ] thus , for the specific case , @xmath58 , hillery s criterion for amplitude squared squeezing can be obtained as @xmath66 where @xmath67 similarly , we can obtain specific criteria of amplitude powered squeezing for other values of @xmath68 further , there exists a set of higher order inseparability criteria . to be precise ,
all the criteria that are used for witnessing the existence of multi - partite ( multi - mode ) entanglement are essentially higher order criteria @xcite as they always uncover some higher order correlation .
interestingly , even in bipartite ( two - mode ) case one can introduce operational criterion for detection of higher order entanglement . for example , hillery - zubairy introduced following two criteria of higher order intermodal entanglement @xcite @xmath69 and@xmath70 here @xmath71 and @xmath56 are non - zero positive integers and the lowest possible values of @xmath71 and @xmath56 are @xmath72 which reduces ( [ hoe - criteria-1 ] ) and ( [ hoe - criteria-2 ] ) to usual hz - i criterion ( i.e. , ( [ hz1 ] ) ) and hz - ii criterion ( i.e. , ( [ hz2 ] ) ) , respectively .
thus these two criteria may be viewed as generalized versions of the well known lower order criteria of hillery and zubairy ( i.e. , the criteria described in ( [ hz1 ] ) and ( [ hz2 ] ) ) . however , these generalized criterion can also be obtained as special cases of more general criterion of shchukin and vogel @xcite . for the convenience of the readers , we refer to ( [ hoe - criteria-1 ] ) and ( [ hoe - criteria-2 ] ) as hz - i criterion and hz - ii criterion respectively in analogy to the lowest order cases . in
what follows , a quantum state will be called ( bipartite ) higher order entangled state if it is found to satisfy ( [ hoe - criteria-1 ] ) and/or ( [ hoe - criteria-2 ] ) for any choice of integers @xmath71 and @xmath56 satisfying @xmath73 existence of higher order entanglement can also be viewed through the criteria of multi - partite entanglement .
for example , li et al .
@xcite proved that a three - mode ( tripartite ) quantum state is not bi - separable in the form @xmath74 ( i.e. , compound mode @xmath75 is entangled with the mode @xmath8 ) if the following inequality holds for the three - mode system    @xmath76    where @xmath77 are positive integers and annihilation operators @xmath78 correspond to the three modes .
a quantum state that satisfy ( [ eq : tripartite ent1 ] ) is referred to as @xmath74 entangled state .
the three - mode inseparability criterion mentioned above is not unique .
there exist various alternative criteria of three - mode entanglement .
for example , an alternative criterion for detection of @xmath74 entangled state is @xcite @xmath79 for @xmath80 this criterion coincides with miranowicz et al.s criterion @xcite . using ( [ eq : tripartite ent1 ] ) and ( [ eq : tripartite ent2 ] ) , we can easily obtain criteria for detection of @xmath81 and @xmath82 entangled states and use them to obtain a simple criterion for detection of fully entangled tripartite state . specifically , using ( [ eq : tripartite ent1 ] ) and ( [ eq : tripartite ent2 ] ) respectively we can conclude that the three modes of our interest are not bi - separable in any form if one of the following two sets of inequalities are satisfied simultaneously @xmath83    @xmath84    further , for a fully separable pure state we always have @xmath85^{\frac{1}{2}}.\label{eq : fully ent2}\ ] ] thus , a three - mode pure state that violates ( [ eq : fully ent2 ] ) ( i.e. , satisfies @xmath86 and simultaneously satisfies either ( [ eq : fully enta 0 ] ) or ( [ eq : fully enta 1 ] ) is a fully entangled state as it is neither fully separable nor bi - separable in any form .
spatial evolution of different operators that are relevant for witnessing nonclassicality can be obtained using the perturbative solutions ( [ 2])-([eq : soln ] ) reported here .
for example , using ( [ 2])-([eq : soln ] ) we can obtain the following closed form expressions for number operators of various field modes @xmath87,\end{array}\label{eq : na}\ ] ] @xmath88,\end{array}\label{eq : nb}\ ] ] @xmath89.\label{eq : nb2}\ ] ] it is now straight forward to compute the average values of the number of photons in different modes with respect to a given initial state . in the present work we consider that the initial state is product of three coherent states : @xmath90 where @xmath91 and @xmath92 are eigen kets of annihilation operators @xmath93 and @xmath8 , respectively .
thus , @xmath94 and @xmath95 are the number of input photons in the field modes @xmath93 and @xmath8 , respectively . for a spontaneous process ,
@xmath96 and @xmath97 whereas , for a stimulated process , the complex amplitudes are not necessarily zero and it seems reasonable to consider @xmath98      using ( [ 2])-([eq : soln ] ) and ( [ eq : condition for squeezing])-([eq : two - mode quadrature ] ) we obtain analytic expressions for variance in single - mode and compound mode quadratures as    @xmath99 & = & \frac{1}{4}\left[1\pm\left\ { \left(f_{1}f_{4}+f_{2}f_{3}\right)\gamma+{\rm c.c}.\right\ } \right],\\ \left[\begin{array}{c } \left(\delta x_{b_{1}}\right)^{2}\\ \left(\delta y_{b_{1}}\right)^{2}\end{array}\right ] & = & \frac{1}{4}\left[1\pm\left\ { \left(g_{1}g_{4}+g_{2}g_{3}\right)\gamma+{\rm c.c}.\right\ } \right],\\ \left[\begin{array}{c } \left(\delta x_{b_{2}}\right)^{2}\\ \left(\delta y_{b_{2}}\right)^{2}\end{array}\right ] & = & \frac{1}{4},\end{array}\label{eq : squeezing}\ ] ]    and    @xmath100 & = & \frac{1}{4}\left[1\pm\frac{1}{2}\left\ { \left(\left(f_{1}+g_{1}\right)\left(f_{4}+g_{4}\right)+\left(f_{2}+g_{2}\right)\left(f_{3}+g_{3}\right)\right)\gamma+{\rm c.c}.\right\ } \right],\\ \left[\begin{array}{c } \left(\delta x_{ab_{2}}\right)^{2}\\ \left(\delta y_{ab_{2}}\right)^{2}\end{array}\right ] & = & \frac{1}{4}\left[1\pm\frac{1}{2}\left\ { \left(f_{1}f_{4}+f_{2}f_{3}\right)\gamma+{\rm c.c}.\right\ } \right]=\frac{1}{2}\left[\begin{array}{c } \left(\delta x_{a}\right)^{2}\\ \left(\delta y_{a}\right)^{2}\end{array}\right]+\frac{1}{8},\\ \left[\begin{array}{c } \left(\delta x_{b_{1}b_{2}}\right)^{2}\\ \left(\delta y_{b_{1}b_{2}}\right)^{2}\end{array}\right ] & = & \frac{1}{4}\left[1\pm\frac{1}{2}\left\ { \left(g_{1}g_{4}+g_{2}g_{3}\right)\gamma+{\rm c.c}.\right\ } \right]=\frac{1}{2}\left[\begin{array}{c } \left(\delta x_{b_{1}}\right)^{2}\\ \left(\delta y_{b_{1}}\right)^{2}\end{array}\right]+\frac{1}{8},\end{array}\label{eq : intermodal squeezing}\ ] ]    respectively . from eq .
( [ eq : squeezing ] ) it is clear that no squeezing is observed in @xmath8 mode .
however , squeezing is possible in @xmath28 mode and @xmath101 mode as illustrated in fig .
[ fig : squeezing - del k - g ] a - b , d - e . further , we observed intermodal squeezing in quadratures @xmath102 and @xmath103 by plotting right hand sides of eq .
( [ eq : intermodal squeezing ] ) in fig .
[ fig : squeezing - del k - g ] c and fig .
[ fig : squeezing - del k - g ] f. variation of amount of squeezing in different modes with phase mismatch @xmath9 and nonlinear coupling constant @xmath104 are shown in fig .
[ fig : squeezing - del k - g ] a - c and fig .
[ fig : squeezing - del k - g ] d - f , respectively .
we have also studied the effect of linear coupling constant @xmath20 on the amount of squeezing , but its effect is negligible in all other quadratures except the quadratures of compound mode @xmath105 in compound mode @xmath75 , after a short distance depth of squeezing is observed to increase with the decrease in the linear coupling constant @xmath20 ( this is not illustrated through figure ) .
intermodal squeezing in compound mode quadratures @xmath106 and @xmath107 can be visualized from the last two rows of eq.([eq : intermodal squeezing ] ) .
specifically , we can see that variance in compound mode quadrature @xmath108 and @xmath109 have _ bijective _ ( both one - to - one and onto ) correspondence with the variance in @xmath110 and @xmath111 , respectively , where @xmath112 to be precise , quadrature squeezing in single - mode @xmath110(@xmath111 ) implies quadrature squeezing in @xmath113 and vice versa .
for example , @xmath114 or , @xmath115 this is why we have not explicitly shown the variance of compound mode quadratures @xmath116 and @xmath117 with different parameters as we have done for the other cases . as we have shown squeezing in @xmath118 and @xmath119 through fig .
[ fig : squeezing - del k - g ] a - b , and c - d this implies the existence of squeezing in quadrature @xmath120 and @xmath121 thus we have observed intermodal squeezing in compound modes involving @xmath8 .
this nonclassical feature was not observed in earlier studies @xcite as in those studies @xmath8 mode was considered classical .
the plots do not show quadrature squeezing in @xmath122 @xmath123 , and @xmath102 ( for some specific values of @xmath104 ) .
however , a suitable choice of phase of the input coherent state would lead to squeezing in these quadratures .
for example , if we replace @xmath124 by @xmath125 ( i.e. , if we chose @xmath126 instead of present choice of @xmath127 ) then we would observe squeezing in all these quadratures , but the squeezing that is observed now with the original choice of @xmath124 would vanish .
this is so as all the expressions of variance of quadrature variables that are @xmath128 have a common functional form : @xmath129 ( c.f . , eqs .
( [ eq : squeezing ] ) and ( [ eq : intermodal squeezing ] ) ) .      after establishing the existence of squeezing in single - modes and compound modes
we now examine the possibility of higher order squeezing using eqs .
( [ 2])-([eq : soln ] ) , ( [ eq : na])-([eq : nb2 ] ) and ( [ eq : criterion - amplitude squared ] ) .
in fact , we obtain @xmath130 & = & \pm\frac{n^{2}}{4}\left[\gamma\left(f_{1}f_{4}+f_{2}f_{3}\right)\left(f_{1}\alpha+f_{2}\beta\right)^{2n-2}+{\rm c.c.}\right],\end{array}\label{eq : asq-1}\ ] ]    @xmath131 & = & \pm\frac{n^{2}}{4}\left[\gamma\left(g_{1}g_{4}+g_{2}g_{3}\right)\left(g_{1}\alpha+g_{2}\beta\right)^{2n-2}+{\rm c.c.}\right],\end{array}\label{eq : as - q-2}\ ] ]    and @xmath132 & = & 0.\end{array}\label{eq : as - q-3}\ ] ] thus , we do not get any signature of amplitude powered squeezing in @xmath8 mode using the present solution . in contrary , mode @xmath28 @xmath133 is found to show amplitude powered squeezing in one of the quadrature variables for any value of interaction length as both @xmath134 and @xmath135 ( @xmath136 and @xmath137 ) can not be positive simultaneously . to study the possibilities of amplitude powered squeezing in further detail we have plotted the spatial variation of @xmath138 and @xmath139 in fig .
[ fig : amplitude - squared - squeezing ] . negative regions of these two plots clearly illustrate the existence of amplitude powered squeezing in both @xmath28 and @xmath101 modes for @xmath58 and @xmath140 . extending our observations in context of single - mode squeezing and intermodal squeezing
we can state that the appearance of amplitude powered squeezing in a particular quadrature can be controlled by suitable choice of phase of input coherent state @xmath124 as the expressions for amplitude powered squeezing reported in ( [ eq : asq-1 ] ) and ( [ eq : as - q-2 ] ) have a common functional form @xmath141 .      the condition of hoa is already provided through the inequality ( [ hoa-1 ] ) . now using this inequality along with eqns .
( [ 2])-([eq : soln ] ) and ( [ eq : na])-([eq : nb2 ] ) we can obtain closed form analytic expressions for @xmath142 for various modes as follows @xmath143 @xmath144    @xmath145    further , using the condition of intermodal antibunching described in ( [ eq : antib2 ] ) and eqns .
( [ 2])-([eq : soln ] ) we obtain following closed form expressions of @xmath146 @xmath147 @xmath148    @xmath149    from the above expressions it is clear that neither the single - mode antibunching nor the intermodal antibunching is obtained involving @xmath8 mode .
as the expressions obtained in the right hand sides of ( [ eq : dan ] ) , ( [ eq : dbn ] ) and ( [ eq : dab1 ] ) are not simple , we plot them to investigate the existence of single - mode and compound mode antibunching .
the plots for usual antibunching and intermodal antibunching are shown in fig .
[ fig : antibunching ] .
existence of antibunching is obtained in single - mode @xmath28 for @xmath127 and the same is illustrated through fig .
[ fig : antibunching ] a. however , in a effort to obtain antibunching in @xmath101 mode , we do not observe any antibunching in @xmath101 mode for @xmath127 .
interestingly , from ( [ eq : dbn ] ) it is clear that if we replace @xmath127 by @xmath150 as before and keep @xmath151 unchanged , then we would observe antibunching for all values of rescaled interaction length @xmath152 this is true in general for all values of @xmath124 . to be precise
, if we observe antibunching ( bunching ) in a mode for @xmath153 we will always observe bunching ( antibunching ) for @xmath154 if we keep @xmath151 unchanged .
this fact is illustrated through fig .
[ fig : antibunching ] b where we plot variation of @xmath155 with @xmath156 and have observed the existence of antibunching . in compound mode
@xmath157 we can observe existence of antibunching for both @xmath127 and @xmath158 .
however , in fig .
[ fig : antibunching ] c we have illustrated the existence of intermodal antibunching in compound mode @xmath75 by plotting variation @xmath159 with rescaled interaction length @xmath156 for @xmath158 as region of nonclassicality is relatively larger ( compared to the case where @xmath127 and @xmath151 are same ) in this case .
now we may extend the discussion to hoa and plot right hand sides of ( [ eq : dan ] ) and ( [ eq : dbn ] ) for various values of @xmath56 .
the plots are shown in fig .
[ fig : hoa ] which clearly illustrates the existence of hoa and also demonstrate that the depth of nonclassicality increases with @xmath68 this is consistent with earlier observations on hoa in other systems @xcite .
we first examine the existence of intermodal entanglement in compound mode @xmath75 using hz - i criterion ( [ hz1 ] ) .
to do so we use eqns .
( [ 2])-([eq : soln ] ) and ( [ eq : na])-([eq : nb2 ] ) and obtain@xmath160 similarly , using hz - ii criterion ( [ hz2 ] ) we obtain    @xmath161.\end{array}\label{eq : hz2-ab1}\ ] ]    it is easy to observe that eqns .
( [ eq : hz-1-ab1 ] ) and ( [ eq : hz2-ab1 ] ) provide us the following simple relation that is valid for the present case : @xmath162 , which implies that for any particular choice of rescaled interaction length @xmath156 either hz - i criterion or hz - ii criterion would show the existence of entanglement in contradirectional asymmetric nonlinear optical coupler as both of them can not be simultaneously positive .
thus the compound mode @xmath75 is always entangled .
the same is explicitly illustrated through fig .
[ fig : ent - hz1-ii ] .
similar investigations using hz - i and hz - ii criteria in the other two compound modes ( i.e. , @xmath163 and @xmath164 failed to obtain any signature of entanglement in these cases .
further , signature of intermodal entanglement was not witnessed using duan et al .
criterion as using the present solution and ( [ duan-1 ] ) we obtain @xmath165 however , it does not ensure separability of these modes as hz - i , hz - ii and duan et al .
inseparability criteria are only sufficient and not necessary .
( color online ) hillery - zubairy criterion i ( solid line ) and criterion ii ( dashed line ) for entanglement are showing intermodal entanglement between modes @xmath28 and @xmath166 here @xmath167 ( solid line ) and @xmath168 ( dashed line ) are plotted with rescaled interaction length @xmath169 for mode @xmath75 with @xmath170 for @xmath171 ( thin blue lines ) and @xmath172 ( thick red lines ) . ]
we may now study the possibilities of existence of higher order entanglement using eqns .
( [ hoe - criteria-1])-([eq : fully ent2 ] ) .
to begin with , we use ( [ 2])-([eq : soln ] ) and ( [ hoe - criteria-1 ] ) to yield @xmath173 similarly , using ( [ 2])-([eq : soln ] ) and ( [ hoe - criteria-2 ] ) we can produce an analytic expression for @xmath174 and observe that @xmath175 from relation ( [ eq : emnprimeab1 ] ) , it is clear that the higher order entanglement between @xmath28 mode and @xmath101 mode would always exist for any choice of @xmath176 @xmath71 and @xmath56 .
this is so because @xmath177 and @xmath174 can not be simultaneously positive . using ( [ eq : emnab1 ] ) and
( [ eq : emnprimeab1 ] ) , it is a straight forward exercise to obtain analytic expressions of @xmath177 and @xmath174 for specific values of @xmath71 and @xmath68 such analytic expressions are not reported here as the existence of higher order entanglement is clearly observed through ( [ eq : emnprimeab1 ] ) .
however , in fig .
[ fig : hoe ] we have illustrated the variation of @xmath178 and @xmath179 with the rescaled interaction length , @xmath156 .
negative parts of the plots shown in fig .
[ fig : hoe ] illustrate the existence of higher order intermodal entanglement in compound mode @xmath105 as expected from ( [ eq : emnprimeab1 ] ) , we observe that for any value of @xmath156 compound mode @xmath75 is higher order entangled .
further , it is observed that hillery - zubairy s higher order entanglement criteria ( [ hoe - criteria-1])-([hoe - criteria-2 ] ) fail to detect any signature of higher order entanglement in compound modes @xmath163 and @xmath180    one can also investigate higher order entanglement using criterion of multi - partite ( multi - mode ) entanglement as all the multi - mode entangled states are essentially higher order entangled .
here we have only three modes in the coupler and thus we can study higher order entanglement by investigating the existence of three - mode entanglement . a three - mode pure state that violates ( [ eq : fully ent2 ] ) ( i.e. , satisfies @xmath86 and simultaneously satisfies either ( [ eq : fully enta 0 ] ) or ( [ eq : fully enta 1 ] ) is a fully entangled state . using ( [ 2])-([eq : soln ] ) and ( [ eq : tripartite ent1])-([eq : fully ent2 ] ) we obtain following set of interesting relations for @xmath181 :    @xmath182    @xmath183    and @xmath184 from ( [ eq : relation1 ] ) , it is easy to observe that three modes of the coupler are not bi - separable in the form @xmath81 and @xmath185 for any choice of rescaled interaction length @xmath186 further , eqn .
( [ eq : relation3 ] ) shows that the three modes of the coupler are not fully separable for @xmath187 ( c.f .
positive regions of plot of @xmath167 shown in fig .
[ fig : ent - hz1-ii ] ) .
however , ( [ eq : relation2 ] ) illustrate that the present solution does not show entanglement between coupled mode @xmath75 and single - mode @xmath188 thus the three modes present here are not found to be fully entangled . specifically , three - mode ( higher order ) entanglement is observed here , but signature of fully entangled three - mode state is not observed .
further , we have observed that in all the figures depth of nonclassicality increases with @xmath189    ( color online ) higher order entanglement is observed using hillery - zubairy criteria .
solid lines show spatial variation of @xmath178 and dashed lines show spatial variation of @xmath179 with @xmath170 for @xmath171 ( thin blue lines ) and @xmath172 ( thick red lines ) .
it is observed that depth of nonclassicality increases with increase in @xmath190 . ]    in the present paper , various lower order and higher order nonclassical phenomena have been observed in contradirectional asymmetric nonlinear optical coupler .
however , so far we have discussed only the stimulated cases as no nonclassical phenomenon is expected to be observed in spontaneous case .
this is so because all the useful non - vanishing expressions for witnessing nonclassicality ( i.e. , eqs .
( [ eq : asq-1])-([eq : relation3 ] ) ) are proportional to @xmath191 thus all these expressions would vanish for @xmath192 .
approximation @xcite & present work +   1 .
& squeezing & @xmath8 & not investigated & not observed +   2 . &
intermodal squeezing & @xmath193 & not observed & observed +   3 . &
amplitude squared squeezing & @xmath93 & not investigated & observed +   4 . & amplitude squared squeezing & @xmath8 & not investigated & not observed +   5 . & lower order and higher order intermodal entanglement & @xmath75 & not investigated & observed +   6 . & lower order and higher order intermodal entanglement & @xmath193 & not investigated & not observed +   7 . &
three - mode ( higher order ) bi - separable entanglement & @xmath81 , @xmath185 & not investigated & observed +   8 . &
three - mode ( higher order ) bi - separable entanglement & @xmath81 & not investigated & not observed +
in the present study we report lower order and higher order nonclassicalities in a contradirectional asymmetric nonlinear optical coupler using a set of criteria of entanglement , single - mode squeezing , intermodal squeezing , antibunching , intermodal antibunching etc .
variation of nonclassicality with various parameters , such as number of input photon in the linear mode , linear coupling constant , nonlinear coupling constant and phase mismatch is also studied and it is observed that amount of nonclassicality can be controlled by controlling these parameters .
the contradirectional asymmetric nonlinear optical coupler studied in the present work was studied earlier using a short - length solution and considering @xmath8 mode as classical @xcite . in contrast , a completely quantum mechanical solution of the equations of motion is obtained here using sen - mandal approach which is not restricted by length .
the use of better solution and completely quantum mechanical treatment led to the identification of several nonclassical characters of a contradirectional asymmetric nonlinear optical coupler that were not reported in earlier studies .
all such nonclassical phenomena that are observed here and were not observed in earlier studies are listed in table [ tab : nonclassicalities - observed - in ] .
further , there exist a large number of nonclassicality criteria that are not studied here and are based on expectation values of moments of annihilation and creation operators ( c.f .
table i and ii of ref .
as we already have compact expressions for the field operators it is a straight forward exercise to extend the present work to investigate other signatures of nonclassicality , such as , photon hyperbunching @xcite , sum and difference squeezing of an - tinh and hillery , inseparability criterion of manicini et al .
@xcite , simon @xcite and miranowicz et al .
@xcite etc .
further , the present work can be extended to investigate the nonclassical phenomena in other types of contradirectional optical couplers ( e.g. , contradirectional parametric coupler , contradirectional raman coupler , etc . ) that are either not studied till date or studied using short - length solution .
recently , allevi et al .
@xcite and avenhaus et al .
@xcite have independently reported that they have experimentally measured @xmath194 which is sufficient to completely characterize bipartite multi - mode states .
it is easy to observe that ability to experimentally measure @xmath194 ensures that we can experimentally detect signatures of nonclassicalities reported here .
it is true that most of the nonclassicalities reported here can be observed in some other bosonic systems , too .
however , the present system has some intrinsic advantages over most of the other systems as it can be used as a component in the integrated waveguide based structures in general and photonic circuits in particular [ @xcite and references therein ] .
thus the nonclassicalities reported in this easily implementable waveguide based system is expected to be observed experimentally .
further , the system studied here is expected to play important role as a source of nonclassical fields in the integrated waveguide based structures .
* acknowledgment : * k. t. and a. p. thank the department of science and technology ( dst ) , india , for support provided through dst project no .
sr / s2/lop-0012/2010 and a. p. also thanks the operational program education for competitiveness - european social fund project cz.1.07/2.3.00/20.0017 of the ministry of education , youth and sports of the czech republic .
j. p. thanks the operational program research and development for innovations - european regional development fund project cz.1.05/2.1.00/03.0058 of the ministry of education , youth and sports of the czech republic .
j. lugani , studies on guided wave devices and ultracold atoms in optical cavity for applications in quantum optics and quantum information  , phd thesis , iit delhi ( 2013 ) , http://web.iitd.ac.in/~sankalpa/jasleen_corrected_thesis.pdf .