arxiv: v1 [quant-ph] 25 Jan 2018

Transcription

1 Self-testing quantum states and measurements in the prepare-and-measure scenario Armin Tavakoli, Jędrzej Kaniewski, Tamás Vértesi, 3 Denis Rosset,, 5 and Nicolas Brunner Département de Physique Appliquée, Université de Genève, CH- Genève, Switzerland QMATH, Department of Mathematical Sciences, University of Copenhagen, Universitetsparken 5, 00 Copenhagen, Denmark 3 Institute for Nuclear Research, Hungarian Academy of Sciences, P.O. Box 5, 00 Debrecen, Hungary Perimeter Institute for Theoretical Physics, 3 Caroline St. N, Waterloo, Ontario, Canada, NL Y5 5 Institute for Quantum Optics and Quantum Information (IQOQI, Boltzmangasse 3, 090 Vienna, Austria The goal of self-testing is to characterize an a priori unknown quantum system based solely on measurement statistics, i.e. using an uncharacterized measurement device. Here we develop self-testing methods for quantum prepare-and-measure experiments, thus not necessarily relying on entanglement and/or violation of a Bell inequality. We present noise-robust techniques for self-testing sets of quantum states and measurements, assuming an upper bound on the Hilbert space dimension. We discuss in detail the case of a random access code with qubits, for which we provide analytically optimal self-tests. The simplicity and noise robustness of our methods should make them directly applicable to experiments. arxiv:0.050v [quant-ph] 5 Jan 0 Introduction. Predicting the results of measurements performed on a given physical system has traditionally been the main concern of physics. However, with the advent of deviceindependent quantum information processing [ 3], the opposite question has become relevant. More specifically, given an initially unknown system and an uncharacterized measurement device, what can be inferred about the physics of the experiment based solely on the observed measurement statistics? Despite the apparent generality of this question, certain cases do allow for a precise characterization of the system. This is referred to as self-testing [, 5]. The possibility to self-test quantum states and measurements usually relies on quantum nonlocality. Consider two distant observers performing local measurements on a shared quantum state. When the resulting statistics leads to violation of a Bell inequality [6], it is necessarily the case that the shared quantum state is entangled, and moreover, that the local quantum measurements are incompatible; see e.g. [7]. Furthermore, for specific Bell inequalities, maximal violation (i.e. the largest possible value in quantum theory implies that the quantum state and the measurements can be uniquely identified (up to local isometries. For instance, a maximal violation of the Clauser-Horne-Shimony-Holt (CHSH Bell inequality [] implies maximally incompatible measurements (two anti-commuting Pauli observables and a shared maximally entangled two-qubit state [9 ]. More recently, it has been demonstrated that all bipartite pure entangled states can be self-tested [3], as well as certain multipartite entangled states [ 6]. Another important progress is the development of self-testing methods robust to noise [7 3]. For instance, given a certain level of violation of a Bell inequality (but not necessarily maximal, the fidelity between the initially unknown state and a given target state can be lowerbounded. Self-testing thus offers promising perspectives for the certification of quantum systems in experiments (see e.g. [], as well as for device-independent quantum information protocols. It is therefore natural to ask whether the concept of self-testing can be applied to more general quantum experiments, beyond those based on entanglement and nonlocality. In the present work, we develop self-testing methods tailored to the prepare-and-measure scenario. This covers a broad class of experiments, where quantum communication schemes (e.g. the BB quantum cryptography protocol are prominent examples. In this setting, a preparation device initially prepares a quantum system in different possible states. The system is then transmitted to a measurement device, which performs different possible measurements on it. While it is still possible in this case to characterize certain physical properties of the system based only on statistics, this requires in general an assumption on the devices. One possibility, which we will follow here, is to assume that the set of quantum states and measurements admit a full description in a Hilbert space of given dimension [5 7]. Intuitively this means that the amount of information communicated from the preparation device to the measurement device is assumed to be upper bounded. Such a scenario considering quantum systems of fixed dimension, but otherwise uncharacterized, is referred to as semi-device-independent, and opens interesting possibilities for quantum information processing [ 3]. Here we demonstrate techniques for robustly self-testing sets of prepared quantum states, as well as sets of quantum measurements. These methods allow one to (i assess the compatibility of given sets of preparations and measurements with the observed statistics, and (ii lower-bound the average fidelity between the unknown preparations (measurements and a set of ideal quantum states (measurements. We discuss in detail a simple prepare-and-measure scenario, namely the random access code (RAC. This allows us to provide analytically optimal self-tests for a pair of anti-commuting Pauli observables, and for a set of four qubit states corresponding to the eigenstates of two anti-commuting Pauli observables. We then generalise these results to other prepareand-measure scenarios. The simplicity and robustness of our methods should make them directly applicable to experiments. We conclude with a number of open questions. Scenario. We consider a quantum prepare-and-measure experiment. Upon receiving input x, a preparation device emits a physical system in a quantum state ρ x. The system is then transmitted to a measurement device, which, upon receiving an input y, performs a quantum measurement returning an outcome b. Formally, the measurement is described by a set of positive operators My, b that equal identity when summed over b. Importantly both the specific states ρ x and measurements

2 M b y are a priori unknown to the observer. The statistics of the experiment is then given by P (b x, y = tr(ρ x M b y. In this setting, any possible probability distribution can be obtained, given that the prepared states ρ x can be taken in a sufficiently large Hilbert space. This is however no longer the case when we limit the Hilbert space dimension; specifically we impose that ρ x L(C d for some given d < x (where x denotes the number of possible inputs x. This also corresponds to limiting the amount of information sent from the preparation device to the measurement device to be at most log d bits. In this case, limits on the set of possible distributions can be captured via inequalities of the form A = x,y,b α xyb P (b x, y Q d, ( where α xyb are real coefficients. These are referred to as dimension witnesses [5] whenever (i the inequality is always satisfied when a d-dimensional representation of the experiment is possible, and (ii there exists an experiment involving systems of dimension larger than d such that the inequality is violated. Moreover, it is also meaningful to find the limitations on the set of distributions P (b x, y given that the preparations admit a classical d-dimensional representation. Formally, this means that there exists a d-dimensional basis such that all states ρ x are diagonal in this basis. We denote by C d the maximal value of the quantity A in this case. Interestingly, for wellchosen quantities A, one finds that C d < Q d. Thus, for a given system dimension d (or equivalently for a limit of log d bits of communication, quantum systems outperform classical ones, in the sense that certain quantum distributions cannot be reproduced classically [5]. This quantum advantage can be viewed as the origin for the possibility of developing self-testing methods for the prepare-and-measure scenario; in analogy to Bell inequality violation being the root for selftesting entangled states. In the following we present robust self-testing techniques based on specific dimension witnesses A, both for states and measurements. That is, based only on the value of A, which is directly accessible from the experiment statistics, we will characterize the (initially unknown prepared states and measurements. Specifically, when the maximal value of the witness is obtained, i.e. A = Q d, then a specific set of pure states ρ x = ψ x ψ x and a specific set of projective measurements My b must have been used (up to unitaries. Moreover, when a non-maximal value A < Q d is obtained, one can still efficiently lower bound the fidelity between the created states and measurements and the ideal (or target states and measurement leading to A = Q d. The random access code. We discuss in detail the problem of self-testing states and measurements in a simple prepare-and-measure experiment. This involves four possible preparations, denoted by x = (x 0, x (where x j {0, }, and two possible binary measurements, y {0, } and b {0, }. The score is given by A = x 0,x,y P (b = x y x 0, x, y. ( This means that, upon receiving input y, the measurement device should return the output b = x y, i.e. the y-th bit of the input bit-string x received by the preparation device. Hence the name of a RAC [33 35]. Note that all inputs are assumed to be chosen uniformly at random. Indeed, this task is nontrivial only when d < ; here we will consider the case d =, i.e. qubits. In this case, one finds the tight bounds C = 3/ and Q = ( + / / 0.5 [33]. The classical bound C can be obtained by simply always sending the bit x 0. The quantum bound Q is obtained via the following ideal strategy. The four qubit preparations correspond to the pure states: ρ ideal 00 = + σ x ρ ideal 0 = σ z ρ ideal 0 = + σ z ρ ideal = σ x. (3 These are simply the eigenstates of the Pauli observables σ x and σ z. Next, the measurements are projective and given by two anti-commuting Pauli observables M ideal y = (M 0 y ideal (M y ideal = σ x + ( y σ z ( These qubit preparations and measurements represent the ideal situation, where the maximal value A = Q is achieved. In the following we will determine what restrictions apply to the possible preparations and measurements, given that a particular value of A is observed. In particular, when the maximal value A = Q is attained, both the states and the measurements must be the ideal ones as given above (up to a unitary. Self-testing preparations. Here we find restrictions on the set of prepared states given an observed value of A. For convenience, we write the qubit preparations as ρ x0x = ( + m x0x σ /, where m x0x denotes the Bloch vector (satisfying m x0x and σ = (σ x, σ y, σ z denotes the vector of Pauli matrices. The first step consists in re-expressing A = + y tr ( M 0 y V y + y tr ( M 0 y V y tr ( M 0 y (5 where V y = x 0,x ( xy ρ x0x. In the second step we have used the Cauchy-Schwarz inequality adapted for operators [3]: for any positive semi-definite operator O and any Hermitian operator R, it holds that tr (OR tr ( OR tr (O. Without loss of generality, we can assume that the measurements are projective rank-one operators, since for any set of preparations the optimal measurements are of this form. Consequently, we have that tr ( My 0 =. Next, we obtain Vy = (β + ( y α, where β = m x0x m 00 m m 0 m 0 x 0,x α = ( m 00 m ( m 0 m 0.

3 3 Finally, we find that Eq. (5 reduces to A + [ β + α + β α ]. (6 This provides a tight self-test of the prepared states (in terms of their Bloch vectors, for any given value of A. Let us start with the case A = Q. In order to evaluate the maximum of the right-hand-side of Eq. (6, notice that β + α + β α = β + β α, which is maximal iff α = 0 and β is maximal. This turns out to be achievable. In order to maximise β, we need (i x 0 x : m x0x =, i.e. that all preparations are pure states, and (ii that m 00 m = m 0 m 0 =, i.e. Âĺ the states correspond to (pairwise antipodal Bloch vectors. We define r 0 = m 00 = m and r = m 0 = m 0. Consequently, we find α = r 0 r. Therefore, in order to have α = 0, we must choose r 0 r = 0. This implies that the right-hand-side of Eq. (6 is upper-bounded by Q. Therefore, we conclude that when observing maximal value A = Q, the set of four prepared states must be equivalent (up to a unitary rotation to the set of four ideal states of Eq. (3. More generally, for any value A, one can find a set of preparations (and corresponding measurements such that the inequality (6 is saturated (see Appendix A. For the case of classical preparations (i.e. diagonal in a given basis, the Bloch vectors can simply be replaced by numbers m x0x [, ], and we get A C. Self-testing measurements. Let us now consider selftesting of measurements. Using that M y = My 0 My, we write A = + tr ( ρ x0x 6 [( x0 M 0 + ( x M ] x 0,x (7 + λ max [( x0 M 0 + ( x M ], ( 6 x 0,x where λ max [X] is the largest eigenvalue of the (Hermitian operator X. Since the upper bound corresponds to choosing the optimal preparations for a fixed pair of observables, it simply quantifies the optimal performance achievable using these observables. If M 0 and M are qubit observables the upper bound can be evaluated exactly (see Appendix A to give A + 6 ( µ + ν η + + µ ν η, (9 where µ = tr ( M0 + M, ν = tr{m0, M } and η ± = tr(m 0 ± M. The right-hand side reaches the optimal value Q iff µ =, η ± = 0 and ν = 0, which implies anticommuting projective observables (i.e. projective measurement operators. In other words, observing A = Q implies that the measurements are unitarily equivalent to the ideal ones given in Eq. (. Moreover, note that inequality (9 is tight; for any value of A one can find measurements (and corresponding states such that inequality is saturated (see Appendix A. It follows that any pair of projective, rank-one observables that is incompatible ( ν < can lead to A > C. Robust self-testing of the preparations. We now discuss the problem of characterizing the fidelity between the realized preparations and the ideal ones, as given in Eq. (3. This will allow us to quantify the distance of the prepared states with respect to the ideal ones. Again, we want to develop selftesting methods which are based only on the value of A. More formally, given an arbitrary set of preparations, we define the average fidelity with the ideal preparations to be S({ρ x0x } = max Λ x 0,x F (ρ ideal x 0x, Λ[ρ x0x ]/, where Λ is a quantum channel i.e. a completely positive tracepreserving map. Here the fidelities F (ρ, σ = tr ( ρσ ρ simplify to F (ρ ideal x 0x, Λ[ρ x0x ] = tr ( Λ[ρ x0x ]ρ ideal x 0x, as the ρ ideal x 0x are pure states. We derive lower bounds on the smallest possible value of S given a value of A, i.e., F (A = min S [{ρ x 0x }]. (0 {ρ x0 x } R(A Note that this involves a minimisation over all sets of four preparations R(A that are compatible with an observed value A. In order to lower bound F, we use an approach inspired by Ref. []. From Eq. (7, we have A = + x 0,x tr (W x0x ρ x0x, where W x0x = 6 y ( xy M y. We define operators corresponding to some suitably chosen channel acting on the ideal preparations: K x0x (M 0, M = Λ (M 0, M [ρ ideal x 0x ], ( where Λ is the channel dual to Λ. We aim to construct operator inequalities of the form K x0x (M 0, M sw x0x + t x0x (M 0, M, ( for all inputs (x 0, x, for any given measurements, where s and t x0x (M 0, M are real coefficients. Finding such inequalities, as well as a suitable channel Λ (which could depend on s allows us to lower bound S as follows: S tr (K x0x ρ x0x s tr (W x0x ρ x0x x 0,x x 0,x + t x0x = s (A / + t x0x, (3 x 0,x x 0,x Applying a minimisation over M 0 and M to the right-handside, the above inequality becomes valid for all preparations. Consequently, F(A s (A / + t L (A, ( where t / min M0,M x 0,x t x0x (M 0, M. Next, we construct a channel and a corresponding operator inequality tailored to our problem. Firstly, since for any given set of preparations, the optimal measurements are projective and rank-one, we express the measurements, without loss of generality, as M y = cos(θσ x + ( y sin(θσ z for θ [0, π/]. Secondly, we choose a dephasing channel: Λ θ [ρ] = + c(θ ρ + c(θ Γ(θρΓ(θ. (5

4 bound the following quantity: F (A = min S ( {My} b, (6 {My b} R (A where R (A represents all sets of measurements compatible with a certain value of A. We first rewrite A = y,b tr(m yz b yb, where Z yb = x 0,x ρ x0x δ b,xy. Next, we construct operator inequalities K yb ({ρ x0x } sz yb + t y ({ρ x0x }, (7 FIG.. Average fidelity F (F for prepared states (measurements, as a function of the observed value of A. The black line is our analytical lower bound of Eq. (. The blue region is accessible via single qubit strategies without shared randomness, as confirmed by strong numerical evidence (see Appendix C. When allowing for shared randomness between the devices, the accessible region (obtained by taking the convex hull of the blue region now also includes the grey area, and our analytic lower bound is tight in general. For θ [0, π/] we have c(θ = min{, s sin θ} and Γ = σ x, while for θ (π/, π/] we have c(θ = min{, s cos θ} and Γ = σ z. In Appendix B we construct the operator inequalities (, leading to s = ( + and t = ( /. Inserting these values into Eq. (, we obtain the lower bound. This provides a robust self-testing for the preparations. A maximal value A = Q implies F =, i.e. the preparations must be the ideal ones (i.e. the states of Eq. (3, up to a unitary. For A = C, i.e. a maximal value given a set of classical states, we get that F 3/. This bound can be attained via the set of pure states ρ x0x = ( + ( x0x σ z / (diagonal in the same basis, hence classical, combined with the measurements M 0 = M = σ z. Therefore, we see that our bound F (A L(A is optimal, as far as linear inequalities are concerned (see Fig.. It is then interesting to consider the intermediate region C < A < Q. First, focusing on strategies involving a single set of states and measurements, we observe numerically that the linear bound F (A L(A cannot be saturated anymore, and conjecture the form of optimal states and measurements; see red curve in Fig. and Appendix C for details. Second, allowing for shared randomness between the preparation and measurement device (such that convex combinations of qubit strategies are now possible, the linear bound becomes tight, a direct consequence of the linearity of F and A in terms of the states and measurements. Robust self-testing of the measurements. Similarly, we can quantify the average fidelity of the measurements with respect to the ideal ones (given in Eq. (: S ( {M b y} = max Λ y,b F ( (M b y ideal, Λ[M b y] /, where Λ must be a unital channel (i.e. mapping the identity to itself, in order to ensure that measurements are mapped to measurements. In analogy with the case of preparations, our goal is to lower given the unital channel K yb = Λ[(My b ideal ]. Similarly to the case of preparations, strong operator inequalities can be derived by choosing carefully the channel; all details are given in Appendix D. Finally, this leads to a lower bound on the average fidelity F (A min {ρ x0 x } tr ( K yb My b L(A. ( y,b That is, we find that F can be lower-bounded by a linear expression in terms of A, which turns out to be the same as for the case of preparations. This provides a robust self-test for the measurements. Observing the maximal value of A = Q implies that F =, hence the measurements are equivalent to the ideal ones (i.e. two anti-commuting Pauli observables as given in Eq. (, up to a rotation. For a value of A = C, we have that F 3/. This lower bound can be attained by choosing M 0 = σ z and M =, with the states ρ 00 = ρ 0 = ( + σ z / and ρ 0 = ρ = ( σ z /. For the intermediate regime C < A < Q, we find numerically that the inequality ( cannot be saturated using a single set of measurements and states (see Fig.. Details, in particular a conjecture for the form of the optimal measurements, are given in Appendix C. Similarly as for the case of states, when allowing for convex combinations of qubit strategies (using shared randomness our linear bound is tight. Generalisations. The above results can be generalised in several directions, by considering more general prepare-andmeasure scenarios. First, a generalisation of the RAC provides a self-testing condition for any pair of incompatible Pauli observables, as shown in Appendix E. Next, we consider the N RAC, where the preparation device receives as input an N-bit string x = (x,..., x N and the measurement device gets input y {,..., N}. The average score is then given by A N = N N P (b = x y x, y. (9 x,y The methods discussed above (for N = can be generalised and lead to self-testing relations for preparations and measurements; details are given in Appendix F. Thus we obtain necessary conditions for a set of preparations (and measurements to be compatible with a given observed value of A N. The case of N = 3 is of particular interest. Here, the best possible score with qubits is A 3 = ( + / 3/; see e.g. [35].

5 5 In this case, our self-testing conditions can certify that (i the eight prepared states correspond to Bloch vectors forming a cube on the Bloch sphere, and (ii the measurements correspond to three mutually unbiased bases (i.e. three pairwise anti-commuting Pauli observables. Finally, it is worth mentioning that our methods can also be extended beyond the case of qubits. In Appendix G, we discuss the case of the RAC with qutrit preparations, i.e. ρ x0x L(C 3, and derive a self-testing condition for projective measurements. Outlook. We presented methods for self-testing quantum states and measurements in the prepare-and-measure scenario. These techniques demonstrate strong robustness to noise, and should therefore be directly amenable to experiments, where they could provide useful certification techniques in a semidevice-independent setting. It would be interesting to develop robust self-testing techniques for more general scenarios, and also consider higher dimensional quantum systems. Another direction would be to consider scenarios beyond prepare-and-measure, for instance adding between the preparation and measurement devices a transformation device [36, 37], and self-test the latter. Finally, while we have focused here on the prepare-andmeasure scenario with an assumption on the dimension of the system, it would be interesting to develop self-testing methods based on different assumptions, such as an upper-bound on the energy [3] or a lower-bound on the overlap [39]. Acknowledgements. This work was supported by the Swiss national science foundation (Starting grant DIAQ, QSIT, and Early Postdoc Mobility fellowship PGEP_6060, the European Union s Horizon 00 research and innovation programme under the Marie Skłodowska-Curie Action ROSETTA (grant no. 7936, the European Research Council (grant no , the Danish Council for Independent Research (Sapere Aude, VILLUM FONDEN via the QMATH Centre of Excellence (grant no. 0059, and the National Research, Development and Innovation Office NKFIH (Grants No. K73, and No. KH5096. [] A. Acín, N. Brunner, N. Gisin, S. Massar, S. Pironio, and V. Scarani. Phys. Rev. Lett., 9, 3 (007. [] R. Colbeck, PhD thesis, University of Cambridge, 007. ArXiv: [3] S. Pironio, A. Acín, S. Massar, A. Boyer De La Giroday, N. D. Matsukevich, P. Maunz, S. Olmschenk, D. Hayes, L. Luo, T. A. Manning, and C. Monroe, Nature 6 0 (00. [] D. Mayers, and A. Yao, Proc. 39th FOCS (Washington, DC: IEEE Computer Society p 503 (99. [5] D. Mayers, and A. Yao, Quantum Information and Computation, (:73-6, July 00 [6] J. S. Bell, Physics, 95 (96. [7] N. Brunner, D. Cavalcanti, S. Pironio, V. Scarani, and S. Wehner, Rev. Mod. Phys. 6, 9 (0. [] J. F. Clauser, M. A. Horne, A. Shimony, and R. A. Holt, Phys. Rev. Lett., 3 (5, 969. [9] S. J. Summers, and R. Werner J. Math. Phys. 0 (97. [0] S. Popescu and D. Rohrlich, Phys. Lett. A, 69(6, 99. [] B. S. Tsirelson, Hadronic J. Suppl.,, 39 (993. [] B. W. Reichardt, F. Unger and U. Vazirani, Nature 96, 56 (03. [3] A. Coladangelo, K. T. Goh, and V. Scarani, Nature Communications, 55 (07. [] M. McKague, Proc. TQC, LNCS, 675, 0. [5] K. F. Pál, T. Vértesi, and M. Navascués, Phys. Rev. A, 90, 030 (0. [6] X. Wu, Y. Cai, T. H. Yang, H. N. Le, J.-D. Bancal, and V. Scarani. Phys. Rev. A, 90, 0339 (0. [7] C.-E. Bardyn, T. C. H. Liew, S. Massar, M. McKague, and V. Scarani. Phys. Rev. A 0, 0637 (009. [] M. McKague, T. H. Yang, and V. Scarani, J. Phys. A: Math. Theor., (0. [9] T. H. Yang, and M. Navascués, Phys. Rev. A, 7, 0500(R (03. [0] C. Bamps and S. Pironio. Phys. Rev. A, 9, 05 (05. [] T. H. Yang, T. Vértesi, J-D. Bancal, V. Scarani, and M. Navascués, Phys. Rev. Lett. 3, 000 (0. [] J. Kaniewski, Phys. Rev. Lett., 7, 0700 (06. [3] J. Kaniewski, Phys. Rev. A 95, 0633 (07. [] T. R. Tan, Y. Wan, S. Erickson, P. Bierhorst, D. Kienzler, S. Glancy, E. Knill, D. Leibfried, and D. J. Wineland, Phys. Rev. Lett., 3003 (07. [5] R. Gallego, N. Brunner, C. Hadley, and A. Acin, Phys.Rev.Lett. 05, 3050, (00. [6] M. Hendrych, R. Gallego, M. Micuda, N. Brunner, A. Acin, and J. P. Torres, Nature Physics, 5 (0. [7] J. Ahrens, P. Badziag, A. Cabello, and M. Bourennane, Nature Physics, 59 (0. [] M. Pawłowski and N. Brunner, Phys. Rev. A, 0030(R (0. [9] H-W Li, M. Pawłowski, Z-Q. Yin, G-C. Guo, and Z-F. Han, Phys. Rev. A, 0330 (0. [30] E. Woodhead, and S. Pironio, Phys. Rev. Lett. 5, 5050 (05. [3] T. Lunghi, J. B. Brask, C. C. W. Lim, Q. Lavigne, J. Bowles, A. Martin, H. Zbinden, and N. Brunner, Phys. Rev. Lett., 5050 (05. [3] P. Mironowicz, A. Tavakoli, A. Hameedi, B. Marques, M. Pawłowski, and M. Bourennane, New J. Phys., (06. [33] A. Ambainis, A. Nayak, A. Ta-Shama, U. Vazirani, in Proceedings of 3st ACM Symposium on Theory of Computing, pp , 999. [3] A. Nayak, in Proceedings of the 0th IEEE Symposium on Foundations of Computer Science (FOCS 99, pp , 999. [35] A. Tavakoli, A. Hameedi, B. Marques, and M. Bourennane, Phys. Rev. Lett., 7050 (05. [36] P. Trojek, C. Schmid, M. Bourennane, C. Brukner, M. Żukowski, and H. Weinfurter, Phys. Rev. A 7, (R (005. [37] J. Bowles, N. Brunner, M. Pawłowski, Phys. Rev. A 9, 035 (05. [3] T. Van Himbeeck, E. Woodhead, N. J. Cerf, R. García-Patron, S. Pironio, Quantum, 33 (07. [39] J. B. Brask, A. Martin, W. Esposito, R. Houlmann, J. Bowles,

6 H. Zbinden, and N. Brunner, Phys. Rev. Applied 7, 050 (07. 6

7 7 Appendix A: Self-testing relations for preparations and measurements In this appendix we provide a simple example of preparations that saturate the bound given in Eq. (6. Moreover, we derive the upper bound given in Eq. (9 and provide a family of measurements that saturates the bound. First, we consider the case of preparations. Consider preparations such that ρ 00 and ρ, and ρ 0 and ρ 0 correspond to antipodal Bloch vectors with a relative angle θ, the maximal quantum value of A, is obtained from the left-hand-side of Eq. (5 by A = + λ max [V y ], where V y = x 0,x ( xy ρ x0x. We represent the preparations on the Bloch sphere as ρ x0x = / ( + m x0x σ, where m 00 = [cos(θ/, 0, sin(θ/] and m 0 = [cos(θ/, 0, sin(θ/], with m = m 00 and m 0 = m 0. This gives V 0 = cos(θ/σ x and V = sin(θ/σ z. The respective largest eigenvalues are λ max [V 0 ] = cos(θ/ and λ max [V ] = sin(θ/, leading to A = + [ ] + cos θ + cos θ. (A It is straightforward to see that this achieves the upper bound of Eq. (6; indeed the above choice of preparations leads to β = and α = cos θ. To derive the upper bound given in Eq. (9 we evaluate [ ] λ max ( x0 M 0 + ( x M (A3 x 0,x for arbitrary qubit observables M 0, M. We take advantage of the fact that y (A λ max [T ] + λ max [ T ] = λ max [T ] λ min [T ], (A which for a matrix can be evaluated analytically. More specifically, if T is a Hermitian matrix with eigenvalues λ 0 λ, let and then χ := tr T = λ 0 + λ, ζ := tr T = λ 0 + λ λ 0 λ = ζ χ. (A5 Evaluating this expression for T = M 0 ± M gives the upper bound given in Eq. (9. In order to show that Eq. (9 can be saturated, we consider a specific family of measurements. Choose the measurements as projective rank-one with associated observables described by Bloch vectors n 0 = [cos (θ/, 0, sin (θ/] and n = [cos (θ/, 0, sin (θ/]. Evaluating Eq. (9, we find that η ± = 0, µ = and ν = cos θ, leading to the upper bound A + ( + cos θ + cos θ. (A6 In order to saturate this bound, we calculate the maximal value of A, as given in Eq. (, for such measurements: Thus, we have saturated the bound of Eq. (9. A = + (cos (θ/ + sin (θ/ = + ( + cos θ + cos θ. (A7 Appendix B: Operator inequalities for robust self-testing of preparations In this section we provide a detailed derivation of the lower bound on the average fidelity F(A. As explained in the main text, our approach relies on deriving operator inequalities of the form K x0x (M 0, M sw x0x + t x0x (M 0, M, (B

8 where W x0x = 6 y ( xy M y and K x0x (M 0, M = Λ (M 0, M [ρ ideal x 0x ], for some suitably chosen channel Λ. We choose a dephasing channel of the form Λ θ (ρ = + c(θ ρ + c(θ Γ(θρΓ(θ, where for 0 θ π/ we use Γ = σ x, while for π/ < θ π/ we use Γ = σ z. The function c(θ [, ] will be specified later. In the interval 0 θ π/, the action of the channel leads to (B K 00 = + σ x K 0 = + c(θσ z whereas in the interval π/ < θ π/, we have K 0 = c(θσ z K = σ x, (B3 K 00 = + c(θσ x K 0 = + σ z K 0 = σ z K = c(θσ x. (B As discussed in the main text, for any given set of preparations, the optimal measurements are projective and rank-one. Furthermore, any two such measurements can be represented on an equator of the Bloch sphere. Due to the freedom of setting the reference frame, we can without loss of generality represent the two measurements in the xz-plane, i.e. We can therefore write W x0x as M 0 = cos θ σ x + sin θ σ z M = cos θ σ x sin θ σ z. W 00 = cos θσ x W 0 = sin θσ z W 0 = sin θσ z W = cos θσ x. (B5 We can reduce the number of operator inequalities (B by exploiting the apparent symmetries in the expressions for W x0x and K x0x : we restrict ourselves so that t o t 0 = t 0 and t e t 00 = t. Thus, we have to consider two operator inequalities in each interval θ [0, π/] and θ (π/, π/]. In the first interval, the two operator inequalities are + σ x s cos θσ x t e 0 + c(θσ z s sin θσ z t o 0. (B6 In the second interval, the two operator inequalities are + c(θσ x s cos θσ x t e 0 + σ z We now focus on the former interval. Solving the two inequalities for t o and t e we obtain s sin θσ z t o 0. (B7 t e s cos θ, t o ( + c(θ s sin θ, (B t e s cos θ, t o ( c(θ + s sin θ. (B9 Any choice of t o and t e satisfying these constraints gives rise to valid operator inequalities. In order to obtain the strongest bound, we choose the largest values of t o and t e consistent with their respective constraints, i.e., t e = min { s cos θ, s } { cos θ t o = min ( + c(θ s sin θ, } ( c(θ + s sin θ. (B0 A similar procedure for the interval θ (π/, π/] leads to { t e = min ( + c(θ s cos θ, } ( c(θ + s cos θ t o = min { s sin θ, s } sin θ. (B It is worth pointing out that the two intervals only differ by exchanging t e t o and sin θ cos θ. Hence, for any given θ, we have constructed operator inequalities of the form (B. From Eq. ( of the main text, we obtain our lower bound on the average fidelity F (A s (A / + min M 0,M t(m 0, M L(A, (B

9 9 where t(m 0, M = (t e + t o /. To compute this quantity we fix the value of s to be s = ( + (B3 and choose the dephasing function as c(θ = min{, s sin θ} whenever θ [0, π/], and c(θ = min{, s cos θ} whenever θ (π/, π/]. It is easy to see that c(θ [0, ], which ensures that Λ θ is a valid quantum channel, and that c(θ is continuous at θ = π/. A simple calculation shows that in this case which gives the lower bound F(A One can check that choosing distinct values of s will not lead to improved lower bounds. t =, (B ( + A 3 L(A. (B5 Appendix C: Tightness of fidelity bounds In the main text, we have derived fidelity bounds for both the preparations and the measurements, based on operator inequalities. Specifically, we obtain a lower bound on the average fidelity F of the prepared states (with respect to the ideal ones given by the linear expression F(A ( + A 3 L(A. (C For measurements, a similar bound is obtained on the average fidelity F with respect to the ideal ones. In the present appendix, we discuss the tightness of these bounds. We start with our bound on the fidelity of the states. As discussed in the main text, obtaining A = Q implies F =, i.e. the states are the ideal ones (up to a unitary. Let us refer to the optimal strategy (with the ideal states as strategy S. Then, for A = C, our bound gives F 3/. This bound is tight and can be obtained via the set of pure states ρ x0x = ( + ( x0x σ z / (diagonal in the same basis, hence classical, combined with the measurements M 0 = M = σ z. Let us refer to this strategy as S. The above shows that our bound (C is tight as far as linear inequalities are concerned. More generally, the bound is in fact tight in general, when shared randomness between the preparation and measurement devices is taken into account. In this case, taking a convex combination between strategies S and S allows us to get any point on the line (i.e. pair of values F and A between S and S. It is also interesting to understand what happens when shared randomness between the devices is not taken into account. In this case, the end points (A = Q, F = and (A = C, F = 3/ can still be obtained. To understand what happens in the intermediate region C < A < Q, we first performed a numerical analysis. Specifically, we chose randomly four qubit states, and compute (i the maximal value of A (optimizing over the measurements, and (ii the average fidelity F (where the optimization over channels is restricted here to unitaries. The resulting points are shown on Fig. (blue circles. This indicates that for C < A < Q, the bound (C cannot be saturated anymore. Moreover, we conjecture that an optimal class of strategies is given by the pure states ψ 00 = 0 ψ = ψ 0 = cos θ 0 + sin θ ψ 0 = cos θ 0 sin θ (C and the measurements M y = cos(ϕσ z + ( y sin(ϕσ x. Straightforward calculations show that, taking tan ϕ = sin θ, leads to A = + + tan (ϕ and F = (3 + tan ϕ. (C3 This gives a parametric curve, as a function of ϕ [0, π/], given by the red curve in Fig.. This curve is in excellent agreement with the numerical results obtained before. Note that this class of strategies interpolates between the strategies S (setting ϕ = 0 and S (setting ϕ = π/. Next we discuss the bound on the average fidelity of measurements. As discussed in the main text, the linear bound F (A L(A is optimal as far as linear inequalities are concerned. Moreover, when allowing for shared randomness the bound is tight in general for C A Q. This is obtained by considering convex combinations of strategy S (defined as the optimal

10 0 FIG.. The black line is the analytic lower bound on the average fidelity F (F for prepared states (measurements, as a function of the observed value of A. To characterize the region accessible via pure qubit strategies (i.e. without shared randomness, we perform numerics generating randomly sets of qubit preparations (blue circles and crosses; here we show the numerical results for the case of states, but similar results are obtained for the case of measurements. In the region C < A < Q, we conjecture that the class of strategies given in the text (corresponding to the red curve are optimal; both for F and F. Finally, the green dashed line is our conjectured upper bound on the average fidelity. strategy S, up to a rotation of π/ around the y axis; see below, and the following strategy (referred to as S 3 : take M 0 = σ z and M =, with the states ρ 00 = ρ 0 = ( + σ z / and ρ 0 = ρ = ( σ z /. Similarly to the case of states, we now consider the situation where shared randomness between the devices is not allowed. Performing a numerical analysis similar to the one described above (except that measurements are no generated randomly, we observe that the accessible region (in terms of F vs A appears to be exactly the same as for the case of states (i.e. the blue region in Fig.. We conjecture that the lower bound is given by the following class of optimal strategies: take the measurements: M 0 = σ z and M = ησ x + ( η (C with the states ψ 00 = cos θ 0 + sin θ, ψ 0 = cos θ 0 sin θ, ψ 0 = cos θ + sin θ 0, and ψ = cos θ sin θ 0. Setting η = tan θ, we get A = cos (θ + + sin (θ cos(θ and F = (3 + tan(θ. (C5 This gives a parametric curve, as a function of θ [0, π/], given by the red curve in Fig.. This curve is in excellent agreement with the numerical results obtained before. Also, this curve turns out to be exactly the same as the curve we obtained above for the case of states. Note that this class of strategies interpolates between the strategies S (setting θ = π/ and S 3 (setting θ = 0. Finally, note that the numerics also suggests that there is a linear upper bound on the average fidelities F (F as a function of A (see Fig. ; specifically F Q Q A 3/ + Q 3/ Q and similarly for F 3/. It would be interesting to provide a proof of these upper bounds. Appendix D: Operator inequalities for robust self-testing of measurements In this section, we account for the detailed derivation of the lower bound on the average fidelity of the measurements F (A. The approach bears significant resemblance to the case of robustly self-testing preparations, as outlined in Appendix B. We aim to derive operator inequalities of the form K yb ({ρ x0x } sz yb + t y ({ρ x0x }, (D where Z yb = x 0,x ρ x0x δ b,xy and K yb ({ρ x0x } = Λ [(My b ideal ]. For the sake of simplicity, we first apply a unitary channel to (My b ideal align these operators with the eigenstates of σ x and σ z. Then, we adopt the same (unital channel Λ as

11 specified in Eq. (5 with the same coefficients as used to robustly self-test the preparations: c(θ = min{, s sin θ} when θ [0, π/] and c(θ = min{, s cos θ} when θ (π/, π/]. It is straightforward to see that, for any given pair of measurements, the optimal choice of preparations are four pure qubit states, such that ρ 00 and ρ, and ρ 0 and ρ 0 respectively, correspond to antipodal vectors on the Bloch sphere. Therefore, we can without loss of generality restrict to such preparations since these impose the weakest constraints on the measurements of our interest. We can therefore parametrise the preparations ρ x0x = / ( + m x0x σ by Bloch vectors m 00 = [cos θ, 0, sin θ] m = [cos θ, 0, sin θ] m 0 = [cos θ, 0, sin θ] m 0 = [ cos θ, 0, sin θ]. (D Expressing Z yb in terms of these preparations gives Z 00 = ( + cos θσ x Z 0 = ( cos θσ x Z 0 = ( + sin θσ z σ z Z = ( sin θσ z. (D3 Due to symmetries, we restrict ourselves so that t o t 0 = t 0 and t e t 00 = t. Thus, we have to consider two operator inequalities in each interval θ [0, π/] and θ (π/, π/]. In the first interval, the two operator inequalities are + σ x s ( + cos θσ x t e 0 + c(θσ z s ( + sin θσ z t o 0. (D In the second interval, the two operator inequalities are + c(θσ x s ( + cos θσ x t e 0 + σ z s ( + sin θσ z t o 0. Just as in Appendix B, we solve these inequalities for t e and t o, and choose the largest value compatible with the solutions. In the first interval, this gives { t e = min ( s s cos θ, s } (cos θ (D5 { t o = min (c(θ s sin θ s +, } ( c(θ + s sin θ s +. (D6 A similar procedure for the interval θ (π/, π/] leads to { t e = min (c(θ s cos θ s +, } ( c(θ + s cos θ s + { s t o = min (sin θ, } ( s s sin θ. (D7 For any choice of θ, we have constructed operator inequalities of the form (D. In order to obtain our lower bound on F, we must minimise the quantity t(θ = (t e + t o / for a specific choice of s. In analogy with the procedure in Appendix D, we choose s = ( +, which returns min θ t(θ = 3/(. Hence, we have obtained the lower bound ( F (A + A 3 = L(A. (D Appendix E: Self-testing all pairs of incompatible Pauli observables Consider a generalisation of the RAC, in which we introduce a bias on the score associated to certain inputs. Specifically, whenever the game is successful, i.e. b = x y, the awarded score is q/ if x 0 x = 0, and ( q/ if x 0 x =, for some q [0, ]. The average score reads A q = x 0,x,y r(x 0, x P (b = x y x 0, x, y, where r(x 0, x = q/ if x 0 x = 0 and r(x 0, x = ( q/ if x 0 x =. Note that for q = /, we recover the standard RAC. Based on the quantity A q, we will now see how to derive a self-testing condition for any pair of incompatible Pauli observables, i.e. any pair of non-commuting projective rank-one qubit measurements. (E

12 We start by expressing A q for a quantum strategy: A q = + r(x 0, x tr (ρ x0x (( x0 M 0 + ( x M + r(x 0, x λ max [( x0 M 0 + ( x M ]. x 0,x x 0,x (E [ ] [ ] Denoting µ k = λ min M0 + ( k M and νk = λ max M0 + ( k M, for k = 0,, we obtain A q + [q (µ 0 ν 0 + ( q (µ ν ]. (E3 Following a derivation analogous to that appearing in Appendix A to obtain Eq. (9 in the main text, we obtain A q + [ q β + α + ( q ] β α, (E where β = tr ( M 0 + M tr (M0 tr (M and α = tr ({M 0, M } tr (M 0 tr (M. Treating α and β as independent variables, we obtain the largest value of the right-hand-side of Eq. (E by demanding that the derivative with respect to α equals zero, and checking that the second derivative is negative at this point. We obtain the optimality constraint α = q β. (E5 q + q Inserting this value back into Eq. (E, we find an upper bound on A q as obtained by independent variables α and β. It turns out that this bound can be saturated by the de facto coupled variables α and β. From Eq. (E, it is clear that a necessary condition for optimality is to maximise β. This amounts to the observables M 0 and M being traceless and such that M0 = M =, leading to β =. This implies that the observables represent projective rank-one measurements. Hence, we can write M y = n y σ where the Bloch vector satisfies n y =. Hence, we have α = n 0 n. Thus, Eq.(E5 becomes n 0 n = q q + q, (E6 which has a solution for any choice of q. Note that setting q = / reduces the above to n 0 n = 0 which we recognize as the optimality constraint for the standard random access code. In conclusion, for any pair of incompatible Pauli observables (characterized by the scalar product n 0 n, we have a game A q (where q is chosen in order to satisfy the above equation, such that the maximal score can only be attained by using that specific pair of Pauli observables. We thus obtain a general class of self-tests for any pair of Pauli observables, corresponding to saturating the maximal quantum value of A q for a given value of q: A q ( + q + q. (E7 Appendix F: Self testing for the N random access code In this appendix, we extend the results presented in the main text to self-test the preparations and measurements in an N RAC. The latter is a straightforward generalisation of the RAC considered in the main text. The input of the preparation device is a random N-bit string x (x,..., x N, while the input of the measuremet device is y {,..., N}. The average score is A N = N N P (b = x y x, y. Considering qubit states ρ x, and measurement observables M y, we get x,y (F A N = + N N+ ( xy tr (ρ x M y. x,y (F

13 3. Compatibility of measurements We determine whether a set of measurements can explain (i.e. are compatible with a given value of A N. Since rank-one projective measurements are optimal for any set of preparations, we choose for simplicity to restrict our consideration to such measurements. However, it is straightforward to consider general measurements using the method outlined in the main text and Appendix A. Specifically, we first write A N = + N N+ tr (ρ x W x + N N+ λ max [W x ], x where W x = y ( xy M y. Note λ max [W x ] = λ min [W x ], where x = ( x,..., x N the bit string obtained from x by flipping all bits. Thus it is sufficient to only calculate eigenvalues for the strings not obtainable from each other under a full bit-flip operation. To this end let z = x... x N, 0 and λ z,0 (λ z, be the largest (smallest eigenvalue of W z. Thus we write x (F3 A N + N N+ [λ z,0 λ z, ]. (F Since λ z,0 and λ z, are eigenvalues of W z, we have λ z,0 + λ z, = tr ( W z, which is equivalent to z N λ z,0 + λ z, = tr ( My + ( z k+z l tr ({M k, M l }. y= k<l (F5 This equation, together with the relation (λ z,0 λ z, ( λ z,0 + λ z,, imply that Eq. (F becomes A N [ + N N N+ tr ( ] / My + ( z k+z l tr ({M k, M l }. (F6 z y= This provides a robust self-testing condition, allowing one to determine whether a given set of measurements is compatible with the observed value of A N. Furthermore, we can derive an upper bound on the maximal value of A N by assuming (incorrectly for N > 3 that there exists N mutually unbiased bases in C. This means that all measurements are maximally incompatible, i.e. that tr ({M k, M l } = 0 for k l. Consequently, Eq. (F6 reduces to A N k<l ( +. (F7 N We emphasize that only three mutually unbiased bases exist in C and hence this bound is only tight for N =, 3. For N =, we recover the result presented in the main text. For N = 3, this implies that a maximal value of A 3 (i.e. achieving the righthand side of the above inequality ensures that the measurements are three mutually unbiased qubit observables, such as the three Pauli observables σ x, σ y and σ z. Going one step further, we can then also self-test the preparations (still assuming maximal value of A 3. Indeed, each preparation ρ x must be pure, and correspond to the eigenvector of W x associated to its largest eigenvalue. Such a set of preparations correspond to a set of Bloch vectors forming a cube on the surface of the Bloch sphere.. Compatibility of preparations We ask whether a given value of A N can be explained by a particular set of preparations. We suitably express (F in a quantum model and subsequently apply the Cauchy-Schwarz inequality for operators to obtain A N = + N N [ N tr y= M 0 y ] ( xy ρ x x + N N N tr y= M 0 y ( ( xy ρ x. (F x

14 In the last expression, the squared operator is evaluated to ( ( xy ρ x = x x ρ x + k<l( ky+ly {ρ k, ρ l }. (F9 If necessary, the anticommutators can be evaluated using Bloch sphere representation with the relation {ρ k, ρ l } = / (( + m k m l + ( m k + m l σ. However, it is more convenient to consider a basis-independent representation. Importantly, note that since an equal number of positive and negative terms appear inside the square, the operator x ( xy ρ x is a linear combination of {σ x, σ y, σ z } and hence its square is proportional to the identity operator. Therefore, when re-inserting Eq. (F9 into Eq. (F, we find A N + N N [ N tr ( / ρ x + ( ky+ly tr ({ρ k, ρ l }]. (F0 y= x This is a self-testing condition for preparations, assessing whether a given set of preparations is compatible with a given value of A N. In particular, a classical strategy in which the preparations are binary messages corresponds to x : tr ( ρ x = and tr ({ρ k, ρ l } = δ E(k,E(l, where E is the specific classical encoding strategy, i.e. a function E : {0, } N {0, }. k<l Appendix G: Self-testing with three-level systems In the main text, we have considered self-testing in the random access code when the physical system transmitted from Alice to Bob is a qubit. Clearly, if that system is allowed to carry two bits of information, the task is trivial since Alice can send both her inputs to Bob. Here, we consider the remaining non-trivial case of Alice communicating a three-level quantum system. To simplify the analysis we restrict ourselves to projective measurements for which all possible arrangements admit a compact characterisation. We show that the optimal quantum value equals A = ( / and find all the optimal arrangements of observables (we argue that the optimal value is achieved only if both measurements are projective. Our argument is robust in the sense that we are able to certify incompatibility of M 0 and M whenever the success probability exceeds the classical bound for three-level systems, which turns out to be A 7/. To obtain a statement which only depends on the observables we employ the bound given in Eq. (, i.e. we are interested in evaluating the sum [ ] λ max ( x0 M 0 + ( x M. (G x 0,x Jordan s lemma states that any two projective observables can be simultaneously diagonalised such that the resulting blocks are or. For observables acting on a qutrit, we only need to consider two cases: (a three -dimensional subspaces or (b one subspace of each type. Case (a corresponds to classical strategies and it is easy to check that these satisfy A 7/. In case (b the observables (up to a unitary can be written as ( ( cos α σx + sin α σ M 0 = z cos α σx sin α σ and M r = z (G s for some angle α [0, π] and r, s {±}. A simple calculation yields λ max [M 0 + M ] = max{ cos α, r + s}, λ max [M 0 M ] = max{ sin α, r s}, λ max [ M 0 + M ] = max{ sin α, r + s}, λ max [ M 0 M ] = max{ cos α, r s} and, therefore, x 0,x λ max [ ( x0 M 0 + ( x M ] = { + sin α + cos α if r = s, + sin α + cos α if r s. (G3 For r = s the right-hand side is maximised for α {c, c +π, c +π, c +π}, where c is the unique solution to tan c = in the interval [0, π/]. Similarly, for r s the right-hand side is maximised for α {c, c + π, c + π, c + π}, where c is the unique solution to tan c = / in the interval [0, π/].

15 While the different optimal arrangements are not unitarily equivalent, they are of similar form. The optimal arrangement characterised by r = s = and α = c yields the following optimal preparations ( ( ( ( 0 ( + σz / ( σz / ( σx / ρ 00 =, ρ 0 =, ρ 0 0 =, ρ 0 =. (G 0 Indeed, it is always the case that one preparation lives in the subspace, whereas the other three occupy the subspace (two of them form a basis to which the last one is unbiased. To see that the optimal winning probability requires projective measurements note that for every set of preparations the optimal observables can be chosen projective. However, all sets of preparations optimal for projective observables are of the form given above and one can check that for these preparations the optimal measurements must be projective (a direct consequence of the fact that the operators ρ 00 + ρ 0 ρ 0 ρ and ρ 00 ρ 0 + ρ 0 ρ are full-rank. It is the presence of multiple inequivalent maximisers that prevents us from writing down a simple self-testing statement. However, Eq. (G3 allows us to deduce the range of α compatible with the observed value of A (note that the conclusion will be stronger if we know whether r = s or r s. In particular, any value exceeding the classical bound of 7/ implies a lower bound on the incompatibility between M 0 and M on the subspace. 5