This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Negativity volume of the generalized Wigner function in gravitating hybrid system

Daisuke Miki,1 Akira Matsumura,1 Kazuhiro Yamamoto1,2,3 1Department of Physics, Kyushu University, 744 Motooka, Nishi-Ku, Fukuoka 819-0395, Japan 2Research Center for Advanced Particle Physics, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395, Japan [email protected],
[email protected],
[email protected]
3International Center for Quantum-field Measurement Systems for Studies of the Universe and Particles (QUP), KEK, Oho 1-1, Tsukuba, Ibaraki 305-0801, Japan
Abstract

We evaluate the gravity-induced negativity volume of the generalized Wigner function in a hybrid system consisting of a particle in a two-localized superposition state and an oscillator. The generalized Wigner function can capture the nonclassicality of the system. The increase in the negativity volume of the generalized Wigner function can be an indicator of the entanglement generation, which is demonstrated in the hybrid system generating the gravity-induced entanglement in various initial states. Moreover, by comparing the behaviors of the negativity volume with the entanglement fidelity, we show that the nonclassical feature of entanglement is properly identified by the criterion based on the negativity volume of the generalized Wigner function when the oscillator is initially in a thermal state.

I INTRODUCTION

Wigner function describes a distribution function in the phase space, which appears in various fields of sciences (e.g., [1]). It is known that the Wigner function defined for a particle in classical mechanics always takes positive values. In contrast to classical mechanics, the Wigner function defined for a continuous variable system in quantum mechanics may have negative values, which is used as an indicator of the nonclassicality of the system [2]. Such a property of the Wigner function can be measured in experiments, which is commonly used in quantum optics as a quantum signature [3, 4, 5, 6].

The Wigner function for qubit systems has been extensively proposed [7, 8, 9, 10, 11, 12], which is defined on the Bloch sphere using the Euler angles. This Wigner function enables us to analyze the discrete systems in the continuous phase space parameterized by the angles. However, unlike the usual Wigner function for continuous variable systems, the Wigner function for qubit systems has negative regions on the Bloch sphere even for the arbitrary pure state of qubits. Hence, the negative region of the Wigner function does not simply characterize the nonclassicality of the discrete system. Thus, the Wigner function as an indicator of quantumness including qubit systems is not trivial. Under such circumstances, Ref. [10] introduced the negativity volume of the generalized Wigner function as an indicator of the nonclassical feature for a hybrid system.

The hybrid system, consisting of the discrete and continuous variable subsystems, is one of the central research in quantum information science; quantum information processing [13], quantum teleportation [14, 15], and quantum entanglement [16, 17]. The generalized Wigner function in the hybrid system is obtained in Ref. [7], but it is demonstrated that the negative values do not necessarily characterize the nonclassicality because of the property of the Wigner function in qubit systems. However, Ref. [10] showed that there is a critical value for the negativity volume of the Wigner function for the separable hybrid system; if the negativity volume is larger than the critical value, nonclassicality appears. We note that this criterion is a sufficient condition, not a necessary condition.

In the present study, we consider the generalized Wigner function of the hybrid system consisting of a particle in a two-state superposition and an oscillator, which are gravitationally interacting. Refs. [18, 20, 21, 19] discussed the quantum entanglement through the gravitational interaction in this hybrid system. Recently, the experimental proposal of the gravity-induced entanglement between two particles has attracted interest [22, 23]. The quantum entanglement is known as a nonlocal quantum correlation, which is not newly generated by classical processes (e.g. [24]). Refs. [22, 23] showed that two particles in a superposition state are entangled if Newtonian gravity obeys quantum mechanics, which was shown as a feasible experimental proposal in the future. This two-particle model has been the subject of several advanced studies [25, 26, 27, 28, 29]. In addition, gravitational entanglement between mechanical oscillators in optomechanical systems has also been investigated [30, 31, 32, 33, 34, 35]; experiments [36, 37, 38] and theories [39, 40, 41, 42, 43] have been developed towards the realization of macroscopic quantum states to verify the quantum nature of gravity. Furthermore, quantum entanglement generation by Newtonian gravity has been discussed from the viewpoint of quantum field theory and gravitons [44, 45, 46, 47, 48, 49, 50, 51, 52, 53].

For the hybrid system considered in the present paper, we evaluate the gravity-induced negativity volume of the generalized Wigner function as a quantity to characterize the quantumness of gravity and discuss its correspondence with the condition for the quantum entanglement based on the generalized Wigner function [10]. We discuss the applicability of the quantum entanglement criterion with the generalized Wigner function in comparison with the entanglement fidelity criterion proposed in Ref. [54, 55, 56]. In particular, we show that the condition based on the generalized Wigner function correctly captures the entanglement even when the oscillator is initially in a thermal state.

The present paper is organized as follows: In Sec. II, we briefly review the entanglement criterion with the generalized Wigner function for a hybrid system, which was proposed in Ref. [10]. In Sec. III, after briefly reviewing the hybrid system to detect the gravity-induced entanglement considered in Refs. [18, 19], we evaluate the evolution of the negativity volume in the hybrid system for three different initial states to determine the entanglement generation. In Sec. IV, we evaluate the entanglement fidelity for the three states introduced in Sec. III, to verify that the criterion with the negativity volume works properly. Sec. V is devoted to summary and conclusions. In Appendix A, we introduce the generalized Wigner function and explain the negativity volume for a hybrid system following Ref. [10].

II WIGNER FUNCTION FOR HYBRID SYSTEM

In this section, we introduce the generalized Wigner function for the hybrid system described by a continuous variable and a discrete variable. We first consider the Wigner function for the continuous variable system defined by

Wb\displaystyle W_{b} =Tr[ρbΔ^b],Δ^b=D^bΠ^bD^b,\displaystyle=\text{Tr}[\rho_{b}\hat{\Delta}_{b}],\quad\hat{\Delta}_{b}=\hat{D}_{b}\hat{\Pi}_{b}\hat{D}_{b}^{\dagger}, (1)

where ρb\rho_{b} is the density matrix, D^b=eβb^βb^\hat{D}_{b}=e^{\beta\hat{b}^{\dagger}-\beta^{*}\hat{b}} is the displacement operator, and Π^b=2πeiπb^b^\hat{\Pi}_{b}=\frac{2}{\pi}e^{i\pi\hat{b}^{\dagger}\hat{b}} is the parity operator, and b^\hat{b} and b^\hat{b}^{\dagger} are the creation and annihilation operators. It is well known that this Wigner function characterizes the nonclassical nature of continuous variable systems and that it has negative values in quantum states such as a Shrödinger-cat state.

On the other hand, the Wigner function for the discrete system is defined to satisfy the Stratonovich-Weyl conditions as follows [7, 8, 9, 10, 11, 12]:

Wq\displaystyle W_{q} =Tr[ρqΔ^q],Δ^q=U^qΠ^qU^q,\displaystyle=\text{Tr}[\rho_{q}\hat{\Delta}_{q}],\quad\hat{\Delta}_{q}=\hat{U}_{q}\hat{\Pi}_{q}\hat{U}_{q}^{\dagger}, (2)

where U^q=ei2ϕσ^zei2θσ^yei2Φσ^z\hat{U}_{q}=e^{\frac{i}{2}\phi\hat{\sigma}_{z}}e^{\frac{i}{2}\theta\hat{\sigma}_{y}}e^{\frac{i}{2}\Phi\hat{\sigma}_{z}} is the unitary operator of SU(2) and Π^q=12(𝟙^3σ^z)\hat{\Pi}_{q}=\frac{1}{2}(\hat{\mathbbm{1}}-\sqrt{3}\hat{\sigma}_{z}) is the parity operator. ϕ,Φ[0,2π],θ[0,π]\phi,\Phi\in[0,2\pi],\theta\in[0,\pi] are the Euler angles and σ^y,σ^z\hat{\sigma}_{y},~{}\hat{\sigma}_{z} are the Pauli operators. The negative value of the Wigner function of the discrete system does not simply characterize the quantum nature of the system. For the general pure qubit state,

