Quantum entanglement in coupled lossy waveguides

We investigate the viability of coupled waveguides as basic units of quantum circuits. We study entanglement when the waveguides are fed in by light produced by a down-converter working either in low gain limit or under large gain. We present explicit analytical results for the measure of entanglement in terms of the logarithmic negativity for a variety of input states. We also address the effect of loss on entanglement dynamics of waveguide modes. Our results indicate that the waveguide structures are reasonably robust against the effect of loss and thus quite appropriate for quantum architectures as well as for the study of coherent phenomena like random walks. Our analysis is based on realistic structures used currently. © 2010 Optical Society of America OCIS codes: (270.5585) Quantum information and processing; (250.5300) Photonic integrated circuits; (230.7370) Waveguides. References and links 1. R. Morandotti, U. Peschel, J. S. Aitchison, H. S. Eisenberg, and Y. Silberberg, “Experimental Observation of Linear and Nonlinear Optical Bloch Oscillations,” Phys. Rev. Lett. 83, 4756–4759 (1999). 2. D. N. Christodoulides, F. Lederer, and Y. Silberberg, “Discretizing light behaviour in linear and nonlinear waveguide lattices,” Nature (London) 424, 817–823 (2003). 3. S. Longhi, “Optical analog of population trapping in the continuum: Classical and quantum interference effects,” Phys. Rev. A 79, 023811 (2009). 4. H. B. Perets, Y. Lahini, F. Pozzi, M. Sorel, R. Morandotti, and Y. Silberberg, “Realization of Quantum Walks with Negligible Decoherence in Waveguide Lattices,” Phys. Rev. Lett. 100, 170506 (2008). 5. Y. Bromberg, Y. Lahini, R. Morandotti, and Y. Silberberg, “Quantum and Classical Correlations in Waveguide Lattices,” Phys. Rev. Lett. 102, 253904 (2009). 6. U. Peschel, T. Pertsch, and F. Lederer, “Optical Bloch oscillations in waveguide arrays,” Opt. Lett. 23, 1701–1703 (1998). 7. T. Pertsch, P. Dannberg, W. Elflein, A. Brauer, and F. Lederer, “Optical Bloch Oscillations in Temperature Tuned Waveguide Arrays,” Phys. Rev. Lett. 83, 4752–4755 (1999). 8. R. Iwanow, D. A. May-Arrioja, D. N. Christodoulides, G. I. Stegeman, Y. Min, and W. Sohler, “Discrete Talbot Effect in Waveguide Arrays,” Phys. Rev. Lett. 95, 053902 (2005). 9. S. Longhi, “Optical Bloch Oscillations and Zener Tunneling with Nonclassical Light,” Phys. Rev. Lett. 101, 193902 (2008). 10. A. Rai, G. S. Agarwal, and J. H. H. Perk, “Transport and quantum walk of nonclassical light in coupled waveguides,” Phys. Rev. A 78, 042304 (2008). (C) 2010 OSA 15 March 2010 / Vol. 18, No. 6 / OPTICS EXPRESS 6241 #121578 $15.00 USD Received 21 Dec 2009; accepted 18 Feb 2010; published 12 Mar 2010 11. A. Rai and G. S. Agarwal, “Possibility of coherent phenomena such as Bloch oscillations with single photons via W states,” Phys. Rev. A 79, 053849 (2009). 12. A. Politi, M. J. Cryan, J. G. Rarity, S. Yu, and J. L. O’Brien, “Silica-on-Silicon Waveguide Quantum Circuits,” Science 320, 646–649 (2008). 13. J. C. F. Matthews, A. Politi, A. Stefanov, and J. L. O’Brien, “Manipulation of multiphoton entanglement in waveguide quantum circuits,” Nature Photonics 3, 346-350 (2009). 14. A. N. Boto, P. Kok, D. S. Abrams, S. L. Braunstein, C. P. Williams, and J. P. Dowling, “Quantum Interferometric Optical Lithography: Exploiting Entanglement to Beat the Diffraction Limit,” Phys. Rev. Lett. 85, 2733–2736 (2000). 15. C. H. Bennett and D. P. DiVincenzo, “Quantum information and computation,” Nature 404, 247–255 (2000). 16. W. H. Zurek, “Decoherence, einselection, and the quantum origins of the classical ,” Rev. Mod. Phys. 75, 715-775 (2003). 17. B. Do, M. L. Stohler, S. Balasubramanian, D. S. Elliott, C. Eash, E. Fischbach, M. A. Fischbach, A. Mills, and B. Zwickl, “Experimental realization of a quantum quincunx by use of linear optical elements,” J. Opt. Soc. Am. B 22, 499–504 (2005). 18. P. K. Pathak and G. S. Agarwal, “Quantum random walk of two photons in separable and entangled states,” Phys. Rev. A 75, 032351 (2007). 19. A. Politi, J. C. F. Matthews, and J. L. O’Brien, “Shors Quantum Factoring Algorithm on a Photonic Chip,” Science 325, 1221 (2009). 20. B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics, 2nd Edition, (Wiley, New York 2007), p. 319. 21. W. K. Lai, V. Bužek, and P. L. Knight, “Nonclassical fields in a linear directional coupler ,” Phys. Rev. A 43, 6323-6336 (1991). 22. R. Iwanow, R. Schiek, G. I. Stegeman, T. Pertsch, F. Lederer, Y. Min, and W. Sohler, “Observation of Discrete Quadratic Solitons ,” Phys. Rev. Lett. 93, 113902 (2004). 23. U. Peschel, R. Morandotti, J. M. Arnold, J. S. Aitchison, H. S. Eisenberg, Y. Silberberg, T. Pertsch, and F. Lederer, “Optical Discrete Solitons in Waveguide Arrays. 2. Dynamic Properties,” J. Opt. Soc. Am. B 19, 2637– 2644 (2002). 24. K. B. Mogensen, F. Eriksson, O. Gustafsson, R. P. H. Nikolajsen, and J. P. Kutter, “Pure-silica optical waveguides, fiber couplers, and high-aspect ratio submicrometer channels for electrokinetic separation devices,” Electrophoresis 25, 3788–3795 (2004). 25. G. Vidal and R. F. Werner, “Computable measure of entanglement,” Phys. Rev. A 65, 032314 (2002). 26. C. K. Hong, Z. Y. Ou, and L. Mandel, “Measurement of subpicosecond time intervals between two photons by interference,” Phys. Rev. Lett. 59, 2044-2046 (1987). 27. G. S. Agarwal, “Entropy, the Wigner Distribution Function, and the Approach to Equilibrium of a System of Coupled Harmonic Oscillators ,” Phys. Rev. A 3, 828-831 (1971). 28. P. J. Dodd and J. J. Halliwell, “Disentanglement and decoherence by open system dynamics ,” Phys. Rev. A 69, 052105 (2004). 29. L. -M . Duan, G. Giedke, J. I. Cirac, and P. Zoller, “Inseparability Criterion for Continuous Variable Systems,” Phys. Rev. Lett. 84, 2722-2725 (2000). 30. R. Simon, “Peres-Horodecki Separability Criterion for Continuous Variable Systems,” Phys. Rev. Lett. 84, 27262729 (2000). 31. G. Adesso, A. Serafini, and F. Illuminati, “Extremal entanglement and mixedness in continuous variable systems,” Phys. Rev. A 70, 022318 (2004). 32. S.M. Barnett and P. Radmore, Methods in Theoretical Quantum Optics, (Oxford University Press, 2002), p. 168. 33. Sumanta Das and G. S. Agarwal, “Bright and dark periods in the entanglement dynamics of interacting qubits in contact with the environment,” J. Phys. B 42, 141003 (2009). 34. Sumanta Das and G. S. Agarwal, “Decoherence effects in interacting qubits under the influence of various environments,” J. Phys. B 42, 205502 (2009). 35. S. Longhi, “Transfer of light waves in optical waveguides via a continuum,” Phys. Rev. A 78, 013815 (2008).


