# Heralded entangled coherent states between spatially separated massive resonators

## Abstract

We put forward an experimentally feasible scheme for heralded entanglement generation between two distant macroscopic mechanical resonators. The protocol exploits a hybrid quantum device, a qubit interacting with a mechanical resonator as well as a cavity mode, for each party. The cavity modes interfere on a beam-splitter followed by suitable heralding detections which post-select a hybrid entangled state with success probability 1/2. Subsequently, by local measurements on the qubits a mechanical entangled coherent state can be achieved. The mechanical entanglement can be further verified via monitoring the entanglement of the qubit pair. The setup is envisioned as a test bench for sensing gravitational effects on the quantum dynamics of gravitationally coupled massive objects. As a concrete example, we illustrate the implementation of our protocol using the current circuit QED architectures.

###### pacs:

03.65.Ud, 42.50.Pq, 85.40.Xx, 03.65.Yz## I Introduction

Experimental creations of nonclassical states of macroscopic continuous variable systems have been motivated from different perspectives. long-term motivations are to explore the limits of the standard quantum mechanics and perhaps observe possible corrections at scales where collapse phenomena Arndt and Hornberger (2014); Chen (2013); Asadian *et al.* (2014); Gittsovich *et al.* (2015) or gravitational effects become significant to fully account for the quantum dynamics in table-top quantum optics experiments Pfister *et al.* (2015).
A particularly interesting case of nonclassical states are entangled states between different modes Palomaki *et al.* (2013).
There are several proposals for entangling two mechanical oscillators either in a system interacting with a common field Mancini *et al.* (2002); Xue *et al.* (2007); Hartmann and Plenio (2008) or two distant resonators without direct interaction Pirandola *et al.* (2006); Abdi *et al.* (2012). Recently, a scheme for generating Einstein-Podolsky-Rosen entangled state of two mechanical systems has been proposed Schnabel (2015).
These schemes mainly deal with creation of bipartite *Gaussian* entanglement demonstrated by the second moments of the phase space quadratures.

Heralded entanglement generation is tailored to entangle two remote particles where the entanglement is produced conditionally based on measurement outcomes without any direct interactions between the particles Barrett and Kok (2005).
This method has been proposed for entangling distant qubits and recently has been implemented in a number of different systems Barz *et al.* (2010); Usmani *et al.* (2012); Bernien *et al.* (2013); Roy *et al.* (2015); Hofmann *et al.* (2012).
It would be intriguing to extend such technique for creating remote entanglement between macroscopic mechanical resonators.
In this work, we propose a protocol for generating entangled states of two well-separated noninteracting macroscopic mechanical resonators. Our heralding technique is useful for preparing entangled coherent states Sanders (1992), i.e. ,
between massive objects (see, Fig. 1) indicating entanglement between first moments (centers of mass) of the resonators.
This type of continuous variable entanglement is non-Gaussian which is—unlike Gaussian entangled states—characterized by negative Wigner function in phase space representation, and thus is fundamentally inconsistent with any classical description.
An interesting scenario here could be to consider the entanglement between two test masses which are only gravitationally interacting with each other and to track a genuine gravitational decoherence on the entanglement dynamics Kafri *et al.* (2014). Furthermore, the setup is useful for Bell test performed on remote macroscopic subsystems Hofer *et al.* (2015); Vivoli *et al.* (2015); Vlastakis *et al.* (2015); Pfister *et al.* (2015).
We show that the proposed protocol is experimentally feasible with the current technology of the so-called circuit QED devices Wallraff *et al.* (2004).

## Ii Model

We use a particular hybrid quantum device for implementing our entangling protocol.
The hybrid device is composed of a mechanical oscillator with frequency , a photonic cavity mode with frequency , and a qubit system made up of ground and excited energy levels.
We also need an extra, auxiliary, excited state for conditionally entangling the qubit with an emitted cavity photon.
Solid-state qubits such as spin of a nitrogen-vacancy center Iakoubovskii *et al.* (2000) or superconducting transmon qubits Houck *et al.* (2007) have such a structure.
The qubit strongly interacts with both the mechanical and the electromagnetic modes, i.e., its coupling rates to these modes are greater than the respective decoherence and damping rates of the system.
The mechanical resonator is pre-cooled to its motional ground state which can be achieved by sideband cooling of the mechanical resonators via its coupling to the qubit Rabl (2010).
Since typically the thermal excitation numbers of the cavity and qubit itself are very small at cryogenic ambient temperatures, this way one could practically prepare the system very close to its ground state.
In order to create a mechanical entangled coherent state we are specifically interested in a qubit–mechanics interaction of the form with and being position of the resonator. This describes a state-dependent force generating a time evolution of state-dependent displacement of the resonator.
In this model the cavity mode resonantly interacts with the qubit described by the Jaynes-Cummings Hamiltonian , where is the cavity field annihilation operator.
This leads to oscillatory transition which can be exploited to entangle the qubit with a traveling photon.

