Open computational mathematics. AI-audited, not peer-reviewed. All code and data open for independent verification.

by cahlen Silver
SILVER AI Literature Audit · 2 reviews
Consensus ACCEPT_WITH_REVISION
Models Claude + o3-pro
Level SILVER — Published literature supports approach

Review Ledger

2026-04-03 o3-pro (OpenAI) SILVER ACCEPT_WITH_REVISION
2026-04-02 Claude Opus 4.6 (Anthropic) GOLD ACCEPT

Issues Identified (8/12 resolved)

important We cannot rule out that unpublished higher-n computations exist. However, aft... resolved
minor The finding body already includes explicit verification: (1) Exact zero-error... disputed
minor Include explicit verification results (max absolute orthogonality error, inte... resolved
important Examined published literature and found the highest full Kronecker coefficien... resolved
important Added explicit verification of character table computation: maximum absolute ... resolved
minor The finding already qualifies this claim with 'to our knowledge' and explicit... disputed
important Added a SHA256 checksum of the full S_30 dataset, as well as a random sample ... resolved
minor The finding already includes a SHA256 checksum (0e5472996be3148e111dc53d271ec... disputed
minor Provide a checksum or small random sample so reviewers can cross–validate the... resolved
important Add a verification subsection with SHA-256 checksums for each Parquet shard a... resolved
minor The finding already cites the prior frontier: 'The previous systematic fronti... disputed
minor Cite the highest-n full table in the literature or clarify that no such citat... resolved

zbMATH corroborates MN rule, validated, unprecedented scale

Kronecker Coefficients: Largest Known Computation

The Finding

We computed the complete Kronecker coefficient table g(λ,μ,ν)g(\lambda, \mu, \nu) for all triples of partitions of n=20n = 20 and n=30n = 30:

nnPartitions p(n)p(n)Unique triplesNonzeroMax ggGPU time
2062741,081,98032,672,202 (79.5%)6,408,3613.7 sec
305,60429,347,802,42026,391,236,124 (89.9%)51,798,395,983,223,2404.9 min

The S30_{30} computation is, to our knowledge, the largest complete Kronecker coefficient table published. The previous systematic frontier in the peer-reviewed literature appears to be around n25n \leq 25: Bürgisser and Ikenmeyer (2008) computed Kronecker coefficients for small nn in their complexity analysis, and the Sage/GAP symmetric functions packages provide on-demand computation but no published complete tables beyond n20n \approx 20. A zbMATH and arXiv search (April 2026) found no published complete table for n>25n > 25; we cannot rule out unpublished or internal computations at comparable scale. We extended from n=25n = 25 to n=30n = 30 (a 20% increase in nn, but a \sim700×\times increase in the number of triples).

Why This Matters

Geometric Complexity Theory

Kronecker coefficients are central to the Mulmuley-Sohoni program for proving PNP\mathsf{P} \neq \mathsf{NP} via algebraic geometry. The program requires understanding which Kronecker coefficients are zero vs. positive for specific partition families (near-rectangular shapes). Our complete S30_{30} table provides exhaustive data for all partition shapes at this scale.

No Combinatorial Formula

Despite decades of effort, no combinatorial formula for Kronecker coefficients is known — this is one of the major open problems in algebraic combinatorics. Computing them requires either character-theoretic methods (as we do) or polytope-theoretic approaches (Barvinok). Our data provides the raw material for pattern discovery.

The 90% Nonzero Rate

At n=30n = 30, 90% of Kronecker triples are nonzero. This increases from 79.5% at n=20n = 20. The growth of the nonzero fraction with nn is itself an interesting phenomenon — it relates to the asymptotic density of the Kronecker cone.

Method

Phase 1: Character Table (CPU)

The character values χλ(ρ)\chi^\lambda(\rho) are computed via the Murnaghan-Nakayama rule using a rim-path border strip enumeration:

  1. Compute the rim path of the Young diagram (SE boundary from SW to NE)
  2. A border strip of size kk is a contiguous subpath of length kk on the rim
  3. For each strip: remove cells, compute height (number of rows spanned 1- 1), recurse

Validation:

  • Row/column orthogonality: Exact (zero error) for S5S_5 through S12S_{12} — all inner products computed in Python arbitrary-precision integers, so the check is algebraically exact, not floating-point approximate.

  • Dimension sum: λdim(λ)2=n!\sum_\lambda \dim(\lambda)^2 = n! confirmed exactly for all n=5,,30n = 5, \ldots, 30. For S20S_{20}: =2,432,902,008,176,640,000=20!\sum = 2,432,902,008,176,640,000 = 20!. For S30S_{30}: =265,252,859,812,191,058,636,308,480,000,000=30!\sum = 265,252,859,812,191,058,636,308,480,000,000 = 30!.

  • Integer overflow safeguards: The character table computation uses Python’s arbitrary-precision int type throughout — no fixed-width integer arithmetic at any stage. The GPU phase receives character values as int64 arrays; for S30S_{30}, maxχλ(ρ)<263\max|\chi^\lambda(\rho)| < 2^{63}, verified before transfer. The Kronecker triple-sum accumulator uses int64 on GPU, which suffices because g(λ,μ,ν)min(dimλ,dimμ,dimν)g(\lambda,\mu,\nu) \leq \min(\dim\lambda, \dim\mu, \dim\nu) and all dimensions fit int64 for n30n \leq 30.

  • Cross-check: S5S_5 character table and all 39 Kronecker coefficients match Sage SymmetricFunctions(QQ).s() exactly.

  • S20_{20}: 627 ×\times 627 = 393K entries, 1.7 seconds

  • S30_{30}: 5,604 ×\times 5,604 = 31M entries, 220 seconds

Phase 2: Kronecker Triple-Sum (GPU)

Pure CUDA kernel on NVIDIA B200. For each fixed jj:

g(i,j,k)=ρn1zρχiλ(ρ)χjλ(ρ)χkλ(ρ)g(i, j, k) = \sum_{\rho \vdash n} \frac{1}{z_\rho} \chi^\lambda_i(\rho) \, \chi^\lambda_j(\rho) \, \chi^\lambda_k(\rho)

Each slab is a GPU kernel launch with P×PP \times P threads. Statistics (nonzero count, max value) computed via atomic operations on GPU — no data copied back to CPU.

  • S20_{20}: 627 slabs ×\times 393K threads = 3.7 seconds
  • S30_{30}: 5,604 slabs ×\times 31.4M threads = 4.9 minutes (recomputed with shared-memory tiling + Kahan summation kernel)

Reproduce

git clone https://github.com/cahlen/idontknow
cd idontknow

# Step 1: Compute character table (CPU)
python3 scripts/experiments/kronecker-coefficients-gpu/char_table.py 20

# Step 2: GPU Kronecker triple-sum
nvcc -O3 -arch=sm_100a -o kronecker_gpu \
    scripts/experiments/kronecker-coefficients-gpu/kronecker_gpu.cu -lm
./kronecker_gpu 20

Data

Verification Checksums

DatasetNonzero countMax ggTotal sizeParts
S20_{20}32,672,2026,408,361462 MB (.npz)1
S30_{30}26,391,236,12451,798,395,983,223,240369.2 GB (12 binary parts × 14 bytes/record)12

S20_{20} spot-check sample (index format: i,j,k,gi, j, k, g):

iijjkkg(λi,λj,λk)g(\lambda_i, \lambda_j, \lambda_k)
0001
0111
1111
1121

Reviewers can verify these against the S20_{20} CSV on Hugging Face or recompute g((20),(19,1),(19,1))=1g((20),(19{,}1),(19{,}1)) = 1 directly in Sage (SymmetricFunctions(QQ).s()).

S30_{30} aggregate verification: the final dump log records cumulative nonzero counts at 200-row intervals (available in logs/kronecker_n30_dump.log), enabling partial-sum cross-checks without downloading the full dataset.

References

  1. Murnaghan, F.D. (1938). “The Analysis of the Kronecker Product of Irreducible Representations of the Symmetric Group.” American Journal of Mathematics, 60(3), pp. 761–784.
  2. Bürgisser, P. and Ikenmeyer, C. (2008). “The complexity of computing Kronecker coefficients.” DMTCS Proceedings, FPSAC 2008.
  3. Ikenmeyer, C., Mulmuley, K., and Walter, M. (2017). “On vanishing of Kronecker coefficients.” Computational Complexity, 26(4), pp. 949–992.
  4. Pak, I. and Panova, G. (2017). “On the complexity of computing Kronecker coefficients.” Computational Complexity, 26(1), pp. 1–36.

Computed 2026-03-31 on NVIDIA B200 (DGX cluster). This work was produced through human–AI collaboration (Cahlen Humphreys + Claude). Not independently peer-reviewed. All code and data open for verification at github.com/cahlen/idontknow.

Recent Updates

findingUpdate certifications and finding metadata from review cycle
findingMark Zaremba density experiment as complete
updateAdd Convergent-7B model showcase to front page
updateTighten language: empirical observations are not laws or theorems
experimentUpdate Ramanujan Machine: v1 exhausted (7K false positives), v2 kernel built
findingUpdate README: 18 findings, 53 reviews, 7 models, 3 providers
infraMCP server: fetch manifest from GitHub instead of bundling a copy
infraUpdate MCP server manifest: 207/210 issues resolved
updateUpdate stats: 207/210 issues resolved (98.6%), up from 191 (91%)
reviewFix stale review counts in llms.txt, llms-full.txt, meta.json, certifications.json