This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Asymptotic implementation of multipartite quantum channels and other quantum instruments using local operations and classical communication

Scott M. Cohen [email protected] Department of Physics, Portland State University, Portland Oregon, USA 97201
Abstract

We prove a necessary condition that a quantum channel on a multipartite system may be approximated arbitrarily closely using local operations and classical communication (LOCC). We then extend those arguments to obtain a condition that applies to all quantum instruments, which range from the most refined case, a generalized measurement, to the most coarse-grained, which is a quantum channel. We illustrate these results by a detailed analysis of a quantum instrument that is known not to be implementable by LOCC, but which can be arbitrarily closely approximated within that framework. As one outgrowth of this analysis, we find a quantum measurement that falls into the same category: it cannot be implemented exactly by LOCC, but can be approximated by LOCC arbitrarily closely. This measurement has an infinite number of outcomes, leaving open the question as to whether or not there exists a measurement within this same category but having only a finite number of outcomes.

pacs:
03.65.Ta, 03.67.Ac

I Introduction

In two recent papers [18, 20], we developed a new approach to the study of quantum measurements implemented by multiple parties utilizing local operations in their respective laboratories and classical communication to their partners to share the results of those local operations, a fundamental paradigm in quantum information processing commonly known as LOCC. We believe our approach provides valuable insights and tools for analyzing what can be accomplished under these circumstances, where quantum systems are not (under practical conditions, often cannot be) transmitted from one party to another. In the first of these papers [18], we identified certain geometric conditions that must hold in order that a measurement can be approximated arbitrarily closely by LOCC, and then in the second paper [20], we refined that approach for use in studying local quantum state discrimination, showing how one can often determine whether success with such a task is impossible, even while knowing only the set of states to be discriminated but nothing about which measurements, if any, might be possible candidates for successfully accomplishing such a task. Here, we generalize our approach to address the question of whether or not a given quantum channel can be closely approximated using LOCC. We then extend the ideas further to apply to any quantum instrument [11].

The paradigm of LOCC has been widely studied [4, 1, 3, 2, 36, 13, 32, 34, 8, 5, 26, 24, 7, 16, 17, 5] since publication of the seminal paper of Peres and Wootters [33]. While much progress has been made, and many important and sometimes surprising results have been obtained, LOCC has also been widely recognized as being quite difficult to characterize. One particularly challenging aspect arises when one is interested in whether or not a task can be accomplished while allowing an arbitrarily small error to be incurred in doing so. In studying such questions, one must allow the parties to share their outcomes back and forth for as many rounds of communication as needed in order to reduce the error as much as possible. With some tasks that can be easily accomplished by acting globally on the entire multipartite system, it is known that when restricted to LOCC, even allowing for an infinite number of rounds is not enough to reduce the error to zero. One such celebrated example involves perfect discrimination of the so-called domino states [3], which demonstrated that separable measurements—ones that involve measurement operators that are simple tensor products on a multipartite system and can be easily implemented globally—cannot always be implemented, or even closely approximated, by LOCC, even with the use of an infinite number of communication rounds. This discovery gave birth to a wide-ranging effort at understanding the surprisingly difficult, and rich, questions concerning the difference between global actions and those restricted to LOCC.

While LOCC has a rich and difficult to characterize structure, one commonly discusses it in terms of a relatively simple picture, that of a rooted tree graph, or LOCC tree. In the case of a finite LOCC tree and starting at the root node, one may follow one of the tree’s many branches to its endpoint, which is known as a leaf node. Each node along any given branch represents an outcome of a (local) measurement performed by one of the parties, and the next node along the branch represents an outcome obtained in a subsequent (local) measurement by a different party. These nodes are connected by edges, and it is easy to envision a path following this branch along these successive edges from the root node all the way to its leaf. We recognized [18] that this path is not just an abstract visualization of the sequence of outcomes, otherwise having no useful meaning, but rather that it corresponds directly to a particular type of path lying within a specific geometric object. Significantly, that object is determined uniquely by the outcomes associated with the leaf nodes appearing at the ends of the collection of all branches in the tree. That is to say the measurement {\cal M} implemented by the LOCC protocol represented by this tree uniquely identifies a geometric object ZZ_{\cal M}, within which lies each and every path associated with the branches in the tree. We have characterized these paths [18] and used that characterization to obtain a powerful condition, reproduced below as Theorem 1, which must hold if it is possible for {\cal M} to be implemented arbitrarily closely by LOCC.

Here, we extend that analysis beyond measurements to encompass the much more general notion of quantum channels and indeed, any quantum instrument [11]. For channels, we find a condition, see Theorem 2 below, which must hold if quantum channel {\cal E} can be closely approximated by LOCC. The condition is a direct generalization of that for measurements in Theorem 1. It also involves characterization of the paths discussed in the preceding paragraph, which in similar fashion must lie within a more general geometric object ZZ_{\cal E}, which itself is uniquely determined by {\cal E}. The difference in structures of ZZ_{\cal E} and ZZ_{\cal M} reflects the fact that in implementing a quantum channel, no classical information about outcomes can be available, in contrast to measurements, where complete (classical) knowledge of those outcomes must be known. Both channels and measurements are special cases of quantum instruments, which make allowance for that classical information to vary from all to none.

We have found that dealing with ZZ_{\cal E} is more difficult as compared to ZZ_{\cal M}. As we shall see, the latter is a special case of the former, and it is this greater generality of ZZ_{\cal E} which makes it harder to analyze. Furthermore, whereas it was quite easy to find examples of previously studied measurements where we could apply, and learn from, our condition for measurements, Theorem 1, examples of quantum channels to which we can fruitfully apply Theorem 2 are much more difficult to find. Perhaps the explanation for this dearth of examples is that quantum channels—in particular, multipartite quantum channels—are just more difficult to study, in general, so comparatively little work has been done in this direction. One may hope that the results presented here will stimulate further efforts in that direction.

A dearth is not the same as a complete absence, however. In fact, there exists a two-qubit example of a quantum instrument for which it is known that it cannot be implemented exactly but can be approximated arbitrarily closely by LOCC [11]. The inspiration for this instrument came from ideas developed in a series of earlier papers on the random distillation of entanglement from a three-qubit system onto two of those qubits [23, 6, 9, 10]. We study the two-qubit example in detail below, first as a quantum channel, then as a general instrument, and finally as a quantum measurement. This analysis leads to insights regarding how the critical property of our paths, that of being piecewise local, see below, disappears in the limit as the number of rounds tends to infinity. It also illustrates other properties of these paths, such as the fact that certain matrices must have a block diagonal form determined by the structure of the given instrument in terms of its constituent completely positive (CP) maps. We also find from this example an interesting result regarding the closure of LOCC measurements.

The relationship between LOCC and its closure, LOCC¯\overline{\textrm{LOCC}}, has long been a point of interest to researchers. An early result along these lines was the demonstration, in the context of random distillation of two-party entanglement starting from the three qubit WW state, that a convergent sequence of transformations is not achievable by LOCC in the limit [10]. Later, utilizing a two-party generalization of those arguments, it was shown that LOCC is also not closed for two-qubit systems [11]. In each of these cases, the parties must rely upon classical information obtained in these protocols, leaving open the question of whether or not LOCC is closed when considering quantum channels, for which no classical information is retained. It has also been an open question whether LOCC is closed as applied to quantum measurements, where the maximum amount of classical information is needed to uniquely identify each and every measurement outcome. Here, in the context of the two-qubit example just mentioned, we discover a measurement which is not in LOCC but is within its closure. As far as we are aware, this is the first example of a measurement exhibiting such a property.

The remainder of the paper is organized as follows. In the next section, we provide an overview of the tools we will use in the subsequent sections, including a discussion of quantum channels and instruments and their mathematical representation, finite-round LOCC and its extension to infinite rounds, and the geometric objects, ZZ_{\cal M} and ZZ_{\cal E}, mentioned above. In Section III, we present our main result, Theorem 2. Section IV provides a proof of Theorem 2, which utilizes a result linking convergence of a sequence of quantum channels to convergence of a sequence of the geometric objects, ZZ_{\cal E}. The latter result is proved in Appendix A. In Section V the reader will find an extensive analysis, in the context of Theorem 2, of the two-qubit example channel from Ref. [11]. Section VI provides a generalization of Theorem 2 to the case of arbitrary quantum instruments, encompassing and interpolating between Theorem 1 and Theorem 2. There, the reader will also find a brief discussion in this context of the two-qubit example of Section IV. The important question of whether or not LOCC is a closed set, a question that may depend on the type of circumstances to which it is being applied, is addressed in Section VII where it is shown that LOCC is not closed when applied to quantum measurements. Section VIII presents a generalization of the example of Section IV to any number of qubits, and then finally, we offer our conclusions in Section IX.

II Preliminaries

In this section, we describe basic ideas underlying the various tools we will use for the analysis presented in this paper. We begin with a discussion of quantum channels and instruments, continue with local operations and classical communication (LOCC), followed by a discussion of infinite-round LOCC, and finally, the geometric objects we will need, known as zonoids.

II.1 Quantum Channels, Measurements, and Instruments

Let {\cal H} be a Hilbert space of finite dimension dd. Denote the set of bounded linear operators acting on {\cal H} as (){\cal B}({\cal H}). A quantum channel {\cal E} is a completely positive, trace-preserving (CPTP) linear map which sends density operators ρ()\rho\in{\cal B}({\cal H}) to density operators (ρ)(o){\cal E}(\rho)\in{\cal B}({\cal H}_{o}), where a density operator ρ\rho is positive semidefinite and has trace equal to unity. Any channel {\cal E} may be represented as follows. Consider a trio of Hilbert spaces, ,o,e{\cal H},{\cal H}_{o},{\cal H}_{e}, and an isometry VV, so that VV=IV^{\dagger}V=I_{{\cal H}}, V:oeV:{\cal H}\to{\cal H}_{o}\otimes{\cal H}_{e}, with II_{{\cal H}} the identity operator on {\cal H}. This maps operator ρ()\rho\in{\cal B}({\cal H}) to operator VρV(o)(e)V\rho V^{\dagger}\in{\cal B}({\cal H}_{o})\otimes{\cal B}({\cal H}_{e}). {\cal H} describes the input to the channel, o{\cal H}_{o} the output, and e{\cal H}_{e} is commonly referred to as the environment, which can be viewed as introducing noise into the evolution of the input system as it evolves to the output. Taking partial traces, we obtain two separate channels,

(ρ)\displaystyle{\cal E}(\rho) =Tre(VρV)\displaystyle=\textrm{Tr}_{e}\left(V\rho V^{\dagger}\right)
(ρ)\displaystyle{\cal F}(\rho) =Tro(VρV).\displaystyle=\textrm{Tr}_{o}\left(V\rho V^{\dagger}\right). (1)

Channel {\cal F}, which maps density operators on (){\cal B}({\cal H}) to those on (e){\cal B}({\cal H}_{e}), is called the complementary channel with reference to {\cal E}, and vice-versa. Note that output and environment spaces for {\cal F} are exchanged relative to those for {\cal E}, while the inputs are the same. One can always write

V=j=1de|jeKj=l=1do|loQl,\displaystyle V=\sum_{j=1}^{d_{e}}|j\rangle_{e}\otimes K_{j}=\sum_{l=1}^{d_{o}}|l\rangle_{o}\otimes Q_{l}, (2)

where |jee|j\rangle_{e}\in{\cal H}_{e}, while KjK_{j} maps states on {\cal H} to states on o{\cal H}_{o} and are referred to as Kraus operators [27] for channel {\cal E}. Similarly, |loo|l\rangle_{o}\in{\cal H}_{o} while Ql:eQ_{l}:{\cal H}\to{\cal H}_{e} and are Kraus operators for the complementary channel, {\cal F}. From Eq. (2), it is easily seen that l|Kj=j|Qleo{}_{o}\langle l|K_{j}=~{}_{e}\langle j|Q_{l}. This relation, along with the representation of {\cal E} and {\cal F} by VV, will be used in Appendix A.

In the main text of this paper, we will use an alternative representation for a quantum channel, which may be obtained from Eq. (II.1) and Eq. (2) in terms of Kraus operators as (ρ)=jKjρKj{\cal E}(\rho)=\sum_{j}K_{j}\rho K_{j}^{\dagger} (and (ρ)=lQlρQl{\cal F}(\rho)=\sum_{l}Q_{l}\rho Q_{l}^{\dagger}). The choice of such a representation is not unique. Rather, with κ\kappa the Kraus rank, which gives the minimum number of operators needed for a Kraus representation of {\cal E}, if {K^m}m=1κ\{\hat{K}_{m}\}_{m=1}^{\kappa} is a minimal representation, then the set {Kj}\{K_{j}\} is also a representation of {\cal E} as long as Kj=jWjmK^mK_{j}=\sum_{j}W_{jm}\hat{K}_{m} for every jj, where WW is any isometry, WW=IκW^{\dagger}W=I_{\kappa} and IκI_{\kappa} is the κ×κ\kappa\times\kappa identity matrix. There are no constraints on individual Kraus operators, but for channels, as a collection they must satisfy the condition jKjKj=VV=I\sum_{j}K_{j}^{\dagger}K_{j}=V^{\dagger}V=I_{{\cal H}}. Any given set of Kraus operators also represents a specific measurement, each KjK_{j} transforming state ρKjρKj\rho\to K_{j}\rho K_{j}^{\dagger} with probability pj=Tr(KjKjρ)p_{j}=\textrm{Tr}\left(K_{j}^{\dagger}K_{j}\rho\right). When there is a classical record providing index jj, we will refer to the set {Kj}\{K_{j}\} as a generalized measurement. The collection of operators {KjKj}\{K_{j}^{\dagger}K_{j}\} constitutes a positive-operator valued measure (POVM), which arises from a measurement after which the final state of the system is discarded, and each Ej=KjKjE_{j}=K_{j}^{\dagger}K_{j} is called a POVM element.

A quantum instrument is a set of completely positive (but generally not trace-preserving) maps, ={𝒥r}{\cal E}=\{{\cal J}_{r}\}, so that 𝒥r(ρ)=jSrKjρKj{\cal J}_{r}(\rho)=\sum_{j\in S_{r}}K_{j}\rho K_{j}^{\dagger}, and SrS_{r} is a subset of the full index set identifying those KjK_{j} that are associated with CP map 𝒥r{\cal J}_{r}. The collection of these CP maps is trace-preserving, rjSrKjKj=I\sum_{r}\sum_{j\in S_{r}}K_{j}^{\dagger}K_{j}=I_{\cal H}. A generalized measurement is equivalent to a fine-grained quantum instrument having one and only one Kraus operator for each CP map, whereas more general instruments may involve multiple Kraus operators for one or more of its CP maps. For example, from a measurement involving three operators, K1,K2,K3K_{1},K_{2},K_{3}, a fine-grained quantum instrument (generalized measurement) is the collection 𝒥r(ρ)=KrρKr,r=1,2,3{\cal J}_{r}(\rho)=K_{r}\rho K_{r}^{\dagger},~{}r=1,2,3. On the other hand, the pair of CP maps, 𝒥1(ρ)=K1ρK1+K2ρK2{\cal J}_{1}^{\prime}(\rho)=K_{1}\rho K_{1}^{\dagger}+K_{2}\rho K_{2}^{\dagger} and 𝒥2(ρ)=K3ρK3{\cal J}_{2}^{\prime}(\rho)=K_{3}\rho K_{3}^{\dagger} together also constitute a quantum instrument, being a coarse-grained version of the former. When the sum over terms is unrestricted, the instrument involves only a single CP map, which is then trace-preserving and is therefore, itself, a quantum channel. On the other hand, every quantum instrument is uniquely associated with a quantum channel through the identification,

(ρ)=r𝒥r(ρ)[r]c.\displaystyle{\cal E}(\rho)=\sum_{r}{\cal J}_{r}(\rho)\otimes[r]_{c}. (3)

The projector [r]c[r]_{c} is to be understood as a register holding the classical information about the various CP maps, since measuring this register in the |rc|r\rangle_{c} basis provides knowledge of which of these CP maps was implemented. These channels are referred to as quantum-classical maps.

II.2 Local Operations and Classical Communication

In this paper, we are concerned with quantum instruments on multipartite systems, where a multipartite system on PP parties is described by a tensor product Hilbert space, =αα{\cal H}=\bigotimes_{\alpha}{\cal H}_{\alpha}, α{A,B,}\alpha\in\{A,B,\ldots\}, and we shall refer to α{\cal H}_{\alpha} as the local Hilbert space of party α\alpha. We wish to study the question of whether or not channel (or instrument) {\cal E}, acting on multipartite system {\cal H} can be closely approximated when the parties are restricted to using local operations and classical communication. That is, can parties A,B,A,B,\ldots closely simulate {\cal E} by making local measurements in their individual laboratories, communicating the outcomes of those measurements to the other parties, continuing in this way for as many rounds of measurement and communication as they wish, and then at the very end, forgetting the outcomes of those local measurements? To analyze these questions, it is common to represent an LOCC protocol by a tree graph. The root of the tree corresponds to the situation before any party has done anything, so we label this node by the identity operator, II_{\cal H}. For each local measurement by any given party, each outcome is represented by a node with an edge extending upward in the tree (in the direction of the root) to the node representing the situation just before that measurement was performed. The latter node is known as the parent of its collection of child nodes. A finite branch stretches from the root node to a leaf node, which has no children, but a protocol may also have infinite branches, which since they do not end, have no leaf nodes. Since these trees represent LOCC protocols, for which only one party measures at a time, the difference between parent and child is local. To see this explicitly suppose, as we have found it productive to do [14, 16, 17, 18, 20], that each node in the tree is labeled by the POVM element corresponding to the action of all parties up to that point in the protocol. That is, if the parent node is labeled as the POVM element 𝒜{\cal A}\otimes{\cal B} and Alice measures next, for example, then the child node will be 𝒜{\cal A}^{\prime}\otimes{\cal B}; only Alice’s local part has changed. Finally, we recall a lemma from [15], which applies to all finite-round protocols where each node is labeled by a positive semidefinite operator, as just described.

Lemma 1.

[15]Each node nn in a finite-round LOCC tree is equal to the sum of the collection of all leaf nodes that are descendant from that node nn.

The foregoing discussion may give the appearance that only POVMs are being considered, but one can show that any generalized measurement consisting of operators Kj=UjKjKjK_{j}^{\prime}=U_{j}\sqrt{K_{j}^{\dagger}K_{j}}, with UjUj=IU_{j}^{\dagger}U_{j}=I_{\cal H}, can be implemented by essentially the same LOCC protocol as that implementing POVM elements KjKjK_{j}^{\dagger}K_{j}, as long as each UjU_{j} is a product operator of the form UjAUjBU_{jA}\otimes U_{jB}\otimes\ldots. This is because product operators are the only ones LOCC protocols are capable of implementing, so since KjK_{j} and KjK_{j}^{\prime} must both be product operators, then UjU_{j} must be one, as well. In order to implement the KjK_{j}^{\prime}, the parties simply perform the same protocol as that implementing the KjK_{j}, after which they add one last step at the end, each party α\alpha performing UjαU_{j\alpha} after KjKj\sqrt{K_{j}^{\dagger}K_{j}} has been obtained (here viewing the original protocol as implementing a generalized measurement rather than a POVM). Since the reverse implication is obvious, that implementation of the KjK_{j} implies that the KjKjK_{j}^{\dagger}K_{j} can also be implemented, we see that as long as all the UjU_{j} are product unitaries, then the generalized measurement can be implemented by LOCC if and only if the corresponding POVM can be so implemented. A couple of other points are worth noting. First, even though a quantum channel can be represented by many alternative sets of Kraus operators, any given LOCC protocol implements a specific such set, so may equally well be viewed as having implemented the particular measurement corresponding to that set. We will use this observation below and will refer to each corresponding measurement {\cal M} as being compatible with the given channel (or instrument) {\cal E}. Secondly, we note that the outcome of each local measurement in an LOCC protocol must be communicated to the other parties, and any given implementation of the protocol will follow a single branch with a string of local outcomes leading to a given final outcome, where these outcomes have all been communicated and therefore, recorded. Technically, however, when the outcome is known, one has failed to implement the channel. This is why, in discussing LOCC implementation of quantum channels, we imagine that at the end of the protocol, the parties forget the collection of outcomes obtained in each of the individual local measurements. Similarly, for implementing quantum instruments, other than the fully fine-grained measurements, one forgets all classical information other than index rr identifying CP map 𝒥r{\cal J}_{r}.

