Connecting processes with indefinite causal order and multi-time quantum states

Recently, the possible existence of quantum processes with indefinite causal order has been extensively discussed, in particular using the formalism of process matrices. Here we give a new perspective on this question, by establishing a direct connection to the theory of multi-time quantum states. Specifically, we show that process matrices are equivalent to a particular class of pre- and post- selected quantum states. This offers a new conceptual point of view to the nature of process matrices. Our results also provide an explicit recipe to experimentally implement any process matrix in a probabilistic way, and allow us to generalize some of the previously known properties of process matrices. Furthermore we raise the issue of the difference between the notions of indefinite temporal order and indefinite causal order, and show that one can have indefinite causal order even with definite temporal order.


I. INTRODUCTION
When describing physical phenomena it is commonly assumed that there exists an underlying causal order. Loosely speaking, later events can be influenced by previous ones, but not the other way around. Recently, however, the idea that physical theories necessarily require a causal order has been challenged, in particular in the context of quantum theory. Indeed, one may imagine that the notion of causal order can be subject to fundamental quantum principles, such as the superposition or uncertainty principle, resulting in indefinite causal structures. Ref. [1] introduced the notions of a "superposition of quantum evolutions" and of what is now known as a "quantum switch" [2]. The question of the existence of quantum correlations with indefinite causal order has recently triggered an intense research effort and subsequently several frameworks for characterizing quantum processes with undefined or dynamical causal structures have been developed, see e.g. [2][3][4][5][6].
In [7], Oreshkov, Costa, and Brukner introduced the framework of 'process matrices' where causal structure can be partially relaxed. Here, operations in local laboratories are described by quantum theory (making use of a local causal structure), however, no global causal structure is assumed. Crucially, this framework captures situations that cannot be explained by any classical causal structure, as witnessed by the violation of so called "causal inequalities" (analogous to Bell inequalities) [7][8][9][10][11][12]. Possible implications for quantum information processing have also been discussed [13,14].
Although the process matrix formalism contains features reminiscent of quantum mechanics, it is not derived from quantum theory. The crucial question of its relation to quantum theory was left open in the original work. It could have been the case that process matrices simply represented a reformulation of some specific situations in quantum mechanics, or that they contained new physics outside of quantum theory.
Here we offer an answer to this question. We show that process matrices correspond to a subclass of "twotime states", which were introduced by Aharonov and his collaborators [15][16][17]. This means that the world described by process matrices is equivalent to a particular case of a quantum world with fundamental post-selection, i.e. a quantum world with independent initial and final boundary conditions, which are both guaranteed to occur (as opposed to ordinary post-selection in which the final state occurs only probabilistically).
In [10] an alternative connection between quantum mechanics and process matrices was found, via a powerful generalisation of the process matrix formalism. In particular, as a corollary, it was shown that any process matrix can be simulated in quantum mechanics using postselection. As any two-time state can be probabilistically implemented experimentally in standard quantum theory via post-selection, we also recover this result in our formalism.
A number of further insights also follow from our results. First, and most importantly, they show that there is a distinction between indefinite temporal order and indefinite causal order, and that one can have indefinite causal order with definite temporal order. Furthermore, they allow us to generalise some of the previously known properties of process matrices, which might be of interest in future work.
Finally, we show that the subset of pre-and postselected states corresponding to process matrices has very special properties, leading to probabilities that are linear functions of the states and of the measurements. This raises a new and important question: Why only this subset? What would be problematic if the probabilities were non-linear functions, as is the case for most pre-and postselected states?

