When a quantum system is distributed to spatially separated parties, it is natural to consider how the system evolves when the parties perform local quantum operations with classical communication (LOCC). However, the structure of LOCC channels is exceedingly complex, leaving many important physical problems unsolved. In this paper, we consider generalized resource theories of entanglement based on different relaxations to the class of LOCC. The behavior of various entanglement measures is studied under non-entangling channels, as well as the newly introduced classes of dually non-entangling and positive partial transpose (PPT)-preserving channels. In an effort to better understand the nature of LOCC bound entanglement, we study the problem of entanglement distillation in these generalized resource theories. We first show that unlike LOCC, general non-entangling maps can be superactivated, in the sense that two copies of the same non-entangling map can, nevertheless, be entangling. On the single-copy level, we demonstrate that every NPT entangled state can be converted into an LOCC-distillable state using channels that are both dually non-entangling and having a PPT Choi representation and that every state can be converted into an LOCC-distillable state using operations belonging to any family of polytopes that approximate LOCC. We then turn to the stochastic convertibility of multipartite pure states and show that any two states can be interconverted by any polytope approximation to the set of separable channels. Finally, as an analog to k-positive maps, we introduce and analyze the set of k-non-entangling channels.

In the distant laboratory paradigm of quantum information theory, a multipartite quantum system is distributed to spatially separated parties. The individual parties are restricted to performing only local quantum operations on their respective subsystems, but they are permitted to coordinate their actions through global classical communication. Quantum channels arising from this scenario are known as LOCC (local operations and classical communication), and they emerge as the natural operational class to consider for distributed quantum information processing tasks. Entanglement becomes a resource under the restriction of LOCC; it is consumed when performing information processing tasks (such as quantum teleportation), and entangled states cannot be generated using LOCC alone. It is, therefore, natural to study entanglement in the framework of quantum resource theory, which characterizes the convertibility of entangled states under LOCC and investigates their applicability in various information-theoretic tasks (see Refs. 1–3 for review).

Despite the conceptually intuitive description of LOCC in the distant-lab setting, mathematically characterizing the set of LOCC operations is notoriously complex.4 One way to overcome this difficulty is by relaxing the constraints of LOCC to encompass a more mathematically simple class of operations. Any task found to be impossible by these more general operations is, therefore, also impossible by LOCC. In fact, LOCC is not the largest class of channels that cannot create entanglement from unentangled, or separable, states. There exists a plethora of operational classes “beyond” LOCC that act invariantly on the set of separable states, and each of these defines a different quantum resource theory of entanglement manipulation.

It is interesting to study what differences arise in these generalized resource theories of entanglement compared to the standard LOCC scenario. Doing so has already proven to be successful in the asymptotic scenario where, under LOCC, certain mixed states demonstrate irreversibility in terms of pure-state entanglement distillation and the reverse process of entanglement dilution.5,6 The most dramatic example of this irreversibility is the phenomenon of bound entanglement,7 which refers to entangled mixed states that require the consumption of pure-state entanglement for their preparation, and yet no pure-state entanglement can be distilled back from them. Remarkably, the theory of entanglement distillation and dilution becomes much more amenable and elegant after enlarging the class of operations to include all asymptotically non-entangling operations. The work in Refs. 8–10 shows that bound entanglement no longer exists under this larger class of operations, and more importantly, the entanglement in any mixed state can be reversibly distilled and diluted at a rate given by the regularized relative entropy of entanglement.

In this paper, we analyze in detail the entanglement manipulation using operations beyond LOCC. Specifically, we study which state transformations are possible and which entanglement measures retain their validity within generalized resource theories of entanglement. For example, separable maps are those whose Choi representations are separable,11 while non-entangling maps are precisely those that do not produce an entangled output state for any separable input. The class of dually non-entangling maps, which we introduce here, consists of channels Λ such that both Λ and its dual map Λ* are non-entangling. All classes of operations that we consider contain LOCC as a proper subset.

One of our motivations is to better understand the nature of bound entanglement within the standard framework of LOCC. While it is known that no pure-state entanglement can be distilled from any state having a positive partial transpose (PPT), one of the major open problems in quantum information theory is to determine if the converse is also true;12 namely, does there exist non-PPT (NPT) bound entanglement? A strategy that would allow us to answer this question in the affirmative would be to identify a class of operations that contains LOCC, which lacks the ability to distill certain NPT states.

The first comprehensive attempt in this direction was carried out in Ref. 13, where it was shown that NPT bound entanglement does not exist when the class of operations is enlarged to contain all so-called PPT operations, as originally introduced in Refs. 14 and 15. This is the class of channels whose Choi representation is PPT (see also Sec. II for definitions). Even stronger, the results of Refs. 8–10 imply that no bound entanglement (not even PPT bound entanglement) exists for the class of all strictly non-entangling channels.

Beyond the question of distillability, we also examine the difference between the class of channels whose Choi operator is PPT and the class of PPT-preserving channels. While the former refers to the well-known class of PPT channels originally defined by Rains,14,15 to our knowledge, the latter has not been previously studied. Finally, analogous to the relation between completely positive (CP) and k-positive maps, we also introduce here the classes of k-non-entangling and completely non-entangling maps. The former refers to maps that are non-entangling when tensored with the identity on a k-dimensional ancilla, while the latter refers to maps that are k-non-entangling for all k.

The rest of the paper is structured as follows. Section II introduces the necessary background and notation, including precise definitions for the classes of channels that we study. Section III investigates how certain entanglement measures behave under different classes of non-entangling channels. The entanglement measures studied in this section include the robustness of entanglement, the Schmidt rank, the Rényi α-entropies of entanglement (as well as the relative Rényi entropies of entanglement), and the negativity. Although the robustness remains monotonic under non-entangling channels, we use the robustness to derive conditions for the convertibility of pure states under non-entangling channels that are independent of the celebrated majorization criterion for LOCC channels.16 We also show that the Rényi α-entropies of entanglement (which are entanglement measures for all α ∈ [0, +]) can be increased arbitrarily under non-entangling channels for α ∈ [0, 1/2). On the other hand, we show that the α-entropies of entanglement in the range α ∈ [1/2, +] coincide with a related relative Rényi entropy of entanglement measure, and therefore, they cannot increase under non-entangling channels. In particular, this gives a closed-form expression for the relative α-entropies of entanglement in certain cases. We also demonstrate that the Schmidt rank can be increased under the more restrictive class of dually non-entangling maps. Finally, we show that the negativity can be increased by an arbitrarily large fractional amount under PPT-preserving maps, thus highlighting the difference between the classes of PPT and PPT-preserving maps.

Entanglement distillation by classes of channels larger than LOCC is investigated in Sec. IV. We demonstrate that all entangled states can be converted to an LOCC-distillable state using dually non-entangling channels. The significance of this result is somewhat tempered by the fact that two copies of a non-entangling channel need not be non-entangling, something that we will also demonstrate in Sec. IV. Consequently, single-copy convertibility of a given state to an LOCC-distillable one by some generalized class of operations does not imply that the original state is distillable within generalized resource theory. Finally, we close the section by showing that any class of operations that generalizes LOCC by placing a finite number of linear constraints on the Choi matrix always allows us to obtain LOCC-distillable states from any quantum state.

In Sec. V, we relax the constraint of trace-preserving (TP) maps and turn to the problem of stochastic convertibility between pure states. For multipartite state spaces, a natural way to categorize entanglement is in terms of stochastic convertibility under LOCC. That is, two states |ψ⟩ and |ψ̂ belong to the same entanglement class if there is invertible LOCC transformation from |ψ⟩ to |ψ̂ that succeeds with some nonzero probability.17 We show that this entanglement structure completely collapses under any class of entangling undetected channels. Connections to the problem of tensor rank calculation are also discussed.

The structure of k-non-entangling maps and k-PPT-preserving maps are analyzed in Sec. VI. We show that complete non-entangling is equivalent to d-non-entangling when the local systems have dimension d, and we show that the structure of k-non-entangling maps is nontrivial. That is, for every k < d, there exists k-non-entangling maps that are not (k + 1)-non-entangling. We also investigate distillability under this class of maps, but we have not been able to answer whether or not they allow for bound entanglement. In particular, we show that previously used constructions are not 3-non-entangling, so the weaker question of whether every entangled state can be converted to an LOCC-distillable one under this class of maps remains open. Finally, a detailed Proof of Theorem 4 is presented in  Appendix A, while a concluding discussion is given in Sec. VII.

Let L(H) denote the space of linear operators on a finite-dimensional Hilbert space H with orthonormal basis denoted by {|1⟩, …, |n⟩} (where n is the dimension of the space). An operator ρL(H) is in the set D(H) of density operators (or states) on H if ρ ≥ 0 and Tr(ρ) = 1. The tensor product of two Hilbert spaces HA and HB is denoted by HAHB or HAB. The identity operator on a Hilbert space HA is denoted as 1A, while the identity map on L(HA) is denoted as idA, and the identity operator on CdCd is denoted by 1d2.

We write ψ to denote the density operator ψ = |ψ⟩ ⟨ψ| associated with a pure state vector |ψH. For bipartite pure states |ψHAB, we typically assume, without loss of generality, that |ψ⟩ is in the Schmidt form

where |ii⟩ is the shorthand for |i⟩ ⊗ |i⟩ and λ = (λ1, …, λd) are the Schmidt coefficients of ψ (in decreasing order) such that i=1dλi=1 and dmindimHA,dimHB. The Schmidt rank of the state is defined to be the number of non-zero entries in λ, which we denote by rank(λ). A positive operator σL(HAB) is said to be separable with respect to A : B if it can be written as σ = ∑iτiωi for some positive operators τiL(HA) and ωiL(HB). The set of separable density operators on HAB will be denoted as S(HA:B), or simply S if the spaces A and B are clear from context.

The partial transpose of an operator XL(HAB), denoted as XΓA, is the operator that results from applying the transposition map AAT (with respect to a given basis of HA) to system A. We simply write XΓ if the system being transposed is clear from context. A positive operator ρL(HAB) is positive under partial transpose (PPT) with respect to A : B if ρΓ ≥ 0, and it is otherwise said to be non-PPT (or NPT) (we note that the definition of PPT is independent of the basis chosen for the respective transposition). It is well-known that every separable state is PPT, with the converse holding true only in dimensions 2 ⊗ 2 and 2 ⊗ 3.18 