II.3 Asymptotic LOCC and convergence of quantum channels

To study what happens with asymptotic LOCC, where an infinite number of rounds is allowed, one considers the infinite-round protocol to be the limit of a given sequence of finite-round protocols. As discussed in Ref. [11], there are two distinct kinds of such sequences. In the first case, each protocol in the sequence is obtained from the protocol immediately preceding it by adding one (or more) additional rounds. Such protocols are still LOCC. Let us now discuss the other case, in which the next protocol in the sequence again has additional rounds but is also allowed to differ from its preceding protocol even in earlier rounds. As such, apart from the assumption of convergence, there need not be a specific relationship between successive protocols of the sequence. With either type of asymptotic LOCC, if the instruments implemented by the sequence of protocols have a limiting instrument, then we say that the latter is in, or can be implemented by, the closure of LOCC, denoted LOCC¯\overline{\textrm{LOCC}}. That is, if protocol 𝒫ν{\cal P}^{\nu} implements instrument ν{\cal E}^{\nu} and limνν=\lim_{\nu\to\infty}{\cal E}^{\nu}={\cal E}, then we shall say that instrument {\cal E} can be implemented by LOCC¯\overline{\textrm{LOCC}}: LOCC¯{\cal E}\in\overline{\textrm{LOCC}}. What these limits mean is that for every real ϵ>0\epsilon>0 there exists natural number μ\mu such that for every natural number ν>μ\nu>\mu,

ν<ϵ,\displaystyle\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}<\epsilon, (4)

where the diamond norm distance is defined as

ν=maxX11()X(ν)X1,\displaystyle\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}=\max_{\left\|X\right\|_{1}\leq 1}\left\|({\cal I}\otimes{\cal E})X-({\cal I}\otimes{\cal E}^{\nu})X\right\|_{1}, (5)

with {\cal I} the identity channel on (){\cal B}({\cal H}), XX is an operator on ()(){\cal B}({\cal H})\otimes{\cal B}({\cal H}), and 1\left\|\cdot\right\|_{1} is the trace norm, X1=Tr(XX)\left\|X\right\|_{1}=\textrm{Tr}\left(\sqrt{X^{\dagger}X}\right). A lower bound for this distance may be obtained in terms of these channels’ Choi-Jamiolkowski operators, Φ\Phi and Φν\Phi^{\nu}, respectively. This is done by choosing a maximally entangled state on ()(){\cal B}({\cal H})\otimes{\cal B}({\cal H}) in Eq. (5), X=|ΩΩ|=ij|ij||ij|/dX=|\Omega\rangle\langle\Omega|=\sum_{ij}|i\rangle\langle j|\otimes|i\rangle\langle j|/d, and noting the definition of the Choi-Jamiolkowski operator for {\cal E} as Φ=()|ΩΩ|\Phi=({\cal I}\otimes{\cal E})|\Omega\rangle\langle\Omega|, with an identical relation between Φν\Phi^{\nu} and ν{\cal E}^{\nu} ({\cal E} and ν{\cal E}^{\nu} must act on the same input Hilbert space if a sequence of the latter is to converge to the former). Inserting this in Eq. (4) yields,

ΦΦν1<ϵ.\displaystyle\left\|\Phi-\Phi^{\nu}\right\|_{1}<\epsilon. (6)

Therefore, convergence of the channels to {\cal E} implies convergence of the Choi-Jamiolkowski operators to Φ\Phi.

II.4 Zonotopes and Zonoids

We will make use of certain geometric objects known as zonoids, which are limits of objects known as zonotopes. A zonotope is a closed, convex set, the Minkowski sum of a finite number, NN, of line segments, [0,Ej][0,E_{j}], written Z=j[0,Ej]Z=\sum_{j}[0,E_{j}]. An alternative definition is Z={z|z=jcjEj,0cj1}Z=\left\{z\left|z=\sum_{j}c_{j}E_{j},0\leq c_{j}\leq 1\right.\right\}. Zonoids are limits of sequences of zonotopes as NN\to\infty. We have utilized zonotopes, ZZ_{\cal M}, generated by POVM elements EjE_{j}\in{\cal M}, so that the EjE_{j} are positive semidefinite operators in a Hilbert space, (){\cal B}({\cal H}), the sum of which is equal to II_{\cal H} [18, 20]. The zonoids, ZZ_{\cal E}, used here will be generated by a minimal Kraus representation, {K^m}\{\hat{K}_{m}\}, of a quantum channel {\cal E} in the sense that

Z={z|z=mmC^mmK^mK^m,IκC^0},\displaystyle Z_{\cal E}=\left\{z\left|z=\sum_{mm^{\prime}}\hat{C}_{mm^{\prime}}\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}},I_{\kappa}\geq\hat{C}\geq 0\right.\right\}, (7)

with C^\hat{C} a positive semidefinite κ×κ\kappa\times\kappa matrix, κ\kappa is the Kraus rank of channel {\cal E}, and IκI_{\kappa} is the κ×κ\kappa\times\kappa identity matrix. Since positive semidefinite matrices can always be diagonalized, we can write C^mm=jWjmcjWjm\hat{C}_{mm^{\prime}}=\sum_{j}W_{jm}^{\ast}c_{j}W_{jm^{\prime}} for some isometry, WW, so that zz in the above definition of ZZ_{\cal E} becomes

z\displaystyle z =jcj(mWjmK^m)(mWjmK^m)=jcjKjKj=jcjEj,\displaystyle=\sum_{j}c_{j}\left(\sum_{m}W_{jm}^{\ast}\hat{K}_{m}^{\dagger}\right)\left(\sum_{m^{\prime}}W_{jm^{\prime}}\hat{K}_{m^{\prime}}\right)=\sum_{j}c_{j}K_{j}^{\dagger}K_{j}=\sum_{j}c_{j}E_{j}, (8)

Since Eq. (8) holds for every set of coefficients, {cj}\{c_{j}\} and every isometry WW, we see that if zZz\in Z_{\cal M}, then zZz\in Z_{\cal E} as well, and we have thus proved the following lemma.

Lemma 2.

Given zonoid ZZ_{{\cal E}} as defined in Eq. (7) and zonotope ZZ_{{\cal M}} for any {\cal M} compatible with {\cal E} in the sense that measurement operators KjK_{j} of {\cal M} constitute a Kraus representation of {\cal E}, then ZZZ_{\cal M}\subseteq Z_{\cal E}.

The following two examples may provide the reader with a bit of intuition about these geometric objects. First, consider a simple quantum channel on a single qubit, d=2d=2, represented by Kraus operators that are rank-11 projectors, {[0],[1]}\{[0],[1]\}, with [ψ]=|ψψ|[\psi]=|\psi\rangle\langle\psi|. Operators Kl=Vl0[0]+Vl1[1]K_{l}=V_{l0}[0]+V_{l1}[1], with l=1,,Nl=1,\ldots,N for any NN, also represent this same quantum channel, as long as VV=I2V^{\dagger}V=I_{2}, the 2×22\times 2 identity matrix. Then, for any zZz\in Z_{\cal E}, we have that

z=C00[0]+C11[1].\displaystyle z=C_{00}[0]+C_{11}[1]. (9)

Since [0][1][0][1] vanishes, there is no contribution from the off-diagonal entries of CC. Thus, ZZ_{\cal E} is generated by a finite number (two) of line segments, these segments stretching from the zero operator to [0][0] and [1][1], respectively. Therefore, ZZ_{\cal E} is in fact a zonotope, rather than a zonoid, and has the shape of a square.

The second example is again on a single qubit, in this case represented by the two Kraus operators |00|,|01||0\rangle\langle 0|,|0\rangle\langle 1|. In this case, any set of operators, |0ψl||0\rangle\langle\psi_{l}|, satisfying l[ψl]=I2\sum_{l}[\psi_{l}]=I_{2} also represents this same channel. Any zZz\in Z_{\cal E} is of the form

z=C00[0]+C01|01|+C10|10|+C11[1],\displaystyle z=C_{00}[0]+C_{01}|0\rangle\langle 1|+C_{10}|1\rangle\langle 0|+C_{11}[1], (10)

showing that similarly to CC, the only constraint on zz is that 0zI0\leq z\leq I_{\cal H}. Then, ZZ_{\cal E} is a four-dimensional cone-like object which, denoting the fourth dimension as r=Tr(z)r=\textrm{Tr}\left(z\right), stretches from 0r20\leq r\leq 2, and every three-dimensional slice through ZZ_{\cal E} and perpendicular to the rr-axis is a (Bloch-like) sphere of radius rr. Given these spherical cross-sections, ZZ_{\cal E} can only be generated by NN line segments in the limit of NN\to\infty, so this object is a true zonoid.

II.5 Notation

A heads-up for the reader about notation: superscripts that are Greek letters, such as μ,ν\mu,\nu, are not exponents, but are labels denoting position in a sequence; see for example, superscript ν\nu in the various equations of Section II.3, above. Those superscripts that are lowercase letters from the Roman/Latin alphabet, such as nn, represent exponents. Subscripts are always just labels, and all lowercase Roman/Latin labels may be assumed to be natural numbers.

III Approximating quantum channels by LOCC

In this section we present our main result, a necessary condition that LOCC¯{\cal E}\in\overline{\textrm{LOCC}}, Theorem 2, the proof of which is given in the subsequent section. We begin by recalling the main theorem of Ref. [18]. There, we proved that if measurement ={Ej}{\cal M}=\{E_{j}\} can be approximated by LOCC arbitrarily closely, then for each jj, there exists a continuous, monotonic path, Πj(s)\Pi_{j}(s), of positive semidefinite product operators from II_{\cal H} to each outcome EjE_{j} of that measurement and lying entirely within the zonotope, ZZ_{{\cal M}}, generated by those POVM elements, {Ej}\{E_{j}\}. Here, monotonicity means that the trace of Πj(s)\Pi_{j}(s) is non-increasing along the path. The proof of that result relied on the observation that every finite-round LOCC protocol implements a measurement for which there exists a continuous, monotonic, piecewise local111The notion of a piecewise local path is borrowed from the piecewise constant curves we all learned about when studying calculus. Instead of being constant, however, the pieces here correspond to local measurement outcomes and correspond to the edges between a parent node and one of its children in an LOCC tree. Along this edge, only one party’s part of the full measurement operator changes, justifying use of the word local. See Ref. [18] for a more detailed explanation of this concept. path of positive semidefinite product operators to each of the outcomes of that measurement. Using this observation, it was then argued that from any sequence of finite-round LOCC protocols implementing measurements ν{\cal M}^{\nu} for which limνν=\lim_{\nu\to\infty}{\cal M}^{\nu}={\cal M}, there exists a sequence of piecewise local paths whose limit is a (not necessarily piecewise local) path lying within ZZ_{{\cal M}} and ending at one of the outcomes of {\cal M}, and that this holds for each such outcome. The resulting theorem reads as follows.

Theorem 1.

[18] If LOCC¯{\cal M}\in\overline{\textrm{LOCC}}, with measurement {\cal M} consisting of POVM elements EjE_{j}, then for each jj, there exists a continuous, monotonic path of product operators from {\cal I}_{\cal H} to a point on the (half-open) line segment (0,Ej](0,E_{j}], and this path lies entirely within zonotope Z=j[0,Ej]Z_{\cal M}=\sum_{j}[0,E_{j}].

As argued above, the theorem holds whether one views {\cal M} as a POVM or as a generalized measurement. In either case, it is the POVM elements EjE_{j} that determine zonotope ZZ_{\cal M} and the endpoints of the consequent paths.

In this section, we are concerned with quantum channels rather than measurements. Since any LOCC protocol implementing channel {\cal E} may just as well be viewed as implementing the generalized measurement consisting of a specific set of Kraus operators which represents {\cal E}, Theorem 1 provides a solid foundation for the following theorem, which is our main result, providing a necessary condition that a quantum channel may be simulated arbitrarily closely by LOCC.

Theorem 2.

Given quantum channel {\cal E} acting on input space (){\cal B}({\cal H}) and represented by the minimal set of κ\kappa Kraus operators, {K^m}\{\hat{K}_{m}\}, then if {\cal E} can be implemented by LOCC with vanishingly small error, the following conditions must hold:

  1. 1.

    There exists a set, {Πj(s)}\{\Pi_{j}(s)\}, of continuous, monotonic paths of positive semidefinite product operators, each of which begins at II_{\cal H}, the identity operator on {\cal H}, and ends at a positive semidefinite product operator of the form, Ej=mm(C^j)mmK^mK^mE_{j}=\sum_{mm^{\prime}}(\hat{C}_{j})_{mm^{\prime}}\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}}.

  2. 2.

    Each C^j\hat{C}_{j} is a positive semidefinite matrix, C^j0\hat{C}_{j}\geq 0, and may be chosen to have rank equal to unity, with the collection satisfying jC^j=Iκ\sum_{j}\hat{C}_{j}=I_{\kappa}, where IκI_{\kappa} is the κ×κ\kappa\times\kappa identity matrix and κ\kappa is the Kraus rank of {\cal E}.

  3. 3.

    Each of these paths of product operators lies entirely within the zonoid, ZZ_{{\cal E}}, consisting of all positive semidefinite linear combinations of the operators, K^mK^m\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}}: Z={z|z=mmC^mmK^mK^m,IκC^0}Z_{{\cal E}}=\left\{z\left|z=\sum_{mm^{\prime}}\hat{C}_{mm^{\prime}}\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}},I_{\kappa}\geq\hat{C}\geq 0\right.\right\}, with 1rank(C^)κ1\leq\textrm{rank}(\hat{C})\leq\kappa.

Note that a positive semidefinite linear combination simply means a linear combination of the form given in the above definition of ZZ_{{\cal E}}, that is, a linear combination where the coefficients constitute a positive semidefinite matrix. In addition, the paths are parametrized by s=Tr(Πj(s))s=\textrm{Tr}\left(\Pi_{j}(s)\right), whose range is ds0d\geq s\geq 0 with dd the dimension of input Hilbert space {\cal H}; see [18] for a more detailed explanation of this parametrization.

The condition of Theorem 2 that {\cal E} can be implemented by LOCC with vanishingly small error means there exists a sequence of finite-round protocols 𝒫ν{\cal P}^{\nu} implementing quantum channel ν{\cal E}^{\nu} such that limνν=\lim_{\nu\to\infty}{\cal E}^{\nu}={\cal E}. Since as discussed in Section II.2, each of these finite-round protocols may equally well be viewed as implementing a measurement ν{\cal M}^{\nu}, which implies the existence of the concomitant paths of product operators in Theorem 1 [18], this (at least) suggests the existence of the set of limiting paths alluded to in Theorem 2. Thus, one should not be surprised to find that the results for approximating measurements by LOCC [18] can be generalized, as we do here, to apply in similar fashion to approximating quantum channels by LOCC. One notable difference was discussed in Section II.4, that while the approximation of measurements involves paths lying entirely within a zonotope, ZZ_{{\cal M}}, for the approximation of channels one must utilize the zonoid, ZZ_{{\cal E}}.

Recall from Section II.4 that in order to implement a quantum channel, say ν{\cal E}^{\nu}, by LOCC, one may implement any measurement ν{\cal M}^{\nu} compatible with ν{\cal E}^{\nu}. This means that simply knowing one can implement a given channel leaves one without knowledge of which compatible measurement(s) can be implemented. Thus, we also do not know the set of Kraus operators implemented, whereas that set is what determines the endpoints of the paths introduced in Theorem 1 [18]. Therefore, we are left with the question, how do we know what the endpoints of all the paths of Theorem 2 are? The answer is, we do not. Instead, as is indicated in Item 2 of Theorem 2, we only know that each of these endpoints corresponds to a rank-11 matrix C^j\hat{C}_{j} and the collection of these matrices constitutes a rank-11 decomposition of the κ×κ\kappa\times\kappa identity matrix. This is the extent to which we are able to characterize those endpoints without knowing in advance what specific set of Kraus operators is to be implemented. Nonetheless, this is sufficient to ensure that the resulting set of Kraus operators are a representation of {\cal E}.

We now turn to the proof of Theorem 2.

IV Proof of Theorem 2

We begin the proof of Theorem 2 with the following lemma showing that convergence of a sequence of quantum channels implies convergence of the corresponding sequence of zonoids.

Lemma 3.

If there exists a sequence of quantum channels, ν{\cal E}^{\nu} of Kraus rank κν\kappa^{\nu} and with minimal Kraus representation {K^mν}\{\hat{K}_{m}^{\nu}\}, converging to channel {\cal E} of Kraus rank κ\kappa and with minimal Kraus representation {K^m}\{\hat{K}_{m}\}limνν=\lim_{\nu\to\infty}{\cal E}^{\nu}={\cal E}—then there also exists a sequence of zonoids, Zν={zν|zν=mmC^mmνK^mνK^mν,IκνC^ν0}Z_{{\cal E}^{\nu}}=\left\{z^{\nu}\left|z^{\nu}=\sum_{mm^{\prime}}\hat{C}^{\nu}_{mm^{\prime}}\hat{K}_{m}^{\nu{\dagger}}\hat{K}_{m^{\prime}}^{\nu},I_{\kappa^{\nu}}\geq\hat{C}^{\nu}\geq 0\right.\right\} associated with channel ν{\cal E}^{\nu}, which converges to zonoid Z={z|z=mmC^mmK^mK^m,IκC^0}Z_{\cal E}=\left\{z\left|z=\sum_{mm^{\prime}}\hat{C}_{mm^{\prime}}\hat{K}_{m}^{{\dagger}}\hat{K}_{m^{\prime}},I_{\kappa}\geq\hat{C}\geq 0\right.\right\} associated with channel {\cal E}. That is, limνZν=Z\lim_{\nu\to\infty}Z_{{\cal E}^{\nu}}=Z_{\cal E}.

The proof is given in Appendix A.

