Nearly maximal violation of the Mermin-Klyshko inequality with multimode entangled coherent states

Entangled coherent states for multiple bosonic modes, also referred to as multimode cat states, not only are of fundamental interest, but also have practical applications. The nonclassical correlation among these modes is well characterized by the violation of the Mermin-Klyshko inequality. We here study Mermin-Klyshko inequality violations for such multi-mode entangled states with rotated quantum-number parity operators. Our results show that the Mermin-Klyshko signal obtained with these operators can approach the maximal value even when the average quantum number in each mode is only 1, and the inequality violation exponentially increases with the number of entangled modes. The correlations among the rotated parities of the entangled bosonic modes are in distinct contrast with those among the displaced parities, with which a nearly maximal Mermin-Klyshko inequality violation requires the size of the cat state to be increased by about 15 times.


I. INTRODUCTION
Quantum mechanics predicts many counterintuitive phenomena in microscopic world, among which nonlocality is particularly appealing. According to quantum mechanics, when the wavefunctions of two objects are entangled, there exist nonlocal correlations between the features of these two objects, irrespective of their spatial distance. Einstein, Podolsky, and Rosen (EPR) rejected this point of view, and thought that the results of measurements performed on two particles should be mutually independent when their separation is sufficiently large [1]. Based on this locality condition, they concluded that quantum mechanical description of nature is incomplete. Despite the debate between Einstein and Bohr [2], it was not until Bell's discovery in 1964 that it was realized that predictions of quantum mechanics conflict with local realism for some entangled states [3]. In the framework of local realism, Bell developed an inequality, which sets an upper bound on the statistical correlations of the results of measurements performed on two distant spin-1/2 particles. According to quantum mechanics, this bound is exceeded with proper settings of the measurement orientations when these two particles are prepared in a highly entangled state before the measurements.
Various forms of Bell inequality have been formulated, among which Clauser, Horne, Shimony and Holt (CHSH) version [4] is most famous and has been widely used for experimental test of quantum nonlocality [5][6][7]. The Bell-CHSH inequality has been generalized in different directions. On one hand, an inequality, referred to as the Mermin-Klyshko (MK) inequality, was derived for * Electronic address: zhirz@fzu.edu.cn † Electronic address: t96034@fzu.edu.cn systems composed of multiple spin-1/2 particles [8,9], which shows that the conflict between local hidden theories and quantum mechanics becomes stronger as the number of entangled particles increases. Experimental violation of the MK inequality was reported with fourphoton entanglement [10], and was analyzed in the context of ion traps [11]. On the other hand, Banaszek and Wödkiewicz constructed a Bell-CHSH-type operator for two harmonic oscillators with continuous variables, and formulated an inequality in phase space [12,13], which was later generalized to maximize the Bell signal [14,15]; however, with this formalism the inequality is not maximally violated by the original EPR state. In contrast, based on the so-called pseudospin operators introduced by Chen et al. [16], the original EPR state can maximally violate the inequality, but the experimental measurements remains challenging. We have constructed a new Bell operator [17], with which the Bell signal for a two-mode cat state of moderate size can approach the upper bound of quantum mechanics. The issue of quantum nonlocality associated with multimode continuous variable states has also been addressed [18][19][20][21][22][23]. In particular, Jeong et al. [18] have discussed the quantum nonlocality associated with entangled coherent states for three bosonic modes using joint quasiprobability distribution function in phase space. With this method, there are 12 variables to be optimized to find the global maximum value of the MK signal, which is a difficult task. In Ref. [18], only some local maximum values were numerically found. For the Greenberger-Horne-Zeilinger-type (GHZ-type) [24] entangled coherent states, also referred to as 3-mode Schrödinger cat states [25][26][27][28], the obtained MK signal increases with the amplitude α of the coherent state components, reaching a maximal value of 3.6 when α → ∞; this maximum is still lower than the quantummechanical upper bound of 4 for a three-partite system by 10%. With this method, the number of the variables to be optimized increases with the number of entangled bosonic modes, which makes it extremely difficult to obtain the optimal MK signal when more bosonic modes are involved. We note that, multimode entangled coherent states are intriguing both from the fundamental view point and for practical applications because the components forming these highly nonclassical states have classical analogs and can be used for fault-tolerant quantum computation [29].
In this manuscript, we analyzed the MK inequality violations for GHZ-type entangled coherent states of n bosonic modes based on effective rotated parity correlations proposed in Ref. [17]. Our results show that the MK signal quickly increases with α, and approaches the quantum-mechanical upper bound when α is only about 1. This is in distinct contrast with the result obtained in Ref. [18], which shows the MK signal for the 3-mode cat state at α = 1 is even significantly below the lower bound for confirming true 3-partite entanglement. We further show that the MK signal obtained with the joint quasiprobability distribution function can also approach the quantum-mechanical upper bound for a sufficiently large cat size. However, with this approach the size of the cat state needs to be increased by about 15 times to obtain a nearly maximal MK inequality violation.
The paper is organized as follows. In Sec. 2, we construct the MK operator for n-mode cat states with effective rotated parity operators. With this formalism, the correlations among the bosonic modes are in almost perfect analogy to those among n qubits prepared in a GHZ state even for α ∼ 1, in contrast with the results obtained in continuous-variable representation. In Sec. 3, We detailedly analyze the MK inequality violation for the 3-mode cat state based on these correlations. Analytical results qualitatively show that the size of the cat state does not need to be large for the MK signal to approach the quantum-mechanical upper bound. We confirm this prediction with numerical simulation, which demonstrates that a nearly maximal inequality violation can be obtained even when α is only about 1. We find that, when using joint quasiprobability distribution function, the average quantum number in each mode needs to be increased by about 15 times to obtain the same degree of violation. In Sec. 4, we investigate the MK signals for the 4-mode and 5-mode cat states. Numerical simulations show that, for these cases the MK signals can also nearly approach the corresponding quantum-mechanical upper bounds even when the average quantum number in each mode is only 1. Conclusions are presented in Sec. 5.

