**Abstract**

Time has a fundamentally different character in quantum mechanics and in general relativity. In quantum theory events unfold in a fixed order while in general relativity temporal order is influenced by the distribution of matter. When matter requires a quantum description, temporal order is expected to become non-classical—a scenario beyond the scope of current theories. Here we provide a direct description of such a scenario. We consider a thought experiment with a massive body in a spatial superposition and show how it leads to entanglement of temporal orders between time-like events. This entanglement enables accomplishing a task, violation of a Bell inequality, that is impossible under local classical temporal order; it means that temporal order cannot be described by any pre-defined local variables. A classical notion of a causal structure is therefore untenable in any framework compatible with the basic principles of quantum mechanics and classical general relativity.

**Introduction**

Quantum mechanics forces us to question the view that physical quantities (such as spin, positions or energy) have predefined values: Bell’s theorem shows that if observable quantities were determined by some locally defined classical variables, it would be impossible to accomplish certain tasks—such as the violation of Bell’s inequalities—whereas such tasks are possible according to quantum mechanics1,2 and have been realised in experiments3,4,5,6. However, the causal relations between events remain fixed in quantum theory: whether an event A is in the past, in the future, or space-like separated from another event B is predefined by the location of such events in space-time7,8. In contrast, in general relativity, space-time itself is dynamical: the presence of massive objects affects local clocks and thus causal relations between events defined with respect to them. Nonetheless, the dynamical causal structure of general relativity is still classically predefined: the causal relation between any pair of events is uniquely determined by the distribution of matter-energy degrees of freedom (DOFs) in their past light cone. In other words, causal relations are always determined by local classical variables. This picture is expected to change if we consider quantum states of gravitating DOFs: if a massive system is prepared in a superposition of two distinct states, each yielding an observably different causal structure for future events, would it be possible to observe causal relations that display genuine quantum features?

A main obstacle in the analysis of macroscopic superpositions of gravitating bodies is that, in the absence of a classical space-time manifold, it becomes unclear how to identify space-like surfaces on which quantum states are defined, or global fields of time-like vectors to define time evolution. Indeed, some models even postulate that such superpositions are simply not valid physical states and must decohere (or collapse) fast enough to preserve a classical description of space-time and dynamical laws9,10,11,12,13. A very different mindset underlies various quantum gravity frameworks14—where quantum features of the metric and therefore of the causal relations are indeed expected. However, to date, none of the quantum gravity frameworks has been applied to analyse such an epitomic example as superpositions of space-times with macroscopically distinct causal structures. Therefore, it is unclear whether there exists any phenomenology unequivocally associated with quantum causal structures, nor whether quantum gravity frameworks can circumvent or directly address the objections against superpositions of manifolds. Independently, quantum formalisms have been recently developed to study quantum causal structures at an abstract level in the context of quantum-information processing8,15,16. However, although quantum features of space-time are among the motivations for these studies, no direct link with quantum gravity has yet been established.

This work provides the first direct analysis of quantum causal relations arising from a spatial superposition of a massive object. We show how the temporal order between time-like events can become superposed or even entangled. We further discuss a thought experiment, an admissible albeit remote physical scenario, where these non-classical causal relations arise among physical events. In order to prove their non-classicality, we formulate a Bell-type theorem for temporal order: We define a task that cannot be accomplished if the time order between the events was predetermined by local variables, while the task becomes possible if the events are in a space-time region affected by the gravitational field of a massive object in an appropriate quantum state. Our approach provides a method to directly describe scenarios so far considered to be out of reach for standard theoretical physics. We show explicitly how to overcome the difficulties with describing superpositions of metrics that motivated collapse models. On the other hand, our result is independent of the high-energy completion of any specific quantum gravity framework—we do not assume any new physics, the results are based entirely on well-established, low-energy general relativity and on quantum mechanics. Our results are therefore robust against particular mathematical approaches to quantising gravity, thus providing a benchmark for specific frameworks. Furthermore, the time and energy scale at which entangled temporal order arises is closer than the Planck scale, typically invoked in this context, and is also far remote from the scale given by the decoherence models—which therefore do not preclude quantum features of space-time to arise. Our results thus reveal that both the above approaches are missing crucial intuition and correct physical understanding of the phenomena associated with causal structures at the interface of quantum and gravitational physics. In turn, our work provides a robust method to quantitatively assess these phenomena, helping to build correct physical intuition for quantum causal structures.

**Results**

Dynamical causal structure in general relativity

In classical general relativity, the causal structure is the structure of light cones of the space-time metric17,18. As the matter-energy DOFs determine the metric through Einstein’s equations, the causal structure of a region of space-time is dynamical: it depends on the state of the matter energy in its past light cone. A major obstacle towards a quantum theory of gravity is that it is not clear how to transpose the mathematical notion of causal relations to scenarios where matter DOFs can be in general quantum states, as such scenarios seem to preclude the use of any underlying space-time manifold with respect to which events, light cones and causal relations could be defined. To overcome this obstacle, our approach is to start from a physical understanding of events and their causal relations. Even in classical general relativity a physical event cannot be directly identified with a point on a space-time manifold, a fundamental aspect of the theory captured mathematically by diffeomorphism invariance19. Although it can be debated whether or not space-time points have an intrinsic physical meaning, a natural way to define diffeomorphism-invariant events is to specify them operationally, relative to physical systems; for example, positions and proper times of physical systems used as clocks20. We adopt this notion of events throughout the work. Causal relations are then understood as the possibility to exchange non-faster-than-light signals—or more generally, physical systems—between operationally defined events.

The presence of massive bodies generally alters the relative rates at which clocks tick. For example, in a weak field limit, a clock in a gravitational potential Φ exchanging signals with an identical clock far away from the source of Φ, where the potential effectively vanishes, will appear to tick slower by a factor 1+2Φc2−−−−−−√. In classical physics, this leads to the well-tested time-dilation21,22 and redshift effects23. When the clocks are described as quantum systems, new effects arise from the combination of quantum and general relativistic theories. For a clock in superposition of different distances to the mass, its time-keeping DOFs become entangled to the clock’s position24,25,26. This entanglement implies a universal decoherence mechanism for generic macroscopic systems under time dilation27,28. The regime of low-energy quantum systems in curved space-time can be described within a framework of general relativistic composite quantum particles29. Here we additionally exploit the fact that only the distance between a clock and a mass has physical significance and due to linearity of quantum theory this must hold also for a superposition of different distances. (There is no difference in the relative ticking rates of two clocks whether we think that the clocks are being positioned at different distances—possibly in a superposition—from the mass, or that the mass is positioned at different distances from the clocks30.)

Consider two agents, a and b, with two initially synchronised clocks, each following a fixed world line. A third agent prepares one of two mass configurations, KA≺B or KB≺A, so as to induce time dilation between the clocks of a and b. If configuration KA≺B is prepared, event A—defined by the clock of agent a showing proper time ta = τ*—will be in the past light cone of the event B, which is defined in an analogous way: by the clock of agent b showing proper time tb = τ*. If configuration KB≺A is prepared, event B will be in the past light cone of event A. To keep the world lines of the agents independent of the mass configuration, their laboratories can be embedded in tight enough trapping potentials, that is, much stronger than the gravitational field (which is feasible since our protocol does not require macroscopic source masses, see Methods). In Supplementary Note 4 we discuss other mass configurations, which have the desired effect on temporal order, but for which the agents a, b can remain inertial.