The next step is to recall that each LOCC protocol 𝒫ν{\cal P}^{\nu}, which implements channel ν{\cal E}^{\nu}, may also be thought of as implementing a POVM, ν{\cal M}^{\nu}, consisting of elements Ejν=KjνKjνE_{j}^{\nu}=K_{j}^{\nu{\dagger}}K_{j}^{\nu}, where with WνW^{\nu} an isometry, Kjν=mWjmνK^mνK_{j}^{\nu}=\sum_{m}W_{jm}^{\nu}\hat{K}_{m}^{\nu} is a Kraus operator in a collection representing channel ν{\cal E}^{\nu}. Therefore, by Theorem 1, we know that there exists a continuous, monotonic path of product operators from II_{\cal H} to a point on the half-open line segment (0,Ejν](0,E_{j}^{\nu}], for each outcome, jj, and that these paths lie entirely within the zonotope ZνZ_{{\cal M}^{\nu}} generated by the collection of these operators, Ejν{E_{j}^{\nu}}. In fact, by Theorem 22 of Ref. [18], these paths are piecewise local. As discussed at the end of the preceding section, we do not know ν{\cal M}^{\nu}. Therefore, we consider all possible ν{\cal M}^{\nu} compatible with channel ν{\cal E}^{\nu}, which means that instead of considering our continuous paths as lying within zonotope ZνZ_{{\cal M}^{\nu}}, we instead view these paths as lying within zonoid ZνZ_{{\cal E}^{\nu}}. According to Lemma 2, ZνZνZ_{{\cal M}^{\nu}}\subseteq Z_{{\cal E}^{\nu}} for every such ν{\cal M}^{\nu}, so this is always possible. We could, of course, restrict to measurements that are implementable by LOCC, so that such paths of product operators actually exist, but it is sufficient to follow the easier route and include all compatible measurements, as we do here.

Thus, we have a sequence of zonoids, {Zν}\{Z_{{\cal E}^{\nu}}\}, which converges to zonoid ZZ_{\cal E}, and within each ZνZ_{{\cal E}^{\nu}}, we have a collection of continuous paths of product operators starting at II_{\cal H}. This situation is very similar to the one encountered in Ref. [18], including that the paths are all Lipschitz continuous with Lipschitz constant equal to unity [18]. By the Arzéla-Ascoli theorem [22], this tells us that there exists a subsequence of these paths which converges uniformly on the range, s[d,0]s\in[d,0], to a continuous limiting path. However, to this point, we have many paths for each ν\nu and have not actually identified a sequence of paths to begin with. We now turn to the task of identifying such a sequence. In fact we will identify, and characterize, many such sequences.

Whereas in our earlier work, we had a sequence of known measurements, ν{\cal M}^{\nu}, whose limit was another known measurement, {\cal M}, here we do not know what the ν{\cal M}^{\nu} or {\cal M} are. Instead, we only know that each ν{\cal M}^{\nu} is compatible with channel ν{\cal E}^{\nu}, the sequence of which has the limit {\cal E}. In the previous case, knowing ν{\cal M}^{\nu} meant that we also knew the final outcomes, EjνE_{j}^{\nu}, of the corresponding LOCC protocol, and also that there was an ordering of these outcomes, provided by labels ν\nu for fixed jj, such that for each jj, the sequence {Ejν}ν\{E_{j}^{\nu}\}_{\nu} satisfied limνEjν=Ej\lim_{\nu\to\infty}E_{j}^{\nu}=E_{j}, with EjE_{j} an outcome of the limiting measurement {\cal M}. This allowed us to identify sequences of paths—i.e., those ending at EjνE_{j}^{\nu} with the aforementioned fixed label, jj—from which, according to Arzéla-Ascoli, it followed that there was a limiting path from II_{\cal H} to a point along (0,Ej](0,E_{j}]. In our present case, with limνν=\lim_{\nu\to\infty}{\cal E}^{\nu}={\cal E}, we do not know EjνE_{j}^{\nu} or the limiting outcomes EjE_{j}. Thus, we need a way to partition the various paths we’ve identified as lying within the ZνZ_{{\cal E}^{\nu}} so that each subset provides a sequence of paths whose limit we can characterize. In particular, we wish to characterize the endpoints of the limiting paths. We will discuss how to partition the paths forthwith. The following lemma concerning the Kraus ranks of the ν{\cal E}^{\nu} will be useful, and it may be of interest in other contexts, as well.

Lemma 4.

Consider the sequence of quantum channels, {ν}\{{\cal E}^{\nu}\}, with the Kraus rank of ν{\cal E}^{\nu} being κν\kappa^{\nu}. If limνν=\lim_{\nu\to\infty}{\cal E}^{\nu}={\cal E} with Kraus rank of {\cal E} equal to κ\kappa, then there exists μ\mu such that κνκ\kappa^{\nu}\geq\kappa for all ν>μ\nu>\mu.

Proof.

We begin the proof by considering the Choi-Jamiolkowski operator Φ\Phi (Φν)(\Phi^{\nu}) for {\cal E} (ν)({\cal E}^{\nu}), defined below Eq. (5). It is well known that the rank of Φ\Phi is equal to the Kraus rank κ\kappa of {\cal E}, and similarly for Φν\Phi^{\nu}. By the Eckart-Young-Mirsky low rank approximation theorem [21], which holds for all unitarily invariant norms [30], we have that when rank(Φν)<κ(\Phi^{\nu})<\kappa, then ΦΦν1σκ\left\|\Phi-\Phi^{\nu}\right\|_{1}\geq\sigma_{\kappa}, where σκ\sigma_{\kappa} is the κ\kappath largest singular value of Φ\Phi. Note that since {\cal E}, and thus Φ\Phi, is known, σκ\sigma_{\kappa} is some fixed, positive number. Choose ϵ\epsilon to satisfy σκ>ϵ>0\sigma_{\kappa}>\epsilon>0 and then identify μ\mu such that for every ν>μ\nu>\mu, ϵ>ν\epsilon>\left\|{\cal E}^{\nu}-{\cal E}\right\|_{\diamond}. Since limνν=\lim_{\nu\to\infty}{\cal E}^{\nu}={\cal E}, this is always possible.

The proof of the lemma now proceeds by contradiction. Suppose for some ν>μ\nu>\mu, κν<κ\kappa^{\nu}<\kappa. Recalling from the discussion around Eqs. (4) and (6) that νΦΦν1\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}\geq\left\|\Phi-\Phi^{\nu}\right\|_{1}, then for every ϵ<σκ\epsilon<\sigma_{\kappa}, we have

σκ>ϵ>νΦΦν1σκ,\displaystyle\sigma_{\kappa}>\epsilon>\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}\geq\left\|\Phi-\Phi^{\nu}\right\|_{1}\geq\sigma_{\kappa}, (11)

a contradiction. Therefore, either ν<μ\nu<\mu or κνκ\kappa^{\nu}\geq\kappa. That is to say, for each ϵ>0\epsilon>0, there exists μ\mu such that κνκ\kappa^{\nu}\geq\kappa for all ν>μ\nu>\mu, and the proof is complete.\blacksquare

The meaning of this result is clear: a sequence of low rank operators cannot become arbitrarily close to an operator of higher rank. On the other hand, it is clearly possible for higher rank operators to become closer and closer to an operator of lower rank.

We now turn to a characterization of endpoints of our paths. To begin, notice that each branch in finite-round LOCC protocol 𝒫ν{\cal P}^{\nu} ends at a leaf node corresponding to POVM element EjνE_{j}^{\nu}, or equivalently at measurement operator, Kjν=VνEjν=lWjlνK^lνK_{j}^{\nu\prime}=V^{\nu}\sqrt{E_{j}^{\nu}}=\sum_{l}W_{jl}^{\nu}\hat{K}_{l}^{\nu}, where VνV^{\nu} and WνW^{\nu} are isometries (acting on different spaces, of course) and the K^lν\hat{K}_{l}^{\nu} constitute a minimal Kraus representation of the quantum channel ν{\cal E}^{\nu} implemented by 𝒫ν{\cal P}^{\nu}. Then,

Ejν\displaystyle E_{j}^{\nu} =llWjlνWjlνK^lνK^lν\displaystyle=\sum_{ll^{\prime}}W_{jl}^{\nu\ast}W_{jl^{\prime}}^{\nu}\hat{K}_{l}^{\nu{\dagger}}\hat{K}_{l^{\prime}}^{\nu}
=ll(C^jν)llK^lνK^lν,\displaystyle=\sum_{ll^{\prime}}(\hat{C}_{j}^{\nu})_{ll^{\prime}}\hat{K}_{l}^{\nu{\dagger}}\hat{K}_{l^{\prime}}^{\nu}, (12)

with C^jν=WjνWjν\hat{C}_{j}^{\nu}=\vec{W}_{j}^{\nu}\vec{W}_{j}^{\nu{\dagger}} a rank-11 matrix, and Wjν\vec{W}_{j}^{\nu} is the jjth column of isometry WνW^{\nu}. Note that from Ref.  [18], we know that each branch corresponds directly to a continuous monotonic path of product operators lying entirely within ZνZ_{{\cal M}^{\nu}}, and thus within ZνZ_{{\cal E}^{\nu}}. We see here that each of these paths ends at a point corresponding to a rank-11 matrix, C^jν\hat{C}_{j}^{\nu}.

Let us consider the collection of these endpoints and the paths Πjν(s)\Pi_{j}^{\nu}(s) to each of them. Following Ref. [18], we parametrize these paths by s=Tr(Πjν(s))s=\textrm{Tr}\left(\Pi_{j}^{\nu}(s)\right). Define NνN^{\nu} to be the number of these paths, equal to the number of leaf nodes, j=1,,Nν,j=1,\ldots,N^{\nu}, for each protocol 𝒫ν{\cal P}^{\nu}. By the results of Ref. [17], we can assume that NνN^{\nu} is finite for all ν<\nu<\infty. Consider the first λ\lambda of these protocols, ν=1,,λ\nu=1,\ldots,\lambda, and set N0(λ)=maxνλNνN_{0}(\lambda)=\max_{\nu\leq\lambda}{N^{\nu}}. For each ν\nu such that Nν<N0(λ)N^{\nu}<N_{0}(\lambda), split up some of the paths into pjνp_{j}^{\nu} duplicates. We can do this splitting at the very final (local) measurement of each path that is split. For example, if the final outcome of Πjν(s)\Pi_{j}^{\nu}(s) is EjνE_{j}^{\nu}, replace that single outcome by pjνp_{j}^{\nu} individual outcomes, each of which is Ejν/pjνE_{j}^{\nu}/p_{j}^{\nu}. This increases the number of outcomes of that final local measurement by pjν1p_{j}^{\nu}-1. It doesn’t matter which final outcomes are duplicated, or how many parts we choose to split each of them into, as long as we end up with a new number of paths such that Nν=N0(λ)N^{\nu\prime}=N_{0}(\lambda) for each ν\nu. Note that each of these paths correspond to a branch in an LOCC protocol, 𝒫ν{\cal P}^{\nu\prime}, which is only slightly modified from a branch in 𝒫ν{\cal P}^{\nu} by splitting a subset of the final outcomes as just described. Redefine each of these final segments of the duplicated paths as Πjν(s)=Πjν(s)/pjν\Pi_{j}^{\nu\prime}(s)=\Pi_{j}^{\nu}(s)/p_{j}^{\nu}, for ss in the appropriate range for that final local measurement (the paths are not altered for the remainder of the range of ss). This maintains the condition that the sum of POVM elements corresponding to the entire collection of leaf nodes for protocol 𝒫ν{\cal P}^{\nu\prime} is equal to II_{\cal H}, as it must be. We now have that each protocol 𝒫ν{\cal P}^{\nu\prime} has N0(λ)N_{0}(\lambda) final outcomes, with a continuous path Πjν(s)\Pi_{j}^{\nu\prime}(s) to each one, and each of these final outcomes corresponds to a positive semidefinite operator satisfying Eq. (IV) with matrix C^jν\hat{C}_{j}^{\nu} also positive semidefinite and having rank equal to unity. Since the sum of all final outcomes is equal to II_{\cal H}, there is always a choice of these matrices C^jν\hat{C}_{j}^{\nu} such that they are a rank-11 decomposition of the identity, IκνI_{\kappa^{\nu}} (if, for example, the KlνKlνK_{l}^{\nu{\dagger}}K_{l^{\prime}}^{\nu} are linearly dependent, then there can be other possible choices for the C^\hat{C}-matrices, but we are specifically looking for those that are rank-11).

For each j=1,,N0(λ)j=1,\ldots,N_{0}(\lambda), we have λ\lambda paths, Πjν(s),ν=1,,λ\Pi_{j}^{\nu}(s),~{}\nu=1,\ldots,\lambda (dropping the primes, for convenience). Since N0(λ)N_{0}(\lambda) is a maximum over the first λ\lambda protocols, then N0(λ)N0(λ)N_{0}(\lambda^{\prime})\geq N_{0}(\lambda) whenever λ>λ\lambda^{\prime}>\lambda; N0(λ)N_{0}(\lambda) is a non-decreasing function of λ\lambda. When incrementing λλ+1\lambda\to\lambda+1, one path is added for each value of jN0(λ)j\leq N_{0}(\lambda). If N0(λ+1)>N0(λ)N_{0}(\lambda+1)>N_{0}(\lambda), additional values of jj are then added (by splitting outcomes as discussed in the preceding paragraph) up to N0(λ+1)N_{0}(\lambda+1), along with λ+1\lambda+1 paths for each new value of jj, one from each of the protocols 𝒫ν{\cal P}^{\nu} for ν=1,,λ+1\nu=1,\ldots,\lambda+1. It is possible that N^0=limλN0(λ)\hat{N}_{0}=\lim_{\lambda\to\infty}N_{0}(\lambda) diverges in the limit λ\lambda\to\infty (see Section V for a detailed analysis of such a case) but in any event, we end up with an infinite number of paths for each j=1,,N^0j=1,\ldots,\hat{N}_{0} in this limit. For each path, we have a rank-11 matrix C^jν\hat{C}_{j}^{\nu} associated with its endpoint. We will now argue that for each jj, the sequence of sets of matrices {C^jν}ν\{\hat{C}_{j}^{\nu}\}_{\nu} has a subsequence with a limit {C^j}\{\hat{C}_{j}\} the collection of which is also a rank-11 decomposition of the identity. This argument is somewhat complicated by the fact that κν\kappa^{\nu} can vary as ν\nu increases.

For each jj, we apply the Arzelà-Ascoli theorem [22] (see Ref. [18] for details applicable to the present context) to obtain a limiting path, but we make these limits for different jj consistent with each other as follows. For j=1j=1 there exists a subsequence, S1={ν11,ν12,}S_{1}=\{\nu_{11},\nu_{12},\ldots\}, of paths Π1ν1j(s)\Pi_{1}^{\nu_{1j}}(s) for which, by Arzelá-Ascoli, a limiting path of product operators exists, denote it Π1(s)\Pi_{1}(s). Since the endpoint of each path in this subsequence of paths is associated with a rank-11 matrix C^1ν1j\hat{C}_{1}^{\nu_{1j}}, then by the Eckart-Young-Mirsky theorem [21], the endpoint of the limiting path must also be represented by a rank-11 matrix C^1\hat{C}_{1}. Next, for j=2j=2, consider only those paths Π2ν(s)\Pi_{2}^{\nu}(s) and matrices C^2ν\hat{C}_{2}^{\nu} such that νS1\nu\in S_{1}. Here, again by Arzelá-Ascoli, there exists a subsequence S2={ν21,ν22,}S_{2}=\{\nu_{21},\nu_{22},\ldots\} of S1S_{1}, for which the corresponding set of paths has a limiting path of product operators, Π2(s)\Pi_{2}(s), and by the same arguments used above for S1S_{1}, the endpoint of this limiting path is again represented by a rank-11 matrix C^2\hat{C}_{2}. Since S2S_{2} is a subsequence of S1S_{1}, we have that limν,νS2Π1ν(s)=Π1(s)\lim_{\nu\to\infty,\nu\in S_{2}}\Pi_{1}^{\nu}(s)=\Pi_{1}(s), and limν,νS2C^1ν=C^1\lim_{\nu\to\infty,\nu\in S_{2}}\hat{C}_{1}^{\nu}=\hat{C}_{1}; these limits remain the same. Continuing this process for each jj, we find a subsequence SjS_{j} of Sj1S_{j-1} along with a limiting path Πj(s)\Pi_{j}(s) and rank-11 matrix C^j\hat{C}_{j}. Furthermore, the limiting paths and matrices obtained with subsequence SjS_{j^{\prime}} continue to be valid limits when restricted to SjS_{j}, for all j<jj^{\prime}<j. Since for each ν\nu, jC^jν=Iκν\sum_{j}\hat{C}_{j}^{\nu}=I_{\kappa^{\nu}}, one might naïvely expect that the set of limiting matrices C^j\hat{C}_{j} satisfies a similar condition, jC^j=Iκ\sum_{j}\hat{C}_{j}=I_{\kappa}. This conclusion is (essentially) correct, but let us make the argument with care.

As discussed above, rank-11 matrices C^jν\hat{C}_{j}^{\nu} obey the condition, limνC^jν=C^j\lim_{\nu\to\infty}\hat{C}_{j}^{\nu}=\hat{C}_{j}, and their sum behaves as limνjC^jν=limνIκν\lim_{\nu\to\infty}\sum_{j}\hat{C}_{j}^{\nu}=\lim_{\nu\to\infty}I_{\kappa^{\nu}}. And what is the limit of a sequence of identity operators? Note that IκIκνI_{\kappa}-I_{\kappa^{\nu}} can never be small unless κν=κ\kappa^{\nu}=\kappa. From the proof of Lemma 4, we know that κν>κ\kappa^{\nu}>\kappa is possible, so this equality may not hold. Nonetheless, we must recognize that this is not the relevant limit for discussing convergence of a sequence of quantum channels implemented by LOCC protocols. Rather, as we have shown above, what is relevant are the operators, Ejν=mm(C^jν)mmK^mνK^mνE_{j}^{\nu}=\sum_{mm^{\prime}}(\hat{C}_{j}^{\nu})_{mm^{\prime}}\hat{K}_{m}^{\nu{\dagger}}\hat{K}_{m^{\prime}}^{\nu}. Now, jEjν=mm(Iκν)mmK^mνK^mν=mK^mνK^mν=I\sum_{j}E_{j}^{\nu}=\sum_{mm^{\prime}}(I_{\kappa^{\nu}})_{mm^{\prime}}\hat{K}_{m}^{\nu{\dagger}}\hat{K}_{m^{\prime}}^{\nu}=\sum_{m}\hat{K}_{m}^{\nu{\dagger}}\hat{K}_{m}^{\nu}=I_{\cal H}, as required. It may happen that as ν\nu becomes large, the norm of K^mν\hat{K}_{m}^{\nu} may become negligible, tending asymptotically to zero. In that case, even though jC^jν=Iκν\sum_{j}\hat{C}_{j}^{\nu}=I_{\kappa^{\nu}}, the only non-negligible parts of the sum for EjνE_{j}^{\nu} (see a few sentences previous to this one) correspond to matrix elements (C^jν)mm(\hat{C}_{j}^{\nu})_{mm^{\prime}} for which K^mν,K^mν\left\|\hat{K}_{m}^{\nu}\right\|,\left\|\hat{K}_{m^{\prime}}^{\nu}\right\| do not tend to zero. Under these circumstances, while jC^jν=Iκν\sum_{j}\hat{C}_{j}^{\nu}=I_{\kappa^{\nu}}, it is only the part of IκνI_{\kappa^{\nu}} lying in that subspace corresponding to the non-negligible K^mν\hat{K}_{m}^{\nu} that is relevant. Since by assumption the limit of this sequence of channels exists and is equal to {\cal E} of Kraus rank κ\kappa, then according to these arguments that relevant subspace must be of size κ\kappa. By Lemma 4, we have κνκ\kappa^{\nu}\geq\kappa for large enough ν\nu, which guarantees the existence of a subspace of size κ\kappa within the κν\kappa^{\nu}-dimensional space upon which IκνI_{\kappa^{\nu}} acts, for large enough ν\nu. Therefore, IκνI_{\kappa^{\nu}} may (essentially) be replaced by the projector onto this subspace, which is (essentially) equivalent to IκI_{\kappa}. With this, the proof of Theorem 1 is complete.\blacksquare