For any integer d ≥ 2, the maximally entangled pure state on CdCd is denoted as

and we write ϕd+=|ϕd+ϕd+| to denote the corresponding density operator. We denote the so-called “flip” operator on CdCd by

which is precisely the partial transpose of ϕd+ scaled by d. Analogously, we write ϕA+ to denote the (non-normalized) maximally entangled operator on the Hilbert space HAHA,

where {|i⟩} is an orthonormal basis for HA. Note that we can write the maximally entangled operator of the four-partite space ϕAB+L(HABHAB) as

Given a linear map Λ:L(HA1)L(HA2), its Choi representation is the linear operator J(Λ)L(HA2A1) defined by

For any XL(HA1), the action of the map Λ on X can be written in terms of its Choi representation as Λ(X)=TrA1(J(Λ)(1A2XT)). Recall that Λ is completely positive (CP) if and only if its Choi representation is a positive semi-definite operator, and Λ is trace preserving (TP) if and only if it holds that TrA2(J(Λ))=1A1. A linear map that is both completely positive and trace preserving (CPTP) is called a quantum channel. Given two density operators ρ1,ρ2D(H) and an operator AL(H), the linear map Λ:L(H)L(H) defined by

(1)

is CPTP as long as 0 ≤ A1A. Many of the maps constructed in this work take this form. Given a linear map Λ:L(H)L(H), its dual map Λ*:L(H)L(H) is the unique linear map such that Tr(Λ(X)Y) = Tr(XΛ*(Y)) holds for all operators XL(H) and YL(H). It follows readily from the definition that the dual of a map of the form in (1) is given by

while the Choi representation of a map of the form in (1) is given by

We now define all the classes of operations that will be explored in this paper. We first recall the definitions of separable and PPT maps. Let Λ:L(HA1B1)L(HA2B2) be a completely positive linear map. Then, Λ is said to be separable if its Choi representation J(Λ)L(HA2B2HA1B1) is separable with respect to A2A1 : B2B1. To keep the bipartition clear, we will often denote HA2B2HA1B1 as HA2A1:B2B1. Note that Λ is separable if and only if it can be written as

for some CP maps Φj:L(HA1)L(HA2) and Ψj:L(HB1)L(HB2). Analogously, the map Λ is a PPT map if the partial transpose of its Choi representation is a positive semi-definite, i.e., if J(Λ)Γ=J(Λ)ΓA2A10. It holds that Λ is PPT if and only if the map ΛΓ defined by

is also completely positive. In particular, note that JΓ) = J(Λ)Γ. This definition of PPT maps is due to Rains.14 

In this paper, we are primarily concerned with the non-entangling maps, which are defined explicitly as follows. We provide here an analogous definition for PPT-preserving maps.

Definition 2.1.

Let Λ:L(HA1B1)L(HA2B2) be a completely positive map. Then, Λ is separable-preserving (or non-entangling) if Λ(σ) is separable (with respect to A2 : B2) for any operator σ that is separable (with respect to A1 : B1). Analogously, Λ is PPT-preserving if Λ(σ) is PPT for any PPT operator σ.

By definition, the non-entangling channels comprise the largest class of channels which cannot create entanglement, while the PPT-preserving channels are exactly those that cannot generate an NPT state from PPT input states. It is clear that every separable (respectively PPT) map is separable-preserving (respectively PPT-preserving). We remark that the classes of separable-preserving and PPT-preserving maps are distinct from the classes separable and PPT maps. Indeed, the bipartite swap operator is clearly both separable- and PPT-preserving, but its corresponding Choi representation is proportional to the maximally entangled state which is neither separable nor PPT.

We will also consider the following set of maps that is even more restricted than the non-entangling maps which nonetheless contains the set of LOCC channels.

Definition 2.2.

A completely positive map Λ:L(HA1B1)L(HA2B2) is said to be dually non-entangling if both Λ and its dual map Λ* are separable-preserving.

Note that a map is completely positive if and only if its dual map is, but the dual of a trace-preserving map is not necessarily trace-preserving. Dually non-entangling channels form a subset of the non-entangling channels which still contains all LOCC. Indeed, the dual of any separable map is clearly separable; hence, separable maps are dually non-entangling (and thus, LOCC channels are as well).

Any quantum channel is a linear map that must not only preserve positivity but also preserve positivity when it acts on the systems tensored with any ancilla space for a linear map to represent a physical manipulation of quantum states. Motivated by this requirement of complete positivity for a quantum channel, in the following, we define a notion of complete non-entangling. The analogous definition is provided for PPT-preserving maps. We note that completely non-entangling maps (respectively, completely PPT-preserving maps) are exactly those maps that are separable (respectively, PPT). This will be discussed and proved in Sec. VI.

Definition 2.3.

Let Λ:L(HA1B1)L(HA2B2) be a completely positive map and let k ≥ 1 be an integer. Then, Λ is k-separable-preserving (or k-non-entangling) if, for all k-dimensional systems HA=HB=Ck, the map Λ ⊗ idAB is separable-preserving (with respect to A1A : B1B and A2A : B2B). If Λ is k-non-entangling for all integers k, then Λ is said to be completely non-entangling.

Definition 2.4.

Let Λ:L(HA1B1)L(HA2B2) be a completely positive map and let k ≥ 1 be an integer. Then, Λ is k-PPT-preserving if, for all k-dimensional systems HA and HB, the map Λ ⊗ idAB is PPT-preserving. If Λ is k-PPT-preserving for all integers k, then it is said to be completely PPT-preserving.

For the final family of channels that we will use, we recall the notion of entanglement witnesses. A Hermitian operator WL(HAB) is an entanglement witness with respect to A : B if it holds that Tr() ≥ 0 for all separable states ρS(HA:B). Since the Choi operator for any separable map is separable, we can use entanglement witnesses on the level of Choi operators to detect non-separable maps. Even stronger, it is known that entanglement witnesses provide a full characterization of separability in the following sense. A positive operator X is separable if and only if Tr(WX) ≥ 0 holds for all entanglement witnesses W.18 

Geometrically, every entanglement witness can be seen as a hyperplane in the space of Hermitian operators that separates the positive convex cone of separable operators (see Fig. 1). Furthermore, the set of separable operators is not a polytope,19,20 and therefore, any finite set of entanglement witnesses is insufficient to determine separability. Therefore, one obtains a strict relaxation of the separability constraint by requiring that Tr(WX) ≥ 0 hold for not all witnesses but instead for just some finite subset of witnesses. This idea motivates the following type of operations.

FIG. 1.

Given a finite number of entanglement witnesses {Wi}i=1n on L(HA2A1:B2B1), the set of entangling undetected channels have Choi representations that are not witnessed by any Wi. This polytope clearly contains all separable channels.

FIG. 1.

Given a finite number of entanglement witnesses {Wi}i=1n on L(HA2A1:B2B1), the set of entangling undetected channels have Choi representations that are not witnessed by any Wi. This polytope clearly contains all separable channels.

Close modal

Definition 2.5.

Let {Wi}i=1nL(HA2B2A1B1) be a finite collection of entanglement witnesses with respect to A2A1 : B2B1. A completely positive map Λ:L(HA1B1)L(HA2B2) is said to be entangling undetected by {Wi}i=1n if it holds that Tr(WiJ(Λ)) ≥ 0 for all i ∈{1, …, n}.

Note that for any convex set in Rd, a containing polytope can be constructed that approximates the set to arbitrary precision.21 Thus, for any ϵ > 0, there exists a collection of entanglement witnesses {Wi}i=1n such that a non-separable Λ is entangling undetected by {Wi}i=1n only if ∥J(Λ) − Ω∥1 < ϵ for all operators ΩL(HA2B1A1B1) that are separable with respect to A2A1 : B2B1, where ∥·∥1 denotes the trace norm.

In this section, we explore what kinds of transformations among states are possible under non-entangling channels but not possible under LOCC. This is done by finding transformations under non-entangling channels that increase some entanglement measure. We examine which entanglement measures retain their monotonicity under non-entangling channels, dually non-entangling channels, and PPT-preserving channels.

We first explore certain entanglement measures that can be increased with non-entangling channels and show how certain other measures can yield conditions for conversion under such maps. Although most of the claims in this section follow from the already known results, it is useful for our purposes to collect information about which entanglement measures remain monotonic under all non-entangling maps and which ones do not.

Consider first the relative entropy of entanglement,22,23

where S(ρσ) ≔ Tr(ρ log ρ) − Tr(ρ log σ) is the relative entropy of ρ and σ. This is an entanglement measure that is clearly also monotonic under non-entangling channels, a fact that follows directly from the joint monotonicity of the relative entropy, i.e., S(Λ(ρ)‖Λ(σ)) ≤ S(ρσ) for every CPTP map Λ.24 Another entanglement measure that retains its monotonicity under non-entangling channels, as already observed in Ref. 9, is the robustness of entanglement,25 which is defined as

(2)

where, for any separable state σS,

(3)

Indeed, if Λ is a non-entangling channel, it is clearly seen that R(Λ(ρ)‖Λ(σ)) ≤ R(ρσ) holds for any separable σ. For pure states, the robustness reduces (see Ref. 25) to

(4)

In particular, it holds that R(ϕk+)=k1 for any integer k ≥ 2.

The robustness yields necessary and sufficient conditions for certain transformations among certain states under non-entangling channels. The following result provides necessary and sufficient conditions for the transformation of a maximally entangled state to any other state purely in terms of the robustness. This is effectively a statement about the exact one-shot entanglement cost under non-entangling channels and, thus, follows straightforwardly from Theorem 2 in Ref. 10 (similar methods are also used in Ref. 9). For the sake of completeness, we provide a proof of this fact as it presents techniques and ideas that will be used later in this work.

Theorem 1

(Ref. 10). For every bipartite stateρand every positive integerk, there exists a non-entangling channel Λ such thatΛ(ϕk+)=ρif and only ifR(ρ)R(ϕk+).

