Moment-Controlled Mutation
Implementation & benchmark of arXiv:2607.02171 — Theory of collective learning in populations of adaptive agents
Source paper: Theory of collective learning in populations of adaptive agents arXiv:2607.02171 ⓘ · analyzed Aug 30, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Use the paper's mean and variance dynamics to control exploration in a population of neural-network adapters. Estimate local reward curvature from the current candidates, then choose mutation strength so selection contracts diversity only when the reward landscape is locally reliable. Increase diffusion when reward noise or selection causes population collapse.
Formulas
Mathematical statement
For a normalized policy distribution \(\phi(P,t)\), define its mean \(\mu=\int P\phi(P)dP\) and variance \(\sigma^2=\int(P-\mu)^2\phi(P)dP\). Multiplying the replicator-mutation equation by \(P\) and \((P-\mu)^2\) gives \(\dot\mu=\kappa\operatorname{Cov}(P,R)\) and \(\dot\sigma^2=\kappa\operatorname{Cov}((P-\mu)^2,R)+2D\), where \(\kappa=\tilde\lambda_T\rho\), \(R(P)\) is reward, and \(D\) is diffusion. Under a locally quadratic reward \(R(P)\approx R(\mu)+R'(\mu)(P-\mu)+\tfrac12R''(\mu)(P-\mu)^2\) and a Gaussian policy population, the paper obtains \(\dot\mu\approx\kappa\sigma^2R'(\mu)\) and \(\dot\sigma^2\approx\kappa R''(\mu)\sigma^4+2D\). Negative curvature contracts variance, positive curvature expands it, and mutation supplies a guaranteed additive diversity contribution. For vector adapter parameters \(z\in\mathbb R^d\), use diagonal covariance and estimate curvature coordinatewise from weighted candidate rewards.
Implementation notes
1. Integration point: modify the mutation or noise line of a population-based LoRA or prompt-tuning optimizer. Keep N candidate vectors z_i, normalized weights w_i, EMA rewards r_i, weighted mean mu=sum_i w_i z_i, and diagonal covariance C=sum_i w_i(z_i-mu)^2. The backbone remains frozen. 2. Pseudocode: after scoring candidates, compute bar_r=sum_i w_i r_i and q_i=(z_i-mu)^2-C elementwise. Estimate local curvature per coordinate with R2_hat=sum_i w_i(r_i-bar_r)q_i/(sum_i w_i q_i^2+epsilon). Estimate reward uncertainty using repeated evaluation on a small probe batch or the weighted residual variance v_r=sum_i w_i(r_i-bar_r)^2. Set D_j=clip(D_min+gamma v_r/(1+abs(R2_hat_j)),D_min,D_max). Apply selection w_i <- softmax(log w_i+kappa(r_i-bar_r)); then mutate z_i,j <- z_i,j+sqrt(2D_j Delta_t)xi_i,j. If C_j falls below a target C_min, enforce D_j <- max(D_j,(C_min-C_j)/(2 Delta_t)), reflecting the +2D term in the variance equation. Reduce D only when R2_hat<0 and effective sample size exceeds N/2. 3. The moment equations and diffusion contribution are taken from the paper; reward, curvature, and noise estimates are empirical. Clip curvature estimates because finite populations make fourth-moment estimates noisy. 4. First experiment: train rank-2 or rank-4 LoRA adapters with N=16 on CIFAR-100 using a frozen ViT-Tiny. Compare fixed Gaussian mutation, AdamW, and moment-controlled mutation at equal candidate evaluations. Log validation accuracy, covariance trace, effective sample size, reward variance, and mutation magnitude. Success is higher final accuracy or faster accuracy growth without population collapse; a secondary success signal is stable covariance near the target rather than abrupt zero diversity or unbounded exploration.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Confirmed ✓
Stage 2 · Practical benchmark result: Significantly worse ✗
Methodology: 8 paired seeds, permutation test (p<0.05) against a learning-rate-tuned baseline. How verification works
Stage 1 — Mechanism check agent confidence 7/10
Built a NumPy MVP implementing weighted selection, coordinatewise curvature/noise-based diffusion, diversity-floor enforcement, and a fixed-diffusion control. The moment equations were numerically verified: mean and variance-selection relative errors were 0.14% and 0.16%, while the diffusion variance increment matched 2DΔt within 0.00002. Across 8 seeds, moment control improved final mean reward (-0.292 vs -0.341) and stayed closer to target covariance (absolute error 0.0205 vs 0.0448), although fixed diffusion had better mean best-candidate reward (-0.120 vs -0.152); this is a promising small toy signal, not evidence of a general optimizer win.
- Agent confidence
- 7/10
- Baseline
- Fixed D=0.018: mean final reward -0.3411±0.0811, best reward -0.1201±0.0511, covariance-target error 0.0448±0.0439, final covariance 0.1866±0.0439.
- Idea
- Moment-controlled diffusion: mean final reward -0.2918±0.0421, best reward -0.1522±0.0371, covariance-target error 0.0205±0.0228, final covariance 0.1523±0.0228, final D 0.0163±0.0049.
Limitations: This is a noisy synthetic quadratic-vector population, not CIFAR-100, LoRA, or a frozen ViT. The supplied curvature estimator was approximately factor-of-two mis-scaled for an exact Gaussian quadratic: estimate -1.136 versus true curvature -2.3. The MVP follows the supplied estimator formula rather than correcting it. No FLOP-normalized neural-network comparison, hyperparameter sweep, or larger-population test was performed.
How to run: python3 moment_controlled_mutation.py
Stage 2 — Benchmark vs. tuned baseline (arXiv:2607.02171)
Significantly WORSE than the tuned baseline (+70.8%, p=0.0081, wins 0 / 8 paired seeds).
- Benchmark
- Tabular regression (Friedman)
- Model
- mlp_tiny
- Paired seeds
- 8
- Baseline mean
- 9.2592
- Idea mean
- 15.8122
- Effect (Δ)
- 6.553 (+70.8%; negative = idea better)
- Wins
- 0 / 8 paired seeds
- p-value
- 0.0081 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±14.9%
- Mechanism
- Confirmed ✓
- Practical effect
- Hurts
- Baseline tuning
- swept over 3 configs
Only the registered tabular track was tested; vision, sequence, and dynamics were not tested. The small 400/400 split and mlp_tiny were used. Candidate-update FLOPs were not independently normalized against AdamW, and the diffusion signature is a coarse checkpoint-level measurement.
How to run: python3 stage2_bench.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Stage 2 — Benchmark latest Failed ✗
Agent confidence: 9/10
Moment-Controlled Mutation was implemented as a population training optimizer on the registered tabular/Friedman regression track with mlp_tiny. The tuned AdamW baseline had test MSE 9.2592, while the best idea setting had MSE 15.8122; the paired permutation test showed the idea was significantly worse (p=0.0081). The NN-scale diffusion signature was directionally confirmed, but this did not produce a task-metric win.
- Baseline
- AdamW: 9.2592 +/- 1.0373 MSE over 8 seeds; best lr=0.01.
- Idea
- Moment-controlled mutation: 15.8122 +/- 1.9494 MSE over 8 seeds; best lr=0.003; paired delta +6.5530 MSE; 0/8 wins; p=0.0081.
Benchmark result
Significantly WORSE than the tuned baseline (+70.8%, p=0.0081, wins 0 / 8 paired seeds).
- Benchmark
- Tabular regression (Friedman)
- Model
- mlp_tiny
- Paired seeds
- 8
- Baseline mean
- 9.2592
- Idea mean
- 15.8122
- Effect (Δ)
- 6.553 (+70.8%; negative = idea better)
- Wins
- 0 / 8 paired seeds
- p-value
- 0.0081 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±14.9%
- Mechanism
- Confirmed ✓
- Practical effect
- Hurts
- Baseline tuning
- swept over 3 configs
Only the registered tabular track was tested; vision, sequence, and dynamics were not tested. The small 400/400 split and mlp_tiny were used. Candidate-update FLOPs were not independently normalized against AdamW, and the diffusion signature is a coarse checkpoint-level measurement.
How to run: python3 stage2_bench.py
Artifacts
- bench_report.json 3.8 KB View Raw
- moment_controlled_mutation.py 4.9 KB View Raw
- report.md 1.9 KB View
- report_bench_2026-09-04T173828.md 5.0 KB View
- results/fixed_repetitions.json 0.7 KB View Raw
- results/moment_repetitions.json 0.8 KB View Raw
- results/summary.json 33.7 KB View Raw
- stage2_bench.py 6.3 KB View Raw