ρqp\displaystyle\rho_{q}^{p} =|qq|,|q=α|0+eiχ1α|1,α[0,1],χ[0,2π],\displaystyle=\ket{q}\bra{q},\quad\ket{q}=\sqrt{\alpha}\ket{0}+e^{i\chi}\sqrt{1-\alpha}\ket{1},\quad\alpha\in[0,1],\quad\chi\in[0,2\pi], (3)

we derive the Wigner function as

Wq[ρqp]\displaystyle W_{q}[\rho_{q}^{p}] =3α(1α)sinθcos(χ+2ϕ)+3(12α)/2cosθ+1/2.\displaystyle=\sqrt{3\alpha(1-\alpha)}\sin\theta\cos(\chi+2\phi)+\sqrt{3}(1-2\alpha)/2\cos\theta+1/2. (4)

Considering the state |q=|0\ket{q}=\ket{0} for α=1\alpha=1, the Wigner function reduces to

Wq[|00|]\displaystyle W_{q}[\ket{0}\bra{0}] =12(13cosθ),\displaystyle=\frac{1}{2}\left(1-\sqrt{3}\cos\theta\right), (5)

which has the negative values on the Bloch sphere. Hence, the Wigner function of the discrete system has negative values even for a pure state.

To characterize the nonclassicality of the hybrid system, the authors of Ref. [10] introduced the negativity volume of the generalized Wigner function as

𝒱[ρ]\displaystyle\mathcal{V}[\rho] =12(𝑑ϕ𝑑θd2βsinθ2π|W|1),\displaystyle=\frac{1}{2}\left(\int d\phi d\theta d^{2}\beta\frac{\sin\theta}{2\pi}|W|-1\right), (6)

where WW is the generalized Wigner function defined by

W[ρ]=\displaystyle W[\rho]= Tr[ρΔ^qΔ^b].\displaystyle\text{Tr}[\rho\hat{\Delta}_{q}\otimes\hat{\Delta}_{b}]. (7)

The Wigner function and the negativity volume for the reduced density matrix for the continuous variable system and for the qubit are respectively defined as

Wb\displaystyle W_{b} =Tr[ρbΔb],𝒱b[ρb]=12(d2β|Wb|1),\displaystyle=\text{Tr}[\rho_{b}\Delta_{b}],\qquad\mathcal{V}_{b}[\rho_{b}]=\frac{1}{2}\left(\int d^{2}\beta|W_{b}|-1\right), (8)
Wq\displaystyle W_{q} =Tr[ρqΔq],𝒱q[ρq]=12(sinθ2π𝑑ϕ𝑑θ|Wq|1),\displaystyle=\text{Tr}[\rho_{q}\Delta_{q}],\qquad\mathcal{V}_{q}[\rho_{q}]=\frac{1}{2}\left(\int\frac{\sin\theta}{2\pi}d\phi d\theta|W_{q}|-1\right), (9)

where ρb=Trq[ρ]\rho_{b}=\text{Tr}_{q}[\rho] and ρq=Trb[ρ]\rho_{q}=\text{Tr}_{b}[\rho]. The authors of Ref. [10] showed that the negativity volume of the arbitrary pure qubit state (3) is the same value

𝒱q[ρqp]\displaystyle\mathcal{V}_{q}[\rho_{q}^{p}] =1312,\displaystyle=\frac{1}{\sqrt{3}}-\frac{1}{2}, (10)

and suggested that this negativity volume is the upper bound for any qubit state. They also expected that the entanglement criterion for the hybrid system is expressed as follows; the state is entangled if the inequality

𝒱[ρ]\displaystyle\mathcal{V}[\rho] 𝒱cr=23ipi𝒱b[ρbi]+1312,\displaystyle\leq\mathcal{V}_{\text{cr}}=\frac{2}{\sqrt{3}}\sum_{i}p_{i}\mathcal{V}_{b}[\rho_{b}^{i}]+\frac{1}{\sqrt{3}}-\frac{1}{2}, (11)

is violated. Here, the critical value 𝒱cr\mathcal{V}_{\text{cr}} corresponds to the maximum value for the separable hybrid system written as ρsep=ipiρqiρbi\rho^{\text{sep}}=\sum_{i}p_{i}\rho_{q}^{i}\otimes\rho_{b}^{i}. Hence, this criterion indicates that when the negativity volume is larger than the critical value, additional nonclassicality appears, namely entanglement.

On the other hand, there is no upper limit on the negativity volume for the reduced continuous variable system 𝒱b\mathcal{V}_{b}. There is the possibility in which a basis transformation can result in different values of the reduced negativity volume, even for the same state. For example, we consider the state |ϕ|0q|αb+|1q|αb\ket{\phi}\propto\ket{0}_{q}\ket{\alpha}_{b}+\ket{1}_{q}\ket{-\alpha}_{b}, where |0q\ket{0}_{q} and |1q\ket{1}_{q} are the states of a qubit, and |±αb\ket{\pm\alpha}_{b} is the coherent state of a continuous variable. Tracing over the qubit system, the reduced state is mixed state ρb|αα|bb+|αα|bb\rho_{b}\propto\ket{\alpha}{}{}_{b}{}_{b}{}\bra{\alpha}+\ket{-\alpha}{}{}_{b}{}_{b}{}\bra{-\alpha} and the reduced negativity volumes 𝒱b[|αα|bb]\mathcal{V}_{b}[\ket{\alpha}{}{}_{b}{}_{b}{}\bra{\alpha}] and 𝒱b[|αα|bb]\mathcal{V}_{b}[\ket{-\alpha}{}{}_{b}{}_{b}{}\bra{-\alpha}] are zero. However, we consider the transformation of the basis as |+q=|0q+|1q\ket{+}_{q}=\ket{0}_{q}+\ket{1}_{q} and |q=|0q|1q\ket{-}_{q}=\ket{0}_{q}-\ket{1}_{q}, and the reduced state is written in the mixed state ρb(|α+b|α)b(α|b+α|b)+(|αb|α)b(α|bα|b)\rho_{b}\propto(\ket{\alpha}{}{}_{b}+\ket{-\alpha}{}{}_{b})({}_{b}{}\bra{\alpha}+{}_{b}{}\bra{-\alpha})+(\ket{\alpha}{}{}_{b}-\ket{-\alpha}{}{}_{b})({}_{b}{}\bra{\alpha}-{}_{b}{}\bra{-\alpha}). Then, the reduced negativity volumes 𝒱b[(|α+b|α)b(α|b+α|b)]\mathcal{V}_{b}[(\ket{\alpha}{}{}_{b}+\ket{-\alpha}{}{}_{b})({}_{b}{}\bra{\alpha}+{}_{b}{}\bra{-\alpha})] and 𝒱b[(|αb|α)b(α|bα|b)]\mathcal{V}_{b}[(\ket{\alpha}{}{}_{b}-\ket{-\alpha}{}{}_{b})({}_{b}{}\bra{\alpha}-{}_{b}{}\bra{-\alpha})] are non-zero. Thus, the treatment of negativity volume for the reduced system of a continuous variable system coupled to a qubit is non-trivial. Despite such a problem, in the present paper, we demonstrate that the criterion based on the negativity volume of the generalized Wigner function is efficient by comparing the fidelity-based condition for various initial states such as superposed qubit-coherent, -thermal, and -cat state. We omit the subscript "b""b" and "q""q" for simplicity.

III negativity volume AND ENTANGLEMENT INDUCED BY GRAVITY IN HYBRID SYSTEM

In this section, we consider the hybrid system consisting of a spin particle with mass mm and an oscillator with mass MM considered in Refs. [18, 19]. Figure 1 shows a schematic plot of the system. The spin particle is prepared as the superposition state of the two spatially localized states separated by distance \ell by using the Stern-Gerlach scheme. Then, the oscillator also becomes the superposition state due to gravitational interaction with the spatially-localized superposed particle. The Hamiltonian of the hybrid system is

