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

Indefinite temporal order without gravity

Kacper Dębski\scalerel* | [email protected] Institute of Theoretical Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland    Magdalena Zych [email protected]    Fabio Costa [email protected] Australian Research Council Centre for Engineered Quantum Systems, School of Mathematics and Physics, The University of Queensland, St Lucia, QLD 4072, Australia    Andrzej Dragan\scalerel* | [email protected] Institute of Theoretical Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, 117543 Singapore, Singapore
Abstract

According to the general theory of relativity, time can flow at different rates depending on the configuration of massive objects, affecting the temporal order of events. Recent research has shown that, combined with quantum theory, this gravitational effect can result in events with an indefinite temporal order, which can be tested through the violation of Bell-type inequalities. According to Einstein, we shall assume physical equivalence of a uniform gravitational field and a corresponding acceleration of a reference system. Here we construct a non-gravitational scenario where accelerating particles interacting with optical cavities result in a violation of the temporal Bell inequalities analogous to the gravitational case. However, we find that the inequalities can also be violated by time-like events, exposing an ambiguity in their use as a theory-independent test of indefinite temporal order.

Introduction

Quantum theory and general relativity rest on very different foundations. In quantum theory, systems do not posses definite physical properties independent of their measurement; however, processes take place on a fixed background spacetime, where the causal relations between events—space-like or time-like—are defined independently of any operation or physical process. On the other hand, in general relativity, causal relations are not fixed a priori, as the geometry depends on the configuration of mass-energy. Therefore, it is expected that a unification of the two theories should result in non-classical, or indefinite, causal structure Butterfield2001 ; hardy2007towards .

However, the physical meaning of such an indefinite causal structure is not clear. Most approaches to quantum gravity attempt to establish a complete theoretical framework oriti2009 ; kiefer2012quantum , but do not offer a direct physical interpretation of non-classical causal relations.

Recent work has proposed a Gedankenexperiment to directly pinpoint the physical meaning of the non-classicality of temporal order, independently of the details of a full quantum theory of gravity Bell4time . The idea is that a superposition of mass configurations will produce a corresponding non-classical time dilation on any system one might use as “clock” to identify spacetime events, for example, producing effectively a superposition of orderings of time-like events. By choosing an appropriate protocol, one can perform a task—the violation of a Bell inequality—that would be impossible if the events were ordered according to a classical causal structure. This would provide a theory independent certification of the non-classicality of causal structure emerging from the incorporation of the superposition principle into general relativity.

As the proposal in Ref. Bell4time relies on time dilation, it is legitimate to ask whether a similar experiment would be possible in a purely special-relativistic setting, where the back-action on spacetime of the involved matter is negligible but where a superposition of the clock states of motion causes the corresponding non-classical time dilation. Furthermore, a flat spacetime version of the argument would allow in principle a complete/fundamental description of a potential practical implementation and thus provide insights into the requirements entering the original, gravitational, argument. Special relativistic implementation of the experiment might also be less demanding and provide a possible outlook for feasible laboratory implementations.

Here we present such an in principle complete scheme where accelerating particles interact with quantum fields according to their own internal clock degrees of freedom (DoFs). By letting the particles evolve in superposition along appropriately accelerating trajectories, the particle-field interaction events reproduce the “entanglement of temporal order” found in the gravitational case. Further measurements on the particles and fields can thus violate the “Bell inequalities for temporal order” presented in Ref. Bell4time .

Surprisingly, we find that a violation of the Bell inequalities persists when the entangled events are space-like, challenging the interpretation that the protocol demonstrates a non-classical temporal order of time-like events. We find that this is due to the failure of one of the assumptions behind the original protocol—namely that the superposed alternatives only differ in the event order, while all local evolutions are trivial. We argue that the failure of this assumption is ubiquitous and would be present in a generic dynamical context, including a gravitational implementation of the protocol. We also propose an interpretation of the Bell inequality violation in scenarios where non-classicality of temporal order cannot possibly explain the results (e.g. for spacelike separated events, as mentioned above).

Our results raise a fundamental question, whether it is possible to formulate an operational scenario that can unambiguously discriminate the non-classicality of the causal structure of spacetime from dynamical effects, necessarily present in its physical implementations, or other laboratory implementations of quantum causal structures Procopio2015 ; Rubinoe1602589 ; rubino2017experimental ; Goswami2018 ; guo2018experimental ; Wei2019 ; taddei2020experimental . We present a discussion on what extensions to the original protocol might be required in order to achieve that.

Throughout this work we use natural units, for which we have =c=1\hbar=c=1.

.1 Gravitational implementation of an indefinite causal order

For later reference. we will briefly review the salient aspects of the protocol introduced in Ref. Bell4time . The goal of the protocol is to realise four events, A1A_{1}, B1B_{1}, A2A_{2}, B2B_{2}, whose pair-wise order is ‘entangled’: A1A_{1} is in the causal past/past lightcone of B1B_{1}, denoted A1B1A_{1}\prec B_{1}, when A2B2A_{2}\prec B_{2}, and A1A_{1} is in the causal future/future lightcone of B1B_{1}, denoted A1B1A_{1}\succ B_{1} when A2B2A_{2}\succ B_{2}; the full scenario is arranged such that this entanglement leads to correlations that do not admit a classical explanation, which is formalised in analogy to Bell-like scenarios for local classical properties. Crucially, an “event” is here understood operationally as “something that happens at a particular time and place” and thus is defined by some physical reference system. In the protocol, the reference systems are four clocks, a1a_{1}, b1b_{1}, a2a_{2}, b2b_{2}, while the events are associated with quantum operations performed on some additional system when the corresponding clock reaches a specified proper time: A1A_{1} takes place (an operation is performed) when a1a_{1} reaches time τ\tau^{*}, B1B_{1} takes place when b1b_{1} reaches time τ\tau^{*}, and so on.

In flat spacetime, if the clocks are initially synchronised (in an arbitrary chosen reference frame), all events are space-like separated. However, introducing a massive body closer to some clocks than others causes a differential time dilation, which can ‘push’ some events into the future light-cone of other events. Thereby, the position of the mass provides a control of the time order of events, see Fig. 1.

Refer to caption
Figure 1: Position of a mass as a control of time order. Two identical clocks, a,ba,b are synchronised (with any source-mass sufficiently far away). Events A,BA,B are defined as the location and fixed proper time τ\tau^{*} of the corresponding clock. In the absence of any source mass, events A,BA,B are space-like. However, if a massive object is initially (just after synchronisation) placed closer to clock aa than to bb, event AA can be in the future light cone of BB (for sufficiently large τ\tau^{*}), denoted BAB\prec A. Analogously, for the mass closer to bb, event AA can be up in the past light cone of BB, ABA\prec B.

If the control mass is prepared in a semiclassical configuration denoted KAK_{A}, event AjA_{j}, j=1,2j=1,2 is in the past of BjB_{j}, while for a different configuration, KBK_{B}, event BjB_{j} is in the past of AjA_{j}. Thus, AjA_{j} is time-like from BjB_{j}, for each mass configuration, but their order is interchanged. Moreover, for both mass configurations the pair A1A_{1}, B1B_{1} is space-like from A2A_{2}, B2B_{2}.

The full protocol features also two of the above mentioned additional systems on which the operations are performed, S1S_{1} and S2S_{2}, referred to as ’targets’. The operations at events A1A_{1} and B1B_{1} are applied only on the target system S1S_{1}, while those at A2A_{2} and B2B_{2} are applied only on S2S_{2}. The two target systems are initially in a product state, |ψS1|ψS2\ket{\psi}^{S_{1}}\ket{\psi}^{S_{2}}, and the considered operations are unitaries: U^A1\hat{U}^{A_{1}} is applied on S1S_{1} at event A1A_{1}, etc. By preparing the control mass in a superposition state, |K=12(|KA+|KB)\ket{K}=\frac{1}{\sqrt{2}}\left(\ket{K_{A}}+\ket{K_{B}}\right), one obtains the final state

|Ψfin=12(|KAU^B1U^A1|ψS1U^B2U^A2|ψS2+|KBU^A1U^B1|ψS1U^A2U^B2|ψS2).\ket{\Psi^{\textrm{fin}}}=\frac{1}{\sqrt{2}}\left(\ket{K_{A}}\hat{U}^{B_{1}}\hat{U}^{A_{1}}\ket{\psi}^{S_{1}}\hat{U}^{B_{2}}\hat{U}^{A_{2}}\ket{\psi}^{S_{2}}\right.\\ +\left.\ket{K_{B}}\hat{U}^{A_{1}}\hat{U}^{B_{1}}\ket{\psi}^{S_{1}}\hat{U}^{A_{2}}\hat{U}^{B_{2}}\ket{\psi}^{S_{2}}\right). (1)

Next, the control mass is measured in the basis |±=12(|KA±|KB)\ket{\pm}=\frac{1}{\sqrt{2}}\left(\ket{K_{A}}\pm\ket{K_{B}}\right), leaving the target system in

|Ψpost=12(|ψAS1|ψAS2±|ψBS1|ψBS2),\ket{\Psi^{\textrm{post}}}=\frac{1}{\sqrt{2}}\left(\ket{\psi_{A}}^{S_{1}}\ket{\psi_{A}}^{S_{2}}\pm\ket{\psi_{B}}^{S_{1}}\ket{\psi_{B}}^{S_{2}}\right), (2)

where |ψA=U^BU^A|ψ\ket{\psi_{A}}=\hat{U}^{B}\hat{U}^{A}\ket{\psi}, |ψB=U^AU^B|ψ\ket{\psi_{B}}=\hat{U}^{A}\hat{U}^{B}\ket{\psi}, and we are using the same unitaries in the two ‘wings’: U^A1=U^A2U^A\hat{U}^{A_{1}}=\hat{U}^{A_{2}}\equiv\hat{U}^{A}, U^B1=U^B2U^B\hat{U}^{B_{1}}=\hat{U}^{B_{2}}\equiv\hat{U}^{B}.

In general, |Ψpost\ket{\Psi^{\textrm{post}}} is an entangled state, unless |ψA=|ψB\ket{\psi_{A}}=\ket{\psi_{B}} (and it is elementary to find examples of unitaries U^A\hat{U}^{A}, U^B\hat{U}^{B} such that this is not the case). In the final step of the protocol, the entangled state is measured in appropriate bases, that lead to a violation of a Bell inequality.

The argument presented in Ref. Bell4time is that, given the initial product state of the target systems, local operations that are performed in a definite order would not be able to produce entanglement, even after conditioning on the control system. Therefore, if a set of conditions is satisfied, a violation of Bell inequalities implies that the operations were not performed in a definite order.

Another assumption of the protocol, which ends up a focal point of the present work, is that any additional evolution of the target systems (including their free evolution) between the events of interest can be neglected.

In the “Methods” and the “Results” sections, we focus on reproducing the gravitational protocol using special relativistic time dilation. We present an operational setting where the dynamics of all the relevant DoFs is incorporated, in particular the DoFs whose interactions realize the four uniatries U^A,U^B\hat{U}^{A},\hat{U}^{B}. In the “Ambiguity in the signature of indefinite temporal order” section we discuss our main result – that entanglement can be generated and Bell inequality for temporal order can be violated even if temporal order is classically defined – and argue that it is due to the failure of the additional assumption mentioned above, that target systems have trivial evolution apart from the unitaries marking the four spacetime events of interest. In the “Conclusion” section we discuss the implications of our result, including the fundamental question of how to isolate quantum features of a causal structure from other non-classical effects.

Methods

Refer to caption
Figure 2: General scheme of our protocol for a violation of Bell inequalities for temporal order. In each wing of the experiment we have a quantum field in a cavity, and two composite ’molecules’. After following entangled pairs of trajectories the molecules interact with the fields at a fixed proper time of their internal clocks. Due to time dilation, entangled state of motion gets transferred to the order of the interaction events (as well as to other degrees of freedom of the molecules and the cavities).

