# Tripartite and bipartite entanglement in continuous-variable tripartite systems

## Abstract

We examine one asymmetric adnd two fully symmetric Gaussian continuous-variable systems in terms of their tripartite and bipartite entanglement properties. We treat pure states and are able to find analytic solutions using the undepleted pump approximation for the Hamiltonian models, and standard beamsplitter relations for a model that mixes the outputs of optical parametric oscillators. Our two symmetric systems exhibit perfect tripartite correlations, but only in the unphysical limit of infinite squeezing. For more realistic squeezing parameters, all three systems exhibit both tripartite and bipartite entanglement. We conclude that none of the outputs are completely analogous to either GHZ or W states, but there are parameter regions where they produce T states introduced by Adesso et al.The qualitative differences in the output states for different interaction parameters indicate that continuous-variable tripartite quantum information systems offer a versatility not found in bipartite systems.

###### pacs:

42.50.-p,42.50.Dv,42.65.Lm,03.65.Ud## I Introduction

In the theory of discrete variable tripartite entanglement, the two most famous states are known as the GHZ (Greenberger-Horne-Zeilinger) state (1) and the W state (2). One essential difference between these two is seen when one mode is traced over and any resulting bipartite entanglement is looked for. The reduced GHZ state then becomes completely separable, while the reduced W state will demonstrate remnant bipartite entanglement. The concept of these two states when transferred to the Gaussian continuous-variable (CV) quantum information (3) has been extensively analysed by Adesso et al. (4); (5), who also introduced the T states, which exhibit tripartite entanglement only.

In this paper we consider three processes which are known to produce CV entanglement and show analytically that they also produce bipartite entanglement in some, but not all, of the parameter space. We detect the bipartite entanglement using the Duan-Simon measure (6); (7), which is both necessary and sufficient for Gaussian systems. The fact that both bipartite and tripartite entanglement are exhibited suggests that they are none of the above states over the entire operating regime. Since the bipartite entanglement is not maximal, they are not proper W states, and because the bipartite entanglement exists, they are neither GHZ nor T states. There are also operating regimes where only tripartite entanglement exists, although this can only be perfect in the limit of infinite squeezing. As this limit is not physical, they necessarily lack one of the characteristics of the true GHZ state. Since no bipartite entanglement exists in these regions, the outputs then qualify as T states.

We write the Hamiltonians in the non-depleted pump approximations for the two nonlinear systems (8); (9) and use standard beamsplitter relations for the other (10). The solutions in terms of expectation values of the second moments of the quadrature operators allow us to find analytical expressions for the Duan-Simon (6); (7) and Reid Einstein-Podolsky-Rosen (EPR) (11); (12) correlations used to denote bipartite entanglement and EPR-steering (13). We also calculate the correlations developed by van Loock and Furusawa (14) for the detection of tripartite entanglement and the three-mode EPR-steering correlations developed by Olsen et al. (15).

Our results, showing that continuous-variable systems for the production of tripartite entanglement will behave qualitatively differently as the input fields are varied, add a dimension to continuous-variable quantum information which is not present in the discrete variable version. We will show that, by simple tuning of the inputs, different classes of output states are produced. This tunability brings a versatility to these systems that is not present in the discrete variable equivalents, and opens new possibilities for quantum information technologies such as quantum key distribution (16).

## Ii Common classes of tripartite state

In the discrete variable regime, there are two classes of tripartite entangled states of three qubits that are commonly considered. The first of these is the GHZ state, introduced by Greenberger et al. (1), commonly represented as

(1) |

This state possesses maximal tripartite entanglement and no remnant bipartite entanglement whatsoever whenever any one of the qubits is traced over. This state also gives a maximal Bell violation.

The canonical representative of the W states is written as

(2) |

which also exhibits maximal tripartite entanglement but also exhibits remnant bipartite entanglement if any one mode is traced over (2). Both of these states as written are pure and symmetric.