## Iii Protocol

Having the above system in two remote sites, the following protocol can be applied for entangling their distant mechanical resonators:

(1) Initialization. The separated parties and prepare their system in the ground state: , where within each ‘ket’ the first argument stands for the qubit, the second for the mechanical mode being in ground state, and the third argument denotes the vacuum mode of the cavity. Immediately afterwards a pulse is applied to each qubit giving where .

(2) Displacing the resonators.
After the state initialization, the qubit and its respective mechanical mode interact for a time duration of .
For a coupling rate the maximum achievable displacement is which is obtained when .
Therefore, the mechanical modes will be conditionally displaced in their phase space.
The distance determines the distinguishability of the two mechanical coherent states.
Thus, one needs to increase it by enhancing the qubit–mechanical coupling rate which is actually very demanding from a technical point of view though notable efforts for attaining it Pirkkalainen *et al.* (2015).
Instead, the long coherence time of the current technology of qubits allows for exploiting the method proposed in Ref. Tian (2005) to increase this distance.
The idea is to apply a sequence of pulses to flip the qubit periodically synchronized with the resonator frequency.
By choosing an odd number of such pulses at half mechanical period time intervals the state prepared at the end of this stage will be
where with the number of pulses.

(3) Conditional photon emission. We now ‘copy’ the qubit excitations into the cavity photons by employing their third state . One first applies a pulse to bring up the qubit from to . Once being back in , due to relaxation or other mechanism, a single photon will be emitted conditioned on the state of the qubit prior to the pulse. This, therefore, results in the following state at each party

(1) |

already indicating a local three-body entangled state between, the qubit, mechanical resonator and cavity photon.

(4) Heralded hybrid entanglement. The cavity modes of the parties and interfere on a 50:50 beam-splitter resulting in

where denotes the state of the output modes after the action of the beam-splitter.
An essential feature appearing here is Hong-Ou-Mandel effect in the second line. This effect has been recently realized in microwave regime Lang *et al.* (2013) This arises due to the indistinguishability of the two input mode photons interfering on the beam-splitter. Therefore, generation of indistinguishable photons at the input ports is crucial for successful heralded entanglement. The first two lines show no entanglement while the third and forth lines contain hybrid entanglements.
By inspecting the above expression, we realize that the projection onto an entangled state in a single shot measurement can be achieved by adopting a suitable photon detection scheme.
Here, we choose it to be the photon-number parity detection at the beam-splitter output modes. This can be done by placing a ‘detector’ qubit at each output port and performing a controlled-phase gate causing a bit flip in the corresponding qubit conditioned on the arrival of odd number (here, only one photon) of photons.
Therefore, we can map the photon parity of the output modes onto the qubit states and serve as parity detectors.
The bit flip can be detected in a Ramsey measurement by inserting between two pulses on the detector qubits. Accordingly,
detection of a bit flip in either output port heralds the corresponding qubit-mechanical entangled state

(2) |

The above state involves superposition of the two resonators being in different relative distances (see, Fig. 1a). Another form of entangled state can be created by simply overturning the direction of the state-dependent force in one of the qubit devices giving superposition of the center-of-mass in different locations (see, Fig. 1b).

(5) Entangled cat state. Finally, in order to disentangle the mechanics from the qubits and obtain a purely mechanical entanglement, each site applies a pulse to its qubit resulting in:

Here, we have defined the state for the mechanical parties which already exhibits entanglement between the first moments of the two distant mechanical resonators. The users now can read out their qubit and post-select with probability according to the measurement outcomes.

## Iv Implementation

Here, we focus on a specific implementation of our protocol.
High controllability and tunability of superconducting qubits makes them a versatile tool for engineering different regimes of interactions and control over photonic and mechanical systems Xiang *et al.* (2013); Kurizki *et al.* (2015).
There are important achievements in fabricating such hybrid devices where a vibrational mode of a mechanical resonator is coupled to a superconducting qubit and at the same time the qubit strongly interacts with a coplanar microwave resonator Pirkkalainen *et al.* (2015); Lecocq *et al.* (2015).
These make the hybrid circuit quantum electrodynamical devices a promising framework for implementing our protocol. We consider, in particular, a recently proposed hybrid electromechanical circuit which can basically be fabricated and employed by the current technology Abdi *et al.* (2015).

The device consists of a transmon qubit capacitively coupled to a microwave coplanar waveguide and a mechanical resonator (see Fig. 1).
The reduced anharmonicity of a transmon qubit makes it possible to access its higher levels, thus giving us the three-level ladder system for producing cavity photons from the qubit excitations.
The Hamiltonian of the system at site A (similarly at site B) reads Koch *et al.* (2007); Abdi *et al.* (2015)

(3) | |||||