V Detailed example of a quantum channel

In this section, we consider the example [11] of a two-qubit state transformation which can be arbitrarily closely approximated, but cannot be implemented exactly by LOCC. In Section VIII, we consider a generalization of this example to any number of qubits. Consider a quantum instrument [11], which consists of the following three CP maps, the sum of which form the quantum channel {\cal E} we consider here, but are not individually trace preserving,

𝒥1(ρ)\displaystyle{\cal J}_{1}(\rho) =[11]ρ[11]\displaystyle=[11]\rho[11]
𝒥2(ρ)\displaystyle{\cal J}_{2}(\rho) =i=12(Ti[0])ρ(Ti[0])\displaystyle=\sum_{i=1}^{2}(T_{i}\otimes[0])\rho(T_{i}^{\dagger}\otimes[0])
𝒥3(ρ)\displaystyle{\cal J}_{3}(\rho) =i=12([0]Ti)ρ([0]Ti),\displaystyle=\sum_{i=1}^{2}([0]\otimes T_{i})\rho([0]\otimes T_{i}^{\dagger}), (13)

where T1=I/3T_{1}=I/\sqrt{3} and T2=([0]+2[1])/6T_{2}=\left([0]+2[1]\right)/\sqrt{6}. This suggests describing channel {\cal E} by the five product Kraus operators,

K1\displaystyle K_{1}^{\prime} =[11]\displaystyle=[11]
K2\displaystyle K_{2}^{\prime} =13([00]+[10])\displaystyle=\sqrt{\frac{1}{3}}\left([00]+[10]\right)
K3\displaystyle K_{3}^{\prime} =16([00]+2[10])\displaystyle=\sqrt{\frac{1}{6}}\left([00]+2[10]\right)
K4\displaystyle K_{4}^{\prime} =13([00]+[01])\displaystyle=\sqrt{\frac{1}{3}}\left([00]+[01]\right)
K5\displaystyle K_{5}^{\prime} =16([00]+2[01]).\displaystyle=\sqrt{\frac{1}{6}}\left([00]+2[01]\right). (14)

It turns out that there is a smaller set of (product) Kraus operators that represent the same channel. For the sake of simplicity, we will use this smaller set, which is

K^1\displaystyle\hat{K}_{1} =[11]\displaystyle=[11]
K^2\displaystyle\hat{K}_{2} =23[00]+[10]\displaystyle=\frac{2}{3}[00]+[10]
K^3\displaystyle\hat{K}_{3} =23[00]+[01]\displaystyle=\frac{2}{3}[00]+[01]
K^4\displaystyle\hat{K}_{4} =13[00].\displaystyle=\frac{1}{3}[00]. (15)

It is easily checked that Kj=jWjlK^lK_{j}^{\prime}=\sum_{j}W_{jl}\hat{K}_{l}, with WW the isometry,

W=(1000013013023016001313002316).\displaystyle W=\begin{pmatrix}1&0&0&0\\ 0&\sqrt{\frac{1}{3}}&0&\sqrt{\frac{1}{3}}\\ 0&\sqrt{\frac{2}{3}}&0&-\sqrt{\frac{1}{6}}\\ 0&0&\sqrt{\frac{1}{3}}&\sqrt{\frac{1}{3}}\\ 0&0&\sqrt{\frac{2}{3}}&-\sqrt{\frac{1}{6}}\end{pmatrix}. (16)

The following LOCC protocol [23, 11], 𝒫ν{\cal P}^{\nu}, implements a quantum channel, ν{\cal E}^{\nu}, which approximates the channel {\cal E} of Eq. (V). Given single qubit (Hermitian) operators, A0=1ϵ[0]A+[1]AA_{0}=\sqrt{1-\epsilon}[0]_{A}+[1]_{A} and A1=ϵ[0]AA_{1}=\sqrt{\epsilon}[0]_{A}, Alice begins by measuring {A0,A1}\{A_{0},A_{1}\} on her system AA, with corresponding POVM elements 𝒜j=AjAj=Aj2,j=0,1{\cal A}_{j}=A_{j}^{\dagger}A_{j}=A_{j}^{2},~{}j=0,1. She shares her outcome with Bob, who follows with the same measurement—having POVM elements 0=(1ϵ)[0]B+[1]B{\cal B}_{0}=(1-\epsilon)[0]_{B}+[1]_{B} and 1=ϵ[1]B{\cal B}_{1}=\epsilon[1]_{B}—on his system BB, sharing his outcome in turn with Alice. They halt at this point unless they both obtained outcome 0, in which case they repeat their local measurements as just described. Then, after a maximum of ν\nu iterations of this procedure, they stop. Noting that A0n=(1ϵ)n/2[0]A+[1]AA_{0}^{n}=(1-\epsilon)^{n/2}[0]_{A}+[1]_{A} and A1A0n=ϵ1/2(1ϵ)n/2[0]AA_{1}A_{0}^{n}=\epsilon^{1/2}(1-\epsilon)^{n/2}[0]_{A}, and similarly for Bob’s measurement operators, then the first few rounds of this LOCC protocol—represented as a tree graph with nodes labeled by the corresponding POVM elements implemented to that point in the protocol—is depicted in Figure 1 (note that the notation is condensed slightly in this figure by defining η=1ϵ\eta=1-\epsilon). In the figure, there is a single main branch, stretching down along the upper right side, which continues (note ellipses at bottom right) for 2ν2\nu segments, while all other branches extend either one or two segments off of that main branch and then terminate. As discussed in [11], the limiting channel is {\cal E}, as desired, as long as we choose ϵ=νc\epsilon=\nu^{-c} with 0<c<10<c<1, which we will do in what follows.

Refer to caption
Figure 1: Tree graph for the first several local measurements performed by Alice and Bob in any one of the sequence of LOCC protocols described in the text. Note that we have defined η=1ϵ\eta=1-\epsilon.

Consider that main branch in each protocol and the sequence of these main branches formed as we increment ν\nu. Each of these main branches describes a path, Πν(s)\Pi^{\nu}(s), of product operators in operator space, as described in Sections III and IV. The sequence of such paths of product operators has a limiting path, which we will now identify. Defining s=Tr(𝒜0n0n)=[(1ϵ)n+1]2s=\textrm{Tr}\left({\cal A}_{0}^{n}\otimes{\cal B}_{0}^{n}\right)=[(1-\epsilon)^{n}+1]^{2}, we have s1=(1ϵ)n\sqrt{s}-1=(1-\epsilon)^{n}, at least at those particular nodes along this main branch. This suggests the quantity s1\sqrt{s}-1 may play a role. Noting that we will be considering the limit as ϵ0\epsilon\to 0, and that in this limit, we will also have the length of the branch, ν\nu\to\infty, care must be taken in considering how the quantity (1ϵ)n(1-\epsilon)^{n} behaves in these limits (since nn ranges up to ν\nu). Nonetheless, the quantity s1\sqrt{s}-1 is perfectly well-behaved in these limits, and in fact, we can prove the following Lemma.

Lemma 5.

The limiting path of product operators obtained from the sequence of paths, Πν(s)\Pi^{\nu}(s), is ΠM(s)=limνΠν(s)=M(s)M(s)\Pi_{M}(s)=\lim_{\nu\to\infty}\Pi^{\nu}(s)=M(s)\otimes M(s), 4s14\geq s\geq 1, with M(s)=(s1)[0]+[1]M(s)=\left(\sqrt{s}-1\right)[0]+[1].

Note that Tr(ΠM(s))=s\textrm{Tr}\left(\Pi_{M}(s)\right)=s, so the path is parametrized by the trace, as has been done in [18].

Proof.

Define tn=(1ϵ)n+1t_{n}=(1-\epsilon)^{n}+1, s2n=Tr(𝒜0n0n)=tn2s_{2n}=\textrm{Tr}\left({\cal A}_{0}^{n}\otimes{\cal B}_{0}^{n}\right)=t_{n}^{2} and s2n+1=Tr(𝒜0n+10n)=tntn+1s_{2n+1}=\textrm{Tr}\left({\cal A}_{0}^{n+1}\otimes{\cal B}_{0}^{n}\right)=t_{n}t_{n+1}. For each ν,ϵ\nu,\epsilon, and 0nν0\leq n\leq\nu, consider a piecewise local segment from 𝒜0n0n{\cal A}_{0}^{n}\otimes{\cal B}_{0}^{n} to 𝒜0n+10n{\cal A}_{0}^{n+1}\otimes{\cal B}_{0}^{n} along the main branch of the ν\nuth protocol, with s2nss2n+1s_{2n}\geq s\geq s_{2n+1} or tn2stntn+1t_{n}^{2}\geq s\geq t_{n}t_{n+1}. This segment along the path may be parametrized as

Πν(s)\displaystyle\Pi^{\nu}(s) =ss2n+1s2ns2n+1𝒜0n0n+s2nss2ns2n+1𝒜0n+10n\displaystyle=\frac{s-s_{2n+1}}{s_{2n}-s_{2n+1}}{\cal A}_{0}^{n}\otimes{\cal B}_{0}^{n}+\frac{s_{2n}-s}{s_{2n}-s_{2n+1}}{\cal A}_{0}^{n+1}\otimes{\cal B}_{0}^{n}
=1tn(tntn+1)[(stntn+1)(tn1001)+(tn2s)(tn+11001)](tn1001)\displaystyle=\frac{1}{t_{n}(t_{n}-t_{n+1})}\left[(s-t_{n}t_{n+1})\begin{pmatrix}t_{n}-1&0\\ 0&1\end{pmatrix}+(t_{n}^{2}-s)\begin{pmatrix}t_{n+1}-1&0\\ 0&1\end{pmatrix}\right]\otimes\begin{pmatrix}t_{n}-1&0\\ 0&1\end{pmatrix}
=1tn(stn00tn)(tn1001).\displaystyle=\frac{1}{t_{n}}\begin{pmatrix}s-t_{n}&0\\ 0&t_{n}\end{pmatrix}\otimes\begin{pmatrix}t_{n}-1&0\\ 0&1\end{pmatrix}. (17)

We need to show that for each ss and for all δ>0\delta>0, there exists μ\mu such that for all ν>μ\nu>\mu, δ>Πν(s)ΠM(s)1\delta>\left\|\Pi^{\nu}(s)-\Pi_{M}(s)\right\|_{1}. Consider

tn00|[Πν(s)ΠM(s)]|00\displaystyle t_{n}\langle 00|\left[\Pi^{\nu}(s)-\Pi_{M}(s)\right]|00\rangle =(stn)(tn1)tn(s1)2\displaystyle=(s-t_{n})(t_{n}-1)-t_{n}(\sqrt{s}-1)^{2}
=2tnstn2s\displaystyle=2t_{n}\sqrt{s}-t_{n}^{2}-s
2tns2ntn2s2n+1=tn2tntn+1=tn(tntn+1).\displaystyle\leq 2t_{n}\sqrt{s_{2n}}-t_{n}^{2}-s_{2n+1}=t_{n}^{2}-t_{n}t_{n+1}=t_{n}(t_{n}-t_{n+1}). (18)

Since tntn+1=ϵ(1ϵ)nt_{n}-t_{n+1}=\epsilon(1-\epsilon)^{n}, this is small. Similarly,

tn01|[ΠM(s)Πν(s)]|01\displaystyle t_{n}\langle 01|\left[\Pi_{M}(s)-\Pi^{\nu}(s)\right]|01\rangle =tn(s1)(stn)\displaystyle=t_{n}(\sqrt{s}-1)-(s-t_{n})
=tnss\displaystyle=t_{n}\sqrt{s}-s
tn(tntn+1),\displaystyle\leq t_{n}(t_{n}-t_{n+1}), (19)

and it is easy to see that 10|[Πν(s)ΠM(s)]|10=tnstntn+1\langle 10|\left[\Pi^{\nu}(s)-\Pi_{M}(s)\right]|10\rangle=t_{n}-\sqrt{s}\leq t_{n}-t_{n+1}, as well. Since the 1111 diagonal element vanishes identically, we have that

Πν(s)ΠM(s)13(tntn+1)=3ϵ(1ϵ)n3ϵ=3νc.\displaystyle\left\|\Pi^{\nu}(s)-\Pi_{M}(s)\right\|_{1}\leq\sqrt{3}(t_{n}-t_{n+1})=\sqrt{3}\epsilon(1-\epsilon)^{n}\leq\sqrt{3}\epsilon=\sqrt{3}\nu^{-c}. (20)

Considering the segments from 𝒜0n+10n{\cal A}_{0}^{n+1}\otimes{\cal B}_{0}^{n} to 𝒜0n+10n+1{\cal A}_{0}^{n+1}\otimes{\cal B}_{0}^{n+1}, corresponding to s2n+1ss2n+2s_{2n+1}\geq s\geq s_{2n+2}, one may follow the same argument as above—the only difference being that tn+1tn+2t_{n+1}\to t_{n+2} and tntn+1t_{n}\to t_{n+1}—yielding the exact same result. Since this encompasses all cases, then for the entire path along that main branch, we may choose μ=(δ/3)1/c\mu=\left(\delta/\sqrt{3}\right)^{-1/c} yielding Πν(s)ΠM(s)1<δ\left\|\Pi^{\nu}(s)-\Pi_{M}(s)\right\|_{1}<\delta for all ν>μ\nu>\mu and for all ss, and this ends the proof.\blacksquare

It is instructive to understand what happens with the other outcomes of our sequence of protocols. Recalling that s=[(1ϵ)n+1]2s=[(1-\epsilon)^{n}+1]^{2} and that ϵ\epsilon decreases as ν\nu increases, the value of ss at which the nnth side-branch emerges from the main branch increases with ν\nu. One may envision this as these side-branches sliding up the main branch toward the root node in Fig. 1 as ν\nu increases, as more and more of these side-branches are added further down the main branch. In the limit ν\nu\to\infty, the gap between each side-branch and its nearest neighbors vanishes, and we have a continuous set of complete separable (generally, not LOCC222Note that since the limiting main branch is not piecewise local—that is, each point along this branch differs from the next in both party’s local operators, not one party’s at a time—these measurements at each ss are not LOCC.) measurements being made, one at each ss, with one outcome being ΠM(sΔs)\Pi_{M}(s-\Delta s) and the other outcomes we will now determine. By taking the limit as Δs0\Delta s\to 0, we see that it is simply neccessary to differentiate the path along the main branch. This gives

dΠM(s)=ds2s([0]M(s)+M(s)[0]),\displaystyle d\Pi_{M}(s)=\frac{ds}{2\sqrt{s}}\left([0]\otimes M(s)+M(s)\otimes[0]\right), (21)

which is not a product operator, as it must be for the limiting channel to be in LOCC¯\overline{\textrm{LOCC}}. However, we can simply split it into two parts, each of which is clearly a product operator, providing two additional outcomes for each of the continuous set of measurements appearing at the various values of ss. Notice the path from ΠM(s)\Pi_{M}(s) to either one of these parts is, in fact, piecewise local, extending from M(s)M(s)M(s)\otimes M(s) to either ds[0]M(s)d\sqrt{s}[0]\otimes M(s) or dsM(s)[0]d\sqrt{s}M(s)\otimes[0]. One of these outcomes clearly corresponds to those of the form ϵ(1ϵ)n[0]0n\epsilon(1-\epsilon)^{n}[0]\otimes{\cal B}_{0}^{n} in each of our sequence of LOCC protocols, and the other to ϵ(1ϵ)n𝒜0n[0]\epsilon(1-\epsilon)^{n}{\cal A}_{0}^{n}\otimes[0]. One may conjecture that the finite-round parts of (convergent) sequences of LOCC protocols always converge to piecewise local paths, as we have now seen happens for this specific example, which also demonstrates how the parts whose number of rounds tends to infinity can converge to paths that are not piecewise local. Note that in the limit ϵ0\epsilon\to 0, there are no paths to outcomes proportional to [00][00], consistent with the observation (see Fig. 1) that in each of the sequence of protocols, these outcomes are proportional to ϵ2\epsilon^{2} and so become negligible. On the other hand, outcomes proportional to ϵ\epsilon are not negligible, since their number diverges as ν=ϵ1/c\nu=\epsilon^{-1/c}, and taken together, they provide a significant contribution. To summarize then, the POVM implemented by the protocol illustrated in Fig. 1 in the limit as ν\nu\to\infty (ϵ0\epsilon\to 0) consists of the following operators,

E1\displaystyle E_{1} =[11]\displaystyle=[11]
E2(s)\displaystyle E_{2}(s) =dsM(s)[0],4s1\displaystyle=d\sqrt{s}M(s)\otimes[0],~{}4\geq s\geq 1
E3(s)\displaystyle E_{3}(s) =ds[0]M(s),4s1.\displaystyle=d\sqrt{s}[0]\otimes M(s),~{}4\geq s\geq 1. (22)

One immediately obtains from this the limiting generalized measurement operators, for example K2(s)=E2(s)K_{2}(s)=\sqrt{E_{2}(s)}, when viewing the sequence of protocols as implementing generalized measurements. That is,

K1\displaystyle K_{1} =[11]\displaystyle=[11]
K2(s)\displaystyle K_{2}(s) =dsM(s)[0],4s1\displaystyle=\sqrt{d\sqrt{s}M(s)\otimes[0]},~{}4\geq s\geq 1
K3(s)\displaystyle K_{3}(s) =ds[0]M(s),4s1.\displaystyle=\sqrt{d\sqrt{s}[0]\otimes M(s)},~{}4\geq s\geq 1. (23)

We now calculate the associated Choi-Jamiolkowski operator to demonstrate that this infinite set of limiting Kraus operators represents the desired channel {\cal E}, see Eq. (V). There is one limiting Kraus operator obtained at the end of the main branch in the limit, which is M(1)M(1)=[11]=K^1\sqrt{M(1)\otimes M(1)}=[11]=\hat{K}_{1}. In addition, within the range 4s14\geq s\geq 1, we have two infinite sets of Kraus operators, K2(s)K_{2}(s) and K3(s)K_{3}(s).

Define σ=s1\sigma=\sqrt{s}-1 so that ds=dσd\sqrt{s}=d\sigma, and let the Choi-Jamiolkowski operator be written as Ω=Ω1+Ω2+Ω3\Omega=\Omega_{1}+\Omega_{2}+\Omega_{3}, with

Ω1=w,x,y,z=01|wxyz|ABK^1|wxyz|ABK^1=[1111],\displaystyle\Omega_{1}=\sum_{w,x,y,z=0}^{1}|wx\rangle\langle yz|^{A^{\prime}B^{\prime}}\otimes\hat{K}_{1}|wx\rangle\langle yz|^{AB}\hat{K}_{1}=[1111], (24)
Ω2\displaystyle\Omega_{2} =w,x,y,z=01|wxyz|AB12𝑑s(M(s)[0])|wxyz|AB(M(s)[0])\displaystyle=\sum_{w,x,y,z=0}^{1}|wx\rangle\langle yz|^{A^{\prime}B^{\prime}}\otimes\int_{1}^{2}d\sqrt{s}\left(\sqrt{M(s)}\otimes[0]\right)|wx\rangle\langle yz|^{AB}\left(\sqrt{M(s)}\otimes[0]\right)
=12𝑑σ(σ|0000+|1010)(σ0000|+1010|)\displaystyle=\int_{1}^{2}d\sigma\left(\sqrt{\sigma}|0000\rangle+|1010\rangle\right)\left(\sqrt{\sigma}\langle 0000|+\langle 1010|\right)
=12[0000]+23(|00001010|+|10100000|)+[1010],\displaystyle=\frac{1}{2}[0000]+\frac{2}{3}\left(|0000\rangle\langle 1010|+|1010\rangle\langle 0000|\right)+[1010], (25)

