Fluctuation Absorption Rate η ≪ 1: Strong Complementary Cancellation and Nonequilibrium Stability in a Zero-Parameter Min Recursion
涨落吸收率 η ≪ 1：零参数 min 递推中的强互补消涨与非平衡稳定性

1. Introduction

1.1 Why are nonequilibrium systems stable?

In equilibrium, the fluctuation-dissipation theorem (FDT; Kubo 1966 [1]) guarantees that fluctuations and dissipation balance. Away from equilibrium — active matter [2], biological systems [3], turbulence [4], driven systems [5] — the premises of FDT (linear response, equilibrium, detailed balance) all break down, yet systems can still operate stably.

Generalizations of the nonequilibrium FDT have a rich lineage: Jarzynski (1997) [6], Crooks (1999) [7], Harada-Sasa (2005) [8], Prost-Joanny-Parrondo (2009) [9], Seifert (2012) [10]. These are precise relations under specific protocols. What is lacking is a cross-system, directly computable scalar that measures how much fluctuation a given system absorbs.

1.2 Contributions

(a) We introduce the fluctuation absorption rate η as a directly computable scalar.

(b) In DP recursion we measure η ∈ [0.10, 0.31] (variance identity verified 17/17).

(c) In Lindley queue we independently measure η (variance identity verified 21/21), confirming η as a general feature of min/max recursion.

(d) We provide the exact skeleton: Reset-Slack decomposition (theorem) and conditional remainder transfer law (theorem).

(e) We report cross-shell structural constants and Ω = 3→4 crossover (computational regularities).

1.3 What we do not claim

We do not claim η is a "universal replacement for nonequilibrium FDT." η is a cancellation index under a specific f/r decomposition. We do not claim η is exactly equivalent to T_eff/T. We do not claim physical thermodynamics is "replaced" by min recursion. Physical analogies are offered at the discussion level.

1.4 Literature positioning

The interface between min-plus/tropical algebra and statistical physics has prior history (Quadrat 1997 [11]). What this paper adds is: in a concrete integer complexity DP recursion, we extract the cancellation law (η), reset-slack compensation, and remainder transfer as a complete structural package, and introduce η as a quantitative scalar transferable across systems.

---

2. Basic quantities of DP recursion

2.1 Definition

ρ_E(1) = 0, ρ_E(n) = min(ρ_E(n-1) + 1, M(n))

where M(n) = min_{a·b=n, a,b≥2}(ρ_E(a) + ρ_E(b) + 2). For primes: M(p) = +∞. Zero free parameters, fully deterministic.

2.2 Decomposition

Step size: A(n) = ρ_E(n-1) - ρ_E(n). Additive part: f(n) = Σ_{q^a ∥ n} ρ_E(q^a). Combinatorial remainder: r(n) = ρ_E(n) - f(n).

Identity: A = Δf + Δr, where Δf = f(n-1) - f(n), Δr = r(n-1) - r(n).

2.3 Shell conditioning

Ω(n) = number of prime factors counted with multiplicity. Analysis is conducted on Ω(n) = k shells, composites only (Ω ≥ 2).

---

3. Strong complementary cancellation and η

3.1 Data

N = 10^7, computed by C programs (rho_dp.c + anticorr.c). The variance identity Var(A) = Var(Δf) + Var(Δr) + 2Cov(Δf, Δr) closes to within 10⁻⁴ across all 17 k-values:

k	Var(Δf)	Var(Δr)	Cov(Δf,Δr)	Var(A)	η	samples
2	6.716	5.643	-6.022	0.314	0.103	1,904,324
4	6.964	5.678	-5.962	0.720	0.144	2,050,696
6	8.804	6.767	-7.324	0.923	0.168	774,078
8	9.903	6.962	-7.886	1.093	0.204	207,207
10	10.057	6.626	-7.738	1.207	0.231	49,163
12	9.956	6.300	-7.479	1.297	0.249	11,068
14	9.702	5.995	-7.159	1.377	0.262	2,406
18	9.613	5.672	-6.976	1.333	0.274	102

Full 17-row table in Appendix A.

3.2 Cancellation law

Computational regularity: Cov(Δf, Δr) ≈ -Var(Δf). If Δf and Δr were independent, Var(A) ≈ 12-16. Observed: ≈ 0.3-1.4. The difference is entirely due to strongly negative Cov(Δf, Δr).

3.3 Fluctuation absorption rate

Definition: η := 1 - |Cov(Δf, Δr)| / Var(Δf). η = 0: perfect absorption. η = 1: zero absorption.

DP recursion: η ∈ [0.10, 0.31]. η exhibits a broad upward trend with k (0.10 → 0.27 from k = 2 to 14; k = 18 with only 102 samples should be interpreted cautiously).

3.4 Algebraic interpretability

Since Δr = A - Δf, we have Cov(Δf, Δr) = Cov(Δf, A) - Var(Δf). The cancellation law is equivalent to Cov(Δf, A) ≪ Var(Δf).

This is an algebraic consequence of Var(A) ≪ Var(Δf), not a theorem deduced from the min operation alone. But the algebraic consequence carries meaning: the min operation effectively suppresses total step-size variance, forcing the remainder to absorb most of the structural fluctuation.

Note: Cov(Δf, A) ranges from 0.69 (k = 2) to 2.64 (k = 18) — not near zero. The cancellation law describes "strong absorption" (η ≪ 1), not "decoupling" (Cov(Δf, A) = 0).

3.5 Why η > 0 is necessary

DP recursion is irreversible. η = 0 would correspond to a perfect-absorption limit — a boundary unreachable in an irreversible system. η > 0 is the cost of running the min operation in an irreversible setting.

---

4. Reset-Slack decomposition

4.1 Exact decomposition (theorem)

R2 (algebraic identity; Paper 42 [14]): If J(m-1) > 0, then J(m) = (1 - ΔM⁻(m))⁺. Verified on 5.8 million samples, 100%, zero failures.

