Documentation

MIPRE.Background.LIDT.MIPStarRE.LDT.MakingMeasurementsProjective.QXPLayerIdentities.LayerAlgebra

Section 5 — Q/X/XHat/P algebraic identities #

Algebraic identities for the Q/X/XHat/P layer, including the restatements of Q_a, P_a, the mixed product, and projectivity of P.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.xa_t {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (a : Outcome) :
Xa data a = Ta data.qLayer a * data.x

X_a = T_a X (lem:xa-t).

theorem MIPStarRE.LDT.MakingMeasurementsProjective.qaRestated {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (a : Outcome) :
Qa data.qLayer a = (Xa data a).conjTranspose * Xa data a Qa data.qLayer a = data.x.conjTranspose * Ta data.qLayer a * data.x Qa data.qLayer a = (Xa data a).conjTranspose * data.x

Q_a restated (lem:qa-restated).

Rewrites the paper's operator Q_a in terms of X_a, X, and T_a.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.xHat_eq_x_of_x_mul_conjTranspose_eq_one_of_mixed_eq_gram {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (hx_left : data.x * data.x.conjTranspose = 1) (hmixed : data.x.conjTranspose * data.xHat = data.x.conjTranspose * data.x) :
data.xHat = data.x

If the original X rows are already coisometric and the mixed product X† XHat agrees with the Gram operator X† X, then the polar replacement XHat is equal to X.

This is the algebraic form of the observation used in the residual-domination route: on a part where the row block already lies in the unit singular subspace, the XHat replacement does not change it.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.xHatA_eq_xa_of_x_mul_conjTranspose_eq_one_of_mixed_eq_gram {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (a : Outcome) (hx_left : data.x * data.x.conjTranspose = 1) (hmixed : data.x.conjTranspose * data.xHat = data.x.conjTranspose * data.x) :
XHatA data a = Xa data a

Row-block form of xHat_eq_x_of_x_mul_conjTranspose_eq_one_of_mixed_eq_gram.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.xHat_eq_x_of_x_mul_conjTranspose_eq_one {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (hx_left : data.x * data.x.conjTranspose = 1) :
data.xHat = data.x

If the original X rows are already coisometric, then the polar replacement XHat is equal to X.

Indeed, X X† = I makes the Gram operator X† X idempotent. Since this Gram operator is positive, its positive square root is itself; the stored mixed identity X† XHat = sqrt (X† X) therefore reduces to the hypothesis of xHat_eq_x_of_x_mul_conjTranspose_eq_one_of_mixed_eq_gram.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.xHatA_eq_xa_of_x_mul_conjTranspose_eq_one {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (a : Outcome) (hx_left : data.x * data.x.conjTranspose = 1) :
XHatA data a = Xa data a

Row-block form of xHat_eq_x_of_x_mul_conjTranspose_eq_one.

The total Q operator in a QXP layer is Hermitian.

This follows from lem:X-squared: the total operator is the right Gram matrix X†X. The statement gives downstream spectral arguments a canonical Hermitian witness for QTotal data.qLayer.

The total Q operator in a QXP layer is positive semidefinite.

This is the positivity companion to qtotal_isHermitian_of_x_squared; after QTotal is identified with the Gram matrix X†X, positivity follows from the standard Gram-matrix argument.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.qtotal_eigenvalues_nonneg_of_x_squared {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (i : ι) :

The spectral eigenvalues of the total Q operator are nonnegative.

This is the scalar form of qtotal_posSemidef_of_x_squared used when the rectangular polar construction separates the positive and zero spectral subspaces.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.exists_unitaryGroup_with_qxp_positive_gram_spectrum_rows {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) :
∃ (e : { i : ι // 0 < .eigenvalues i } data.qLayer.auxSpace.carrier) (U : (Matrix.unitaryGroup data.qLayer.auxSpace.carrier )), ∀ (i : { i : ι // 0 < .eigenvalues i }) (r : data.qLayer.auxSpace.carrier), U (e i) r = positiveGramSpectrumImageRows data.x (QTotal data.qLayer) i r

QXP-layer form of the positive-Gram row extension theorem.

For the matrix X stored in a QXP layer, the positive spectral part of Q = X†X determines normalized image rows. These rows can be placed in distinct auxiliary coordinates and then completed to a square unitary group element on the auxiliary Hilbert space.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.xExpressionToQExpression {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (a : Outcome) :
(Xa data a).conjTranspose * ((data.x * data.x.conjTranspose - 1) * (data.x * data.x.conjTranspose - 1)) * Xa data a = Qa data.qLayer a * QTotal data.qLayer * Qa data.qLayer a - Qa data.qLayer a

X-expression to Q-expression (lem:X-expression-to-Q-expression).

Converts the quadratic error term in X X† - I to the corresponding Q_a Q Q_a - Q_a expression.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.paRestated {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (a : Outcome) :
Pa data a = data.xHat.conjTranspose * Ta data.qLayer a * data.xHat Pa data a = (XHatA data a).conjTranspose * data.xHat

P_a restated (lem:pa-restated).

Rewrites P_a in terms of XHat, XHat_a, and T_a.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.pa_eq_xHatA_adjoint_mul_xHatA {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (a : Outcome) :
Pa data a = (XHatA data a).conjTranspose * XHatA data a

The P_a operator is the right Gram matrix of the corresponding XHat_a row block.

This is the row-block form of paRestated: the auxiliary projector T_a is idempotent, so inserting the second copy of T_a does not change the operator.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.qa_eq_pa_of_xHatA_eq_xa {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (a : Outcome) (hrow : XHatA data a = Xa data a) :
Qa data.qLayer a = Pa data a

If the XHat construction preserves one row block, then the corresponding Q and P outcomes agree.

This lemma isolates the row-block identity used for the fresh outcome in the option-completed orthonormalization step.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.fresh_outcome_le_of_xHatA_eq_xa {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData (Option Outcome) ι) (hrow : XHatA data none = Xa data none) :

Fresh-outcome domination follows from preservation of the fresh row block in the option-completed QXP layer.

This is only a QXP-internal comparison: it upgrades fresh-row preservation to Q_none ≤ P_none. The former generic RestrictSome monotone-total route would also have needed a source-to-Q comparison (optionCompletion A).outcome none ≤ Q_none; that comparison is not part of this lemma, and the present formal development no longer uses that generic route. See docs/reports/issue-1642-restrictsome-residual-domination-obstruction.md.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.xHat_adjoint_mul_x_eq_sqrt {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) :

The adjoint mixed product Xhat† X equals the positive square root of Q.

This is the adjoint form of the stored identity X† Xhat = sqrt Q. It is used to identify the operator Y = X Xhat† in the proof of lem:squared-difference.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.x_mul_xHat_adjoint_isHermitian {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) :

The operator X Xhat† is Hermitian in a QXP layer.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.x_mul_xHat_adjoint_sq {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) :
data.x * data.xHat.conjTranspose * (data.x * data.xHat.conjTranspose) = data.x * data.x.conjTranspose

The square of X Xhat† is X X† in a QXP layer.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.x_mul_xHat_adjoint_nonneg {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) :
0 data.x * data.xHat.conjTranspose

The operator X Xhat† is positive semidefinite in a QXP layer.

Spectral form of the paper's identity X * Xhat† = U * Σ * U† in lem:X-times-X-hat.

The unitary is the Mathlib eigenvector unitary of the already proved Hermitian operator X * Xhat†; the diagonal matrix is the corresponding real spectrum, embedded in . This is the unconditional form used by the later lem:squared-difference argument.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.squaredDifference {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) :
(data.x - data.xHat) * (data.x - data.xHat).conjTranspose (data.x * data.x.conjTranspose - 1) * (data.x * data.x.conjTranspose - 1)

Squared difference (lem:squared-difference).

Bounds the defect between X and XHat by the squared defect of X X† from the auxiliary identity.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.sum_pa_eq_xHat_adjoint_mul_xHat {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) :
a : Outcome, Pa data a = data.xHat.conjTranspose * data.xHat

The sum of the QXP P-operators is the Gram operator XHat† XHat.

This is the total-mass identity for the canonical projective submeasurement produced from the Q/X/XHat/P layer.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.pProjectivity {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) :
∃ (P : ProjSubMeas Outcome ι), ∀ (a : Outcome), P.outcome a = Pa data a

Projectivity of P (lem:P-projectivity).

The family P_a built from XHat and T_a is a projective submeasurement.

noncomputable def MIPStarRE.LDT.MakingMeasurementsProjective.qxpProjSubMeas {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) :
ProjSubMeas Outcome ι

The canonical projective submeasurement obtained from the Q/X/XHat/P layer. Its outcomes are the paper's operators P_a.

Equations
Instances For
    @[simp]
    theorem MIPStarRE.LDT.MakingMeasurementsProjective.qxpProjSubMeas_outcome {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (a : Outcome) :
    (qxpProjSubMeas data).outcome a = Pa data a

    The total of the canonical QXP projective submeasurement is XHat† XHat.

    This exposes the total-mass identity implicit in pProjectivity, which is needed when comparing the repaired projective family to the source submeasurement in later monotonicity arguments.

    theorem MIPStarRE.LDT.MakingMeasurementsProjective.qxpProjSubMeas_total_ev_eq_sum_pa_ev {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (ψ : QuantumState ι) :
    ev ψ (qxpProjSubMeas data).total = a : Outcome, ev ψ (Pa data a)

    The expectation of the QXP projective total is the sum of the expectations of the paper projectors P_a.

    theorem MIPStarRE.LDT.MakingMeasurementsProjective.qxpProjSubMeas_total_le_of_outcome_le {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (A : SubMeas Outcome ι) (hpoint : ∀ (a : Outcome), Pa data a A.outcome a) :

    Outcomewise domination of the QXP projectors implies domination of the total operator.

    This is the summation form of the monotone-total invariant needed downstream: once the concrete repair proves P_a ≤ A_a for every outcome, the canonical QXP projective total is bounded by the source submeasurement total.

    structure MIPStarRE.LDT.MakingMeasurementsProjective.QXPLayerTotalDomination {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) (A : SubMeas Outcome ι) :

    Total-domination invariant for the QXP repair.

    This proposition is the construction-level operator comparison required by the paper-tight monotone-total route: the canonical projective family obtained from the QXP layer has total operator bounded by the source submeasurement total. It is deliberately stronger than state-dependent-distance closeness and should be proved from the concrete repair, not inferred from the orthonormalization error estimate alone.

    Instances For
      theorem MIPStarRE.LDT.MakingMeasurementsProjective.QXPLayerTotalDomination.of_outcome_le {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] {data : QXPLayerData Outcome ι} {A : SubMeas Outcome ι} (hpoint : ∀ (a : Outcome), Pa data a A.outcome a) :

      Outcomewise operator domination is a sufficient way to prove the QXP total-domination invariant.

      theorem MIPStarRE.LDT.MakingMeasurementsProjective.QXPLayerTotalDomination.rightTensor_total_ev_le {Outcome : Type u_1} {ιLeft : Type u_2} {ι : Type u_3} [Fintype ιLeft] [DecidableEq ιLeft] [Fintype ι] [DecidableEq ι] [Fintype Outcome] {ψ : QuantumState (ιLeft × ι)} {data : QXPLayerData Outcome ι} {A : SubMeas Outcome ι} (hdom : QXPLayerTotalDomination data A) :

      A QXP total-domination witness gives the scalar right-register comparison used by the final-fields transport.