II. PRE-AND POST-SELECTION: EXPERIMENTAL VERSUS FUNDAMENTAL
The standard procedure for collecting statistics in quantum mechanics is to fix an initial state for a system and accumulate experimental results based on this chosen state, or perhaps on an evolution of it later in time. As pointed out by Aharonov et al. in [15], in addition to fixing an initial state for a system, one may also specify an independent final state. One way to realise such a specification is with post-selection, which can be understood by considering the following scenario.
At the initial time t 1 , Alice prepares a quantum system in the state |ψ . In the time interval between the initial and a final time t 2 she performs some arbitrary experiments and records their results. At the final time t 2 , she measures an observable O, one of whose non-degenerate eigenstates is |φ , the desired final state (which is arbitrary and independent of |ψ ). Alice considers her experiment to be successful if the measurement of O yields the eigenvalue corresponding to the eigenstate |φ ; otherwise she discards the experiment.
This way, if Alice repeats her experiment on an ensemble of particles, all prepared in the initial state |ψ , she ends up with a sub-ensemble, which we call a preand post-selected ensemble, characterised by the initial state |ψ and the final state |φ . In this sub-ensemble the statistics of the results of the intermediate measurements is, in general, different from the statistics over the entire ensemble. Importantly, the procedure described above is purely quantum mechanical, albeit not the one that is usually considered (the standard paradigm considers only preselected ensembles). In this procedure, one cannot guarantee a priori that the final measurement will yield the eigenvalue corresponding to |φ . Thus, at no intermediate point can one know if the post-selection will be successful or not (Alice may even decide not to measure O). It follows that there is also no intermediate time at which Alice can know that she is in the desired pre-and post-selected sub-ensemble. Only after the final time, t 2 , will she know which events are successful and must be kept, and which have to be discarded. It is only then, looking back at the records of her intermediate measurement results, that she is be able to find out the statistics corresponding to a specific pre-and post-selected ensemble.
Contrary to this experimental realisation of a pre-and post-selected ensemble, the quantum mechanical formalism also allows one to fix a guaranteed final state. Unlike the case of experimental post-selection, here Alice would already see the statistics corresponding to the pre-and post-selected ensemble at the intermediate times t 1 < t < t 2 . We call this a situation with 'fundamental' post-selection, as opposed to the ordinary, measurement based, probabilistic post-selection.
There is no evidence that fundamental post-selection exists in nature, however, some authors have suggested that it could be present in some exotic situations, for example having a final state of the universe [18][19][20], or a final state at the singularity of a black-hole [21].