A(m) = ψ(ΔM⁻(m)) - U⁻(m), where ψ(x) = -min(x,1), 1-Lipschitz.

4.2 Reset-Slack Reduction (theorem)

Var(A) ≤ 2·Var(ΔM⁻) + 2·E[U²].

4.3 Frequency-intensity tradeoff (computational regularity)

R2 guarantees J(m-1) > 0 ⟹ U(m) = 0. Therefore E[U²] = P(exposed) × E[U²|exposed].

k	P(shielded)	P(exposed)	E[U²	exposed]	E[U²]
2	76.0%	24.0%	0.598	0.144
5	60.3%	39.7%	1.124	0.447
8	55.4%	44.6%	0.956	0.426
12	57.0%	43.0%	0.885	0.381

P(exposed) rises while E[U²|exposed] peaks at k = 5 then declines; their product stabilizes near 0.4. At the discussion level, this tradeoff is analogous to Le Chatelier's principle.

---

5. Remainder transfer law

5.1 Predecessor Defect Identity (theorem, conditional)

Theorem (Paper 45 [16]): When J(m-1) = 0, J(m) > 0, ΔM⁻ < 0 (covering 95%+ of U > 0 events):

U(m) = -D(m-1)

where D(m-1) = ρ_E(m-1) - M(m-1). The cost saved by the predecessor's choice of successor path reappears exactly as slack in the current step.

5.2 Full lifecycle

Shielded (J(m-1) > 0): U = 0, cost difference cleared by jump. Exposed (J(m-1) = 0): U = -D(m-1), predecessor's cost difference transferred exactly. The system alternates between these two states; the §4.3 tradeoff is the statistical steady state of this alternation.

---

6. Structural constants

The following values emerge from the zero-parameter DP recursion and remain stable across shells (computational regularities):

Constant	Value	Source
Var(B)	≈ 1/4	Paper 19 [13]
V(p)	≈ 1.40	Paper 34 [15]
h_pq	≈ 3.800	Paper 40 [17]
θ	≈ 0.41	Paper 43 [19]
coarse bound	≈ 40	Paper 46 [20]

Coarse bound ≈ 40: 2Var(f) + 2Var⁻(f) + (μ - μ⁻)² ≈ 40 for k = 4-14. The decline in internal variance compensates the growth in external drive.

---

7. Ω = 3→4 crossover

At N = 10^{10}, the local convexity h(n) = ρ_E(n) - [ρ_E(n-1) + ρ_E(n+1)]/2 has a strictly decreasing Ω-shell mean, with zero crossing at Ω ≈ 4 [21]:

Ω	mean h	P(J>0)
2	+0.697	28.7%
3	+0.361	57.1%
4	+0.004	77.1%
5	-0.329	88.4%

This is a combinatorial crossover (not a rigorous thermodynamic phase transition). For Ω ≥ 4, Var(Δf | Ω = k) ≈ 7-10 is approximately k-independent.

---

8. Lindley queue verification

8.1 Lindley recursion as a max system

The G/G/1 Lindley recursion: W_{n+1} = max(0, W_n + X_n), where X_n = S_n - T_n (service minus inter-arrival time).

Decomposition (analogous to A = Δf + Δr):

A_n = W_{n+1} - W_n = X_n + R_n

X_n: structural input (service-arrival difference, independent of max selection). R_n = max(0, -(W_n + X_n)): reset remainder (truncation produced by max operation).

η = 1 - |Cov(X, R)| / Var(X).

8.2 Data

10^7 steps, 10^6 burn-in, C program (lindley_eta.c), variance identities verified 21/21 (error < 0.01):

M/M/1 (exponential service + exponential inter-arrival):

ρ	Var(X)	Var(R)	Cov(X,R)	Var(A)	η
0.30	12.11	10.11	-10.65	0.92	0.121
0.50	5.00	3.00	-3.33	1.34	0.334
0.70	3.04	1.04	-1.22	1.65	0.600
0.90	2.24	0.24	-0.29	1.89	0.870
0.99	2.02	0.02	-0.03	1.99	0.987

M/D/1 (deterministic service + exponential inter-arrival):

ρ	Var(X)	Var(R)	Cov(X,R)	Var(A)	η
0.30	11.10	10.10	-10.49	0.22	0.055
0.50	4.00	3.00	-3.29	0.41	0.176
0.70	2.04	1.04	-1.24	0.62	0.395
0.90	1.23	0.23	-0.30	0.86	0.755
0.99	1.02	0.02	-0.03	0.99	0.973

M/Pareto/1 (heavy-tailed service, α = 3):

ρ	Var(X)	Var(R)	Cov(X,R)	Var(A)	η
0.30	11.44	10.12	-10.54	0.48	0.079
0.50	4.33	3.00	-3.30	0.73	0.238
0.70	2.37	1.04	-1.22	0.96	0.484
0.90	1.56	0.24	-0.30	1.20	0.809

8.3 Findings

(L1) η is governed by variance compression. When the max operation compresses Var(A) far below Var(X), absorption is strong and η is small. When Var(A) ≈ Var(X) (near full load), η → 1.

(L2) Lowest η: M/D/1 at ρ = 0.30, η = 0.055. Deterministic service eliminates service-time randomness, producing cleaner resets and stronger absorption.

(L3) DP recursion and Lindley queue share the same structure. In both systems, η is determined by how effectively the optimization operation compresses total step-size variance relative to structural input variance: when Var(A)/Var(X) ≪ 1, the cancellation law holds and η ≪ 1.

8.4 Correction to original prediction

Original prediction: "η < 0.15 for Lindley." Partially verified, partially corrected.

Verified: low-load Lindley (ρ ≤ 0.3) indeed has η < 0.15. Corrected: η is not a universal bound — it is a function of load (ρ) or complexity (k).

