alphaFeedback

WorksPhysics and Astronomy

Correlated Catalysis Removes the Strong Converse for Entanglement Distillation

AI contributions
  • ChatGPT 6 Astra
View details

Abstract

Entanglement distillation converts noisy shared quantum states into maximally entangled pairs using local operations and classical communication. The distillable entanglement EDE_D determines the maximum asymptotic rate at which this process can be performed reliably. For several classes of quantum states, distillation above this rate obeys a strong converse: the fidelity with the target state vanishes asymptotically. Here we show that correlated catalysis removes the strong converse for every distillable bipartite state. Specifically, at any target rate R>EDR>E_D, the optimal asymptotic fidelity is given by ED/RE_D/R. Thus, although catalysis does not increase the maximum reliable distillation rate, it allows nonvanishing fidelity at every finite rate. Remarkably, this advantage persists even when the catalyst is returned exactly at every finite blocklength and the correlations between the output and catalyst vanish asymptotically. We further characterize the catalyst size required to sustain nonvanishing fidelity above the reliable distillation rate.

Introduction

Quantum entanglement [1] is a fundamental type of correlation between quantum particles, enabling communication protocols that surpass classical limits. As a uniquely quantum phenomenon, entanglement acts as a vital resource for tasks such as quantum teleportation [2], superdense coding [3], and quantum key distribution [4, 5]. These applications motivate the development of methods for preparing high-quality entangled states.

Entanglement is often studied in the context of resource theory. In this setting, two parties, say Alice and Bob, can perform only local quantum channels, and communicate through a classical channel. These operations are called local operations and classical communication (LOCC) [6].

Entanglement distillation enables the extraction of nearly maximally entangled states from noisy entangled input using local operations and classical communication [7, 8, 9, 10, 11, 12, 13]. For a finite number of input copies, both deterministic protocols based on quantum error correction [8, 13] and probabilistic protocols based on postselection [7, 9] have been studied. In the asymptotic regime, the central objective is to determine the optimal rate at which maximally entangled pairs can be extracted as the number of input copies tends to infinity [8, 14, 15].

In the asymptotic regime, a distillation protocol Dn\mathcal{D}_n transforms nn copies of a bipartite state ρ\rho into an output state approximating mnm_n maximally entangled qubit pairs. The distillation rate is the number of target pairs per input copy, R=lim⁡n→∞mn/nR = \lim_{n\to\infty} m_n/n. The distillable entanglement ED(ρ)E_D(\rho) is the supremum of rates achievable such that the fidelity between the output state and [Rn][Rn] maximally entangled states tends to unity [16]:

ED(ρ)=sup⁡{R≥0:lim sup⁡n→∞Tr⁡[Φ2⊗[Rn]Dn(ρ⊗n)]=1}.(1)E_D(\rho)=\sup\left\{R\geq0:\limsup_{n\to\infty}\operatorname{Tr}\left[\Phi_2^{\otimes[Rn]}\mathcal{D}_n\left(\rho^{\otimes n}\right)\right]=1\right\}. \tag*{(1)}

where the supremum is taken over rates RR for which there exists a sequence of LOCC channels {Dn}n∈N\{\mathcal{D}_n\}_{n\in\mathbb{N}}, and Φ2\Phi_2 denotes the density operator of a maximally entangled qubit pair. Operationally, ED(ρ)E_D(\rho) quantifies the asymptotic number of Bell pairs that can be extracted per input copy with vanishing error [17].

The definition of distillable entanglement rules out asymptotically perfect distillation at rates R>ED(ρ)R > E_D(\rho), but does not determine how the fidelity behaves in this regime. A stronger requirement is the strong converse property: for every rate R>ED(ρ)R > E_D(\rho), the fidelity of the output with the target state Φ2⊗[Rn]\Phi_2^{\otimes[Rn]} must tend to zero as n→∞n\to\infty, irrespective of the sequence of LOCC protocols employed [18]. Thus, when the strong converse holds, allowing a fixed nonzero error does not increase the asymptotically achievable distillation rate.

It is known that the strong converse holds for bipartite pure states [19] and for maximally correlated states [20]. More generally, the Rains relative entropy provides a strong-converse upper bound on the distillation rate: whenever the rate exceeds this bound, the fidelity vanishes exponentially with the number of input copies [21, 22]. Whether the strong converse holds for unassisted distillation from arbitrary bipartite states remains open. In this paper, we investigate whether the strong converse also holds for entanglement distillation protocols assisted by a catalyst.

Entanglement catalysis provides another way to extend the possibilities of quantum state manipulation under LOCC [23, 24, 25, 26, 27]. In its original formulation [23], Alice and Bob have access to an auxiliary entangled state ωC\omega_C, called a catalyst, which assists a transformation while being returned unchanged and uncorrelated with the output. A catalytic transformation therefore takes the form

ΛLOCC(ρ⊗ωC)=σ⊗ωC.(2)\Lambda_{\mathrm{LOCC}}(\rho\otimes\omega_C)=\sigma\otimes\omega_C. \tag*{(2)}

Remarkably, such a transformation may be possible even when ρ\rho cannot be converted into σ\sigma by LOCC alone. This framework was subsequently extended to correlated catalysis [28, 29, 30], in which the catalyst is returned exactly in its marginal state but may become correlated with the output system [31, 32, 33]. The final state ΩSC\Omega_{SC} is such that

ΩSC=ΛLOCC(ρ⊗ωC),ΩS=σ,ΩC=ωC,(3)\Omega_{SC}=\Lambda_{\mathrm{LOCC}}(\rho\otimes\omega_C),\qquad\Omega_S=\sigma,\qquad\Omega_C=\omega_C, \tag*{(3)}

instead of ΩSC=σ⊗ωC\Omega_{SC}=\sigma\otimes\omega_C. Allowing such residual correlations can enlarge the set of achievable state transformations [31].

Entangled states with positive partial transpose remain undistillable under correlated catalytic local operations and classical communication, while retaining a positive entanglement cost [33]. Thus, correlated catalysis does not make entanglement theory asymptotically reversible.

In this work, we show that correlated catalysis removes the strong converse for every distillable bipartite state, even when the correlations between the output and the catalyst vanish asymptotically. More precisely, at any rate R>ED(ρ)R>E_D(\rho), the optimal asymptotic fidelity equals ED(ρ)/RE_D(\rho)/R, whereas fidelity approaching unity is achievable for R≤ED(ρ)R\le E_D(\rho). These limits remain achievable under any prescribed positive sequence of upper bounds on the output–catalyst mutual information tending to zero.

We also determine the catalytic resources needed to sustain non-vanishing fidelity above the reliable distillation rate. At rates where non-catalytic distillation obeys an exponential strong converse, the catalyst’s max-relative entropy of entanglement must grow at least linearly with the input blocklength nn. Requiring the output–catalyst mutual information to satisfy I(S:C)≤δnI(S:C)\le\delta_n, for any prescribed positive sequence δn→0\delta_n\to0, strengthens this lower bound to Ω(nlog⁡2(1/δn))\Omega(n\log_2(1/\delta_n)), and hence imposes the same lower bound on the logarithm of each local catalyst dimension. Conversely, for every distillable source, every R>ED(ρ)R>E_D(\rho), and every fixed 0<f0<ED(ρ)/R0<f_0<E_D(\rho)/R, we construct correlated catalytic protocols achieving fidelity at least f0f_0 with log⁡2dim⁡CA,log⁡2dim⁡CB=O(nlog⁡2(1/δn))\log_2\dim C_A,\log_2\dim C_B=O(n\log_2(1/\delta_n)). Thus, whenever the exponential-converse assumption holds, the optimal catalyst dimension satisfies log⁡2Dnopt=Θ(nlog⁡2(1/δn))\log_2 D_n^{\mathrm{opt}}=\Theta(n\log_2(1/\delta_n)). In particular, this conclusion applies to pure and maximally correlated source states at every rate above their distillable entanglement.

Definitions and problem statement

Let

∣Φd⟩=1d∑x=0d−1∣x⟩A∣x⟩B,Φd=∣Φd⟩⟨Φd∣(4)\lvert\Phi_d\rangle=\frac{1}{\sqrt{d}}\sum_{x=0}^{d-1}\lvert x\rangle_A\lvert x\rangle_B,\qquad\Phi_d=\lvert\Phi_d\rangle\langle\Phi_d\rvert \tag*{(4)}

be a maximally entangled state with local dimensions dd. The mm-fold tensor product of a maximally entangled qubit pair is therefore Φ2⊗m=Φ2m\Phi_2^{\otimes m}=\Phi_{2^m}.

As discussed in the introduction, catalysis can enable transformations that are otherwise impossible under LOCC alone. Here, we study entanglement distillation assisted by a bipartite catalyst. We define a distillation protocol with catalyst as follows [31, 32]:

Definition 1. Correlated catalytic distillation protocol

Fix a bipartite source state ρ\rho and a target of mm Bell pairs. A correlated catalytic distillation protocol is a deterministic LOCC channel Ddist\mathcal{D}_{\mathrm{dist}} that takes nn copies of the state ρ\rho and a bipartite catalyst ωC\omega_C and outputs

ΩSC=Ddist(ρ⊗n⊗ωC),ΩC=ωC.(5)\Omega_{SC}=\mathcal{D}_{\mathrm{dist}}(\rho^{\otimes n}\otimes\omega_C),\qquad\Omega_C=\omega_C. \tag*{(5)}

We quantify the performance of a catalytic distillation protocol by the overlap of its output marginal with the target maximally entangled state:

f(ΩSC,m):=Tr⁡[(Φ2m⊗IC)ΩSC]=Tr⁡(Φ2mΩS).(6)f(\Omega_{SC},m):=\operatorname{Tr}\left[(\Phi_{2^m}\otimes I_C)\Omega_{SC}\right]=\operatorname{Tr}(\Phi_{2^m}\Omega_S). \tag*{(6)}

The correlations between the output and the catalyst can be quantified by the quantum mutual information [34]

I(S:C)Ω=S(ΩS)+S(ωC)−S(ΩSC),(7)I(S:C)_\Omega=S(\Omega_S)+S(\omega_C)-S(\Omega_{SC}), \tag*{(7)}

where S(ρ)=−Tr⁡(ρlog⁡(ρ))S(\rho)=-\operatorname{Tr}(\rho\log(\rho)) denotes the von Neumann entropy [27]. Additionally, we can bound the trace distance to uncorrelated catalysis by the quantum Pinsker inequality with quantum mutual information [35],

12∥ΩSC−ΩS⊗ωC∥1≤ln⁡22I(S:C)Ω.(8)\frac{1}{2}\lVert\Omega_{SC}-\Omega_S\otimes\omega_C\rVert_1\le\sqrt{\frac{\ln2}{2}I(S:C)_\Omega}. \tag*{(8)}

In particular, I(S:C)Ω=0I(S:C)_\Omega=0 implies ΩSC=ΩS⊗ωC\Omega_{SC}=\Omega_S\otimes\omega_C. Moreover, vanishing mutual information implies that the joint state approaches the product of its marginals in trace distance.

Throughout this work, we consider finite-dimensional catalysts with unrestricted dimension and entanglement, and investigate protocols in which the residual correlations between the system ΩS\Omega_S and catalyst ΩC\Omega_C can be made arbitrarily small while preserving the exact catalyst marginal.

Suppose that we have nn copies of a state ρ\rho and want to produce mm copies of Φ2\Phi_2 using correlated catalytic entanglement distillation. The maximum achievable fidelity with a target state Φ2m\Phi_{2^m} is:

Fρ→Φ2cc(n,m;δ)=sup⁡ωC,DdistTr⁡(Φ2m⊗mΩS),(9)F_{\rho\to\Phi_2}^{\mathrm{cc}}(n,m;\delta)=\sup_{\omega_C,\mathcal{D}_{\mathrm{dist}}}\operatorname{Tr}\left(\Phi_{2^m}^{\otimes m}\Omega_S\right), \tag*{(9)}

where the supremum is taken over all finite-dimensional catalysts and deterministic LOCC channels satisfying ΩC=ωC\Omega_C=\omega_C and I(S:C)Ω≤δI(S:C)_\Omega\le\delta. See Eq. (5) for the definition of ΩSC\Omega_{SC}. We use δ=∞\delta=\infty when no restriction on the output–catalyst correlations is imposed.

For a fixed distillation rate R>0R>0, let mn=⌊Rn⌋m_n=\lfloor Rn\rfloor. The reliable catalytic distillation rate is defined as

Dcc(ρ)=sup⁡{R:lim sup⁡n→∞Fρ→Φ2cc(n,mn;∞)=1}.(10)D^{\mathrm{cc}}(\rho)=\sup\left\{R:\limsup_{n\to\infty}F_{\rho\to\Phi_2}^{\mathrm{cc}}(n,m_n;\infty)=1\right\}. \tag*{(10)}

This quantity coincides with the correlated-catalytic distillable entanglement Ed,LOCCcc(ρ)E_{d,\mathrm{LOCCcc}}(\rho) defined in [33], since vanishing trace-distance error and fidelity tending to one are equivalent for the pure target state. The corresponding strong-converse threshold is

D~cc(ρ)=sup⁡{R:lim sup⁡n→∞Fρ→Φ2cc(n,mn;∞)>0}.(11)\widetilde{D}^{\mathrm{cc}}(\rho)=\sup\left\{R:\limsup_{n\to\infty}F_{\rho\to\Phi_2}^{\mathrm{cc}}(n,m_n;\infty)>0\right\}. \tag*{(11)}

The strong converse holds if every rate exceeding the reliable distillation rate necessarily leads to vanishing global fidelity in the limit n→∞n\to\infty. In other words, for an input state ρ\rho, the strong converse holds iff Dcc(ρ)=D~cc(ρ)D^{\mathrm{cc}}(\rho)=\widetilde{D}^{\mathrm{cc}}(\rho). In this paper we show that strong converse does not hold for catalytic entanglement distillation protocols.

Throughout this work, all logarithms are base two unless stated otherwise. We use standard asymptotic notation: for nonnegative sequences ana_n and bnb_n, we write an=O(bn)a_n=O(b_n) if there exists a constant C>0C>0 such that an≤Cbna_n\le Cb_n for all sufficiently large nn. We write an=Ω(bn)a_n=\Omega(b_n) if there exists a constant c>0c>0 such that an≥cbna_n\ge cb_n for all sufficiently large nn. Finally, we write an=Θ(bn)a_n=\Theta(b_n) if there exist constants c,C>0c,C>0 such that cbn≤an≤Cbncb_n\le a_n\le Cb_n for all sufficiently large nn.

Throughout this work, fidelity refers to the squared fidelity; for a pure target Φ\Phi, it is given by Tr⁡(Φρ)\operatorname{Tr}(\Phi\rho).

Main results

We first show that exact catalysis with vanishing correlations removes the strong converse for every distillable source state. At any rate R>ED(ρ)R>E_D(\rho), the optimal asymptotic fidelity is ED(ρ)/RE_D(\rho)/R, even when output–catalyst correlations vanish asymptotically.

Theorem 2. Let ρ\rho be a finite-dimensional bipartite state with ED(ρ)>0E_D(\rho)>0. For every prescribed positive sequence δn→0\delta_n\to0 and every R>0R>0,

lim⁡n→∞Fρ→Φ2cc(n,⌊Rn⌋;δn)=min⁡{1,ED(ρ)R}.(12)\lim_{n\to\infty}F_{\rho\to\Phi_2}^{\mathrm{cc}}(n,\lfloor Rn\rfloor;\delta_n)=\min\left\{1,\frac{E_D(\rho)}{R}\right\}. \tag*{(12)}

Consequently,

D~cc(ρ)=∞.(13)\widetilde{D}^{\mathrm{cc}}(\rho)=\infty. \tag*{(13)}

The proof is given in § 1.3. In particular, if ρ=Φ2\rho=\Phi_2 we obtain

lim⁡n→∞FΦ2→Φ2cc(n,⌊Rn⌋;δn)=min⁡{1,1R}.(14)\lim_{n\to\infty}F_{\Phi_2\to\Phi_2}^{\mathrm{cc}}(n,\lfloor Rn\rfloor;\delta_n)=\min\left\{1,\frac{1}{R}\right\}. \tag*{(14)}
D~cc(Φ2)=∞.(15)\widetilde{D}^{\mathrm{cc}}(\Phi_2)=\infty. \tag*{(15)}

Note that the result Dcc(ρ)=ED(ρ)D_{\mathrm{cc}}(\rho)=E_D(\rho) was previously established for distillable states in Ref. [32].

The construction used in the proof of Theorem 2 achieves the optimal asymptotic fidelity ED(ρ)/RE_D(\rho)/R at rates R>ED(ρ)R>E_D(\rho), but requires a large catalyst memory.

In the following, we derive lower bounds on the catalyst resources needed to maintain a fixed positive fidelity at rates where non-catalytic distillation obeys an exponential decay of fidelity.

Let

Fρ→Φ2not(n,m)=sup⁡Λ∈LOCCTr⁡(Φ2mΛ(ρ⊗n))(16)F_{\rho\to\Phi_2}^{\mathrm{not}}(n,m)=\sup_{\Lambda\in\mathrm{LOCC}}\operatorname{Tr}\left(\Phi_2^m\Lambda(\rho^{\otimes n})\right) \tag*{(16)}

denote the optimal fidelity without a catalyst. For a bipartite catalyst state ωC\omega_C, define its max-relative entropy of entanglement [36]

Emax⁡(ωC)=inf⁡τ∈SEPDmax⁡(ωC∥τ),(17)E_{\max}(\omega_C)=\inf_{\tau\in\mathrm{SEP}}D_{\max}(\omega_C\mathbin\Vert\tau), \tag*{(17)}
Dmax⁡(x∥y)=inf⁡{λ:x≤2λy}.(18)D_{\max}(x\mathbin\Vert y)=\inf\{\lambda:x\le2^\lambda y\}. \tag*{(18)}

The Rains relative entropy Rent(ρ)R_{\mathrm{ent}}(\rho), defined as [37]

Rent(ρ):=inf⁡σ≥0∥σTB∥1≤1D(ρ∥σ),(19)R_{\mathrm{ent}}(\rho):=\inf_{\sigma\ge0\atop\|\sigma^{T_B}\|_1\le1}D(\rho\mathbin\Vert\sigma), \tag*{(19)}

provides an exponential converse-bound for non-catalytic distillation of any finite-dimensional bipartite state [21, 22], i.e. every rate RR satisfying R>Rent(ρ)R>R_{\mathrm{ent}}(\rho) leads to exponential decay of Fρ→Φ2notF_{\rho\to\Phi_2}^{\mathrm{not}}. In Eq. (19), TBT_B denotes partial transposition on Bob’s subsystem [38, 39] and D(ρ∥σ)=Tr⁡[ρ(log⁡2ρ−log⁡2σ)]D(\rho\mathbin\Vert\sigma)=\operatorname{Tr}[\rho(\log_2\rho-\log_2\sigma)] is the quantum relative entropy, with value +∞+\infty when supp⁡ρ⊄supp⁡σ\operatorname{supp}\rho\not\subset\operatorname{supp}\sigma.

For pure and maximally correlated states, ED(ρ)=Rent(ρ)E_D(\rho)=R_{\mathrm{ent}}(\rho) [20], so for these states every rate above distillable entanglement leads to exponential decay of fidelity with target state Φ2\Phi_2. Our catalyst-resource lower bounds require an exponential fidelity bound at a threshold strictly below the target rate. We formulate this assumption below.

Definition 3 (Exponential Converse Bound). Let ρ\rho be a finite-dimensional bipartite state. We say that non-catalytic distillation from ρ\rho satisfies an exponential-converse bound with parameters rsc≥0r_{\mathrm{sc}}\ge0, γ>0\gamma>0, and N0∈NN_0\in\mathbb{N} if

Fρ→Φ2not(N,L)≤2−γN(20)F_{\rho\to\Phi_2}^{\mathrm{not}}(N,L)\le2^{-\gamma N} \tag*{(20)}

for every integer N≥N0N\ge N_0 and every integer L≥rscNL\ge r_{\mathrm{sc}}N. The parameters are independent of NN and $L.

We next derive lower bounds on the catalyst entanglement required to maintain nonvanishing fidelity at rates where non-catalytic distillation obeys an exponential converse bound. The following theorem applies at every target rate R>Rent(ρ)R>R_{\mathrm{ent}}(\rho), since its exponential-converse assumption is satisfied by choosing a threshold rate rscr_{\mathrm{sc}} such that Rent(ρ)<rsc<RR_{\mathrm{ent}}(\rho)<r_{\mathrm{sc}}<R.

Theorem 4 (Necessary catalyst entanglement). Let ρ\rho be a finite-dimensional bipartite state satisfying the exponential-converse bound of Definition 3 with parameters rscr_{\mathrm{sc}}, γ\gamma, and N0N_0. Consider a deterministic catalytic LOCC protocol with ρ⊗n\rho^{\otimes n} as input, ωC\omega_C as catalyst such

that it achieves the rate R>rscR > r_{\mathrm{sc}} and the fidelity of this protocol with a target Bell state is greater than a constant f0∈(0,1)f_0 \in(0,1).

Then, for all sufficiently large nn and with no restriction on catalyst correlation, the catalyst’s max-relative entropy of entanglement satisfies

Emax⁡(ωC)≥γn−log⁡2(1/f0).(21)E_{\max}(\omega_C) \ge\gamma n - \log_2(1/f_0). \tag*{(21)}

Furthermore, if the correlations of the joint output state ΩSC\Omega_{SC} are restricted by

I(S:C)Ω≤δn,0<δn⟶0,(22)I(S:C)_{\Omega} \le\delta_n,\qquad0 < \delta_n \longrightarrow0, \tag*{(22)}

then

Emax⁡(ωC)=Ω(nlog⁡21δn).(23)E_{\max}(\omega_C) = \Omega\left(n\log_2\frac{1}{\delta_n}\right). \tag*{(23)}

The proof is given in § 1.4.

We next translate the lower bounds on Emax⁡(ωC)E_{\max}(\omega_C) into lower bounds on the local catalyst dimensions. Let

Dn:=min⁡{dim⁡CA,dim⁡CB}(24)D_n := \min\{\dim C_A,\dim C_B\} \tag*{(24)}

denote the smaller local catalyst dimension. We show the following corollary:

Corollary 5 (Necessary catalyst dimension). Under the assumptions of Theorem 4 and without a correlation constraint, the following lower bound on the catalyst dimension holds:

Dn≥f02γn(25)D_n \ge f_0 2^{\gamma n} \tag*{(25)}

for all sufficiently large nn.

Moreover, if we also assume the correlation constraint I(Sn:Cn)≤δn→0I(S_n:C_n)\le\delta_n\to0, we obtain

log⁡Dn=Ω(nlog⁡1δn).(26)\log D_n = \Omega\left(n\log\frac{1}{\delta_n}\right). \tag*{(26)}

Equivalently, there exists a constant c>0c>0 such that

Dn≥(1δn)cn(27)D_n \ge\left(\frac{1}{\delta_n}\right)^{cn} \tag*{(27)}

for all sufficiently large $n.

Proof. The max-relative entropy of entanglement satisfies the standard dimension bound [36, 40]

Emax⁡(ωC)≤log⁡2min⁡{dim⁡CA,dim⁡CB}.(28)E_{\max}(\omega_C) \le\log_2\min\{\dim C_A,\dim C_B\}. \tag*{(28)}

Consequently,

log⁡2Dn≥Emax⁡(ωC).(29)\log_2 D_n \ge E_{\max}(\omega_C). \tag*{(29)}

Combining this inequality with Theorem 4 gives

Dn≥2γn−log⁡2(1/f0)=f02γn(30)D_n \ge2^{\gamma n-\log_2(1/f_0)}=f_0 2^{\gamma n} \tag*{(30)}

for protocols without a constraint on output–catalyst correlations. Under the additional correlation assumptions, Theorem 4 also yields

log⁡2Dn=Ω(nlog⁡21δn),(31)\log_2 D_n = \Omega\left(n\log_2\frac{1}{\delta_n}\right), \tag*{(31)}

which proves the remaining claims. □

Although the above corollary provides a lower bound on the dimension of the catalyst, it does not provide a protocol that uses such a catalyst. The following theorem constructs a protocol which reaches the scaling specified by the lower bound on catalyst dimensions (26) of Corollary 5.

Theorem 6 (Achievable catalyst dimension). Let ρ\rho be a finite-dimensional bipartite state with ED(ρ)>0E_D(\rho)>0. Fix R>ED(ρ)R>E_D(\rho), 0<f0<ED(ρ)R0<f_0<\frac{E_D(\rho)}{R}, and let δn>0\delta_n>0 be any prescribed sequence with δn→0\delta_n\to0.

Then, for all sufficiently large nn, there exists a deterministic correlated catalytic LOCC protocol with input ρ⊗n\rho^{\otimes n} and target Φ2mn\Phi_{2^{m_n}}, where mn=⌊Rn⌋m_n=\lfloor Rn\rfloor, such that

ΩC=ωC,fn:=Tr⁡(Φ2mnΩS)≥f0,I(S:C)Ω≤δn.(32)\Omega_C=\omega_C,\quad f_n:=\operatorname{Tr}(\Phi_{2^{m_n}}\Omega_S)\ge f_0,\quad I(S:C)_{\Omega}\le\delta_n. \tag*{(32)}

Moreover, both local catalyst dimensions satisfy

log⁡2dim⁡CA, log⁡2dim⁡CB=O(nlog⁡21δn).(33)\log_2\dim C_A,\ \log_2\dim C_B=O\left(n\log_2\frac{1}{\delta_n}\right). \tag*{(33)}

The proof is given in Section VI.

Now, let DnoptD_n^{\mathrm{opt}} be the minimum value of min⁡{dim⁡CA,dim⁡CB}\min\{\dim C_A,\dim C_B\} over all deterministic correlated catalytic LOCC protocols with input ρ⊗n\rho^{\otimes n}, target Φ2⌊Rn⌋\Phi_{2^{\lfloor Rn\rfloor}}, fidelity at least f0f_0, and output–catalyst mutual information at most δn\delta_n. Using above results we show the following corollary.

Corollary 7 (Optimal order of catalyst dimension). Under the assumptions of Theorem 6, suppose additionally that the exponential-converse hypothesis of Theorem 4 holds for some threshold strictly below RR. Assume that the output–catalyst mutual information is at most δn\delta_n. Then

log⁡2Dnopt=Θ(nlog⁡21δn).(34)\log_2 D_n^{\mathrm{opt}}=\Theta\left(n\log_2\frac{1}{\delta_n}\right). \tag*{(34)}

Proof. Corollary 5 gives

log⁡2Dnopt=Ω(nlog⁡21δn).(35)\log_2 D_n^{\mathrm{opt}}=\Omega\left(n\log_2\frac{1}{\delta_n}\right). \tag*{(35)}

The matching upper bound follows from Theorem 6. □

The constants in these bounds may depend on ρ\rho, RR, f0f_0, and the fixed exponential-converse parameters, but not on nn or δn\delta_n. The strict inequality f0<ED(ρ)/Rf_0<E_D(\rho)/R provides the fixed positive gap required for our construction. Thus, under the assumptions of Corollary 7, the optimal catalyst dimension satisfies log⁡2Dnopt=Θ(nlog⁡2(1/δn))\log_2 D_n^{\mathrm{opt}}=\Theta(n\log_2(1/\delta_n)) for every fixed fidelity requirement strictly below the optimal asymptotic fidelity. This does not determine the optimal dimension scaling when the fidelity itself converges to ED(ρ)/RE_D(\rho)/R.

Operationally, this corresponds to an optimal local catalyst memory cost of Θ(nlog⁡2(1/δn))\Theta(n\log_{2}(1/\delta_n)) qubits: every admissible protocol requires Ω(nlog⁡2(1/δn))\Omega(n\log_{2}(1/\delta_n)) qubits on each side, and our construction achieves the matching upper bound on both sides. This counts the memory needed to represent the catalyst, not the number of Bell pairs it contains.

Proof of Theorem 2

In this section, we detail the proof of theorem 2. Within the proof, we will use several lemmas whose full statements and proofs can be found in Section VII.

To begin, assume that we have a state ρ\rho such that ED(ρ)>0E_D(\rho)>0. First, we find an upper bound on lim sup⁡n→∞Fρ→Φ2cc(n,mn;∞)\limsup_{n\to\infty}F_{\rho\to\Phi_2}^{\mathrm{cc}}(n,m_n;\infty).

Let a catalytic protocol take ρ⊗n\rho^{\otimes n} as input and produce an mnm_n-pair output of fidelity f(ΩSC,mn)f(\Omega_{SC},m_n). Lemmas 9 and 10 give, for mn≥8m_n\geq8,

(mn−⌈2log⁡mn⌉)f(ΩSC,mn)−2≤ED(ΩS)≤nED(ρ).(36)\left(m_n-\lceil2\log m_n\rceil\right)f(\Omega_{SC},m_n)-2\leq E_D(\Omega_S)\leq nE_D(\rho). \tag*{(36)}

Hence

f(ΩSC,mn)≤nED(ρ)+2mn−⌈2log⁡mn⌉.(37)f(\Omega_{SC},m_n)\leq\frac{nE_D(\rho)+2}{m_n-\lceil2\log m_n\rceil}. \tag*{(37)}

For mn=⌊Rn⌋m_n=\lfloor Rn\rfloor this yields

lim sup⁡n→∞Fρ→Φ2cc(n,mn;∞)≤min⁡{1,ED(ρ)R}.(38)\limsup_{n\to\infty}F_{\rho\to\Phi_2}^{\mathrm{cc}}(n,m_n;\infty)\leq\min\left\{1,\frac{E_D(\rho)}{R}\right\}. \tag*{(38)}

The bound is independent of catalyst dimension, catalyst entanglement, and output–catalyst correlation.

We now construct a protocol that asymptotically attains this upper bound. The protocol consists of two stages. First, correlated catalytic amplification produces a mixture of m~>n\widetilde{m}>n copies of the source state and a separable state. Second, ordinary LOCC distillation with a rate rr is applied to obtain mnm_n Bell pairs. In the following paragraph we describe the first stage.

Fix a rate rr with 0<r<ED(ρ)0<r<E_D(\rho) and a requested rate R>rR>r. For all sufficiently large nn set

mn=⌊Rn⌋,m~=⌈mnr⌉.(39)m_n=\lfloor Rn\rfloor,\qquad\widetilde{m}=\left\lceil\frac{m_n}{r}\right\rceil. \tag*{(39)}

Since r<ED(ρ)r<E_D(\rho) and ⌊rm~⌋≥mn\lfloor r\widetilde{m}\rfloor\geq m_n, there are LOCC distillation channels on m~\widetilde{m} source copies producing mnm_n Bell pairs with fidelity tending to one, after discarding any excess output pairs.

Let τ\tau be a separable bipartite state on the same Hilbert space as ρ\rho. Choose an integer N≥2N\geq2 such that

nNm~∈N,log⁡2(N+1)N≤δn.(40)\frac{nN}{\widetilde{m}}\in\mathbb{N},\qquad\frac{\log_2(N+1)}{N}\leq\delta_n. \tag*{(40)}

Such a finite NN exists because log⁡2(N+1)/N→0\log_2(N+1)/N\to0 along the multiples of m~/gcd⁡(n,m~)\widetilde{m}/\gcd(n,\widetilde{m}). By Proposition 12, there exists a correlated catalytic LOCC protocol with input ρ⊗n\rho^{\otimes n} and output marginal

ΩS=nm~ρ⊗m~+(1−nm~)τ⊗m~,(41)\Omega_S=\frac{n}{\widetilde{m}}\rho^{\otimes\widetilde{m}}+\left(1-\frac{n}{\widetilde{m}}\right)\tau^{\otimes\widetilde{m}}, \tag*{(41)}

satisfying

ΩC=ωC,I(S:C)Ω≤log⁡2(N+1)N≤δn.(42)\Omega_C=\omega_C,\qquad I(S:C)_\Omega\leq\frac{\log_2(N+1)}{N}\leq\delta_n. \tag*{(42)}

Now we proceed to the second stage. We apply the ordinary m~\widetilde{m}-copy distillation channel only to the system output. Let us denote it by DdD_d. This does not change the catalyst marginal and cannot increase I(S:C)I(S:C). Additionally, we have

Tr⁡[Φ2⊗mnDd(ΩS)]=nm~Tr⁡[Φ2⊗mnDd(ρ⊗m~)]+(1−nm~)Tr⁡[Φ2⊗mnDd(τ⊗m~)].(43)\operatorname{Tr}\left[\Phi_2^{\otimes m_n}D_d(\Omega_S)\right]=\frac{n}{\widetilde{m}}\operatorname{Tr}\left[\Phi_2^{\otimes m_n}D_d(\rho^{\otimes\widetilde{m}})\right]+\left(1-\frac{n}{\widetilde{m}}\right)\operatorname{Tr}\left[\Phi_2^{\otimes m_n}D_d(\tau^{\otimes\widetilde{m}})\right]. \tag*{(43)}

Since r<ED(ρ)r<E_D(\rho), ordinary distillation achieves

Tr⁡[Φ2⊗mnDd(ρ⊗m~)]=1−o(1).(44)\operatorname{Tr}\left[\Phi_2^{\otimes m_n}D_d(\rho^{\otimes\widetilde{m}})\right]=1-o(1). \tag*{(44)}

Moreover, since τ\tau is separable and DdD_d is LOCC, the state Dd(τ⊗m~)D_d(\tau^{\otimes\widetilde{m}}) is separable. Thus, by Lemma 8,

0≤Tr⁡[Φ2⊗mnDd(τ⊗m~)]≤2−mn=o(1).(45)0\leq\operatorname{Tr}\left[\Phi_2^{\otimes m_n}D_d(\tau^{\otimes\widetilde{m}})\right]\leq2^{-m_n}=o(1). \tag*{(45)}

Consequently, the achieved fidelity satisfies

Tr⁡[Φ2⊗mnDd(ΩS)]=nm~(1−o(1))+o(1)⟶rR,(46)\operatorname{Tr}\left[\Phi_2^{\otimes m_n}D_d(\Omega_S)\right]=\frac{n}{\widetilde{m}}(1-o(1))+o(1) \longrightarrow\frac{r}{R}, \tag*{(46)}

where we used Eq. (39).

For R>ED(ρ)R>E_D(\rho), this construction therefore gives

lim inf⁡n→∞Fρ→Φ2cc(n,mn;δn)≥rR.(47)\liminf_{n\to\infty}F_{\rho\to\Phi_2}^{\mathrm{cc}}(n,m_n;\delta_n)\geq\frac{r}{R}. \tag*{(47)}

Since this holds for every fixed 0<r<ED(ρ)0<r<E_D(\rho), letting r↑ED(ρ)r\uparrow E_D(\rho) yields

lim inf⁡n→∞Fρ→Φ2cc(n,mn;δn)≥ED(ρ)R.(48)\liminf_{n\to\infty}F_{\rho\to\Phi_2}^{\mathrm{cc}}(n,m_n;\delta_n)\geq\frac{E_D(\rho)}{R}. \tag*{(48)}

For R<ED(ρ)R<E_D(\rho), ordinary unassisted distillation already achieves fidelity tending to one. At R=ED(ρ)R=E_D(\rho), the same catalytic construction gives an asymptotic fidelity of at least r/ED(ρ)r/E_D(\rho) for every fixed r<ED(ρ)r<E_D(\rho). Letting r↑ED(ρ)r\uparrow E_D(\rho) therefore yields a lower bound of one. Combining all three cases, we obtain

lim inf⁡n→∞Fρ→Φ2cc(n,mn;δn)≥min⁡{1,ED(ρ)R}.(49)\liminf_{n\to\infty}F_{\rho\to\Phi_2}^{\mathrm{cc}}(n,m_n;\delta_n)\geq\min\left\{1,\frac{E_D(\rho)}{R}\right\}. \tag*{(49)}

For every R<ED(ρ)R<E_D(\rho), the limiting fidelity is one; for every R>ED(ρ)R>E_D(\rho), it is ED(ρ)/R<1E_D(\rho)/R<1. Hence the reliable rate is ED(ρ)E_D(\rho). For every finite RR, the limiting fidelity is positive because ED(ρ)>0E_D(\rho)>0. Therefore the nonvanishing-fidelity threshold is infinite, proving (13).

Proof of Theorem 4

In this section, we detail the proof of theorem 4. Let fnf_n denote the fidelity of the catalytic protocol with the target state Φ2mn\Phi_{2^{m_n}}. By assumption, we have fn>f0f_n > f_0. Since mn/n→R>rscm_n/n \to R > r_{\mathrm{sc}}, for all sufficiently large nn, we have n≥N0n \ge N_0 and mn≥rscnm_n \ge r_{\mathrm{sc}}n. Lemma 13 and Eq. (20) therefore give

f0≤fn≤2Emax⁡(ωC)Fρ→Φ2not(n,mn)≤2Emax⁡(ωC)−γn.(50)f_0 \le f_n \le2^{E_{\max}(\omega_C)}F_{\rho\to\Phi_2^{\mathrm{not}}}(n,m_n) \le2^{E_{\max}(\omega_C)-\gamma n}. \tag*{(50)}

Taking logarithms proves Eq. (21).

We now repeat the protocol kk times, each time using a fresh input block ρ⊗n\rho^{\otimes n} and the same catalyst. Since the catalyst marginal is returned exactly after each use, the resulting process is a catalytic LOCC protocol with input ρ⊗kn\rho^{\otimes kn}, target kmnkm_n Bell pairs, and catalyst ωC\omega_C. Since

kmn≥rsckn,(51)km_n \ge r_{\mathrm{sc}}kn, \tag*{(51)}

our exponential-converse assumption implies

Fρ→Φ2not(kn,kmn)≤2−γkn.(52)F_{\rho\to\Phi_2^{\mathrm{not}}}(kn,km_n) \le2^{-\gamma kn}. \tag*{(52)}

Now, we impose the additional constraints that the correlations within the joint output state ΩSC\Omega_{SC} are restricted by

I(S:C)Ω≤δn,0<δn⟶0.(53)I(S:C)_{\Omega} \le\delta_n,\qquad0 < \delta_n \longrightarrow0. \tag*{(53)}

Then, by applying Lemma 14, we obtain

f0k≤fnk≤2Emax⁡(ωC)−γkn+ln⁡22(k−1)δn.(54)f_0^k \le f_n^k \le2^{E_{\max}(\omega_C)-\gamma kn}+\sqrt{\frac{\ln2}{2}(k-1)\delta_n}. \tag*{(54)}

Set

kn=⌈ℓn4b⌉,b=log⁡2(1/f0)>0,ℓn=log⁡2(1/δn).(55)k_n=\left\lceil\frac{\ell_n}{4b}\right\rceil,\qquad b=\log_2(1/f_0)>0,\qquad\ell_n=\log_2(1/\delta_n). \tag*{(55)}

Since δn→0\delta_n\to0, eventually kn≥2k_n\ge2. Moreover,

f0kn=2−bkn≥2−ℓn/4=δn1/4.(56)f_0^{k_n}=2^{-bk_n}\ge2^{-\ell_n/4}=\delta_n^{1/4}. \tag*{(56)}

Consequently,

ln⁡22(kn−1)δnf0kn≤ln⁡22knδn1/2≤ln⁡28bℓn2−ℓn/4⟶0,(57)\frac{\sqrt{\frac{\ln2}{2}(k_n-1)\delta_n}}{f_0^{k_n}} \le\sqrt{\frac{\ln2}{2}k_n\delta_n^{1/2}} \le\sqrt{\frac{\ln2}{8b}\ell_n2^{-\ell_n/4}}\longrightarrow0, \tag*{(57)}

where we used the fact that exponential decay dominates polynomial growth. Thus, for all sufficiently large nn,

ln⁡22(kn−1)δn≤12f0kn.(58)\sqrt{\frac{\ln2}{2}(k_n-1)\delta_n}\le\frac{1}{2}f_0^{k_n}. \tag*{(58)}

Substituting k=knk=k_n into Eq. (54) gives

2Emax⁡(ωC)−γknn≥12f0kn.(59)2^{E_{\max}(\omega_C)-\gamma k_n n}\ge\frac{1}{2}f_0^{k_n}. \tag*{(59)}

Taking logarithms yields

Emax⁡(ωC)≥kn(γn−b)−1.(60)E_{\max}(\omega_C)\ge k_n(\gamma n-b)-1. \tag*{(60)}

For sufficiently large nn, kn≥ℓn8bk_n\ge\frac{\ell_n}{8b} and n≥2bγn\ge\frac{2b}{\gamma}. Consequently,

Emax⁡(ωC)≥γ16bnℓn−1.(61)E_{\max}(\omega_C)\ge\frac{\gamma}{16b}n\ell_n-1. \tag*{(61)}

The threshold in nn and the constants in this bound depend only on the fixed parameters and the prescribed sequences (mn)(m_n) and (δn)(\delta_n), not on the protocol or catalyst. Thus, the bound holds uniformly over all admissible protocols. Eq. (23) follows.

Proof of Theorem 6

In this section, we detail the proof of theorem 6. As in the proof of Theorem 2, the protocol consists of two stages. First, correlated catalysis transforms nn copies of ρ\rho into a mixture of ρ⊗m~\rho^{\otimes\tilde m}, with m~>n\tilde m>n, and a separable state. Second, ordinary LOCC distillation is applied to this output to produce mnm_n approximate Bell pairs. To reduce the catalyst memory cost, the first stage uses a truncated binomial distribution instead of the fixed-weight distribution used in Proposition 12.

Choose fixed constants p,rp,r such that

f0<p<rR<ED(ρ)R.(62)f_0<p<\frac{r}{R}<\frac{E_D(\rho)}{R}. \tag*{(62)}

For all sufficiently large nn, set

mn=⌊Rn⌋,m~=⌈mnr⌉>n.(63)m_n=\lfloor Rn\rfloor,\qquad\tilde m=\left\lceil\frac{m_n}{r}\right\rceil>n. \tag*{(63)}

Then

nm~⟶rR>p.(64)\frac{n}{\tilde m}\longrightarrow\frac{r}{R}>p. \tag*{(64)}

Since r<ED(ρ)r<E_D(\rho) and ⌊rm~⌋≥mn\lfloor r\tilde m\rfloor\ge m_n, there exist deterministic LOCC distillation channels Dd\mathcal{D}_d such that

Tr⁡[Φ2⊗mnDd(ρ⊗m~)]=1−εn,εn→0.(65)\operatorname{Tr}\left[\Phi_2^{\otimes m_n}\mathcal{D}_d(\rho^{\otimes\tilde m})\right]=1-\varepsilon_n,\qquad\varepsilon_n\to0. \tag*{(65)}

We now construct the first stage of the protocol. Let τ\tau be a separable state on the same bipartite Hilbert space as ρ\rho, and define

α=ρ⊗n⊗τ⊗(m~−n).(66)\alpha=\rho^{\otimes n}\otimes\tau^{\otimes(\tilde m-n)}. \tag*{(66)}

The additional copies of τ\tau can be prepared freely by LOCC.

For an integer N≥2N\ge2, to be chosen below, set

K=⌊nNm~⌋.(67)K=\left\lfloor\frac{nN}{\tilde m}\right\rfloor. \tag*{(67)}

Let QpQ_p be the probability distribution on {0,1}N\{0,1\}^N given by

Qp(x)=p∣x∣(1−p)N−∣x∣,(68)Q_p(x)=p^{|x|}(1-p)^{N-|x|}, \tag*{(68)}

where ∣x∣=∑i=1Nxi|x| = \sum_{i=1}^{N} x_i is the Hamming weight of the binary string x=(x1,…,xN)x=(x_1,\ldots,x_N), i.e., the number of entries equal to one. Define

Z=∑∣x∣≤KQp(x),P(x)=Qp(x)1{∣x∣≤K}Z,(69)Z=\sum_{|x|\le K}Q_p(x),\qquad P(x)=\frac{Q_p(x)\mathbf{1}_{\{|x|\le K\}}}{Z}, \tag*{(69)}

where

1{∣x∣≤K}:={1,∣x∣≤K,0,∣x∣>K.(70)\mathbf{1}_{\{|x|\le K\}}:=\begin{cases}1,&|x|\le K,\\0,&|x|>K.\end{cases} \tag*{(70)}

Thus, P(x)P(x) is the probability of sampling the string xx. Consider the permutation-invariant state

ΓN=∑x∈{0,1}NP(x)θx1⊗⋯⊗θxN,θ1=ρ⊗m~,θ0=τ⊗m~.(71)\Gamma_N=\sum_{x\in\{0,1\}^N}P(x)\theta_{x_1}\otimes\cdots\otimes\theta_{x_N},\qquad\theta_1=\rho^{\otimes\widetilde{m}},\quad\theta_0=\tau^{\otimes\widetilde{m}}. \tag*{(71)}

The transformation a⊗N→ΓNa^{\otimes N}\to\Gamma_N is deterministic LOCC. Indeed, Alice samples xx according to PP and communicates it to Bob. Every string in the support of PP requires ∣x∣m~≤nN|x|\widetilde{m}\le nN copies of ρ\rho. The parties discard the unused copies of ρ\rho, replace them with copies of τ\tau, and perform matching local permutations to prepare the corresponding tensor product.

Let Γj\Gamma_j denote the jj-block marginal of ΓN\Gamma_N, with Γ0=1\Gamma_0=1 (it is a scalar, see the comment in Lemma 11) on the trivial system. By Lemma 11, the state

ωC=1N∑j=0N−1∣j⟩⟨j∣J⊗a⊗j⊗ΓN−1−j(72)\omega_C=\frac{1}{N}\sum_{j=0}^{N-1}|j\rangle\langle j|_J\otimes a^{\otimes j}\otimes\Gamma_{N-1-j} \tag*{(72)}

enables a correlated catalytic conversion of one fresh block aa into an output with marginal Γ1\Gamma_1. The complete catalyst marginal is returned exactly, and the output–catalyst mutual information at this stage equals

I(S:C)=1ND(ΓN∥Γ1⊗N).(73)I(S:C)=\frac{1}{N}D(\Gamma_N\|\Gamma_1^{\otimes N}). \tag*{(73)}

To bound this quantity, define

σp=pρ⊗m~+(1−p)τ⊗m~.(74)\sigma_p=p\rho^{\otimes\widetilde{m}}+(1-p)\tau^{\otimes\widetilde{m}}. \tag*{(74)}

Let qp(0)=1−pq_p(0)=1-p and qp(1)=pq_p(1)=p, so that Qp(x)=∏i=1Nqp(xi)Q_p(x)=\prod_{i=1}^{N}q_p(x_i). Denote by E\mathcal{E} the classical-to-quantum channel that prepares θx1⊗⋯⊗θxN\theta_{x_1}\otimes\cdots\otimes\theta_{x_N} conditional on the classical string xx. It maps the distribution PP to ΓN\Gamma_N and QpQ_p to σp⊗N\sigma_p^{\otimes N}. Indeed,

E(Qp)=∑x1,…,xN∈{0,1}(∏i=1Nqp(xi))θx1⊗⋯⊗θxN=⨂i=1N(∑xi=01qp(xi)θxi)=((1−p)θ0+pθ1)⊗N=σp⊗N.(75)\begin{aligned} \mathcal{E}(Q_p) &=\sum_{x_1,\ldots,x_N\in\{0,1\}}\left(\prod_{i=1}^{N}q_p(x_i)\right)\theta_{x_1}\otimes\cdots\otimes\theta_{x_N}\\ &=\bigotimes_{i=1}^{N}\left(\sum_{x_i=0}^{1}q_p(x_i)\theta_{x_i}\right)=\bigl((1-p)\theta_0+p\theta_1\bigr)^{\otimes N}\\ &=\sigma_p^{\otimes N}. \tag*{(75)} \end{aligned}

Data processing therefore gives

D(ΓN∥σp⊗N)≤D(P∥Qp)=−log⁡2Z.(76)D(\Gamma_N\|\sigma_p^{\otimes N})\le D(P\|Q_p)=-\log_2 Z. \tag*{(76)}

Furthermore, since every one-block marginal of ΓN\Gamma_N equals Γ1\Gamma_1,

D(ΓN∥σp⊗N)=D(ΓN∥Γ1⊗N)+ND(Γ1∥σp).(77)D(\Gamma_N\|\sigma_p^{\otimes N})=D(\Gamma_N\|\Gamma_1^{\otimes N})+ND(\Gamma_1\|\sigma_p). \tag*{(77)}

Consequently,

1ND(ΓN∥Γ1⊗N)≤−1Nlog⁡2Z.(78)\frac{1}{N}D(\Gamma_N\|\Gamma_1^{\otimes N})\le-\frac{1}{N}\log_2 Z. \tag*{(78)}

Choose a fixed η>0\eta>0 such that, for all sufficiently large nn,

nm~≥p+η.(79)\frac{n}{\widetilde{m}}\ge p+\eta. \tag*{(79)}

For X∼Bin⁡(N,p)X\sim\operatorname{Bin}(N,p), Hoeffding’s inequality yields [41]

1−Z=Pr⁡(X>K)=Pr⁡(XN>nm~)≤e−2η2N.(80)\begin{aligned} 1-Z&=\Pr(X>K)=\Pr\left(\frac{X}{N}>\frac{n}{\widetilde{m}}\right)\\ &\le e^{-2\eta^2N}. \tag*{(80)} \end{aligned}

Whenever 1−Z≤1/21-Z\le1/2, we also have

−log⁡2Z≤2ln⁡2(1−Z).(81)-\log_2Z\le\frac{2}{\ln2}(1-Z). \tag*{(81)}

Combining these inequalities with Eq. (78) gives

1ND(ΓN∥Γ1⊗N)≤2Nln⁡2e−2η2N.(82)\frac{1}{N}D(\Gamma_N\|\Gamma_1^{\otimes N})\le\frac{2}{N\ln2}e^{-2\eta^2N}. \tag*{(82)}

The condition Z≥1/2Z\ge1/2 means that the admissible configurations carry at least half of the probability under the original distribution QpQ_p. Since 1−Z≤e−2η2N1-Z\le e^{-2\eta^2N}, it is sufficient to take

N≥ln⁡22η2.(83)N\ge\frac{\ln2}{2\eta^2}. \tag*{(83)}

As η>0\eta>0 is fixed, this condition is automatically satisfied for sufficiently large NN. It imposes no additional assumption on the source state. Note that ZZ is a normalization factor, not the success probability of the catalytic protocol: the protocol samples directly from PP, whose support contains only admissible configurations.

For all sufficiently large nn, we have δn<1\delta_n<1. We may then choose

N=max⁡{2,⌈12η2ln⁡(2δnln⁡2)⌉}.(84)N=\max\left\{2,\left\lceil\frac{1}{2\eta^2}\ln\left(\frac{2}{\delta_n\ln2}\right)\right\rceil\right\}. \tag*{(84)}

Here ln⁡\ln denotes the natural logarithm. This choice ensures

1ND(ΓN∥Γ1⊗N)≤δn,N=O(log⁡21δn).(85)\frac{1}{N}D(\Gamma_N\|\Gamma_1^{\otimes N})\le\delta_n,\qquad N=O\left(\log_2\frac{1}{\delta_n}\right). \tag*{(85)}

and N→∞N \to\infty.

By permutation invariance,

Γ1=pNρ⊗m~+(1−pN)τ⊗m~,pN=1NE[X∣X≤K].(86)\Gamma_1 = p_N \rho^{\otimes\tilde{m}} + (1-p_N)\tau^{\otimes\tilde{m}}, \qquad p_N = \frac{1}{N}\mathbb{E}[X \mid X \le K]. \tag*{(86)}

Since conditioning on X≤KX \le K removes only an upper tail,

0≤p−pN≤1−Z≤e−2n2N.(87)0 \le p-p_N \le1-Z \le e^{-2n^2N}. \tag*{(87)}

Thus pN→pp_N \to p.

Finally, apply Dd\mathcal{D}_d (recall that Dd\mathcal{D}_d is ordinary distillation protocol) only to the released output block. This preserves the catalyst marginal and cannot increase the output–catalyst mutual information. The final state therefore satisfies

ΩC=ωC,I(S:C)Ω≤δn,ΩS=Dd(Γ1).(88)\Omega_C = \omega_C,\qquad I(S:C)_\Omega\le\delta_n,\qquad\Omega_S = \mathcal{D}_d(\Gamma_1). \tag*{(88)}

Because Dd(τ⊗m~)\mathcal{D}_d(\tau^{\otimes\tilde{m}}) is separable,

0≤Tr⁡[Φ2mnDd(τ⊗m~)]≤2−mn.(89)0 \le\operatorname{Tr}\left[\Phi_{2^{m_n}}\mathcal{D}_d(\tau^{\otimes\tilde{m}})\right] \le2^{-m_n}. \tag*{(89)}

Hence

pN(1−ϵn)≤fn≤pN+(1−pN)2−mn,(90)p_N(1-\epsilon_n) \le f_n \le p_N + (1-p_N)2^{-m_n}, \tag*{(90)}

which implies

fn⟶p>f0.(91)f_n \longrightarrow p > f_0. \tag*{(91)}

In particular, fn≥f0f_n \ge f_0 for all sufficiently large nn. No comparison between ϵn\epsilon_n and δn\delta_n is required.

It remains to bound the catalyst dimensions. Let dA,dBd_A,d_B be the local dimensions of ρ\rho. Each catalyst block contains m~\tilde{m} copies of the source Hilbert space, and the classical clock has NN values. Including the clock, the local dimensions can be chosen to satisfy

log⁡2dim⁡CA≤m~(N−1)log⁡2dA+log⁡2N,(92)\log_2 \dim C_A \le\tilde{m}(N-1)\log_2 d_A+\log_2 N, \tag*{(92)}
log⁡2dim⁡CB≤m~(N−1)log⁡2dB+log⁡2N.(93)\log_2 \dim C_B \le\tilde{m}(N-1)\log_2 d_B+\log_2 N. \tag*{(93)}

Since m~=O(n)\tilde{m}=O(n) and N=O(log⁡2(1/δn))N=O(\log_2(1/\delta_n)), both quantities are O(nlog⁡2(1/δn))O(n\log_2(1/\delta_n)), as claimed.

AUXILIARY LEMMAS AND PROPOSITIONS

Lemma 8 (Separable target overlap [22]). For every separable state σAB\sigma_{AB} with local dimension DD,

Tr⁡(ΦDσ)≤1D.(94)\operatorname{Tr}(\Phi_D\sigma) \le\frac{1}{D}. \tag*{(94)}

Lemma 9 (Distillable entanglement of a finite catalytic output). Assume that ED(ρ)>0E_D(\rho)>0. Suppose a finite-block exact-return catalytic LOCC protocol maps ρ⊗n\rho^{\otimes n} to an output marginal ΩS\Omega_S. Then

ED(ΩS)≤nED(ρ).(95)E_D(\Omega_S) \le nE_D(\rho). \tag*{(95)}

Proof. Fix nn and initial catalyst ωC\omega_C and assume instead that ED(ΩS)>nED(ρ)E_D(\Omega_S)>nE_D(\rho). Consequently, there exists rr such that nED(ρ)<r<ED(ΩS)nE_D(\rho)<r<E_D(\Omega_S). Run tt independent copies of the finite protocol on ρ⊗n\rho^{\otimes n}, using independent copies of its finite catalyst. The system output is exactly ΩS⊗t\Omega_S^{\otimes t} and the complete catalyst marginal is exactly ωC⊗t\omega_C^{\otimes t}. Apply an ordinary distillation protocol only to ΩS⊗t\Omega_S^{\otimes t}. Such protocol does not change the catalyst marginal. Since r<ED(ΩS)r<E_D(\Omega_S) there exist protocol which distill ⌊rt⌋\lfloor rt\rfloor copies of Φ2\Phi_2 state as t→∞t\to\infty. Thus, we constructed a protocol which consumes ntnt copies of a state ρ\rho and distills rtrt copies of Φ2\Phi_2. The rate of this protocol is ⌊rt⌋nt=rn>ED(ρ)\frac{\lfloor rt\rfloor}{nt}=\frac{r}{n}>E_D(\rho) in the limit t→∞t\to\infty. This contradicts the catalytic rate-invariance theorem (Theorem 2 of Ref. [32]). Therefore, (95) holds.

Lemma 10 (Bell-product overlap lower-bounds EDE_D). Let σ\sigma be a state on m≥8m\ge8 qubit pairs and let

f=Tr⁡(Φ2mσ).(96)f=\operatorname{Tr}(\Phi_{2^m}\sigma). \tag*{(96)}

Set

t=⌈2log⁡m⌉,ℓ=m−t>0.(97)t=\lceil2\log m\rceil,\qquad\ell=m-t>0. \tag*{(97)}

Then

ED(σ)≥fℓ−2=f(m−⌈2log⁡m⌉)−2.(98)E_D(\sigma)\ge f\ell-2=f\left(m-\lceil2\log m\rceil\right)-2. \tag*{(98)}

Proof. Let D=2mD=2^m. Apply the isotropic twirl U⊗U‾U\otimes\overline{U}, which is LOCC and preserves ff. Since distillable entanglement cannot increase under LOCC, it is enough to lower-bound the distillable entanglement of

σf=fΦD+u(I−ΦD),u=1−fD2−1.(99)\sigma_f=f\Phi_D+u(I-\Phi_D),\qquad u=\frac{1-f}{D^2-1}. \tag*{(99)}

Decompose each local DD-dimensional space into two subsystems of dimensions d1d_1 and d2d_2 such that:

d1=2t,d2=2ℓ,D=d1d2.(100)d_1=2^t,\qquad d_2=2^\ell,\qquad D=d_1d_2. \tag*{(100)}

The maximally entangled state factors as ΦD=Φd1⊗Φd2\Phi_D=\Phi_{d_1}\otimes\Phi_{d_2}.

Alice and Bob measure their d1d_1-dimensional subsystems in the computational basis. They keep the d2d_2-dimensional subsystems exactly when their outcomes agree and forget the common outcome. The ideal component ΦD\Phi_D always passes and leaves Φd2\Phi_{d_2}. The identity operator has d1d_1 matching outcomes, so the unnormalized retained state is

X=fΦd2+u(d1Id2−Φd2).(101)X=f\Phi_{d_2}+u(d_1I_{d_2}-\Phi_{d_2}). \tag*{(101)}

Its success probability and unnormalized Bell overlap are

s=Tr⁡X=f+u(d1d22−1),(102)s=\operatorname{Tr}X=f+u(d_1d_2^2-1), \tag*{(102)}
h=Tr⁡(Φd2X)=f+u(d1−1).(103)h=\operatorname{Tr}(\Phi_{d_2}X)=f+u(d_1-1). \tag*{(103)}

The normalized successful state is isotropic with local dimension d2d_2, with Bell fidelity

α=hs.(104)\alpha=\frac{h}{s}. \tag*{(104)}

@references 久久精品 Its local marginals are maximally mixed, and therefore its coherent information is

I(A⟩B)X/s=ℓ−h2(α)−(1−α)log⁡(d22−1)(105)I(A\rangle B)_{X/s} = \ell- h_2(\alpha) - (1-\alpha)\log(d_2^2-1) \tag*{(105)}
≥ℓ−1−2ℓ(1−α)(106)\ge\ell- 1 - 2\ell(1-\alpha) \tag*{(106)}
=(2α−1)ℓ−1,(107)= (2\alpha-1)\ell-1, \tag*{(107)}

where h2h_2 denotes the binary entropy. Here, AA and BB denote Alice’s and Bob’s retained subsystems, respectively, each of dimension d2d_2. The hashing inequality states that the positive part of coherent information is an achievable one-way distillation rate [15]. Applying the filter independently to many copies and hashing the successful outputs gives

ED(σ)≥ED(σf)≥smax⁡{0,(2α−1)ℓ−1}≥(2h−s)ℓ−s.(108)\begin{aligned} E_D(\sigma)&\ge E_D(\sigma_f)\\ &\ge s\max\{0,(2\alpha-1)\ell-1\}\ge(2h-s)\ell-s. \tag*{(108)} \end{aligned}

The final inequality is valid even when the expression is negative, because max⁡{0,x}≥x\max\{0,x\}\ge x.

We now estimate the last expression. First, h≥fh\ge f, see Eq. (103). Second,

u(d1d22−1)=(1−f)(d1d22−1)d12d22−1≤1d1=2−t.(109)u(d_1d_2^2-1)=\frac{(1-f)(d_1d_2^2-1)}{d_1^2d_2^2-1}\le\frac{1}{d_1}=2^{-t}. \tag*{(109)}

Plugging it into Eq. (102), we obtain that s≤f+2−ts\le f+2^{-t}. Since also s≤1s\le1,

(2h−s)ℓ−s≥(f−2−t)ℓ−1(110)(2h-s)\ell-s\ge(f-2^{-t})\ell-1 \tag*{(110)}
≥fℓ−2,(111)\ge f\ell-2, \tag*{(111)}

where the second inequality uses the fact that 2−tℓ≤m−2m≤12^{-t}\ell\le m^{-2}m\le1. Recalling (97) completes the proof of (98).

The following construction is a variant of the catalysts introduced in [29, 42].

Lemma 11 (Batch catalyst and entropy identity). Let aa be an arbitrary bipartite state, and suppose that a deterministic LOCC channel transforms a⊗Na^{\otimes N} exactly into a permutation-invariant NN-block state ΓN\Gamma_N. Denote by Γj\Gamma_j its jj-block marginal and set Γ0=1\Gamma_0=1 (it is a scalar, since all the subsystems are traced out).

Define a catalyst of N−1N-1 blocks and a classical clock by

ωC=1N∑j=0N−1∣j⟩⟨j∣J⊗a⊗j⊗ΓN−1−j.(112)\omega_C=\frac{1}{N}\sum_{j=0}^{N-1}|j\rangle\langle j|_J\otimes a^{\otimes j}\otimes\Gamma_{N-1-j}. \tag*{(112)}

Then, there is an exact catalytic conversion of one fresh block aa to an output with marginal Γ1\Gamma_1. Moreover,

I(S:C)Ω=S(Γ1)−1NS(ΓN)=1ND(ΓN∥Γ1⊗N),(113)I(S:C)_\Omega=S(\Gamma_1)-\frac{1}{N}S(\Gamma_N)=\frac{1}{N}D(\Gamma_N\Vert\Gamma_1^{\otimes N}), \tag*{(113)}

where D(ρ∥σ)D(\rho\Vert\sigma) is quantum relative entropy.

Proof. We construct the LOCC protocol that implements the claimed transformation. Apply the following channel. For clock value j<N−1j<N-1, regard one block of the ΓN−1−j\Gamma_{N-1-j} part as the output, place the fresh aa into the bank, and increment the clock. The catalyst marginal in the new clock branch is

a⊗(j+1)⊗ΓN−2−j,(114)a^{\otimes(j+1)}\otimes\Gamma_{N-2-j}, \tag*{(114)}

which is precisely the term of (112) corresponding to j+1j+1. Correlations between the released block and the retained Γ\Gamma blocks are kept; only the marginal of the retained blocks is ΓN−2−j\Gamma_{N-2-j}.

For clock value j=N−1j=N-1, the bank and fresh input contain a⊗Na^{\otimes N}. Apply the assumed exact channel to obtain ΓN\Gamma_N, release one block, retain the remaining N−1N-1 blocks, and reset the clock to zero. The retained marginal is ΓN−1\Gamma_{N-1}, again the required next catalyst branch. Since the clock is uniformly distributed and shifted cyclically, the complete catalyst marginal returns exactly. Every released block has marginal Γ1\Gamma_1 by permutation invariance.

It remains to compute the correlation. Put S(Γ0)=0S(\Gamma_0)=0. Orthogonality of the classical clock labels and standard properties of von Neumann entropy give:

S(ωC)=log⁡N+1N∑j=0N−1(jS(a)+S(ΓN−1−j)),(115)S(\omega_C)=\log N+\frac{1}{N}\sum_{j=0}^{N-1}\bigl(jS(a)+S(\Gamma_{N-1-j})\bigr), \tag*{(115)}
S(ΩSC)=log⁡N+1N∑j=0N−1(jS(a)+S(ΓN−j)).(116)S(\Omega_{SC})=\log N+\frac{1}{N}\sum_{j=0}^{N-1}\bigl(jS(a)+S(\Gamma_{N-j})\bigr). \tag*{(116)}

The second expression is indexed by the final clock value; at clock zero its conditional joint state is ΓN\Gamma_N, and at clock j>0j>0 it is a⊗j⊗ΓN−ja^{\otimes j}\otimes\Gamma_{N-j} with one Γ\Gamma block designated as SS. Since the catalyst returns exactly and ΩS=Γ1\Omega_S=\Gamma_1,

I(S:C)Ω=S(Γ1)+S(ωC)−S(ΩSC)=S(Γ1)−S(ΓN)N,(117)\begin{aligned} I(S:C)_\Omega&=S(\Gamma_1)+S(\omega_C)-S(\Omega_{SC})\\ &=S(\Gamma_1)-\frac{S(\Gamma_N)}{N}, \tag*{(117)} \end{aligned}

where the sums telescope. Finally, all one-block marginals of ΓN\Gamma_N are Γ1\Gamma_1, so

D(ΓN∥Γ1⊗N)=−S(ΓN)−Tr⁡ΓNlog⁡(Γ1⊗N)=NS(Γ1)−S(ΓN).(118)D(\Gamma_N\Vert\Gamma_1^{\otimes N})=-S(\Gamma_N)-\operatorname{Tr}\Gamma_N\log(\Gamma_1^{\otimes N})=NS(\Gamma_1)-S(\Gamma_N). \tag*{(118)}

Dividing by NN and comparing with Eq. (117) proves (113).

Proposition 12

Proposition 12 (Decoupled source-state amplifier). Let ρ\rho be a finite-dimensional bipartite state and let τ\tau be a separable state on the same Hilbert space. Let n,m~,Nn,\widetilde m,N be positive integers satisfying m~>n\widetilde m>n, N≥2N\ge2, and K=nN/m~∈NK=nN/\widetilde m\in\mathbb{N}. There exists a deterministic correlated catalytic LOCC protocol with input ρ⊗n\rho^{\otimes n} and output marginal

ΩS=nm~ρ⊗m~+(1−nm~)τ⊗m~.(119)\Omega_S=\frac{n}{\widetilde m}\rho^{\otimes\widetilde m}+\left(1-\frac{n}{\widetilde m}\right)\tau^{\otimes\widetilde m}. \tag*{(119)}

The catalyst is returned exactly, ΩC=ωC\Omega_C=\omega_C, and the residual correlations satisfy

I(S:C)Ω≤h2(nm~)−1Nlog⁡2(NK)≤log⁡2(N+1)N.(120)I(S:C)_{\Omega}\le h_2\left(\frac{n}{\widetilde{m}}\right)-\frac{1}{N}\log_2\binom{N}{K} \le\frac{\log_2(N+1)}{N}. \tag*{(120)}

Proof. Set

a1=ρ⊗n⊗τ⊗(m~−n),K=nNm~.(121)a_1=\rho^{\otimes n}\otimes\tau^{\otimes(\widetilde{m}-n)},\qquad K=\frac{nN}{\widetilde{m}}. \tag*{(121)}

Since τ\tau is separable, the additional copies of τ\tau can be prepared freely by LOCC. Define

θ1=ρ⊗m~,θ0=τ⊗m~,(122)\theta_1=\rho^{\otimes\widetilde{m}},\qquad\theta_0=\tau^{\otimes\widetilde{m}}, \tag*{(122)}

and the permutation-invariant state

ΓN=1(NK)∑x∈{0,1}N∣x∣=Kθx1⊗⋯⊗θxN,(123)\Gamma_N=\frac{1}{\binom{N}{K}}\sum_{\substack{x\in\{0,1\}^N\\ |x|=K}}\theta_{x_1}\otimes\cdots\otimes\theta_{x_N}, \tag*{(123)}

where ∣x∣=∑i=1Nxi|x|=\sum_{i=1}^{N}x_i. The state a1⊗Na_1^{\otimes N} contains exactly nNnN copies of ρ\rho and (m~−n)N(\widetilde{m}-n)N copies of τ\tau. Every tensor product term in Eq. (123) contains exactly Km~=nNK\widetilde{m}=nN copies of ρ\rho and (N−K)m~=(m~−n)N(N-K)\widetilde{m}=(\widetilde{m}-n)N copies of τ\tau.

Alice chooses a binary string of weight KK uniformly at random and communicates it to Bob. Both parties then perform matching local permutations to arrange the copies of ρ\rho and τ\tau into the corresponding tensor product term in Eq. (123). Consequently, a1⊗N→ΓNa_1^{\otimes N}\to\Gamma_N is a deterministic LOCC transformation.

Since ΓN\Gamma_N is permutation invariant, every one-block marginal is identical. The probability that a given block contains ρ⊗m~\rho^{\otimes\widetilde{m}} equals K/N=n/m~K/N=n/\widetilde{m}. Therefore,

Γ1=nm~ρ⊗m~+(1−nm~)τ⊗m~.(124)\Gamma_1=\frac{n}{\widetilde{m}}\rho^{\otimes\widetilde{m}}+\left(1-\frac{n}{\widetilde{m}}\right)\tau^{\otimes\widetilde{m}}. \tag*{(124)}

By Lemma 11, the transformation a1⊗N→ΓNa_1^{\otimes N}\to\Gamma_N yields a correlated catalytic conversion of a1a_1 into an output with marginal ΩS=Γ1\Omega_S=\Gamma_1. Since the additional copies of τ\tau can be prepared freely, this gives the claimed protocol with input ρ⊗n\rho^{\otimes n} and proves Eq. (119).

By Lemma 11, the residual correlation is

I(S:C)Ω=1ND(ΓN∥Γ1⊗N).(125)I(S:C)_{\Omega}=\frac{1}{N}D\left(\Gamma_N\middle\|\Gamma_1^{\otimes N}\right). \tag*{(125)}

Let p=K/N=n/m~p=K/N=n/\widetilde{m}, and let PN,KP_{N,K} denote the uniform probability distribution over binary strings of length NN and weight KK. Let Qp=Bern⁡(p)⊗NQ_p=\operatorname{Bern}(p)^{\otimes N} be the distribution of NN independent Bernoulli random variables. Consider the classical-to-quantum channel E\mathcal{E} that prepares θx1⊗⋯⊗θxN\theta_{x_1}\otimes\cdots\otimes\theta_{x_N} conditioned on a classical string x=(x1,…,xN)x=(x_1,\ldots,x_N). By construction,

E(PN,K)=ΓN,E(Qp)=Γ1⊗N.(126)\mathcal{E}(P_{N,K})=\Gamma_N,\qquad\mathcal{E}(Q_p)=\Gamma_1^{\otimes N}. \tag*{(126)}

The data-processing inequality for quantum relative entropy therefore gives

D(ΓN∥Γ1⊗N)≤D(PN,K∥Qp)=−log⁡2[(NK)pK(1−p)N−K].(127)D\left(\Gamma_N\middle\|\Gamma_1^{\otimes N}\right)\le D(P_{N,K}\|Q_p) =-\log_2\left[\binom{N}{K}p^K(1-p)^{N-K}\right]. \tag*{(127)}

Finally, define

qK=(NK)pK(1−p)N−K.(128)q_K=\binom{N}{K}p^K(1-p)^{N-K}. \tag*{(128)}

Since K=NpK=Np is a mode of the binomial distribution, and the distribution has N+1N+1 possible outcomes,

qK≥1N+1.(129)q_K\ge\frac{1}{N+1}. \tag*{(129)}

Consequently, by Eq. (127),

I(S:C)Ω≤−1Nlog⁡2qK≤log⁡2(N+1)N,(130)I(S:C)_{\Omega}\le-\frac{1}{N}\log_2q_K\le\frac{\log_2(N+1)}{N}, \tag*{(130)}

which proves Eq. (120).

Auxiliary-state bound

Lemma 13 (Auxiliary-state bound). Let Λ\Lambda be a deterministic LOCC channel acting on ρ⊗n⊗ωC\rho^{\otimes n}\otimes\omega_C, where ωC\omega_C is a finite-dimensional bipartite auxiliary state. Write

ΩSC=Λ(ρ⊗n⊗ωC),f=Tr⁡(Φ2mΩS).(131)\Omega_{SC}=\Lambda(\rho^{\otimes n}\otimes\omega_C),\qquad f=\operatorname{Tr}(\Phi_{2^m}\Omega_S). \tag*{(131)}

Then

f≤min⁡{1,2Emax⁡(ωC)Fρ→Φ2not(n,m)}.(132)f\le\min\left\{1,2^{E_{\max}(\omega_C)}F_{\rho\to\Phi_2}^{\mathrm{not}}(n,m)\right\}. \tag*{(132)}

No return condition on the auxiliary state is required.

Proof. Fix η>0\eta>0. By the definition of Emax⁡E_{\max} (see Eq. (17)), there exists a normalized separable state τC\tau_C such that [36]

ωC≤2Emax⁡(ωC)+ητC.(133)\omega_C\le2^{E_{\max}(\omega_C)+\eta}\tau_C. \tag*{(133)}

Tensoring with ρ⊗n\rho^{\otimes n} and applying the positive map Λ\Lambda gives

Λ(ρ⊗n⊗ωC)≤2Emax⁡(ωC)+ηΛ(ρ⊗n⊗τC).(134)\Lambda(\rho^{\otimes n}\otimes\omega_C)\le2^{E_{\max}(\omega_C)+\eta}\Lambda(\rho^{\otimes n}\otimes\tau_C). \tag*{(134)}

Consequently,

f≤2Emax⁡(ωC)+ηTr⁡[(Φ2m⊗IC)Λ(ρ⊗n⊗τC)]≤2Emax⁡(ωC)+ηFρ→Φ2not(n,m).(135)f\le2^{E_{\max}(\omega_C)+\eta}\operatorname{Tr}\left[(\Phi_{2^m}\otimes I_C)\Lambda(\rho^{\otimes n}\otimes\tau_C)\right] \le2^{E_{\max}(\omega_C)+\eta}F_{\rho\to\Phi_2}^{\mathrm{not}}(n,m). \tag*{(135)}

Indeed, preparing the separable state τC\tau_C, applying Λ\Lambda, and discarding CC defines a non-catalytic LOCC protocol. Taking η↓0\eta\downarrow0 and using f≤1f\le1 proves the claim. □

Lemma 14 (Repeated-use bound). Under the assumptions of Lemma 13, suppose additionally that

ΩC=ωC,I(S:C)Ω≤δ.(136)\Omega_C=\omega_C,\qquad I(S:C)_{\Omega}\le\delta. \tag*{(136)}

Then, for every integer k≥2k\ge2,

fk≤2Emax⁡(ωC)Fρ→Φ2not(kn,mk)+ln⁡22(k−1)δ.(137)f^k\le2^{E_{\max}(\omega_C)}F_{\rho\to\Phi_2}^{\mathrm{not}}(kn,mk)+\sqrt{\frac{\ln2}{2}(k-1)\delta}. \tag*{(137)}

Proof. Apply the same protocol kk times, using fresh input blocks ρ⊗n\rho^{\otimes n} and the same catalyst, without acting again on previous outputs. Let SiS_i denote the entire output block of use ii, including both parties’ systems, and let ΘS1⋯Sk\Theta_{S_1\cdots S_k} denote the final joint state of all output blocks. Each use produces the same output marginal, so ΘSi=ΩS\Theta_{S_i}=\Omega_S. Note that the total output might be correlated and in general ΘS1⋯Sk≠ΩS⊗k\Theta_{S_1\cdots S_k}\ne\Omega_S^{\otimes k}.

Before each use, the catalyst marginal is ωC\omega_C, and the fresh input is independent of the catalyst and all previous outputs. The reduced input to each use is therefore ρ⊗n⊗ωC\rho^{\otimes n}\otimes\omega_C. Consequently, immediately after each use, the joint marginal of the new output and the catalyst is

ΩSiC=Λ(ρ⊗n⊗ωC)=ΩSC,(138)\Omega_{S_iC}=\Lambda(\rho^{\otimes n}\otimes\omega_C)=\Omega_{SC}, \tag*{(138)}

up to relabeling SS as SiS_i. In particular,

ΩC=ωC,I(Si:C)Ω=I(S:C)Ω≤δ.(139)\Omega_C=\omega_C,\qquad I(S_i:C)_\Omega=I(S:C)_\Omega\le\delta. \tag*{(139)}

Here ΩSiC\Omega_{S_iC} refers to the state immediately after use ii. Since subsequent uses do not act on SiS_i, its final marginal also satisfies ΘSi=ΩS\Theta_{S_i}=\Omega_S.

Fix i<ki<k. Starting immediately after use ii, all future outputs are obtained by processing the catalyst together with fresh independent input blocks. Including these fixed input states in the channel definition, the remaining protocol followed by discarding the final catalyst is a quantum channel

Ti:C⟶Si+1⋯Sk.(140)\mathcal{T}_i:C\longrightarrow S_{i+1}\cdots S_k. \tag*{(140)}

Thus

ΘSi⋯Sk=(id⁡Si⊗Ti)(ΩSiC).(141)\Theta_{S_i\cdots S_k}=(\operatorname{id}_{S_i}\otimes\mathcal{T}_i)(\Omega_{S_iC}). \tag*{(141)}

Data processing gives

I(Si:Si+1⋯Sk)Θ≤I(Si:C)Ω≤δ,(142)I(S_i:S_{i+1}\cdots S_k)_\Theta\le I(S_i:C)_\Omega\le\delta, \tag*{(142)}

where each mutual information is evaluated on the corresponding joint state.

Using ΘSi=ΩS\Theta_{S_i}=\Omega_S and the entropy chain identity, we obtain (recall that DD is quantum relative entropy.)

D(ΘS1⋯Sk∥ΩS⊗k)=∑i=1kS(ΘSi)−S(ΘS1⋯Sk)=∑i=1k−1I(Si:Si+1⋯Sk)Θ≤(k−1)δ.(143)\begin{aligned} D(\Theta_{S_1\cdots S_k}\Vert\Omega_S^{\otimes k}) &=\sum_{i=1}^{k}S(\Theta_{S_i})-S(\Theta_{S_1\cdots S_k})\\ &=\sum_{i=1}^{k-1}I(S_i:S_{i+1}\cdots S_k)_\Theta\le(k-1)\delta. \tag*{(143)} \end{aligned}

In the last step we used Eq. (142). Quantum Pinsker’s inequality therefore implies [35]

12∥Θ−ΩS⊗k∥1≤ln⁡22(k−1)δ.(144)\frac{1}{2}\lVert\Theta-\Omega_S^{\otimes k}\rVert_1\le\sqrt{\frac{\ln2}{2}(k-1)\delta}. \tag*{(144)}

Eq. (131) gives Tr⁡(Φ2m⊗kΩS⊗k)=fk\operatorname{Tr}(\Phi_{2^m}^{\otimes k}\Omega_S^{\otimes k})=f^k. Since 0≤Φ2m⊗k≤I0\le\Phi_{2^m}^{\otimes k}\le I, the trace-distance bound yields

∣fk−Tr⁡(Φ2m⊗kΘS1⋯Sk)∣=∣Tr⁡(Φ2m⊗kΩS⊗k)−Tr⁡(Φ2m⊗kΘS1⋯Sk)∣≤12∥ΩS⊗k−ΘS1⋯Sk∥1(145)\left|f^k-\operatorname{Tr}(\Phi_{2^m}^{\otimes k}\Theta_{S_1\cdots S_k})\right| =\left|\operatorname{Tr}(\Phi_{2^m}^{\otimes k}\Omega_S^{\otimes k})-\operatorname{Tr}(\Phi_{2^m}^{\otimes k}\Theta_{S_1\cdots S_k})\right| \le\frac{1}{2}\lVert\Omega_S^{\otimes k}-\Theta_{S_1\cdots S_k}\rVert_1 \tag*{(145)}

Using the fact that ∣a−b∣≤c⇒a≤b+c|a-b|\le c\Rightarrow a\le b+c we obtain:

fk≤Tr⁡(Φ2m⊗kΘS1⋯Sk)+12∥ΩS⊗k−ΘS1⋯Sk∥1(146)f^k\le\operatorname{Tr}(\Phi_{2^m}^{\otimes k}\Theta_{S_1\cdots S_k})+\frac{1}{2}\lVert\Omega_S^{\otimes k}-\Theta_{S_1\cdots S_k}\rVert_1 \tag*{(146)}

The complete repeated-use process is a deterministic LOCC protocol acting on ρ⊗kn⊗ωC\rho^{\otimes kn}\otimes\omega_C, with the catalyst supplied only once. Since Φ2m⊗k\Phi_{2^m}^{\otimes k} is locally unitarily equivalent to Φ2mk\Phi_{2^{mk}}, the auxiliary-state bound gives

Tr⁡(Φ2m⊗kΘ)≤2Emax⁡(ωC)Fρ→Φ2not(kn,mk).(147)\operatorname{Tr}(\Phi_{2^m}^{\otimes k}\Theta)\le2^{E_{\max}(\omega_C)}F^{\mathrm{not}}_{\rho\to\Phi_2}(kn,mk). \tag*{(147)}

Combining Eqs. (147), (146) and (144), we obtain:

fk≤2Emax⁡(ωC)Fρ→Φ2not(kn,mk)+ln⁡22(k−1)δ.(148)f^k\le2^{E_{\max}(\omega_C)}F^{\mathrm{not}}_{\rho\to\Phi_2}(kn,mk)+\sqrt{\frac{\ln2}{2}(k-1)\delta}. \tag*{(148)}

Conclusions

In this work, we have shown that correlated catalysis removes the strong converse for entanglement distillation from every finite-dimensional distillable bipartite state. At any target rate R>0R>0, the optimal asymptotic fidelity is min⁡{1,ED(ρ)/R}\min\{1,E_D(\rho)/R\}, even under an arbitrarily small positive bound on the output–catalyst mutual information tending to zero. The catalyst is returned exactly for every input blocklength. Consequently, catalysis preserves the reliable distillation rate while allowing nonvanishing fidelity at every finite rate.

An analogous catalytic enhancement occurs in coherence manipulation under time-translation covariant operations [43]. Under suitable conditions on the coherent modes, correlated catalysts can increase the reliable distillation rate from zero to infinity [43]. Our result identifies a different effect for entanglement: for every distillable source, the reliable rate remains ED(ρ)E_D(\rho), while the strong-converse threshold becomes infinite, even with asymptotically vanishing output–catalyst mutual information.

The construction used in the proof of Theorem 2 establishes the optimal asymptotic fidelity, but requires a larger catalyst memory. This motivates our second construction, which reduces memory cost while retaining any fixed fidelity requirement 0<f0<ED(ρ)/R0<f_0<E_D(\rho)/R at rates R>ED(ρ)R>E_D(\rho). For every distillable source, it requires only O(nlog⁡2(1/δn))O(n\log_2(1/\delta_n)) qubits of local catalyst memory while ensuring that the fidelity fnf_n with the target state satisfies fn≥f0f_n\ge f_0 and that I(S:C)Ω≤δn→0I(S:C)_\Omega\le\delta_n\to0. Whenever noncatalytic distillation obeys an exponential strong converse at a threshold below RR, a matching lower bound holds uniformly over all admissible protocols. Under this assumption, the optimal local memory cost for such fixed fidelity requirements is therefore Θ(nlog⁡2(1/δn))\Theta(n \log_2(1/\delta_n)).

The fixed fidelity requirement can be chosen arbitrarily close to ED(ρ)/RE_D(\rho)/R, but the constants in the memory bound and the required input blocklength may depend on this choice. Thus, although the construction of Theorem 2 establishes the optimal fidelity limit, our matching memory bounds do not determine the optimal memory cost when fn→ED(ρ)/Rf_n \to E_D(\rho)/R. Characterizing the joint dependence of this cost on the fidelity gap and the allowed output–catalyst correlations remains an open question.

Acknowledgment

This work was supported by the National Science Centre Poland (Grant No. 2022/46/E/ST2/00115 and 2024/55/B/ST2/01590).

ChatGPT 6 Astra was used during the preparation of this manuscript. The authors checked and verified all the results.

References

References

  1. [1]R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, Quantum entanglement, Rev. Mod. Phys. 81, 865 (2009).DOI
  2. [2]C. H. Bennett, G. Brassard, C. Crépeau, 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).pubmed.ncbi.nlm.nih.gov/10053414
  3. [3]C. H. Bennett and S. J. Wiesner, Communication via one- and two-particle operators on Einstein-Podolsky-Rosen states, Phys. Rev. Lett. 69, 2881 (1992).pubmed.ncbi.nlm.nih.gov/10046665
  4. [4]C. H. Bennett and G. Brassard, Quantum cryptography: Public key distribution and coin tossing, in Proceedings of IEEE International Conference on Computers, Systems and Signal Processing (Bangalore, India, 1984) pp. 175–179.
  5. [5]A. K. Ekert, Quantum cryptography based on Bell’s theorem, Phys. Rev. Lett. 67, 661 (1991).pubmed.ncbi.nlm.nih.gov/10044956
  6. [6]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), Communications in Mathematical Physics 328, 387 (2014).
  7. [7]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).pubmed.ncbi.nlm.nih.gov/10061534
  8. [8]C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. Wootters, Mixed-state entanglement and quantum error correction, Phys. Rev. A 54, 3824 (1996).DOI
  9. [9]D. Deutsch, A. Ekert, R. Jozsa, C. Macchiavello, S. Popescu, and A. Sanpera, Quantum privacy amplification and the security of quantum cryptography over noisy channels, Phys. Rev. Lett. 77, 2818 (1996).DOI
  10. [10]J. Dehaene, M. Van den Nest, B. De Moor, and F. Verstraete, Local permutations of products of Bell states and entanglement distillation, Phys. Rev. A 67, 022310 (2003).DOI
  11. [11]M. Horodecki and P. Horodecki, Reduction criterion of separability and limits for a class of distillation protocols, Phys. Rev. A 59, 4206 (1999).
  12. [12]W. Dür and H. J. Briegel, Entanglement purification and quantum error correction, Reports on Progress in Physics 70, 1381 (2007).
  13. [13]R. Matsumoto, Conversion of a general quantum stabilizer code to an entanglement distillation protocol, Journal of Physics A: Mathematical and General 36, 8113 (2003).
  14. [14]C. H. Bennett, H. J. Bernstein, S. Popescu, and B. Schumacher, Concentrating partial entanglement by local operations, Phys. Rev. A 53, 2046 (1996).pubmed.ncbi.nlm.nih.gov/9913106
  15. [15]I. Devetak and A. Winter, Distillation of secret key and entanglement from quantum states, Proceedings of the Royal Society A 461, 207 (2005).
  16. [16]E. M. Rains, Rigorous treatment of distillable entanglement, Phys. Rev. A 60, 173 (1999).DOI
  17. [17]M. B. Plenio and S. Virmani, An introduction to entanglement measures, Quantum Information and Computation 7, 1 (2007).arxiv.org/abs/quant-ph/0504163
  18. [18]L. Lami and B. Regula, No second law of entanglement manipulation after all, Nature Physics 19, 184 (2023).DOI
  19. [19]M. Hayashi, M. Koashi, K. Matsumoto, F. Morikoshi, and A. Winter, Error exponents for entanglement concentration, Journal of Physics A: Mathematical and General 36, 527 (2003).DOI
  20. [20]K. Fang, X. Wang, M. Tomamichel, and R. Duan, Non-asymptotic entanglement distillation, IEEE Transactions on Information Theory 65, 6454 (2019).
  21. [21]M. Tomamichel, M. M. Wilde, and A. Winter, Strong converse rates for quantum communication, IEEE Transactions on Information Theory 63, 715 (2017).
  22. [22]S. Khatri and M. M. Wilde, Principles of quantum communication theory: A modern approach (2024).
  23. [23]D. Jonathan and M. B. Plenio, Entanglement-assisted local manipulation of pure quantum states, Phys. Rev. Lett. 83, 3566 (1999).DOI
  24. [24]C. Datta, T. V. Kondra, M. Miller, and A. Streltsov, Catalysis of entanglement and other quantum resources, Reports on Progress in Physics 86, 116002 (2023).DOI
  25. [25]P. Lipka-Bartosik, H. Wilming, and N. H. Y. Ng, Catalysis in quantum information theory, Rev. Mod. Phys. 96, 025005 (2024).
  26. [26]J. Eisert and M. Wilkens, Catalysis of entanglement manipulation for mixed states, Phys. Rev. Lett. 85, 437 (2000).pubmed.ncbi.nlm.nih.gov/10991302
  27. [27]M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information: 10th Anniversary Edition (Cambridge University Press, 2010).
  28. [28]H. Wilming, R. Gallego, and J. Eisert, Axiomatic characterization of the quantum relative entropy and free energy, Entropy 19, 241 (2017).arxiv.org/abs/1702.08473
  29. [29]N. Shiraishi and T. Sagawa, Quantum thermodynamics of correlated-catalytic state conversion at small scale, Phys. Rev. Lett. 126, 150502 (2021).DOI
  30. [30]H. Wilmington, Entropy and reversible catalysis, Phys. Rev. Lett. 127, 260402 (2021).
  31. [31]T. V. Kondra, C. Datta, and A. Streltsov, Catalytic transformations of pure entangled states, Phys. Rev. Lett. 127, 150503 (2021).DOI
  32. [32]R. Ganardi, T. V. Kondra, and A. Streltsov, Catalytic and asymptotic equivalence for quantum entanglement, Phys. Rev. Lett. 133, 250201 (2024).DOI
  33. [33]L. Lami, B. Regula, and A. Streltsov, No-go theorem for entanglement distillation using catalysis, Phys. Rev. A 109, L050401 (2024).DOI
  34. [34]M. M. Wilde, Quantum Information Theory (Cambridge University Press, 2013).
  35. [35]F. Hiai, M. Ohya, and M. Tsukada, Sufficiency, KMS condition and relative entropy in von Neumann algebras, Pacific Journal of Mathematics 96, 99 (1981).
  36. [36]N. Datta, Min- and max-relative entropies and a new entanglement monotone, IEEE Transactions on Information Theory 55, 2816 (2009).DOI
  37. [37]E. Rains, A semidefinite program for distillable entanglement, IEEE Transactions on Information Theory 47, 2921 (2001).DOI
  38. [38]M. Horodecki, P. Horodecki, and R. Horodecki, Separability of mixed states: necessary and sufficient conditions, Physics Letters A 223, 1 (1996).DOI
  39. [39]A. Peres, Separability criterion for density matrices, Phys. Rev. Lett. 77, 1413 (1996).DOI
  40. [40]N. Datta, Max-relative entropy of entanglement, alias log robustness, International Journal of Quantum Information 07, 475 (2009).arxiv.org/abs/0807.2536
  41. [41]W. Hoeffding, Probability inequalities for sums of bounded random variables, Journal of the American Statistical Association 58, 13 (1963).
  42. [42]R. Duan, Y. Feng, X. L. Li, and M. Ying, Multiple-copy entanglement transformation and entanglement catalysis, Phys. Rev. A 71, 042319 (2005).DOI
  43. [43]N. Shiraishi and R. Takagi, Arbitrary amplification of quantum coherence in asymptotic and catalytic transformation, Phys. Rev. Lett. 132, 180202 (2024).DOI

Paper details

vlnc:W_k32mfvzq
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Physics and Astronomy

ChatGPT 6 Astra · View statement →

Looking for related papers…
Contents