where, and are Josephson and charging energy of the superconducting qubit, respectively. Here, is the annihilation operator of the microwave photons inside cavity , while and are respectively the normalized mechanical position and momentum operators with commutation relation . Also, and are the superconducting charge number and phase operators, satisfying the commutation relation and is the Lamb-Dicke parameter. Since the charging energy of the qubit depends on the position of the mechanical resonator, one arrives at for the transmon–mechanical mode interaction featuring state-dependent force on the mechanical mode. Here, where is the number of exchanged Cooper pairs with corresponding eigenstate . By writing this in energy eigenstates of the transmon Hamiltonian one retrieves the interaction given above. The influence of the number of transferred cooper pairs by the cavity field also results in a Jaynes-Cummings interaction between the qubit transitions and the cavity mode with Rabi frequency [see Appendix A]. This can be employed for coherently converting the qubit excitations into the cavity photons as explained below.

In such devices, an external magnetic field applied to the Cooper pair box tunes the transition frequencies of the qubit.
Thus, bringing either of the transmon transitions into resonance or taking them off-resonance from the cavity mode frequency.
In the third step of the protocol, we specifically are interested in the situation where frequency of the cavity matches qubit’s first to second excited state transition: .
To make a travelling microwave photon conditioned on the state of the qubit, one first applies a -pulse which flips the qubit from to then takes the state transfer interaction between transmon and cavity into resonance for half period of a Rabi oscillation .
This brings the qubit into its first excited state accompanied with emission of a microwave photon.
We remind that high fidelity single-qubit gate operations can be performed by properly shaped dispersive microwave pulses which can be used for controlled flipping of the qubit in either of its transitions Motzoi *et al.* (2009); Chow *et al.* (2010).

Emitting indistinguishable photons from each cavity is essential for a faithful projection onto an entangled state. Indistinguishability of the photons can indeed be guaranteed by employing two low finesse microwave resonators on the sites. This leads to a broad wave-packet which increases overlap of the photonic wavepackets, and therefore, their indistinguishability. Its other consequence is lowering emission time of the photons which basically makes them travelling photons. Therefore, reducing the protocol run time.

Moreover, the controlled phase gate performed on detector qubits discussed in the fourth step of the protocol can also be realized using superconducting qubits dispersively coupled to a coplanar microwave resonator at the output of the beam-splitter. The Hamiltonian of this part of the system is that of a qubit dispersively coupled to a cavity with frequency , the same frequency of the site cavities (see below)

(4) |

where is the qubit-cavity detuning.
Because of this dispersive coupling the unitary time evolution operator of the detector qubits in a frame rotating at the qubit and cavity frequencies is .
Our goal is to have a conditional -phase shift given that the photon number in the cavity is odd.
Therefore, the interaction time of the qubit and the photons must satisfy .
By storing the incident photons in a high finesse cavity the interaction time between photons and the detector qubits can be increased sufficiently large.
Therefore, after arrival of the photons to these secondary cavities one waits for seconds then measures the parity by performing a Ramsey pulse sequence on detector qubits Sun *et al.* (2014).
Today’s experiments are able to measure the parity of the storage microwave cavity with fidelities above 90% via another read out cavity coupled to a transmon superconducting qubit Vlastakis *et al.* (2015).

It worth mentioning here that to maximize the absorption of these photons (after mixing at the beam-splitter) into the detecting cavities two conditions must be met: First, the resonance frequency of the cavities must match, their linewidth must be close to each other.
The former condition, in principle, is easily met by employing equal frequency cavities, while the latter sounds contradictory.
On the one hand, in the detection parts, we need to employ high finesse cavities in order to store the arrived photons and give them enough time to rotate the qubit states.
On the other hand, as discussed above, site cavities must have high decay rates to ensure indistinguishability of the outgoing photons.
This can be resolved by employing a cavity with tunable decay rate in the detecting parts Pierre *et al.* (2014).
These cavities are designed such that can be tuned in situ to different decay rates for the purpose of maximal capture, storing, and retrieving microwave photons.
In our case, the detecting cavities can be first tuned to a linewidth which allows for maximal capture of the incoming photon(s), then the photon(s) can be stored for performing the parity measurement.

## V Discussion

The question of how to verify the nonlocal coherence involving the mechanical parts can be approached by looking into the dynamics of the two qubits coupled to their respective resonators.
Since in Eq. (2) is periodic in time it becomes zero at some points at which the mechanical part is disentangled from the qubit part, and thus leaving the qubit pair in a maximally entangled state. This is manifested in the entanglement revival in the qubit pair which can be considered as a signature of the entangled superposition state between the masses and qubits, i.e., Eq. (2), at intermediate times. A similar consideration has been used previously for probing macroscopic superposition states in Refs. Marshall *et al.* (2003); Armour *et al.* (2002). In our case this feature can point to even a stronger indication of nonclassicality. Monotonicity of the entanglement allows the entanglement revival if and only if the qubit pair has access to a global coherent operation (e.g. direct interaction) inducing this entanglement. Since there is no direct interaction between the qubits this revival has to be provided by the dynamics of the resonators being in an entangled superposition with the pair prior to the revival time. This is an unambiguous way of verifying entangled superposition including the mechanical parts. Because the action of separable resonators on the qubits fall under the local operation and classical communication (LOCC) operations which cannot increase the entanglement. For this purpose, one should use an entanglement monotone such as concurrence Wootters (1998) or negativity Vidal and Werner (2002) for monitoring the entanglement. Concurrence has already been measured experimentally on superconducting qubits with high fidelity using state tomography DiCarlo *et al.* (2009); Ristè *et al.* (2013).
In Fig. 2 we plot time evolution of the concurrence of the two qubits for two different displacement amplitudes and the effect of decoherence on this evilution.
The detailed analytical and numerical analyses are presented in Appendix B.