III. THE TWO-TIME STATE FORMALISM
Regardless of whether a two-time state arises from experimental or fundamental post-selection, the statistics generated by a given pre-and post-selected ensemble are the same. First introduced for pure states [15,16], the formalism was recently extended to general mixtures of two-time states [17]. In this section we will review this formalism following the convenient notations of Ref. [17].
Suppose, once more, that Alice starts at t 1 with her system in the state |ψ , measures the operator O and selects the cases in which this final measurement yields the eigenvalue corresponding to |φ at time t 2 . For simplicity we assume that the Hamiltonian between t 1 and t 2 is zero and that in the intermediate time Alice performs a detailed measurement, described by the set of Kraus operators {Ê a = k,l β a,kl |k l|}. By "detailed measurement" we mean a measurement where each outcome a corresponds to a single Kraus operator, such that the normalisation condition that the operators obey is aÊ † aÊa = I. The probability to obtain the outcome a given the pre-and post-selection is then given by (1) In order to use the two-time formalism to its full advantage, we transfer to the two-time language. Formally, we take the state space of Alice to be the tensor Hilbert space H A2 ⊗ H A1 , where H A1 is the Hilbert space of the preselected states, which are denoted by ket vectors (with raised labels) and that evolve forward in time, and H A2 is the Hilbert space of the post-selected states, which are denoted by bra vectors (with lowered labels) and that evolve backward in time. The structure is made explicit by defining a two-time state and a two-time version of the Kraus operator in the following way where we differentiate between the usual Kraus operator and its two-time version by the presence or absence of a hat. This notation places states and measurement operators on an equal footing such that they are dual to one another, which allows us to re-write the amplitude φ|Ê a |ψ appearing in (1) as Ψ • E a . The operation (•) applies to vectors belonging to different Hilbert spaces and is a combined composition/contraction: it contracts FIG. 1. Representation of how to experimentally prepare an entangled two-time state. An experimenter prepares the state ij αij |j A 1 ⊗ |i C 1 of Alice's system and an ancilla (denoted C). After Alice has performed her operation, the experimenter then post-selects the final state of Alice and the ancilla to be the maximally entangled state ( where d indicates the dimension of the space on A2 and also that of the ancilla. Note that this procedure can be intuitively understood as entanglement swapping: the entanglement between the system and ancilla is effectively swapped to the two-time state of the system. One can verify that the statistics of the outcomes of Alice's operation in this case obey (4).
any dual vector pairs (i.e. bra and ket pairs) in Hilbert spaces with the same labels to generate a scalar 1 , e.g.
A2 φ| • |i A2 = φ|i ; and it performs the tensor product between any unpaired vectors. In the present example all the vectors are paired, so • effectively results in a scalar product; the tensoring will be useful in later examples. 2 The probability to obtain the outcome a, as seen in (1), now takes on the alternative form The advantage of viewing pre-and post-selections as two-time states in the tensor product Hilbert space H A2 ⊗ H A1 is that one can take superpositions of simple "direct product" two-time states and obtain general pure twotime states in which the pre-selection is entangled with the post-selection. An arbitrary pure two-time state has the form This general two-time state can be understood as a fundamentally post-selected state, or one that can be implemented via experimental post-selection using entangled ancillas [22], the procedure for which is detailed in Fig. 1. Generalizing further, one may also consider mixtures of pure two-time states [17]. Consider that the experimenter follows a similar procedure as above (see Fig. 1).
However, instead of preparing a pure state for A 1 and the ancilla C, the experimenter prepares any possible mixed state, i.e. the ensemble { ij α r,ij |j A1 ⊗ |i C } r with associated probabilities p r .
In this case, the statistics of Alice's operation obey Analogous to the case of ordinary quantum mechanics where the density operator captures the information of a mixed state at a single time, one can also construct a similar object for multi-time states. The density vector of a pure multi-time state Ψ is given by Ψ⊗Ψ † , where the Hilbert spaces pertaining to Ψ † are differentiated by dagger labels, such that every ket (bra) in Ψ is transformed into a daggered bra (ket) in Ψ † . 3 For the pure state in (5), We will use the shorthand The density vector of the ensemble is the convex combination of the pure density vectors, Similar to the construction of a density vector for states, one can also construct a 'Kraus density vector' J a = E a ⊗ E † a for any Kraus operatorÊ a . If the measurement is not detailed with respect to the outcome a (in other words multiple Kraus operators correspond to a single outcome, a so-called 'coarse-grained' measurement), then the density vector corresponding to the outcome is the sum of all the Kraus density vectors corresponding to that outcome. IfÊ µ a denotes the Kraus operator, where a is the outcome and µ is an index running over all the operators corresponding to that outcome, It is important to note that J a contains all the information about the dynamics induced by the measurement, and not only information about the outcome probabilities. That is, it is not equivalent to the POVM element µÊ µ † aÊ µ a , which contains no information about the post-measurement state.
The normalisation condition in the standard formalism, a,µÊ µ † aÊ µ a = I, associated to the fact that the Kraus operators {Ê µ a } a,µ form a completely-positive trace-preserving channel, expressed in our notation is 4 where J = a J a and we have introduced the 'identity' vector I A := I A † A ∈ H A given by Considering the form of (10), we can think of J that satisfy this condition as future identity preserving.
The generalization to mixed states and non-detailed measurements preserves the duality between states and measurements, that is now between the density vectors corresponding to two-time states and the density vectors corresponding to a measurement outcome. The statistics of Alice's operation given by (6) can now be concisely written as The statistics are thus fully captured by the density vector of the two-time ensemble (η), and the density vectors corresponding to each outcome of the measurement (J a ). Note that the denominator, which ensures that the probabilities sum to one, generally depends upon the choice of measurement. In such cases it is impossible to remove it by merely normalising η. 5 Note that because of this, the probabilities are in general non-linear with respect to the state η, as well as the measurement J a . Finally, we note that we defined density vectors η as being mixtures of pure density vectors. We can easily show that this condition is equivalent to asking that the state η is "positive", in the sense that for any vector v⊗v † we have η • (v ⊗ v † ) ≥ 0 (which is the two-time version of the usual definition of a positive operator: A is positive if, for any |ψ , ψ|A|ψ ≥ 0). Hence every density vector η is positive and every positive η is a density vector. Furthermore, note that the J, as defined in (9) are positive in the same sense. In particular, this ensures that η produces positive probabilities via (12). 4 This representation is the analogue of the condition that if Λ(·) = a,µÊ µ a (·)Ê µ † a is a trace-preserving quantum channel, then the conjugate channel Λ † (·) = a,µÊ µ † a (·)Ê µ a is unital, i.e. satisifes Λ † (I) = I. 5 This also means that two density vectors which are identical up to normalization represent the same physical state.
We will be particularly interested in this paper in bipartite two-time states (shared between Alice and Bob), which can be characterised by Fig. 3 for the preparation of such a state via post-selection). Given that Alice measures J a and Bob measures K b , the joint probability to obtain outcomes a and b is given by

