Tests of multimode quantum non-locality with homodyne measurements

# Tests of multimode quantum non-locality with homodyne measurements

Antonio Acín ICFO-Institut de Ciencies Fotoniques, Mediterranean Technology Park, 08860 Castelldefels (Barcelona), Spain    Nicolas J. Cerf QuIC, Ecole Polytechnique, CP 165, Université Libre de Bruxelles, 1050 Brussels, Belgium    Alessandro Ferraro ICFO-Institut de Ciencies Fotoniques, Mediterranean Technology Park, 08860 Castelldefels (Barcelona), Spain    Julien Niset QuIC, Ecole Polytechnique, CP 165, Université Libre de Bruxelles, 1050 Brussels, Belgium
July 11, 2019
###### Abstract

We investigate the violation of local realism in Bell tests involving homodyne measurements performed on multimode continuous-variable states. By binning the measurement outcomes in an appropriate way, we prove that the Mermin-Klyshko inequality can be violated by an amount that grows exponentially with the number of modes. Furthermore, the maximum violation allowed by quantum mechanics can be attained for any number of modes, albeit requiring a quantum state that is rather unrealistic. Interestingly, this exponential increase of the violation holds true even for simpler states, such as multipartite GHZ states. The resulting benefit of using more modes is shown to be significant in practical multipartite Bell tests by analyzing the increase of the robustness to noise with the number of modes. In view of the high efficiency achievable with homodyne detection, our results thus open a possible way to feasible loophole-free Bell tests that are robust to experimental imperfections. We provide an explicit example of a three-mode state (a superposition of coherent states) which results in a significantly high violation of the Mermin-Klyshko inequality (around ) with homodyne measurements.

## I Introduction

The non-local nature of quantum mechanics is unarguably one of its most counterintuitive aspects, marking a fundamental departure from the classical picture of physical systems. First recognized by Einstein, Podolsky, and Rosen EPR (), the intrinsic non-locality of the quantum theory was put into a testable form via the celebrated Bell inequalities Bell (). Since the first experimental violation of a Bell inequality Aspect (), a great variety of tests have been performed. In all of them, a local realistic model is shown to be incompatible with the actual experimental observations. Nonetheless, a fully conclusive evidence for rejecting local realism has not been achieved yet, for all the experiments performed so far suffer from loopholes, forcing us to rely on supplementary assumptions in order to reject local realism. Specifically, a conclusive experiment should close two main loopholes: the locality loophole (namely, the measured correlations must be collected with spacelike separated events) and the detection-efficiency loophole (namely, the proportion between detected and undetected events has to be high enough for the data to be fully representative of the whole ensemble without the need to assume a “fair sampling”).

There is a large variety of quantum systems for which a test of local realism may be envisaged. However, the quest for a loophole-free Bell test has recently focused the research towards experiments involving propagating light modes measured with homodyne detectors Munro (); Wenger (); Patron (); Nha (). The advantage of such a “continuous-variable” approach is twofold. First, light modes can be easily sent to spacelike separated detectors suffering only a tolerable degree of decoherence. Second, the current technology of homodyne detectors achieves a degree of detection efficiency high enough to potentially close the detection efficiency loophole. On the other hand, such an approach involves also drawbacks whose resolution is still challenging. The main issue comes from the fact that the results of homodyne measurements can be described by means of the Wigner function. Thus, in order to avoid a local hidden variable description of the measured correlations, one should necessarily perform the test with a state endowed with a non-positive Wigner function (otherwise, the Wigner function is a genuine probability distribution, which provides an explicit local realistic model of the data). Even if single-mode states of traveling light have already been generated with a non-positive Wigner function NGaus_exp1 (); NGaus_exp2 (), such (non-Gaussian) states are hard to master in the laboratory. Among the above mentioned proposals of Bell tests with homodyne detection, those nearer to an experimental realization only give a small violation of Bell inequalities Patron (); Nha (), whereas higher violations involve states whose actual generation seems rather unrealistic Wenger (). In this context, the search for Bell tests involving feasible resources and giving, at the same time, violations that are high enough to be robust against experimental noise is very desirable. This is the thrust of the present paper.