A possible way to realise configuration KA≺B is to place an approximately point-like body of mass M closer to b than to a, see Fig. 1. The light-cone structure of the resulting space-time is fully determined by the metric tensor gμν, for which we adopt the sign convention (−, +, +, +). In isotropic coordinates in the first-order post-Newtonian expansion the metric components are31 g00(r)=−(1+2Φ(r)c2) and gij(r)=δij(1+2Φ(r)c2)−1, i, j = 1, 2, 3, where Φ(r)=−GMr is the gravitational potential and r is the spatial distance between the mass and the event where the metric is evaluated. For an event with a spatial coordinate Ra and the mass at a spatial coordinate rM (where the spatial coordinates are defined, for example, by a far-away agent as in Fig. 1), we have r ≡ |Ra − rM|. Note that we use a common coordinate system to describe the different mass configurations and the associated space-time metrics. Operationally, we can associate such coordinates with the far-away agent, whose local clocks are not affected by the change in the matter distribution. However, this is only a convenient interpretation, we can always think of the coordinates in analogy to gauge fixing—any physical prediction regarding proper times of the clocks and exchange of the signals will not depend on the choice of coordinates.

We consider that a and b remain at fixed coordinate distances from the mass, ra and rb = ra − h, respectively, and find the parameters for which event A ends up in the past light cone of B for KA≺B (and vice versa for KB≺A). An infinitesimal proper time element along a world line at a distance r from the mass is given by dτ(r)=−g00(r)−−−−−−√dt, where t is the coordinate time, and a photon travelling in the radial direction from ra reaches rb after a coordinate time Tc=1c∫rarbdr′−grr(r′)g00(r′)−−−−−−√. Therefore, if the photon is emitted at the local time ta = τ*, it reaches rb when b’s local time is t¯b=−g00(rb)−−−−−−−√(τ∗−g00(ra)√+Tc), assuming that the local clocks are synchronised so that ta = 0 and tb = 0 coincide with the coordinate time t = 0. For

τ∗>Tc−g00(rb)−−−−−−−√1−g00(rb)g00(ra)−−−−−√,

(1)

we have t¯b≤τ∗, which means that there is enough time for a not-faster-than-light signal emitted at event A (defined by ta = τ*) to travel the distance h and reach agent b at event B (defined by tb = τ*). This means that event A is in the causal past of event B as required. For example, for h≪ra condition (1) is satisfied for τ∗>2r2acGM. Configuration KB≺A can be arranged analogously, by placing the mass closer to a than to b. Then, the condition τ∗>2r2bcGM, for h≪rb, ensures that B is in the causal past of A. Note that with the above conditions on τ* the events A and B are always time-like separated, but have different time orders for the two mass configurations—these conditions guarantee that the time order between A and B is swapped in all reference frames.

The example above simply illustrates that in general relativity causal structure is dynamical and depends on the stress-energy tensor of the matter DOFs: preparing different matter distributions on a space-like hypersurface can result in different causal relations between events in its causal future.

Quantum control of temporal order

When A is in the past light cone of B, a physical system can in principle be transferred from A to B. Consider a quantum system S initially prepared in state |ψ〉S, which undergoes a unitary UA at event A (at the space-time location where the clock of agent a marks proper time τ*) and a unitary UB at event B. Such ordered events can therefore result in the following state of S:

∣∣ψ~1⟩S=UBUA|ψ⟩S.

(2)

If B is before A, and S is prepared in the same initial state, the final state of S is

∣∣ψ~2⟩S=UAUB|ψ⟩S.

(3)

A situation can therefore be arranged such that state (2) is produced for configuration KA≺B and (3) is produced for KB≺A. (We ignore a possible additional time evolution between the two events for simplicity.) Different mass configurations can result in different temporal orders of local operations, which holds in quantum as well as in classical theory. Let us make the following assumptions:

(a)

Macroscopically distinguishable states of physical systems can be assigned orthogonal quantum states.

(b)

Gravitational time dilation in a classical limit reduces to that predicted by general relativity.

(c)

The quantum superposition principle holds (regardless of the mass or nature of the involved system).

Even though the above assumptions hold in the standard quantum and general relativistic frameworks, it is not known if a fundamental theory of quantum gravity satisfies them. Our aim is to investigate their consequences for the notion of temporal order.

The coordinates introduced in the previous section define a foliation of space-time into equal-time slices. As long as no horizons are present in any of the considered configurations, such slices define space-like hypersurfaces. With each hypersurface one can associate a Hilbert space, containing the quantum states of interest at the given time. The time coordinate corresponds to the time t in Fig. 1 and is operationally defined as the time measured by the local clock of the far-away agent (not affected by the mass configurations). These quantum states can be understood operationally as states assigned by the far-away agent. However, as discussed in the previous section, such an interpretation is not strictly necessary but is merely a convenient way to define the relevant mathematical objects and to carry out the calculations.

The two mass configurations KA≺B,KB≺A can thus be assigned quantum states |KA≺B⟩M, |KB≺A⟩M. By assumption (a) these states are orthogonal. Since each state individually satisfies the classical limit (mass is sufficiently localised around a single world line), following assumption (b), the system S will evolve as in Eqs. (2) or (3) depending whether the mass is in state |KA≺B⟩M or |KB≺A⟩M, respectively. Finally, by assumption (c), a superposition |K+⟩M:=12√(|KA≺B⟩M+|KB≺A⟩M) is a physically allowed mass configuration, and will yield the following final state of the joint system:

|ψsup⟩MS=12–√(|KA≺B⟩MUBUA|ψ⟩S+|KB≺A⟩MUAUB|ψ⟩S).

(4)

An explicit calculation showing how this state arises is presented in Methods. We note that not only classical gravity but also semi-classical14 and stochastic gravity32 theories would not yield Eq. (4) since these frameworks describe gravitational interactions in terms of classical, possibly stochastic, variables, thus violating assumption (c).

Note that, given a specific physical system used as a clock, it is possible to simulate its time dilation using non-gravitational interactions. For example, an electric field can shift atomic energy levels and thus “time dilate” a clock based on atomic transitions. Therefore, one can produce a state analogous to (4) without using gravity. However, only gravity can alter the relative ordering of events independently of the nature of the systems and interactions used as clocks, due to the universality of time dilation: the preparation and manipulation of the massive object can be carried out without any knowledge of other aspects of the protocol. Such a universality underpins a fundamental distinction between our gravitational protocol and other, non-gravitational, methods to control causal relations between operationally defined events33,34,35,36,37,38,39. (See also Supplementary Note 4 for further discussion.)

Finally, the state (4) is the result of a process wherein the order of operations on a target system (S) is determined by the quantum state of a control system (position of the massive body). Such a process is known as a quantum switch15 and has been studied as a possible quantum-information resource40,41,42,43,44. The state |ψsup〉MS is a superposition of two amplitudes corresponding to different predefined, classical orders between events A and B. Note that, if the control system is discarded, the reduced state of S is

12(|ψ~1⟩⟨ψ~1|S+|ψ~2⟩⟨ψ~2|S),

(5)

which is indistinguishable from a probabilistic mixture of |ψ~1⟩ and |ψ~2⟩. The state in Eq. (5) can be interpreted as arising from events A and B with a classical, albeit unknown, temporal order. Therefore, any protocol aimed at testing operationally quantum features of temporal order necessarily requires a measurement of the control system.

Bell’s theorem for temporal order

The above argument shows that superpositions of massive objects can in principle result in a coherent quantum control of temporal order between events. However, one might question whether such a conclusion has a direct physical meaning or whether it relies on a particular interpretation of state (4). Furthermore, the state assignment is defined in terms of a given coordinate system, while we would like to base our conclusions on coordinate-independent physical events. Since the very meaning of quantum states and measurements might be put into question in the absence of a classical space-time, a proof of non-classical causal relations should not rely on the validity of the quantum formalism. In the following we show that it is possible to probe the nature of temporal order irrespective of the validity of quantum theory. We formulate a theory-independent argument—which does not rely on the quantum framework and provides means to exclude the very possibility of explaining data from a hypothetical experiment in terms of a classical temporal order (which can be stochastic and dynamical) within a broad class of probabilistic theories, not limited to quantum mechanics. Our formulation is analogous to Bell’s theorem for local hidden variables1,2 (see Methods) and we thus refer to the theorem below as Bell’s theorem for temporal order of events. The core of the argument is simple: given a bipartite system prepared in a separable state, it is not possible to violate any bipartite Bell inequality by performing local operations (transformations and measurements) on the two parts, as long as the local operations are applied in a definite order.