and in entirely similar fashion, we find that

Ω3\displaystyle\Omega_{3} =12[0000]+23(|00000101|+|01010000|)+[0101].\displaystyle=\frac{1}{2}[0000]+\frac{2}{3}\left(|0000\rangle\langle 0101|+|0101\rangle\langle 0000|\right)+[0101]. (26)

Adding these together yields

Ω\displaystyle\Omega =I+23(|00000101|+|01010000|+|00001010|+|10100000|)`,\displaystyle=I_{{\cal H}\otimes{\cal H}}+\frac{2}{3}\left(|0000\rangle\langle 0101|+|0101\rangle\langle 0000|+|0000\rangle\langle 1010|+|1010\rangle\langle 0000|\right)`, (27)

in agreement with the results of [11] for the Choi-Jamiolkowski operator of {\cal E}.333There seem to be typographical errors in [11] in their calculation of this Choi-Jamiolkowski operator—we believe their AAA^{\prime}A matrices should be 4×44\times 4 instead of 2×22\times 2—but their intended meaning is clear.

Another way to demonstrate that these limiting Kraus operators represent the desired channel {\cal E} is by showing that they are related to the set in Eq. (V) by an isometry. In this case, the isometry will have an infinite number of rows. We have that one of these Kraus operators is K^1\hat{K}_{1}, so that the first row of our isometry is [1000][1~{}0~{}0~{}0]. The other rows are obtained from K2(s)K_{2}(s) and K3(s)K_{3}(s), written as

K2(s)\displaystyle K_{2}(s) =dσM(s)[0]=dσ(K^2+(3σ2)K^4),\displaystyle=\sqrt{d\sigma M(s)\otimes[0]}=\sqrt{d\sigma}\left(\hat{K}_{2}+(3\sqrt{\sigma}-2)\hat{K}_{4}\right),
K3(s)\displaystyle K_{3}(s) =dσ[0]M(s)=dσ(K^3+(3σ2)K^4).\displaystyle=\sqrt{d\sigma[0]\otimes M(s)}=\sqrt{d\sigma}\left(\hat{K}_{3}+(3\sqrt{\sigma}-2)\hat{K}_{4}\right). (28)

That is, we obtain a continuous set of additional rows from each of these: dσ[010(3σ2)]\sqrt{d\sigma}[0~{}1~{}0~{}(3\sqrt{\sigma}-2)] and dσ[001(3σ2)]\sqrt{d\sigma}[0~{}0~{}1~{}(3\sqrt{\sigma}-2)]. As an isometry, all columns must be orthonormal. The first column is [1000]T[1~{}0~{}0~{}\ldots~{}0]^{T}, so this column is normalized as well as orthogonal to the other columns, since each of their first entries is zero. Normalization of the second and third columns follows from the fact that 12𝑑σ=1\int_{1}^{2}d\sigma=1. Orthogonality of this pair of columns is clear, since the non-zero entries in one are always zero in the other. Normalization of the last column is

212𝑑σ(3σ2)2\displaystyle 2\int_{1}^{2}d\sigma(3\sqrt{\sigma}-2)^{2} =201𝑑σ(9σ12σ+4)=1,\displaystyle=2\int_{0}^{1}d\sigma(9\sigma-12\sqrt{\sigma}+4)=1, (29)

where the factor of 22 is included because there are two entries in this column for each value of ss, one from K2(s)K_{2}(s) and the other from K3(s)K_{3}(s), see Eq. (V). Thus, we see that all columns are normalized. Orthogonality of the last column to either of the middle two is

01𝑑σ(3σ2)=0\displaystyle\int_{0}^{1}d\sigma(3\sqrt{\sigma}-2)=0 (30)

Thus, we do indeed see that this continuous set of Kraus operators is related to those in Eq. (V) by an isometry.

As noted in Item 2 of Theorem 2, the endpoints of our limiting paths should be described by positive semidefinite matrices which have rank equal to unity. The current example illustrates this aspect, as well. We have that dσM(s)[0]=mnC2(s)mnK^mK^nd\sigma M(s)\otimes[0]=\sum_{mn}C_{2}(s)_{mn}\hat{K}_{m}^{\dagger}\hat{K}_{n}, with the K^m\hat{K}_{m} given in Eq. (V) and

C2(s)=dσ(00000103σ2000003σ20(3σ2)2).\displaystyle C_{2}(s)=d\sigma\begin{pmatrix}0&0&0&0\\ 0&1&0&3\sqrt{\sigma}-2\\ 0&0&0&0\\ 0&3\sqrt{\sigma}-2&0&\left(3\sqrt{\sigma}-2\right)^{2}\end{pmatrix}. (31)

Similarly, dσ[0]M(s)=mnC3(s)mnK^mK^nd\sigma[0]\otimes M(s)=\sum_{mn}C_{3}(s)_{mn}\hat{K}_{m}^{\dagger}\hat{K}_{n}, with

C3(s)=dσ(000000000013σ2003σ2(3σ2)2).\displaystyle C_{3}(s)=d\sigma\begin{pmatrix}0&0&0&0\\ 0&0&0&0\\ 0&0&1&3\sqrt{\sigma}-2\\ 0&0&3\sqrt{\sigma}-2&\left(3\sqrt{\sigma}-2\right)^{2}\end{pmatrix}. (32)

Both C2(s)C_{2}(s) and C3(s)C_{3}(s) are rank-11 matrices for every ss. They also provide an easy way to find a matrix C1(s)C_{1}(s) from which ΠM(s)=mnC1(s)mnK^mK^n\Pi_{M}(s)=\sum_{mn}C_{1}(s)_{mn}\hat{K}_{m}^{\dagger}\hat{K}_{n}, the path along the (limit of the) main branch in Fig. 1, may be obtained. That is, since the collection of these matrices must sum to IκI_{\kappa}, we have

C1(s)\displaystyle C_{1}(s) =I4s2(C2(s)+C3(s))\displaystyle=I_{4}-\int_{\sqrt{s}}^{2}\left(C_{2}(s)+C_{3}(s)\right)
=(1000000000000000)+σ(00000102(σ1)0012(σ1)02(σ1)2(σ1)9s116σ).\displaystyle=\begin{pmatrix}1&0&0&0\\ 0&0&0&0\\ 0&0&0&0\\ 0&0&0&0\end{pmatrix}+\sigma\begin{pmatrix}0&0&0&0\\ 0&1&0&2\left(\sqrt{\sigma}-1\right)\\ 0&0&1&2\left(\sqrt{\sigma}-1\right)\\ 0&2\left(\sqrt{\sigma}-1\right)&2\left(\sqrt{\sigma}-1\right)&9\sqrt{s}-1-16\sqrt{\sigma}\end{pmatrix}. (33)

This is equal to the identity matrix at σ=1\sigma=1, of course, and since the second term in the second line of Eq. (V) vanishes at σ=0\sigma=0, it leads to the endpoint of the limiting main branch being described by ΠM(1)=K^1K^1\Pi_{M}(1)=\hat{K}_{1}^{\dagger}\hat{K}_{1}, with rank(C1(1))=1\textrm{rank}\left(C_{1}(1)\right)=1, as required. Furthermore, one calculates using the set of operators in Eq. (V) that mn(C1(s))mnK^mK^n=M(s)M(s)\sum_{mn}(C_{1}(s))_{mn}\hat{K}_{m}^{\dagger}\hat{K}_{n}=M(s)\otimes M(s) for all ss, as it must.

Note that C2(s)C_{2}(s) and C3(s)C_{3}(s) correspond to the endpoints of side branches that begin at position ss of the main branch. It is easily seen that the entire length of these side branches may be represented by matrices C2(s,x)=(1x)C1(s)+xC2(s)C_{2}(s,x)=(1-x)C_{1}(s)+xC_{2}(s) and C3(s,x)=(1x)C1(s)+xC3(x)C_{3}(s,x)=(1-x)C_{1}(s)+xC_{3}(x), 0x10\leq x\leq 1. Since these are convex combinations of positive semidefinite matrices, they are themselves positive semidefinite for all ss and xx.

According to the definition of ZZ_{\cal E} in Item 3 of Theorem 1, the fact that C1(s)C_{1}(s) and Cj(s,x),j=2,3,C_{j}(s,x),~{}j=2,3, are all positive semidefinite matrices shows that the paths associated with all these limiting measurement outcomes lie entirely within ZZ_{\cal E}. It is not difficult to see that most of the length of the paths associated with the main branch of the individual LOCC protocols in our approximating sequence do not lie in ZZ_{\cal E}, but only approach it in the limit. While the endpoints of those approximating side branches do lie within ZZ_{\cal E}, they begin at the main branch so are at least partly outside of ZZ_{\cal E}, as well.

VI Quantum instruments

Generalized measurements and quantum channels share similar characteristics, both involving transformations of state ρ\rho to KjρKjK_{j}\rho K_{j}^{\dagger}, with the channel involving an unrestricted sum over the latter quantities, there being no accessible classical information about the outcome jj. Both are special cases of the more general notion of a a quantum instrument, which covers circumstances where some classical information is available but possibly less than is required for a full generalized measurement. The proof of Theorem 2 for channels in Section IV is also capable of encompassing all quantum instruments as long as we make allowances for this required classical information, which can be embedded in the form of the positive semidefinite CC-matrices introduced in Theorem 2, as we now discuss.

In the definition we have given here for zonoid ZZ_{\cal E}, sums over terms of the form C^mmK^mK^m\hat{C}_{mm^{\prime}}\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}} appear. On the other hand, in our analysis of quantum measurements [18, 20], the zonotopes used there involved simple, positive linear combinations of KmKmK_{m}^{\dagger}K_{m}. In other words, for measurements, the CC-matrices are diagonal, C^mm=cmδmm\hat{C}_{mm^{\prime}}=c_{m}\delta_{mm^{\prime}}. We will now argue that a simple generalization to block diagonal matrices will suffice for all quantum instruments.

Recall that to account for all possible Kraus representations of {\cal E}, we started with a positive linear combination, z=jcjKjKjz=\sum_{j}c_{j}K_{j}^{\dagger}K_{j} with {Kj}\{K_{j}\} constituting an arbitrary such representation, and then utilized isometric freedom, Kj=mWjmK^mK_{j}=\sum_{m}W_{jm}\hat{K}_{m}, to obtain z=mmC^mmK^mK^mz=\sum_{mm^{\prime}}\hat{C}_{mm^{\prime}}\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}}. Suppose we have a quantum instrument {\cal E} consisting of CP maps, 𝒥r{\cal J}_{r}, for some set of indices, rr. To understand what happens to isometric freedom in this case, consider as illustration instrument ={𝒥1,𝒥2}{\cal E}=\{{\cal J}_{1},{\cal J}_{2}\}, with 𝒥1(ρ)=K1ρK1+K2ρK2{\cal J}_{1}(\rho)=K_{1}\rho K_{1}^{\dagger}+K_{2}\rho K_{2}^{\dagger} and 𝒥2(ρ)=K3ρK3+K4ρK4{\cal J}_{2}(\rho)=K_{3}\rho K_{3}^{\dagger}+K_{4}\rho K_{4}^{\dagger}. If an LOCC protocol has a final outcome K=μ1K1+μ3K3K=\mu_{1}K_{1}+\mu_{3}K_{3}, this will mix the two CP maps, in general. This is so unless K1=αK3+βK4K_{1}=\alpha K_{3}+\beta K_{4} (or K3=αK1+βK2K_{3}=\alpha^{\prime}K_{1}+\beta^{\prime}K_{2}), in which case KK can be written as a linear combination of only those Kraus operators involved in expressing 𝒥2{\cal J}_{2} (or 𝒥1{\cal J}_{1}), and then the parties can obtain enough classical information to identify KK as contributing to one of the CP maps and not the other. If such a separation between the CP maps is impossible for one or more outcomes of the protocol, then there is no way to coarse-grain the collection of outcomes such that each subset exactly represents one of those CP maps. Thus for a successful implementation, it must be possible to express each of the final outcomes in a way that does not mix the different CP maps of the desired instrument {\cal E}. As the above example illustrates, any linear dependence of the full set of Kraus operators allows for multiple ways of expressing any given outcome, but nonetheless, there must be a choice of these expressions which keeps the different CP maps of {\cal E} separate.

Thus, we see that the full isometric freedom of quantum channels is no longer in effect in this case, but it is still allowed within each 𝒥r{\cal J}_{r}. That is, for example, if ={𝒥r}{\cal E}=\{{\cal J}_{r}\} and 𝒥r(ρ)=mSrK^mρK^m{\cal J}_{r}(\rho)=\sum_{m\in S_{r}}\hat{K}_{m}\rho\hat{K}_{m}, with {K^m}mSr\{\hat{K}_{m}\}_{m\in S_{r}} a minimal set of Kraus operators for 𝒥r{\cal J}_{r}, we must be able to write any outcome KjK_{j} in terms of one, and only one, of these sets. That is for some fixed rr, Kj=mSr(Wr)jmK^mK_{j}=\sum_{m\in S_{r}}(W_{r})_{jm}\hat{K}_{m} with WrWr=IκrW_{r}^{\dagger}W_{r}=I_{\kappa_{r}} the κr×κr\kappa_{r}\times\kappa_{r} identity matrix and κr\kappa_{r} is the Kraus rank of 𝒥r{\cal J}_{r}. Note that this leads to KjKj=m,mSrC^mmK^mK^mK_{j}^{\dagger}K_{j}=\sum_{m,m^{\prime}\in S_{r}}\hat{C}_{mm^{\prime}}\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}} with C^mm=(Wr)jm(Wr)jm\hat{C}_{mm^{\prime}}=(W_{r})_{jm}^{\ast}(W_{r})_{jm^{\prime}}, so that CC is a rank-11 matrix, as required by Theorem 2. In addition, since these outcomes are to be coarse-grained (collected together) to represent 𝒥r{\cal J}_{r}, it must be that

jSrKjρKj=mSrK^mρK^m=𝒥r(ρ)\displaystyle\sum_{j\in S_{r}^{\prime}}K_{j}\rho K_{j}=\sum_{m\in S_{r}}\hat{K}_{m}\rho\hat{K}_{m}={\cal J}_{r}(\rho) (34)

for some fixed rr and index set SrS_{r}^{\prime}. This equality follows directly from our identification of WrW_{r} as an isometry acting on a κr\kappa_{r}-dimensional space.

Defining the direct sum, W=rWrW=\oplus_{r}W_{r}, we have that WW is a block-diagonal isometry mapping the collection of minimal (within each 𝒥r{\cal J}_{r}) Kraus sets to a complete set of Kraus operators {Kj}\{K_{j}\} for the full protocol implementing instrument {\cal E}. Since by Lemma 1 we have that each node in a finite-round LOCC protocol is a positive linear combination of the KjKjK_{j}^{\dagger}K_{j}, then each of these nodes—and therefore, every point along our piecewise local paths identified from these protocols—can be written as a positive semidefinite combination of the K^mK^m\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}} with the corresponding CC-matrix a block diagonal matrix, being itself a sum of block diagonal (rank-11) matrices. Thus, the proof of Theorem 2 holds for the case of quantum instruments, except that there is this additional constraint that the CC-matrices be block diagonal. These matrices keep the different CP maps separate from each other, precluding any mixing between the Kraus operators for 𝒥r{\cal J}_{r} and those for 𝒥r{\cal J}_{r^{\prime}} when rrr\neq r^{\prime}.

These intuitive ideas, which are formalized in Appendix B, lead directly to a generalization of Theorem 2, encompassing the most general notion of quantum instruments. Note that to generalize the proof of Lemma 3 for this case, see Appendix B for details, we simply add this block diagonal constraint on isometries W,WνW,W^{\nu} introduced just before—and matrices C,CνC,C^{\nu} appearing in—Eq. (A), and the rest of the proof of this lemma goes through as before.

Theorem 3.

Given quantum instrument ={𝒥r}{\cal E}=\{{\cal J}_{r}\} acting on input space (){\cal B}({\cal H}), with each 𝒥r{\cal J}_{r} represented by the minimal set of κr\kappa_{r} Kraus operators, {K^m}mSr\{\hat{K}_{m}\}_{m\in S_{r}}, then if {\cal E} can be implemented by LOCC with vanishingly small error, the following conditions must hold:

  1. 1.

    There exists a set, {Πj(s)}\{\Pi_{j}(s)\}, of continuous, monotonic paths of positive semidefinite product operators, each of which begins at II_{\cal H}, the identity operator on {\cal H}, and ends at a positive semidefinite operator of the form, Ej=mm(C^j)mmK^mK^mE_{j}=\sum_{mm^{\prime}}(\hat{C}_{j})_{mm^{\prime}}\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}}.

  2. 2.

    Each C^j\hat{C}_{j} is a positive semidefinite matrix, C^j0\hat{C}_{j}\geq 0, and may be chosen to have rank equal to unity, with the collection satisfying jC^j=Iκ\sum_{j}\hat{C}_{j}=I_{\kappa}, where IκI_{\kappa} is the κ×κ\kappa\times\kappa identity matrix and κ=rκr\kappa=\sum_{r}\kappa_{r}.

  3. 3.

    Each of these paths of product operators lies entirely within the zonoid, ZZ_{{\cal E}}, consisting of all positive semidefinite linear combinations of the operators, K^mK^m\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}}: Z={z|z=mmC^mmK^mK^m,IκC0}Z_{{\cal E}}=\left\{z\left|z=\sum_{mm^{\prime}}\hat{C}_{mm^{\prime}}\hat{K}_{m}^{\dagger}\hat{K}_{m^{\prime}},I_{\kappa}\geq C\geq 0\right.\right\}, with 1rank(C)κ1\leq\textrm{rank}(C)\leq\kappa and CC is block diagonal, the rrth diagonal block having size κr×κr\kappa_{r}\times\kappa_{r}.

  4. 4.

    There exists a partition of the index set jj into subsets SrS_{r}^{\prime} such that jSrC^j=Pκr\sum_{j\in S_{r}^{\prime}}\hat{C}_{j}=P_{\kappa_{r}}, with PκrP_{\kappa_{r}} a κ×κ\kappa\times\kappa matrix having the κr×κr\kappa_{r}\times\kappa_{r} identity matrix embedded within its rrth block and everything else being zero. That is, every matrix element of PκrP_{\kappa_{r}} vanishes except that the diagonal elements within the rrth block are equal to unity.

To illustrate this theorem, we return to the example of the preceding section. The Kraus operators K1,K2(s),K3(s)K_{1},K_{2}(s),K_{3}(s) of Eq. (V) can be written in terms of those of Eq. (V) as