When we consider continuous-variable tripartite states, the situation becomes somewhat different. The systems demonstrated by Smithers (8) and Aoki (10) will exhibit the perfect correlations of a GHZ state in the regime where they become perfect quadrature eigenstates. Unfortunately, this limit is not obtainable in practice. In the Aoki scheme, for example, it would require perfect squeezing, which is unphysical. In the version of the Smithers scheme that we will analyse here, it would require infinite interaction strength, which is also not a physically relevant concept. However, Adesso and Illuminati (4) have addressed this problem, dividing CV tripartite states into five separate classes. These range from states which are totally inseparable under any of the three possible bipartitions to those which are separable under all possible bipartitions. In terms of the entanglement properties, they add a type to the GHZ-types and W-types commonly in use. For a state which exhibits tripartite entanglement only, but without being in the GHZ limit, they introduce the nomenclature T state. We will use this nomenclature in this article.

## Iii The inequalities

### iii.1 Bipartite measures

We will evaluate bipartite entanglement and EPR-steering in terms of the functions of quadrature operators developed by Duan et al. (6), Simon (7), and Reid (12), which are appropriate measures for two-mode optical systems. The quadrature operators are defined as and . This allows us to define the Duan-Simon inequalities as

(3) |

with and being mode indices. Violation of either of these is a demonstration of bipartite entanglement. We will call the first of these combined variance sums and the second . Because all the systems we consider are Gaussian and pure, these entanglement correlations are both necessary and sufficient for this demonstration (17).

The EPR paradox is detected by the well-known criteria developed by Reid (12), in terms of inferred quadrature variances,

(4) |

The inferred variances are defined, with the value of being inferred from measurements of (and similarly for ), as,

(5) |

from which we immediately see that there is an implied asymmetry since we can equally define , swapping the roles of the people measuring each mode. In some circumstances the values measured at can be inferred from measurements on mode , but not vice-versa. This was first predicted, in sum frequency generation, by Olsen and Bradley (18) and has recently been further analysed by Ji et al. (19). In what follows we will label the product , inferred from and , as .

### iii.2 Tripartite

The van Loock-Furusawa conditions (14) give a set of inequalities

(6) |

for which the violation of any two demonstrates tripartite entanglement. The , which are arbitrary and real, can be optimised (20), using the variances and covariances, as

(7) |

which is the process we follow with the results presented below.

Another set of inequalities was also presented by van Loock and Furusawa, the violation of any one of which is sufficient to prove tripartite entanglement,

(8) |

The generalisation of the Reid EPR inequalities to three modes by Olsen et al. (15) involves using either one mode to infer combined properties of the other two, or combined properties of two of the modes to infer properties of the third mode. We define

(9) |

with a demonstration of the paradox requiring

(10) |

When this is satisfied for , and , we have established tripartite entanglement. We can also use the inferred variances of the combined modes

(11) |

with a demonstration of the paradox when

(12) |

In the interests of brevity we will label these two correlations (Eq. 10) and (Eq. 12). As above, a demonstration for the three possible combinations establishes tripartite entanglement. In the language of EPR-steering introduced by Wiseman et al. (13), a demonstration via means that two of the participants have combined to steer the third. A demonstration via means that one participant can steer the combined properties measured by the other two, without steering either of them individually.

## Iv A symmetric model from a single Optical Parametric Amplifier

This model consists of triply concurrent downconversion (21), with the intracavity version being analysed by Bradley et al. (9), where it was noted that the state created tended towards a GHZ state in the limit of infinite squeezing, but was analogous to a W state for finite squeezing. The interaction Hamiltonian in the undepleted pump approximation is written as

(13) |

where represents the product of the optical nonlinearity and the corresponding pump fields. Note that we have set the two pump fields as equal. Setting

(14) |

the Heisenberg equations of motion for the quadrature operators are solved as (15)

(15) |

from which we can calculate the second order moments necessary for the Duan-Simon, Reid EPR, and van Loock-Furusawa (VLF) correlations. We find