H^\displaystyle\hat{H} =ωa^a^+H^g,\displaystyle=\omega\hat{a}^{\dagger}\hat{a}+\hat{H}_{g}, (12)

where a^,a^\hat{a},~{}\hat{a}^{\dagger} are the creation and annihilation operators and ω\omega is the frequency of the oscillator. Here, we ignore the Larmor precession for simplicity. Hamiltonian H^g\hat{H}_{g} describes the gravitational interaction given by

Hg\displaystyle H_{g} =GMmL2+(q^+σ^z/2)2GMmL2+2/43σ^zq^+const.,\displaystyle=-\frac{GMm}{\sqrt{L^{2}+(\hat{q}^{\prime}+\hat{\sigma}_{z}\ell/2)^{2}}}\simeq\frac{GMm\ell}{\sqrt{L^{2}+\ell^{2}/4}^{3}}\hat{\sigma}_{z}\hat{q}^{\prime}+\text{const.}, (13)

where q^\hat{q} is the position operator of the oscillator, LL is the distance between the spin particle and the oscillator, \ell is the distance between the positions of the superposed particle and GG is the gravitational constant. Here, we define the dimensionless position operator q^=Mωq^\hat{q}=\sqrt{M\omega}\hat{q}^{\prime} and introduce the gravitational coupling gg as

g\displaystyle g =GMmlL2+2/4312Mω.\displaystyle=\frac{GMml}{\sqrt{L^{2}+\ell^{2}/4}^{3}}\frac{1}{\sqrt{2M\omega}}. (14)

Hamiltonian HgH_{g} is rewritten as H^g=2gσ^zq^\hat{H}_{g}=\sqrt{2}g\hat{\sigma}_{z}\hat{q}^{\prime} and the unitary evolution is

U^(t)\displaystyle\hat{U}(t) =eiH^t=eiωa^a^teσ^z(αta^αta^)+iCt2,\displaystyle=e^{-i\hat{H}t}=e^{-i\omega\hat{a}\hat{a}^{\dagger}t}e^{-\hat{\sigma}_{z}(\alpha^{*}_{t}\hat{a}^{\dagger}-\alpha_{t}\hat{a})+iC^{2}_{t}}, (15)

where we defined

αt\displaystyle\alpha_{t} =λ(eiωt1),Ct=λ2(ωtsin[ωt]),withλ=gω.\displaystyle=\lambda(e^{-i\omega t}-1),\quad C_{t}=\lambda^{2}(\omega t-\sin[\omega t]),\quad\text{with}\quad\lambda=\frac{g}{\omega}. (16)
Refer to caption
Figure 1: Schematic plot of the model as a hybrid system. This figure assumes γ=0\gamma=0 in Eq. (18).

III.1 INITIAL QUBIT-COHENRENT STATE

We assume that the particle and the oscillator are initially in a superposed state and a coherent state, respectively. The state of the total system is

|ϕ\displaystyle\ket{\phi} =12(|0+|1)|γ,\displaystyle=\frac{1}{\sqrt{2}}(\ket{0}+\ket{1})\otimes\ket{\gamma}, (17)

which evolves by the Hamiltonian (12) as

|ϕ(t)\displaystyle\ket{\phi(t)} =12(e(αtγαtγ)/2|0|αt+γt+e(αtγαtγ)/2|1|αt+γt),\displaystyle=\frac{1}{\sqrt{2}}\left(e^{-(\alpha^{*}_{t}\gamma^{*}-\alpha_{t}\gamma)/2}\ket{0}\ket{\alpha_{t}+\gamma_{t}}+e^{(\alpha^{*}_{t}\gamma^{*}-\alpha_{t}\gamma)/2}\ket{1}\ket{-\alpha_{t}+\gamma_{t}}\right), (18)

where γt=eiωtγ\gamma_{t}=e^{-i\omega t}\gamma and we ignore the global phase. The generalized Wigner function of the state ρ(t)=|ϕ(t)ϕ(t)|\rho(t)=\ket{\phi(t)}\bra{\phi(t)} is given by

W[ρ(t)]\displaystyle W[\rho(t)] =12π((13cosθ)e2|βαtγt|2+(1+3cosθ)e2|β+αtγt|2\displaystyle=\frac{1}{2\pi}\left((1-\sqrt{3}\cos\theta)e^{-2|\beta-\alpha_{t}-\gamma_{t}|^{2}}+(1+\sqrt{3}\cos\theta)e^{-2|\beta+\alpha_{t}-\gamma_{t}|^{2}}\right.
+23sinθe2|βγt|2cos(ϕ+4Im[αtβ])).\displaystyle\left.\qquad+2\sqrt{3}\sin\theta e^{-2|\beta-\gamma_{t}|^{2}}\cos(\phi+4\text{Im}[\alpha_{t}^{*}\beta])\right). (19)

Here, the reduced density matrix of the oscillator ρb(t)=Trq[ρ(t)]\rho_{b}(t)=\text{Tr}_{q}[\rho(t)] is

ρb(t)\displaystyle\rho_{b}(t) =12(ρb1(t)+ρb2(t)),\displaystyle=\frac{1}{2}(\rho_{b1}(t)+\rho_{b2}(t)), (20)
ρb1(t)\displaystyle\rho_{b1}(t) =|αt+γtαt+γt|,ρb2(t)=|αt+γtαt+γt|,\displaystyle=\ket{\alpha_{t}+\gamma_{t}}\bra{\alpha_{t}+\gamma_{t}},\quad\rho_{b2}(t)=\ket{-\alpha_{t}+\gamma_{t}}\bra{-\alpha_{t}+\gamma_{t}}, (21)

and the corresponding Wigner functions are

Wb[ρb1(t)]\displaystyle W_{b}[\rho_{b1}(t)] =12πe2|βαtγt|2,Wb[ρb2(t)]=12πe2|β+αtγt|2.\displaystyle=\frac{1}{2\pi}e^{-2|\beta-\alpha_{t}-\gamma_{t}|^{2}},\quad W_{b}[\rho_{b2}(t)]=\frac{1}{2\pi}e^{-2|\beta+\alpha_{t}-\gamma_{t}|^{2}}. (22)

Because this state is the mixed state of two Gaussian states, the negativity volume of the reduced oscillator’s state is zero. Hence, if the inequality

𝒱[ρ(t)]\displaystyle\mathcal{V}[\rho(t)] 1312,\displaystyle\leq\frac{1}{\sqrt{3}}-\frac{1}{2}, (23)

is violated, the state is expected to be entangled. Fig. 2 shows the behavior of the negativity volume for γ=0\gamma=0. We here use λ=101\lambda=10^{-1} to easily understand the behavior of the negativity volume of the Wigner function, although λ\lambda is much smaller in realistic situations. The red dashed line represents the critical value 𝒱cr\mathcal{V}_{\text{cr}} in the right-hand side of the inequality (23). The negativity volume at the initial time in the separable state is the same as the critical value. The result is consistent with the criterion that the state is entangled when the negativity volume is larger than the critical value 1/31/21/\sqrt{3}-1/2 in this case.

Refer to caption
Figure 2: The behavior of the negativity volume as a function of dimensionless time ωt/π\omega t/\pi for λ=101,γ=0\lambda=10^{-1},\gamma=0. The blue solid curve and the red dashed line represent the negativity volume 𝒱[ρ(t)]\mathcal{V}[\rho(t)] and the critical value 𝒱cr=1/31/2\mathcal{V}_{cr}=1/\sqrt{3}-1/2, respectively.

III.2 INITIAL QUBIT-THERMAL STATE