Information transmission about mass’ positions due to, for instance, the interaction with the environment suppresses the magnitude of the entanglement revival and the state evolves towards statistical mixture with the same rate regardless of which form of mechanical entanglement (a) or (b) shown in figure 1 was created [see Appendix C]. However, if we think of an unconventional situation in which the two masses gravitationally interact with each other, then the masses’ configuration in the mechanical parts of the entanglement becomes important. One model of gravitational decoherence is presented in Ref. Diosi (2011) which is equivalent to a different formulation proposed in Ref. Kafri *et al.* (2014). In these models the decoherence rate is completely determined by the gradient
of the gravitational force between the masses. The gradient of the gravitational forces in the two types of entanglement shown in Fig. 1(a) and (b) are different. Therefore, adding this new source of noise may lead to a detectable gap in the decay of entanglement between the two cases undergoing identical environmental decoherence, yet different gravitational decoherence rates. This gap then can be attributed to a genuine gravitational effect. This is an intriguing strategy. Because, the main challenge in observation of gravitational effects is that it is hard to distinguish the intrinsic gravitational field effects from those of environmental (conventional) decoherence in quantum dynamics as both these sources of noise leads to similar reduction of the superposition states.

## Vi Conclusion

In this work we proposed an experimentally feasible protocol for creating heralded entangled cat states between spatially separated mechanical resonators. The scheme can be implemented in currently available circuit QED architectures. We expect that the setup will provide a useful platform for experimentally probing the interface between gravity and quantum physics where the mechanical resonators in the present scenario serve as test masses undergoing gravitational decoherence. Finally, we think that going beyond a single particle superposition, and monitoring the dynamics of ‘gravitationally’ different types of entangled states may provide a more visible test of gravitational decoherence models in a new regime of gravitational sensing.

###### Acknowledgements.

The authors are grateful of O. Gühne and G. J. Milburn for useful comments. A.A. acknowledges the support by Erwin Schrödinger Stipendium No. J3653-N27 and M.A. by the Alexander von Humboldt Foundation via a postdoctoral fellowship.## Appendix A Implementation Hamiltonian

The Hamiltonian describing the device we are considering in this paper is composed of the transmon qubit coupled to the mechanical resonator via its charging energy, that is ()

(5) | |||||

where is the circuit QED Lamb-Dicke parameter and is the induced dc gate charge. Here, is the externally applied magnetic flux through the superconducting loop and is the superconducting flux quantum. The operators and denote the number of Cooper pairs transferred between the islands and the gauge-invariant phase difference between the superconductors, respectively.

Since in transmon qubits , the nonlinearity of the qubit is reduced such that the higher levels play a role in its dynamics. We simplify the above Hamiltonian by first Taylor expanding the and keeping only to the first order in . Then applying rotating wave approximations and truncating the transmon Hilbert space to its first three levels. This brings us to

(6) | |||||

with . Here, is the annihilation operator of the microwave photons inside the cavity, while and are respectively the normalized mechanical position and momentum operators with commutation relation . The Hamiltonian (6), already features a state-dependent force on the mechanics via the qubit which is crucial in our protocol. The second line of the Hamiltonian is the generalized Jaynes-Cummings transmon-cavity interaction.

The probe qubits discussed in the fourth step of the protocol could be realized by superconducting qubits dispersively coupled to the coplanar microwave transmission line at the output of the beam-splitter. Here, the goal is to have a parity flip in the qubits conditioned on the odd incident photon numbers. Therefore, the interaction time of the qubit and the photons must satisfy . The pulse duration of the incident single-photon states is roughly the same as the time it has taken to be emitted in the sites, i.e. . Therefore, one needs to fulfil for getting a half rotation about the -axis for every photon. Since we have from the dispersive coupling regime, this can be achieved only for . This looks, however, impractical for a waveguide because of the finite coherence time of the qubits. By storing the incident photons in a high finesse superconducting resonator the interaction time between photons and the probe qubits can be significantly increased.

## Appendix B Monitoring the entanglement dynamics under the influence of decoherence

