Pith. sign in

REVIEW 2 major objections 6 minor 79 references

Segregation, ordering, and precipitation in WTaV-based concentrated refractory alloys

T0 review · 2 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash

Pith's one-line read CrV precipitates in WTaCrV are semicoherent, not coherent: simulations and atom probe data agree.

desk verdict A solid, well-validated tabGAP study that likely explains the CrV/TaW vs CrTa/WV discrepancy via semicoherent interfaces; the main weakness is that the crossover curve leans on tabGAP-only semicoherent interfacial energies. read the letter →

arxiv 2412.13750 v2 pith:IWDBMQS5 submitted 2024-12-18 cond-mat.mtrl-sci

classification cond-mat.mtrl-sci
keywords refractoryhigh-entropyalloysmachine-learnedinteratomicpotentialtabGAPshort-rangeordersegregationsemicoherentinterfacesprecipitationatomprobetomography
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper argues that the CrV-rich precipitates observed in irradiated WTaCrV refractory alloys are not coherently embedded in the matrix; instead they form semicoherent bcc-to-bcc interfaces carrying networks of misfit dislocations. The authors build a machine-learned interatomic potential for the W-Ta-Cr-V system, validate it against density functional theory, and use hybrid Monte Carlo/molecular dynamics simulations plus interface-energy calculations to show that coherent CrV precipitates in a TaW matrix are thermodynamically unstable at all sizes. Above roughly 2.5 to 3 nm in slab thickness or radius, a semicoherent CrV/TaW configuration becomes the lowest-energy structure, matching the 2 to 5 nm precipitates seen in atom probe tomography. This resolves a discrepancy in which coherent-lattice simulations predicted CrTa/WV short-range order while experiments and an earlier cluster-expansion model see CrV precipitates.

What carries the argument

The load-bearing machinery is the tabGAP machine-learned interatomic potential for the full W-Ta-Cr-V composition space, trained on about 13,483 density functional theory structures. To explain precipitation, the authors combine the interfacial superlattice formation energy E_f/N = 2Sσ/N + ζ + $ΔE_f^{{AB}}$/N with coherency strain energies computed from Eq. 10, obtaining coherent interface energies with both DFT and tabGAP and semicoherent interface energies with tabGAP for the mismatched CrV/TaW system. The crossover curves in Fig. 14 are generated from these ingredients and are what place semicoherent CrV/TaW below coherent CrTa/VW above 2.5 to 3 nm.

What would settle it

A DFT calculation of a semicoherent CrV/TaW interface with a misfit dislocation network, or direct high-resolution transmission electron microscopy of the interfaces in the irradiated W38Ta36Cr15V11 films, would settle whether the precipitates are semicoherent. If DFT placed semicoherent CrV/TaW above coherent CrTa/VW at all sizes up to 5 nm, or if imaging showed no misfit dislocation network while precipitates were still present, the central claim would fail.

Watch

Extended reading notes

Core claim

The central claim is that the experimentally observed CrV-rich precipitates in WTaCrV form semicoherent bcc-to-bcc interfaces with the surrounding TaW-rich matrix, because coherent precipitates are unstable due to excessive lattice mismatch. The evidence is a crossover analysis: coherent CrV three-dimensional precipitates in a TaW matrix have positive formation energy of about 0.17 eV/atom at large radius, while semicoherent CrV/TaW becomes the lowest-energy configuration above 2.5 to 3 nm in slab thickness or 3 nm in radius, in the size range of the 2 to 5 nm precipitates observed by atom probe tomography. A Monte Carlo/molecular dynamics simulation starting from a random W38Ta36Cr15V11 solution with a compressed embedded precipitate produces a semicoherent precipitate with roughly 70% Cr and 30% V composition and a network of misfit dislocations, quantitatively matching the earlier experiments. The paper therefore concludes that earlier coherent-lattice simulations predicted the wrong short-range order because they excluded semicoherent interfaces.

Load-bearing premise

The conclusion that the observed precipitates must be semicoherent rests on semicoherent interfacial energies computed only with the machine-learned potential, without a DFT cross-check, and on model systems that are ideal equiatomic B2 binaries rather than the real roughly 70% Cr and 30% V precipitates.

Editorial extensions

If this is right

  • If the central claim is correct, the CrV-rich precipitates seen in irradiated WTaCrV are surrounded by misfit dislocation networks, and any simulation that fixes atoms on a single coherent lattice will miss the observed short-range order.
  • The crossover at 2.5 to 3 nm implies that small coherent CrV clusters are not thermodynamically stable, explaining why the experimental precipitates are 2 to 5 nm and why no smaller precipitates are observed.
  • The uniform segregation rule that small atoms (Cr, V) prefer compressed regions and large atoms (Ta, Nb) prefer tensile regions extends across WTaV, WTaCrV, and MoNbTaVW, offering a simple size-based guideline for other refractory alloys.
  • Grain boundary segregation in WTaV approaches TaV2 composition, suggesting grain boundaries can act as nucleation sites for Laves phases.
  • In WTaCrV, Cr segregation to defects is accompanied by V to form stable CrV mixtures, making it the most segregation-resistant of the three alloys studied.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the crossover depends on the equiatomic B2 idealization, then for the real roughly 70% Cr and 30% V precipitates the crossover size and interfacial energy may shift; extending the calculation to non-equiatomic compositions would be a direct test.
  • The semicoherent interfacial energies are currently computed only with the machine-learned potential; a DFT calculation of a misfit-dislocation-bearing CrV/TaW supercell would either confirm or revise the 2.5 to 3 nm crossover.
  • The paper leaves open how semicoherent precipitates form; the proposed pathway through segregation to interstitial dislocation loops could be tested by kinetic simulations of CrV nucleation at interstitial clusters.
  • If semicoherent interfaces are indeed the stable form, mechanical properties of irradiated WTaCrV such as hardening may be controlled more by misfit dislocation networks at precipitate interfaces than by short-range order alone, which would affect alloy design considerations.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 6 minor

