ZFCρ Series · Thermodynamic Interface

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

Abstract

We discover and quantify a strong fluctuation absorption mechanism in a zero-parameter arithmetic recursion — the integer complexity dynamic programming (DP) ρ_E(n) = min(ρ_E(n-1)+1, M(n)).

Decomposing the step size A = ρ_E(n-1) - ρ_E(n) into structural input Δf and optimization remainder Δr, the system satisfies a strong complementary cancellation law: Cov(Δf, Δr) ≈ -Var(Δf). Two quantities with individual variances ≈ 7-10 cancel precisely to produce a total variance ≈ 0.3-1.3. We define the fluctuation absorption rate:

η := 1 - |Cov(Δf, Δr)| / Var(Δf)

At N = 10^7 across Ω = 2-18 shells: η ∈ [0.10, 0.31], meaning 69-90% of structural fluctuations are absorbed by the optimization remainder.

The exact skeleton of the cancellation law rests on two theorems: (1) Reset-Slack decomposition A = ψ(ΔM⁻) - U⁻, yielding Var(A) ≤ 2Var(ΔM⁻) + 2E[U²]; (2) conditional remainder transfer law U(m) = -D(m-1) (covering 95%+ of U > 0 events).

First external verification: Lindley queue (G/G/1). For M/M/1, M/D/1, and M/Pareto/1 queues, η ranges from 0.055 (M/D/1 at ρ = 0.3) to 0.987 (M/M/1 at ρ = 0.99), with variance identities verified 21/21. η is determined by how effectively the max operation compresses step-size variance relative to structural input variance. DP recursion and Lindley queue share the same cancellation structure — η is not a peculiarity of DP arithmetic, but a general feature of min/max recursion systems.

Keywords: fluctuation absorption, nonequilibrium stability, min recursion, max recursion, integer complexity, DP recursion, Lindley queue, strong cancellation, reset-slack decomposition, remainder transfer law.

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:

kVar(Δf)Var(Δr)Cov(Δf,Δr)Var(A)ηsamples
26.7165.643-6.0220.3140.1031,904,324
46.9645.678-5.9620.7200.1442,050,696
68.8046.767-7.3240.9230.168774,078
89.9036.962-7.8861.0930.204207,207
1010.0576.626-7.7381.2070.23149,163
129.9566.300-7.4791.2970.24911,068
149.7025.995-7.1591.3770.2622,406
189.6135.672-6.9761.3330.274102

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].

kP(shielded)P(exposed)E[U²|exposed]E[U²]
276.0%24.0%0.5980.144
560.3%39.7%1.1240.447
855.4%44.6%0.9560.426
1257.0%43.0%0.8850.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):

ConstantValueSource
Var(B)≈ 1/4Paper 19 [13]
V(p)≈ 1.40Paper 34 [15]
h_pq≈ 3.800Paper 40 [17]
θ≈ 0.41Paper 43 [19]
coarse bound≈ 40Paper 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 hP(J>0)
2+0.69728.7%
3+0.36157.1%
4+0.00477.1%
5-0.32988.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.3012.1110.11-10.650.920.121
0.505.003.00-3.331.340.334
0.703.041.04-1.221.650.600
0.902.240.24-0.291.890.870
0.992.020.02-0.031.990.987

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

ρVar(X)Var(R)Cov(X,R)Var(A)η
0.3011.1010.10-10.490.220.055
0.504.003.00-3.290.410.176
0.702.041.04-1.240.620.395
0.901.230.23-0.300.860.755
0.991.020.02-0.030.990.973

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

ρVar(X)Var(R)Cov(X,R)Var(A)η
0.3011.4410.12-10.540.480.079
0.504.333.00-3.300.730.238
0.702.371.04-1.220.960.484
0.901.560.24-0.301.200.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η rangeCondition
DP recursion0.10-0.31k = 2-18
M/D/1 queue0.055-0.97ρ = 0.3-0.99
M/M/1 queue0.12-0.99ρ = 0.3-0.99
M/Pareto/1 queue0.079-0.98ρ = 0.3-0.99
Weakly driven colloids~0.1-0.5T_eff/T conversion (approximate)
Active matter~0.5-0.99T_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)

kVar(Δf)Var(Δr)Cov(Δf,Δr)Var(A)ηCov(Δf,A)samples
26.7165.643-6.0220.3140.1030.6941,904,324
36.5395.389-5.6850.5580.1310.8542,444,359
46.9645.678-5.9620.7200.1441.0032,050,696
57.7346.166-6.5330.8330.1551.2001,349,779
68.8046.767-7.3240.9230.1681.480774,078
79.3666.869-7.6151.0050.1871.751409,849
89.9036.962-7.8861.0930.2042.017207,207
99.8726.689-7.7001.1620.2202.173101,787
1010.0576.626-7.7381.2070.2312.31949,163
119.8076.258-7.4071.2500.2452.39923,448
129.9566.300-7.4791.2970.2492.47711,068
139.2315.655-6.7921.3010.2642.4385,210
149.7025.995-7.1591.3770.2622.5422,406
158.7915.267-6.3741.3100.2752.4171,124
169.7835.841-7.0851.4530.2762.698510
177.4144.155-5.0961.3780.3132.318233
189.6135.672-6.9761.3330.2742.637102

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.