K1\displaystyle K_{1} =K1\displaystyle=K_{1}^{\prime}
K2(s)\displaystyle K_{2}(s) =dσM(s)[0]=3dσ[(2σ1)K2+2(1σ)K3],\displaystyle=\sqrt{d\sigma M(s)\otimes[0]}=\sqrt{3d\sigma}\left[\left(2\sqrt{\sigma}-1\right)K_{2}^{\prime}+\sqrt{2}\left(1-\sqrt{\sigma}\right)K_{3}^{\prime}\right],
K3(s)\displaystyle K_{3}(s) =dσ[0]M(s)=3dσ[(2σ1)K4+2(1σ)K5].\displaystyle=\sqrt{d\sigma[0]\otimes M(s)}=\sqrt{3d\sigma}\left[\left(2\sqrt{\sigma}-1\right)K_{4}^{\prime}+\sqrt{2}\left(1-\sqrt{\sigma}\right)K_{5}^{\prime}\right]. (35)

We see that these Kraus operators, obtained from the limit of our sequence of protocols, are indeed related to those of the quantum instrument, see Eq. (V), by a block diagonal isometry. The first diagonal block has size 1×11\times 1, whereas the next two blocks are each infinite in length and two columns wide, and the five columns form an orthonormal set. This implies that the CC-matrices describing our paths are also block diagonal, as required by Theorem 3. From another point of view, we may coarse-grain all outcomes K2(s)K_{2}(s) together on the one hand, and all outcomes K3(s)K_{3}(s) on the other to find

K2(s)ρK2(s)\displaystyle\int K_{2}(s)\rho K_{2}(s) =12𝑑σM(s)[0]ρM(s)[0]=K2ρK2+K3ρK3=𝒥2(ρ)\displaystyle=\int_{1}^{2}d\sigma\sqrt{M(s)}\otimes[0]\rho\sqrt{M(s)}\otimes[0]=K_{2}^{\prime}\rho K_{2}^{\prime}+K_{3}^{\prime}\rho K_{3}^{\prime}={\cal J}_{2}(\rho)
K2(s)ρK3(s)\displaystyle\int K_{2}(s)\rho K_{3}(s) =12𝑑σ[0]M(s)ρ[0]M(s)=K4ρK4+K5ρK5=𝒥3(ρ),\displaystyle=\int_{1}^{2}d\sigma[0]\otimes\sqrt{M(s)}\rho[0]\otimes\sqrt{M(s)}=K_{4}^{\prime}\rho K_{4}^{\prime}+K_{5}^{\prime}\rho K_{5}^{\prime}={\cal J}_{3}(\rho), (36)

where we have used Eq. (VI). Thus, recalling that 𝒥1(ρ)=K1ρK1{\cal J}_{1}(\rho)=K_{1}^{\prime}\rho K_{1}^{\prime}, then by retaining the classical information about which outcome is obtained and when, one is able to coarse-grain the outcomes to separately obtain each of the three CP maps of the full quantum instrument.

VII LOCCLOCC¯\neq\overline{\textrm{LOCC}}

In this section, we consider the difference between LOCC and LOCC¯\overline{\textrm{LOCC}} for quantum channels and quantum measurements, the extremes of the range of quantum instruments. For quantum channels, we are unable to shed much light on this question. Certainly as we have seen in Section V, while the side-branches of the limit of the LOCC protocols implementing channels ν{\cal E}^{\nu} are all piecewise local, the limiting main branch, ΠM(s)\Pi_{M}(s), is not. At first blush, one may wonder if this demonstrates that channel {\cal E} cannot be exactly implemented by LOCC, which requires a piecewise local path to each outcome, even if one or more paths may contain an infinite number of pieces. However, it is possible that there is a different set of paths, possibly to the same Kraus representation of {\cal E} or possibly to a different one, all of which are piecewise local. We suspect that this is not the case but have not found a proof.

While we are presently unable to answer this question for quantum channels, we are able to show, we believe for the first time, that LOCC is not closed when considering implementation of quantum measurements, whether the measurements are viewed as POVMs or as generalized measurements. That is, we have the following theorem.

Theorem 4.

The set of LOCC measurements is not closed.

Proof.

We continue with the example of [10, 11], for which we have obtained the limiting POVM, see Eq. (V), of the sequence of protocols illustrated in Fig. 1, and the associated limiting generalized measurement, see Eq. (V). Our measurement operators for the limiting generalized measurement are [11][11], dσ[0]M(s)\sqrt{d\sigma}[0]\otimes\sqrt{M(s)}, and dσM(s)[0]\sqrt{d\sigma}\sqrt{M(s)}\otimes[0], and recall that 0σ10\leq\sigma\leq 1 and σ=s1\sigma=\sqrt{s}-1. We wish to show that this measurement cannot be implemented by LOCC, even though it is the limit of a sequence of LOCC measurements. To show this, we use an argument very similar to the one used in [11] for the same sequence of LOCC protocols, but there coarse-grained to the quantum instrument of Eq. (V). That is, we will consider the action of this generalized measurement—each element of which is tensored with the identity operator on a third party, CC—on the WW state, |W=(|001+|010+|100)/3|W\rangle=\left(|001\rangle+|010\rangle+|100\rangle\right)/\sqrt{3}. Then, it is not difficult to see that

(K1IC)|W\displaystyle(K_{1}\otimes I_{C})|W\rangle =([11]IC)|W=0\displaystyle=\left([11]\otimes I_{C}\right)|W\rangle=0
(K2(s)IC)|W\displaystyle(K_{2}(s)\otimes I_{C})|W\rangle =dσ([0]AM(s)BIC)|W=dσ3[0]A(σ|01+|10)BC\displaystyle=\sqrt{d\sigma}\left([0]_{A}\otimes\sqrt{M(s)_{B}}\otimes I_{C}\right)|W\rangle=\sqrt{\frac{d\sigma}{3}}[0]_{A}\otimes\left(\sqrt{\sigma}|01\rangle+|10\rangle\right)_{BC}
(K3(s)IC)|W\displaystyle(K_{3}(s)\otimes I_{C})|W\rangle =dσ([0]BM(s)AIC)|W=dσ3[0]B(σ|01+|10)AC.\displaystyle=\sqrt{d\sigma}\left([0]_{B}\otimes\sqrt{M(s)_{A}}\otimes I_{C}\right)|W\rangle=\sqrt{\frac{d\sigma}{3}}[0]_{B}\otimes\left(\sqrt{\sigma}|01\rangle+|10\rangle\right)_{AC}. (37)

The latter pair are a continuum of elements, and what we’re interested in are outcomes in the form (K2(s)IC)[W](K2(s)IC)\left(K_{2}(s)\otimes I_{C}\right)[W]\left(K_{2}(s)\otimes I_{C}\right)^{\dagger}, which we may sum up by integration over σ\sigma. This yields

01(K2(s)IC)[W](K2(s)IC)=[0]A([01]/3+4|0110|/9+4|1001|/9+2[10]/3)BC\displaystyle\int_{0}^{1}\left(K_{2}(s)\otimes I_{C}\right)[W]\left(K_{2}(s)\otimes I_{C}\right)^{\dagger}=[0]_{A}\otimes\left([01]/3+4|01\rangle\langle 10|/9+4|10\rangle\langle 01|/9+2[10]/3\right)_{BC} (38)

with probability equal to one-half, and the exact same thing with the same probability but with parties AA and BB interchanged for the term involving K3(s)K_{3}(s). Given that Eq. (38) is equal to (𝒥2C)([W])\left({\cal J}_{2}\otimes{\cal I}_{C}\right)([W]), compare Eq. (VI) (and similarly for (𝒥3C)([W])\left({\cal J}_{3}\otimes{\cal I}_{C}\right)([W]); C{\cal I}_{C} is the identity channel on the Hilbert space, C{\cal H}_{C}, representing system CC), this is in agreement with the result obtained in Section 4.34.3 of [11], from which they showed by a concurrence argument [37] that this could not be achieved by LOCC. That argument applies here, as well, proving that this is an example of a generalized measurement that is in LOCC¯\overline{\textrm{LOCC}} but not in LOCC. To put it in simple terms, if this measurement could be implemented by LOCC, then by coarse-graining, so could the instrument of Eq. (V), but this is known to be impossible [11], so the result follows directly.

Note that if the corresponding POVM could be implemented within LOCC, then so could the generalized measurement. Since we have just argued that the latter is impossible, we have also now provided an example of a POVM that is in LOCC¯\overline{\textrm{LOCC}} but not in LOCC, and this ends the proof.\blacksquare

It is perhaps of interest that this POVM involves an infinite number of outcomes, as does the corresponding generalized measurement. Therefore, one may ask whether or not LOCC is closed for measurements that are restricted to have a finite number of outcomes. If the latter question can be answered in the affirmative, it would show an interesting relationship between the closure properties of LOCC and the number of outcomes allowed in a measurement. The same relationship would then apparently also hold when considering the size of Kraus representation needed to implement a quantum channel or quantum instrument.

Note that these arguments require that one knows which outcomes are in the set, K2(s)K_{2}(s), and which are in K3(s)K_{3}(s), treating the two sets separately while coarse-graining over the range of ss. Thus, this example also requires classical information about the outcomes, so it does not apply to the full quantum channel, {\cal E}. Therefore, the question of whether or not LOCC is closed for implementing quantum channels remains open.

VIII Generalization of {\cal E} to any number of qubits

Consider a PP-qubit system with the parties performing the same local measurements as described in Section V. That is, each of their local POVM elements are R=(1ϵ)[0]+[1]R=(1-\epsilon)[0]+[1] and IRI-R. The parties measure one after the other, and if after any cycle of their PP local measurements, any one of them has gotten the second outcome, they stop. Only if they all got RR do they continue to the next cycle, and we shall refer to the continued sequence of these cycles as the main branch of the protocol and denote it as Πν(s)\Pi^{\nu}(s), as has been done in Section V for the two-qubit case. The parties always stop after ν\nu rounds, with ϵ=νc\epsilon=\nu^{-c} for some 0<c<10<c<1, even if they’ve always gotten the first outcome, having thus implemented channel ν{\cal E}^{\nu}. The limiting paths have a form that are a simple generalization of those for the two-qubit case.

Lemma 6.

In the limit ν\nu\to\infty (equivalently, ϵ0\epsilon\to 0), the main branch of this protocol has a limit given by ΠM(s)=limνΠν(s)=M(s)P\Pi_{M}(s)=\lim_{\nu\to\infty}\Pi^{\nu}(s)=M(s)^{\otimes P}, with M(s)=(sP1)[0]+[1]M(s)=(\sqrt[P]{s}-1)[0]+[1].

The proof, which requires one to show that limνΠM(s)Πν(s)1=0\lim_{\nu\to\infty}\left\|\Pi_{M}(s)-\Pi^{\nu}(s)\right\|_{1}=0, is straightforward and is left as an exercise for the reader. We share a few hints. Note that after nn complete cycles followed by ll of the parties having completed their part of the (n+1)(n+1)st cycle, the segment of the (piecewise local) path from this point to the point after the (l+1)(l+1)st party performs the next measurement may be written as

Πν(s)=stn+1l+1tnPl1tn+1ltnPl1(tntn+1)(Rn+1)l(Rn)(Pl)+tn+1ltnPlstn+1ltnPl1(tntn+1)(Rn+1)l+1(Rn)(Pl1),\displaystyle\Pi^{\nu}(s)=\frac{s-t_{n+1}^{l+1}t_{n}^{P-l-1}}{t_{n+1}^{l}t_{n}^{P-l-1}(t_{n}-t_{n+1})}\left(R^{n+1}\right)^{\otimes l}\left(R^{n}\right)^{\otimes(P-l)}+\frac{t_{n+1}^{l}t_{n}^{P-l}-s}{t_{n+1}^{l}t_{n}^{P-l-1}(t_{n}-t_{n+1})}\left(R^{n+1}\right)^{\otimes l+1}\left(R^{n}\right)^{\otimes(P-l-1)}, (39)

with tn+1ltnPlstn+1l+1tnPl1t_{n+1}^{l}t_{n}^{P-l}\geq s\geq t_{n+1}^{l+1}t_{n}^{P-l-1}, and tn=Tr(Rn)=(1ϵ)n+1t_{n}=\textrm{Tr}\left(R^{n}\right)=\left(1-\epsilon\right)^{n}+1. Since all operators are diagonal, the proof of this lemma simply requires one to show that |N|(ΠM(s)Πν(s))|N|ϵP\left|\langle N|\left(\Pi_{M}(s)-\Pi^{\nu}(s)\right)|N\rangle\right|\leq\epsilon P for every 0N2P10\leq N\leq 2^{P}-1, where using the binary expansion, N=bP1b1b0N=b_{P-1}\ldots b_{1}b_{0} with bj=0b_{j}=0 or 11, we define

|N=|b0|b1|bP1.\displaystyle|N\rangle=|b_{0}\rangle|b_{1}\rangle\ldots|b_{P-1}\rangle. (40)

Recall that all non-zero local matrix elements are equal to unity unless that party’s part of the binary expansion of NN is equal to 0, in which case one obtains a factor of a power of (1ϵ)\left(1-\epsilon\right); for example, 0|Rn|0=(1ϵ)n=tn1\langle 0|R^{n}|0\rangle=(1-\epsilon)^{n}=t_{n}-1.

As done in the case of P=2P=2 above, we may differentiate ΠM(s)\Pi_{M}(s) to obtain other possible Kraus operators for our limiting channel, P{\cal E}_{P}. This gives

dΠM(s)=dsP([0]M(s)(P1)+M(s)[0]M(s)(P2)++M(s)(P1)[0]),\displaystyle d\Pi_{M}(s)=d\!\!\sqrt[P]{s}\left([0]\otimes M(s)^{\otimes(P-1)}+M(s)\otimes[0]\otimes M(s)^{\otimes(P-2)}+\ldots+M(s)^{\otimes(P-1)}\otimes[0]\right), (41)

suggesting a separable Kraus representation as consisting of [111][11\ldots 1] and each of the individual terms appearing in the preceding equation. In fact, we can show this is correct by calculating the Choi-Jamiolkowski operator for this conjectured channel, as well as the limit of those for the approximating sequence of channels described above in terms of local operators RR. In the following, |i,|j|i\rangle,|j\rangle are to be understood as a binary expansion, in similar fashion to that for |N|N\rangle in Eq. (40), lil_{i} is the number of 0’s appearing in that binary expansion of ii, and mijm_{ij} counts the number of positions where both ii and jj each have a zero. In other words, li=h(¬i)l_{i}=h(\neg i), the Hamming weight of ¬i\neg{i} with ¬\neg denoting the bitwise NOT function, and mij=h(¬(i|j))m_{ij}=h\left(\neg(i|j)\right) with || the bitwise OR function. Then, we have

Lemma 7.

The Choi-Jamiolkowski operator for channel P{\cal E}_{P} and that for the ν\nu\to\infty limit of the above-described approximating channels, ν{\cal E}^{\nu}, are both given by

ΦP=2P[2P1][2P1]+2Pi,j=02P2|ij|2mijli+lj|ij|,\displaystyle\Phi_{P}=2^{-P}[2^{P}-1]\otimes[2^{P}-1]+2^{-P}\sum_{i,j=0}^{2^{P}-2}|i\rangle\langle j|\otimes\frac{2m_{ij}}{l_{i}+l_{j}}|i\rangle\langle j|, (42)

and since this is equal to 2Pi,j=02P1|ij|P(|ij|)2^{-P}\sum_{i,j=0}^{2^{P}-1}|i\rangle\langle j|\otimes{\cal E}_{P}(|i\rangle\langle j|), the corresponding channels are given by

P(ρ)=2P1|ρ|2P1[2P1]+i,j=02P22mijρijli+lj|ij|,\displaystyle{\cal E}_{P}(\rho)=\langle 2^{P}-1|\rho|2^{P}-1\rangle[2^{P}-1]+\sum_{i,j=0}^{2^{P}-2}\frac{2m_{ij}\rho_{ij}}{l_{i}+l_{j}}|i\rangle\langle j|, (43)

with ρij=i|ρ|j\rho_{ij}=\langle i|\rho|j\rangle.

The proof is given in Appendix C.

One may consider the action of P{\cal E}_{P} on the state [W][111][W]\otimes[11\ldots 1] where the first tensor factor is a projector onto the WW-state of three qubits and the second one is the tensor product of P3P-3 local (qubit) projectors onto the |1|1\rangle state. Since the action of M(s)M(s) on |1|1\rangle is identical to the action of the identity operator, as is the action of [1][1], whereas [0][0] annihilates |1|1\rangle, the last P3P-3 parties’ actions are irrelevant to the evolution of the state of the other three parties. It has been shown that the transformation 3([W]){\cal E}_{3}([W]) cannot be achieved exactly by LOCC [10]. Therefore, it is also the case that P([W]){\cal E}_{P}([W]) cannot be implemented exactly by LOCC—since if it could, then so could 3([W]){\cal E}_{3}([W])—even though, as we’ve demonstrated above, LOCC can be used to approximate P{\cal E}_{P} arbitrarily closely. Given the results of [11, 10], these last points are perhaps unsurprising, but the limiting measurement operators in Eq. (41) may be of some interest.

IX Conclusions

In this paper, we have extended our earlier result providing a necessary condition that a quantum measurement can be closely approximated by LOCC [18, 20], to a similar condition which now applies to the full range of quantum instruments from the most refined up to the fully coarse-grained case of quantum channels. The proof of this new result includes a demonstration that if a sequence of quantum instruments converges to some limiting instrument, then there also exists a sequence of geometric objects known as zonoids, which converges to a limiting zonoid. Each such zonoid is uniquely determined by the instrument from which it is derived, and the limiting zonoid plays a critical role in determining whether or not the limiting instrument can be closely approximated by LOCC. When this approximation is possible, there must exist certain paths of product operators lying within the limiting zonoid, see Theorems 2 and 3.

We have studied in detail a particular quantum instrument on two-qubits, {\cal E}, which has been shown [11] to be the limit of a sequence of instruments, each of which can be implemented by a finite-round LOCC protocol, but that {\cal E} itself cannot be exactly achieved using LOCC. Our analysis shows why this is so. For each of the approximating instruments in the given sequence, there exists a set of piecewise local paths to each of a set of outcomes, where that set of outcomes can be coarse-grained to represent the given instrument. When one considers the limit of these paths, however, one finds that while they correspond to a set of outcomes representing {\cal E}, one of the limiting paths fails to remain piecewise local in the limit. Since for LOCC implementation, a crucial characteristic of these paths is that they must be piecewise local, see the discussion above Theorem 1, this failure is consistent with the finding that {\cal E} cannot be exactly implemented by LOCC [11].

We initially treated this example as a quantum channel, see Section V, and found the Kraus representation of {\cal E} corresponding to the limit of the Kraus representations for the approximating channels ν{\cal E}^{\nu}. This allowed us to identify the isometry which transforms the minimal representation of {\cal E} to the representation implemented in the limit, which contains an infinite number of Kraus operators. Then after generalizing Theorem 2 to the case of general quantum instruments, Theorem 3, we treated this same example as a (less than fully coarse-grained) quantum instrument, see Section VI. For the latter case, it was shown [11] that this instrument cannot be exactly implemented by LOCC, but when considered as a channel, a proof of this statement is still lacking. Thus, it remains an open question whether LOCC=LOCC¯=\overline{\textrm{LOCC}} in the case of quantum channels. On the other hand, we have succeeded in answering this same question for the case of quantum measurements, which as far as we are aware, had not been previously answered. That is, we have given an example of a quantum measurement taken directly from the limiting Kraus representation for the above-mentioned {\cal E}, which can be arbitrarily closely approximated, but not implemented exactly, by LOCC. This measurement consists of an infinite number of outcomes.