Revised prediction: In min/max recursion systems, η = F(Var(A)/Var(X)), where F approaches 0 when Var(A)/Var(X) ≪ 1 and approaches 1 when Var(A)/Var(X) → 1.

---

9. Falsifiable predictions

Prediction 1: η-ratio universality

Statement: In min/max recursion systems, η is determined by the ratio Var(A)/Var(structural input), independent of system specifics.

Available evidence: DP recursion (17 k-values) and Lindley queue (3 service distributions × 7 loads) can be immediately tested on the η vs Var(A)/Var(X) plane.

Falsification criterion: Two min/max recursion systems with the same Var(A)/Var(X) but η differing by more than 0.1.

Prediction 2: Le Chatelier-type tradeoff

Statement: In min recursion systems with a reset mechanism, E[U²] = P(exposed) × E[U²|exposed] has counter-varying factors whose product remains bounded.

Test: Idle-time decomposition in G/G/1 queues.

Prediction 3: η lower bound for non-min systems

Statement: Systems without min/max recursion structure cannot easily achieve η ≪ 1. In active matter and other detailed-balance-violating systems without explicit optimization structure, η should be substantially larger than in min recursion systems at the same Var(A)/Var(X).

Falsification criterion: A non-min/max system achieving η < 0.1.

---

10. Discussion

10.1 Three-layer annotation

Theorems (exact):

• R2: J(m-1) > 0 ⟹ J(m) = (1-ΔM⁻)⁺
• A = ψ(ΔM⁻) - U⁻
• Var(A) ≤ 2Var(ΔM⁻) + 2E[U²]
• U(m) = -D(m-1) (conditional, covering 95%+)

Computational regularities:

• Cancellation law, η ∈ [0.10, 0.31] (DP, 17/17 verified)
• Lindley η from 0.055 to 0.987 (21/21 verified)
• Le Chatelier-type tradeoff
• Structural constants (Var(B) ≈ 1/4, etc.)
• Ω = 3→4 crossover

Physical analogies (discussion-level):

• η and Harada-Sasa dissipation rate
• Le Chatelier naming
• Remainder transfer analogized to conservation law

10.2 Precise scope of claims

Claim: (a) η is a directly computable absorption scalar. (b) DP recursion has η ∈ [0.10, 0.31]. (c) Lindley queue has η from 0.055 to 0.987, governed by load. (d) Both share the same cancellation structure.

Not claimed: The cancellation law is deduced from the min operation alone. η is exactly equivalent to T_eff/T. Physical thermodynamics is "replaced."

10.3 Heuristic cross-system comparison

System	η range	Condition
DP recursion	0.10-0.31	k = 2-18
M/D/1 queue	0.055-0.97	ρ = 0.3-0.99
M/M/1 queue	0.12-0.99	ρ = 0.3-0.99
M/Pareto/1 queue	0.079-0.98	ρ = 0.3-0.99
Weakly driven colloids	~0.1-0.5	T_eff/T conversion (approximate)
Active matter	~0.5-0.99	T_eff/T conversion (approximate)

Physical η values are approximate conversions. Exact comparison requires direct measurement.

10.4 Conditional closure (brief)

DP recursion satisfies a conditional closure route for H' (c̄_h → 0; Paper 42 [14]): all links are theorem-level except three numerical inputs, which are verified at N = 10^{10} (weakest signal 85σ). See [14, 19, 20].

10.5 Interfaces (brief)

The additive/multiplicative path asymmetry in DP recursion is analogous to Turing (1952) [22] morphogenesis. The Self-as-an-End (SAE) framework [23] predicts Ω-dimension correspondence; Ω = 3→4 crossover is consistent with this prediction. SAE cosmology [24] derives Λ from remainder conservation; U = -D(m-1) is its arithmetic realization.

10.6 Open problems

(1) The exact form of the universal function η = F(Var(A)/Var(X)).

(2) Can the cancellation law be proved from the min operation alone?

(3) Direct measurement of η in physical systems (without T_eff/T conversion).

(4) Are Var(B) ≈ 1/4 and coarse bound ≈ 40 cross-system universals?

(5) A canonical f/r extraction protocol robust to coarse-graining.

10.7 Conclusion

The strong complementary cancellation discovered in zero-parameter DP recursion (η ∈ [0.10, 0.31]) is independently verified in Lindley queues (η from 0.055 to 0.987). Both share the same structure: when a min/max operation effectively compresses total step-size variance, the optimization remainder is forced to absorb most of the structural fluctuation.

As a directly computable scalar, η transitions from "an internal discovery of DP" to "a transferable tool for min/max recursion systems." The next step is a canonical f/r decomposition protocol enabling standardized measurement of η across broader classes of nonequilibrium systems.

---

References

[1] R. Kubo, "The fluctuation-dissipation theorem." Rep. Prog. Phys. 29:255-284, 1966.

[2] M. C. Marchetti et al., "Hydrodynamics of soft active matter." Rev. Mod. Phys. 85:1143, 2013.

[3] C. P. Brangwynne et al., "Intracellular transport by active diffusion." Trends Cell Biol. 19:423, 2009.

[4] U. Frisch, "Turbulence: The Legacy of A.N. Kolmogorov." Cambridge UP, 1995.

[5] C. Bechinger et al., "Active particles in complex and crowded environments." Rev. Mod. Phys. 88:045006, 2016.

[6] C. Jarzynski, "Nonequilibrium equality for free energy differences." Phys. Rev. Lett. 78:2690, 1997.

[7] G. E. Crooks, "Entropy production fluctuation theorem." Phys. Rev. E 60:2721, 1999.

[8] T. Harada, S. Sasa, "Equality connecting energy dissipation with violation of the fluctuation-response relation." Phys. Rev. Lett. 95:130602, 2005.