IV. THE PROCESS MATRIX FORMALISM
One of the key motivating factors behind the introduction of process matrices in [7] was the ability of the formalism to capture "indefinite causal orders". We will describe the formalism for two parties; the extension to more parties is straightforward and can be found in [11].
The set-up is summarised as follows: Alice and Bob reside in spatially separated laboratories, which are sealed off from the outside world. The doors of Alice's (Bob's) lab may open once to let a system in, and once to let a system out. Within their respective labs Alice and Bob may perform local, quantum mechanical measurements (with the aid of local ancillas, should they require them). Finally, Alice and Bob's labs may be connected in some way. For example it could be the case that on receiving a system, processing it, and then releasing it, Alice passes a quantum system to Bob, who opens his lab doors to let it in (or vice versa). However, they may be connected in completely different ways, corresponding to indefinite causal order, which may not be realisable in quantum mechanics.
The novel element in this formalism is the so-called "process matrix" W which specifies how the two labs are connected, and loosely speaking, plays the role of a joint state for Alice and Bob. More precisely, we associate two Hilbert spaces to Alice, one to the system that enters her laboratory, and one to the system that leaves it, and similarly for Bob. The process matrix is a mathematical object playing the role of a generalised "density operator" acting on these four Hilbert spaces.
On the other hand, measurements in this formalism are usual quantum measurements. Suppose Alice performs a measurement whereÊ µ a denotes the (possibly non-detailed) Kraus operators corresponding to the outcome a; the index µ runs over the various operators corresponding to this individual outcome. In the standard quantum formalism, if the state of Alice's system is ρ, the unnormalised state after the measurement, given that the outcome a was obtained, is ρ a = M a (ρ) = µÊ µ a ρÊ µ † a , with M a : L(H A1 ) → L(H A2 ) denoting the linear map from linear operators on the input Hilbert space to linear operators on the output space, thus describing the measurement. We similarly denote Bob's Kraus operators bŷ In the process-matrix formalism the measurements of Alice are represented (via a Choi-Jamio lkowski transformation) by [7,9] where |Φ + = i |i |i ∈ H A1 ⊗H A 1 is the unnormalized maximally entangled state on two copies of Alice's input Hilbert space, and . Having defined the measurements in this framework the next step is to define the object on which they act. This object is the process matrix W A1A2B1B2 ∈ L(H A1 ⊗ H A2 ⊗ H B1 ⊗ H B2 ), and is an operator acting on the input and output Hilbert spaces of Alice and Bob. The probability rule for the measurement outputs is then given in analogy to the Born rule of quantum mechanics, For W to be a valid process matrix, the probabilities P W (a, b) must be positive and normalised, ab P W (a, b) = 1, in addition to W itself being a positive operator. These are the only requirements to be a valid process matrix, and the set of valid process matrices is any operator satisfying these constraints.

V. CONNECTING PROCESS MATRICES WITH MULTI-TIME STATES
We are now in position to present our main result, namely to connect process matrices and multi-time states. In particular, we will show that to every process matrix we can associate a two-time state that will produce the same probabilities for all measurements. This shows that every process matrix can be realised within quantum theory, if both pre-and post-selection are allowed.
We will present the mapping in two ways; first in terms of matrix elements and second via an isomorphism between Hilbert spaces. Starting with the former, let us denote an arbitrary process matrix as where this decomposition is in the same basis as the |Φ + used to define M A1A2 a and M B1B2 b . Then, to this process matrix we associate the bipartite two-time state Henceforth, we denote a two-time state which has been mapped from a process matrix W as η W . One can explicitly check that W and η W lead to the same probabilities for all measurements. We show in appendix A1 that where J a = µ E µ a ⊗ E µ † a are the Kraus density vectors of Alice's measurement, and K b = ν F ν b ⊗F ν † b of Bob's. Since all process matrices satisfy ab P W (a, b) = 1, the above shows that every η W satisfies and hence which completes the claim, and shows that W and η W lead to the same probabilities. From this reasoning we also deduce that every η W has the special property that the probability rule becomes linear.
In particular, whereas the standard probability rule for two-time states in (12) is a non-linear function with respect to the density vector η and Kraus density vectors J a and K b , for those η W which arise through (17) we have P η W (a, b) = η W •(J a ⊗K b ), which is now a linear function of the two-time state, and measurements. We will refer to any state that satisfies the property as a linear two-time state. Note that this is equivalent to This also highlights the fact that process matrices only map into a subset of two-time states, see Fig. 2.
We have shown above that any process matrix maps into a linear two-time state. Conversely, every linear twotime state corresponds to some process matrix W . In particular, if one writes η as in (17) and then considers the corresponding matrix W given by (16) then the positivity of W follows from the positivity of η, and the normalisation of probabilities for any measurement follows from inserting the linearity condition (21) in (18), and noting that P η (a, b) is normalised. Thus the matrix W is a process matrix. We end by giving a second way in which one can express the connection between process matrices and twotime states via a mapping between operators that act on the Hilbert spaces in the process matrix formalism, and vectors in the Hilbert spaces of the two-time state formalism. Recall that W A1A2B1B2 ∈ L(H A1 ⊗ H A2 ⊗ H B1 ⊗ H B2 ) is an operator acting on a four-party Hilbert space. We introduce the following mappings Applying this set of transformations to any process matrix W as in (16) leads to the two-time state η W as in (17). Note that the second and fourth relations, which map the spaces L(H A2 ) and L(H B2 ) from the process matrix formalism to the bra-vector spaces and their conjugates, H A2 ≡ H A2 ⊗ H A † 2 and H B2 ≡ H B2 ⊗ H B † 2 respectively, are basis-dependent mappings, and must be done in the same basis in which the process matrix is represented (i.e. in this case in the same basis as used in (16)).

VI. THE SET OF TWO-TIME STATES CORRESPONDING TO PROCESS MATRICES
In [7] it is shown that the condition that process matrices give normalised probabilities (that are also positive) for all possible measurements is in fact equivalent to a finite set of conditions, which provide an explicit and compact characterisation of the set of valid process matrices. Here we give a translation of these conditions in the two-time state formalism.
In the bipartite case the set of process matrices is specified by the following five necessary and sufficient condi-tions: where d A2 (d B2 ) is the dimension of Alice's (Bob's) output Hilbert space; the notation in the last three equations is defined using the 'trace-and-replace' operation Returning to two-time states, we first need to introduce two operations defined in terms of the identity vector (11). The first corresponds to a 'throw-away-and-replace' operation, replacing a pre-selected system A by the maximally mixed (pre-selected) state on A , and is analogous to the 'traceand-replace' operation from (24) for process matrices. The second corresponds to the 'do-nothing' operation, taking any state of A to the same state of A . Note that in both of these operations we introduce new primed spaces. Depending on the context, these primed spaces should either be thought of as corresponding to a time t 1 just after time t 1 , or to a time t 2 just before time t 2 . For example, T is a two-time state between the times t 1 and t 2 .
With these operations, a translation of the conditions in (23) is 6 η W ≥ 0, (28) 6 Strictly speaking, the third condition from (23) translated into two-time notation reads which is the same as the third condition in (28), with maximally mixed states on Bob's input and output spaces tensored in on both sides of the equation, which does not change the condition. A similar argument holds for the fourth conditions from (23) and (28). 3. Experimental preparation of a mixed bipartite twotime state η. At time t1, the experimenter prepares the mixed state ρ = ijkl,pqrs η ijkl,pqrs |i |j |kl p| q| rs| ∈ H A 1 ⊗ H B 1 ⊗ H C 1 between Alice, Bob and an ancilla (C1) of dimension dA 2 dB 2 . After Alice and Bob have performed their operations, at a later time t2, the experimenter post-selects on the pure state mn m| ⊗ n| ⊗ mn| ∈ HA 2 ⊗ HB 2 ⊗ HC 2 . The effect of this is to have created the bipartite two-time state η between times t1 and t2 with the equivalent statistics given by the operators Ja = Ea ⊗ E † a and

FIG.
The double arrow on the ancilla indicates that the dimension is greater than that of Alice or Bob's systems.

VII. GENERALISED PROPERTIES OF PROCESS MATRICES
In this section we show that the properties (28) of process matrices, that are expressed in terms of the 'throwaway-and-replace' and the 'do-nothing' channels (T A A and I A A ), can be generalized to arbitrary channels. This demonstrates the usefulness of connecting the process matrix formalism to the formalism of two-time states, and might be useful in future research on process matrices, as it provides further structural information about them.
In Appendix B we prove the following result concerning linear bipartite two-time states (i.e. about η W states): Theorem Given any η W , i.e. any bipartite two time state that satisfies 7 for all J A 2 A 1 and K B 2 B 1 corresponding to completely positive trace preserving maps (i.e. of the form (9) and such that ) then the following 7 Since the overall normalisation of a two-time state is unphysical, given a two-time state η that satisfies (J B 1 , we can always re-scale η → 1 λ η such that this condition holds. properties hold: corresponding to completely positive trace preserving maps. 8 We can translate the above conditions into the process matrix formalism, and present a new set of conditions that are satisfied by all process matrices. To do so, let us introduce arbitrary Choi-Jamio lkowski operators corresponding to completely positive and trace-preserving (CPTP) maps R A 2 A2 ,R A 2 A2 and M A1A2 for Alice, and similarly S B 2 B2 ,S B 2 B2 and N B1B2 for Bob. Then, the conditions (30) can be re-expressed as (31) 8 Note that the final condition of (30) can be expressed more sym- )] • ηW = 0, but we adopt the above form to emphasise its similarity with (23).
These conditions are not directly implied by the conditions (23), and may prove of independent interest for proving new results regarding process matrices.

VIII. INDEFINITE CAUSAL STRUCTURE VERSUS INDEFINITE TEMPORAL STRUCTURE.
In the usual causal scenario an event A can be a cause of an event B only if A takes place before B, and B can be a cause of A only if B takes place first. The most general causal structure is a simple mixture of these cases, by probabilistically sometimes letting A be before B and sometimes B before A. In quantum mechanics one can envisage a more interesting situation in which one can arrange for a superposition of these two time orderings, resulting in an "indefinite" causal order. The process matrix formalism is inspired by this idea.
Note that in the above the causal structure is constrained by the temporal structure, and indefinite causal order is induced by indefinite temporal order. On the other hand, given the correspondence between two-time states and process matrices presented here, we have implicitly shown that in quantum mechanics with postselection (and presumably in every post-selected probabilistic theory) temporal order and causal order are independent notions. In particular one can have indefinite causal order even with definite temporal order.
In the process matrix formalism the actual time at which Alice's and Bob's actions take place is left undefined. The only ordering we have is a local time order inside Alice and Bobs labs, and nothing is said about the relative timing between Alice and Bob. The possibility is left open for Alice to pass her output state to Bob as his input state, and for Bob to pass his output state to Alice as her input state, in some sort of potential superposition, though how exactly this is encoded in the formalism is less clear.
On the other hand, in the two-state formalism the temporal order is perfectly well-defined: pre-selection takes place at time t 1 , post-selection at time t 2 , and the actions of both Alice and Bob take place between these times. Since their initial times are synchronised, in the usual causal scenario, there would be no way for Alice to send her final state to Bob, or vice-versa. However, using post-selection "effective communication" can be realised using the entanglement present in the joint twotime state, which enables, for example, Alice's final state to be teleported backward in time to Bobs initial time. Thus we can realise a causal order from Alice to Bob, or vice versa, or even an indefinite causal order (when the joint two-time state corresponds to a process matrix with indefinite causal order), despite having definite temporal order. 9 9 Note that in fact any fixed timings for the pre-and post-selections

IX. DISCUSSION
In our paper we addressed the question of the relation between the process matrix formalism and quantum theory: does the process matrix formalism imply new phenomena, not present in quantum theory or not? We showed that correlations predicted by the process matrix formalism are identical to those obtainable in a subset of pre-and post-selected quantum states. In this sense, the process matrix formalism is contained within quantum theory. The subtle question is that of the probability with which we can prepare situations described by process matrices. In ordinary quantum theory pre-and post-selected states cannot be prepared with certainty, unless they are trivial (i.e. no post-selection): as their name suggests, sometimes we don't succeed in obtaining them, so, at the end of the experiment we need to select the cases when we were successful and reject the other cases. Hence, if Nature would turn out to be such that all process matrices can be prepared with certainty, then Nature may not be described by quantum theory, but something new. (Here we say "may not" rather than "would not", since we have not proven that trivial post-selections are not enough, or that quantum mechanics does not allow simulating W matrices by some other means than by pre-and post-selection.) On the other hand, if one enlarges the scope of quantum theory to quantum theory + fundamental post-selection, in which nature provides the "post-selected" state with certainty (such as by giving a final state of the universe, in addition to and independent of the initial state), then the situations described by the process matrix formalism are completely contained within it.
As noted above, the correlations described by the process matrix formalism are the same as those arising from a subset of pre-and post-selected states. This raises a new and important question: Why only this subset? Why not other pre-and post-selected states? As we have seen, the limitation to this set of states stems from the requirement that the probabilities are linear functions of the states and of the measurements. What would be problematic if the probabilities were non-linear functions? After all, in the context of pre-and post selection the probabilities are non-linear in general. We will address this question in forthcoming work.
We generalised the defining conditions of the process matrices which are given by relations that are obeyed when the systems are subjected to simple "trace and replace channels to arbitrary channels. This may be useful in proving further results, or to allow us to gain a deeper understanding of process matrices.
We also discussed the issues of indefinite temporal and of Alice and Bob, not necessarily the synchronised timings considered above, allows for indefinite causal order with definite temporal order.
causal order and argued that they are actually two independent concepts. Finally, we note a technical issue which stems from a deeper conceptual issue and shows an advantage of viewing process matrix situations within the two-time state formalism. The process matrix formalism, which is based on the Choi-Jamio lkowski representation, and in which process-matrices are like density operators on standard ket vector spaces, is basis dependent. In particular, this leads to the rather complicated relation (14) between the Kraus operators that describe Alice and Bobs measurements and the corresponding operators M A1A2 a and N B1B2 b that are used for calculating probabilities. On the other hand, the two-time formalism, in which a ket vector and a bra vector space are used to represent the forward and backward in time propagating states respectively (and vice versa for the daggered spaces), is basis independent. Indeed, the state i |i A1 |i A2 is basis dependent, while the two-time equivalent i |i A1 ⊗ A2 i| is basis independent. Correspondingly, the relation (9) between the Kraus operators E a and the operator J a used for calculating probabilities is straightforward.

APPENDIX Appendix A: Proof of the connection between process matrices and two-time states
Here we give the proof of (18), which connects P W (a, b) and η W .
Appendix B: Proof of the generalized properties of process matrices In this appendix we give the proof of the generalized properties of process matrices given in the main text (30). We begin by restating the theorem: Theorem Given any η W , i.e. any bipartite two time state that satisfies for all J A 2 A 1 and K B 2 B 1 corresponding to completely positive trace preserving maps (i.e. of the form (9) and such that I A 2 • J A 2 A 1 = I A 1 and I B 2 • K B 2 B 1 = I B 1 ) then the following properties hold: corresponding to completely positive trace preserving maps. Proof: Consider the following CPTP map for Alicẽ where T A 2 A 1 is the throw-away-and-replace operation (c.f. (25)), C A 2 correspond to two arbitrary CPTP maps, A 1 ∈ H A1 ⊗ H A 2 is an arbitrary vector 10 , and > 0 is a positive constant, taken to be sufficiently small such that J A 2 A 1 is positive. Similarly, for Bob,K (B4) 10 The only requirement on X are two arbitrary CPTP maps, Y B 2 B 1 ∈ H B1 ⊗H B 2 is an arbitrary vector, and δ > 0 is a sufficiently small positive constant so thatK B 2 B 1 is positive. Now, for linear two-time states it follows that for all CPTP maps K B 2 B 1 . Using (B3) to expand outJ A 2 A 1 , this is seen to be equivalent to where we have used the fact that > 0 to cancel it. The only way that this can hold for an arbitrary vector X which is the first condition of the theorem. An identical argument can be made, by starting from (J A 2 is an arbitrary CPTP map, which leads to the second result of the theorem, Finally, by using the relation it follows that having used the two results already proved. Since this must hold for all X which completes the proof.