Pith. sign in

REVIEW 4 major objections 3 minor 9 references

This paper establishes that the faster lithium diffusion in defect-rich MoNb12O33 comes from a change in which lithium sites fill first, not from faster hopping in any single channel. In the ordered block structure, lithium first occupies s

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-03 22:35 UTC pith:FG3L4MEX

load-bearing objection Strong experimental/computational pairing that plausibly explains defect-enhanced Li diffusion in MoNb12O33, with the main caveat that the Wadsley-defect models rely on a charge-compensation assumption the paper's own EXAFS/STEM do not support for the real MNO-800 sample. the 4 major comments →

arxiv 2511.09521 v1 pith:FG3L4MEX submitted 2025-11-12 cond-mat.mtrl-sci

Role of Wadsley Defects and Cation Disorder to Enhance MoNb12O33 Diffusion

classification cond-mat.mtrl-sci
keywords Wadsley–Roth niobatesMoNb12O33lithium-ion diffusioncation disorderWadsley defectscluster expansionmachine-learned interatomic potentialanode materials
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper establishes that the faster lithium diffusion in defect-rich MoNb12O33 comes from a change in which lithium sites fill first, not from faster hopping in any single channel. In the ordered block structure, lithium first occupies slow pocket sites and only reaches fast window sites late in lithiation. In computational models containing either cation disorder or enlarged Wadsley-defect blocks, the fast window sites become occupied at much lower lithium content, matching the earlier onset of rapid diffusion seen experimentally. The result matters because it turns defect engineering into a transferable design strategy: defects that reorder site occupancy can speed up intercalation materials without changing their intrinsic per-site mobility.

Core claim

In the ordered T[3×4] block structure, lithium first fills the slow five-coordinate pocket sites and only populates the fast four-coordinate window sites (horizontal Wh and vertical Wv) late in lithiation. In models with either cation disorder (Mo displaced from tetrahedral sites into the block) or enlarged Wadsley-defect blocks (with charge balance restored by oxygen vacancies), the fast Wh and Wv sites become occupied at substantially lower lithium concentrations. Since per-site hopping barriers change only modestly, the main effect of both defect types is to reorder the sequence of site occupancy rather than to speed up any individual channel.

What carries the argument

The analysis rests on classifying lithium sites in Wadsley–Roth blocks into slow five-coordinate pocket sites and fast four-coordinate window sites (horizontal Wh and vertical Wv), then computing two things with a machine-learned interatomic potential trained to density functional theory: (1) per-site hopping barriers and diffusion coefficients, and (2) the ground-state sequence of lithium site occupancy at increasing lithiation via cluster expansion and Monte Carlo. Comparing pristine T[3×4], cation-disordered T[3×4], enlarged T[3×5], and mixed T[3×4–3×5] models isolates the two defect effects.

Load-bearing premise

The computational models of Wadsley defects keep Mo only in tetrahedral sites and restore charge balance with oxygen vacancies, but the experiments show Mo in the defect-rich sample is not uniformly tetrahedral; if the real defect-rich structure balances charge differently, the predicted early activation of fast lithium sites may not apply.

What would settle it

Measure the occupancy of fast window sites at low lithiation (e.g., x≈3–5) in MNO-800 using neutron diffraction or 7Li NMR; if fast Wh/Wv sites are not occupied early, the proposed mechanism is wrong. Alternatively, compute the defect-rich model with Mo placed in octahedral sites and no oxygen vacancies; if the fast sites no longer activate at low x, the conclusion is tied to the assumed charge-balance mode.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • MNO-800's higher capacity at 0.1C (307 vs ~293 mAh/g) and its 200 mAh/g at 10C follow from faster solid-state diffusion, which ICI overpotential analysis identifies as the dominant rate-limiting step.
  • The ~3x capacity-weighted diffusivity of MNO-800 is attributed to early activation of fast window sites; the measured D(x) curves converge at high lithiation, matching model predictions at ~75% Li.
  • Wadsley defects and cation disorder are predicted to be independently beneficial, so a material containing both should retain the advantage of each.
  • The mechanism provides a design principle: defects that lower the energy of fast window sites relative to slow pocket sites can accelerate rate capability across Wadsley–Roth anode materials.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • A testable extension is to partially prelithiate the ordered MNO-900 to x≈3–5 before high-rate cycling: if the site-occupancy sequence is the controlling factor
  • this could emulate the defect-rich early activation and recover some of MNO-800's rate performance.
  • The computational separation of defect types rests on charge balancing enlarged blocks with oxygen vacancies
  • but the paper's own EXAFS shows Mo in MNO-800 is not uniformly tetrahedral
  • if the real defect-rich structure instead places Mo in octahedral sites
  • the predicted early-activation effect would need to be re-evaluated for that charge-balance mode.
  • The same site-occupancy-reordering mechanism may operate in other intercalation hosts with multiple lithium site types of differing mobility
  • where defect engineering could be screened by computing only relative site energies rather than full diffusivities.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 3 minor