Introduction
Discrete optical systems like coupled waveguides are known to be extremely efficient in manipulating the flow of light and have been investigated extensively in the last two decades [1,2,3,4,5,6,7,8,9,10].Many key quantum effects like quantum interference, entanglement and quantum walk has been investigated in these systems [3,4,5,10].For example, using coherent beam Peretes et al. [4] have observed quantum walk effects in a system consisting of large number of waveguides.In another experiment, Bromberg et al. investigated the quantum correlations in GaAs waveguide arrays [5] using two-photon input states.In particular, they considered both the separable and entangled two photon state and observed various features as-sociated with quantum interference.In addition, the coupled waveguide arrays have been used to study the discrete analogue of the Talbot effect [8].The entanglement between the waveguide modes and behavior of nonclassical light in coupled waveguides has been studied [9,11].In a recent experiment Politi et.al. [12] have shown how a CNOT gate can be implemented on a single Silicon chip using coupled silica waveguides, thus showing possible application of waveguides in quantum computation.They also observed Hong-Ou -Mandel two photon interference in these coupled waveguides.In a following experiment [13] coupled silica waveguides was used to generate a multimode interferometer on an integrated chip.It was further shown that these interferometers can be used to generate arbitrary quantum circuits.They also showed that two and four photon entangled states similar to NOON states [14] can be generated on the silicon chip.All these studies have hence given a new impetus to the field of quantum information processing and quantum optics with waveguides.In particular, for effective use of these waveguide circuits in quantum computation and communication tasks sustainability of generated entanglement is very important [15].In light of this, it is imperative to study entanglement in waveguides using quantitative measures for entanglement.This is the main purpose of the present study.Moreover, in practice the waveguides are not completely lossless.Thus an immediate question of interest would be how does this loss affects the entanglement in the waveguide modes ?It is well known that entanglement is quite susceptible to decoherence [16] and thus the above question bears immense interest in context to quantum information processing using waveguides.Further it is important to understand the role of loss in coherent phenomena like quantum random walk [4,17,18].
In this paper we investigate these in a simple system of two single mode waveguides, which are coupled through the overlap of evanescent fields.This simple system serves as a unit or the basic element for constructing a quantum circuit [19].The input light to the coupled waveguide system is usually produced by the parametric down-conversion process.At high gain the parametric down-conversion process produces a squeezed state of light while single photon states are produced at low gain.Behavior of photon number states such as the single photon state and the NOON state have also been investigated in these systems [12,13,19].We thus consider a variety of nonclassical input states like squeezed states and photon number states which have been extensively investigated in couple waveguide system and study their respective entanglement dynamics.We quantify the evolution of entanglement in terms of logarithmic negativity and present explicit analytical results for both squeezed and number state inputs.We further investigate the question of possible effects of loss on the entanglement dynamics in waveguides by considering lossy waveguide modes.We find that in this case, for both number state inputs as well as squeezed state inputs, entanglement shows considerable robustness against loss.
The organization of the paper is as follows: In section 2, we describe the model and derive analytical result for the field modes of the coupled waveguide system.In section 3, we study the evolution of entanglement for two classes of photon number states, (A) separable single photon state |1, 1 and (B) entangled two-photon NOON state.We quantify the degree of entanglement of these states by using the logarithmic negativity.In section 4 , we then study the time evolution of entanglement by evaluating the logarithmic negativity for two classes of squeezed input states (A) separable two mode squeezed state and (B) entangled two mode squeezed state.The effect of loss in waveguides on the entanglement dynamics is then discussed in section 5. Finally we summarize our results in section VI with a future outlook.

