Entropic entanglement criteria in phase space Stefan Floerchinger,1, ∗ Martin Gärttner,1, 2, 3, y Tobias Haas,1, z and Oliver R. Stockdale3, x 1Institut für Theoretische Physik, Universität Heidelberg, Philosophenweg 16, 69120 Heidelberg, Germany 2Physikalisches Institut, Universität Heidelberg, Im Neuenheimer Feld 226, 69120 Heidelberg, Germany 3Kirchhoff-Institut für Physik, Universität Heidelberg, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany We derive inseparability criteria for the phase space representation of quantum states in terms of variants of Wehrl’s entropy. In contrast to entropic criteria involving differential entropies of marginal phase space distributions, our criteria are based on the Husimi Q-distribution. This is experimentally accessible through the heterodyne detection scheme, avoiding costly tomographic measurements. We apply our entropic criteria to Gaussian states and show that they imply a pair of second-order criteria for moments. We exemplify the strengths of our entropic approach by considering several classes of non-Gaussian states where second-order criteria fail. We show that our criteria certify entanglement in previously undetectable regions highlighting the strength of using the Husimi Q-distribution for entanglement detection. Introduction—Entanglement is the distinguishing fea- ferent approach and characterize the inseparability of a ture of quantum systems and will be a key resource in given bipartite state not by its marginal distributions, but future quantum technologies [1]. Fundamental problems by its joint probability distribution. This distribution, such as the thermalization of isolated quantum systems known as the Husimi Q-distribution, is a quasiprobability and characterizing quantum phase transitions, strongly distribution that contains the full information about the rely on a deep understanding of how entanglement man- state [33–37]. As it is non-negative and normalized, it has ifests and evolves [2, 3]. However a central problem in an associated entropy known as the Wehrl entropy [38, 39] studying entanglement is being able to derive experimen- and is bound from below per the Wehrl-Lieb inequal- tally accessible witnesses to detect it [4, 5]. The complex- ity [40, 41]. Crucially, the Husimi Q-function (and hence ity of such a problem depends upon the Hilbert space size, the Wehrl entropy) is an experimentally accessible quan- such that continuous variable systems with an infinite tity that can be measured via a heterodyne detection dimensional Hilbert space pose a particular challenge. protocol [12]. This is well-established within quantum For continuous quantum variables, many entanglement optics and has been realized in quantum optics experi- criteria rely on measuring the second-order moments of ments [42–44]. Furthermore, it was shown recently that two marginal distributions [6–11]. These criteria are most the Husimi Q-distribution can be measured in ultracold powerful whenever the state is Gaussian but are often atomic gases [45, 46]. insensitive elsewhere [12, 13]. This can pose significant The use of the the Wehrl entropy in the context of problems as there are many important classes of highly entanglement theory has only been discussed to define an entangled non-Gaussian states that cannot be witnessed entanglement monotone [47]. Here, we make use of the by second-order criteria [14–16]. Wehrl-Lieb inequality to derive entanglement criteria in To capture higher-order moments of measured distri- terms of entropies of the Husimi distribution. We show butions [17, 18], one can use differential entropies of that these criteria are stronger than previously known measured marginal distributions [19, 20]. Differential ones for certain classes of states within the non-Gaussian entropies reach beyond the scope of second-order criteria regime. arXiv:2106.08788v1 [quant-ph] 16 Jun 2021 since they are a functional of the full probability density Notation—We set ~ = 1 and disregard operator hats. function. Examples include criteria that rely on entropic We use capital letters for quantum operators Xj and uncertainty relations [21–25] as well as the complexity small letters for their corresponding eigenvalues xj and based criterion [26]. Other approaches are predicated on eigenvectors jxji. entropic uncertainty relations with (quantum) memory Husimi distribution—We consider a set of continuous formulated in terms of conditional differential entropies quantum variables Xj and Pj that fulfil the bosonic com- [27, 28] or as entropic steering inequalities [29, 30], see mutation relation [Xj;Pk] = iδjk where the subindices also [31]. denote the two subsystems j; k 2 f1; 2g. One disadvantage to these methods is that measuring A canonical transformation allows to do rotations in marginals of a distribution is often costly and impractical the local phase spaces by angles #j, as it requires angle tomography. This is particularly diffi- cult in ultracold quantum gas experiments where statistics R cos # sin # X j = j j j : (1) are limited [32]. In this work, we take a conceptually dif- Sj − sin #j cos #j Pj 2 To a set of (possibly rotated) position and momentum [50]. One can associate to them an entropy SM(Q±) that operators Rj andpSj, we define annihilation operators is defined in analogy to Wehrl’s entropy in Eq. (3) even Aj = (Rj + iSj)= 2 such that coherent states are their though s∓ is not the conjugate momentum of r±. eigenstates Aj jαi = αj jαi. The complexp eigenvalues αj Inseparability criteria—We first derive criteria for pure are parameterized as αj = (rj + isj)= 2. One can as- states and show they generalize to mixed states. We sociate a positive operator-valued measure (POVM) to consider pure separable states, for which the density op- pure coherent state projectors of the form Eα = jαi hαj. erator is a product ρ = ρ1 ⊗ ρ2. Here, the global Husimi This POVM is experimentally accessible through the het- Q-distribution factorizes erodyne detection protocol commonly used in quantum optical systems [12, 34, 35], which corresponds to mea- Q(r1; s1; r2; s2) = Q1(r1; s1) Q2(r2; s2); (7) suring the state ρ in the pure coherent state basis. This where Q (r ; s ) denotes the marginals of the global allows for full access to the phase space via a single mea- j j j Husimi Q-distribution. surement setting. Inserting Eq. (7) in Eq. (6) yields The distribution associated to the measurement out- comes for a set of conjugate quantum variables is the (±) Q±(r±; s∓) = (Q1 ∗ Q )(r±; s∓); (8) global Husimi Q-distribution 2 where ∗ denotes a convolution and Q(±) ≡ Q (±r; ∓s). In- Q(r ; s ; r ; s ) = TrfρE g = hαjρjαi : (2) 2 2 1 1 2 2 α voking the two-dimensional entropy power inequality [51– Due to the non-orthogonality of the coherent states, the 53] Husimi Q-distribution is a quasi-probability distribution eS(QA∗QB ) ≥ eS(QA) + eS(QB ); (9) in phase space. The Husimi Q-distribution is bounded, 0 ≤ Q ≤ 1, and normalized to unity with respect to the for any two two-dimensional Husimi Q-distributions Q Q A phase space measure j drjdsj=(2π). and QB, as well as the invariance of the Wehrl entropy Hence, an entropy associated to the Husimi Q- under mirror reflections in phase space allows to write distribution exists, known as the Wehrl entropy [38, 39] eSM(Q±) ≥ eSW(Q1) + eSW(Q2): (10) Z Y drj dsj S (Q) = − Q ln Q: (3) W 2π Thus, we find the pair of inequalities j SW(Q1) SW(Q2) The Wehrl entropy is strictly monotonous under partial SM(Q±) ≥ ln e + e ; (11) trace and can therefore be regarded as a classical (differ- which provide a state-dependent lower bound on the en- ential) entropy [48]. tropies S (Q ) obeyed by all pure product states. Hence, The Wehrl-Lieb inequality bounds the Wehrl entropy M ± pure states for which S (Q ) violates this bound are nec- from below [40, 41] (for 2N dimensional phase space) M ± essarily entangled. We call Eq. (11) the strong criteria. To obtain a state-independent bound, we apply the SW(Q) ≥ N; (4) Wehrl-Lieb inequality Eq. (4) to both subsystems, leaving and can therefore be understood as an entropic uncer- us with the criteria tainty relation in the quantum mechanical phase space. In contrast to many common uncertainty relations, it is SM(Q±) ≥ 1 + ln 2: (12) not invariant under one-mode squeezing but is instead As the latter relations are in general less tight than preserved under rotations in phase space. Eq. (11), we call them the weak criteria. Mixed marginalized distributions—It is useful to make We can generalize the weak criteria to mixed separable a linear transformation to non-local EPR-type vari- P i i P states ρ = i pi ρ1 ⊗ ρ2 , where pi ≥ 0 and i pi = 1. ables [49], On the level of the global Husimi Q-distributions, one has an analogous decomposition, leading, via Eq. (6), to r± = r1 ± r2; s± = s1 ± s2: (5) X i By integrating over half of these variables one obtains Q±(r±; s∓) = pi Q±(r±; s∓): (13) from (2) marginalized distributions i Using concavity of the entropy S (Q ) [54], we find Z dr ds M ± Q (r ; s ) = Q(r; s; ∓r ± r ; ±s ∓ s ): (6) ± ± ∓ 2π ± ∓ X i SM(Q±) ≥ piSM(Q±) ≥ 1 + ln 2; (14) These are normalized to unity with respect to the mea- i sure dr±ds∓=(2π) but, because of the twisted assignment where we have employed the strong pure state criteria (r±; s∓), they are not themselves Husimi distributions Eq. (11) and then the Wehrl-Lieb inequality Eq. (4) in 3 the last step. Therefore, the weak Wehrl entropic criteria To validate these second-order criteria, we consider how for pure product states Eq. (12) generalize identically to symplectic transformations in the ±-variables affect the mixed separable states [55]. The violation of inequality Husimi Q-distribution.