Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Approximate recovery with locality and symmetry constraints

Definition
QAlgorithms_Papers_Independent_Intakethird100_P44

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

prove2quantum

The objects the formalized statements of Approximate recovery with locality and symmetry constraints (arXiv:1806.10324) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • ObservableAlgebra: PDF p.3: unital dagger algebra with its canonical Hilbert--Schmidt orthogonal projection channel. Range and orthogonality uniquely specify it.
  • IsPhysicalChannel: PDF p.3, Definition 1; no condition is imposed on forbidden output coherences.
  • FixesObservables: PDF pp.7--8: fixing an observable algebra is a Heisenberg condition.
  • fixedFidelity: Root entanglement fidelity on the specified normalized input, retaining the identity reference. PDF p.4, Eq.(5).
  • Comparison: The two comparison domains in Eqs.(5)--(6), PDF pp.3--4. An encoding is an isometry; no surjectivity or fixed-charge assumption occurs.
  • comparisonRootFidelity: PDF p.4, Eq.(6): infimum over ALL normalized logical densities. For channels the value set is nonempty and contained in [0,1]; compactness makes this infimum a minimum. The fixed branch uses ROOT fidelity as well.
  • SameMaximum: Equality of attained maxima, avoiding totalized real suprema on empty sets. This includes the existence of both optimizers stated on PDF pp.5,8--9.
  • IsChannelDilation: A true Stinespring dilation of N, on every matrix. PDF pp.4,8.
  • krausComplement: Ordinary complementary channel of a normalized Kraus family. The entry (i,j) is Tr(E_j† E_i rho), PDF p.4, Eq.(9).
  • IsGeneralizedComplement: PDF p.4: complements include channels reversibly postprocessable to a Stinespring complement, with BOTH postprocessings CPTP.
  • localComplement: PDF p.8, Definition 3, Eq.(33): the local environment retains B TOGETHER with the dilation environment, rather than the environment alone.
  • Sectors: PDF p.5, Section III B: nonzero orthogonal sector projectors, with pairwise distinct real charge eigenvalues.
  • sectorProject: PDF p.5, Eq.(14), the canonical blind charge measurement.
  • sectorChargeObservable: PDF p.5: the charge is computed from the sector data.
  • PhysicallyCorrectable: Exact PHYSICAL correction to P on the encoded space, PDF pp.3,5. Equality on all logical matrices is equivalent by linearity to equality on all logical density matrices; reference extensions are thereby retained.
  • ChargeStateCondition: PDF p.5, Eqs.(18),(20): each conditional output is normalized only AFTER summing its classical output sectors. Individual blocks need not have trace 1.
  • ChargeCoefficients: PDF p.5, Eq.(21): one coefficient family, uniform in all subsequent indices.
  • constantChannel: PDF p.9, Eq.(48): constant-output channel, on ALL matrices.
  • TensorCorrectable: PDF pp.9--10: tensor-local exact correction on an arbitrary JOINT code.
  • codeSlice: PDF p.10, Eq.(51): W_n = (I tensor <n|) W.
  • Fermions: PDF pp.2--3: a finite spinless-fermion CAR representation on its actual 2^n-dimensional Fock space. Minimal dimension prevents extra spectator sectors.
  • totalParity: PDF pp.2--3: total parity, up to the immaterial global sign convention. i^n times the ordered product of 2n Majoranas is a Hermitian involution. Printed, PDF p.2 Eq.(1): i^(|omega|-1) for Majorana labels. Changed: i^n on 2n ordered modes. Why: the printed four-mode value is anti-Hermitian and squares to -I; parity must be a Hermitian involution.
  • parityProject: PDF p.3: blind parity measurement.
  • DefiniteParity: PDF p.3, Proposition 1: both preserving and flipping parity are allowed.
  • evenRegion: PDF pp.2--3,10: the even algebra supported on a set of Majorana indices. Repeated indices are harmless under CAR; the empty word supplies the identity.
  • StrongLocal: PDF p.7, Definition 2: strong locality retains BOTH clauses.
  • siteModes: The two Majoranas for each fermionic SITE in omega, PDF pp.2--3,10.
  • regionParity: PDF p.10: regional parity for site regions, including the empty region. Printed, PDF p.2 Eq.(1): i^(|omega|-1) for Majorana labels. Changed: i^m on the 2m ordered modes of m sites, up to overall sign. Why: the printed four-mode value is anti-Hermitian and squares to -I, contrary to the source's parity involution and sector-projector conventions.
  • paritySector: PDF pp.3,10: plus sector indexed by false, minus by true.
  • StrongCorrectable: PDF pp.9--10: exact strong-local correction of the algebra A vee A'. For the even algebra of a SITE region this target is regional parity dephasing.
  • RegionalStateCondition: PDF p.10, Eq.(53); one normalized environmental state per regional sector. The order X_ij E_i† B E_j is copied from the rendered page.
  • RegionalCoefficients: PDF p.10, Eq.(54): the coefficient family is UNIFORM in the commutant observable Y. Sector compression is on both sides, without division by its weight.
  • ringGap: PDF pp.6,10: number of paired Majoranas strictly between adjacent free modes.
  • RingLayout: PDF pp.6--7,10--11: ordered free modes and the alternating nearest-neighbor paired intervals between them. Even gaps give the alternating pairing types.
  • modeAt: Cyclic mode at a specified offset.
  • pairingProjector: PDF pp.6--7,10--11: ACTUAL pairing-ground projector. Each paired nearest-neighbor fermion is initialized in its vacuum. The ordered bond (a,b) has vacuum projector (I-i w_a w_b)/2; reversing bond orientation changes its sign. The free modes receive no projector, so ALL their states are encoded.
  • RingCode: PDF pp.6--7: an isometry ONTO the computed pairing ground code, not an arbitrary invariant subspace. Logical dimension is exactly 2^k.
  • minimumGap: PDF p.10: shortest paired interval. Initial bound 2n dominates every gap; the nonempty list has length 2k. Empty paired intervals therefore give d=0.
  • HasIntervalSupport: PDF p.10: a WHOLE Kraus operator is an even polynomial supported in one cyclic interval containing length modes; length=0 permits only scalar operators.
  • poisoning: PDF pp.6--7: equally weighted jumps on the first two unpaired Majoranas.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_QAlgorithms_Defs_QuantumInformation
import Definitions.Def_QAlgorithms_Defs_Measurement
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P60
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P67
import Definitions.Def_QAlgorithms_Papers_Independent_Claimed50_P48
import Definitions.Def_QAlgorithms_Papers_Independent_Renew46_P43

/-!
Statements from Bény, Zimborás and Pastawski, arXiv:1806.10324v2.
All spaces are finite matrix spaces. Nonempty ambient and logical index types
are explicit theorem assumptions. New interfaces are paper-local.
The independent-paper override permits imports from the pinned accepted library;
exactly matching matrix, tensor, adjoint and purification interfaces are reused.
The four source-ambiguity nodes are recorded in HARD.md, not asserted here.
-/

noncomputable section

namespace QAlgorithms.Papers.Independent.Intakethird100.P44

open scoped Matrix ComplexConjugate ComplexOrder

open QAlgorithms.Papers.Independent.Information100.P60 (Mat)
open QAlgorithms.Papers.Independent.Information100.P67 (mapFirst)
open QAlgorithms.Papers.Independent.Claimed50.P48 (operatorTraceDual)
open QAlgorithms.Papers.Independent.Renew46.P43 (purification)

variable {H K L J E F B : Type*}
variable [Fintype H] [DecidableEq H] [Fintype K] [DecidableEq K]
variable [Fintype L] [DecidableEq L] [Fintype J] [DecidableEq J]
variable [Fintype E] [DecidableEq E] [Fintype F] [DecidableEq F]
variable [Fintype B] [DecidableEq B]

/-- PDF p.3: unital dagger algebra with its canonical Hilbert--Schmidt
orthogonal projection channel. Range and orthogonality uniquely specify it. -/
structure ObservableAlgebra (H : Type*) [Fintype H] [DecidableEq H] where
  algebra : StarSubalgebra ℂ (Mat H)
  project : Superop H H
  channel : IsChannel project
  range_mem : ∀ X, project X ∈ algebra
  orthogonal : ∀ X Y, Y ∈ algebra → ((X - project X)ᴴ * Y).trace = 0

/-- PDF p.3, Definition 1; no condition is imposed on forbidden output coherences. -/
def IsPhysicalChannel (P : Superop H H) (Q : Superop K K) (N : Superop H K) : Prop :=
  ∀ X, Q (N (P X)) = Q (N X)

/-- PDF pp.7--8: fixing an observable algebra is a Heisenberg condition. -/
def FixesObservables (A : Set (Mat H)) (N : Superop H H) : Prop :=
  ∀ X ∈ A, operatorTraceDual N X = X

/-- Root entanglement fidelity on the specified normalized input, retaining
the identity reference. PDF p.4, Eq.(5). -/
def fixedFidelity (rho : Mat H) (N M : Superop H K) : ℝ :=
  rootFidelity (idTensorSuperop N (purification rho))
    (idTensorSuperop M (purification rho))

/-- The two comparison domains in Eqs.(5)--(6), PDF pp.3--4.
An encoding is an isometry; no surjectivity or fixed-charge assumption occurs. -/
inductive Comparison (H : Type*) [Fintype H] [DecidableEq H] where
  | fixed (rho : Mat H) (normalized : IsDensityMatrix rho)
  | code (k : ℕ) (positive : 0 < k) (W : Matrix H (Fin k) ℂ)
      (isometry : Wᴴ * W = 1)

