pith. sign in

arxiv: 2605.23612 · v2 · pith:V76QV5W5new · submitted 2026-05-22 · ❄️ cond-mat.mtrl-sci

First-principles transition-state tensorial cluster expansion of vacancy diffusion in Ta-W beyond the kinetically-resolved activation approximation

Jacob Jeffries , Brianna Sebastian-Olazabal , Enrique Martinez This is my paper

Pith reviewed 2026-05-25 03:49 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords vacancy diffusionTa-W alloystensorial cluster expansionkinetic Monte Carlomigration barriersfirst-principlescomposition dependenceactivation energy
0
0 comments X

The pith

A tensorial cluster expansion learns full environment-dependent migration barriers for vacancy diffusion in Ta-W without the KRA approximation.

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

The paper develops a first-principles framework that represents migration barriers via a tensorial cluster expansion trained on DFT and nudged-elastic-band data, avoiding the kinetically resolved activation approximation. The barriers are then used directly in on-lattice kinetic Monte Carlo simulations of vacancy diffusion. Applied to the Ta-W binary system, the simulations produce a maximum in apparent activation energy at intermediate compositions. This maximum arises from a crossover between regimes of solute trapping and percolated low-barrier pathways. The work supplies a scalable method for incorporating explicit transition-state energetics into mesoscale kinetic models of chemically complex alloys.

Core claim

By training a tensorial cluster expansion on transition-state energies computed from density-functional theory and nudged-elastic-band calculations, the method parameterizes the full local-environment dependence of vacancy migration barriers; when these barriers drive kinetic Monte Carlo simulations of Ta-W, they produce nontrivial composition dependence in which apparent activation energy reaches a maximum near intermediate compositions because of the crossover between solute-trapping and percolated low-barrier transport pathways.

What carries the argument

Tensorial cluster expansion of transition states that encodes the full local atomic environment dependence of migration barriers for direct use in kinetic Monte Carlo.

If this is right

  • The framework supplies a general route for integrating first-principles transition-state energetics into mesoscale kinetic simulations.
  • It enables predictive multiscale modeling of diffusion in chemically complex materials.
  • It reveals emergent transport phenomena such as composition-dependent maxima in activation energy.
  • Reliable extrapolation of barriers to unseen environments becomes possible in KMC without reduced kinetic approximations.

Where Pith is reading between the lines

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

  • The same expansion approach could be applied to other binary refractory alloys to test whether similar trapping-to-percolation crossovers appear.
  • Alloy design for high-temperature service might deliberately target or avoid the intermediate-composition regime where the activation-energy maximum occurs.
  • Extension to ternary or higher-order systems could show whether the tensorial representation remains tractable when more chemical species are present.

Load-bearing premise

The tensorial cluster expansion can accurately reproduce the full environment dependence of DFT-computed migration barriers without the KRA approximation.

What would settle it

Experimental measurements of vacancy diffusion coefficients in Ta-W alloys that show no maximum in apparent activation energy near intermediate compositions would falsify the predicted crossover behavior.

Figures

Figures reproduced from arXiv: 2605.23612 by Brianna Sebastian-Olazabal, Enrique Martinez, Jacob Jeffries.