The scenario involves a bipartite system with subsystems S1 and S2 and a system M that can influence the temporal order of events. For j = 1, 2, each system Sj undergoes two transformations, TAj and TBj, at space-time events Aj, Bj, respectively. Each system is then measured at an event Cj according to some measurement setting ij, producing a measurement outcome oj. Additionally, M is measured at an event D, space-like separated from both C1 and C2, producing an outcome z, see Fig. 2. We now define the notion of classical order between events:

Definition 1: A set of events is classically ordered if, for each pair of events A and B, there exists a space-like surface and a classical variable λ defined on it that determines the causal relation between A and B: for each given λ, either A≼B (A in the past causal cone of B), B≼A (A in the past causal cone of B) or A||B (A and B space-like separated).

Classically ordered events do not necessarily form a partially ordered set: classical order can be dynamical (the order between two events can depend on some operation performed in the past, i.e. some agent can prepare λ) and stochastic (λ might be distributed according to some probability, and not specified deterministically)45,46.

Bell’s theorem for temporal order. No states, set of transformations and measurements which obey assumptions 1–5 below can result in a violation of the Bell inequalities:

1.

Local state: The initial state ω of S1, S2 and M is separable (as defined in Methods).

2.

Local operations: All transformations performed on the systems are local (as defined in Methods).

3.

Classical order: The events at which operations (transformations and measurements) are performed are classically ordered.

4.

Space-like separation: Events (A1, B1) are space-like separated from events (A2, B2); C1, C2, and D are pair-wise space-like separated.

5.

Free-choice: The measurement choices in the Bell measurement are independent of the rest of the experiment. (This is a standard assumption necessary in Bell-like theorems.)

More formally, let us denote by T=(TA1,TB1,TA2,TB2) the set of all local transformations irrespective of their order. The thesis of the theorem can be rephrased as: the conditional probability

P(o1,o2|i1,i2,z,T,ω)

(6)

produced under assumptions 1–5 does not violate Bell’s inequalities for any value of z. The proof of the theorem is presented in Methods.

Violation of Bell inequalities for temporal order

Here we show how the gravitational quantum control of temporal order from the first section can result in events whose temporal order is entangled: a bipartite quantum system, initially in a product state |ψ1⟩S1|ψ2⟩S2, is sent to two different regions of space such that a1, b1 and c1 only interact with S1, while a2, b2 and c2 only interact with S2. Agents a1, a2 perform, respectively, the unitaries UA1, UA2 at the events A1, A2, while agents b1, b2, perform the unitaries UB1, UB2 at the events B1, B2. Finally, c1 and c2 measure S1 and S2 at events C1 and C2, respectively, see Fig. 3. Assume that a massive system can be prepared in two configurations, KA≺B and KB≺A, such that A1≺B1≺C1 (A1 in the past light cone of B1, etc.) and A2≺B2≺C2 for KA≺B, while B1≺A1≺C1 and B2≺A2≺C2 for KB≺A, and such that the events are space-like separated as per assumption 4, which can always be achieved by having the groups sufficiently separated. If the mass is prepared in superposition 12√(|KA≺B⟩M+|KB≺A⟩M), the joint state of the mass and the systems after the application of the unitaries is

12–√⎛⎝⎜⎜|KA≺B⟩MUB1UA1|ψ1⟩S1UB2UA2|ψ2⟩S2+|KB≺A⟩MUA1UB1|ψ1⟩S1UA2UB2|ψ2⟩S2⎞⎠⎟⎟.

(7)

Agent d at the event D measures the mass in the superposition basis |±⟩=12√(|KA≺B⟩±|KB≺A⟩). Conditioned on the outcome, the joint state of S1 and S2 reads

12–√(UB1UA1|ψ1⟩S1UB2UA2|ψ2⟩S2±UA1UB1|ψ1⟩S1UA2UB2|ψ2⟩S2).

(8)

If the states UB1UA1|ψ1⟩S1, UB2UA2|ψ2⟩S2 are orthogonal to UA1UB1|ψ1⟩S1, UA2UB2|ψ2⟩S2, respectively, then the state (8) is maximally entangled. Local measurements can thus be performed on subsystems S1, S2 whose outcomes will violate Bell inequalities, conditioned on the measurement outcome at D (see Supplementary Note 2 for an example).

The above thought experiment can in principle be realised in a scenario where it is meaningful to argue that assumptions 1, 2 and 4, 5 are satisfied. Violation of the Bell’s inequality would then imply that assumption 3 does not hold, proving non-classicality of temporal order. In order to maximally violate the inequality, the time-dilated clocks of the agents need to decorrelate from the systems Si. In the Methods section we present a particular scenario using photons that satisfies also this requirement. In Supplementary Note 3 we present two concrete examples of our thougth experiment, using as the systems Si polarisation states of photons, depicted in Supplementary Fig. 1, or spatial modes of a quantum field, depicted in Supplementary Fig. 2.

Discussion

The non-classical causal structures discussed in this work arise in a semi-classical, albeit non-perturbative, regime where no explicit quantisation of the gravitational field is needed (which is complementary to the regime of most quantum gravity frameworks14). Our approach shows that general relativity and standard quantum mechanics are sufficient to analyse scenarios involving superpositions of macroscopically different classical backgrounds. Not only is there no tension between the two frameworks, but there is also no ambiguity in the prediction of physical effects that arise: for each probability amplitude, the time-dilation effects introduced by the mass can be treated classically. The considered processes involve a simple superposition of such amplitudes and the final probability amplitude is given by the usual Feynman sum. Note that, even though no explicit quantisation of the metric is used, the amplitudes in the Feynman sum do correspond to macroscopically distinct space-time metrics: this is because each of these amplitudes contains a different causal structure, which determines the metric up to a conformal factor17,18. Quantisation of the metric is therefore implicit in our result, in a similar way as in recently considered witnesses for quantum gravity in interferometric scenarios47,48,49.

A practical realisation of the Bell test for time order would be extremely challenging, even in light of current efforts to prepare superposition states of massive objects and test their gravitational interactions50,51,52,53,54. However, there would be far reaching consequences if a such a test were fundamentally impossible: this would imply that time order, and thus time itself, can be described with a classical parameter even in space-times originating from a quantum state of a massive object—with no need to invoke any other mechanism, such as refs. 9,10,11,12,13, that would decohere these states (see also Supplementary Note 5 for further discussion). On the other hand, since these mechanisms postulate a specific decoherence time of spatial superpositions, one could think that they preclude the preparation of non-classical causal structures. This is not the case: the time required to complete our protocol can be shorter than the decoherence time postulated by these models (see Methods). Thus, contrary to some motivations11,13, these models do not enforce fundamentally classical space-time with a fixed causal structure (i.e. there is a parameter regime where entangled causal structures could form but decoherence postulated by these models is negligible). Finally, classical temporal order could not be excluded also in a scenario where massive bodies can be prepared in quantum states but one (or more) of the assumptions 1, 2, 4, and 5 cannot be satisfied for some fundamental reason. We note that in particular the notion of locality may be fundamentally limited in the context of quantum gravity55,56.

We should note that proof-of-principle realisations of indefinite causal order, analogous to the examples discussed here, have been realised in the laboratory. However, such realisations cannot be interpreted as proofs of non-classical space-time in the sense of general relativity, see Supplementary Note 4 for a discussion of the key differences between the gravitational and other methods for a quantum control of temporal order. The full extent of the relation between gravitational and non-gravitational realisations of quantum causal structures merits an in-depth study on its own.