Instead of using a massive object to control the spacetime geometry—and thus causal order via gravitational time dilation—here we want to control the trajectories of particles so to induce special-relativistic time dilation. The protocol involves several degrees of freedom, which can be best thought of as multiple particles ‘glued’ together (up to the moment when we need to break them apart, as detailed later). We will refer to a bunch of joined particles as a ‘molecule’, although the details of what binds the particles together are irrelevant to the discussion.

Our protocol involves four molecules going through two optical cavities (two molecules per cavity). Each molecule is composed of three particles: a ‘clock’, a ‘detector’, and a ‘control’. Figure 2 presents a general scheme of our protocol.

The clock is simply a particle with some time-evolving internal state; if the molecule evolves along a classical trajectory, the internal state evolves at a rate proportional to the trajectory’s proper time. The role of the clock is to trigger an interaction between the detector and the cavity at the desired proper time. Thanks to the universality of time dilation, the protocol does not depend on the particular mechanism by which the clocks evolve—all we need is that the clock reaches two different orthogonal states depending on the two proper times involved in the protocol (see, e.g., Refs. Zych2011 ; pikovski2017time ).

The detector is a particle with two internal energy levels that (at the proper time specified by the clock), interacts with a quantum field confined in a cavity. We use the Unruh-DeWitt detector model for the interaction; see the “Unruh-DeWitt coupling” section below for the details of the coupling.

Finally, the control is a spin-12\frac{1}{2} particle, whose two orthogonal spin states serve to define the molecules’ trajectories (see, e.g., Ref. Margalit2021 for a realisation of coherent spin-dependent trajectories which could be used here). Although each detector interacts with the cavity at the same proper time of its local clock, special-relativistic time dilation implies that the interactions take place at different coordinate times depending on the molecule’s trajectory (which, in turn, depends on the spin).

In the protocol, the detectors in each molecule are prepared in their ground state and the clocks are synchronised at a reference starting time, while the spins of different molecules are prepared in an appropriate entangled state (defined below). Each molecule is sent to a cavity along a trajectory that depends on the spin. The clock triggers an interaction between the detector in the molecule and the field, creating entanglement between the two. (The proper time at which the interaction happens is chosen such that the molecule is in the cavity for either trajectory). At this point, the molecule is ‘broken apart’: the detector and the clock stay next to the corresponding cavity, while all the controls are brought together in a middle location. A joint measurement on the controls prepares an entangled state of the remaining systems, which can then be used to violate a Bell inequality.

In our setup, the ‘target’ system on each side, say S1S_{1}, comprises the field in the cavity as well as the two detectors that go through that cavity. Crucially, the field-detector interaction leaves the clock and control unaffected. This is important to ensure that each of the operations UA,UBU^{A},U^{B} acts only on the target system. Without this assumption, one could simply entangle each target system with an additional DoF, e.g. an extra particle, bring the two extra particles together and, by measuring them, induce entanglement on the two target systems. This would be a form of entanglement swapping Yurke1992 that does not require any control of time ordering (nor any control system for that matter).

Furthermore, note that since the control DoF goes through the cavity together with the rest of the molecule, we effectively need to trust the involved devices when we assume that the local operations leave the control untouched. This is the reason such a test for indefinite temporal order cannot be formulated in a device-independent way (just as in the gravitational protocol of Ref. Bell4time , see also discussion therein).

Gravity Cavity
Control system massive body spin-12\frac{1}{2} particles (one per molecule, two on each side)
Target system a single 2-level system (e.g., a spin-12\frac{1}{2} particle) Optical cavity mode and two detectors
Local operations unitaries on each system Interaction between cavity and detectors
Table 1: Comparison between the degrees of freedom involved in the gravitational scheme and ours. Here we only consider the main scheme from Ref. Bell4time (variations of the scheme involve multiple control systems or different target DoFs.)
Refer to caption
Figure 3: Trajectories of the molecules in one cavity giving rise to the required time dilation. x1x_{1}, x2x_{2} are initial positions of the molecules. The trajectory of each molecule depends on the spin of its particle called the ‘control’. For spin |\ket{\uparrow}, the trajectory is described by the proper acceleration 𝒜\mathcal{A}_{\uparrow}, and analogously for spin |\ket{\downarrow}. Dots along the trajectories divide each curve into four geometrically identical hyperbolic segments. The trajectories of the molecules in the second cavity are fully analogous.

.2 Trajectories

Each molecule has two possible trajectories, depending on the spin state. The specific trajectories we propose are shown in Fig. 3. The involved trajectories are analogous between the two wings of the experiment: They start and end at a common event but the proper time elapsing along each trajectory is different. Both trajectories can be constructed by joining four identical hyperbolic segments, characterised by proper accelerations 𝒜\mathcal{A}_{\uparrow}, 𝒜\mathcal{A}_{\downarrow} for spin up, |\ket{\uparrow}, or spin down, |\ket{\downarrow}, respectively. These segments have to be rotated or flipped according to Fig. 3, so that the acceleration for each trajectory switches sign three times. Let us assume that the acceleration for the spin |\ket{\uparrow} is bigger than an acceleration for the spin |\ket{\downarrow}. Knowing the value of proper time along the generic hyperbolic trajectory (see e.g. Ref. Dragan2021 ), we can find the difference between proper times measured at the common end of one such pair of trajectories:

Δτ=4𝒜asinh(𝒜T𝒜4)τ4𝒜asinh(𝒜T𝒜4)τ,\displaystyle\Delta\tau=\underbrace{\frac{4}{\mathcal{A}_{\downarrow}}\mathrm{asinh}(\frac{\mathcal{A}_{\downarrow}T_{\mathcal{A}}}{4})}_{\tau_{\downarrow}}-\underbrace{\frac{4}{\mathcal{A}_{\uparrow}}\mathrm{asinh}(\frac{\mathcal{A}_{\uparrow}T_{\mathcal{A}}}{4})}_{\tau_{\uparrow}}, (3)

where τ/\tau_{\uparrow/\downarrow} is the proper time measured along trajectory of a /\uparrow/\downarrow spin and T𝒜T_{\mathcal{A}} is the total travel time in a common inertial frame of reference, e.g. the frame used to initially synchronize the clocks.

Note that it is possible to achieve a large value of Δτ\Delta\tau using a small value of proper acceleration. For a large value of T𝒜T_{\mathcal{A}} we have:

Δτ=2(log(𝒜24)𝒜log(𝒜24)𝒜)\displaystyle\Delta\tau=2\left(\frac{\log(\frac{\mathcal{A}_{\downarrow}^{2}}{4})}{\mathcal{A}_{\downarrow}}-\frac{\log(\frac{\mathcal{A}_{\uparrow}^{2}}{4})}{\mathcal{A}_{\uparrow}}\right) (4)
+4(1𝒜1𝒜)log(T𝒜)+𝒪(1T𝒜2).\displaystyle+4\left(\frac{1}{\mathcal{A}_{\downarrow}}-\frac{1}{\mathcal{A}_{\uparrow}}\right)\log(T_{\mathcal{A}})+\mathcal{O}\left(\frac{1}{T_{\mathcal{A}}^{2}}\right). (5)

It means that we can achieve the required Δτ\Delta\tau (to make the events defined by such time dilated pair of clocks timelike) by accelerating for long enough time, even for arbitrarily small accelerations. Finally, taking the initial state of the spin of the two molecules to be entangled, one of the two clocks in each cavity will “be older” than another in a correspondingly correlated manner. Thus, the order of operations, controlled by the clocks, will likewise be “entagled”, as a direct consequence of the initial spin entanglement and of time dilation. In this scenario, the joint spin state of the two molecules plays the role of the control, which is played by the position of a massive object in the gravitational case Bell4time . Table 1 summarises the differences between the gravitational and our cavity-based implementation of the protocol, in terms of control systems, target(s), and local operations.

.3 Unruh-DeWitt coupling

.3.1 Interaction between a detector and a cavity

In this section we define the interaction between a two-level system, the ‘detector’, and the scalar field inside the cavity via the pointlike Unruh-DeWitt Hamiltonian Unruh1976 ; Unruh1984 . We first discuss key properties of the field operators.

We consider a scalar field of a mass mm governed by the Klein-Gordon equation111For generality, we write everything for arbitrary mm, although in the numerical calculations below we set m=0m=0.,

(+m2)ϕ=0,\left(\Box+m^{2}\right)\phi=0, (6)

in a cavity of length LL fulfilling Dirichlet boundary conditions, ϕ(x=0)=ϕ(x=L)=0\phi(x=0)=\phi(x=L)=0. The field has the following mode solutions:

un(x,t)=1ωnLsin((knx))eiωntun(x)eiωnt,u_{n}(x,t)=\frac{1}{\sqrt{\omega_{n}L}}\sin{\left(k_{n}x\right)}e^{-i\omega_{n}t}\equiv u_{n}(x)e^{-i\omega_{n}t}, (7)

where ωn=kn2+m2\omega_{n}=\sqrt{k_{n}^{2}+m^{2}}, kn=nπLk_{n}=\frac{n\pi}{L}, nn\in\mathbb{N}. Using these modes, the field operator ϕ^\hat{\phi} can be decomposed as:

ϕ^(x)=n[a^nun(x)+a^nun(x)],\hat{\phi}(x)=\sum_{n}\left[\hat{a}_{n}^{\dagger}u_{n}(x)+\hat{a}_{n}u_{n}(x)\right], (8)

where a^n\hat{a}_{n} and a^n\hat{a}_{n}^{\dagger} are annihilation and creation bosonic operators satisfying the canonical commutation relations, [a^n,a^k]=δnk\left[\hat{a}_{n},\hat{a}_{k}^{\dagger}\right]=\delta_{nk} and [a^n,a^k]=[a^n,a^k]=0\left[\hat{a}_{n},\hat{a}_{k}\right]=\left[\hat{a}_{n}^{\dagger},\hat{a}_{k}^{\dagger}\right]=0.

For the detector we consider a two-level system, the simplest model of an atom, with an energy gap Ω\Omega, and position parameter denoted xdx_{d}, (where the subscript dd hereafter stands for the ’detector’). The full Hamiltonian consists of the free Hamiltonians of the scalar field and the detector, and an interaction Hamiltonian. One of the simplest choices of the interaction between a scalar field and a two-level system is the pointlike Unruh-DeWitt (UDW) Hamiltonian which in the Schrödinger picture has the following form:

H^UDW=λχd(t)μ^Sϕ^(xd),\hat{H}_{\mathrm{UDW}}=\lambda~{}\chi_{d}(t)~{}\hat{\mu}_{\mathrm{S}}~{}\hat{\phi}(x_{d}), (9)

where λ\lambda is a dimensionless coupling constant; the real function χd(t)\chi_{d}(t) is equal to 0 when the detector does not interact and 11 for any other time and is commonly referred to as the switching function; μ^S\hat{\mu}_{\mathrm{S}} is the monopole operator μ^S=σ^++σ^=|ge|+|eg|\hat{\mu}_{\mathrm{S}}=\hat{\sigma}^{+}+\hat{\sigma}^{-}=\ket{g}\bra{e}+\ket{e}\bra{g}, where |g\ket{g} is the ground state of the two-level system and |e\ket{e} is its excited state. The Hilbert space spanned by |g\ket{g} and |e\ket{e} will be called the internal Hilbert space of the detector. Finally, ϕ^(xd)\hat{\phi}(x_{d}) is the field operator evaluated at the position of the detector.

As mentioned above, the full Hamiltonian includes also the time-independent free Hamiltonian of the field and of the two-level system, which reads H^0=nωna^na^n𝟙+𝟙Ωσ^+σ^\hat{H}_{0}=\sum_{n}\omega_{n}\hat{a}_{n}^{\dagger}\hat{a}_{n}\otimes\mathds{1}+\mathds{1}\otimes\Omega\hat{\sigma}^{+}\hat{\sigma}^{-}. Thus, the evolution of the state of the full system in the Dirac picture (also called the interaction picture) is here given by the unitary of the form:

U^=𝒯exp(idtH^UDW(D)(t)),\hat{U}=\mathcal{T}\exp{-i\int_{-\infty}^{\infty}\mathrm{d}t\hat{H}_{\mathrm{UDW}}^{(D)}(t)}, (10)

where (D)(D) stands for the Dirac picture and 𝒯\mathcal{T} is the time-ordering operator. It can be shown that RevModPhys.80.787 :

H^UDW(D)(t)=χ(t)λμ^(D)ϕ^(D),\displaystyle\hat{H}_{\mathrm{UDW}}^{(D)}(t)=\chi(t)~{}\lambda~{}\hat{\mu}^{(D)}~{}\hat{\phi}^{(D)}, (11)

where:

μ^(D)\displaystyle\hat{\mu}^{(D)} =\displaystyle= (eiΩtσ^++eiΩtσ^),\displaystyle\left(e^{i\Omega t}\hat{\sigma}^{+}+e^{-i\Omega t}\hat{\sigma}^{-}\right), (12)
ϕ^(D)(xd)\displaystyle\hat{\phi}^{(D)}(x_{d}) =\displaystyle= n(a^nun(xd)eiωnt+H.c.).\displaystyle\sum_{n}\left(\hat{a}_{n}^{\dagger}u_{n}(x_{d})e^{i\omega_{n}t}+H.c.\right). (13)

The evolution operator (10) can be expanded into the Dyson series. For small value of the coupling constant λ\lambda we limit this series only to the first order, the second order expansion does not change our results as shown in Supplementary Note 1. Thus, in the above approximation:

U^=𝟙iλdtχd(t)(eiΩtσ^++eiΩtσ^)\displaystyle\hat{U}=\mathds{1}-i\lambda\int_{-\infty}^{\infty}\mathrm{d}t~{}\chi_{d}(t)\left(e^{i\Omega t}\hat{\sigma}^{+}+e^{-i\Omega t}\hat{\sigma}^{-}\right)
×n(a^nun(xd)eiωnt+a^nun(xd)eiωnt).\displaystyle\times\sum_{n}\left(\hat{a}_{n}^{\dagger}u_{n}(x_{d})e^{i\omega_{n}t}+\hat{a}_{n}u_{n}(x_{d})e^{-i\omega_{n}t}\right). (14)

.3.2 Two molecules in the cavity

Before moving to the full protocol, let us consider the interaction between one cavity and two molecules. Each molecule contains a two-level detector, which interacts with the cavity via the UDW Hamiltonian (9), plus a clock and a control spin. At the beginning of the protocol, the cavity field is in the vacuum, denoted |0\ket{0}, and the two detectors are in their respective ground states. We synchronize the clocks, namely we prepare them in the same state |τ0\ket{\tau_{0}}, and we prepare the two control spins in the entangled state 12(|+|)\frac{1}{\sqrt{2}}(\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}+\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow}). After the spin-dependent accelerations described in the “Trajectories” section, the joint state of cavity and molecules is

12(|ττ|+|ττ|)|gg|0,\displaystyle\frac{1}{\sqrt{2}}(\ket{\tau_{\uparrow}\tau_{\downarrow}}\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}+\ket{\tau_{\downarrow}\tau_{\uparrow}}\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow})\ket{gg}\ket{0}, (15)

where |τ/\ket{\tau_{\uparrow/\downarrow}} is the state of the clock after it evolved for time τ/\tau_{\uparrow/\downarrow}. Note that for later convenience we have grouped the degrees of freedom by type (clock, spin, detector), rather than by molecule. For example, instead of

|τ||g1stmolecule|τ||g2ndmolecule,\displaystyle\underbrace{\ket{\tau_{\uparrow}}\ket{\uparrow}\ket{g}}_{1^{st}~{}molecule}\otimes\underbrace{\ket{\tau_{\downarrow}}\ket{\downarrow}\ket{g}}_{2^{nd}~{}molecule}, (16)

we write |ττ||gg\ket{\tau_{\uparrow}\tau_{\downarrow}}\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}\ket{gg}, where inside each ‘ket’ symbol we write first the state relative to molecule 1 and then the state relative to molecule 2.

We see from Eq. (15) that, after the spin-dependent accelerations, the clocks become entangled with the remaining systems. This can be problematic because, in the final protocol, we aim to observe entanglement in the target systems (detectors+cavities) after projecting the control in an appropriate state. Anything that correlates to the target, including the clocks, would effectively degrade entanglement. We can circumvent this by re-synchronising the clocks after they had passed through the cavity. This can be achieved, for example, by flipping the molecules’ spins and imparting the accelerations identical to those in the first phase, whereby at the end, all the trajectories accrue equal proper times (this is analogous to the decorrelation of clocks in the gravitational scenario, discussed in Ref. Bell4time ). This procedure allows us to ignore the clocks in the final state. Therefore, we simplify the notation and map |τ/|/|/\ket{\tau_{\uparrow/\downarrow}}\ket{\uparrow/\downarrow}\mapsto\ket{\uparrow/\downarrow} which allows us to replace state (15) with

|Ψ01=12(|+|)|gg|0,\displaystyle\ket{\Psi_{0}^{\mathrm{1}}}=\frac{1}{\sqrt{2}}(\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}+\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow})\ket{gg}\ket{0}, (17)

where the subscript 0 informs that this is the initial state of the system (before detectors and cavity interact) and the index 11 in the superscript informs us that this is the state of the field and two molecules in cavity number 11 (see Fig. 2).

Note that according to the scheme shown in the “Interaction between a detector and a cavity” section, the state |\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow} corresponds to the case when the detector at position x2x_{2} is the first to interact with the cavity (the spin of the molecule that interacts earlier is \downarrow). Similarly, for state |\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow}, the detector at x1x_{1} interacts earlier than detector at x2x_{2}. In further calculations, we assume that the duration TT of each detector’s interaction with the cavity is smaller than the time dilation between the two trajectories, i.e., TΔτT\leq\Delta\tau, to ensure that the two detectors do not interact simultaneously from the perspective of the cavity’s reference frame.

As the interaction between each detector and the cavity is described by Eq. (10) and knowing that the order of interactions is given by the molecule’s spin, we can write the final state in the form:

|Ψ1=12(|U^1U^2|g|g|0+|U^2U^1|g|g|0).\displaystyle\ket{\Psi^{\mathrm{1}}}=\frac{1}{\sqrt{2}}\left(\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}\hat{U}_{1}\hat{U}_{2}\ket{g}\ket{g}\ket{0}+\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow}\hat{U}_{2}\hat{U}_{1}\ket{g}\ket{g}\ket{0}\right). (18)

Here, operators U^1\hat{U}_{1} and U^2\hat{U}_{2} act on the internal states of the first (left) and the second (right) detector respectively. For ease of notation we define the state |ψR:=U^1U^2|g|g|0\ket{\psi_{R}}:=\hat{U}_{1}\hat{U}_{2}\ket{g}\ket{g}\ket{0}, where the subscript RR denotes that the right detector interacts before the left one, and analogously define |ψL:=U^2U^1|g|g|0\ket{\psi_{L}}:=\hat{U}_{2}\hat{U}_{1}\ket{g}\ket{g}\ket{0}. We thus have:

|Ψ1=12||ψR+12||ψL.\displaystyle\ket{\Psi^{\mathrm{1}}}=\frac{1}{\sqrt{2}}\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}\ket{\psi_{R}}+\frac{1}{\sqrt{2}}\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow}\ket{\psi_{L}}. (19)

Using the evolution operator (10) we find an explicit form of |ψR\ket{\psi_{R}}, |ψL\ket{\psi_{L}}, up to leading order in the interaction parameter λ\lambda using Eq. (14) (see Supplementary Note 2 for details). The results of this calculation are

|ψR=|gg|0+|ge|ϕgeR+|eg|ϕegR+𝒪(λ2),\displaystyle\ket{\psi_{R}}=\ket{gg}\ket{0}+\ket{ge}\ket{\phi_{ge}^{R}}+\ket{eg}\ket{\phi_{eg}^{R}}+\mathcal{O}(\lambda^{2}), (20)
|ψL=|gg|0+|ge|ϕgeL+|eg|ϕegL+𝒪(λ2),\displaystyle\ket{\psi_{L}}=\ket{gg}\ket{0}+\ket{ge}\ket{\phi_{ge}^{L}}+\ket{eg}\ket{\phi_{eg}^{L}}+\mathcal{O}(\lambda^{2}), (21)

where |ϕgeL/R\ket{\phi_{ge}^{L/R}} and |ϕegL/R\ket{\phi_{eg}^{L/R}} are first order in λ\lambda, describing field states containing a single excitation, see Supplementary Note 2 for their explicit expression.

I Results

I.1 Protocol – two cavities, four molecules

We now move to the full protocol, whose aim is to demonstrate indefinite causal order in a fully operational setting. We consider two ‘wings’ of the experiment, where in each wing two operations take place. The goal is to correlate the order of each pair of operations with a control system and, after measuring the control, produce an entangled state between the two wings, which would not have been possible had the operations been realised in a definite order.

Recall that we want to entangle the ‘target’ systems which comprise the fields in both cavities as well as the detectors (which are two per cavity), while the four spin-12\frac{1}{2} particles play the role of the control. The full protocol involves four molecules (two in each wing), where each molecule contains a clock, a detector, and a spin-12\frac{1}{2} particle. Initially, the state of these molecules reads:

|τ0τ0τ0τ012(|+|)|gggg,\displaystyle\ket{\tau_{0}\tau_{0}\tau_{0}\tau_{0}}\frac{1}{\sqrt{2}}(\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}+\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow\downarrow\mathrel{\mspace{-1.0mu}}\uparrow})\ket{gggg}, (22)

where |τ0\ket{\tau_{0}} is the initial state of one clock (and thus all four clocks are initially decorrelated and synchronised) and each detector is in its ground state, while the spins are entangled.

Note that this state is constructed in such a way that the four molecules can be divided into two identical pairs. These pairs are then accelerated as explained in the “Trajectories” section. After the process of acceleration, the state of all four molecules is:

12(|ττττ|+|ττττ|)|gggg.\displaystyle\frac{1}{\sqrt{2}}(\ket{\tau_{\uparrow}\tau_{\downarrow}\tau_{\uparrow}\tau_{\downarrow}}\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}+\ket{\tau_{\downarrow}\tau_{\uparrow}\tau_{\downarrow}\tau_{\uparrow}}\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow\downarrow\mathrel{\mspace{-1.0mu}}\uparrow})\ket{gggg}. (23)

Due to time dilation, after the acceleration takes place, one clock from each pair of molecules will be older than the other from the same pair. We cannot determine which one because the initial state of spins is a superposition of two different possibilities. Next, each pair of molecules enters a cavity. The first pair is placed in cavity 11 at the respective positions x1x_{1} and x2x_{2} (determined relative to the boundary of this cavity at x=0x=0). The second pair of molecules is placed in cavity 22 at the same positions relative to that cavity. Recall that the field in each cavity is initially in the vacuum state. As discussed in the “Two molecules in the cavity” section, we can decorrelate the clock DoF by re-synchronising the clocks. This is done after the molecules interact with the cavities. To keep track of the effect of the clocks on the detector-field interactions, it is enough to note that a detector belonging to a molecule with spin \downarrow interacts first (and with spin \uparrow interacts second). To sum up the protocol, the initial state of the total system including the two cavities is:

|Ψ0tot=12(|+|)|gg|0S1|gg|0S2,\displaystyle\ket{\Psi_{0}^{\mathrm{tot}}}=\frac{1}{\sqrt{2}}(\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}+\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow\downarrow\mathrel{\mspace{-1.0mu}}\uparrow})\underbrace{\ket{gg}\ket{0}}_{\in S_{1}}\otimes\underbrace{\ket{gg}\ket{0}}_{\in S_{2}}, (24)

