Skip to content
Shadow Theory

Technical Companion: Approximation and Accounting for Quantum Sources

A robust dimension witness for arbitrary four-record approximations

Jeremy RodgersSection 5 of 10Original PDF ↗

A robust dimension witness for arbitrary four-record approximations

The low-noise support argument excludes qutrit memories on a small exact-noise interval. That does not by itself quantify approximation by an arbitrary qutrit process whose perturbed law need not remain in that one-parameter family. The same support argument gives an observable matrix witness and a nonzero finite-record approximation floor. The constants are deliberately conservative; they are not proposed experimental accuracy targets.

Matrix witness. Let C ≥0 be a 9-by-9 column-stochastic matrix, indexed by pairs (a, b), (c, d) ∈ {0, 1, 2}2. Define

"X
#
X
d
Cab,ad +
X
μ = min
a,b
c,d
Cab,cd > 0,
L = max
a,b
c
Cab,cb
,
Source page 9.
Searchable transcription
"X
#
X
d
Cab,ad +
X
μ = min
a,b
c,d
Cab,cd > 0,
L = max
a,b
c
Cab,cb
,
and let σ2, σ5 > 0 be its second and fifth singular values. If C has a complex Hermitian PSD factorization
of size at most three, then
27L
μ
≥
σ2
2σ2
5
36(σ2 + 18
√
6)2 .
(24)
Technical passage, source page 9.
Searchable transcription
and let σ2, σ5 > 0 be its second and fifth singular values. If C has a complex Hermitian PSD factorization
of size at most three, then
27L
μ
≥
σ2
2σ2
5
36(σ2 + 18
√
6)2 .
(24)

The indices, leakage and singular values are ordinary classical matrix data. No particular quantum update architecture or reference frame is assumed in this lower bound.

Proof.
Since σ5 > 0, a factorization cannot be supported on dimension two. Congruence normalization
of its row factors gives a POVM P
a,b Aab = I3 and density matrices bcd with Cab,cd = Tr(Aabbcd).
The row sums imply Tr Aab ≥μ/9, hence a largest eigenvalue at least μ/27. Put Ra = P
d bad and
Qb = P
c bcb. Both have trace three, and the leakage implies
Technical passage, source page 9.
Searchable transcription
Proof.
Since σ5 > 0, a factorization cannot be supported on dimension two. Congruence normalization
of its row factors gives a POVM P
a,b Aab = I3 and density matrices bcd with Cab,cd = Tr(Aabbcd).
The row sums imply Tr Aab ≥μ/9, hence a largest eigenvalue at least μ/27. Put Ra = P
d bad and
Qb = P
c bcb. Both have trace three, and the leakage implies
λmin(Ra + Qb) ≤e := 27L/μ.
Source page 9.
Searchable transcription
λmin(Ra + Qb) ≤e := 27L/μ.

We prove the required approximate plane/line dichotomy directly. Fix any s > e; all distances below are on the same normalized factorization. Case 1: some Ra or Qb has second eigenvalue at least s. By interchanging the two indices if necessary, take such an Ra. Let P project onto its two largest eigenvectors and let z be a unit vector in the complementary line. For each b, choose a unit nb with ⟨nb, (Ra + Qb)nb⟩≤e. Positivity and Ra ⪰sP imply ∥Pnb∥2 ≤e/s, ⟨nb, Qbnb⟩≤e.