A crucial aspect of Bell’s theorem for temporal order is that it provides a theory independent result—it applies to any framework where causal relations are described classically, such as classical, semi-classical14 and stochastic gravity32 theories. Moreover, joint validity of the quantum superposition principle and gravitational time dilation, assumptions (a)–(c), suffice for a maximal possible violation of the bound. Therefore, a classical notion of temporal order is untenable in any theory compatible with these basic principles. Finally, the way in which a non-classical causal structure can be engineered exploiting time dilation from a massive body in a quantum state reveals a close connection between the information-theoretic framework of quantum combs/process matrices and joint effects of quantum mechanics and general relativity.

Methods

Quantum gravitational control of temporal order

According to the Einstein equations, a massive object gives rise to a space-time metric gμν, μ, ν = 0, …, 3, which in isotropic coordinates and a post-Newtonian expansion reads31: g00(r)=−(1+2Φ(r)c2), gij(r)=δij(1−2Φ(r)c2), i, j = , 2, 3, where r denotes the distance to the location of the mass. In other words, if a test mass or a clock is positioned at a spatial coordinate Ra as described by a far-away agent (as in Fig. 1) and the massive object is at a coordinate rM, then r = |Ra − rM|, which for clarity we denote below by Ra − rM. It is important to note that the same coordinates describe scenarios where the mass is placed at different locations at a finite distance from rM, as long as it remains far from an asymptotic region so that the spatial and temporal coordinates of the far-away agent remain unaffected (i.e. are those of flat Minkowski space-time). In these coordinates, the Hamiltonian of a clock—a particle with internal DOFs—reads

Ha=−g00(Ra−rM)(Ω2a+c2gij(Ra−rM)PiPj)−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−√,

(9)

(see e.g. refs. 57,58,59) where Pi, i = 1, 2, 3 are the components of the momentum operator, and Ωa is the internal Hamiltonian, describing the local time evolution of the internal DOFs. Note that we can restrict ourself to an effectively one-dimensional scenario, so only one of the spatial coordinates has been kept in the above expression. In the first post-Newtonian expansion and considering that both the mass and the clock follow fixed world lines at constant Ra and rM, respectively, Eq. (9) becomes

Ha≈Ωa(1+Φ(Ra−rM)c2).

(10)

The asymptotic time coordinate t defines space-like hypersurfaces that are independent of the location of the mass and on which one can define states of all the involved systems (the clocks, the target systems and the mass itself) and Hamiltonian (10) describes their time evolution of with respect to t. Due to the interactions between the mass and the clocks—effected by the space-time metric, which contains the potential Φ(Ra − rM)—the time evolution of the clocks depends on their relative distance Ra − rM to the mass. Crucially, by the definition of t and the Hamiltonian our description includes both considered different mass configurations: the mass can be semi-classically localised around a single spatial coordinate r or in superposition of different spatial coordinates and the associated states belong to the same Hilbert space associated with a space-like hypersurface labelled by t. We thus have all the tools to analyse time evolution in the presence of a superposition state of the mass, even though it leads to a quantifiably non-classical causal structure.

With respect to t and the associated foliation of space-time, the evolution of the clock, which at t = 0, is in an internal state |sa(τ0)〉, where τ0 denotes the clock’s proper time at t = 0, reads

e−iΩat(1+Φ(Ra−rM)c2)|Ra⟩|sa(τ0)⟩=|Ra⟩|sa(τ0+τ(Ra−rM,t))⟩,

(11)

where τ(Ra−rM,t):=t(1+Φ(Ra−rM)c2) is the proper time elapsing for the clock at a radial distance |Ra − rM| from the mass when the elapsed coordinate time is t; and for clarity we set ħ = 1.

Before continuing on to the gravitational quantum control, we give an example of an internal Hamiltonian, state, and evolution. Let us take Ωa = E0|0〉〈0| + E1|1〉〈1| and |sa(τ0=0)⟩=12√(|0⟩+|1⟩), which describe, for example, an atom in an equal superposition of some two electronic energy levels |0〉,|1〉 with energies E0, E1, respectively. Under Ha from Eq. (10) internal state |sa(0)〉 from Eq. (11) evolves as

e−iΩat(1+Φ(Ra−rM)c2)|sa(0)⟩=≡12√e−iE0t(1+Φ(Ra−rM)c2)|0⟩+12√e−iE1t(1+Φ(Ra−rM)c2)|1⟩12√e−iE0τ(Ra−rM,t)|0⟩+12√e−iE1τ(Ra−rM,t)|1⟩,

(12)

which is simply |sa(τ(Ra − rM, t))〉.

We now use the above to show how the quantum superposition principle and general relativity lead to the prediction that quantised matter acts as a quantum control of temporal order. To this end, we assume conditions (a)–(c) from the Results section and consider two clocks positioned at RA and RB, respectively. The Hamiltonian of clock a is thus Eq. (10) and fully analogously for b, Hb≈Ωb(1+Φ(Rb−rM)c2). The clocks are initially synchronised with each other and with a clock of the distant agent so that at t0 = 0 both clocks are at τ0 = 0. We further consider a target system, for example, a mode of the electromagnetic field, initially in a state |ψ〉S, on which an operation OA is performed at an event A = (Ra, τa = τ*) and an operation OB at an event B = (Rb, τb = τ*), where τa, τb refer to the proper times of the clock A, B, respectively. We effectively represent these operations as OA=δ(τa−τ∗,r−Ra)OA, where δ(τA − τ*, r − Ra) is a Dirac delta distribution and OA is an operator (e.g. describing rotation of the polarisation of an electromagnetic field mode by a particular half-wave plate) independent of time and location. The total Hamiltonian reads

Htot=Ha+Hb+OA+OB,

(13)

which for simplicity assumes trivial time evolution of the mass and of the target system between the application of the operations. We furthermore consider the following initial (at t0 = 0) state of the mass, clocks and the target system:

|ψ(0)⟩MSab=|Ra⟩|Rb⟩|sa(τ0=0)⟩|sb(τ0=0)⟩|ψ⟩S(|rL⟩M+|rR⟩M),

(14)

where positions rL, rR of the mass refer to the configurations in the left and the right panel of Fig. 1, respectively, that is, they realise configurations KA≺B and KB≺A: for |rL〉 the mass is at a distance ra = rL − Ra from clock a and at rb = ra − h from b, while for |rR〉 the relative distances are swapped and the mass is at a distance ra − h from a and at ra from b. After coordinate time t such that τ(ra, t) > τ* (where τ∗>2r2bcGM, see main text) the state evolves to

|ψ(t)⟩MSab=|Ra⟩|Rb⟩(|sa(τ(ra,t))⟩|sb(τ(ra−h,t))⟩e−iOBe−iOA|ψ⟩|rL⟩M+|sa(τ(ra−h,t))⟩|sb(τ(ra,t))⟩e−iOAe−iOB|ψ⟩S|rR⟩M).

(15)

The order of applying unitary transformations UA=e−iOA and UB=e−iOB to the target system is controlled by the position of the mass, which due to time-dilation changes causal relations between events A and B. Swapping the mass distribution: |rL〉 → |rR〉, |rR〉 → |rL〉 and letting the state evolve for another time interval t results in the final state where the clocks become synchronised again

|ψ(t)⟩MSab=|Ra⟩|Rb⟩|sa(τf)⟩|sb(τf)⟩(UBUA|ψ⟩S|rR⟩M+UAUB|ψ⟩S|rL⟩M),

(16)

where τf = τ(ra, t) + τ(ra − h, t). Measuring the mass in a superposition basis |rL〉M ± |rR〉M prepares the target system in the corresponding superposition state UBUA|ψ〉S ± UAUB|ψ〉S.

The above example demonstrates that under very conservative assumptions a spatial superposition of a mass generates a quantum-controlled application of unitary operations. More fundamentally, this effect stems from the superposition of different causal structures associated with the superposed states of the mass.