Summary. The manuscript develops a tabulated Gaussian approximation potential (tabGAP) for the W–Ta–Cr–V system, validates it extensively against DFT for binary formation energies, energy-volume curves, point defects, surfaces, and grain boundaries, and then uses hybrid MC/MD simulations to study short-range order and segregation in WTaV, WTaCrV, and MoNbTaVW. Atom probe tomography of W38Ta36Cr15V11 thin films is used to characterize compositional striations in the as-deposited state and globular CrV-rich precipitates after irradiation. To explain these observations, the authors compute coherent and semicoherent interface energetics between B2-ordered CrV/TaW and CrTa/VW binaries, combining DFT coherency-strain energies with tabGAP interface energies. The resulting crossover curves in Fig. 14 indicate that coherent CrV precipitates in TaW are unstable at all sizes, and that semicoherent CrV/TaW becomes competitive with coherent CrTa/VW at thicknesses or radii around 2.5–3 nm, from which the authors conclude that the experimentally observed CrV-rich precipitates are semicoherent. The paper also draws general conclusions about segregation trends linked to pure-element surface, grain-boundary, and atomic-size properties.

Significance. If the central claim holds, the paper resolves a real discrepancy: coherent-lattice MC/MD simulations and cluster-expansion models disagree on whether CrV or CrTa ordering is preferred in WTaCrV, and the paper explains this by showing that the coherent lattice constraint selects CrTa/VW while the experimentally relevant state is semicoherent CrV/TaW. The work is strengthened by the breadth of DFT validation of the tabGAP, the iterative training on short-range-order-relevant structures, the cross-validation of two independent tabGAPs for WTaV, and the fact that the APT observations were not used as fitting input for any parameter. The conclusion that coherent CrV 3D precipitates are unstable is robust because the DFT coherency-strain energy (0.248 eV/atom) alone exceeds the combined B2 formation-energy gain of CrV and TaW (−0.08 eV/atom). The main uncertainty is quantitative: the semicoherent branch of Fig. 14, which is load-bearing for the semicoherency conclusion, rests on tabGAP-only interface energies for misfit-dislocation-bearing interfaces, and the crossover is computed for ideal equiatomic B2 binaries rather than the measured precipitate composition.

major comments (2)
  1. [III D 2, Fig. 14, Table IV] The semicoherent branch of Fig. 14 is the load-bearing input for the statement that CrV-rich precipitates in WTaCrV must be semicoherent, and it rests entirely on tabGAP interfacial energies: the text states that 'Semicoherent interfacial energies are always from tabGAP.' The tabGAP is validated against DFT for bulk phases, point defects, surfaces, and grain boundaries, but not for the 1/2<111> misfit dislocation networks that define these semicoherent interfaces. In Eq. (9) with ζ = 0, the interface term for a spherical precipitate is 3vσ/r; with v ≈ 12.2 Å^3 and σ in the 0.05–0.15 eV/Å^2 range reported, an uncertainty of only ±0.03 eV/Å^2 changes this term by roughly 0.04 eV/atom at r = 3 nm, which is larger than the 0.03 eV/atom asymptotic energy difference between the CrV/TaW and CrTa/VW combinations. Such an uncertainty can move the crossover outside the observed 2–5 nm precipitate-size range and thereby undermine the 'must be semicoherent' inference. I request either a DFT benchmark for at least one semicoherent interface (for example, the smallest commensurate 9:10 CrV/TaW supercell) or a systematic sensitivity analysis over σ demonstrating that the crossover location is stable to ±0.03 eV/Å^2.
  2. [III D 2, Figs. 13–14] The crossover calculation is performed for ideal equiatomic B2 CrV/TaW and CrTa/VW binaries, as the paper itself acknowledges with 'the theoretical crossovers are for ideal equiatomic B2 binaries.' The APT-measured precipitates, however, are approximately 70% Cr and 30% V, and the matrix is W38Ta36Cr15V11 rather than equiatomic. Off-equiatomic compositions will change both the bulk formation-energy input to Eq. (9) (the −0.08 vs −0.05 eV/atom values) and the interfacial energies. A calculation at the measured precipitate composition, or at least an explicit demonstration that the crossover is insensitive to composition, is needed before the observed precipitates can be quantitatively identified with the calculated crossover.