[9] J. Prost, J.-F. Joanny, J. M. R. Parrondo, "Generalized fluctuation-dissipation theorem for steady-state systems." Phys. Rev. Lett. 103:090601, 2009.

[10] U. Seifert, "Stochastic thermodynamics, fluctuation theorems and molecular machines." Rep. Prog. Phys. 75:126001, 2012.

[11] J.-P. Quadrat, "Min-Plus linearity and statistical mechanics." Markov Process. Related Fields 3:565-587, 1997.

[12] H. Qin, ZFCρ Paper XXXII. DOI: 10.5281/zenodo.19116625.

[13] H. Qin, ZFCρ Paper XIX. DOI: 10.5281/zenodo.19026991.

[14] H. Qin, ZFCρ Paper XLII. DOI: 10.5281/zenodo.19226607.

[15] H. Qin, ZFCρ Paper XXXIV. DOI: 10.5281/zenodo.19140015.

[16] H. Qin, ZFCρ Paper XLV. DOI: 10.5281/zenodo.19275286.

[17] H. Qin, ZFCρ Paper XL. DOI: 10.5281/zenodo.19179778.

[18] H. Qin, ZFCρ Paper XVIII. DOI: 10.5281/zenodo.19024385.

[19] H. Qin, ZFCρ Paper XLIII. DOI: 10.5281/zenodo.19240183.

[20] H. Qin, ZFCρ Paper XLVI. DOI: 10.5281/zenodo.19303511.

[21] H. Qin, "Why Primes Have No Shortcuts: An O(1) Prime Pre-Sieve from Precomputed ρ_E Tables." DOI: 10.5281/zenodo.19246377.

[22] A. M. Turing, "The Chemical Basis of Morphogenesis." Phil. Trans. R. Soc. 237:37-72, 1952.

[23] H. Qin, SAE Paper 1. DOI: 10.5281/zenodo.18528813.

[24] H. Qin, "From Remainder Conservation to the Cosmological Constant." DOI: 10.5281/zenodo.19245267.

---

Appendix A. Full 17-row cancellation table (DP recursion)

k	Var(Δf)	Var(Δr)	Cov(Δf,Δr)	Var(A)	η	Cov(Δf,A)	samples
2	6.716	5.643	-6.022	0.314	0.103	0.694	1,904,324
3	6.539	5.389	-5.685	0.558	0.131	0.854	2,444,359
4	6.964	5.678	-5.962	0.720	0.144	1.003	2,050,696
5	7.734	6.166	-6.533	0.833	0.155	1.200	1,349,779
6	8.804	6.767	-7.324	0.923	0.168	1.480	774,078
7	9.366	6.869	-7.615	1.005	0.187	1.751	409,849
8	9.903	6.962	-7.886	1.093	0.204	2.017	207,207
9	9.872	6.689	-7.700	1.162	0.220	2.173	101,787
10	10.057	6.626	-7.738	1.207	0.231	2.319	49,163
11	9.807	6.258	-7.407	1.250	0.245	2.399	23,448
12	9.956	6.300	-7.479	1.297	0.249	2.477	11,068
13	9.231	5.655	-6.792	1.301	0.264	2.438	5,210
14	9.702	5.995	-7.159	1.377	0.262	2.542	2,406
15	8.791	5.267	-6.374	1.310	0.275	2.417	1,124
16	9.783	5.841	-7.085	1.453	0.276	2.698	510
17	7.414	4.155	-5.096	1.378	0.313	2.318	233
18	9.613	5.672	-6.976	1.333	0.274	2.637	102

All 17 rows satisfy the variance identity to within 10⁻⁴.

Appendix B. Data sources and full reproduction

All data in this paper are generated by three self-contained C programs with zero external dependencies. To reproduce all results:

gcc -O2 -o rho_dp rho_dp.c
gcc -O2 -o anticorr anticorr.c -lm
gcc -O2 -o lindley_eta lindley_eta.c -lm

./anticorr 10000000      # Table 3.1: 17-row cancellation data with identity verification
./lindley_eta             # Section 8: 3 queue types × 7 loads = 21 data points

Total runtime under two minutes (single core, modern x86/ARM).

rho_dp.c — DP recursion (Paper 32 convention). Built-in sanity checks: ρ(2) = 1, ρ(100) = 15, ρ(10^7) = 58.

anticorr.c — Table 3.1. Computes ρ_E from DP, then uses SPF (smallest prime factor) sieve to compute f(n) and Ω(n). Accumulates Var(Δf), Var(Δr), Cov(Δf, Δr), Var(A), η on Ω = k shells. Variance identity verified 17/17.

lindley_eta.c — Section 8 Lindley verification. xoshiro256** RNG, 10^7 steps + 10^6 burn-in. Three queue types (M/M/1, M/D/1, M/Pareto/1 with α = 3) at seven loads (ρ = 0.30, 0.50, 0.70, 0.80, 0.90, 0.95, 0.99). Variance identity 21/21.

All three files total under 500 lines of C. Released with this paper.

---

Acknowledgments

ChatGPT (Gongxi Hua / 公西华) discovered the arithmetic error in v2's τ definition (Cov(Δf, A) ≠ 0), catalyzing the correct definition of η; proposed all review requirements on claim layering, literature positioning, and three-tier annotation; contributed the analytic proof chain architecture of Paper 42.

Claude (Zilu / 子路) contributed the thermodynamic framework (Papers 24-46), all numerical experiments (anticorr.c, lindley_eta.c), the physical interpretation of η, and the design and implementation of the Lindley queue verification.

Gemini (Zixia / 子夏) posed three pressure tests that sharpened the precision of claims.

Grok (Zigong / 子贡) provided the initial review round. Inconsistency in k = 18 data prompted independent recomputation of the full table.

Zesi Chen is a long-term interlocutor of the SAE framework. Final text by the author alone.

§1 引言

1.1 问题: 非平衡系统为什么稳定