We aim to monitor the entanglement dynamics of the two-qubit system in which the qubits and their respective resonator evolve under the influence of the interaction with the environment. Here, we use the Wooter’s concurrence, which is defined as where are the eigenvalues of the matrix and is the maximum eignevalue.

### b.1 Environmental decoherence: Analytic treatment

As an illustration, in this section we treat a simple scenario of probing the entanglement dynamics under decoherence in which the qubits coupled with constant coupling strengths to their respective mechanical resonators. We evaluate the effect of qubit decoherence and mechanical decoherence due to a weak coupling of the resonator to a finite temperature bath where the Markovian master equation of a Lindblad form is applied. Therefore, we model these decoherence processes by a master equation of the form

(7) |

where

(8) |

describes qubit dephasing acting locally with dephasing time and

(9) |

describes the mechanical dissipation within a single resonator, where is the mechanical damping rate for mechanical resonator with quality factor and is the equilibrium occupation number (identical resonators and dissipation is assumed). In the high temperature limit we obtain as the relevant mechanical decoherence rate.

Let us now consider the effect of decoherence on the entanglement dynamics of the qubit pair after the preparation of the hybrid entangled state obtained by the heralded detection. The state at later time after a period of obtained from Liouville superoperator of the time evolution which is the solution to Eq. (7). In particular we are interested in the dynamics of the off-diagonal term of the qubit pair’s reduced state

Now, let us analyze a scenario in which after a fast preparation of the hybrid entangled state at time the system undergos decoherence processes according to Eq.(7) at the later time () For the expectation value at we need to solve

(10) |

According to the master equation (7) and Hamiltonian (1) in the main text this operator evolves as

(11) |

We have assumed for both of the systems.

From above we can obtain the time evolution of an entanglement monotone for the two-qubit system. In the special case we consider, concurrence simplifies to

(12) |

Note that, the concurrence is completely characterized by the off-diagonal term of the qubit pair’s state in which the qubits undergo only dephasing and the mechanical resonators interact with a finite temperature bath. The plot is shown in Fig. 2(a). In the absence of the resonator decoherence

The displacement amplitude is which is periodic in time, and thus the concurrence demonstrates collapse and revival. Here, we have taken . Concurrence is the same regardless of which state

(13) |

or

(14) |

was prepared at time . For entangled state (14) the exerted force from the second qubit is reverted with respect to the -axis. Therefore, the associate configuration indicate superposition between two different relative distance between the resonators.

Equivalently, we can define the Wigner characteristic function and and write

(15) |

The evolution of the characteristic function is given by the Fokker-Planck equation

(16) | ||||

where . We solve this equation in three steps. First, we make the ansatz

(17) |

where , and . For the remaining equation for we have

(18) |

We now make the second ansatz

(19) |

where , and . This leaves us with the remaining equation for , which is given by

(20) |

This equation is solved by a function of the form

(21) |

and the specific expression for is determined by the initial conditions

(22) |

Ã For an initial thermal state and therefore

(23) |

In the limit where the dissipative part of the master equation is valid we obtain

(24) |

This shows that for the signal of a single measurement decays with a total decoherence rate . For we obtain . Similar conclusions are obtain, when starting from a precooled state or for , when a -pulse sequence is obtained to amplify the displacement amplitude.

We should remark that in the decoherence model it is assumed that the oscillator damping is very small (). We use this model of damping in this section to provide a simple illustrative evaluation of the effect of the resonator’s dissipation on the qubit pair’s entanglement dynamics. In the next section we use quantum Brownian motion master equation giving more accurate results in higher damping rate and reduce to above master equation using rotating wave approximation. In the next section, we treat this scenario numerically in a more involved decoherence processes.

### b.2 Numerical simulation of the full dynamics of the system

The system dynamics is composed of free evolution and the dissipations. The free evolution part is captured by the system Hamiltonians in the mechanical sites given by (5) and (3) in the main text, and the dispersive Hamiltonian describing the dynamics of the probe sites given by the Eq. (4) in the main text. And the dissipations and decoherence stemming from dephasing of the transmon qubits and relaxations in the cavity modes, and finally thermalization in the mechanical modes. Here, we assume that the cavity decay rates is majorly due to its coupling the the transmission lines. This will not lead to serious restriction as it is almost the practical case. Moreover, note that the cavities are irreversibly coupled to each other via transmission lines mediated by a beam-splitter. These are the major sources of imperfection and one includes them all in a single master equation to study full dynamics of the whole system.

(25) |

where is the total Hamiltonian which is composed of the both mechanical sites and the detection facilities denoted by ‘p’ subscript. We have also decomposed the dissipator into three part each of which for a distinct form of energy: photonic, charge, and phonoic. The mechanical resonators damping rate is and because of the low frequency its diffusion is largely affected by the environmental thermal phonons . The mechanical dissipator is already introduced in Eq. (9).