II. CONSTRUCTION OF THE MK OPERATOR FOR N-MODE CAT STATES
The n-mode Schrödinger cat state under consideration is defined as the equal superposition of all the n bosonic modes being in the coherent state |α and all in |−α , where N n = 2 + 2e −2n|α| 2 −1/2 and the subscripts j = 1 to n label these modes. When |α| 2 ≫ 1, the two coherent states |α and |−α are approximately orthogonal and can be considered two logic states of a qubit. With this encoding, the n-mode Schrödinger cat state is in analogy with the GHZ state of n qubits (2) According to quantum mechanics, when n qubits are prepared in a GHZ state, there are nonlocal correlations among the outcomes of measurements individually performed on them. The conflict between the predictions of quantum mechanics and the results allowed by local realism becomes stronger as the number of entangled qubits increases, as evidenced by the violations of the MK inequality. The MK operator for a multipartite entangled state is recursively defined as Here σ k,a k and σ k,a ′ k represent two-valued observables for the kth party, with a k and a ′ k denoting the corresponding measurement settings, and O ′ k is obtained from O k by exchanging all a j and a ′ j (j ≤ k). For a multi-qubit system, a k is a unit vector and σ k,a k corresponds to the Pauli operator along the direction a k . In the framework of local hidden variable (LHV) theories, one can assign a value of 1 or −1 to σ k,a k and σ k, As a result, the outcome of each measurement on O n is 2 (n−1)/2 or −2 (n−1)/2 , which implies that the MK signal within a LHV model, defined as S l n = | O n | LHV , satisfies S l n ≤ D n = 2 (n−1)/2 . To show the violation of the MK inequality by the n-qubit GHZ state, we set all of a k and a ′ k are on the xy-plane, and for simplicity, use the notation σ k (φ k ) to replace σ k,a k , with φ k being the angle between a k and the x axis. Set Then we have where H.c. denotes the Hermitian conjugate. According to quantum mechanics, the expectation value of O n for the ideal n-qubit GHZ state is S n = O n QM = 2 n−1 , which exceeds the bound imposed by the local realism by a factor of 2 (n−1)/2 , indicating the violation of the MK inequality exponentially grows with the number of entangled qubits. We first analyze the MK inequality violation for the n-mode cat state based on the effective cat state qubit rotation operators. For the kth cat state qubit, the quantum-number parity operator P k = (−1) a † k a k acts as σ k,x , where a † k and a k is the quantum-number rising and lowering operators. We note the other components of the Pauli operators in the xy-plane for the cat state qubits can be operationally constructed by combining the parity operator and the effective rotation operator, defined as [17] where is the phase gate, and is effectively equivalent to a rotation of the kth cat state qubit around the z axis of the Bloch sphere by an angle φ k . The observables needed for test of the MK inequality can be expressed in terms of this kind of rotation operators and the parity operator The effective rotated parity operator σ k (φ k ) corresponds to the Pauli operator along the axis with an angle φ k to the x axis of the cat state qubit. For a finite value of α, the coherent states |α j and |−α j are not strictly orthogonal, so that the n-mode cat state is not perfectly equivalent to the n-qubit GHZ state. As a consequence, the MK signal S n is smaller than 2 (n−1)/2 . However, the overlap between |α j and |−α j , e −4|α| 2 , is quite small even for a moderate value of α. For example, for α = √ 2 this overlap is only about 2.5 × 10 −4 . Therefore, the MK signal S n for the n-mode cat state of a moderate size, obtained with the effective cat state qubit rotation operators, can approximate the result for an n-qubit GHZ state.
The test of quantum nonlocality in phase space is based on the measurement of displaced quantum-number parity observables, defined as In this formalism, P k (β k ) and P k (β ′ k ) corresponds to σ k,a k and σ k,a ′ k , and the recursive definition of the MK operator can be written as The expectation value of each term of Q n is proportional to the joint Wigner function (quasiprobability distribution function) for these n bosonic modes at the corresponding point in phase space. For a limited value of α, the displacement operation D k (β k ) moves the corresponding bosonic mode out of the logic space of the corresponding cat state qubit, and is not equivalent to a qubit rotation. Therefore, the obtained MK signal is expected to be smaller than that based on the effective cat qubit rotation operators.