minor comments (6)
  1. [II F] Equation (9) is derived for a periodic superlattice with two interfaces, but the 3D spherical-cluster application is described only as 'we use Eq. 9.' Please state explicitly how the 3D formula is obtained (one interface, S = 4πr^2, N = V/v) and whether the factor 2 in Eq. (9) is dropped for the spherical case.
  2. [Table IV] Several tabulated interfacial energies are negative (for example, σ(100) = −0.055 eV/Å^2 for CrV/TaW tabGAP). A sentence explaining that these negative values reflect favorable chemical mixing across the interface rather than an unphysical negative interface energy would prevent misinterpretation.
  3. [Fig. 14 caption] The label 'DFT/tabGAP pred.' is used in the figure legend but is defined only in the main text; please define it in the caption as well, and state which quantities are from DFT and which are from tabGAP.
  4. [III C 3] The phrase 'volemetric strains' should read 'volumetric strains.'
  5. [III D 2] The wording 'must be semicoherent' is stronger than the evidence supports, since APT cannot resolve the misfit dislocation network and the paper notes that no high-resolution imaging was possible. 'Are predicted to be semicoherent' would better match the computational nature of the evidence.
  6. [References] Reference [43] is cited as 'W-Ta-Cr-V tabGAP: Potential files, training data, and input (2025)' without a repository URL or DOI; please provide the full citation so the potential and training data are actually accessible.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the central semicoherent-precipitate claim is a computed prediction compared against APT, not a re-statement of fitted inputs.

full rationale

The derivation chain is not circular. The WTaCrV tabGAP is fitted to DFT energies, forces, and virials for 13,483 structures, and the paper validates it against independent DFT data such as MC/MD frame energies (RMSE 2.54 meV/atom), vacancy migration and formation energies, binary formation energies, and grain-boundary energies. The semicoherent interfacial energies that set the crossover in Fig. 14 are predictions from this DFT-trained potential, not parameters fitted to the APT-observed precipitates; the APT data are used only as an external comparison. The conclusion that CrV-rich precipitates in WTaCrV must be semicoherent follows from the computed free-energy crossover, with the paper explicitly flagging that semicoherent interfacial energies are tabGAP-only and that the model uses ideal equiatomic B2 binaries. Those are accuracy and robustness limitations, not circular reductions. Self-citations to Refs. [11] and [19] support the tabGAP methodology and the MoNbTaVW potential, but they are not load-bearing for the new WTaCrV central claim, which is independently validated against DFT and experiment.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new physical entities. Its central conclusion depends on the accuracy of the ML potential and on the idealization of interfacial model systems. The tabGAP is fitted to DFT data, and the semicoherent interface energies are not independently checked against DFT. No ad hoc parameters were introduced to force the semicoherent result; the crossover arises from computed strain and interfacial energies.

free parameters (3)
  • tabGAP hyperparameters (sparse points, cutoffs, regularization) = M_2b=20, M_EAM=20, M_3b=300; r_cut=5 Å (2b/EAM) and 4.1 Å (3b); regularization 2 to 5 meV/atom energies, 0.1 eV/Å forces
    Chosen by hand, identical to prior MoNbTaVW work [19]; they control the ML potential used for all segregation and interface energies, but are not fitted to the experimental observation being explained.
  • Chemical potentials for W, Ta, Cr, V in equiatomic WTaCrV = mu_Cr=-9.388 eV, mu_Ta=-11.724 eV, mu_V=-9.033 eV, mu_W=-13.029 eV
    Obtained by least-squares fit to 4000 substitution energies (Eq. 7); used for vacancy and interstitial formation energies, not for the semicoherent precipitate conclusion.
  • ZBL repulsive potential parameters zeta_i and eta_i = Six parameters per element pair
    Prefitted to all-electron DFT dimer data; part of the base potential, not central to the semicoherent claim.
assumptions (5)
  • domain assumption Neglect of spin polarization in DFT training data and in the tabGAP is acceptable for WTaCrV alloys and for CrV/TaW interfaces relevant to the conclusions.
    Appendix A shows magnetic moments vanish in B2 CrV except at very large volumes, and spin-polarized 3D coherency strain energy is only slightly lower (0.230 vs 0.248 eV/atom), preserving conclusions. However, the potential does not reproduce the magnetic Cr ground state, and stretched Cr-rich environments could be affected.
  • domain assumption Equiatomic B2-ordered binary alloys CrV, TaW, CrTa, and VW are adequate model systems for the experimentally observed non-equiatomic CrV-rich precipitates and the TaW-rich matrix.
    The paper states 'the theoretical crossovers are for ideal equiatomic B2 binaries, while the concentration profiles in Fig. 13 show more chemical variations.' The crossover sizes (2 to 3 nm) are therefore approximate.
  • domain assumption The tabGAP extrapolates reliably to semicoherent interfaces with misfit dislocations, for which no DFT energies are computed.
    Semicoherent interfacial energies in Table IV and Fig. 14 are tabGAP-only because the cells are too large for DFT. The potential is validated for coherent interfaces and bulk properties, but not for dislocation-bearing interfaces.
  • domain assumption Hybrid MC/MD simulations at 300 K for 100,000 MD steps produce representative segregation states around defects, even though full equilibrium is not reached.
    The paper says 'the systems do not reach equilibrium and the chemical order in bulk regions far from the defect are not fully optimised'; the segregation conclusions assume near-defect regions are converged.
  • domain assumption PBE-GGA DFT with PAW potentials is a sufficient reference for the alloy thermodynamics and interface energies.
    Standard for this field, but PBE has known errors for magnetic Cr; the paper partially addresses this in Appendix A.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Segregation, ordering, and precipitation in WTaV-based concentrated refractory alloys." pith.science (2026). https://pith.science/paper/IWDBMQS5

