Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Simulation of non-Pauli Channels

Definition
QAlgorithms_Papers_Independent_Intakethird100_P13

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

prove2quantum

The objects the formalized statements of Simulation of non-Pauli Channels (arXiv:1706.05384) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • QubitOp: PDF p.3: qubit operators in the fixed Bool computational basis.
  • ResourceOp: PDF pp.3,5: two-qubit operators, Alice/reference first and Bob/output second.
  • QubitMap: PDF p.5: the imported all-operator superoperator type specialized to qubits.
  • BlochCoordinates: PDF p.3, Eqs.(14),(23): three real Bloch/correlation coordinates, not algebraic points.
  • teleportationPauli: PDF pp.2--3: the ordered basis (I,X,Y,Z), with the library's computational basis.
  • FourOutcomeLaw: PDF pp.3,5, Eqs.(9),(36): a normalized four-letter probability law.
  • ClassicalTeleportationLaw: PDF p.5, Eq.(36): p(l|k), output index FIRST; no row-sum condition is imposed.
  • teleportationBellAmplitude: PDF p.3, Eqs.(10)--(13): Bell vectors Phi+,Psi+,Psi-,Phi-. The Y-generated third vector differs from Psi- by an irrelevant phase.
  • teleportationBellProjector: PDF p.3, Eqs.(10)--(13): normalized Bell projectors, phase independent.
  • teleportationBellMixture: PDF p.3, Eq.(24): the actual Bell-diagonal resource with weights w.
  • qubitBellBranch: PDF pp.2,5: unnormalized Bob state after Bell outcome k on C,A. Every branch is retained; there is no division by an outcome probability.
  • noisyTeleport: PDF p.5: the physical noisy Bell/Pauli protocol, averaged over both classical labels. The resource tau and the conditional law p are fixed before rho.
  • Simulable: PDF p.5: simulability by exactly this protocol on every operator, hence with references.
  • augmentedMoment: PDF p.3, Eq.(8): augmented t_(j,i); j=0 is Bob's local Bloch vector.
  • bobVector: PDF p.3, Eq.(8): Bob's local Bloch vector, not Alice's.
  • resourcePauliCorrelation: PDF p.3, Eq.(8): correlations in Alice-row, Bob-column order.
  • labelKroneckerDelta: PDF p.5, Eq.(38): Kronecker delta in the exponent.
  • sCoordinate: PDF p.5, Eq.(38): the six deltas and common 1/4 are essential. Appendix B's summation instruction omits that factor and is not used.
  • fCoordinate: PDF p.5, Eqs.(37)--(39): f_ij = tprime_(j,i) S_ij.
  • bellVertex: PDF p.3, Eq.(23): source correlation vertices; actual Y multipliers have opposite sign.
  • pauliTetrahedron: PDF p.3, Eq.(23): the closed Bell-correlation tetrahedron, not the multiplier tetrahedron.
  • pauliMixture: PDF p.3, Eq.(9): an actual probability mixture of Pauli conjugations.
  • IsPauli: PDF p.3, Eq.(9): Pauli channels always carry a normalized probability law.
  • pauliPointChannel: PDF pp.3,6, Eqs.(15),(52): q is a correlation point, so the Y multiplier is -q2.
  • blochAmplitudeDamping: PDF pp.6--7, Eqs.(47),(51): amplitude damping on all operators. Its physical parameter range is supplied by each theorem.
  • normalizedQubitChoi: PDF pp.3,6, Eqs.(26),(48): normalized Choi resource, reference first.
  • pauliDamping: PDF p.7, Eq.(53): Bool false/true encodes u=0/1; this asserts no uniqueness.
  • attainableVertices: PDF p.7, Eq.(59): eight vertices with paired signs, regular regime only.
  • attainableRegion: PDF pp.7--8, Eq.(59): complete convex attainable q region for 0<gamma<1.
  • axialReplacement: PDF p.8: axial replacement, used to state gamma=1 without 0/0.
  • channelTraceDistance: PDF p.9, Eq.(68): full trace norm supremum over normalized qubit states. For the channel arguments in the statements the domain is nonempty and the image bounded.
  • channelDiamondDistance: PDF p.9, Eq.(72): full norm with a qubit reference, on normalized density inputs.
  • closestPoint: PDF pp.9,14: explicit closest Pauli correlation point, translation deleted.
  • closestPauli: PDF pp.9,14: concrete minimizing channel, with the source's Y-coordinate convention.
  • discardClassical: PDF p.10: forgetting Alice's outcome and always applying I retains all branches.
  • rescaling: PDF p.12, Lemma 20: regular quotient, used only with a strictly positive denominator.
  • contractPoint: PDF p.12, Lemma 22: roots apply to alpha, never to coordinates.
  • regularPauliPoint: PDF pp.12--13, Eq.(C3): q is constructed from S, rather than freely supplied.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_QAlgorithms_Defs_QubitBlochMaps
import Definitions.Def_QAlgorithms_Defs_MatrixDistances

/-! Statements from Cope--Hetzel--Banchi--Pirandola, arXiv:1706.05384v2.
All local constructions use the fixed computational basis false=0, true=1.
Source ambiguities and the operational-capacity interface gap are recorded in HARD.md.
Compilation is pending the mechanical batch handoff. -/

noncomputable section

open scoped BigOperators Matrix ComplexConjugate ComplexOrder

namespace QAlgorithms.Papers.Independent.Intakethird100.P13

attribute [local instance] Classical.propDecidable Classical.decEq

/-- PDF p.3: qubit operators in the fixed Bool computational basis. -/
abbrev QubitOp := Matrix Bool Bool ℂ

/-- PDF pp.3,5: two-qubit operators, Alice/reference first and Bob/output second. -/
abbrev ResourceOp := Matrix (Bool × Bool) (Bool × Bool) ℂ

/-- PDF p.5: the imported all-operator superoperator type specialized to qubits. -/
abbrev QubitMap := QAlgorithms.Superop Bool Bool

/-- PDF p.3, Eqs.(14),(23): three real Bloch/correlation coordinates, not algebraic points. -/
abbrev BlochCoordinates := Fin 3 → ℝ

/-- PDF pp.2--3: the ordered basis (I,X,Y,Z), with the library's computational basis. -/
def teleportationPauli : Fin 4 → QubitOp := ![1, QAlgorithms.QubitChannels.sigma 0,
  QAlgorithms.QubitChannels.sigma 1, QAlgorithms.QubitChannels.sigma 2]

/-- PDF pp.3,5, Eqs.(9),(36): a normalized four-letter probability law. -/
def FourOutcomeLaw (w : Fin 4 → ℝ) : Prop :=
  (∀ k, 0 ≤ w k) ∧ ∑ k, w k = 1

/-- PDF p.5, Eq.(36): p(l|k), output index FIRST; no row-sum condition is imposed. -/
def ClassicalTeleportationLaw (p : Fin 4 → Fin 4 → ℝ) : Prop :=
  ∀ k, FourOutcomeLaw (fun l => p l k)

/-- PDF p.3, Eqs.(10)--(13): Bell vectors Phi+,Psi+,Psi-,Phi-.
The Y-generated third vector differs from Psi- by an irrelevant phase. -/
def teleportationBellAmplitude (k : Fin 4) (c a : Bool) : ℂ :=
  (Real.sqrt 2 : ℂ)⁻¹ * teleportationPauli k a c

/-- PDF p.3, Eqs.(10)--(13): normalized Bell projectors, phase independent. -/
def teleportationBellProjector (k : Fin 4) : ResourceOp :=
  fun x y => teleportationBellAmplitude k x.1 x.2 * conj (teleportationBellAmplitude k y.1 y.2)

/-- PDF p.3, Eq.(24): the actual Bell-diagonal resource with weights w. -/
def teleportationBellMixture (w : Fin 4 → ℝ) : ResourceOp :=
  ∑ k, (w k : ℂ) • teleportationBellProjector k