For the qubits, the dissipator is can be divided into the four local dissipators (for four qubits) . This assumption is true so long as there is no dissipative coupling between the qubits, and keeps to be true in our case, since the qubits are well separated from each other and any potential coupling between them will happen coherently via the coupled cavities (see below). The dissipator of every qubit must include both relaxation and its pure dephasing. For example is given by

(26) | |||||

where with are the three lowest transmon states. The first of the above dissipator and the second line correspond respectively to the relaxation and the pure dephasing of qubit . The relaxation happens with rate and the pure dephasing rate is such that the total dephasing time of the qubit is . The same arguments hold for the remaining three qubits.

Finally, for the cavities, the photons experience both decay into the coupled transmission line and loss (absorption and decay into the free space), leading to the total decay rate . However, the current technology superconducting resonators have negligible loss and the microwave photons mostly escape to the coupled transmission line . We therefore, have a dissipator in the Lindblad form with the decay rates and for the site and probe cavities, respectively. We also include the unidirectional coupling of the detecting cavities and the site cavities. Actually, this coupling is not direct and it happens via the beam-splitter. Thus, the detecting cavities are irreversibly fed by both of the site cavities. Since we are considering a 50:50 beam-splitter, the effective mode which would couple to the and cavities are . By taking this into account, the following Liouvillian holds for the cavity modes

(27) | |||||

where is the Lindblad dissipator. Here, the last two lines express the irreversible coupling between the probe and site cavities including the beam-splitter mixing effect and is the efficiency of the transmission channels which takes values Gardiner (1993); Carmichael (1993). The parameter contains the waveguide losses, which are typically negligible, and reflection of the microwave photons at the port of the detecting cavities and the beam-splitter. As we have discussed in the main text the reflection effects can be minimized by appropriate choice of devices, therefore, giving very close to one.

### b.3 Simulation’s parameters

One numerically solves the full master equation (25) with the above Liouvillians and by a post-selection simulated by a projective measurement one of the states (2) will be obtained. However, this cannot be done by the available computational resources because of the very large system size. Instead, we turn to simulate the protocol step by step. Therefore, we first solve the master equation for each local site coupled the waveguide at the output which at the end of the third step of the protocol gives with , then we merge their photon parties by assuming a perfect 50:50 beam-splitter. In the next step, the dynamics in the probe cavities are simulated with the initial separable quadripartite qubit-cavity state , where is the photonic party of the state mixed at the beam-splitter. Here expresses the beam-splitters unitary operation, while and are the outputs of the first stage of the simulations. The measurements are simulated as perfect projections. Therefore, the post-selected state according to the parity of the detecting qubits and state of the site qubits gives us the final state which turns out to be an entangled coherent state.

In Fig. 3 we plot the evolution of the probability of simultaneously finding the and qubits in the ground state and excited state, respectively. The parameters used here are feasible with the current technology [see Table 1]. In the detection section, a microwave resonator with kHz is considered. The detector qubits are transmon qubits with the same properties as the site qubits and operated at which makes them well away from the cavity resonance (). This means having giving enough time for changing the qubit parity. The ambient temperature is taken to be mK.

It worth to mention here that in the detection parts, we need to employ high-Q cavities in order to store the arrived photons and give them enough time to rotate the qubit states and thus perform the parity measurement.
On the other hand, the photons leaving the and cavities must have broader band to ensure indistinguishability, the key point for Hong-Ou-Mandel effect.
To maximize the absorption of these photons (after mixing at the beam-splitter) into the detecting cavities two conditions must be met: First, the resonance frequency of the cavities must match, their linewidth must be close to each other.
The former condition, in principle, is easily met by employing equal frequency cavities.
However, the second condition requires , which is in contradiction with our requirements of the parity measurement.
This can be resolved by employing a cavity with tunable decay rate in the detecting parts Pierre *et al.* (2014).
These are designed such that can be tuned in situ to different decay rates with three orders of magnitude difference.
In our case the difference between and is only one order of magnitude: s and s within the reported cavity lifetimes in Ref. Pierre *et al.* (2014).

Quantity | Symbol | Value |
---|---|---|

Mechanical mass | 3 pg | |

Mechanical frequency | 1 MHz | |

Mechanical quality factor | ||

Transmon-mechanics coupling rate | 50 kHz | |

Josephson energy | 35–55 GHz | |

Charging energy | 0.5 GHz | |

Transmon relaxation rate | 5 kHz | |

Transmon pure dephasing rate | 20 kHz | |

Transmon-cavity coupling rate | 45 MHz | |

Cavity frequency | 11 GHz | |

Cavity decay rate | 200 kHz |

## Appendix C Gravitational decoherence of collective modes of two gravitationally coupled harmonic oscillators

For the sake of gaining more insight into the idea discussed in the main text let us focus on a simple and ideal example.
An interesting scenario happens if the two mechanical resonators are gravitationally coupled to each other for sufficiently large gravitational coupling strength. In this case dynamics of the two resonators is described by two independent collective modes of oscillations with coordinates , for the center of mass mode, and , for the breathing mode, at corresponding frequency . For further detail on the model see Ref. Kafri *et al.* (2014).