where we used \otimes to separate states from the two cavities and showed which degrees of freedom comprise the targets S1S_{1} and S2S_{2} introduced in the “Methods” Section.

After the interactions between the detectors and the cavity fields, the state of the system reads:

|Ψtot=12(||ψR|ψR+||ψL|ψL),\displaystyle\ket{\Psi^{\mathrm{tot}}}=\frac{1}{\sqrt{2}}\left(\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}\ket{\psi_{R}}\otimes\ket{\psi_{R}}+\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow\downarrow\mathrel{\mspace{-1.0mu}}\uparrow}\ket{\psi_{L}}\otimes\ket{\psi_{L}}\right), (25)

where the detectors and cavity states, |ψL/R\ket{\psi_{L/R}}, are those in Eqs (20), (21).

Next the molecules are ‘broken apart’. All spins are send to a common location where, at an event labelled DD, they are jointly measured in the basis 12(|±|)\frac{1}{\sqrt{2}}(\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}\pm\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow\downarrow\mathrel{\mspace{-1.0mu}}\uparrow}). This prepares the remaining systems— cavities and detectors (which stay next to their cavities)—in the state

|Ψ±=12(|ψR|ψR±|ψL|ψL),\displaystyle\ket{\Psi^{\pm}}=\frac{1}{\sqrt{2}}\left(\ket{\psi_{R}}\otimes\ket{\psi_{R}}\pm\ket{\psi_{L}}\otimes\ket{\psi_{L}}\right), (26)

where the sign ±\pm depends on the measurement outcome. This state is entangled as long as |ψL|ψR\ket{\psi_{L}}\neq\ket{\psi_{R}}. This is typically shown by finding the scalar product between |ψL\ket{\psi_{L}} and |ψR\ket{\psi_{R}} which requires the Dyson series up to second order in λ\lambda.

However, there is an equivalent but technically much simpler method to show that the final state is in general entangled, which is to consider another measurement performed on each detector pair in a basis containing the vector 1/2(|ge+|eg)1/\sqrt{2}\left(\ket{ge}+\ket{eg}\right). The resulting conditional states of the fields (arising from Eqs (20) and (21)) read

12(ge|+eg|)|ψR=12(|ϕgeR+|ϕegR)+𝒪(λ2),\displaystyle\frac{1}{\sqrt{2}}(\bra{ge}+\bra{eg})\ket{\psi_{R}}=\frac{1}{\sqrt{2}}\left(\ket{\phi_{ge}^{R}}+\ket{\phi_{eg}^{R}}\right)+\mathcal{O}(\lambda^{2}), (27)
12(ge|+eg|)|ψL=12(|ϕgeL+|ϕegL)+𝒪(λ2).\displaystyle\frac{1}{\sqrt{2}}(\bra{ge}+\bra{eg})\ket{\psi_{L}}=\frac{1}{\sqrt{2}}\left(\ket{\phi_{ge}^{L}}+\ket{\phi_{eg}^{L}}\right)+\mathcal{O}(\lambda^{2}). (28)

Upon measuring the control in the above defined entangled basis, 12(|±|)\frac{1}{\sqrt{2}}(\ket{\uparrow\mathrel{\mspace{-1.0mu}}\downarrow\uparrow\mathrel{\mspace{-1.0mu}}\downarrow}\pm\ket{\downarrow\mathrel{\mspace{-1.0mu}}\uparrow\downarrow\mathrel{\mspace{-1.0mu}}\uparrow}), the joint state of two cavities fields up to the leading order is:

|Ψ~±(|ϕgeR+|ϕegR)|ΦR(|ϕgeR+|ϕegR)|ΦR±(|ϕgeL+|ϕegL)|ΦL(|ϕgeL+|ϕegL)|ΦL.\displaystyle\ket{\tilde{\Psi}^{\pm}}\propto\underbrace{\left(\ket{\phi_{ge}^{R}}+\ket{\phi_{eg}^{R}}\right)}_{\ket{\Phi_{R}}}\otimes\underbrace{\left(\ket{\phi_{ge}^{R}}+\ket{\phi_{eg}^{R}}\right)}_{\ket{\Phi_{R}}}\pm\underbrace{\left(\ket{\phi_{ge}^{L}}+\ket{\phi_{eg}^{L}}\right)}_{\ket{\Phi_{L}}}\otimes\underbrace{\left(\ket{\phi_{ge}^{L}}+\ket{\phi_{eg}^{L}}\right)}_{\ket{\Phi_{L}}}. (29)

Note that the states |ΦL\ket{\Phi_{L}} and |ΦR\ket{\Phi_{R}} describe only the fields in the cavities. The conditional state |Ψ~±\ket{\tilde{\Psi}^{\pm}} is entangled as long as |ΦR|ΦL\ket{\Phi_{R}}\neq\ket{\Phi_{L}} which is easier to verify than the condition |ψR|ψL\ket{\psi_{R}}\neq\ket{\psi_{L}}. We further note that the measurement on the detectors has different possible outcomes, not all of which result in entangled field states. Crucially however, for a product state across S1S_{1} and S2S_{2}, all local detector measurement would result in product states of the fields. Therefore, obtaining entangled state for one of the outcomes is sufficient to prove that the original state was entangled.

To sum up, we have shown that using two cavities, four molecules, and the basic effects of special relativity, we can obtain the state given in Eq. (29). If this state is entangled, in the next steps of the protocol this entanglement can be used to violate Bell’s inequalities with the goal to prove that the detector–field interaction events were not classically ordered.

I.2 Entanglement of the final state

We now proceed to prove that the state in Eq. (29) can indeed be entangled by showing that there exist parameters for which |ΦR|ΦL\ket{\Phi_{R}}\neq\ket{\Phi_{L}}. The scalar product between these states reads:

ΦR|ΦL=keiΔτ(ωk+Ω)(ωk+Ω)2[1cos(T(ωk+Ω))](uk(x1)+eiΔτ(ωk+Ω)uk(x2))22nsin2T(ωn+Ω)2(ωn+Ω)2[un(x2)2+2un(x2)un(x1)cos(Δτ(ωn+Ω))+un(x1)2].\displaystyle\innerproduct{\Phi_{R}}{\Phi_{L}}=\frac{\sum_{k}\frac{e^{-i\Delta\tau(\omega_{k}+\Omega)}}{(\omega_{k}+\Omega)^{2}}\left[1-\cos\left(T(\omega_{k}+\Omega)\right)\right]\left(u_{k}(x_{1})+e^{i\Delta\tau(\omega_{k}+\Omega)}u_{k}(x_{2})\right)^{2}}{2\sum_{n}\frac{\sin^{2}\frac{T(\omega_{n}+\Omega)}{2}}{(\omega_{n}+\Omega)^{2}}\left[u_{n}(x_{2})^{2}+2u_{n}(x_{2})u_{n}(x_{1})\cos\left(\Delta\tau(\omega_{n}+\Omega)\right)+u_{n}(x_{1})^{2}\right]}\,. (30)

Explicit derivation is presented in Supplementary Note 3. Therein we further show how to choose parameters for which the states are not just different but orthogonal. To summarise this procedure, the parameters that have to be chosen are:

  1. 1.

    The length LL of each cavity and the positions x1,x2x_{1},x_{2} of the molecules relative to their respective cavity.

  2. 2.

    The energy gap Ω\Omega of the detectors.

  3. 3.

    The duration TT of the interaction between each detector and the cavity field.

  4. 4.

    The time dilation Δτ\Delta\tau between clocks within each pair of molecules (within each cavity), arising from the different accelerations 𝒜/\mathcal{A}_{\uparrow/\downarrow}.

Remarkably, it is also possible to find parameters such that |ΦR|ΦL\ket{\Phi_{R}}\neq\ket{\Phi_{L}} even when the two field-detector interactions within each cavity are space-like separated, see Fig. 4. This means that it is possible to obtain an entangled state also in the case when, in some reference frame, the relevant operations are performed in the same temporal order. Although from the perspective of such a reference frame the operations would take place over four different coordinate regions, their order would be the same in each of the amplitudes. Obviously then, entanglement generated in this scheme cannot be simply attributed to non-classical temporal order. We discuss implications and argue for generality of this result in the next sections.

Refer to caption
Figure 4: Spacetime diagrams of interacting detectors. Spacetime diagram showing the regions where the detector-field interactions lead to an entangled state, Eq. (29). The diagrams are for one cavity – the interaction regions are idenatically defined for the second cavity. The first row shows spacelike separated regions which nevertheless yield entanglement. The second row shows timelike separation which yields maximal entanglement. The columns correspond to the two amplitudes of the process that are superposed using the control (spin) state. See Supplementary Note 3 for the supporting calculations.

I.3 Ambiguity in the signature of indefinite temporal order

We have explicitly modelled a special-relativistic version of a protocol where gravitational time dilation and quantum superposition lead to an indefinite temporal order of events from Ref. Bell4time . This protocol was formulated in terms of a ‘Bell inequality for temporal order’. The idea was to formulate a protocol where operations in a definite order cannot produce an entangled state, if an appropriate set of assumptions is satisfied. (The final step of the protocol requires measurements on the state in order to violate a Bell inequality and verify the entanglement.) The motivation was to find a test of temporal order that is theory independent, which was based on the observation that a violation of a Bell inequality would prove indefinite temporal order without assuming that the final state is described by quantum mechanics. In this section, we reexamine the assumptions made in Ref. Bell4time , showing that the theory-independent nature of the argument is problematic.

We first remark that the issue we have identified does not arise if we believe that quantum mechanics is valid, i.e., if we wish to provide experimental evidence for non-classical temporal order assuming that the involved states and transformations are faithfully described by quantum theory. In such a case, already the one-cavity part of our protocol, the “Methods” section, or the gravitational equivalent, is sufficient. The reason is that at an abstract level, these protocols implement a “quantum switch” chiribella09 – a scenario where two local operations, represented by unitary operators U^A\hat{U}^{A} U^B\hat{U}^{B}, act on a target system in an order determined by a control system which is prepared in a superposition, thus producing a final state of the form

|ψfin=12(|0U^AV^0U^B+|1U^BV^1U^A)|ψ.\ket{\psi_{\textrm{fin}}}=\\ \frac{1}{\sqrt{2}}\left(\ket{0}\hat{U}^{A}\hat{V}_{0}\hat{U}^{B}+\ket{1}\hat{U}^{B}\hat{V}_{1}\hat{U}^{A}\right)\ket{\psi}. (31)

Here, |0\ket{0}, |1\ket{1} are two basis states of the control, and |ψ\ket{\psi} is the initial state of the target system. V^0\hat{V}_{0} and V^1\hat{V}_{1} are two arbitrary unitary operators, representing the evolution of the target between the two operations. Most presentations of the switch do not include the intermediate evolution, but we will see shortly that this is important in our context222For generality, one can also include a control-dependent initial state, but this is not necessary for our analysis.. Given an implementation of the switch, by making appropriate final measurements for a set of suitably chosen operations U^A\hat{U}^{A}, U^B\hat{U}^{B} it is possible to prove that the operations are not performed in a definite order. This procedure is known as measuring a causal witness araujo15 ; Branciard2016a and it could in principle be incorporated within a single-cavity variant of our protocol to demonstrate that, if the quantum description of the experiment is correct, the cavity-detector interactions do not take place in a definite order.

Returning to the Bell inequality approach, it essentially is an ‘entangled’ version of the switch, with the final state of the form

|ψfin=12(|0U^A1V^0U^B1|ψU^A2V^0U^B2|ψ+|1U^B1V^1U^A1|ψU^B2V^1U^A2|ψ).\ket{\psi_{\textrm{fin}}}=\\ \frac{1}{\sqrt{2}}\left(\ket{0}\hat{U}^{A_{1}}\hat{V}_{0}\hat{U}^{B_{1}}\ket{\psi}\otimes\hat{U}^{A_{2}}\hat{V}_{0}\hat{U}^{B_{2}}\ket{\psi}\right.\\ +\left.\ket{1}\hat{U}^{B_{1}}\hat{V}_{1}\hat{U}^{A_{1}}\ket{\psi}\otimes\hat{U}^{B_{2}}\hat{V}_{1}\hat{U}^{A_{2}}\ket{\psi}\right). (32)