/-- PDF p.4, Eq.(6): infimum over ALL normalized logical densities.
For channels the value set is nonempty and contained in [0,1]; compactness
makes this infimum a minimum. The fixed branch uses ROOT fidelity as well. -/
def comparisonRootFidelity (c : Comparison H) (N M : Superop H K) : ℝ :=
  match c with
  | .fixed rho _ => fixedFidelity rho N M
  | .code k _ W _ =>
      sInf {t : ℝ | ∃ rho : Mat (Fin k), IsDensityMatrix rho ∧
        t = fixedFidelity rho (N ∘ conjSuperop W) (M ∘ conjSuperop W)}

/-- Equality of attained maxima, avoiding totalized real suprema on empty sets.
This includes the existence of both optimizers stated on PDF pp.5,8--9. -/
def SameMaximum {X Y : Type*} (allowedX : X → Prop) (f : X → ℝ)
    (allowedY : Y → Prop) (g : Y → ℝ) : Prop :=
  ∃ x y, allowedX x ∧ allowedY y ∧ f x = g y ∧
    (∀ x2, allowedX x2 → f x2 ≤ f x) ∧
    (∀ y2, allowedY y2 → g y2 ≤ g y)

/-- A true Stinespring dilation of N, on every matrix. PDF pp.4,8. -/
def IsChannelDilation (N : Superop H K) (V : Matrix (K × E) H ℂ) : Prop :=
  Vᴴ * V = 1 ∧ ∀ X, N X = traceRight (conjSuperop V X)

/-- Ordinary complementary channel of a normalized Kraus family. The entry
(i,j) is Tr(E_j† E_i rho), PDF p.4, Eq.(9). -/
def krausComplement (A : J → Matrix K H ℂ) : Superop H J :=
  fun X i j => (A i * X * (A j)ᴴ).trace

/-- PDF p.4: complements include channels reversibly postprocessable to a
Stinespring complement, with BOTH postprocessings CPTP. -/
def IsGeneralizedComplement (N : Superop H K) (C : Superop H E) : Prop :=
  ∃ (m : ℕ) (A : Fin m → Matrix K H ℂ), IsKrausFamily A ∧ N = krausApply A ∧
    ∃ (R : Superop (Fin m) E) (S : Superop E (Fin m)),
      IsChannel R ∧ IsChannel S ∧ C = R ∘ krausComplement A ∧
      krausComplement A = S ∘ C

/-- PDF p.8, Definition 3, Eq.(33): the local environment retains B TOGETHER
with the dilation environment, rather than the environment alone. -/
def localComplement (A : ObservableAlgebra H) (V : Matrix (H × E) H ℂ) :
    Superop H (H × E) :=
  mapFirst A.project ∘ conjSuperop V

/-- PDF p.5, Section III B: nonzero orthogonal sector projectors, with
pairwise distinct real charge eigenvalues. -/
structure Sectors (H : Type*) [Fintype H] [DecidableEq H] (r : ℕ) where
  P : Fin r → Mat H
  hermitian : ∀ i, (P i)ᴴ = P i
  nonzero : ∀ i, P i ≠ 0
  orthogonal : ∀ i j, P i * P j = if i = j then P i else 0
  complete : ∑ i, P i = 1
  values : Fin r → ℝ
  distinct : Function.Injective values

/-- PDF p.5, Eq.(14), the canonical blind charge measurement. -/
def sectorProject {r : ℕ} (s : Sectors H r) : Superop H H :=
  fun X => ∑ i, s.P i * X * s.P i

/-- PDF p.5: the charge is computed from the sector data. -/
def sectorChargeObservable {r : ℕ} (s : Sectors H r) : Mat H :=
  ∑ i, (s.values i : ℂ) • s.P i

/-- Exact PHYSICAL correction to P on the encoded space, PDF pp.3,5.
Equality on all logical matrices is equivalent by linearity to equality on
all logical density matrices; reference extensions are thereby retained. -/
def PhysicallyCorrectable (P : Superop H H) (N : Superop H H) (W : Matrix H L ℂ) : Prop :=
  ∃ R : Superop H H, IsChannel R ∧ IsPhysicalChannel P P R ∧
    ∀ X : Mat L, R (N (conjSuperop W X)) = P (conjSuperop W X)

/-- PDF p.5, Eqs.(18),(20): each conditional output is normalized only
AFTER summing its classical output sectors. Individual blocks need not have trace 1. -/
def ChargeStateCondition {r : ℕ} (s : Sectors H r) (A : J → Mat H)
    (W : Matrix H L ℂ) : Prop :=
  ∃ sigma : Fin r → Fin r → Mat J,
    (∀ i j, (sigma i j).PosSemidef) ∧
    (∀ i, ∑ j, (sigma i j).trace = 1) ∧
    ∀ j n m, Wᴴ * (A n)ᴴ * s.P j * A m * W =
      ∑ i, (sigma i j n m) • (Wᴴ * s.P i * W)

/-- PDF p.5, Eq.(21): one coefficient family, uniform in all subsequent indices. -/
def ChargeCoefficients {r : ℕ} (s : Sectors H r) (A : J → Mat H)
    (W : Matrix H L ℂ) : Prop :=
  ∃ c : Fin r → Fin r → J → J → ℂ,
    ∀ j n m, Wᴴ * (A n)ᴴ * s.P j * A m * W =
      ∑ i, c i j n m • (Wᴴ * s.P i * W)

/-- PDF p.9, Eq.(48): constant-output channel, on ALL matrices. -/
def constantChannel (sigma : Mat K) : Superop H K :=
  fun X => X.trace • sigma

/-- PDF pp.9--10: tensor-local exact correction on an arbitrary JOINT code. -/
def TensorCorrectable (N : Superop H H) (W : Matrix (H × B) L ℂ) : Prop :=
  ∃ R : Superop H H, IsChannel R ∧ ∀ X : Mat L,
    mapFirst (R ∘ N) (conjSuperop W X) = conjSuperop W X

/-- PDF p.10, Eq.(51): W_n = (I tensor <n|) W. -/
def codeSlice (W : Matrix (H × B) L ℂ) (n : B) : Matrix H L ℂ :=
  fun a l => W (a, n) l

/-- PDF pp.2--3: a finite spinless-fermion CAR representation on its actual
2^n-dimensional Fock space. Minimal dimension prevents extra spectator sectors. -/
structure Fermions (n : ℕ) where
  w : Fin (2 * n) → Mat (Fin (2 ^ n))
  hermitian : ∀ j, (w j)ᴴ = w j
  car : ∀ i j, w i * w j + w j * w i =
    if i = j then (2 : ℂ) • 1 else 0

/-- PDF pp.2--3: total parity, up to the immaterial global sign convention.
i^n times the ordered product of 2n Majoranas is a Hermitian involution.
Printed, PDF p.2 Eq.(1): i^(|omega|-1) for Majorana labels.
Changed: i^n on 2n ordered modes. Why: the printed four-mode value is
anti-Hermitian and squares to -I; parity must be a Hermitian involution. -/
def totalParity {n : ℕ} (f : Fermions n) : Mat (Fin (2 ^ n)) :=
  (Complex.I ^ n) • (List.ofFn f.w).prod

/-- PDF p.3: blind parity measurement. -/
def parityProject (C : Mat H) : Superop H H :=
  fun X => (1 / 2 : ℂ) • (X + C * X * C)

/-- PDF p.3, Proposition 1: both preserving and flipping parity are allowed. -/
def DefiniteParity (C X : Mat H) : Prop :=
  X * C = C * X ∨ X * C = -(C * X)

/-- PDF pp.2--3,10: the even algebra supported on a set of Majorana indices.
Repeated indices are harmless under CAR; the empty word supplies the identity. -/
def evenRegion {n : ℕ} (f : Fermions n) (S : Finset (Fin (2 * n))) :
    Submodule ℂ (Mat (Fin (2 ^ n))) :=
  Submodule.span ℂ {X | ∃ word : List (Fin (2 * n)),
    (∀ j ∈ word, j ∈ S) ∧ Even word.length ∧ X = (word.map f.w).prod}

/-- PDF p.7, Definition 2: strong locality retains BOTH clauses. -/
def StrongLocal (A : Set (Mat H)) (N : Superop H H) : Prop :=
  IsChannel N ∧ (∀ X ∈ A, operatorTraceDual N X ∈ A) ∧
    FixesObservables (Set.centralizer A) N

/-- The two Majoranas for each fermionic SITE in omega, PDF pp.2--3,10. -/
def siteModes {n : ℕ} (omega : Finset (Fin n)) : Finset (Fin (2 * n)) :=
  Finset.univ.filter fun j => ∃ s ∈ omega, j.val / 2 = s.val

/-- PDF p.10: regional parity for site regions, including the empty region.
Printed, PDF p.2 Eq.(1): i^(|omega|-1) for Majorana labels.
Changed: i^m on the 2m ordered modes of m sites, up to overall sign.
Why: the printed four-mode value is anti-Hermitian and squares to -I,
contrary to the source's parity involution and sector-projector conventions. -/
def regionParity {n : ℕ} (f : Fermions n) (omega : Finset (Fin n)) :
    Mat (Fin (2 ^ n)) :=
  (Complex.I ^ omega.card) •
    (((List.finRange (2 * n)).filter fun j => j ∈ siteModes omega).map f.w).prod

/-- PDF pp.3,10: plus sector indexed by false, minus by true. -/
def paritySector (C : Mat H) (k : Bool) : Mat H :=
  (1 / 2 : ℂ) • (1 + (if k then -C else C))