We conjecture that LOCC is closed for measurements that are restricted to have a finite number of outcomes. It is certainly possible this is false, but our reasoning is as follows: That LOCC is not closed would mean that within ZZ_{\cal M}, the zonotope associated with measurement {\cal M}, any sequence of protocols whose limiting measurement is {\cal M} leads to a limiting path that is not piecewise local, since otherwise it would be LOCC. But as shown in our example of Section V, obtaining a path that isn’t piecewise local appears to require that the sequence of approximating LOCC protocols involves paths that branch off at closer and closer intervals, leaving a continuum of branches in the limit. This continuum of branches implies a continuum of leaf nodes at the ends of those branches, and this suggests an infinite number of outcomes. However, it is possible to have multiple leaf nodes that are equal to each other, including possibly an infinite number of them. Therefore, it is not entirely clear whether or not the conjecture is correct. Of course, in our example we find that each leaf node is unique, and it is at least conceivable that this must be true for any case demonstrating this property of sequences of piecewise local paths having a limiting path that is not piecewise local.

Finally, we would like to encourage further study of multipartite quantum channels, for which there appear to be relatively few examples of physical interest that have been studied in the literature. Hopefully, the results presented here will stimulate new interest in exploring such avenues of research.

Acknowledgments — We would like to thank Vlad Gheorghiu and Fei Shi for comments on an early version of this manuscript.

Appendix A Convergence of channels implies convergence of the associated zonoids

Suppose we have a sequence of LOCC protocols 𝒫ν{\cal P}^{\nu}, each implementing a quantum channel ν{\cal E}^{\nu}, and limνν=\lim_{\nu\to\infty}{\cal E}^{\nu}={\cal E}, where :()(o){\cal E}:{\cal B}({\cal H})\to{\cal B}({\cal H}_{o}) and ν:()(oν){\cal E}^{\nu}:{\cal B}({\cal H})\to{\cal B}({\cal H}_{o}^{\nu}), having recognized that the input Hilbert space (){\cal B}({\cal H}) must be the same for all these channels. (Recall that superscript ν\nu denotes the position in the sequence of quantum channels and is not to be read as an exponent.) That is, for every δ>0\delta>0, there exists natural number μ\mu such that for every ν>μ\nu>\mu

δ/16>ν,\displaystyle\delta/16>\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}, (44)

where \left\|\cdot\right\|_{\diamond} is the diamond norm, see Eq. (47) below. From the results of Ref. [29, 28], we have that

δ/4>νinfV,VνVVνν+ν.\displaystyle\sqrt{\delta}/4>\sqrt{\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}}\geq\inf_{V,V^{\nu}}\left\|V-V^{\nu}\right\|_{\infty}\geq\frac{\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}}{\sqrt{\left\|{\cal E}\right\|_{\diamond}}+\sqrt{\left\|{\cal E}^{\nu}\right\|_{\diamond}}}. (45)

where the operator norm, X\left\|X\right\|_{\infty}, is equal to the largest singular value of XX, and the infimum appearing here is over all Stinespring dilations [35], V:oeV:{\cal H}\to{\cal H}_{o}\otimes{\cal H}_{e} and Vν:oνeνV^{\nu}:{\cal H}\to{\cal H}_{o}^{\nu}\otimes{\cal H}_{e}^{\nu}, of channels {\cal E} and ν{\cal E}^{\nu}, respectively. Below, we will make use of channels, :e{\cal F}:{\cal H}\to{\cal H}_{e} complementary to {\cal E}, and ν:eν{\cal F}^{\nu}:{\cal H}\to{\cal H}_{e}^{\nu} complementary to ν{\cal E}^{\nu}, and without loss of generality, we will consider minimal representations of {\cal E} (ν)({\cal E}^{\nu}), so that the dimension of e{\cal H}_{e} (eν)({\cal H}_{e}^{\nu}) is equal to the Kraus rank κ\kappa of {\cal E} (κν\kappa^{\nu} of ν{\cal E}^{\nu}). The dimensions of ,o{\cal H},{\cal H}_{o} and oν{\cal H}_{o}^{\nu} are d,dod,d_{o} and doνd_{o}^{\nu}, respectively.

The trace norm of a quantum channel TT is defined as

T1=maxX11T(X)1,\displaystyle\left\|T\right\|_{1}=\max_{\left\|X\right\|_{1}\leq 1}\left\|T(X)\right\|_{1}, (46)

where X1=Tr(XX)\left\|X\right\|_{1}=\textrm{Tr}\left(\sqrt{X^{\dagger}X}\right) is the trace norm of operator XX. For X0X\geq 0, we have that X1=Tr(X)\left\|X\right\|_{1}=\textrm{Tr}\left(X\right). With {\cal I}_{\cal H} the identity channel on input Hilbert space {\cal H}, recall the definition of the diamond norm as

T=T1=maxX11(T)X1.\displaystyle\left\|T\right\|_{\diamond}=\left\|T\otimes{\cal I}_{\cal H}\right\|_{1}=\max_{\left\|X\right\|_{1}\leq 1}\left\|(T\otimes{\cal I}_{\cal H})X\right\|_{1}. (47)

It is not difficult to show (by use of Eq. (77) in Ref. [31], for example) that for any quantum channel {\cal E},

=1.\displaystyle\left\|{\cal E}\right\|_{\diamond}=1. (48)

Then, Eq. (45) becomes

δ/4>νinfV,VνVVν12ν.\displaystyle\sqrt{\delta}/4>\sqrt{\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}}\geq\inf_{V,V^{\nu}}\left\|V-V^{\nu}\right\|_{\infty}\geq\frac{1}{2}\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}. (49)

From Eq. (44) we have an (arbitrarily small) upper bound proportional to δ\delta on ν\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}, but we are here seeking to find a similar bound on the Hausdorff distance between corresponding zonotopes, Z,ZνZ_{\cal E},Z_{{\cal E}^{\nu}}. The following arguments achieve that goal.

While the infimum in Eq. (49) may not be achieved, it is always arbitrarily close to at least one point in the set over which it is taken. Therefore, for every δ>0\delta>0 there exists Stinespring dilations V~,V~ν\tilde{V},\tilde{V}^{\nu} of ,ν{\cal E},{\cal E}^{\nu}, respectively, such that δ/2>infV,VνVVν+δ/4V~V~νinfV,VνVVν\sqrt{\delta}/2>\inf_{V,V^{\nu}}\left\|V-V^{\nu}\right\|_{\infty}+\sqrt{\delta}/4\geq\left\|\tilde{V}-\tilde{V}^{\nu}\right\|_{\infty}\geq\inf_{V,V^{\nu}}\left\|V-V^{\nu}\right\|_{\infty}. Therefore,

δ>2V~V~νν.\displaystyle\sqrt{\delta}>2\left\|\tilde{V}-\tilde{V}^{\nu}\right\|_{\infty}\geq\left\|{\cal E}-{\cal E}^{\nu}\right\|_{\diamond}. (50)

Since each of these isometries V~\tilde{V} (V~ν)(\tilde{V}^{\nu}) is a dilation not only of channel {\cal E} (ν)({\cal E}^{\nu}) but also of a corresponding complementary channel {\cal F} (ν)({\cal F}^{\nu}), it also holds that [28]

δ>2V~V~νν.\displaystyle\sqrt{\delta}>2\left\|\tilde{V}-\tilde{V}^{\nu}\right\|_{\infty}\geq\left\|{\cal F}-{\cal F}^{\nu}\right\|_{\diamond}. (51)

Since the diamond norm is a maximum over inputs XX, see Eq. (47), we get a lower bound on this norm by evaluating ν{\cal F}-{\cal F}^{\nu} acting on any density operator XX, so choose X=|ΩΩ|X=|\Omega\rangle\langle\Omega|, where |Ω=i|i|i/d|\Omega\rangle=\sum_{i}|i\rangle\otimes|i\rangle/\sqrt{d} is a maximally entangled state on input Hilbert space {\cal H}\otimes{\cal H}. Then, we obtain

δ>2V~V~νΦΦν1,\displaystyle\sqrt{\delta}>2\left\|\tilde{V}-\tilde{V}^{\nu}\right\|_{\infty}\geq\left\|\Phi-\Phi^{\nu}\right\|_{1}, (52)

where Φ=(I)(|ΩΩ|)\Phi=(I_{\cal H}\otimes{\cal F})(|\Omega\rangle\langle\Omega|) is the Choi-Jamiolkowski operator [12, 25] for channel {\cal F}, and similarly for Φν\Phi^{\nu} and ν{\cal F}^{\nu}.

Consider the two Kraus representations {Km}\{K_{m}\} of {\cal E} with ded_{e} members, and {Kmν}\{K_{m}^{\nu}\} of ν{\cal E}^{\nu} with deνd_{e}^{\nu} members. Define d~e=max(de,deν)\tilde{d}_{e}=\max(d_{e},d_{e}^{\nu}), expand the smaller of e,eν{\cal H}_{e},{\cal H}_{e}^{\nu} to dimension d~e\tilde{d}_{e}, and add extra zero operators to the smaller of the two representations, so that these each have d~e\tilde{d}_{e} operators. For the smaller of these two Hilbert spaces, {\cal F} or ν{\cal F}^{\nu} still maps into the same unexpanded space, e{\cal H}_{e} or eν{\cal H}_{e}^{\nu}, which is now a subspace of the new, expanded space. This can be seen, in the case that de<deνd_{e}<d_{e}^{\nu} for example, from the relation noted below Eq. (2) of the main text, l|Km=m|Ql\langle l|K_{m}=\langle m|Q_{l}, implying that for those mm such that KmK_{m} is the null operator, the added kets correspond to a subspace of e{\cal H}_{e} which {\cal F} does not map into. With these considerations, we have

Φ\displaystyle\Phi =1di,j=1dl=1do|ij|Ql|ij|Ql\displaystyle=\frac{1}{d}\sum_{i,j=1}^{d}\sum_{l=1}^{d_{o}}|i\rangle\langle j|\otimes Q_{l}|i\rangle\langle j|Q_{l}^{\dagger}
=1di,j=1dm,m=1d~el=1do|ij||mm|Ql|ij|Ql|mm|\displaystyle=\frac{1}{d}\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}\sum_{l=1}^{d_{o}}|i\rangle\langle j|\otimes|m^{\prime}\rangle\langle m^{\prime}|Q_{l}|i\rangle\langle j|Q_{l}^{\dagger}|m\rangle\langle m|
=1di,j=1dm,m=1d~el=1do|ij|l|Km|ij|Km|l|mm|\displaystyle=\frac{1}{d}\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}\sum_{l=1}^{d_{o}}|i\rangle\langle j|\otimes\langle l|K_{m^{\prime}}|i\rangle\langle j|K_{m}^{\dagger}|l\rangle|m^{\prime}\rangle\langle m|
=1di,j=1dm,m=1d~el=1do|ij|j|Km|ll|Km|i|mm|\displaystyle=\frac{1}{d}\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}\sum_{l=1}^{d_{o}}|i\rangle\langle j|\otimes\langle j|K_{m}^{\dagger}|l\rangle\langle l|K_{m^{\prime}}|i\rangle|m^{\prime}\rangle\langle m|
=1di,j=1dm,m=1d~ej|KmKm|i|ij||mm|,\displaystyle=\frac{1}{d}\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}\langle j|K_{m}^{\dagger}K_{m^{\prime}}|i\rangle|i\rangle\langle j|\otimes|m^{\prime}\rangle\langle m|, (53)

and similarly for Φν\Phi^{\nu}, the only difference being the replacement of KmK_{m} by KmνK_{m}^{\nu}. In going from the second to the third line of Eq. (A), we have used the relationship, m|Ql=l|Km\langle m|Q_{l}=\langle l|K_{m}, see just below Eq. (2) of the main text. Note that matrices Φ\Phi and Φν\Phi^{\nu} are of the same size, and we have

ΦΦν1=1di,j=1dm,m=1d~ej|(KmKmKmνKmν)|i|ij||mm|1.\displaystyle\left\|\Phi-\Phi^{\nu}\right\|_{1}=\frac{1}{d}\left\|\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}\langle j|\left(K_{m}^{\dagger}K_{m^{\prime}}-K_{m}^{\nu{\dagger}}K_{m^{\prime}}^{\nu}\right)|i\rangle|i\rangle\langle j|\otimes|m^{\prime}\rangle\langle m|\right\|_{1}. (54)

Let K^n\hat{K}_{n} (K^nν)\hat{K}_{n}^{\nu}) be one of κ\kappa (κν)(\kappa^{\nu}) members of a minimal Kraus representation for {\cal E} (ν)({\cal E}^{\nu}). Then, Km=nWmnK^nK_{m}=\sum_{n}W_{mn}\hat{K}_{n} with WW a d~e×κ\tilde{d}_{e}\times\kappa isometry so that WW=IκW^{\dagger}W=I_{\kappa}, and Kmν=nWmnνK^nνK_{m}^{\nu}=\sum_{n}W_{mn}^{\nu}\hat{K}_{n}^{\nu} with WνW^{\nu} a d~e×κν\tilde{d}_{e}\times\kappa^{\nu} isometry so that WνWν=IκνW^{\nu{\dagger}}W^{\nu}=I_{\kappa^{\nu}}. Recalling from Lemma 4 that κνκ\kappa^{\nu}\geq\kappa, add zero operators if needed to {K^n}\{\hat{K}_{n}\}, so that this set also has κν\kappa^{\nu} members, along with zero columns to WW so that WW=PκW^{\dagger}W=P_{\kappa}, a κν×κν\kappa^{\nu}\times\kappa^{\nu} projector onto a κ\kappa-dimensional subspace. Then, from Eqs. (52) and (54), we can write

δ\displaystyle\sqrt{\delta} >1di,j=1dm,m=1d~ej|(KmKmKmνKmν)|i|ij||mm|1\displaystyle>\frac{1}{d}\left\|\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}\langle j|\left(K_{m}^{\dagger}K_{m^{\prime}}-K_{m}^{\nu{\dagger}}K_{m^{\prime}}^{\nu}\right)|i\rangle|i\rangle\langle j|\otimes|m^{\prime}\rangle\langle m|\right\|_{1}
=1d(i,j=1dm,m=1d~e|j|(KmKmKmνKmν)|i|2)1/2\displaystyle=\frac{1}{d}\left(\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}\left|\langle j|\left(K_{m}^{\dagger}K_{m^{\prime}}-K_{m}^{\nu{\dagger}}K_{m^{\prime}}^{\nu}\right)|i\rangle\right|^{2}\right)^{1/2}
1dsup0CId~e(i,j=1dm,m=1d~e|Cmmj|(KmKmKmνKmν)|i|2)1/2\displaystyle\geq\frac{1}{d}\sup_{0\leq C\leq I_{\tilde{d}_{e}}}\left(\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}\left|C_{mm^{\prime}}\langle j|\left(K_{m}^{\dagger}K_{m^{\prime}}-K_{m}^{\nu{\dagger}}K_{m^{\prime}}^{\nu}\right)|i\rangle\right|^{2}\right)^{1/2}
1dsup0CId~e(i,j=1d|m,m=1d~eCmmj|(KmKmKmνKmν)|i|2)1/2\displaystyle\geq\frac{1}{d}\sup_{0\leq C\leq I_{\tilde{d}_{e}}}\left(\sum_{i,j=1}^{d}\left|\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}C_{mm^{\prime}}\langle j|\left(K_{m}^{\dagger}K_{m^{\prime}}-K_{m}^{\nu{\dagger}}K_{m^{\prime}}^{\nu}\right)|i\rangle\right|^{2}\right)^{1/2}
=1dsup0CId~einf0CνId~e(i,j=1d|m,m=1d~en,n=1κνj|(WmnCmmWmnK^nK^nWmnνCmmνWmnνK^nνK^nν)|i|2)1/2\displaystyle=\frac{1}{d}\sup_{0\leq C\leq I_{\tilde{d}_{e}}}\inf_{0\leq C^{\nu}\leq I_{\tilde{d}_{e}}}\left(\sum_{i,j=1}^{d}\left|\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}\sum_{n,n^{\prime}=1}^{\kappa^{\nu}}\langle j|\left(W_{mn}^{\ast}C_{mm^{\prime}}W_{m^{\prime}n^{\prime}}\hat{K}_{n}^{\dagger}\hat{K}_{n^{\prime}}-W_{mn}^{\nu\ast}C_{mm^{\prime}}^{\nu}W_{m^{\prime}n^{\prime}}^{\nu}\hat{K}_{n}^{\nu{\dagger}}\hat{K}_{n^{\prime}}^{\nu}\right)|i\rangle\right|^{2}\right)^{1/2}
=1dsup0C^Pκinf0C^νIκν(i,j=1d|n,n=1κνj|(C^nnK^nK^nC^nnνK^nνK^nν)|i|2)1/2\displaystyle=\frac{1}{d}\sup_{0\leq\hat{C}\leq P_{\kappa}}\inf_{0\leq\hat{C}^{\nu}\leq I_{\kappa^{\nu}}}\left(\sum_{i,j=1}^{d}\left|\sum_{n,n^{\prime}=1}^{\kappa^{\nu}}\langle j|\left(\hat{C}_{nn^{\prime}}\hat{K}_{n}^{\dagger}\hat{K}_{n^{\prime}}-\hat{C}_{nn^{\prime}}^{\nu}\hat{K}_{n}^{\nu{\dagger}}\hat{K}_{n^{\prime}}^{\nu}\right)|i\rangle\right|^{2}\right)^{1/2}
=1dsup0C^Pκinf0C^νIκνn,n=1κν(C^nnK^nK^nC^nnνK^nνK^nν)1,\displaystyle=\frac{1}{d}\sup_{0\leq\hat{C}\leq P_{\kappa}}\inf_{0\leq\hat{C}^{\nu}\leq I_{\kappa^{\nu}}}\left\|\sum_{n,n^{\prime}=1}^{\kappa^{\nu}}\left(\hat{C}_{nn^{\prime}}\hat{K}_{n}^{\dagger}\hat{K}_{n^{\prime}}-\hat{C}_{nn^{\prime}}^{\nu}\hat{K}_{n}^{\nu{\dagger}}\hat{K}_{n^{\prime}}^{\nu}\right)\right\|_{1}, (55)

with C^=WCW\hat{C}=W^{\dagger}CW and C^ν=WνCνWν\hat{C}^{\nu}=W^{\nu{\dagger}}C^{\nu}W^{\nu}. To obtain the third line, we inserted |Cmm|1\left|C_{mm^{\prime}}\right|\leq 1, which will hold for all m,mm,m^{\prime} when 0CId~e0\leq C\leq I_{\tilde{d}_{e}}, and then we used the triangle inequality in the line after that. In entirely analogous fashion, one finds that

δ>1dsup0C^νIκνinf0C^Pκn,n=1κν(C^nnK^nK^nC^nnνK^nνK^nν)1.\displaystyle\sqrt{\delta}>\frac{1}{d}\sup_{0\leq\hat{C}^{\nu}\leq I_{\kappa^{\nu}}}\inf_{0\leq\hat{C}\leq P_{\kappa}}\left\|\sum_{n,n^{\prime}=1}^{\kappa^{\nu}}\left(\hat{C}_{nn^{\prime}}\hat{K}_{n}^{\dagger}\hat{K}_{n^{\prime}}-\hat{C}_{nn^{\prime}}^{\nu}\hat{K}_{n}^{\nu{\dagger}}\hat{K}_{n^{\prime}}^{\nu}\right)\right\|_{1}. (56)