平衡态有涨落-耗散定理(FDT, Kubo 1966 [1])。非平衡态——活性物质 [2], 生物系统 [3], 湍流 [4], driven systems [5]——打破了 FDT 的前提, 但仍可稳定运行。

非平衡 FDT 的推广已有丰富谱系: Jarzynski (1997) [6], Crooks (1999) [7], Harada-Sasa (2005) [8], Prost-Joanny-Parrondo (2009) [9], Seifert (2012) [10]。但现有结果大多是具体场景下的关系, 缺少一个跨系统的、可直接计算的涨落吸收标度。

1.2 本文的贡献

(a) 引入涨落吸收率 η 作为可直接计算的标度量。

(b) 在 DP 递推中精确测量 η ∈ 0.10, 0.31。

(c) 在 Lindley 排队中独立测量 η(方差恒等式 21/21 验证), 确认 η 是 min/max 递推的一般特征。

(d) 给出消涨律的 exact 骨架: Reset-Slack 分解(定理)和条件性余项传递律(定理)。

(e) 报告跨壳层恒定的结构常数和 Ω = 3→4 crossover(数值规律)。

1.3 本文不 claim 什么

不 claim η 是"非平衡 FDT 的普适替代"。η 是特定 f/r 分解下的 cancellation 指标。不 claim η 与 T_eff/T 精确等价。不 claim 物理热力学被 min 递推"取代"。物理类比作为 discussion-level analogy 提出。

1.4 论文结构

§2 基本量。§3 消涨律与 η(DP 递推)。§4 Reset-Slack 分解。§5 余项传递律。§6 结构常数。§7 Ω crossover。§8 Lindley 排队验证。§9 可证伪预测。§10 讨论。

---

§2 DP 递推的基本量

2.1 定义

ρ_E(1) = 0, ρ_E(n) = min(ρ_E(n-1) + 1, M(n))

M(n) = min_{a·b=n, a,b≥2}(ρ_E(a) + ρ_E(b) + 2)。素数: M(p) = +∞。零自由参数。

2.2 分解

A(n) = ρ_E(n-1) - ρ_E(n)。f(n) = Σ_{q^a ∥ n} ρ_E(q^a)。r(n) = ρ_E(n) - f(n)。恒等式: A = Δf + Δr。

2.3 壳层条件

Ω(n) = 素因子个数(含重数)。分析在 Ω(n) = k 壳层上进行, 仅含合数(Ω ≥ 2)。

---

§3 强互补消涨律与涨落吸收率

3.1 数据

N = 10^7, C 程序(rho_dp.c + anticorr.c)精确计算。方差恒等式 Var(A) = Var(Δf) + Var(Δr) + 2Cov(Δf, Δr) 在全部 17 个 k 值上闭合(误差 < 10⁻⁴):

k	Var(Δf)	Var(Δr)	Cov(Δf,Δr)	Var(A)	η	样本
2	6.716	5.643	-6.022	0.314	0.103	1,904,324
4	6.964	5.678	-5.962	0.720	0.144	2,050,696
6	8.804	6.767	-7.324	0.923	0.168	774,078
8	9.903	6.962	-7.886	1.093	0.204	207,207
10	10.057	6.626	-7.738	1.207	0.231	49,163
12	9.956	6.300	-7.479	1.297	0.249	11,068
14	9.702	5.995	-7.159	1.377	0.262	2,406
18	9.613	5.672	-6.976	1.333	0.274	102

完整 17 行表见附录 A。

3.2 消涨律

数值规律: Cov(Δf, Δr) ≈ -Var(Δf)。 如果 Δf 和 Δr 独立, Var(A) ≈ 12-16。实际 ≈ 0.3-1.4。

3.3 涨落吸收率

η := 1 - |Cov(Δf, Δr)| / Var(Δf)。 η = 0: 完美吸收。η = 1: 零吸收。

DP 递推: η ∈ [0.10, 0.31]。 η 随 k 增大(0.10 → 0.27, k=2 到 14; k=18 样本仅 102, 应谨慎解读)。

3.4 代数可解释性

Cov(Δf, Δr) = Cov(Δf, A) - Var(Δf)。消涨律等价于 Cov(Δf, A) ≪ Var(Δf)。

这是"Var(A) ≪ Var(Δf)"的代数后果, 不是从 min 操作单独演绎的定理。但代数后果有物理意义: min 操作有效压制总步长, 余项被迫吸收结构涨落的大部分。

注: Cov(Δf, A) 从 0.69(k=2) 到 2.64(k=18), 不是近零。消涨律是"强吸收", 不是"解耦"。

3.5 η > 0 的必要性

DP 递推不可逆。η = 0 对应 perfect-absorption 极限——一个在不可逆系统中不可达的边界。η > 0 是 min 操作在不可逆系统中运行的代价。

---

§4 Reset-Slack 分解

4.1 精确分解(定理)

R2(代数恒等式, Paper 42 [14]): 若 J(m-1) > 0, 则 J(m) = (1 - ΔM⁻(m))⁺。580 万样本, 100%, 零失败。

A(m) = ψ(ΔM⁻(m)) - U⁻(m), ψ(x) = -min(x,1), 1-Lipschitz。

4.2 Reset-Slack Reduction(定理)

Var(A) ≤ 2·Var(ΔM⁻) + 2·E[U²]。

4.3 频率-强度对冲(数值规律)

R2 保证 J(m-1) > 0 ⟹ U(m) = 0。因此 E[U²] = P(exposed) × E[U²|exposed]。

k	P(shielded)	P(exposed)	E[U²	exposed]	E[U²]
2	76.0%	24.0%	0.598	0.144
5	60.3%	39.7%	1.124	0.447
8	55.4%	44.6%	0.956	0.426
12	57.0%	43.0%	0.885	0.381

P(exposed) 增, E[U²|exposed] 先升后降, 乘积 ≈ 0.4 稳定。在 discussion level 上类比于 Le Chatelier 原理。