/-- PDF pp.9--10: exact strong-local correction of the algebra A vee A'.
For the even algebra of a SITE region this target is regional parity dephasing. -/
def StrongCorrectable (A : Set (Mat H)) (C : Mat H)
    (N : Superop H H) (W : Matrix H L ℂ) : Prop :=
  ∃ R : Superop H H, StrongLocal A R ∧ ∀ X : Mat L,
    R (N (conjSuperop W X)) = parityProject C (conjSuperop W X)

/-- PDF p.10, Eq.(53); one normalized environmental state per regional sector.
The order X_ij E_i† B E_j is copied from the rendered page. -/
def RegionalStateCondition (A : Set (Mat H)) (C : Mat H)
    (noise : J → Mat H) (W : Matrix H L ℂ) : Prop :=
  ∃ rho : Bool → Mat J, (∀ k, IsDensityMatrix (rho k)) ∧
    ∀ (Y : Mat H), Y ∈ Set.centralizer A → ∀ X : Mat J,
      (∑ i, ∑ j, X i j • (Wᴴ * (noise i)ᴴ * Y * noise j * W)) =
      ∑ k : Bool, (rho k * X).trace •
        (Wᴴ * paritySector C k * Y * paritySector C k * W)

/-- PDF p.10, Eq.(54): the coefficient family is UNIFORM in the commutant
observable Y. Sector compression is on both sides, without division by its weight. -/
def RegionalCoefficients (A : Set (Mat H)) (C : Mat H)
    (noise : J → Mat H) (W : Matrix H L ℂ) : Prop :=
  ∃ lambda : J → J → Bool → ℂ, ∀ i j k (Y : Mat H), Y ∈ Set.centralizer A →
    Wᴴ * (noise i)ᴴ * noise j * (paritySector C k * Y * paritySector C k) * W =
      lambda i j k • (Wᴴ * (paritySector C k * Y * paritySector C k) * W)

/-- PDF pp.6,10: number of paired Majoranas strictly between adjacent free modes. -/
def ringGap {n k : ℕ} (free : Fin (2 * k) → Fin (2 * n))
    (j : Fin (2 * k)) : ℕ :=
  ((free (finRotate (2 * k) j)).val + 2 * n - (free j).val - 1) % (2 * n)

/-- PDF pp.6--7,10--11: ordered free modes and the alternating nearest-neighbor
paired intervals between them. Even gaps give the alternating pairing types. -/
structure RingLayout (n k : ℕ) where
  sites_pos : 0 < n
  logical_pos : 0 < k
  free : Fin (2 * k) → Fin (2 * n)
  ordered : StrictMono free
  even_gaps : ∀ j, Even (ringGap free j)

/-- Cyclic mode at a specified offset. -/
def modeAt {n : ℕ} (hn : 0 < n) (j : Fin (2 * n)) (offset : ℕ) : Fin (2 * n) :=
  ⟨(j.val + offset) % (2 * n), Nat.mod_lt _ (by omega)⟩

/-- PDF pp.6--7,10--11: ACTUAL pairing-ground projector. Each paired
nearest-neighbor fermion is initialized in its vacuum. The ordered bond (a,b)
has vacuum projector (I-i w_a w_b)/2; reversing bond orientation changes its sign.
The free modes receive no projector, so ALL their states are encoded. -/
def pairingProjector {n k : ℕ} (f : Fermions n) (r : RingLayout n k) :
    Mat (Fin (2 ^ n)) :=
  (((List.finRange (2 * k)).flatMap fun j =>
    (List.range (ringGap r.free j / 2)).map fun t =>
      (1 / 2 : ℂ) • (1 - Complex.I •
        (f.w (modeAt r.sites_pos (r.free j) (1 + 2 * t)) *
         f.w (modeAt r.sites_pos (r.free j) (2 + 2 * t)))))).prod

/-- PDF pp.6--7: an isometry ONTO the computed pairing ground code, not an
arbitrary invariant subspace. Logical dimension is exactly 2^k. -/
structure RingCode {n k : ℕ} (f : Fermions n) (r : RingLayout n k) where
  W : Matrix (Fin (2 ^ n)) (Fin (2 ^ k)) ℂ
  isometry : Wᴴ * W = 1
  onto_pairing_ground : W * Wᴴ = pairingProjector f r

/-- PDF p.10: shortest paired interval. Initial bound 2n dominates every gap;
the nonempty list has length 2k. Empty paired intervals therefore give d=0. -/
def minimumGap {n k : ℕ} (r : RingLayout n k) : ℕ :=
  (List.ofFn (ringGap r.free)).foldl min (2 * n)

/-- PDF p.10: a WHOLE Kraus operator is an even polynomial supported in one
cyclic interval containing length modes; length=0 permits only scalar operators. -/
def HasIntervalSupport {n : ℕ} (f : Fermions n)
    (X : Mat (Fin (2 ^ n))) (length : ℕ) : Prop :=
  ∃ start : Fin (2 * n), X ∈ evenRegion f
    (Finset.univ.filter fun j =>
      (j.val + 2 * n - start.val) % (2 * n) < length)

/-- PDF pp.6--7: equally weighted jumps on the first two unpaired Majoranas. -/
def poisoning {n k : ℕ} (f : Fermions n) (r : RingLayout n k) :
    Superop (Fin (2 ^ n)) (Fin (2 ^ n)) :=
  let a : Fin (2 * k) := ⟨0, by have := r.logical_pos; omega⟩
  let b : Fin (2 * k) := ⟨1, by have := r.logical_pos; omega⟩
  fun X => (1 / 2 : ℂ) •
    (conjSuperop (f.w (r.free a)) X + conjSuperop (f.w (r.free b)) X)

/-- Sanity check for the reused trace adjoint on matrix units. -/
example (N : Superop H K) (X : Mat K) (i j : H) :
    operatorTraceDual N X i j = (N (Matrix.single j i 1) * X).trace := by
  rfl

example {n k : ℕ} (f : Fermions n) (r : RingLayout n k)
    (hz : ∀ j, ringGap r.free j = 0) : pairingProjector f r = 1 := by
  simp only [pairingProjector, hz, Nat.zero_div, List.range_zero, List.map_nil]
  have hempty : ∀ js : List (Fin (2 * k)),
      js.flatMap (fun _ => ([] : List (Mat (Fin (2 ^ n))))) = [] := by
    intro js
    induction js with
    | nil => rfl
    | cons j js ih => simpa only [List.flatMap_cons, List.nil_append] using ih
  rw [hempty]
  rfl

example (N : Superop H K) : IsPhysicalChannel id id N := by
  intro X
  rfl

example (N : Superop H K) (X : Mat (H × B)) (i j : K) (b b2 : B) :
    mapFirst N X (i, b) (j, b2) = N (fun a a2 => X (a, b) (a2, b2)) i j := by
  rfl

example (sigma : Mat K) : constantChannel (H := H) sigma (0 : Mat H) = 0 := by
  simp [constantChannel]

example (W : Matrix (H × B) L ℂ) (a : H) (b : B) (l : L) :
    codeSlice W b a l = W (a, b) l := by
  rfl

example (f : Fermions 0) : totalParity f = 1 := by
  simp [totalParity]

example {n : ℕ} (f : Fermions n) : regionParity f ∅ = 1 := by
  simp [regionParity, siteModes]

example : paritySector (1 : Mat H) true = 0 := by
  simp [paritySector]

example (C : Mat H) : DefiniteParity C 1 := by
  exact Or.inl (by simp)

end QAlgorithms.Papers.Independent.Intakethird100.P44
Source
Approximate recovery with locality and symmetry constraints, arXiv:1806.10324v2, https://arxiv.org/abs/1806.10324v2
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P44.ChargeCoefficients

Let HHH be a finite index type with decidable equality, J,LJ,LJ,L index types, and rrr a natural number. Let Pi∈Mat⁡(H)P_i\in\operatorname{Mat}(H)Pi​∈Mat(H), 0≤i<r0\le i<r0≤i<r, be nonzero Hermitian projections, mutually orthogonal and summing to III, with pairwise distinct real values viv_ivi​. Given matrices An∈Mat⁡(H)A_n\in\operatorname{Mat}(H)An​∈Mat(H), indexed by n∈Jn\in Jn∈J, and a complex H×LH\times LH×L matrix WWW, the coefficient condition means that there exist complex numbers cijnmc_{ijnm}cijnm​, indexed by i,j<ri,j<ri,j<r and n,m∈Jn,m\in Jn,m∈J, such that

W†An†PjAmW=∑i=0r−1cijnm W†PiWfor every j<r, n,m∈J.W^\dagger A_n^\dagger P_jA_mW =\sum_{i=0}^{r-1}c_{ijnm}\,W^\dagger P_iW \quad\text{for every }j<r,\ n,m\in J.W†An†​Pj​Am​W=i=0∑r−1​cijnm​W†Pi​Wfor every j<r, n,m∈J.

The coefficients have no positivity, Hermitian symmetry, or normalization requirements. No isometry condition on WWW or Kraus normalization of the AnA_nAn​ is imposed. If JJJ is empty, the equalities are vacuous; r=0r=0r=0 is also permitted when HHH is empty.

QAlgorithms.Papers.Independent.Intakethird100.P44.ChargeStateCondition

