## Abstract

The standard formalism of quantum theory treats space and time in fundamentally different ways. In particular, a composite system at a given time is represented by a joint state, but the formalism does not prescribe a joint state for a composite of systems at different times. If there were a way of defining such a joint state, this would potentially permit a more even-handed treatment of space and time, and would strengthen the existing analogy between quantum states and classical probability distributions. Under the assumption that the joint state over time is an operator on the tensor product of single-time Hilbert spaces, we analyse various proposals for such a joint state, including one due to Leifer and Spekkens, one due to Fitzsimons, Jones and Vedral, and another based on discrete Wigner functions. Finding various problems with each, we identify five criteria for a quantum joint state over time to satisfy if it is to play a role similar to the standard joint state for a composite system: that it is a Hermitian operator on the tensor product of the single-time Hilbert spaces; that it represents probabilistic mixing appropriately; that it has the appropriate classical limit; that it has the appropriate single-time marginals; that composing over multiple time steps is associative. We show that no construction satisfies all these requirements. If Hermiticity is dropped, then there is an essentially unique construction that satisfies the remaining four criteria.

### 1. Introduction

Quantum theory, as usually formalized, contains a fundamental asymmetry between space and time. This is evident when considering the description of composite quantum systems. If a composite system consists of several components existing at a given time, then the formalism specifies that the joint state is given by a density matrix acting on the tensor product of the Hilbert spaces associated with the components. (This is the case regardless of whether the components are spatially separated systems or different degrees of freedom of a single system.) But in principle, there is another sort of composite system besides these, namely one where the components are time-like separated. In this case, there is no standard prescription for the joint state of the composite. Rather, states are defined only at a single time, and evolve over time under the action of a Hamiltonian. But is this asymmetry fundamental or merely an artefact of a particular formalism? Relativity theory has revealed that many distinctions between space and time previously thought to be fundamental are not. While spatial and temporal dimensions do have opposite signs in the metric and causal ordering is preserved under Lorentz transformations, the distinction between spatial and temporal intervals becomes a merely observer-dependent decomposition of the fundamental concept of a spatio-temporal interval. This naturally prompts the question of whether the distinction in quantum theory between how one models composite systems at a single time and how one models composite systems over time is like the distinction between spatial and temporal intervals, and can be eliminated, or whether it is like the sign distinction in the metric and cannot. The asymmetry is also in marked contrast with classical probability theory, where joint probabilities can be defined for sets of events whatever their spatio-temporal relationships. Hence, although there is a well-developed analogy between density matrices and classical probability distributions, the analogy is limited in scope when it comes to time-like separated systems.

If the apparent asymmetry between space and time in quantum theory is fundamental, then the understanding of time given by quantum theory must be different from that suggested by a combination of relativity and classical probability theory. If, on the other hand, it is not fundamental, then it should be removable. One way of removing the asymmetry would be to construct quantum states for composites over time. States would then be defined across both space and time, without a separate formalism describing evolution over time.

There have been a variety of proposals for expressing quantum theory in a manner that treats space and time in a more even-handed fashion. The sum-over-histories approach to quantum dynamics [1,2] and the programme of consistent/decoherent histories [3–5] are examples. An alternative family of proposals includes the multi-time formalism (see [6,7] and follow-up work [8,9]), quantum combs (see [10], follow-up work [11,12] and related formalisms [13,14]), process matrices (see [15] and follow-up work [16]) and the causaloid formalism [17,18]. These latter proposals, although differing from one another in notation and in the types of problems that have been addressed, share the feature that a quantum system at a localized region of space–time is associated with two Hilbert spaces, one carrying an incoming state, and one an outgoing state. Part of the reason for this is that each proposal allows for the possibility that an agent situated at that region of space–time can intervene upon the system, with the intervention corresponding to a quantum instrument: the quantum instrument is a set of trace non-increasing completely positive maps, one for each classical outcome of the intervention, that mediate the incoming and outgoing Hilbert spaces.

In such approaches, the main criterion of success is whether the formalism can be used to compute the joint probability distribution over the outcome variables for a given set of interventions. Our project, by contrast, is also motivated by the goal of providing a *causal account* of the operational predictions. To do so, it is critical to have a formulation of quantum theory that makes a clean distinction between the aspects of the formalism that are about causal influence and those that are about Bayesian inference. For instance, rather than merely predicting a correlation between the outcomes of two measurements, such a formulation specifies whether this correlation is due to a common cause of the two variables or rather a cause–effect relationship between them [19]. One proposal for how to achieve this separation takes joint quantum states to be inferential objects, the quantum analogues of joint probability distributions [20]. Given that a joint probability distribution is the appropriate way to describe an agent's incomplete information of a composite system, *regardless* of the spatio-temporal relations that hold among its components, it is natural to ask whether there is a notion of a joint quantum state that is similarly applicable to an arbitrary composite and which coincides with the standard notion when the composite consists of a set of systems considered at a given time. We focus on a special case of this question: whether there is a notion of a joint quantum state for a composite over time which mirrors that of a joint quantum state for a composite at a single time.

If such a project could succeed, not only would it extend and strengthen the analogy between quantum states and classical probability distributions, it would also lend support to the view that a quantum state can be thought of as information directly about a system as opposed to a form of information that can only be expressed in terms of outcomes of potential measurements on the system. That is, it would lend support to the view that the quantum state represents a state of knowledge of an underlying reality associated to the system as opposed to merely representing a state of knowledge of the outcomes of measurements that one might implement on the system. It would also shed light on the subject of the quantum-classical correspondence: if a subset of the systems in a network describe measurement pointers (and are therefore suitably macroscopic and decohered), then the joint quantum state for this subset would encode the joint probability distribution over the pointer variables. Such a formulation could be of practical use for the analysis of quantum information-processing protocols. It could also help to pinpoint where any remaining asymmetry between space and time arises in the quantum formalism.

We therefore undertake to investigate the possibility of defining quantum joint states over time in a manner that closely matches the standard quantum treatment of composite systems at a given time. In particular, we assume that a temporally localized *d*-dimensional quantum system is associated with a single *d*-dimensional Hilbert space (rather than a pair of such Hilbert spaces, as is assumed in [6–18,21]), and that the state of a composite system over time is an operator on the tensor product of the Hilbert spaces associated with the temporally localized components. We explore a number of definitions along these lines, including a proposal of Leifer & Spekkens (LS) [20] (inspired by the view that quantum theory is a generalization of classical probability theory), a proposal of Fitzsimons, Jones and Vedral (FJV) [22], and a novel proposal based on discrete Wigner functions. Following a close examination of these different options, we distil five criteria that the state of a composite over time should satisfy if it is to have the essential qualities of the state of a composite at a single time. Our main theorem states that there is no way to construct a state over time that satisfies all of these criteria. We also show that if one of the criteria is dropped—namely, the criterion that the state over time is a Hermitian operator—then the only construction satisfying the remaining four criteria corresponds to straightforward matrix multiplication of the initial state and the operators representing the quantum channels that evolve from one time slice to the next. Finally, we show that if Hermiticity is retained, but an assumption of associativity is dropped, then the FJV construction satisfies the remaining criteria.

