Null Controls Are Necessary but Not Sufficient: Evaluating Acceptance Guards for Scientific Discovery and Interpretability
Organizations: Project Aether
Abstract
Automated discovery systems propose scientific explanations and test which to accept. Does passing a test establish the intended claim? We develop four executable operations for changing physical descriptions: identifying redundancies (L1), coarse-graining dynamics (L2), transferring structure (L3), and inferring missing quantities from conservation imbalances (L4). We evaluate these operations and three interpretability settings through null acceptance (NAR), positive acceptance (PAR), and consequence agreement (CAR). A reduction mapping every state to zero predicts its own future perfectly. It passes closure on 96% of 129 public chaotic systems despite retaining no dynamical variation. An independent observable-preservation test rejects it throughout. On controlled benchmarks, the combined criteria accept 24/24 genuine reductions and reject 72/72 nulls. The other operations distinguish symmetries from redundancies, validate transfers, and reject spurious conservation imbalances. In language and ESM-2 protein models, linear CKA remains unchanged when a proposed correspondence is randomly rotated, whereas localized interventions distinguish genuine mappings from tested nulls. In sparse autoencoders, enrichment rejects all 16 tested nulls, yet independent ablation supports 45% of evaluable enriched features, compared with 0.5% of matched controls. These results reveal two acceptance gaps: descriptive tests can accept scientifically empty candidates, and accepted candidates can lack independently verified consequences. Reliable automated discovery requires evaluating both what a test accepts and the scientific consequences it establishes. Code and data: https://real-project-aether.github.io/project-aether/.
Figures & tables
| Candidate revision | Descriptive ambiguity | Independent consequence | |
|---|---|---|---|
| L1 | Remove unobservable directions | Global symmetry resembles redundancy | Locality separates 1 redundancy from 4 conserved symmetries |
| L2 | Reduce dynamical variables | Empty states satisfy closure | Preserve observables; 0/72 nulls and 24/24 true subspaces accepted |
| L3 | Transfer structure across systems | 2/6 broken controls match structurally | Target dynamics validate 6/6 transfers, reject 6/6 controls |
| L4 | Infer an unobserved quantity | Calibration offsets appear significant | Noise-scale test rejects offsets; recovers vs. |
| Arm | Paired | Selected CAR | Control CAR | CAR |
|---|---|---|---|---|
| Pythia-160m L4 | 93 | 0.24 | 0.00 | |
| Pythia-160m L8 | 111 | 0.54 | 0.00 | |
| Pythia-410m L16 | 102 | 0.58 | 0.00 | |
| Pythia-410m L8 | 99 | 0.42 | 0.02 | |
| pooled | 405 | 0.45 (CI – ) | 0.005 |
Appendix figures & tables14 assets
Supplementary material from the paper’s appendix.
Appendix
| invariant point by point | redundancy: quotient it [L1] |
| invariant only uniformly | symmetry: keep it, read off a conserved charge [L4] |
| control | closes | pays | refused by |
|---|---|---|---|
| constant reduction (every state to zero) | 0.000 | guard 2 only | |
| random subspace | 2.167 | guard 1 | |
| the fast modes instead of the slow ones | 0.000 | guard 2 only |
| structure, in-sample | structure, held out ( ) | charge test ( ) | |
|---|---|---|---|
| real | median — 6/6 | median , worst — 6/6 pass | median , worst — 6/6 pass |
| control | median — 5/6 | median , best — 2/6 pass | median , best — 0/6 pass |
| type | preserves | destroys |
|---|---|---|
| degenerate | nothing | all structure (constant reduction, zero map) |
| randomised | dimensionality | any relation to the source |
| correspondence-breaking | marginal statistics | the pairing between source and target |
| structure-preserving | the invariant itself | only its redundancy status |
| task-confounded | a trained model’s geometry | the semantic content of the input |
| sub-threshold | an effect’s significance | its size against the noise scale |
| a reduction that keeps everything | must score exactly |
| a guard | must never rank a null candidate above the real one |
| a structural score | must be reported on data the candidate was not fitted to |
| the guard accepts… | seen in | |
|---|---|---|
| null acceptance (vacuity) | candidates carrying nothing the operation seeks | L1–L4, X1, X3 |
| unsupported acceptance | candidates an independent test does not support | X2 ( paired features) |
| episode | named | blind | bare |
|---|---|---|---|
| dark energy [ Riess and others, 1998 , Perlmutter and others, 1999 ] vs. grey dust [ Aguirre, 1999 ] | 67% | 42% | 50% |
| relativistic mass vs. Lorentz’s aether [ Lorentz, 1904 , Einstein, 1905 ] | 100% | 100% | 100% |
| quasicrystals [ Shechtman et al., 1984 ] vs. twinning [ Pauling, 1985 ] | 100% | 50% | 42% |
| two constructed controls | 83% | 25% | 38% |
| pooled, 60 trials | 87% | 48% | 53% |
| candidate | type | CKA | guard 1 | causal | AUROC | accepted | |
| X1: Pythia pairs, configurations, interventions each | |||||||
| the correspondence under test | real | 0.82 | 0.61 | 4/4 | +0.80 | – | 4/4 |
| pairings permuted | corr.-breaking | -0.04 | 0.61 | 3/4 | -0.00 | 0.995 | 0/4 |
| the same map, randomly rotated | corr.-breaking | -0.96 | 0.61 | 3/4 | -0.01 | 0.996 | 0/4 |
| covariance-matched Gaussian target | randomised | -0.04 | 0.61 | 3/4 | +0.00 | 0.995 | 0/4 |
| fitted on disjoint text | task-confounded | 0.65 | 0.61 | 3/4 | +0.01 | 0.996 | 0/4 |
| stage | 160m L4 | 160m L8 | 410m L16 | 410m L8 |
|---|---|---|---|---|
| all SAE features | 6144 | 6144 | 8192 | 8192 |
| firing rate in range | 3772 | 1248 | 2515 | 3540 |
| concept frequent enough | 3772 | 1248 | 2508 | 3540 |
| enrichment | 1263 | 605 | 1169 | 1453 |
| selected for intervention | 150 | 150 | 150 | 150 |
| arm | defined | detected | wrong sign |
|---|---|---|---|
| Pythia-160m L4 | 94 | 0.79 | 0.00 |
| Pythia-160m L8 | 110 | 1.00 | 0.00 |
| Pythia-410m L16 | 102 | 0.55 | 0.19 |
| Pythia-410m L8 | 99 | 0.07 | 0.93 |
| system | what it does | rung | why not higher |
|---|---|---|---|
| BACON [ Langley, 1981 ] | rediscovers laws from tables | L0 | searches a fixed term language |
| Eureqa [ Schmidt and Lipson, 2009 ] | symbolic regression on motion | L0 | operator set given in advance |
| AI Feynman [ Udrescu and Tegmark, 2020 ] | physics-informed symbolic regr. | L0 | recovers known forms; invents no term |
| PySR [ Cranmer, 2023 ] | evolutionary symbolic regression | L0 | fixed primitive set |
| LaSR [ Grayeli et al., 2024 ] | symbolic regr. with learned library | L0/L2 | abstractions are shortcuts, not retypings |
| DreamCoder [ Ellis et al., 2021 ] | wake–sleep library learning | L2 | names structure; no theory to reinterpret |
| candidate | accepted by the pair | ||
|---|---|---|---|
| true subspace, known by construction | 24 | 24 / 24 | recall |
| the proposer’s estimate of it | 24 | 17 / 24 | recall |
| fast, random, constant | 72 | 0 / 72 | precision |
| closure test alone | 72 | 47 / 72 |
| reduction | closes | median closure err. | max pay score | accepted by both |
|---|---|---|---|---|
| slow modes (the real candidate) | 0 / 129 | 1.810 | 0 | |
| control constant (all to zero) | 124 / 129 | 0.000 | 0 | |
| control fast modes | 0 / 129 | 1.465 | 0 | |
| control random subspace | 0 / 129 | 1.902 | 0 |
Scientific Discovery
Scientific discovery uses computational methods to find hypotheses, laws, mechanisms, and materials from data, literature, and experiments.
See all 85 papers