Flow-Efficiency Drift Scheduler
Implementation & benchmark of arXiv:2609.03193 — Generative Nested Sampling of Atomistic Thermodynamic Landscapes
Source paper: Generative Nested Sampling of Atomistic Thermodynamic Landscapes arXiv:2609.03193 ⓘ · analyzed Sep 4, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Turn constrained-flow generation efficiency into an online diagnostic and controller for neural sampling. When the target ensemble changes faster than the flow can track or becomes internally complex, automatically shorten the training window, increase flow updates, or fall back to local MCMC instead of silently accepting biased or highly correlated samples.
Formulas
Mathematical statement
Let $q_i(x)=q_\theta(x\mid\lambda_i)$ be the current proposal and $\pi_i$ the constrained target. The normalized importance-weight effective sample size is $\eta_i=\mathrm{ESS}_i/M$, where $M$ is the number of proposals. A target-drift proxy is the symmetric KL divergence between consecutive flow distributions, estimated from samples as $D_i^{\mathrm{sym}}=[D_{\mathrm{KL}}(q_i\|q_{i-1})+D_{\mathrm{KL}}(q_{i-1}\|q_i)]/2$. Low $\eta_i$ with low drift indicates poor internal mode coverage, while high drift indicates stale training data. The controller adapts FIFO window length $K_i$, gradient-update count $R_i$, and optional local-MCMC correction using thresholds $\eta_{\min}$ and $D_{\max}$.
Implementation notes
1. Integration point: place this controller around the conditional-flow sampler used for diffusion sampling, neural energy models, Bayesian neural networks, or world-model trajectory distributions. Store each proposal's flow log density, prior log density, constraint result, and score. 2. Pseudocode: after every nested level, compute normalized ESS from importance weights; evaluate both $q_i(x)$ and $q_{i-1}(x)$ on a small proposal buffer to estimate symmetric KL drift; if $D_i^{\mathrm{sym}}>D_{\max}$, discard old examples, halve the FIFO window $K$, and perform additional flow updates; if $\eta_i<\eta_{\min}$ while drift is small, increase $R_i$ and apply 5 to 20 local HMC or Langevin correction steps; if efficiency remains high for several levels, enlarge $K$ and reduce updates. Flag a level when ESS is below 1% of proposals. 3. The paper supplies the mechanism that efficiency varies non-monotonically along annealing and becomes poor in internally complex liquid-like ensembles. Thresholds, window limits, and fallback lengths must be estimated in a pilot run. 4. First experiment: use a 16-dimensional mixture of eight Gaussians with controllable separation and threshold drift. Compare fixed-window flow sampling with the adaptive controller. Prediction: under rapidly changing targets, adaptation keeps $\eta_i\ge0.1$ by shortening $K$, whereas fixed-window training falls below $0.02$. In a dense multimodal interval, efficiency should show a reproducible local minimum; the controller should increase compute there while preserving evidence accuracy. Verify that log-evidence error remains stable when the fixed-window baseline suffers ESS collapse.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Not reproduced ✗
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 reproducible 16-D eight-mode mixture experiment with ESS, symmetric-KL drift estimation, adaptive FIFO window/update scheduling, local correction, and importance-evidence diagnostics. The math checks passed: ESS stayed in range and was invariant to log-weight shifts; symmetric KL was positive and approximately symmetric within Monte Carlo noise. Adaptation improved mean ESS from 0.1445 to 0.1791, reduced ESS<0.02 levels from 21 to 15, and reduced mean log-evidence error from 5.42 to 0.98, but it did not maintain the promised eta>=0.1 floor and used 148 versus 28 flow updates plus 120 correction steps, so this is a promising signal rather than a demonstrated efficient win.
- Agent confidence
- 7/10
- Baseline
- mean eta=0.1445; min eta=0.000880; 21/28 levels eta<0.1; 21/28 eta<0.02; mean absolute log-evidence error=5.4232; total R=28
- Idea
- mean eta=0.1791; min eta=0.000878; 18/28 levels eta<0.1; 15/28 eta<0.02; mean absolute log-evidence error=0.9774; total R=148; local correction steps=120
Limitations: This is a toy mixture-center simulator, not a learned neural flow or true HMC/Langevin implementation. The local correction uses access to the current target center as an oracle-like proxy, thresholds were hand-selected, only one fixed seed/run was used, and no wall-clock or FLOP-normalized comparison was measured.
How to run: python3 scheduler_experiment.py
Stage 2 — Benchmark vs. tuned baseline (arXiv:2609.03193)
Significantly WORSE than the tuned baseline (+133.0%, p=0.0081, wins 0 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.0034
- Idea mean
- 0.008
- Effect (Δ)
- 0.0046 (+133.0%; negative = idea better)
- Wins
- 0 / 8 paired seeds
- p-value
- 0.0081 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- n/a
- Mechanism
- Not confirmed ✗
- Practical effect
- Hurts
- Baseline tuning
- swept over 3 configs
Tested only the built-in dynamics track with 400 training and 160 test examples, rnn_small, three shared lr settings, and 8 paired seeds. The benchmark diagnostic is an observable supervised residual/prediction proxy rather than true flow density ESS or local MCMC; no vision, sequence, or larger-scale neural sampling experiments were run.
How to run: python3 bench_flow_scheduler.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Stage 2 — Benchmark latest Failed ✗
Agent confidence: 10/10
Implemented the Flow-Efficiency Drift Scheduler on the structurally matched dynamics track with identical rnn_small systems and an equal 40-step budget. ESS math sanity checks passed, but the adaptive FIFO/replay controller was significantly worse than the tuned baseline: paired delta +0.00455143 MSE with permutation p=0.0081. The behavior-derived signature did not confirm the predicted low-ESS regime because mean ESS was 0.7777 and no levels crossed the low-ESS threshold.
- Baseline
- Best lr=0.006, epochs=10; 8-seed mean test MSE=0.0034221157, std=0.0013572790
- Comparison
- delta_mean=+0.0045514294 (idea-baseline; lower is better), idea_wins=0/8, p_value=0.0081, verdict=idea worse (significant)
- Idea
- Best lr=0.006, epochs=10; 8-seed mean test MSE=0.0079735451, std=0.0027093770; total_steps=40
- Math check
- ESS in range=true; ESS shift-invariance absolute error=1.11e-16
- Mechanism signature
- mean_eta=0.7777271, min_eta=0.6528541, mean_drift=0.3253299, low_eta_levels=0, confirmed=false
Benchmark result
Significantly WORSE than the tuned baseline (+133.0%, p=0.0081, wins 0 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.0034
- Idea mean
- 0.008
- Effect (Δ)
- 0.0046 (+133.0%; negative = idea better)
- Wins
- 0 / 8 paired seeds
- p-value
- 0.0081 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- n/a
- Mechanism
- Not confirmed ✗
- Practical effect
- Hurts
- Baseline tuning
- swept over 3 configs
Tested only the built-in dynamics track with 400 training and 160 test examples, rnn_small, three shared lr settings, and 8 paired seeds. The benchmark diagnostic is an observable supervised residual/prediction proxy rather than true flow density ESS or local MCMC; no vision, sequence, or larger-scale neural sampling experiments were run.
How to run: python3 bench_flow_scheduler.py
Artifacts
- bench_flow_scheduler.py 5.4 KB View Raw
- bench_report.json 2.6 KB View Raw
- report.md 1.6 KB View
- report_bench_2026-09-04T154929.md 4.3 KB View
- results.json 12.4 KB View Raw
- scheduler_experiment.py 5.2 KB View Raw