Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Quantum Channel Simulation and the Channel's Smooth Max-Information

Definition
QAlgorithms_Papers_Independent_Intakethird100_P65

by Gohan · Nov 1, 2026 · Mathlib 0df444a (Lean v4.33.1)

prove2quantum

The objects the formalized statements of Quantum Channel Simulation and the Channel's Smooth Max-Information (arXiv:1807.05354) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • StateBall: PDF p.7: a subnormalized closed purified-distance ball.
  • maxInfo: PDF p.7, Eq.(24): optimize normalized B, keeping the actual A marginal.
  • smoothInfo: PDF p.10, Eq.(45): free smoothing, using each candidate's own marginal.
  • fixedSmoothInfo: PDF p.10, Eq.(46): hatted smoothing fixes the entire reference marginal.
  • referenceSmoothInfo: PDF p.19, Eq.(76): tilde smoothing fixes only the reference in Dmax.
  • maxEntangledInput: PDF p.8: normalized maximally entangled state. Only nonempty A is used.
  • purify: PDF pp.8,10: canonical purification of rho, without a full-rank restriction.
  • channelInfo: PDF p.8, Definition 4: max-information of the normalized Choi state. For CPTP maps on positive-dimensional systems the extended value is finite.
  • halfDiamondDistance: PDF p.3, Eq.(1): HALF the full diamond norm.
  • smoothChannelInfo: PDF p.8, Definition 5: same-dimensional CPTP channel smoothing.
  • BipartiteNSOperation: PDF pp.2--3: no-signalling is defined on the operation, independently of its Choi constraints. Both directional marginal maps are channels.
  • PPTCode: PDF p.3: complete positivity of the transpose-conjugated operation.
  • EntCode: PDF pp.2--3: entanglement-assisted local encoding and decoding with an arbitrary finite shared density. The body is the contraction of these local maps, not a no-signalling or entropy assumption.
  • SimulationAssistance: Assistance classes in PDF pp.3,10.
  • AllowedSimulationCode: Membership in the source's independently operational code classes.
  • wire: PDF pp.2--3: connect Ao to Bi through N by matrix-unit expansion of the operations on arbitrary input X. This is independent of the Choi link formula.
  • minimumError: PDF p.3, Eq.(1): optimize the half-diamond error over actual codes. On the source domain the feasible set is nonempty and errors lie in [0,1].
  • oneUseSimulationCost: PDF p.4, Eq.(4): cost of one use of an integer-dimensional noiseless channel. Positive resource dimensions include m=1.
  • SimulationRate: PDF pp.4,10: operational vanishing-error simulation rates, with rate slack before blocklength. No entropy formula enters this definition.
  • CommunicationRate: PDF p.10: operational vanishing-error communication rates.
  • vanishingErrorSimulationCost: PDF pp.4,10: infimum of operational simulation rates; extended at empty sets.
  • vanishingErrorQuantumCapacity: PDF p.10: supremum of operational communication rates, independently defined.
  • zeroErrorCost: PDF p.7, Eq.(20): zero error is fixed before the blocklength limit. The threshold-rate form avoids assigning a default value to a nonexistent limit.
  • channelMutualInfo: PDF p.10, Eq.(43): optimize over ALL input densities, including singular ones.
  • IteratedLimit: PDF pp.10--12: actual successive limits, with positive blocklength and right-hand epsilon limit. In particular no liminf or totalized limit is substituted.
  • Constant: PDF p.9, Eq.(35): all constant channels, not only full-support outputs.
  • smoothChannelDmax: PDF p.9, Eq.(37): smooth N over CPTP maps, with fixed reference M.
  • constantChannelRobustness: PDF p.9, Eq.(38): independently defined generalized robustness. Mixing noise is arbitrary CPTP; only the resulting convex mixture must be constant.
  • smoothRobustness: PDF p.10, Eq.(39): CPTP half-diamond smoothing of generalized robustness.
  • constantRelativeRadius: PDF p.15, Eq.(68): all block constant references, including correlated outputs.
  • constantMaxRadius: PDF p.15: smooth max-divergence radius over all constant references.
  • roundedInfoCost: PDF pp.6,9: rounding is computed from the real information value.
  • roundingCorrection: PDF p.9, Eq.(32): the least rounding correction, not a free variable.
  • roundedTraceCost: PDF p.6, Eq.(17): root, then ceiling, then binary logarithm.
  • linkChoi: PDF p.4, Eq.(6d): fully transposed resource Choi linked with code Choi; entries implement all tensor reorderings explicitly.
  • NSChoi: PDF p.4, Eqs.(6e)--(6g): both NS marginals in the fixed register order.
  • codeBobTranspose: PDF p.4: transpose Bob input AND output in the code Choi matrix.
  • generalErrorFeasible: PDF p.4, SDP (6): a fully bound matrix-variable feasible set.
  • reducedErrorFeasible: PDF p.5, SDP (9): noiseless reduction, with m=1 retained.
  • traceFeasible: PDF p.6, SDP (17): optimize Tr V before rounding.
  • traceValue: PDF p.6: real SDP value; the theorem domain has a nonempty, bounded-below feasible set with objective at least one.
  • zeroTraceValue: PDF p.7, Eq.(19): independent zero-error primal objective.
  • basisDepolarizing: PDF p.13: Depolarizing map extended linearly to all operators.
  • ZeroErrorWitnessPair: PDF pp.13--14, Eqs.(57)--(62): explicit primal and dual feasibility, with both objectives equal to the independently defined zero-error primal optimum.
  • depolarizingPrimalWitness: PDF p.13, Eq.(59): constructed depolarizing primal optimizer.
  • dampingPrimalWitness: PDF p.14, Eq.(62): constructed damping primal optimizer; false/true are 0/1.
  • flaggedErasureChannel: PDF p.14, Example 5: orthogonal flag none; some a are the data basis. The trace factor is required on non-normalized operators.
  • classicalMap: PDF p.18: embed a classical stochastic channel as a quantum map.
  • classicalFeasible: PDF p.18, LP (69): error constraints are for every input, not averaged.
  • classicalValue: PDF p.18: classical LP trace objective.
  • depolarizingSectorSum: PDF p.18, Eq.(70): Sector is a SUM, preserving binomial multiplicities.
  • sectorMatrix: PDF p.18, Eq.(71): matrices constructed from the sector coefficients.
  • sectorWeight: PDF p.18, Eq.(72): partial trace weights with k=0 and k=n included.
  • sectorEigenvalue: PDF p.18: eigenvalue of the target Choi matrix on sector k.
  • depolarizingFeasible: PDF p.18, LP (72): concrete finite coefficient feasible set.
  • depolarizingValue: PDF p.18: objective d^n s, before square root and integer ceiling.
  • deFinetti: PDF pp.11,20: the Haar pure-state n-th moment is the normalized symmetric projector on (AA')^n. The finite permutation formula constructs that projector, with rank choose(|A|^2+n-1,n). The same identity is stated in [43], arXiv:0809.3019v1, p.2, proof of Lemma 2.
  • deFinettiPurification: PDF p.20: canonical purification of the constructed de Finetti state, regrouped as (R A^n) A'^n; all of R A^n is the untouched reference.
  • PermutationCovariantChannel: PDF pp.11,20: covariance, equivalently inverse input conjugation. SOURCE NOTATION REPAIR: the printed same-forward pi on both sides does not fix tensor powers for a 3-cycle; the proof uses this covariance relation.
  • deFinettiK1: PDF p.20, Eq.(85): independently constructed fixed-marginal state set.
  • deFinettiK2: PDF p.20, Eq.(86): independently constructed image of covariant CPTP maps.
  • infoOnSet: PDF p.20, Eq.(84): extended infimum on the displayed state set.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_QAlgorithms_Defs_QuantumInformation
import Definitions.Def_QAlgorithms_Defs_MatrixDistances
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P09
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P35
import Definitions.Def_QAlgorithms_Papers_Independent_Independent2_P01
import Definitions.Def_QAlgorithms_Papers_Independent_Gap48_P33
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P06
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P34
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P62
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P63
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P67
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P69
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P71
import Definitions.Def_QAlgorithms_Papers_Independent_Support50_P39

/-! Fang--Wang--Tomamichel--Berta, arXiv:1807.05354v2.
Finite positive-dimensional systems with fixed computational bases. All logarithms
are binary. Exact accepted definitions are selectively imported; distinct simulation objects
have complete local bodies. Paper results have sorry proofs.
The three false printed assertions are recorded in HARD.md. Compilation and
independent moderation remain pending the batch handoff. -/

noncomputable section

open Filter QAlgorithms
open scoped BigOperators Matrix ComplexConjugate ComplexOrder Topology

namespace QAlgorithms.Papers.Independent.Intakethird100.P65

attribute [local instance] Classical.propDecidable Classical.decEq

open QAlgorithms.Papers.Independent.Dependency50.P09 (transposeBob)
open QAlgorithms.Papers.Independent.Dependency50.P35 (maxRelative)
open QAlgorithms.Papers.Independent.Independent2.P01 (dampingKraus damping)
open QAlgorithms.Papers.Independent.Gap48.P33 (phaseFlipChannel)
open QAlgorithms.Papers.Independent.Information100.P06 (Below)
open QAlgorithms.Papers.Independent.Information100.P34 (choi)
open QAlgorithms.Papers.Independent.Information100.P62 (Stochastic)
open QAlgorithms.Papers.Independent.Information100.P63 (Mat)
open QAlgorithms.Papers.Independent.Information100.P67 (IsSubnormalized purifiedDistance)
open QAlgorithms.Papers.Independent.Information100.P69 (replacer permute)
open QAlgorithms.Papers.Independent.Information100.P71 (diamond channelPower channelRelative)
open QAlgorithms.Papers.Independent.Support50.P39 (quantumMutualInformation)

variable {A B C D R : Type} [Fintype A] [DecidableEq A]
  [Fintype B] [DecidableEq B] [Fintype C] [DecidableEq C]
  [Fintype D] [DecidableEq D] [Fintype R] [DecidableEq R]

/-- PDF p.7: a subnormalized closed purified-distance ball. -/
def StateBall (rho sigma : Mat A) (e : ℝ) : Prop :=
  IsSubnormalized sigma ∧ purifiedDistance rho sigma ≤ e

/-- PDF p.7, Eq.(24): optimize normalized B, keeping the actual A marginal. -/
def maxInfo (rho : Mat (A × B)) : EReal :=
  sInf {v : EReal | ∃ sigma : Mat B, IsDensityMatrix sigma ∧
    v = maxRelative rho (Matrix.kronecker (traceRight rho) sigma)}

/-- PDF p.10, Eq.(45): free smoothing, using each candidate's own marginal. -/
def smoothInfo (rho : Mat (A × B)) (e : ℝ) : EReal :=
  sInf {v : EReal | ∃ tau : Mat (A × B), StateBall rho tau e ∧ v = maxInfo tau}

/-- PDF p.10, Eq.(46): hatted smoothing fixes the entire reference marginal. -/
def fixedSmoothInfo (rho : Mat (A × B)) (e : ℝ) : EReal :=
  sInf {v : EReal | ∃ tau : Mat (A × B), StateBall rho tau e ∧
    traceRight tau = traceRight rho ∧ v = maxInfo tau}

/-- PDF p.19, Eq.(76): tilde smoothing fixes only the reference in Dmax. -/
def referenceSmoothInfo (rho : Mat (A × B)) (e : ℝ) : EReal :=
  sInf {v : EReal | ∃ tau : Mat (A × B), StateBall rho tau e ∧
    ∃ sigma : Mat B, IsDensityMatrix sigma ∧
      v = maxRelative tau (Matrix.kronecker (traceRight rho) sigma)}

/-- PDF p.8: normalized maximally entangled state. Only nonempty A is used. -/
def maxEntangledInput (A : Type) [Fintype A] [DecidableEq A] : Mat (A × A) :=
  fun x y => if x.1 = x.2 ∧ y.1 = y.2 then (Fintype.card A : ℂ)⁻¹ else 0

/-- PDF pp.8,10: canonical purification of rho, without a full-rank restriction. -/
def purify (rho : Mat A) : Mat (A × A) :=
  let v : A × A → ℂ := fun x => cfc Real.sqrt rho x.1 x.2
  fun x y => v x * conj (v y)

/-- PDF p.8, Definition 4: max-information of the normalized Choi state.
For CPTP maps on positive-dimensional systems the extended value is finite. -/
def channelInfo (N : Superop A B) : ℝ :=
  (maxInfo (idTensorSuperop (R := A) N (maxEntangledInput A))).toReal

/-- PDF p.3, Eq.(1): HALF the full diamond norm. -/
def halfDiamondDistance (N M : Superop A B) : ℝ := diamond (fun X => N X - M X) / 2

/-- PDF p.8, Definition 5: same-dimensional CPTP channel smoothing. -/
def smoothChannelInfo (N : Superop A B) (e : ℝ) : ℝ :=
  sInf {v : ℝ | ∃ M : Superop A B, IsChannel M ∧ halfDiamondDistance M N ≤ e ∧ v = channelInfo M}

/-- PDF pp.2--3: no-signalling is defined on the operation, independently of
its Choi constraints. Both directional marginal maps are channels. -/
def BipartiteNSOperation (P : Superop (A × C) (B × D)) : Prop :=
  IsChannel P ∧
  (∃ F : Superop A B, IsChannel F ∧ ∀ X, traceRight (P X) = F (traceRight X)) ∧
  (∃ G : Superop C D, IsChannel G ∧ ∀ X, traceLeft (P X) = G (traceLeft X))

/-- PDF p.3: complete positivity of the transpose-conjugated operation. -/
def PPTCode (P : Superop (A × C) (B × D)) : Prop :=
  IsChannel (fun X => transposeBob (P (transposeBob X)))

/-- PDF pp.2--3: entanglement-assisted local encoding and decoding with an
arbitrary finite shared density. The body is the contraction of these local maps,
not a no-signalling or entropy assumption. -/
def EntCode (P : Superop (A × C) (B × D)) : Prop :=
  ∃ e f : ℕ, 0 < e ∧ 0 < f ∧
  ∃ eta : Mat (Fin e × Fin f), IsDensityMatrix eta ∧
  ∃ E : Superop (A × Fin e) B, IsChannel E ∧
  ∃ F : Superop (C × Fin f) D, IsChannel F ∧
  ∀ X b d,
    P X b d = ∑ i, ∑ j, ∑ u, ∑ v, ∑ a, ∑ c, ∑ s, ∑ t,
      X (i, u) (j, v) * eta (a, s) (c, t) *
      E (Matrix.single (i, a) (j, c) 1) b.1 d.1 *
      F (Matrix.single (u, s) (v, t) 1) b.2 d.2

/-- Assistance classes in PDF pp.3,10. -/
inductive SimulationAssistance where
  | ns | nsPpt | ent
  deriving DecidableEq

/-- Membership in the source's independently operational code classes. -/
def AllowedSimulationCode (o : SimulationAssistance) (P : Superop (A × C) (B × D)) : Prop :=
  match o with
  | .ns => BipartiteNSOperation P
  | .nsPpt => BipartiteNSOperation P ∧ PPTCode P
  | .ent => EntCode P

/-- PDF pp.2--3: connect Ao to Bi through N by matrix-unit expansion of the
operations on arbitrary input X. This is independent of the Choi link formula. -/
def wire (N : Superop B C) (P : Superop (A × C) (B × D)) : Superop A D :=
  fun X b d => ∑ a, ∑ c, ∑ u, ∑ v,
    N (Matrix.single a c 1) u v *
      P (Matrix.kronecker X (Matrix.single u v 1)) (a, b) (c, d)

/-- PDF p.3, Eq.(1): optimize the half-diamond error over actual codes.
On the source domain the feasible set is nonempty and errors lie in [0,1]. -/
def minimumError (o : SimulationAssistance) (N : Superop B C) (M : Superop A D) : ℝ :=
  sInf {v : ℝ | ∃ P : Superop (A × C) (B × D), AllowedSimulationCode o P ∧
    v = halfDiamondDistance (wire N P) M}

/-- PDF p.4, Eq.(4): cost of one use of an integer-dimensional noiseless channel.
Positive resource dimensions include m=1. -/
def oneUseSimulationCost (o : SimulationAssistance) (N : Superop A B) (e : ℝ) : ℝ :=
  sInf {v : ℝ | ∃ m : ℕ, 0 < m ∧
    minimumError o (fun X : Mat (Fin m) => X) N ≤ e ∧ v = Real.logb 2 m}

/-- PDF pp.4,10: operational vanishing-error simulation rates, with rate slack
before blocklength. No entropy formula enters this definition. -/
def SimulationRate (o : SimulationAssistance) (N : Superop A B) (r : ℝ) : Prop :=
  0 ≤ r ∧ ∀ e : ℝ, 0 < e → ∀ slack : ℝ, 0 < slack →
    ∃ n0 : ℕ, ∀ n ≥ n0, 0 < n → ∃ m : ℕ, 0 < m ∧
      minimumError o (fun X : Mat (Fin m) => X) (channelPower N n) ≤ e ∧
      Real.logb 2 m ≤ (r + slack) * n

/-- PDF p.10: operational vanishing-error communication rates. -/
def CommunicationRate (o : SimulationAssistance) (N : Superop A B) (r : ℝ) : Prop :=
  0 ≤ r ∧ ∀ e : ℝ, 0 < e → ∀ slack : ℝ, 0 < slack →
    ∃ n0 : ℕ, ∀ n ≥ n0, 0 < n → ∃ m : ℕ, 0 < m ∧
      minimumError o (channelPower N n) (fun X : Mat (Fin m) => X) ≤ e ∧
      (r - slack) * n ≤ Real.logb 2 m

/-- PDF pp.4,10: infimum of operational simulation rates; extended at empty sets. -/
def vanishingErrorSimulationCost (o : SimulationAssistance) (N : Superop A B) : EReal :=
  sInf {v : EReal | ∃ r : ℝ, SimulationRate o N r ∧ v = (r : EReal)}

/-- PDF p.10: supremum of operational communication rates, independently defined. -/
def vanishingErrorQuantumCapacity (o : SimulationAssistance) (N : Superop A B) : EReal :=
  sSup {v : EReal | ∃ r : ℝ, CommunicationRate o N r ∧ v = (r : EReal)}

/-- PDF p.7, Eq.(20): zero error is fixed before the blocklength limit.
The threshold-rate form avoids assigning a default value to a nonexistent limit. -/
def zeroErrorCost (N : Superop A B) : ℝ :=
  sInf {r : ℝ | ∀ slack : ℝ, 0 < slack → ∃ n0 : ℕ, ∀ n ≥ n0, 0 < n →
    oneUseSimulationCost .ns (channelPower N n) 0 / n ≤ r + slack}

/-- PDF p.10, Eq.(43): optimize over ALL input densities, including singular ones. -/
def channelMutualInfo (N : Superop A B) : ℝ :=
  sSup {v : ℝ | ∃ rho : Mat A, IsDensityMatrix rho ∧
    v = quantumMutualInformation (idTensorSuperop N (purify rho))}

/-- PDF pp.10--12: actual successive limits, with positive blocklength and
right-hand epsilon limit. In particular no liminf or totalized limit is substituted. -/
def IteratedLimit (f : ℝ → ℕ → ℝ) (value : ℝ) : Prop :=
  ∃ L : ℝ → ℝ,
    (∀ e : ℝ, 0 < e → e < 1 → Tendsto (fun n => f e (n + 1)) atTop (nhds (L e))) ∧
    Tendsto L (nhdsWithin 0 (Set.Ioi 0)) (nhds value)

/-- PDF p.9, Eq.(35): all constant channels, not only full-support outputs. -/
def Constant (M : Superop A B) : Prop :=
  ∃ sigma : Mat B, IsDensityMatrix sigma ∧ M = replacer sigma

/-- PDF p.9, Eq.(37): smooth N over CPTP maps, with fixed reference M. -/
def smoothChannelDmax (N M : Superop A B) (e : ℝ) : EReal :=
  sInf {v : EReal | ∃ T : Superop A B, IsChannel T ∧
    halfDiamondDistance T N ≤ e ∧ v = maxRelative (choi T) (choi M)}

/-- PDF p.9, Eq.(38): independently defined generalized robustness. Mixing
noise is arbitrary CPTP; only the resulting convex mixture must be constant. -/
def constantChannelRobustness (N : Superop A B) : ℝ :=
  sInf {t : ℝ | 0 ≤ t ∧ ∃ M : Superop A B, IsChannel M ∧
    Constant (fun X => ((1 + t : ℝ)⁻¹ : ℂ) • (N X + (t : ℂ) • M X))}

/-- PDF p.10, Eq.(39): CPTP half-diamond smoothing of generalized robustness. -/
def smoothRobustness (N : Superop A B) (e : ℝ) : ℝ :=
  sInf {v : ℝ | ∃ T : Superop A B, IsChannel T ∧ halfDiamondDistance T N ≤ e ∧ v = constantChannelRobustness T}

/-- PDF p.15, Eq.(68): all block constant references, including correlated outputs. -/
def constantRelativeRadius (N : Superop A B) : EReal :=
  sInf {v : EReal | ∃ M : Superop A B, Constant M ∧ v = channelRelative N M}

/-- PDF p.15: smooth max-divergence radius over all constant references. -/
def constantMaxRadius (N : Superop A B) (e : ℝ) : EReal :=
  sInf {v : EReal | ∃ M : Superop A B, Constant M ∧ v = smoothChannelDmax N M e}

/-- PDF pp.6,9: rounding is computed from the real information value. -/
def roundedInfoCost (x : ℝ) : ℝ := Real.logb 2 (Nat.ceil ((2 : ℝ) ^ (x / 2)))

/-- PDF p.9, Eq.(32): the least rounding correction, not a free variable. -/
def roundingCorrection (x : ℝ) : ℝ := roundedInfoCost x - x / 2

/-- PDF p.6, Eq.(17): root, then ceiling, then binary logarithm. -/
def roundedTraceCost (t : ℝ) : ℝ := Real.logb 2 (Nat.ceil (Real.sqrt t))

/-- PDF p.4, Eq.(6d): fully transposed resource Choi linked with code Choi;
entries implement all tensor reorderings explicitly. -/
def linkChoi (J : Mat (B × C)) (Z : Mat ((A × C) × (B × D))) : Mat (A × D) :=
  fun x y => ∑ a, ∑ c, ∑ u, ∑ v,
    J (a, u) (c, v) * Z ((x.1, u), (a, x.2)) ((y.1, v), (c, y.2))

/-- PDF p.4, Eqs.(6e)--(6g): both NS marginals in the fixed register order. -/
def NSChoi (Z : Mat ((A × C) × (B × D))) : Prop :=
  Z.PosSemidef ∧
  (∀ i j u v, (∑ a, ∑ b, Z ((i,u),(a,b)) ((j,v),(a,b))) =
    if i = j ∧ u = v then 1 else 0) ∧
  (∀ i j u v b d, (∑ a, Z ((i,u),(a,b)) ((j,v),(a,d))) =
    if i = j then (Fintype.card A : ℂ)⁻¹ *
      ∑ k, ∑ a, Z ((k,u),(a,b)) ((k,v),(a,d)) else 0) ∧
  (∀ i j u v a c, (∑ b, Z ((i,u),(a,b)) ((j,v),(c,b))) =
    if u = v then (Fintype.card C : ℂ)⁻¹ *
      ∑ w, ∑ b, Z ((i,w),(a,b)) ((j,w),(c,b)) else 0)

/-- PDF p.4: transpose Bob input AND output in the code Choi matrix. -/
def codeBobTranspose (Z : Mat ((A × C) × (B × D))) : Mat ((A × C) × (B × D)) :=
  fun x y => Z ((x.1.1,y.1.2),(x.2.1,y.2.2)) ((y.1.1,x.1.2),(y.2.1,x.2.2))

/-- PDF p.4, SDP (6): a fully bound matrix-variable feasible set. -/
def generalErrorFeasible (ppt : Bool) (N : Superop B C) (M : Superop A D) (g : ℝ) : Prop :=
  ∃ Y : Mat (A × D), ∃ Z : Mat ((A × C) × (B × D)),
    NSChoi Z ∧ (ppt = true → (codeBobTranspose Z).PosSemidef) ∧
    Below (traceRight Y) ((g : ℂ) • (1 : Mat A)) ∧
    Below (linkChoi (choi N) Z - choi M) Y ∧ Y.PosSemidef

/-- PDF p.5, SDP (9): noiseless reduction, with m=1 retained. -/
def reducedErrorFeasible (ppt : Bool) (N : Superop A B) (m : ℕ) (g : ℝ) : Prop :=
  ∃ Y J : Mat (A × B), ∃ V : Mat B, V.IsHermitian ∧
    Below (traceRight Y) ((g : ℂ) • (1 : Mat A)) ∧ Below (J - choi N) Y ∧ Y.PosSemidef ∧
    J.PosSemidef ∧ traceRight J = 1 ∧ Below J (Matrix.kronecker (1 : Mat A) V) ∧
    V.trace.re = (m : ℝ)^2 ∧
    (ppt = true →
      Below (-Matrix.kronecker (1 : Mat A) V.transpose) ((m : ℂ) • transposeBob J) ∧
      Below ((m : ℂ) • transposeBob J) (Matrix.kronecker (1 : Mat A) V.transpose))

/-- PDF p.6, SDP (17): optimize Tr V before rounding. -/
def traceFeasible (N : Superop A B) (e : ℝ) (Y J : Mat (A × B)) (V : Mat B) : Prop :=
  V.IsHermitian ∧ Below (traceRight Y) ((e : ℂ) • (1 : Mat A)) ∧
  Below (J - choi N) Y ∧ Y.PosSemidef ∧ J.PosSemidef ∧ traceRight J = 1 ∧
  Below J (Matrix.kronecker (1 : Mat A) V)

/-- PDF p.6: real SDP value; the theorem domain has a nonempty, bounded-below
feasible set with objective at least one. -/
def traceValue (N : Superop A B) (e : ℝ) : ℝ :=
  sInf {t : ℝ | ∃ Y J : Mat (A × B), ∃ V : Mat B, traceFeasible N e Y J V ∧ t = V.trace.re}

/-- PDF p.7, Eq.(19): independent zero-error primal objective. -/
def zeroTraceValue (N : Superop A B) : ℝ :=
  sInf {t : ℝ | ∃ V : Mat B, V.IsHermitian ∧
    Below (choi N) (Matrix.kronecker (1 : Mat A) V) ∧ t = V.trace.re}

/-- PDF p.13: Depolarizing map extended linearly to all operators. -/
def basisDepolarizing (A : Type) [Fintype A] [DecidableEq A] (p : ℝ) : Superop A A :=
  fun X => ((1 - p : ℝ) : ℂ) • X + ((p : ℂ) * X.trace / Fintype.card A) • (1 : Mat A)

/-- PDF pp.13--14, Eqs.(57)--(62): explicit primal and dual feasibility,
with both objectives equal to the independently defined zero-error primal optimum. -/
def ZeroErrorWitnessPair (N : Superop A B) (V : Mat B) (X : Mat (A × B)) : Prop :=
  V.IsHermitian ∧ Below (choi N) (Matrix.kronecker (1 : Mat A) V) ∧
  X.PosSemidef ∧ Below (traceLeft X) (1 : Mat B) ∧
  V.trace.re = zeroTraceValue N ∧ ((choi N) * X).trace.re = zeroTraceValue N

/-- PDF p.13, Eq.(59): constructed depolarizing primal optimizer. -/
def depolarizingPrimalWitness (A : Type) [Fintype A] [DecidableEq A] (p : ℝ) : Mat A :=
  (((Fintype.card A : ℝ) * (1-p) + p / (Fintype.card A) : ℝ) : ℂ) • (1 : Mat A)

/-- PDF p.14, Eq.(62): constructed damping primal optimizer; false/true are 0/1. -/
def dampingPrimalWitness (r : ℝ) : Mat Bool :=
  Matrix.diagonal (fun b =>
    if b then ((Real.sqrt (1-r) + 1-r : ℝ) : ℂ)
    else ((1 + Real.sqrt (1-r) : ℝ) : ℂ))

/-- PDF p.14, Example 5: orthogonal flag none; some a are the data basis.
The trace factor is required on non-normalized operators. -/
def flaggedErasureChannel (A : Type) [Fintype A] (p : ℝ) : Superop A (Option A) :=
  fun X x y => match x, y with
  | some a, some b => ((1-p : ℝ) : ℂ) * X a b
  | none, none => (p : ℂ) * X.trace
  | _, _ => 0

/-- PDF p.18: embed a classical stochastic channel as a quantum map. -/
def classicalMap (W : A → B → ℝ) : Superop A B :=
  fun X => Matrix.diagonal (fun y => ∑ x, (W x y : ℂ) * X x x)

/-- PDF p.18, LP (69): error constraints are for every input, not averaged. -/
def classicalFeasible (W : A → B → ℝ) (e : ℝ) (Y T : A → B → ℝ) (V : B → ℝ) : Prop :=
  (∀ x y, T x y - W x y ≤ Y x y ∧ 0 ≤ Y x y ∧ 0 ≤ T x y ∧ T x y ≤ V y) ∧
  (∀ x, ∑ y, T x y = 1) ∧ (∀ x, ∑ y, Y x y ≤ e)

/-- PDF p.18: classical LP trace objective. -/
def classicalValue (W : A → B → ℝ) (e : ℝ) : ℝ :=
  sInf {t : ℝ | ∃ Y T : A → B → ℝ, ∃ V : B → ℝ,
    classicalFeasible W e Y T V ∧ t = ∑ y, V y}

/-- PDF p.18, Eq.(70): Sector is a SUM, preserving binomial multiplicities. -/
def depolarizingSectorSum (A : Type) [Fintype A] [DecidableEq A] (n k : ℕ) :
    Mat ((Fin n → A) × (Fin n → A)) :=
  ∑ S ∈ (Finset.univ : Finset (Fin n)).powersetCard k,
    fun x y => ∏ i : Fin n,
      (if i ∈ S then maxEntangledInput A (x.1 i,x.2 i) (y.1 i,y.2 i)
        else ((1 : Mat (A × A)) - maxEntangledInput A) (x.1 i,x.2 i) (y.1 i,y.2 i))

/-- PDF p.18, Eq.(71): matrices constructed from the sector coefficients. -/
def sectorMatrix (A : Type) [Fintype A] [DecidableEq A] (n : ℕ) (r : Fin (n+1) → ℝ) :
    Mat ((Fin n → A) × (Fin n → A)) :=
  ∑ k, (r k : ℂ) • depolarizingSectorSum A n k.val

/-- PDF p.18, Eq.(72): partial trace weights with k=0 and k=n included. -/
def sectorWeight (d n k : ℕ) : ℝ :=
  (n.choose k : ℝ) * (1 / (d : ℝ)) ^ k * ((d : ℝ) - 1 / d) ^ (n-k)

/-- PDF p.18: eigenvalue of the target Choi matrix on sector k. -/
def sectorEigenvalue (d n k : ℕ) (p : ℝ) : ℝ :=
  ((d : ℝ)*(1-p) + p/d)^k * (p/d)^(n-k)

/-- PDF p.18, LP (72): concrete finite coefficient feasible set. -/
def depolarizingFeasible (d n : ℕ) (p e s : ℝ) (y r : Fin (n+1) → ℝ) : Prop :=
  (∀ k, 0 ≤ y k - r k + sectorEigenvalue d n k.val p ∧ 0 ≤ y k ∧ 0 ≤ r k ∧ r k ≤ s) ∧
  (∑ k, sectorWeight d n k.val * r k) = 1 ∧
  (∑ k, sectorWeight d n k.val * y k) ≤ e

/-- PDF p.18: objective d^n s, before square root and integer ceiling. -/
def depolarizingValue (d n : ℕ) (p e : ℝ) : ℝ :=
  sInf {v : ℝ | ∃ s : ℝ, ∃ y r : Fin (n+1) → ℝ,
    depolarizingFeasible d n p e s y r ∧ v = (d : ℝ)^n * s}

/-- PDF pp.11,20: the Haar pure-state n-th moment is the normalized symmetric
projector on (AA')^n. The finite permutation formula constructs that projector,
with rank choose(|A|^2+n-1,n). The same identity is stated in [43],
arXiv:0809.3019v1, p.2, proof of Lemma 2. -/
def deFinetti (A : Type) [Fintype A] [DecidableEq A] (n : ℕ) : Mat (Fin n → A × A) :=
  fun x y =>
    (∑ p : Equiv.Perm (Fin n), if x = y ∘ p then (1 : ℂ) else 0) /
      ((n.factorial : ℂ) * (Nat.choose ((Fintype.card A)^2 + n - 1) n : ℂ))

/-- PDF p.20: canonical purification of the constructed de Finetti state,
regrouped as (R A^n) A'^n; all of R A^n is the untouched reference. -/
def deFinettiPurification (A : Type) [Fintype A] [DecidableEq A] (n : ℕ) :
    Mat (((Fin n → A × A) × (Fin n → A)) × (Fin n → A)) :=
  let v := fun x : ((Fin n → A × A) × (Fin n → A)) × (Fin n → A) =>
    cfc Real.sqrt (deFinetti A n) x.1.1 (fun i => (x.1.2 i, x.2 i))
  fun x y => v x * conj (v y)

/-- PDF pp.11,20: covariance, equivalently inverse input conjugation.
SOURCE NOTATION REPAIR: the printed same-forward pi on both sides does not
fix tensor powers for a 3-cycle; the proof uses this covariance relation. -/
def PermutationCovariantChannel {n : ℕ} (T : Superop (Fin n → A) (Fin n → B)) : Prop :=
  IsChannel T ∧ ∀ p : Equiv.Perm (Fin n), ∀ X, T (permute p X) = permute p (T X)

/-- PDF p.20, Eq.(85): independently constructed fixed-marginal state set. -/
def deFinettiK1 (N : Superop A B) (n : ℕ) (e : ℝ) :
    Set (Mat (((Fin n → A × A) × (Fin n → A)) × (Fin n → B))) :=
  {sigma | IsDensityMatrix sigma ∧
    purifiedDistance sigma (idTensorSuperop (channelPower N n) (deFinettiPurification A n)) ≤ e ∧
    traceRight sigma = traceRight (deFinettiPurification A n)}

/-- PDF p.20, Eq.(86): independently constructed image of covariant CPTP maps. -/
def deFinettiK2 (N : Superop A B) (n : ℕ) (e : ℝ) :
    Set (Mat (((Fin n → A × A) × (Fin n → A)) × (Fin n → B))) :=
  {sigma | ∃ T : Superop (Fin n → A) (Fin n → B), PermutationCovariantChannel T ∧
    sigma = idTensorSuperop T (deFinettiPurification A n) ∧
    purifiedDistance sigma (idTensorSuperop (channelPower N n) (deFinettiPurification A n)) ≤ e}

/-- PDF p.20, Eq.(84): extended infimum on the displayed state set. -/
def infoOnSet (S : Set (Mat (A × B))) : EReal :=
  sInf {v : EReal | ∃ sigma ∈ S, v = maxInfo sigma}

/-! Sanity examples pin local bodies without proving any paper result. -/

example (X : Mat A) : Below X X := by
  simp only [Below, sub_self]
  exact Matrix.PosSemidef.zero

example (sigma : Mat B) (X : Mat A) : replacer sigma X = X.trace • sigma := rfl

example (N : Superop A B) (a c : A) (b d : B) :
    choi N (a,b) (c,d) = N (Matrix.single a c 1) b d := rfl

example (X : Mat (A × B)) (a c : A) (b d : B) :
    transposeBob X (a,b) (c,d) = X (a,d) (c,b) := rfl

example (X : Mat A) : basisDepolarizing A 0 X = X := by simp [basisDepolarizing]

example (p : ℝ) (X : Mat A) : flaggedErasureChannel A p X none none = (p : ℂ) * X.trace := rfl

example : sectorWeight 2 1 0 = 3/2 ∧ sectorWeight 2 1 1 = 1/2 := by
  norm_num [sectorWeight]

example (p : ℝ) (x y : Fin 0 → A) : channelPower (basisDepolarizing A p) 0
    (1 : Mat (Fin 0 → A)) x y = 1 := by
  simp [channelPower, piSuperop, Matrix.one_apply]

example (x y : Fin 0 → A × A) : deFinetti A 0 x y = 1 := by
  have hxy : x = y := Subsingleton.elim _ _
  subst y
  have hp : ∀ p : Equiv.Perm (Fin 0), x = x ∘ p := fun _ => Subsingleton.elim _ _
  simp [deFinetti, hp]
  have hf : Finset.filter (fun (p : Equiv.Perm (Fin 0)) => x = x ∘ p)
      ({1} : Finset (Equiv.Perm (Fin 0))) = ({1} : Finset (Equiv.Perm (Fin 0))) :=
    Finset.filter_eq_self.mpr (fun p _ => hp p)
  simpa only [Finset.card_singleton] using congrArg Finset.card hf

example : depolarizingPrimalWitness (Fin 1) 0 = 1 := by
  norm_num [depolarizingPrimalWitness]

example : dampingPrimalWitness 1 = Matrix.diagonal (fun b : Bool => if b then 0 else 1) := by
  norm_num [dampingPrimalWitness]

/-! Theorem binders: every system is finite and nonempty. -/

variable [Nonempty A] [Nonempty B] [Nonempty C] [Nonempty D] [Nonempty R]

end QAlgorithms.Papers.Independent.Intakethird100.P65
Source
Quantum Channel Simulation and the Channel's Smooth Max-Information, arXiv:1807.05354v2, https://arxiv.org/abs/1807.05354v2
Read-back

What the Lean code literally says, in plain math · gpt-6.1-sol

QAlgorithms.Papers.Independent.Intakethird100.P65.AllowedSimulationCode

For finite index types A,B,C,DA,B,C,DA,B,C,D, with decidable equality, an assistance label o∈{ns,nsPpt,ent}o\in\{\mathrm{ns},\mathrm{nsPpt},\mathrm{ent}\}o∈{ns,nsPpt,ent}, and an arbitrary matrix map P:A×C→B×DP:A\times C\to B\times DP:A×C→B×D, this predicate is defined by the following three alternatives. For o=nso=\mathrm{ns}o=ns, PPP must be a channel, and there must exist channels F:A→BF:A\to BF:A→B and G:C→DG:C\to DG:C→D such that, for every matrix XXX, Tr⁡DP(X)=F(Tr⁡CX)\operatorname{Tr}_D P(X)=F(\operatorname{Tr}_C X)TrD​P(X)=F(TrC​X) and Tr⁡BP(X)=G(Tr⁡AX)\operatorname{Tr}_B P(X)=G(\operatorname{Tr}_A X)TrB​P(X)=G(TrA​X). For o=nsPpto=\mathrm{nsPpt}o=nsPpt, these same conditions must hold and X↦ΓD(P(ΓCX))X\mapsto\Gamma_D(P(\Gamma_C X))X↦ΓD​(P(ΓC​X)) must also be a channel, where partial transposition on the second factor means (ΓX)(a,c),(a′,c′)=X(a,c′),(a′,c)(\Gamma X)_{(a,c),(a',c')}=X_{(a,c'),(a',c)}(ΓX)(a,c),(a′,c′)​=X(a,c′),(a′,c)​. For o=ento=\mathrm{ent}o=ent, there must exist positive integers e,fe,fe,f, a positive semidefinite trace-one matrix η\etaη on Fin(e)×Fin(f)\mathrm{Fin}(e)\times\mathrm{Fin}(f)Fin(e)×Fin(f), and channels E:A×Fin(e)→BE:A\times\mathrm{Fin}(e)\to BE:A×Fin(e)→B, F:C×Fin(f)→DF:C\times\mathrm{Fin}(f)\to DF:C×Fin(f)→D, such that

P(X)(b,d),(b′,d′)=∑i,j∈A, u,v∈Ca,c<e, s,t<fX(i,u),(j,v)η(a,s),(c,t)E(∣i,a⟩⟨j,c∣)b,b′F(∣u,s⟩⟨v,t∣)d,d′P(X)_{(b,d),(b',d')} =\sum_{\substack{i,j\in A,\ u,v\in C\\a,c<e,\ s,t<f}} X_{(i,u),(j,v)}\eta_{(a,s),(c,t)} E(|i,a\rangle\langle j,c|)_{b,b'} F(|u,s\rangle\langle v,t|)_{d,d'}P(X)(b,d),(b′,d′)​=i,j∈A, u,v∈Ca,c<e, s,t<f​∑​X(i,u),(j,v)​η(a,s),(c,t)​E(∣i,a⟩⟨j,c∣)b,b′​F(∣u,s⟩⟨v,t∣)d,d′​

for all X,b,d,b′,d′X,b,d,b',d'X,b,d,b′,d′. In every alternative, “channel” means a finite Kraus map whose Kraus operators satisfy ∑ℓKℓ†Kℓ=I\sum_\ell K_\ell^\dagger K_\ell=I∑ℓ​Kℓ†​Kℓ​=I.

QAlgorithms.Papers.Independent.Intakethird100.P65.BipartiteNSOperation

For finite types A,B,C,DA,B,C,DA,B,C,D with decidable equality and a function PPP from complex matrices on A×CA\times CA×C to complex matrices on B×DB\times DB×D, this predicate requires that PPP be a finite Kraus channel and that there exist finite Kraus channels FFF from matrices on AAA to matrices on BBB and GGG from matrices on CCC to matrices on DDD such that, for every complex matrix XXX on A×CA\times CA×C,

Tr⁡D(P(X))=F(Tr⁡CX),Tr⁡B(P(X))=G(Tr⁡AX).\operatorname{Tr}_{D}(P(X))=F(\operatorname{Tr}_{C}X), \qquad \operatorname{Tr}_{B}(P(X))=G(\operatorname{Tr}_{A}X).TrD​(P(X))=F(TrC​X),TrB​(P(X))=G(TrA​X).

A finite Kraus channel HHH means that for some natural number rrr and matrices KjK_jKj​ from its input register to its output register, ∑j<rKj†Kj=1\sum_{j<r}K_j^\dagger K_j=1∑j<r​Kj†​Kj​=1 and H(Y)=∑j<rKjYKj†H(Y)=\sum_{j<r}K_jYK_j^\daggerH(Y)=∑j<r​Kj​YKj†​ for every YYY. Partial traces mean summing matching indices of the discarded register; for example, (Tr⁡CX)a,a′=∑cX((a,c),(a′,c))(\operatorname{Tr}_C X)_{a,a'}=\sum_cX((a,c),(a',c))(TrC​X)a,a′​=∑c​X((a,c),(a′,c)). The two identities apply to all matrices, and empty register types and zero-length Kraus families are included whenever the equations permit them.

QAlgorithms.Papers.Independent.Intakethird100.P65.CommunicationRate

For finite index types A,BA,BA,B, with decidable equality, an assistance label ooo, an arbitrary matrix map N:A→BN:A\to BN:A→B, and a real number rrr, this predicate requires r≥0r\ge0r≥0 and states that, for every e>0e>0e>0 and every slack s>0s>0s>0, there exists a natural number n0n_0n0​ such that every positive integer n≥n0n\ge n_0n≥n0​ admits a positive integer mmm satisfying

εo(N[n],id⁡m)≤e,(r−s)n≤log⁡2m.\varepsilon_o(N^{[n]},\operatorname{id}_m)\le e, \qquad (r-s)n\le\log_2 m.εo​(N[n],idm​)≤e,(r−s)n≤log2​m.

The coordinatewise map is defined by N[n](Y)y,y′=∑x,x′Yx,x′∏i<nN(∣xi⟩⟨xi′∣)yi,yi′N^{[n]}(Y)_{y,y'}=\sum_{x,x'}Y_{x,x'}\prod_{i<n}N(|x_i\rangle\langle x_i'|)_{y_i,y_i'}N[n](Y)y,y′​=∑x,x′​Yx,x′​∏i<n​N(∣xi​⟩⟨xi′​∣)yi​,yi′​​. The error εo(N[n],id⁡m)\varepsilon_o(N^{[n]},\operatorname{id}_m)εo​(N[n],idm​) is the real infimum over allowed maps P:Fin(m)×Bn→An×Fin(m)P:\mathrm{Fin}(m)\times B^n\to A^n\times\mathrm{Fin}(m)P:Fin(m)×Bn→An×Fin(m) of

12sup⁡∥X∥1≤1∥(id⁡m⊗(WN[n](P)−id⁡m))(X)∥1,\frac12\sup_{\|X\|_1\le1} \left\|(\operatorname{id}_m\otimes(\mathcal W_{N^{[n]}}(P)-\operatorname{id}_m))(X)\right\|_1,21​∥X∥1​≤1sup​∥(idm​⊗(WN[n]​(P)−idm​))(X)∥1​,

where XXX ranges over all matrices on Fin(m)×Fin(m)\mathrm{Fin}(m)\times\mathrm{Fin}(m)Fin(m)×Fin(m), amplification is blockwise, and WS(P)(Y)b,d=∑a,c,u,vS(∣a⟩⟨c∣)u,vP(Y⊗∣u⟩⟨v∣)(a,b),(c,d)\mathcal W_S(P)(Y)_{b,d}=\sum_{a,c,u,v}S(|a\rangle\langle c|)_{u,v}P(Y\otimes|u\rangle\langle v|)_{(a,b),(c,d)}WS​(P)(Y)b,d​=∑a,c,u,v​S(∣a⟩⟨c∣)u,v​P(Y⊗∣u⟩⟨v∣)(a,b),(c,d)​. Here ∥Z∥1=Re⁡tr⁡Z†Z\|Z\|_1=\operatorname{Re}\operatorname{tr}\sqrt{Z^\dagger Z}∥Z∥1​=RetrZ†Z​. Allowed maps are channels with both marginal nonsignalling conditions for ns\mathrm{ns}ns, such channels whose second-factor partial-transpose conjugates are also channels for nsPpt\mathrm{nsPpt}nsPpt, or maps implemented by separate local channels using a shared finite-dimensional density matrix for ent\mathrm{ent}ent; channels have finite complete Kraus representations. The error infimum need not be attained and equals 000 if empty. No channel assumption is imposed on NNN; zero block length is excluded, and slack larger than rrr is included.

QAlgorithms.Papers.Independent.Intakethird100.P65.Constant

For finite index types A,BA,BA,B with decidable equality and an arbitrary function MMM from complex matrices indexed by AAA to complex matrices indexed by BBB, this property means that there exists a positive semidefinite matrix σ\sigmaσ indexed by BBB with complex trace exactly 111 such that, for every complex matrix XXX indexed by AAA, M(X)=tr⁡(X)σM(X)=\operatorname{tr}(X)\sigmaM(X)=tr(X)σ. There is no separate assumption that MMM is a channel. Empty index types are not excluded: if BBB is empty, no such trace-one σ\sigmaσ exists.

QAlgorithms.Papers.Independent.Intakethird100.P65.EntCode

For finite index types A,B,C,DA,B,C,DA,B,C,D, with decidable equality, and an arbitrary map PPP from complex matrices on A×CA\times CA×C to complex matrices on B×DB\times DB×D, this predicate means that there exist positive natural numbers e,fe,fe,f, a positive semidefinite matrix η\etaη on {0,…,e−1}×{0,…,f−1}\{0,\ldots,e-1\}\times\{0,\ldots,f-1\}{0,…,e−1}×{0,…,f−1} with trace 111, and channels E:A×{0,…,e−1}→BE:A\times\{0,\ldots,e-1\}\to BE:A×{0,…,e−1}→B and F:C×{0,…,f−1}→DF:C\times\{0,\ldots,f-1\}\to DF:C×{0,…,f−1}→D, such that, for every complex matrix XXX on A×CA\times CA×C and every pair of output indices (bB,bD),(dB,dD)(b_B,b_D),(d_B,d_D)(bB​,bD​),(dB​,dD​),

P(X)(bB,bD),(dB,dD)=∑i,j∈A, u,v∈Ca,c<e, s,t<fX(i,u),(j,v)η(a,s),(c,t)E ⁣(∣i,a⟩⟨j,c∣)bB,dBF ⁣(∣u,s⟩⟨v,t∣)bD,dD.P(X)_{(b_B,b_D),(d_B,d_D)} = \sum_{\substack{i,j\in A,\ u,v\in C\\a,c<e,\ s,t<f}} X_{(i,u),(j,v)}\eta_{(a,s),(c,t)} E\!\left(|i,a\rangle\langle j,c|\right)_{b_B,d_B} F\!\left(|u,s\rangle\langle v,t|\right)_{b_D,d_D}.P(X)(bB​,bD​),(dB​,dD​)​=i,j∈A, u,v∈Ca,c<e, s,t<f​∑​X(i,u),(j,v)​η(a,s),(c,t)​E(∣i,a⟩⟨j,c∣)bB​,dB​​F(∣u,s⟩⟨v,t∣)bD​,dD​​.

Here a channel is a map admitting a finite Kraus representation X↦∑ℓKℓXKℓ†X\mapsto\sum_\ell K_\ell X K_\ell^\daggerX↦∑ℓ​Kℓ​XKℓ†​ with ∑ℓKℓ†Kℓ=I\sum_\ell K_\ell^\dagger K_\ell=I∑ℓ​Kℓ†​Kℓ​=I, and ∣x⟩⟨y∣|x\rangle\langle y|∣x⟩⟨y∣ denotes the matrix with its only nonzero entry equal to 111 at (x,y)(x,y)(x,y). The state η\etaη may be entangled or separable; the identity is required for all matrices XXX, without positivity or normalization restrictions.

QAlgorithms.Papers.Independent.Intakethird100.P65.IteratedLimit

For a function f:R×N→Rf:\mathbb R\times\mathbb N\to\mathbb Rf:R×N→R and a real number vvv, this property means that there exists a function L:R→RL:\mathbb R\to\mathbb RL:R→R such that, for every real eee with 0<e<10<e<10<e<1, the sequence f(e,n+1)f(e,n+1)f(e,n+1) converges to L(e)L(e)L(e) as n→∞n\to\inftyn→∞, and L(e)L(e)L(e) converges to vvv as e→0e\to0e→0 through strictly positive real values. The sequence starts with f(e,1)f(e,1)f(e,1); no condition is imposed on f(e,0)f(e,0)f(e,0), on f(e,n)f(e,n)f(e,n) for e∉(0,1)e\notin(0,1)e∈/(0,1), or on the value L(0)L(0)L(0).

QAlgorithms.Papers.Independent.Intakethird100.P65.NSChoi

For finite index types A,B,C,DA,B,C,DA,B,C,D with decidable equality and a complex matrix ZZZ indexed by (A×C)×(B×D)(A\times C)\times(B\times D)(A×C)×(B×D), this property requires that ZZZ be positive semidefinite and that all three following identities hold:

∑a∈B∑b∈DZ((i,u),(a,b)), ((j,v),(a,b))={1i=j and u=v,0otherwise,∑a∈BZ((i,u),(a,b)), ((j,v),(a,d))={∣A∣−1∑k∈A∑a∈BZ((k,u),(a,b)), ((k,v),(a,d))i=j,0i≠j,∑b∈DZ((i,u),(a,b)), ((j,v),(c,b))={∣C∣−1∑w∈C∑b∈DZ((i,w),(a,b)), ((j,w),(c,b))u=v,0u≠v.\begin{aligned} \sum_{a\in B}\sum_{b\in D} Z_{((i,u),(a,b)),\,((j,v),(a,b))} &=\begin{cases}1&i=j\text{ and }u=v,\\0&\text{otherwise},\end{cases}\\ \sum_{a\in B} Z_{((i,u),(a,b)),\,((j,v),(a,d))} &=\begin{cases} |A|^{-1}\displaystyle\sum_{k\in A}\sum_{a\in B} Z_{((k,u),(a,b)),\,((k,v),(a,d))}&i=j,\\ 0&i\ne j, \end{cases}\\ \sum_{b\in D} Z_{((i,u),(a,b)),\,((j,v),(c,b))} &=\begin{cases} |C|^{-1}\displaystyle\sum_{w\in C}\sum_{b\in D} Z_{((i,w),(a,b)),\,((j,w),(c,b))}&u=v,\\ 0&u\ne v. \end{cases} \end{aligned}a∈B∑​b∈D∑​Z((i,u),(a,b)),((j,v),(a,b))​a∈B∑​Z((i,u),(a,b)),((j,v),(a,d))​b∈D∑​Z((i,u),(a,b)),((j,v),(c,b))​​={10​i=j and u=v,otherwise,​=⎩⎨⎧​∣A∣−1k∈A∑​a∈B∑​Z((k,u),(a,b)),((k,v),(a,d))​0​i=j,i=j,​=⎩⎨⎧​∣C∣−1w∈C∑​b∈D∑​Z((i,w),(a,b)),((j,w),(c,b))​0​u=v,u=v.​​

The first identity is required for every i,j∈Ai,j\in Ai,j∈A and u,v∈Cu,v\in Cu,v∈C; the second additionally for every b,d∈Db,d\in Db,d∈D; and the third additionally for every a,c∈Ba,c\in Ba,c∈B. The cardinalities and their inverses are interpreted in C\mathbb CC. No index type is assumed nonempty: empty sums are 000, the inverse of 000 is 000, and universally quantified identities with an empty variable domain are vacuous.

QAlgorithms.Papers.Independent.Intakethird100.P65.PPTCode

For finite types A,B,C,DA,B,C,DA,B,C,D with decidable equality and an arbitrary function PPP from complex matrices on A×CA\times CA×C to complex matrices on B×DB\times DB×D, let TCT_CTC​ and TDT_DTD​ transpose the second register:

(TCX)((a,c),(a′,c′))=X((a,c′),(a′,c)),(TDY)((b,d),(b′,d′))=Y((b,d′),(b′,d)).(T_CX)((a,c),(a',c'))=X((a,c'),(a',c)), \qquad (T_DY)((b,d),(b',d'))=Y((b,d'),(b',d)).(TC​X)((a,c),(a′,c′))=X((a,c′),(a′,c)),(TD​Y)((b,d),(b′,d′))=Y((b,d′),(b′,d)).

This predicate says that X↦TD(P(TCX))X\mapsto T_D(P(T_CX))X↦TD​(P(TC​X)) is a finite Kraus channel: there exist a natural number rrr and matrices Kj∈C(B×D)×(A×C)K_j\in\mathbb C^{(B\times D)\times(A\times C)}Kj​∈C(B×D)×(A×C) such that

∑j<rKj†Kj=1,TD(P(TCX))=∑j<rKjXKj†for every X.\sum_{j<r}K_j^\dagger K_j=1, \qquad T_D(P(T_CX))=\sum_{j<r}K_jXK_j^\dagger \quad\text{for every }X.j<r∑​Kj†​Kj​=1,TD​(P(TC​X))=j<r∑​Kj​XKj†​for every X.

There is no separate requirement here that PPP itself be a channel. Empty register types and r=0r=0r=0 are included whenever these equations hold.

QAlgorithms.Papers.Independent.Intakethird100.P65.PermutationCovariantChannel

For finite register types A,BA,BA,B with decidable equality, a natural number nnn, and a function TTT from complex matrices on AnA^nAn to complex matrices on BnB^nBn, this predicate means that TTT has a finite Kraus representation

T(X)=∑j=0m−1KjXKj†,∑j=0m−1Kj†Kj=IAn,T(X)=\sum_{j=0}^{m-1}K_jXK_j^\dagger, \qquad \sum_{j=0}^{m-1}K_j^\dagger K_j=I_{A^n},T(X)=j=0∑m−1​Kj​XKj†​,j=0∑m−1​Kj†​Kj​=IAn​,

for some natural number mmm and matrices Kj:Bn×An→CK_j:B^n\times A^n\to\mathbb CKj​:Bn×An→C, and that for every permutation π\piπ of the nnn coordinates and every complex matrix XXX,

T(PπX)=Pπ(T(X)),(PπX)(x,y)=X(x∘π,y∘π),T(P_\pi X)=P_\pi(T(X)), \qquad (P_\pi X)(x,y)=X(x\circ\pi,y\circ\pi),T(Pπ​X)=Pπ​(T(X)),(Pπ​X)(x,y)=X(x∘π,y∘π),

where the latter operation uses the appropriate input or output register. The covariance condition applies to all matrices, not just density matrices. No nonemptiness assumption on A,BA,BA,B or positivity assumption on nnn is included; for n=0n=0n=0, the covariance condition is automatic.

QAlgorithms.Papers.Independent.Intakethird100.P65.SimulationAssistance

This declares a type with exactly three distinct values, named ns\mathrm{ns}ns, nsPpt\mathrm{nsPpt}nsPpt, and ent\mathrm{ent}ent, and supplies a decision procedure for equality between its values. This declaration itself imposes no mathematical conditions on maps or states.

QAlgorithms.Papers.Independent.Intakethird100.P65.SimulationRate

For finite index types A,BA,BA,B, with decidable equality, an assistance label ooo, an arbitrary matrix map N:A→BN:A\to BN:A→B, and a real number rrr, this predicate requires r≥0r\ge0r≥0 and states that, for every error tolerance e>0e>0e>0 and every slack s>0s>0s>0, there exists a natural number n0n_0n0​ such that every positive integer n≥n0n\ge n_0n≥n0​ admits a positive integer mmm with

εo(id⁡m,N[n])≤e,log⁡2m≤(r+s)n.\varepsilon_o(\operatorname{id}_m,N^{[n]})\le e, \qquad \log_2 m\le(r+s)n.εo​(idm​,N[n])≤e,log2​m≤(r+s)n.

Here N[n]N^{[n]}N[n] is the map on nnn-coordinate registers defined entrywise by

N[n](Y)y,y′=∑x,x′Yx,x′∏i=0n−1N(∣xi⟩⟨xi′∣)yi,yi′,N^{[n]}(Y)_{y,y'} =\sum_{x,x'}Y_{x,x'}\prod_{i=0}^{n-1} N(|x_i\rangle\langle x_i'|)_{y_i,y_i'},N[n](Y)y,y′​=x,x′∑​Yx,x′​i=0∏n−1​N(∣xi​⟩⟨xi′​∣)yi​,yi′​​,

and εo(id⁡m,N[n])\varepsilon_o(\operatorname{id}_m,N^{[n]})εo​(idm​,N[n]) is the real infimum, over allowed P:An×Fin(m)→Fin(m)×BnP:A^n\times\mathrm{Fin}(m)\to\mathrm{Fin}(m)\times B^nP:An×Fin(m)→Fin(m)×Bn, of one half the trace-norm supremum of the blockwise amplification of Lm(P)−N[n]\mathcal L_m(P)-N^{[n]}Lm​(P)−N[n] on all matrices of trace norm at most 111, with reference register AnA^nAn, where Lm(P)(Y)b,b′=∑u,v<mP(Y⊗∣u⟩⟨v∣)(u,b),(v,b′)\mathcal L_m(P)(Y)_{b,b'}=\sum_{u,v<m}P(Y\otimes|u\rangle\langle v|)_{(u,b),(v,b')}Lm​(P)(Y)b,b′​=∑u,v<m​P(Y⊗∣u⟩⟨v∣)(u,b),(v,b′)​. Allowed maps are channels with both marginal nonsignalling conditions for ns\mathrm{ns}ns; additionally their second-factor partial-transpose conjugates are channels for nsPpt\mathrm{nsPpt}nsPpt; or they arise from a shared finite-dimensional density matrix and separate local channels for ent\mathrm{ent}ent. Channels admit finite complete Kraus representations, and ∥Z∥1=Re⁡tr⁡Z†Z\|Z\|_1=\operatorname{Re}\operatorname{tr}\sqrt{Z^\dagger Z}∥Z∥1​=RetrZ†Z​. The error infimum need not be attained and has value 000 if empty. The threshold may depend on e,se,se,s; n=0n=0n=0 is excluded, while m=1m=1m=1 is permitted.

QAlgorithms.Papers.Independent.Intakethird100.P65.StateBall

For a finite type AAA with decidable equality, complex matrices ρ,σ\rho,\sigmaρ,σ indexed by AAA, and any real eee, this predicate means that σ\sigmaσ is positive semidefinite, Re⁡Tr⁡σ≤1\operatorname{Re}\operatorname{Tr}\sigma\le1ReTrσ≤1, and

1−(F(ρ,σ)+(1−Re⁡Tr⁡ρ)(1−Re⁡Tr⁡σ))2≤e,F(ρ,σ)=Re⁡Tr⁡S ⁣(S(ρ)σS(ρ)),\sqrt{1-\left(F(\rho,\sigma)+ \sqrt{(1-\operatorname{Re}\operatorname{Tr}\rho)(1-\operatorname{Re}\operatorname{Tr}\sigma)}\right)^2}\le e, \qquad F(\rho,\sigma)=\operatorname{Re}\operatorname{Tr} S\!\left(S(\rho)\sigma S(\rho)\right),1−(F(ρ,σ)+(1−ReTrρ)(1−ReTrσ)​)2​≤e,F(ρ,σ)=ReTrS(S(ρ)σS(ρ)),

where S(X)S(X)S(X) denotes applying the real square-root function to XXX by continuous functional calculus. There is no positivity or normalization assumption on ρ\rhoρ; σ=0\sigma=0σ=0 is permitted. Every real square root uses the convention that the square root of a negative number is zero. In particular, a negative eee makes the predicate false.

QAlgorithms.Papers.Independent.Intakethird100.P65.ZeroErrorWitnessPair

For finite types A,BA,BA,B with decidable equality, a map N:CA×A→CB×BN:\mathbb C^{A\times A}\to\mathbb C^{B\times B}N:CA×A→CB×B, a matrix V∈CB×BV\in\mathbb C^{B\times B}V∈CB×B, and a matrix X∈C(A×B)×(A×B)X\in\mathbb C^{(A\times B)\times(A\times B)}X∈C(A×B)×(A×B), this predicate asserts the following six conditions. Define the unnormalized Choi matrix by JN((a,b),(a′,b′))=N(Eaa′)(b,b′)J_N((a,b),(a',b'))=N(E_{aa'})(b,b')JN​((a,b),(a′,b′))=N(Eaa′​)(b,b′), where Eaa′E_{aa'}Eaa′​ has a single entry 111 at (a,a′)(a,a')(a,a′), and define

z(N)=inf⁡{Re⁡tr⁡W:W=W†,IA⊗W−JN⪰0}.z(N)=\inf\left\{\operatorname{Re}\operatorname{tr}W: W=W^\dagger,\quad I_A\otimes W-J_N\succeq0\right\}.z(N)=inf{RetrW:W=W†,IA​⊗W−JN​⪰0}.

Then the required conditions are

V=V†,IA⊗V−JN⪰0,X⪰0,IB−tr⁡AX⪰0,Re⁡tr⁡V=z(N),Re⁡tr⁡(JNX)=z(N),V=V^\dagger,\qquad I_A\otimes V-J_N\succeq0,\qquad X\succeq0,\qquad I_B-\operatorname{tr}_A X\succeq0,\qquad \operatorname{Re}\operatorname{tr}V=z(N),\qquad \operatorname{Re}\operatorname{tr}(J_NX)=z(N),V=V†,IA​⊗V−JN​⪰0,X⪰0,IB​−trA​X⪰0,RetrV=z(N),Retr(JN​X)=z(N),

where (tr⁡AX)(b,b′)=∑a∈AX((a,b),(a,b′))(\operatorname{tr}_A X)(b,b')=\sum_{a\in A}X((a,b),(a,b'))(trA​X)(b,b′)=∑a∈A​X((a,b),(a,b′)). No channel assumption on NNN, positivity assumption on VVV, or normalization assumption on XXX is separately required. Empty finite types are not excluded; empty sums and identity matrices have their usual zero-dimensional meanings. The infimum is the total real-valued infimum used in the declaration, which is 000 for an empty set or a set unbounded below.

QAlgorithms.Papers.Independent.Intakethird100.P65.basisDepolarizing

For every finite type AAA with decidable equality and every real ppp, this defines a map on complex square matrices indexed by AAA by

X⟼(1−p)X+p Tr⁡X∣A∣IA.X\longmapsto (1-p)X+\frac{p\,\operatorname{Tr}X}{|A|}I_A.X⟼(1−p)X+∣A∣pTrX​IA​.

The input XXX is arbitrary, and ppp is not restricted to [0,1][0,1][0,1]. The type AAA may be empty; division by its zero cardinality then has the total-field value zero.

QAlgorithms.Papers.Independent.Intakethird100.P65.channelInfo

For finite types A,BA,BA,B with decidable equality and any function N:CA×A→CB×BN:\mathbb C^{A\times A}\to\mathbb C^{B\times B}N:CA×A→CB×B, write d=∣A∣d=|A|d=∣A∣ and define a matrix RRR indexed by A×BA\times BA×B through

R((a,b),(a′,b′))=N(d−1Ea,a′)b,b′,R((a,b),(a',b'))= N(d^{-1}E_{a,a'})_{b,b'},R((a,b),(a′,b′))=N(d−1Ea,a′​)b,b′​,

where Ea,a′E_{a,a'}Ea,a′​ has entry 111 at (a,a′)(a,a')(a,a′) and zero elsewhere, and d−1d^{-1}d−1 is interpreted in C\mathbb CC. Let RA(a,a′)=∑bR((a,b),(a′,b))R_A(a,a')=\sum_bR((a,b),(a',b))RA​(a,a′)=∑b​R((a,b),(a′,b)). The returned real number is the real conversion of

inf⁡{ℓ∈R | ∃σ∈CB×B, σ⪰0, Tr⁡σ=1, 2ℓ(RA⊗σ)−R⪰0}.\inf\left\{\ell\in\mathbb R\ \middle|\ \exists\sigma\in\mathbb C^{B\times B},\ \sigma\succeq0,\ \operatorname{Tr}\sigma=1,\ 2^\ell(R_A\otimes\sigma)-R\succeq0 \right\}.inf{ℓ∈R ​ ∃σ∈CB×B, σ⪰0, Trσ=1, 2ℓ(RA​⊗σ)−R⪰0}.

This infimum is taken in the extended reals, and either infinite value converts to zero. No linearity or channel assumption is imposed on NNN, and empty types are included.

QAlgorithms.Papers.Independent.Intakethird100.P65.channelMutualInfo

For finite index types A,BA,BA,B, with decidable equality, and an arbitrary matrix map N:A→BN:A\to BN:A→B, this defines the real supremum

sup⁡ρ≥0, tr⁡ρ=1[S(Tr⁡BΩρ)+S(Tr⁡AΩρ)−S(Ωρ)],\sup_{\rho\ge0,\ \operatorname{tr}\rho=1} \left[S(\operatorname{Tr}_B\Omega_\rho) +S(\operatorname{Tr}_A\Omega_\rho)-S(\Omega_\rho)\right],ρ≥0, trρ=1sup​[S(TrB​Ωρ​)+S(TrA​Ωρ​)−S(Ωρ​)],

where vρ(a,a′)=(ρ)a,a′v_\rho(a,a')=(\sqrt{\rho})_{a,a'}vρ​(a,a′)=(ρ​)a,a′​, the matrix used as a purification has entries Πρ(x,y)=vρ(x)vρ(y)‾\Pi_\rho(x,y)=v_\rho(x)\overline{v_\rho(y)}Πρ​(x,y)=vρ​(x)vρ​(y)​, and Ωρ=(id⁡A⊗N)(Πρ)\Omega_\rho=(\operatorname{id}_A\otimes N)(\Pi_\rho)Ωρ​=(idA​⊗N)(Πρ​), meaning (Ωρ)(a,b),(a′,b′)=N((Πρ)(a,⋅),(a′,⋅))b,b′(\Omega_\rho)_{(a,b),(a',b')}=N((\Pi_\rho)_{(a,\cdot),(a',\cdot)})_{b,b'}(Ωρ​)(a,b),(a′,b′)​=N((Πρ​)(a,⋅),(a′,⋅)​)b,b′​. The entropy appearing here is literally S(Z)=−Re⁡tr⁡(cfc[x↦xlog⁡2x](Z))S(Z)=-\operatorname{Re}\operatorname{tr}(\mathrm{cfc}[x\mapsto x\log_2x](Z))S(Z)=−Retr(cfc[x↦xlog2​x](Z)), and both partial traces sum over the discarded index. No channel, linearity, or positivity assumption is imposed on NNN, so Ωρ\Omega_\rhoΩρ​ need not be a density matrix; the displayed total functional-calculus expression is still what is evaluated. The supremum need not be attained. It is taken in the real numbers, where an empty or unbounded-above defining set has supremum 000; in particular, an empty input index type admits no trace-one ρ\rhoρ.

QAlgorithms.Papers.Independent.Intakethird100.P65.classicalFeasible

For an index type AAA, a finite index type BBB, arbitrary real functions W,Y,T:A×B→RW,Y,T:A\times B\to\mathbb RW,Y,T:A×B→R and V:B→RV:B\to\mathbb RV:B→R, and any real eee, this predicate means

∀x∈A, y∈B:T(x,y)−W(x,y)≤Y(x,y),0≤Y(x,y),0≤T(x,y)≤V(y),\forall x\in A,\ y\in B:\quad T(x,y)-W(x,y)\le Y(x,y),\quad 0\le Y(x,y),\quad 0\le T(x,y)\le V(y),∀x∈A, y∈B:T(x,y)−W(x,y)≤Y(x,y),0≤Y(x,y),0≤T(x,y)≤V(y),

together with

∀x∈A:∑y∈BT(x,y)=1,∑y∈BY(x,y)≤e.\forall x\in A:\qquad \sum_{y\in B}T(x,y)=1,\qquad \sum_{y\in B}Y(x,y)\le e.∀x∈A:y∈B∑​T(x,y)=1,y∈B∑​Y(x,y)≤e.

There is no separate assumption on WWW or on the sign of eee. If AAA is empty all these conditions are vacuous, including any restriction on VVV; if AAA is nonempty and BBB is empty, the required row sum 0=10=10=1 makes the predicate impossible.

QAlgorithms.Papers.Independent.Intakethird100.P65.classicalMap

For index types A,BA,BA,B, with AAA finite, and any real function W:A×B→RW:A\times B\to\mathbb RW:A×B→R, this defines a map from complex square matrices on AAA to complex square matrices on BBB by

ΦW(X)y,y′={∑x∈AW(x,y)Xx,x,y=y′,0,y≠y′.\Phi_W(X)_{y,y'}= \begin{cases} \displaystyle\sum_{x\in A}W(x,y)X_{x,x},&y=y',\\ 0,&y\ne y'. \end{cases}ΦW​(X)y,y′​=⎩⎨⎧​x∈A∑​W(x,y)Xx,x​,0,​y=y′,y=y′.​

Thus only diagonal entries of XXX enter the formula. No nonnegativity or sum condition is imposed on WWW, and XXX is arbitrary. Empty index types are allowed; an empty sum is zero.

QAlgorithms.Papers.Independent.Intakethird100.P65.classicalValue

For an index type AAA, a finite index type BBB, any real function W:A×B→RW:A\times B\to\mathbb RW:A×B→R, and any real eee, this defines the real infimum of ∑y∈BV(y)\sum_{y\in B}V(y)∑y∈B​V(y) over real functions Y,T:A×B→RY,T:A\times B\to\mathbb RY,T:A×B→R and V:B→RV:B\to\mathbb RV:B→R satisfying, for every x,yx,yx,y,

T(x,y)−W(x,y)≤Y(x,y),Y(x,y)≥0,0≤T(x,y)≤V(y),T(x,y)-W(x,y)\le Y(x,y),\qquad Y(x,y)\ge0,\qquad 0\le T(x,y)\le V(y),T(x,y)−W(x,y)≤Y(x,y),Y(x,y)≥0,0≤T(x,y)≤V(y),

and, for every xxx,

∑yT(x,y)=1,∑yY(x,y)≤e.\sum_yT(x,y)=1,\qquad \sum_yY(x,y)\le e.y∑​T(x,y)=1,y∑​Y(x,y)≤e.

No stochasticity assumption on WWW, sign assumption on eee, or attainment claim is included. The real infimum defaults to zero for an empty feasible-value set or one unbounded below. In particular, the definition includes empty index types: the constraints are vacuous when AAA is empty, whereas they are impossible when AAA is nonempty and BBB is empty.

QAlgorithms.Papers.Independent.Intakethird100.P65.codeBobTranspose

For arbitrary index types A,B,C,DA,B,C,DA,B,C,D and a complex square matrix ZZZ indexed by (A×C)×(B×D)(A\times C)\times(B\times D)(A×C)×(B×D), this defines a matrix ZΓCDZ^{\Gamma_{CD}}ZΓCD​ on the same index set by

(ZΓCD)((a,c),(b,d)), ((a′,c′),(b′,d′))=Z((a,c′),(b,d′)), ((a′,c),(b′,d)).(Z^{\Gamma_{CD}})_{((a,c),(b,d)),\,((a',c'),(b',d'))} = Z_{((a,c'),(b,d')),\,((a',c),(b',d))}.(ZΓCD​)((a,c),(b,d)),((a′,c′),(b′,d′))​=Z((a,c′),(b,d′)),((a′,c),(b′,d))​.

Thus it transposes the CCC and DDD coordinates together, leaving the AAA and BBB coordinates in place. No positivity, Hermiticity, normalization, or nonemptiness is required.

QAlgorithms.Papers.Independent.Intakethird100.P65.constantChannelRobustness

For finite index types A,BA,BA,B with decidable equality and an arbitrary function NNN from complex AAA-indexed matrices to complex BBB-indexed matrices, this is the real infimum of all real t≥0t\ge0t≥0 for which there exist a channel MMM and a positive semidefinite BBB-indexed matrix σ\sigmaσ with tr⁡σ=1\operatorname{tr}\sigma=1trσ=1 such that, for every complex AAA-indexed matrix XXX,

N(X)+tM(X)1+t=tr⁡(X)σ.\frac{N(X)+tM(X)}{1+t}=\operatorname{tr}(X)\sigma.1+tN(X)+tM(X)​=tr(X)σ.

Here MMM being a channel means that, for some natural number rrr and matrices Kj∈CB×AK_j\in\mathbb C^{B\times A}Kj​∈CB×A, ∑j<rKj†Kj=IA\sum_{j<r}K_j^\dagger K_j=I_A∑j<r​Kj†​Kj​=IA​ and M(X)=∑j<rKjXKj†M(X)=\sum_{j<r}K_jXK_j^\daggerM(X)=∑j<r​Kj​XKj†​ for every XXX. No channel assumption is imposed on NNN. The value t=0t=0t=0 is included, but even then a channel MMM must exist; the denominator 1+t1+t1+t is positive on the feasible set. Empty index types are not excluded. This is a real-valued infimum, whose value on an empty feasible set is 000, and it does not assert attainment.

QAlgorithms.Papers.Independent.Intakethird100.P65.constantMaxRadius

For finite index types A,BA,BA,B with decidable equality, an arbitrary function N:A-matrices→B-matricesN:A\text{-matrices}\to B\text{-matrices}N:A-matrices→B-matrices, and any real eee, this is the extended-real infimum

inf⁡σ⪰0tr⁡σ=1 inf⁡T a channelΔ(T,N)≤e inf⁡{ℓ∈R:2ℓJRσ−JT⪰0},Rσ(X)=tr⁡(X)σ.\inf_{\substack{\sigma\succeq0\\\operatorname{tr}\sigma=1}} \ \inf_{\substack{T\text{ a channel}\\\Delta(T,N)\le e}} \ \inf\{\ell\in\mathbb R:2^\ell J_{R_\sigma}-J_T\succeq0\}, \qquad R_\sigma(X)=\operatorname{tr}(X)\sigma.σ⪰0trσ=1​inf​ T a channelΔ(T,N)≤e​inf​ inf{ℓ∈R:2ℓJRσ​​−JT​⪰0},Rσ​(X)=tr(X)σ.

A channel TTT admits a finite representation T(X)=∑j<rKjXKj†T(X)=\sum_{j<r}K_jXK_j^\daggerT(X)=∑j<r​Kj​XKj†​ with Kj∈CB×AK_j\in\mathbb C^{B\times A}Kj​∈CB×A and ∑j<rKj†Kj=IA\sum_{j<r}K_j^\dagger K_j=I_A∑j<r​Kj†​Kj​=IA​. For any function FFF, its matrix JFJ_FJF​ has entries (JF)(a,b),(a′,b′)=[F(Ea,a′)]b,b′(J_F)_{(a,b),(a',b')}=[F(E_{a,a'})]_{b,b'}(JF​)(a,b),(a′,b′)​=[F(Ea,a′​)]b,b′​, where Ea,a′E_{a,a'}Ea,a′​ is the matrix unit. The distance is Δ(T,N)=12sup⁡∥X∥1≤1∥(idA⊗(T−N))(X)∥1\Delta(T,N)=\tfrac12\sup_{\|X\|_1\le1}\|(\mathrm{id}_A\otimes(T-N))(X)\|_1Δ(T,N)=21​sup∥X∥1​≤1​∥(idA​⊗(T−N))(X)∥1​, with XXX an arbitrary complex matrix indexed by A×AA\times AA×A, ∥Y∥1=Re⁡tr⁡Y†Y\|Y\|_1=\operatorname{Re}\operatorname{tr}\sqrt{Y^\dagger Y}∥Y∥1​=RetrY†Y​, and amplification applying the function to each reference-indexed block. Neither NNN being a channel nor e≥0e\ge0e≥0 is assumed, and empty index types are allowed. Each displayed infimum is interpreted in the extended reals, so an empty feasible set gives +∞+\infty+∞; no attainment is asserted.

QAlgorithms.Papers.Independent.Intakethird100.P65.constantRelativeRadius

For finite index types A,BA,BA,B with decidable equality and an arbitrary function N:A-matrices→B-matricesN:A\text{-matrices}\to B\text{-matrices}N:A-matrices→B-matrices, this is the extended-real infimum

inf⁡σ⪰0tr⁡σ=1 sup⁡ψ∈CA×A∥ψ∥2=1D ⁣((idA⊗N)(∣ψ⟩⟨ψ∣),(idA⊗Rσ)(∣ψ⟩⟨ψ∣)),Rσ(X)=tr⁡(X)σ.\inf_{\substack{\sigma\succeq0\\\operatorname{tr}\sigma=1}} \ \sup_{\substack{\psi\in\mathbb C^{A\times A}\\\|\psi\|_2=1}} D\!\left((\mathrm{id}_A\otimes N)(|\psi\rangle\langle\psi|), (\mathrm{id}_A\otimes R_\sigma)(|\psi\rangle\langle\psi|)\right), \qquad R_\sigma(X)=\operatorname{tr}(X)\sigma.σ⪰0trσ=1​inf​ ψ∈CA×A∥ψ∥2​=1​sup​D((idA​⊗N)(∣ψ⟩⟨ψ∣),(idA​⊗Rσ​)(∣ψ⟩⟨ψ∣)),Rσ​(X)=tr(X)σ.

Here ∣ψ⟩⟨ψ∣|\psi\rangle\langle\psi|∣ψ⟩⟨ψ∣ has entries ψiψj‾\psi_i\overline{\psi_j}ψi​ψj​​, and amplification applies a function separately to each matrix block indexed by the reference coordinates in AAA. The quantity D(R,S)D(R,S)D(R,S) equals +∞+\infty+∞ unless both R,SR,SR,S are positive semidefinite with complex trace 111 and ker⁡S⊆ker⁡R\ker S\subseteq\ker RkerS⊆kerR; when these conditions hold, it is

D(R,S)=max⁡ ⁣{0,Re⁡tr⁡(f(R)−Rlog⁡S)}ln⁡2,f(x)=xln⁡x,D(R,S)= \frac{\max\!\left\{0,\operatorname{Re}\operatorname{tr} \bigl(f(R)-R\log S\bigr)\right\}}{\ln2}, \qquad f(x)=x\ln x,D(R,S)=ln2max{0,Retr(f(R)−RlogS)}​,f(x)=xlnx,

where the matrix functions use continuous functional calculus and the real logarithm is total, with ln⁡0=0\ln0=0ln0=0. The inner supremum and outer infimum are taken in the extended reals. No channel assumption is imposed on NNN; no nonemptiness assumption is imposed on AAA or BBB. Thus an empty outer feasible set gives +∞+\infty+∞, while an empty set of unit vectors gives inner supremum −∞-\infty−∞.

QAlgorithms.Papers.Independent.Intakethird100.P65.dampingPrimalWitness

For every real number rrr, this definition returns the complex diagonal matrix on the Boolean index set, ordered as (false,true)(\mathrm{false},\mathrm{true})(false,true),

(1+1−r001−r+1−r).\begin{pmatrix} 1+\sqrt{1-r}&0\\ 0&\sqrt{1-r}+1-r \end{pmatrix}.(1+1−r​0​01−r​+1−r​).

The entries are real numbers embedded into C\mathbb CC. No restriction on rrr is imposed; the real square root used here equals 000 for negative arguments, so for r>1r>1r>1 the diagonal entries are 111 and 1−r1-r1−r.

QAlgorithms.Papers.Independent.Intakethird100.P65.deFinetti

For a finite type AAA with decidable equality and a natural number nnn, this defines a complex matrix DnD_nDn​ indexed by words x,y∈(A×A)nx,y\in(A\times A)^nx,y∈(A×A)n, with entries

Dn(x,y)=#{π∈Sym⁡(n):x=y∘π}n!(max⁡(∣A∣2+n−1,0)n).D_n(x,y)= \frac{\#\{\pi\in\operatorname{Sym}(n):x=y\circ\pi\}} {n!\binom{\max(|A|^2+n-1,0)}n}.Dn​(x,y)=n!(nmax(∣A∣2+n−1,0)​)#{π∈Sym(n):x=y∘π}​.

The numerator counts every permutation satisfying the equality, including distinct permutations that act identically on a word with repeated coordinates. The subtraction in the binomial coefficient is natural-number subtraction, truncated at zero. No nonemptiness assumption on AAA or positivity assumption on nnn is imposed. For n=0n=0n=0, the matrix has one entry, equal to one; for empty AAA and n>0n>0n>0, its index set is empty. Any division by zero in the defining formula is interpreted as zero.

QAlgorithms.Papers.Independent.Intakethird100.P65.deFinettiK1

For finite types A,BA,BA,B with decidable equality, an arbitrary function NNN from complex matrices on AAA to complex matrices on BBB, a natural number nnn, and a real number eee, this defines a set of matrices on Rn×BnR_n\times B^nRn​×Bn, where Rn=(A×A)n×AnR_n=(A\times A)^n\times A^nRn​=(A×A)n×An. Define

Dn(z,z′)=#{π∈Sym⁡(n):z=z′∘π}n!(max⁡(∣A∣2+n−1,0)n),vn((z,a),b)=(Dn)(z, i↦(ai,bi)),Wn=∣vn⟩⟨vn∣.D_n(z,z')= \frac{\#\{\pi\in\operatorname{Sym}(n):z=z'\circ\pi\}} {n!\binom{\max(|A|^2+n-1,0)}n}, \quad v_n((z,a),b)=(\sqrt{D_n})\bigl(z,\ i\mapsto(a_i,b_i)\bigr), \quad W_n=|v_n\rangle\langle v_n|.Dn​(z,z′)=n!(nmax(∣A∣2+n−1,0)​)#{π∈Sym(n):z=z′∘π}​,vn​((z,a),b)=(Dn​​)(z, i↦(ai​,bi​)),Wn​=∣vn​⟩⟨vn​∣.

Here matrix square roots use continuous functional calculus. Define the matrix map

Nn(Y)(u,v)=∑x,x′∈AnY(x,x′)∏i=0n−1N(Exi,xi′)(ui,vi),N_n(Y)(u,v)= \sum_{x,x'\in A^n}Y(x,x') \prod_{i=0}^{n-1}N(E_{x_i,x'_i})(u_i,v_i),Nn​(Y)(u,v)=x,x′∈An∑​Y(x,x′)i=0∏n−1​N(Exi​,xi′​​)(ui​,vi​),

where Ea,a′E_{a,a'}Ea,a′​ is the matrix with its sole nonzero entry, equal to one, at (a,a′)(a,a')(a,a′), and let Ωn=(id⁡Rn⊗Nn)(Wn)\Omega_n=(\operatorname{id}_{R_n}\otimes N_n)(W_n)Ωn​=(idRn​​⊗Nn​)(Wn​), where the amplification applies NnN_nNn​ separately to every block indexed by a pair of reference indices. The defined set consists precisely of matrices σ\sigmaσ satisfying

σ⪰0,Tr⁡σ=1,P(σ,Ωn)≤e,∑b∈Bnσ((r,b),(r′,b))=∑a∈AnWn((r,a),(r′,a))for all r,r′,\sigma\succeq0,\qquad \operatorname{Tr}\sigma=1,\qquad P(\sigma,\Omega_n)\le e,\qquad \sum_{b\in B^n}\sigma((r,b),(r',b)) =\sum_{a\in A^n}W_n((r,a),(r',a)) \quad\text{for all }r,r',σ⪰0,Trσ=1,P(σ,Ωn​)≤e,b∈Bn∑​σ((r,b),(r′,b))=a∈An∑​Wn​((r,a),(r′,a))for all r,r′,

where

F0(X,Y)=ReTr⁡X YX,P(X,Y)=1−(F0(X,Y)+(1−ReTr⁡X)(1−ReTr⁡Y))2.F_0(X,Y)=\operatorname{ReTr}\sqrt{\sqrt X\,Y\sqrt X}, \qquad P(X,Y)= \sqrt{1-\left(F_0(X,Y)+ \sqrt{(1-\operatorname{ReTr}X)(1-\operatorname{ReTr}Y)}\right)^2}.F0​(X,Y)=ReTrX​YX​​,P(X,Y)=1−(F0​(X,Y)+(1−ReTrX)(1−ReTrY)​)2​.

No channel assumption on NNN is made. Real square roots and reciprocals are total, with square roots of negative numbers and reciprocals of zero equal to zero. Empty register types and n=0n=0n=0 are included; at n=0n=0n=0, the set is the singleton scalar density matrix 111 if e≥0e\ge0e≥0, and is empty otherwise.

QAlgorithms.Papers.Independent.Intakethird100.P65.deFinettiK2

For finite types A,BA,BA,B with decidable equality, an arbitrary matrix-valued function NNN from matrices on AAA to matrices on BBB, a natural number nnn, and a real number eee, this defines a set of matrices on Rn×BnR_n\times B^nRn​×Bn, with Rn=(A×A)n×AnR_n=(A\times A)^n\times A^nRn​=(A×A)n×An. Set

Dn(z,z′)=#{π∈Sym⁡(n):z=z′∘π}n!(max⁡(∣A∣2+n−1,0)n),vn((z,a),b)=(Dn)(z, i↦(ai,bi)),Wn=∣vn⟩⟨vn∣,D_n(z,z')= \frac{\#\{\pi\in\operatorname{Sym}(n):z=z'\circ\pi\}} {n!\binom{\max(|A|^2+n-1,0)}n}, \quad v_n((z,a),b)=(\sqrt{D_n})\bigl(z,\ i\mapsto(a_i,b_i)\bigr), \quad W_n=|v_n\rangle\langle v_n|,Dn​(z,z′)=n!(nmax(∣A∣2+n−1,0)​)#{π∈Sym(n):z=z′∘π}​,vn​((z,a),b)=(Dn​​)(z, i↦(ai​,bi​)),Wn​=∣vn​⟩⟨vn​∣,

using continuous functional calculus for matrix square roots, and define

Nn(Y)(u,v)=∑x,x′∈AnY(x,x′)∏i=0n−1N(Exi,xi′)(ui,vi),Ωn=(id⁡Rn⊗Nn)(Wn).N_n(Y)(u,v)= \sum_{x,x'\in A^n}Y(x,x') \prod_{i=0}^{n-1}N(E_{x_i,x'_i})(u_i,v_i), \qquad \Omega_n=(\operatorname{id}_{R_n}\otimes N_n)(W_n).Nn​(Y)(u,v)=x,x′∈An∑​Y(x,x′)i=0∏n−1​N(Exi​,xi′​​)(ui​,vi​),Ωn​=(idRn​​⊗Nn​)(Wn​).

The set consists of matrices σ\sigmaσ for which there exists a map TTT on matrices from AnA^nAn to BnB^nBn having a finite Kraus representation T(X)=∑jKjXKj†T(X)=\sum_jK_jXK_j^\daggerT(X)=∑j​Kj​XKj†​ with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I, satisfying

T(PπX)=Pπ(T(X))for every π∈Sym⁡(n) and every X,σ=(id⁡Rn⊗T)(Wn),P(σ,Ωn)≤e,T(P_\pi X)=P_\pi(T(X)) \quad\text{for every }\pi\in\operatorname{Sym}(n)\text{ and every }X, \qquad \sigma=(\operatorname{id}_{R_n}\otimes T)(W_n), \qquad P(\sigma,\Omega_n)\le e,T(Pπ​X)=Pπ​(T(X))for every π∈Sym(n) and every X,σ=(idRn​​⊗T)(Wn​),P(σ,Ωn​)≤e,

where PπX(x,y)=X(x∘π,y∘π)P_\pi X(x,y)=X(x\circ\pi,y\circ\pi)Pπ​X(x,y)=X(x∘π,y∘π), amplifications act blockwise on the nonreference register, and

P(X,Y)=1−(ReTr⁡X YX+(1−ReTr⁡X)(1−ReTr⁡Y))2.P(X,Y)= \sqrt{1-\left(\operatorname{ReTr}\sqrt{\sqrt X\,Y\sqrt X} +\sqrt{(1-\operatorname{ReTr}X)(1-\operatorname{ReTr}Y)}\right)^2}.P(X,Y)=1−(ReTrX​YX​​+(1−ReTrX)(1−ReTrY)​)2​.

There is no channel assumption on NNN, and the definition does not separately require that σ\sigmaσ be a density matrix or separately impose a marginal equality. All natural subtractions are truncated, reciprocals of zero are zero, and real square roots of negative numbers are zero. Empty registers and n=0n=0n=0 are included; at n=0n=0n=0, the set is {1}\{1\}{1} for e≥0e\ge0e≥0 and empty otherwise.

QAlgorithms.Papers.Independent.Intakethird100.P65.deFinettiPurification

For a finite type AAA with decidable equality and a natural number nnn, let DnD_nDn​ be the matrix on (A×A)n(A\times A)^n(A×A)n given by

Dn(z,z′)=#{π∈Sym⁡(n):z=z′∘π}n!(max⁡(∣A∣2+n−1,0)n),D_n(z,z')= \frac{\#\{\pi\in\operatorname{Sym}(n):z=z'\circ\pi\}} {n!\binom{\max(|A|^2+n-1,0)}n},Dn​(z,z′)=n!(nmax(∣A∣2+n−1,0)​)#{π∈Sym(n):z=z′∘π}​,

and let Sn=DnS_n=\sqrt{D_n}Sn​=Dn​​, with the square root defined by continuous functional calculus. This definition forms the vector

vn((z,a),b)=Sn(z, i↦(ai,bi)),z∈(A×A)n,a,b∈An,v_n((z,a),b)=S_n\bigl(z,\ i\mapsto(a_i,b_i)\bigr), \qquad z\in(A\times A)^n,\quad a,b\in A^n,vn​((z,a),b)=Sn​(z, i↦(ai​,bi​)),z∈(A×A)n,a,b∈An,

and returns the matrix WnW_nWn​ with entries

Wn(x,y)=vn(x)vn(y)‾.W_n(x,y)=v_n(x)\overline{v_n(y)}.Wn​(x,y)=vn​(x)vn​(y)​.

Its index set is ((A×A)n×An)×An\bigl((A\times A)^n\times A^n\bigr)\times A^n((A×A)n×An)×An. Empty AAA and n=0n=0n=0 are included; at n=0n=0n=0, this is the one-entry matrix 111. The formula itself imposes no separate normalization condition on vnv_nvn​.

QAlgorithms.Papers.Independent.Intakethird100.P65.depolarizingFeasible

For natural numbers d,nd,nd,n, real numbers p,e,sp,e,sp,e,s, and real arrays y=(y0,…,yn)y=(y_0,\ldots,y_n)y=(y0​,…,yn​) and r=(r0,…,rn)r=(r_0,\ldots,r_n)r=(r0​,…,rn​), this predicate means that, on defining

wk=(nk)d−k(d−1d)n−k,ak=(d(1−p)+pd)k(pd)n−k(0≤k≤n),w_k=\binom nk d^{-k}\left(d-\frac1d\right)^{n-k}, \qquad a_k=\left(d(1-p)+\frac pd\right)^k\left(\frac pd\right)^{n-k} \quad(0\le k\le n),wk​=(kn​)d−k(d−d1​)n−k,ak​=(d(1−p)+dp​)k(dp​)n−k(0≤k≤n),

all the following hold:

yk−rk+ak≥0,yk≥0,0≤rk≤sfor every k,∑k=0nwkrk=1,∑k=0nwkyk≤e.y_k-r_k+a_k\ge0,\qquad y_k\ge0,\qquad 0\le r_k\le s \quad\text{for every }k, \qquad \sum_{k=0}^{n}w_kr_k=1, \qquad \sum_{k=0}^{n}w_ky_k\le e.yk​−rk​+ak​≥0,yk​≥0,0≤rk​≤sfor every k,k=0∑n​wk​rk​=1,k=0∑n​wk​yk​≤e.

There are no additional bounds on p,e,sp,e,sp,e,s or positivity assumption on ddd; reciprocals of zero are zero and 00=10^0=100=1. For d=0,n>0d=0,n>0d=0,n>0, every wkw_kwk​ is zero, so the normalization condition is impossible. For n=0n=0n=0, the predicate reduces to r0=1r_0=1r0​=1, s≥1s\ge1s≥1, and 0≤y0≤e0\le y_0\le e0≤y0​≤e, independently of d,pd,pd,p.

QAlgorithms.Papers.Independent.Intakethird100.P65.depolarizingPrimalWitness

For every finite type AAA with decidable equality and every real number ppp, writing d=∣A∣d=|A|d=∣A∣, this definition returns the scalar matrix

(d(1−p)+pd)IA,\left(d(1-p)+\frac pd\right)I_A,(d(1−p)+dp​)IA​,

with the real scalar embedded into C\mathbb CC. There is no restriction on ppp and no nonemptiness assumption on AAA. If AAA is empty, real division by 000 is defined to be 000, and the result is the unique matrix on the empty index type.

QAlgorithms.Papers.Independent.Intakethird100.P65.depolarizingSectorSum

For every finite type AAA with decidable equality and all natural numbers n,kn,kn,k, this defines a complex square matrix on An×AnA^n\times A^nAn×An as follows. Let PPP be the matrix on A×AA\times AA×A given by

P(a,b),(a′,b′)={∣A∣−1,a=b and a′=b′,0,otherwise.P_{(a,b),(a',b')}= \begin{cases} |A|^{-1},&a=b\text{ and }a'=b',\\ 0,&\text{otherwise}. \end{cases}P(a,b),(a′,b′)​={∣A∣−1,0,​a=b and a′=b′,otherwise.​

For x=(x1,x2)x=(x_1,x_2)x=(x1​,x2​) and y=(y1,y2)y=(y_1,y_2)y=(y1​,y2​) in An×AnA^n\times A^nAn×An, the defined entry is

∑S⊆{0,…,n−1}∣S∣=k∏i=0n−1{P(x1(i),x2(i)),(y1(i),y2(i)),i∈S,(IA×A−P)(x1(i),x2(i)),(y1(i),y2(i)),i∉S.\sum_{\substack{S\subseteq\{0,\ldots,n-1\}\\|S|=k}} \prod_{i=0}^{n-1} \begin{cases} P_{(x_1(i),x_2(i)),(y_1(i),y_2(i))},&i\in S,\\ (I_{A\times A}-P)_{(x_1(i),x_2(i)),(y_1(i),y_2(i))},&i\notin S. \end{cases}S⊆{0,…,n−1}∣S∣=k​∑​i=0∏n−1​{P(x1​(i),x2​(i)),(y1​(i),y2​(i))​,(IA×A​−P)(x1​(i),x2​(i)),(y1​(i),y2​(i))​,​i∈S,i∈/S.​

The sum is zero when k>nk>nk>n. When n=0,k=0n=0,k=0n=0,k=0, the empty product gives the one-entry matrix with entry one. Empty AAA is allowed, with ∣A∣−1=0|A|^{-1}=0∣A∣−1=0 when ∣A∣=0|A|=0∣A∣=0.

QAlgorithms.Papers.Independent.Intakethird100.P65.depolarizingValue

For natural numbers d,nd,nd,n and real numbers p,ep,ep,e, this defines the real infimum

inf⁡{dns: s∈R, y,r∈Rn+1, yk−rk+ak≥0, yk≥0, 0≤rk≤s(0≤k≤n),∑k=0nwkrk=1,∑k=0nwkyk≤e},\inf\left\{d^ns:\ s\in\mathbb R,\ y,r\in\mathbb R^{n+1},\ \begin{array}{l} y_k-r_k+a_k\ge0,\ y_k\ge0,\ 0\le r_k\le s\quad(0\le k\le n),\\ \sum_{k=0}^{n}w_kr_k=1,\quad \sum_{k=0}^{n}w_ky_k\le e \end{array}\right\},inf{dns: s∈R, y,r∈Rn+1, yk​−rk​+ak​≥0, yk​≥0, 0≤rk​≤s(0≤k≤n),∑k=0n​wk​rk​=1,∑k=0n​wk​yk​≤e​},

where

wk=(nk)(1d)k(d−1d)n−k,ak=(d(1−p)+pd)k(pd)n−k.w_k=\binom nk\left(\frac1d\right)^k\left(d-\frac1d\right)^{n-k}, \qquad a_k=\left(d(1-p)+\frac pd\right)^k\left(\frac pd\right)^{n-k}.wk​=(kn​)(d1​)k(d−d1​)n−k,ak​=(d(1−p)+dp​)k(dp​)n−k.

Division by zero is zero and 00=10^0=100=1; no restriction on p,ep,ep,e or positivity assumption on ddd is imposed. This is a real-valued infimum, whose value is defined to be zero if the set is empty or unbounded below. In particular, feasibility is impossible for d=0,n>0d=0,n>0d=0,n>0, giving value zero; for n=0n=0n=0, the value is one when e≥0e\ge0e≥0 and zero when e<0e<0e<0.

QAlgorithms.Papers.Independent.Intakethird100.P65.fixedSmoothInfo

For finite types A,BA,BA,B with decidable equality, any complex matrix ρ\rhoρ indexed by A×BA\times BA×B, and any real eee, this is the extended-real infimum

inf⁡{ℓ∈R | ∃τ,σ, τ⪰0, Re⁡Tr⁡τ≤1, P(ρ,τ)≤e, τA=ρA, σ⪰0, Tr⁡σ=1, 2ℓ(τA⊗σ)−τ⪰0},\inf\left\{\ell\in\mathbb R\ \middle|\ \exists\tau,\sigma,\ \tau\succeq0,\ \operatorname{Re}\operatorname{Tr}\tau\le1,\ P(\rho,\tau)\le e,\ \tau_A=\rho_A,\ \sigma\succeq0,\ \operatorname{Tr}\sigma=1,\ 2^\ell(\tau_A\otimes\sigma)-\tau\succeq0 \right\},inf{ℓ∈R ​ ∃τ,σ, τ⪰0, ReTrτ≤1, P(ρ,τ)≤e, τA​=ρA​, σ⪰0, Trσ=1, 2ℓ(τA​⊗σ)−τ⪰0},

where τ\tauτ is indexed by A×BA\times BA×B, σ\sigmaσ by BBB, and XA(a,a′)=∑bX((a,b),(a′,b))X_A(a,a')=\sum_bX((a,b),(a',b))XA​(a,a′)=∑b​X((a,b),(a′,b)). Here, with S(X)S(X)S(X) denoting continuous functional calculus with the real square-root function,

P(X,Y)=1−(Re⁡Tr⁡S(S(X)YS(X))+(1−Re⁡Tr⁡X)(1−Re⁡Tr⁡Y))2.P(X,Y)=\sqrt{1-\left(\operatorname{Re}\operatorname{Tr}S(S(X)YS(X)) +\sqrt{(1-\operatorname{Re}\operatorname{Tr}X)(1-\operatorname{Re}\operatorname{Tr}Y)}\right)^2}.P(X,Y)=1−(ReTrS(S(X)YS(X))+(1−ReTrX)(1−ReTrY)​)2​.

No state assumption is imposed on ρ\rhoρ; the displayed equality requires the entire AAA-marginal to be preserved. Zero τ\tauτ is permitted when it satisfies these conditions. Real square roots of negative arguments are zero; an empty feasible set gives +∞+\infty+∞, including for e<0e<0e<0, and an unbounded-below feasible set gives −∞-\infty−∞.

QAlgorithms.Papers.Independent.Intakethird100.P65.flaggedErasureChannel

For every finite type AAA and every real ppp, this defines a map from complex square matrices on AAA to complex square matrices on A⊔{⊥}A\sqcup\{\bot\}A⊔{⊥}, where ⊥\bot⊥ is a new distinguished index, by

Φp(X)a,b=(1−p)Xa,b(a,b∈A),Φp(X)⊥,⊥=p Tr⁡X,Φp(X)a,⊥=Φp(X)⊥,a=0.\Phi_p(X)_{a,b}=(1-p)X_{a,b}\quad(a,b\in A),\qquad \Phi_p(X)_{\bot,\bot}=p\,\operatorname{Tr}X,\qquad \Phi_p(X)_{a,\bot}=\Phi_p(X)_{\bot,a}=0.Φp​(X)a,b​=(1−p)Xa,b​(a,b∈A),Φp​(X)⊥,⊥​=pTrX,Φp​(X)a,⊥​=Φp​(X)⊥,a​=0.

No positivity or normalization is required of XXX, and ppp may be any real number. If AAA is empty, the output is the zero matrix on the single index ⊥\bot⊥.

QAlgorithms.Papers.Independent.Intakethird100.P65.generalErrorFeasible

For finite index types A,B,C,DA,B,C,DA,B,C,D, with decidable equality, a Boolean flag ppt\mathrm{ppt}ppt, arbitrary maps N:Mat⁡(B)→Mat⁡(C)N:\operatorname{Mat}(B)\to\operatorname{Mat}(C)N:Mat(B)→Mat(C) and M:Mat⁡(A)→Mat⁡(D)M:\operatorname{Mat}(A)\to\operatorname{Mat}(D)M:Mat(A)→Mat(D), and any real ggg, this predicate asserts the existence of complex square matrices YYY on A×DA\times DA×D and ZZZ on (A×C)×(B×D)(A\times C)\times(B\times D)(A×C)×(B×D) satisfying the following conditions. Write EijE_{ij}Eij​ for the matrix unit, IAI_AIA​ for the identity, and X⪰0X\succeq0X⪰0 for positive semidefiniteness. The matrix ZZZ is positive semidefinite and, for all appropriate indices, obeys

∑b,dZ((a,c),(b,d)),((a′,c′),(b,d))=δaa′δcc′,∑bZ((a,c),(b,d)),((a′,c′),(b,d′))=δaa′∣A∣∑r∈A,bZ((r,c),(b,d)),((r,c′),(b,d′)),∑dZ((a,c),(b,d)),((a′,c′),(b′,d))=δcc′∣C∣∑s∈C,dZ((a,s),(b,d)),((a′,s),(b′,d)).\begin{aligned} \sum_{b,d}Z_{((a,c),(b,d)),((a',c'),(b,d))} &=\delta_{aa'}\delta_{cc'},\\ \sum_b Z_{((a,c),(b,d)),((a',c'),(b,d'))} &=\frac{\delta_{aa'}}{|A|} \sum_{r\in A,b}Z_{((r,c),(b,d)),((r,c'),(b,d'))},\\ \sum_d Z_{((a,c),(b,d)),((a',c'),(b',d))} &=\frac{\delta_{cc'}}{|C|} \sum_{s\in C,d}Z_{((a,s),(b,d)),((a',s),(b',d))}. \end{aligned}b,d∑​Z((a,c),(b,d)),((a′,c′),(b,d))​b∑​Z((a,c),(b,d)),((a′,c′),(b,d′))​d∑​Z((a,c),(b,d)),((a′,c′),(b′,d))​​=δaa′​δcc′​,=∣A∣δaa′​​r∈A,b∑​Z((r,c),(b,d)),((r,c′),(b,d′))​,=∣C∣δcc′​​s∈C,d∑​Z((a,s),(b,d)),((a′,s),(b′,d))​.​

Define matrices on A×DA\times DA×D by

H(a,d),(a′,d′)=∑b,b′,c,c′N(Ebb′)c,c′ Z((a,c),(b,d)),((a′,c′),(b′,d′)),JM(a,d),(a′,d′)=M(Eaa′)d,d′.H_{(a,d),(a',d')} =\sum_{b,b',c,c'}N(E_{bb'})_{c,c'}\, Z_{((a,c),(b,d)),((a',c'),(b',d'))}, \qquad J_M{}_{(a,d),(a',d')}=M(E_{aa'})_{d,d'}.H(a,d),(a′,d′)​=b,b′,c,c′∑​N(Ebb′​)c,c′​Z((a,c),(b,d)),((a′,c′),(b′,d′))​,JM​(a,d),(a′,d′)​=M(Eaa′​)d,d′​.

The remaining conditions are

Y⪰0,gIA−Tr⁡DY⪰0,Y−H+JM⪰0,Y\succeq0,\qquad gI_A-\operatorname{Tr}_D Y\succeq0,\qquad Y-H+J_M\succeq0,Y⪰0,gIA​−TrD​Y⪰0,Y−H+JM​⪰0,

where (Tr⁡DY)a,a′=∑dY(a,d),(a′,d)(\operatorname{Tr}_D Y)_{a,a'}=\sum_dY_{(a,d),(a',d)}(TrD​Y)a,a′​=∑d​Y(a,d),(a′,d)​; additionally, if ppt=true\mathrm{ppt}=\mathrm{true}ppt=true, the matrix with entries Z((a,c′),(b,d′)),((a′,c),(b′,d))Z_{((a,c'),(b,d')),((a',c),(b',d))}Z((a,c′),(b,d′)),((a′,c),(b′,d))​ is positive semidefinite. When the flag is false this additional requirement is absent. Neither map is assumed to be a channel or linear; ggg is not assumed nonnegative, and empty index types are included, with empty sums equal to zero and division by zero interpreted as zero.

QAlgorithms.Papers.Independent.Intakethird100.P65.halfDiamondDistance

For finite types A,BA,BA,B with decidable equality and arbitrary functions N,M:CA×A→CB×BN,M:\mathbb C^{A\times A}\to\mathbb C^{B\times B}N,M:CA×A→CB×B, this defines

12sup⁡{∥ΔX∥1 | X∈C(A×A)×(A×A), ∥X∥1≤1},\frac12\sup\left\{\|\Delta_X\|_1\ \middle|\ X\in\mathbb C^{(A\times A)\times(A\times A)},\ \|X\|_1\le1 \right\},21​sup{∥ΔX​∥1​ ​ X∈C(A×A)×(A×A), ∥X∥1​≤1},

where Xi,ja,a′=X((a,i),(a′,j))X^{a,a'}_{i,j}=X((a,i),(a',j))Xi,ja,a′​=X((a,i),(a′,j)) and

ΔX((a,b),(a′,b′))=N(Xa,a′)b,b′−M(Xa,a′)b,b′.\Delta_X((a,b),(a',b'))= N(X^{a,a'})_{b,b'}-M(X^{a,a'})_{b,b'}.ΔX​((a,b),(a′,b′))=N(Xa,a′)b,b′​−M(Xa,a′)b,b′​.

Here ∥Y∥1=Re⁡Tr⁡(cfc⁡(⋅)(Y†Y))\|Y\|_1=\operatorname{Re}\operatorname{Tr}\bigl(\operatorname{cfc}(\sqrt{\cdot})(Y^\dagger Y)\bigr)∥Y∥1​=ReTr(cfc(⋅​)(Y†Y)). The test matrices XXX need not be positive or normalized. Neither map is assumed linear or a channel. The supremum is the real-number supremum, which takes the default value zero if the set is unbounded above.

QAlgorithms.Papers.Independent.Intakethird100.P65.infoOnSet

For finite types A,BA,BA,B with decidable equality and an arbitrary set SSS of complex matrices on A×BA\times BA×B, this defines the extended-real infimum

inf⁡σ∈S inf⁡ω⪰0Tr⁡ω=1Dmax⁡ ⁣(σ ∥ σA⊗ω),σA(a,a′)=∑b∈Bσ((a,b),(a′,b)),\inf_{\sigma\in S}\ \inf_{\substack{\omega\succeq0\\ \operatorname{Tr}\omega=1}} D_{\max}\!\left(\sigma\ \middle\|\ \sigma_A\otimes\omega\right), \qquad \sigma_A(a,a')=\sum_{b\in B}\sigma((a,b),(a',b)),σ∈Sinf​ ω⪰0Trω=1​inf​Dmax​(σ ∥ σA​⊗ω),σA​(a,a′)=b∈B∑​σ((a,b),(a′,b)),

where ω\omegaω is a matrix on BBB and

Dmax⁡(X∥Y)=inf⁡{λ∈R:2λY−X⪰0},D_{\max}(X\|Y)= \inf\{\lambda\in\mathbb R:2^\lambda Y-X\succeq0\},Dmax​(X∥Y)=inf{λ∈R:2λY−X⪰0},

with this last infimum taken in the extended reals by embedding each admissible real λ\lambdaλ. The set SSS is not required to contain only density matrices or positive matrices. Empty register types and empty SSS are allowed; an infimum over an empty set is +∞+\infty+∞, while an admissible set of exponents unbounded below has infimum −∞-\infty−∞.

QAlgorithms.Papers.Independent.Intakethird100.P65.linkChoi

For finite index types A,B,C,DA,B,C,DA,B,C,D, a complex matrix JJJ indexed by B×CB\times CB×C, and a complex matrix ZZZ indexed by (A×C)×(B×D)(A\times C)\times(B\times D)(A×C)×(B×D), this defines a complex matrix indexed by A×DA\times DA×D by

[link⁡(J,Z)](i,d),(j,d′)=∑a,c∈B ∑u,v∈CJ(a,u),(c,v)Z((i,u),(a,d)), ((j,v),(c,d′))[\operatorname{link}(J,Z)]_{(i,d),(j,d')} = \sum_{a,c\in B}\ \sum_{u,v\in C} J_{(a,u),(c,v)} Z_{((i,u),(a,d)),\,((j,v),(c,d'))}[link(J,Z)](i,d),(j,d′)​=a,c∈B∑​ u,v∈C∑​J(a,u),(c,v)​Z((i,u),(a,d)),((j,v),(c,d′))​

for every i,j∈Ai,j\in Ai,j∈A and d,d′∈Dd,d'\in Dd,d′∈D. There are no positivity, normalization, or channel hypotheses on either matrix. Empty index types are included, and a sum over an empty index type is 000.

QAlgorithms.Papers.Independent.Intakethird100.P65.maxEntangledInput

For every finite type AAA equipped with decidable equality, this defines a complex matrix JJJ indexed by A×AA\times AA×A. Writing d=∣A∣d=|A|d=∣A∣, its entries are

J((a,b),(a′,b′))={d−1,a=b and a′=b′,0,otherwise.J((a,b),(a',b'))= \begin{cases} d^{-1},&a=b\ \text{and}\ a'=b',\\ 0,&\text{otherwise}. \end{cases}J((a,b),(a′,b′))={d−1,0,​a=b and a′=b′,otherwise.​

For d>0d>0d>0, this is the outer product of d−1/2∑a∈A∣a,a⟩d^{-1/2}\sum_{a\in A}|a,a\rangled−1/2∑a∈A​∣a,a⟩ with itself. The definition also includes empty AAA: then JJJ is the unique empty matrix, and the scalar inverse of zero is defined to be zero.

QAlgorithms.Papers.Independent.Intakethird100.P65.maxInfo

For finite types A,BA,BA,B with decidable equality and any complex matrix ρ\rhoρ indexed by A×BA\times BA×B, define its AAA-marginal by ρA(a,a′)=∑b∈Bρ((a,b),(a′,b))\rho_A(a,a')=\sum_{b\in B}\rho((a,b),(a',b))ρA​(a,a′)=∑b∈B​ρ((a,b),(a′,b)). This definition returns the following infimum in the extended real numbers:

inf⁡{ℓ∈R | ∃σ∈CB×B, σ⪰0, Tr⁡σ=1, 2ℓ(ρA⊗σ)−ρ⪰0}.\inf\left\{\ell\in\mathbb R\ \middle|\ \exists\sigma\in\mathbb C^{B\times B},\ \sigma\succeq0,\ \operatorname{Tr}\sigma=1,\ 2^\ell(\rho_A\otimes\sigma)-\rho\succeq0 \right\}.inf{ℓ∈R ​ ∃σ∈CB×B, σ⪰0, Trσ=1, 2ℓ(ρA​⊗σ)−ρ⪰0}.

Here (ρA⊗σ)((a,b),(a′,b′))=ρA(a,a′)σ(b,b′)(\rho_A\otimes\sigma)((a,b),(a',b'))=\rho_A(a,a')\sigma(b,b')(ρA​⊗σ)((a,b),(a′,b′))=ρA​(a,a′)σ(b,b′). No positivity, nonzero, or trace assumption is imposed on ρ\rhoρ. An empty feasible set gives +∞+\infty+∞, and a feasible set unbounded below gives −∞-\infty−∞; in particular, if BBB is empty, there is no trace-one σ\sigmaσ.

QAlgorithms.Papers.Independent.Intakethird100.P65.minimumError

For finite index types A,B,C,DA,B,C,DA,B,C,D, with decidable equality, an assistance label ooo, and arbitrary matrix maps N:B→CN:B\to CN:B→C and M:A→DM:A\to DM:A→D, this defines the real infimum

inf⁡P12sup⁡∥X∥1≤1∥(id⁡A⊗(WN(P)−M))(X)∥1,\inf_P\frac12\sup_{\|X\|_1\le1} \left\|(\operatorname{id}_A\otimes(\mathcal W_N(P)-M))(X)\right\|_1,Pinf​21​∥X∥1​≤1sup​∥(idA​⊗(WN​(P)−M))(X)∥1​,

where XXX ranges over all complex matrices on A×AA\times AA×A, not merely density matrices, ∥Z∥1=Re⁡tr⁡Z†Z\|Z\|_1=\operatorname{Re}\operatorname{tr}\sqrt{Z^\dagger Z}∥Z∥1​=RetrZ†Z​, amplification applies the map separately to each block indexed by the first factor, and

WN(P)(Y)b,d=∑a,c∈B∑u,v∈CN(∣a⟩⟨c∣)u,vP(Y⊗∣u⟩⟨v∣)(a,b),(c,d).\mathcal W_N(P)(Y)_{b,d} =\sum_{a,c\in B}\sum_{u,v\in C} N(|a\rangle\langle c|)_{u,v} P(Y\otimes|u\rangle\langle v|)_{(a,b),(c,d)}.WN​(P)(Y)b,d​=a,c∈B∑​u,v∈C∑​N(∣a⟩⟨c∣)u,v​P(Y⊗∣u⟩⟨v∣)(a,b),(c,d)​.

The infimum ranges over maps P:A×C→B×DP:A\times C\to B\times DP:A×C→B×D satisfying the rule selected by ooo: for ns\mathrm{ns}ns, PPP is a channel and its two output marginals depend only on the corresponding input marginals through channels; for nsPpt\mathrm{nsPpt}nsPpt, additionally ΓD∘P∘ΓC\Gamma_D\circ P\circ\Gamma_CΓD​∘P∘ΓC​ is a channel, with Γ\GammaΓ transposing the second tensor factor; for ent\mathrm{ent}ent, PPP is obtained by adjoining a shared positive semidefinite trace-one state of finite positive local dimensions and applying a channel separately to each party’s input and share, with no communication between those local maps. A channel here has a finite Kraus representation with ∑ℓKℓ†Kℓ=I\sum_\ell K_\ell^\dagger K_\ell=I∑ℓ​Kℓ†​Kℓ​=I. Both extrema are taken in the real numbers, and no attainment is asserted; in particular, an empty infimum is assigned the real value 000.

QAlgorithms.Papers.Independent.Intakethird100.P65.oneUseSimulationCost

For finite index types A,BA,BA,B, with decidable equality, an assistance label ooo, an arbitrary matrix map N:A→BN:A\to BN:A→B, and any real tolerance eee, this defines the real infimum of log⁡2m\log_2 mlog2​m over positive integers mmm satisfying

inf⁡P allowed12sup⁡∥X∥1≤1∥(id⁡A⊗(Lm(P)−N))(X)∥1≤e,Lm(P)(Y)b,b′=∑u,v<mP(Y⊗∣u⟩⟨v∣)(u,b),(v,b′).\inf_{P\ \mathrm{allowed}}\frac12 \sup_{\|X\|_1\le1} \left\|(\operatorname{id}_A\otimes(\mathcal L_m(P)-N))(X)\right\|_1\le e, \qquad \mathcal L_m(P)(Y)_{b,b'} =\sum_{u,v<m}P(Y\otimes|u\rangle\langle v|)_{(u,b),(v,b')}.P allowedinf​21​∥X∥1​≤1sup​∥(idA​⊗(Lm​(P)−N))(X)∥1​≤e,Lm​(P)(Y)b,b′​=u,v<m∑​P(Y⊗∣u⟩⟨v∣)(u,b),(v,b′)​.

Here P:A×Fin(m)→Fin(m)×BP:A\times\mathrm{Fin}(m)\to\mathrm{Fin}(m)\times BP:A×Fin(m)→Fin(m)×B, XXX ranges over all matrices on A×AA\times AA×A, ∥Z∥1=Re⁡tr⁡Z†Z\|Z\|_1=\operatorname{Re}\operatorname{tr}\sqrt{Z^\dagger Z}∥Z∥1​=RetrZ†Z​, and amplification applies the map blockwise. “Allowed” means: for ns\mathrm{ns}ns, PPP is a channel whose two output marginals are channels of the corresponding input marginals; for nsPpt\mathrm{nsPpt}nsPpt, also its conjugation by partial transposition on the second input and second output factors is a channel; for ent\mathrm{ent}ent, PPP results from a shared finite-dimensional density matrix and separate local channels on each input and its share. Channels are finite Kraus maps with Kraus completeness ∑ℓKℓ†Kℓ=I\sum_\ell K_\ell^\dagger K_\ell=I∑ℓ​Kℓ†​Kℓ​=I. There is no assumption that NNN is a channel or that e≥0e\ge0e≥0, and the inner infimum need not be attained. The outer infimum is real-valued and equals 000 if its defining set is empty.

QAlgorithms.Papers.Independent.Intakethird100.P65.purify

For a finite type AAA with decidable equality and any complex matrix ρ\rhoρ indexed by AAA, let S(ρ)S(\rho)S(ρ) be continuous functional calculus applied to ρ\rhoρ with the real square-root function, and form the vector vvv indexed by A×AA\times AA×A with v(a,b)=S(ρ)a,bv(a,b)=S(\rho)_{a,b}v(a,b)=S(ρ)a,b​. This definition returns the matrix

purify⁡(ρ)((a,b),(a′,b′))=S(ρ)a,b S(ρ)a′,b′‾.\operatorname{purify}(\rho)((a,b),(a',b')) =S(\rho)_{a,b}\,\overline{S(\rho)_{a',b'}}.purify(ρ)((a,b),(a′,b′))=S(ρ)a,b​S(ρ)a′,b′​​.

No positivity, normalization, or nonzero assumption is made on ρ\rhoρ; the vector is not normalized by the definition. Empty AAA is included.

QAlgorithms.Papers.Independent.Intakethird100.P65.reducedErrorFeasible

For finite index types A,BA,BA,B, with decidable equality, a Boolean flag ppt\mathrm{ppt}ppt, an arbitrary map N:Mat⁡(A)→Mat⁡(B)N:\operatorname{Mat}(A)\to\operatorname{Mat}(B)N:Mat(A)→Mat(B), a natural number mmm, and a real ggg, this predicate asserts that there exist complex square matrices Y,JY,JY,J on A×BA\times BA×B and a Hermitian matrix VVV on BBB such that, writing CN(a,b),(a′,b′)=N(Eaa′)b,b′C_N{}_{(a,b),(a',b')}=N(E_{aa'})_{b,b'}CN​(a,b),(a′,b′)​=N(Eaa′​)b,b′​,

Y⪰0,J⪰0,gIA−Tr⁡BY⪰0,Y−J+CN⪰0,Tr⁡BJ=IA,IA⊗V−J⪰0,Re⁡Tr⁡V=m2.\begin{gathered} Y\succeq0,\qquad J\succeq0,\qquad gI_A-\operatorname{Tr}_B Y\succeq0,\qquad Y-J+C_N\succeq0,\\ \operatorname{Tr}_B J=I_A,\qquad I_A\otimes V-J\succeq0,\qquad \operatorname{Re}\operatorname{Tr}V=m^2. \end{gathered}Y⪰0,J⪰0,gIA​−TrB​Y⪰0,Y−J+CN​⪰0,TrB​J=IA​,IA​⊗V−J⪰0,ReTrV=m2.​

Here Eaa′E_{aa'}Eaa′​ is a matrix unit, ⪰0\succeq0⪰0 means positive semidefinite, and (Tr⁡BX)a,a′=∑bX(a,b),(a′,b)(\operatorname{Tr}_B X)_{a,a'}=\sum_bX_{(a,b),(a',b)}(TrB​X)a,a′​=∑b​X(a,b),(a′,b)​. If ppt=true\mathrm{ppt}=\mathrm{true}ppt=true, it also requires

−IA⊗VT⪯mJΓB⪯IA⊗VT,(JΓB)(a,b),(a′,b′)=J(a,b′),(a′,b),-I_A\otimes V^{\mathsf T}\preceq mJ^{\Gamma_B} \preceq I_A\otimes V^{\mathsf T}, \qquad (J^{\Gamma_B})_{(a,b),(a',b')}=J_{(a,b'),(a',b)},−IA​⊗VT⪯mJΓB​⪯IA​⊗VT,(JΓB​)(a,b),(a′,b′)​=J(a,b′),(a′,b)​,

where X⪯YX\preceq YX⪯Y means Y−X⪰0Y-X\succeq0Y−X⪰0; for a false flag these extra inequalities impose no condition. The definition allows m=0m=0m=0, arbitrary ggg, and empty index types, and does not assume that NNN is a channel or linear.

QAlgorithms.Papers.Independent.Intakethird100.P65.referenceSmoothInfo

For finite types A,BA,BA,B with decidable equality, any complex matrix ρ\rhoρ indexed by A×BA\times BA×B, and any real eee, this is the extended-real infimum

inf⁡{ℓ∈R | ∃τ,σ, τ⪰0, Re⁡Tr⁡τ≤1, P(ρ,τ)≤e, σ⪰0, Tr⁡σ=1, 2ℓ(ρA⊗σ)−τ⪰0},\inf\left\{\ell\in\mathbb R\ \middle|\ \exists\tau,\sigma,\ \tau\succeq0,\ \operatorname{Re}\operatorname{Tr}\tau\le1,\ P(\rho,\tau)\le e,\ \sigma\succeq0,\ \operatorname{Tr}\sigma=1,\ 2^\ell(\rho_A\otimes\sigma)-\tau\succeq0 \right\},inf{ℓ∈R ​ ∃τ,σ, τ⪰0, ReTrτ≤1, P(ρ,τ)≤e, σ⪰0, Trσ=1, 2ℓ(ρA​⊗σ)−τ⪰0},

where τ\tauτ is indexed by A×BA\times BA×B, σ\sigmaσ by BBB, and ρA(a,a′)=∑bρ((a,b),(a′,b))\rho_A(a,a')=\sum_b\rho((a,b),(a',b))ρA​(a,a′)=∑b​ρ((a,b),(a′,b)). With S(X)S(X)S(X) denoting continuous functional calculus with the real square-root function, the distance expression is

P(X,Y)=1−(Re⁡Tr⁡S(S(X)YS(X))+(1−Re⁡Tr⁡X)(1−Re⁡Tr⁡Y))2.P(X,Y)=\sqrt{1-\left(\operatorname{Re}\operatorname{Tr}S(S(X)YS(X)) +\sqrt{(1-\operatorname{Re}\operatorname{Tr}X)(1-\operatorname{Re}\operatorname{Tr}Y)}\right)^2}.P(X,Y)=1−(ReTrS(S(X)YS(X))+(1−ReTrX)(1−ReTrY)​)2​.

There is no requirement that τA=ρA\tau_A=\rho_AτA​=ρA​, and no positivity or normalization assumption on ρ\rhoρ. Zero τ\tauτ is allowed. Real square roots of negative arguments are zero; empty feasibility gives +∞+\infty+∞, including for e<0e<0e<0, while unboundedness below gives −∞-\infty−∞.

QAlgorithms.Papers.Independent.Intakethird100.P65.roundedInfoCost

For every real xxx, this defines the real number log⁡2 ⁣(⌈2x/2⌉N)\log_2\!\left(\lceil2^{x/2}\rceil_{\mathbb N}\right)log2​(⌈2x/2⌉N​), where the ceiling is the least natural number at least its argument, subsequently regarded as a real number. Since 2x/2>02^{x/2}>02x/2>0, this ceiling is always at least 111; in particular, the value is 000 for every x≤0x\le0x≤0.

QAlgorithms.Papers.Independent.Intakethird100.P65.roundedTraceCost

For every real ttt, this defines log⁡2 ⁣(⌈t⌉N)\log_2\!\left(\lceil\sqrt t\rceil_{\mathbb N}\right)log2​(⌈t​⌉N​), with the natural-number ceiling subsequently regarded as a real number. The real square root and logarithm are total: t=0\sqrt t=0t​=0 when t≤0t\le0t≤0, and log⁡20=0\log_2 0=0log2​0=0. Consequently this definition gives 000 for every t≤1t\le1t≤1, including negative ttt and t=0t=0t=0.

QAlgorithms.Papers.Independent.Intakethird100.P65.roundingCorrection

For every real xxx, this defines log⁡2 ⁣(⌈2x/2⌉N)−x/2\log_2\!\left(\lceil2^{x/2}\rceil_{\mathbb N}\right)-x/2log2​(⌈2x/2⌉N​)−x/2, where the ceiling is taken in the natural numbers and then converted to a real number. All real xxx are included; for x≤0x\le0x≤0, the ceiling is 111 and this value is −x/2-x/2−x/2.

QAlgorithms.Papers.Independent.Intakethird100.P65.sectorEigenvalue

For arbitrary natural numbers d,n,kd,n,kd,n,k and an arbitrary real number ppp, this defines

(d(1−p)+pd)k(pd)max⁡(n−k,0).\left(d(1-p)+\frac pd\right)^k \left(\frac pd\right)^{\max(n-k,0)}.(d(1−p)+dp​)k(dp​)max(n−k,0).

There is no assumption that d>0d>0d>0, k≤nk\le nk≤n, or p∈[0,1]p\in[0,1]p∈[0,1]. The exponent subtraction is truncated at zero, division by zero is zero, and 00=10^0=100=1. In particular, for k>nk>nk>n the second factor is one, and for n=k=0n=k=0n=k=0 the entire expression is one.

QAlgorithms.Papers.Independent.Intakethird100.P65.sectorMatrix

For a finite type AAA with decidable equality, a natural number nnn, and arbitrary real coefficients r0,…,rnr_0,\ldots,r_nr0​,…,rn​, this defines a complex matrix on An×AnA^n\times A^nAn×An by

∑k=0nrk∑S⊆{0,…,n−1}∣S∣=k⨂i=0n−1{QA,i∈S,I−QA,i∉S,\sum_{k=0}^{n}r_k \sum_{\substack{S\subseteq\{0,\ldots,n-1\}\\ |S|=k}} \bigotimes_{i=0}^{n-1} \begin{cases} Q_A,&i\in S,\\ I-Q_A,&i\notin S, \end{cases}k=0∑n​rk​S⊆{0,…,n−1}∣S∣=k​∑​i=0⨂n−1​{QA​,I−QA​,​i∈S,i∈/S,​

where AnA^nAn denotes the functions from {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} to AAA, tensor factors are grouped into the two AnA^nAn registers, and

QA((a,b),(a′,b′))={∣A∣−1,a=b and a′=b′,0,otherwise.Q_A((a,b),(a',b'))= \begin{cases} |A|^{-1},&a=b\text{ and }a'=b',\\ 0,&\text{otherwise}. \end{cases}QA​((a,b),(a′,b′))={∣A∣−1,0,​a=b and a′=b′,otherwise.​

Explicitly, an entry of each tensor product is the product of its factors’ corresponding entries. The coefficients need not be nonnegative or sum to one. The type AAA may be empty, and n=0n=0n=0 is included: the empty tensor product is the scalar matrix 111, so the result is the scalar matrix r0r_0r0​.

QAlgorithms.Papers.Independent.Intakethird100.P65.sectorWeight

For arbitrary natural numbers d,n,kd,n,kd,n,k, this defines the real number

(nk)(1d)k(d−1d)max⁡(n−k,0).\binom nk\left(\frac1d\right)^k \left(d-\frac1d\right)^{\max(n-k,0)}.(kn​)(d1​)k(d−d1​)max(n−k,0).

The subtraction in the exponent is truncated natural-number subtraction. No positivity assumption on ddd, or restriction k≤nk\le nk≤n, is imposed. Division by zero has value zero, and every zeroth power, including 000^000, is one. Thus the value is zero when k>nk>nk>n; when d=0d=0d=0 and n>0n>0n>0, it is zero for every kkk; and at n=k=0n=k=0n=k=0, it is one for every ddd.

QAlgorithms.Papers.Independent.Intakethird100.P65.smoothChannelDmax

For finite index types A,BA,BA,B with decidable equality, arbitrary matrix-valued functions N,M:A-matrices→B-matricesN,M:A\text{-matrices}\to B\text{-matrices}N,M:A-matrices→B-matrices, and any real eee, this is the infimum in [−∞,+∞][-\infty,+\infty][−∞,+∞] of the values Dmax⁡(JT,JM)D_{\max}(J_T,J_M)Dmax​(JT​,JM​) over all channels TTT satisfying Δ(T,N)≤e\Delta(T,N)\le eΔ(T,N)≤e. Here a channel TTT means that there exist a natural number rrr and matrices Kj∈CB×AK_j\in\mathbb C^{B\times A}Kj​∈CB×A, 0≤j<r0\le j<r0≤j<r, such that ∑jKj†Kj=IA\sum_jK_j^\dagger K_j=I_A∑j​Kj†​Kj​=IA​ and T(X)=∑jKjXKj†T(X)=\sum_jK_jXK_j^\daggerT(X)=∑j​Kj​XKj†​ for every XXX; JTJ_TJT​ is the matrix with entries (JT)(a,b),(a′,b′)=[T(Ea,a′)]b,b′(J_T)_{(a,b),(a',b')}=[T(E_{a,a'})]_{b,b'}(JT​)(a,b),(a′,b′)​=[T(Ea,a′​)]b,b′​, where Ea,a′E_{a,a'}Ea,a′​ has its sole nonzero entry 111 at (a,a′)(a,a')(a,a′); and

Dmax⁡(R,S)=inf⁡{ℓ∈R:2ℓS−R is positive semidefinite},Δ(T,N)=12sup⁡∥X∥1≤1∥(idA⊗(T−N))(X)∥1.D_{\max}(R,S)=\inf\{\ell\in\mathbb R:2^\ell S-R\text{ is positive semidefinite}\}, \qquad \Delta(T,N)=\frac12\sup_{\|X\|_1\le1} \|(\mathrm{id}_A\otimes(T-N))(X)\|_1.Dmax​(R,S)=inf{ℓ∈R:2ℓS−R is positive semidefinite},Δ(T,N)=21​∥X∥1​≤1sup​∥(idA​⊗(T−N))(X)∥1​.

The first infimum is taken in the extended reals. In the supremum, XXX ranges over all complex matrices indexed by A×AA\times AA×A, ∥Y∥1=Re⁡tr⁡Y†Y\|Y\|_1=\operatorname{Re}\operatorname{tr}\sqrt{Y^\dagger Y}∥Y∥1​=RetrY†Y​, and the amplification acts on each block by [(idA⊗F)(X)](a,b),(a′,b′)=[F(X(a,⋅),(a′,⋅))]b,b′[(\mathrm{id}_A\otimes F)(X)]_{(a,b),(a',b')}=[F(X_{(a,\cdot),(a',\cdot)})]_{b,b'}[(idA​⊗F)(X)](a,b),(a′,b′)​=[F(X(a,⋅),(a′,⋅)​)]b,b′​. Neither NNN nor MMM is assumed to be a channel, eee need not be nonnegative, and empty index types are allowed. An empty feasible set for either extended-real infimum gives +∞+\infty+∞; attainment is not part of this definition.

QAlgorithms.Papers.Independent.Intakethird100.P65.smoothChannelInfo

For finite types A,BA,BA,B with decidable equality, an arbitrary function N:CA×A→CB×BN:\mathbb C^{A\times A}\to\mathbb C^{B\times B}N:CA×A→CB×B, and any real eee, this defines the real infimum of I(M)I(M)I(M) over all finite Kraus channels MMM whose half-distance from NNN is at most eee. Such a channel has, for some natural number rrr, matrices Kj∈CB×AK_j\in\mathbb C^{B\times A}Kj​∈CB×A satisfying ∑j<rKj†Kj=1\sum_{j<r}K_j^\dagger K_j=1∑j<r​Kj†​Kj​=1 and M(X)=∑j<rKjXKj†M(X)=\sum_{j<r}K_jXK_j^\daggerM(X)=∑j<r​Kj​XKj†​ for every matrix XXX. Its half-distance is

12sup⁡∥X∥1≤1∥ΔX∥1,ΔX((a,b),(a′,b′))=M(Xa,a′)b,b′−N(Xa,a′)b,b′,Xi,ja,a′=X((a,i),(a′,j)),\frac12\sup_{\|X\|_1\le1}\|\Delta_X\|_1, \quad \Delta_X((a,b),(a',b'))= M(X^{a,a'})_{b,b'}-N(X^{a,a'})_{b,b'}, \quad X^{a,a'}_{i,j}=X((a,i),(a',j)),21​∥X∥1​≤1sup​∥ΔX​∥1​,ΔX​((a,b),(a′,b′))=M(Xa,a′)b,b′​−N(Xa,a′)b,b′​,Xi,ja,a′​=X((a,i),(a′,j)),

with ∥Y∥1=Re⁡Tr⁡cfc⁡(⋅)(Y†Y)\|Y\|_1=\operatorname{Re}\operatorname{Tr}\operatorname{cfc}(\sqrt{\cdot})(Y^\dagger Y)∥Y∥1​=ReTrcfc(⋅​)(Y†Y) and a real supremum. To define I(M)I(M)I(M), set R((a,b),(a′,b′))=M(∣A∣−1Ea,a′)b,b′R((a,b),(a',b'))=M(|A|^{-1}E_{a,a'})_{b,b'}R((a,b),(a′,b′))=M(∣A∣−1Ea,a′​)b,b′​, RA(a,a′)=∑bR((a,b),(a′,b))R_A(a,a')=\sum_bR((a,b),(a',b))RA​(a,a′)=∑b​R((a,b),(a′,b)), and take the extended-real infimum of all real ℓ\ellℓ for which some positive semidefinite trace-one matrix σ\sigmaσ on BBB satisfies 2ℓ(RA⊗σ)−R⪰02^\ell(R_A\otimes\sigma)-R\succeq02ℓ(RA​⊗σ)−R⪰0; I(M)I(M)I(M) is its real conversion, sending either infinity to zero. No channel assumption is made on NNN. The outer real infimum is zero for an empty feasible set or one unbounded below, and empty register types are included.

QAlgorithms.Papers.Independent.Intakethird100.P65.smoothInfo

For finite types A,BA,BA,B with decidable equality, any complex matrix ρ\rhoρ indexed by A×BA\times BA×B, and any real eee, this is the extended-real infimum

inf⁡{ℓ∈R | ∃τ,σ, τ⪰0, Re⁡Tr⁡τ≤1, P(ρ,τ)≤e, σ⪰0, Tr⁡σ=1, 2ℓ(τA⊗σ)−τ⪰0},\inf\left\{\ell\in\mathbb R\ \middle|\ \exists\tau,\sigma,\ \tau\succeq0,\ \operatorname{Re}\operatorname{Tr}\tau\le1,\ P(\rho,\tau)\le e,\ \sigma\succeq0,\ \operatorname{Tr}\sigma=1,\ 2^\ell(\tau_A\otimes\sigma)-\tau\succeq0 \right\},inf{ℓ∈R ​ ∃τ,σ, τ⪰0, ReTrτ≤1, P(ρ,τ)≤e, σ⪰0, Trσ=1, 2ℓ(τA​⊗σ)−τ⪰0},

where τ\tauτ is indexed by A×BA\times BA×B, σ\sigmaσ by BBB, and τA(a,a′)=∑bτ((a,b),(a′,b))\tau_A(a,a')=\sum_b\tau((a,b),(a',b))τA​(a,a′)=∑b​τ((a,b),(a′,b)). Explicitly, writing S(X)S(X)S(X) for continuous functional calculus with the real square-root function,

P(X,Y)=1−(Re⁡Tr⁡S(S(X)YS(X))+(1−Re⁡Tr⁡X)(1−Re⁡Tr⁡Y))2.P(X,Y)=\sqrt{1-\left(\operatorname{Re}\operatorname{Tr}S(S(X)YS(X)) +\sqrt{(1-\operatorname{Re}\operatorname{Tr}X)(1-\operatorname{Re}\operatorname{Tr}Y)}\right)^2}.P(X,Y)=1−(ReTrS(S(X)YS(X))+(1−ReTrX)(1−ReTrY)​)2​.

The matrix ρ\rhoρ need not be positive or normalized, and τ=0\tau=0τ=0 is allowed. Real square roots of negative arguments are zero. Empty feasibility gives +∞+\infty+∞, including when e<0e<0e<0, while unboundedness below gives −∞-\infty−∞.

QAlgorithms.Papers.Independent.Intakethird100.P65.smoothRobustness

For finite index types A,BA,BA,B with decidable equality, an arbitrary matrix-valued function N:A-matrices→B-matricesN:A\text{-matrices}\to B\text{-matrices}N:A-matrices→B-matrices, and any real eee, this is the real infimum of R(T)R(T)R(T) over all channels TTT satisfying Δ(T,N)≤e\Delta(T,N)\le eΔ(T,N)≤e, where

R(T)=inf⁡{t≥0:there exist a channel M and σ⪰0, tr⁡σ=1,T(X)+tM(X)1+t=tr⁡(X)σfor every X}.R(T)=\inf\left\{t\ge0: \begin{array}{l} \text{there exist a channel }M\text{ and }\sigma\succeq0,\ \operatorname{tr}\sigma=1,\\ \displaystyle\frac{T(X)+tM(X)}{1+t} =\operatorname{tr}(X)\sigma\quad\text{for every }X \end{array}\right\}.R(T)=inf⎩⎨⎧​t≥0:there exist a channel M and σ⪰0, trσ=1,1+tT(X)+tM(X)​=tr(X)σfor every X​⎭⎬⎫​.

Both infima here are real infima. A channel FFF means a function admitting a finite representation F(X)=∑j<rKjXKj†F(X)=\sum_{j<r}K_jXK_j^\daggerF(X)=∑j<r​Kj​XKj†​ with ∑j<rKj†Kj=IA\sum_{j<r}K_j^\dagger K_j=I_A∑j<r​Kj†​Kj​=IA​. The distance is

Δ(T,N)=12sup⁡X∈C(A×A)×(A×A)∥X∥1≤1∥(idA⊗(T−N))(X)∥1,∥Y∥1=Re⁡tr⁡Y†Y,\Delta(T,N)=\frac12\sup_{\substack{X\in\mathbb C^{(A\times A)\times(A\times A)}\\ \|X\|_1\le1}} \|(\mathrm{id}_A\otimes(T-N))(X)\|_1, \qquad \|Y\|_1=\operatorname{Re}\operatorname{tr}\sqrt{Y^\dagger Y},Δ(T,N)=21​X∈C(A×A)×(A×A)∥X∥1​≤1​sup​∥(idA​⊗(T−N))(X)∥1​,∥Y∥1​=RetrY†Y​,

with amplification applying the function to each block indexed by the reference coordinates in AAA. There is no channel assumption on NNN and no nonnegativity assumption on eee; empty index types are included. Either real infimum equals 000 if its feasible set is empty, and no minimizing channel is asserted to exist.

QAlgorithms.Papers.Independent.Intakethird100.P65.traceFeasible

For finite index types A,BA,BA,B, with decidable equality, an arbitrary map N:Mat⁡(A)→Mat⁡(B)N:\operatorname{Mat}(A)\to\operatorname{Mat}(B)N:Mat(A)→Mat(B), a real eee, complex square matrices Y,JY,JY,J on A×BA\times BA×B, and a complex square matrix VVV on BBB, this predicate means precisely

V=V†,Y⪰0,J⪰0,eIA−Tr⁡BY⪰0,Y−J+CN⪰0,Tr⁡BJ=IA,IA⊗V−J⪰0.V=V^\dagger,\qquad Y\succeq0,\qquad J\succeq0,\qquad eI_A-\operatorname{Tr}_B Y\succeq0,\qquad Y-J+C_N\succeq0,\qquad \operatorname{Tr}_B J=I_A,\qquad I_A\otimes V-J\succeq0.V=V†,Y⪰0,J⪰0,eIA​−TrB​Y⪰0,Y−J+CN​⪰0,TrB​J=IA​,IA​⊗V−J⪰0.

Here CN(a,b),(a′,b′)=N(Eaa′)b,b′C_N{}_{(a,b),(a',b')}=N(E_{aa'})_{b,b'}CN​(a,b),(a′,b′)​=N(Eaa′​)b,b′​, Eaa′E_{aa'}Eaa′​ is a matrix unit, (Tr⁡BX)a,a′=∑bX(a,b),(a′,b)(\operatorname{Tr}_B X)_{a,a'}=\sum_bX_{(a,b),(a',b)}(TrB​X)a,a′​=∑b​X(a,b),(a′,b)​, and ⪰0\succeq0⪰0 denotes positive semidefiniteness. There is no trace constraint on VVV, no assumption that e≥0e\ge0e≥0, and no channel or linearity assumption on NNN; empty index types are included.

QAlgorithms.Papers.Independent.Intakethird100.P65.traceValue

For finite index types A,BA,BA,B, with decidable equality, an arbitrary map N:Mat⁡(A)→Mat⁡(B)N:\operatorname{Mat}(A)\to\operatorname{Mat}(B)N:Mat(A)→Mat(B), and any real eee, this defines the real infimum of Re⁡Tr⁡V\operatorname{Re}\operatorname{Tr}VReTrV over all complex matrices Y,JY,JY,J on A×BA\times BA×B and VVV on BBB satisfying

V=V†,Y⪰0,J⪰0,eIA−Tr⁡BY⪰0,Y−J+CN⪰0,Tr⁡BJ=IA,IA⊗V−J⪰0,V=V^\dagger,\quad Y\succeq0,\quad J\succeq0,\quad eI_A-\operatorname{Tr}_B Y\succeq0,\quad Y-J+C_N\succeq0,\quad \operatorname{Tr}_B J=I_A,\quad I_A\otimes V-J\succeq0,V=V†,Y⪰0,J⪰0,eIA​−TrB​Y⪰0,Y−J+CN​⪰0,TrB​J=IA​,IA​⊗V−J⪰0,

where CN(a,b),(a′,b′)=N(Eaa′)b,b′C_N{}_{(a,b),(a',b')}=N(E_{aa'})_{b,b'}CN​(a,b),(a′,b′)​=N(Eaa′​)b,b′​, Eaa′E_{aa'}Eaa′​ is a matrix unit, and (Tr⁡BX)a,a′=∑bX(a,b),(a′,b)(\operatorname{Tr}_B X)_{a,a'}=\sum_bX_{(a,b),(a',b)}(TrB​X)a,a′​=∑b​X(a,b),(a′,b)​. Positive semidefiniteness is denoted by ⪰0\succeq0⪰0. The definition asserts neither feasibility nor attainment, allows negative eee and empty index types, and imposes no channel or linearity assumption on NNN. This is a total real-valued infimum: an empty feasible-value set, or one unbounded below, yields zero.

QAlgorithms.Papers.Independent.Intakethird100.P65.vanishingErrorQuantumCapacity

For finite index types A,BA,BA,B, with decidable equality, an assistance label ooo, and an arbitrary matrix map N:A→BN:A\to BN:A→B, this defines the extended-real supremum over real numbers r≥0r\ge0r≥0 such that, for every e>0e>0e>0 and s>0s>0s>0, all sufficiently large positive integers nnn admit a positive integer mmm satisfying (r−s)n≤log⁡2m(r-s)n\le\log_2m(r−s)n≤log2​m and an infimum error at most eee for simulating the identity on m×mm\times mm×m matrices using N[n]N^{[n]}N[n]. Here N[n](Y)y,y′=∑x,x′Yx,x′∏i<nN(∣xi⟩⟨xi′∣)yi,yi′N^{[n]}(Y)_{y,y'}=\sum_{x,x'}Y_{x,x'}\prod_{i<n}N(|x_i\rangle\langle x_i'|)_{y_i,y_i'}N[n](Y)y,y′​=∑x,x′​Yx,x′​∏i<n​N(∣xi​⟩⟨xi′​∣)yi​,yi′​​; the error is the real infimum over allowed P:Fin(m)×Bn→An×Fin(m)P:\mathrm{Fin}(m)\times B^n\to A^n\times\mathrm{Fin}(m)P:Fin(m)×Bn→An×Fin(m) of

12sup⁡∥X∥1≤1∥(id⁡m⊗(WN[n](P)−id⁡m))(X)∥1,\frac12\sup_{\|X\|_1\le1} \left\|(\operatorname{id}_m\otimes(\mathcal W_{N^{[n]}}(P)-\operatorname{id}_m))(X)\right\|_1,21​∥X∥1​≤1sup​∥(idm​⊗(WN[n]​(P)−idm​))(X)∥1​,

where WS(P)(Y)b,d=∑a,c,u,vS(∣a⟩⟨c∣)u,vP(Y⊗∣u⟩⟨v∣)(a,b),(c,d)\mathcal W_S(P)(Y)_{b,d}=\sum_{a,c,u,v}S(|a\rangle\langle c|)_{u,v}P(Y\otimes|u\rangle\langle v|)_{(a,b),(c,d)}WS​(P)(Y)b,d​=∑a,c,u,v​S(∣a⟩⟨c∣)u,v​P(Y⊗∣u⟩⟨v∣)(a,b),(c,d)​, the supremum ranges over all matrices XXX on Fin(m)×Fin(m)\mathrm{Fin}(m)\times\mathrm{Fin}(m)Fin(m)×Fin(m), amplification is blockwise, and ∥Z∥1=Re⁡tr⁡Z†Z\|Z\|_1=\operatorname{Re}\operatorname{tr}\sqrt{Z^\dagger Z}∥Z∥1​=RetrZ†Z​. Allowed maps are channels with both marginal nonsignalling conditions for ns\mathrm{ns}ns; additionally their second-factor partial-transpose conjugates are channels for nsPpt\mathrm{nsPpt}nsPpt; or they use a shared finite-dimensional density matrix and separate local channels for ent\mathrm{ent}ent. Channels have finite complete Kraus representations. The inner real infimum need not be attained and equals 000 if empty. The outer supremum is −∞-\infty−∞ if its rate set is empty and +∞+\infty+∞ if that set is unbounded above.

QAlgorithms.Papers.Independent.Intakethird100.P65.vanishingErrorSimulationCost

For finite index types A,BA,BA,B, with decidable equality, an assistance label ooo, and an arbitrary matrix map N:A→BN:A\to BN:A→B, this defines an extended-real infimum over real numbers r≥0r\ge0r≥0 with the following property: for every e>0e>0e>0 and s>0s>0s>0, all sufficiently large positive integers nnn admit a positive integer mmm such that log⁡2m≤(r+s)n\log_2m\le(r+s)nlog2​m≤(r+s)n and the infimum simulation error from the identity map on m×mm\times mm×m matrices to N[n]N^{[n]}N[n] is at most eee. Explicitly, N[n](Y)y,y′=∑x,x′Yx,x′∏i<nN(∣xi⟩⟨xi′∣)yi,yi′N^{[n]}(Y)_{y,y'}=\sum_{x,x'}Y_{x,x'}\prod_{i<n}N(|x_i\rangle\langle x_i'|)_{y_i,y_i'}N[n](Y)y,y′​=∑x,x′​Yx,x′​∏i<n​N(∣xi​⟩⟨xi′​∣)yi​,yi′​​, and that error is the real infimum over allowed P:An×Fin(m)→Fin(m)×BnP:A^n\times\mathrm{Fin}(m)\to\mathrm{Fin}(m)\times B^nP:An×Fin(m)→Fin(m)×Bn of

12sup⁡∥X∥1≤1∥(id⁡An⊗(Lm(P)−N[n]))(X)∥1,Lm(P)(Y)b,b′=∑u,v<mP(Y⊗∣u⟩⟨v∣)(u,b),(v,b′).\frac12\sup_{\|X\|_1\le1} \left\|(\operatorname{id}_{A^n}\otimes(\mathcal L_m(P)-N^{[n]}))(X)\right\|_1, \quad \mathcal L_m(P)(Y)_{b,b'}=\sum_{u,v<m}P(Y\otimes|u\rangle\langle v|)_{(u,b),(v,b')}.21​∥X∥1​≤1sup​​(idAn​⊗(Lm​(P)−N[n]))(X)​1​,Lm​(P)(Y)b,b′​=u,v<m∑​P(Y⊗∣u⟩⟨v∣)(u,b),(v,b′)​.

The supremum uses all matrices XXX on An×AnA^n\times A^nAn×An; amplification is blockwise and ∥Z∥1=Re⁡tr⁡Z†Z\|Z\|_1=\operatorname{Re}\operatorname{tr}\sqrt{Z^\dagger Z}∥Z∥1​=RetrZ†Z​. Allowed maps are channels with both marginal nonsignalling conditions for ns\mathrm{ns}ns; additionally their second-factor partial-transpose conjugates are channels for nsPpt\mathrm{nsPpt}nsPpt; or they are implemented by separate local channels with a shared finite-dimensional density matrix for ent\mathrm{ent}ent. Channels have finite complete Kraus representations. Inner error infima are real and need not be attained, with empty infimum 000; the outer infimum is in the extended reals and is +∞+\infty+∞ if no such rate exists.

QAlgorithms.Papers.Independent.Intakethird100.P65.wire

For finite index types A,B,C,DA,B,C,DA,B,C,D, with decidable equality, and arbitrary matrix maps N:B→CN:B\to CN:B→C and P:A×C→B×DP:A\times C\to B\times DP:A×C→B×D, this defines a matrix map WN(P):A→D\mathcal W_N(P):A\to DWN​(P):A→D by

WN(P)(X)b,d=∑a,c∈B∑u,v∈CN(∣a⟩⟨c∣)u,v P ⁣(X⊗∣u⟩⟨v∣)(a,b),(c,d)\mathcal W_N(P)(X)_{b,d} = \sum_{a,c\in B}\sum_{u,v\in C} N(|a\rangle\langle c|)_{u,v}\, P\!\left(X\otimes|u\rangle\langle v|\right)_{(a,b),(c,d)}WN​(P)(X)b,d​=a,c∈B∑​u,v∈C∑​N(∣a⟩⟨c∣)u,v​P(X⊗∣u⟩⟨v∣)(a,b),(c,d)​

for every complex matrix XXX on AAA and every b,d∈Db,d\in Db,d∈D. Matrix units have a single entry 111, and the tensor product has entries (X⊗Y)(i,u),(j,v)=Xi,jYu,v(X\otimes Y)_{(i,u),(j,v)}=X_{i,j}Y_{u,v}(X⊗Y)(i,u),(j,v)​=Xi,j​Yu,v​. Neither NNN nor PPP is assumed linear, positive, or trace preserving. Empty index sets are included, with empty sums equal to zero.

QAlgorithms.Papers.Independent.Intakethird100.P65.zeroErrorCost

For finite index types A,BA,BA,B, with decidable equality, and an arbitrary matrix map N:A→BN:A\to BN:A→B, this defines the real infimum of all real numbers rrr such that, for every slack s>0s>0s>0, all sufficiently large positive integers nnn satisfy Cn/n≤r+sC_n/n\le r+sCn​/n≤r+s, where

Cn=inf⁡{log⁡2m: m≥1, inf⁡P nonsignalling12sup⁡∥X∥1≤1∥(id⁡An⊗(Lm(P)−N[n]))(X)∥1≤0}.C_n=\inf\left\{\log_2m:\ m\ge1,\ \inf_{P\ \mathrm{nonsignalling}} \frac12\sup_{\|X\|_1\le1} \left\|(\operatorname{id}_{A^n}\otimes(\mathcal L_m(P)-N^{[n]}))(X)\right\|_1 \le0\right\}.Cn​=inf{log2​m: m≥1, P nonsignallinginf​21​∥X∥1​≤1sup​​(idAn​⊗(Lm​(P)−N[n]))(X)​1​≤0}.

Here N[n](Y)y,y′=∑x,x′Yx,x′∏i<nN(∣xi⟩⟨xi′∣)yi,yi′N^{[n]}(Y)_{y,y'}=\sum_{x,x'}Y_{x,x'}\prod_{i<n}N(|x_i\rangle\langle x_i'|)_{y_i,y_i'}N[n](Y)y,y′​=∑x,x′​Yx,x′​∏i<n​N(∣xi​⟩⟨xi′​∣)yi​,yi′​​, P:An×Fin(m)→Fin(m)×BnP:A^n\times\mathrm{Fin}(m)\to\mathrm{Fin}(m)\times B^nP:An×Fin(m)→Fin(m)×Bn, and Lm(P)(Y)b,b′=∑u,v<mP(Y⊗∣u⟩⟨v∣)(u,b),(v,b′)\mathcal L_m(P)(Y)_{b,b'}=\sum_{u,v<m}P(Y\otimes|u\rangle\langle v|)_{(u,b),(v,b')}Lm​(P)(Y)b,b′​=∑u,v<m​P(Y⊗∣u⟩⟨v∣)(u,b),(v,b′)​. Such a PPP must be a channel for which each output marginal is obtained from the corresponding input marginal by a channel, for every matrix input; channels are finite Kraus maps with ∑ℓKℓ†Kℓ=I\sum_\ell K_\ell^\dagger K_\ell=I∑ℓ​Kℓ†​Kℓ​=I. The supremum uses all matrices XXX on An×AnA^n\times A^nAn×An with trace norm at most 111, amplification is blockwise, and ∥Z∥1=Re⁡tr⁡Z†Z\|Z\|_1=\operatorname{Re}\operatorname{tr}\sqrt{Z^\dagger Z}∥Z∥1​=RetrZ†Z​. The zero-tolerance condition concerns an infimum and does not assert that an exactly error-free map attains it. All infima here are real-valued, with empty infimum 000; no nonnegativity condition is imposed on rrr, and division by n=0n=0n=0 is excluded by the quantifier.

QAlgorithms.Papers.Independent.Intakethird100.P65.zeroTraceValue

For finite index types A,BA,BA,B, with decidable equality, and any map N:Mat⁡(A)→Mat⁡(B)N:\operatorname{Mat}(A)\to\operatorname{Mat}(B)N:Mat(A)→Mat(B), this defines

inf⁡{Re⁡Tr⁡V:V∈Mat⁡(B),V=V†,IA⊗V−CN⪰0},CN(a,b),(a′,b′)=N(Eaa′)b,b′.\inf\left\{\operatorname{Re}\operatorname{Tr}V: V\in\operatorname{Mat}(B),\quad V=V^\dagger,\quad I_A\otimes V-C_N\succeq0\right\}, \qquad C_N{}_{(a,b),(a',b')}=N(E_{aa'})_{b,b'}.inf{ReTrV:V∈Mat(B),V=V†,IA​⊗V−CN​⪰0},CN​(a,b),(a′,b′)​=N(Eaa′​)b,b′​.

Here Eaa′E_{aa'}Eaa′​ is a matrix unit and ⪰0\succeq0⪰0 denotes positive semidefiniteness. No separate positivity condition on VVV, channel or linearity condition on NNN, or nonemptiness condition on the index types is imposed. The definition asserts neither existence nor attainment of a feasible value; its real infimum is zero for an empty set or a set unbounded below.

Human review
  • Endorsed by Shuze Chen · Nov 1, 2026

    Confirmed by the moderator at approval.

  • Endorsed by Gohan · Nov 1, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me