Proof.
Since we have already pointed out that R cannot increase under non-entangling channels, it remains to explicitly construct a channel achieving the transformation whenever R(ρ)R(ϕk+). This is done by the channel defined as
(5)
where σ is any separable state for which R(ρ) = R(ρσ). To see that this channel is, indeed, non-entangling, note that for all other separable states σS, one has that Tr(ϕk+σ)1/k and, thus, there exists a value p ≤ 1/k satisfying Λ(σ′) = + (1 − p)σ. However, one has that 1/k=1/(1+R(ϕk+)) and that pρ+(1p)σS for any p ≤ 1/(1 + R(ρ)) by construction. Hence, the above channel is non-entangling whenever R(ρ)R(ϕk+).□

Whereas pure state transformation under LOCC is governed strictly by majorization of Schmidt coefficients,16 the result of Theorem 1 indicates that in certain cases, the Schmidt rank of a pure state (which cannot be increased under LOCC or separable channels26) can be increased under non-entangling channels. Indeed, for any positive integers k and d with k < d, there exist pure states with the Schmidt rank k that can be transformed to states with the Schmidt rank d. To see this, note that there exist pure states ψ with the Schmidt rank d but with robustness R(ψ) < 1/k. As R(ϕk+)=1/k, Theorem 1 implies that ϕk+ may be transformed into such pure states ψ by non-entangling operations. This fact is formalized in the following corollary.

Corollary 2.

For all integerskanddwith 2 ≤ k < d, there exists a non-entangling channel Λ such thatΛ(ϕk+)is a pure state with the Schmidt rankd.

The above result demonstrates that the Schmidt rank is not monotonic under non-entangling channels. It is natural to ask then which other entanglement measures remain monotonic under non-entangling channels and which ones do not. We now investigate other entanglement measures that can be increased under non-entangling channels. The Rényi entropy of entanglement of order α is defined in terms of the Schmidt coefficients of pure states as

for any α ∈ [0, +], where Hα are the Rényi entropies

These quantities are the Schur concave as functions of the Schmidt coefficients. Hence, the function Eα is an entanglement measure for pure states for any α. Note that the Schmidt rank of a pure state is in direct correspondence with the Rényi entropy of order zero. Since we have seen that E0 can be increased under non-entangling channels, it is natural to consider whether Eα is not monotonic under non-entangling channels for other values of α as well. The Rényi entropy of order 1 reduces to the standard entropy of entanglement

which is known to be monotonic under non-entangling channels since it coincides with the relative entropy of entanglement ER for pure states.23 Furthermore, note from (4) that the Rényi entropy of entanglement of order 1/2 is a monotonic function of the robustness since

so E1/2 is also monotonic under non-entangling channels by monotonicity of R.

In the following, we show that Eα can, in fact, be increased under non-entangling channels for any α ∈ [0, 1/2), while it cannot be increased under non-entangling channels for any α ≥ 1/2. The proof of the former is actually a simple corollary of Theorem 1. To show this, it suffices to find a pure state ψ with E1/2(ψ)=E1/2(ϕ2+) but that Eα(ψ)>Eα(ϕ2+) for α < 1/2. Since E1/2(ϕ2+)=1, Theorem 1 indicates that the transformation ϕ2+ψ is possible under non-entangling channels as long as E1/2(ψ) ≤ 1.

Corollary 3.

The measuresEαcan be increased by non-entangling channels for eachα ∈ [0, 1/2).

Proof.
Let ψ be a pure state with the Schmidt rank greater than 2 such that E1/2(ψ) = 1. States with this property certainly exist since we may take, for instance, the state with Schmidt coefficients λ = (1 − 2ϵ, ϵ, ϵ) and ϵ = 1/18. Note that Eα(ψ) is a continuous and non-increasing function of α and that
for this choice of ψ. It follows that Eα(ψ) > 1 for all α < 1/2. Note, however, that Eα(ϕ2+)=1 holds for all α. In particular, E1/2(ϕ2+)=1 so the transformation ϕ2+ψ is possible with non-entangling channels by Theorem 9. This completes the proof.□

To prove that Eα is monotonic under non-entangling channels for all α ≥ 1/2, it suffices to show that the value of Eα coincides with that of another entanglement measure that is known to be monotonic under non-entangling channels. For every α ∈ [0, +) and density operators ρ and σ, the Rényiα-relative entropy of ρ with respect to σ is defined as

where supp denotes the support of an operator. Note that S1 = S is the standard relative entropy. The Rényiα-relative entropy of entanglement of a state ρ is

where the minimization is taken over all separable states σ. The Rényi α-relative entropies are known to be contractive under CPTP maps for all α ∈ [0, 2], i.e.,

for all CPTP maps Λ.27 It readily follows that ER,α is monotonic under non-entangling channels for α ∈ [0, 2]. Further properties of ER,α have been studied in Ref. 28. In the case α = 1, it is known that ER,1 reduces to the well-known relative entropy of entanglement ER. Remarkably, for all pure states ψ, one has that ER,α(ψ) = E1/α(ψ) for all α ∈ [0, 2] (see Corollary 2 in Ref. 28).

Theorem 4

(Ref. 28). For pure statesψ, it holds thatER,α(ψ) = E1/α(ψ) for allα ∈ [0, 2].

The result of Theorem 4 has been recently proved in Ref. 28 by relating the resource theories of coherence and entanglement. However, we have independently obtained an alternative proof of this fact that uses similar methods to those in Ref. 23, where it is proven that E(ψ) = ER(ψ). We include our Proof of Theorem 4 in  Appendix A for the interested reader.

Corollary 5.

For allα ∈ [1/2, +], the measuresEαare monotonic for pure states under non-entangling channels.

While majorization provides necessary and sufficient conditions for conversions of pure states under LOCC, whether one can find similar necessary and sufficient conditions for convertibility of arbitrary pure states under non-entangling channels remains unknown. In the following, we provide a sufficient condition for convertibility under non-entangling channels of arbitrary pure states that is independent of the majorization criterion. We first provide the following Lemma (which was first proved in Ref. 29).

Lemma 6.

For every pure stateψ, it holds thatmaxσSTr(ψσ)=λ1, whereλ1is the largest Schmidt coefficient ofψ.

Proof.
We may suppose without loss of generality that ψ is in the Schmidt form |ψ=i=1nλi|ii. We first show that Tr(ψσ) ≤ λ1 holds for all σS. It will suffice to show that Tr(ψϕ) ≤ λ1 for all separable pure states |ϕ⟩ = ∑i,juivj|ij⟩ with ∑i|ui|2 = ∑j|vj|2 = 1. Note that Tr(ψϕ) = |⟨ψ|ϕ⟩|2 and
as desired. Since |11⟩ is separable and |⟨ψ|11⟩|2 = λ1, the maximum is achieved.□

Theorem 7.

Letψandϕbe pure states with Schmidt coefficientsλandμ, respectively. If it holds that 1 + R(ϕ) ≤ 1/λ1, then there exists a non-entangling channel Λ such that Λ(ψ) = ϕ.

Proof.
The proof is very similar to that of Theorem 1. The desired channel is given by
(6)
where σS is a separable state satisfying R(ϕ) = R(ϕσ). It is evident that Λ(ψ) = ϕ, so it remains to show that Λ is non-entangling. By construction, it holds that pϕ+(1p)σS whenever p ≤ 1/(1 + R(ϕ)). By Lemma 6, it holds that Tr(ψρ) ≤ λ1 for any separable state ρ. If the condition that 1 + R(ϕ) ≤ 1/λ1 is met, it follows that Λ(ρ)S whenever ρ is separable and, thus, Λ is non-entangling.□

We now explore what transformations are possible under the more restrictive class of dually non-entangling channels. Recall that these are the channels Λ such that both it and its dual Λ* do not generate entanglement (although Λ* is not necessarily a channel). In particular, we show that the Schmidt rank can be increased by this more restrictive class of transformations.

We first recall some conditions for certain operators to be separable. For a pure state ψ with Schmidt coefficients λ = (λ1, …, λd) (in decreasing order), it holds that R(ψ1d21d2)=d2λ1λ2 (see Appendix B of Ref. 25). Hence, operators of the form ψ+s1d2 for s ≥ 0 are separable if and only if sλ1λ2. Additionally, for any Hermitian operator Δ on a bipartite space with local dimensions d1 and d2, the operator 1d1d2+Δ is always separable whenever ‖Δ∥ ≤ 1, where ∥·∥ is the Frobenius norm (see Theorem 1 of Ref. 30). We now state a useful condition for certain types of channels to be dually non-entangling.

Lemma 8.
Letψandϕbe pure states with Schmidt coefficients (in decreasing order)λandμ, respectively, and consider the CPTP mapΛ:L(CkCk)L(CdCd)defined by
(7)
If it holds that
(8)
then Λ is non-entangling. If it further holds that
(9)
then Λ is also dually non-entangling.

Proof.
Outputs of the channel in (7) are of the form Λ(ρ)=pϕ+(1p)1d2/d2 for any state ρ, where p = Tr(ψρ). From the observation above, this is separable if and only if p1/(1+d2μ1μ2) since ∥ϕ∥ = 1. Since Tr(ψρ) ≤ λ1 for any separable state ρ (see Lemma 6), we obtain the desired result that Λ is non-entangling whenever the condition in (8) holds. On the other hand, for any density operator ρ, its output under the dual map of Λ is given by
Let ρ be a separable state. If it holds that Tr(ϕρ) ≤ 1/d2, then the operator (d2Tr(ϕρ)1)ψ+1k2 is separable by Theorem 1 of30 since ∥ψ∥ = 1 and |d2Tr(ϕρ) − 1| ≤ 1. We may, therefore, assume that Tr(ϕρ) > 1/d2, in which case Λ*(ρ) is separable if and only if d2Tr(ϕρ)1+1/λ1λ2. Since Tr(ϕρ) ≤ μ1 holds for any separable state ρ, it follows that Λ*(ρ) is separable for all separable states ρ, as long as the condition in (9) holds.□

Now that we have conditions for a channel of the form in (7) to be dually non-entangling, we may construct such dually non-entangling channels that increase the Schmidt rank of pure states arbitrarily. The trick is to construct pure states ψ and ϕ with Schmidt coefficients λ and μ that satisfy both (8) and (9). This is done in the Proof of Theorem 9.

Theorem 9.