Bell tests relying on homodyne detection typically involve the following scenario. Two parties perform spacelike separated homodyne measurements by randomly choosing each between two settings, thus measuring two quadratures of the incoming electromagnetic fields. The violation of the Bell Clauser-Horne-Shimony-Holt (CHSH) inequality CHSH () is then tested. Since this inequality is devised for two-outcome measurements, the collected data, which are distributed (ideally) in a continuous way, have to be discretized; such a procedure is referred to as binning process. In this work, we will generalize this scenario to more than two parties, considering the Mermin-Klyshko (MK) type of Bell inequalities, which involves parties, two measurement settings, and two outcomes Mermin_Klyshko (). Our main motivation is that higher violations of local realism are expected as the number of parties increases GBP98 (), which may help in the search for feasible proposals of loophole-free Bell tests able to tolerate the experimental noise.

The earlier investigations of Bell tests in infinite-dimensional spaces with multipartite settings have considered several scenarios. In Ref. VLoock (), a test based on the measurement of the light field parity was envisaged. It was found that the violation does not increase exponentially, as one would have expected, but this may be due to the fact that no optimization over the possible measurement settings was performed job (). In Ref. pseudospin () instead, a maximal violation of Mermin-Klyshko inequality was found for continuous variable states, considering the measurement of a different class of operators which can be seen as the continuous-variable analogue of the spin operators. However, both these approaches deal with non-Gaussian measurements described by a non-positive Wigner function, which are far from the reach of current detection technology. Recently, another approach has been introduced that does not rely on the use of inequalities with discrete outcomes, thus avoiding the need of a binning procedure Cavalcanti (). There, quadrature measurements via homodyne detection are considered, and the measurement outcomes are used to directly test a novel Bell inequality for continuous variables. Interestingly, the authors find that the violation of the local realistic bound is exponential in the number of parties. However, a possible experimental implementation is still very challenging as it would require at least ten space-like separated homodyne measurements. It is then unclear whether such a novel approach can give advantages, from a practical perspective, over non-locality tests involving binning strategies.

In this paper, we will stick with the use of discrete-variable Bell inequalities in the tests of non-locality for continuous-variable states, thus using some binning procedure. We will address the following question: Is it possible to have an exponential increase of the violation of local realism in a test involving quadrature measurements of modes and considering the Mermin-Klyshko inequalities? We will answer this question by the affirmative, providing specific examples of states exhibiting such a behavior. These states belong to the class of photon-number correlated states, and generalize to an arbitrary number of modes the approach of Ref. Munro (). Remarkably, we will show that it is even possible to reach the maximal violation of the Mermin-Klyshko inequalities for modes by properly choosing the binning procedure and the quantum state, thereby extending the results obtained in Ref. Wenger () for two modes.

The theoretical interest of our results resides in the fact that we are not considering a direct mapping from continuous- to discrete-variable states. Thus, the possibility to obtain a maximal violation of Mermin-Klyshko inequality for any number of parties is a priori non-trivial. Furthermore, we may have anticipated that the binning procedure causes an irreparable loss of information, preventing the possibility to reach a maximum violation. From an experimental perspective, our results imply that various strategies allow for an increase of the violation of locality in homodyne-based Bell tests. This, in turn, gives specific insights in the search of multipartite states appropriate for an actual experiment.

The paper is organized as follows. In Sec. II, we analyze a Bell test involving multimode states with perfect correlations in the number of photons in all modes. We analytically show that the violation of local realism can grow exponentially with the number of modes involved. This is true despite the fact that a simple sign binning strategy is used. Nonetheless, we see that the resulting violation does not reach the maximum value achievable within quantum mechanics. However, we show in Sec. III that a different binning strategy, properly tailored for some specific states, allows attaining the maximal violation of locality achievable within quantum mechanics. Then, in Sec. IV, we address the robustness of these tests of nonlocality. We see that, at least for the analyzed noise model, the tolerable noise increases with the number of modes, making the violation better testable in practice. In Sec. V, we exploit all these results from a practical oriented perspective. In particular, we give an explicit example of a class of three-mode states which, on the one hand, exhibit a gain in the violation by considering more than two modes, and, on the other hand, may be implemented with the near future technology. Finally, we close the paper in Sec.VI with some concluding remarks.

## Ii Correlated photon number states

We show in this Section that in order to obtain an exponential violation of local realism with quadrature measurements, it is possible to rely on a simple binning strategy and the standard Mermin-Klyshko Bell inequalities. In achieving this goal we generalize to parties the work of Ref. Munro (), where the bipartite scenario was analyzed.