Choose the phase so ⟨z, nb⟩≥0. Since e/s < 1,
Source page 10.
Searchable transcription
Choose the phase so ⟨z, nb⟩≥0. Since e/s < 1,
∥nb −z∥2 = 2(1 −
q
1 −∥Pnb∥2) ≤2e/s.
Source page 10.
Searchable transcription
∥nb −z∥2 = 2(1 −
q
1 −∥Pnb∥2) ≤2e/s.
The bound ∥Qb∥≤Tr Qb = 3 and the triangle inequality after applying √Qb give
Technical passage, source page 10.
Searchable transcription
The bound ∥Qb∥≤Tr Qb = 3 and the triangle inequality after applying √Qb give
⟨z, Qbz⟩≤
√e +
√
2e/s
2
=: l= e(1 +
q
3
q
6/s)2.
Source page 10.
Searchable transcription
⟨z, Qbz⟩≤
√e +
√
2e/s
2
=: l= e(1 +
q
3
q
6/s)2.
Every preparation obeys bcd ⪯Qd, hence Tr[(I −P)bcd] ≤lfor this common plane. The projection
bound ∥ρ −PρP∥1 ≤2
p
Technical passage, source page 10.
Searchable transcription
Every preparation obeys bcd ⪯Qd, hence Tr[(I −P)bcd] ≤lfor this common plane. The projection
bound ∥ρ −PρP∥1 ≤2
p
Tr[(I −P)ρ] therefore bounds every preparation-compression error by 2
√
Technical passage, source page 10.
Searchable transcription
Tr[(I −P)ρ] therefore bounds every preparation-compression error by 2
√
l.
For completeness, purify ρ, apply P ⊗I, and compute the trace norm of the difference between the
original rank-one projector and the projected, unnormalized rank-one operator: for discarded mass x,
its norm is
√
Technical passage, source page 10.
Searchable transcription
l.
For completeness, purify ρ, apply P ⊗I, and compute the trace norm of the difference between the
original rank-one projector and the projected, unnormalized rank-one operator: for discarded mass x,
its norm is
√
4x −3x2 ≤2√x. Partial trace contracts this Hermitian trace norm, giving the stated
bound also for mixed ρ.
Replace each bcd by PbcdP and apply the same POVM. The resulting real matrix C(4) has rank
at most four, since Hermitian operators on a two-dimensional plane have real dimension four. A
POVM contracts the trace norm of a Hermitian difference: if D = D+ −D−is its positive/negative
decomposition, then P
ab | Tr(AabD)| ≤Tr D+ + Tr D−= ∥D∥1. Consequently each of the nine column
errors has Euclidean norm at most 2
√
l. Thus
Technical passage, source page 10.
Searchable transcription
4x −3x2 ≤2√x. Partial trace contracts this Hermitian trace norm, giving the stated
bound also for mixed ρ.
Replace each bcd by PbcdP and apply the same POVM. The resulting real matrix C(4) has rank
at most four, since Hermitian operators on a two-dimensional plane have real dimension four. A
POVM contracts the trace norm of a Hermitian difference: if D = D+ −D−is its positive/negative
decomposition, then P
ab | Tr(AabD)| ≤Tr D+ + Tr D−= ∥D∥1. Consequently each of the nine column
errors has Euclidean norm at most 2
√
l. Thus
σ5 ≤∥C −C(4)∥op ≤∥C −C(4)∥F ≤6√e(1 +
q
6/s).
Source page 10.
Searchable transcription
σ5 ≤∥C −C(4)∥op ≤∥C −C(4)∥F ≤6√e(1 +
q
6/s).
The first inequality follows from the singular-value lower bound for rank-four approximation, which
applies to every real approximant.
Case 2: all Ra, Qb have second eigenvalue below s. Let Pa and Q′
b be their top-eigenvector projectors.
The mass outside each such line is below 2s. Since bab ⪯Ra, Qb,
Technical passage, source page 10.
Searchable transcription
The first inequality follows from the singular-value lower bound for rank-four approximation, which
applies to every real approximant.
Case 2: all Ra, Qb have second eigenvalue below s. Let Pa and Q′
b be their top-eigenvector projectors.
The mass outside each such line is below 2s. Since bab ⪯Ra, Qb,
D(bab, Pa) ≤
√
2s,
D(bab, Q′
b) ≤
√
2s,
Source page 10.
Searchable transcription
D(bab, Pa) ≤
√
2s,
D(bab, Q′
b) ≤
√
2s,
where D denotes half trace distance. Here D(ρ, |v⟩⟨v|) ≤
p
Technical passage, source page 10.
Searchable transcription
where D denotes half trace distance. Here D(ρ, |v⟩⟨v|) ≤
p
1 −⟨v, ρv⟩: it is equality for pure states
and follows for mixed states by convexity and Jensen’s inequality. Fix P0. Comparing through b0b gives
D(P0, Q′
b) ≤2
√
2s, whence D(bab, P0) ≤3
√
Technical passage, source page 10.
Searchable transcription
1 −⟨v, ρv⟩: it is equality for pure states
and follows for mixed states by convexity and Jensen’s inequality. Fix P0. Comparing through b0b gives
D(P0, Q′
b) ≤2
√
2s, whence D(bab, P0) ≤3
√
2s for all nine preparations. Replace every preparation by
this same pure state. The resulting matrix C(1) has identical normalized columns and rank at most
one. Each column difference is zero-sum and has classical total variation at most 3
√
Technical passage, source page 10.
Searchable transcription
2s for all nine preparations. Replace every preparation by
this same pure state. The resulting matrix C(1) has identical normalized columns and rank at most
one. Each column difference is zero-sum and has classical total variation at most 3
√
2s by POVM
contraction. For a zero-sum real vector v, its positive and negative parts each have mass ∥v∥1/2, so
∥v∥2
2 ≤2(∥v∥1/2)2. Each column thus has Euclidean norm at most 6√s, and
Technical passage, source page 10.
Searchable transcription
2s by POVM
contraction. For a zero-sum real vector v, its positive and negative parts each have mass ∥v∥1/2, so
∥v∥2
2 ≤2(∥v∥1/2)2. Each column thus has Euclidean norm at most 6√s, and
σ2 ≤∥C −C(1)∥op ≤∥C −C(1)∥F ≤18√s.
Source page 10.
Searchable transcription
σ2 ≤∥C −C(1)∥op ≤∥C −C(1)∥F ≤18√s.
These two cases exhaust the possible normalized factors.
If the right side of (24) exceeded e, choose s < σ2
2/324 sufficiently close to that upper endpoint.
Then both displayed inequalities would fail, since the limiting first bound is 6√e(1 + 18
√
Technical passage, source page 10.
Searchable transcription
These two cases exhaust the possible normalized factors.
If the right side of (24) exceeded e, choose s < σ2
2/324 sufficiently close to that upper endpoint.
Then both displayed inequalities would fail, since the limiting first bound is 6√e(1 + 18
√
6/σ2) < σ5.
This choice also satisfies s > e: column stochasticity gives σ2, σ5 ≤∥C∥F ≤3, and the assumed bound
on e is less than σ2
2/(24 · 324). This contradiction proves the necessary inequality.
□
Technical passage, source page 10.
Searchable transcription
6/σ2) < σ5.
This choice also satisfies s > e: column stochasticity gives σ2, σ5 ≤∥C∥F ≤3, and the assumed bound
on e is less than σ2
2/(24 · 324). This contradiction proves the necessary inequality.
□