Summary. This manuscript reports a comparative study of two MoNb12O33 powders from a single sol-gel batch calcined at 800 °C and 900 °C. The defect-rich MNO-800 sample shows larger XRD microstrain, more Wadsley-type block-size variations and mixed T/E interfaces in STEM, higher and less well-defined Mo coordination in EXAFS, higher capacity and rate capability, ~3× larger capacity-weighted diffusivity from ICI, and lower volume expansion on lithiation. To explain the transport enhancement, the authors train a MACE-based MLIP on DFT and use NEB, cluster-expansion/Monte Carlo, and MD to study pristine T[3×4], cation-disordered T[3×4] variants, and enlarged-block T[3×5] and T[3×4-3×5] models with O vacancies. The computations predict that both cation disorder and Wadsley-type block enlargement cause the fast Wh/Wv window sites to become occupied at lower Li content than in the pristine model, matching the earlier onset of high diffusivity in MNO-800. The paper concludes that both defect types contribute independently to the enhanced performance.

Significance. The experimental dataset is rich and internally consistent: same-batch synthesis, multiple structural probes, operando WAXS, and a defensible ICI/SAXS diffusivity protocol. The computational work is mostly independent of the measured transport—no diffusion coefficient is fitted to experiment—and the MLIP is checked against DFT-NEB with ~10–15 meV MAE. If the mechanistic conclusion holds, the paper would provide a design principle for defect-enhanced WR anodes and would be of broad interest. However, the computational representation of the defect structures is the weakest link: the Wadsley-defect models are not the defect chemistry measured in MNO-800, and the cation-disorder conclusion rests on one selected motif rather than the SQS representative. These issues are load-bearing because the final attribution of the ~3× diffusivity enhancement to 'both defect types' depends on them.

major comments (4)
  1. [Results – Computational Analysis (Type-3 models, Fig. 5c,d)] The Wadsley-defect leg of the central claim is carried by T[3×5] and T[3×4-3×5] models with Mo only in tetrahedral sites and charge balance restored by O vacancies. This conflicts with the paper's own structural evidence for MNO-800: Mo K-edge EXAFS could not be fit to a model and indicates a markedly higher average Mo coordination than ~4.2; STEM-EDS maps Mo throughout [3×8] blocks and along shear planes. The Type-3 formulas also change the Mo:Nb ratio (1:15 and 2:27 vs 1:12). Because the O vacancy itself creates a 528 meV barrier and slows middle windows, the 'block size' effect on early window-site activation is entangled with the assumed charge-compensation mechanism. Please test enlarged-block models with Mo in octahedral/corner sites at the measured composition, or restrict the claim to O-vacancy-bearing block enlargements.
  2. [Results – Computational Analysis (Type-2 models, T[3×4-dis-5])] The cation-disorder mechanism is demonstrated with T[3×4-dis-1] (Mo at the block center), but the SQS model T[3×4-dis-5], described as representative of the disordered configuration, does not enhance diffusion relative to pristine T[3×4] and is not used in the site-occupation comparison. The paper should justify why dis-1 is the appropriate representation of disorder in MNO-800 (e.g., by matching the Mo-in-block distribution seen by STEM-EDS) or qualify the conclusion that cation disorder generally lowers the x threshold for activation of fast sites.
  3. [Conclusion and abstract] The abstract and conclusion attribute the improvement to 'both defect types,' but no model combines Wadsley block-size variation with cation disorder in a single simulation, and no model represents the actual integrated defect structure of MNO-800. The separate Type-2 and Type-3 results can show that each isolated motif is capable of shifting site occupancies; they do not establish that the effects are independent or additive in the real material. The manuscript itself notes after the STEM-EDS analysis that 'there remains ambiguity as to how each of these defects alters lithium diffusion.' A combined-defect calculation or an explicit discussion of interaction effects is needed before the concluding attribution can be accepted.
  4. [Methods – Molecular dynamics; Table 1] Table 1 is used to state that D is higher at 25% and 50% Li (order: dis-1 > 3×4-3×5 > 3×5 > pristine). These values are 300 K extrapolations from 2 ns MLIP-MD runs at 500–900 K, but no statistical uncertainty, trajectory convergence, or Arrhenius-fit quality is reported, and each concentration uses only three configurations. Please report error bars or confidence intervals so the ordering and factor differences in Table 1 can be evaluated.