The model
We consider a system with two single mode waveguides, coupled through nearest-neighbor interaction as shown in Fig. 1.Let a and b be the field operators for the modes in each waveguide.These obey bosonic commutation relations [a, a † ] = 1; (a → b).The Hamiltonian describing the evanescent coupling between the waveguide mode in such a system of two coupled waveguides can be derived using the coupled mode theory [20].The coupling among the waveguides is incorporated in this framework by treating it as a perturbation to the mode amplitudes.It is assumed that the presence of the second waveguide perturbs the medium outside the first waveguide.This creates a source of polarization outside the first waveguide, which thereby leads to modification of the amplitude of the mode in it.Further, the amplitude of the modes in each waveguide is assumed to be a slowly varying function of the propagation distance.Moreover, in this perturbative approach the coupling does not effect the propagation constant or transverse spatial distribution of the waveguide modes.The field of the first waveguide has a similar effect on the second waveguide.Under these assumptions, the field mode of the composite structure are governed by the Helmholtz equation which gives two coupled first order differential equations which can be solved to obtain the time evolution of field modes in the coupled waveguide structure.The corresponding description for the nonclassical light can be studied by quantizing the field amplitudes as has been done in the work of Lai et.al. [21].Following an approach similar to that developed by Lai et.al., we can write the corresponding quantum mechanical Hamiltonian for the coupled waveguide as where the first two terms correspond to the free energy of the waveguide modes and the last two terms account for the evanescent coupling between the waveguide modes with J as the coupling strength.The coupling J depend on the distance between the waveguides.The input to the coupled waveguide system can be in a separable or an entangled state.Let γ be the loss rates of the modes a and b.The loss γ arises from the loss in the material of the waveguide.1. Approximate values of some of the parameters used in waveguide structures [22,23,24].The loss, usually quoted in dB/cm, for different waveguides is converted to frequency units used in this paper by using the formula, 10 Log( P out P in ) ≡ 10 Log(e −2γ/c ), where P in is the input power, P out is the power after traveling unit length.