We next consider the thermal state as the initial state of the oscillator. In a realistic situation assuming an experiment, the thermal effect should be taken into account, and the generation of the gravity-induced entanglement is non-trivial. The entanglement in the hybrid system under the thermal effect has not been explicitly estimated [18, 19]. We perform the estimation through the negativity volume of the generalized Wigner function of the hybrid system.

The initial density matrix of the total system is

ρth(0)\displaystyle\rho^{th}(0) =12(|0+|1)(0|+1|)1πn¯d2γe|γ|2/n¯|γγ|,\displaystyle=\frac{1}{2}(\ket{0}+\ket{1})(\bra{0}+\bra{1})\otimes\frac{1}{\pi\bar{n}}\int d^{2}\gamma e^{-|\gamma|^{2}/\bar{n}}\ket{\gamma}\bra{\gamma}, (24)

where n¯=kBT/2ω\bar{n}=k_{B}T/2\omega with the temperature TT is thermal phonon occupancy number. The evolution of the density matrix by the Hamiltonian (12) is

ρth(t)\displaystyle\rho^{th}(t) =12πn¯d2γe|γ|2/n¯{|0|αt+γt0|αt+γt|+e(αtγαtγ)|0|αt+γt1|αt+γt|\displaystyle=\frac{1}{2\pi\bar{n}}\int d^{2}\gamma e^{-|\gamma|^{2}/\bar{n}}\left\{\ket{0}\ket{\alpha_{t}+\gamma_{t}}\bra{0}\bra{\alpha_{t}+\gamma_{t}}+e^{-(\alpha_{t}^{*}\gamma^{*}-\alpha_{t}\gamma)}\ket{0}\ket{\alpha_{t}+\gamma_{t}}\bra{1}\bra{-\alpha_{t}+\gamma_{t}}\right.
+eαtγαtγ|1|αt+γt0|αt+γt|+|1|αt+γt1|αt+γt|},\displaystyle\qquad\qquad\qquad\left.+e^{\alpha_{t}^{*}\gamma^{*}-\alpha_{t}\gamma}\ket{1}\ket{-\alpha_{t}+\gamma_{t}}\bra{0}\bra{\alpha_{t}+\gamma_{t}}+\ket{1}\ket{-\alpha_{t}+\gamma_{t}}\bra{1}\bra{-\alpha_{t}+\gamma_{t}}\right\}, (25)

and the generalized Wigner function is

W[ρth(t)]=12π(2n¯+1)\displaystyle W[\rho^{th}(t)]=\frac{1}{2\pi(2\bar{n}+1)} {(13cosθ)e22n¯+1|βαt|2+(1+3cosθ)e22n¯+1|β+αt|2\displaystyle\left\{(1-\sqrt{3}\cos\theta)e^{-\frac{2}{2\bar{n}+1}|\beta-\alpha_{t}|^{2}}+(1+\sqrt{3}\cos\theta)e^{-\frac{2}{2\bar{n}+1}|\beta+\alpha_{t}|^{2}}\right.
+23e22n¯+1|β|2sinθcos[ϕ+4Im[αtβ]]}.\displaystyle\left.+2\sqrt{3}e^{-\frac{2}{2\bar{n}+1}|\beta|^{2}}\sin\theta\cos[\phi+4\text{Im}[\alpha_{t}^{*}\beta]]\right\}. (26)

Then, the reduced density matrix of the oscillator ρbth\rho^{th}_{b} is

ρbth(t)\displaystyle\rho^{th}_{b}(t) =12(ρb1th(t)+ρb2th(t)),\displaystyle=\frac{1}{2}(\rho^{th}_{b1}(t)+\rho^{th}_{b2}(t)),
ρb1th(t)\displaystyle\rho^{th}_{b1}(t) =1πn¯d2γe|γ|2/n¯|αt+γtαt+γt|,ρb2th(t)=1πn¯d2γe|γ|2/n¯|αt+γtαt+γt|,\displaystyle=\frac{1}{\pi\bar{n}}\int d^{2}\gamma e^{-|\gamma|^{2}/\bar{n}}\ket{\alpha_{t}+\gamma_{t}}\bra{\alpha_{t}+\gamma_{t}},\quad\rho^{th}_{b2}(t)=\frac{1}{\pi\bar{n}}\int d^{2}\gamma e^{-|\gamma|^{2}/\bar{n}}\ket{-\alpha_{t}+\gamma_{t}}\bra{-\alpha_{t}+\gamma_{t}},

and the Wigner function is obtained by

Wb[ρb1th]\displaystyle W_{b}[\rho_{b1}^{th}] =12π(2n¯+1)e12n¯+1|βαt|2,Wb[ρb2th]=12π(2n¯+1)e12n¯+1|β+αt|2.\displaystyle=\frac{1}{2\pi(2\bar{n}+1)}e^{-\frac{1}{2\bar{n}+1}|\beta-\alpha_{t}|^{2}},\quad W_{b}[\rho_{b2}^{th}]=\frac{1}{2\pi(2\bar{n}+1)}e^{-\frac{1}{2\bar{n}+1}|\beta+\alpha_{t}|^{2}}. (27)

These Wigner functions of the reduced state are always non-negative and the negativity volume is zero, 𝒱[ρb1th]=𝒱[ρb2th]=0\mathcal{V}[\rho_{b1}^{th}]=\mathcal{V}[\rho_{b2}^{th}]=0. Then, the inequality (11) reduces to the inequality (23). Hence, the violation of the inequality (23) is the criterion for the entanglement generation due to gravity. Figure 3 shows the behavior of the negativity volume and critical value for λ=101,n¯=3\lambda=10^{-1},\bar{n}=3. The entanglement behavior for the initial qubit-thermal state is similar to that for the initial qubit-coherent state, but the maximum value of the negativity volume is smaller due to the thermal effects. However, the criterion based on the negativity volume of the generalized Wigner function detects the entanglement even for the mixed state in the hybrid system.

Refer to caption
Figure 3: The behavior of the negativity volume (blue curve) as a function of dimensionless time ωt/π\omega t/\pi for the initial qubit-thermal state with λ=101,n¯=3\lambda=10^{-1},\bar{n}=3. The negativity volume exhibits the same periodic behavior as that in the initial qubit-coherent state of Fig. 2, though the peak value is smaller. The red dashed line is the critical value 1/31/21/\sqrt{3}-1/2 in this case.
Refer to caption
Figure 4: The negativity volume for the initial qubit-cat state with λ=101,γ=1\lambda=10^{-1},\gamma=1. The red dashed curve is the right-hand side of Eq. (36).

III.3 INITIAL QUBIT-CAT STATE

The results in the previous subsections are the case that the Wigner function of the reduced state of the oscillator is non-negative. We here consider the cat state of the oscillator, where its Wigner function has negative values and the negativity volume of the Wigner function is non-zero. We investigate this model to verify that the inequality (11) can be used as an entanglement criterion of a hybrid system.

The initial density matrix of the total system is

ρcat(0)\displaystyle\rho^{cat}(0) =|ϕcat(0)ϕcat(0)|,|ϕcat(0)=12N(|0+|1)(|γ+|γ),\displaystyle=\ket{\phi^{cat}(0)}\bra{\phi^{cat}(0)},\qquad\ket{\phi^{cat}(0)}=\frac{1}{\sqrt{2N}}(\ket{0}+\ket{1})\otimes(\ket{\gamma}+\ket{-\gamma}), (28)

where N=2+2e2|γ|2N=2+2e^{-2|\gamma|^{2}}. The evolved state by unitary operator (15) is given by

