Definitions for Quantum Channel Simulation and the Channel's Smooth Max-Information
DefinitionQAlgorithms_Papers_Independent_Intakethird100_P65The 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.
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
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 , with decidable equality, an assistance label , and an arbitrary matrix map , this predicate is defined by the following three alternatives. For , must be a channel, and there must exist channels and such that, for every matrix , and . For , these same conditions must hold and must also be a channel, where partial transposition on the second factor means . For , there must exist positive integers , a positive semidefinite trace-one matrix on , and channels , , such that
for all . In every alternative, “channel” means a finite Kraus map whose Kraus operators satisfy .
QAlgorithms.Papers.Independent.Intakethird100.P65.BipartiteNSOperation
For finite types with decidable equality and a function from complex matrices on to complex matrices on , this predicate requires that be a finite Kraus channel and that there exist finite Kraus channels from matrices on to matrices on and from matrices on to matrices on such that, for every complex matrix on ,
A finite Kraus channel means that for some natural number and matrices from its input register to its output register, and for every . Partial traces mean summing matching indices of the discarded register; for example, . 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 , with decidable equality, an assistance label , an arbitrary matrix map , and a real number , this predicate requires and states that, for every and every slack , there exists a natural number such that every positive integer admits a positive integer satisfying
The coordinatewise map is defined by . The error is the real infimum over allowed maps of
where ranges over all matrices on , amplification is blockwise, and . Here . Allowed maps are channels with both marginal nonsignalling conditions for , such channels whose second-factor partial-transpose conjugates are also channels for , or maps implemented by separate local channels using a shared finite-dimensional density matrix for ; channels have finite complete Kraus representations. The error infimum need not be attained and equals if empty. No channel assumption is imposed on ; zero block length is excluded, and slack larger than is included.
QAlgorithms.Papers.Independent.Intakethird100.P65.Constant
For finite index types with decidable equality and an arbitrary function from complex matrices indexed by to complex matrices indexed by , this property means that there exists a positive semidefinite matrix indexed by with complex trace exactly such that, for every complex matrix indexed by , . There is no separate assumption that is a channel. Empty index types are not excluded: if is empty, no such trace-one exists.
QAlgorithms.Papers.Independent.Intakethird100.P65.EntCode
For finite index types , with decidable equality, and an arbitrary map from complex matrices on to complex matrices on , this predicate means that there exist positive natural numbers , a positive semidefinite matrix on with trace , and channels and , such that, for every complex matrix on and every pair of output indices ,
Here a channel is a map admitting a finite Kraus representation with , and denotes the matrix with its only nonzero entry equal to at . The state may be entangled or separable; the identity is required for all matrices , without positivity or normalization restrictions.
QAlgorithms.Papers.Independent.Intakethird100.P65.IteratedLimit
For a function and a real number , this property means that there exists a function such that, for every real with , the sequence converges to as , and converges to as through strictly positive real values. The sequence starts with ; no condition is imposed on , on for , or on the value .
QAlgorithms.Papers.Independent.Intakethird100.P65.NSChoi
For finite index types with decidable equality and a complex matrix indexed by , this property requires that be positive semidefinite and that all three following identities hold:
The first identity is required for every and ; the second additionally for every ; and the third additionally for every . The cardinalities and their inverses are interpreted in . No index type is assumed nonempty: empty sums are , the inverse of is , and universally quantified identities with an empty variable domain are vacuous.
QAlgorithms.Papers.Independent.Intakethird100.P65.PPTCode
For finite types with decidable equality and an arbitrary function from complex matrices on to complex matrices on , let and transpose the second register:
This predicate says that is a finite Kraus channel: there exist a natural number and matrices such that
There is no separate requirement here that itself be a channel. Empty register types and are included whenever these equations hold.
QAlgorithms.Papers.Independent.Intakethird100.P65.PermutationCovariantChannel
For finite register types with decidable equality, a natural number , and a function from complex matrices on to complex matrices on , this predicate means that has a finite Kraus representation
for some natural number and matrices , and that for every permutation of the coordinates and every complex matrix ,
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 or positivity assumption on is included; for , the covariance condition is automatic.
QAlgorithms.Papers.Independent.Intakethird100.P65.SimulationAssistance
This declares a type with exactly three distinct values, named , , and , 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 , with decidable equality, an assistance label , an arbitrary matrix map , and a real number , this predicate requires and states that, for every error tolerance and every slack , there exists a natural number such that every positive integer admits a positive integer with
Here is the map on -coordinate registers defined entrywise by
and is the real infimum, over allowed , of one half the trace-norm supremum of the blockwise amplification of on all matrices of trace norm at most , with reference register , where . Allowed maps are channels with both marginal nonsignalling conditions for ; additionally their second-factor partial-transpose conjugates are channels for ; or they arise from a shared finite-dimensional density matrix and separate local channels for . Channels admit finite complete Kraus representations, and . The error infimum need not be attained and has value if empty. The threshold may depend on ; is excluded, while is permitted.
QAlgorithms.Papers.Independent.Intakethird100.P65.StateBall
For a finite type with decidable equality, complex matrices indexed by , and any real , this predicate means that is positive semidefinite, , and
where denotes applying the real square-root function to by continuous functional calculus. There is no positivity or normalization assumption on ; is permitted. Every real square root uses the convention that the square root of a negative number is zero. In particular, a negative makes the predicate false.
QAlgorithms.Papers.Independent.Intakethird100.P65.ZeroErrorWitnessPair
For finite types with decidable equality, a map , a matrix , and a matrix , this predicate asserts the following six conditions. Define the unnormalized Choi matrix by , where has a single entry at , and define
Then the required conditions are
where . No channel assumption on , positivity assumption on , or normalization assumption on 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 for an empty set or a set unbounded below.
QAlgorithms.Papers.Independent.Intakethird100.P65.basisDepolarizing
For every finite type with decidable equality and every real , this defines a map on complex square matrices indexed by by
The input is arbitrary, and is not restricted to . The type may be empty; division by its zero cardinality then has the total-field value zero.
QAlgorithms.Papers.Independent.Intakethird100.P65.channelInfo
For finite types with decidable equality and any function , write and define a matrix indexed by through
where has entry at and zero elsewhere, and is interpreted in . Let . The returned real number is the real conversion of
This infimum is taken in the extended reals, and either infinite value converts to zero. No linearity or channel assumption is imposed on , and empty types are included.
QAlgorithms.Papers.Independent.Intakethird100.P65.channelMutualInfo
For finite index types , with decidable equality, and an arbitrary matrix map , this defines the real supremum
where , the matrix used as a purification has entries , and , meaning . The entropy appearing here is literally , and both partial traces sum over the discarded index. No channel, linearity, or positivity assumption is imposed on , so 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 ; in particular, an empty input index type admits no trace-one .
QAlgorithms.Papers.Independent.Intakethird100.P65.classicalFeasible
For an index type , a finite index type , arbitrary real functions and , and any real , this predicate means
together with
There is no separate assumption on or on the sign of . If is empty all these conditions are vacuous, including any restriction on ; if is nonempty and is empty, the required row sum makes the predicate impossible.
QAlgorithms.Papers.Independent.Intakethird100.P65.classicalMap
For index types , with finite, and any real function , this defines a map from complex square matrices on to complex square matrices on by
Thus only diagonal entries of enter the formula. No nonnegativity or sum condition is imposed on , and is arbitrary. Empty index types are allowed; an empty sum is zero.
QAlgorithms.Papers.Independent.Intakethird100.P65.classicalValue
For an index type , a finite index type , any real function , and any real , this defines the real infimum of over real functions and satisfying, for every ,
and, for every ,
No stochasticity assumption on , sign assumption on , 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 is empty, whereas they are impossible when is nonempty and is empty.
QAlgorithms.Papers.Independent.Intakethird100.P65.codeBobTranspose
For arbitrary index types and a complex square matrix indexed by , this defines a matrix on the same index set by
Thus it transposes the and coordinates together, leaving the and coordinates in place. No positivity, Hermiticity, normalization, or nonemptiness is required.
QAlgorithms.Papers.Independent.Intakethird100.P65.constantChannelRobustness
For finite index types with decidable equality and an arbitrary function from complex -indexed matrices to complex -indexed matrices, this is the real infimum of all real for which there exist a channel and a positive semidefinite -indexed matrix with such that, for every complex -indexed matrix ,
Here being a channel means that, for some natural number and matrices , and for every . No channel assumption is imposed on . The value is included, but even then a channel must exist; the denominator 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 , and it does not assert attainment.
QAlgorithms.Papers.Independent.Intakethird100.P65.constantMaxRadius
For finite index types with decidable equality, an arbitrary function , and any real , this is the extended-real infimum
A channel admits a finite representation with and . For any function , its matrix has entries , where is the matrix unit. The distance is , with an arbitrary complex matrix indexed by , , and amplification applying the function to each reference-indexed block. Neither being a channel nor is assumed, and empty index types are allowed. Each displayed infimum is interpreted in the extended reals, so an empty feasible set gives ; no attainment is asserted.
QAlgorithms.Papers.Independent.Intakethird100.P65.constantRelativeRadius
For finite index types with decidable equality and an arbitrary function , this is the extended-real infimum
Here has entries , and amplification applies a function separately to each matrix block indexed by the reference coordinates in . The quantity equals unless both are positive semidefinite with complex trace and ; when these conditions hold, it is
where the matrix functions use continuous functional calculus and the real logarithm is total, with . The inner supremum and outer infimum are taken in the extended reals. No channel assumption is imposed on ; no nonemptiness assumption is imposed on or . Thus an empty outer feasible set gives , while an empty set of unit vectors gives inner supremum .
QAlgorithms.Papers.Independent.Intakethird100.P65.dampingPrimalWitness
For every real number , this definition returns the complex diagonal matrix on the Boolean index set, ordered as ,
The entries are real numbers embedded into . No restriction on is imposed; the real square root used here equals for negative arguments, so for the diagonal entries are and .
QAlgorithms.Papers.Independent.Intakethird100.P65.deFinetti
For a finite type with decidable equality and a natural number , this defines a complex matrix indexed by words , with entries
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 or positivity assumption on is imposed. For , the matrix has one entry, equal to one; for empty and , 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 with decidable equality, an arbitrary function from complex matrices on to complex matrices on , a natural number , and a real number , this defines a set of matrices on , where . Define
Here matrix square roots use continuous functional calculus. Define the matrix map
where is the matrix with its sole nonzero entry, equal to one, at , and let , where the amplification applies separately to every block indexed by a pair of reference indices. The defined set consists precisely of matrices satisfying
where
No channel assumption on 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 are included; at , the set is the singleton scalar density matrix if , and is empty otherwise.
QAlgorithms.Papers.Independent.Intakethird100.P65.deFinettiK2
For finite types with decidable equality, an arbitrary matrix-valued function from matrices on to matrices on , a natural number , and a real number , this defines a set of matrices on , with . Set
using continuous functional calculus for matrix square roots, and define
The set consists of matrices for which there exists a map on matrices from to having a finite Kraus representation with , satisfying
where , amplifications act blockwise on the nonreference register, and
There is no channel assumption on , and the definition does not separately require that 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 are included; at , the set is for and empty otherwise.
QAlgorithms.Papers.Independent.Intakethird100.P65.deFinettiPurification
For a finite type with decidable equality and a natural number , let be the matrix on given by
and let , with the square root defined by continuous functional calculus. This definition forms the vector
and returns the matrix with entries
Its index set is . Empty and are included; at , this is the one-entry matrix . The formula itself imposes no separate normalization condition on .
QAlgorithms.Papers.Independent.Intakethird100.P65.depolarizingFeasible
For natural numbers , real numbers , and real arrays and , this predicate means that, on defining
all the following hold:
There are no additional bounds on or positivity assumption on ; reciprocals of zero are zero and . For , every is zero, so the normalization condition is impossible. For , the predicate reduces to , , and , independently of .
QAlgorithms.Papers.Independent.Intakethird100.P65.depolarizingPrimalWitness
For every finite type with decidable equality and every real number , writing , this definition returns the scalar matrix
with the real scalar embedded into . There is no restriction on and no nonemptiness assumption on . If is empty, real division by is defined to be , and the result is the unique matrix on the empty index type.
QAlgorithms.Papers.Independent.Intakethird100.P65.depolarizingSectorSum
For every finite type with decidable equality and all natural numbers , this defines a complex square matrix on as follows. Let be the matrix on given by
For and in , the defined entry is
The sum is zero when . When , the empty product gives the one-entry matrix with entry one. Empty is allowed, with when .
QAlgorithms.Papers.Independent.Intakethird100.P65.depolarizingValue
For natural numbers and real numbers , this defines the real infimum
where
Division by zero is zero and ; no restriction on or positivity assumption on 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 , giving value zero; for , the value is one when and zero when .
QAlgorithms.Papers.Independent.Intakethird100.P65.fixedSmoothInfo
For finite types with decidable equality, any complex matrix indexed by , and any real , this is the extended-real infimum
where is indexed by , by , and . Here, with denoting continuous functional calculus with the real square-root function,
No state assumption is imposed on ; the displayed equality requires the entire -marginal to be preserved. Zero is permitted when it satisfies these conditions. Real square roots of negative arguments are zero; an empty feasible set gives , including for , and an unbounded-below feasible set gives .
QAlgorithms.Papers.Independent.Intakethird100.P65.flaggedErasureChannel
For every finite type and every real , this defines a map from complex square matrices on to complex square matrices on , where is a new distinguished index, by
No positivity or normalization is required of , and may be any real number. If is empty, the output is the zero matrix on the single index .
QAlgorithms.Papers.Independent.Intakethird100.P65.generalErrorFeasible
For finite index types , with decidable equality, a Boolean flag , arbitrary maps and , and any real , this predicate asserts the existence of complex square matrices on and on satisfying the following conditions. Write for the matrix unit, for the identity, and for positive semidefiniteness. The matrix is positive semidefinite and, for all appropriate indices, obeys
Define matrices on by
The remaining conditions are
where ; additionally, if , the matrix with entries is positive semidefinite. When the flag is false this additional requirement is absent. Neither map is assumed to be a channel or linear; 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 with decidable equality and arbitrary functions , this defines
where and
Here . The test matrices 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 with decidable equality and an arbitrary set of complex matrices on , this defines the extended-real infimum
where is a matrix on and
with this last infimum taken in the extended reals by embedding each admissible real . The set is not required to contain only density matrices or positive matrices. Empty register types and empty are allowed; an infimum over an empty set is , while an admissible set of exponents unbounded below has infimum .
QAlgorithms.Papers.Independent.Intakethird100.P65.linkChoi
For finite index types , a complex matrix indexed by , and a complex matrix indexed by , this defines a complex matrix indexed by by
for every and . 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 .
QAlgorithms.Papers.Independent.Intakethird100.P65.maxEntangledInput
For every finite type equipped with decidable equality, this defines a complex matrix indexed by . Writing , its entries are
For , this is the outer product of with itself. The definition also includes empty : then 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 with decidable equality and any complex matrix indexed by , define its -marginal by . This definition returns the following infimum in the extended real numbers:
Here . No positivity, nonzero, or trace assumption is imposed on . An empty feasible set gives , and a feasible set unbounded below gives ; in particular, if is empty, there is no trace-one .
QAlgorithms.Papers.Independent.Intakethird100.P65.minimumError
For finite index types , with decidable equality, an assistance label , and arbitrary matrix maps and , this defines the real infimum
where ranges over all complex matrices on , not merely density matrices, , amplification applies the map separately to each block indexed by the first factor, and
The infimum ranges over maps satisfying the rule selected by : for , is a channel and its two output marginals depend only on the corresponding input marginals through channels; for , additionally is a channel, with transposing the second tensor factor; for , 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 . Both extrema are taken in the real numbers, and no attainment is asserted; in particular, an empty infimum is assigned the real value .
QAlgorithms.Papers.Independent.Intakethird100.P65.oneUseSimulationCost
For finite index types , with decidable equality, an assistance label , an arbitrary matrix map , and any real tolerance , this defines the real infimum of over positive integers satisfying
Here , ranges over all matrices on , , and amplification applies the map blockwise. “Allowed” means: for , is a channel whose two output marginals are channels of the corresponding input marginals; for , also its conjugation by partial transposition on the second input and second output factors is a channel; for , 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 . There is no assumption that is a channel or that , and the inner infimum need not be attained. The outer infimum is real-valued and equals if its defining set is empty.
QAlgorithms.Papers.Independent.Intakethird100.P65.purify
For a finite type with decidable equality and any complex matrix indexed by , let be continuous functional calculus applied to with the real square-root function, and form the vector indexed by with . This definition returns the matrix
No positivity, normalization, or nonzero assumption is made on ; the vector is not normalized by the definition. Empty is included.
QAlgorithms.Papers.Independent.Intakethird100.P65.reducedErrorFeasible
For finite index types , with decidable equality, a Boolean flag , an arbitrary map , a natural number , and a real , this predicate asserts that there exist complex square matrices on and a Hermitian matrix on such that, writing ,
Here is a matrix unit, means positive semidefinite, and . If , it also requires
where means ; for a false flag these extra inequalities impose no condition. The definition allows , arbitrary , and empty index types, and does not assume that is a channel or linear.
QAlgorithms.Papers.Independent.Intakethird100.P65.referenceSmoothInfo
For finite types with decidable equality, any complex matrix indexed by , and any real , this is the extended-real infimum
where is indexed by , by , and . With denoting continuous functional calculus with the real square-root function, the distance expression is
There is no requirement that , and no positivity or normalization assumption on . Zero is allowed. Real square roots of negative arguments are zero; empty feasibility gives , including for , while unboundedness below gives .
QAlgorithms.Papers.Independent.Intakethird100.P65.roundedInfoCost
For every real , this defines the real number , where the ceiling is the least natural number at least its argument, subsequently regarded as a real number. Since , this ceiling is always at least ; in particular, the value is for every .
QAlgorithms.Papers.Independent.Intakethird100.P65.roundedTraceCost
For every real , this defines , with the natural-number ceiling subsequently regarded as a real number. The real square root and logarithm are total: when , and . Consequently this definition gives for every , including negative and .
QAlgorithms.Papers.Independent.Intakethird100.P65.roundingCorrection
For every real , this defines , where the ceiling is taken in the natural numbers and then converted to a real number. All real are included; for , the ceiling is and this value is .
QAlgorithms.Papers.Independent.Intakethird100.P65.sectorEigenvalue
For arbitrary natural numbers and an arbitrary real number , this defines
There is no assumption that , , or . The exponent subtraction is truncated at zero, division by zero is zero, and . In particular, for the second factor is one, and for the entire expression is one.
QAlgorithms.Papers.Independent.Intakethird100.P65.sectorMatrix
For a finite type with decidable equality, a natural number , and arbitrary real coefficients , this defines a complex matrix on by
where denotes the functions from to , tensor factors are grouped into the two registers, and
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 may be empty, and is included: the empty tensor product is the scalar matrix , so the result is the scalar matrix .
QAlgorithms.Papers.Independent.Intakethird100.P65.sectorWeight
For arbitrary natural numbers , this defines the real number
The subtraction in the exponent is truncated natural-number subtraction. No positivity assumption on , or restriction , is imposed. Division by zero has value zero, and every zeroth power, including , is one. Thus the value is zero when ; when and , it is zero for every ; and at , it is one for every .
QAlgorithms.Papers.Independent.Intakethird100.P65.smoothChannelDmax
For finite index types with decidable equality, arbitrary matrix-valued functions , and any real , this is the infimum in of the values over all channels satisfying . Here a channel means that there exist a natural number and matrices , , such that and for every ; is the matrix with entries , where has its sole nonzero entry at ; and
The first infimum is taken in the extended reals. In the supremum, ranges over all complex matrices indexed by , , and the amplification acts on each block by . Neither nor is assumed to be a channel, need not be nonnegative, and empty index types are allowed. An empty feasible set for either extended-real infimum gives ; attainment is not part of this definition.
QAlgorithms.Papers.Independent.Intakethird100.P65.smoothChannelInfo
For finite types with decidable equality, an arbitrary function , and any real , this defines the real infimum of over all finite Kraus channels whose half-distance from is at most . Such a channel has, for some natural number , matrices satisfying and for every matrix . Its half-distance is
with and a real supremum. To define , set , , and take the extended-real infimum of all real for which some positive semidefinite trace-one matrix on satisfies ; is its real conversion, sending either infinity to zero. No channel assumption is made on . 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 with decidable equality, any complex matrix indexed by , and any real , this is the extended-real infimum
where is indexed by , by , and . Explicitly, writing for continuous functional calculus with the real square-root function,
The matrix need not be positive or normalized, and is allowed. Real square roots of negative arguments are zero. Empty feasibility gives , including when , while unboundedness below gives .
QAlgorithms.Papers.Independent.Intakethird100.P65.smoothRobustness
For finite index types with decidable equality, an arbitrary matrix-valued function , and any real , this is the real infimum of over all channels satisfying , where
Both infima here are real infima. A channel means a function admitting a finite representation with . The distance is
with amplification applying the function to each block indexed by the reference coordinates in . There is no channel assumption on and no nonnegativity assumption on ; empty index types are included. Either real infimum equals if its feasible set is empty, and no minimizing channel is asserted to exist.
QAlgorithms.Papers.Independent.Intakethird100.P65.traceFeasible
For finite index types , with decidable equality, an arbitrary map , a real , complex square matrices on , and a complex square matrix on , this predicate means precisely
Here , is a matrix unit, , and denotes positive semidefiniteness. There is no trace constraint on , no assumption that , and no channel or linearity assumption on ; empty index types are included.
QAlgorithms.Papers.Independent.Intakethird100.P65.traceValue
For finite index types , with decidable equality, an arbitrary map , and any real , this defines the real infimum of over all complex matrices on and on satisfying
where , is a matrix unit, and . Positive semidefiniteness is denoted by . The definition asserts neither feasibility nor attainment, allows negative and empty index types, and imposes no channel or linearity assumption on . 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 , with decidable equality, an assistance label , and an arbitrary matrix map , this defines the extended-real supremum over real numbers such that, for every and , all sufficiently large positive integers admit a positive integer satisfying and an infimum error at most for simulating the identity on matrices using . Here ; the error is the real infimum over allowed of
where , the supremum ranges over all matrices on , amplification is blockwise, and . Allowed maps are channels with both marginal nonsignalling conditions for ; additionally their second-factor partial-transpose conjugates are channels for ; or they use a shared finite-dimensional density matrix and separate local channels for . Channels have finite complete Kraus representations. The inner real infimum need not be attained and equals if empty. The outer supremum is if its rate set is empty and if that set is unbounded above.
QAlgorithms.Papers.Independent.Intakethird100.P65.vanishingErrorSimulationCost
For finite index types , with decidable equality, an assistance label , and an arbitrary matrix map , this defines an extended-real infimum over real numbers with the following property: for every and , all sufficiently large positive integers admit a positive integer such that and the infimum simulation error from the identity map on matrices to is at most . Explicitly, , and that error is the real infimum over allowed of
The supremum uses all matrices on ; amplification is blockwise and . Allowed maps are channels with both marginal nonsignalling conditions for ; additionally their second-factor partial-transpose conjugates are channels for ; or they are implemented by separate local channels with a shared finite-dimensional density matrix for . Channels have finite complete Kraus representations. Inner error infima are real and need not be attained, with empty infimum ; the outer infimum is in the extended reals and is if no such rate exists.
QAlgorithms.Papers.Independent.Intakethird100.P65.wire
For finite index types , with decidable equality, and arbitrary matrix maps and , this defines a matrix map by
for every complex matrix on and every . Matrix units have a single entry , and the tensor product has entries . Neither nor 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 , with decidable equality, and an arbitrary matrix map , this defines the real infimum of all real numbers such that, for every slack , all sufficiently large positive integers satisfy , where
Here , , and . Such a 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 . The supremum uses all matrices on with trace norm at most , amplification is blockwise, and . 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 ; no nonnegativity condition is imposed on , and division by is excluded by the quantifier.
QAlgorithms.Papers.Independent.Intakethird100.P65.zeroTraceValue
For finite index types , with decidable equality, and any map , this defines
Here is a matrix unit and denotes positive semidefiniteness. No separate positivity condition on , channel or linearity condition on , 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.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.