For all integersd, k ≥ 2, there exists pure statesψandϕwith Schmidt ranks equal tokandd, respectively, and a dually non-entangling channel Λ such that Λ(ψ) = ϕ.

Proof.
Let d, k ≥ 2 be integers and choose positive real numbers δ and ϵ small enough so that
(10)
(11)
are both satisfied. This can be done since, for example, the values δ = d−4 and ϵ = d−12 satisfy both (10) and (11) (see  Appendix B for proof). Let ψ and ϕ be pure states with Schmidt coefficients λ and μ given by
(12)
(13)
such that ψ and ϕ have the Schmidt rank k and d, respectively. The desired channel is given by
(14)
which is dually non-entangling by Lemma 8 and performs the transformation Λ(ψ) = ϕ.□

There is nothing remarkable about the statement in Theorem 9 in the case when kd. However, if one chooses 2 ≤ k < d, Theorem 9 indicates that it is possible to increase the Schmidt rank of a pure state from k to d by a non-entangling channel. We remark that it does not seem possible to strengthen the above results by considering maps whose outputs are mixtures of ϕ and some separable state σ other than the identity, as we did in Theorems 1 and 7 for the case of non-entangling channels. The reason is that σ can be rank-deficient30 so that there exists some other separable state ρS such that Tr(σρ) = 0. With this state, the output of the dual map Λ*(ρ) is entangled and thus, Λ cannot be dually non-entangling.

Note that a completely positive map Λ is separable if and only if its dual map Λ* is separable. As we have seen here, the same is not necessarily true of separability preserving maps, a fact that lead us to explore the class of maps that are dually non-entangling. Similarly, a map Λ is PPT if and only if its dual map is PPT. Indeed, for any map Λ, it is clear that (Λ*)Γ=(ΛΓ)*. Since a map is completely positive if and only if its dual is, it follows that ΛΓ is completely positive if and only if (Λ*)Γ is. The condition that a map is PPT-preserving, however, is not equivalent to the condition that its dual map be PPT-preserving. The class of “dually PPT-preserving” maps may also be studied, in principle, but are not investigated in this paper.

We now investigate the difference between the class of PPT maps (whose Choi representation is PPT) defined by Rains14 and the class of PPT-preserving maps defined in this paper. Despite some confusion in the literature,2,31 these classes are not equivalent49 (this distinction goes away, however, if we require that a map be PPT-preserving when tensored with the identity on any ancilla space). We have seen that these two classes of maps are, indeed, distinct as, for example, the swap operation is clearly PPT-preserving but its Choi representation is not PPT. To further distinguish between these classes of maps, we present in this section the examples of transformations under PPT-preserving channels that are not possible under PPT channels. This will be done by investigating the negativity, a quantity defined as

for bipartite states ρ. It is well-known that the negativity is monotonic under PPT channels (and thus, under LOCC as well), so the negativity is an entanglement measure. Here, however, we show that the negativity can be increased by PPT-preserving channels. This indicates that the class of PPT-preserving channels is far more powerful for entanglement manipulation than the class of channels whose Choi representations are PPT.

Before proceeding, we first note that states of the form a1d2/d2+(1a)ϕd+ are PPT if and only if ad/(d + 1). Indeed, note that

which is positive exactly when ad/(d + 1) since the eigenvalues of Fd are ±1. We now show that the negativity can be increased by PPT-preserving channels.

Theorem 10.

The negativity is not a monotone under PPT-preserving channels.

Proof.
Let d > 3 be an integer. We construct a PPT-preserving channel Λ satisfying N(Λ(ρ)) > N(ρ) for some bipartite state ρD(CdCd). Consider the operator
and the corresponding channel Λ:L(CdCd)L(CdCd) defined by
(15)
By the observation above, it holds that Λ(ρ) is PPT whenever Tr() ≥ d/(d + 1), which certainly holds for any PPT state ρ. We conclude that Λ is PPT-preserving. Now consider the state
which is the most entangled Werner state and has negativity equal to N(ρ) = 1/d. However,
which has negativity
As d > 3, it follows that N(Λ(ρ)) > N(ρ) for this state.□

Note that for ρ and Λ as in the proof above, we have N(Λ(ρ))/N(ρ) = (d − 1)/2. Therefore, this ratio can be made arbitrarily big by taking a large enough d. Moreover, a PPT-preserving channel that is not PPT can be constructed from any NPT state ρ and a corresponding pure state vector |ψ⟩ corresponding to a negative eigenvalue of ρΓ. Indeed, for such ρ and |ψ⟩, define the operator

The corresponding channel Λ of the form in (15) is clearly PPT-preserving since it maps every state to a PPT state, which results from the fact that Tr() ≥ d/(d + 1) for every PPT state σ. Nonetheless, this channel is not PPT since the operator J(Λ)Γ is not a positive semi-definite. Indeed, note that (1d2Fd)ρT0 but that

since Tr(Γ) < d/(d + 1).

In order to explore the possible existence of NPT bound entanglement under LOCC, we address the question of distillable entanglement under classes of channels that are strictly larger than LOCC. As mentioned in the Introduction, it has been shown that all NPT states are distillable by PPT channels,13 but non-entangling channels are independent from this class of maps. Reference 9 has shown that all entangled states are distillable not only under asymptotically non-entangling channels but also under non-entangling channels.

A state ρ is said to be distillable under some class of channels if there exists a sequence of channels {Λn} in this class such that limnΛn(ρn)ϕ2+1=0, where ∥·∥1 denotes the trace norm. It has been shown that every entangled two-qubit state is distillable by LOCC channels.32 It is tempting to argue that if some class of operations containing LOCC can convert ρ to a two-qubit entangled state, then ρ is distillable using these operations. However, this argument only holds if the class of operations is closed under tensor products. More precisely, one requires that Λn belongs to the class whenever Λ does. Indeed, if Λ is in a class of operations that is closed under tensor products and it is the case that Λ(ρ) is a two-qubit entangled state, then one can use this class of operations to generate the n-copy state Λ(ρ)n=Λn(ρn) for any n. For large enough n, this can further be converted to ϕ2+ with arbitrary precision using the LOCC distillation protocol of Ref. 32. It is true that Λn is LOCC or PPT whenever Λ is LOCC or PPT. However, for other classes of mappings, this will not be true (see below). We must, therefore, distinguish between the ability of a class of operations to transform ρ into an LOCC-distillable state on the single-copy level vs its ability to distill ρ into a maximally entangled state in the asymptotic limit.

Although obtaining an LOCC-distillable state is weaker than obtaining the desired state ϕ2+, we believe that studying this question is operationally justified. LOCC operations are generally regarded as being less resource intensive since they do not require the transmission of quantum information. Thus, our approach can be seen as partitioning the distillation of ϕ2+ into two steps: the LOCC part and the non-LOCC part. In this way, we can minimize the quantum communication resources needed to perform the distillation by completing as much of the task as possible using LOCC. An analogous rationale is applied in other quantum resource theories such as magic-state quantum computation, in which a quantum computation is divided into two parts: the generation of the so-called “magic states” and then the manipulation of these states using much simpler stabilizer operations.33 

It is easy to see that not all non-entangling maps are completely non-entangling. A simple example is given by the flip map Λ(⋅) = Fd ⋅ Fd (see also Sec. VI). However, this does not forbid in principle that Λn belongs to the class of non-entangling maps whenever Λ does. We start this section by answering this question providing an explicit example of a non-entangling map that becomes entangling when taken in a two-copy form.

Theorem 11.

Non-entangling maps can be superactivated. That is, there exists non-entangling maps Λ such that Λ⊗2is not non-entangling.

Proof.
We provide an explicit example. Consider the map Λ:L(C2C2)L(C2C2) given by
(16)
where ψ2±=(|01±|10)/2 are elements of the Bell basis. This map is non-entangling since Tr[ϕ2+X]12 for all separable X, and states of the form pψ2++(1p)12(ψ2++ψ2) are separable whenever p ≤ 1/2 as an entangled Bell-diagonal state requires one eigenvalue greater than 1/2.34,35 However, the map Λ⊗2 is entangling. To see this, consider the so-called Smolin state,36 
(17)
Despite its appearance, this state is separable in the cut AA′ : BB′.36 The action of Λ⊗2 on this state is
(18)
Taking a partial transpose on system BB′ yields
(19)
Note that F2 has an eigenvalue −1 with eigenvector |ψ2. Hence,
(20)
which implies that Λ⊗2(ρ) is entangled.□

Given the above result, we now turn to the problem of obtaining LOCC-distillable states on the single-copy level. Our techniques will involve maps whose outputs are always mixtures of some operator with the maximally entangled state of two qubits. Note that the state p14/4+(1p)ϕ2+ is separable if and only if p ≥ 2/3.37 We first show how any entangled state can be transformed to an LOCC-distillable state using dually non-entangling channels.

Lemma 12.
LetWbe an entanglement witnessWsuch thatW2 ≤ 1 (where ∥·∥2is the Frobenius norm), and letA = (W + 21)/3. Then, the channel Λ defined by
(21)
is dually non-entangling.

Proof.
First note that ∥W2 ≤ 1 implies that −1W1 and, thus, 0 ≤ A1, so Λ is CPTP. Since W is an entanglement witness, it holds that Tr() ≥ 0 for all separable states σ. Hence, Tr() ≥ 2/3 holds for any separable state σ, and therefore, as discussed above, the state Λ(σ) is separable whenever σ is separable. Thus, Λ is non-entangling. To conclude that Λ is also dually non-entangling, we show that Λ*(σ) is separable for all two-qubit states σ. Note that
and that 0Tr(ϕ2+σ) holds for any two-qubit state σ. Since
and ∥W2 ≤ 1, it holds that Λ*(σ) is separable by Theorem 1 of Ref. 30.□

Theorem 13.

All entangled states may be converted by dually non-entangling channels to a state that is distillable under LOCC.

Proof.

Let ρ be an arbitrary entangled state. There exists an entanglement witness W such that Tr() < 0 but that Tr() ≥ 0 for all separable states σS. We may suppose without loss of generality ∥W2 ≤ 1. Indeed, we may otherwise take the witness W′ = W/∥W2. Consider the channel Λ in (21) with A = (W + 21)/3, which is dually non-entangling by Lemma 12 and satisfies Tr(Λ(ρ)ϕ2+)>1/2. It follows from Lemma 6 that Λ(ρ) is an entangled state of two qubits and is, therefore, distillable under LOCC.□