Waveguide Type
Coupling parameter J (sec −1 ) Loss γ (sec As known the silica waveguides have very little intrinsic loss and should be preferable in many applications.Nevertheless the loss is to be included as this could be detrimental in long propagation for example in the study of quantum random walks.Since the two waveguides are identical, we have taken the loss rate of both the modes to be the same.We can model the loss in waveguides in the framework of system-reservoir interaction well known in quantum optics and is given by, where ρ is the density operator corresponding to the system consisting of fields in the modes a and b.The dynamical evolution of any measurable O in the coupled waveguide system is then governed by the quantum-Louiville equation of motion given by, where Ȯ = Tr{O ρ}, the commutator gives the unitary time evolution of the system under the influence of coupling and the last term account for the loss.Note that in absence of loss (lossless waveguides) the time evolution of the field operators can be evaluated using the Heisenberg equation of motion and is given by, Next we will study the entanglement characteristics of photon number and squeezed input states as they propagate through the waveguides.To keep the analysis simple in the next few sections we consider the case of lossless waveguide modes (γ = 0).We defer the discussion of loss on entanglement to Sec V.

Evolution of entanglement for input states at single photon level
In this section we study the dynamics of entanglement for photon number input state.We quantify the entanglement of the system by studying the time evolution for the logarithmic negativity [25].For a bipartite system described by the density matrix ρ the logarithmic negativity is where ρ T is the partial transpose of ρ and the symbol denotes the trace norm.Also N(ρ) is the absolute value of the sum of all the negative eigenvalues of the partial transpose of ρ.The log negativity is a non-negative quantity and a non-zero value of E N would mean that the state is entangled.

Separable photon number state as an input
We first consider the case when there is no loss and hence we set γ = 0. We assume that the input is in a separable state.Further, for studying the entanglement dynamics for photon number states we first consider the case of a single photon input in each waveguide.Thus the initial state is |ψ(0) = |1, 1 .Using Eq. ( 4) we can show that a single photon input state given by |ψ(0) evolves into a state : The coefficients α 1 , β 1 and δ 1 are given by : The density matrix corresponding to the state |ψ(t) can be written as : ρ = |ψ(t) ψ(t)| .Thus, using Eq. ( 5), we can write the log negativity E N for the state |ψ(t) as: We show the time evolution of E N for the single photon input state |1, 1 in the red curve of Fig. 2 (a).We would like to emphasize that the values of θ studied here are very similar to the ones employed in the recent experiments [5,12].At time t = 0, we begin with a separable input state and thus the value of log negativity is E N = 0.The entanglement quantified by the log negativity increases with time and attains a maximum value of 1.58 for θ 0.15.In this case the single photon state evolves into a maximally entangled state given by: Further, for θ = 1/4, we get an analog of the well known Hong-Ou-Mandel interference [26].Note that in this case the logarithmic negativity E N attains a value of 1 which is less than the corresponding value of E N for the maximally entangled state |ψ m .In addition, for θ = 1/2, we find that E N vanishes and the state at this point is e iπ |1, 1 .At later times, we see a periodic behavior which can be attributed to the inter-waveguide coupling J.We next consider the case where we have two photons in one waveguide and none in the other input.Thus the initial state can be written as :|ϕ(0) = |2, 0 .Again using Eq. ( 4) we find that the |ϕ(0) evolves into a state given by :|ϕ(t) → α 2 |2, 0 + β 2 |1, 1 + δ 2 |0, 2 .The coefficients α 2 , β 2 and δ 2 are given by : Using a similar procedure as discussed above we can evaluate the log negativity E N for the state |ϕ(t) .We show the result for the log negativity in black curve of Fig. 2 (a).In this case we find that the log negativity increases and attains a maximum value of 1.54.After reaching the maximum value the log negativity decreases and eventually becomes equal to zero.Thus the state becomes disentangled at this point of time.At later times we see a periodic behavior and the system gets entangled and disentangled periodically.Clearly the entanglement dynamics of the states |ψ(0) and |ϕ(0) are different.Unlike the earlier case for the |1, 1 input state, we don't see any interference effects in this case [5].

NOON state as an input
Next we consider the entangled state prepared in a N photon NOON state [14] as our initial state : Using Eq. ( 4) we can show that the input state given by |ψ in evolves into a state : where a(t) and b(t) are given by Eq. ( 4).Using Eq. ( 4) in the above equation, we get (11)   where C(N, k) is the Binomial coefficient given by: The density matrix corresponding to the state |ψ out can be written as : 2 is a diagonal matrix and the eigenvalues of (ρ T out ) 2 is of the form: Thus the negative eigenvalues of ρ T out are of the form |β k ||β m | (k = m) and the log negativity E N can be written as: We can use the above equation to study entanglement dynamics for the multi-photon NOON state.We first consider the entangled state prepared in a two photon NOON state as our initial state.As shown in the black curve of Fig. 2 (b), the value of E N at time t = 0 is equal to 1 which indicates entanglement.The log negativity E N in the black curve of Fig. 2 (b) shows a behavior that is similar to the result for the |1, 1 state shown in the red curve of Fig. 2 (a).Also note the shift of π/4 between the corresponding results in Figs. 2 (a) and 2 (b).As in the case of single photon input state |1, 1 , the initial state evolves into a maximally entangled state corresponding to a value of E N which is equal to 1.58.In addition, for θ = 1/2, we again see a signature of quantum interference such that the probability of getting the single photons at each of the output port vanishes [5].The logarithmic negativity E N at this point is equal to 1.At later times the entanglement shows an oscillatory behavior and the system gets periodically entangled and disentangled.The red curve in Fig. 2 (b) shows the negativity (E N − 1) for the four photon NOON state.As earlier, the value of E N at time t = 0 is equal to 1 which indicates entanglement.The curve for four photon NOON state also shows quantum interference effect.Further, the logarithmic negativity never becomes zero in this case and hence the initially entangled state remains entangled for later times.

Separable two mode squeezed state as an input
We next study the generation and evolution of entanglement for the case of squeezed input states.For this purpose we first consider a separable squeezed input state coupled to the modes a and b of the waveguide given by:|ζ = |ζ a ⊗|ζ b , where |ζ a (|ζ b ) are single mode squeezed states defined as, where r is taken to be real.Here r is proportional to the gain of the down-converter.Such a separable two mode state can be obtained as the field at the output ports of a 50-50 beam splitter whose input ports are fed by light from a down-converter.It is well known that a two mode squeezed state like |ζ can be completely characterized by its first and second statistical moments given by the first moment : ( x 1 , p 1 , x 2 , p 2 ) and the covariance matrix σ .The squeezed vacuum state falls under the class of Gaussian states.It is to be noted that evolution of Gaussian states has been studied for many different model Hamiltonians [27,28].We focus on the practical case of propagation of light produced by a down converter in coupled waveguides which currently are used in quantum architectures and quantum random walks.Note that since the first statistical moments can be arbitrarily adjusted by local unitary operations, it does not affect any property related to entanglement or mixedness and thus the behavior of the covariance matrix σ is all important for the study of entanglement.The measure of entanglement for a Gaussian state is best characterized by the logarithmic negativity E N , a quantity evaluated in terms of the symplectic eigenvalues of the covariance matrix σ [29,30].The elements of the covariance matrix σ are given in terms of conjugate observables, x and p in the form, where α, β and μ are 2 × 2 matrices given by x 1 p 2 +p 2 x 1 2 (13) Here x 1 , x 2 and p 1 , p 2 are given in terms of the normalized bosonic annihilation (creation) operators a(a † ), b(b † ) associated with the modes a and b respectively, The observables, x j , p j satisfy the cannonical commutation relation [x k , p j ] = iδ k j .The condition for entanglement of a Gaussian state like |ζ is derived from the PPT criterion [30], according to which the smallest symplectic eigenvalue ν< of the transpose of matrix σ should satisfy, ν< < 1 where ν± is given by, where Δ(σ ) = Δ( σ ) = Det(α)+Det(β )−2Det(μ).Thus according to the condition (15) when ν< ≥ 1/2 a Gaussian state become separable.The corresponding quantification of entanglement is given by the logarithmic negativity E N [25,31] defined as, which constitute an upper bound to the distillable entanglement of any Gaussian state [31].On evaluating the covariance matrix σ for the state ζ for γ = 0 (no loss), using equation ( 3), ( 4) and ( 14) we find, where c, d, e are given by c = 1 2 {cosh(2r) + sinh(2r) cos(2Jt)} ( 19) The corresponding symplectic eigenvalues ν± are then given by One can clearly see from equations ( 17), ( 19), ( 20) and ( 21) the dependence of logarithmic negativity E N on coupling strength J between the waveguides and the squeezing parameter r.In figure 3. we plot the logarithmic negativity as a function of scaled time, θ = Jt for the state |ζ .
Here t is related to the length l of the waveguide and its refractive index n by t = nl/v, v being the velocity of light.We see from figure (3) that as |ζ is separable at t = 0, E N = 0 initially but as Jt increases, it oscillates periodically between a non-zero and zero value.Thus the initially separable state |ζ becomes periodically entangled and disentangled as its propagates through the waveguide.We attribute this periodic generation of entanglement to the coupling J among the waveguides.We further find that ν< = 1/2 at certain points along the waveguide given by 2θ = (k + 1)π, k = 0, 1, 2, 3, ......Note that at this points E N vanishes and |ζ becomes separable.At all other points the state |ζ = |ζ a ⊗ |ζ b .We see that E N is maximum and has a value equal to the amount of squeezing 2r at the points given by 2θ = (k + 1)π/2.Hence at this points the initial seperable state |ζ becomes maximally entangled and is given by :|ζ ≡ exp{e iπ r(a † b † + ab)}|00 .

Entangled two mode squeezed state as an input
Let us now study the dynamical evolution of a two mode squeezed state |ξ = exp[r(a † b † − ab)]|00 as an input to the waveguide.As before we consider r to be real.To quantify the entanglement of the state |ξ we need to evaluate the logarithmic negativity E N .Thus we first evaluate the covariance matrix σ for the state |ξ using equations ( 3) with γ = 0, ( 4) and ( 14).We find σ to be where f , g and h are given by, The corresponding symplectic eigenvalues ν± is then given by, The logarithmic negativity E N can then be evaluated using equations ( 15), ( 17) and (24).From equations ( 23) and ( 24) the dependence of E N on the squeezing r and the coupling J between the waveguides is clearly visible .From equation ( 24) we find that E N = 0 i.e. entanglement become zero when, 2θ = (k + 1)π/2 as then ν< = 1/2 and thus the initially entangled state |ξ becomes separable, i.e |ξ = exp{ r 2 e iπ (a †2 + a 2 )}|0 ⊗ exp{ r 2 e iπ (b †2 + b 2 )}|0 .The behavior of time evolution of logarithmic negativity for the initial two mode squeezed state is found to be similar to figure (3).In this case though the points of zero entanglement is shifted by π/4 with respect to that for the initial separable state |ζ .The oscillatory behavior of entanglement is as discussed before, due to the coupling J among the waveguides.Each time the states get separable the presence of coupling leads to interaction among the modes of the waveguides and creates back the entanglement.Further we find that the logarithmic negativity E N reaches maximum at times given by the points 2θ = (k + 1)π and is equal to 2r.Thus at this points the state |ξ regains its initial form.

Lossy waveguides
In this section we study the entanglement dynamics of lossy waveguides (γ = 0).The loss γ arises from the loss in the material of the waveguide.In this case the dynamical evolution of the waveguide modes is governed by the full quantum-Louiville equation (3).We next consider the cases of both photon number state and squeezed states at the input of the waveguide and discuss the influence of the loss on their respective entanglement evolution.

Effect of loss on entanglement for photon number states
As discussed above, we first study the effect of loss on the entanglement dynamics of the waveguide modes for photon number input states.For this purpose we consider a single photon input state |1, 1 as the initial state.In this case we can analytically solve the quantum-Louiville equation described in (3).To proceed further, we work in the interaction picture in which the density matrix is ρ(t) = e iJt(a † b+b † a) ρ(t)e −iJt(a † b+b † a) .Then in the interaction picture we can where ã and b are given by Using Eq. ( 26), we can rewrite Eq. ( 25) as For the separable input state |1, 1 , the solution for the density matrix ( 27) can be written as [32] : Further, we can write ρ(t) in terms ρ(t) using the following equation: The decay rates of the modes are given by γ/J = 0.1 (solid black), γ/J = 0.2 (broken black) and γ/J = 0.3 (red).
The above equation gives the time evolution of the density matrix corresponding to the single photon state |1, 1 .Following a similar approach as discussed in section 3, we can evaluate the log negativity for the lossy waveguide case.But the resulting expressions are lengthy and do not exhibit a simple structure.Thus we only give the numerical results for the lossy waveguide case.In Fig. (4), we show the decay of entanglement, as a function of scaled time for the state |1, 1 .
Note that the range of γ/J values studied here are similar to the numerical values used in the experiments [22,23].For example, the coupling parameter J for the lithium niobate waveguide lie between 1.83 × 10 10 sec −1 and 4.92 × 10 10 sec −1 .The loss parameter for these waveguides is close to 3 × 10 9 sec −1 [22] which corresponds to a value of γ/J between 1/7 and 1/20.For AlGaAs waveguides the loss γ is close to 2.7 × 10 10 sec −1 [23].The coupling parameter J for these waveguides is about 2.46 × 10 11 sec −1 .Thus the γ/J value for these waveguides is of the order of 1/10.It is worth mentioning that the γ/J value for silica waveguides is significantly lower than the corresponding values for the lithium niobate and AlGaAs waveguides.
This means that even a small loss would add up to a significant decoherence in these complex quantum systems.From Fig. ( 4) we find that for the lossy waveguide case the entanglement between the waveguide modes decrease with time.In addition, we find that increasing the value of γ/J makes the waveguide modes more fragile, as is evident from Fig. (4).However, we find that the decrease in entanglement is not substantial.Our results indicate that the waveguide system can sustain the entanglement even for the higher decay rates.Thus the coupled waveguide system can be used as an efficient tool for the study of basic quantum optical effects.In addition, the persistence of entanglement suggests that the coupled waveguide system can be used effectively for various applications in quantum information processing [13].For example, the single photon entanglement described here is a key step for the successful implementation of the CNOT gate [12].We also studied the behavior of log negativity for the entangled initial state in the form of two photon NOON state.In this case also we found that the entanglement quantified by E N shows a considerable robustness against the decoherence effect.

Effect of loss on entanglement for squeezed input states
For the initially separable two mode squeezed state |ζ , we find that elements of the covariance matrix σ in presence of loss become dependent on the decay rate γ and is given by, The corresponding symplectic eigenvalue ν of the covariance matrix is then found to be ν± = √ c d ± e .On substituting ν± in equations (15) and using (17) we get the logarithmic negativity-for lossy waveguides.
To study the dependence of entanglement on loss of the waveguide modes we plot the logarithmic negativity E N for different decay rates γ/J in figure (5).As for the case of single photon states we focus on the range of θ important from the experiment point of view.We see new features in the entanglement dynamics as an effect of the loss.We see from figure (5) that in presence of loss the maximum value of entanglement for the state |ζ reduces in comparison to the case of lossless waveguides.However it is important to note that this decrease is not substantial.We further find that with increase in decay rate, the entanglement maximum decreases but does not show considerable reduction (the maximum changes by only 0.4 as the decay rate becomes three times).Thus we see that entanglement is quite robust against decoherence in this coupled waveguide systems.The robustness of entanglement dynamics is an artifact of coherent coupling among the waveguide modes.These findings hence suggest that coupled waveguide can be used as an effective quantum circuit for use in quantum information computations.Further we see another new feature in entanglement in figure (5).We find that there exist an interval of θ during which the state |ζ remains separable.Note that in absence of loss the state |ζ becomes separable momentarily and entanglement starts to build up instantaneously once it becomes zero (see figure 3 decay rate of γ/J = 0.1 and squeezing parameter r = 0.9.We see that entanglement decays slowly with increasing θ as the magnitude of E N diminish successively with every oscillations.In addition periods of disentanglement arises repeatedly in its oscillations.We find that the length of this periods increases with increasing θ .It is worth mentioning here that this kind of behavior has been predicted earlier for two qubit entanglement [33,34].
Next we study the effect of the decay of waveguide mode on the entanglement dynamics of the initial entangled squeezed state |ξ .We find in this case the covariance matrix to be, where f , g , h are given by, f = 1 2 + e −2γt sinh 2 (r); g = − 1 2 e −2γt sinh(2r) sin(2Jt); h = 1 2 e −2γt sinh(2r) cos(2Jt).(34) The symplectic eigenvalues ν± are found to be dependent on the decay rate of the waveguide modes and is given by: The corresponding measure of entanglement given by the logarithmic negativity E N can then be calculated by using equation ( 35), (15) and (17).The behavior of the time dependent logarithmic negativity E N for the state |ξ in presence of loss is found to be similar to that for the separable state |ζ .For brevity we do not show the plot here.Thus as for the separable states, in case of initial entangled input states entanglement is found to be quite robust in the face of loss.
The loss in waveguides that we discussed in this section arises due to material properties like change in refractive index and absorption.On the other hand there can be decay of the waveguide modes in the form of leakage to its surrounding also.It should be noted that leakage is inherently different from the evanescent coupling as the former can arise due to scattering and refraction due to refractive index difference at the waveguide boundaries.Thus the analysis of this section is also valid when the leakage is important as for example is the case when one couples channel waveguides to slab waveguides [3,35].

Conclusion
To conclude, we investigated the time evolution of entanglement in a coupled waveguide system.We quantified the degree of entanglement between the waveguide modes in terms of logarithmic negativity.We have given explicit analytical results for logarithmic negativity in case of initially separable single photon states and for separable as well as entangled squeezed states.We have also addressed the question of decoherence in coupled waveguide systems by considering loss of waveguide modes.For the lossy waveguides we found that the entanglement shows considerable robustness even for substantial loss.Note that our results are based on experimental parameters and thus should be relevant for applications of waveguides in quantum information sciences.Our results serve as guide for experiments dealing with entanglement in waveguide structures.For efficient use of these waveguides, one should choose the waveguide parameter like θ so that one is away from values where the entanglement is minimum.
We are grateful to NSF grant no CCF-0829860 for supporting this research.

Fig. 1 .
Fig. 1. (Color online) Schematic diagram of a coupled waveguide system.The parameter J gives the coupling between the waveguide modes.

Fig. 2 .
Fig. 2. (a) Time evolution of log negatively for the separable input state.The red curve shows the result for |1, 1 state while the black curve shows the result for the |2, 0 state.(b) The behavior of log negatively for the NOON state.The black curve shows the result for two photon NOON state while the red curve shows the result (E N − 1) for the four photon NOON state

Fig. 3 .
Fig.3.Plot of the time dependent logarithmic negativity E N for the state |ζ .Here amount of squeezing is taken to be r = 0.9.

Fig. 4 .
Fig.4.Time evolution of the logarithmic negativity E N in presence of loss of the waveguide modes for the initial separable input state |1, 1 .The decay rates of the modes are given by γ/J = 0.1 (solid black), γ/J = 0.2 (broken black) and γ/J = 0.3 (red).

Fig. 5 .
Fig. 5. (a)Time evolution of the logarithmic negativity E N in presence of loss of the waveguide modes for the input state |ζ .The decay rates of the modes are given by γ/J = 0.1 (solid black), γ/J = 0.2 (broken black) and γ/J = 0.3 (red).Here the squeezing is taken to be r = 0.9.The loss leads to new behavior in the entanglement.(b) Long time behavior of the logarithmic negativity in presence of loss of the waveguide modes for the the state |ζ .