### 2. States over time

We begin with the simplest situation in which to discuss the construction of a quantum state over time: a single system considered at two different times, as depicted in figure 1. Although Hilbert spaces are conventionally associated with systems that persist over time, we shall here associate Hilbert spaces to quantum systems that are considered to be localized to a given region of space–time (similar in spirit to ‘events’ in relativity). The system at the first time is denoted by *A*, and that at the second time by *B*. Assuming that the system is initially uncorrelated with its environment, the most general possible evolution of the system corresponds to a quantum channel, that is, a trace-preserving completely positive map, denoted by ${\mathcal{E}}_{B|A}$. Conventionally, the quantum state of *A* is represented by an operator acting on ${\mathcal{H}}_{A}$, denoted by *ρ*_{A}, and the quantum state of *B* is represented by an operator acting on ${\mathcal{H}}_{B}$, denoted *ρ*_{B}. The state of *B* can be determined by the state of *A* and the channel, via ${\rho}_{B}={\mathcal{E}}_{B|A}({\rho}_{A})$.

It is natural to consider the possibility of a quantum state for the composite system *AB*, represented by an operator *ρ*_{AB} acting on ${\mathcal{H}}_{A}\otimes {\mathcal{H}}_{B}$. The scenario of figure 1 is fully defined by the state *ρ*_{A} and the channel ${\mathcal{E}}_{B|A}$ (as these determine the output state *ρ*_{B}), so we assume that the state over time is a function of these:

*ρ*

_{AB}can be represented by a square matrix, but we do not require that

*ρ*

_{AB}be a standard density matrix. The properties that should be required of the function

*f*, and the possible forms

*f*can take, are the main questions addressed in this paper.

#### (a) Three possible definitions of a state over time

It will be useful in the following to represent the channel ${\mathcal{E}}_{B|A}$ by an operator acting on ${\mathcal{H}}_{A}\otimes {\mathcal{H}}_{B}$, denoted by *E*_{B|A}. As it is also useful if the representation of the channel is basis-independent, we choose it to be the operator that is Jamiołkowski-isomorphic [23] to ${\mathcal{E}}_{B|A}$, defined as

*i*〉} is any orthonormal basis of ${\mathcal{H}}_{A}$. (Up to normalization, this operator is the partial transpose of the operator that is Choi-isomorphic [24] to ${\mathcal{E}}_{B|A}$.) The operator

*E*

_{B|A}is independent of the basis used in its definition, and the state of the output of the channel can be expressed in terms of it as [20]

*E*

_{B|A}the

*channel state*.

The joint state *ρ*_{AB} can be understood as being given by a binary operation ⋆ that takes a pair of operators on ${\mathcal{H}}_{A}\otimes {\mathcal{H}}_{B}$ to an operator on ${\mathcal{H}}_{A}\otimes {\mathcal{H}}_{B}$:

We call this operation a *star product*. For notational convenience, we will in the following often suppress tensor products with identity, writing for example *ρ*_{A} instead of ${\rho}_{A}\otimes {\mathrm{\U0001d7d9}}_{B}$, so that

There have been a number of previous attempts to define a joint state *ρ*_{AB} along these lines. The proposal due to LS constructs joint states over time in the context of a programme of reframing quantum theory as a theory of Bayesian inference [20], which builds on [25,26]. Just as a density matrix on a composite of *A* and *B*, *ρ*_{AB}, can be understood to be a quantum analogue of a joint probability distribution, *P*(*XY*), LS observe that a channel state *E*_{B|A} can be understood as a quantum analogue of a conditional probability distribution, *P*(*Y* |*X*), which describes a classical channel with input *X* and output *Y* . For example, the analogue of the normalization condition for conditional probabilities, $\sum _{Y}P(Y|X)=1$, is ${\mathrm{Tr}}_{B}({E}_{B|A})={\mathrm{\U0001d7d9}}_{A}$, and the analogue of the fact that the output of a classical channel satisfies $P(Y)=\sum _{X}P(Y|X)P(X)$ is the fact that the channel state satisfies *ρ*_{B}=Tr_{A}(*E*_{B|A}*ρ*_{A}). In the classical case, there is no obstacle to defining a joint probability distribution for the pair of systems constituting the input and the output of a channel via *P*(*XY*)=*P*(*Y* |*X*)*P*(*X*). In the quantum case, however, the operators *E*_{B|A} and ${\rho}_{A}\otimes {\mathrm{\U0001d7d9}}_{B}$ do not in general commute, so it is not obvious what should replace this equation. LS considered the following option:

*ρ*

^{(LS)}

_{AB}is Hermitian, and although it is not in general positive, it is locally positive, meaning that 〈

*ab*|

*ρ*

^{(LS)}

_{AB}|

*ab*〉≥0 for all |

*a*〉 in ${\mathcal{H}}_{A}$ and |

*b*〉 in ${\mathcal{H}}_{B}$. Significant problems with this proposal were already noted by LS in Sec. VII of [20].

An alternative proposal for a state over time is that of Fitzsimons *et al.* [22]. As it stands, the proposal is for multi-qubit systems, rather than systems of arbitrary dimension. In this work, it will suffice to consider the FJV construction for the special case in which *A* and *B* are qubits. In order to define the FJV joint state, imagine that a Pauli measurement *σ*_{i} is performed on *A*, with outcome ±1 (here, the Pauli operators are defined as ${\sigma}_{0}=\mathrm{\U0001d7d9},\hspace{0.17em}{\sigma}_{1}={\sigma}_{x},\hspace{0.17em}{\sigma}_{2}={\sigma}_{y},\hspace{0.17em}{\sigma}_{3}={\sigma}_{z}$). After this measurement, the state of *A* is imagined to be updated according to the projection postulate, before evolving according to the channel ${\mathcal{E}}_{B|A}$, and a Pauli measurement *σ*_{j} is then performed on *B*. Let 〈*σ*_{i}⊗*σ*_{j}〉 denote the expectation value of the product of the outcomes of the two measurements. The FJV joint state is given by

*A*and on

*B*, the joint state is intended as a description of the system

*AB*in the case where no measurements are made on either, but rather

*A*simply evolves into

*B*according to the given channel. The matrix

*ρ*

^{(FJV)}

_{AB}is Hermitian, but not necessarily positive (nor even necessarily locally positive). It is clear that

*ρ*

^{(FJV)}

_{AB}depends only on

*E*

_{B|A}and

*ρ*

