Experiments / V2.41
V2.41
Deep Numerical Tests COMPLETE

V2.41 - Integrated Convergence Pipeline — Report

V2.41: Integrated Convergence Pipeline — Report

Status: PARTIAL (2/5 checks PASS, 8/8 tests pass)

Objective

Wire V2.19’s factor-of-4 corrected QFI and V2.20’s discrete entanglement entropy into V2.17’s 8-step pipeline, then run a convergence study from N=500 to N=2000 with multiple seeds to establish whether Gamma* converges toward 1.0 in the large-N limit.

Why This Matters

V2.17 established the 8-step pipeline (C -> SJ -> capacity -> T -> S -> Clausius -> Einstein -> metric), but used an uncorrected QFI (missing factor of 4, identified in V2.19) and fell back to the analytic entropy formula S = piT/6 when discrete entropy failed (identified in V2.20). V2.41 integrates both corrections and runs the first multi-seed, multi-N convergence study with all corrections applied.

If Gamma* converges to 1.0 as N -> infinity, it means the channel capacity of the SJ vacuum reproduces the Unruh temperature exactly, and the entire pipeline from quantum information to Einstein’s equations becomes self-consistent in the continuum limit.

Corrections Applied

StepV2.17 (old)V2.41 (new)
3. QFI capacityF = max(omega^2 * |F(omega)|)F = max(4 * omega^2 * |F(omega)|) (V2.19 Braunstein-Caves)
5. EntropyAnalytic fallback S = pi*T/6Discrete symplectic eigenvalues only (V2.20)
2. SJ vacuumDense eigh onlyOptional sparse eigsh for N > 500

Results

Phase 1: Integration Validation (N=200, dense SJ)

MetricValue
Pipeline completeTrue
Gamma* (from run_integrated_pipeline)-0.27
Gamma* (from V2.19 compute_gamma_star)Computed separately
Clausius residual2.51 (251%)
All 8 steps succeedTrue
QFI factor-of-4 appliedTrue
Discrete entropy (no analytic)True

Note: The negative Gamma* from the direct pipeline is an artifact of the slope-law fit at small N=200. The V2.19 compute_gamma_star method (used in the convergence study) gives positive, meaningful values.

Phase 2: Sparse SJ Validation (N=500)

MethodGamma*Time
Dense (reference)-0.1060.3s
Sparse (k=50)-0.14112.7s
  • Relative difference: 32.8%
  • Speedup: 0.02x (sparse is 40x slower)
  • [FAIL] Agreement < 15%: 32.8%

The sparse eigsh with k=50 modes does not agree well with the full dense eigendecomposition at N=500. This is because k=50 captures only the top 50 modes of iDelta, discarding information that matters for the capacity profile. The speedup is also negative because scipy.sparse.linalg.eigsh has significant overhead that only pays off at N >> 2000.

Phase 3: Convergence Study (N=500, 1000, 2000; 3 seeds each) — SUPERSEDED

Note: Phase 3 used the broken sparse SJ (k=50 eigsh for all N > 100). Results below are included for historical reference but are unreliable. Phase 5 replaces this with correct dense SJ.

NGamma* (median)Gamma* (std)c/3 (median)Clausius residualSuccess
5001.6600.6923.291992%3/3
10000.4370.1521.751127%3/3
20001.0760.5030.99304%3/3

Phase 4: Convergence Trend — SUPERSEDED

Power-law fit from Phase 3 (broken sparse): gamma* = 1 + 10554 * N^(-1.558). See Phase 5 for correct results.

Phase 5: Large-N Dense Convergence (N=1000, 2000, 5000; 3 seeds each)

Bug fix applied: use_sparse = N > 100 changed to use_sparse = False. All runs now use exact dense np.linalg.eigh with the memory-efficient sj_wightman_efficient(C) construction.

NGamma* (median)Gamma* (std)c/3 (median)Clausius residualTime/run
10000.5360.0648.41305%0.7s
20000.9500.91719.02756%4.0s
50001.1760.10854.58917%52.9s

Gamma trend:* 0.536 → 0.950 → 1.176 (approaching 1.0 from below, slight overshoot at N=5000).

Phase 5: Convergence Trend

The trend fitter reports “already near target” at N=5000: Gamma* = 1.176.

Unlike the superseded Phase 3 fit, the Phase 5 data shows Gamma* approaching 1.0 from below and slightly overshooting, not monotonically decreasing from above. This is consistent with a non-monotone approach to the continuum limit.