Finally, while it is known that all non-PPT entangled states are distillable under PPT channels, we next show that all such states can be converted to a state that is distillable under LOCC by only using channels in the more restrictive class of maps that are both PPT and dually non-entangling.

Theorem 14.

All NPT states may be converted by channels that are both dually non-entangling and PPT to a state that is distillable under LOCC.

Proof.
Let ρ denote an arbitrary d1 × d2 NPT state and let |η⟩ denote a normalized eigenvector of ρΓ corresponding to a negative eigenvalue. Then, W = |η⟩ ⟨η|Γ is an entanglement witness that detects ρ since Tr() < 0 but Tr() ≥ 0 for all separable states σ. Furthermore, note that the Frobenius norm of this witness is ∥W2 =  ∥|η⟩ ⟨η|∥2 = 1. It follows from Lemma 12 that the channel Λ in (21) is dually non-entangling. To show that this channel is also PPT, note that the Choi representation of this map is
where A=(W+21d1d2)/3. The partial transpose of this operator is
where F2=2(ϕ2+)Γ is the flip operator. This is clearly positive since −14F2. Finally, note that Λ(ρ) is an entangled two-qubit state, so it is distillable under LOCC.□

Finally, we show that all states (even separable ones) can be converted to an LOCC-distillable state under any class of entangling undetected channels.

Theorem 15.

Let{Wi}i=1nL(HA2B2A2B1)be a finite number of entanglement witnesses with respect toA2A1 : B2B1. Every state can be converted to an LOCC-distillable state using channels that are entangling undetected by{Wi}i=1n.

Proof.
Let Π=i,j=01|ijij| be the projector onto an arbitrary two-qubit subspace of HA2B2. We first note that the projected partial traces of the witnesses ViΠTrA1B1(Wi)Π are again entanglement witnesses for two-qubit states with support on the span of {|ij}i,j=01HA2:B2. Indeed, let σ be a separable state with support on span{|ij}i,j=01. Then, σ1A1B1 is separable with respect to A2A1 : B2B1, and thus,
Next consider the set of two-qubit density operators which are not detected by Vi. This is the intersection of a convex set (D, the set of all two-qubit density operators) with a polyhedron (P, the set of all Hermitian operators satisfying Tr(ViX) ≥ 0 for all i = 1, …n). If D is contained entirely in the polyhedron, then there obviously exists an entangled two-qubit state ρ not detected by any of Vi. If D is not contained entirely in the polyhedron, then DP has a facet.38 However, the set of two-qubit separable states does not have a facet,20 and thus, an entangled state ρ exists in DP. We then define the map Λ:L(HA1B1)L(HA2B2) by Λ(X) = Tr(X)ρ. This map has Choi representation J(Λ)=ρ1A1B1, which clearly satisfies Tr(WiJ(Λ)) ≥ 0 for all i. The desired result follows since any two-qubit entangled state ρ is distillable under LOCC channels.□

In this section, we relax the condition that the quantum operation is trace-preserving. Let Λ1 be an arbitrary non-trace-preserving CP map, which by appropriate scaling, we can assume is trace non-increasing, i.e., Tr(X) ≥Tr(Λ1(X)) ≥ 0 for all X ≥ 0. The rescaled map Λ1 can always be “completed” to a channel by another CP map Λ2 so that Λ(⋅) = Λ1(⋅) ⊗|1⟩ ⟨1| + Λ2(⋅) ⊗|2⟩ ⟨2| is the CPTP map. After performing the channel Λ on ρ, the classical register |i⟩ ⟨i| can then be measured so that the post-measurement state Λi(ρ)/Tr[Λi(ρ)] is obtained with probability Tr[Λi(ρ)]. Thus, any non-trace-preserving CP map Λ1 can be stochastically implemented in this way.

While all the discussion thus far has focused on bipartite entanglement, here we consider multipartite systems. To fix notation, let S1SN be N systems with joint state space HS1SN. Convertibility between any two states by stochastic LOCC (SLOCC) provides a natural way to classify entanglement. Under this partitioning of states, two states ρ and ρ̂ are said to be equivalent in terms of their entanglement iff ρρ̂ and ρ̂ρ by SLOCC.17 When considering pure states, it is well-known that two states |ψ⟩ and |ψ̂ belong to the same entanglement class iff there exists invertible linear operators Ai such that |ψ̂=A1A2AN|ψ. A paradigmatic example of two inequivalent states is the three-qubit GHZ state, |GHZ=1/2(|000+|111), and the W state, |W=1/3(|100+|010+|001), which cannot be converted from one to the other using SLOCC.17 In fact, for systems with more than three parties, or for tripartite systems with dimensions greater than 2 ⊗ 3 ⊗ 6, there are an infinite number of inequivalent SLOCC entanglement classes.39 

We will now show that the situation is dramatically different if any family of entangling undetected operations is considered. This is similar in spirit to Theorem 15, except our proof is simpler in that it does not rely on the convex structure of the set of separable states. For a CP map Λ:L(HS1SN)L(HŜ1ŜN), its Choi operator is an N-partite positive semi-definite operator

where ϕS1SN+=ϕS1+ϕSN+. Conversely, for any positive operator ΩL(HŜ1S1::ŜNSN), the map ΛΩ(X)=TrS1SNΩ[1Ŝ1ŜNXT] is CP. Direct calculation shows that J(Ω) is a fully separable operator whenever Λ is a separable map, and conversely ΛΩ is a separable map whenever Ω is fully separable. By fully separable, it is meant that Ω can be expressed as a convex combination of product operators on HŜ1S1::ŜNSN. Hence, we have a one-to-one correspondence between separable operators and separable maps.

Like in the bipartite case, we say that a Hermitian operator W on HŜ1S1::ŜNSN is an entanglement witness if Tr() ≥ 0 for all fully separable N-partite states ρS(HŜ1S1::ŜNSN).

Theorem 16.

Let{Wi}i=1nbe a finite number of entanglement witnesses for operators onHŜ1S1::ŜNSN, and let|ψHS1SNand|ψ̂HŜ1ŜNbe any two states. Then, there exists a CP map Λ that is entangling undetected by{Wi}i=1nand such thatΛ|ψψ|=|ψ̂ψ̂|.

Proof.
Let |ψ*⊥⟩ = |a1⟩⋯|aN⟩ be any product state that is orthogonal to |ψ*⟩. Define the Hermitian operators
(22)
which by construction are entanglement witnesses on HŜ1ŜN. There must exist a separable state ω such that Tr(ωRi) > 0 for all i = 1, …, n. Indeed, such a state can be constructed as follows. For each Ri, let ρi be a separable state such that Tr(Riρi) > 0; since the set of separable states span the space of Hermitian operators, we are assured that such a ρi exists. Then, take ω=i=1nρi, and the fact that each Ri is an entanglement witness guarantees that Tr(Riω) > 0 for all i. Define
(23)
and the positive semi-definite operator
(24)
It can be directly seen that Tr(WiΩ) ≥ 0 for all i, and the CP map given by
satisfies Λ|ψψ|=|ψ̂ψ̂|.□

Theorem 16 has interesting consequences for the problem of tensor rank calculation. The tensor rank of an N-partite state |ψ⟩ is the minimum number of product states whose linear span contains |ψ⟩, i.e.,

(25)

It is easy to see that the tensor rank of a state can be equivalently characterized as an SLOCC convertibility problem,

(26)

where |GHZr(N)=1/ri=1r|iiS1SN.40 Given the one-to-one correspondence between separable maps and separable operators, Eq. (26) can be expressed in terms of entanglement witnesses as

(27)

where Φr(N)=|GHZr(N)GHZr(N)| and WŜ1S1::ŜNSN is the collection of all entanglement witnesses on HŜ1S1:,ŜNSN.

From (25), the tensor rank can be seen as a multipartite generalization of the Schmidt rank. However, unlike the Schmidt rank, the tensor rank is in general very difficult to compute; in fact already in tripartite systems, the problem is NP-complete.41 One way to tackle the tensor rank problem is to start from the characterization (27) and move “beyond” SLOCC. Specifically, instead of considering all entanglement witnesses in the minimization of (27), one could relax the problem and work with some finite subset. A lower bound on this relaxed problem would also be a lower bound on the tensor rank of |ψ⟩. However, Theorem 16 implies that such a strategy will fail since the minimum for any finite set of witnesses will always be one. This finding reflects the general difficulty even in computing non-trivial lower bounds for the tensor rank.

We now examine the classes of k-non-entangling maps and k-PPT-preserving maps. As Proposition 17 indicates, the class of completely non-entangling maps coincides with the class of separable maps (a fact that was first observed in Ref. 42). In particular, a map Λ is completely non-entangling if and only if it is d non-entangling, where d is the dimension of Alice’s and Bob’s input spaces. Since a map is separable if and only if its Choi representation is a separable operator, this parallels the notion that complete positivity of a map is equivalent to positivity of its Choi representation. The same holds for maps that are PPT and completely PPT-preserving. We provide a proof of this fact for the sake of completeness.

Proposition 17

(Ref. 42). LetHA1andHB1bed-dimensional systems and consider a completely positive mapΛ:L(HA1B1)L(HA2B2). The following are equivalent:

  • Λ is separable (respectively, PPT),

  • Λ isd-non-entangling (respectively, d-PPT-preserving), and

  • Λ is completely non-entangling (respectively, completely PPT-preserving).

Proof.
The implications (iii) ⇒ (ii) ⇒ (i) follow from the definitions. It remains to prove the implication (i) ⇒ (iii), which follows from the following observations. Note that the tensor product of separable maps is separable; since if Λ:L(HA1B1)L(HA2B2) and Λ:L(HA1B1)L(HA2B2) are both separable, then ΛΛ:L(HA1A1HB1B1)L(HA2A2HB2B2) can be written as
for some CP maps
such that Λ = ∑iΦi ⊗ Ψi and Λ′ = ∑jΦ′j ⊗ Ψ′j, and thus, Λ ⊗ Λ′ is separable. Analogously, the tensor product of any PPT maps is PPT, which can be clearly seen from the fact that (Λ ⊗ Λ′)Γ = ΛΓ ⊗ Λ′Γ. In particular, the identity map idAB = idA ⊗ idB is separable (respectively, PPT) for any systems A and B; hence, the map Λ ⊗ idAB is separable (respectively, PPT) whenever Λ is separable (respectively, PPT). The desired result follows from the fact that every separable map is non-entangling (and every PPT map is PPT-preserving).□