The inter-mode interaction is , where is the coupling rate with the universal gravitational constant and the distance between the mechanical resonators. One then diagonalizes the total Hamiltonian by applying the unitary operator . The effect of Gravitational coupling will appear in the free evolution part of Eq. (11) and after applying the above unitary transformation the equation reads

(28) |

where are the normal annihilation operators and .

The initial conditions (the coupling signs to the qubits and the initial positions) are fixed such that only one normal mode is excited and put into a superposition via heralded technique. In realistic situation the normal mode splitting is expected to be small, and thus is approximated to be Kafri *et al.* (2014). In the first scenario we imagine that only the center-of-mass mode with frequency is excited by the proper constant coupling strength to the qubit: . This corresponds to Eq.(14) and gives

Here, . While for the second scenario where only breathing mode is excited with the proper coupling strength to the qubits: . This corresponds to Eq.. Therefore, we have

oscillating with normal mode frequency .

Therefore, in this representation we can think of a single massive system being in a harmonic potential oscillating with frequency . Environmental noise as described by (9) is coupled locally to each resonator’s coordinate, and therefore qubit pair entanglement, according to the definition (12), undergoes the same decoherence rate regardless of the collective mode. In the presence of large enough normal mode spitting , due to gravitational force between the mechanical resonators, the collective mechanical modes undergo different gravitational decoherence rates. The difference in the decoherence generates a gap in the amount of entanglement revival per single oscillation period. The effect of normal mode frequency splitting, i.e. the gap in the entanglement revival, becomes more pronounced for larger . The decoherence channel opens up by the gravitational interaction may captured by

(29) |

This form of master equation leads to the suppression of superposition in the position coordinate of a collective mode. This effect can also be easily seen by rewriting the double commutator in the position space,

describing in particular the spatial collapse of the oscillating mode due to gravitational decoherence.
The gravitational decoherence is proportional to the inverse of the normal mode frequency, i.e., Kafri *et al.* (2014). Therefore, unlike environmental decoherence gravitational decoherence attribute different decoherence rates for the two collective modes which can be manifested in the gap in the amount of the entanglement revival between the two scenarios. This is a preliminary illustration of how non-Gaussian entangled state might be useful for probing quantum dynamical phenomena which are of purely gravitational effect. This leads us to ask, can non-Gaussian mechanical entanglement provide us a sensitive detection of the genuine effect of gravitational decoherence? To estimate this effect in real experimental situation one need to take into account all the practical limitations which is a challenging issue but not impossible.

### References