minor comments (3)
  1. [Methods – Molecular dynamics / Results Table 1 paragraph] The text says concentration-dependent MD were 'performed at 300 K' while Methods states D at 300 K was extracted by Arrhenius extrapolation from 500–900 K. Please clarify which values are direct and which are extrapolated.
  2. [Fig. 4e] The three parabolic D(x) regions are described by widths and positions but no fitting procedure is given. State whether these are fits or guides and report the fit function/criteria.
  3. [General] Typos and notation: 'lowers' should be 'lower'; 'intensity of the peak reducing with annealing temperature' is ambiguous; non-integer oxygen stoichiometry in MoNb15O40.5 and Mo2Nb27O73.5 should be explicitly defined (per formula unit vs supercell).

Circularity Check

0 steps flagged

No significant circularity; central mechanism rests on DFT-trained MLIP plus independent ICI/electrochemistry; self-citations are methodological only.

full rationale

The central claim—that Wadsley defects and cation disorder independently advance Li occupation of fast Wh/Wv window sites to lower lithiation extents—is derived from CE/MC site-occupancy calculations using an MLIP trained on DFT geometries and energies, with MLIP-NEB validated against DFT-NEB (MAE 10.3/14.9 meV). No diffusion coefficient or site-occupation sequence is fitted to the experimental D(x) from ICI. The experimental 3X capacity-weighted diffusivity is measured from ICI and SAXS surface area, not inferred from the models. The paper's self-citations (refs 6, 13, 46, 51) concern nomenclature, a comparison metric, and prior applications of the ICI/parabola analysis; they are methodological precedents rather than load-bearing proofs, and the ICI method traces to an external source. The O-vacancy charge-compensation in the Type-3 Wadsley-defect models is an acknowledged modeling assumption that may not match MNO-800's measured Mo coordination, but this is a correctness/transferability limitation, not a circular reduction: the model's output is not defined in terms of the experimental result it explains.

Axiom & Free-Parameter Ledger

1 free parameters · 5 axioms · 1 invented entities

The central mechanistic claim rests on DFT/MLIP calculations. The only genuine fitted numeric is the Mo Hubbard U; the more consequential burden is a set of modeling assumptions: the MLIP's transferability, Arrhenius extrapolation, high-temperature CE occupation order, and the O-vacancy representation of Wadsley defects. The experimental half carries no free parameters of its own.

free parameters (1)
  • Hubbard U for Mo = 4.38 eV
    PBE+U parameter taken from MPRelaxSet; it controls redox and migration energetics in all DFT/MLIP calculations but is calibrated to Materials Project, not to the measured diffusivities of this paper.
axioms (5)
  • domain assumption PBE+U with U(Mo)=4.38 eV gives reliable formation and Li-migration energetics for MoNb12O33.
    All DFT, NEB, and MLIP training energies inherit this functional/parameter choice; no benchmark against other functionals or experimental migration barriers is provided.
  • domain assumption The MACE MLIP fine-tuned from MACE-MP0 reproduces DFT barriers well enough (MAE 10.3/14.9 meV on NEB paths) for MD diffusivities.
    MLIP-NEB validation is limited to T[3x4] and T[3x5]; MD diffusion and CE occupation use the same MLIP, including configurations not in the validation set.
  • domain assumption Li diffusivity follows Arrhenius behavior from 500-900 K down to 300 K.
    D values at 300 K are obtained by extrapolation of 500-900 K MD runs; no 300 K MD statistics or error estimates are reported.
  • domain assumption CE/MC at 2000 K finds ground-state Li site occupations that determine room-temperature occupation order.
    Site-occupation sequence is computed from high-temperature Monte Carlo and used to interpret the experimental D(x) parabola positions.
  • ad hoc to paper Wadsley-defect models T[3x5] and T[3x4-3x5] with Mo only in tetrahedral sites and compensating O vacancies are representative of defect-rich MNO-800.
    This stoichiometry/charge-balance choice is made to fix composition, but EXAFS and STEM suggest Mo occupies non-tetrahedral sites in MNO-800; the O-vacancy representation is not experimentally validated.