/-- PDF pp.2,5: unnormalized Bob state after Bell outcome k on C,A.
Every branch is retained; there is no division by an outcome probability. -/
def qubitBellBranch (tau : ResourceOp) (k : Fin 4) (rho : QubitOp) : QubitOp :=
  fun b b2 => ∑ c : Bool, ∑ a : Bool, ∑ d : Bool, ∑ a2 : Bool,
    conj (teleportationBellAmplitude k c a) * rho c d * tau (a,b) (a2,b2) * teleportationBellAmplitude k d a2

/-- PDF p.5: the physical noisy Bell/Pauli protocol, averaged over both classical labels.
The resource tau and the conditional law p are fixed before rho. -/
def noisyTeleport (tau : ResourceOp) (p : Fin 4 → Fin 4 → ℝ) : QubitMap :=
  fun rho => ∑ k, ∑ l, (p l k : ℂ) •
    (teleportationPauli l * qubitBellBranch tau k rho * (teleportationPauli l)ᴴ)

/-- PDF p.5: simulability by exactly this protocol on every operator, hence with references. -/
def Simulable (tau : ResourceOp) (E : QubitMap) : Prop :=
  ∃ p, ClassicalTeleportationLaw p ∧ noisyTeleport tau p = E

/-- PDF p.3, Eq.(8): augmented t_(j,i); j=0 is Bob's local Bloch vector. -/
def augmentedMoment (tau : ResourceOp) (j : Fin 4) (i : Fin 3) : ℝ :=
  ((Matrix.kronecker (teleportationPauli j) (teleportationPauli i.succ) * tau).trace).re

/-- PDF p.3, Eq.(8): Bob's local Bloch vector, not Alice's. -/
def bobVector (tau : ResourceOp) : BlochCoordinates := fun i => augmentedMoment tau 0 i

/-- PDF p.3, Eq.(8): correlations in Alice-row, Bob-column order. -/
def resourcePauliCorrelation (tau : ResourceOp) : Matrix (Fin 3) (Fin 3) ℝ :=
  fun j i => augmentedMoment tau j.succ i

/-- PDF p.5, Eq.(38): Kronecker delta in the exponent. -/
def labelKroneckerDelta (a b : Fin 4) : ℕ := if a = b then 1 else 0

/-- PDF p.5, Eq.(38): the six deltas and common 1/4 are essential.
Appendix B's summation instruction omits that factor and is not used. -/
def sCoordinate (p : Fin 4 → Fin 4 → ℝ) (i : Fin 3) (j : Fin 4) : ℝ :=
  (1 / 4 : ℝ) * ∑ k, ∑ l,
    (-1 : ℝ) ^ (labelKroneckerDelta k 0 + labelKroneckerDelta j 2 + labelKroneckerDelta j 0 + labelKroneckerDelta k j +
      labelKroneckerDelta i.succ l + labelKroneckerDelta 0 l) * p l k

/-- PDF p.5, Eqs.(37)--(39): f_ij = tprime_(j,i) S_ij. -/
def fCoordinate (tau : ResourceOp) (p : Fin 4 → Fin 4 → ℝ)
    (i : Fin 3) (j : Fin 4) : ℝ :=
  augmentedMoment tau j i * sCoordinate p i j

/-- PDF p.3, Eq.(23): source correlation vertices; actual Y multipliers have opposite sign. -/
def bellVertex : Fin 4 → BlochCoordinates :=
  ![![1,-1,1], ![1,1,-1], ![-1,-1,-1], ![-1,1,1]]

/-- PDF p.3, Eq.(23): the closed Bell-correlation tetrahedron, not the multiplier tetrahedron. -/
def pauliTetrahedron : Set BlochCoordinates := convexHull ℝ (Set.range bellVertex)

/-- PDF p.3, Eq.(9): an actual probability mixture of Pauli conjugations. -/
def pauliMixture (w : Fin 4 → ℝ) : QubitMap :=
  fun rho => ∑ k, (w k : ℂ) • (teleportationPauli k * rho * (teleportationPauli k)ᴴ)

/-- PDF p.3, Eq.(9): Pauli channels always carry a normalized probability law. -/
def IsPauli (E : QubitMap) : Prop :=
  ∃ w, FourOutcomeLaw w ∧ E = pauliMixture w

/-- PDF pp.3,6, Eqs.(15),(52): q is a correlation point, so the Y multiplier is -q2. -/
def pauliPointChannel (q : BlochCoordinates) : QubitMap :=
  QAlgorithms.QubitChannels.blochAffineMap 0 (Matrix.diagonal ![q 0, -q 1, q 2])

/-- PDF pp.6--7, Eqs.(47),(51): amplitude damping on all operators.
Its physical parameter range is supplied by each theorem. -/
def blochAmplitudeDamping (gamma : ℝ) : QubitMap :=
  QAlgorithms.QubitChannels.blochAffineMap ![0,0,gamma]
    (Matrix.diagonal ![Real.sqrt (1-gamma), Real.sqrt (1-gamma), 1-gamma])

/-- PDF pp.3,6, Eqs.(26),(48): normalized Choi resource, reference first. -/
def normalizedQubitChoi (E : QubitMap) : ResourceOp :=
  idTensorSuperop (R := Bool) E (teleportationBellProjector 0)

/-- PDF p.7, Eq.(53): Bool false/true encodes u=0/1; this asserts no uniqueness. -/
def pauliDamping (u : Bool) (eta : ℝ) (q : BlochCoordinates) : QubitMap :=
  fun rho => let Y := blochAmplitudeDamping eta (pauliPointChannel q rho)
    if u then pauliX * Y * pauliXᴴ else Y

/-- PDF p.7, Eq.(59): eight vertices with paired signs, regular regime only. -/
def attainableVertices (gamma eta : ℝ) : Fin 8 → BlochCoordinates :=
  let a := Real.sqrt ((1-gamma)/(1-eta))
  let b := 1-eta/gamma
  ![![a,a*b,-a^2*b], ![a,-a*b,a^2*b],
    ![a*b,a,-a^2*b], ![-a*b,a,a^2*b],
    ![-a,a*b,a^2*b], ![-a,-a*b,-a^2*b],
    ![a*b,-a,a^2*b], ![-a*b,-a,-a^2*b]]

/-- PDF pp.7--8, Eq.(59): complete convex attainable q region for 0<gamma<1. -/
def attainableRegion (gamma eta : ℝ) : Set BlochCoordinates :=
  convexHull ℝ (Set.range (attainableVertices gamma eta))

/-- PDF p.8: axial replacement, used to state gamma=1 without 0/0. -/
def axialReplacement (z : ℝ) : QubitMap :=
  QAlgorithms.QubitChannels.blochAffineMap ![0,0,z] 0

/-- PDF p.9, Eq.(68): full trace norm supremum over normalized qubit states.
For the channel arguments in the statements the domain is nonempty and the image bounded. -/
def channelTraceDistance (E F : QubitMap) : ℝ :=
  sSup {v | ∃ rho : QubitOp, IsDensityMatrix rho ∧ v = traceNorm (E rho - F rho)}

/-- PDF p.9, Eq.(72): full norm with a qubit reference, on normalized density inputs. -/
def channelDiamondDistance (E F : QubitMap) : ℝ :=
  sSup {v | ∃ rho : ResourceOp, IsDensityMatrix rho ∧
    v = traceNorm (idTensorSuperop (R := Bool) E rho -
      idTensorSuperop (R := Bool) F rho)}

/-- PDF pp.9,14: explicit closest Pauli correlation point, translation deleted. -/
def closestPoint (u : Bool) (eta : ℝ) (q : BlochCoordinates) : BlochCoordinates :=
  let a := Real.sqrt (1-eta)
  if u then ![a*q 0, -a*q 1, -(1-eta)*q 2]
  else ![a*q 0, a*q 1, (1-eta)*q 2]

/-- PDF pp.9,14: concrete minimizing channel, with the source's Y-coordinate convention. -/
def closestPauli (u : Bool) (eta : ℝ) (q : BlochCoordinates) : QubitMap :=
  pauliPointChannel (closestPoint u eta q)