Note that 1-non-entangling is the same as non-entangling. (Similarly, 1-PPT-preserving is the same as PPT-preserving.) While it is evident that non-entangling is a strictly larger set than completely non-entangling channels, it is not evident from the definitions that there must exist maps that are k-non-entangling for some k < d that are not also completely non-entangling. In the following, we show that there exist k-non-entangling maps which are not (k + 1)-non-entangling for every k < d. This shows that the structure of these maps is as rich as it could be. Nonetheless, it remains an open question as to whether there exist transformations among states that are possible with k-non-entangling but not possible with (k + 1)-non-entangling channels.

We first examine a useful condition for determining when a map on L(CdCd) is k-non-entangling. Note that any pure state |uCdCk can be obtained from |ϕd+ by some rank-k operator X such that

A map Λ on L(CdCd) is k-non-entangling if and only if

is separable for all pure states |u,|vCdCk. To show that Λ is k-non-entangling, it suffices to check only when X is a k-dimensional projection. A k-dimensional projection is a linear operator P:CdCk such that PP = 1k. In particular, for such an operator, it holds that P=Tr(1k)=k, where ∥·∥2 is the Frobenius norm. The above observations allow us to state the following characterization of k-non-entangling maps in Lemma 18.

Lemma 18.
A map Λ onL(CdCd)isk-non-entangling if and only if the operator
(28)
is separable for allk-dimensional projectionsPandQ.

The construction of k-non-entangling maps that are not (k + 1)-non-entangling is based on the well known Werner states. The Werner states are a family of symmetric states on CdCd that are defined by

(29)

for −(d + 1) ≤ βd − 1, where Fd=d(ϕd+)Γ. These states are entangled for β > 0, and furthermore, they are PPT if and only if they are separable. The following lemma presents a fact about Werner states on CdCd that will be employed in our construction of k-non-entangling maps in Theorem 20.

Lemma 19.

Letkanddbe integers with 2 ≤ k < dand let −(d + 1) ≤ βd − 1. The operator (PQ)ρd(β)(PQ)is separable for allk-dimensional projectionsPandQif and only if it holds thatβ ≤ (dk)/k.

Proof.
First assume that β ≤ (dk)/k and let P and Q be k-dimensional projections. Let P¯ denote the matrix whose entries are complex conjugate of those of P, and let |ψ=P¯Q|ϕd+/c denote the normalized pure state with Schmidt rank k and normalization constant c > 0 given by
where the inequality follows from the fact that QPPQ1k with equality if and only if P = Q. Hence, ck/d with equality if and only if P = Q. Since ∥ψΓ2 = 1, note that
(30)
is separable as long as c(β + 1) ≤ 1 (by Theorem 1 of Ref. 30). Then, it is clearly separable since β ≤ (dk)/k by assumption. On the other hand, if β > (dk)/k, we may choose P = Q such that
(31)
We conclude that (PQ)ρd(β)(PQ) is NPT and, therefore, entangled.□

We are now ready to present the construction of k-non-entangling maps in Theorem 20. The maps presented here are trivially non-entangling since they map every state to a separable state. Nonetheless, the maps defined in Theorem 20 are not completely non-entangling and showcase the hierarchy of k-non-entangling maps.

Theorem 20.

For all integerskanddwith 2 ≤ k < d, there exists ak-non-entangling mapΛ:L(CdCd)L(C2C2)that is not (k + 1)-non-entangling.

Proof.
Let k and d be integers with 2 ≤ k < d, and let β > 0 such that
(32)
Consider the map Λβ defined by
(33)
Note that this map is trivially non-entangling since every output is separable. The corresponding Choi operator of this map is given by
(34)
From Lemma 18, we have that Λβ is k-non-entangling but not (k + 1)-non-entangling iff
(35)
is separable for all k-dimensional projections P and Q on Cd but entangled for some (k + 1)-dimensional projections. Note that the operator 1d2ρd(β) is always separable for any β (by Ref. 30), and moreover, the operator in Eq. (35) is separable iff so is (PQ)ρd(β)(PQ). Thus, the use of Lemma 19 together with the assumption in Eq. (32) completes the proof.□

Having seen that k-non-entangling maps for 1 < k < d is a strict superset of separable maps and a strict subset of non-entangling maps, this leads one to consider whether these classes of maps allow for bound entanglement. This question remains unanswered as does the version as to whether all entangled states can be transformed to an LOCC-distillable state under these classes of maps. This is due to the interesting fact that the non-entangling maps constructed used for this task in Sec. IV are not necessarily k-non-entangling for every k < d.

Proposition 21.

The dually non-entangling and PPT maps used in the Proof of Theorem 14 that transform the statesρd(β) with (d − 3)/3 < β ≤ (d − 2)/2 into an LOCC-distillable state are not 3-non-entangling.

Proof.
Lemma 19 tells us that for the values of β mentioned above, there exist projections on a 3-dimensional subspace P and Q such that (PQ)ρd(β)(PQ) is entangled (and NPT). This immediately implies that there exists a pure state |η⟩ with the Schmidt rank equals to 3 such that Tr(ρd(β)W) < 0, where W = |η⟩⟨η|Γ is an entanglement witness. Thus, as argued in Theorem 14, the map Λ of Eq. (21) with A = (W + 21)/3 has the property that Λ(ρd(β)) is an entangled 2-qubit state. We are now going to show that this map is not 3-non-entangling. Since for (d − 3)/3 < β ≤ (d − 2)/2 the states ρd(β) are 1-undistillable,43 the above map cannot be separable and, therefore, its Choi operator
(36)
even though PPT, must be entangled. Now, by Lemma 18, the map Λ can be seen not to be 3-non-entangling if J(Λ)=14(PQ)J(Λ)14(PQ) is still entangled for some choice of projections on a 3-dimensional subspace P and Q. Choosing P and Q to project on the Schmidt bases of |η¯, we obtain that
(37)
which amounts to
(38)
Now, since the operator corresponding to the second term in the right-hand side of the above equation is separable, if J′(Λ) were separable too, this would imply that J(Λ) is separable, which is a contradiction. Thus, J′(Λ) must be entangled.□

Although rather mildly, this suggests that this class of maps could lead to obstructions to distillation allowing for the existence of (NPT) bound entanglement. Thus, we believe that further investigation of the properties of this class deserves future attention.

Although LOCC maps provide the most physically meaningful choice of free operations for the resource theory of entanglement, from a resource-theoretic point of view, other choices of operations still provide consistent and well-defined theories. Furthermore, insight into LOCC is gained by studying more general resource theories since impossibility results for the latter imply impossibility results for the former. We have examined many different classes of non-entangling channels that are still larger than the class of LOCC (and separable) channels. While many entanglement measures are still monotonic under such non-entangling channels, we have presented examples of entanglement measures whose monotonicity fails to hold. This allowed us to find transformations among entangled states that are possible under non-entangling maps but not possible under LOCC (or separable maps). PPT channels have also been widely studied in the context of entanglement theory and, in a similar vein, here we have also shown that, perhaps surprisingly, the negativity is no longer a monotone when this class is extended to the set of PPT-preserving channels.

Although we believe that resource theories of entanglement under different classes of operations is an interesting research topic in itself, one of the motivations behind this present work is to understand more deeply the phenomenon of bound entanglement. The question of whether NPT bound entanglement exists is one of the most challenging open problems in the field. If NPT bound entanglement does exist, this could be proven by finding some NPT entangled state that still remains undistillable by a larger class of operations than LOCC. So far, only PPT operations had been considered in this context.13 Here, we have systematically studied the extraction of entanglement under hierarchies of operational classes that go beyond LOCC. Although this does not necessarily exclude the phenomenon of bound entanglement, we have proven that in most of these classes every entangled state can be converted to a state that is LOCC-distillable. The only class for which we have been unable to prove that all NPT states can be converted to an LOCC-distillable state is that of k-non-entangling maps. Thus, further investigation into the properties of this class could be a promising route to answer the long-standing question of NPT bound entanglement. In particular, even though we have proved that the class of k-non-entangling maps is a strict superset of LOCC and SEP maps, this alone does not necessarily imply that there exist pure-state transformations achievable by the former class of maps, which are impossible by the latter. It would be interesting to study if such examples exist and, if so, which entanglement measures loose their monotonicity under k-non-entangling maps. This would also make it possible to find differences in state manipulation under k-non-entangling maps for different values of k.

The authors would like to thank an anonymous referee who pointed out several similarities and overlaps of our work with the previously known results.

E.C. is supported by the National Science Foundation (NSF) Award No. 1914440. J.I.d.V. acknowledges support from the Spanish MINECO through Grants MTM2017-84098-P and MTM2017-88385-P and from the Comunidad de Madrid through grant QUITEMAD-CM P2018/TCS-4342. M.G. is supported by the Natural Sciences and Engineering Research Council of Canada (NSERC), the Canadian Institute for Advanced Research (CIFAR), and through funding provided to IQC by the Government of Canada. M.G. acknowledges further support from an Izaak Walton Killam Memorial Scholarship and an Alberta Innovates–Technology Futures (AITF) Graduate Student Scholarship. G.G. acknowledges support from NSERC.

This section presents an alternative proof of the result of Theorem 4, which shows the equivalence of the α-entropy of entanglement and the (1/α)-relative entropy of entanglement for α ∈ [0, 2] on pure states. Although this result is already known (see Corollary 2 in Ref. 28), the new proof presented in this appendix is quite different than the one in Ref. 28. Here, we use methods that generalize the methods used for proving the result in the case when α = 1 that were first used in Ref. 23.

We first state the background necessary for proving this result. Recall that the Rényi α-relative entropies Sα are monotonic under CPTP maps for α ∈ [0, 2], and the corresponding α-relative entropies of entanglement are defined by

When ρ = ψ is a pure state, the α-relative entropy reduces (for α ≠ 1) to