invented entities (1)
  • Oxygen vacancies in enlarged Wadsley blocks no independent evidence
    purpose: Charge-balance device to keep Nb/Mo ratio ~12 when block size increases while Mo remains tetrahedral
    No direct O-vacancy measurement; XPS shows only fully oxidized metals and EXAFS/STEM indicate Mo is not confined to tetrahedra, so charge compensation in the real material is unresolved.

pith-pipeline@v1.3.0-alltime-deepseek · 21865 in / 15559 out tokens · 165610 ms · 2026-08-03T22:35:34.650340+00:00 · methodology

0 comments
read the original abstract

Wadsley-Roth (WR) niobates have emerged as high-rate anode materials that can combine rapid ionic diffusion with good electronic conductivity. WR compounds have been defect-enhanced by limited annealing, however, such materials often contain multiple types of defects. In particular, both Wadsley defects (variable block size) and transition metal disorder have the potential to modify transport rates, however the corresponding effects are not well understood mechanistically. Here, MoNb12O33 (MNO) was calcined at two different temperatures to compare a defect-rich condition (MNO-800) with a proximal order-rich condition (MNO-900) as assessed through XRD, XANES, EXAFS, and STEM characterizations. Galvanostatically cycled lithium half cells of MNO-800 exhibited additional capacity (307 mAh/g at 0.1C, 4.66% higher) and improved high-rate capacity of 200 mAhg-1 at 10C. ICI-based overpotential analysis identified solid state diffusion as the dominant rate limiting process where MNO-800 correspondingly exhibited ~3X faster capacity-weighted diffusivity. A machine-learning interatomic potential was trained to density functional theory and then applied with molecular dynamics (MLIP-MD) to examine the possible roles of Wadsley defects and transition metal disorder. For both defect-types, Li was found to populate and activate fast diffusion paths from window sites at lower extents of lithiation as compared to the order-rich model.

Figures

Figures reproduced from arXiv: 2511.09521 by Aaron Hegler, Christopher Sutton, CJ Sturgill, Coby S. Collins, Hsin-Yun Joy Chao, Iva Milisavljevic, Manish Kumar, Morgan Stefik, Nima Karimitari, Santosh Kiran Balijepalli, Scott Misture.

