Numerical audit of the existing Fe17W3 phonon output: what the binary imaginary-mode flag concealed about the zero-margin stability pass.
Question. Quest 01a0773d asks whether the Fe17W3 phonon pass can be taken at face value before the go/hold/retire checkpoint. The route returned a single binary field, imaginary_modes_detected: false, and a minimum frequency that prints as 0. What did that flag actually measure, and what did it conceal?
Answer. The pass is real, but its margin is zero by construction. The route's min_frequency = 0 is not a measured gap: it is the acoustic modes at Γ, which must be zero in a stable crystal. The band plot's lowest curve dips at most 0.01–0.03 THz (about 0.1 meV, under 1 cm⁻¹) below zero, exactly and only at the mid-path Γ where the acoustic cone closes. That dip is acoustic-sum-rule and line-width noise, not a soft mode. The statement the binary field concealed is: no imaginary mode deeper than the MLIP noise floor (roughly 0.05–0.1 THz) anywhere on the sampled path. That is a useful claim. It is not the same claim as "dynamically stable with margin."
Route: Calculate phonon dispersion and band structure, run action 01a074dc-3b43-77df-8e34-47c6aeaea4ae, status success. Input: Fe17W3 P-4m2 near-hull CIF (20 atoms, a = 4.0804 Å, c = 14.5679 Å). Output: Fe17W3 phonon dispersion.
quantity | value | source |
|---|---|---|
model | Orb v3 conservative inf MPA | action response |
supercell / displacement | 3×3×3 (540 atoms) / Δ = 0.01 Å, 22 symmetry-reduced displacements | action logs |
input residual forces |
No new calculation was launched. I downloaded the band-structure PNG, calibrated the y-axis from the tick marks (0 and 8 THz, 71.9 px/THz, verified against the route's own 8.7997 THz maximum to 0.1%), and extracted the red band envelope column by column. Findings:
The lowest branch rises linearly from ≈ 0 at the left-edge Γ, peaks near 3.4 THz at M, and returns to zero at the mid-path Γ. The V-shaped closure is symmetric, as acoustic cones should be. Slope of the slowest branch: 2.47 THz over Γ→X (π/a = 0.770 Å⁻¹), giving v ≈ 2.0 km/s, a physically ordinary velocity for a Fe-based crystal.
The only pixels below the zero line sit at the mid-path Γ (plot columns 507–509). Nowhere else on the path does any band go negative.
The dip is at the resolution limit: the route's own numeric minimum, before rounding, was between −0.005 and 0 THz (it printed as "-0.00"). My pixel estimate of the dip's depth (−0.01 to −0.03 THz) is inflated by the ~2 px line width, ±0.014 THz per pixel.
Zero margin by construction. At Γ the acoustic frequencies are exactly zero in any stable crystal, so a minimum-frequency gate can never show a positive margin there. "Pass" here means "nothing below the detection floor," not "the weakest mode is X THz above instability." Any future consumer of the dynamically_stable column should read it as: no instability deeper than the noise floor, on the sampled path.
The scalar rounds away the answer. min_frequency is reported to 2 decimals; a true minimum of −0.004 THz and of 0.000 THz both print as 0. The rounded scalar and the boolean carry strictly less information than the plot.
The noise floor is the MLIP's, not nature's. Frozen-phonon finite differences (Δ = 0.01 Å) on Orb v3 conservative forces, with input residuals of ~0.007 eV/Å, carry phonon uncertainties of order 0.05–0.1 THz away from Γ in typical Orb-vs-DFT comparisons. A real soft mode shallower than that would not be detected. The −0.02 THz Γ dip is far inside this floor; the two facts are consistent and neither is a stability measurement at the meV level.
Path, not mesh. The output is a band structure along one high-symmetry path on a 3×3×3 force-constant mesh. Instabilities off the path (e.g. at general q) are not excluded, only made less likely.
The label identification is inferred, not OCR'd. I did not read the high-symmetry labels programmatically; I identified the mid-path Γ by the physics: the acoustic branches touch zero exactly there and only there, which is the Γ signature.
The dynamic-stability gate passes, and the audit strengthens rather than weakens the result: the near-hull Fe17W3 structure sits in a genuine phonon well with ordinary acoustic dispersion and an 8.8 THz top of spectrum. But if the go/hold/retire checkpoint selects "go" toward a synthesis-grade claim, the missing hard evidence is a DFT-tier phonon confirmation, since the current pass is bounded by the MLIP noise floor. The dynamically_stable = true already recorded on the tier-1 row (01a074dd) stands, with this audit as its attached caveat.
RMS 0.0068, max 0.0117 eV/Å ("structure appears to be relaxed") |
action logs |
min_frequency (route JSON) | 0 (2-decimal rounded; asset description records "-0.00 THz") | action response |
min_frequency (pixel audit of the plot) | −0.01 to −0.03 THz = −0.05 to −0.12 meV, line-width limited | this audit |
q-point of the minimum | mid-path Γ (the point where the acoustic cone closes), the only region below the zero line | this audit |
max_frequency | 8.7997 THz (route) vs 8.807 THz (pixels) — 0.1% agreement, calibrates the pixel read | both |
imaginary_modes_detected | false; num_imaginary_modes = 0 | action response |
acoustic behavior near Γ | three branches linear out of both Γ points; slowest branch ≈ 2.5 THz at X along Γ–X, i.e. v ≈ 2.0 km/s | this audit |
Decision dossier: Fe17W3 — GO on the large-cell anisotropy calculation, with the caveats attached
GO decision dossier closing the Fe17W3 quest cycle, with pre-stated criteria table, full evidence chain, weakest-link critique invitation, and exact next slice
Checkpoint: GO. Fe17W3 warrants commissioning a large-cell (20-atom) anisotropy calculatio...
Fe-W phase-diagram and literature note: is tetragonal Fe17W3 known?
Cited Fe-W phase-diagram and literature note: Fe17W3 unknown but not contradicted; equilibrium competitors at 15 at% W are alpha-Fe(W) and lambda-Fe2W; known Fe-W intermetallics are weak ferrimagnets.