- M. Arndt and K. Hornberger, Nature Phys. 10, 271 (2014).
- Y. Chen, J. Phys. B: At. Mol. Opt. Phys. 46, 104001 (2013).
- A. Asadian, C. Brukner, and P. Rabl, Phys. Rev. Lett. 112, 190402 (2014).
- O. Gittsovich, T. Moroder, A. Asadian, O. Gühne, and P. Rabl, Phys. Rev. A 91, 022114 (2015).
- C. Pfister, J. Kaniewski, M. Tomamichel, A. Mantri, R. Schmucker, N. McMahon, G. Milburn, and S. Wehner, arXiv:1503.00577 [quant-ph] (2015).
- T. A. Palomaki, J. D. Teufel, R. W. Simmonds, and K. W. Lehnert, Science 342, 710 (2013).
- S. Mancini, V. Giovannetti, D. Vitali, and P. Tombesi, Phys. Rev. Lett. 88, 120401 (2002).
- F. Xue, Y. xi Liu, C. Sun, and F. Nori, Phys. Rev. B 76, 064305 (2007).
- M. J. Hartmann and M. B. Plenio, Phys. Rev. Lett. 101, 200503 (2008).
- S. Pirandola, D. Vitali, P. Tombesi, and S. Lloyd, Phys. Rev. Lett. 97, 150403 (2006).
- M. Abdi, S. Pirandola, P. Tombesi, and D. Vitali, Phys. Rev. Lett. 109, 143601 (2012).
- R. Schnabel, Phys. Rev. A 92, 012126 (2015).
- S. D. Barrett and P. Kok, Phys. Rev. A 71, 060310(R) (2005).
- S. Barz, G. Cronenberg, A. Zeilinger, and P. Walther, Nat. Photon. 4, 553 (2010).
- I. Usmani, C. Clausen, F. Bussieres, N. Sangouard, M. Afzelius, and N. Gisin, Nat. Photon. 6, 234 (2012).
- H. Bernien, B. Hensen, W. Pfaff, G. Koolstra, M. S. Blok, L. Robledo, T. H. Taminiau, M. Markham, D. J. Twitchen, L. Childress, and R. Hanson, Nature 497, 86 (2013).
- A. Roy, L. Jiang, A. D. Stone, and M. Devoret, arXiv:1505.01178 [quant-ph] (2015).
- J. Hofmann, M. Krug, N. Ortegel, L. Gérard, M. Weber, W. Rosenfeld, and H. Weinfurter, Science 337, 72 (2012).
- B. C. Sanders, Phys. Rev. A 45, 6811 (1992).
- D. Kafri, J. M. Taylor, and G. J. Milburn, New J. Phys. 16, 065020 (2014).
- S. G. Hofer, K. W. Lehnert, and K. Hammerer, arXiv:1506.08097 [quant-ph] (2015).
- V. C. Vivoli, T. Barnea, C. Galland, and N. Sangouard, arXiv:1506.06116 [quant-ph] (2015).
- B. Vlastakis, A. Petrenko, N. Ofek, L. Sun, Z. Leghtas, K. Sliwa, Y. Liu, M. Hatridge, J. Blumoff, L. Frunzio, M. Mirrahimi, L. Jiang, M. H. Devoret, and R. J. Schoelkopf, arXiv:1504.02512 [quant-ph] (2015).
- A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Nature 431, 162 (2004).
- K. Iakoubovskii, G. J. Adriaenssens, and M. Nesladek, J. Phys.: Condens. Matter 12, 189 (2000).
- A. A. Houck, D. I. Schuster, J. M. Gambetta, J. A. Schreier, B. R. Johnson, J. M. Chow, L. Frunzio, J. Majer, M. H. Devoret, S. M. Girvin, and R. J. Schoelkopf, Nature 449, 328 (2007).
- P. Rabl, Phys. Rev. B 82, 165320 (2010).
- J.-M. Pirkkalainen, S. U. Cho, F. Massel, J. Tuorila, T. T. Heikkilä, P. J. Hakonen, and M. A. Sillanpää, Nat. Commun. 6, 6981 (2015).
- L. Tian, Phys. Rev. B 72, 195411 (2005).
- C. Lang, C. Eichler, L. Steffen, J. M. Fink, M. J. Woolley, A. Blais, and A. Wallraff, Nature Phys. 9, 345 (2013).
- Z.-L. Xiang, S. Ashhab, J. Q. You, and F. Nori, Rev. Mod. Phys. 85, 623 (2013).
- G. Kurizki, P. Bertet, Y. Kubo, K. Mølmer, D. Petrosyan, P. Rabl, and J. Schmiedmayer, Proc. Natl. Acad. Sci. 112, 3866 (2015).
- F. Lecocq, J. D. Teufel, J. Aumentado, and R. W. Simmonds, Nature Phys. 11, 635 (2015).
- M. Abdi, M. Pernpeintner, R. Gross, H. Huebl, and M. J. Hartmann, Phys. Rev. Lett. 114, 173602 (2015).
- J. Koch, T. Yu, J. Gambetta, A. Houck, D. Schuster, J. Majer, A. Blais, M. Devoret, S. Girvin, and R. Schoelkopf, Phys. Rev. A 76, 042319 (2007).
- F. Motzoi, J. M. Gambetta, P. Rebentrost, and F. K. Wilhelm, Phys. Rev. Lett. 103, 110501 (2009).
- J. M. Chow, L. DiCarlo, J. M. Gambetta, F. Motzoi, L. Frunzio, S. M. Girvin, and R. J. Schoelkopf, Phys. Rev. A 82, 040305 (2010).
- L. Sun, A. Petrenko, Z. Leghtas, B. Vlastakis, G. Kirchmair, K. M. Sliwa, A. Narla, M. Hatridge, S. Shankar, J. Blumoff, L. Frunzio, M. Mirrahimi, M. H. Devoret, and R. J. Schoelkopf, Nature 511, 444 (2014).
- M. Pierre, I. M. Svensson, S. Raman Sathyamoorthy, G. Johansson, and P. Delsing, App. Phys. Lett. 104, 232604 (2014).
- W. Marshall, C. Simon, R. Penrose, and D. Bouwmeester, Phys. Rev. Lett. 91, 130401 (2003).
- A. D. Armour, M. P. Blencowe, and K. C. Schwab, Phys. Rev. Lett. 88, 148301 (2002).
- W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998).
- G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 (2002).
- L. DiCarlo, J. M. Chow, J. M. Gambetta, L. S. Bishop, B. R. Johnson, D. I. Schuster, A. B. J. Majer, L. Frunzio, S. M. Girvin, and R. J. Schoelkopf, Nature 460, 240 (2009).
- D. Ristè, M. Dukalski, C. A. Watson, G. de Lange, M. J. Tiggelman, Y. M. Blanter, K. W. Lehnert, R. N. Schouten, and L. DiCarlo, Nature 502, 350 (2013).
- L. Diosi, J. Phys.: Conf. Ser 306, 012006 (2011).
- C. W. Gardiner, Phys. Rev. Lett. 70, 2269 (1993).
- H. J. Carmichael, Phys. Rev. Lett. 70, 2273 (1993).