_{A}, hence can in principle be written in the form of equation (2.3), for a suitable definition of the star product. In fact, as we show in appendix A, the FJV construction for two qubits corresponds to the Jordan product:

Finally, one might also hope to construct a state over time in a manner that is analogous to the definition of a joint probability distribution, but where probabilities are replaced by the *quasi-probabilities* appearing in a discrete Wigner representation of the states and channels [27,28]. A discrete Wigner (W) representation for a *d*-dimensional system *A* is defined by a set ${\mathrm{\Omega}}^{A}=\{{K}_{i}^{A}\}$ of *phase-point operators* which form a basis for the space of operators on ${\mathcal{H}}_{A}$, satisfying $\mathrm{Tr}({K}_{i}^{A}{K}_{j}^{A})=d{\delta}_{ij}$ and $\sum _{i}{K}_{i}^{A}=d{\mathrm{\U0001d7d9}}_{A}$, hence $\mathrm{Tr}({K}_{i}^{A})=1$. A density matrix *ρ*_{A} can be written ${\rho}_{A}=\sum _{i}{r}^{A}(i){K}_{i}^{A}$, where *r*^{A}(*i*) is a real-valued function on *Ω*^{A}, with $\sum _{i}{r}^{A}(i)=1$. Hence a system *A* can be described by the function *r*^{A}(*i*)∈[−1,1], which has the form of a quasi-probability distribution on a discrete phase space. Similarly, the operator *E*_{B|A} representing a channel can be written ${E}_{B|A}=\frac{1}{d}\sum _{ij}{r}^{B|A}(j|i){K}_{i}^{A}\otimes {K}_{j}^{B}$, where *r*^{B|A}(*j*|*i*) is a real-valued function on *Ω*^{A}×*Ω*^{B}. The extra factor of 1/*d* is introduced so that *r*^{B|A}(*j*|*i*) satisfies $\sum _{j}{r}^{B|A}(j|i)=1$; these can therefore be thought of as conditional quasi-probabilities. The natural definition for the W representation of the state of the composite *AB* is then *r*^{AB}(*ij*)≡*r*^{B|A}(*j*|*i*)*r*^{A}(*i*). This implies that the state over time is represented by the operator

*ρ*

_{A}and

*E*

_{B|A}, has the form

There are various ways, then, in which one might aim to define quantum states over time. This leads us to pose the question: can we isolate specific axioms or desiderata that a state over time ought to satisfy? Furthermore, given those axioms, is it possible to find a construction that satisfies them all?

### 3. Five criteria for a star product

We introduce five basic axioms for a star product, motivating each by a corresponding desideratum for the properties of a state over time.

#### (a) Hermiticity

We assume that, for a finite-dimensional system considered at *k* distinct time steps, the state over time, *ρ*_{A1…Ak}, is a Hermitian operator on the *k*-fold tensor product of copies of its Hilbert space, ${\mathcal{H}}_{{A}_{1}}\otimes \cdots \otimes {\mathcal{H}}_{{A}_{k}}$. We denote the space of such operators by ${\mathbb{H}}_{{A}_{1}\dots {A}_{k}}$, so that ${\rho}_{{A}_{1}\dots {A}_{k}}\in {\mathbb{H}}_{{A}_{1}\dots {A}_{k}}$. In particular, if *A* and *B* denote a system at two times, then ${\rho}_{AB}\in {\mathbb{H}}_{AB}$.

Given that ${\rho}_{A}\otimes {\mathrm{\U0001d7d9}}_{B}\in {\mathbb{H}}_{AB}$ and ${E}_{B\hspace{0.17em}|\hspace{0.17em}A}\in {\mathbb{H}}_{AB}$, it follows that we are assuming that the star product of equation (2.2) is a map of the form

#### (b) Preservation of probabilistic mixtures

The output state of a quantum channel for a given input state is a linear function of both the input state and the channel. This is necessary to ensure that a probabilistic mixture of input states maps under the channel to the corresponding mixture of output states, and that a probabilistic mixture of channels acting on a given input state yields the corresponding mixture of output states. By analogy, we would like the joint state over time for a probabilistic mixture of input states to be the corresponding mixture of the joint states for each input state, and similarly for probabilistic mixtures of channels.

Consider as an example the situation in figure 1, but where the quantum system *A* is conditioned on a classical variable *x*. For instance, it could be that *x* is the outcome of a fair coin toss, and a qubit *A* is prepared in a state that depends on this outcome: *ρ*_{A,x=h}≡|0〉〈0| for heads, and *ρ*_{A,x=t}≡|1〉〈1| for tails. If the channel between *A* and *B* is the identity channel, then ${E}_{B|A}^{(\mathrm{id})}\equiv 2|{\varphi}^{+}\u27e9\u27e8{\varphi}^{+}{|}^{{T}_{B}}$ where $|{\varphi}^{+}\u27e9\equiv (1/\sqrt{2})(|0\u27e9|0\u27e9+|1\u27e9|1\u27e9)$ and where *T*_{B} denotes the partial transpose on *B*. In this case, we would like the star product to satisfy

In general, we assume the following:

While the FJV and W formulations are both straightforwardly convex-bilinear, there is a failure of this axiom in the case of the LS star product, which does not preserve convex combinations in *ρ*_{A}. For example, using the classical coin toss example, it is easy to show that the LS star product yields the mixture $\frac{1}{2}|00\u27e9\u27e800|+\frac{1}{2}|11\u27e9\u27e811|$ for the right-hand side of (3.1), but the non-separable |*ϕ*^{+}〉〈*ϕ*^{+}|^{TB} for the left-hand side.

#### (c) Preservation of classical limit

In classical probability theory, defining joint probabilities for events at different times is unproblematic. This gives a constraint on the form of quantum temporal joint states, if it is assumed that classical joint probabilities must be reproduced when a quantum system is behaving entirely classically. For the situation of figure 1, this is the case if both the input quantum state and the operator corresponding to the channel are diagonal in the same basis. Consider, for example, a channel that is completely dephasing in some basis, and an input state that is diagonal in the same basis. For some probabilities *p*(*i*) and conditional probabilities *p*(*j*|*i*), the input state can be written as ${\rho}_{A}=\sum _{i}p(i)|i\u27e9\u27e8i{|}_{A}$, and the behaviour of the channel is given by

*p*(

*i*,

*j*)=

*p*(

*i*)

*p*(

*j*|

*i*), which, when encoded in a matrix in the same way, are given by

In order to reproduce classical probability theory, therefore, when states and channel states are simultaneously diagonalizable the joint state over time should be given by their matrix product. This is ensured if

*x*.

It is easy to see that the LS joint state satisfies *ρ*^{LS}_{AB}=*E*_{B|A}*ρ*_{A} when *E*_{B|A} and *ρ*_{A} commute. This is also the case for the FJV joint state for a qubit at two times, as is evident from the form of equation (2.6). However, the W construction does not satisfy this axiom. Consider, for example, an input state that is a pure state of a qubit, and a completely dephasing channel that preserves this state. Direct calculation shows that unless the input state is an eigenstate of one of the Pauli operators, the Wigner joint state ${\rho}_{AB}^{(W)}\ne {E}_{B|A}{\rho}_{A}$.

#### (d) Preservation of marginal states

An obvious desideratum for a joint state *ρ*_{AB} is that it returns the correct marginal states for *A* and *B*, i.e. Tr_{B}*ρ*_{AB}=*ρ*_{A} and Tr_{A}*ρ*_{AB}=*ρ*_{B}. Given that *ρ*_{AB}=*E*_{B|A}⋆*ρ*_{A} by definition, and that *ρ*_{B}=Tr_{A}(*ρ*_{A}*E*_{B|A}), as noted earlier, it follows that we require Tr_{A}(*E*_{B|A}⋆*ρ*_{A})=Tr_{A}(*ρ*_{A}*E*_{B|A}), hence also

#### (e) Compositionality

The final axiom for the star product concerns how it composes when we wish to describe a system at more than two distinct times. Consider the situation of figure 2, with three localized regions, labelled *A*, *B* and *C*, joined by two channels, ${\mathcal{E}}_{B|A}$ and ${\mathcal{E}}_{C|B}$.

There are two ways in which one can imagine forming the joint state *ρ*_{ABC} on the composite of the three regions. The first way proceeds as follows: form the joint state over the first two regions, *ρ*_{AB}; define a novel kind of channel, denoted by ${\mathcal{E}}_{C|AB}$, whose input is the composite of the first two regions and whose output is the third region; finally, take the star product of *ρ*_{AB} with *E*_{C|AB}, the operator corresponding to ${\mathcal{E}}_{C|AB}$:

*ρ*

_{A}with

*E*

_{BC|A}, the operator corresponding to this channel:

We now consider the consequences, for both approaches, of the fact that the only dependence that *C* has on *A* is mediated by *B*.

In the first approach, this implies that the channel ${\mathcal{E}}_{C|AB}$ is simply ${\mathcal{E}}_{C|B}$, so that ${E}_{C|AB}={E}_{C|B}\otimes {\mathrm{\U0001d7d9}}_{A}$. Using this, together with the fact that *ρ*_{AB}=*E*_{B|A}⋆*ρ*_{A}, equation (3.9) becomes

*E*

_{BC|A}≡

*E*

_{C|B}⋆

*E*

_{B|A}. This can be understood as a quantum analogue of the Markov condition for classical causal networks [29].

We have here allowed a channel state to appear as the second argument of the star product. This is, however, no more objectionable than allowing a state over time to appear as the second argument, as we did in equation (3.9). This is because, for every channel state *E*_{B|A}, there is a state over time *ρ*_{AB} such that *E*_{B|A}=*d*_{A}*ρ*_{AB} (where *d*_{A} is the dimension of the Hilbert space of *A*). To see that this is correct, we note that, for channel state *E*_{B|A} and input state $(1/{d}_{A}){\mathrm{\U0001d7d9}}_{A}$, the state over time is ${\rho}_{AB}={E}_{B|A}\star (1/{d}_{A}){\mathrm{\U0001d7d9}}_{A}$. Therefore, by our assumption of being product on commuting pairs, this reduces to *ρ*_{AB}=(1/*d*_{A})*E*_{B|A}. Invoking convex-bilinearity, we can therefore write *E*_{C|B}⋆*E*_{B|A} as *d*_{A}(*E*_{C|B}⋆*ρ*_{AB}).

Using *E*_{BC|A}≡*E*_{C|B}⋆*E*_{B|A}, equation (3.10) then becomes

Our last condition on the star product is motivated by the idea that the two techniques just described for defining a joint state over three times should give the same result,

The W formulation defines an associative star product. The LS star product, however, is not associative, as noted by LS in Sec. VII of [20]. In particular, the second argument of the LS star product is required to be a positive operator—for the square root of the operator to be uniquely defined—and yet non-positive operators arise as joint states over time. Indeed, the fact that the LS star product does not support composition across time was highlighted by LS as one of its most significant shortcomings, motivating the search for alternative approaches.

The FJV construction is interesting with regard to associativity. One can verify that the FJV construction for a single qubit at three times gives a joint state:

_{FJV}is the Jordan product. Although the Jordan product is non-associative in general, it is associative for triples of operators of the form ${E}_{C|B}\otimes {\mathrm{\U0001d7d9}}_{A},{E}_{B|A}\otimes {\mathrm{\U0001d7d9}}_{B}$ and ${\rho}_{A}\otimes {\mathrm{\U0001d7d9}}_{B}\otimes {\mathrm{\U0001d7d9}}_{C}$, that is (suppressing identities as usual)

*is*associative when restricted to operators of the correct form to describe a quantum system passing through a sequence of channels. Hence a unique joint state is defined without privileging any particular way of grouping the different time steps. When one considers an arbitrary number of qubits, the situation becomes more subtle: an FJV star product can again be inferred, but in this case it is not the Jordan product of equation (2.6). A similar situation arises for systems whose dimension is a power of a prime. For other systems, the proposal is ambiguous. We omit the details as these issues are non-essential for our purposes.

### 4. Main theorem

We now have five criteria for states across time, expressed as axioms for a star product: Hermiticity, convex-bilinearity, product when commuting, product when traced, and associativity. The satisfaction or otherwise of these axioms, for each of LS, FJV and the W construction, are summarized in table 1. None of these constructions satisfies all of the axioms.

convex- | product | product | |||
---|---|---|---|---|---|

Hermitian? | bilinear? | on commuting pairs? | when traced? | associative? | |

LS | ✓ | × | ✓ | ✓ | × |

FJV | ✓ | ✓ | ✓ | ✓ | ×^{a} |

W | ✓ | ✓ | × | ✓ | ✓ |

### Theorem 4.1

*Let* ${\mathbb{H}}_{n}$ *be the set of n-by-n Hermitian matrices. There is no function* $\star :{\mathbb{H}}_{n}\times {\mathbb{H}}_{n}\to {\mathbb{H}}_{n}$ *satisfying convex-bilinearity, product on commuting pairs, product when traced, and associativity.*

In other words, there is no definition of a state over time that satisfies the five axioms given in this paper. We provide the proof of this result below.

It is natural to consider what happens when one or more of the axioms is dropped.

Suppose that Hermiticity is not assumed, and that the star product is assumed to be a map $\star :{\mathbb{M}}_{n}\times {\mathbb{M}}_{n}\to {\mathbb{M}}_{n}$, where ${\mathbb{M}}_{n}$ is the complex vector space of all *n*-by-*n* matrices. In this case, as we also show below, the remaining four axioms force the following form of the star product.

### Theorem 4.2

*Let* ${\mathbb{M}}_{n}$ *be the complex vector space of n-by-n matrices. The only functions* $\star :{\mathbb{M}}_{n}\times {\mathbb{M}}_{n}\to {\mathbb{M}}_{n}$ *satisfying convex-bilinearity, product on commuting pairs, product when traced and associativity are x⋆y=xy and x⋆y=yx.*

In this case, there are only two possibilities for the star product, corresponding to straightforward (left and right) matrix multiplication. The possibility of taking the star product to be matrix multiplication was noted in [20] as one way to overcome the shortcomings of the LS star product. We comment on some of the consequences of such a choice in the Conclusion section.

If associativity is dropped, on the other hand, with the other four axioms (including Hermiticity) maintained, then the FJV construction provides an interesting example of a means to construct states over time.

Finally, in appendix B, we also consider the possibility that the result of the star product is a matrix in a larger (or smaller) state space, with

*product on commuting pairs*and

*associativity*, the conclusion of the theorem remains essentially the same.

#### (a) Derivation of main theorems

Key to the derivation of theorems 4.1 and 4.2 is the following result [30], Theorem 1. Using conditions on a function $\star :{\mathbb{M}}_{n}\times {\mathbb{M}}_{n}\to {\mathbb{M}}_{n}$, given as

(A) ‘Associativity’:

*x*⋆(*y*⋆*z*)=(*x*⋆*y*)⋆*z*for all $x,y,z\in {\mathbb{M}}_{n}$;(B) ‘Bilinearity’: (

*λx*)⋆*y*=*λ*(*x*⋆*y*)=*x*⋆(*λy*), (*x*+*y*)⋆*z*=(*x*⋆*z*)+(*y*⋆*z*), and*x*⋆(*y*+*z*)=(*x*⋆*y*)+(*x*⋆*z*) for scalars*λ*and $x,y,z\in {\mathbb{M}}_{n}$;(I) ‘Identity’:

*x*⋆1=*x*for all $x\in {\mathbb{M}}_{n}$, where $1\in {\mathbb{M}}_{n}$ is the identity matrix.(T) ‘Trace’: Tr(

*x*⋆*y*)=Tr(*xy*) for all $x,y\in {\mathbb{M}}_{n}$.(O) ‘Orthogonality’:

*x*⋆*y*=0 when*xy*=*yx*=0,*xx*=*x*, and*yy*=*y*,for $x,y\in {\mathbb{M}}_{n}$ that have rank one.

the result is as follows:

The function $\star :{\mathbb{M}}_{n}\times {\mathbb{M}}_{n}\to {\mathbb{M}}_{n}$ satisfies (A), (B), (I), (T) and (O) if and only if either *x*⋆*y*=*xy* for all $x,y\in {\mathbb{M}}_{n}$, or *x*⋆*y*=*yx* for all $x,y\in {\mathbb{M}}_{n}$.

Here, condition (A) is precisely our associativity condition and (T) is the product when traced axiom. Condition (I) is the identity composition, equation (3.8) above, and condition (O) is that

#### Lemma 4.3

*Let* ${\mathbb{H}}_{n}$ *be the set of n-by-n Hermitian matrices, and* ${\mathbb{M}}_{n}$ *be the set of all n-by-n complex matrices. Any convex-bilinear function* $\star :{\mathbb{H}}_{n}\times {\mathbb{H}}_{n}\to {\mathbb{H}}_{n}$ *satisfying product on commuting pairs, product when traced, and associativity, extends uniquely to a complex bilinear function* $\star :{\mathbb{M}}_{n}\times {\mathbb{M}}_{n}\to {\mathbb{M}}_{n}$, *which also satisfies product on commuting pairs, product when traced and associativity*.

Theorem 4.2 follows from [30], Theorem 1, because if the star product is merely convex-bilinear and satisfies the other requirements, then it follows from the lemma that it must in fact be complex bilinear and still satisfy the other requirements. Theorem 4.1 follows because if there were a function on ${\mathbb{H}}_{n}\times {\mathbb{H}}_{n}$ satisfying the five axioms, then by the lemma it would extend to a function on ${\mathbb{M}}_{n}\times {\mathbb{M}}_{n}$ that satisfies (A), (B), (I), (T) and (O), hence could only be left or right matrix multiplication: but matrix multiplication does not preserve Hermiticity.

We now provide the proof of lemma 4.3.

#### Proof.

We first show that any convex-bilinear function $\star :{\mathbb{H}}_{n}\times {\mathbb{H}}_{n}\to {\mathbb{H}}_{n}$ satisfying the other axioms is in fact bilinear over the reals. It follows from the product on commuting pairs axiom that *x*⋆0=0⋆*x*=0, where 0 here denotes the *n*-by-*n* matrix all of whose components are 0. By considering convex combinations of the form *x*= *rx*′+(1−*r*)0 (for 0≤*r*≤1), or *x*=(1/*r*)*x*′+(1−1/*r*)0 (for *r*≥1), it is then easy to show that (*rx*)⋆*y*=*r*(*x*⋆*y*) for any *r*≥0. Recall that a *conic combination* of ${x}_{i}\in {\mathbb{H}}_{n}$ is $\sum _{i}{r}_{i}{x}_{i}$ for non-negative real weights *r*_{i}. Writing $r=\sum _{i}{r}_{i}$ yields $(\sum _{i}{r}_{i}{x}_{i})\star y=r((1/r\sum _{i}{r}_{i}{x}_{i})\star y)=\sum _{i}{r}_{i}({x}_{i}\star y)$, where the last equality follows from convex-bilinearity. Hence the function ⋆ is linear over conic combinations. Finally, consider a linear combination $x=\sum _{i\in I}{s}_{i}{x}_{i}$ with ${s}_{i}\in \mathbb{R}$. Chop the index set *I* into *I*_{−}={*i*∈*I*|*s*_{i}<0} and *I*_{+}={*i*∈*I*|*s*_{i}≥0}. Then $x+\sum _{i\in {I}_{-}}|{s}_{i}|{x}_{i}=\sum _{i\in {I}_{+}}{s}_{i}{x}_{i}$, where each side of this equation is a conic combination of elements of ${\mathbb{H}}_{n}$. Taking the star product with *y* yields $(x+\sum _{i\in {I}_{-}}|{s}_{i}|{x}_{i})\star y=(\sum _{i\in {I}_{+}}{s}_{i}{x}_{i})\star y$, which after using conic linearity and rearranging gives $x\star y=\sum _{i}{s}_{i}({x}_{i}\star y)$. Similar reasoning applied to the second argument yields bilinearity over the reals, as claimed.

We now prove that any product $\star :{\mathbb{H}}_{n}\times {\mathbb{H}}_{n}\to {\mathbb{H}}_{n}$, which is bilinear over the reals, and satisfies the remaining axioms, extends to a product on ${\mathbb{M}}_{n}$, which is bilinear over the complex numbers, and still satisfies the remaining axioms. Recall that a matrix $x\in {\mathbb{M}}_{n}$ is Hermitian if *x*^{†}=*x*, and anti-Hermitian if *x*^{†}=−*x*; there is a bijection *x*↦*ix* between Hermitian and anti-Hermitian matrices. Any matrix $x\in {\mathbb{M}}_{n}$ equals *x*=*x*_{h}+*x*_{a} for a unique Hermitian ${x}_{h}\in {\mathbb{H}}_{n}$ and anti-Hermitian *x*_{a}, given by ${x}_{h}={\textstyle \frac{1}{2}}(x+{x}^{\u2020})$ and ${x}_{a}={\textstyle \frac{1}{2}}(x-{x}^{\u2020})$. Define $\star :{\mathbb{M}}_{n}\times {\mathbb{M}}_{n}\to {\mathbb{M}}_{n}$ by

*λ*=

*α*+

*iβ*for $\alpha ,\beta \in \mathbb{R}$. Observe that (

*λx*)

_{h}=

*αx*

_{h}+

*βix*

_{a}, and

*i*(

*λx*)

_{a}=

*αix*

_{a}−

*βx*

_{h}. Then

*x*⋆(

*λy*)=

*λ*(

*x*⋆

*y*). Using (

*x*+

*y*)

_{h}=

*x*

_{h}+

*y*

_{h}and (

*x*+

*y*)

_{a}=

*x*

_{a}+

*y*

_{a}gives (

*x*+

*y*)⋆

*z*=

*x*⋆

*z*+

*y*⋆

*z*similarly easily.

Finally, consider the condition of being product on commuting pairs.

First, consider the case where $x\in {\mathbb{M}}_{n}$ is normal, that is, *xx*^{†}=*x*^{†}*x*. By assumption *xy*=*yx*, and so *x*^{†}*y*^{†}=(*yx*)^{†}=(*xy*)^{†}=*y*^{†}*x*^{†}, too. By Fuglede's theorem [31], also *xy*^{†}=*y*^{†}*x*, and so *x*^{†}*y*=*yx*^{†}. Therefore,

*x*

_{h}

*y*

_{a}=

*y*

_{a}

*x*

_{h},

*x*

_{a}

*y*

_{h}=

*y*

_{h}

*x*

_{a}and

*x*

_{a}

*y*

_{a}=

*y*

_{a}

*x*

_{a}. It now follows that

*x*

_{h}⋆

*y*

_{h}=

*x*

_{h}

*y*

_{h},

*x*

_{h}⋆

*iy*

_{a}=

*ix*

_{h}

*y*

_{a},

*ix*

_{a}⋆

*y*

_{h}=

*ix*

_{a}

*y*

_{h}and

*ix*

_{a}⋆

*iy*

_{a}=−

*x*

_{a}

*y*

_{a}. But this, together with the bilinearity of $\star :{\mathbb{M}}_{n}\times {\mathbb{M}}_{n}\to {\mathbb{M}}_{n}$, implies

Now, we consider the case of a general $x\in {\mathbb{M}}_{n}$ (not necessarily normal). Because both hermitian and anti-hermitian matrices are normal, an arbitrary matrix $x\in {\mathbb{M}}_{n}$ can be written as a sum of two normal matrices *x*=*x*_{h}+*x*_{a}. Combining the result for normal matrices with bilinearity of the extended function $\star :{\mathbb{M}}_{n}\times {\mathbb{M}}_{n}\to {\mathbb{M}}_{n}$, we infer that *x*⋆*y*=*xy* for any $y\in {\mathbb{M}}_{n}$.

This concludes the proof of lemma 4.3. □

### 5. Conclusion

We have considered methods for constructing quantum states over time, under the assumption that the state over time is an operator on the tensor product of Hilbert spaces associated with the temporally localized subsystems. The motivation was to explore the possibilities for a quantum treatment of time-like separated systems that is as close as possible to the standard quantum treatment of composite systems at a single time. In particular, given the well-developed analogy between density matrices and probability distributions—and the fact that classical joint probabilities can be defined regardless of the spatio-temporal relationships between a set of events—it is interesting to see whether this analogy can be extended, or whether it is fundamentally limited in scope by an asymmetric treatment of space and time by quantum theory. We focused on simple cases of spatio-temporally localized quantum systems evolving through a sequence of quantum channels. We identified a set of physically motivated criteria for such a construction, and theorem 4.1 shows that there is no construction that satisfies all of the criteria.

One possible conclusion that may be drawn from theorem 4.1, then, is that there is something fundamentally misguided in attempting to treat quantum systems over time in the same manner as composite systems at a single time. Interpreted thus, theorem 4.1 might be seen, for example, as lending support to the view that a quantum system localized in time should be described with two Hilbert spaces, representing an incoming and outgoing state, as per the multi-time [6–9], quantum combs [10–12] and process matrices [15,16,21] formalisms.

But one should also ask if there is a way that the apparent no-go result of theorem 4.1 can be avoided. The most obvious way is to question whether the axioms that we have identified are themselves reasonable. The axioms of *convex-bilinearity* and of *product when traced* seem to be essential if the joint state is to play a role that is analogous to the role played by the joint state of composite systems at a given time. Reproducing classical joint probabilities when the initial state and the channel states all commute also seems essential, at least if the quantum states over time are to be seen as generalizations of classical probability distributions in the way that standard quantum density matrices are seen as generalizations of classical probability distributions. However, it is reasonable to judge a pairing of initial state and channel state to be classical only when they are jointly diagonalizable *in a basis that is a tensor product of a fixed basis for each system* rather than simply jointly diagonalizable. Therefore, the assumption that the joint state should be equivalent to a joint probability distribution for such pairs only warrants that the star product should reduce to the matrix product on certain pairs of commuting operators, rather than all pairs of commuting operators, as the condition asserts. Hence there may be some wiggle room in constructing joint states that reproduce classical joint probabilities in the desired sense but which do not satisfy the axiom.

As noted in the previous section, if the state over time is not assumed to be Hermitian, then theorem 4.2 shows that matrix multiplication is the only construction that satisfies the remaining axioms. We also noted that if one does not assume associativity, then the FJV construction satisfies the other four axioms. Furthermore, although the FJV star product is not associative for arbitrary triples of operators, we saw that on the restricted domain of triples of operators corresponding to a sequence of channel states, it *is* associative. A disadvantage of the FJV construction as stated, however, is that it is defined only for collections of qubits. It can be applied to systems whose dimension is a power of 2 by regarding them as composite systems made up of qubits, but then there is an arbitrariness in the choice of how to factorize the larger Hilbert space into component systems. There is also a natural generalization of the proposal from a qubit to a system of prime dimension, but how to generalize the proposal to systems of arbitrary dimension remains ambiguous. Note that there is an interesting relationship between the FJV construction and that corresponding to matrix multiplication, at least for a single qubit considered at two times: the joint state *ρ*^{(FJV)}_{AB} is a Hermitian operator whose eigenspectrum coincides with the positive and real eigenspectrum of *ρ*^{(MM)}_{AB}=*E*_{B|A}*ρ*_{A}; that is, *ρ*^{(FJV)}_{AB} can be obtained from *ρ*^{(MM)}_{AB} by removing the latter's imaginary eigenvalues.

In all of the proposals for defining a state over time that we have considered, the set of operators that can describe a joint state of a system considered at two times is distinct from the set of operators that can describe a state of two copies of the system at a given time. This contrasts with classical probability theory, where the possible joint probability distributions on a pair of systems is the same in the two cases. This suggests that even if a sensible notion of a state over time can be defined, quantum theory still manifests an asymmetry between space and time that is not present in classical probability theory (analogous, perhaps, to the fact that relativity theory still distinguishes space and time through the signature of the metric). Indeed, this conclusion is supported by the results of [19], which demonstrate that in quantum theory, unlike classical theories, the set of correlations that can be generated by two systems that are related as cause and effect is different from the set that can be generated by two systems that are both effects of a common cause. Further investigations into operationally motivated definitions of states over time are likely to shed additional light on the precise nature of this asymmetry.

Finally, in seeking to define a state for a composite over time, we have here considered to what extent the formal properties of such a state could resemble those of a state for a composite at a given time. It would also be interesting to approach the problem from a purely operational point of view, where the aim is to define a joint state by reference to predictions for experimental interventions on systems at different times. For instance, it would be interesting to consider the relationships between the various constructions considered here and the predictions for the outcomes of weak measurements performed across multiple time steps [32].

### Data accessibility

This work does not have any experimental data.

### Authors' contributions

The initial question came out of discussions between J.B. and R.W.S. The original axioms resulted from discussions between J.B. and M.F.P., and were further developed by R.W.S. and D.H. C.H. proved the main theorems. D.H. and M.F.P. developed the FJV construction as a star product. D.H. wrote most of the article.

### Competing interests

We declare we have no competing interests.

### Funding

This research was supported in part by Perimeter Institute for Theoretical Physics. Research at Perimeter Institute is supported by the Government of Canada through the Department of Innovation, Science and Economic Development Canada and by the Province of Ontario through the Ministry of Research, Innovation and Science. D.H. and J.B. acknowledge support from an FQXi Large Grant, and from the NQIT Quantum Hub. D.H. was also supported by EPSRC grant no. EP/L022303/1. This project/publication was made possible in part through the support of a grant from the John Templeton Foundation. The opinions expressed in this publication are those of the author(s) and do not necessarily reflect the views of the John Templeton Foundation. C.H. was supported by EPSRC Fellowship EP/L002388/1.

## Appendix A. The FJV construction as a star product

In this section, we show that for a qubit at two times, the Fitzsimons–Jones–Vedral construction (2.5) corresponds to the Jordan product, ${\rho}_{AB}=\frac{1}{2}({\rho}_{A}{E}_{B|A}+{E}_{B|A}{\rho}_{A})$. Recall that the FJV construction defines a state of two systems as

*σ*

_{i}⊗

*σ*

_{j}〉 is read for a state over time as the expectation value of the product of outcomes of a projective measurement of

*σ*

_{i}followed by a measurement of

*σ*

_{j}(where, after the first measurement,

*A*is in the eigenstate of

*σ*

_{i}corresponding to the outcome of that measurement).

Using the fact that the state at the output of the channel is given by Tr_{A}(*E*_{B|A}*ρ*_{A}), we have

*A*after projection onto the +1 eigenspace of

*σ*

_{i}, and similarly for ${\rho}_{A}^{i-}$, i.e. ${\rho}_{A}^{i\pm}\equiv {P}^{i\pm}{\rho}_{A}{P}^{i\pm}$ where ${P}^{i\pm}\equiv (\mathrm{\U0001d7d9}\pm {\sigma}_{i})/2$.

Noting that

Equation (A 1) implies that

## Appendix B. Larger state spaces

One could consider functions $\star :{\mathbb{M}}_{n}\times {\mathbb{M}}_{n}\to {\mathbb{M}}_{m}$ for *m* possibly different from *n*. The condition of convex-bilinearity then still makes sense, but associativity and being product on commuting pairs do not. The most natural solution is to require that there exists a linear function $p:{\mathbb{M}}_{m}\to {\mathbb{M}}_{n}$ and adapt the criteria accordingly:

–

*associativity*:*x*⋆*p*(*y*⋆*z*)=*p*(*x*⋆*y*)⋆*z*;–

*product on commuting pairs*: if*xy*=*yx*, then*p*(*x*⋆*y*)=*xy*;–

*product when traced*: Tr(*p*(*x*⋆*y*))=Tr(*xy*).

But then one easily sees that $p\circ \star :{\mathbb{M}}_{n}\times {\mathbb{M}}_{n}\to {\mathbb{M}}_{n}$ satisfies the original criteria for a star product. So, by theorem 4.2, either *p*(*x*⋆*y*)=*xy* for all $x,y\in {\mathbb{M}}_{n}$ or *p*(*x*⋆*y*)=*yx*. Hence $p(x\star \mathrm{\U0001d7d9})=x$ and it follows that *rank*(*p*)≥*n*^{2} and so *m*≥*n*. Thus, the function $j:{\mathbb{M}}_{n}\to {\mathbb{M}}_{m}$ given by $j(x)=x\star \mathrm{\U0001d7d9}$ satisfies *p*(*j*(*x*))=*x*. But then Tr(*x*⋆*y*)=Tr(*j*(*xy*)), so that *x*⋆*y* and *j*(*xy*) result in the same probabilistic predictions when measuring a state. In other words, there is nothing to be gained by moving from ${\mathbb{M}}_{n}$ to ${\mathbb{M}}_{m}$ for *m*>*n*.

### References

- 1
Feynman RP . 1948 Space-time approach to non-relativistic quantum mechanics.**Rev. Mod. Phys.**, 367. (doi:10.1103/RevModPhys.20.367) Crossref, ISI, Google Scholar**20** - 2
Sorkin RD . 1997 Forks in the road, on the way to quantum gravity.**Int. J. Theor. Phys.**, 2759–2781. (doi:10.1007/BF02435709) Crossref, ISI, Google Scholar**36** - 3
Griffiths RB . 1984 Consistent histories and the interpretation of quantum mechanics.**J. Stat. Phys.**, 219–272. (doi:10.1007/BF01015734) Crossref, ISI, Google Scholar**36** - 4
Omnes R . 1988 Logical reformulation of quantum mechanics. I. Foundations.**J. Stat. Phys.**, 893–932. Crossref, ISI, Google Scholar**53** - 5
Gell-Mann M, Hartle JB . 1990 Quantum mechanics in the light of quantum cosmology. In*Complexity, entropy and the physics of information*(ed. WH Zurek), vol. 8. Redwood City, CA: Addison-Wesley. Google Scholar - 6
Aharonov Y, Vaidman L . 2007 The two-state vector formalism: an updated review. In*Time in quantum mechanics*(eds J Muga, R Mayato, I Egusquiza), pp. 399–447. Berlin, Germany: Springer. Google Scholar - 7
Aharonov Y, Popescu S, Tollaksen J, Vaidman L . 2009 Multiple-time states and multiple-time measurements in quantum mechanics.**Phys. Rev. A**, 052110. (doi:10.1103/PhysRevA.79.052110) Crossref, ISI, Google Scholar**79** - 8
Aharonov Y, Popescu S, Tollaksen J . 2013 Each instant of time a new universe. In*Quantum theory: a two-time success story*(eds DC Struppa, JM Tollaksen), p. 21. Yakir Aharonov Festschrift. Milan, Italy: Springer. Google Scholar - 9
Silva R, Guryanova Y, Brunner N, Linden N, Short AJ, Popescu S . 2014 Pre- and postselected quantum states: density matrices, tomography, and kraus operators.**Phys. Rev. A**, 012121. (doi:10.1103/PhysRevA.89.012121) Crossref, ISI, Google Scholar**89** - 10
Chiribella G, D'Ariano GM, Perinotti P . 2009 Theoretical framework for quantum networks.**Phys. Rev. A**, 022339. (doi:10.1103/PhysRevA.80.022339) Crossref, ISI, Google Scholar**80** - 11
Chiribella G . 2012 Perfect discrimination of no-signalling channels via quantum superposition of causal structures.**Phys. Rev. A**, 040301. (doi:10.1103/PhysRevA.86.040301) Crossref, ISI, Google Scholar**86** - 12
Chiribella G, D'Ariano GM, Perinotti P, Valiron B . 2013 Quantum computations without definite causal structure.**Phys. Rev. A**, 022318. (doi:10.1103/PhysRevA.88.022318) Crossref, ISI, Google Scholar**88** - 13
Gutoski G, Watrous. J . 2007 Toward a general theory of quantum games. In*STOC '07 Proc. 39th Annual ACM Symp. on Theory of Computing*, pp. 565–574. New York, NY: ACM. Google Scholar - 14
Oreshkov O, Cerf NJ . 2016 Operational quantum theory without predefined time.**New. J. Phys.**, 073037. (doi:10.1088/1367-2630/18/7/073037) Crossref, ISI, Google Scholar**18** - 15
Oreshkov O, Costa F, Brukner C . 2012 Quantum correlations with no causal order.**Nat. Commun.**, 1092. (doi:10.1038/ncomms2076) Crossref, PubMed, ISI, Google Scholar**3** - 16
Costa F, Shrapnel S . 2016 Quantum causal modelling.**New. J. Phys.**, 063032. (doi:10.1088/1367-2630/18/6/063032) Crossref, ISI, Google Scholar**18** - 17
Hardy L . 2007 Towards quantum gravity: a framework for probabilistic theories with non-fixed causal structure.**J. Phys. A: Math. Theor.**, 3081. (doi:10.1088/1751-8113/40/12/S12) Crossref, ISI, Google Scholar**40** - 18
Hardy L . 2012 The operator tensor formulation of quantum theory.**Phil. Trans. R. Soc. A**, 3385–3417. (doi:10.1098/rsta.2011.0326) Link, ISI, Google Scholar**370** - 19
Ried K, Agnew M, Vermeyden L, Janzing D, Spekkens RW, Resch KJ . 2015 A quantum advantage for inferring causal structure.**Nat. Phys.**, 414–420. (doi:10.1038/nphys3266) Crossref, ISI, Google Scholar**11** - 20
Leifer MS, Spekkens RW . 2013 Towards a formulation of quantum theory as a causality neutral theory of Bayesian inference.**Phys. Rev. A**, 052130. (doi:10.1103/PhysRevA.88.052130) Crossref, ISI, Google Scholar**88** - 21
Araújo M, Costa F, Brukner Č . 2014 Computational advantage from quantum-controlled ordering of gates.**Phys. Rev. Lett.**, 250402. (doi:10.1103/PhysRevLett.113.250402) Crossref, PubMed, ISI, Google Scholar**113** - 22
Fizsimons J, Jones JA, Vedral V . 2013 Quantum correlations which imply causation. (http://arxiv.org/abs/1302.2731) Google Scholar - 23
Jamiołkowski A . 1972 Linear transformations which preserve trace and positive semidefiniteness of operators.**Rep. Math. Phys.**, 275–278. (doi:10.1016/0034-4877(72)90011-0) Crossref, Google Scholar**3** - 24
Choi M-D . 1975 Completely positive linear maps on complex matrices.**Linear. Algebra. Appl.**, 285–290. (doi:10.1016/0024-3795(75)90075-0) Crossref, ISI, Google Scholar**10** - 25
Leifer MS . 2006 Quantum dynamics as an analog of conditional probability.**Phys. Rev. A**, 042310. (doi:10.1103/PhysRevA.74.042310) Crossref, ISI, Google Scholar**74** - 26
Leifer MS, Poulin D . 2008 Quantum graphical models and belief propagation.**Ann. Phys. (N. Y)**, 1899–1946. (doi:10.1016/j.aop.2007.10.001) Crossref, ISI, Google Scholar**323** - 27
Wooters WK . 1987 A Wigner-function formulation of finite-state quantum mechanics.**Ann. Phys. (N. Y)**, 1–21. (doi:10.1016/0003-4916(87)90176-X) Crossref, ISI, Google Scholar**176** - 28
Gross D . 2006 Hudson's theorem for finite-dimensional quantum systems.**J. Math. Phys.**, 122107. (doi:10.1063/1.2393152) Crossref, ISI, Google Scholar**47** - 29
Pearl J . 2000**Causality: models, reasoning, and inference**. New York, NY: Cambridge University Press. Google Scholar - 30
Heunen C, Horsman C . 2013 Matrix multiplication is determined by orthogonality and trace.**Linear Algebra Appl.**, 4130–4134. (doi:10.1016/j.laa.2013.09.039) Crossref, ISI, Google Scholar**439** - 31
Fuglede B . 1950 A commutativity theorem for normal operators.**Proc. Natl Acad. Sci. USA**, 35–40. (doi:10.1073/pnas.36.1.35) Crossref, PubMed, ISI, Google Scholar**36** - 32
Marcovitch S, Reznik B . 2011 Structural unification of space and time correlations in quantum theory. (http://arxiv.org/abs/1103.2557) Google Scholar