Before proceeding, let us recall the general form of the Mermin-Klyshko Bell inequalities. Consider two dichotomic observables and for each party . The Mermin-Klyshko inequalities are based on the recursive definition of the Bell operator

 Bt≡Bt−12⊗(Ot+O′t)+B′t−12⊗(Ot−O′t), (1)

where , , and denotes the same expression as but with all ’s and ’s exchanged GBP98 (). The Mermin-Klyshko inequality for parties then reads

 Bm≡|⟨Bm⟩|≤2. (2)

In the case of two parties (), this expression reduces to the well-known CHSH inequality, for which

 B2=O1⊗O2+O1⊗O′2+O′1⊗O2−O′1⊗O′2 (3)

Let us now consider a generic photon-number correlated state of bosonic modes

 |Ψ⟩=∞∑r=0cr|r⟩1|r⟩2...|r⟩m, (4)

with . The two observables and to be measured are chosen as the two quadratures of the electromagnetic field and (corresponding to angles and ). These can be measured by applying an homodyne detection on each mode . In the following, we use the notation for the quadrature , and for the outcome of its measurement (and similarly for primed quantities).

The joint probability to obtain the outcomes when measuring the quadratures is given by

 P(x1,...,xm)= |1⟨x1|...m⟨xm|Ψ⟩|2=∞∑r,s=0crc∗s ×eiϕ(r−s)(π2r+sr!s!)m/2m∏t=1e−x2tHr(xt)Hs(xt), (5)

where , is the Hermite polynomial of degree , and are the eigenvectors of the quadrature operator .

Consider the following simple binning strategy, known as sign binning: when the result of a quadrature measurement falls in the domain the value is associated to it. Now, let us start by calculating the probability that a “” result is observed in all the measuring sites:

 P+1,...,+1 = ∫∞0dx1⋯∫∞0dxmP(x1,...,xm) (6) = ∞∑r,s=0crc∗seiϕ(r−s)(π2r+sr!s!)m/2 ×m∏t=1∫∞0dxte−x2tHr(xt)Hs(xt).

The integrals above can be evaluated by recalling the following properties of Hermite polynomials for ,

 ∫∞0dxe−x2Hr(x)Hs(x)=π2r+sr−s[F(r,s)−F(s,r)], (7)

where we defined as

 F(r,s)−1=Γ(12−12r)Γ(−12s), (8)

with being the Gamma function. For , one has instead

 ∫∞0dxe−x2H2r(x)=2r−1r!√π, (9)

Defining the functions

 G(ϕ,m) = 2∑r>sRe(crc∗s)gr,s(ϕ,m), gr,s(ϕ,m) = (π2r+sr!s!)m/2[F(r,s)−F(s,r)r−s]m (10) ×cos[ϕ(r−s)],

one obtains that

 P+1,...,+1=12m+G(ϕ,m). (11)

The other probabilities can be obtained in a similar way. Let us define the multi-index , with denoting the measurement outcome obtained for mode after binning. Then, the joint probability for a generic collection of measurement outcomes will be indicated by . It can be obtained similarly as by recalling that an Hermite polynomial of even (odd) degree is an even (odd) function, namely

 Pd=12m+σ(d)G(ϕ,m), (12)

where . Now, we are in the position to calculate the generic correlation function between the measurement results (note that it only depends on the sum of the angles ). By definition, we have

 E(ϕ,m)=∑dσ(d)Pd, (13)

where the sum goes over all possible collections of measurement outcomes. Since the number of possible collections for which is equal to that for which , one finally has, by substituting Eq. (12) into Eq. (13), that

 E(ϕ,m)=2mG(ϕ,m). (14)

Let us show that we can reach an exponential violation with a simple analytically manageable example. For the simple case of three parties (), the Mermin-Klyshko inequality then reads , with

 B3= O1⊗O2⊗O′3+O1⊗O′2⊗O3+O′1⊗O2⊗O3 −O′1⊗O′2⊗O′3 (15)

Considering a tripartite GHZ state,

 |GHZ3⟩=1√2(|000⟩+|111⟩), (16)

we have and , so that and Eq. (14) becomes

 E(ϕ,3)=23G(ϕ,3)=(2/π)3/2cos(ϕ). (17)

The GHZ-like angles (, , , and ), give the maximum violation of the inequality, namely

 B3=|3E(π,3)−E(0,3)|=4(2π)3/2≃2.032. (18)