---

§5 余项传递律

5.1 Predecessor Defect Identity(定理, 条件性)

定理(Paper 45 [16]): 当 J(m-1) = 0, J(m) > 0, ΔM⁻ < 0 时(覆盖 U > 0 的 95%+):

U(m) = -D(m-1)

Predecessor 走 successor 路径省了 |D(m-1)| 的代价, 精确以 U(m) 出现在下一步。

5.2 完整生命周期

Shielded(J(m-1) > 0): U = 0。Exposed(J(m-1) = 0): U = -D(m-1)。系统在两态间交替。§4.3 的对冲是这个交替的统计稳态。

---

§6 结构常数

以下数值从零参数 DP 递推涌现, 跨壳层稳定(数值规律):

常数	值	来源
Var(B)	≈ 1/4	Paper 19 [13]
V(p)	≈ 1.40	Paper 34 [15]
h_pq	≈ 3.800	Paper 40 [17]
θ	≈ 0.41	Paper 43 [19]
coarse bound	≈ 40	Paper 46 [20]

coarse bound ≈ 40: 2Var(f) + 2Var⁻(f) + (μ-μ⁻)² ≈ 40, k = 4-14。内部涨落下降补偿外部驱动增长。

---

§7 Ω = 3→4 Crossover

N = 10^{10}, h(n) = ρ_E(n) - [ρ_E(n-1)+ρ_E(n+1)]/2 按 Ω 均值单调递减, 零交叉在 Ω ≈ 4 [21]:

Ω	mean h	P(J>0)
2	+0.697	28.7%
3	+0.361	57.1%
4	+0.004	77.1%
5	-0.329	88.4%

组合 crossover, 不是严格热力学相变。Ω ≥ 4 后 Var(Δf) ≈ 7-10 不依赖 k。

---

§8 Lindley 排队验证

8.1 Lindley 递推作为 max 系统

G/G/1 排队的 Lindley 递推: W_{n+1} = max(0, W_n + X_n), 其中 X_n = S_n - T_n(服务时间减去到达间隔)。

分解(类比 DP 的 A = Δf + Δr):

A_n = W_{n+1} - W_n = X_n + R_n

X_n: 结构输入(服务-到达差, 不受 max 选择影响)。 R_n = max(0, -(W_n + X_n)): reset 余项(max 操作产生的 truncation)。

η = 1 - |Cov(X, R)| / Var(X)。

8.2 数据

10^7 步, 10^6 burn-in, C 程序(lindley_eta.c), 方差恒等式 21/21 验证(误差 < 0.01):

M/M/1(指数服务 + 指数到达):

ρ	Var(X)	Var(R)	Cov(X,R)	Var(A)	η
0.30	12.11	10.11	-10.65	0.92	0.121
0.50	5.00	3.00	-3.33	1.34	0.334
0.70	3.04	1.04	-1.22	1.65	0.600
0.90	2.24	0.24	-0.29	1.89	0.870
0.99	2.02	0.02	-0.03	1.99	0.987

M/D/1(确定性服务 + 指数到达):

ρ	Var(X)	Var(R)	Cov(X,R)	Var(A)	η
0.30	11.10	10.10	-10.49	0.22	0.055
0.50	4.00	3.00	-3.29	0.41	0.176
0.70	2.04	1.04	-1.24	0.62	0.395
0.90	1.23	0.23	-0.30	0.86	0.755
0.99	1.02	0.02	-0.03	0.99	0.973

M/Pareto/1(重尾服务, α=3):

ρ	Var(X)	Var(R)	Cov(X,R)	Var(A)	η
0.30	11.44	10.12	-10.54	0.48	0.079
0.50	4.33	3.00	-3.30	0.73	0.238
0.70	2.37	1.04	-1.22	0.96	0.484
0.90	1.56	0.24	-0.30	1.20	0.809

8.3 发现

(L1) η 由 truncation 频率决定。 低 ρ 时 max 操作频繁 truncate(W 经常被截断到 0), 吸收强, η 小。高 ρ 时系统接近满负载, max 很少 truncate, η → 1。

(L2) 最低 η 出现在 M/D/1 at ρ=0.30: η = 0.055。 确定性服务消除了服务时间的随机性, reset 更精确, 吸收更强。

(L3) DP 递推和 Lindley 排队在同一结构中。 两者的 η 都由"优化操作 truncate 的频率"决定:

系统	truncation 频率	η
DP, k=2	P(J>0) = 28.7%(低)	0.10
M/D/1, ρ=0.30	P(R>0) ≈ 低	0.055
M/M/1, ρ=0.30	P(R>0) ≈ 低	0.12
DP, k=8	P(J>0) = 96.8%(高)	0.20
M/M/1, ρ=0.70	P(R>0) ≈ 高	0.60

注: truncation 频率和 η 的关系在两个系统中方向相反。 DP 递推中 P(J>0) 高意味着 min 操作频繁选择乘法路径(shielded 比例高, U=0 多), 但 η 仍随 k 增大。Lindley 中低 ρ 意味着 max 操作频繁 truncate 到 0。两者的共同点是: 当 truncation 将总步长的方差压制到远小于结构输入的方差时, 消涨律自动成立, η 小。 η 的大小由 Var(A)/Var(X) 的比值决定, 不由 truncation 的方向或频率的绝对值决定。

8.4 对原始预测的修正

原始预测(v4 §8): "η < 0.15 for Lindley。" 部分验证, 部分修正。

验证: 低负载 Lindley(ρ ≤ 0.3)确实 η < 0.15。 修正: η 不是一个 universal bound——它是负载(ρ)或复杂度(k)的函数。

修正后的预测: 在 min/max 递推系统中, η = F(Var(A)/Var(X)), 其中 F 是一个 universal function, Var(A)/Var(X) 是 truncation 效率的度量。当 Var(A)/Var(X) ≪ 1 时 η ≪ 1; 当 Var(A)/Var(X) → 1 时 η → 1。