Proof of Bell’s theorem for temporal order

Bell’s theorem in general asserts that, under certain assumptions, the correlations between the outcomes of independent measurements on two subsystems must satisfy a class of inequalities. The two measuring parties are referred to as Alice and Bob. In every experimental run, each of them measures one of two properties of the subsystem they receive. For each of the properties, one of two outcomes is obtained, for convenience chosen to be ±1. Bell’s inequalities follow from the conjunction of the following assumptions: (1) measurement results are determined by properties that exist prior to and independent of the experiment (hidden variables); (2) results obtained at one location are independent of any measurements or actions performed at space-like separation (locality); (3) any process that leads to the choice of which measurement will be carried out is independent from other processes in the experiment (free choice). The outcomes of Alice A(i, λ) and Bob B(i, λ) thus only depend on their own choice of setting, index i, and on the property of the system, variable λ. The correlation between outcomes A(i, λ) and B(i, λ) for the measurement choices i, j is described by E(Ai,Bj)=∫dλP(λ)A(i,λ)B(j,λ), where P(λ) is the probability distribution over the properties of the systems. It is straightforward to check that one possible inequality satisfied by the correlations E(Ai, Bj) is the so-called Clauser–Horne–Shimony–Holt inequality: |E(A1, B1) + E(A1, B2) + E(A2, B1) − E(A2, B2)| ≤ 2. Crucially, quantum theory allows for the left-hand side of this inequality to reach a value >2, and experimental measurements of this (and other inequalities) have confirmed such a violation3,4,5,6. The significance of the violations of Bell’s inequalities is in showing that neither nature nor quantum mechanics obey all three assumptions mentioned above.

The assumption of classical order is sufficient to derive Causal Inequalities16,60: tasks that, without any further assumptions, cannot be performed on a classical causal structure. However, it is not possible to violate causal inequalities using quantum control of order45,61, this is why additional assumptions were required in the present context. It is an open question whether a gravitational implementation of a scenario that does allow for a violation of causal inequalities is possible.

The theorem we have formulated is theory independent, but not fully device-independent, as it requires the notions of a physical state and a physical transformation (in addition to the measured probability distributions), which we introduce below and then proceed to the proof. Discussion of the present work in the context of the theory-dependent framework of causally non-separable quantum processes16,45,61 and the fully theory- and device-independent approach of causal inequalities16,60 is presented in Supplementary Note 1.

We consider a sufficiently broad framework to describe physical systems that can undergo transformations and measurements, similar to generalised probabilistic theories62,63,64. This framework is more general than quantum or classical theory and we thus need to define key notions required in the proof. In this framework, a state ω is a complete specification of the probabilities P(o|i, ω) for observing outcome o given that a measurement with setting i is performed on the system. We are interested in situations where a system can be split up in subsystems, say S1 and S2, with space-like separated agents performing independent operations on S1 and S2. We say ω is a product state, and write ω = ω1 ⊗ ω2, if probabilities for local measurements factorise as P(o1, o2|i1, i2, ω) = P(o1|i1, ω1)P(o2|i2, ω2). If state ωf1 is prepared for system S1 and state ωf2 is prepared for system S2, according to a probability distribution P(f) for some variable f, we write ω=∫dfP(f)ωf1⊗ωf2 and say the state is separable. Probabilities are then given by the corresponding mixture: P(o1,o2|i1,i2,ω)=∫dfP(o1|i1,ωf1)P(o2|i2,ωf2)P(f). Note that for such a decomposition Bell inequalities cannot be violated1,65.

A physical transformation of the system is represented by a function ω↦T(ω). To make our arguments precise we need a notion of local transformations, namely, realised at the time and location defined by a local clock. If S1 is the subsystem on which a local transformation T1 acts, and S2 labels the DOFs space-like separated from T1, then, by definition, T1 transforms product states as ω1⊗ω2↦T1(ω1)⊗ω2 and separable states by convex extension. How local operations act on general, non-separable states can depend on the particular physical theory; however, action on separable states will suffice for our purposes. We further need to define how the local transformations combine. This depends on their relative spatio-temporal locations: if transformations T1, T2 are space-like separated they combine as (T1 ⊗ T2)(ω1 ⊗ ω2) = T1(ω1) ⊗ T2(ω2), which follows from the definition above; if T1 is in the future of T2, we define their combination as T1 ο T2(ω) = T1(T2(ω)). (For simplicity, we omit possible additional transformations taking place between the specified events, as they are of no consequence for our argument.)

Proof Assumption (1) says that there is a random variable f determining the local states ωf1, ωf2 of systems S1, S2, respectively. Assumption (3) says there is a random variable λ that determines the order of events. In general, the two variables can be correlated by some joint probability distribution P(λ, f). By assumption (4), events labelled A1, B1 are space-like separated from events A2, B2 and the order between events within each set (Aj, Bj), j = 1, 2 can be defined by a permutation σj. Most generally, there is a probability P(σj|λ) that the permutation σj is realised for a given λ. By assumption (2), for each given order the system undergoes a transformation Tσ1⊗Tσ2, where Tσ1 is the transformation obtained by composing TA1 and TB1 in the order corresponding to the permutation σ1 and similarly for Tσ2. (For example, if σ1 corresponds to the order A1≺B1, then Tσ1=TB1∘TA1.) Furthermore, at event D an outcome z is obtained with a probability P(z|λ, f, σ1, σ2). Finally, using assumption (1), we write the probabilities for all outcomes as

P(o1,o2,z|i1,i2,T,ω)=∑σ1σ2∫dλdfP(o1|i1,Tσ1(ωf1))P(o2|i2,Tσ2(ωf2))P(σ1|λ)P(σ2|λ)P(z|λ,f,σ1,σ2)P(λ,f).

(17)

A simple Bayesian inversion P(σ1|λ)P(σ2|λ)P(z|λ, f, σ1, σ2)P(λ, f) = P(λ, f, σ1, σ2|z)P(z), where we used P(σj|λ) = P(σj|λ, f), gives the desired probabilities

P(o1,o2|i1,i2,z,T,ω)==∑σ1σ2∫dλdfP(o1|i1,Tσ1(ωf1))P(o2|i2,Tσ2(ωf2))P(λ,f,σ1,σ2|z)∫df~P(o1|i1,Tσ1)P(o2|i2,Tσ2)P(f~|z),

(18)

where f~ is a short-hand for the variables λ, f, σ1, and σ2. The above probability distribution satisfies the hypothesis of Bell’s theorem and thus cannot violate any Bell inequality.

Exemplary scenario realising Bell test for temporal order of events

The protocol allowing for the violation of Bell’s inequalities for temporal order exploits correlations between the clocks of the agents a1, b1 and the agents a2, b2, created due to time dilation induced by the mass. It should be noted that the protocol allows maximal violation of the Bell inequality if the joint state of the systems S1 and S2 is pure (and maximally entangled) when the Bell measurements are realised. Thus, for a maximal violation, the clocks need to decorrelate from the mass after the application of the unitaries. Below we sketch a scenario that can achieve this.

The space-time arrangement of the mass and the agents in this example is presented in Fig. 4. It can be realised in one spatial dimension: agents acting on the system S1 are located at distance h from each other, and the mass is placed at distance r (configuration KB≺A) or r + L (configuration KA≺B) from agent a1. Agents acting on system S2 are placed symmetrically on the opposite side of the mass, such that the mass is at a distance r + L from a2 in configuration KB≺A and r in configuration KA≺B. Here, events Aj are defined by the local time τa that differs from the local time τb defining Bj, j = 1, 2. In such a case, even though the mass is always closer to aj than to bj, the two mass configurations can lead to different event orders—as they induce different relative time dilations. (Equivalently, one can introduce an initial offset in the synchronisation of the clocks.) Note that the time orders between the two groups are here “anti-correlated”: A1≺B1 and B2≺A2 for KA≺B, and vice versa for KB≺A. Since otherwise the scenario is the same for S1 and S2, we focus on the operations performed on S1. The key observation is that swapping the mass distribution, as depicted in Fig. 4, will eventually disentangle the clocks from the mass, and since the clocks must be suitably time dilated when the operations are performed, the operations must not take place in the future light cone of the swapped mass state.