Writing Declaration: All text in this paper, both English and Chinese, is generated by the author Han Qin alone. The paper represents the author's sole responsibility.

摘要

我们在一个零参数的纯算术递推——整数复杂度的动态规划(DP) ρ_E(n) = min(ρ_E(n-1)+1, M(n))——中发现并量化了一条极强的涨落吸收机制。

将步长 A = ρ_E(n-1) - ρ_E(n) 分解为结构输入 Δf 和优化余项 Δr 后, 系统满足强互补消涨律: Cov(Δf, Δr) ≈ -Var(Δf)。两个各自方差约 7-10 的量精确对消到方差约 0.3-1.3。定义涨落吸收率:

η := 1 - |Cov(Δf, Δr)|/Var(Δf)

在 N = 10^7, Ω = 2-18 壳层上: η ∈ [0.10, 0.31], 即结构涨落的 69-90% 被优化余项吸收。

消涨律的 exact 骨架由两个定理支撑: (1) Reset-Slack 分解 A = ψ(ΔM⁻) - U⁻, 给出 Var(A) ≤ 2Var(ΔM⁻) + 2E[U²]; (2) 条件性余项传递律 U(m) = -D(m-1)(覆盖 U > 0 的 95%+)。

首个外部验证: Lindley 排队(G/G/1)。 对 M/M/1, M/D/1, M/Pareto/1 三种排队, η 从低负载时的 0.06(M/D/1, ρ=0.3) 到高负载时的 0.99(M/M/1, ρ=0.99), 方差恒等式 21/21 验证通过。η 由 max 操作的 truncation 频率决定: truncation 越频繁, 吸收越强, η 越小。DP 递推和 Lindley 排队在同一条"η vs truncation 频率"曲线上——η 不是 DP 的特殊算术性质, 而是 min/max 递推的一般特征。

关键词: 涨落吸收, 非平衡稳定性, min 递推, max 递推, 整数复杂度, DP 递推, Lindley 排队, 强互补消涨, Reset-Slack 分解, 余项传递律。

§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](方差恒等式 17/17 验证)。

(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⁻⁴):

kVar(Δf)Var(Δr)Cov(Δf,Δr)Var(A)η样本
26.7165.643-6.0220.3140.1031,904,324
46.9645.678-5.9620.7200.1442,050,696
68.8046.767-7.3240.9230.168774,078
89.9036.962-7.8861.0930.204207,207
1010.0576.626-7.7381.2070.23149,163
129.9566.300-7.4791.2970.24911,068
149.7025.995-7.1591.3770.2622,406
189.6135.672-6.9761.3330.274102

完整 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]。

kP(shielded)P(exposed)E[U²|exposed]E[U²]
276.0%24.0%0.5980.144
560.3%39.7%1.1240.447
855.4%44.6%0.9560.426
1257.0%43.0%0.8850.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/4Paper 19 [13]
V(p)≈ 1.40Paper 34 [15]
h_pq≈ 3.800Paper 40 [17]
θ≈ 0.41Paper 43 [19]
coarse bound≈ 40Paper 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 hP(J>0)
2+0.69728.7%
3+0.36157.1%
4+0.00477.1%
5-0.32988.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.3012.1110.11-10.650.920.121
0.505.003.00-3.331.340.334
0.703.041.04-1.221.650.600
0.902.240.24-0.291.890.870
0.992.020.02-0.031.990.987

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

ρVar(X)Var(R)Cov(X,R)Var(A)η
0.3011.1010.10-10.490.220.055
0.504.003.00-3.290.410.176
0.702.041.04-1.240.620.395
0.901.230.23-0.300.860.755
0.991.020.02-0.030.990.973

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

ρVar(X)Var(R)Cov(X,R)Var(A)η
0.3011.4410.12-10.540.480.079
0.504.333.00-3.300.730.238
0.702.371.04-1.220.960.484
0.901.560.24-0.301.200.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 将总步长的方差压制到远小于结构输入的方差时, 消涨律自动成立, η 小。η 的大小由 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.31k=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.5T_eff/T 换算(近似)
活性物质~0.5-0.99T_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 行消涨表

kVar(Δf)Var(Δr)Cov(Δf,Δr)Var(A)ηCov(Δf,A)样本
26.7165.643-6.0220.3140.1030.6941,904,324
36.5395.389-5.6850.5580.1310.8542,444,359
46.9645.678-5.9620.7200.1441.0032,050,696
57.7346.166-6.5330.8330.1551.2001,349,779
68.8046.767-7.3240.9230.1681.480774,078
79.3666.869-7.6151.0050.1871.751409,849
89.9036.962-7.8861.0930.2042.017207,207
99.8726.689-7.7001.1620.2202.173101,787
1010.0576.626-7.7381.2070.2312.31949,163
119.8076.258-7.4071.2500.2452.39923,448
129.9566.300-7.4791.2970.2492.47711,068
139.2315.655-6.7921.3010.2642.4385,210
149.7025.995-7.1591.3770.2622.5422,406
158.7915.267-6.3741.3100.2752.4171,124
169.7835.841-7.0851.4530.2762.698510
177.4144.155-5.0961.3780.3132.318233
189.6135.672-6.9761.3330.2742.637102

全部 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 框架的长期对话者。最终文本由作者独立完成。

写作声明: 本论文所有文本, 无论中英文, 均由作者秦汉独立完成。论文代表作者唯一的责任。