While Sα is jointly convex for all α ∈ [0, 1], it is also known to be convex in the second argument for α ∈ [0, 2] (see Ref. 44). This allows us to state a useful minimization criterion for the α-relative entropies of entanglement.45 Indeed, for any state ρ and any separable state σ, it holds that

if and only if it holds that

(A1)

for all other separable states σS. Furthermore, we can limit our consideration only to pure separable states σ′ = |ϕ⟩ ⟨ϕ|, where ϕ is separable.

To show that ER,α(ψ) = E1/α(ψ) holds for all pure states ψ and all α ∈ [0, 2], we must first find a separable state σ such that Sα(ψσ) = E1/α(ψ). We then show that (A1) holds for all other separable pure states σ′. The necessary background for computing the derivatives in (A1) is given in Appendix A.1. Finally, the Proof of Theorem 4 is given in Appendix A.2.

1. Directional derivatives of matrix trace functionals

Any real-valued function f:(0,)R can be extended to the set of positive definite Hermitian matrices by means of the spectral theorem. In particular, for any n × n Hermitian matrix A with spectral decomposition A=i=1nαi|uiui|, one defines

Moreover, for any n × n positive semi-definite matrix P, one may define a function fP on positive definite n × n matrices as

for all positive definite n × n matrices A. Such functions may be extended to certain non-positive-definite matrices in the following way. If A is a n × n positive semi-definite matrix, one defines fP(A) as

(A2)

whenever supp(P) ⊆ supp(A), where the sum is taken only over the indices corresponding to nonzero eigenvalues of A=i=1nλi|uiui|. By extending the function fP to be defined on certain positive semi-definite matrices in this way, it is possible to compute the directional derivatives

A derivation and proof of these directional derivatives is provided in Ref. 46. We provide the details here without proof for clarity.

For any differentiable function f:(0,)R, its (first-order) divided differences are defined as

for all x, y ∈ (0, ). Let f:(0,)R be a differentiable function. For every positive semi-definite n × n matrix A, we will define a linear mapping Φf,A on the space of n × n Hermitian matrices as follows. If A = diag(λ1, …, λn) is diagonal, we define an matrix n × n Hermitian matrix Df,A of the corresponding matrix of divided differences (restricted to the strictly positive eigenvalues) as

(A3)

and define the linear mapping Φf,A as

for all n × n Hermitian matrices B, where AB denotes the entry wise product of the two matrices with matrix entries (AB)i,j = Ai,jBi,j for any n × n Hermitian matrices A and B. Thus, Φf,A(B) has matrix elements

If A is not diagonal, there is an n × n diagonalizing unitary matrix U such that the matrix UAU is diagonal, and one defines the linear mapping Φf,A as

(A4)

for all n × n Hermitian matrices B, and this mapping is independent of the choice of diagonalizing unitary U.

The main result of Ref. 46 (Theorem 1) states that directional derivatives of functions of the form A ↦ Tr(Pf(A)) can be computed at certain non-positive-definite matrices as follows. Suppose that f:(0,)R is differentiable and satisfies

(A5)

Let P be a positive semi-definite n × n matrix, and let A be another positive semi-definite n × n matrix satisfying supp(P) ⊆ supp(A) such that fP(A) may be defined as in (A2). Let B be any n × n Hermitian matrix and suppose there exists a positive value ɛ > 0 such that A + tB is a positive semi-definite for all t ∈ [0, ɛ). Then, one may compute the directional derivative of fP at A in the direction B as

(A6)

In the case when A is a positive definite, the mapping Φf,A is precisely the Fréchet derivative of the function f at A (see, e.g., Ref. 47, Sec. X.4 and Ref. 48, Theorem 3.25)

Note that the function f(x) = log(x) satisfies the condition in (A5), as well as the functions fα(x) = x1−α for all α < 2. However, although the condition is not satisfied by the function f(x) = x−1, in Ref. 46, Sec. V, it is shown that the expression in (A6) is still a lower bound for the directional derivative. That is, if the matrices P, A, and B satisfy all the criteria as given above, then

(A7)

and this inequality is strict in general if A is not a positive definite.

Before providing the Proof of Theorem 4, we make a few more remarks regarding directional derivatives the linear mappings Φf,A given above as to how they will be applied here. Let f:(0,)R be a differentiable function satisfying (A5), and let ρ,σ,σD(H) be states on a finite-dimensional Hilbert space H satisfying supp(ρ) ⊆ supp(σ). It holds that

is a positive semi-definite for all t ∈ [0, 1]; hence, the directional derivative of fρ at σ is computed as

Moreover, as the mapping Φf,σ [as defined in (A4)] is a linear, it can be seen that

where f′(σ) is defined as in (A2).

2. Proof of Theorem

We now present the proof that ER,α(ψ) = E1/α(ψ) for all pure states ψ and all α ∈ [0, 2]. This result is already known when α = 1,23 so here we prove it only for α ≠ 1. We will consider the case when α = 0 separately.

Here, we will make use of the functions fα for α ∈ [0, 1) ∪ (1, 2] that are defined by fα(x) = x1−α and whose divided differences are given by

(A8)

for all x, y ∈ (0, ). For these functions, it is straightforward to check that

holds for any positive real constant c > 0 and any x, y > 0. We now prove Theorem 4.

As remarked in Appendix A.1, the function fα satisfies condition (A5) for all α < 2. This allows us to compute the directional derivatives of the α-relative entropies Sα(ρσ). For α = 2, we may compute a lower bound to the derivative as in (A7).

Proof of Theorem 4.
Let ψ=iλi|ii be a pure state with Schmidt coefficients λ = (λ1, …, λd). First suppose that α ≠ 0 and define the following separable density operator
(A9)
where λ1/α=iλi1/αα. Note that the α-relative entropy of ψ and σ is
from which it follows that Sα(ψσ) = E1/α(ψ). The matrix σ has eigenvalues λi1/α/λ1/α1/α, and the corresponding matrix Dfα,σ of divided differences as defined in (A3) has non-zero matrix elements given by
(A10)
such that the linear mapping Φfα,σ is defined as
for all B. We first suppose that α ∈ (0, 1) ∪ (1, 2) (the case α = 2 must be considered separately). For any other separable state σ′, we may use (A6) to compute the necessary derivative
where the final equality follows from the fact that ⟨ψ|σ1−α|ψ⟩ = ∥λ1/α, the linearity of the mapping Φfα,σ, and the fact that Φfα,σ(σ)=(1α)σ1α. Thus,
When σ′ = |ϕ⟩ ⟨ϕ| is a separable pure state for some |ϕ⟩ = ∑i,juivj|ij⟩ with ∑i|ui|2 = ∑j|vj|2 = 1, it holds that
Note that fα[1](x,y)/(1α)0 for any x, y > 0. Hence,
where the second inequality is due to Lemma 22 (see below) and the final inequality follows from the Cauchy–Schwarz inequality. This completes the proof since we have found that the inequality
(A11)
holds for all pure separable states σ′.
To prove the claim in the case when α = 2, we note that only a lower bound the directional derivative can be provided for the function f2(x) = x−1. In this case, we have
and using the same arguments as above one finds
which is non-negative (by Lemma 22 and the same arguments as above).
Finally, to prove the claim when α = 0, set σ = |11⟩ ⟨11|. Note that S0(ψσ) = −log λ1, where λ1 is the largest Schmidt coefficient of ψ. For any other separable state σS, it holds that ⟨ψ|σ′|ψ⟩ ≤ λ1 (by Lemma 6) and thus,
from which it follows that ER,0(ψ) = −log λ1. Since E+(ψ) = −log λ1, the result follows.□

The Proof of Theorem 4 above in the case when α ≠ 0 relies on the following lemma, which we prove below.

Lemma 22.
Letα ∈ (0, 1) ∪ (1, 2]. For allp, q ∈ (0, 1], it holds that
(A12)
wherefα:(0,)Ris the function defined byfα(x) = x1−α.

The proof is trivial in the case when α = 2. Indeed, f2(x) = x−1 and thus, for pq, it holds that

The proof in the case when α ≠ 2 is a bit more technical and requires the following pair of results.

Lemma 23.
For allx, y > 0 withxyandr ∈ (−1, 0) ∪ (0, 1), it holds that
(A13)

Proof.
For x, y > 0 and r ∈ (−1, 0) ∪ (0, 1), we have the following integral representations:
(A14)
(A15)
Let x, y > 0 with xy. Note that x+y2xy and thus,
holds for all t ≥ 0. Let r ∈ (−1, 0) ∪ (0, 1) and note that r sin() ≥ 0. Thus,
where we use the representations in (A14) and (A15).□

Lemma 24.
For allp, q ∈ (0, 1] withpq, and allα ∈ (0, 1) ∪ (1, 2), it holds that
(A16)

Proof.
Setting r = 1 − α, x = p1/α, and y = q1/α, an application of of Lemma 23 yields
as desired.□

A Proof of Lemma 22 in the case α ≠ 2 now follows.

Proof of Lemma 22.
As the statement has been shown to be true for α = 2, we may assume that α ≠ 2. If pq, we can apply Lemma 24. Otherwise, we have
and thus, p1αfα[1](p1/α,p1/α)=1, which completes the proof.□

Here, we show that the values δ = d−4 and ϵ = d−12 satisfy both (10) and (11). Indeed, for these values of ϵ and δ, it is straightforward to check that

and thus,

which implies that (10) holds. On the other hand, it is clear that

Furthermore, it holds that 1k1 since k ≥ 2. It follows that

which proves that (11) holds as desired.