Check Summary

CheckStatus
[PASS] Gamma* near 1.0 at N=50001.176 (within 18%)
[FAIL] c/3 converging to 1/3Growing: 8.4 → 19.0 → 54.5
[PASS] Convergence trendNear target at N=5000
[FAIL] Clausius residual < 20%8917% at N=5000 (diverging)
[FAIL] Sparse/dense agreementN/A (sparse disabled)

Key Findings

  1. Gamma approaches 1.0.* With exact dense SJ: 0.536 → 0.950 → 1.176. At N=5000, Gamma* = 1.176 is within 18% of the target. The non-monotone approach (below at N=1000, above at N=5000) is consistent with oscillatory convergence.

  2. The sparse SJ bug corrupted all previous convergence data. The Phase 3 results (Clausius 1992% → 304%) were measured with broken rank-50 eigsh, which had 32.8% disagreement with dense. The “improving” Clausius trend was an artifact of error cancellation in the sparse approximation.

  3. Clausius residual DIVERGES with correct SJ. With exact dense SJ: 1305% → 2756% → 8917%. This is the opposite of what was expected. The discrete entropy S grows much faster than the continuum S = piT/6 as N increases.

  4. c/3 diverges (UV catastrophe). The central charge proxy grows as c/3 ~ N (8.4 → 19 → 54.5), indicating the discrete entanglement entropy has an unregularized UV divergence. In the continuum CFT, S ~ (c/3) ln(L/epsilon); on the causal set, the UV cutoff epsilon ~ 1/sqrt(rho) shrinks as N grows, and the discrete symplectic eigenvalue method doesn’t properly regulate this divergence.

  5. Low seed variance at N=1000 and N=5000. At N=1000: std = 0.064 (12%). At N=5000: std = 0.108 (9%). The outlier at N=2000 (seed 179, Gamma* = 2.86) drives the high variance there.

  6. Dense SJ scales well. N=5000 completes in 53s with the memory-efficient construction. N=10000 is feasible (~150s, 4 GB peak).

Non-Circularity Audit

StepMetric-Free?Input Source
1. Causal sprinklingConformal onlyCoordinates (flat 1+1D)
2. SJ vacuumYesCausal matrix C
3. QFI capacity (V2.19)YesWightman W
4. TemperatureYesCapacity profile
5. Entropy (V2.20)YesWightman W, Pauli-Jordan Delta
6. ClausiusYesT, S from above
7. EinsteinYesClausius result
8. Metric recoveryYesEinstein result

Score: 7/8 metric-free (Step 1 uses conformal structure for sprinkling).

Limitations

  1. Discrete entropy has UV divergence. The symplectic eigenvalue method gives S ~ N, not S ~ ln(N). This is the dominant failure mode — Clausius, c/3, and all entropy-derived quantities diverge.
  2. Only 3 seeds per N. N=2000 shows a large outlier (Gamma* = 2.86), suggesting 5-10 seeds are needed.
  3. Non-monotone Gamma convergence.* Gamma* approaches 1.0 from below (N=1000: 0.54) then overshoots (N=5000: 1.18). This could be oscillatory convergence or systematic bias.
  4. Gamma from direct pipeline is negative.* The slope-law alpha from extract_temperature can be negative at small N, giving meaningless Gamma*. The V2.19 compute_gamma_star method is more robust.

Path Forward

  1. Fix entropy UV divergence (CRITICAL). The discrete symplectic entropy needs proper UV regularization. Options: (a) subtract the vacuum contribution (mutual information instead of entropy), (b) use the SJ-specific regularization from Sorkin & Yazdi (2018), (c) compute entropy density dS/dA directly instead of total S.
  2. Push to N=8000-10000 for tighter Gamma* statistics (dense SJ feasible at ~4 GB).
  3. Increase to 5-10 seeds to reduce variance, especially at N=2000.
  4. Separate Gamma from Clausius.* The capacity → temperature → Gamma* pipeline works independently of entropy. Focus the convergence claim on Gamma* → 1.0, which doesn’t depend on the broken entropy.
  5. Feed corrected Gamma into V2.42* (forward Jacobson) with eta derived from regularized entropy.

Test Coverage

8/8 tests pass: QFI correction validation, entropy discreteness check, end-to-end pipeline (3 sub-tests), Gamma* computability, convergence study smoke test, non-circularity audit.