Let HHH be a finite index type with decidable equality, JJJ a finite index type, LLL an index type, and rrr a natural number. Let Pi∈Mat⁡(H)P_i\in\operatorname{Mat}(H)Pi​∈Mat(H), 0≤i<r0\le i<r0≤i<r, be nonzero Hermitian projections with PiPj=0P_iP_j=0Pi​Pj​=0 for i≠ji\ne ji=j and ∑iPi=I\sum_iP_i=I∑i​Pi​=I, carrying pairwise distinct real values viv_ivi​. Given matrices An∈Mat⁡(H)A_n\in\operatorname{Mat}(H)An​∈Mat(H), indexed by n∈Jn\in Jn∈J, and a complex H×LH\times LH×L matrix WWW, the condition means that there exist matrices σij∈Mat⁡(J)\sigma_{ij}\in\operatorname{Mat}(J)σij​∈Mat(J), for every i,j<ri,j<ri,j<r, such that each σij\sigma_{ij}σij​ is positive semidefinite, ∑jtr⁡(σij)=1\sum_j\operatorname{tr}(\sigma_{ij})=1∑j​tr(σij​)=1 for each iii, and

W†An†PjAmW=∑i=0r−1(σij)nm W†PiWfor every j<r, n,m∈J.W^\dagger A_n^\dagger P_jA_mW =\sum_{i=0}^{r-1}(\sigma_{ij})_{nm}\,W^\dagger P_iW \quad\text{for every }j<r,\ n,m\in J.W†An†​Pj​Am​W=i=0∑r−1​(σij​)nm​W†Pi​Wfor every j<r, n,m∈J.

No isometry condition on WWW or Kraus normalization of AnA_nAn​ is imposed by this definition. If JJJ is empty and r>0r>0r>0, the trace normalization is impossible; if r=0r=0r=0, all the displayed requirements are vacuous, with the sector data forcing HHH to be empty.

QAlgorithms.Papers.Independent.Intakethird100.P44.Comparison

For a finite type HHH with decidable equality, a comparison is data of exactly one of two forms: either a square complex matrix ρ\rhoρ indexed by HHH, together with proofs that ρ\rhoρ is positive semidefinite and Tr⁡ρ=1\operatorname{Tr}\rho=1Trρ=1; or a natural number k>0k>0k>0, together with a complex matrix WWW having rows indexed by HHH and columns indexed by {0,…,k−1}\{0,\ldots,k-1\}{0,…,k−1}, satisfying W†W=IkW^\dagger W=I_kW†W=Ik​. Thus the second form specifies an isometric embedding of a positive-dimensional space. If HHH is empty, neither form can be supplied; k=0k=0k=0 is excluded.

QAlgorithms.Papers.Independent.Intakethird100.P44.DefiniteParity

For complex square matrices C,XC,XC,X on an index type HHH, this predicate means

XC=CXorXC=−CX.XC=CX\quad\text{or}\quad XC=-CX.XC=CXorXC=−CX.

The disjunction is inclusive: both equations may hold, as happens when X=0X=0X=0 or C=0C=0C=0. No normalization, Hermiticity, or invertibility condition is imposed on either matrix.

QAlgorithms.Papers.Independent.Intakethird100.P44.Fermions

For every natural number nnn, this structure consists of 2n2n2n complex 2n×2n2^n\times2^n2n×2n matrices w0,…,w2n−1w_0,\ldots,w_{2n-1}w0​,…,w2n−1​, together with the requirements that every wjw_jwj​ is Hermitian and that, for all i,ji,ji,j,