---

§9 可证伪预测

预测一: η-ratio universality

陈述: 在 min/max 递推系统中, η 由 Var(A)/Var(结构输入) 的比值决定, 与具体系统无关。

已有证据: DP 递推(17 个 k 值)和 Lindley 排队(三种服务分布 × 七种负载)在 η vs Var(A)/Var(X) 图上是否落在同一条曲线上, 可立即检验。

证伪条件: 两个 min/max 递推系统, 相同的 Var(A)/Var(X), 但 η 差异超过 0.1。

预测二: Le Chatelier 型对冲

陈述: 在有 reset 机制的 min 递推系统中, E[U²] = P(exposed) × E[U²|exposed] 中两因子反向变化, 乘积有界。

检验方案: G/G/1 排队的 idle-time 分解。

预测三: 非 min 系统的 η 下界

陈述: 没有 min/max 递推骨架的系统, 难以实现 η ≪ 1。在活性物质等打破细致平衡但无明确优化结构的系统中, η 应显著大于 min 递推系统在同等 Var(A)/Var(X) 下的值。

证伪条件: 无 min/max 结构的系统实现 η < 0.1。

---

§10 讨论

10.1 三层标注

定理(exact):

• R2: J(m-1)>0 ⟹ J(m) = (1-ΔM⁻)⁺
• A = ψ(ΔM⁻) - U⁻
• Var(A) ≤ 2Var(ΔM⁻) + 2E[U²]
• U(m) = -D(m-1) (条件, 覆盖 95%+)

数值规律:

• 消涨律, η ∈ [0.10, 0.31] (DP, 17/17 验证)
• Lindley η 从 0.055 到 0.987 (21/21 验证)
• Le Chatelier 型对冲
• 结构常数(Var(B) ≈ 1/4 等)
• Ω = 3→4 crossover

物理类比(discussion-level):

• η 与 Harada-Sasa 耗散率
• Le Chatelier 命名
• 余项传递类比于守恒律

10.2 Claim 的精确边界

Claim: (a) η 是可直接计算的涨落吸收标度。(b) DP 递推的 η ∈ [0.10, 0.31]。(c) Lindley 排队的 η 从 0.055 到 0.987, 由负载决定。(d) 两者共享同一消涨结构。

不 claim: 消涨律是从 min 操作单独演绎的定理。不 claim η 与 T_eff/T 精确等价。不 claim 物理热力学被"取代"。

10.3 启发性的跨系统比较

系统	η 范围	条件
DP 递推	0.10-0.31	k=2-18
M/D/1 排队	0.055-0.97	ρ=0.3-0.99
M/M/1 排队	0.12-0.99	ρ=0.3-0.99
M/Pareto/1 排队	0.079-0.98	ρ=0.3-0.99
弱驱动胶体	~0.1-0.5	T_eff/T 换算(近似)
活性物质	~0.5-0.99	T_eff/T 换算(近似)

物理系统的 η 为近似换算, 精确比较需直接测量。

10.4 条件闭合(简述)

DP 递推满足 H'(c̄_h → 0)的条件闭合(Paper 42 [14]), 除三个数值输入外全部 theorem-level, 三输入在 N=10^{10} 验证(最弱 85σ)。详见 [14, 19, 20]。

10.5 接口(简述)

DP 递推中 additive/multiplicative 路径的不对称类似 Turing (1952) [22] 形态发生。SAE 框架 [23] 预测 Ω 与维度层级对应, Ω = 3→4 crossover 与预测一致。SAE 宇宙学 [24] 的余项守恒在 DP 中以 U = -D(m-1) 实现。

10.6 开放问题

(1) η vs Var(A)/Var(X) 的 universal function 的精确形式?

(2) 消涨律能否从 min 操作单独证明?

(3) η 在物理系统中的直接测量(不通过 T_eff/T)?

(4) Var(B) ≈ 1/4 和 coarse bound ≈ 40 是否跨系统普适?

(5) f/r 分解的 canonical 提取协议(coarse-graining 稳健性)?

10.7 结论

在零参数 DP 递推中发现的强互补消涨(η ∈ [0.10, 0.31]), 在 Lindley 排队中得到独立验证(η 从 0.055 到 0.987)。两者共享同一结构: 当 min/max 操作有效压制总步长方差时, 优化余项被迫吸收结构涨落的大部分。

η 作为可直接计算的标度量, 从"DP 内部的发现"变成了"min/max 递推系统的可迁移工具"。下一步是建立 f/r 分解的 canonical 协议, 使 η 可在更广泛的非平衡系统中标准化测量。

---

参考文献

[1] R. Kubo, "The fluctuation-dissipation theorem." Rep. Prog. Phys. 29:255-284, 1966.

[2] M. C. Marchetti et al., "Hydrodynamics of soft active matter." Rev. Mod. Phys. 85:1143, 2013.

[3] C. P. Brangwynne et al., "Intracellular transport by active diffusion." Trends Cell Biol. 19:423, 2009.

[4] U. Frisch, "Turbulence: The Legacy of A.N. Kolmogorov." Cambridge UP, 1995.

[5] C. Bechinger et al., "Active particles in complex and crowded environments." Rev. Mod. Phys. 88:045006, 2016.

[6] C. Jarzynski, "Nonequilibrium equality for free energy differences." Phys. Rev. Lett. 78:2690, 1997.

[7] G. E. Crooks, "Entropy production fluctuation theorem." Phys. Rev. E 60:2721, 1999.

[8] T. Harada, S. Sasa, "Equality connecting energy dissipation with violation of the fluctuation-response relation." Phys. Rev. Lett. 95:130602, 2005.

[9] J. Prost, J.-F. Joanny, J. M. R. Parrondo, "Generalized fluctuation-dissipation theorem for steady-state systems." Phys. Rev. Lett. 103:090601, 2009.