In order to use such a state to disprove classical temporal order among U^A\hat{U}^{A}, U^B\hat{U}^{B}, the state has to arise in a scenario satisfying all assumptions used to derive the Bell inequality for temporal order, apart from the assumption of classical order itself. All assumptions used in Ref Bell4time for the derivation of the inequality are: the initial state of the target systems S1,S2S_{1},S_{2} is separable; transformations performed on the targets are local (i.e. the operations on target SjS_{j} act as identity on other degrees of freedom of the system); events/spacetime regions at which transformations and measurements take place are suitably separated: both interaction events in one wing are spacelike separated from both interaction events in the other wing, and event D\mathrm{D} is spacelike from the events at which the Bell measurements are performed (which as, usual in Bell inequalities, are assumed to be spacelike from each other); the choices of bases for Bell measurements are independent of all other aspects of the experiment (often referred to as ‘free choice’ assumption); and, finally, of course the assumption that events at which transformations and measurements are performed are classically ordered.

The surprising result identified in the “Entanglement of the final state” section is that entanglement is produced while all the above assumptions are met, including the assumption of classical order. Clearly, some other assumption was made to derive the inequality and is violated in our implementation. Indeed, as already mentioned, the additional assumption made in Ref. Bell4time is that the target systems do not have non-trivial evolution apart from the transformations U^A\hat{U}^{A} and U^B\hat{U}^{B}. Below we explain why this assumption is violated in the present implementation and in the “Ambiguity in the signature of indefinite temporal order” section we argue that this will remain true in a generic dynamical implementation.

In order to see where the assumption of no free evolution enters and why it is the culprit it is again sufficient to look at one wing of our setup. Comparing a single-cavity scenario, Eq. (18), and a generic quantum switch, Eq. (31), one finds that the evolution operators U^1\hat{U}_{1}, U^2\hat{U}_{2} in Eq. (18) are not directly representing the local operations U^A\hat{U}^{A}, U^B\hat{U}^{B}. The reason is that U^1\hat{U}_{1}, U^2\hat{U}_{2} are written in the Dirac (interaction) picture, which necessarily includes time evolution with respect to the free Hamiltonian (starting from some initial time established in a common reference frame). By unravelling this time evolution, one finds precisely a state of the form (31), where V^0\hat{V}_{0} and V^1\hat{V}_{1} represent the free evolution of the targets (i.e. cavity and detectors) between the interactions. Note that the time intervals between the events, and thus intervals of free evolution, are equal in the reference frame of the cavity, which means that in that frame V^0=V^1\hat{V}_{0}=\hat{V}_{1}. On the other hand, U^A\hat{U}^{A}, U^B\hat{U}^{B} describe only the field-cavity interactions in the Schrödinger picture333Strictly speaking, we should add free-evolution operators also before and after the two local unitaries—not only in between. However, free evolution acts trivially on the vacuum state (our initial state), while the final evolution can be reabsorbed in the definition of the measurement basis..

The reason this is relevant here is that if the free Hamiltonian does not commute with the interaction, the free evolution does not commute with U^A,U^B\hat{U}^{A},\hat{U}^{B}, and this is how the presence of entanglement can be explained in the state in Eq. (32) in a frame where V^0=V^1\hat{V}_{0}=\hat{V}_{1} and the events are spacelike separated. Note that in that case there is a reference frame where AjBjA_{j}\prec B_{j} for both states of the control, however, as we mentioned above, in that frame necessarily V^0V^1\hat{V}_{0}\neq\hat{V}_{1}, since the time intervals of free evolution along the worldlines of the molecule 1 and 2 are in such a frame necessarily different. In such a frame the final state in the two wing scenario becomes

|ψfin=12(|0U^B1V^0U^A1|ψU^B2V^0U^A2|ψ+|1U^B1V^1U^A1|ψU^B2V^1U^A2|ψ),\ket{\psi_{\textrm{fin}}}=\\ \frac{1}{\sqrt{2}}\left(\ket{0}\hat{U}^{B_{1}}\hat{V}_{0}\hat{U}^{A_{1}}\ket{\psi}\otimes\hat{U}^{B_{2}}\hat{V}_{0}\hat{U}^{A_{2}}\ket{\psi}\right.\\ +\left.\ket{1}\hat{U}^{B_{1}}\hat{V}_{1}\hat{U}^{A_{1}}\ket{\psi}\otimes\hat{U}^{B_{2}}\hat{V}_{1}\hat{U}^{A_{2}}\ket{\psi}\right), (33)

and the presence of entanglement is thus interpreted as due to overall different dynamics depending on the control, U^BiV^0U^AiU^BiV^1U^Ai\hat{U}^{B_{i}}\hat{V}_{0}\hat{U}^{A_{i}}\neq\hat{U}^{B_{i}}\hat{V}_{1}\hat{U}^{A_{i}} i.e., in that frame while the order of operations is common, when they take place relative to periods of free dynamics, and thus the overll evolution, depends on the control. Crucially, in this case temporal order among events cannot be even defined as changes depending on the reference frame. A further discussion of the role played by free evolution in this protocol is presented in Supplementary Note 4.

In fact, if free evolution and the applied operations do not commute, even simpler scenarios can illustrate the issue. Consider that one of the operations is trivial, say U^B=𝕀^\hat{U}^{B}=\hat{\mathbb{I}}, and so in fact only one operation is applied. The non-commutativity between U^A\hat{U}_{A} and V^=V^0=V^1\hat{V}=\hat{V}_{0}=\hat{V}_{1} (we are in the reference frame of the cavity) would again result in different final states depending on when U^A\hat{U}^{A} is applied relative to V^\hat{V}. This would again lead to an entangled final state in a two-wing scenario, even though in this case there is even no time order of events to speak of.

Finally, we note that in all quantum switch scenarios, including the entangled switch, there is an assumption that each local operation is performed ‘only once’. This condition is, however, naturally satisfied in the relativistic implementations such as ours, as each operation is performed at a specific time of a local clock.

Conclusion

In this paper, we constructed a non-gravitational scenario where accelerating particles, interacting with quantum fields, according to their own internal clock degrees of freedom, can lead to a violation of the temporal Bell inequalities analogous to the gravitational case. In the “Methods” section we introduced the formalism required to reproduce the gravitational protocol using special relativistic time dilation. We defined the kinematics of all particles involved in the protocol and the appropriate coupling between them and the quantum field. In the “Results” section, we proposed the full protocol, that demonstrated indefinite causal order in a fully operational setting. We described the procedure than would lead us to the violation of Bell’s inequalities for the proposed system. We found that a violation of Bell inequalities occurs also when the entangled events are space-like what we interpreted as an ambiguity in the signature of indefinite temporal order. We finally found that this surprising conclusion is the result of failure of the additional assumption, requiring that target systems have no other evolution except, the one governed by unitaries applied in a specific time order. A detailed discussion of this problem we presented in the “Ambiguity in the signature of indefinite temporal order” section of our work.

We have focused so far on a particular realisation of the entangled switch protocol – with two-level ‘detectors’ and cavity-confined quantum field modes as targets, and the position DoFs of the detectors as the control. However, our main result and its explanation applies to any physical realisation of the protocol. Indeed, we have shown that entanglement can be generated for spacelike separated operations and identified that this is due to free evolution of the targets. For any physical implementation of the protocol, if the applied operations do not commute with the free evolution, the final state will in general be entangled regardless of the commutation relations between the operations themselves, and thus also regardless of their temporal order. Entanglement can arise even if only one operation is applied, as discussed in the previous section.

It is of course possible to identify implementations where the assumption of no free evolution of the targets does hold (e.g. with polarisation or angular momentum DoFs of photons) or where the operations are performed within a degenerate subspace of the Hamiltonian of the targets (and thus commute with free evolution), so that only non-commutativity between the two local operations is relevant. Such strategies would indeed lead to an entangled final state only if the local operations are timelike and applied in a non-classical order. However, identifying such implementations does require a theoretical description of the states and dynamics of the involved systems. In other words, one needs theory-dependent assumptions to interpret a violation of the final Bell inequalities as a signature of indefinite temporal order.

The observations above do not depend on the special-relativistic setting studied in this work and hold equally for gravitational protocols. Indeed, nothing in the argument depends on how the time dilation of the clocks is achieved. For example, the position of a massive body can determine – through gravitational time dilation – the time at which a single operation takes place (relative to some mass-independent coordinates). This can lead to the same situation as described earlier: creation of entanglement (in a two-wing scenario) even with a single operation per wing. Again, theory-dependent assumptions would be required to ensure that entanglement can only arise as a result of an indefinite order of events.

Being able to describe and experimentally test entangled temporal order (i.e. even working fully within quantum mechanics) is of interest in its own right. Our result, however, opens the question whether it is possible to formulate a stronger, theory independent, test of temporal order. The insight from the present study is that it is problematic to separate out the effect of the free dynamics of the system from that of the local operations. A possible avenue to circumvent this is to consider more general operations than the fixed unitaries discussed thus far: they can involve a measurement of the system, producing a classical variable as outcome. Furthermore, a setting variable for each party can model a choice among different operations. In this fashion, one can consider directly the causal relations between parties – understood operationally as correlations between settings and outcomes – without relying on a theory-dependent description of the transformations. We leave further investigation of this possibility to future work.

II acknowledgements

K.D. is financially supported by the (Polish) National Science Center Grant 2021/41/N/ST2/01901. F.C. acknowledges support through the Australian Research Council (ARC) DECRA grant DE170100712, M.Z. acknowledges support through ARC Future Fellowship grant FT210100675, F.C.and M.Z. acknowledge support through ARC Centre of Excellence EQuS CE170100009. The University of Queensland (UQ) acknowledges the Traditional Owners and their custodianship of the lands on which UQ operates.