/-- PDF p.10: forgetting Alice's outcome and always applying I retains all branches. -/
def discardClassical : Fin 4 → Fin 4 → ℝ := fun l _ => if l = 0 then 1 else 0

/-- PDF p.12, Lemma 20: regular quotient, used only with a strictly positive denominator. -/
def rescaling (gamma s30 : ℝ) : ℝ := (1-gamma)/(1-gamma*s30)

/-- PDF p.12, Lemma 22: roots apply to alpha, never to coordinates. -/
def contractPoint (alpha : ℝ) (q : BlochCoordinates) : BlochCoordinates :=
  ![Real.sqrt alpha*q 0, Real.sqrt alpha*q 1, alpha*q 2]

/-- PDF pp.12--13, Eq.(C3): q is constructed from S, rather than freely supplied. -/
def regularPauliPoint (gamma s30 : ℝ) (s : BlochCoordinates) : BlochCoordinates :=
  contractPoint (rescaling gamma s30) s

/-- Definition sanity: source's first resourcePauliCorrelation vertex. -/
example : bellVertex 0 = ![1,-1,1] := rfl

/-- Definition sanity: normalized classical reset, including all four input columns. -/
example : ClassicalTeleportationLaw discardClassical := by
  intro k
  constructor
  · intro l; by_cases h : l = 0 <;> simp [discardClassical, h]
  · simp [discardClassical]

/-- Definition sanity: at alpha=0 the contraction is the origin. -/
example (q : BlochCoordinates) : contractPoint 0 q = 0 := by
  funext i
  fin_cases i <;> simp [contractPoint]

/-- Definition sanity: no singular quotient is needed at zero damping. -/
example (s : ℝ) : rescaling 0 s = 1 := by simp [rescaling]

/-- Definition sanity: the first coordinate of the regular vertex is computed directly. -/
example (gamma : ℝ) :
    attainableVertices gamma gamma 0 0 = Real.sqrt ((1-gamma)/(1-gamma)) := rfl

/-- Definition sanity: eta=gamma makes b=0 in both remaining coordinates. -/
example (gamma : ℝ) (h : gamma ≠ 0) :
    attainableVertices gamma gamma 0 1 = 0 ∧
    attainableVertices gamma gamma 0 2 = 0 := by
  simp [attainableVertices, h]

/-- Definition sanity: the first Bell vector has positive normalized 00 amplitude. -/
example : teleportationBellAmplitude 0 false false = (Real.sqrt 2 : ℂ)⁻¹ := by
  simp [teleportationBellAmplitude, teleportationPauli, QAlgorithms.QubitChannels.sigma]

/-- Definition sanity: the third Bell vector has the allowed global Y phase. -/
example : teleportationBellAmplitude 2 false true = (Real.sqrt 2 : ℂ)⁻¹ * Complex.I := by
  simp [teleportationBellAmplitude, teleportationPauli, QAlgorithms.QubitChannels.sigma, pauliY]

/-- Definition sanity: Bell-branch averaging is performed before choosing an input. -/
example (tau : ResourceOp) (p : Fin 4 → Fin 4 → ℝ) (rho : QubitOp) :
    noisyTeleport tau p rho = ∑ k, ∑ l, (p l k : ℂ) •
      (teleportationPauli l * qubitBellBranch tau k rho * (teleportationPauli l)ᴴ) := rfl

end QAlgorithms.Papers.Independent.Intakethird100.P13
Source
Simulation of non-Pauli Channels, arXiv:1706.05384v2, https://arxiv.org/abs/1706.05384v2
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P13.BlochCoordinates

Bloch coordinates are a function from {0,1,2}\{0,1,2\}{0,1,2} to R\mathbb RR, equivalently an arbitrary ordered triple of real numbers. No norm or other constraint is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P13.ClassicalTeleportationLaw

For any real-valued function ppp of two indices in {0,1,2,3}\{0,1,2,3\}{0,1,2,3}, this predicate means that, for each fixed second index kkk, its values over the first index form a probability distribution:

plk≥0for every l,k∈{0,1,2,3},∑l=03plk=1for every k.p_{lk}\ge0\quad\text{for every }l,k\in\{0,1,2,3\}, \qquad \sum_{l=0}^{3}p_{lk}=1\quad\text{for every }k.plk​≥0for every l,k∈{0,1,2,3},l=0∑3​plk​=1for every k.

Zero entries are allowed; no condition requiring a sum over kkk to equal one is included.

QAlgorithms.Papers.Independent.Intakethird100.P13.FourOutcomeLaw

For any real-valued function www on {0,1,2,3}\{0,1,2,3\}{0,1,2,3}, this predicate means that every value is nonnegative and the four values sum to one:

wk≥0for every k∈{0,1,2,3},∑k=03wk=1.w_k\ge0\quad\text{for every }k\in\{0,1,2,3\}, \qquad \sum_{k=0}^{3}w_k=1.wk​≥0for every k∈{0,1,2,3},k=0∑3​wk​=1.

Individual values may be zero.

QAlgorithms.Papers.Independent.Intakethird100.P13.IsPauli

For any function EEE from complex 2×22\times22×2 matrices to complex 2×22\times22×2 matrices, IsPauli⁡(E)\operatorname{IsPauli}(E)IsPauli(E) means that there exist real weights w0,w1,w2,w3w_0,w_1,w_2,w_3w0​,w1​,w2​,w3​ satisfying wk≥0w_k\ge0wk​≥0 for every kkk and ∑k=03wk=1\sum_{k=0}^3w_k=1∑k=03​wk​=1, such that, for every complex 2×22\times22×2 matrix ρ\rhoρ,

E(ρ)=∑k=03wkPkρPk†,E(\rho)=\sum_{k=0}^3w_kP_k\rho P_k^\dagger,E(ρ)=k=0∑3​wk​Pk​ρPk†​,

where (P0,P1,P2,P3)=(I2,X,Y,Z)(P_0,P_1,P_2,P_3)=(I_2,X,Y,Z)(P0​,P1​,P2​,P3​)=(I2​,X,Y,Z) are the identity and standard Pauli matrices in the ordered Boolean basis. Zero weights are permitted, and the equality is required on all input matrices.

QAlgorithms.Papers.Independent.Intakethird100.P13.QubitMap

A qubit map is an arbitrary function from complex 2×22\times22×2 matrices indexed by Boolean values to matrices of the same kind. Linearity, positivity, and trace preservation are not required.

QAlgorithms.Papers.Independent.Intakethird100.P13.QubitOp

A qubit operator is an arbitrary complex matrix whose rows and columns are indexed by the two Boolean values, equivalently a 2×22\times22×2 complex matrix. No Hermiticity, positivity, trace, or unitarity condition is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P13.ResourceOp

A resource operator is an arbitrary complex matrix whose rows and columns are indexed by ordered pairs of Boolean values, equivalently a 4×44\times44×4 complex matrix. No Hermiticity, positivity, trace, or unitarity condition is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P13.Simulable

For any complex matrix τ\tauτ indexed by {0,1}2\{0,1\}^2{0,1}2 and any function EEE from complex 2×22\times22×2 matrices to complex 2×22\times22×2 matrices, Simulable⁡(τ,E)\operatorname{Simulable}(\tau,E)Simulable(τ,E) means that there exist real numbers plkp_{lk}plk​, indexed by k,l∈{0,1,2,3}k,l\in\{0,1,2,3\}k,l∈{0,1,2,3}, such that plk≥0p_{lk}\ge0plk​≥0, ∑l=03plk=1\sum_{l=0}^3p_{lk}=1∑l=03​plk​=1 for every kkk, and, for every complex 2×22\times22×2 matrix ρ\rhoρ,

E(ρ)=∑k=03∑l=03plkPlBk(τ,ρ)Pl†.E(\rho)=\sum_{k=0}^3\sum_{l=0}^3p_{lk}P_lB_k(\tau,\rho)P_l^\dagger.E(ρ)=k=0∑3​l=0∑3​plk​Pl​Bk​(τ,ρ)Pl†​.