Figure 2
Figure 2. Figure 2: Vacancy migration on the bcc multi-lattice with the ini [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Example minimum energy path computed from NEB and [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Transition energy E ‡ as a function of average energy (E ′ + E ◦ )/2 (left), and energy barrier E ‡ − E ◦ as a function of energy difference E ′ − E ◦ (right) from DFT + NEB calculations of vacancy hopping in Ta-W. This means that, at least for our dataset, the energy barrier cannot be written as a function of solely E′ − E◦ , regardless of how those energies are computed, therefore the KRA barrier express… view at source ↗
Figure 6
Figure 6. Figure 6: Pairwise interaction coefficients ε (top) and three-body interaction coefficients ζ (bottom). For the pairwise coefficients, each heatmap is labeled with an interaction order n, with each cell corresponding to a given α-β pair. For the three-body coefficients, each heatmap is labeled with an interaction order n as well as an atom type γ, with each cell corresponding to a given α-β pair. In both cases, the … view at source ↗
Figure 7
Figure 7. Figure 7: Parity plot of energy barrier predictions. [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Distributions of traversed barriers within the KMC trajectories for each tested composition and temperature. [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Predicted vacancy diffusivity as a function of composi [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Predicted activation energy for vacancy diffusion as a [PITH_FULL_IMAGE:figures/full_fig_p010_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Ta excess near vacancy ΓTa throughout the KMC sim￾ulations. Each cell represents a simulation at a fixed temperature and nominal composition. The purple shades represent attraction of the vacancy to Ta, i.e. excess Ta, while the green shades rep￾resent repulsion, i.e. depleted Ta, with shade corresponding to magnitude of excess/depletion. The results for excess Ta near the vacancy are con￾sistent with the… view at source ↗
read the original abstract

Predicting diffusion in chemically complex alloys remains challenging due to the strong dependence of migration barriers on local atomic environments. Migration barriers computed using density functional theory and nudged elastic band calculations are represented via a tensorial cluster expansion including transition states and deployed in on-lattice kinetic Monte Carlo simulations. Applied to the Ta-W system, the framework captures nontrivial composition-dependent diffusion behavior arising from a crossover between solute trapping and percolated low-barrier transport pathways, yielding a maximum in the apparent activation energy near intermediate compositions. This approach establishes a general and scalable route for integrating first-principles transition-state energetics into mesoscale kinetic simulations, enabling predictive multiscale modeling of diffusion in chemically complex materials and providing a pathway for uncovering emergent transport phenomena.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a first-principles framework that represents environment-dependent vacancy migration barriers in Ta-W via a tensorial cluster expansion fitted to DFT/NEB data, deploys the expansion directly in on-lattice kinetic Monte Carlo without the kinetically resolved activation (KRA) approximation, and reports that the resulting composition-dependent diffusivities exhibit a maximum in apparent activation energy at intermediate compositions arising from a crossover between solute trapping and percolated low-barrier pathways.

Significance. If the tensorial cluster expansion faithfully reproduces the full set of DFT barriers across local environments, the work supplies a scalable route for embedding first-principles transition-state energetics into mesoscale simulations of diffusion in chemically complex alloys. The explicit avoidance of the KRA approximation and the demonstration of emergent composition dependence constitute a concrete advance over reduced models, with potential applicability to other alloy systems.

major comments (2)
  1. [Methods/Results] The central claim that the tensorial cluster expansion captures the full environment dependence of migration barriers (thereby enabling reliable KMC extrapolation) requires quantitative validation; the manuscript should report RMSE, cross-validation scores, or parity plots comparing CE predictions to held-out DFT/NEB barriers for representative local environments (Methods or Results section).
  2. [Results] The reported maximum in apparent activation energy at intermediate compositions is load-bearing for the crossover interpretation; the KMC results should include statistical uncertainties from multiple runs or sensitivity tests to the CE fit parameters to establish that the maximum is robust rather than an artifact of finite sampling or fitting noise.
minor comments (2)
  1. [Methods] Notation for the tensorial cluster expansion (e.g., how the transition-state tensors are indexed and symmetrized) should be defined explicitly with an equation or table early in the Methods section for reproducibility.
  2. [Figures] Figure captions for the composition-dependent diffusivity and activation-energy plots should state the temperature range, vacancy concentration, and number of KMC steps used.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and insightful comments, which have helped us strengthen the validation and statistical rigor of the manuscript. We address each major comment below and have revised the manuscript accordingly.

read point-by-point responses

  1. Referee: [Methods/Results] The central claim that the tensorial cluster expansion captures the full environment dependence of migration barriers (thereby enabling reliable KMC extrapolation) requires quantitative validation; the manuscript should report RMSE, cross-validation scores, or parity plots comparing CE predictions to held-out DFT/NEB barriers for representative local environments (Methods or Results section).

    Authors: We agree that explicit quantitative validation is required to support the central claim. In the revised manuscript we have added a dedicated paragraph in the Methods section that reports the root-mean-square error of the tensorial CE fit to the full DFT/NEB training set, the leave-one-out cross-validation score, and parity plots for a held-out subset of barriers spanning representative local environments at several compositions. These metrics demonstrate that the expansion reproduces the DFT barriers to within 15 meV on average, thereby justifying direct use in KMC without the KRA approximation. revision: yes

  2. Referee: [Results] The reported maximum in apparent activation energy at intermediate compositions is load-bearing for the crossover interpretation; the KMC results should include statistical uncertainties from multiple runs or sensitivity tests to the CE fit parameters to establish that the maximum is robust rather than an artifact of finite sampling or fitting noise.

    Authors: We acknowledge the need to demonstrate robustness of the composition-dependent maximum. The revised Results section now includes error bars obtained from five independent KMC trajectories per composition (different random seeds) and a sensitivity analysis in which the CE coefficients were perturbed within their fitting uncertainties before re-running the KMC. Both the location and magnitude of the maximum remain statistically significant across these tests, confirming that the feature arises from the physical crossover between solute trapping and percolated low-barrier pathways rather than sampling or fitting artifacts. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper computes migration barriers directly from DFT+NEB calculations for local environments in Ta-W, fits a tensorial cluster expansion to those computed values, and deploys the resulting model in on-lattice KMC to obtain composition-dependent diffusivities. The reported maximum in apparent activation energy and crossover between trapping and percolated pathways are emergent outputs of the KMC simulations on the fitted landscape, not definitions or statistical re-statements of the input DFT data. No load-bearing self-citation, uniqueness theorem, or ansatz smuggling is present in the abstract or described framework; the method is explicitly constructed to avoid the KRA approximation by direct environment-dependent fitting. The chain is therefore self-contained against external first-principles benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No specific free parameters, axioms, or invented entities are identifiable from the abstract alone.

pith-pipeline@v0.9.0 · 5698 in / 1035 out tokens · 54876 ms · 2026-05-25T03:49:49.430919+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.