Setting δ=(ϵ/d)2\delta=(\epsilon/d)^{2}, it immediately follows from Eqs. 44, A and 56 that for every ϵ>0\epsilon>0 there exists μ\mu such that for every ν>μ\nu>\mu

ϵ>dH(Z,Zν),\displaystyle\epsilon>d_{H}\left(Z_{\cal E},Z_{{\cal E}^{\nu}}\right), (57)

where dH(X,Y)=max{supxXinfyY|xy|,supyYinfxX|xy|}d_{H}(X,Y)=\max\{\sup_{x\in X}\inf_{y\in Y}\left|x-y\right|,\sup_{y\in Y}\inf_{x\in X}\left|x-y\right|\} is the Hausdorff distance between sets X,YX,Y, with |xy|\left|x-y\right| the distance between elements x,yx,y, which for present purposes is the operator trace norm, 1\left\|\cdot\right\|_{1}.

We therefore have that the sequence of zonoids, ZνZ_{{\cal E}^{\nu}}, associated with quantum channels ν{\cal E}^{\nu}, converges to the zonoid ZZ_{\cal E}, associated with the limiting quantum channel {\cal E}. This completes the proof of Lemma 3.

Appendix B Generalization of ZZ_{\cal E} for the case of quantum instruments

Assume that the sequence of instruments, ν{\cal E}^{\nu}, has as its limit instrument ={𝒥r}{\cal E}=\{{\cal J}_{r}\}. To formalize the intuitive arguments leading up to the statement of Theorem 3 in Section VI, let us consider qc{\cal E}_{qc} the quantum-classical channel, ρr𝒥r(ρ)[r]c\rho\to\sum_{r}{\cal J}_{r}(\rho)\otimes[r]_{c}, uniquely associated with instrument {\cal E} [11], and qcν{\cal E}_{qc}^{\nu} associated with each ν{\cal E}^{\nu}, as well. Then, limνqcν=qc\lim_{\nu\to\infty}{\cal E}_{qc}^{\nu}={\cal E}_{qc}. Here, subscript cc on [r]c[r]_{c} denotes an ancillary output containing the classical information as to which CP map, 𝒥r{\cal J}_{r}, was obtained. As these are quantum channels, the analysis in Appendix A holds, but let us take care in identifying the Kraus operators representing these channels. Recognizing that a set of Kraus operators for qc{\cal E}_{qc} is Km|f(m)cK_{m}\otimes|f(m)\rangle_{c}, where f(m)=rf(m)=r for KmK_{m} a Kraus operator representing 𝒥r{\cal J}_{r}, then returning to Eq. (A) and noting that the first two lines there remain unchanged, we may follow through the same steps (but with minor re-arrangement of factors to make clear how the classical register is to be handled) to derive the Choi-Jamiolkowski operator for the channel complementary to qc{\cal E}_{qc}. We find

Φ\displaystyle\Phi =1di,j=1dm,m=1d~el=1do|ij||mm|j|(Kmf(m)|c)(|ll|Ic)(Km|f(m)c)|i\displaystyle=\frac{1}{d}\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}\sum_{l=1}^{d_{o}}|i\rangle\langle j|\otimes|m^{\prime}\rangle\langle m|\langle j|\left(K_{m}^{\dagger}\otimes\langle f(m)|_{c}\right)\left(|l\rangle\langle l|\otimes I_{c}\right)\left(K_{m^{\prime}}\otimes|f(m^{\prime})\rangle_{c}\right)|i\rangle
=1di,j=1dm,m=1d~e|ij||mm|j|KmKm|if(m)|f(m)c\displaystyle=\frac{1}{d}\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}|i\rangle\langle j|\otimes|m\rangle\langle m^{\prime}|\langle j|K_{m}^{\dagger}K_{m^{\prime}}|i\rangle\langle f(m)|f(m^{\prime})\rangle_{c}
=1di,j=1dm,m=1d~e|ij||mm|j|KmKm|iδf(m),f(m),\displaystyle=\frac{1}{d}\sum_{i,j=1}^{d}\sum_{m,m^{\prime}=1}^{\tilde{d}_{e}}|i\rangle\langle j|\otimes|m\rangle\langle m^{\prime}|\langle j|K_{m}^{\dagger}K_{m^{\prime}}|i\rangle\delta_{f(m),f(m^{\prime})}, (58)

and similarly for Φν\Phi^{\nu} associated with each qcν{\cal E}_{qc}^{\nu} but with a different coarse-graining map, fν(m)f^{\nu}(m). Due to the appearance of the Kronecker delta in the last line, this is a block-diagonal operator, each block associated with one of the CP maps 𝒥r{\cal J}_{r}. Furthermore, each CP map has its own set of minimal Kraus operators, the collection of which (over all the CP maps in instrument {\cal E}) we denote as {K^n}\{\hat{K}_{n}\}. For the associated quantum-classical channel, these Kraus operators become K^n|g(n)c\hat{K}_{n}\otimes|g(n)\rangle_{c}, where g(n)g(n) is another coarse-graining map, possibly having a different domain than, but the same range as, f(n)f(n). As before, g(n)=rg(n)=r for those operators representing 𝒥r{\cal J}_{r}. Then, isometric freedom of this channel implies

Km|f(m)c\displaystyle K_{m}\otimes|f(m)\rangle_{c} =nWmnK^n|g(n)c\displaystyle=\sum_{n}W_{mn}\hat{K}_{n}\otimes|g(n)\rangle_{c}
Kmν|fν(m)c\displaystyle K_{m}^{\nu}\otimes|f^{\nu}(m)\rangle_{c} =nWmnνK^nν|gν(n)c.\displaystyle=\sum_{n}W_{mn}^{\nu}\hat{K}_{n}^{\nu}\otimes|g^{\nu}(n)\rangle_{c}. (59)

Clearly, since the kets appearing here represent classical states and are therefore mutually orthogonal (and we can always take them to be normalized by absorbing a factor into KmK_{m} or K^n\hat{K}_{n}, for example), there should be no contributions on the right-hand sides of these expressions except for those nn satisfying g(n)=f(m)g(n)=f(m) or gν(n)=fν(m)g^{\nu}(n)=f^{\nu}(m). That is to say, isometries W(Wν)W~{}(W^{\nu}), introduced on the fifth line of Eq. (A), are here block diagonal with blocks matching those appearing in C(Cν)C~{}(C^{\nu}), also introduced in Eq. (A). By arguments completely analogous to those in the proof of Lemma 4, the blocks defined by fν(m)f^{\nu}(m) must contain those blocks defined by f(m)f(m), at least for large enough ν\nu. This means that a CP map of ν{\cal E}^{\nu} cannot have smaller Kraus rank than the corresponding CP map of {\cal E} (or possibly CP maps, plural, since the blocks defined by fν(m)f^{\nu}(m) may contain more than one of those defined by f(m)f(m)).

Thus, by continuing with the arguments following Eq. (A) in Appendix A, one obtains the exact same result as previously obtained in Eqs. (56) and (57), but with the additional condition that matrices C^(C^ν)\hat{C}~{}(\hat{C}^{\nu}) appearing in the definition of Z(Zν)Z_{\cal E}~{}(Z_{{\cal E}^{\nu}}) must be block diagonal. The remainder of the proof of Theorem 3 follows precisely the arguments in the proof of Theorem 2.

Appendix C Proof of Lemma 7

Define

Kα(s)=dsPM(s)(α1)[0]M(s)(Pα),\displaystyle K_{\alpha}(s)=\sqrt{d\!\!\sqrt[P]{s}M(s)^{\otimes(\alpha-1)}\otimes[0]\otimes M(s)^{\otimes(P-\alpha)}}, (60)

where the [0][0] appearing here is a local operator for party α\alpha. Then, defining σ=sP1\sigma=\sqrt[P]{s}-1, we have that

01𝑑σKα(s)|ij|Kα=01𝑑σσ(li+lj2)|ij|=2li+lj|ij|,\displaystyle\int_{0}^{1}d\sigma K_{\alpha}(s)|i\rangle\langle j|K_{\alpha}^{\dagger}=\int_{0}^{1}d\sigma\sqrt{\sigma}^{(l_{i}+l_{j}-2)}|i\rangle\langle j|=\frac{2}{l_{i}+l_{j}}|i\rangle\langle j|, (61)

when the binary representation of ii and jj each have a 0 in the α\alphath position, whereas this integral vanishes otherwise. Summing over all α\alpha gives an additional factor of mijm_{ij}. This accounts for all the Kraus operators except for the one equal to [2P1][2^{P}-1], which gives [2P1][2^{P}-1] when acting on itself, and zero otherwise. Therefore, adding all the terms together, we see that Eq. (42) and Eq. (43) in the main text are correct for channel P{\cal E}_{P}.

For the approximating channel ν{\cal E}^{\nu}, the Kraus operators are KMν=RνK_{M}^{\nu}=\sqrt{R}^{\otimes\nu} at the end of the main branch, which asymptotically approaches [2P1][2^{P}-1] in the limit ν\nu\to\infty (equivalently, ϵ0\epsilon\to 0). In addition, we have Kαnν=ϵ(1ϵ)(n1)Rn(α1)[0]Rn(Pα1)K_{\alpha n}^{\nu}=\sqrt{\epsilon(1-\epsilon)^{(n-1)}}\sqrt{R^{n}}^{\otimes(\alpha-1)}\otimes[0]\otimes\sqrt{R^{n}}^{\otimes(P-\alpha-1)}, while all other Kraus operators are 𝒪(ϵ)\cal{O}(\epsilon) and will yield vanishing contributions in the limit. Now, Kαnν|ij|Kαnν=ϵ(1ϵ)(n1)(1ϵ)n(li+lj2)/2|ij|K_{\alpha n}^{\nu}|i\rangle\langle j|K_{\alpha n}^{\nu{\dagger}}=\epsilon(1-\epsilon)^{(n-1)}\left(1-\epsilon\right)^{n(l_{i}+l_{j}-2)/2}|i\rangle\langle j| when i,ji,j each have a 0 appearing at position α\alpha in their binary expansions. Summing over α\alpha gives a factor of mijm_{ij}, and then the sum over nn from 11 to ν\nu yields factor 2/(li+lj)+𝒪(ϵ)2/(l_{i}+l_{j})+\cal{O}(\epsilon). Therefore, the limiting channel also agrees with Eq. (42) and Eq. (43), and the proof is complete.\blacksquare

References

  • Bennett et al. [1993] C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, and W. K. Wootters. Teleporting an unknown quantum state via dual classical and Einstein-Podolsky-Rosen channels. Phys. Rev. Lett., 70:1895, 1993. doi: https://doi.org/10.1103/PhysRevLett.70.1895.
  • Bennett et al. [1996] C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher, J. A. Smolin, and W. K. Wootters. Purification of noisy entanglement and faithful teleportation via noisy channels. Phys. Rev. Lett., 76:722, 1996. doi: https://doi.org/10.1103/PhysRevLett.76.722.
  • Bennett et al. [1999] C. H. Bennett, D. P. DiVincenzo, C. A. Fuchs, T. Mor, E. Rains, P. W. Shor, J. A. Smolin, and W. K. Wootters. Quantum nonlocality without entanglement. Phys. Rev. A, 59:1070, 1999. doi: https://doi.org/10.1103/PhysRevA.59.1070.
  • Bergou [2007] János A. Bergou. Quantum state discrimination and selected applications. Journal of Physics: Conference Series, 84:012001, 2007.
  • Childs et al. [2013] Andrew M. Childs, Debbie Leung, Laura Manc̆inska, and Maris Ozols. A framework for bounding nonlocality of state discrimination. Communications in Mathematical Physics, 323:1121–1153, 2013. doi: https://doi.org/10.1007/s00220-013-1784-0.
  • Chitambar [2011] Eric Chitambar. Local quantum transformations requiring infinite rounds of classical communication. Phys. Rev. Lett., 107:190502, 2011. doi: https://doi.org/10.1103/PhysRevLett.107.190502.
  • Chitambar and Hsieh [2013] Eric Chitambar and Min-Hsiu Hsieh. Revisiting the optimal detection of quantum information. Phys. Rev. A, 88:020302(R), 2013. doi: https://doi.org/10.1103/PhysRevA.88.020302.
  • Chitambar and Hsieh [2014] Eric Chitambar and Min-Hsiu Hsieh. Asymptotic state discrimination and a strict hierarchy in distinguishability norms. Journal of Mathematical Physics, 55:112204, 2014. doi: https://doi.org/10.1063/1.4902027.
  • Chitambar et al. [2012a] Eric Chitambar, Wei Cui, and Hoi-Kwong Lo. Entanglement monotones for W-type states. Phys. Rev. A, 85:062316, 2012a. doi: https://doi.org/10.1103/PhysRevA.85.062316.
  • Chitambar et al. [2012b] Eric Chitambar, Wei Cui, and Hoi-Kwong Lo. Increasing entanglement monotones by separable operations. Phys. Rev. Lett., 108:240504, 2012b. doi: https://doi.org/10.1103/PhysRevLett.108.240504.
  • Chitambar et al. [2014] Eric Chitambar, Debbie Leung, Laura Manc̆inska, Maris Ozols, and Andreas Winter. Everything you always wanted to know about LOCC (but were afraid to ask). Commun. Math. Phys., 328:303, 2014. doi: https://doi.org/10.1007/s00220-014-1953-9.
  • Choi [1975] M-D Choi. Completely positive linear maps on complex matrices. Lin. Alg. Appl., 10:285, 1975. doi: https://doi.org/10.1016/0024-3795(75)90075-0.
  • Cirac et al. [1999] J. I. Cirac, A. K. Ekert, S. F. Huelga, and C. Macchiavello. Distributed quantum computation over noisy channels. Phys. Rev. A, 59:4249, 1999. doi: https://doi.org/10.1103/PhysRevA.59.4249.
  • Cohen [2011] S. M. Cohen. When a quantum measurement can be implemented locally, and when it cannot. Phys. Rev. A, 84:052322, 2011. doi: https://doi.org/10.1103/PhysRevA.84.052322.
  • Cohen [2016] S. M. Cohen. Devising local protocols for multipartite quantum measurements. Phys. Rev. A, 93:042334, 2016. doi: https://doi.org/10.1103/PhysRevA.93.042334.
  • Cohen [2017] S. M. Cohen. General approach to quantum channel impossibility by local operations and classical communication. Phys. Rev. Lett., 118:020501, 2017. doi: https://doi.org/10.1103/PhysRevLett.118.020501.
  • Cohen [2019] S. M. Cohen. Strong bounds on required resources for quantum channels by local operations and classical communication. Phys. Rev. A, 100:042324, 2019. doi: https://doi.org/10.1103/PhysRevA.100.042324.
  • Cohen [2022] S. M. Cohen. Local approximation of multipartite quantum measurements. Phys. Rev. A, 105:022207, 2022. doi: https://doi.org/10.1103/PhysRevA.105.022207.
  • Cohen [2023a] S. M. Cohen. Erratum: Local approximation for perfect discrimination of quantum states. Phys. Rev. A, 108:029902, 2023a. doi: https://doi.org/10.1103/PhysRevA.108.029902.
  • Cohen [2023b] S. M. Cohen. Local approximation for perfect discrimination of quantum states. Phys. Rev. A, 107:012401, 2023b. doi: https://doi.org/10.1103/PhysRevA.107.012401. See also Erratum [19].
  • Eckart and Young [1936] C. Eckart and G. Young. The approximation of one matrix by another of lower rank. Psychometrika, 1:211, 1936. doi: https://doi.org/10.1007/BF02288367.
  • Folland [1999] Gerald B. Folland, editor. Real Analysis: Modern Techniques and Their Applications. John Wiley and Sons, Inc., New York, 2nd edition, 1999.
  • Fortescue and Lo [2007] B. Fortescue and H.-K. Lo. Random bipartite entanglement from W and W-like states. Phys. Rev. Lett., 98:260501, 2007. doi: https://doi.org/10.1103/PhysRevLett.98.260501.
  • Fu et al. [2014] Honghao Fu, Debbie Leung, and Laura Mančinska. When the asymptotic limit offers no advantage in the local-operations-and-classical-communication paradigm. Phys. Rev. A, 89:052310, 2014. doi: https://doi.org/10.1103/PhysRevA.89.052310.
  • Jamiolkowski [1972] A Jamiolkowski. Linear transformations which preserve trace and positive semi-definiteness of operators. Rep. Math. Phys., 3:275, 1972. doi: https://doi.org/10.1016/0034-4877(72)90011-0.
  • Kleinmann et al. [2011] Matthias Kleinmann, Hermann Kampermann, and Dagmar Bruß. Asymptotically perfect discrimination in the local-operation-and-classical-communication paradigm. Phys. Rev. A, 84:042326, 2011. doi: https://doi.org/10.1103/PhysRevA.84.042326.
  • Kraus [1983] K. Kraus. States, Effects, and Operations. Spring-Verlag, Berlin, 1983.
  • Kretschmann et al. [2008a] Dennis Kretschmann, Dirk Schlingemann, and Reinhard F. Werner. The information-disturbance tradeoff and the continuity of Stinespring’s representation. IEEE Transactions on Information Theory, 54:1708, 2008a. doi: https://doi.org/10.1109/TIT.2008.917696.
  • Kretschmann et al. [2008b] Dennis Kretschmann, Dirk Schlingemann, and Reinhard F. Werner. A continuity theorem for Stinespring’s dilation. Journal of Functional Analysis, 255:1889, 2008b. doi: https://doi.org/10.1016/j.jfa.2008.07.023.
  • Li and Strang [2020] Chi-Kwong Li and Gilbert Strang. An elementary proof of Mirsky’s low rank approximation theorem. The Electronic Journal of Linear Algebra (ELA), 36:694, 2020.
  • Nechita et al. [2018] Ion Nechita, Zbigniew Puchała, Łukasz Pawela, and Karol Życzkowski. Almost all quantum channels are equidistant. J. Math. Phys., 59:052201, 2018. doi: https://doi.org/10.1063/1.5019322.
  • Nielsen [1999] M. A. Nielsen. Conditions for a class of entanglement transformations. Phys. Rev. Lett., 83:436, 1999. doi: https://doi.org/10.1103/PhysRevLett.83.436.
  • Peres and Wootters [1991] A. Peres and W. K. Wootters. Optimal detection of quantum information. Phys. Rev. Lett., 66:1119, 1991. doi: https://doi.org/10.1103/PhysRevLett.66.1119.
  • Spee et al. [2017] C. Spee, J. I. de Vicente, D. Sauerwein, and B. Kraus. Entangled pure state transformations via local operations assisted by finitely many rounds of classical communication. Phys. Rev. Lett., 118:040503, 2017. doi: https://doi.org/10.1103/PhysRevLett.118.040503.
  • Stinespring [1955] W. F. Stinespring. Positive functions on C*-algebras. Proc. Amer. Math. Soc., 6:211, 1955. doi: https://doi.org/10.1090/S0002-9939-1955-0069403-4.
  • Walgate et al. [2000] J. Walgate, A. J. Short, L. Hardy, and V. Vedral. Local distinguishability of multipartite orthogonal quantum states. Phys. Rev. Lett., 85:4972, 2000. doi: https://doi.org/10.1103/PhysRevLett.85.4972.
  • Wootters [1998] W. K. Wootters. Entanglement of formation of an arbitrary state of two qubits. Phys. Rev. Lett., 80:2245, 1998. doi: https://doi.org/10.1103/PhysRevLett.80.2245.