Now, consider the multipartite generalization of the GHZ state

 |GHZm⟩=1√2(|0...0⟩+|1...1⟩), (19)

Eq. (14) becomes

 E(ϕ,m)=(2/π)m/2cosϕ. (20)

The dependence on the angle of the above correlations is the same as the one appearing in a standard spin-like test for a multipartite GHZ state, namely . In that case, it is known that the choice of GHZ-like angles and , gives the highest value of the Bell factor, namely (see e.g. Ref. GBP98 ()). Therefore, using the same angles, the corresponding Bell factor reads

 Bm=√2(4/π)m/2, (21)

giving rise to an exponential violation of local realism.

Apart from this simple analytical example, one can use a numerical approach to show the exponential violation of local realism as formula (14) is easily amenable to perform numerical calculations for a fixed number of parties . In order to find the state (coefficients ’s) that maximally violates the Mermin-Klyshko inequality, one has to evaluate the corresponding Bell factor for a given configuration of measuring angles. The Bell factors are expressed in general by a linear combination of correlation functions given each by Eq. (14), with the prescription given in Eq. (1).

Let us search the state which maximizes the violation of the Mermin-Klyshko inequality for the particular GHZ-like choice of angles. For three modes (), defining the (infinite dimensional) real symmetric matrix as

 [B3]r,s=23(3gr,s(π,3)−gr,s(0,3)), (22)

with the diagonal elements being set to zero, we note that the Bell factor can be re-expressed as , where the elements of the vector are given by the coefficients of the input state, i.e. . Consequently, the maximal violation of the Mermin-Klyshko inequality is simply given by the maximal eigenvalue of the matrix , while the optimal input state is determined by its corresponding eigenvector. In order to perform a numerical analysis, one has to truncate the Hilbert space dimension of to some arbitrary . For example, for the optimal choice turns out to be the GHZ state (16), giving a violation of . By increasing the dimension , the asymptotic violation is given by . In Fig. 1, we show the coefficients for the optimal state in the case , for which the Bell factor is .

The same procedure can be applied for any number of parties. In the case of two parties (), one recovers the results given by Munro in Ref. Munro () provided that the constraint is taken into account, namely . Interestingly, a higher violation can be achieved if we consider negative coefficients for . As an example, we report in Fig. 2 the coefficients for the optimal state in the case , for which the Bell factor raises up to (recall that the Bell factor can be written in this case as , where we have chosen , as in Munro ()).

## Iii Maximal violation

Let us now see if we can find a class of states and a binning strategy that allows for a maximal violation of the Mermin-Klyshko inequality for any number of parties and using quadrature measurements. We recall that the maximal quantum violation of these inequalities is given by GBP98 (). Inspired by the results of Ref. Wenger (), we introduce the state

 |Ψ⟩=1√2(|f⟩⊗m+eiθ|g⟩⊗m), (23)

where is a real and even function of some quadrature , while is real and odd (thus, and are orthogonal). Both functions are also normalized to unity. Note that because is real and even, it has a real and even Fourier transform , while has an imaginary Fourier transform , with real and odd.

Each party chooses to measure one of two conjugated quadrature via homodyne detection, either or , and obtains a continuous variable or depending on the choice of measurement setting. Suppose the first parties measure the quadrature, while the remaining measure the conjugate quadrature . The joint probability that they obtain the results is given by

 P(x1,... ,xk,pk+1,...,pm)=|⟨x1|...⟨xk|⟨pk+1|...⟨pm|Ψ⟩|2 = 12(f2(x1)...f2(xk)~f2(pk+1)...~f2(pm) +g2(x1)...g2(xk)~h2(pk+1)...~h2(pm) +2cos(θ+(m−k)π2)×f(x1)g(x1)...f(xk)g(xk) ×~f(pk+1)~h(pk+1)...~f(pm)~h(pm)). (24)

To exploit the parity properties of and , we introduce the root binning defined in Ref. Wenger (). This binning depends on the roots of the known functions and . If party measures the quadrature, the result will be interpreted as a “” if the measured value lies in the interval where and have the same sign, and “” if their signs are opposite, i.e. we consider the following domains

 D+x ={x∈R|f(x)g(x)≥0} D−x ={x∈R|f(x)g(x)<0}. (25)

