Documentation

LeanPool.LowWeightPauliDynamics.Pauli.Discard

The operator discarded by a Pauli truncation #

The component Õ^{(d)}_{≥w*+1} of apd:eq:step_component is an operator difference, (1 - Π_{≤ w*}) A, where A is the evolved observable immediately before the truncation that ends step d. The scalar highNorm w A measures the high-weight coefficients of A, but that definition alone does not identify it with the norm of the discarded operator. This file proves the identification.

sub_truncOp proves the general operator identity A - Π_S A = Π_{Sᶜ} A, using completeness of the Pauli expansion. pauliNorm_sub_truncOp_highSet_compl then identifies the norm of the operator discarded by a weight cut with highNorm. This is the norm identification in the last equality of apd:thm:triangle, not that proposition's telescoping or expectation bound. No trajectory, truncation schedule, or entanglement hypothesis enters here; trajectories and schedules are in Pauli/TrotterTruncate, and the telescoping bound is in Pauli/TruncationError.

All identities hold for arbitrary complex matrices indexed by bit strings, with Pauli strings represented by the entrywise model toMatrix.

Main results #

theorem Lean4LPD.PauliString.coeff_sub {n : ℕ} (A B : Matrix (Bits n) (Bits n) ℂ) (p : PauliIndex n) :
coeff (A - B) p = coeff A p - coeff B p

Coefficients preserve operator subtraction, an algebraic helper for the discarded component (1 - Π_{≤ w*}) A in apd:eq:step_component.

Coefficient vectors preserve subtraction, so the operator difference of apd:eq:step_component can be compared in coefficient space.

theorem Lean4LPD.PauliString.coeff_ext {n : ℕ} {A B : Matrix (Bits n) (Bits n) ℂ} (h : ∀ (p : PauliIndex n), coeff A p = coeff B p) :
A = B

Equality of every Pauli coefficient determines an operator. This is the completeness helper used to identify the discarded operator in apd:eq:step_component, rather than merely matching a scalar norm.

theorem Lean4LPD.PauliString.sub_truncOp {n : ℕ} (S : Finset (PauliIndex n)) (O : Matrix (Bits n) (Bits n) ℂ) :
O - truncOp S O = truncOp Sᶜ O

The discarded operator is the complementary projection. This generalizes the 1 - Π_{≤ w*} in apd:eq:step_component to any retained set S.

The complementary retained set discards exactly the S component; at S = highSet n w* this is the operator identity in apd:eq:step_component.

The coefficient vector of the discarded operator O - Π_S O is the coefficient vector of O restricted to the complement Sᶜ; a coefficient-space form of apd:eq:step_component.

A projected operator's Pauli norm is its restricted coefficient norm, the norm bridge used for the discarded component in apd:thm:triangle.

The norm of the operator discarded by retaining S is the norm on Sᶜ. This generalizes the final norm identification of apd:thm:triangle.

The norm of the operator discarded by retaining Sᶜ is the norm on S, the general-set form of the final norm identification in apd:thm:triangle.

The high-weight projection has norm highNorm, identifying the high-weight operator norm appearing in apd:thm:triangle.

The high-weight scalar is the norm of the discarded operator. Retaining the classes of weight at most w discards an operator whose Pauli norm is exactly highNorm w O. This connects the component Õ^{(d)}_{≥w*+1} of apd:eq:step_component with the high-weight norm apd:eq:def_high_weight_norm that the ladder bounds (the identification of the truncated mass at the start of the proof of apd:thm:one_step_truncation_error), and it holds for an arbitrary operator O.