III. VIOLATIONS OF THE MK INEQUALITY WITH THE 3-MODE CAT STATE
To quantitatively show the violations of MK inequality based on the above-mentioned rotated parity correlations, we will derive analytic expressions for the quantum mechanical expectation values of the MK observable O n and then perform numerical simulations. We first consider the 3-mode cat state, for which the MK operator based on the effective rotated parity correlations among the three modes is given by The MK signal, defined as the absolute value of the expectation value of O 3 , is S 3 = |E 3 (0, −π/4, π/4) + E 3 (0, π/4, −π/4) + E 3 (π/2, −π/4, −π/4) − E 3 (π/2, π/4, π/4)|, (12) where with When e −2|α| 2 ≪ 1, the dominant term of . The ratio of this term to the corresponding correlation for a GHZ state is 2N 2 3 , which approximates 1 even for a very moderate value of α. For example, it is about 0.9975 for α = 1. This allows a nearly maximal violation of the MK inequality with the average quantum number in each mode being only about 1.
To confirm the validity of the rotated parity operators for revealing the quantum nonlocality among the three bosonic modes, we first perform numerical simulations of the four correlations E 3 (0, −π/4, π/4), E 3 (0, π/4, −π/4), E 3 (π/2, −π/4, −π/4), and E 3 (π/2, π/4, π/4) as functions of α. As shown in Eqs. (13) and (14), these correlations are independent of the argument of the complex amplitude α, which is taken to be a positive real number for simplicity. The results are respectively shown in Fig. 1(a), (b), (c), and (d), where the solid lines denote the correlations of the 3-mode cat states, while the dotted lines represent the corresponding correlations for the equally-weighted classical mixture of the two components |α 1 |α 2 |α 3 and |−α 1 |−α 2 |−α 3 , which are given by For α = 0, each mode is in the vacuum state |0 , which is not affected by the rotation operator R † k,z (φ k ). Consequently, the rotated parity operator σ k (φ k ) is equivalent to the parity operator P k , which has a definite value of 1, so that each correlation is 1. In this case, all correlations are completely classical as there is no entanglement. As α increases, for both |α and |−α the difference between the occupations between the even and odd quantumnumber states quickly decreases, and consequently the correlations E m 3 (φ 1 , φ 2 , φ 3 ) for the classical mixture drop fast, and approach 0 when e −2|α| 2 ≪ 1. In contrast, for this case the 3-mode cat state exhibits nearly perfect correlations, like the three-qubits GHZ state |ψ GHZ , which is the eigenstate of each of the four terms in O 3 , with the eigenvalues corresponding to the first three terms being 1 and that associated with the last term being −1. This is in distinct contrast with the two-mode case, where the correlations used to construct the Bell signal are not perfect even when α → ∞; the absolute value of each correlation tends to √ 2/2 in this limit [17]. These reuslts unambiguously demonstrate that, when e −2|α| 2 ≪ 1, the correlations associated with the 3-mode cat state almost completely arise from the quantum coherence between |α 1 |α 2 |α 3 and |−α 1 |−α 2 |−α 3 , which is responsible for the entanglement among the three modes.
The MK signal (S 3 ) for the 3-mode cat state as a function of α is shown in Fig. 2, with the blue solid line representing the simulated result while the black dotted line denoting the bound imposed by local hidden variable theories. The result shows that S 3 first drops below 2, and then quickly increases with α. This non-monotonous behavior is similar to the two-mode case [17]. When α < 0.266, E 3 (π/2, −π/4, −π/4) and E 3 (π/2, π/4, π/4) almost have the same value and thus cancel each other out, so that S 3 is approximately equal to the sum of the other two correlations, each of which slightly decreases as α increases in this regime. Consequently, S 3 drops as α increases in this regime. When α > 0.266, the difference between E 3 (π/2, −π/4, −π/4) and E 3 (π/2, π/4, π/4) increases quickly, whose contribution makes S 3 rise quickly. (a), E3(0, −π/4, π/4); (b), E3(0, π/4, −π/4); (c), E3(π/2, −π/4, −π/4); (d), E3(π/2, π/4, π/4). The solid and dotted lines corrrespond to the correlations for the cat state and the equally-weighted classical mixture of the two components |α 1 |α 2 |α 3 and |−α 1 |−α 2 |−α 3 , respectively. and exceeds the lower bound 2 √ 2 for true 3-partite entanglement [30] when α > 0.581. When α further increases, it quickly approaches the quantum-mechanical upper bound 4, which implies that a nearly maximal MK inequality violation can be obtained with a small cat state. For example, when α = 1, S 3 ≃ 3.916, very close to upper bound 4. This can be interpreted as follows. The two coherent states |α and |−α of each bosonic mode are approximately orthogonal to each other and can act as the basis states of a logic qubit even for a moderate value of α. The bosonic mode is well restricted in the logic space of the cat state qubit under the application of the effective cat state qubit rotation operator R k,z (φ k ). Consequently, the rotated parity operator σ k (φ k ) is approximate to the corresponding Pauli operator, and the 3-mode cat state is approximately an eigenstate of each of the four correlation operators in Eq. (11).
As noted in Ref. [18], for the MK signal based on the joint Wigner function, there are 6 complex displacement parameters, β k and β ′ k with k = 1 to 3, corresponding to 12 variables to be optimized to obtain the maximal MK signal for each value of α, which is a difficult task. Instead of the global maximum value, a local maximum value for each value of α was numerically obtained in Ref. [18], which monotonously increases with α, and approaches 3.6, instead of 4, when α → ∞. We note that the upper bound 4 can be reached in this limit with the choice β 1 = 0, β ′ 1 = iπ/8α, β k = −iπ/16α, β ′ k = iπ/16α, with k > 1. This is due to the fact that the displaced parity operators are equivalent to the corresponding Pauli operators of the cat state qubits when α → ∞, so that the MK inequality can be maximally violated in this limit. Based on these displaced parity correlations, the MK signal at According to the result of Ref. [18], the MK signal at α = 1 is only about 2.5, which is significantly below the lower bound for confirming true 3-partite entanglement. This is in stark contrast with the result obtained with our framework, where the upper bound allowed by 3-partite entanglement can be approached for α ≃ 1.

V. CONCLUSIONS
In conclusion, we have analyzed quantum nonlocality for n bosonic modes in an entangled coherent state using rotated parity operators. Our results show that the correlations obtained with this approach are close to those for the n-qubit GHZ state, and the MK inequality can be nearly maximally violated even when the average quan-tum number in each mode is only about 1. This is in distinct contrast with the previous study based on displaced parity operators, where the obtained MK signal for the 3-mode cat state even does not exeed the lower bound for genuine 3-partite entanglement for this size. With this framework, a nearly maximal inequality violation requires the quantum number in each mode to be improved by about 15 times compared to the present result. Since the decoherence rate of a cat state increases with the size, our results are useful for experimentally investigating quantum nonlocality for entanglement of quasiclassical states of multiple bosonic modes.