We can similarly define the domains and associated to the measurement of the quadrature . For the choice of measurement settings defined above, we can thus calculate probability corresponding to the observation of a given collection of binary results. For example, the probability that each party observes a “” result reads

 P+1,...,+1=∫D+x dx1...∫D+xdxk∫D+pdpk+1...∫D+pdpm ×P(x1,...,xk,pk+1,...,pm). (26)

We are now in the position to calculate the correlation function . Note that since and ( and ) are even and odd respectively, and ( and ) are even functions. Hence the first two terms of the right hand side of (III) are even functions also, and their contribution to the correlation function will vanish. We thus obtain the remarkably simple expression

 E(X1,⋯Xk,Pk+1,⋯Pm) =VkWm−k ×cos[θ+(m−k)π2] (27)

where

 V =∫∞−∞dx|f(x)g(x)| W =∫∞−∞dp|~f(p)~h(p)| (28)

Interestingly, the correlation function (III) only depends on the number of sites where and are measured. All correlation functions corresponding to measurements of the quadrature, and measurements of the quadrature are therefore equal. We will denote them to emphasize this property.

Let us illustrate the power of this compact notation with an example. For , the Bell factor reads

 B3= |E(X1,X2,P3)+E(P1,X2,X3)+E(X1,P2,X3) −E(P1,P2,P3)| = |3E(2,1)−E(0,3)| = |3V2Wcos(θ+π2)−W3cos(θ+3π2)| (29)

We see that the maximal violation, i.e. , can be reached with a state such that and . Although such a state is quite unrealistic, one can define a family of physical states that approximates it arbitrarily well. The corresponding and functions are trains of gaussians, and as the number of peaks goes to infinity. We refer the reader to Ref. Wenger () for their exact analytical expression.

Let us now generalize this result for an arbitrary . First note that the Bell factor (2) can be written as

 Bm=12|⟨XmBm−1⟩+⟨PmBm−1⟩ +⟨XmB′m−1⟩ −⟨PmB′m−1⟩| (30)

with and . In order to benefit from our compact notation, we explicitely expand the expectation values of and in terms of correlation functions

 ⟨Bm−1⟩ =m−1∑k=0αkE(k,m−1−k) ⟨B′m−1⟩ =m−1∑k=0αkE(m−1−k,k) (31)

where the ’s are some known coefficients. When for example, we have , , and . As the correlation function only depends on the number of and measurements, the average values of the four operators of (III) can be easily calculated from and . Suppose has a term proportional to the operator, which leads to the correlation function . The operator will thus have a term proportional to leading to the correlation function , i.e. at the level of corelation functions we only need to replace by as the quadrature is measured at one additional site. A similar argument for the expectation values of , and leads to

 Bm= 12∣∣m−1∑k=0αk[E(k+1,m−1−k)+E(k,m−k) + E(m−k,k)−E(m−k−1,k+1)]∣∣ (32)

To maximize this expression, we note that for two and three parties the maximum violation is reached for a state with . We thus make the reasonable hypothesis that it remains true for an arbitrary . Recall that we know how to choose and such as to reach these values. When , we have

 E(k,m−k) =cos[θ+(m−k)π2] (33)

and Eq. (III) becomes

 ⟨Bm−1⟩ =m−1∑k=0αkcos[θ+(m−1−k)π2] ⟨B′m−1⟩ =m−1∑k=0αkcos[θ+kπ2] (34)

Introducing Equation (33) into Eq. (III) and using some well known trigonometric formulae, the Bell factor simplifies to

 Bm= ∣∣∣cos(θ+mπ4)+sin(θ+mπ4)∣∣∣ ×∣∣ ∣∣m−1∑k=0αkcos[(m−2k)π4]∣∣ ∣∣ (35)

Maximizing the violation of local realism boils down to finding the optimal phase such that the first factor of the right hand side is maximum. This term achieves its maximum of for a value of the phase

 θm=(1−m)π4 (36)

We also note that

 (m−2k)π4=θm−1+(m−1−k)π2 (37)

hence the maximum value of the Bell factor can be finally written as

 Bmaxm =√2∣∣ ∣∣m−1∑k=0αkcos[θm−1+(m−1−k)π2]∣∣ ∣∣ (38) =√2Bmaxm−1 (39)