wiwj+wjwi={2I,i=j,0,i≠j.w_iw_j+w_jw_i= \begin{cases} 2I,&i=j,\\ 0,&i\ne j. \end{cases}wi​wj​+wj​wi​={2I,0,​i=j,i=j.​

Thus each wi2=Iw_i^2=Iwi2​=I, and distinct matrices anticommute. The case n=0n=0n=0 is included: the underlying matrices have size 1×11\times11×1, but there are no wjw_jwj​, so both requirements are vacuous.

QAlgorithms.Papers.Independent.Intakethird100.P44.FixesObservables

For a finite index type HHH with decidable equality, a set A⊆Mat⁡(H)\mathcal A\subseteq\operatorname{Mat}(H)A⊆Mat(H), and an arbitrary map N:Mat⁡(H)→Mat⁡(H)N:\operatorname{Mat}(H)\to\operatorname{Mat}(H)N:Mat(H)→Mat(H), this predicate means that for every X∈AX\in\mathcal AX∈A and every a,b∈Ha,b\in Ha,b∈H,

Tr⁡(N(Eba)X)=Xab,\operatorname{Tr}\bigl(N(E_{ba})X\bigr)=X_{ab},Tr(N(Eba​)X)=Xab​,

where EbaE_{ba}Eba​ is the matrix with entry 111 at row bbb, column aaa, and zero elsewhere. No channel or linearity assumption on NNN is part of the predicate. The set A\mathcal AA may be empty, in which case the condition is vacuous.

QAlgorithms.Papers.Independent.Intakethird100.P44.HasIntervalSupport

For every natural number nnn, a family fff of Hermitian complex 2n×2n2^n\times2^n2n×2n matrices w0,…,w2n−1w_0,\ldots,w_{2n-1}w0​,…,w2n−1​ satisfying wiwj+wjwi=2Iw_iw_j+w_jw_i=2Iwi​wj​+wj​wi​=2I when i=ji=ji=j and 000 otherwise, a complex 2n×2n2^n\times2^n2n×2n matrix XXX, and a natural number ℓ\ellℓ, this predicate says that there exists a starting index s∈{0,…,2n−1}s\in\{0,\ldots,2n-1\}s∈{0,…,2n−1} such that XXX belongs to the complex linear span of all products

wj1⋯wj2hwith(ja+2n−s) mod 2n<ℓ for every a.w_{j_1}\cdots w_{j_{2h}} \quad\text{with}\quad (j_a+2n-s)\bmod 2n<\ell \ \text{for every }a.wj1​​⋯wj2h​​with(ja​+2n−s)mod2n<ℓ for every a.

Here hhh is any natural number, repeated indices and arbitrary ordering are permitted, and the empty product is III. Thus ℓ=0\ell=0ℓ=0 permits precisely scalar multiples of III when a starting index exists; for n>0n>0n>0 and ℓ≥2n\ell\ge 2nℓ≥2n, every mode index is permitted. For n=0n=0n=0, the predicate is false because there is no starting index.

QAlgorithms.Papers.Independent.Intakethird100.P44.IsChannelDilation

For finite index types H,K,EH,K,EH,K,E with decidable equality, an arbitrary map N:Mat⁡(H)→Mat⁡(K)N:\operatorname{Mat}(H)\to\operatorname{Mat}(K)N:Mat(H)→Mat(K), and a complex matrix VVV with rows indexed by K×EK\times EK×E and columns indexed by HHH, this predicate means

V†V=IHand(N(X))ab=∑e∈E(VXV†)(a,e),(b,e)V^\dagger V=I_H \quad\text{and}\quad (N(X))_{ab} =\sum_{e\in E}(VXV^\dagger)_{(a,e),(b,e)}V†V=IH​and(N(X))ab​=e∈E∑​(VXV†)(a,e),(b,e)​

for every matrix X∈Mat⁡(H)X\in\operatorname{Mat}(H)X∈Mat(H) and all a,b∈Ka,b\in Ka,b∈K. The equality concerns every input matrix. Empty index types are permitted, subject to the displayed isometry equation.

QAlgorithms.Papers.Independent.Intakethird100.P44.IsGeneralizedComplement

For finite index types H,K,EH,K,EH,K,E with decidable equality and maps N:Mat⁡(H)→Mat⁡(K)N:\operatorname{Mat}(H)\to\operatorname{Mat}(K)N:Mat(H)→Mat(K), C:Mat⁡(H)→Mat⁡(E)C:\operatorname{Mat}(H)\to\operatorname{Mat}(E)C:Mat(H)→Mat(E), this predicate asserts the existence of a natural number mmm, matrices Aj:H→KA_j:H\to KAj​:H→K for 0≤j<m0\le j<m0≤j<m, and maps R:Mat⁡(Fin m)→Mat⁡(E)R:\operatorname{Mat}(\mathrm{Fin}\,m)\to\operatorname{Mat}(E)R:Mat(Finm)→Mat(E), S:Mat⁡(E)→Mat⁡(Fin m)S:\operatorname{Mat}(E)\to\operatorname{Mat}(\mathrm{Fin}\,m)S:Mat(E)→Mat(Finm), such that

∑j<mAj†Aj=IH,N(X)=∑j<mAjXAj†,C=R∘CA,CA=S∘C,\sum_{j<m}A_j^\dagger A_j=I_H,\qquad N(X)=\sum_{j<m}A_jXA_j^\dagger,\qquad C=R\circ C_A,\qquad C_A=S\circ C,j<m∑​Aj†​Aj​=IH​,N(X)=j<m∑​Aj​XAj†​,C=R∘CA​,CA​=S∘C,

where CA(X)ij=Tr⁡(AiXAj†)C_A(X)_{ij}=\operatorname{Tr}(A_iXA_j^\dagger)CA​(X)ij​=Tr(Ai​XAj†​). Both RRR and SSS must be channels: each has a finite Kraus representation T(X)=∑ℓBℓXBℓ†T(X)=\sum_\ell B_\ell XB_\ell^\daggerT(X)=∑ℓ​Bℓ​XBℓ†​ with ∑ℓBℓ†Bℓ=I\sum_\ell B_\ell^\dagger B_\ell=I∑ℓ​Bℓ†​Bℓ​=I on its input space. All map equalities hold on every input matrix. The definition allows m=0m=0m=0 and empty ambient index types whenever these conditions can hold.

QAlgorithms.Papers.Independent.Intakethird100.P44.IsPhysicalChannel

For finite index types H,KH,KH,K with decidable equality and arbitrary maps P:Mat⁡(H)→Mat⁡(H)P:\operatorname{Mat}(H)\to\operatorname{Mat}(H)P:Mat(H)→Mat(H), Q:Mat⁡(K)→Mat⁡(K)Q:\operatorname{Mat}(K)\to\operatorname{Mat}(K)Q:Mat(K)→Mat(K), and N:Mat⁡(H)→Mat⁡(K)N:\operatorname{Mat}(H)\to\operatorname{Mat}(K)N:Mat(H)→Mat(K), this predicate means

Q(N(P(X)))=Q(N(X))for every X∈Mat⁡(H),Q(N(P(X)))=Q(N(X)) \quad\text{for every }X\in\operatorname{Mat}(H),Q(N(P(X)))=Q(N(X))for every X∈Mat(H),

where Mat⁡(I)\operatorname{Mat}(I)Mat(I) denotes the square complex matrices indexed by III. It imposes no channel, linearity, positivity, or normalization requirement on these maps, and XXX ranges over all matrices, not merely density matrices; empty index types are included.

QAlgorithms.Papers.Independent.Intakethird100.P44.ObservableAlgebra

For a finite type HHH with decidable equality, an observable algebra consists of a unital complex star subalgebra A\mathcal AA of the square complex matrices indexed by HHH, together with a map PPP on those matrices satisfying three conditions: PPP is a channel, meaning that there are m∈Nm\in\mathbb Nm∈N and matrices E1,…,EmE_1,\ldots,E_mE1​,…,Em​ with ∑jEj†Ej=I\sum_j E_j^\dagger E_j=I∑j​Ej†​Ej​=I and P(X)=∑jEjXEj†P(X)=\sum_j E_jXE_j^\daggerP(X)=∑j​Ej​XEj†​ for every matrix XXX; P(X)∈AP(X)\in\mathcal AP(X)∈A for every XXX; and

Tr⁡((X−P(X))†Y)=0for every matrix X and every Y∈A.\operatorname{Tr}\bigl((X-P(X))^\dagger Y\bigr)=0 \quad\text{for every matrix }X\text{ and every }Y\in\mathcal A.Tr((X−P(X))†Y)=0for every matrix X and every Y∈A.

The quantification over HHH permits the empty type, and the Kraus family may be empty when its completeness equation holds.

QAlgorithms.Papers.Independent.Intakethird100.P44.PhysicallyCorrectable

Let HHH be a finite index type with decidable equality, LLL a finite index type, P,N:Mat⁡(H)→Mat⁡(H)P,N:\operatorname{Mat}(H)\to\operatorname{Mat}(H)P,N:Mat(H)→Mat(H) arbitrary maps, and WWW an arbitrary complex matrix with rows indexed by HHH and columns by LLL. Physical correctability means that there exists a map R:Mat⁡(H)→Mat⁡(H)R:\operatorname{Mat}(H)\to\operatorname{Mat}(H)R:Mat(H)→Mat(H) admitting a finite Kraus representation R(Y)=∑aVaYVa†R(Y)=\sum_aV_aYV_a^\daggerR(Y)=∑a​Va​YVa†​ with ∑aVa†Va=I\sum_aV_a^\dagger V_a=I∑a​Va†​Va​=I, such that

P(R(P(Y)))=P(R(Y))for every Y∈Mat⁡(H),R(N(WXW†))=P(WXW†)for every X∈Mat⁡(L).P(R(P(Y)))=P(R(Y))\quad\text{for every }Y\in\operatorname{Mat}(H), \qquad R(N(WXW^\dagger))=P(WXW^\dagger)\quad\text{for every }X\in\operatorname{Mat}(L).P(R(P(Y)))=P(R(Y))for every Y∈Mat(H),R(N(WXW†))=P(WXW†)for every X∈Mat(L).

There is no requirement here that WWW be an isometry or that P,NP,NP,N themselves be channels; the second equality is required for all complex matrices XXX, including matrices that are not density matrices. Empty index types are permitted.

QAlgorithms.Papers.Independent.Intakethird100.P44.RegionalCoefficients

For finite index types H,L,JH,L,JH,L,J with decidable equality, a set AAA of complex H×HH\times HH×H matrices, matrices C∈CH×HC\in\mathbb C^{H\times H}C∈CH×H and W∈CH×LW\in\mathbb C^{H\times L}W∈CH×L, and a family (Ni)i∈J(N_i)_{i\in J}(Ni​)i∈J​ of complex H×HH\times HH×H matrices, this condition means that there exist complex scalars λijk\lambda_{ijk}λijk​, indexed by i,j∈Ji,j\in Ji,j∈J and Boolean kkk, such that, for every such i,j,ki,j,ki,j,k and every complex H×HH\times HH×H matrix YYY commuting with every element of AAA,

W†Ni†Nj(PkYPk)W=λijk W†(PkYPk)W,Pfalse=I+C2,Ptrue=I−C2.W^\dagger N_i^\dagger N_j(P_kYP_k)W =\lambda_{ijk}\,W^\dagger(P_kYP_k)W, \qquad P_{\mathrm{false}}=\frac{I+C}{2},\quad P_{\mathrm{true}}=\frac{I-C}{2}.W†Ni†​Nj​(Pk​YPk​)W=λijk​W†(Pk​YPk​)W,Pfalse​=2I+C​,Ptrue​=2I−C​.

No positivity, normalization, or compatibility conditions on the scalars are required. In particular, when JJJ is empty, the condition holds vacuously.

QAlgorithms.Papers.Independent.Intakethird100.P44.RegionalStateCondition

For finite index types H,L,JH,L,JH,L,J with decidable equality, a set AAA of complex H×HH\times HH×H matrices, matrices C∈CH×HC\in\mathbb C^{H\times H}C∈CH×H and W∈CH×LW\in\mathbb C^{H\times L}W∈CH×L, and a family (Ni)i∈J(N_i)_{i\in J}(Ni​)i∈J​ of complex H×HH\times HH×H matrices, this condition means that there exist two positive semidefinite J×JJ\times JJ×J matrices ρfalse,ρtrue\rho_{\mathrm{false}},\rho_{\mathrm{true}}ρfalse​,ρtrue​, each of trace 111, such that, for every H×HH\times HH×H matrix YYY commuting with every element of AAA and every complex J×JJ\times JJ×J matrix XXX,

∑i,j∈JXij W†Ni†YNjW=∑k∈{false,true}tr⁡(ρkX) W†PkYPkW,Pfalse=I+C2,Ptrue=I−C2.\sum_{i,j\in J}X_{ij}\,W^\dagger N_i^\dagger YN_jW = \sum_{k\in\{\mathrm{false},\mathrm{true}\}} \operatorname{tr}(\rho_kX)\,W^\dagger P_kYP_kW, \qquad P_{\mathrm{false}}=\frac{I+C}{2},\quad P_{\mathrm{true}}=\frac{I-C}{2}.i,j∈J∑​Xij​W†Ni†​YNj​W=k∈{false,true}∑​tr(ρk​X)W†Pk​YPk​W,Pfalse​=2I+C​,Ptrue​=2I−C​.

The matrices X,YX,YX,Y need not be states or Hermitian. There are no normalization assumptions on the noise family or WWW; if JJJ is empty, the required trace-one matrices cannot exist.

QAlgorithms.Papers.Independent.Intakethird100.P44.RingCode

For natural numbers n,kn,kn,k, a family (wa)a<2n(w_a)_{a<2n}(wa​)a<2n​ of Hermitian complex 2n×2n2^n\times2^n2n×2n matrices satisfying wawb+wbwa=2Iw_aw_b+w_bw_a=2Iwa​wb​+wb​wa​=2I for a=ba=ba=b and 000 otherwise, and a ring layout with n>0, k>0n>0,\ k>0n>0, k>0, a strictly increasing function f:{0,…,2k−1}→{0,…,2n−1}f:\{0,\ldots,2k-1\}\to\{0,\ldots,2n-1\}f:{0,…,2k−1}→{0,…,2n−1}, and even gaps gj=((f(j+)+2n)\dotminusf(j)\dotminus1) mod (2n)g_j=((f(j^+)+2n)\mathbin{\dotminus}f(j)\mathbin{\dotminus}1)\bmod(2n)gj​=((f(j+)+2n)\dotminusf(j)\dotminus1)mod(2n), a ring code consists of a complex 2n×2k2^n\times2^k2n×2k matrix WWW and proofs that

W†W=I,WW†=∏j=02k−1 ∏t=0gj/2−112(I−i w(f(j)+1+2t) mod (2n)w(f(j)+2+2t) mod (2n)).W^\dagger W=I,\qquad WW^\dagger= \prod_{j=0}^{2k-1}\ \prod_{t=0}^{g_j/2-1} \frac12\left( I-\mathrm{i}\, w_{(f(j)+1+2t)\bmod(2n)} w_{(f(j)+2+2t)\bmod(2n)} \right).W†W=I,WW†=j=0∏2k−1​ t=0∏gj​/2−1​21​(I−iw(f(j)+1+2t)mod(2n)​w(f(j)+2+2t)mod(2n)​).

Here j+j^+j+ is the cyclic successor, \dotminus\dotminus\dotminus denotes subtraction truncated at zero, and the product follows increasing ttt within increasing jjj, with empty products equal to III. Thus the required matrix is an isometry whose range projection equals exactly the displayed ordered product.

QAlgorithms.Papers.Independent.Intakethird100.P44.RingLayout

For natural numbers n,kn,kn,k, a ring layout consists of proofs that n>0n>0n>0 and k>0k>0k>0, a strictly increasing function f:{0,…,2k−1}→{0,…,2n−1}f:\{0,\ldots,2k-1\}\to\{0,\ldots,2n-1\}f:{0,…,2k−1}→{0,…,2n−1}, and a proof that every number

gf(j)=((f(j+)+2n)\dotminusf(j)\dotminus1) mod (2n)g_f(j)=\bigl((f(j^+)+2n)\mathbin{\dotminus}f(j)\mathbin{\dotminus}1\bigr)\bmod(2n)gf​(j)=((f(j+)+2n)\dotminusf(j)\dotminus1)mod(2n)

is even, where j+j^+j+ is the cyclic successor and \dotminus\dotminus\dotminus denotes subtraction truncated at zero. Strict increase implies injectivity and hence k≤nk\le nk≤n. There are no layouts with n=0n=0n=0, k=0k=0k=0, or k>nk>nk>n.

QAlgorithms.Papers.Independent.Intakethird100.P44.SameMaximum

For arbitrary types X,YX,YX,Y, predicates P:X→PropP:X\to\mathrm{Prop}P:X→Prop, Q:Y→PropQ:Y\to\mathrm{Prop}Q:Y→Prop, and real-valued functions f:X→Rf:X\to\mathbb Rf:X→R, g:Y→Rg:Y\to\mathbb Rg:Y→R, this predicate asserts that there exist x∈Xx\in Xx∈X and y∈Yy\in Yy∈Y such that P(x)P(x)P(x), Q(y)Q(y)Q(y), f(x)=g(y)f(x)=g(y)f(x)=g(y), every x′x'x′ satisfying P(x′)P(x')P(x′) has f(x′)≤f(x)f(x')\le f(x)f(x′)≤f(x), and every y′y'y′ satisfying Q(y′)Q(y')Q(y′) has g(y′)≤g(y)g(y')\le g(y)g(y′)≤g(y). Thus both maxima must be attained and equal; uniqueness is not required. If either allowed domain is empty, the predicate is false.

QAlgorithms.Papers.Independent.Intakethird100.P44.Sectors

For a finite index type HHH with decidable equality and a natural number rrr, sector data consist of complex matrices Pi∈Mat⁡(H)P_i\in\operatorname{Mat}(H)Pi​∈Mat(H) and real numbers viv_ivi​, indexed by 0≤i<r0\le i<r0≤i<r, such that every PiP_iPi​ is Hermitian and nonzero,

PiPj={Pi,i=j,0,i≠j,∑i=0r−1Pi=I,P_iP_j= \begin{cases} P_i,&i=j,\\ 0,&i\ne j, \end{cases} \qquad \sum_{i=0}^{r-1}P_i=I,Pi​Pj​={Pi​,0,​i=j,i=j,​i=0∑r−1​Pi​=I,

and the values viv_ivi​ are pairwise distinct. In particular, the matrices are mutually orthogonal projections. The definition permits r=0r=0r=0; then completeness requires 0=I0=I0=I, which holds precisely when HHH is empty. If HHH is empty, the nonzero requirement prevents r>0r>0r>0.

QAlgorithms.Papers.Independent.Intakethird100.P44.StrongCorrectable

For finite index types H,LH,LH,L with decidable equality, a set AAA of complex H×HH\times HH×H matrices, a complex H×HH\times HH×H matrix CCC, a function NNN from complex H×HH\times HH×H matrices to themselves, and a complex H×LH\times LH×L matrix WWW, this condition means that there exists a function RRR on complex H×HH\times HH×H matrices with a finite Kraus representation R(Z)=∑aEaZEa†R(Z)=\sum_a E_aZE_a^\daggerR(Z)=∑a​Ea​ZEa†​, where ∑aEa†Ea=I\sum_aE_a^\dagger E_a=I∑a​Ea†​Ea​=I, such that its trace dual

[R∨(Y)]ab=tr⁡ ⁣(R(Eba)Y)[R^\vee(Y)]_{ab}=\operatorname{tr}\!\bigl(R(E_{ba})Y\bigr)[R∨(Y)]ab​=tr(R(Eba​)Y)

maps every Y∈AY\in AY∈A into AAA and fixes every matrix commuting with all elements of AAA, and such that, for every complex L×LL\times LL×L matrix XXX,

R ⁣(N(WXW†))=12(WXW†+C WXW†C).R\!\left(N(WXW^\dagger)\right) =\frac12\left(WXW^\dagger+C\,WXW^\dagger C\right).R(N(WXW†))=21​(WXW†+CWXW†C).

Here EbaE_{ba}Eba​ denotes the matrix unit. The quantifier over XXX includes arbitrary matrices, and the definition imposes no channel assumption on NNN, no isometry assumption on WWW, and no Hermiticity or involution assumption on CCC.

QAlgorithms.Papers.Independent.Intakethird100.P44.StrongLocal

For a finite index type HHH with decidable equality, an arbitrary set AAA of complex matrices on HHH, and a map NNN on those matrices, this predicate requires that NNN be a channel, that its trace dual

N∨(X)ab=tr⁡ ⁣(N(Eba)X)N^\vee(X)_{ab}=\operatorname{tr}\!\bigl(N(E_{ba})X\bigr)N∨(X)ab​=tr(N(Eba​)X)

belong to AAA whenever X∈AX\in AX∈A, and that N∨(Y)=YN^\vee(Y)=YN∨(Y)=Y for every matrix YYY commuting with every member of AAA. Here EbaE_{ba}Eba​ is the matrix unit, and being a channel means admitting a finite representation N(X)=∑jKjXKj†N(X)=\sum_jK_jXK_j^\daggerN(X)=∑j​Kj​XKj†​ with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I. The set AAA is not required to be an algebra or nonempty; if A=∅A=\varnothingA=∅, the preservation condition is vacuous and the fixed-observable condition applies to every matrix.

QAlgorithms.Papers.Independent.Intakethird100.P44.TensorCorrectable

For finite index types H,B,LH,B,LH,B,L, with decidable equality, a map N:CH×H→CH×HN:\mathbb C^{H\times H}\to\mathbb C^{H\times H}N:CH×H→CH×H and a complex matrix W:CL→CH⊗CBW:\mathbb C^L\to\mathbb C^H\otimes\mathbb C^BW:CL→CH⊗CB, this predicate means that there exists a channel RRR on HHH such that, for every complex L×LL\times LL×L matrix XXX,

((R∘N)⊗idB)(WXW†)=WXW†.((R\circ N)\otimes\mathrm{id}_B)(WXW^\dagger)=WXW^\dagger.((R∘N)⊗idB​)(WXW†)=WXW†.

Here applying a map FFF to the first factor means that the ((h,b),(h′,b′))((h,b),(h',b'))((h,b),(h′,b′)) entry of its result on ZZZ is F((Z(a,b),(a′,b′))a,a′)h,h′F((Z_{(a,b),(a',b')})_{a,a'})_{h,h'}F((Z(a,b),(a′,b′)​)a,a′​)h,h′​. A channel means a map admitting a finite representation R(Y)=∑jEjYEj†R(Y)=\sum_j E_jYE_j^\daggerR(Y)=∑j​Ej​YEj†​ with ∑jEj†Ej=I\sum_jE_j^\dagger E_j=I∑j​Ej†​Ej​=I. The predicate imposes no channel assumption on NNN, no isometry assumption on WWW, and no nonemptiness assumption on the index types; its equality concerns all XXX, including matrices that are not density matrices.

QAlgorithms.Papers.Independent.Intakethird100.P44.codeSlice

For index types H,B,LH,B,LH,B,L, a complex matrix WWW with rows indexed by H×BH\times BH×B and columns indexed by LLL, and any n∈Bn\in Bn∈B, the slice at nnn is the matrix with rows indexed by HHH and columns indexed by LLL whose entries are

(Wn)aℓ=W(a,n),ℓ.(W_n)_{a\ell}=W_{(a,n),\ell}.(Wn​)aℓ​=W(a,n),ℓ​.

No normalization, positivity, or isometry condition on WWW is required.

QAlgorithms.Papers.Independent.Intakethird100.P44.comparisonRootFidelity

For finite index types H,KH,KH,K with decidable equality, a comparison ccc on HHH, and arbitrary maps N,M:Mat⁡(H)→Mat⁡(K)N,M:\operatorname{Mat}(H)\to\operatorname{Mat}(K)N,M:Mat(H)→Mat(K), let S(Z)S(Z)S(Z) denote square-root continuous functional calculus, put (Πρ)(a,b),(c,d)=S(ρ)baS(ρ)dc‾(\Pi_\rho)_{(a,b),(c,d)}=S(\rho)_{ba}\overline{S(\rho)_{dc}}(Πρ​)(a,b),(c,d)​=S(ρ)ba​S(ρ)dc​​, and define Fρ(T,U)=Re⁡Tr⁡S(S(A)BS(A))F_\rho(T,U)=\operatorname{Re}\operatorname{Tr}S(S(A)BS(A))Fρ​(T,U)=ReTrS(S(A)BS(A)), where A=(id⊗T)(Πρ)A=(\mathrm{id}\otimes T)(\Pi_\rho)A=(id⊗T)(Πρ​), B=(id⊗U)(Πρ)B=(\mathrm{id}\otimes U)(\Pi_\rho)B=(id⊗U)(Πρ​), and the reference extension applies the map separately to each matrix block indexed by the reference coordinates. If ccc is the fixed form with positive semidefinite trace-one matrix ρ\rhoρ, the declared value is Fρ(N,M)F_\rho(N,M)Fρ​(N,M). If ccc is the code form with k>0k>0k>0 and W†W=IkW^\dagger W=I_kW†W=Ik​, its value is the real infimum

inf⁡ρ∈Mat⁡({0,…,k−1})ρ⪰0, Tr⁡ρ=1Fρ ⁣(N∘(X↦WXW†),M∘(X↦WXW†)).\inf_{\substack{\rho\in\operatorname{Mat}(\{0,\ldots,k-1\})\\ \rho\succeq0,\ \operatorname{Tr}\rho=1}} F_\rho\!\left(N\circ(X\mapsto WXW^\dagger), M\circ(X\mapsto WXW^\dagger)\right).ρ∈Mat({0,…,k−1})ρ⪰0, Trρ=1​inf​Fρ​(N∘(X↦WXW†),M∘(X↦WXW†)).

No channel assumption on N,MN,MN,M, and no attainment of this infimum, is included.

QAlgorithms.Papers.Independent.Intakethird100.P44.constantChannel

For a finite input index type HHH, an output index type KKK, and any complex matrix σ∈Mat⁡(K)\sigma\in\operatorname{Mat}(K)σ∈Mat(K), this defines the map

X⟼tr⁡(X) σ(X∈Mat⁡(H)).X\longmapsto\operatorname{tr}(X)\,\sigma \qquad(X\in\operatorname{Mat}(H)).X⟼tr(X)σ(X∈Mat(H)).

No positivity or trace normalization of σ\sigmaσ is assumed, and no channel property is asserted. If HHH is empty, its trace is zero and the map is identically zero.

QAlgorithms.Papers.Independent.Intakethird100.P44.evenRegion

For every natural number nnn, a family of Hermitian 2n×2n2^n\times2^n2n×2n matrices wjw_jwj​, 0≤j<2n0\leq j<2n0≤j<2n, satisfying wiwj+wjwi=2δijIw_iw_j+w_jw_i=2\delta_{ij}Iwi​wj​+wj​wi​=2δij​I, and a finite subset S⊆{0,…,2n−1}S\subseteq\{0,\ldots,2n-1\}S⊆{0,…,2n−1}, this definition gives the complex linear span of all matrices

wj1wj2⋯wj2r,r≥0,j1,…,j2r∈S.w_{j_1}w_{j_2}\cdots w_{j_{2r}}, \qquad r\geq0,\quad j_1,\ldots,j_{2r}\in S.wj1​​wj2​​⋯wj2r​​,r≥0,j1​,…,j2r​∈S.

The indices may repeat, and their order is retained. The empty word is allowed and contributes III; consequently, when S=∅S=\varnothingS=∅, the defined space is the span of III. This also describes the case n=0n=0n=0.

QAlgorithms.Papers.Independent.Intakethird100.P44.fixedFidelity

For finite index types H,KH,KH,K with decidable equality, any matrix ρ∈Mat⁡(H)\rho\in\operatorname{Mat}(H)ρ∈Mat(H), and arbitrary maps N,M:Mat⁡(H)→Mat⁡(K)N,M:\operatorname{Mat}(H)\to\operatorname{Mat}(K)N,M:Mat(H)→Mat(K), define S(Z)S(Z)S(Z) to be the continuous functional calculus of the real square-root function applied to ZZZ, and define the matrix Πρ\Pi_\rhoΠρ​ on H×HH\times HH×H by

(Πρ)(a,b),(c,d)=S(ρ)ba S(ρ)dc‾.(\Pi_\rho)_{(a,b),(c,d)} =S(\rho)_{ba}\,\overline{S(\rho)_{dc}}.(Πρ​)(a,b),(c,d)​=S(ρ)ba​S(ρ)dc​​.

For any map TTT, its reference extension is defined entrywise by

[(idH⊗T)(X)](a,u),(b,v)=[T ⁣((X(a,i),(b,j))i,j∈H)]uv.[(\mathrm{id}_H\otimes T)(X)]_{(a,u),(b,v)} =\left[T\!\left((X_{(a,i),(b,j)})_{i,j\in H}\right)\right]_{uv}.[(idH​⊗T)(X)](a,u),(b,v)​=[T((X(a,i),(b,j)​)i,j∈H​)]uv​.

The declared real number is

Re⁡Tr⁡S ⁣(S(A)BS(A)),A=(idH⊗N)(Πρ),B=(idH⊗M)(Πρ).\operatorname{Re}\operatorname{Tr} S\!\left(S(A)BS(A)\right), \qquad A=(\mathrm{id}_H\otimes N)(\Pi_\rho),\quad B=(\mathrm{id}_H\otimes M)(\Pi_\rho).ReTrS(S(A)BS(A)),A=(idH​⊗N)(Πρ​),B=(idH​⊗M)(Πρ​).

There is no density-matrix assumption on ρ\rhoρ, no channel assumption on N,MN,MN,M, and no nonemptiness assumption on the index types; the functional calculus is used with its total definition for all matrices.

QAlgorithms.Papers.Independent.Intakethird100.P44.krausComplement

For finite index types H,K,JH,K,JH,K,J with decidable equality and any family of complex matrices AjA_jAj​ with rows indexed by KKK and columns indexed by HHH, this definition gives the map C:Mat⁡(H)→Mat⁡(J)C:\operatorname{Mat}(H)\to\operatorname{Mat}(J)C:Mat(H)→Mat(J) whose entries are

C(X)ij=Tr⁡(AiXAj†).C(X)_{ij}=\operatorname{Tr}(A_iXA_j^\dagger).C(X)ij​=Tr(Ai​XAj†​).

No equation ∑jAj†Aj=IH\sum_jA_j^\dagger A_j=I_H∑j​Aj†​Aj​=IH​, channel assumption, or restriction on the input matrix XXX is imposed. The indexing family may be empty.

QAlgorithms.Papers.Independent.Intakethird100.P44.localComplement

For finite index types H,EH,EH,E with decidable equality, let AAA consist of a unital complex star subalgebra A⊆Mat⁡(H)\mathcal A\subseteq\operatorname{Mat}(H)A⊆Mat(H) and a map PPP admitting a finite Kraus representation P(X)=∑jEjXEj†P(X)=\sum_jE_jXE_j^\daggerP(X)=∑j​Ej​XEj†​, ∑jEj†Ej=IH\sum_jE_j^\dagger E_j=I_H∑j​Ej†​Ej​=IH​, such that P(X)∈AP(X)\in\mathcal AP(X)∈A and Tr⁡((X−P(X))†Y)=0\operatorname{Tr}((X-P(X))^\dagger Y)=0Tr((X−P(X))†Y)=0 for every matrix XXX and every Y∈AY\in\mathcal AY∈A. For any matrix VVV with rows indexed by H×EH\times EH×E and columns indexed by HHH, the declared map LA,V:Mat⁡(H)→Mat⁡(H×E)L_{A,V}:\operatorname{Mat}(H)\to\operatorname{Mat}(H\times E)LA,V​:Mat(H)→Mat(H×E) is

LA,V(X)(a,e),(b,f)=[P ⁣(((VXV†)(i,e),(j,f))i,j∈H)]ab.L_{A,V}(X)_{(a,e),(b,f)} = \left[P\!\left(((VXV^\dagger)_{(i,e),(j,f)})_{i,j\in H}\right)\right]_{ab}.LA,V​(X)(a,e),(b,f)​=[P(((VXV†)(i,e),(j,f)​)i,j∈H​)]ab​.

No isometry condition on VVV is imposed here, and empty index types are included.

QAlgorithms.Papers.Independent.Intakethird100.P44.minimumGap

For natural numbers n,kn,kn,k and a ring layout rrr, let uj=r.free(j)u_j=r.\mathrm{free}(j)uj​=r.free(j) for j=0,…,2k−1j=0,\ldots,2k-1j=0,…,2k−1. A ring layout requires n>0n>0n>0, k>0k>0k>0, the indices uj∈{0,…,2n−1}u_j\in\{0,\ldots,2n-1\}uj​∈{0,…,2n−1} to be strictly increasing, and every cyclic gap

gj=(u(j+1) mod 2k+2n−uj−1) mod 2ng_j=(u_{(j+1)\bmod 2k}+2n-u_j-1)\bmod 2ngj​=(u(j+1)mod2k​+2n−uj​−1)mod2n

to be even. The defined natural number is

min⁡(2n,g0,…,g2k−1),\min(2n,g_0,\ldots,g_{2k-1}),min(2n,g0​,…,g2k−1​),

computed by starting at 2n2n2n and successively taking the minimum with each gap in increasing index order. Zero gaps are allowed; the gap list is nonempty because k>0k>0k>0.

QAlgorithms.Papers.Independent.Intakethird100.P44.modeAt

For every natural number nnn with a supplied proof that n>0n>0n>0, every index j∈{0,…,2n−1}j\in\{0,\ldots,2n-1\}j∈{0,…,2n−1}, and every natural-number offset aaa, this defines the index

(j+a) mod (2n)∈{0,…,2n−1}.(j+a)\bmod(2n)\in\{0,\ldots,2n-1\}.(j+a)mod(2n)∈{0,…,2n−1}.

Offsets may be zero or arbitrarily large; they wrap around modulo 2n2n2n. The supplied positivity proof ensures the modulus is nonzero and the result is a valid index.

QAlgorithms.Papers.Independent.Intakethird100.P44.pairingProjector

For natural numbers n,kn,kn,k, a family (wa)a<2n(w_a)_{a<2n}(wa​)a<2n​ of Hermitian complex 2n×2n2^n\times2^n2n×2n matrices satisfying wawb+wbwa=2Iw_aw_b+w_bw_a=2Iwa​wb​+wb​wa​=2I for a=ba=ba=b and 000 otherwise, and a ring layout consisting of n>0, k>0n>0,\ k>0n>0, k>0, a strictly increasing function f:{0,…,2k−1}→{0,…,2n−1}f:\{0,\ldots,2k-1\}\to\{0,\ldots,2n-1\}f:{0,…,2k−1}→{0,…,2n−1}, and even gaps gj=((f(j+)+2n)\dotminusf(j)\dotminus1) mod (2n)g_j=((f(j^+)+2n)\mathbin{\dotminus}f(j)\mathbin{\dotminus}1)\bmod(2n)gj​=((f(j+)+2n)\dotminusf(j)\dotminus1)mod(2n), this defines the matrix

∏j=02k−1 ∏t=0gj/2−112(I−i w(f(j)+1+2t) mod (2n)w(f(j)+2+2t) mod (2n)).\prod_{j=0}^{2k-1}\ \prod_{t=0}^{g_j/2-1} \frac12\left( I-\mathrm{i}\, w_{(f(j)+1+2t)\bmod(2n)} w_{(f(j)+2+2t)\bmod(2n)} \right).j=0∏2k−1​ t=0∏gj​/2−1​21​(I−iw(f(j)+1+2t)mod(2n)​w(f(j)+2+2t)mod(2n)​).

Here \dotminus\dotminus\dotminus is natural subtraction truncated at zero, j+j^+j+ is the cyclic successor, and the factors are multiplied first in increasing ttt within each gap and then in increasing jjj. A zero gap contributes no factors, and if all gaps are zero, the result is III.

QAlgorithms.Papers.Independent.Intakethird100.P44.parityProject

For a finite index type HHH with decidable equality and any complex matrix CCC on HHH, this definition assigns the map

PC(X)=12(X+CXC)P_C(X)=\frac12(X+CXC)PC​(X)=21​(X+CXC)

on all complex matrices XXX on HHH. It assumes no Hermiticity, unitarity, or involution property of CCC; the rightmost factor is CCC itself.

QAlgorithms.Papers.Independent.Intakethird100.P44.paritySector

For a finite index type HHH with decidable equality, every complex matrix CCC indexed by H×HH\times HH×H, and every Boolean kkk, this defines

Pk=12(I+(−1)kC),P_k=\frac12\bigl(I+(-1)^kC\bigr),Pk​=21​(I+(−1)kC),

where false corresponds to k=0k=0k=0 and true to k=1k=1k=1. Thus Pfalse=(I+C)/2P_{\mathrm{false}}=(I+C)/2Pfalse​=(I+C)/2 and Ptrue=(I−C)/2P_{\mathrm{true}}=(I-C)/2Ptrue​=(I−C)/2. No Hermiticity or involution assumption on CCC is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P44.poisoning

For natural numbers n,kn,kn,k, let fff consist of Hermitian complex 2n×2n2^n\times2^n2n×2n matrices wjw_jwj​, indexed by 0≤j<2n0\le j<2n0≤j<2n, with wiwj+wjwi=2δijIw_iw_j+w_jw_i=2\delta_{ij}Iwi​wj​+wj​wi​=2δij​I, and let rrr be a ring layout: n,k>0n,k>0n,k>0, its free indices uj=r.free(j)u_j=r.\mathrm{free}(j)uj​=r.free(j) are strictly increasing for 0≤j<2k0\le j<2k0≤j<2k, and all cyclic gaps (u(j+1) mod 2k+2n−uj−1) mod 2n(u_{(j+1)\bmod 2k}+2n-u_j-1)\bmod 2n(u(j+1)mod2k​+2n−uj​−1)mod2n are even. The defined map sends every complex 2n×2n2^n\times2^n2n×2n matrix XXX to

N(X)=12(wu0Xwu0†+wu1Xwu1†).\mathcal N(X)=\frac12\left(w_{u_0}Xw_{u_0}^{\dagger} +w_{u_1}Xw_{u_1}^{\dagger}\right).N(X)=21​(wu0​​Xwu0​†​+wu1​​Xwu1​†​).

The two selected free indices have zero-based positions 000 and 111; their existence follows from k>0k>0k>0. No positivity or normalization condition is imposed on the argument XXX.

QAlgorithms.Papers.Independent.Intakethird100.P44.regionParity

For every natural number nnn, every family f=(wj)0≤j<2nf=(w_j)_{0\le j<2n}f=(wj​)0≤j<2n​ of complex 2n×2n2^n\times2^n2n×2n matrices satisfying wj†=wjw_j^\dagger=w_jwj†​=wj​ and wiwj+wjwi=2Iw_iw_j+w_jw_i=2Iwi​wj​+wj​wi​=2I when i=ji=ji=j, and 000 otherwise, and every finite set ω⊆{0,…,n−1}\omega\subseteq\{0,\ldots,n-1\}ω⊆{0,…,n−1}, this defines

Cω=i∣ω∣∏0≤j<2n∃s∈ω, ⌊j/2⌋=sin increasing order of jwj.C_\omega=\mathrm{i}^{|\omega|} \prod_{\substack{0\le j<2n\\ \exists s\in\omega,\ \lfloor j/2\rfloor=s}}^{\text{in increasing order of }j}w_j.Cω​=i∣ω∣0≤j<2n∃s∈ω, ⌊j/2⌋=s​∏in increasing order of j​wj​.

The multiplication order is the increasing order of the selected mode indices. The empty product is III, so ω=∅\omega=\varnothingω=∅, including the case n=0n=0n=0, gives Cω=IC_\omega=ICω​=I.

QAlgorithms.Papers.Independent.Intakethird100.P44.ringGap

For natural numbers n,kn,kn,k, any function f:{0,…,2k−1}→{0,…,2n−1}f:\{0,\ldots,2k-1\}\to\{0,\ldots,2n-1\}f:{0,…,2k−1}→{0,…,2n−1}, and any j<2kj<2kj<2k, this defines

gf(j)=((f(j+)+2n)\dotminusf(j)\dotminus1) mod (2n),g_f(j)=\bigl((f(j^+)+2n)\mathbin{\dotminus}f(j)\mathbin{\dotminus}1\bigr)\bmod(2n),gf​(j)=((f(j+)+2n)\dotminusf(j)\dotminus1)mod(2n),

where j+j^+j+ is the cyclic successor of jjj, and a\dotminusb=max⁡(a−b,0)a\mathbin{\dotminus}b=\max(a-b,0)a\dotminusb=max(a−b,0) is natural-number subtraction, performed in the displayed order. No injectivity, ordering, or positivity hypotheses are imposed. When k=0k=0k=0, there is no index jjj; when n=0n=0n=0, a function fff together with such a jjj cannot exist.

QAlgorithms.Papers.Independent.Intakethird100.P44.sectorChargeObservable

Let HHH be a finite index type with decidable equality, rrr a natural number, and suppose the sector data comprise nonzero Hermitian matrices Pi∈Mat⁡(H)P_i\in\operatorname{Mat}(H)Pi​∈Mat(H), mutually orthogonal with Pi2=PiP_i^2=P_iPi2​=Pi​ and ∑iPi=I\sum_iP_i=I∑i​Pi​=I, and pairwise distinct real numbers viv_ivi​, for 0≤i<r0\le i<r0≤i<r. The associated charge observable is the complex matrix

C=∑i=0r−1viPi.C=\sum_{i=0}^{r-1}v_iP_i.C=i=0∑r−1​vi​Pi​.

For r=0r=0r=0, this is the empty sum 000, and the sector assumptions require HHH to be empty.

QAlgorithms.Papers.Independent.Intakethird100.P44.sectorProject

Let HHH be a finite index type with decidable equality, rrr a natural number, and let Pi∈Mat⁡(H)P_i\in\operatorname{Mat}(H)Pi​∈Mat(H) be nonzero Hermitian matrices satisfying PiPj=PiP_iP_j=P_iPi​Pj​=Pi​ when i=ji=ji=j, PiPj=0P_iP_j=0Pi​Pj​=0 otherwise, and ∑iPi=I\sum_iP_i=I∑i​Pi​=I, together with pairwise distinct real values viv_ivi​. The associated map on all complex matrices on HHH is

X⟼∑i=0r−1PiXPi.X\longmapsto\sum_{i=0}^{r-1}P_iXP_i.X⟼i=0∑r−1​Pi​XPi​.

The values viv_ivi​ do not enter this formula. For r=0r=0r=0, the sum is zero and the sector assumptions force HHH to be empty.

QAlgorithms.Papers.Independent.Intakethird100.P44.siteModes

For every natural number nnn, including n=0n=0n=0, and every subset ω\omegaω of {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1}, this defines the finite set

Sω={j∈{0,…,2n−1}:∃s∈ω, ⌊j/2⌋=s}.S_\omega=\{j\in\{0,\ldots,2n-1\}:\exists s\in\omega,\ \lfloor j/2\rfloor=s\}.Sω​={j∈{0,…,2n−1}:∃s∈ω, ⌊j/2⌋=s}.

Thus each site s∈ωs\in\omegas∈ω contributes the two indices 2s2s2s and 2s+12s+12s+1. If ω\omegaω is empty, including necessarily when n=0n=0n=0, the result is empty.

QAlgorithms.Papers.Independent.Intakethird100.P44.totalParity

For any natural number nnn and any family fff of Hermitian 2n×2n2^n\times2^n2n×2n matrices wjw_jwj​, indexed by 0≤j<2n0\leq j<2n0≤j<2n, satisfying wiwj+wjwi=2δijIw_iw_j+w_jw_i=2\delta_{ij}Iwi​wj​+wj​wi​=2δij​I, this definition assigns the matrix

Cf=i nw0w1⋯w2n−1,C_f=\mathrm{i}^{\,n}w_0w_1\cdots w_{2n-1},Cf​=inw0​w1​⋯w2n−1​,

with the factors multiplied in increasing index order. For n=0n=0n=0, the product is empty and Cf=IC_f=ICf​=I on the one-dimensional space.

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