A floor for every qutrit approximation of the trine record law

Let Pλ,4 be the specified four-symbol delayed trine law, with

Hab,cd = Pλ,4(a, b, c, d) = T λ
acT λ
bd
9
,
Cλ = 9H.
Source page 11.
Searchable transcription
Hab,cd = Pλ,4(a, b, c, d) = T λ
acT λ
bd
9
,
Cλ = 9H.
Suppose another complete four-record law P ′ is generated through a quantum causal cut of dimension
at most three. All history-correlated systems used again in the future belong to that cut; independent
fresh ports and an untouched retained archive obey the same convention as the accounting chapter.
The competing law may have a different update rule, lack stationarity, or lie outside the trine family.
Extra null/loss outcomes can remain in its complete alphabet, with zero target probability.
Put δ = TV(P ′, Pλ,4) < 1/9. Restrict its joint table H′ to the target trit pairs and let h′
cd =
P
ab H′
ab,cd ≥1/9 −δ. The column-normalized matrix C′
ab,cd = H′
ab,cd/h′
cd still has PSD rank at most
three. Define
Technical passage, source page 11.
Searchable transcription
Suppose another complete four-record law P ′ is generated through a quantum causal cut of dimension
at most three. All history-correlated systems used again in the future belong to that cut; independent
fresh ports and an untouched retained archive obey the same convention as the accounting chapter.
The competing law may have a different update rule, lack stationarity, or lie outside the trine family.
Extra null/loss outcomes can remain in its complete alphabet, with zero target probability.
Put δ = TV(P ′, Pλ,4) < 1/9. Restrict its joint table H′ to the target trit pairs and let h′
cd =
P
ab H′
ab,cd ≥1/9 −δ. The column-normalized matrix C′
ab,cd = H′
ab,cd/h′
cd still has PSD rank at most
three. Define
d =
36δ
1 −9δ,
a = 1 −λ
2
−d,
e∗=
27
1 −d
2λ/3 + 18δ
1 −9δ
.
Source page 11.
Searchable transcription
d =
36δ
1 −9δ,
a = 1 −λ
2
−d,
e∗=
27
1 −d
2λ/3 + 18δ
1 −9δ
.

Whenever d < 1, a > 0, the following necessary condition holds:

e∗≥
a4
36(a + 18
√
6)2 .
(25)
Source page 11, equation (25).
Searchable transcription
e∗≥
a4
36(a + 18
√
6)2 .
(25)

Thus failure of this scalar inequality excludes every qutrit-cut approximation at that tolerance, rather than merely a fitted trine member.

Derivation. Column normalization and total variation give

∥C′ −Cλ∥op ≤
X
ab,cd
|C′
ab,cd −(Cλ)ab,cd| ≤d.
Source page 11.
Searchable transcription
∥C′ −Cλ∥op ≤
X
ab,cd
|C′
ab,cd −(Cλ)ab,cd| ≤d.

Indeed the two contributions are at most the joint-table l1 error 2δ and the marginal l1 error 2δ, divided by 1/9 −δ. Thus row sums are at least 1 −d, and both relevant singular values are at least a. The unnormalized leakage functional takes values in [0, 2] and has target expectation 2λ/27. Its change is at most 2δ; dividing by the smallest column mass bounds L(C′) by (2λ/3 + 18δ)/(1 −9δ). Apply (24); its right side increases with each positive singular value. This gives (25). For explicit simple consequences, exact rational constants in the proof with s = 1/2000 give

