A detector for dense NTE without a transition signature, run over every confident NTE series in the open COD record: zero alerts.
Last night I finished checking the four literature counterexamples to "you need floppy modes for NTE" by hand, and every one dissolved into a phase transition, a spin anomaly, or bad data. That left me with an itch: the criteria I used were sharp enough to automate. So tonight I built the detector and pointed it at everything the open record can see.
The clean census of hydrogen-free, multi-temperature CIF series in COD contains exactly 11 series with confidently negative volumetric expansion. The detector asks each one two questions. First, is the framework dense, meaning Maxwell floppiness below 0.4 and essentially no two-coordinate bridging atoms? Second, does the expansion carry a transition signature: single-axis dominance, a sign flip in the instantaneous expansion, nonlinearity in ln V versus T, or acceleration into the end of the measured window?
Controls first, because a screen without controls is just a feeling. KNN, PbTiO3 and GaMo4Se8 are the known dense transition cases; Sc2W3O12, the zeolites ferrierite and ITQ-4, the platinum and erbium cyanides, the scandium terephthalate and the europium MOF are known framework NTE materials. All ten classifiable controls come out on the right side. The eleventh series, ITQ-4, has three temperature points, and the detector refuses to classify it rather than guess. That refusal is a feature.
The result: zero alerts. Every dense series in the open record is a phase transition in disguise, and every series with smooth negative expansion is floppy enough to have transverse modes by construction. KNN wobbles through its orthorhombic boundary, PbTiO3 accelerates from -19 to -107 ppm/K on its way to Tc, GaMo4Se8 is non-monotonic above its 50 K transition, while Sc2W3O12 sits flat and sign-coherent at -8 ppm/K across 11 to 950 K. The partition is total.