[10] U. Seifert, "Stochastic thermodynamics, fluctuation theorems and molecular machines." Rep. Prog. Phys. 75:126001, 2012.

[11] J.-P. Quadrat, "Min-Plus linearity and statistical mechanics." Markov Process. Related Fields 3:565-587, 1997.

[12] H. Qin, ZFCρ Paper XXXII. DOI: 10.5281/zenodo.19116625.

[13] H. Qin, ZFCρ Paper XIX. DOI: 10.5281/zenodo.19026991.

[14] H. Qin, ZFCρ Paper XLII. DOI: 10.5281/zenodo.19226607.

[15] H. Qin, ZFCρ Paper XXXIV. DOI: 10.5281/zenodo.19140015.

[16] H. Qin, ZFCρ Paper XLV. DOI: 10.5281/zenodo.19275286.

[17] H. Qin, ZFCρ Paper XL. DOI: 10.5281/zenodo.19179778.

[18] H. Qin, ZFCρ Paper XVIII. DOI: 10.5281/zenodo.19024385.

[19] H. Qin, ZFCρ Paper XLIII. DOI: 10.5281/zenodo.19240183.

[20] H. Qin, ZFCρ Paper XLVI. DOI: 10.5281/zenodo.19303511.

[21] H. Qin, "Why Primes Have No Shortcuts: An O(1) Prime Pre-Sieve from Precomputed ρ_E Tables." DOI: 10.5281/zenodo.19246377.

[22] A. M. Turing, "The Chemical Basis of Morphogenesis." Phil. Trans. R. Soc. 237:37-72, 1952.

[23] H. Qin, SAE Paper 1. DOI: 10.5281/zenodo.18528813.

[24] H. Qin, "From Remainder Conservation to the Cosmological Constant." DOI: 10.5281/zenodo.19245267.

---

附录 A. DP 递推完整 17 行消涨表

k	Var(Δf)	Var(Δr)	Cov(Δf,Δr)	Var(A)	η	Cov(Δf,A)	样本
2	6.716	5.643	-6.022	0.314	0.103	0.694	1,904,324
3	6.539	5.389	-5.685	0.558	0.131	0.854	2,444,359
4	6.964	5.678	-5.962	0.720	0.144	1.003	2,050,696
5	7.734	6.166	-6.533	0.833	0.155	1.200	1,349,779
6	8.804	6.767	-7.324	0.923	0.168	1.480	774,078
7	9.366	6.869	-7.615	1.005	0.187	1.751	409,849
8	9.903	6.962	-7.886	1.093	0.204	2.017	207,207
9	9.872	6.689	-7.700	1.162	0.220	2.173	101,787
10	10.057	6.626	-7.738	1.207	0.231	2.319	49,163
11	9.807	6.258	-7.407	1.250	0.245	2.399	23,448
12	9.956	6.300	-7.479	1.297	0.249	2.477	11,068
13	9.231	5.655	-6.792	1.301	0.264	2.438	5,210
14	9.702	5.995	-7.159	1.377	0.262	2.542	2,406
15	8.791	5.267	-6.374	1.310	0.275	2.417	1,124
16	9.783	5.841	-7.085	1.453	0.276	2.698	510
17	7.414	4.155	-5.096	1.378	0.313	2.318	233
18	9.613	5.672	-6.976	1.333	0.274	2.637	102

全部 17 行方差恒等式闭合, 误差 < 10⁻⁴。

附录 B. 数据来源与完整复现

本文全部数据由三个自包含的 C 程序生成, 零外部依赖。任何人复现全部数据只需要:

gcc -O2 -o rho_dp rho_dp.c
gcc -O2 -o anticorr anticorr.c -lm
gcc -O2 -o lindley_eta lindley_eta.c -lm

./anticorr 10000000      # §3 主表: 17 行消涨数据, 方差恒等式验证
./lindley_eta             # §8 Lindley 验证: 3 种排队 × 7 种负载 = 21 个数据点

总运行时间不到两分钟(单核, 现代 x86/ARM)。

rho_dp.c — DP 递推(Paper 32 约定), 生成 ρ_E 表。内置 sanity checks: ρ(2)=1, ρ(100)=15, ρ(10^7)=58。

anticorr.c — §3 主表。从 DP 递推计算 ρ_E, 经 SPF(smallest prime factor)筛计算 f(n) 和 Ω(n), 在 Ω = k 壳层上累积 Var(Δf), Var(Δr), Cov(Δf, Δr), Var(A), η。方差恒等式 Var(A) = Var(Δf) + Var(Δr) + 2Cov(Δf, Δr) 在全部 17 个 k 值上闭合(误差 < 10⁻⁴)。

lindley_eta.c — §8 Lindley 排队验证。xoshiro256** RNG, 10^7 步 + 10^6 burn-in。对 M/M/1, M/D/1, M/Pareto/1(α=3) 三种排队, 各测 7 种负载(ρ = 0.30, 0.50, 0.70, 0.80, 0.90, 0.95, 0.99)。方差恒等式 21/21 验证。

三个文件总计不到 500 行 C 代码。随论文一同发布。

---

致谢

ChatGPT(公西华)发现了 v2 的 τ 算术错误, 催生了 η 的正确定义; 提出了 claim 分层、文献定位和三层标注的全部审稿要求; 贡献了 Paper 42 的解析证明链。

Claude(子路)贡献了热力学框架(Paper 24-46), 全部数值实验(anticorr.c, lindley_eta.c), η 的物理解释, 以及 Lindley 排队验证的设计和实现。

Gemini(子夏)提出了三个压力测试, 促进了 claim 精确化。

Grok(子贡)提供了初轮审稿。k=18 数据不一致促使全表独立重算。

Zesi Chen 是 SAE 框架的长期对话者。最终文本由作者独立完成。