A scheme for performing strong and weak sequential measurements of non-commuting observables

Quantum systems usually travel a multitude of different paths when evolving through time from an initial to a final state. In general, the possible paths will depend on the future and past boundary conditions, as well as the system’s dynamics. We present a gedanken experiment where a single system apparently follows mutually exclusive paths simultaneously, each with probability one, depending on which measurement was performed. This experiment involves the measurement of observables that do not correspond to Hermitian operators. Our main result is a scheme for measuring these operators. The scheme is based on the erasure protocol [Brodutch and Cohen (Phys. Rev. Lett. 116:070404, 2016)] and allows a wide range of sequential measurements at both the weak and strong limits. At the weak limit the back action of the measurement cannot be used to account for the surprising behaviour and the resulting weak values provide a consistent yet strange account of the system’s past.


Introduction
The joint measurement of sequential observables in quantum mechanics is not canonically defined. A sequential measurement of two observables at different times can lead to incompatible results that depend on how the experiment is carried out [1][2][3]. When the observables are compatible (i.e., they commute) the results of a measurement can depend on the context of the measurement, leading to the well known contextuality and non locality theorems [4][5][6][7][8][9]. When the observables are incompatible, the results of measurements can be correlated in a non-trivial manner that leads to Leggett-Garg inequalities [10][11][12][13][14][15][16][17][18] for standard invasive measurements and to other predictions when the measurements are less invasive [2,19].
Questions regarding sequential measurements are even more significant in systems with both past and future boundary conditions. The observable quantities that depend on past (pre-selected) and future (post-selected) bound- Fig. 1 The erasure scheme for the sequential operator B A To make a variable strength measurement of the operator B A we first use an ancilla to couple strongly to A at time t 1 and then use the meter to couple to both ancilla and the system to measure the joint observable. The measurement is completed by erasing the result recorded in the ancilla, a probabilistic procedure. The strength of the measurement can be selected by changing the strength of the interaction at t 2

Measurements in the TSVF
The TSVF is a time symmetric formulation of quantum mechanics, that can be used to calculate probabilities for measurements between a pre-selection |ψ and a post-selection φ|. It usually helps in acquiring better understanding of a system between two projective measurements.
The following setup is relevant for the rest of the paper: At time t 0 a system S is prepared in the state |ψ . Later, at t f the system undergoes a projective rank-1 measurement and is found to be in the state |φ . For clarity of notation we represent this post-selected state as φ|. Since quantum mechanics is time-symmetric, it is possible to think of the state |ψ as evolving forward in time and on the state φ| as evolving backward in time. We assume the free evolution is the identity, so if there are no interactions at intermediate times, the state of the system is best described by the two-state vector φ| |ψ [21].
We want to predict the results of measurements made at the intermediate time interval t 1 to t 2 , where t 0 < t 1 < t 2 < t f . We consider two types of measurements, a strong (projective) measurement and a weak measurement. In this section we introduce the standard formalism for these two types of measurements. In the following sections we will extend them to the case of sequential measurements using two techniques: Modular measurements [1,46] and the erasure method [32].

Strong measurements
We use the term strong measurement to refer to a measurement that gives an unambiguous result and projects S onto the corresponding subspace. Let A be a measurement operator characterized by the projectors A k such that A = k a k A k . Given the two-state vector φ| |ψ , the probability that an intermediate strong measurement of A yields a result corresponding to A k is given by the ABL formula [20,21]: We usually think of A as an Hermitian operator and A k are projectors onto the degenerate subspaces of this operator. However, as we will see below (Sect. 3.4), the formula can be extended to a more general case where A is a quantum channel and A k are the corresponding Kraus operators (see also [47]).