The proper time τa that has to elapse for the clock of a1 such that the order of events is A1≺B1 for |KA≺B⟩ and B1≺A1 for |KB≺A⟩ for the present case reads

τa=−g00(r)−−−−−−√Tc(r,h)+Tc(r+L,h)g00(r+L+h)g00(r+h)−−−−−−−−√1−g00(r)g00(r+L+h)g00(r+h)g00(r+L)−−−−−−−−−−−√,

(19)

where Tc(r, L/2) is the coordinate travel time of light between radial distances r and r + L/2 from the mass. The coordinate time corresponding to τa is Ta=τa/−g00(r)−−−−−−√. The proper time of event B1 is then defined as:

τb=−g00(r+L+h)−−−−−−−−−−−−−√(τa−g00(r+L)−−−−−−−−−−√+Tc(r+L,h)).

(20)

It can directly be checked that when the mass is placed in configuration KA≺B—at a distance r + L from a1—the event A1 defined by local clock of a1 showing proper time τa from Eq. (18) is in the past light cone of event B1, which is defined by the local clock of b1 showing proper time τb from Eq. (19). When the mass is placed in configuration KB≺A, event B1 ends up in the past of the event A1. The coordinate time required for the application of the operations can be estimated as twice the travel time of light between the agents, To = 2Tc(r + L/2, h).

The world lines of the mass can be arranged such that: (a) the mass is moving slow so that the two amplitudes of the mass are swapped in a time interval longer than To; (b) during the application of the operations the distance of each agent to the mass is approximately the same for both mass configurations (as in Fig. 4). The first guarantees that there is enough time to apply the operations after the clocks get correlated, the second—that the slow-down of light in curved space-time, the Shapiro delay66,67, can be neglected.

The coordinate-time duration of the entire protocol can be estimated as Tp = 2Ta + 4L/2c, where L/2c is the minimal time required to put the mass in superposition of amplitudes separated by the distance L/2. Taking as an example M ~ 0.1 μg, L = h ~ 0.1 μm, r ~ 1 fm, the protocol in Fig. 4 takes Tp ~10 h. Furthermore, we note that a quantum treatment of the local clocks is central to our protocol since the application of the operations on the target systems is conditioned on the states of the clocks. The time-energy uncertainty68,69 thus poses a limitation to a single-shot precision with which space-time events can be defined with physical clocks. The optimal clock state in this context—evolving the fastest—is a balanced superposition of energy eigenstates; for an energy gap ħ ⋅ 2πνc, where νc is the clock frequency, the smallest time that can be resolved by a single quantum system is the so-called orthogonalisation time70,71,72 t⊥ = 1/2νc. For the values of parameters quoted above, the coordinate-time difference between the superposed locations of the events Ai, i = 1, 2 is ~10−15 s, and we thus need a system with frequency νc ≥1015 Hz such as a clock based on optical transitions in ytterbium73 or mercury74, which both give t⊥ ~10−16 s. While this ideal limit is not reached with practical systems, the resolution of current atomic clocks based on such atoms far exceeds this theoretical bound due to averaging over many atoms, with 2.5 × 10−19 uncertainty of the clock frequency recently demonstrated in ref. 75. We further note that by using n entangled atoms, the orthogonalisation time of the entire system becomes t⊥/n and can thus be even a few orders of magnitude smaller76 than required. Finally, such atoms have masses ~10−25 kg and their back action on the metric produced by M ~ 10−7 kg would thus be negligible. Since the mass difference between the atom in the two involved energy levels is 2πħνc/c2 ~10−35 kg also quantum effects from the clocks’ mutual gravitational interactions58 can be neglected.

We conclude that it is in principle possible to achieve the required entanglement of orders, swap the mass distribution so as to finally disentangle the clocks from the mass, and satisfy the locality conditions on the events. Although a direct experiment in such a regime is not practical, the above example surprisingly shows that the regime where entangled temporal order arises is in no way related to the Planck scale. It is usually assumed that the Planckscale marks the regime where quantum gravity effects become relevant (first discussed in this context by Bronstein77), but this is not the case for the superposition of temporal order. In terms of a potential experiment, one could also take a different (theory-dependent) approach and explore possible witnesses of entangled temporal order61, in analogy to witnesses of entanglement in quantum-information theory78. A witness would probe the quantum nature of temporal order indirectly and under further assumptions, but in a relaxed parameter range. Such an approach may lead to more feasible experiments, which will be explored in a future study.

A spatial superposition state of a mass such as used in our protocol is postulated to decohere in various gravity-inspired collapse models9,10,11,12,13 (which thus violate assumption (c) in the first section). However, even if endorsed, these models do not immediately preclude realisation of our protocol: the decoherence time scale in those models is the Diosi–Penrose time10,11 TDP=2δ3ℏG(ML)2, where δ is a free parameter. For every value of δ one can find the mass and the relevant distances (M, r, L, h) so that the duration of our entire protocol is shorter than TDP. For example, following the recent ref. 79 and taking δ = 10−7 m, for r = 1010RSch, L = 5r, h = r and M = 1 g, where RSch ≈10−30 m, the protocol from Fig. 4 takes Tp ≈7 × 10−18 s, while TDP ≈0.5 s. Taking instead the originally proposed value δ = 10−15 m10, the desired regime is achieved, for example, for M = 10−7 kg, r = 107RSch, L = 5 × 105r, h = 105r; with Tp ~10−23 s and TDP ~10−13 s. Thus, the above models in principle still allow for events with entangled temporal order, and do not enforce the classicality of the causal structure of space time.

Data availability

The data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request.

References

1.

Bell, J. S. On the Einstein–Poldolsky–Rosen paradox. Physics 1, 195–200 (1964).

Article

Google Scholar

2.

Clauser, J. F., Horne, M. A., Shimony, A. & Holt, R. A. Proposed experiment to test local hidden-variable theories. Phys. Rev. Lett. 23, 880–884 (1969).

ADS

Article

Google Scholar

3.

Freedman, S. J. & Clauser, J. F. Experimental test of local hidden-variable theories. Phys. Rev. Lett. 28, 938–941 (1972).

ADS

CAS

Article

Google Scholar

4.

Hensen, B. et al. Loophole-free Bell inequality violation using electron spins separated by 1.3 kilometres. Nature 526, 682–686 (2015).

ADS

CAS

Article

Google Scholar

5.

Giustina, M. et al. Significant-loophole-free test of bell’s theorem with entangled photons. Phys. Rev. Lett. 115, 250401 (2015).

ADS

Article

Google Scholar

6.

Shalm, L. K. et al. Strong loophole-free test of local realism. Phys. Rev. Lett. 115, 250402 (2015).

ADS

Article

Google Scholar

7.

Hardy, L. Probability theories with dynamic causal structure: a new framework for quantum gravity. Preprint at https://arxiv.org/abs/gr-qc/0509120 (2005).

8.

Hardy, L. Towards quantum gravity: a framework for probabilistic theories with non-fixed causal structure. J. Phys. A 40, 3081–3099 (2007).

ADS

MathSciNet

Article

Google Scholar

9.

Karolyhazy, F. Gravitation and quantum mechanics of macroscopic objects. Il Nuovo Cim. A 42, 390–402 (1966).

ADS

Article

Google Scholar

10.

Diosi, L. Models for universal reduction of macroscopic quantum fluctuations. Phys. Rev. A 40, 1165–1174 (1989).

ADS

CAS