Figure 1
Figure 1. Figure 1: XRD and Rietveld refinement results for (a) MNO [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: HAADF STEM data for MNO-800(a-f, left side) and MNO-900(g-l, right side). Each region is presented with block overlays (a,c,e,g,i,k) indicating the blocks colored according to size and without overlays (b,d,f,h,j,l). Specific bonds (m-q) were associated with the radial distribution functions (r,s). Single-image RDFs (r) from calculated from the indicated insets. Average RDFs for each sample type (s) were c… view at source ↗
Figure 5
Figure 5. Figure 5: Example model structures analyzed, including T[3 [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

9 extracted references · 3 canonical work pages

  1. [1]

    # !$# !"#$%& %!!&#%'!!()# &

    XRD and Rietveld refinement results for (a) MNO-800 and (b) MNO-900 samples (excluded sample holder peak at ~22.5). (c) Nb K edge and (d) Mo K edge radial distribution function of FT EXAFS oscillations, all phase-corrected. Kinetically limited crystallization is well known to result in transient (non-equilibrium) defect-rich states. The early 1970’s studi...

  2. [2]

    carbon black,

    Diffusion coefficient (in m2s-1) at different concentrations of Li for the considered structural models. Li Concentrations Pristine T[3×4] T[3×4-dis-1] T[3×5] T[3×4-3×5] 25% 2.00x10-12 3.07x10-10 3.08x10-11 3.26x10-11 50% 3.12x10-13 3.01x10-11 1.30x10-11 1.57x10-11 75% 3.17x10-12 5.93x10-13 2.40x10-12 8.95x10-12 The faster diffusion is confirmed with MD s...

  3. [4]

    Ganesha Operating Manual

    (61) Joensen, K. Ganesha Operating Manual. https://smif.pratt.duke.edu/capabilities/55527 (accessed 2025-06-17). (62) NIST X-Ray Form Factor, Atten. Scatt. Tables Form Page. https://physics.nist.gov/PhysRefData/FFast/html/form.html (accessed 2025-06-17). (63) Scattering Length Density Calculator. https://www.ncnr.nist.gov/resources/sldcalc.html (accessed ...

  4. [516]

    (43) Xiang, W.; Chen, M.; Zhou, X.; Chen, J.; Huang, H.; Sun, Z.; Lu, Y.; Zhang, G.; Wen, X.; Li, W

    https://doi.org/10.1038/ncomms1527. (43) Xiang, W.; Chen, M.; Zhou, X.; Chen, J.; Huang, H.; Sun, Z.; Lu, Y.; Zhang, G.; Wen, X.; Li, W. Highly Enforced Rate Capability of a Graphite Anode via Interphase Chemistry Tailoring Based on an Electrolyte Additive. J. Phys. Chem. Lett. 2022, 13 (23), 5151–5159. https://doi.org/10.1021/acs.jpclett.2c01183. (44) Ch...

  5. [592]

    (31) Driscoll, E

    https://doi.org/10.1038/s41467-019-08490-w. (31) Driscoll, E. H.; Green, A.; Fortes, D.; Howard, C.; Driscoll, L. L.; Kendrick, E.; Greaves, C.; Slater, P. R. A High-Power 4 × 4: Crystallographic and Electrochemical Insights into a Novel Wadsley–Roth Anode Nb9Ti1.5W1.5O30. Chem. Commun. 2024, 60 (73), 10001–10004. https://doi.org/10.1039/D4CC02844A. (32) ...

  6. [2024]

    (52) Zunger, A.; Wei, S.-H.; Ferreira, L

    https://doi.org/10.1039/D4TA06041E. (52) Zunger, A.; Wei, S.-H.; Ferreira, L. G.; Bernard, J. E. Special Quasirandom Structures. Phys. Rev. Lett. 1990, 65 (3), 353–356. https://doi.org/10.1103/PhysRevLett.65.353. (53) van de Walle, A.; Tiwary, P.; de Jong, M.; Olmsted, D. L.; Asta, M.; Dick, A.; Shin, D.; Wang, Y.; Chen, L.-Q.; Liu, Z.-K. Efficient Stocha...

  7. [2025]

    (75) Waters, M

    https://doi.org/10.48550/arXiv.2401.00096. (75) Waters, M. J.; Rondinelli, J. M. Energy Contour Exploration with Potentiostatic Kinematics. J. Phys. Condens. Matter 2021, 33 (44), 445901. https://doi.org/10.1088/1361-648X/ac1af0. (76) Sanchez, J. M.; Ducastelle, F.; Gratias, D. Generalized Cluster Description of Multicomponent Systems. Phys. Stat. Mech. I...

  8. [2289]

    (45) Geng, Z.; Thiringer, T.; Lacey, M

    https://doi.org/10.1038/s41467-023-37989-6. (45) Geng, Z.; Thiringer, T.; Lacey, M. J. Intermittent Current Interruption Method for Commercial Lithium-Ion Batteries Aging Characterization. IEEE Trans. Transp. Electrification 2022, 8 (2), 2985–2995. https://doi.org/10.1109/TTE.2021.3125418. (46) Kumar, M.; Muhit, M. A. A.; Sturgill, C.; Karimitari, N.; Bar...

  9. [6507]

    (34) Iijima, S

    https://doi.org/10.1038/s41467-025-61646-9. (34) Iijima, S. Direct Observation of Lattice Defects in H -Nb 2 O 5 by High Resolution Electron Microscopy. Acta Crystallogr. Sect. A 1973, 29 (1), 18–24. https://doi.org/10.1107/S0567739473000045. (35) Anderson, J. S. On Infinitely Adaptive Structures. J. Chem. Soc. Dalton Trans. 1973, No. 10, 1107–1115. https...