(16) |

all of which assume vacuum in these modes at .

In Fig. 1 we show the results for the various correlations. Interestingly, despite the fact that we know this system displays full tripartite entanglement, we see that the DS detect a degree of bipartite entanglement at early times. If the state produced were analogous to a continuous-variable GHZ state, there would be no possibility of bipartite entanglement (2). In fact, if we operate on the three-mode vacuum state with the interaction Hamiltonian, we find to first order

(17) |

which has more in common with the W states (2); (22), so that a degree of bipartite entanglement is therefore not ruled out. In fact we can see that tracing over any one mode leaves bipartite entanglement for both and . We found no evidence of bipartite EPR-steering in this system, with the correlation never being less than one. Fig. 1 also shows the results for the two van Loock-Furusawa correlations and the EPR-steering correlation . All of these are equal under any exchange of indices for this system and is an identical shape to , but begins at a value of one. We see that the system begins to exhibit tripartite entanglement and EPR-steering as soon as it is turned on, but that the correlation fails to detect this after some time. In the region where DS, the output is properly caled a T state (4); (5). It becomes more analogous to a GHZ state (which would require ) as increases, but is only totally equivalent in the unphysical limit of infinite interaction time .

## V Three OPOs and a two beamsplitters

One of the pioneering results for continuous variable tripartite entanglement came from van Loock and Braunstein (23), and was implemented experimentally by Aoki et al. (10), who mixed three squeezed states on two beamsplitters to obtain three entangled output beams, as shown in This setup was subsequently analysed in terms of tripartite entanglement and EPR-steering in both the time and frequency domains by Olsen et al. (15) and is a subset of the systems recently analysed by Wang et al. (24). Bipartite measures were not analysed in these articles. The system uses three optical parametric oscillators (OPO), with the first, OPO, producing a state squeezed in the quadrature, while the other two produce squeezed states. The annihilation operator represents the output of OPO. The output of OPO and OPO are mixed on the first beamsplitter, BS, to produce outputs represented by and . The field corresponding to is then mixed with on BS. The outputs of BS are represented by and . With the squeezed inputs, tripartite entanglement is found between the three outputs. In Ref. (15), the van Loock Furusawa correlations were calculated analytically, but without optimisation. The tripartite EPR-steering correlations were also calculated. We will now calculate the correlations with optimisation, and examine both bipartite and tripartite entanglement in this system.

Assigning BS a reflectivity of and BS a reflectivity of , we find the solutions for the in terms of the inputs as

(18) |

which allow us to find all the correlations we require for the bipartite and tripartite correlations we wish to calculate. We note here that, although these expressions appear asymmetric, they become fully symmetric for and , and these are the values we use in our final results. The required variances are

(19) |

and the covariances are

(20) |

from which we have all that is necessary to calculate the Duan-Simon, van Loock Furusawa and EPR-steering correlations. For and as in Aoki et al. (10), the and three-mode EPR correlations are given in Ref. (15). However, possible bipartite entanglement was not analysed in that work, nor were the optimised using the , so we will give these results here.

For a squeezing parameter , equal for each OPO, we may assume minimum uncertainty squeezed states and set

(21) |

which leads to the bipartite correlations

(22) |

of which falls below over a range of . For these parameters, we do not see a demonstration of bipartite EPR-steering.

The optimised are found as

(23) |

and the are

(24) |

with these not changing under permutations of the indices. Note that, with optimisation, the begin at , rather than at the non-optimised value of found in Ref. (15). For completeness we note that

(25) |

so that the two types of tripartite EPR-steering become available as soon as is greater than zero.