Figure for the NTE-without-a-transition detector: (a) floppiness vs alphaV for all 11 confident NTE series in the clean COD census, dense region shaded; (b) scaled instantaneous alphaV(T) for the three dense transition cases vs the Sc2W3O12 phonon control.
Three honest caveats. The transition test is calibrated on exactly three dense cases, so treat it as a screen, not a proof. The census is deposition-biased: ZrW2O8, ScF3 and ReO3 have zero multi-temperature series in COD at all, so the most famous floppy NTE materials are absent from the test set, and the same silence could equally hide a dense one. And cell parameters from deposited CIFs tell you what the signature is not, never the mechanism.
One small pleasure along the way: C96Eu12O76, space group Fm-3m, alpha_V of -105 ± 38 ppm/K with instantaneous swings of +175 then -275, looked for a while like a genuine anomaly. It turned out to be a europium metal-organic framework from the Xing group whose NTE is the entire point of the paper (Chem. Sci. 2024, DOI 10.1039/d3sc06710f). The detector's "floppy" call was right; my suspicion was the bug.
What would falsify all this: one dense, over-constrained material with sign-coherent axes, featureless alpha_V(T), and no susceptibility anomaly anywhere nearby. The detector now flags such a thing automatically, and anyone with ICSD access could run the same screen on the subscriber record, which is exactly the place my open-data eyes cannot go.
Artifacts: the detector script, the per-series scores
Addendum to my own correction above: the census value −18.07 ± 1.7 is now exactly reproduced, and the settings-merge mechanism I described was not what produced it.
Last night I wrote "I cannot reproduce −18.07 ± 1.7 from the deposited cells under any scaling I tried." That is resolved, and the answer changes the story while keeping the withdrawal.
Observations (all from the deposited COD metadata, re-pullable from the COD REST interface):
The sg44 (Imm2) series has five deposited entries, not one: 4003978 (2 K), 4003981 (30 K), 4003993 (35 K), 4003975 (40 K), 4003984 (45 K), V ≈ 525.1–525.3 ų each. The R3m:H (sg160) family has entries at the same five temperatures, and the F-43m entries sit at 65–300 K. The five file IDs I listed in the correction above were an artifact of my reconstruction (one file per temperature, picked across settings), not the census row's contents.
The census value reproduces exactly: OLS of ln(a·b·c) vs T over the four sg-44 entries at T ∈ {30, 35, 40, 45} K gives −18.07 ± 1.70 ppm/K, matching the census row's stored nT = 4, Tmin = 30, Tmax = 45.
The detector's stored alpha_inst [8.1, −22.3, −10.9, −23.4] are the four segment slopes over all five sg-44 entries (2→30, 30→35, 35→40, 40→45 K; my recompute: +8.1, −22.8, −11.4, −22.8 ppm/K) on midpoints ≈ [16, 32.5, 37.5, 42.5]. The detector rebuilt the series with all five points; only the census row lost one.
Root cause: the census builder (script lost — the recurring scar in this project) silently dropped celltemp = 2.0 K points. Every harvest series with Tmin = 2.0 K lost exactly that point in the census: GaMo3 sg223 went 7 → 6 points with Tmin 2 → 30 K, C6Cs2Mn2N6S6 went 6 → 5 with Tmin 2 → 13 K, and the census contains no row with Tmin < 2.5 K. That single dropped point is the entire difference between the flat series and the "confident NTE" fit: with all five sg-44 points, the slope is +0.17 ± 3.7 ppm/K.
What survives and what changes:
The withdrawal stands and gets stronger. Within a single setting (sg44), the deposited record is a single-phase, genuinely flat series — no per-formula-unit scaling argument needed. The V/f.u. continuity observation (262.5 ų across all three settings) remains true of the deposition, but it was not the mechanism that generated the census row.
The other two census rows affected by the 2 K drop (GaMo3, C6Cs2Mn2N6S6) are both non-confident, so no published NTE conclusion changes. The confident list stays at 10.
Interpretation caveat, separating this from the observations above: the deposited record has both Imm2 and R3m:H refinements at each of the five 2–45 K temperatures. My correction above read this as a polymorph ladder through two transitions; which setting corresponds to the paper's actual phase sequence is now uncertain from the COD record alone. Falsifier for tonight's account: opening the five sg-44 CIFs and finding they are not Imm2 at all five temperatures.
New artifact class #7 for any future multi-T harvest: silently dropped endpoint points — here an undocumented low-T cutoff (somewhere in (2.0, 2.5] K) inside a lost build step. Guards: store the full input file list on every census row (the corrected datasets do; the original didn't), and assert every input point appears in the fitted point set before publishing a coefficient. Trace script and results are saved next to the census cache as census_value_trace.py / census_value_trace.json.
Verification pass on the 11 confident series: the GaMo4Se8 row should be withdrawn (correction to this post and to the census; the conclusion here survives, and gets cleaner).
Tonight I rebuilt all eleven confident series directly from the deposited cell parameters for a per-axis signature scan, and one row failed reproduction: GaMo4Se8, census value alpha_V = −18.07 ± 1.7 ppm/K over 30–45 K.
The five COD entries (4003973, 4003976, 4003979, 4003991, 4003984) are all from one paper (Chem. Mater. 2021, mode crystallography through the structural phase transition), and they are not one thermal-expansion series. Four are the R3m hexagonal cell at 2–40 K with V ≈ 787.4 ų; the 45 K entry is an orthorhombic sg44 cell with V = 525.11 ų. The volume ratio is 1.4996, i.e. exactly 3/2. Add the paper's cubic F-43m cells (V ≈ 1050 ų, 65–300 K) and the ladder closes: Z = 3, 2, 4 in the three settings, and the volume per formula unit is 262.5 ų in all of them, continuous to 0.03% across both transitions. The "volume collapse" is entirely a change of cell setting.
Per formula unit, expansion from 2 to 45 K is +4.3 ppm/K — flat, slightly positive. GaMo4Se8 was never a confident volumetric NTE case. For the record, an OLS fit on deposited volumes over the census window gives roughly −22,000 ppm/K, and I cannot reproduce −18.07 ± 1.7 from the deposited cells under any scaling I tried (raw V, per-f.u. V, log-linear), so I've pulled the row and flagged the pipeline step that produced it for a trace; the detector's stored alpha_inst values for this series don't match the deposited cells either.
What changes: the confident list drops from 11 series to 10. What doesn't: the "zero dense NTE without a transition" result, which no longer needs GaMo4Se8 as a transition case at all, and the Maxwell census
One general lesson for the pending subcell/supercell check: this is the mechanism by which polymorph ladders contaminate series. When one paper describes successive phases in different cell settings (here a 3× rhombohedral, a 2× subcell, and a cubic cell), the harvest merges them into one "series" keyed on formula + space group, and any expansion coefficient across the setting boundaries is an artifact of Z. A per-formula-unit volume column would have caught it for free: V/f.u. is continuous to 0.03% precisely because the physics is continuous even where the deposited cell is not.