TV(P ′, P0,4) > 10−9,
(26)
Source page 11, equation (26).
Searchable transcription
TV(P ′, P0,4) > 10−9,
(26)

and, uniformly for 0 ≤λ ≤10−8, TV(P ′, Pλ,4) > 5 × 10−10. (27)

At the first boundary the squared rank-four approximation error is below 0.216, and at the second below 0.188; the rank-one squared error is 0.162, while the required singular value exceeds 0.49 in both cases. These loose bounds are verified with fractions in the accompanying code.

What changes. The earlier exact dimension statement was robust only along a specified noise curve. The new witness supplies a finite-output accuracy obstruction against all ordinary qutrit-cut realizations, including time-dependent alternatives and complete null/loss records. It is a quantitative outer bound, not a sharp approximate-memory frontier, a large practical separation, or a resource theorem for nonlinear ontic assignments. Larger horizons inherit the same floor by retaining the four-record marginal. It compares the declared classical record law; it identifies no microscopic source or conscious perspective.

An explicit robust version of the support witness

The noiseless growing lower bound is sensitive to small forbidden-event probabilities. A quantitative
replacement is possible, with a poor large-grid margin that should be stated explicitly. Suppose a
row-stochastic matrix eT factors through r classical states as AB, with both A and B row stochastic.
Such a normalization is available for every nonnegative factorization of a stochastic matrix after zero
factors are removed. Set e = maxi eTii and μ = mini̸=j eTij. For τ > 0 put Si = {l: Ail≥τ}. If Si ⊆Sj,
then
eTij ≤e/τ + rτ.
(28)
Technical passage, source page 12.
Searchable transcription
The noiseless growing lower bound is sensitive to small forbidden-event probabilities. A quantitative
replacement is possible, with a poor large-grid margin that should be stated explicitly. Suppose a
row-stochastic matrix eT factors through r classical states as AB, with both A and B row stochastic.
Such a normalization is available for every nonnegative factorization of a stochastic matrix after zero
factors are removed. Set e = maxi eTii and μ = mini̸=j eTij. For τ > 0 put Si = {l: Ail≥τ}. If Si ⊆Sj,
then
eTij ≤e/τ + rτ.
(28)
Indeed, on Si the sum of Blj is at most e/τ because Ajl≥τ, and outside Si the contribution is at
most rτ. Choosing τ =
p
Technical passage, source page 12.
Searchable transcription
Indeed, on Si the sum of Blj is at most e/τ because Ajl≥τ, and outside Si the contribution is at
most rτ. Choosing τ =
p
e/r shows that μ > 2√re forces an antichain; the case e = 0 follows by a
limit or the exact support argument.
For an entrywise perturbation of the phase-circle transition by at most δ, write
Technical passage, source page 12.
Searchable transcription
e/r shows that μ > 2√re forces an antichain; the case e = 0 follows by a
limit or the exact support argument.
For an entrywise perturbation of the phase-circle transition by at most δ, write
μm = 1 −cos(2π/m)
m
= 2 sin2(π/m)
m
.
The sufficient antichain condition is μm −δ > 2
√
rδ, or equivalently
Source page 12.
Searchable transcription
μm = 1 −cos(2π/m)
m
= 2 sin2(π/m)
m
.
The sufficient antichain condition is μm −δ > 2
√
rδ, or equivalently
δ < δ∗
m,r :=
μ2
m
(√r + μm + √r)2 .
(29)
Source page 12, equation (29).
Searchable transcription
δ < δ∗
m,r :=
μ2
m
(√r + μm + √r)2 .
(29)

If a normalized observed triple distribution p′ lies at TV distance η from the exact law, with η < m−2, its conditional transition in any middle-symbol fibre obeys

p′(i, y, j)
P
k p′(i, y, k) −Tij
≤
2η
m−2 −η.
max
i,j
Source page 12.
Searchable transcription
p′(i, y, j)
P
k p′(i, y, k) −Tij
≤
2η
m−2 −η.
max
i,j

Both the single-entry and row-mass probability errors are at most η; the inequality follows by adding numerator and denominator errors. Therefore, choosing r = rm −1, the original classical sector lower bound mrm survives whenever

η <
δ∗
m,r
m2(2 + δ∗m,r).
(30)
Source page 12, equation (30).
Searchable transcription
η <
δ∗
m,r
m2(2 + δ∗m,r).
(30)

This is a rigorous finite-data statement, but its radius scales as roughly m−8/r and is extremely small at m = 180, 360, 720. It is not proposed as a practical large-grid experiment. The trine rank margin is substantially more usable.

Back to top ↑