1.
R.
Horodecki
,
P.
Horodecki
,
M.
Horodecki
, and
K.
Horodecki
, “
Quantum entanglement
,”
Rev. Mod. Phys.
81
(
2
),
865
942
(
2009
); arXiv:quant-ph/0702225.
2.
M. B.
Plenio
and
S.
Virmani
, “
An introduction to entanglement measures
,”
Quantum Inf. Comput.
7
(
1
),
1
51
(
2007
); arXiv:quant-ph/0504163.
3.
E.
Chitambar
and
G.
Gour
, “
Quantum resource theories
,”
Rev. Mod. Phys.
91
,
025001
(
2019
); arXiv:1806.06107.
4.
E.
Chitambar
,
D.
Leung
,
L.
Mančinska
,
M.
Ozols
, and
A.
Winter
, “
Everything you always wanted to know about LOCC (but were afraid to ask)
,”
Commun. Math. Phys.
328
(
1
),
303
326
(
2014
); arXiv:1210.4583.
5.
G.
Vidal
and
J. I.
Cirac
, “
Irreversibility in asymptotic manipulations of entanglement
,”
Phys. Rev. Lett.
86
(
25
),
5803
5806
(
2001
); arXiv:quant-ph/0102036.
6.
D.
Yang
,
M.
Horodecki
,
R.
Horodecki
, and
B.
Synak-Radtke
, “
Irreversibility for all bound entangled states
,”
Phys. Rev. Lett.
95
(
19
),
190501
(
2005
); arXiv:quant-ph/0506138.
7.
M.
Horodecki
,
P.
Horodecki
, and
R.
Horodecki
, “
Mixed-state entanglement and distillation: Is there a “bound” entanglement in nature?
,”
Phys. Rev. Lett.
80
(
24
),
5239
5242
(
1998
); arXiv:quant-ph/9801069.
8.
F. G. S. L.
Brandão
and
M. B.
Plenio
, “
Entanglement theory and the second law of thermodynamics
,”
Nat. Phys.
4
(
11
),
873
877
(
2008
); arXiv:0810.2319.
9.
F. G. S. L.
Brandão
and
M. B.
Plenio
, “
A reversible theory of entanglement and its relation to the second law
,”
Commun. Math. Phys.
295
(
3
),
829
851
(
2010
); arXiv:0710.5827.
10.
F. G. S. L.
Brandão
,
M.
Christandl
, and
J.
Yard
, “
Faithful squashed entanglement
,”
Commun. Math. Phys.
306
(
3
),
805
830
(
2011
); arXiv:1010.1750.
11.
E. M.
Rains
, “
Entanglement purification via separable superoperators
,” arXiv:quant-ph/9707002 (
1997
).
12.
Ł.
Pankowski
,
M.
Piani
,
M.
Horodecki
, and
P.
Horodecki
, “
A few steps more towards NPT bound entanglement
,”
IEEE Trans. Inf. Theory
56
(
8
),
4085
4100
(
2010
); arXiv:0711.2613.
13.
T.
Eggeling
,
K. G. H.
Vollbrecht
,
R. F.
Werner
, and
M.
Wolf
, “
Distillability via protocols respecting the positivity of partial transpose
,”
Phys. Rev. Lett.
87
(
25
),
257902
(
2001
); arXiv:quant-ph/0104095.
14.
E. M.
Rains
, “
Rigorous treatment of distillable entanglement
,”
Phys. Rev. A
60
(
1
),
173
178
(
1999
); arXiv:quant-ph/9809078.
15.
E. M.
Rains
, “
Erratum: Bound on distillable entanglement [Phys. Rev. A 60, 179 (1999)]
,”
Phys. Rev. A
63
(
1
),
019902
(
2000
).
16.
M. A.
Nielsen
, “
Conditions for a class of entanglement transformations
,”
Phys. Rev. Lett.
83
(
2
),
436
439
(
1999
).
17.
W.
Dür
,
G.
Vidal
, and
J. I.
Cirac
, “
Three qubits can be entangled in two inequivalent ways
,”
Phys. Rev. A
62
(
6
),
062314
(
2000
); arXiv:quant-ph/0005115.
18.
M.
Horodecki
,
P.
Horodecki
, and
R.
Horodecki
, “
Separability of mixed states: Necessary and sufficient conditions
,”
Phys. Lett. A
223
(
1-2
),
1
8
(
1996
); arXiv:quant-ph/9605038.
19.
L. M.
Ioannou
and
B. C.
Travaglione
, “
Quantum separability and entanglement detection via entanglement-witness search and global optimization
,”
Phys. Rev. A
73
(
5
),
052314
(
2006
); arXiv:quant-ph/0602223.
20.
O.
Gühne
and
N.
Lütkenhaus
, “
Nonlinear entanglement witnesses, covariance matrices and the geometry of separable states
,”
J. Phys. Conf. Ser.
67
,
012004
(
2007
); arXiv:quant-ph/0612108.
21.
E. M.
Bronstein
, “
Approximation of convex sets by polytopes
,”
J. Math. Sci.
153
(
6
),
727
762
(
2008
).
22.
V.
Vedral
,
M. B.
Plenio
,
M. A.
Rippin
, and
P. L.
Knight
, “
Quantifying entanglement
,”
Phys. Rev. Lett.
78
(
12
),
2275
2279
(
1997
); arXiv:quant-ph/9702027.
23.
V.
Vedral
and
M. B.
Plenio
, “
Entanglement measures and purification procedures
,”
Phys. Rev. A
57
(
3
),
1619
1633
(
1998
).
24.
G.
Lindblad
, “
Completely positive maps and entropy inequalities
,”
Commun. Math. Phys.
40
(
2
),
147
151
(
1975
).
25.
G.
Vidal
and
R.
Tarrach
, “
Robustness of entanglement
,”
Phys. Rev. A
59
(
1
),
141
155
(
1999
); arXiv:quant-ph/9806094.
26.
B. M.
Terhal
and
P.
Horodecki
, “
A schmidt number for density matrices
,”
Phys. Rev. A
61
,
040301
(
2000
); arXiv:quant-ph/9911117.
27.
D.
Petz
, “
Quasi-entropies for finite quantum systems
,”
Rep. Math. Phys.
23
(
1
),
57
65
(
1986
).
28.
H.
Zhu
,
M.
Hayashi
, and
L.
Chen
, “
Coherence and entanglement measures based on Rényi relative entropies
,”
J. Phys. A: Math. Theor.
50
(
47
),
475303
(
2017
).
29.
A.
Shimony
, “
Degree of entanglement
,”
Ann. N. Y. Acad. Sci.
755
,
675
679
(
1995
).
30.
L.
Gurvits
and
H.
Barnum
, “
Largest separable balls around the maximally mixed bipartite quantum state
,”
Phys. Rev. A
66
(
6
),
062311
(
2002
); arXiv:quant-ph/0204159.
31.
K.
Audenaert
,
M. B.
Plenio
, and
J.
Eisert
, “
Entanglement cost under positive-partial-transpose-preserving operations
,”
Phys. Rev. Lett.
90
(
2
),
027901
(
2003
); arXiv:quant-ph/0207146v1.
32.
M.
Horodecki
,
P.
Horodecki
, and
R.
Horodecki
, “
Inseparable two spin- 12 density matrices can be distilled to a singlet form
,”
Phys. Rev. Lett.
78
(
4
),
574
577
(
1997
).
33.
S.
Bravyi
and
A.
Kitaev
, “
Universal quantum computation with ideal clifford gates and noisy ancillas
,”
Phys. Rev. A
71
,
022316
(
2005
).
34.
C. H.
Bennett
,
D. P.
DiVincenzo
,
J. A.
Smolin
, and
W. K.
Wootters
, “
Mixed-state entanglement and quantum error correction
,”
Phys. Rev. A
54
(
5
),
3824
3851
(
1996
); arXiv:quant-ph/9604024v2.
35.
R.
Horodecki
and
M.
Horodecki
, “
Information-theoretic aspects of inseparability of mixed states
,”
Phys. Rev. A
54
(
3
),
1838
1843
(
1996
); arXiv:quant-ph/9607007v1.
36.
J. A.
Smolin
, “
Four-party unlockable bound entangled state
,”
Phys. Rev. A
63
(
3
),
032306
(
2001
); arXiv:quant-ph/0001001v6.
37.
G.
Vidal
, “
Entanglement of pure states for a single copy
,”
Phys. Rev. Lett.
83
(
5
),
1046
1049
(
1999
); arXiv:quant-ph/9902033.
38.
A.
Barvinok
,
A Course in Convexity
(
American Mathematical Society
,
2002
).
39.
E.
Chitambar
,
C. A.
Miller
, and
Y.
Shi
, “
Matrix pencils and entanglement classification
,”
J. Math. Phys.
51
(
7
),
072205
(
2010
); arXiv:0911.4058v1.
40.
E.
Chitambar
,
R.
Duan
, and
Y.
Shi
, “
Tripartite entanglement transformations and tensor rank
,”
Phys. Rev. Lett.
101
(
14
),
140502
(
2008
); arXiv:0805.2977.
41.
J.
Håstad
, “
Tensor rank is NP-complete
,”
J. Algorithms
11
(
4
),
644
654
(
1990
).
42.
J. I.
Cirac
,
W.
Dür
,
B.
Kraus
, and
M.
Lewenstein
, “
Entangling operations and their implementation using a small amount of entanglement
,”
Phys. Rev. Lett.
86
(
3
),
544
547
(
2001
); arXiv:quant-ph/0007057.
43.
W.
Dür
,
J.
Cirac
,
M.
Lewenstein
, and
D.
Bruß
, “
Distillability and partial transposition in bipartite systems
,”
Phys. Rev. A
61
(
6
),
062313
(
2000
); arXiv:quant-ph/9910022.
44.
M.
Mosonyi
and
F.
Hiai
, “
On the quantum Rényi relative entropies and related capacity formulas
,”
IEEE Trans. Inf. Theory
57
(
4
),
2474
2487
(
2011
); arXiv:0912.1286.
45.
M. W.
Girard
,
G.
Gour
, and
S.
Friedland
, “
On convex optimization problems in quantum information theory
,”
J. Phys. A: Math. Theor.
47
(
50
),
505302
(
2014
); arXiv:1402.0034.
46.
M. W.
Girard
, “
On directional derivatives of trace functionals of the form Tr(Pf(A))
,”
Linear Algebra Appl.
569
,
62
77
(
2019
); arXiv:1810.05641.
47.
R.
Bhatia
.
Matrix Analysis
, Graduate Texts in Mathematics Vol. 169 (
Springer
,
1997
).
48.
F.
Hiai
and
D.
Petz
,
Introduction to Matrix Analysis and Applications
(
Springer
,
2014
).
49.

Note that a comparison of the distillable entanglement to the entanglement cost under PPT channels has been carried out in Ref. 31, whose authors use the terminology “operations preserving the positivity of partial transpose” for what we call PPT maps.