|ϕcat(t)\displaystyle\ket{\phi^{cat}(t)} =12N{|0(e(αtγαtγ)/2|αt+γt+e(αtγαtγ)/2|αtγt)\displaystyle=\frac{1}{\sqrt{2N}}\left\{\ket{0}(e^{-(\alpha_{t}^{*}\gamma^{*}-\alpha_{t}\gamma)/2}\ket{\alpha_{t}+\gamma_{t}}+e^{(\alpha_{t}^{*}\gamma^{*}-\alpha_{t}\gamma)/2}\ket{\alpha_{t}-\gamma_{t}})\right.
+|1(e(αtγαtγ)/2|αt+γt+e(αtγαtγ)/2|αtγt)},\displaystyle\left.\qquad\quad+\ket{1}(e^{(\alpha_{t}^{*}\gamma^{*}-\alpha_{t}\gamma)/2}\ket{-\alpha_{t}+\gamma_{t}}+e^{-(\alpha_{t}^{*}\gamma^{*}-\alpha_{t}\gamma)/2}\ket{-\alpha_{t}-\gamma_{t}})\right\}, (29)

and the generalized Wigner function is

W[ρcat(t)]\displaystyle W[\rho^{cat}(t)] =1Nπ[13cosθ2(e2|βαtγt|2+e2|βαt+γt|2+2e2|βαt|2cos[4Im[αtγ+γtβ]])\displaystyle=\frac{1}{N\pi}\left[\frac{1-\sqrt{3}\cos\theta}{2}\left(e^{-2|\beta-\alpha_{t}-\gamma_{t}|^{2}}+e^{-2|\beta-\alpha_{t}+\gamma_{t}|^{2}}+2e^{-2|\beta-\alpha_{t}|^{2}}\cos[4\text{Im}[\alpha_{t}^{*}\gamma^{*}+\gamma_{t}^{*}\beta]]\right)\right.
+1+3cosθ2(e2|β+αtγt|2+e2|β+αt+γt|2+2e2|β+αt|2cos[4Im[αtγγtβ]])\displaystyle\left.\qquad\quad~{}+\frac{1+\sqrt{3}\cos\theta}{2}\left(e^{-2|\beta+\alpha_{t}-\gamma_{t}|^{2}}+e^{-2|\beta+\alpha_{t}+\gamma_{t}|^{2}}+2e^{-2|\beta+\alpha_{t}|^{2}}\cos[4\text{Im}[\alpha_{t}^{*}\gamma^{*}-\gamma_{t}^{*}\beta]]\right)\right.
+3sinθcos[4Im[αtβ+ϕ]](e2|βγt|2+e2|β+γt|2+2e2|β|2cos[4Im[γtβ]])].\displaystyle\left.\qquad\quad+\sqrt{3}\sin\theta\cos[4\text{Im}[\alpha_{t}^{*}\beta+\phi]]\left(e^{-2|\beta-\gamma_{t}|^{2}}+e^{-2|\beta+\gamma_{t}|^{2}}+2e^{-2|\beta|^{2}}\cos[4\text{Im}[\gamma_{t}^{*}\beta]]\right)\right]. (30)

Tracing over the qubit system, the reduced density matrix of the oscillator is derived by

ρbcat(t)\displaystyle\rho^{cat}_{b}(t) =12(ρb1cat(t)+ρb2cat(t)),\displaystyle=\frac{1}{2}(\rho^{cat}_{b1}(t)+\rho^{cat}_{b2}(t)), (31)
ρb1cat(t)\displaystyle\rho^{cat}_{b1}(t) =1N(|αt+γtαt+γt|+e(αtγαtγ)|αt+γtαtγt|\displaystyle=\frac{1}{N}\left(\ket{\alpha_{t}+\gamma_{t}}\bra{\alpha_{t}+\gamma_{t}}+e^{-(\alpha_{t}^{*}\gamma^{*}-\alpha_{t}\gamma)}\ket{\alpha_{t}+\gamma_{t}}\bra{\alpha_{t}-\gamma_{t}}\right.
+eαtγαtγ|αtγtαt+γt|+|αtγtαtγt|),\displaystyle\left.\qquad~{}+e^{\alpha_{t}^{*}\gamma^{*}-\alpha_{t}\gamma}\ket{\alpha_{t}-\gamma_{t}}\bra{\alpha_{t}+\gamma_{t}}+\ket{\alpha_{t}-\gamma_{t}}\bra{\alpha_{t}-\gamma_{t}}\right), (32)
ρb2cat(t)\displaystyle\rho^{cat}_{b2}(t) =1N(|αt+γtαt+γt|+eαtγαtγ|αt+γtαtγt|\displaystyle=\frac{1}{N}\left(\ket{-\alpha_{t}+\gamma_{t}}\bra{-\alpha_{t}+\gamma_{t}}+e^{\alpha_{t}^{*}\gamma^{*}-\alpha_{t}\gamma}\ket{-\alpha_{t}+\gamma_{t}}\bra{-\alpha_{t}-\gamma_{t}}\right.
+e(αtγαtγ)|αtγtαt+γt|+|αtγtαtγt|),\displaystyle\left.\qquad~{}+e^{-(\alpha_{t}^{*}\gamma^{*}-\alpha_{t}\gamma)}\ket{-\alpha_{t}-\gamma_{t}}\bra{-\alpha_{t}+\gamma_{t}}+\ket{-\alpha_{t}-\gamma_{t}}\bra{-\alpha_{t}-\gamma_{t}}\right), (33)

and the Wigner function is

Wb[ρb1cat(t)]\displaystyle W_{b}[\rho^{cat}_{b1}(t)] =2Nπ(e2|βαtγt|2+e2|βαt+γt|2+2e2|βαt|2cos[4Im[αtγ+γtβ]]),\displaystyle=\frac{2}{N\pi}\left(e^{-2|\beta-\alpha_{t}-\gamma_{t}|^{2}}+e^{-2|\beta-\alpha_{t}+\gamma_{t}|^{2}}+2e^{-2|\beta-\alpha_{t}|^{2}}\cos[4\text{Im}[\alpha_{t}^{*}\gamma^{*}+\gamma_{t}^{*}\beta]]\right), (34)
Wb[ρb2cat(t)]\displaystyle W_{b}[\rho^{cat}_{b2}(t)] =2Nπ(e2|β+αtγt|2+e2|β+αt+γt|2+2e2|β+αt|2cos[4Im[αtγγtβ]]).\displaystyle=\frac{2}{N\pi}\left(e^{-2|\beta+\alpha_{t}-\gamma_{t}|^{2}}+e^{-2|\beta+\alpha_{t}+\gamma_{t}|^{2}}+2e^{-2|\beta+\alpha_{t}|^{2}}\cos[4\text{Im}[\alpha_{t}^{*}\gamma^{*}-\gamma_{t}^{*}\beta]]\right). (35)

Here, the Wigner function of each reduced state has negative regions because the reduced state is the non-Gaussian state. Hence, the inequality for testing the entanglement is

𝒱[ρcat(t)]\displaystyle\mathcal{V}[\rho^{cat}(t)] 23i=1212𝒱[ρbicat(t)]+1312.\displaystyle\leq\frac{2}{\sqrt{3}}\sum_{i=1}^{2}\frac{1}{2}\mathcal{V}[\rho^{cat}_{bi}(t)]+\frac{1}{\sqrt{3}}-\frac{1}{2}. (36)

Figure 4 shows the negativity volume (blue curve) and critical value (red dashed curve) for λ=101\lambda=10^{-1} and γ=1\gamma=1. The critical value (red dashed curve) is a function of time, but it is almost constant for the parameters. The negativity volume is larger than the critical value except for ωt/π=2n(n=0,1,2,)\omega t/\pi=2n\,(n=0,1,2,\dots). Because the total system is in a pure state, the entanglement should be generated. The result indicates that the entanglement generation due to the gravitational interaction is captured by the negativity volume of the generalized Wigner function.

Refer to caption
Refer to caption
Refer to caption
Figure 5: The behaviors of the fidelity for the initial qubit-coherent state (upper left), qubit-thermal state (upper right), and qubit-cat state (bottom), adopting the same parameters as Figs. 2, 3, and 4, respectively. The state is entangled if and only if the fidelity is not equal to 1.

IV ENTANGLEMENT FIDELITY

In this section, we discuss how the evaluation of the entanglement by the Wigner function is consistent with that by the fidelity proposed in Ref. [54, 55, 56]. We here rewrite the Hamiltonian (12) as

H^\displaystyle\hat{H} =ωa^a^+|00|g(a^+a^)+|11|(g(a^+a^)).\displaystyle=\omega\hat{a}^{\dagger}\hat{a}+\ket{0}\bra{0}\otimes g(\hat{a}+\hat{a}^{\dagger})+\ket{1}\bra{1}\otimes(-g(\hat{a}+\hat{a}^{\dagger})). (37)

According to the entanglement criterion Ref. [54, 55, 56], the state is separable if and only if the equality

[ρb(0),ω^(t)]\displaystyle[\rho_{b}(0),\hat{\omega}(t)] =0\displaystyle=0 (38)

is satisfied, where the operator ω^(t)\hat{\omega}(t) is given by

ω^(t)\displaystyle\hat{\omega}(t) =ei(ωa^a^+g(a^+a^))tei(ωa^a^g(a^+a^))t,\displaystyle=e^{i(\omega\hat{a}^{\dagger}\hat{a}+g(\hat{a}+\hat{a}^{\dagger}))t}e^{-i(\omega\hat{a}^{\dagger}\hat{a}-g(\hat{a}+\hat{a}^{\dagger}))t}, (39)

and we define the initial state as ρ(0)=ρqp(0)ρb(0)\rho(0)=\rho_{q}^{p}(0)\otimes\rho_{b}(0). This criterion is applicable to a initial product state of the pure qubit system and any other system, and is applicable to all the systems analyzed in the previous section. Eq. (38) is rewritten as

ei(ωa^a^+g(a^+a^))tρb(0)ei(ωa^a^+g(a^+a^))t\displaystyle e^{-i(\omega\hat{a}^{\dagger}\hat{a}+g(\hat{a}+\hat{a}^{\dagger}))t}\rho_{b}(0)e^{i(\omega\hat{a}^{\dagger}\hat{a}+g(\hat{a}+\hat{a}^{\dagger}))t} =ei(ωa^a^g(a^+a^))tρb(0)ei(ωa^a^g(a^+a^))t.\displaystyle=e^{-i(\omega\hat{a}^{\dagger}\hat{a}-g(\hat{a}+\hat{a}^{\dagger}))t}\rho_{b}(0)e^{i(\omega\hat{a}^{\dagger}\hat{a}-g(\hat{a}+\hat{a}^{\dagger}))t}. (40)

and this equation holds when the following operators σ0\sigma_{0} and σ1\sigma_{1} are equal:

σ0\displaystyle\sigma_{0} =eiωa^a^tei(αta^αta^)ρb(0)ei(αta^αta^)eiωa^a^t,\displaystyle=e^{-i\omega\hat{a}^{\dagger}\hat{a}t}e^{-i(\alpha_{t}^{*}\hat{a}^{\dagger}-\alpha_{t}\hat{a})}\rho_{b}(0)e^{i(\alpha_{t}^{*}\hat{a}^{\dagger}-\alpha_{t}\hat{a})}e^{i\omega\hat{a}^{\dagger}\hat{a}t}, (41)
σ1\displaystyle\sigma_{1} =eiωa^a^tei(αta^αta^)ρb(0)ei(αta^αta^)eiωa^a^t.\displaystyle=e^{-i\omega\hat{a}^{\dagger}\hat{a}t}e^{i(\alpha_{t}^{*}\hat{a}^{\dagger}-\alpha_{t}\hat{a})}\rho_{b}(0)e^{-i(\alpha_{t}^{*}\hat{a}^{\dagger}-\alpha_{t}\hat{a})}e^{i\omega\hat{a}^{\dagger}\hat{a}t}. (42)

To evaluate the identification of the operators σ0\sigma_{0} and σ1\sigma_{1}, we introduce the fidelity as

F(σ0,σ1)\displaystyle F(\sigma_{0},\sigma_{1}) =(Tr[σ1σ0σ1])2,\displaystyle=\left(\text{Tr}\left[\sqrt{\sqrt{\sigma_{1}}\sigma_{0}\sqrt{\sigma_{1}}}\right]\right)^{2}, (43)

and the state is entangled if and only if Eq. (40) is not satisfied, namely, F1F\neq 1.

For the initial qubit-coherent state (17) in Sec. III-A, we derive the fidelity as

F(σ0,σ1)\displaystyle F(\sigma_{0},\sigma_{1}) =e2|αt|2,\displaystyle=e^{-2|\alpha_{t}|^{2}}, (44)

which is not dependent on γ\gamma. Similarly, the fidelity for the initial qubit-thermal state (24) is obtained [57] as

F(σ0th,σ1th)\displaystyle F(\sigma_{0}^{th},\sigma_{1}^{th}) =e22n¯+1|αt|2,\displaystyle=e^{-\frac{2}{2\bar{n}+1}|\alpha_{t}|^{2}}, (45)

and the fidelity for the initial qubit-cat state is

F(σ0cat(t),σ1cat(t))\displaystyle F(\sigma_{0}^{cat}(t),\sigma_{1}^{cat}(t)) =1N|e2|αtγt|2+e2|αt+γt|2+2e2|αt|2cos(4Im[αtγ])|,\displaystyle=\frac{1}{N}\left|e^{-2|\alpha_{t}-\gamma_{t}|^{2}}+e^{-2|\alpha_{t}+\gamma_{t}|^{2}}+2e^{-2|\alpha_{t}|^{2}}\cos\left(4\text{Im}[\alpha_{t}^{*}\gamma^{*}]\right)\right|, (46)

where N=2+2e2|γ|2N=2+2e^{-2|\gamma|^{2}}.

Each panel of Fig. 5 plots the behavior of the fidelity (44), (45), and (46), as a function of time ωt/π\omega t/\pi. We find that the region where F(σ0,σ1)1F(\sigma_{0},\sigma_{1})\neq 1 is the same as the region where the negativity volume is larger than the critical value in Figs. 234. Moreover, the fidelity for the qubit-thermal state is closer to 11 than that of the qubit-coherent state, which corresponds to the fact that the negativity volume of the qubit-thermal state is smaller than that of the initial qubit-coherent state. Thus the comparison of the conditions for quantum entanglement demonstrated the consistency between the method with the negativity volume of the generalized Wigner function and the method with fidelity, both of which properly detect the quantum entanglement in the hybrid system. In all these models, the signature of entanglement is completely captured by the negativity volume of the generalized Wigner function.

V SUMMARY AND CONCLUSION

We examined the quantum nature of gravity in terms of the generalized Wigner function in a hybrid system consisting of a two-level particle and an oscillator. The two-level particle means a spatially localized superposed state at two different points and the total system evolves into nonclassical entangled states through gravitational interaction. We discussed the entanglement due to gravity by using the quantum entanglement criterion, which states that an increase in the negativity volume of the generalized Wigner function corresponds to quantum entanglement, as proposed in Ref. [10]. Although this criterion is not a necessary condition, we demonstrated that the entanglement is adequately captured by comparing it with the condition using the fidelity. Our results suggest that the criterion with the generalized Wigner function is valid for detecting the gravity-induced entanglement, though the increase in the negativity volume is small in a realistic situation because gravity is weak.

Acknowledgements.
D. Miki thanks Prof. Kae Nemoto for her hospitality during his stay at OIST and for useful discussions on the topic of the present paper. D. M. is supported by JSPS KAKENHI, Grant No. 22J21267, K.Y. is supported by JSPS KAKENHI, Grant No. JP22H05263 and No. JP23H01175.

Appendix A ENTANGLEMENT CRITERION BASED ON THE GENERALIZED WIGNER FUNCTION

We explain the entanglement criterion proposed in Ref. [10]. Here, we consider the separable state of a hybrid system ρsep=ipiρqiρbi\rho^{\text{sep}}=\sum_{i}p_{i}\rho_{q}^{i}\otimes\rho_{b}^{i}, where pip_{i} is positive values satisfying ipi=1\sum_{i}p_{i}=1. The generalized Wigner function of the separable system is

W[ρsep]\displaystyle W[\rho^{\text{sep}}] =ipitr[ρqiΔq^]tr[ρbiΔb^]=ipiWq[ρqi]Wb[ρbi],\displaystyle=\sum_{i}p_{i}\text{tr}[\rho_{q}^{i}\hat{\Delta_{q}}]\text{tr}[\rho_{b}^{i}\hat{\Delta_{b}}]=\sum_{i}p_{i}W_{q}[\rho_{q}^{i}]W_{b}[\rho_{b}^{i}], (47)

and the negativity volume is

𝒱[ρsep]\displaystyle\mathcal{V}[\rho^{\text{sep}}] =12(𝑑ϕ𝑑θd2βsinθ2π|W[ρsep]|1)\displaystyle=\frac{1}{2}\left(\int d\phi d\theta d^{2}\beta\frac{\sin\theta}{2\pi}|W[\rho^{\text{sep}}]|-1\right)
12(ipi𝑑ϕ𝑑θd2βsinθ2π|Wq[ρqi]Wb[ρbi]|1)\displaystyle\leq\frac{1}{2}\left(\sum_{i}p_{i}\int d\phi d\theta d^{2}\beta\frac{\sin\theta}{2\pi}|W_{q}[\rho_{q}^{i}]W_{b}[\rho_{b}^{i}]|-1\right)
=12ipi{(dϕdθsinθ2π|Wq[ρqi]|1)(dβ|Wb[ρbi]|1)\displaystyle=\frac{1}{2}\sum_{i}p_{i}\left\{\left(\int d\phi d\theta\frac{\sin\theta}{2\pi}|W_{q}[\rho_{q}^{i}]|-1\right)\left(\int d\beta|W_{b}[\rho_{b}^{i}]|-1\right)\right.
+(dϕdθsinθ2π|Wq[ρqi]|1)+(dβ|Wb[ρbi]|1)}\displaystyle\left.\quad+\left(\int d\phi d\theta\frac{\sin\theta}{2\pi}|W_{q}[\rho_{q}^{i}]|-1\right)+\left(\int d\beta|W_{b}[\rho_{b}^{i}]|-1\right)\right\}
=ipi(2𝒱q[ρqi]𝒱b[ρbi]+𝒱q[ρqi]+𝒱b[ρbi]).\displaystyle=\sum_{i}p_{i}\left(2\mathcal{V}_{q}[\rho_{q}^{i}]\mathcal{V}_{b}[\rho_{b}^{i}]+\mathcal{V}_{q}[\rho_{q}^{i}]+\mathcal{V}_{b}[\rho_{b}^{i}]\right). (48)

Since the negativity volume of the Wigner function for the qubit system is upper bounded by the pure state, 1/31/21/\sqrt{3}-1/2, then the negativity volume of the separable state of the hybrid system satisfies the following inequality:

𝒱[ρsep]\displaystyle\mathcal{V}[\rho^{\text{sep}}] 23ipi𝒱[ρbi]+1312.\displaystyle\leq\frac{2}{\sqrt{3}}\sum_{i}p_{i}\mathcal{V}[\rho_{b}^{i}]+\frac{1}{\sqrt{3}}-\frac{1}{2}. (49)

Hence, if the inequality is violated, it means the nonclassicality of the hybrid system. We note that this criterion is the only sufficient criterion for the entanglement generation and there are cases where the inequality (49) is satisfied even though the state is quantum entangled, as shown in Ref. [10]. Also, there is no restriction of the negativity volume for the reduced continuous variable system 𝒱b[ρbi]\mathcal{V}_{b}[\rho_{b}^{i}].