@misc{pith2026241213750,
  author       = {Pith},
  title        = {Pith review of: Segregation, ordering, and precipitation in WTaV-based concentrated refractory alloys},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/IWDBMQS5}},
  note         = {Machine review of arXiv:2412.13750}
}
read the original abstract

Tungsten-based low-activation high-entropy alloys are possible candidates for next-generation fusion reactors due to their exceptional tolerance to irradiation, thermal loads, and stress. We develop an accurate and efficient machine-learned interatomic potential for the W-Ta-Cr-V system and use it in hybrid Monte Carlo molecular dynamics simulations of ordering and segregation to all common types of defects in WTaCrV. The predictions are compared to atom probe tomography analysis of segregation and precipitation in WTaCrV thin films. By also considering two other alloys, WTaV and MoNbTaVW, we are able to draw general conclusions about preferred segregation in refractory alloys and the reasons behind it, guiding future alloy design and elucidating experimental observations. We show that the experimentally observed CrV precipitates in WTaCrV form semicoherent bcc-to-bcc interfaces with the surrounding matrix, as coherent precipitates are not thermodynamically stable due to excessive lattice mismatch. The predictions from simulations align well with our atom probe tomography analysis as well as previous experimental observations.

Figures

Figures reproduced from arXiv: 2412.13750 by the authors.