We see from Fig. 2 that the V once again fail to detect entanglement in a parameter regime where it is found by the V. In the Aoki experiment (10), a V of approximately was measured. The fact that we found no bipartite EPR-steering is consistent with the result of Wang et al. (24), who showed that, in an -mode system of this type, with beamsplitters and sources, at least participants must combine to steer any single participant. Fig. 2 shows that the experimental result of Ref. (10) had not entered the T state regime, which begins for V. In the experiment, both bipartite and tripartite entanglement would have been available, with an increase in squeezing needed for the bipartite entanglement to disappear. As with the scheme of section IV, the system tends towards a state with true GHZ properties only in the limit of large squeezing. In fact, the behaviours of the two systems in terms of quantum correlations are very similar, which is as expected since both are examples of fully symmetric Gaussian systems.

## Vi An asymmetric model

Our asymmetric system, which combines downconversion with sum-frequency generation, was first proposed by Smithers and Lu (8), and theoretically analysed in both travelling wave (25) and in an intracavity configuration by Yu et al. (26). The configuration was subsequently analysed in more depth by Pennarun et al. (27), who investigated the stability properties and predicted tripartite entanglement in different regimes. It consists of a nonlinear medium pumped at frequency . The downconversion part of the process, denoted by the effective nonlinearity , generates two fields at and , where . The pump field at can then combine with the field at in a sum frequency generation process (18), to produce a further field at , with effective nonlinearity . In Yu et al. (26) the nonlinear medium is a quasiperiodic superlattice. We will use the annihilation operators to describe the fields at for . If we consider that the pump field is intense and classical so that depletion does not become important, we may write the interaction Hamiltonian as

(26) |

In this case, the represent products of the actual nonlinearity multiplied by the amplitude of the pump field.

We find the Heisenberg equations of motion for the annihilation operators as

(27) |

with those for the creation operators being the Hermitian conjugates. These may be solved analytically. Setting for , we find

(28) |

We note here that these are different to the solutions given by Ferraro et al. (25), who worked in the regime where . They have been given previously by Olsen and Bradley (28), who also calculated the non-optimised VLF measures of Eq. 6, but did not investigate bipartite entanglement. This immediately allows us to write solutions for the quadrature operators, which then allows us to find expressions for all the entanglement and EPR-steering correlations of section III. Setting

(29) |

we find the moments required for the variances and covariances as

(30) |

In Fig. 3 we show that both tripartite and bipartite entanglement are predicted over a range of interaction strength, , for . The Duan-Simon measure, DS, shows that modes and are entangled, while either of V or V, in the region where they are less than , demonstrate tripartite entanglement. Fig. 4 shows the bipartite EPR-steering correlations between modes and , demonstrating that these two modes are able to steer each other. There was no violation of either the bipartite entanglement or EPR-steering inequalities for the pairs and .

When we investigate the three-mode EPR-steering correlations, shown in Fig. 5 and Fig. 6, we see that only fails to violate the inequality, although the violation by is minimal. This means that, for these parameters, modes and cannot be used to steer mode , although all the other combinations are possible over some range of interaction strengths. This is quite different from the symmetric case in the previous section, where the correlations were equivalent under any change of indices. We also draw attention to the fact that the system does not enter the T state regime within the range of interaction strength shown here, although it will for longer interaction times.

## Vii Conclusions

We have investigated two symmetric and one asymmetric systems known to produce tripartite entangled outputs, in terms of both the bipartite and tripartite entanglement available. We found that all three systems produce both tripartite and bipartite entanglement in some operating regimes. In these regimes the outputs may be thought of as W-type states, even though they do not satisfy the criteria completely. Only the asymmetric system was found to produce bipartite EPR-steering. All three systems have operating regimes where only tripartite entangled outputs are introduced, where the outputs qualify as T states. In the limit of large interaction strength or squeezing, the two symmetric systems produce states which may truly be thought of as having GHZ properties. However, this limit is unphysical because either the undepleted pump approximation or energy conservation breaks down long before it is reached. Although the labelling of these states as GHZ or GHZ type is common, some care should be taken with this since they do not satisfy all the criteria. For example, in the physically attainable regimes they do not produce eigenstates of quadrature combinations and thus will not give a yes or no reply to the question of whether tripartite entanglement is present. The answer they give arises from a statistical violation of the entanglement criteria.