Weak measurements
For a system described by the two-state vector φ| |ψ , the weak measurement of an operator A produces an outcome {A} w = φ|A|ψ φ|ψ called the weak value of A [21,23]. The standard method for implementing this weak measurement is to couple the system to a meter with momentum P and position Q via the von Neumann interaction Hamiltonian H I = f (t)AP. This way, the shift in the pointer's position will be proportional to the weak value. The weakness (or strength) of the measurement is usually modified by varying the strength of the interaction g = f (t)dt and/or by changing the variance σ 2 of the initial state of the measuring device. In general, the measurement process between t 1 and t 2 is an interaction between the system S and a meter M, such that at the end of the process the (change in the) meter's state corresponds to the measurement outcome. While the state of the system can be arbitrary, the state of the meter is specified according to the desired properties of the measurement. To avoid confusion, we use the subscript w for a meter in a weak measurement. The meter's initial state is denoted by |0 M w .
The weak measurement is a physical process parameterized by the parameter g. It can be modeled as a map taking a pure product system-meter state ρ SM w 0 = |ψ ψ| ⊗ |0 0| to a joint state W g (ρ SM w 0 ), such that the following properties hold at g → 0 [32].
The probability of post-selecting a state φ|, given by is unaffected by the measurement up to terms of order g 2 (2) Property 2 (Weak potential) After post-selection the meter state |0 M w is shifted by a value proportional to the 'weak value' The measurement must also have a strong limit at g = 1 with the following property: ) is a strong (Lüders) measurement.
While most weak measurements have been realized via a coupling to continuous variables, weak values appearing in couplings to qubit meters have been given increasing attention recently [46,[48][49][50]. In the current work we shall analyze the discrete, qubit meter case. Properties 1, 2 are valid also in the case of coupling to a discrete variable, as long as the coupling is weak enough [49]. Property 3 is easily achieved when the corresponding observables are dichotomic (as in the examples that follow), however, one has to be careful in more general cases.
The case of dichotomic observables is of special interest and has an interesting property in the case of weak measurements [44].
This property follows simply from Property 2 and the ABL rule, Eq. (1).

Sequential measurements
As the first example of non-trivial sequential measurement Aharonov and Albert [1] described an experiment where, using a modular meter, it is possible to perform a correlation-type measurement of a two-time operator σ t 2 x σ t 1 z (that Fig. 2 The basic setup. a Two sequential measurements in the Z and X bases are performed on a pre-/post-selected system | |ψ . b The four possible measurements outcomes are translated into four paths is, a measurement of σ z at time t 1 followed by a measurement of σ x at time t 2 ). Such a measurement can give qualitatively and quantitatively different results than correlations between a sequence of two measurements made on two devices, the first for σ z at time t 1 and the second for σ x at time t 2 . In a more recent approach [27], the path of a photon through a double interferometer was shown to exhibit strange behavior when questions about its location at a given point in time were asked. These were resolved using weak values of sequential observables. The weak values were, however, related to an indirect measurement procedure. If the observable at time t 1 is A 1 and the observable at time t 2 is A 2 , then it is possible to extract the sequential weak value {A 2 A 1 } w by performing a weak measurements of A 1 using a meter q 1 and of A 2 using q 2 . Using an equality initially derived by Resch and Steinberg [51], it can be shown that So the weak value can be extracted although the procedure is not a weak measurement, i.e. it does not have Properties 1-3 above [32,46,52,53].
Below we present a similar example of sequential measurements and show a scheme for performing these measurements directly. We will use this example to illustrate some of the features of weak and strong sequential measurements. To avoid confusion we use the notation P(A) to denote the probability of outcome corresponding to the operator A in a strong measurement and {A} w to denote the weak value of the operator A.

Setup
Our example is based on the following scenario (see Fig. 2a): a system S is prepared in the state and post-selected in the state where we use the standard notation σ z |0 = |0 , σ z |1 = − |1 and σ x |± = ± |± .
These correspond to the following four operators: that add up to the identity A + B + C + D = I. We also note that The four transition amplitudes are

A sequence of two strong measurements
If we use a sequence of two projective measurements (i.e. we measure in the Z basis, read out the result and then measure in the X basis) we get the probabilities corresponding to the four possible paths in Fig. 2b. This type of measurement represents the question 'which path did the system follow?'. We now turn our attention to two other ways of measuring the observables A, B, C, D . The first of these (Fig. 3a) is a measurement that can only distinguish between pairs of the results in Eq. (8), by coupling to a single meter. In the final setup ( Fig.  Fig. 3 Different measurement schemes. a Sequential measurement of the modular values M E = A + D and , each measurement distinguishes one result from the other three, marking a distinct path from |ψ to |φ . Our choice of pre-and post-selection is inspired by the three box paradox [44] and produces a situation where we can deterministically observe the system going through two mutually exclusive paths.

Sequential measurement using modular values
In the spirit of [1] we will now study a measurement of a sequential observable σ Z X = σ t 2 x σ t 1 z . A result of +1 corresponds to both (0, +) and (1, −) in the Z , X measurements, while a −1 result corresponds to (0, −) and (1, +). This is the same as measuring the operators M E = A + D and M O = B + C, respectively (see Fig. 3a).

The probabilities for a strong measurement
The probabilities for a strong measurement are The difference between this measurement and a scheme using two separate measurements is apparent when we consider deterministic outcomes. For Eq. (19) we need to set θ = 0, π/2 and φ = 0, π/2. In contrast, for the modular measurement, we need to set either θ = −φ to obtain b + c = 0 or θ = π/2 − φ to obtain a + d = 0.

The corresponding weak values
To calculate the weak values we can use the fact that weak values are additive, hence We can now see the difference between the four outcome measurement of Eq. (19) and the two outcome measurement of Eq. (20). To experimentally measure the weak and strong values for this setup it is possible to use the modular values method of Kedem and Vaidman [46], or the erasure method [32] (see Sect. 3.4.3 below).

The deterministic path paradox
In the second example the measurement device is set in such a way (see Fig. 1 above and Sect. 3.4.3 below) that it clicks only for one of the four possible outcomes, say C (see Fig. 3b). In this case we have four different measurement settings: A, where a result of 1 ( or 'click') corresponds to path A while a 0 (or 'no-click') corresponds to not going through path A and likewise for B, C and D. Using the ABL rule we arrive at the following probabilities for a 'click' in each experiment.
These probabilities are very different from the sequence of individual measurements described in Eq. (19).

Derivation of the 'paradox'
We want to find θ and φ such that B and C click with certainty, if the paths are projectively measured (i.e with measurements B and C, respectively). Using Eq. (24) the conditions for P B (B) = P C (C) = 1 are Subtracting these we get c = b, which translates to cot(θ ) = cot(φ). Using this, and dividing the first equation by sin θ sin φ = 0, we get and the solution is The apparent paradox can be described as: 1. The paths B and C are mutually exclusive. 2. The system traveled through path B with certainty. 3. The system traveled through path C with certainty.
However, statements 2 and 3 are counterfactual, we cannot ascertain them simultaneously. Or can we?

The weak value of a path
As in the three box paradox [44], weak measurements allow us to make sense of counterfactual statements. While the measurements B and C cannot be performed simultaneously, their weak counterparts can. The results are consistent with the apparent 'paradox'. For general θ and φ the weak values are: For the deterministic case, using the solution (27), we have As expected the total number of particles traversing the four paths is 1. Moreover, the particle traversed both paths B and C with certainty. These results could have been predicted using Property 4. Furthermore, we could have used this fact to arrive at Eq. (26).
This result raises a question regrading the 'past of a quantum particle' [36]. If the particle is understood to be at every place it left a 'weak trace', then we can conclude it traveled through all four paths. However, if we understand the weak value as an effective weak potential [26] then a weak coupling to the particle in B and C will result in the effective H I = +1, while coupling to the particle in A/D will be negative.
This result accords well with the past results of the three box paradox [44,54,55], the negative pressure paradox [21], the Cheshire cat [56,57], Hardy's paradox [5,6,58] and others, where weak measurements reveal a curious behavior of pre-and post-selected systems.

The erasure method for measuring a distinct path
Weak values provide an interesting perspective for the situation described above. However, their physical meaning is lost, if there is no corresponding weak measurement. To perform the desired sequential weak measurement we use the erasure method introduced in [32]. The erasure method uses two basic principles.
1. It is possible to undo (erase) the effect of a von Neumann coupling by making an appropriate measurement on the meter. 2. It is possible to make a von Neumann measurement at any strength using the following procedure • Coupling S to a meter M s with the standard (strong) von Neumann coupling.
• Coupling M s to a second meter M w using a von Neumann coupling of the desired strength.
• Erasing the result on M s .
To show how this is done in our example (see also Fig. 1) we define three types of unitary operators C N OT , C i j R(g) and R kl i j (g) as follows: where the index i takes the values +, − and j takes the values 0, 1. The measurement procedure requires an ancilla and a meter. It works in the following way: first, we perform a C N OT on the system and ancilla to measure the system in the Z basis. Next, we perform C i j R(g) with (i, j) set to be the inverted sequential operator we want, for example (i, j) = (+, 0) for the measurement A of A. This rotates the meter by g around σ x when the system follows the measured path. Finally, we erase the first measurement by post-selecting the ancilla in the |+ state. Failing this final step is the same as a unitary operation on the initial state. It is possible to undo this unitary, but the cost is a change of the measurement operator. We will follow the circuit for an arbitrary pre-and post-selection. |ψ = α |0 + β |1 , |φ = γ |+ + δ |− 1. Initial state: 1-System; 2-Ancilla; 3-Meter 2. CNOT on 1, 2 3. C i j R(g) | 2 4. Erasure: we measure σ x on 2 and discard this subsystem. If the result is |+ we ('succeed' and) get the unnormalized state If the result is |− we ('fail' and) get the unnormalized state We will continue the derivation for successful erasure (Eq. 34) and only later return to the failed case.
Returning to the erasure step, failure would mean the post-erasure state is proportional to Eq. (35). This state corresponds to Eq. (34) with the change β → −β which can be interpreted as a σ z operation on the initial state.
We can 'undo' this operation by applying σ z to the system to get Comparing with Eq. (34) we have the same state up to a change in the measurement operator, switching A ↔ C and B ↔ D. For a weak measurement we can apply A, B, C and D simultaneously, so that the change in measurement operator can be corrected in post-processing.
Similarly, the above method can be used for performing non-local weak measurements of the peculiar weak values discussed in [59].
Finally, it is possible to think of other interesting measurement operators that can create stronger versions of the paradox above. However, it is unclear if we can assign physical meaning to those operators using the erasure scheme or other more complex schemes.

Conclusions and outlook
We presented an extension of the erasure scheme [32] that can be used to perform special types of sequential measurements. To demonstrate the potential of this scheme we described three gedanken experiments that showcase the difference between a sequence of measurements and sequential measurements. First, we ask 'which path did the system go through?', second, we ask 'did the system go through A/D or through B/C?'. Finally, for each path we ask 'did the system go through this path?'. The answers to these three questions do not agree, and moreover, the answer to the last question suggests that the system exhibits surprising behavior by deterministically traversing two routes at the same time. It is possible to attribute the strange behavior to the fact that the questions are mutually exclusive, however, the results are consistent even when weak measurements are used. Our scheme shows that these theoretical predictions can be tested experimentally.
Optical implementations of weak measurement experiments have become common in the last decade [13,36,54,58,60]. It is convenient to think of an optical version of the paradox using a double interferometer as in [27]. The advantage of this setup is that the path described by the Z , X measurements corresponds to an actual path along the arms of the interferometer. However, from an experimental point of view an optical test of this experiment may be hard to realize due to the large number of qubits required (i.e. 3-5 qubits 1 ). The mathematical description of our protocol in terms of Pauli operators (rather than optical paths) was deliberately chosen with the view that these experiments are more likely to be performed in other platforms such as liquid state NMR [50] or superconducting qubits (where the sequential weak measurements on independent meters are performed routinely) [16,17,61]. References [16,17] experimentally demonstrate a qubit meter in a superconducting system with explicitly variable strength, with the latter experiment being realized on a chip with 9 independently controllable qubits, thus possibly allowing the proposed circuit to be implemented in a similar system.
The experiment described here demonstrates only one potential use of the erasure scheme in the context of sequential or non-Hermitian observables. The study of these types of observables might have been overlooked in the past due to the lack of a method to make sequential measurements and partly due to the ambiguity in their definition as described in the three different experiments above. The erasure scheme opens up the possibility for making a large class of such sequential measurements, as well as the possibility for further exploration of their significance. They may lead to two complementary research directions, one related to the foundations of quantum mechanics and one related to development of practical methods in quantum information processing and quantum control.
On the foundational side, sequential measurements of the type described in Sect. 3.4 offer a different view of the evolution of a quantum system and provide a unique approach for probing the causal structure of quantum mechanics. They can be combined with the non-local measurements described in [32] to probe general space-time observables (e.g. in the cases of quantum steering [62,63] or summoning [64,65]) in a way that has not been explored yet in much detail (see also [31] for one implication of these sequential measurements). General space-time observables are interesting at both the weak and strong regimes. It is, however, unclear if the erasure scheme can be somehow modified to perform a continuous sequential measurement of this type which could provide answers to more general questions.
Regarding quantum control, we have already described how the non-local version of the erasure scheme can be used for performing a variable strength Toffoli gate [32]. Examples of this kind are easy to find when discussing non-local interactions since they are very natural in quantum protocols. The role of sequential observables in quantum information processing and quantum control methods have not been explored much, but we believe that our experimentally feasible method is capable of opening up new opportunities in this direction.