Figure 1
Figure 1. FIG. 1: Formation energies of relaxed binary alloys from DFT compared with the W–Ta–Cr–V tabGAP. [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Energy evolution during a MC/MD simulation [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. FIG. 4: Distributions of (a) single-vacancy migration [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figures from the paper (12 more)
Figure 5
Figure 5. Figure 5: FIG. 5: Relaxed grain boundary energies of all six pure [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6: Short-range order parameters for first and second-nearest neighbour shells, 1NN and 2NN, as functions of [PITH_FULL_IMAGE:figures/full_fig_p009_6.png]
Figure 8
Figure 8. Figure 8: FIG. 8: Surface segregation illustrated as [PITH_FULL_IMAGE:figures/full_fig_p010_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9: Radial concentrations and snapshots of the [PITH_FULL_IMAGE:figures/full_fig_p011_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10: Segregation around edge dislocation lines in the three alloys. (a) shows the edge dislocation dipole [PITH_FULL_IMAGE:figures/full_fig_p012_10.png]
Figure 11
Figure 11. Figure 11: FIG. 11: Segregation around screw dislocation cores in [PITH_FULL_IMAGE:figures/full_fig_p012_11.png]
Figure 12
Figure 12. Figure 12: FIG. 12: Radial concentration profiles showing the [PITH_FULL_IMAGE:figures/full_fig_p012_12.png]
Figure 13
Figure 13. Figure 13: FIG. 13: APT analyses of (a, b) as-deposited W [PITH_FULL_IMAGE:figures/full_fig_p013_13.png]
Figure 14
Figure 14. Figure 14: FIG. 14: Formation energies per atom of coherent and semicoherent 3D and 2D periodic interfacial structures of [PITH_FULL_IMAGE:figures/full_fig_p014_14.png]
Figure 15
Figure 15. Figure 15: FIG. 15: Semicoherent interfaces between CrV and WTa after MC/MD relaxation, (a): (1 0 0)/(1 0 0) and (b): [PITH_FULL_IMAGE:figures/full_fig_p016_15.png]
Figure 16
Figure 16. Figure 16: FIG. 16: Energy versus volume for pure Cr and B2 [PITH_FULL_IMAGE:figures/full_fig_p018_16.png]
Figure 17
Figure 17. Figure 17: FIG. 17: Cross-validation of the new W–Ta–Cr–V [PITH_FULL_IMAGE:figures/full_fig_p019_17.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

79 extracted references · 66 canonical work pages

  1. [1]

    [13]) as well as the A15 and Laves C14, C15, and C36 intermetallic phases

    Binary alloys Figure 1 shows formation energies of binary alloys, in- cluding BCC-like random structures and ordered phases (from Ref. [13]) as well as the A15 and Laves C14, C15, and C36 intermetallic phases. The formation energies of random alloys are the averages of three different 1024- atom systems. All structures are fully relaxed (positions and cel...

  2. [2]

    Figure 3 shows the energy as a function of number of MC swap attempts for a tabGAP MC/MD simulation at 300 K of a 128-atom initially random WTaCrV system

    WTaCrV alloys We used a hybrid MC/MD simulation as a test case for how accurate the tabGAP is for short-range ordered structures of WTaCrV. Figure 3 shows the energy as a function of number of MC swap attempts for a tabGAP MC/MD simulation at 300 K of a 128-atom initially random WTaCrV system. The sample contains one va- cancy to trigger some local segreg...

  3. [3]

    Pure-element properties Even though the primary purpose of the WTaCrV tabGAP is to simulate alloys, it can also be used in simu- lations of the pure metals. More importantly, some fun- damental pure-element properties can be used to either 8 0.0 0.5 1.0 1.5 2.0 2.5 3.0 Vacancy migration energy (eV) 0 200 400Count Moving element:(a) Cr T a V W 1.5 2.0 2.5 ...

  4. [4]

    8: Surface segregation illustrated as layer-by-layer concentration profiles for the three alloys, with snapshots of each surface after MC/MD relaxation shown as insets

    Grain boundaries Figure 7 shows snapshots of the nanocrystalline WTaV, WTaCrV, and MoNbTaVW systems after the 0.2 0.4 0.6 0.8 1.0 WTaV (100) WTaV (110) WTaV (111)T a V W 0.2 0.4 0.6 0.8 1.0 WTaCrV (100) WTaCrV (110) WTaCrV (111)Cr T a V W 1 3 5 7 9 0.2 0.4 0.6 0.8 1.0 MoNbTaVW (100) 1 3 5 7 9 MoNbTaVW (110) 1 3 5 7 9 MoNbTaVW (111)Mo Nb T a V W Subsurface...

  5. [5]

    Based on Tab

    Surfaces and voids The size difference and significant variation in surface energies of the elements in the alloys is expected to lead to segregation to surfaces and voids in order to release stress and minimise the energy. Based on Tab. II, pure Nb has the lowest surface energy followed by V and Ta. Fig. 8 shows the results of the MC/MD simulations of 11...

  6. [6]

    10 shows the results of MC/MD relaxation of 1/2⟨1 1 1⟩(1 1 0) edge dislocation dipole structures of the three alloys

    Edge and screw dislocations Fig. 10 shows the results of MC/MD relaxation of 1/2⟨1 1 1⟩(1 1 0) edge dislocation dipole structures of the three alloys. We observe that the dislocation cores at- tract significant segregation, which can be correlated with the strain field induced by the edge dislocations. Fig. 10 also shows the volemetric strains around the ...

  7. [7]

    12 shows the results from MC/MD relaxation of both interstitial- and vacancy-type 1/2⟨1 1 1⟩prismatic dislocation loops

    Dislocation loops Fig. 12 shows the results from MC/MD relaxation of both interstitial- and vacancy-type 1/2⟨1 1 1⟩prismatic dislocation loops. The results are summarised in plots of the radial concentration profiles from the centre of the loop. The same trends as for infinite edge dislocation lines are apparent: V (and Cr in WTaCrV) segregate to the comp...

  8. [8]

    6 and in Ref

    Atom probe tomography analysis There is a remarkable discrepancy between the pre- dicted short-range order in single crystal WTaCrV by ML potentials (both here in Fig. 6 and in Ref. [20]) and that produced by the CE model in Ref. [5]. The latter predicts formation of almost pure binary clusters of Cr (∼65%) and V (∼30%) in a surrounding WTa-rich ma- trix,...

Show all 79 references
  1. [9]

    DFT/tabGAP pred

    Coherent and semicoherent interfaces Fig. 13 shows that both the CrV-rich precipitates in the irradiated sample and the CrV- and TaW-rich layers in the as-deposited sample are only 2–5 nm in length. To analyse the interface energetics, we choose as a model system equiatomic Cr...

  2. [10]

    O. N. Senkov, D. B. Miracle, K. J. Chaput, and J.- P. Couzinie, Development and exploration of refractory high entropy alloys—A review, Journal of Materials Re- search33, 3092 (2018)

  3. [11]

    Cheng, J

    Z. Cheng, J. Sun, X. Gao, Y. Wang, J. Cui, T. Wang, and H. Chang, Irradiation effects in high-entropy al- loys and their applications, Journal of Alloys and Com- pounds930, 166768 (2023)

  4. [12]

    O. N. Senkov, G. B. Wilks, D. B. Miracle, C. P. Chuang, and P. K. Liaw, Refractory high-entropy alloys, Inter- metallics18, 1758 (2010)

  5. [13]

    O. N. Senkov, J. M. Scott, S. V. Senkova, D. B. Miracle, and C. F. Woodward, Microstructure and room temper- ature properties of a high-entropy TaNbHfZrTi alloy, Journal of Alloys and Compounds509, 6043 (2011)

  6. [14]

    El-Atwani, N

    O. El-Atwani, N. Li, M. Li, A. Devaraj, J. K. S. Baldwin, M. M. Schneider, D. Sobieraj, J. S. Wr´ obel, D. Nguyen-Manh, S. A. Maloy, and E. Martinez, Out- standing radiation resistance of tungsten-based high- entropy alloys, Science Advances5, eaav2002 (2019)

  7. [15]

    El Atwani, H

    O. El Atwani, H. T. Vo, M. A. Tunes, C. Lee, A. Al- varado, N. Krienke, J. D. Poplawsky, A. A. Kohnert, J. Gigax, W.-Y. Chen, M. Li, Y. Q. Wang, J. S. Wr´ obel, D. Nguyen-Manh, J. K. S. Baldwin, O. U. Tukac, E. Aydogan, S. Fensin, and E. Martinez, A quinary WTaCrVHf nanocrysta...

  8. [16]

    G. Wei, J. Byggm¨ astar, J. Cui, K. Nordlund, J. Ren, and F. Djurabekova, Revealing the critical role of vana- dium in radiation damage of tungsten-based alloys, Acta Materialia274, 119991 (2024)

  9. [17]

    S. Yin, J. Ding, M. Asta, and R. O. Ritchie, Ab initio modeling of the energy landscape for screw dislocations in body-centered cubic high-entropy alloys, npj Comput Mater6, 1 (2020)

  10. [18]

    Kostiuchenko, F

    T. Kostiuchenko, F. K¨ ormann, J. Neugebauer, and A. Shapeev, Impact of lattice relaxations on phase tran- sitions in a high-entropy alloy studied by machine- learning potentials, npj Comput Mater5, 1 (2019)

  11. [19]

    X.-G. Li, C. Chen, H. Zheng, Y. Zuo, and S. P. Ong, Complex strengthening mechanisms in the NbMoTaW multi-principal element alloy, npj Comput Mater6, 1 (2020)

  12. [20]

    Byggm¨ astar, K

    J. Byggm¨ astar, K. Nordlund, and F. Djurabekova, Modeling refractory high-entropy alloys with effi- cient machine-learned interatomic potentials: Defects and segregation, Phys. Rev. B104, 104101 (2021), arXiv:2106.03369

  13. [21]

    Fern´ andez-Caballero, J

    A. Fern´ andez-Caballero, J. S. Wr´ obel, P. M. Mummery, and D. Nguyen-Manh, Short-Range Order in High En- tropy Alloys: Theoretical Formulation and Application to Mo-Nb-Ta-V-W System, J. Phase Equilib. Diffus.38, 391 (2017)

  14. [22]

    Sobieraj, J

    D. Sobieraj, J. S. Wr´ obel, T. Rygier, K. J. Kurzyd lowski, O. E. Atwani, A. Devaraj, E. M. Saez, and D. Nguyen-Manh, Chemical short-range order in derivative Cr–Ta–Ti–V–W high entropy alloys from the first-principles thermodynamic study, Phys. Chem. Chem. Phys.22, 23929 (2020)

  15. [23]

    N. C. Smith, T.-c. Liu, Y. Xia, and C. Wolverton, Competition between long- and short-range order in size-mismatched medium-entropy alloys, Acta Materi- alia277, 120199 (2024)

  16. [24]

    Fantin, A

    A. Fantin, A. M. Manzoni, H. Springer, R. D. Ka- machali, and R. Maaß, Local lattice distortions and chemical short-range order in MoNbTaW, Materials Re- search Letters12, 346 (2024)

  17. [25]

    Mishin, Machine-learning interatomic potentials for materials science, Acta Materialia214, 116980 (2021)

    Y. Mishin, Machine-learning interatomic potentials for materials science, Acta Materialia214, 116980 (2021). 19 0 1000 2000 3000 MC swap attempts 1429 1428 1427 1426 1425 1424 Energy (eV) RMSE: 1.52 meV/atom tabGAP-WTaCrV DFT tabGAP-MoNbT aVW 0 1000 2000 3000 MC swap attempts ...

  18. [26]

    V. L. Deringer, A. P. Bart´ ok, N. Bernstein, D. M. Wilkins, M. Ceriotti, and G. Cs´ anyi, Gaussian Process Regression for Materials and Molecules, Chem. Rev. 121, 10073 (2021)

  19. [27]

    Behler, Perspective: Machine learning potentials for atomistic simulations, J

    J. Behler, Perspective: Machine learning potentials for atomistic simulations, J. Chem. Phys.145, 170901 (2016)

  20. [28]

    Byggm¨ astar, K

    J. Byggm¨ astar, K. Nordlund, and F. Djurabekova, Sim- ple machine-learned interatomic potentials for complex alloys, Phys. Rev. Materials6, 083801 (2022)

  21. [29]

    S. Lyu, W. Li, Y. Xia, Y. Chen, and A. H. W. Ngan, Effects of chemical randomness on strength contribu- tors and dislocation behaviors in a bcc multiprincipal element alloy, Phys. Rev. Mater.7, 073602 (2023)

  22. [30]

    Kresse and J

    G. Kresse and J. Hafner, Ab initio molecular dynamics for liquid metals, Phys. Rev. B47, 558 (1993)

  23. [31]

    Kresse and J

    G. Kresse and J. Hafner, Ab initio molecular- dynamics simulation of the liquid-metal–amorphous- semiconductor transition in germanium, Phys. Rev. B 49, 14251 (1994)

  24. [32]

    Kresse and J

    G. Kresse and J. Furthm¨ uller, Efficiency of ab-initio to- tal energy calculations for metals and semiconductors using a plane-wave basis set, Computational Materials Science6, 15 (1996)

  25. [33]

    Kresse and J

    G. Kresse and J. Furthm¨ uller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B54, 11169 (1996)

  26. [34]

    P. E. Bl¨ ochl, Projector augmented-wave method, Phys. Rev. B50, 17953 (1994)

  27. [35]

    Kresse and D

    G. Kresse and D. Joubert, From ultrasoft pseudopoten- tials to the projector augmented-wave method, Physical Review B59, 1758 (1999)

  28. [36]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized Gradient Approximation Made Simple, Phys. Rev. Lett. 77, 3865 (1996)

  29. [37]

    H. J. Monkhorst and J. D. Pack, Special points for Brillouin-zone integrations, Phys. Rev. B13, 5188 (1976)

  30. [38]

    Methfessel and A

    M. Methfessel and A. T. Paxton, High-precision sam- pling for Brillouin-zone integration in metals, Phys. Rev. B40, 3616 (1989)

  31. [39]

    Zhao, Defect properties in a VTaCrW equiatomic high entropy alloy (HEA) with the body centered cubic (bcc) structure, Journal of Materials Science & Tech- nology44, 133 (2020)

    S. Zhao, Defect properties in a VTaCrW equiatomic high entropy alloy (HEA) with the body centered cubic (bcc) structure, Journal of Materials Science & Tech- nology44, 133 (2020)

  32. [40]

    El-Atwani, A

    O. El-Atwani, A. Alvarado, K. Unal, S. Fensin, J. Hinks, G. Greaves, J. K. S. Baldwin, S. Maloy, and E. Martinez, Helium Implantation Damage Resistance in Nanocrys- talline W-Ta-V-Cr High Entropy Alloys, Materials To- day Energy , 100599 (2020)

  33. [41]

    Hafner, D

    R. Hafner, D. Spiˇ s´ ak, R. Lorenz, and J. Hafner, Mag- netic ground state of Cr in density-functional theory, Phys. Rev. B65, 184432 (2002)

  34. [42]

    A. P. Bart´ ok, M. C. Payne, R. Kondor, and G. Cs´ anyi, Gaussian Approximation Potentials: The Accuracy of Quantum Mechanics, without the Electrons, Physical Review Letters104, 10.1103/PhysRevLett.104.136403 (2010)

  35. [43]

    Glielmo, C

    A. Glielmo, C. Zeni, and A. De Vita, Efficient nonpara- metric n -body force fields from machine learning, Phys- ical Review B97, 10.1103/PhysRevB.97.184307 (2018)

  36. [44]

    Vandermause, S

    J. Vandermause, S. B. Torrisi, S. Batzner, Y. Xie, L. Sun, A. M. Kolpak, and B. Kozinsky, On-the-fly ac- tive learning of interpretable Bayesian force fields for atomistic rare events, npj Computational Materials6, 1 (2020)

  37. [45]

    A. P. Bart´ ok and G. Cs´ anyi, Gaussian approximation potentials: A brief tutorial introduction, International Journal of Quantum Chemistry115, 1051 (2015)

  38. [46]

    M. S. Daw and M. I. Baskes, Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals, Phys. Rev. B29, 6443 (1984)

  39. [47]

    J. F. Ziegler, J. P. Biersack, and U. Littmarck, The Stopping and Range of Ions in Matter, inTreatise on Heavy-Ion Science(Pergamon, New York, 1985) pp. 93– 129

  40. [48]

    Nordlund, N

    K. Nordlund, N. Runeberg, and D. Sundholm, Repul- sive interatomic potentials calculated using Hartree- Fock and density-functional theory methods, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms132, 45 (1997)

  41. [49]

    Byggm¨ astar, A

    J. Byggm¨ astar, A. Hamedani, K. Nordlund, and F. Djurabekova, Machine-learning interatomic poten- tial for radiation damage and defects in tungsten, Phys. Rev. B100, 144105 (2019)

  42. [50]

    Byggm¨ astar, K

    J. Byggm¨ astar, K. Nordlund, and F. Djurabekova, Gaussian approximation potentials for body-centered- cubic transition metals, Phys. Rev. Materials4, 093802 (2020)

  43. [51]

    Starikov, P

    S. Starikov, P. Grigorev, and P. A. T. Olsson, Angular- dependent interatomic potential for large-scale atom- istic simulation of W-Mo-Nb ternary alloys, Computa- tional Materials Science233, 112734 (2024). 20

  44. [52]

    Byggm¨ astar, W-Ta-Cr-V tabGAP: Potential files, training data, and input (2025)

    J. Byggm¨ astar, W-Ta-Cr-V tabGAP: Potential files, training data, and input (2025)

  45. [53]

    Jeffries, F

    J. Jeffries, F. Abdeljawad, S. Mathaudhu, E. Marquis, and E. Martinez, Prediction of defect properties in con- centrated solid solutions using a Langmuir-like model, Phys. Rev. Mater.9, 033803 (2025)

  46. [54]

    J. S. Wrobel, A. Zhong, A. Goryaeva, D. Nguyen-Manh, A. E. Poisvert, M. Athenes, and M. C. Marinica, Data- driven sampling for predicting composition-dependent thermoelasticity, melting and defect anharmonic contri- butions in high-entropy alloys, Phys. Rev. B , submitted (2025)

  47. [55]

    Schuler, M

    T. Schuler, M. Nastar, K. Li, and C.-C. Fu, Towards ac- curate thermodynamics from random energy sampling, Acta Materialia276, 120074 (2024)

  48. [56]

    K. Li, T. Schuler, C.-C. Fu, and M. Nastar, Vacancy formation free energy in concentrated alloys: Equilib- rium vs. random sampling, Acta Materialia281, 120355 (2024)

  49. [57]

    Stukowski, Visualization and analysis of atomistic simulation data with OVITO–the Open Visualization Tool, Modelling Simul

    A. Stukowski, Visualization and analysis of atomistic simulation data with OVITO–the Open Visualization Tool, Modelling Simul. Mater. Sci. Eng.18, 015012 (2010)

  50. [58]

    Henkelman, B

    G. Henkelman, B. P. Uberuaga, and H. J´ onsson, A climbing image nudged elastic band method for find- ing saddle points and minimum energy paths, J. Chem. Phys.113, 9901 (2000)

  51. [59]

    Zheng, X.-G

    H. Zheng, X.-G. Li, R. Tran, C. Chen, M. Horton, D. Winston, K. A. Persson, and S. P. Ong, Grain bound- ary properties of elemental metals, Acta Materialia186, 40 (2020)

  52. [60]

    A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolin- tineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimp- ton, LAMMPS - a flexible simulation tool for particle- based ...

  53. [61]

    Hirel, Atomsk: A tool for manipulating and convert- ing atomic data files, Computer Physics Communica- tions197, 212 (2015)

    P. Hirel, Atomsk: A tool for manipulating and convert- ing atomic data files, Computer Physics Communica- tions197, 212 (2015)

  54. [62]

    A. H. Larsen, J. J. Mortensen, J. Blomqvist, I. E. Castelli, R. Christensen, Marcin Du lak, J. Friis, M. N. Groves, B. Hammer, C. Hargus, E. D. Hermes, P. C. Jennings, P. B. Jensen, J. Kermode, J. R. Kitchin, E. L. Kolsbjerg, J. Kubal, Kristen Kaasbjerg, S. Lysgaard, J. B. Mar...

  55. [63]

    Thompson, D

    K. Thompson, D. Lawrence, D. J. Larson, J. D. Ol- son, T. F. Kelly, and B. Gorman, In situ site-specific specimen preparation for atom probe tomography, Ul- tramicroscopy107, 131 (2007)

  56. [64]

    O. C. Hellman, J. A. Vandenbroucke, J. R¨ using, D. Isheim, and D. N. Seidman, Analysis of Three- dimensional Atom-probe Data by the Proximity His- togram, Microscopy and Microanalysis6, 437 (2000)

  57. [65]

    Ozoli¸ nˇ s, C

    V. Ozoli¸ nˇ s, C. Wolverton, and A. Zunger, Effects of an- harmonic strain on the phase stability of epitaxial films and superlattices: Applications to noble metals, Phys. Rev. B57, 4816 (1998)

  58. [66]

    Wolverton, V

    C. Wolverton, V. Ozolins, and A. Zunger, Short-range- order types in binary alloys: A reflection of coherent phase stability, J. Phys.: Condens. Matter12, 2749 (2000)

  59. [67]

    R. Tran, Z. Xu, B. Radhakrishnan, D. Winston, W. Sun, K. A. Persson, and S. P. Ong, Surface energies of ele- mental crystals, Sci Data3, 160080 (2016)

  60. [68]

    Singh, A

    P. Singh, A. V. Smirnov, and D. D. Johnson, Atomic short-range order and incipient long-range order in high- entropy alloys, Phys. Rev. B91, 224204 (2015)

  61. [69]

    Toda-Caraballo, J

    I. Toda-Caraballo, J. S. Wr´ obel, D. Nguyen-Manh, P. P´ erez, and P. E. J. Rivera-D ´ ıaz-del-Castillo, Simu- lation and Modeling in High Entropy Alloys, JOM69, 2137 (2017)

  62. [70]

    M. Koskenniemi, Vacancy diffusion in high-entropy al- loys by a kinetic Monte Carlo methodology based on collective variable-driven hyperdynamics and machine learning, University of Helsinki (2023), mSc thesis, http://hdl.handle.net/10138/567555

  63. [71]

    Ma, Unusual dislocation behavior in high-entropy al- loys, Scripta Materialia181, 127 (2020)

    E. Ma, Unusual dislocation behavior in high-entropy al- loys, Scripta Materialia181, 127 (2020)

  64. [72]

    Maresca and W

    F. Maresca and W. A. Curtin, Mechanistic origin of high strength in refractory BCC high entropy alloys up to 1900K, Acta Materialia182, 235 (2020)

  65. [73]

    D. Utt, S. Lee, Y. Xing, H. Jeong, A. Stukowski, S. H. Oh, G. Dehm, and K. Albe, The origin of jerky dislo- cation motion in high-entropy alloys, Nat Commun13, 4777 (2022)

  66. [74]

    Z. T. Kloenne, J.-P. Couzini´ e, M. Heczko, R. Gr¨ oger, G. B. Viswanathan, W. A. T. Clark, and H. L. Fraser, On thebcc/B2 interface structure in a refractory high entropy alloy, Scripta Materialia223, 115071 (2023)

  67. [75]

    Aksoy, M

    D. Aksoy, M. J. McCarthy, I. Geiger, D. Apelian, H. Hahn, E. J. Lavernia, J. Luo, H. Xin, and T. J. Rupert, Chemical order transitions within extended in- terfacial segregation zones in NbMoTaW, Journal of Ap- plied Physics132, 235302 (2022)

  68. [76]

    Y. Shi, Z. Jiang, T. Xia, Z. Wang, J. Wu, X. Cao, and K. Zhu, Deuterium retention and desorption behavior of W-Ta-Cr-V high entropy alloy, Journal of Nuclear Materials568, 153897 (2022)

  69. [77]

    Y. Shi, Z. Jiang, T. Xia, W. Zhang, P. Yang, X. Ren, M. Wang, L. Liang, X. Cao, and K. Zhu, Helium dif- fusion and bubble evolution in single-phase tungsten- based W-Ta-Cr-V complex concentrated alloy, Journal of Nuclear Materials578, 154335 (2023)

  70. [78]

    Kalita, I

    D. Kalita, I. J´ o´ zwik, L. Kurpaska, Y. Zhang, K. Mulewska, W. Chrominski, J. O’Connell, Y. Ge, W. L. Boldman, P. D. Rack, Y. Wang, W. J. Weber, and J. Jagielski, The microstructure and He+ ion ir- radiation behavior of novel low-activation W-Ta-Cr-V refractory high entropy ...

  71. [79]

    D. A. Pankhurst, D. Nguyen-Manh, and D. G. Petti- for, Electronic origin of structural trends across early transition-metal disilicides: Anomalous behavior of CrSi2, Phys. Rev. B69, 075113 (2004)

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.