On a final note, the production of qualitatively different quantum states from the same apparatus, obtained by changing the operating parameters, may be of advantage to quantum information experimentalists. There may be situations where changing the input laser intensities, for example, and moving from a W type state to a T state is advantageous. As was shown in Ref. (16), there are operating regimes of the asymmetric system where two of the participants can practise one-sided device independent quantum key distribution which cannot involve the third. We expect that there will be other applications which take advantage of the flexibility we have demonstrated in this article.

## Acknowledgments

This research was supported by the Australian Research Council under the Future Fellowships Program (Grant ID: FT100100515).

### References

- D.M. Greenberger, M.A. Horne, A. Shimony, and A. Zeilinger, Am. J. Phys. 58, 1131 (1990).
- W. Dur, G. Vidal, and J.I. Cirac, Phys. Rev. A62, 062314 (2000).
- S.L. Braunstein and P. van Loock, Rev. Mod. Phys. 77, 513 (2005).
- G. Adesso, A. Serafini, and F. Illuminati, Phys. Rev. A73,032345 (2006).
- G. Adesso and F. Illuminati, New. J. Phys. 8, 15 (2006).
- L.-M. Duan, G. Giedke, J.I. Cirac and P. Zoller, Phys. Rev. Lett. 84, 2722 (2000).
- R. Simon, Phys. Rev. Lett. 84, 2726 (2000).
- M.E. Smithers and E.Y.C. Lu, Phys. Rev. A10, 1874 (1974).
- A.S. Bradley, M.K. Olsen, O. Pfister and R.C. Pooser, Phys. Rev. A72, 053805 (2005).
- T. Aoki, N. Takei, H. Yonezawa, K. Wakui, T. Hiraoka, A. Furusawa, and P. van Loock, Phys. Rev. Lett. 91, 080404 (2003).
- A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 777 (1937).
- M.D. Reid, Phys. Rev. A40, 913 (1989).
- H.M. Wiseman, S.J. Jones and A.C. Doherty, Phys. Rev. Lett. 98, 140402 (2007).
- P. van Loock and A. Furusawa, Phys. Rev. A67, 052315 (2003).
- M.K. Olsen, A.S. Bradley and M.D. Reid, J. Phys. B 39 2515 (2006).
- M.K. Olsen and E.G. Cavalcanti, Phys. Rev. Lett. submitted.
- R.Y. Teh and M.D. Reid, Phys. Rev. A90, 062337 (2014).
- M.K. Olsen and A.S. Bradley, Phys. Rev. A77, 023813 (2008).
- S.-W. Li, M.S. Kim, and H. Nha, J. Phys. A 48, 135301 (2015).
- M.K. Olsen and A.S. Bradley, Phys. Rev. A74, 063809 (2006).
- O. Pfister, S. Feng, G. Jennings, R. Pooser, and D. Xie, Phys. Rev. A70, 020302 (2004).
- B. Fortescue and H.-K. Lo, Phys. Rev. Lett. 98, 260501 (2007).
- P. van Loock and S.L. Braunstein, Phys. Rev. Lett. 84, 3482 (2000).
- M. Wang, Y. Xiang, Q. He, and Q. Gong, Phys. Rev. A91, 012112 (2015).
- A. Ferraro, M.G.A. Paris, M. Bondani, A. Allevi, E. Puddu, and A. Andreoni, J. Opt. Soc. Am. B6, 1241 (2004).
- Y.B. Yu, Z.D. Xie, X.Q. Yu, H.X. Li, P. Xu, H.M. Yao, and S.N. Zhu, Phys. Rev. A74, 042332 (2006).
- C. Pennarun, A.S. Bradley, and M.K. Olsen, Phys. Rev. A76, 063812 (2007).
- M.K. Olsen and A.S. Bradley, J. Phys. B 39, 127 (2006).