Article

Google Scholar

11.

Penrose, R. On gravity’s role in quantum state reduction. Gen. Relativ. Gravit. 28, 581–600 (1996).

ADS

MathSciNet

Article

Google Scholar

12.

Stamp, P. C. E. Environmental decoherence versus intrinsic decoherence. Philos. Trans. R. Soc. Ser. A 370, 4429–4453 (2012).

ADS

CAS

Article

Google Scholar

13.

Penrose, R. On the gravitization of quantum mechanics 1: quantum state reduction. Found. Phys. 44, 557–575 (2014).

ADS

MathSciNet

Article

Google Scholar

14.

Kiefer, C. Quantum Gravity 3rd edn. International Series of Monographs on Physics (OUP, Oxford, 2012).

15.

Chiribella, G., D’Ariano, G. M., Perinotti, P. & Valiron, B. Quantum computations without definite causal structure. Phys. Rev. A 88, 022318 (2013).

ADS

Article

Google Scholar

16.

Oreshkov, O., Costa, F. M. & Brukner, C. Quantum correlations with no causal order. Nat. Commun. 3, 1092 (2012).

ADS

Article

Google Scholar

17.

Hawking, S. W., King, A. R. & McCarthy, P. J. A new topology for curved space?time which incorporates the causal, differential, and conformal structures. J. Math. Phys. 17, 174–181 (1976).

ADS

MathSciNet

Article

Google Scholar

18.

Malament, D. B. The class of continuous timelike curves determines the topology of spacetime. J. Math. Phys. 18, 1399–1404 (1977).

ADS

MathSciNet

Article

Google Scholar

19.

Stachel, J. The hole argument and some physical and philosophical implications. Living Rev. Relativ. 17, 1 (2014).

ADS

Article

Google Scholar

20.

Rovelli, C. What is observable in classical and quantum gravity? Class. Quantum Gravity 8, 297 (1991).

ADS

MathSciNet

Article

Google Scholar

21.

Hafele, J. C. & Keating, R. E. Around-the-world atomic clocks: observed relativistic time gains. Science 177, 168–170 (1972).

ADS

CAS

Article

Google Scholar

22.

Chou, C.-W., Hume, D., Rosenband, T. & Wineland, D. Optical clocks and relativity. Science 329, 1630–1633 (2010).

ADS

CAS

Article

Google Scholar

23.

Pound, R. & Rebka, G. Apparent weight of photons. Phys. Rev. Lett. 4, 337–341 (1960).

ADS

Article

Google Scholar

24.

Zych, M., Costa, F., Pikovski, I. & Brukner, C. Quantum interferometric visibility as a witness of general relativistic proper time. Nat. Commun. 2, 505 (2011).

ADS

Article

Google Scholar

25.

Zych, M., Costa, F., Pikovski, I., Ralph, T. C. & Brukner, C. General relativistic effects in quantum interference of photons. Class. Quantum Gravity 29, 224010 (2012).

ADS

MathSciNet

Article

Google Scholar

26.

Zych, M., Pikovski, I., Costa, F. & Brukner, Č. General relativistic effects in quantum interference of “clocks”. J. Phys. 723, 012044 (2016).

Google Scholar

27.

Pikovski, I., Zych, M., Costa, F. & Brukner, Č. Universal decoherence due to gravitational time dilation. Nat. Phys. 11, 668–672 (2015).

CAS

Article

Google Scholar

28.

Pikovski, I., Zych, M., Costa, F. & Brukner, Č. Time dilation in quantum systems and decoherence. N. J. Phys. 19, 025011 (2017).

Article

Google Scholar

29.

Zych, M. Quantum Systems under Gravitational Time Dilation. Springer theses (Springer International Publishing, Switzerland, 2017).

Google Scholar

30.

Zych, M., Costa, F. & Ralph, T. C. Relativity of quantum superpositions. Preprint at https://arxiv.org/abs/1809.04999 (2018).

31.

Weinberg, S. Gravitation and Cosmology: Principle and Applications of General Theory of Relativity (Wiley, New York, 1972).

Google Scholar

32.

Hu, B. L. & Verdaguer, E. Stochastic gravity: a primer with applications. Class. Quantum Gravity 20, R1 (2003).

ADS

MathSciNet

Article

Google Scholar

33.

Procopio, L. M. et al. Experimental superposition of orders of quantum gates. Nat. Commun. 6, 7913 (2015).

ADS

CAS

Article

Google Scholar

34.

Rubino, G. et al. Experimental verification of an indefinite causal order. Sci. Adv. 3, https://doi.org/10.1126/sciadv.1602589 (2017).

ADS

Article

Google Scholar

35.

Rubino, G. et al. Experimental entanglement of temporal orders. Preprint at http://arxiv.org/abs/1712.06884 (2017).

36.

Goswami, K. et al. Indefinite causal order in a quantum switch. Phys. Rev. Lett. 121, 090503 (2018).

ADS

CAS

Article

Google Scholar

37.

Goswami, K., Romero, J. & White, A. Communicating via ignorance. Preprint at https://arxiv.org/abs/1807.07383 (2018).

38.

Wei, K. et al. Experimental quantum switching for exponentially superior quantum communication complexity. Phys. Rev. Lett. 122, 120504 (2019).

ADS

CAS

Article

Google Scholar

39.

Guo, Y. et al. Experimental investigating communication in a superposition of causal orders. Preprint at http://arxiv.org/abs/1811.07526 (2018).

40.

Chiribella, G. Perfect discrimination of no-signalling channels via quantum superposition of causal structures. Phys. Rev. A 86, 040301 (2012).

ADS

Article

Google Scholar

41.

Colnaghi, T., D’Ariano, G. M., Facchini, S. & Perinotti, P. Quantum computation with programmable connections between gates. Phys. Lett. A 376, 2940–2943 (2012).

ADS

CAS

Article

Google Scholar

42.

Araújo, M., Costa, F. & Brukner, C. Computational advantage from quantum-controlled ordering of gates. Phys. Rev. Lett. 113, 250402 (2014).

ADS

Article

Google Scholar

43.

Feix, A., Araújo, M. & Brukner, Č. Quantum superposition of the order of parties as a communication resource. Phys. Rev. A 92, 052326 (2015).

ADS

Article

Google Scholar

44.

Guérin, P. A., Feix, A., Araújo, M. & Brukner, Č. Exponential communication complexity advantage from quantum superposition of the direction of communication. Phys. Rev. Lett. 117, 100502 (2016).

ADS

Article

Google Scholar

45.

Oreshkov, O. & Giarmatzi, C. Causal and causally separable processes. N.J. Phys. 18, 093020 (2016).

Article

Google Scholar

46.

Abbott, A. A., Giarmatzi, C., Costa, F. & Branciard, C. Multipartite causal correlations: polytopes and inequalities. Phys. Rev. A 94, 032131 (2016).

ADS

Article

Google Scholar

47.

Bose, S. et al. Spin entanglement witness for quantum gravity. Phys. Rev. Lett. 119, 240401 (2017).

ADS

MathSciNet

Article

Google Scholar

48.

Marletto, C. & Vedral, V. Gravitationally induced entanglement between two massive particles is sufficient evidence of quantum effects in gravity. Phys. Rev. Lett. 119, 240402 (2017).

ADS

CAS

Article

Google Scholar

49.

Belenchia, A. et al. Quantum superposition of massive objects and the quantization of gravity. Phys. Rev. D 98, 126009 (2018).

ADS

CAS

Article

Google Scholar

50.

Bose, S., Jacobs, K. & Knight, P. L. Scheme to probe the decoherence of a macroscopic object. Phys. Rev. A 59, 3204–3210 (1999).

ADS

CAS

Article

Google Scholar

51.

Marshall, W., Simon, C., Penrose, R. & Bouwmeester, D. Towards quantum superpositions of a mirror. Phys. Rev. Lett. 91, 130401 (2003).