Here P0=I2P_0=I_2P0​=I2​, P1=XP_1=XP1​=X, P2=YP_2=YP2​=Y, and P3=ZP_3=ZP3​=Z, with X,Y,ZX,Y,ZX,Y,Z the standard Pauli matrices in the ordered Boolean basis (false,true)(\mathrm{false},\mathrm{true})(false,true), identified with (0,1)(0,1)(0,1); writing βk(c,a)=(Pk)ac/2\beta_k(c,a)=(P_k)_{ac}/\sqrt2βk​(c,a)=(Pk​)ac​/2​, the matrix BkB_kBk​ is defined entrywise by

[Bk(τ,ρ)]bb′=∑c,a,d,a′∈{0,1}βk(c,a)‾ ρcd τ(a,b),(a′,b′) βk(d,a′).[B_k(\tau,\rho)]_{bb'} =\sum_{c,a,d,a'\in\{0,1\}} \overline{\beta_k(c,a)}\,\rho_{cd}\, \tau_{(a,b),(a',b')}\,\beta_k(d,a').[Bk​(τ,ρ)]bb′​=c,a,d,a′∈{0,1}∑​βk​(c,a)​ρcd​τ(a,b),(a′,b′)​βk​(d,a′).

No positivity or trace condition on τ\tauτ, and no separate channel assumption on EEE, is included in this definition.

QAlgorithms.Papers.Independent.Intakethird100.P13.attainableRegion

For every pair of real numbers γ,η\gamma,\etaγ,η, this defines the real convex hull of the following eight triples:

(a,ab,−a2b), (a,−ab,a2b), (ab,a,−a2b), (−ab,a,a2b), (−a,ab,a2b), (−a,−ab,−a2b), (ab,−a,a2b), (−ab,−a,−a2b),(a,ab,-a^2b),\ (a,-ab,a^2b),\ (ab,a,-a^2b),\ (-ab,a,a^2b),\ (-a,ab,a^2b),\ (-a,-ab,-a^2b),\ (ab,-a,a^2b),\ (-ab,-a,-a^2b),(a,ab,−a2b), (a,−ab,a2b), (ab,a,−a2b), (−ab,a,a2b), (−a,ab,a2b), (−a,−ab,−a2b), (ab,−a,a2b), (−ab,−a,−a2b),

where a=(1−γ)/(1−η)a=\sqrt{(1-\gamma)/(1-\eta)}a=(1−γ)/(1−η)​ and b=1−η/γb=1-\eta/\gammab=1−η/γ. Equivalently, the resulting set consists of all sums ∑j=07tjvj\sum_{j=0}^7 t_jv_j∑j=07​tj​vj​ of these triples with real coefficients tj≥0t_j\ge0tj​≥0 and ∑j=07tj=1\sum_{j=0}^7t_j=1∑j=07​tj​=1. No restrictions are placed on the parameters; division by zero gives zero and the real square root of a nonpositive number is zero, so repeated or collapsed triples are included.

QAlgorithms.Papers.Independent.Intakethird100.P13.attainableVertices

For every pair of real numbers γ,η\gamma,\etaγ,η, define

a=1−γ1−η,b=1−ηγ.a=\sqrt{\frac{1-\gamma}{1-\eta}},\qquad b=1-\frac{\eta}{\gamma}.a=1−η1−γ​​,b=1−γη​.

This declaration defines a function from {0,…,7}\{0,\ldots,7\}{0,…,7} to real triples, assigning, in that order,

0↦(a,ab,−a2b),1↦(a,−ab,a2b),2↦(ab,a,−a2b),3↦(−ab,a,a2b),4↦(−a,ab,a2b),5↦(−a,−ab,−a2b),6↦(ab,−a,a2b),7↦(−ab,−a,−a2b).\begin{aligned} 0&\mapsto(a,ab,-a^2b),& 1&\mapsto(a,-ab,a^2b),\\ 2&\mapsto(ab,a,-a^2b),& 3&\mapsto(-ab,a,a^2b),\\ 4&\mapsto(-a,ab,a^2b),& 5&\mapsto(-a,-ab,-a^2b),\\ 6&\mapsto(ab,-a,a^2b),& 7&\mapsto(-ab,-a,-a^2b). \end{aligned}0246​↦(a,ab,−a2b),↦(ab,a,−a2b),↦(−a,ab,a2b),↦(ab,−a,a2b),​1357​↦(a,−ab,a2b),↦(−ab,a,a2b),↦(−a,−ab,−a2b),↦(−ab,−a,−a2b).​

The parameters are unrestricted, division by zero is defined to be zero, and the real square root is zero for nonpositive arguments. In particular, η=1\eta=1η=1 gives a=0a=0a=0, and γ=0\gamma=0γ=0 gives b=1b=1b=1; the eight listed triples need not be distinct.

QAlgorithms.Papers.Independent.Intakethird100.P13.augmentedMoment

For any complex matrix τ\tauτ indexed by {0,1}2\{0,1\}^2{0,1}2, any j∈{0,1,2,3}j\in\{0,1,2,3\}j∈{0,1,2,3}, and any i∈{0,1,2}i\in\{0,1,2\}i∈{0,1,2}, the real number called the augmented moment is

augmentedMoment⁡(τ,j,i)=Re⁡Tr⁡ ⁣[(Pj⊗Pi+1)τ],\operatorname{augmentedMoment}(\tau,j,i) =\operatorname{Re}\operatorname{Tr}\!\left[(P_j\otimes P_{i+1})\tau\right],augmentedMoment(τ,j,i)=ReTr[(Pj​⊗Pi+1​)τ],

where (P0,P1,P2,P3)=(I2,X,Y,Z)(P_0,P_1,P_2,P_3)=(I_2,X,Y,Z)(P0​,P1​,P2​,P3​)=(I2​,X,Y,Z), with X,Y,ZX,Y,ZX,Y,Z the standard Pauli matrices in the Boolean basis ordered as (false,true)(\mathrm{false},\mathrm{true})(false,true). The matrix τ\tauτ is arbitrary; it is not required to be a density matrix.

QAlgorithms.Papers.Independent.Intakethird100.P13.axialReplacement

For every real number zzz, this defines a function on all complex 2×22\times22×2 matrices by

ρ⟼tr⁡ρ2(1+z001−z),\rho\longmapsto\frac{\operatorname{tr}\rho}{2} \begin{pmatrix}1+z&0\\0&1-z\end{pmatrix},ρ⟼2trρ​(1+z0​01−z​),

where rows and columns are ordered false, true. There is no restriction on zzz or on the input matrix ρ\rhoρ.

QAlgorithms.Papers.Independent.Intakethird100.P13.bellVertex

The function bellVertex⁡\operatorname{bellVertex}bellVertex assigns a real three-component vector, with coordinates ordered 0,1,20,1,20,1,2, to each label in {0,1,2,3}\{0,1,2,3\}{0,1,2,3} as follows:

bellVertex⁡(0)=(1,−1,1),bellVertex⁡(1)=(1,1,−1),bellVertex⁡(2)=(−1,−1,−1),bellVertex⁡(3)=(−1,1,1).\operatorname{bellVertex}(0)=(1,-1,1),\qquad \operatorname{bellVertex}(1)=(1,1,-1),\qquad \operatorname{bellVertex}(2)=(-1,-1,-1),\qquad \operatorname{bellVertex}(3)=(-1,1,1).bellVertex(0)=(1,−1,1),bellVertex(1)=(1,1,−1),bellVertex(2)=(−1,−1,−1),bellVertex(3)=(−1,1,1).

QAlgorithms.Papers.Independent.Intakethird100.P13.blochAmplitudeDamping

For every real number γ\gammaγ, this defines a function Dγ\mathcal D_\gammaDγ​ from complex 2×22\times22×2 matrices to complex 2×22\times22×2 matrices by

Dγ(ρ)=tr⁡ρ2(I+γZ)+1−γ2tr⁡(Xρ)X+1−γ2tr⁡(Yρ)Y+1−γ2tr⁡(Zρ)Z,\mathcal D_\gamma(\rho)=\frac{\operatorname{tr}\rho}{2}(I+\gamma Z) +\frac{\sqrt{1-\gamma}}{2}\operatorname{tr}(X\rho)X +\frac{\sqrt{1-\gamma}}{2}\operatorname{tr}(Y\rho)Y +\frac{1-\gamma}{2}\operatorname{tr}(Z\rho)Z,Dγ​(ρ)=2trρ​(I+γZ)+21−γ​​tr(Xρ)X+21−γ​​tr(Yρ)Y+21−γ​tr(Zρ)Z,

where, with rows and columns ordered as false, true,

X=(0110),Y=(0−ii0),Z=(100−1),X=\begin{pmatrix}0&1\\1&0\end{pmatrix},\qquad Y=\begin{pmatrix}0&-i\\i&0\end{pmatrix},\qquad Z=\begin{pmatrix}1&0\\0&-1\end{pmatrix},X=(01​10​),Y=(0i​−i0​),Z=(10​0−1​),

and III is the identity matrix. There is no restriction on γ\gammaγ or ρ\rhoρ; the real square root is zero for nonpositive arguments, so 1−γ=0\sqrt{1-\gamma}=01−γ​=0 when γ≥1\gamma\ge1γ≥1.

QAlgorithms.Papers.Independent.Intakethird100.P13.bobVector

For any complex matrix τ\tauτ indexed by {0,1}2\{0,1\}^2{0,1}2, its defined vector bobVector⁡(τ)\operatorname{bobVector}(\tau)bobVector(τ) is the real three-component vector whose coordinates are

[bobVector⁡(τ)]i=Re⁡Tr⁡ ⁣[(I2⊗σi)τ],i∈{0,1,2},[\operatorname{bobVector}(\tau)]_i =\operatorname{Re}\operatorname{Tr}\!\left[(I_2\otimes \sigma_i)\tau\right], \qquad i\in\{0,1,2\},[bobVector(τ)]i​=ReTr[(I2​⊗σi​)τ],i∈{0,1,2},

where (σ0,σ1,σ2)=(X,Y,Z)(\sigma_0,\sigma_1,\sigma_2)=(X,Y,Z)(σ0​,σ1​,σ2​)=(X,Y,Z) are the standard Pauli matrices in the ordered Boolean basis. No positivity or normalization condition on τ\tauτ is assumed.

QAlgorithms.Papers.Independent.Intakethird100.P13.channelDiamondDistance

For any two functions E,FE,FE,F from complex 2×22\times22×2 matrices to complex 2×22\times22×2 matrices, define their blockwise extensions to matrices indexed by Boolean pairs by

[E^(ρ)](a,b),(c,d)=[E ⁣((ρ(a,i),(c,j))i,j)]bd,[\widehat E(\rho)]_{(a,b),(c,d)} =\left[E\!\left((\rho_{(a,i),(c,j)})_{i,j}\right)\right]_{bd},[E(ρ)](a,b),(c,d)​=[E((ρ(a,i),(c,j)​)i,j​)]bd​,

and define F^\widehat FF in the same way; every index here is Boolean. This declaration defines

sup⁡ρ⪰0tr⁡ρ=1Re⁡tr⁡(E^(ρ)−F^(ρ))†(E^(ρ)−F^(ρ)),\sup_{\substack{\rho\succeq0\\\operatorname{tr}\rho=1}} \operatorname{Re}\operatorname{tr} \sqrt{(\widehat E(\rho)-\widehat F(\rho))^\dagger (\widehat E(\rho)-\widehat F(\rho))},ρ⪰0trρ=1​sup​Retr(E(ρ)−F(ρ))†(E(ρ)−F(ρ))​,

where ρ\rhoρ ranges over positive semidefinite complex 4×44\times44×4 matrices of trace one, †\dagger† is conjugate transpose, and the square root is the positive semidefinite matrix square root obtained by functional calculus. No linearity or channel assumptions are imposed on E,FE,FE,F; the auxiliary index has exactly two values and there is no factor of 1/21/21/2. The real supremum convention returns zero if the set of values is unbounded above.

QAlgorithms.Papers.Independent.Intakethird100.P13.channelTraceDistance

For any two functions E,FE,FE,F from complex 2×22\times22×2 matrices to complex 2×22\times22×2 matrices, this defines the real supremum

sup⁡ρ⪰0tr⁡ρ=1Re⁡tr⁡(E(ρ)−F(ρ))†(E(ρ)−F(ρ)),\sup_{\substack{\rho\succeq0\\\operatorname{tr}\rho=1}} \operatorname{Re}\operatorname{tr} \sqrt{(E(\rho)-F(\rho))^\dagger(E(\rho)-F(\rho))},ρ⪰0trρ=1​sup​Retr(E(ρ)−F(ρ))†(E(ρ)−F(ρ))​,

where ρ\rhoρ ranges over all positive semidefinite complex 2×22\times22×2 matrices of trace one, †\dagger† denotes conjugate transpose, and the square root is the positive semidefinite matrix square root obtained by functional calculus. No linearity or channel assumptions are imposed on EEE or FFF, and there is no factor of 1/21/21/2. The real supremum convention returns zero if the set of values is unbounded above.

QAlgorithms.Papers.Independent.Intakethird100.P13.closestPauli

For every Boolean uuu, every real number η\etaη, and every real triple q=(q0,q1,q2)q=(q_0,q_1,q_2)q=(q0​,q1​,q2​), set a=1−ηa=\sqrt{1-\eta}a=1−η​ and let s=−1s=-1s=−1 if uuu is true and s=1s=1s=1 if uuu is false. This defines the function on complex 2×22\times22×2 matrices

ρ⟼tr⁡ρ2I+12(aq0tr⁡(Xρ)X−saq1tr⁡(Yρ)Y+s(1−η)q2tr⁡(Zρ)Z),\rho\longmapsto \frac{\operatorname{tr}\rho}{2}I+ \frac12\left( aq_0\operatorname{tr}(X\rho)X -sa q_1\operatorname{tr}(Y\rho)Y +s(1-\eta)q_2\operatorname{tr}(Z\rho)Z \right),ρ⟼2trρ​I+21​(aq0​tr(Xρ)X−saq1​tr(Yρ)Y+s(1−η)q2​tr(Zρ)Z),

where III is the identity and, in false, true index order,

X=(0110),Y=(0−ii0),Z=(100−1).X=\begin{pmatrix}0&1\\1&0\end{pmatrix},\qquad Y=\begin{pmatrix}0&-i\\i&0\end{pmatrix},\qquad Z=\begin{pmatrix}1&0\\0&-1\end{pmatrix}.X=(01​10​),Y=(0i​−i0​),Z=(10​0−1​).

There are no restrictions on η,q,ρ\eta,q,\rhoη,q,ρ, and the real square root is zero for nonpositive arguments.

QAlgorithms.Papers.Independent.Intakethird100.P13.closestPoint

For every Boolean uuu, every real number η\etaη, and every real triple q=(q0,q1,q2)q=(q_0,q_1,q_2)q=(q0​,q1​,q2​), let a=1−ηa=\sqrt{1-\eta}a=1−η​. This defines the real triple

{(aq0,−aq1,−(1−η)q2),u=true,(aq0,aq1,(1−η)q2),u=false.\begin{cases} (aq_0,-aq_1,-(1-\eta)q_2),&u=\mathrm{true},\\ (aq_0,aq_1,(1-\eta)q_2),&u=\mathrm{false}. \end{cases}{(aq0​,−aq1​,−(1−η)q2​),(aq0​,aq1​,(1−η)q2​),​u=true,u=false.​

The inputs are unrestricted; the real square root is zero for nonpositive arguments, so a=0a=0a=0 whenever η≥1\eta\ge1η≥1.

QAlgorithms.Papers.Independent.Intakethird100.P13.contractPoint

For every real number α\alphaα and every vector q=(q0,q1,q2)∈R3q=(q_0,q_1,q_2)\in\mathbb R^3q=(q0​,q1​,q2​)∈R3, the definition gives

contractPoint⁡(α,q)=(α q0,α q1,αq2).\operatorname{contractPoint}(\alpha,q) =(\sqrt{\alpha}\,q_0,\sqrt{\alpha}\,q_1,\alpha q_2).contractPoint(α,q)=(α​q0​,α​q1​,αq2​).

There is no restriction on α\alphaα or qqq; the real square root used here is 000 when α<0\alpha<0α<0.

QAlgorithms.Papers.Independent.Intakethird100.P13.discardClassical

This defines a real-valued function of two indices l,k∈{0,1,2,3}l,k\in\{0,1,2,3\}l,k∈{0,1,2,3} by

p(l,k)={1,l=0,0,l≠0.p(l,k)= \begin{cases} 1,&l=0,\\ 0,&l\ne0. \end{cases}p(l,k)={1,0,​l=0,l=0.​

Its value is independent of the second index kkk.

QAlgorithms.Papers.Independent.Intakethird100.P13.fCoordinate

For any complex matrix τ\tauτ indexed by {0,1}2\{0,1\}^2{0,1}2, any real array p:{0,1,2,3}2→Rp:\{0,1,2,3\}^2\to\mathbb Rp:{0,1,2,3}2→R, any i∈{0,1,2}i\in\{0,1,2\}i∈{0,1,2}, and any j∈{0,1,2,3}j\in\{0,1,2,3\}j∈{0,1,2,3}, the defined real coordinate is

fij(τ,p)=Re⁡Tr⁡ ⁣[(Pj⊗Pi+1)τ]⋅14∑k=03∑l=03(−1)δk0+δj2+δj0+δkj+δi+1,l+δ0l plk,f_{ij}(\tau,p)= \operatorname{Re}\operatorname{Tr}\!\left[(P_j\otimes P_{i+1})\tau\right] \cdot\frac14\sum_{k=0}^3\sum_{l=0}^3 (-1)^{\delta_{k0}+\delta_{j2}+\delta_{j0}+\delta_{kj} +\delta_{i+1,l}+\delta_{0l}}\,p_{lk},fij​(τ,p)=ReTr[(Pj​⊗Pi+1​)τ]⋅41​k=0∑3​l=0∑3​(−1)δk0​+δj2​+δj0​+δkj​+δi+1,l​+δ0l​plk​,

where (P0,P1,P2,P3)=(I2,X,Y,Z)(P_0,P_1,P_2,P_3)=(I_2,X,Y,Z)(P0​,P1​,P2​,P3​)=(I2​,X,Y,Z), with the standard Pauli matrices in the ordered Boolean basis, and δab\delta_{ab}δab​ is 111 if a=ba=ba=b and 000 otherwise. No positivity or normalization condition on τ\tauτ or ppp is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P13.labelKroneckerDelta

For every pair of labels a,b∈{0,1,2,3}a,b\in\{0,1,2,3\}a,b∈{0,1,2,3}, labelKroneckerDelta⁡(a,b)\operatorname{labelKroneckerDelta}(a,b)labelKroneckerDelta(a,b) is the natural number

δab={1,a=b,0,a≠b.\delta_{ab}= \begin{cases} 1,&a=b,\\ 0,&a\ne b. \end{cases}δab​={1,0,​a=b,a=b.​

QAlgorithms.Papers.Independent.Intakethird100.P13.noisyTeleport

For any complex matrix τ\tauτ indexed by pairs of Boolean values and any real-valued function ppp on {0,1,2,3}2\{0,1,2,3\}^2{0,1,2,3}2, this defines a function Eτ,pE_{\tau,p}Eτ,p​ from Boolean-indexed complex matrices to matrices of the same kind by

Eτ,p(ρ)=∑k=03∑l=03plk PlCτ,k(ρ)Pl†,[Cτ,k(ρ)]b,b′=12∑c,a,d,a′∈{false,true}(Pk)a,c‾ ρc,d τ(a,b),(a′,b′) (Pk)a′,d.E_{\tau,p}(\rho)= \sum_{k=0}^{3}\sum_{l=0}^{3}p_{lk}\, P_l C_{\tau,k}(\rho)P_l^\dagger, \qquad [C_{\tau,k}(\rho)]_{b,b'}= \frac12\sum_{c,a,d,a'\in\{\mathrm{false},\mathrm{true}\}} \overline{(P_k)_{a,c}}\, \rho_{c,d}\, \tau_{(a,b),(a',b')}\, (P_k)_{a',d}.Eτ,p​(ρ)=k=0∑3​l=0∑3​plk​Pl​Cτ,k​(ρ)Pl†​,[Cτ,k​(ρ)]b,b′​=21​c,a,d,a′∈{false,true}∑​(Pk​)a,c​​ρc,d​τ(a,b),(a′,b′)​(Pk​)a′,d​.

Here ρ\rhoρ is any Boolean-indexed complex matrix, b,b′b,b'b,b′ are Boolean values, a bar denotes complex conjugation, †\dagger† denotes conjugate transpose, and, in Boolean order false,true\mathrm{false},\mathrm{true}false,true,

P0=(1001),P1=(0110),P2=(0−ii0),P3=(100−1).P_0=\begin{pmatrix}1&0\\0&1\end{pmatrix},\quad P_1=\begin{pmatrix}0&1\\1&0\end{pmatrix},\quad P_2=\begin{pmatrix}0&-\mathrm i\\ \mathrm i&0\end{pmatrix},\quad P_3=\begin{pmatrix}1&0\\0&-1\end{pmatrix}.P0​=(10​01​),P1​=(01​10​),P2​=(0i​−i0​),P3​=(10​0−1​).

No density-matrix condition on τ\tauτ or ρ\rhoρ, and no nonnegativity or normalization condition on ppp, is assumed.

QAlgorithms.Papers.Independent.Intakethird100.P13.normalizedQubitChoi

For every function EEE from complex 2×22\times22×2 matrices to complex 2×22\times22×2 matrices, this defines a complex 4×44\times44×4 matrix CEC_ECE​, indexed by pairs of Boolean values, by

(CE)(a,b),(c,d)=[E ⁣(12eac)]bd(a,b,c,d∈{false,true}),(C_E)_{(a,b),(c,d)} =\left[E\!\left(\tfrac12 e_{ac}\right)\right]_{bd} \qquad(a,b,c,d\in\{\mathrm{false},\mathrm{true}\}),(CE​)(a,b),(c,d)​=[E(21​eac​)]bd​(a,b,c,d∈{false,true}),

where eace_{ac}eac​ is the 2×22\times22×2 matrix whose only nonzero entry is a 111 in row aaa, column ccc. This is the blockwise application of EEE to the second factor of the matrix ∣Φ⟩⟨Φ∣|\Phi\rangle\langle\Phi|∣Φ⟩⟨Φ∣, where ∣Φ⟩=(∣false,false⟩+∣true,true⟩)/2|\Phi\rangle=(|\mathrm{false},\mathrm{false}\rangle+|\mathrm{true},\mathrm{true}\rangle)/\sqrt2∣Φ⟩=(∣false,false⟩+∣true,true⟩)/2​. No linearity, positivity, or trace preservation is assumed of EEE.

QAlgorithms.Papers.Independent.Intakethird100.P13.pauliDamping

For every Boolean uuu, every real number η\etaη, and every real triple q=(q0,q1,q2)q=(q_0,q_1,q_2)q=(q0​,q1​,q2​), this defines a function on all complex 2×22\times22×2 matrices ρ\rhoρ. Write

X=(0110),Y=(0−ii0),Z=(100−1),X=\begin{pmatrix}0&1\\1&0\end{pmatrix},\quad Y=\begin{pmatrix}0&-i\\i&0\end{pmatrix},\quad Z=\begin{pmatrix}1&0\\0&-1\end{pmatrix},X=(01​10​),Y=(0i​−i0​),Z=(10​0−1​),

with Boolean indices ordered false, true, and set

Aq(ρ)=tr⁡ρ2I+12(q0tr⁡(Xρ)X−q1tr⁡(Yρ)Y+q2tr⁡(Zρ)Z).A_q(\rho)=\frac{\operatorname{tr}\rho}{2}I+ \frac12\bigl(q_0\operatorname{tr}(X\rho)X-q_1\operatorname{tr}(Y\rho)Y +q_2\operatorname{tr}(Z\rho)Z\bigr).Aq​(ρ)=2trρ​I+21​(q0​tr(Xρ)X−q1​tr(Yρ)Y+q2​tr(Zρ)Z).

For any matrix BBB, set

Dη(B)=tr⁡B2(I+ηZ)+1−η2(tr⁡(XB)X+tr⁡(YB)Y)+1−η2tr⁡(ZB)Z.D_\eta(B)=\frac{\operatorname{tr}B}{2}(I+\eta Z) +\frac{\sqrt{1-\eta}}2\bigl(\operatorname{tr}(XB)X+\operatorname{tr}(YB)Y\bigr) +\frac{1-\eta}{2}\operatorname{tr}(ZB)Z.Dη​(B)=2trB​(I+ηZ)+21−η​​(tr(XB)X+tr(YB)Y)+21−η​tr(ZB)Z.

The output is XDη(Aq(ρ))X†XD_\eta(A_q(\rho))X^\daggerXDη​(Aq​(ρ))X† if uuu is true, and Dη(Aq(ρ))D_\eta(A_q(\rho))Dη​(Aq​(ρ)) if uuu is false, where †\dagger† denotes conjugate transpose. There are no restrictions on η\etaη or qqq; the real square root is zero when its argument is nonpositive.

QAlgorithms.Papers.Independent.Intakethird100.P13.pauliMixture

For any function w:{0,1,2,3}→Rw:\{0,1,2,3\}\to\mathbb Rw:{0,1,2,3}→R, pauliMixture⁡(w)\operatorname{pauliMixture}(w)pauliMixture(w) is the function on all complex 2×22\times22×2 matrices given by

ρ⟼∑k=03wkPkρPk†,\rho\longmapsto\sum_{k=0}^3w_kP_k\rho P_k^\dagger,ρ⟼k=0∑3​wk​Pk​ρPk†​,

where (P0,P1,P2,P3)=(I2,X,Y,Z)(P_0,P_1,P_2,P_3)=(I_2,X,Y,Z)(P0​,P1​,P2​,P3​)=(I2​,X,Y,Z), with X,Y,ZX,Y,ZX,Y,Z the standard Pauli matrices in the Boolean basis ordered as (false,true)(\mathrm{false},\mathrm{true})(false,true), and †{}^\dagger† denotes conjugate transpose. No nonnegativity or normalization condition on www, and no density-matrix condition on ρ\rhoρ, is required.

QAlgorithms.Papers.Independent.Intakethird100.P13.pauliPointChannel

For any real three-component vector q=(q0,q1,q2)q=(q_0,q_1,q_2)q=(q0​,q1​,q2​), pauliPointChannel⁡(q)\operatorname{pauliPointChannel}(q)pauliPointChannel(q) is the function on all complex 2×22\times22×2 matrices defined by

ρ⟼Tr⁡ρ2I2+12(q0Tr⁡(Xρ)X−q1Tr⁡(Yρ)Y+q2Tr⁡(Zρ)Z),\rho\longmapsto \frac{\operatorname{Tr}\rho}{2}I_2+ \frac12\left( q_0\operatorname{Tr}(X\rho)X -q_1\operatorname{Tr}(Y\rho)Y +q_2\operatorname{Tr}(Z\rho)Z \right),ρ⟼2Trρ​I2​+21​(q0​Tr(Xρ)X−q1​Tr(Yρ)Y+q2​Tr(Zρ)Z),

where X,Y,ZX,Y,ZX,Y,Z are the standard Pauli matrices in the ordered Boolean basis and the traces are complex-valued. The vector qqq is arbitrary; membership in any specified region or a channel-positivity condition is not part of this definition.

QAlgorithms.Papers.Independent.Intakethird100.P13.pauliTetrahedron

The defined set pauliTetrahedron⁡\operatorname{pauliTetrahedron}pauliTetrahedron consists of all real three-component vectors expressible as

λ0(1,−1,1)+λ1(1,1,−1)+λ2(−1,−1,−1)+λ3(−1,1,1),λk≥0,∑k=03λk=1.\lambda_0(1,-1,1)+\lambda_1(1,1,-1) +\lambda_2(-1,-1,-1)+\lambda_3(-1,1,1), \qquad \lambda_k\ge0,\quad \sum_{k=0}^3\lambda_k=1.λ0​(1,−1,1)+λ1​(1,1,−1)+λ2​(−1,−1,−1)+λ3​(−1,1,1),λk​≥0,k=0∑3​λk​=1.

Thus it is the real convex hull of these four explicitly listed vectors, including their boundary points.

QAlgorithms.Papers.Independent.Intakethird100.P13.qubitBellBranch

For any complex matrix τ\tauτ indexed by pairs of Boolean values, any k∈{0,1,2,3}k\in\{0,1,2,3\}k∈{0,1,2,3}, and any complex matrix ρ\rhoρ indexed by Boolean values, this defines a Boolean-indexed complex matrix Cτ,k(ρ)C_{\tau,k}(\rho)Cτ,k​(ρ) whose entries, for Boolean b,b′b,b'b,b′, are

[Cτ,k(ρ)]b,b′=12∑c,a,d,a′∈{false,true}(Pk)a,c‾ ρc,d τ(a,b),(a′,b′) (Pk)a′,d,[C_{\tau,k}(\rho)]_{b,b'}= \frac12\sum_{c,a,d,a'\in\{\mathrm{false},\mathrm{true}\}} \overline{(P_k)_{a,c}}\, \rho_{c,d}\, \tau_{(a,b),(a',b')}\, (P_k)_{a',d},[Cτ,k​(ρ)]b,b′​=21​c,a,d,a′∈{false,true}∑​(Pk​)a,c​​ρc,d​τ(a,b),(a′,b′)​(Pk​)a′,d​,

where a bar denotes complex conjugation and, in Boolean order false,true\mathrm{false},\mathrm{true}false,true,

P0=(1001),P1=(0110),P2=(0−ii0),P3=(100−1).P_0=\begin{pmatrix}1&0\\0&1\end{pmatrix},\quad P_1=\begin{pmatrix}0&1\\1&0\end{pmatrix},\quad P_2=\begin{pmatrix}0&-\mathrm i\\ \mathrm i&0\end{pmatrix},\quad P_3=\begin{pmatrix}1&0\\0&-1\end{pmatrix}.P0​=(10​01​),P1​=(01​10​),P2​=(0i​−i0​),P3​=(10​0−1​).

Neither τ\tauτ nor ρ\rhoρ is required to be a density matrix.

QAlgorithms.Papers.Independent.Intakethird100.P13.regularPauliPoint

For all real numbers γ,s30\gamma,s_{30}γ,s30​ and all vectors s=(s0,s1,s2)∈R3s=(s_0,s_1,s_2)\in\mathbb R^3s=(s0​,s1​,s2​)∈R3, set r=(1−γ)/(1−γs30)r=(1-\gamma)/(1-\gamma s_{30})r=(1−γ)/(1−γs30​). The definition gives

regularPauliPoint⁡(γ,s30,s)=(r s0,r s1,rs2).\operatorname{regularPauliPoint}(\gamma,s_{30},s) =(\sqrt r\,s_0,\sqrt r\,s_1,rs_2).regularPauliPoint(γ,s30​,s)=(r​s0​,r​s1​,rs2​).

No bounds or nonzero-denominator hypothesis are imposed: division by zero has value 000, so the resulting vector is zero when 1−γs30=01-\gamma s_{30}=01−γs30​=0, and r=0\sqrt r=0r​=0 when r<0r<0r<0.

QAlgorithms.Papers.Independent.Intakethird100.P13.rescaling

For every pair of real numbers γ,s30\gamma,s_{30}γ,s30​, this defines

1−γ1−γs30.\frac{1-\gamma}{1-\gamma s_{30}}.1−γs30​1−γ​.

No nonzero-denominator or parameter-range hypothesis is imposed; when 1−γs30=01-\gamma s_{30}=01−γs30​=0, the value is zero by the convention for real division by zero.

QAlgorithms.Papers.Independent.Intakethird100.P13.resourcePauliCorrelation

For any complex matrix τ\tauτ indexed by {0,1}2\{0,1\}^2{0,1}2, the defined real 3×33\times33×3 matrix resourcePauliCorrelation⁡(τ)\operatorname{resourcePauliCorrelation}(\tau)resourcePauliCorrelation(τ) has row jjj, column iii entry

[resourcePauliCorrelation⁡(τ)]ji=Re⁡Tr⁡ ⁣[(σj⊗σi)τ],j,i∈{0,1,2},[\operatorname{resourcePauliCorrelation}(\tau)]_{ji} =\operatorname{Re}\operatorname{Tr}\!\left[(\sigma_j\otimes\sigma_i)\tau\right], \qquad j,i\in\{0,1,2\},[resourcePauliCorrelation(τ)]ji​=ReTr[(σj​⊗σi​)τ],j,i∈{0,1,2},

where (σ0,σ1,σ2)=(X,Y,Z)(\sigma_0,\sigma_1,\sigma_2)=(X,Y,Z)(σ0​,σ1​,σ2​)=(X,Y,Z) are the standard Pauli matrices in the ordered Boolean basis. No density-matrix assumption on τ\tauτ is included.

QAlgorithms.Papers.Independent.Intakethird100.P13.sCoordinate

For any real array p:{0,1,2,3}2→Rp:\{0,1,2,3\}^2\to\mathbb Rp:{0,1,2,3}2→R, any i∈{0,1,2}i\in\{0,1,2\}i∈{0,1,2}, and any j∈{0,1,2,3}j\in\{0,1,2,3\}j∈{0,1,2,3}, the defined real coordinate is

sij(p)=14∑k=03∑l=03(−1)δk0+δj2+δj0+δkj+δi+1,l+δ0l plk,s_{ij}(p)=\frac14\sum_{k=0}^3\sum_{l=0}^3 (-1)^{\delta_{k0}+\delta_{j2}+\delta_{j0}+\delta_{kj} +\delta_{i+1,l}+\delta_{0l}}\,p_{lk},sij​(p)=41​k=0∑3​l=0∑3​(−1)δk0​+δj2​+δj0​+δkj​+δi+1,l​+δ0l​plk​,

where δab\delta_{ab}δab​ is 111 when a=ba=ba=b and 000 otherwise, and the exponent is a sum of natural numbers. The array ppp is arbitrary; nonnegativity and column-sum conditions are not required.

QAlgorithms.Papers.Independent.Intakethird100.P13.teleportationBellAmplitude

For every k∈{0,1,2,3}k\in\{0,1,2,3\}k∈{0,1,2,3} and Boolean values c,ac,ac,a, this defines the complex number

βk(c,a)=12(Pk)a,c,\beta_k(c,a)=\frac{1}{\sqrt2}(P_k)_{a,c},βk​(c,a)=2​1​(Pk​)a,c​,

where, in Boolean order false,true\mathrm{false},\mathrm{true}false,true,

P0=(1001),P1=(0110),P2=(0−ii0),P3=(100−1).P_0=\begin{pmatrix}1&0\\0&1\end{pmatrix},\quad P_1=\begin{pmatrix}0&1\\1&0\end{pmatrix},\quad P_2=\begin{pmatrix}0&-\mathrm i\\ \mathrm i&0\end{pmatrix},\quad P_3=\begin{pmatrix}1&0\\0&-1\end{pmatrix}.P0​=(10​01​),P1​=(01​10​),P2​=(0i​−i0​),P3​=(10​0−1​).

Thus aaa is the row index and ccc is the column index of the selected matrix.

QAlgorithms.Papers.Independent.Intakethird100.P13.teleportationBellMixture

For any real-valued function www on {0,1,2,3}\{0,1,2,3\}{0,1,2,3}, this defines the complex matrix TwT_wTw​ indexed by pairs of Boolean values through

(Tw)(c,a),(d,a′)=12∑k=03wk(Pk)a,c (Pk)a′,d‾,(T_w)_{(c,a),(d,a')}= \frac12\sum_{k=0}^{3}w_k(P_k)_{a,c}\,\overline{(P_k)_{a',d}},(Tw​)(c,a),(d,a′)​=21​k=0∑3​wk​(Pk​)a,c​(Pk​)a′,d​​,

for all Boolean c,a,d,a′c,a,d,a'c,a,d,a′, where complex conjugation is denoted by a bar and, in Boolean order false,true\mathrm{false},\mathrm{true}false,true,

P0=(1001),P1=(0110),P2=(0−ii0),P3=(100−1).P_0=\begin{pmatrix}1&0\\0&1\end{pmatrix},\quad P_1=\begin{pmatrix}0&1\\1&0\end{pmatrix},\quad P_2=\begin{pmatrix}0&-\mathrm i\\ \mathrm i&0\end{pmatrix},\quad P_3=\begin{pmatrix}1&0\\0&-1\end{pmatrix}.P0​=(10​01​),P1​=(01​10​),P2​=(0i​−i0​),P3​=(10​0−1​).

The weights are arbitrary real numbers: nonnegativity and a sum of one are not assumed.

QAlgorithms.Papers.Independent.Intakethird100.P13.teleportationBellProjector

For every k∈{0,1,2,3}k\in\{0,1,2,3\}k∈{0,1,2,3}, this defines a complex matrix BkB_kBk​ indexed by pairs of Boolean values: for x=(c,a)x=(c,a)x=(c,a) and y=(d,a′)y=(d,a')y=(d,a′),

(Bk)(c,a),(d,a′)=12(Pk)a,c (Pk)a′,d‾,(B_k)_{(c,a),(d,a')}=\frac12(P_k)_{a,c}\,\overline{(P_k)_{a',d}},(Bk​)(c,a),(d,a′)​=21​(Pk​)a,c​(Pk​)a′,d​​,

where the bar denotes complex conjugation and, in Boolean order false,true\mathrm{false},\mathrm{true}false,true,

P0=(1001),P1=(0110),P2=(0−ii0),P3=(100−1).P_0=\begin{pmatrix}1&0\\0&1\end{pmatrix},\quad P_1=\begin{pmatrix}0&1\\1&0\end{pmatrix},\quad P_2=\begin{pmatrix}0&-\mathrm i\\ \mathrm i&0\end{pmatrix},\quad P_3=\begin{pmatrix}1&0\\0&-1\end{pmatrix}.P0​=(10​01​),P1​=(01​10​),P2​=(0i​−i0​),P3​=(10​0−1​).

QAlgorithms.Papers.Independent.Intakethird100.P13.teleportationPauli

This defines four complex 2×22\times22×2 matrices PkP_kPk​, indexed by k∈{0,1,2,3}k\in\{0,1,2,3\}k∈{0,1,2,3}. With rows and columns ordered as false,true\mathrm{false},\mathrm{true}false,true, they are

P0=(1001),P1=(0110),P2=(0−ii0),P3=(100−1),P_0=\begin{pmatrix}1&0\\0&1\end{pmatrix},\qquad P_1=\begin{pmatrix}0&1\\1&0\end{pmatrix},\qquad P_2=\begin{pmatrix}0&-\mathrm i\\ \mathrm i&0\end{pmatrix},\qquad P_3=\begin{pmatrix}1&0\\0&-1\end{pmatrix},P0​=(10​01​),P1​=(01​10​),P2​=(0i​−i0​),P3​=(10​0−1​),

where i2=−1\mathrm i^2=-1i2=−1.

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