where we have identified the summation of Eq. (38) with Eq. (III) at the optimal angle . Introducing now the maximum value obtained in Ref. Wenger () for the two party case, , we obtain by recursion

 Bmaxm=2(m+1)/2 (40)

which is the known maximal bound imposed by quantum mechanics. Remarkably, the state defined in (23) combined with homodyne detection and a binning strategy called root binning allows for a maximal violation of the MK inequality. This result shows that even if the binning process discretizing the result of the homodyne detection discards some information, it does not prevent to maximally violate tests of local realism based on discrete variables.

## Iv Effect of noise

In Sections II and III, we proved that the search of loophole-free Bell tests might benefit from an increased number of parties involved in the experiment. The signature of this improvement lies in the exponential increase of the Bell factor with the number of parties . However, what makes a Bell test challenging in practice is not the magnitude of the violation, but rather the inevitable noise associated with any real experiment. In many cases, these imperfections are sufficient to hide the non-local correlations that one tries to observe. When the number of parties involved in a Bell test increases, so does the fragility of the state used in the experiment. The risk is thus to rescale the violation in such a way that no benefit of a larger is witnessed in practice. One can therefore correctly argue that an increased violation of local realism is only significative if accompanied by a comparable improvement of the robustness to noise of the test. After all, Bell tests have to be verified in a lab, not on paper.

In a discrete variable setting, the question of the tolerance to noise of a Bell test is often investigated introducing the noise fraction Krasz (). The noise fraction quantifies the maximum amount of depolarizing noise one can add to an entangled state and still detect non-local correlations. The depolarizing noise is characterized by the state , where is the dimension of the Hilbert space. However, in the continuous variable regime, the noise model underlying the noise fraction is irrelevant. Even if we deal with a finite number of photons, such as with the truncated photon number correlated states of Section II, the Hilbert space spanned by the eigenstates of the quadrature operators that are measured remains infinite dimensional. Hence, the appropriate basis is the infinite photon number bases and operators proportional to the identity have no physical meaning. To adopt an objective measure of the magnitude of the violation of local realism, we must therefore introduce a relevant noise model. In this Section, we will consider a probabilistic erasure: with probability , the mode of a random party is erased; otherwise it is untouched. This noise acts independently on each mode and transforms an initial state into

 ρ =(1−p)m|Ψ⟩⟨Ψ| +p(1−p)m−1{m∑t=1Trt(|Ψ⟩⟨Ψ|)⊗|0⟩t⟨0|} (41) +...+pm|0⟩1⟨0|⊗...⊗|0⟩m⟨0| (42)

Such a probabilistic erasure is known to appear in, e.g. atmospheric transmissions, and has recently been studied in Refs. wittmann (); niset ().

Let us first consider the photon number correlated states of Section II, and concentrate on the second term of Eq. (IV). Each element of the sum corresponds to the erasure of one of the subsystems. So, suppose for example that the state of subsystem has been erased and replaced by vacuum while distributing . The corresponding state shared between the parties reads

 ρ1,m =Trm(|Ψ⟩⟨Ψ|)⊗|0⟩m⟨0| (43) =∞∑r=0|cr|2|r⟩1⟨r|⊗...⊗|r⟩m−1⟨r|⊗|0⟩m⟨0|

This state is diagonal in the photon number bases, hence the results of all possible measurements are equiprobable, i.e., , , and all correlation coefficients vanish. This was to be expected from photon correlated states as their entanglement is truly -partite; tracing out one subsystem makes the state become separable. Thus, this property also holds for the other noisy terms of (IV), so that only the erasure-free term will contribute to the Bell factor. We obtain

 Bρ=(1−p)mBm (44)

To illustrate this result, consider the -partite GHZ state, Eq. (19). The noisy Bell factor reads , hence the maximum probability of erasure such that non-local correlations can be detected is

 pmax =1−√π221/2m (45)

This value exponentially tends towards as goes to infinity, and we observe the desired increased robustness to noise as becomes large.

Finally, consider now the states of Section III. First note that , hence in the bases, these states look like GHZ states. We thus expect their robustness to noise to behave like the photon correlated number states. Indeed, going through the calculation, one finds that for every noisy term of (IV), and for every choice of measurement setting the probability is an even function of the results, e.g. if party looses its mode, the probability to obtain given that is measured is an even function. As a result, none of the noisy terms contribute to the correlation coefficients, and the noisy Bell factor is again given by (44) as expected. As a conclusion, with respect to probabilistic erasure, more parties means more robustness.