ADS

MathSciNet

Article

Google Scholar

52.

Kleckner, D. et al. Creating and verifying a quantum superposition in a micro-optomechanical system. N. J. Phys. 10, 095020 (2008).

Article

Google Scholar

53.

Wan, C. et al. Free nano-object ramsey interferometry for large quantum superpositions. Phys. Rev. Lett. 117, 143003 (2016).

ADS

CAS

Article

Google Scholar

54.

Schmöle, J., Dragosits, M., Hepach, H. & Aspelmeyer, M. A micromechanical proof-of-principle experiment for measuring the gravitational force of milligram masses. Class. Quantum Gravity 33, 125031 (2016).

ADS

Article

Google Scholar

55.

Giddings, S. B., Marolf, D. & Hartle, J. B. Observables in effective gravity. Phys. Rev. D 74, 064018 (2006).

ADS

MathSciNet

Article

Google Scholar

56.

Donnelly, W. & Giddings, S. B. Diffeomorphism-invariant observables and their nonlocal algebra. Phys. Rev. D 93, 024030 (2016).

ADS

MathSciNet

Article

Google Scholar

57.

Zych, M. & Brukner, Č. Quantum formulation of the Einstein equivalence principle. Nat. Phys. 14, 1027–1031 (2018).

CAS

Article

Google Scholar

58.

Castro Ruiz, E., Giacomini, F. & Brukner, Č. Entanglement of quantum clocks through gravity. Proc. Natl. Acad. Sci. USA 114, E2303–E2309 (2017).

ADS

MathSciNet

Article

Google Scholar

59.

Zych, M., Rudnicki, L. & Pikovski, I. Gravitational mass of composite systems. Phys. Rev. D 99, 104029 (2019).

ADS

CAS

Article

Google Scholar

60.

Branciard, C., Araújo, M., Feix, A., Costa, F. & Brukner, Č. The simplest causal inequalities and their violation. N. J. Phys. 18, 013008 (2016).

MathSciNet

Article

Google Scholar

61.

Araújo, M. et al. Witnessing causal nonseparability. N. J. Phys. 17, 102001 (2015).

Article

Google Scholar

62.

Hardy, L. Quantum theory from five reasonable axioms. Preprint at https://arxiv.org/abs/quant-ph/0101012 (2001).

63.

Barrett, J. Information processing in generalized probabilistic theories. Phys. Rev. A 75, 032304 (2007).

ADS

Article

Google Scholar

64.

Chiribella, G., D’Ariano, G. M. & Perinotti, P. Probabilistic theories with purification. Phys. Rev. A 81, 062348 (2010).

ADS

Article

Google Scholar

65.

Fine, A. Hidden variables, joint probability, and the bell inequalities. Phys. Rev. Lett. 48, 291–295 (1982).

ADS

MathSciNet

Article

Google Scholar

66.

Shapiro, I. I. Fourth test of general relativity. Phys. Rev. Lett. 13, 789–791 (1964).

ADS

MathSciNet

Article

Google Scholar

67.

Shapiro, I. I. et al. Fourth test of general relativity: new radar result. Phys. Rev. Lett. 26, 1132–1135 (1971).

ADS

Article

Google Scholar

68.

Fleming, G. A unitarity bound on the evolution of nonstationary states. Il Nuovo Cim. A 16, 232–240 (1973).

ADS

Article

Google Scholar

69.

Mandelstam, L. & Tamm, I. The Uncertainty Relation Between Energy and Time in Non-relativistic Quantum Mechanics, in Selected Papers, 115–123 (Springer, Berlin Heidelberg, 1991).

70.

Margolus, N. & Levitin, L. B. The maximum speed of dynamical evolution. Phys. D 120, 188–195 (1998).

Article

Google Scholar

71.

Kosiński, P. & Zych, M. Elementary proof of the bound on the speed of quantum evolution. Phys. Rev. A 73, 024303 (2006).

ADS

Article

Google Scholar

72.

Zieliński, B. & Zych, M. Generalization of the margolus-levitin bound. Phys. Rev. A 74, 034301 (2006).

ADS

Article

Google Scholar

73.

Pizzocaro, M. et al. Absolute frequency measurement of the 1 S 0−3 P 0 transition of 171Yb. Metrologia 54, 102–112 (2017).

ADS

CAS

Article

Google Scholar

74.

Hoyt, C. W. et al. Observation and absolute frequency measurements of the 1 S 0−3 P 0 optical clock transition in neutral ytterbium. Phys. Rev. Lett. 95, 083003 (2005).

ADS

CAS

Article

Google Scholar

75.

Marti, G. E. et al. Imaging optical frequencies with 100 μHz precision and 1.1 μm resolution. Phys. Rev. Lett. 120, 103201 (2018).

ADS

CAS

Article

Google Scholar

76.

Kómár, P. et al. Quantum network of atom clocks: a possible implementation with neutral atoms. Phys. Rev. Lett. 117, 060506 (2016).

ADS

Article

Google Scholar

77.

Gorelik, G. in Studies in the History of General Relativity (eds Eisenstaedt, J. & Kox, A. J.) Vol. 4, 364–379 (Birkhaeuser, Boston, 1992).

78.

Terhal, B. M. Bell inequalities and the separability criterion. Phys. Lett. A 271, 319–326 (2000).

ADS

MathSciNet

CAS

Article

Google Scholar

79.

Bahrami, M., Smirne, A. & Bassi, A. Role of gravity in the collapse of a wave function: a probe into the Diósi–Penrose model. Phys. Rev. A 90, 062105 (2014).

ADS

Article

Google Scholar

Download references

Acknowledgements

We thank G. Chiribella, G. Milburn, H. Wiseman and M. Vojinovic for feedback. M.Z. and F.C. acknowledge support through the Australian Research Council (ARC) Centre of Excellence for Engineered Quantum Systems (CE 110001013), Discovery Early Career Researcher Awards DE180101443, DE170100712 and the Templeton World Charity Foundation (TWCF 0064/AB38). I.P. acknowledges support of the NSF through a grant to ITAMP and the Branco Weiss Fellowship—Society in Science, administered by the ETH Zürich. Č.B. acknowledges the support of the Austrian Science Fund (FWF) through the Doctoral Programme CoQuS, the project I-2526-N27 and I-2906, the research platform TURIS and the ÖAW Innovationsfond “Quantum Regime of Gravitational Source Masses”. This publication was made possible through the support of a grant from the John Templeton Foundation and from the Foundational Questions Institute (FQXi) Fund. The opinions expressed in this publication are those of the authors and do not necessarily reflect the views of the John Templeton Foundation. F.C. and M.Z. acknowledge the traditional owners of the land on which the University of Queensland is situated, the Turrbal and Jagera people.

Author information

Affiliations

Centre for Engineered Quantum Systems, School of Mathematics and Physics, The University of Queensland, St. Lucia, QLD, 4072, Australia

Magdalena Zych & Fabio Costa

ITAMP, Harvard-Smithsonian Center for Astrophysics, Cambridge, MA, 02138, USA

Igor Pikovski

Department of Physics, Harvard University, Cambridge, MA, 02138, USA

Igor Pikovski

Stevens Institute of Technology, Hoboken, NJ, 07030, USA

Igor Pikovski

Vienna Center for Quantum Science and Technology (VCQ), Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Vienna, Austria

Časlav Brukner

Institute for Quantum Optics and Quantum Information (IQOQI), Austrian Academy of Sciences, Boltzmanngasse 3, 1090, Vienna, Austria

Časlav Brukner

Contributions

M.Z., F.C., I.P and Č.B. contributed to all aspects of the research, with the leading input from M.Z.

Corresponding author

Correspondence to Magdalena Zych.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information: Nature Communications thanks Rodolfo Gambini and other anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Peer Review File

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

**Source :**