nowadays cross - kerr nonlinearity is considered as a promising tool for quantum computation and non - classical state generation @xcite .
an entanglement arising between modes participating in the cross - kerr interaction can be used for establishing an interface between matter qubits and ` flying ' light photonic @xcite , for generation of non - gaussian states and schrdinger - cat states , and for performing quantum gate operations @xcite .
an interest to cross - kerr nonlinear interactions is heated up by both developing ways to implement effectively even very weak non - linearities ( which one commonly expects to have in practice ) , and by discovery of methods to produce sufficiently large self - kerr and cross - kerr nonlinearities ( such as implementation of the electromagnetically induced transparency ( eit ) and photonic crystals @xcite ) .
decoherence is a main practical obstacle to implementations of schemes using kerr and cross - kerr nonlinearities .
a genuine example of quantum state degradation due to losses is a decoherence of a quantum superposition state , the effect drastically enhanced if speaking of a macroscopic superposition state of the schrdinger cat type . already more than 20 years ago
the self - kerr nonlinearity was proposed as a tool for generating such a schrdinger - cat state ( more precisely , a superposition of two coherent states with the same amplitudes but opposite phases ) @xcite . however , photon losses turn this superposition into statistical mixture of two coherent states with the rate proportional to the square modulus of the amplitude of these states
. modifications of the scheme for the cross - kerr nonlinearities or four - wave mixing brought no advantage with respect to photon loss @xcite .
this unfortunate circumstance made one look for the ways to circumvent the problem of decoherence that inevitably accompanies kerr nonlinearity .
recent suggestions in this direction are based on the conditional preparation of desired states ( which brings into consideration an additional problem of the finite detection efficiency ) , and are aimed to exploit weak nonlinearities .
recently , even a way to produce cat - states ` on demand ' was suggested using a source of single - photons @xcite .    in our work
we want to discuss an aspect of the decoherence which has been seldom discussed when considering an influence of losses on states generated via kerr nonlinearity .
namely , we address losses arising _ in the process _ of generation and not due to propagation of the generated state via lossy channels .
we concentrate our attention on a feature that might be quite significantly pronounced in modern schemes of generating large kerr nonlinearity : the modal loss can be strongly correlated .
indeed , the modes occupy the same volume and interact with the same physical systems which form the reservoirs .
also , if the kerr - nonlinearity scheme implies a sufficiently strong dispersive coupling of light modes to emitters , then coupling of these emitters to dissipative reservoirs might also appear to be quite strong . as a result , this would mean strongly correlated modal losses .
for example , in photonic crystals high density of states in the vicinity of a modal frequencies and emitter s transition frequency can cause the strong emitter - field coupling ; but it would also imply higher population loss of emitters due to coupling to radiative reservoirs . dephasing losses of emitters would as well invoke a correlated modal dephasing .
coupling to correlated reservoirs can drastically change state dynamics in comparison with loss to uncorrelated reservoirs .
for example , it was demonstrated that coupling to the common reservoir preserves entanglement of a two - mode state @xcite .
moreover , coupling to the common reservoir is capable of creating an entanglement between states of initially unentangled modes even in absence of any direct interaction between them .    in our work
we demonstrate both how the correlated loss arises via kerr nonlinear process , and how it affects the generated states . for this purpose
we derive analytic solution generalizing a powerful and illustrative method of chaturvedi and srinivasan @xcite . on a number of examples
we show how the correlated loss enhances and creates intermodal correlations and even entanglement , and might lead to generation of entangled states quite different from those generated in the same scheme without loss .
correlated loss can result in the significantly enhanced robustness of the generation scheme .
the outline of the paper is as follows . in sec .
ii and in the related appendices we describe how cross - correlation terms emerge via correlations of markovian reservoirs ; we consider an example of the emitter - field interaction schemes producing correlated modal losses in the section iii . then in sec .
iv we describe the method for obtaining exact solutions of the cross - kerr nonlinear interaction between modes in presence of losses to uncorrelated reservoirs and give generalization of the method for some cases of correlated losses . in the sec .
v we analyze influence of losses in the nonclassical state generation process for the case of uncorrelated loss .
some examples of correlated losses are considered in sec .
to illustrate clearly an influence of correlated and uncorrelated losses , we restrict ourselves to the markovian loss accounting for photon losses and dephasing of interacting modes .
we start from the general effective hamiltonian @xmath0 describing both self- and cross - interaction ( for the moment we refrain from detailing it ) and interaction of modes with reservoirs responsible for losses @xmath1 : @xmath2 here we use the interaction picture with respect to the free hamiltonians of reservoirs and the modes participating in the interaction process .
these modes are described by usual bosonic creation and annihilation operators satisfying @xmath3=[a_2,a_2^{\dagger}]=1 , \quad [ a_1,a_2^{\dagger}]=[a_1,a_2]=0.\ ] ] the operators @xmath4 and @xmath5 describe reservoirs responsible , correspondingly , to the photon losses in modes @xmath6 and @xmath7 , and to dephasing of these modes .
they may include also stochastic variables describing different realizations of reservoirs .
it should be emphasized that reservoir operators @xmath4 and @xmath5 ( together with the initial state of the reservoir ) completely describe the reservoir properties with respect to the interaction with the modes .
these operators are built on the basis of underlying microscopic model and account for all relevant physical parameters .
for example , if the photon loss reservoir of the first mode is composed of electromagnetic field modes with frequencies @xmath8 , described by the creation and annihilation operators @xmath9 , then @xmath10 where each @xmath11 is the constant of interaction of the mode @xmath6 with the @xmath12th mode of the reservoir .
sets of frequencies @xmath8 and interaction constants @xmath11 describe completely physical properties of the reservoir . in particular , if the reservoir is the set of electromagnetic modes of a non - absorbing dielectric structure , they are found from the eigensolutions of maxwell s equations for this structure @xcite .
we make a number of standard assumptions in considering decoherence : coupling with the reservoirs is weak , the initial state of the modes @xmath13 are uncorrelated with initial state of reservoirs , and correlation times of reservoirs are small enough to enable an implementation of the born - markov approximation .
then , using , for example , a time - convolutionless projection operator technique @xcite , one can obtain the master equation , the liouville equation of the following form : @xmath14 + { 1\over 2 } [ ( \gamma_1-\gamma_{12})\mathcal{l}(a_1)+ ( \gamma_2-\gamma_{12})\mathcal{l}(a_2)+ ( d_1-d_{12})\mathcal{l}(a_1^{\dagger}a_1)+(d_2-d_{12})\mathcal{l}(a_2^{\dagger}a_2)]\rho(t)+\\ { 1\over 2 } [ \gamma_{12}\mathcal{l}(a_1+a_2)+d_{12}\mathcal{l}(a_1^{\dagger}a_1+a_2^{\dagger}a_2)]\rho(t ) , \label{master_markov}\end{aligned}\ ] ] where the superoperator @xmath15 acts on the density matrix as @xmath16 here we set @xmath17 for simplicity . for more details of the derivation see appendix a. by construction , the master equation ( [ master_markov ] ) provides for non - negative definite @xmath18 for arbitrary @xmath19 .
note that we call reservoirs `` correlated '' , if their integrated cross - correlation function is non - zero , for example , @xmath20 and , correspondingly , the coefficient @xmath21 defined in eq .
( [ gamma12 ] ) is non - zero . of course
, one can always transform eq .
( [ master_markov ] ) to the diagonal form .
however , in this case lindblad operators ( i.e. operators like @xmath22 in the diagonal form ( [ diagonal form ] ) ) will be the linear superpositions of the former lindblad operators , and the transformed equation will be still describing a coupling between physical objects represented by these original lindblad operators ( modes @xmath23 in our case ) .
for example , the possibility that both modes are coupled to the same reservoir ( either the photon loss reservoir or the dephasing one ) corresponds to an equality in relations ( [ rates_unequalities ] ) in appendix a. in this case the operators describing the reservoir are proportional to each other , say , @xmath24 and @xmath25 . equation ( [ master_markov ] ) then reduces to @xmath14 + { 1\over 2 } [ \gamma_{1}\mathcal{l}(a_1+xa_2)+d_{1}\mathcal{l}(a_1^{\dagger}a_1+ya_2^{\dagger}a_2)]\rho(t ) .
\label{master_markov_corr}\end{aligned}\ ] ] the most important point here is that in this equation both modes behave like a single object with respect to relaxation .
such a ` decoherence ' is able to induce enduring entanglement between modes @xmath6 and @xmath7 ( an example is shown in the appendix b ) .
furthermore , interacting light modes might experience correlated losses due to fact that they both interact with the same atoms .
one can illustrate the mechanism of such a loss by the following qualitative consideration .
let us consider a general hamiltonian describing the light - atom interaction @xmath26 plus terms describing interaction of atoms with the dissipative reservoirs @xmath27 where @xmath28 and @xmath29 are operators describing the atoms and dissipative reservoirs , respectively .
the hamiltonian of this form describes a general interaction between field modes and emitters of some type as the correlated loss can occur in different physical systems .
provided atom - reservoir interactions in ( [ ham - loss ] ) is sufficiently weak as not to perturb much interaction between light modes and atoms , the following approximation can be used for the atomic operators : @xmath30 where @xmath31 , and t denotes the time - ordering operator . after averaging out the atomic variable ,
the terms that describe reservoir - mode coupling in the effective interaction hamiltonian will take the form of @xmath32 . here
@xmath33 and @xmath34 denotes the averaging over atomic states ( in general , over emitter s states ) .
thus , one can see that coupling of light modes to the same atom ( emitter ) interacting with the dissipative reservoir under the condition of adiabatic elimination of emitter s variables leads directly to the mode - reservoir interaction terms in the resulting effective hamiltonian .
note , that these terms remain linear in reservoir operators @xmath29 .
hence one can derive a master equation averaging over the reservoir in a standard way . in the appendix
b we give examples of the derivation of effective hamiltonians and the corresponding master equations discussing the simplest two - mode jaynes - cummings system ( fig .
[ fig1](a ) ) .
the described scenario of how correlated loss emerge is quite general , and can take place for a wide range of schemes involving light - shift - induced photonic nonlinearities .
a three - level @xmath35-system interacting with classical driving and quantum fields represent a more realistic example of the scheme with the correlated loss .
[ lambda - system ] consider the large cross - kerr nonlinearity generation suggested in ref .
@xcite and depicted in fig .
[ fig1](b ) .
there two modes are coupled to the transition between 1 and 3 levels of the @xmath35-system in presence of two classical driving fields on transitions 1 - 2 and 1 - 3 . in the rotating - wave approximation and in the interaction picture with respect to free hamiltonian of the reservoir ,
the problem is described by the following interaction hamiltonian : @xmath36\sigma_{13}+\omega\exp\{{{\rm i}}\delta_2t\}\sigma_{23}+\mbox{h . c. } , \label{hamiltonian_lambda_system}\end{aligned}\ ] ] where @xmath37 is the rabi frequency of the driving fields ; @xmath38 , @xmath39 ; and @xmath40 are the interaction constants for the coupling of a light mode to an emitter ( atom ) . here for simplicity we have taken into account only losses on the transition 1 - 3 .
the setup depicted in fig .
[ fig1](b ) can be realized in toroidal microcavities , where @xmath6 and @xmath7 correspond to the clockwise and counter - clockwise propagating modes @xcite .
as usually , we assume the markovian reservoir and the following conditions hold : @xmath41 from the hamiltonian ( [ hamiltonian_lambda_system ] ) , the dynamics is governed by @xmath42(\sigma_{33}-\sigma_{11})+{{\rm i}}\omega\exp\{{{\rm i}}\delta_2t\}\sigma_{21}. \label{31equation_lambda_system}\end{aligned}\ ] ] using the approach described in ref .
@xcite , we assume that level 3 remains practically unpopulated , @xmath43 and @xmath44 , as well as @xmath45 .
thus the @xmath35-system is prepared in the superposition of the metastable levels 1 and 2 ( namely , in the state @xmath46 ) .
then , neglecting small and rapidly oscillating terms , one obtains from eqs.([hamiltonian_lambda_system ] , [ 31equation_lambda_system ] ) the following master equation @xmath47+{1\over2}\left[(\gamma_1-\gamma_{12 } ) \mathcal{l}(a_1)+(\gamma_2-\gamma_{12})\mathcal{l}(a_2)\right]\rho(t)+ \\ { \gamma_{12}\over2}\mathcal{l}(a_1+a_2)\rho(t ) , \label{master_lambda_model}\end{aligned}\ ] ] where @xmath48    the master equation ( [ master_lambda_model ] ) describes the cross - kerr and self - kerr interactions of two light modes plus their coupling to the correlated reservoirs . as follows from the model depicted in fig .
[ fig1](b ) ) , it is actually the same reservoir : one can see that @xmath49 .
naturally , loss rates @xmath50 are much less than the loss rate of the emitter .
however , if one deals with input modes in a coherent state of rather large amplitude for generating large cross - kerr nonlinearities , then even comparatively small losses can strongly influence mode dynamics .
below we consider examples of such an influence , e.g. , an example of a single - mode schdinger - cat state .    the occurrence of a correlated modal loss due to presence of the emitter ( atom in our case ) has been noticed in ref .
they have also pointed out that the loss rate should be proportional to the square of the ratio of the mode - emitter interaction constant and the detuning .
however , we would like to emphasize that this loss rate does not depend on the population of level 3 , as it can be seen from eq .
( [ 31equation_lambda_system ] ) .
now let us turn to the specific nonlinear interaction described by the same type of the hamiltonian as in the examples above and in the appendix b. we will consider the effective hamiltonian in the master equation  ( [ master_markov ] ) @xmath51 in the following general form @xmath52 it describes the cross - kerr and self - kerr interaction with nonlinear coefficients @xmath53 of two modes ( or mode superpositions ) @xmath6 and @xmath7 . to solve the problem described by the master equation ( [ master_markov ] ) with the hamiltonian ( [ h_0_kerr ] ) , we adopt a simple and illustrative ` thermofield ' notation @xcite .
essentially , instead of a density matrix acting on some space @xmath54 , say , @xmath55 , where @xmath56 is the fock state with @xmath57 photons in @xmath54 , we consider a state vector @xmath58 in an extended space @xmath59 , where @xmath60 is the fock state with @xmath61 photons in @xmath62 .
so when the mode operators @xmath63 and @xmath64 in the master equation ( [ master_markov ] ) act on the density matrix from the left , one introduces operators @xmath65 and @xmath66 is such a manner that @xmath67 obviously , operators @xmath63 and @xmath64 commute with @xmath68 and @xmath66 .
an action of the superoperator @xmath69 on the density matrix can be represented in the thermofield notation as @xmath70 also , the commutator of any function of the operators @xmath63 and @xmath64 , e.g. , the hamiltonian @xmath71 is re - written as : @xmath72\longrightarrow \mathbf{h}_0(t)|\rho\rangle=\left(h_0(t;a , a^{\dagger } )     -h_0(t;{\tilde a},{\tilde a}^{\dagger})\right)|\rho\rangle.\ ] ] with help of these notations eq .
( [ master_markov ] ) can be represented in the ` hamiltonian ' form as @xmath73 its solution is then of the form @xmath74    the advantage of using thermofield notation over more traditional algebraic manipulation with superoperators is that in many situations ( and , particularly , ones of our interest ) it enables to simplify , make more illustrative and less cumbersome finding the solution ( [ general solution ] ) and estimation of time - dependent matrix elements . in particular , it allows to represent in a simple form a factorization of the superoperator @xmath75 into multipliers with easily estimated actions on the number states @xcite .    to illustrate this ,
let us consider a simple problem of modal loss in a single mode described by the equation @xmath76 the key to solving this equation lies in the observation that the operators @xmath77 generate the su(1,1 ) algebra with the casimir invariant @xmath78 .
using the disentangling theorem @xcite for this group , we arrive at the simple result @xmath79    \exp[(1-{{\rm e}}^{-\gamma_at})a_-]|\rho(0)\rangle .
\label{solution example}\end{aligned}\ ] ] generally , the dynamics described by the solution ( [ solution example ] ) leads to transforming initially pure states into mixtures .
however , for a coherent initial state with the amplitude @xmath80 , @xmath81 equation  ( [ solution example ] ) gives the following result : @xmath82 returning now to the effective hamiltonian ( [ h_0_kerr ] ) that describes cross - kerr and self - kerr interaction of two modes ( or modal superpositions ) @xmath6 and @xmath7 , we write it in the thermofield notation as follows : @xmath83 operators in eq .
( [ h_0_kerr_general ] ) are @xmath84 , @xmath85 .    for completely uncorrelated reservoirs of different modes ( i.e. @xmath86 ) the master equation ( [ master_markov_term ] ) with the hamiltonian
@xmath26 given by eq .
( [ h_0_kerr_general ] ) can be solved exactly using the following factorization @xcite : @xmath87t\right\ }
\\ \times    \exp[({{\rm i}}p_1-\gamma_1)a_3^{(1)}t ] \exp[({{\rm i}}p_2-\gamma_2)a_3^{(2)}t ]    \exp[\gamma_{1-}a_-^{(1 ) } ] \exp[\gamma_{2-}a_-^{(2 ) } ] \label{ham}\end{aligned}\ ] ] where superoperators @xmath88 are defined similarly as in eq .
( [ algebra definition ] ) and @xmath89}{{{\rm i}}p_k-\gamma_k } , \quad p_k=2\sum\limits_{l=1}^2\chi_{kl}a_0^{(l)}.\end{aligned}\ ] ] it is useful to note that in eq .
( [ ham ] ) all multipliers apart from two last ones are diagonal in the number - state basis . also , it is easy to see that @xmath90=0 $ ] , operators @xmath91 are diagonal in the number - state basis , and operators @xmath92 are simply products of annihilation operators .
thus , eq . ( [ ham ] ) provides for simple analytic solutions both for coherent initial states of interacting modes .
the solution ( [ ham ] ) for uncorrelated reservoirs can be straightforwardly generalized for some special cases of correlated reservoirs and the hamiltonian ( [ h_0_kerr ] ) .
for example , let us consider the problem without dephasing , @xmath93 , and introduce the rotated mode operators @xmath94 as @xmath95 where @xmath96 .
then for the non - unitary part of the master equation ( [ master_markov_term ] ) one has @xmath97 where @xmath98 , @xmath99 .
obviously , if the transformation ( [ rotation ] ) leaves the form of the hamiltonian ( [ h_0_kerr ] ) invariant , one can derive an exact solution in the way described in this section .
in this section we consider specific effects of the uncorrelated losses in the process of the cross - kerr nonlinear interaction .
notably , the losses in such a nonlinear process can lead to loss - mediated correlations between the modes , as seen from eq .
( [ ham ] ) .
these intermodal correlations modify the effect of losses on the quantum state generated in the cross - kerr interaction with respect to the result , which one would intuitively expect treating the generation and loss separately , i.e. subjecting to loss a state that has been produced without losses .
consider a particular problem of generating an entangled two - mode state from initially uncorrelated coherent states .
we assume that the nonlinearity is given purely by the cross - kerr interaction , i.e. we put @xmath100 in the hamiltonian ( [ h_0_kerr ] ) . producing entangled states this way is important in a number of schemes of quantum computation and communication using continuous variables .
we assume modes @xmath6 and @xmath7 to be initially in coherent states with amplitudes @xmath101 and @xmath102 , respectively . as was pointed in the previous section , for this choice of initial states , the solution given by eq .
( [ ham ] ) has a simple form : @xmath103t}2\right]\exp\left\{{{\rm i}}\chi t    ( a_1^\dagger a_1a_2^\dagger a_2-\tilde a_1^\dagger\tilde a_1\tilde a_2^\dagger\tilde a_2 )     -\frac{\gamma_1t}2(a_1^\dagger a_1+\tilde a_1^\dagger\tilde
a_1 )     -\frac{\gamma_2t}2(a_2^\dagger a_2+\tilde a_2^\dagger\tilde a_2)\right\}\\   \times \exp\left[\frac{\gamma_1\left({{\rm e}}^{{{\rm i}}\chi t(a_2^\dagger a_2-\tilde    a_2^\dagger\tilde a_2)-\gamma_1t}-1\right)}{{{\rm i}}\chi ( a_2^\dagger a_2-\tilde    a_2^\dagger\tilde a_2)-\gamma_1}\,|\alpha_1|^2\right ] \exp\left[\frac{\gamma_2\left({{\rm e}}^{{{\rm i}}\chi t(a_1^\dagger a_1-\tilde    a_1^\dagger\tilde a_1)-\gamma_2t}-1\right)}{{{\rm i}}\chi ( a_1^\dagger a_1-\tilde    a_1^\dagger\tilde a_1)-\gamma_2}|\alpha_2|^2\right ]    |\alpha_1\rangle|\tilde\alpha_1^*\rangle|\alpha_2\rangle|\tilde\alpha_2^*\rangle .    \quad \label{evolution}\end{aligned}\ ] ] a distinctive feature of the solution ( [ evolution ] ) is that the exponential terms contain only operators diagonal in the number - state basis .
thus , in this basis the solution ( [ evolution ] ) turns into @xmath104t\right\ } } \\ & & \times\exp{\left\{{{\rm i}}\chi t ( km - ln ) -\frac{\gamma_1t}{2}(k+l)-\frac{\gamma_2t}2(m+n)+f^{(2)}_{mn}(t)+f^{(1)}_{kl}(t)\right\ } }    \label{evolution_fock}\end{aligned}\ ] ] where @xmath105      now let us analyze the solution ( [ evolution_fock ] ) in more detail and consider for the moment the case of no dephasing ( @xmath106 ) . in this case
the expression for the purity of the state given by the solution ( [ evolution_fock ] ) is quite similar in structure to this solution itself : @xmath107    to describe effects of simultaneous damping ( photon loss ) and cross - kerr nonlinearity on the generated state , consider first a ` short - time regime ' where the following relations hold for all the photon numbers @xmath108 ( @xmath109 ) that have a non - negligible probability of occurring in the state @xmath110 ( @xmath111 ) : @xmath112 then expressions ( [ f - functions ] ) can be expanded as @xmath113^j\over(j+1 ) ! } , \quad f^{(1)}_{kl}(t)= \gamma_2t|\alpha_2|^2\sum\limits_{j=0}^{\infty}{[{{\rm i}}\chi t(k - l)-\gamma_2t]^j\over(j+1)!}. \label{f - functions_expand}\end{aligned}\ ] ] the first three terms of these expansions already describe quite well typical effects produced by simultaneous damping and cross - kerr nonlinearity . for example , for @xmath114 in the limit of small times one has @xmath115-{{{\rm
i}}\over6}\gamma_1\chi t^2(m - n)\right).\ ] ] obviously , first three terms in the round brackets of this expressions do not lead to breaking of the purity of the state . retaining only them in expansion
renders unity value for the right - hand side of eq.([evolution_purity ] ) . indeed , assuming @xmath116 one obtains the time - dependent density matrix formally coinciding with the result for no photon loss @xcite : @xmath117^k[\alpha_1^*(t)]^l}{\sqrt{k!l ! } }    \overline{|\alpha_2^*(t){{\rm e}}^{{{\rm i}}\chi lt}}\rangle . \label{rho_no_dephasing}\end{aligned}\ ] ] here time - dependent amplitudes do not depend on the numbers @xmath108 : @xmath118 the state given by eqs .
( [ evolution_fock ] ) remains negligibly affected by losses if @xmath119 .
thus , the considered scheme of non - classical state generation is quite robust with respect to photon loss ( in drastic difference with propagation losses of already generated cat - state where the off - diagonal terms @xmath120 will decay with the rates proportional to @xmath121 ) . in addition
, it is interesting to note , that one might be able to satisfy conditions @xmath119 in schemes involving dispersive atom - field interactions in qed where it is possible to restrict losses to the photon loss of cavity modes @xcite .
moreover , further in this work we consider ways to circumvent an influence of losses by making them correlated .    remarkably , a purity of the generated bimodal state can be preserved not only in the case of small losses , but also for large loss .
indeed , in the limits of large losses one can consider @xmath122 and @xmath123 as small quantities and expand functions ( [ f - functions ] ) in the following manner : @xmath124 + \chi^2(k - l)^2\left[\left({t^2\over2}+{t\over\gamma_2}+{1\over\gamma_2 ^ 2}\right){{\rm e}}^{-\gamma_2t}- { 1\over\gamma_2 ^ 2}\right ] , \\
\nonumber |\alpha_1|^{-2}f^{(2)}_{mn}(t)\approx ( 1-{{\rm e}}^{-\gamma_1t})-{{\rm i}}\chi(m - n)\left[\left(t+{1\over\gamma_1}\right){{\rm e}}^{-\gamma_1t}-{1\over\gamma_1}\right]+ \chi^2(m - n)^2\left[\left({t^2\over2}+{t\over\gamma_1}+{1\over\gamma_1 ^ 2}\right){{\rm e}}^{-\gamma_1t}- { 1\over\gamma_1 ^ 2}\right ] . \label{tailor}\end{aligned}\ ] ] note that this approximation holds for arbitrary interaction times . if the interaction time is sufficiently large to fulfill the conditions @xmath125 then the state given by eqs .
( [ evolution_fock ] ) is practically pure . under the conditions ( [ long_time_conditions ] )
this state is of the form described by eq .
( [ rho_no_dephasing ] ) with the time - dependent amplitudes given by @xmath126\right\ } , \\ \alpha_2(t)=\alpha_2\exp\left\{-\frac{\gamma_2t}{2}-{{\rm i}}\chi|\alpha_1|^2
\left[t{{\rm e}}^{-\gamma_1t}-\gamma_1^{-1}(1-{{\rm e}}^{-\gamma_1t})\right]\right\}. \label{functions}\end{aligned}\ ] ] clearly , purity of the resulting state in the long time - limit is precisely a consequence of a strong photon loss . in this way
strong photon loss paradoxically suppresses state mixing predicted by the general solution ( [ evolution_fock ] ) .
there is another interesting feature that distinguishes the losses occurring in the process of kerr interaction from the losses that take place after the interaction . in particular , the non - gaussian state generated by the cross - kerr interaction in the scheme discussed in @xcite
can retain its non - classical features even for the loss level , which would completely eliminate any such features in case of free propagation of the state . to be specific ,
a typical signature of non - classicality of a quantum state is the fact that its wigner function is negative in some regions of the phase space @xcite .
we will show that for a 50% photon loss occurring in the scheme @xcite during the cross - kerr interaction , the wigner function of the output state retains its negativity while the same photon loss occurring after the interaction would make the wigner function necessarily positive .
of the non - classical non - gaussian state of ref .
@xcite exhibits strong negativity .
wigner function is viewed along the imaginary axis of the @xmath127 plane .
, width=377 ]     for about @xmath128 photon loss in mode @xmath63 corresponding to @xmath129 .
the negative region of the wigner function is now negligible , but still present , and the amplitude of the state is clearly damped .
, width=377 ]    in the scheme @xcite for generating non - classical states , two coherent states in modes @xmath130 ( @xmath131 in notations of @xcite ) interact via cross - kerr effect in a non - linear medium and subsequently the @xmath132-quadrature of mode @xmath7 is measured .
the resulting state of mode @xmath6 exhibits wigner function with negative regions ( see fig .
[ wigner0 ] ) and a characteristic crescent ( or banana ) shape .
suppose mode @xmath6 is subject to losses during the cross - kerr interaction , i.e. , @xmath133 in eq .
( [ evolution ] ) while we assume @xmath134 .
the photon loss in mode @xmath6 is given by the reduction of the coherent amplitude described by eq .
( [ functions1 ] ) .
consider the situation when the mean photon loss is @xmath128 .
this corresponds to @xmath135 .
the plot of the corresponding wigner function is shown in fig .
[ wigner07 ] which shows clearly that although the negative region of the wigner function is surpressed , it is still present .    now compare this with the situation when losses are introduced to mode @xmath6 after the lossless cross - kerr interaction has taken place .
such losses are equivalent to mixing mode @xmath6 with the vacuum state on a beam splitter ( bs ) and discarding one bs output .
a 50% loss corresponds to a 50/50 bs .
it is known  @xcite that for such a balanced beam splitter the wigner function of one bs output can be expressed as a scaled husimi q - function of the input state : @xmath136 the q - function of a state @xmath137 is defined as @xmath138 and is clearly non - negative for all coherent state amplitudes @xmath80 , i.e. , in the whole phase space of the mode
. therefore negative regions of the wigner function can not survive losses larger than 50% if these occur during propagation of the generated state .
hence losses that take place in the process of state generation via the cross - kerr interaction are less harmful to the non - classicality of the output state than the same level of loss after the interaction .
as can be seen from the solution ( [ evolution_fock ] ) , an influence of dephasing into independent reservoirs can be profoundly destructive .
dephasing leads to diminishing of the off - diagonal elements in the number - state basis with rates proportional to the difference of these numbers .
so for the coherent state with the amplitude @xmath139 , a condition @xmath140 should be fulfilled for the interaction time @xmath141 to consider the influence of dephasing negligible . in the recently discussed qed schemes for generating non - linearity ( including eit - like ones ) ,
dephasing is usually disregarded without being estimated ( see , for example , @xcite ) .
however , since the rate of losses increases with increasing intensity of the coherent states used in the generation process , more caution is required with respect to the dephasing . for schemes involving large cross - kerr nonlinearity in the solid - state structures ( such as , for example , photonic crystals )
emitter - mediated dephasing could be a major source of state decoherence and a reason for failure of schemes involving initial coherent states with large number of photons .
however , remarkably , if the cross - kerr nonlinear interaction scheme is designed in such a way that losses due to dephasing are correlated , then it may be possible to avoid their destructive effects .
this is the subject of the next section .
it is well established , that the states of quantum systems can be correlated and even entangled through interaction with the common reservoir .
this phenomenon can occur even in absence of any direct interaction between systems . in in appendix b3
we give an example of such a phenomenon for a scheme of generating the cross - kerr nonlinearity via dispersive interaction of modes with emitters .
there a beam - splitting action of the common reservoir is considered , and it is shown how such a reservoir can produce a stationary entangled state of two modes .    here we focus our attention on another important possibility : namely , on a way to neutralize a destructive influence of losses in the process of generation by rendering these losses correlated and exploit a correlating effect of the reservoir .
let us consider an example of the realistic scheme to produce the large cross - kerr nonlinearity described in section  [ lambda - system ] ( see eq .  ( [ hamiltonian_lambda_system ] ) and the text thereafter ) .
we consider the case of @xmath142 so the master equation ( [ master_lambda_model ] ) now transforms as @xmath143+{1\over2}\left[(\gamma_1-\gamma_{12 } ) \mathcal{l}(a_1)+(\gamma_2-\gamma_{12})\mathcal{l}(a_2)\right]\rho(t)+ { \gamma_{12}\over2}\mathcal{l}(a_1+a_2)\rho(t ) .
\label{master_lambda_model1}\end{aligned}\ ] ] in absence of decoherence ( i.e. , @xmath144 ) , the scheme described by eq .
( [ master_lambda_model ] ) is able to generate entangled superpositions of schrdinger - cat states from initial coherent states of modes @xmath6 and @xmath7 ( for similar schemes see , for example , ) .
it is easy to see that for @xmath145 the scheme produces an entangled superposition of coherent states from a pair of initially uncorrelated coherent states , which reads as @xmath146 in the presence of correlated reservoirs one can find a solution of eq .
( [ master_lambda_model1 ] ) noticing that under the rotation ( [ rotation ] ) the hamiltonian part of eq .
( [ master_lambda_model1 ] ) remains invariant .
thus , performing the rotation one obtains @xmath147+{1\over2}\left[{\bar\gamma}_1 \mathcal{l}(b_1)+{\bar\gamma}_2\mathcal{l}(b_2)\right]\rho(t ) \label{master_lambda_model2}\end{aligned}\ ] ] with the exact solution given by eq .
( [ ham ] ) . also , one immediately sees that for the completely correlated reservoirs ( i.e. @xmath49 ) , the mode @xmath148 is not affected by the loss , as @xmath149 then .
naturally , a cat state can be generated in this mode .
thus , in the limit of large loss , @xmath150 , and for @xmath151 it follows from eq .
( [ master_lambda_model2 ] ) that @xmath152 where @xmath153 , @xmath154 .
we also assume that @xmath155 . from eq .
( [ lambda_dissipation_exact ] ) , we derive a conclusion that the only effect of completely correlated loss is a reduction of amplitudes of the coherent states forming the superposition ( [ lambda_nodissipation_exact ] ) .
thus , we have seen that by making losses completely correlated one can completely avoid decoherence caused by these losses .
of course , in practice one can hardly have completely correlated reservoirs due to presence of additional uncorrelated loss ( such as modal losses due to coupling to additional reservoirs etc . ) .
nevertheless , designing the scheme such as to have predominantly correlated losses might greatly enhance its robustness in production of non - classical states .
we illustrate this with the simple example of the conditioned cat - state generation from the solution of eq .
( [ master_lambda_model ] ) .
if the rotated mode @xmath156 impinges on the detector , in the case of no signal on the detector the rotated mode @xmath148 is ( up to the normalization factor ) in the state : @xmath157 , @xmath158 , and @xmath159    in fig .
[ fig2 ] one can see examples of the wigner function of the state ( [ evolution_fock2 ] ) . for no loss ( fig .
[ fig2](a ) ) and @xmath145 the state ( [ evolution_fock2 ] ) is a usual schrdinger cat state with the pronounced oscillations near the origin .
large correlated loss ( @xmath160 ) changes the size of the cat and rotates it ( fig .
[ fig2](b ) ) , but otherwise leaves it intact .
lower correlated loss distorts the cat ( fig .
[ fig2](d ) ) due to influence of additional mixing between modes in the interaction process ( as it follows from eq .
( [ evolution_fock2 ] ) ) . nevertheless , the state is strongly non - classical .
uncorrelated loss with the same rate eliminates the non - classicality outright ( fig .
[ fig2](c ) ) .
correlated loss allows for the non - classicality to survive ( fig .
[ fig2](f ) ) even if the uncorrelated loss with the rate equal to difference between individual rates and the correlation rate ( e.g. , @xmath161 ) destroys the non - classicality completely ( fig .
[ fig2](e ) ) .
in what presented here , we followed the quest to find ways to impair the decoherence processes in quantum state generation and manipulation . for the particular class of nonlinear interaction processes ,
we have found two striking examples of loss dynamics , for which the losses themselves counteract decoherence : loss - mediated correlations between the interacting modes and losses to the correlated reservoirs .
the latter result in strongly correlated modal loss .
these correlated losses influence the dynamics of the modes undergoing the cross - kerr nonlinear interaction in a completely different way than losses into independent reservoirs , the aspect of quantum nonlinear dynamics to large extent unexplored so far .
thus , remarkably , designing the schemes for the generation of the kerr nonlinearity in such a way that losses in this process are correlated , one can greatly diminish their destructive impact and even exploit them for entanglement generation ( see also appendix b ) .
as to the origin of this effect , if both modes , for example , interact with the same emitter transition , emitter losses are likely to lead to the correlated loss of both of these modes .
note , that in this case a significant modal loss might occur even in the case when the emitter subject to losses stays in superposition of metastable levels with negligibly small probability to occupy higher , decaying levels . for more details and examples on the origin of the correlated losses and
their entangling effect see appendix b.    turning to the other aforementioned unexpected aspect of the quantum dynamics , we have demonstrated that losses in the nonlinear process of state generation affect the quantum state in quite a different way to the propagation losses of already generated nonclassical state .
this is mainly due to the loss - mediated correlations between the modes participating in the cross - kerr interaction .
in addition , losses through coupling to the correlated reservoirs ( correlated loss ) further enhance the difference between the decoherence processes _ during _ and _ after _ the state generation
. in particular , non - classicality seems to be more robust with respect to the generation loss than to the propagation loss .
we discussed an example of generation loss exceeding @xmath128 with negative values of the wigner function preserved , whereas the propagation loss exceeding @xmath128 renders the wigner function completely positive .
the authors gratefully acknowledge funding from the eu project fp7 - 212008 compas ( n. k. and t.t . ) , equind project of 6fp ist-034368 , ffr of belarus project f08p-131 and fapesp , brazil , project 2008/57657 - 3 ( d.m . ) .
also , d.m . would like to thank the school of physics and astronomy , university of st .
andrews , for their hospitality .
we start from the general effective hamiltonian @xmath0 of eq .
( [ ham_ger ] , [ vloss ] ) describing both self- and cross - interaction and interaction of modes with reservoirs responsible for losses ( see section  [ master_section ] ) : @xmath162 setting @xmath17 for sake of simplicity , and using a time - convolutionless projection operator technique , one obtains the following master equation @xcite : @xmath163 - \\ & & \int\limits_0^{t}{\,{\rm d}}\tau \left\{\langle v_{\rm loss}(t)v_{\rm loss}(\tau)\rangle_{\rm r}\rho(t)+\rho(t)\langle v_{\rm loss}(\tau)v_{\rm loss}(t)\rangle_{\rm r } - \langle v_{\rm loss}(\tau)\rho(t)v_{\rm loss}(t)\rangle_{\rm r}- \langle v_{\rm loss}(t)\rho(t)v_{\rm loss}(\tau)\rangle_{\rm r}\right\}\qquad \label{master_gener}\end{aligned}\ ] ] where @xmath18 denotes the density matrix of the system averaged over states of all reservoirs ( and over possible stochastic variables , too ) , and @xmath164 denotes averaging over all reservoirs .
reservoir modes are assumed to be initially in the vacuum states .
further , we consider reservoirs of different types to be independent ; correlation functions @xmath165 , @xmath166 are taken to be zero .
we consider the reservoirs of the same kind to be mutually correlated , i.e. , for markovian dephasing reservoirs we assume @xmath167 here @xmath168 is the delta - function and the rates @xmath169 , are real and non - negative . the cross - correlation parameter @xmath170 is taken to be real .
for the photon loss reservoir , we assume that the term @xmath171 in the hamiltonian ( [ ham_ger ] ) preserves the total number of photons , i.e. only non - zero correlation functions are @xmath172 also here , for simplicity , the rates @xmath173 are assumed to be real and non - negative and the cross - correlation parameter @xmath174 to be real . for self - loss rates , @xmath175 , @xmath176 , and ` cross'-loss rates , @xmath174 , @xmath170 , the following relations hold @xmath177    under the assumptions made above , the master equation of eq .
( [ master_markov ] ) is obtained from eq .
( [ master_gener ] ) .
we will illustrate the process of appearance of a _ correlated modal photon loss _ in a process of off - resonant interaction between a mode and an emitter with a simple example .
consider two modes of the same frequency interacting off - resonantly with just a single two - level system ( tls ) , see fig .
[ fig1](a ) . in the rotating - wave approximation and in the interaction picture with respect to the free hamiltonian of the reservoir , in the frame rotating with the tls transition frequency @xmath178 ,
one has the following hamiltonian describing the problem : @xmath179 here @xmath180 and @xmath181 is the mode frequency ; @xmath40 are interaction constants for the corresponding modes ; @xmath182 and @xmath183 are pauli operators for the tls , @xmath184 , @xmath185 , @xmath186 ; vectors @xmath56 , @xmath187 describe the lower and the upper tls levels , correspondingly .
the reservoir operator @xmath188 describes the tls energy loss .
we assume this reservoir to be markovian and the following relations hold : @xmath189 where brackets @xmath190 denote an averaging over the reservoir . here
we are assuming that losses are weak ( i.e. , @xmath191 ) .
we adopt the usual conditions for an adiabatic elimination of the emitter , i.e. the tls starts at the lower level , and the detuning @xmath192 between the mode frequency and the tls transition frequency is much larger than @xmath40 .
thus , the tls upper level remains practically unpopulated .
changing to the interaction picture with respect to the part of the hamiltonian ( [ jkmodel_ham_spon ] ) corresponding to the absence of the tls - field interaction , we get the following interaction hamiltonian @xmath193 with the bosonic annihilation operator for the collective mode @xmath194 and @xmath195 . a formal solution for @xmath196 without losses up to the third - order terms
can be approximated as @xmath197x^{\dagger}(t_1 ) , \label{jkmodel_solution_sigma_spon}\end{aligned}\ ] ] where @xmath198 and @xmath199 , if one takes into account the fact that the modal dynamics is very slow on the scale of the tls dynamics and @xmath200 can be considered as practically constant . from eq .
( [ jkmodel_solution_sigma_spon ] ) , neglecting small and rapidly oscillating terms , after averaging over the atomic variables one arrives to the following effective interaction hamiltonian @xmath201 . \label{jkmodel_ham_v_app_spon}\end{aligned}\ ] ] deriving the master equation in the standard manner , one obtains an equation describing the correlated photon losses @xmath202+{\bar\gamma } \mathcal{l}(c)\rho(t ) , \label{master_jkmodel_spon}\end{aligned}\ ] ] where @xmath203 , @xmath204 and @xmath205 .    effectively , eq .  ( [ master_jkmodel_spon ] )
describe both coupling between modes and their interaction with the same reservoir .
both these interaction might lead to the entanglement between modes . as seen in section v , even in the absence of direct intermodal coupling ( i.e. for @xmath206 , @xmath207 ) an interaction of modes with the reservoir entangles these modes .
the analysis made above can be readily generalized to different schemes of cross - kerr nonlinearity generation through interaction of two modes with the same ensemble of emitters . in subsection  [ lambda - system ]
we devise the procedure for the scheme of the giant cross - kerr nonlinearity generation suggested in ref .
@xcite .
[ [ an - example - dispersive - two - mode - jaynes - cummings - model - with - dephasing ] ] an example : dispersive two - mode jaynes - cummings model with dephasing ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~    here we illustrate an appearance of _ correlated modal dephasing _ with the example of the jaynes - cummings model ( fig .  [ fig1 ] ) considered in the previous subsection .
we model an influence of the dephasing reservoir as a stochastic fluctuation of the tls transition frequency . in the rotating - wave approximation the problem is described by the following hamiltonian : @xmath208\sigma_z+ \left[\sigma^+(g_1a_1+g_2a_2)+(g_1a_1^{\dagger}+g_2a_2^{\dagger})\sigma^-\right ] , \label{jkmodel_ham}\end{aligned}\ ] ] where @xmath209 is random process describing a small rapid stochastic modulation due to non - radiative interaction with surroundings . for simplicity
we take @xmath209 to be just a white noise satisfying the following relations @xmath210 where @xmath211 denotes classical averaging .
we consider the case of the weak loss , @xmath212 .    as before
, we assume that the conditions for adiabatic elimination of the emitter hold . in the interaction picture with respect to the part of the hamiltonian ( [ jkmodel_ham ] ) corresponding to the absence of the tls - field interaction
, we have the following interaction hamiltonian @xmath213 a formal solution for @xmath196 in this interaction picture can be approximated as @xmath214c^{\dagger}(t_1)f^*(t_1 ) , \label{jkmodel_solution_sigma_def}\end{aligned}\ ] ] where @xmath215      averaging over dephasing noise and neglecting small terms , we get the following master equation using the standard technique implemented to derive eq .
( [ master_gener ] ) : @xmath218+{\bar d } \mathcal{l}(c^{\dagger}c)\rho(t ) , \label{master_jkmodel}\end{aligned}\ ] ] where @xmath219 the calculations of the coefficients in ( [ master_jkmodel_coefficients ] ) are carried out using the following property : @xmath220 and the fact that the detuning @xmath192 is assumed to be large , @xmath221 .
so , one can see that dephasing of the atom leads to appearance of the correlated modal dephasing practically in the same manner as atomic population losses lead to the correlated modal loss considered in the previous subsection .
also , modal dephasing occurs notwithstanding the fact that the upper atomic level remains practically unpopulated .
finally , we demonstrate that the scheme described by the master equation ( [ master_jkmodel_spon ] ) can effectively produce entanglement between the modes .
in fact , this scheme can act as a kind of lossy beam - splitter even in the absence of intermodal interaction in unitary part of eq .
( [ master_jkmodel_spon ] ) .
indeed , let us consider a completely uncorrelated single - photon initial state @xmath222 , and the initial density matrix @xmath223 . in the zero- and single - photon subspaces one can assume the following orthonormal basis : @xmath224 is not affected by the losses described by eq .
( [ master_jkmodel_spon ] ) because @xmath225 . also , the following relations are satisfied : @xmath226 thus the system of equations for the density matrix elements can be obtained from eq .
( [ master_jkmodel_spon ] ) : @xmath227\rho_{+-}(t ) , \\
\rho_{-v}(t)=\rho_{-v}(0 ) , \quad { { \,{\rm d}}\over { \,{\rm d}}t}\rho_{+v}(t)=-[{\bar\gamma}+{{\rm i}}(\delta w+\chi)]\rho_{+v}(t ) . \label{spon_system}\end{aligned}\ ] ]    the solution ( [ spon_system ] ) describes an emergence of entanglement form the initially uncorrelated state of both modes ( single photon in the mode @xmath228 and vacuum of the mode @xmath229 ) .
figure  [ fig5 ] depicts a measure of entanglement , a negativity as given in ref .
@xcite @xmath230 where @xmath231 is the density matrix @xmath137 partially transposed with respect to the first mode .
non - zero value of the negativity means that the state is the entangled one .
it can be seen that the decay rate into the common reservoir does not affect the finally reached entanglement .
this rate affect only time during which a stationary state is reached .
it follows from this system of equations that the initial uncorrelated state @xmath232 under the action of the correlated modal loss asymptotically turns into @xmath233 the state ( [ spon_entangled_state ] ) is entangled for an arbitrary @xmath234 .
however , the maximal degree of asymptotic entanglement is reached when @xmath235 , and with increasing of difference between @xmath236 and @xmath237 the asymptotic entanglement decreases ( figure  [ fig5 ] ) .
note that the state ( [ spon_entangled_state ] ) is influenced neither by the cross - kerr interaction of the modes nor by the linear excitation exchange , and one can set both @xmath238 .
the same type of state is produced by a correlated dephasing described by eq .
( [ master_jkmodel ] ) .
effectively , the correlation of reservoirs allows for existence of decoherence - free subspaces to which the two - mode state eventually evolves @xcite .
entangling through the common reservoir with appearance of the long - living state similar to the one described by eq .
( [ spon_entangled_state ] ) might occur for emitters and collective reservoir modes near the band - edge in photonic crystals @xcite .
d. k. armani , t. j. kippenberg , s. m. spillane and k. j. vahala , nature * 421 * , 925 ( 2003 ) ; t. j. kippenberg , s.m .
spillane , and k. j. vahala , phys .
rev . lett . * 93 * 083904 ( 2004 ) ; s. m. spillane , t. j. kippenberg , and k. j. vahala , k. w. goh , e. wilcut , and h. j. kimble , phys .
rev . a*71 * 013817 ( 2005 ) .