## V A candidate 3-mode state

Let us now see how the results of Sec. III can be exploited in order to prepare multipartite states such that, on the one hand, they exhibit a significantly high violation of local realism and, on the other hand, they may be generated with near future technology. In particular, we will focus on a class of three-party states whose generation involves four optical “Schrodinger-cat” states, i.e., four single mode superpositions of coherent states NGaus_exp1 ().

Before proceeding let us make a remark. The main goal, as said, is now to exploit the gain expected in the amount of violation of locality when the number of involved parties is increased. To this aim, a natural candidate may have been the generalization of the photon-subtracted state discussed in Refs. Patron (); Nha (). Unfortunately, in the case of three parties (and, actually, in the case of any odd number of parties) the generalization of the strategy adopted in Refs. Patron (); Nha () is not effective. This is due to symmetry reasons and the use of the sign binning. Specifically, on the one hand, the joint probability distribution of any number of quadratures is an even function under the exchange of the arguments, on the other hand, the sign binning introduces an odd function (for odd number of parties) in the integration that has to be performed in order to obtain any correlation function. As a consequence, all the correlators are zero giving trivially no violation.

Let us now see how the increase in the number of parties involved in a Bell test can be exploited using a different approach, based on the results given in Sec. III. In Ref. Wenger (), the authors propose to use a superposition of Gaussians to implement the functions and of Eq. (23). In particular, in the case of only two Gaussians they considered the family of states defined by

 f(x) =⟨x|[c+(|α⟩+|−α⟩)], (46a) g(x) =⟨x|[c−(|α⟩−|−α⟩)], (46b)

where . One can then calculate the corresponding and coefficients using Eq. (III), which for large amplitudes gives and . As noticed in Ref. Wenger (), in the case of two parties no violation is possible, i.e . However, as can be seen in Eq. (III), one can achieve a violation of already for three parties. The corresponding state reads

 |Ψ3⟩=1√2[c3−(|α⟩−|−α⟩)⊗3+c3+(|α⟩+|−α⟩)⊗3]. (47)

where we have put .

Since the maximum violation is achieved for large amplitudes, we can consider the following simpler state:

 ∣∣Ψ′3⟩=c′(|α,α,α⟩+|α,−α,−α⟩+|−α,α,−α⟩+|−α,−α,α⟩), (48)

which coincides with for [where we defined ]. In order to obtain the Bell factor corresponding to such a state, we calculated the probabilities of the binned outcomes using Eq. (46) to define the roots, in combination with Eqs. (III) and (III). Specifically, we considered the domains and inherited by the state :

 D+x ={x∈R|x≥0} (49a) D−x ={x∈R|x<0} (49b) D+p ={p∈R|−cos(pα)sin(pα)≥0} (49c) D−p ={p∈R|−cos(pα)sin(pα)<0}. (49d)

The Bell coefficient calculated with such a procedure is shown in Fig. 3. One can see that for amplitudes as small as , the state already gives values above the local bound. We note that in this regime of small amplitudes, and the domains defined in Eqs. (49) might be non-optimal. As is increased, a violation around of the MK inequality is rapidly achieved.

Now let us describe how the state may be conditionally generated by using linear optics and superpositions of coherent states (SCS) of this form:

 (50)

Consider the scheme depicted in Fig. 4.

Two copies of the state in mode and are mixed in a balanced beam splitter. The same action is performed on modes and . Then modes and , as well as and are respectively mixed by means of two other beam splitters. As a last step mode is measured via a homodyne detector. It is straightforward to show that, when the measurement outcome is obtained, then the conditional state of the remaining modes coincides approximately with .

The main issue in implementing the scheme above concerns the generation of the four states , whose preparation is experimentally demanding. However, their generation in traveling light modes has been recently reported by several groups NGaus_exp1 (). Thus, one may envisage that, in the future, the implementation of the whole scheme of Fig. 4 will be possible.

## Vi Conclusions

In this paper, we have investigated the violation of local realism in Bell tests based on homodyne measurements performed on multipartite continuous-variable states. We have proven that the Mermin-Klyshko inequality supplemented with sign binning can be violated by an amount which grows exponentially with the number of parties. Furthermore, we have shown that it is possible to attain the maximal violation allowed by quantum mechanics by tuning the state and the binning strategy appropriately. The benefit of this multipartite approach was then shown to be effective in practice by analyzing the increased robustness to noise of the Bell tests with the number of modes.