References

  • 1 Butterfield, J. & Isham, C. Spacetime and the philosophical challenge of quantum gravity. In Callender, C. & Huggett, N. (eds.) Physics meets philosophy at the Planck scale, 33 (Cambridge University Press, 2001).
  • 2 Hardy, L. Towards quantum gravity: a framework for probabilistic theories with non-fixed causal structure. Journal of Physics A: Mathematical and Theoretical 40, 3081 (2007). URL https://iopscience.iop.org/article/10.1088/1751-8113/40/12/S12https://iopscience.iop.org/article/10.1088/1751-8113/40/12/S12/meta.
  • 3 Oriti, D. (ed.) Approaches to Quantum Gravity:
    Toward a New Understanding of Space, Time and Matter
    (Cambridge University Press, 2009).
  • 4 Kiefer, C. Quantum Gravity: Third Edition. International Series of Monographs on Physics (OUP Oxford, 2012).
  • 5 Zych, M., Costa, F., Pikovski, I. & Časlav Brukner. Bell’s theorem for temporal order. Nature Communications 2019 10:1 10, 1–10 (2019). URL https://www.nature.com/articles/s41467-019-11579-x.
  • 6 Procopio, L. M. et al. Experimental superposition of orders of quantum gates. Nature Communications 2015 6:1 6, 1–6 (2015). URL https://www.nature.com/articles/ncomms8913.
  • 7 Rubino, G. et al. Experimental verification of an indefinite causal order. Science Advances 3 (2017). URL https://www.science.org/doi/abs/10.1126/sciadv.1602589.
  • 8 Rubino, G. et al. Experimental entanglement of temporal orders. Optics InfoBase Conference Papers Part F165-QIM 2019 (2017). URL https://arxiv.org/abs/1712.06884v4.
  • 9 Goswami, K. et al. Indefinite causal order in a quantum switch. Physical Review Letters 121, 090503 (2018). URL https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.121.090503.
  • 10 Guo, Y. et al. Experimental investigating communication in a superposition of causal orders (2018). eprint 1811.07526.
  • 11 Wei, K. et al. Experimental quantum switching for exponentially superior quantum communication complexity. Physical Review Letters 122, 120504 (2019). URL https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.122.120504.
  • 12 Taddei, M. M. et al. Computational advantage from the quantum superposition of multiple temporal orders of photonic gates. PRX Quantum 2, 010320 (2021). URL https://journals.aps.org/prxquantum/abstract/10.1103/PRXQuantum.2.010320.
  • 13 Zych, M., Costa, F., Pikovski, I. & Brukner, Č. Quantum interferometric visibility as a witness of general relativistic proper time. Nature Communications 2011 2:1 2, 1–7 (2011). URL https://www.nature.com/articles/ncomms1498.
  • 14 Pikovski, I., Zych, M., Costa, F. & Brukner, C. Time dilation in quantum systems and decoherence. New Journal of Physics 19 (2017).
  • 15 Margalit, Y. et al. Realization of a complete stern-gerlach interferometer: Toward a test of quantum gravity. Science Advances 7, eabg2879 (2021). eprint https://www.science.org/doi/pdf/10.1126/sciadv.abg2879.
  • 16 Yurke, B. & Stoler, D. Bell’s-inequality experiments using independent-particle sources. Phys. Rev. A 46, 2229–2234 (1992). URL https://link.aps.org/doi/10.1103/PhysRevA.46.2229.
  • 17 Dragan, A. Unusually Special Relativity (WORLD SCIENTIFIC (EUROPE), 2021).
  • 18 Unruh, W. G. Notes on black-hole evaporation. Physical Review D 14, 870 (1976). URL https://journals.aps.org/prd/abstract/10.1103/PhysRevD.14.870.
  • 19 Unruh, W. G. & Wald, R. M. What happens when an accelerating observer detects a rindler particle. Physical Review D 29, 1047–1056 (1984). URL https://journals.aps.org/prd/abstract/10.1103/PhysRevD.29.1047.
  • 20 Crispino, L. C. B., Higuchi, A. & Matsas, G. E. A. The unruh effect and its applications. Rev. Mod. Phys. 80, 787–838 (2008).
  • 21 Chiribella, G., D’Ariano, G. M., Perinotti, P. & Valiron, B. Quantum computations without definite causal structure. Phys. Rev. A 88, 022318 (2013).
  • 22 Araújo, M. et al. Witnessing causal nonseparability. New J. Phys. 17, 102001 (2015).
  • 23 Branciard, C. Witnesses of causal nonseparability: an introduction and a few case studies. Scientific Reports 6 (2016).

Supplementary Information: Indefinite temporal order – cavity model

XXX

XXX

XXX

XXX

Supplementary Note 1: Second order of the Dyson series

In this section, We will show that second-order Dyson expansion does not affect on |ΦR\ket{\Phi_{R}} and |ΦL\ket{\Phi_{L}}.

One can ask the question why (47) depends on λ2\lambda^{2}. Previously we limited our calculation only to the first-order expansion of the Dyson series. We have to verify that second-order expansion does not produce lambda square terms too. Otherwise, there is a possibility that results from (47) will cancel out with these additional terms. Let us write the general evolution operator in the following form:

U^=𝟙iλdtχd(t)(eiΩtσ^++eiΩtσ^)n(a^nun(xd)eiωnt+a^nun(xd)eiωnt)\displaystyle\hat{U}=\mathds{1}-i\lambda\int_{-\infty}^{\infty}\mathrm{d}t~{}\chi_{d}(t)\left(e^{i\Omega t}\hat{\sigma}^{+}+e^{-i\Omega t}\hat{\sigma}^{-}\right)\sum_{n}\left(\hat{a}_{n}^{\dagger}u_{n}(x_{d})e^{i\omega_{n}t}+\hat{a}_{n}u_{n}(x_{d})e^{-i\omega_{n}t}\right) (34)
λ2dt2t2dt1χd(t2)χd(t1)(eiΩt2σ^++eiΩt2σ^)(eiΩt1σ^++eiΩt1σ^)\displaystyle-\lambda^{2}\int_{-\infty}^{\infty}\mathrm{d}t_{2}\int_{-\infty}^{t_{2}}\mathrm{d}t_{1}~{}\chi_{d}(t_{2})\chi_{d}(t_{1})\left(e^{i\Omega t_{2}}\hat{\sigma}^{+}+e^{-i\Omega t_{2}}\hat{\sigma}^{-}\right)\left(e^{i\Omega t_{1}}\hat{\sigma}^{+}+e^{-i\Omega t_{1}}\hat{\sigma}^{-}\right)
×n(a^nun(xd)eiωnt2+a^nun(xd)eiωnt2)m(a^mum(xd)eiωmt1+a^mum(xd)eiωmt1),\displaystyle\times\sum_{n}\left(\hat{a}_{n}^{\dagger}u_{n}(x_{d})e^{i\omega_{n}t_{2}}+\hat{a}_{n}u_{n}(x_{d})e^{-i\omega_{n}t_{2}}\right)\sum_{m}\left(\hat{a}_{m}^{\dagger}u_{m}(x_{d})e^{i\omega_{m}t_{1}}+\hat{a}_{m}u_{m}(x_{d})e^{-i\omega_{m}t_{1}}\right),

where xdx_{d} is a position of a detector and σ^±\hat{\sigma}^{\pm} acts on the internal state of this detector. In our case we have two detectors and two Hilbert spaces of internal degrees of freedom. We can simplify the notation to write operators of the evolution as:

U1^=𝟙iλn(U^1n+σ^1++U^1nσ^1)λ2n,m(U^1nm++σ^1+σ^1++U^1nm+σ^1+σ^1+U^1nm+σ^1σ^1++U^1nmσ^1σ^1),\displaystyle\hat{U_{1}}=\mathds{1}-i\lambda\sum_{n}\left(\hat{U}^{+}_{1n}\hat{\sigma}^{+}_{1}+\hat{U}^{-}_{1n}\hat{\sigma}^{-}_{1}\right)-\lambda^{2}\sum_{n,m}\left(\hat{U}^{++}_{1nm}\hat{\sigma}^{+}_{1}\hat{\sigma}^{+}_{1}+\hat{U}^{+-}_{1nm}\hat{\sigma}^{+}_{1}\hat{\sigma}^{-}_{1}+\hat{U}^{-+}_{1nm}\hat{\sigma}^{-}_{1}\hat{\sigma}^{+}_{1}+\hat{U}^{--}_{1nm}\hat{\sigma}^{-}_{1}\hat{\sigma}^{-}_{1}\right), (35)
U2^=𝟙iλn(U^2n+σ^2++U^2nσ^2)λ2n,m(U^2nm++σ^2+σ^2++U^2nm+σ^2+σ^2+U^2nm+σ^2σ^2++U^2nmσ^2σ^2),\displaystyle\hat{U_{2}}=\mathds{1}-i\lambda\sum_{n}\left(\hat{U}^{+}_{2n}\hat{\sigma}^{+}_{2}+\hat{U}^{-}_{2n}\hat{\sigma}^{-}_{2}\right)-\lambda^{2}\sum_{n,m}\left(\hat{U}^{++}_{2nm}\hat{\sigma}^{+}_{2}\hat{\sigma}^{+}_{2}+\hat{U}^{+-}_{2nm}\hat{\sigma}^{+}_{2}\hat{\sigma}^{-}_{2}+\hat{U}^{-+}_{2nm}\hat{\sigma}^{-}_{2}\hat{\sigma}^{+}_{2}+\hat{U}^{--}_{2nm}\hat{\sigma}^{-}_{2}\hat{\sigma}^{-}_{2}\right), (36)

where U1n±U^{\pm}_{1n}, U2n±U^{\pm}_{2n} are operators from the first-order order expansion and U1nm±±U^{\pm\pm}_{1nm} and U2nm±±U^{\pm\pm}_{2nm} are operators from the second-order expansion of the Dyson series. σ^1±\hat{\sigma}^{\pm}_{1} and σ^2±\hat{\sigma}^{\pm}_{2} are operators σ^±\hat{\sigma}^{\pm} acting on the internal space of the first or the second detector respectivly. To find contributions proportional to the λ2\lambda^{2} to the value of (47), we have to find new terms in the |ψR\ket{\psi_{R}} and |ψL\ket{\psi_{L}} proportional to |ge\ket{ge} or |eg\ket{eg}. Knowing that |ψR=U^1U^2|gg|0\ket{\psi_{R}}=\hat{U}_{1}\hat{U}_{2}\ket{gg}\ket{0}, |ψL=U^2U^1|gg|0\ket{\psi_{L}}=\hat{U}_{2}\hat{U}_{1}\ket{gg}\ket{0}, σ^|g=0\hat{\sigma}^{-}\ket{g}=0 and σ^+|e=0\hat{\sigma}^{+}\ket{e}=0 we can write:

|ψR=[𝟙iλnU^1n+σ^1+λ2n,mU^1nm+σ^1σ^1+][𝟙iλnU^2n+σ^2+λ2n,mU^2nm+σ^2σ^2+]|gg|0,\displaystyle\ket{\psi_{R}}=\left[\mathds{1}-i\lambda\sum_{n}\hat{U}^{+}_{1n}\hat{\sigma}^{+}_{1}-\lambda^{2}\sum_{n,m}\hat{U}^{-+}_{1nm}\hat{\sigma}^{-}_{1}\hat{\sigma}^{+}_{1}\right]\left[\mathds{1}-i\lambda\sum_{n}\hat{U}^{+}_{2n}\hat{\sigma}^{+}_{2}-\lambda^{2}\sum_{n,m}\hat{U}^{-+}_{2nm}\hat{\sigma}^{-}_{2}\hat{\sigma}^{+}_{2}\right]\ket{gg}\ket{0}, (37)
|ψL=[𝟙iλnU^2n+σ^2+λ2n,mU^2nm+σ^2σ^2+][𝟙iλnU^1n+σ^1+λ2n,mU^1nm+σ^1σ^1+]|gg|0.\displaystyle\ket{\psi_{L}}=\left[\mathds{1}-i\lambda\sum_{n}\hat{U}^{+}_{2n}\hat{\sigma}^{+}_{2}-\lambda^{2}\sum_{n,m}\hat{U}^{-+}_{2nm}\hat{\sigma}^{-}_{2}\hat{\sigma}^{+}_{2}\right]\left[\mathds{1}-i\lambda\sum_{n}\hat{U}^{+}_{1n}\hat{\sigma}^{+}_{1}-\lambda^{2}\sum_{n,m}\hat{U}^{-+}_{1nm}\hat{\sigma}^{-}_{1}\hat{\sigma}^{+}_{1}\right]\ket{gg}\ket{0}. (38)

We can observe that second-order contributions from the same detector do not change the internal state, i.e. σ^σ^+|g=|g\hat{\sigma}^{-}\hat{\sigma}^{+}\ket{g}=\ket{g}, while the product of first order terms from both detectors gives σ^1+σ^2+|gg=σ^2+σ^1+|gg=|ee\hat{\sigma}_{1}^{+}\hat{\sigma}_{2}^{+}\ket{gg}=\hat{\sigma}_{2}^{+}\hat{\sigma}_{1}^{+}\ket{gg}=\ket{ee}. Thus, the second order terms do not have support on the subspace spanned by the states |ge,|eg\ket{ge},\ket{eg}. Note that the states |ge\ket{ge} or |eg\ket{eg} appear in this expansion at order λ3\lambda^{3} or higher. This means that second-order terms from the Dyson series do not affect the states |ΦR\ket{\Phi_{R}} and |ΦL\ket{\Phi_{L}} in our approximation.

Supplementary Note 2: Details of calculations of the final state

Using the form of evolution operator (10) we can find that:

|ψR=U^1U^2|g|g|0=|g|g|0iλdt2χ2R(t2)k1ωkLei(ωk+Ω)t2sin(kπLx2)|g|ea^k|0\displaystyle\ket{\psi_{R}}=\hat{U}_{1}\hat{U}_{2}\ket{g}\ket{g}\ket{0}=\ket{g}\ket{g}\ket{0}-i\lambda\int\mathrm{d}t_{2}\chi_{2R}(t_{2})\sum_{k}\frac{1}{\sqrt{\omega_{k}L}}e^{i(\omega_{k}+\Omega)t_{2}}\sin(\frac{k\pi}{L}x_{2})\ket{g}\ket{e}\hat{a}_{k}^{\dagger}\ket{0}
iλdt1χ1R(t1)k1ωkLei(ωk+Ω)t1sin(kπLx1)|e|ga^k|0+𝒪(λ2),\displaystyle-i\lambda\int\mathrm{d}t_{1}\chi_{1R}(t_{1})\sum_{k}\frac{1}{\sqrt{\omega_{k}L}}e^{i(\omega_{k}+\Omega)t_{1}}\sin(\frac{k\pi}{L}x_{1})\ket{e}\ket{g}\hat{a}_{k}^{\dagger}\ket{0}+\mathcal{O}(\lambda^{2}), (39)

where: χ1R\chi_{1R}–switching function for the first detector in the case that right detector interacts before the left one, χ2R\chi_{2R}–switching function for the second detector in the case that right detector interacts earlier. We assume that the interaction starts and ends rapidly, so that χ2R(t)=1\chi_{2R}(t)=1 for t(0,T)t\in(0,T) and χ2R(t)=0\chi_{2R}(t)=0 for any other time. Similarly, χ1R(t)=1\chi_{1R}(t)=1 for t(Δτ,Δτ+T)t\in(\Delta\tau,\Delta\tau+T) and χ1R(t)=0\chi_{1R}(t)=0 for any other time. We can proceed with the same calculation for |ψL\ket{\psi_{L}}:

|ψL=U^2U^1|g|g|0=|g|g|0iλdt2χ2L(t2)k1ωkLei(ωk+Ω)t2sin(kπLx2)|g|ea^k|0\displaystyle\ket{\psi_{L}}=\hat{U}_{2}\hat{U}_{1}\ket{g}\ket{g}\ket{0}=\ket{g}\ket{g}\ket{0}-i\lambda\int\mathrm{d}t_{2}\chi_{2L}(t_{2})\sum_{k}\frac{1}{\sqrt{\omega_{k}L}}e^{i(\omega_{k}+\Omega)t_{2}}\sin(\frac{k\pi}{L}x_{2})\ket{g}\ket{e}\hat{a}_{k}^{\dagger}\ket{0}
iλdt1χ1L(t1)k1ωkLei(ωk+Ω)t1sin(kπLx1)|e|ga^k|0+𝒪(λ2),\displaystyle-i\lambda\int\mathrm{d}t_{1}\chi_{1L}(t_{1})\sum_{k}\frac{1}{\sqrt{\omega_{k}L}}e^{i(\omega_{k}+\Omega)t_{1}}\sin(\frac{k\pi}{L}x_{1})\ket{e}\ket{g}\hat{a}_{k}^{\dagger}\ket{0}+\mathcal{O}(\lambda^{2}), (40)

where: χ1L\chi_{1L}–switching function for the first detector in the case that right detector interacts before the left one, χ2L\chi_{2L}–switching function for the second detector in the case that right detector interacts earlier. In this case we have that χ1L(t)=1\chi_{1L}(t)=1 for t(0,T)t\in(0,T) and χ1L(t)=0\chi_{1L}(t)=0 for any other time. Similarly, χ2L(t)=1\chi_{2L}(t)=1 for t(Δτ,Δτ+T)t\in(\Delta\tau,\Delta\tau+T) and χ2L(t)=0\chi_{2L}(t)=0 for any other time. It worth noticing that χ1L=χ2R\chi_{1L}=\chi_{2R} and χ2L=χ1R\chi_{2L}=\chi_{1R}. After using this property describing relation between switching functions, we have:

|ψR=U^1U^2|g|g|0=|g|g|0iλdt2χ2R(t2)k1ωkLei(ωk+Ω)t2sin(kπLx2)|g|ea^k|0\displaystyle\ket{\psi_{R}}=\hat{U}_{1}\hat{U}_{2}\ket{g}\ket{g}\ket{0}=\ket{g}\ket{g}\ket{0}-i\lambda\int\mathrm{d}t_{2}\chi_{2R}(t_{2})\sum_{k}\frac{1}{\sqrt{\omega_{k}L}}e^{i(\omega_{k}+\Omega)t_{2}}\sin(\frac{k\pi}{L}x_{2})\ket{g}\ket{e}\hat{a}_{k}^{\dagger}\ket{0}
iλdt1χ1R(t1)k1ωkLei(ωk+Ω)t1sin(kπLx1)|e|ga^k|0+𝒪(λ2),\displaystyle-i\lambda\int\mathrm{d}t_{1}\chi_{1R}(t_{1})\sum_{k}\frac{1}{\sqrt{\omega_{k}L}}e^{i(\omega_{k}+\Omega)t_{1}}\sin(\frac{k\pi}{L}x_{1})\ket{e}\ket{g}\hat{a}_{k}^{\dagger}\ket{0}+\mathcal{O}(\lambda^{2}), (41)
|ψL=U^2U^1|g|g|0=|g|g|0iλdt2χ2R(t2)k1ωkLei(ωk+Ω)t2sin(kπLx1)|e|ga^k|0\displaystyle\ket{\psi_{L}}=\hat{U}_{2}\hat{U}_{1}\ket{g}\ket{g}\ket{0}=\ket{g}\ket{g}\ket{0}-i\lambda\int\mathrm{d}t_{2}\chi_{2R}(t_{2})\sum_{k}\frac{1}{\sqrt{\omega_{k}L}}e^{i(\omega_{k}+\Omega)t_{2}}\sin(\frac{k\pi}{L}x_{1})\ket{e}\ket{g}\hat{a}_{k}^{\dagger}\ket{0}
iλdt1χ1R(t1)k1ωkLei(ωk+Ω)t1sin(kπLx2)|g|ea^k|0+𝒪(λ2).\displaystyle-i\lambda\int\mathrm{d}t_{1}\chi_{1R}(t_{1})\sum_{k}\frac{1}{\sqrt{\omega_{k}L}}e^{i(\omega_{k}+\Omega)t_{1}}\sin(\frac{k\pi}{L}x_{2})\ket{g}\ket{e}\hat{a}_{k}^{\dagger}\ket{0}+\mathcal{O}(\lambda^{2}). (42)

For simplicity of further calculation, let us introduce the following notation:

|ψR=U^1U^2|g|g|0=|gg|0+|ge|ϕgeR+|eg|ϕegR+𝒪(λ2),\displaystyle\ket{\psi_{R}}=\hat{U}_{1}\hat{U}_{2}\ket{g}\ket{g}\ket{0}=\ket{gg}\ket{0}+\ket{ge}\ket{\phi_{ge}^{R}}+\ket{eg}\ket{\phi_{eg}^{R}}+\mathcal{O}(\lambda^{2}), (43)
|ψL=U^2U^1|g|g|0=|gg|0+|ge|ϕgeL+|eg|ϕegL+𝒪(λ2),\displaystyle\ket{\psi_{L}}=\hat{U}_{2}\hat{U}_{1}\ket{g}\ket{g}\ket{0}=\ket{gg}\ket{0}+\ket{ge}\ket{\phi_{ge}^{L}}+\ket{eg}\ket{\phi_{eg}^{L}}+\mathcal{O}(\lambda^{2}), (44)

where: |ϕgeL\ket{\phi_{ge}^{L}}–a state of a field when the left detector interacts first but the right detector is excited, |ϕegL\ket{\phi_{eg}^{L}}–a state of a field when the left detector interacts first and the left detector is excited, |ϕgeR\ket{\phi_{ge}^{R}}–a state of a field when the right detector interacts first and the right detector is excited, |ϕegR\ket{\phi_{eg}^{R}}–a state of a field when the right detector interacts first but the left detector is excited as a consequence of the interaction between atoms and the field.

Supplementary Note 3: Method of finding appropriate parameters

In this section we will find an appropriate set of parameters that orthogonalize |ΦR\ket{\Phi_{R}} and |ΦL\ket{\Phi_{L}} vectors.

|ΦR\displaystyle\ket{\Phi_{R}} =iλk1ωkLdt(χ2R(t)sin(kπLx2)+χ1R(t)sin(kπLx1))ei(ωk+Ω)ta^k|0\displaystyle=-i\lambda\sum_{k}\frac{1}{\sqrt{\omega_{k}L}}\int\mathrm{d}t\left(\chi_{2R}(t)\sin(\frac{k\pi}{L}x_{2})+\chi_{1R}(t)\sin(\frac{k\pi}{L}x_{1})\right)e^{i(\omega_{k}+\Omega)t}\hat{a}_{k}^{\dagger}\ket{0}
=iλkdt(χ2R(t)uk(x2)+χ1R(t)uk(x1))ei(ωk+Ω)ta^k|0\displaystyle=-i\lambda\sum_{k}\int\mathrm{d}t\left(\chi_{2R}(t)u_{k}(x_{2})+\chi_{1R}(t)u_{k}(x_{1})\right)e^{i(\omega_{k}+\Omega)t}\hat{a}_{k}^{\dagger}\ket{0}
=λkeiT(ωk+Ω)1ωk+Ω(eiΔτ(ωk+Ω)uk(x1)+uk(x2))a^k|0\displaystyle=-\lambda\sum_{k}\frac{e^{iT(\omega_{k}+\Omega)}-1}{\omega_{k}+\Omega}\left(e^{i\Delta\tau(\omega_{k}+\Omega)}u_{k}(x_{1})+u_{k}(x_{2})\right)\hat{a}_{k}^{\dagger}\ket{0} (45)
|ΦL\displaystyle\ket{\Phi_{L}} =iλk1ωkLdt(χ2R(t)sin(kπLx1)+χ1R(t)sin(kπLx2))ei(ωk+Ω)ta^k|0\displaystyle=-i\lambda\sum_{k}\frac{1}{\sqrt{\omega_{k}L}}\int\mathrm{d}t\left(\chi_{2R}(t)\sin(\frac{k\pi}{L}x_{1})+\chi_{1R}(t)\sin(\frac{k\pi}{L}x_{2})\right)e^{i(\omega_{k}+\Omega)t}\hat{a}_{k}^{\dagger}\ket{0}
=λkeiT(ωk+Ω)1ωk+Ω(eiΔτ(ωk+Ω)uk(x2)+uk(x1))a^k|0\displaystyle=-\lambda\sum_{k}\frac{e^{iT(\omega_{k}+\Omega)}-1}{\omega_{k}+\Omega}\left(e^{i\Delta\tau(\omega_{k}+\Omega)}u_{k}(x_{2})+u_{k}(x_{1})\right)\hat{a}_{k}^{\dagger}\ket{0} (46)
Refer to caption
Figure 5: The analysis of an orthogonality between |ΦL\ket{\Phi_{L}} and |ΦL\ket{\Phi_{L}}.The absolute value of the scalar product as a function of the energy gap of the detector Ω\Omega.

Where the form of the state |ΦL\ket{\Phi_{L}} we get by changing x1x2x_{1}\longleftrightarrow x_{2}. Now it is easy to see that the scalar product can be written as:

ΦR|ΦL\displaystyle\innerproduct{\Phi_{R}}{\Phi_{L}} λ2k|eiT(ωk+Ω)1|2(ωk+Ω)2(eiΔτ(ωk+Ω)uk(x1)+uk(x2))(eiΔτ(ωk+Ω)uk(x2)+uk(x1))\displaystyle\cong\lambda^{2}\sum_{k}\frac{\absolutevalue{e^{iT(\omega_{k}+\Omega)}-1}^{2}}{\left(\omega_{k}+\Omega\right)^{2}}\left(e^{-i\Delta\tau(\omega_{k}+\Omega)}u_{k}(x_{1})+u_{k}(x_{2})\right)\left(e^{i\Delta\tau(\omega_{k}+\Omega)}u_{k}(x_{2})+u_{k}(x_{1})\right)
2λ2keiΔτ(ωk+Ω)(ωk+Ω)2[1cos(T(ωk+Ω))](uk(x1)+eiΔτ(ωk+Ω)uk(x2))2\displaystyle\cong 2\lambda^{2}\sum_{k}\frac{e^{-i\Delta\tau(\omega_{k}+\Omega)}}{(\omega_{k}+\Omega)^{2}}\left[1-\cos\left(T(\omega_{k}+\Omega)\right)\right]\left(u_{k}(x_{1})+e^{i\Delta\tau(\omega_{k}+\Omega)}u_{k}(x_{2})\right)^{2} (47)

We have to remember that states |ΦL\ket{\Phi_{L}} and |ΦR\ket{\Phi_{R}} were not properly normalized so now we can find the norm:

|ΦR\displaystyle\norm{\ket{\Phi_{R}}} =ΦR|ΦR=λ2k|eiT(ωk+Ω)1|2(ωk+Ω)2|eiΔτ(ωk+Ω)uk(x1)+uk(x2)|2\displaystyle=\sqrt{\innerproduct{\Phi_{R}}{\Phi_{R}}}=\sqrt{\lambda^{2}\sum_{k}\frac{\absolutevalue{e^{iT(\omega_{k}+\Omega)}-1}^{2}}{\left(\omega_{k}+\Omega\right)^{2}}\absolutevalue{e^{i\Delta\tau(\omega_{k}+\Omega)}u_{k}(x_{1})+u_{k}(x_{2})}^{2}}
=2λksin2T(ωk+Ω)2(ωk+Ω)2[uk(x1)2+2uk(x1)uk(x2)cos(Δτ(ωk+Ω))+uk(x2)2]\displaystyle=2\lambda\sqrt{\sum_{k}\frac{\sin^{2}\frac{T(\omega_{k}+\Omega)}{2}}{(\omega_{k}+\Omega)^{2}}\left[u_{k}(x_{1})^{2}+2u_{k}(x_{1})u_{k}(x_{2})\cos\left(\Delta\tau(\omega_{k}+\Omega)\right)+u_{k}(x_{2})^{2}\right]} (48)
|ΦL\displaystyle\norm{\ket{\Phi_{L}}} =ΦL|ΦL=2λksin2T(ωk+Ω)2(ωk+Ω)2[uk(x2)2+2uk(x2)uk(x1)cos(Δτ(ωk+Ω))+uk(x1)2]\displaystyle=\sqrt{\innerproduct{\Phi_{L}}{\Phi_{L}}}=2\lambda\sqrt{\sum_{k}\frac{\sin^{2}\frac{T(\omega_{k}+\Omega)}{2}}{(\omega_{k}+\Omega)^{2}}\left[u_{k}(x_{2})^{2}+2u_{k}(x_{2})u_{k}(x_{1})\cos\left(\Delta\tau(\omega_{k}+\Omega)\right)+u_{k}(x_{1})^{2}\right]} (49)

We can notice that |ΦR=|ΦL\norm{\ket{\Phi_{R}}}=\norm{\ket{\Phi_{L}}} and finally the scalar product has the following form:

ΦR|ΦL=ΦR|ΦL|ΦR|ΦL=keiΔτ(ωk+Ω)(ωk+Ω)2[1cos(T(ωk+Ω))](uk(x1)+eiΔτ(ωk+Ω)uk(x2))22nsin2T(ωn+Ω)2(ωn+Ω)2[un(x2)2+2un(x2)un(x1)cos(Δτ(ωn+Ω))+un(x1)2]\displaystyle\innerproduct{\Phi_{R}}{\Phi_{L}}=\frac{\innerproduct{\Phi_{R}}{\Phi_{L}}}{\norm{\ket{\Phi_{R}}}\norm{\ket{\Phi_{L}}}}=\frac{\sum_{k}\frac{e^{-i\Delta\tau(\omega_{k}+\Omega)}}{(\omega_{k}+\Omega)^{2}}\left[1-\cos\left(T(\omega_{k}+\Omega)\right)\right]\left(u_{k}(x_{1})+e^{i\Delta\tau(\omega_{k}+\Omega)}u_{k}(x_{2})\right)^{2}}{2\sum_{n}\frac{\sin^{2}\frac{T(\omega_{n}+\Omega)}{2}}{(\omega_{n}+\Omega)^{2}}\left[u_{n}(x_{2})^{2}+2u_{n}(x_{2})u_{n}(x_{1})\cos\left(\Delta\tau(\omega_{n}+\Omega)\right)+u_{n}(x_{1})^{2}\right]} (50)

This function can be estimated as a finite sum. Let us denote the upper limits of these two sums as NnN_{n} for the sum over nn and NkN_{k} for the sum over kk.

Refer to caption
Refer to caption
Refer to caption
Figure 6: The analysis of an orthogonality between |ΦL\ket{\Phi_{L}} and |ΦL\ket{\Phi_{L}} for different parameters describing interaction.
An absolute value of the scalar product for two detectors standing at x1=L/4x_{1}=L/4 and x2=3L/4x_{2}=3L/4. The energy gap Ω\Omega is chosen as one of the cavity frequencies. Scalar product approximated as a finite sum of 3030 modes i.e. Nk=Nn=30N_{k}=N_{n}=30.

Fig. 5 shows the absolute value of the scalar product as a function of the energy gap Ω\Omega plotted for different numbers of modes NnN_{n} and NkN_{k}. We can see that we can produce entanglement for very short times Δτ=T=L/10\Delta\tau=T=L/10. Thus, it proves that there is an entanglement between two cavities also for spacelike separated events of interaction.

Fig. 5 look quite random. Let us analyze the scalar product for the case of two detectors standing in the positions x1=L/4x_{1}=L/4 and x2=3L/4x_{2}=3L/4. Fig. 6 show the absolute value of the scalar product for different parameters Δτ\Delta\tau and TΔτT\leq\Delta\tau for different values of the energy gap Ω\Omega. We can notice that there are many parameters minimizing the scalar product between two states.

Based on Fig. 6 we can conjecture that point (Δτ,T)=(3L,2L)\left(\Delta\tau,T\right)=\left(3L,2L\right) is a good candidate for the orthogonalization of the states |ΦR\ket{\Phi_{R}} and |ΦL\ket{\Phi_{L}}. To verify this hypothesis let us consider the following calculation: Let L=1L=1, x1=1/4x_{1}=1/4, x2=3/4x_{2}=3/4, Ω=π\Omega=\pi, Δτ=3\Delta\tau=3, T=2+ϵT=2+\epsilon, where ϵ+\epsilon\in\mathbb{R}_{+} is a small parameter. Then

ΦR|ΦL=ke3ikπk(k+1)2[cos((1+k)(2+ϵ)π)1](sinkπ4e3ikπsin3kπ4)22nsin2(n+1)πϵ22n(n+1)2[(1+2(1)n)cosnπ2+cos3nπ24].\displaystyle\innerproduct{\Phi_{R}}{\Phi_{L}}=-\frac{\sum_{k}\frac{e^{-3ik\pi}}{k(k+1)^{2}}\left[\cos\left((1+k)(2+\epsilon)\pi\right)-1\right]\left(\sin\frac{k\pi}{4}-e^{3ik\pi}\sin\frac{3k\pi}{4}\right)^{2}}{2\sum_{n}\frac{\sin^{2}\frac{(n+1)\pi\epsilon}{2}}{2n(n+1)^{2}}\left[(1+2(-1)^{n})\cos\frac{n\pi}{2}+\cos\frac{3n\pi}{2}-4\right]}. (51)

Each of these sums can be done analytically for arbitrary ϵ\epsilon. Then we expand the numerator and denominator around ϵ=0\epsilon=0, to get:

|ΦR|ΦL|=|12π(2log20log(1i)log(1+i))ϵ2+𝒪(ϵ3)16π(9+log86logπ6logϵ)ϵ2+𝒪(ϵ3)|log89+log86logπϵ.\displaystyle\absolutevalue{\innerproduct{\Phi_{R}}{\Phi_{L}}}=\absolutevalue{\frac{\frac{1}{2\pi}\left(2\log 20\log(1-i)-\log(1+i)\right)\epsilon^{2}+\mathcal{O}(\epsilon^{3})}{\frac{1}{6\pi}\left(9+\log 8-6\log\pi-6\log\epsilon\right)\epsilon^{2}+\mathcal{O}(\epsilon^{3})}}\approx\frac{\log 8}{9+\log 8-6\log\pi\epsilon}. (52)

And we see that:

limϵ0+|ΦR|ΦL|=0\displaystyle\lim_{\epsilon\to 0^{+}}\absolutevalue{\innerproduct{\Phi_{R}}{\Phi_{L}}}=0 (53)

Supplementary Note 4: Indefinite order for spacelike separation of events

Results from the main text of the “Entanglement of the final state” section show that also for spacelike separation of events of interactions Bell’s inequality for temporal order can be violated. This statement seems paradoxical because for spacelike separation of events we cannot define their temporal order - it would depend on the reference frame. Here we look more closely at the compatibility of this fact with the locality of evolution.

Refer to caption
Figure 7: Spacetime diagram of a cavity interacting with two detectors. Operators U^1\hat{U}_{1} and U^3\hat{U}_{3} describes the evolution of the system according to the interaction, and U^2\hat{U}_{2} is an operator of the free evolution that occurs between interactions. FLF_{L} and FRF_{R} are parts of the cavity that can interacts with the detector in some finite time. FF is a middle segment of the cavity that evolves only due to the free evolution operator between interactions.

Let us consider the cavity and two molecules placed near its boundaries. We know that entnaglement at spacelike separation appears for a short time of interaction and small time dilation (see Fig. 4). By definition, in such scenarios information about the interaction with the field can not be transferred from the first to the second detector before the latter interacts with the field. To notice the problem with this situation we can split the whole cavity into three parts (see Fig. 7). Let us denote these as: FLF_{L}, FF, and FRF_{R}. The field FLF_{L} refers only to the part of the cavity nearby the left detector and similarly FRF_{R} is describing just the right part of the cavity. The field FF refers to the middle of the whole cavity. The lengths of each part are chosen to make sure that information about interaction with the left/right molecule can be localized only within the left/right part of the cavity. Because of this division of the cavity into three parts, we can focus on two alternative orders of events. Let us focus on the situation when the left molecule interacts first (LRL\prec R). The whole evolution U^\hat{U} of the system containing two molecules and three parts of the cavity can be presented as U^=U^3U^2U^1\hat{U}=\hat{U}_{3}~{}\hat{U}_{2}~{}\hat{U}_{1}, where:

  • U^1\hat{U}_{1} – interaction of the left molecule with FLF_{L} and free evolution of remaining parts.

  • U^2\hat{U}_{2} – free evolution of the whole cavity.

  • U^3\hat{U}_{3} – interaction of the right molecule with the FRF_{R} and free evolution of remaining parts.

Alternatively, for the right molecule interacting first (RLR\prec L), we would have:

  • U^1\hat{U}_{1} – interaction of the right molecule with the FRF_{R} and free evolution of remaining parts.

  • U^2\hat{U}_{2} – free evolution of the whole cavity.

  • U^3\hat{U}_{3} – interaction of the left molecule with the FLF_{L} and free evolution of remaining parts.

We can notice that these two scenarios differ by not only the order of events but also by the free evolution parts. It means that for such a case we can not argue that violation of Bell’s inequalities imply indefinite order of events because we violate one of the assumptions. We are considering a scenario in which not only the order of events are different but also there is an additional part of the evolution of the system that causes a different final state for each order.