References

  • [1] E. A Calzetta and B-L. B. Hu, Nonequilibrium quantum field theory, (Cambridge University Press, New York, 2008)
  • [2] A. Kenfact and K. Życzkowski, Negativity of the Wigner function as an indicator of non-classicality, J. Opt. B: Quantum Semiclass. Opt. 6 (2004) 396-404
  • [3] K. Banaszek, C. Radzewicz, K. Wódkiewicz, and J. S. Krasiński Direct measurement of the Wigner function by photon counting Phys. Rev. A 𝟔𝟎\bm{60}, 674 (1999)
  • [4] A. I. Lvovsky, H. Hansen, T. Aichele, O. Benson, J. Mlynek, and S. Schiller Quantum State Reconstruction of the Single-Photon Fock State Phys. Rev. Lett. 𝟖𝟕\bm{87}, 050402 (2001)
  • [5] Alessandro Zavatta, Silvia Viciani, and Marco Bellini, Tomographic reconstruction of the single-photon Fock state by high-frequency homodyne detection, Phys. Rev. A 𝟕𝟎\bm{70}, 053821 (2004)
  • [6] A. I. Lvovsky and M. G. Raymer, Continuous-variable optical quantum-state tomography, Rev. Mod. Phys. 𝟖𝟏\bm{81}, 299 (2009)
  • [7] T. Tilma, M. J. Everitt, J. H. Samson, W. J. Munro, and K. Nemoto, Wigner functions for arbitrary quantum system, Phys. Rev. Lett. 117, 180401 (2016)
  • [8] R. P. Rundle, P. W. Mills, T. Tilma, J. H. Samson, and M. J. Everitt, Simple procedure for phase-space measurement and entanglement validation, Phys. Rev. A 96, 022117 (2017)
  • [9] Y. Tian, Z. Wang, P. Zhang, G. Li, J. Li, and T. Zhang, Measurement of complete and continuous Wigner functions for discrete atomic systems, Phys. Rev. A 97, 013840 (2018)
  • [10] I. I. Arkhipov, A. Barasiński, and J. Svozilik, Negativity volume of the generalized Wigner function as an entanglement witness for hybrid bipartite states, Sci. Rep. 8, 16955 (2018)
  • [11] R. P. Rundle, T. Tilma, J. H. Samson, V. M. Dwyer, R. F. Bishop, and M. J. Everitt, General approach to quantum mechanics as a statistical theory, Phys. Rev. A 99, 012115 (2019)
  • [12] B. I. Davies, R. P. Rundle, V. M. Dwyer, J. H. Samson, T. Tilma, and M. J. Everitt, Visualizing spin degrees of freedom in atoms and molecules, Phys. Rev. A 100, 042102 (2019)
  • [13] U. L. Andersen, J. S. Neergaard-Nielsen, P. van Loock, and A. Furusawa, Hybrid discrete- and continuous-variable quantum information, Nature Phys. 𝟏𝟏\bm{11}, 713 (2015)
  • [14] N. Gisin and R. Thew, Quantum communication, Nature Photon. 𝟏\bm{1}, 165 (2007)
  • [15] A. E. Ulanov, D. Sychev, A. A. Pushkina, I. A. Fedorov, and A. I. Lvovsky, Quantum teleportation between discrete and continuous encodings of an optical qubit, Phys. Rev. Lett. 𝟏𝟏𝟖\bm{118}, 160501 (2017)
  • [16] K. Kreis and P. van Loock, Classifying, quantifying, and witnessing qudit-qumode hybrid entanglement, Phys. Rev. A 𝟖𝟓\bm{85}, 032307 (2012)
  • [17] H. Jeong, A. Zavatta, M. Kang, Seung-Woo Lee, L. S. Costanzo, S. Grandi, T. C. Ralph, and M. Bellini, Generation of hybrid entanglement of light, Nature Photon. 𝟖\bm{8}, 564 (2014)
  • [18] D. Carney, H. Müller, and J. M. Taylor, Using an Atom Interferometer to Infer Gravitational Entanglement Generation, PRX Quantum 2, 030330 (2021)
  • [19] A. Matsumura and K. Yamamoto, Leggett-Garg inequality for testing quantumness of gravity, Phys. Rev. A 106, 012214 (2022)
  • [20] K. Strltsov, J. S. Pedernales, and M. B. Plenio, On the Significance of Interferometric Revivals for the Fundamental Description of Gravity, Universe 2022, 8(2) 58
  • [21] Yue Ma, Thomas Guff, Gavin Morley, Igor Pikovski, and M. S. Kim, Limits on inference of gravitational entanglement, Phys. Rev. Research 4, 013024 (2022)
  • [22] S. Bose, A. Mazumdar, G. W. Morley, H. Ulbricht, M. Toroš, M. Paternostro, A. A. Geraci, P. F. Barker, M. S. Kim, and G. Milburn, Spin Entanglement Witness for Quantum Gravity, Phys. Rev. Lett. 𝟏𝟏𝟗\bm{119}, 240401 (2017)
  • [23] C. Marletto and V. Vedral, Gravitationally Induced Entanglement between two Massive Particles is Sufficient Evidence of Quantum Effects in Gravity, Phys. Rev. Lett. 𝟏𝟏𝟗\bm{119}, 240402 (2017)
  • [24] R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, Quantum Entanglement, Rev, Mod. Phys. 81, 865 (2009)
  • [25] H. Chau Nguyen and F. Bernards, Entanglement dynamics of two mesoscopic objects with gravitational interaction, Eur. Phys. J. D 74, 69 (2020)
  • [26] D. Miki, A. Matsumura, and K. Yamamoto, Entanglement and decoherence of massive particles due to gravity, Phys. Rev. D 103, 026017 (2021)
  • [27] A. Matsumura, Field-induced entanglement in spatially superposed objects, Phys. Rev. D 104, 046001 (2021)
  • [28] Y. Sugiyama, A. Matsumura, and K. Yamamoto, Effects of photon field on entanglement generation in charged particles, Phys. Rev. D 104, 046001 (2021)
  • [29] T. Feng and V. Vedral, Amplification of gravitationally induced entanglement, Phys. Rev. D 106, 066013 (2022)
  • [30] H. Miao, D. Martynov, H. Yang, and A. Datta, Quantum correlations of light mediated by gravity, Phys. Rev. A 𝟏𝟎𝟏\bm{101}, 063804 (2020)
  • [31] T. Krisnanda, G. Y. Tham, M. Paternostro, and T. Paterek, Observable quantum entanglement due to gravity, Quantum Inf. 𝟔\bm{6}, 12 (2020)
  • [32] A. Matsumura and K. Yamamoto, Gravity-induced entanglement in optomechanical systems, Phys. Rev. D 𝟏𝟎𝟐\bm{102}, 106021 (2020)
  • [33] A. Datta and H. Miao, Signatures of the quantum nature of gravity in the differential motion of two masses, Quantum Sci. Technol. 𝟔\bm{6}, 045014 (2021)
  • [34] D. Miki, A. Matsumura and K. Yamamoto, Non-Gaussian entanglement in gravitating masses: The role of cumulants, Phys. Rev. D 𝟏𝟎𝟓\bm{105}, 026011 (2022)
  • [35] A. D. K. Plato, D. Rätzel, and C. Wan, Enhanced Gravitational Entanglement via Modulated Optomechanics, Quantum 7, 1177 (2023)
  • [36] N. Matsumoto, S. B. Cataño-Lopez, M. Sugawara, S. Suzuki, N. Abe, K. Komori, Y. Michimura, Y. Aso, and K. Edamatsu, Demonstration of Displacement Sensing of a mg-Scale Pendulum for mm- and mg-Scale Gravity Measurements, Phys. Rev. Lett. 122, 071101 (2019)
  • [37] S. B. Cataño-Lopez, J. G. Santiago-Condori, K. Edamatsu, and N. Matsumoto, High-Q Milligram-Scale Monolithic Pendulum for Quantum-Limited Gravity Measurements, Phys. Rev. Lett. 124, 221102 (2020)
  • [38] N. Matsumoto and N. Yamamoto, Preparing mechanical squeezing of a macroscopic pendulum near quantum regimes, arXiv:2008.10848
  • [39] Y. Chen, Macroscopic Quantum Mechanics: Theory and Experimental Concepts of Optomechanics, J. Phys. B: At. Mol. Opt. Phys. 46, 104001 (2013)
  • [40] W. P. Bowen, G. J. Bilburn, Quantum Optomechanics, (CRC Press, 2015)
  • [41] H. Rudolph, U. Delić, M. Aspelmeyer, K. Hornberger, ans B. A. Stickler, Force-Gradient Sensing and Entanglement via Feedback Cooling of Interacting Nanoparticles, Phys. Rev. Lett. 129, 193602 (2022)
  • [42] D. Miki, N. Matsumoto, A. Matsumura, T. Shichijo, Y. Sugiyama, K. Yamamoto, and N. Yamamoto, Generating quantum entanglement between macroscopic objects with continuous measurement and feedback control, Phys. Rev. A 107, 032410 (2023)
  • [43] D. Miki, A. Matsumura and K. Yamamoto, Quantum signature of gravity in optomechanical systems with conditional measurement, arXiv:2311.00563
  • [44] A. Belenchia, R. M. Wald, F. Giacomini, E. Castro-Ruiz, C. Brukner, and M. Aspelmeyer, Quantum superposition of massive objects and the quantization of gravity, Phys. Rev. D 98, 126009 (2018)
  • [45] R.J. Marshman, A. Mazumdar, and S. Bose, Locality and entanglement in table-top testing of the quantum nature of linearized gravity, Phys. Rev. A 101, 052110 (2020)
  • [46] D. Carney, Newton, entanglement, and the graviton, Phys. Rev. D 105, 024029 (2022)
  • [47] D. L. Danielson, G. Satishchandran, and R. M. Wald, Gravitationally Mediated Entanglement: Newtonian Field vs. Gravitons, Phys. Rev. D 105, 086001 (2022)
  • [48] S. Bose, A. Mazumdar, M. Schut, and M. Toroš, Mechanism for the quantum natured gravitons to entangle masses, Phys. Rev. D 105, 106028 (2022)
  • [49] Y. Hidaka, S. Iso, and K. Shimada, Complementarity and causal propagation of decoherence by measurement in relativistic quantum field theories, Phys. Rev. D 106, 076018 (2022)
  • [50] Y. Sugiyama, A. Matsumura, and K. Yamamoto, Consistency between causality and complementarity guaranteed by Robertson inequality in quantum field theory, Phys. Rev. D 106, 125002 (2022)
  • [51] D. Biswas, S. Bose, A. Mazumdar, and M. Toroš, Gravitational Optomechanics: Photon-Matter Entanglement via Graviton Exchange, Phys. Rev. D 108, 064023 (2023)
  • [52] Y. Hidaka, S. Iso, and K. Shimada, Entanglement Generation and Decoherence in a Two-Qubit System Mediated by Relativistic Quantum Field, Phys. Rev. D 107, 085003 (2023)
  • [53] Y. Sugiyama, A. Matsumura, and K. Yamamoto, Quantum uncertainty of gravitational field and entanglement in superposed massive particles, arXiv:2308.03093
  • [54] K. Roszak and Ł. Cywiński, Characterization and measurement of qubit-environment-entanglement generation during pure dephasing, Phys. Rev. A 92, 032310 (2015)
  • [55] K. Roszak, Criteria for system-environment entanglement generation for systems of any size in pure-dephasing evolutions, Phys. Rev. A 98, 052344 (2018)
  • [56] M. Strzałka, D. Kwiatkowski, Ł. Cywiński, and K. Roszak, Qubit-environment negativity versus fidelity of conditional environmental states for a nitrogen-vacancy-center spin qubit interacting with a nuclear environment, Phys. Rev. A 102, 042602 (2020)
  • [57] P. Marian, optimal purifications and fidelity for displaced thermal states, Phys. Rev. A 76, 054307 (2007)