It is worthwhile noting that our results are not based on a direct mapping from continuous to discrete variables. Such a procedure would in fact force a particular measurement to be used in the locality test, as the analysis of Ref. pseudospin () shows. We instead consider homodyne measurement because of its high detection efficiency. In this case, the possibility to obtain a maximal violation of the Mermin-Klyshko inequality is non-trivial, since (i) homodyne measurement represents a small subset of all possible measurements and (ii) the binning procedure may cause an irreversible loss of information. As it turns out, however, such a loss is not detrimental if suitable binning procedures are used, properly adapted to the states under investigation.

The experimental significance of our findings lies in that one can effectively benefit from the increased violation of locality in the considered multipartite continuous-variable scenario. In general, this allows for a greater freedom in the search for feasible multipartite states that can be used to test nonlocality in an actual experiment suffering from noise. As an illustration, inspired by the analysis of the states that maximally violate the local bound, we have proposed an approximation of the latter state which may be experimentally realized in a reasonable future. This state involves three parties and gives a violation of around of the local bound. As a perspective of this work, it is expected that other states of light with similarly high violation and suitably adapted to experimental constraints might be found, inspired by our results.

## Vii Acknowledgments

We gratefully acknowledge the financial support from the EU Projects QAP and COMPAS. A.A. and A.F. acknowledge the Spanish MEC projects FIS2007-60182 and Consolider Ingenio 2010 QOIT, and the Juan de la Cierva grant. N.J.C. and J.N. acknowledge financial support from the IUAP programme of the Belgian government under project Photonics@be. Finally, J.N. acknowledges support from the Belgian FRIA foundation.

## References

• (1) A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 777 (1935).
• (2) J.S. Bell, Physics (Long Island City, N.Y.) 1, 195 (1964).
• (3) A. Aspect, P. Grangier, and G. Roger, Phys. Rev. Lett. 47, 460 (1981).
• (4) W.J. Munro, Phys.  Rev. A 59, 4197 (1999).
• (5) J. Wenger et al., Phys. Rev. A 67, 012105 (2003).
• (6) R. García-Patrón Sánchez et al., Phys. Rev. Lett. 93, 130409 (2004); R. García-Patrón, J. Fiurasek and N. J. Cerf, Phys. Rev. A 71, 022105 (2005).
• (7) H. Nha, and H. J. Carmichael, Phys. Rev. Lett. 93, 020401 (2004).
• (8) A. Ourjoumtsev et al., Science 312, 83 (2006); J.S. Neergaard-Nielsen et al., Phys. Rev. Lett. 97, 083604 (2006); K. Wakui et al., Opt. Express 15, 3568 (2007); A. Ourjoumtsev et al., Nature 448, 784 (2007).
• (9) V. Parigi et al., Science 317, 1890 (2006); A. Zavatta, S. Viciani, M. Bellini, Science 306, 660 (2004).
• (10) J.F. Clauser, M.A. Horne, A. Shimony and R.A. Holt, Phys. Rev. Lett. 23, 880 (1969).
• (11) N.D. Mermin, Phys. Rev. Lett. 65, 1838 (1990); D.N. Klyshko, Phys. Lett. A 172, 399 (1993).
• (12) N. Gisin and H. Bechmann-Pasquinucci, Phys. Lett. A 246, 1 (1998).
• (13) P. van Loock and S.L. Braunstein, Phys. Rev. A 63, 0222106 (2001).
• (14) A. Ferraro and M.G.A. Paris, J. Opt. B 7, 174 (2005).
• (15) Z.B. Chen, J.W. Pan, G. Hou, and Y.D. Zhang, Phys. Rev. Lett. 88, 040406 (2002).
• (16) E.G. Cavalcanti et al., Phys. Rev. Lett. 99, 210405 (2007).
• (17) D. Kraszlikowski, P. Gnacinski, M. Zukowski, W. Miklaszewski and A. Zeilinger, Phys. Rev. Lett. 85, 4418 (2000).
• (18) C. Wittmann, D. Elser, U.L. Andersen, R. Filip, P. Marek and G. Leuchs, quant-ph/07041918 (2007).
• (19) J. Niset, U.L. Andersen and N.J. Cerf, arxiv:0710.4858 (2007).
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters