× Updated! Potentials that share interactions are now listed as related models.

2023--Sikorski-E-L-Cusentino-M-A-McCarthy-M-J-et-al--W-Zr-C

Citation: E.L. Sikorski, M.A. Cusentino, M.J. McCarthy, J. Tranchida, M.A. Wood, and A.P. Thompson (2023), "Machine learned interatomic potential for dispersion strengthened plasma facing components", The Journal of Chemical Physics, 158(11), 114101. DOI: 10.1063/5.0135269.
Abstract: Tungsten (W) is a material of choice for the divertor material due to its high melting temperature, thermal conductivity, and sputtering threshold. However, W has a very high brittle-to-ductile transition temperature, and at fusion reactor temperatures (≥1000 K), it may undergo recrystallization and grain growth. Dispersion-strengthening W with zirconium carbide (ZrC) can improve ductility and limit grain growth, but much of the effects of the dispersoids on microstructural evolution and thermomechanical properties at high temperatures are still unknown. We present a machine learned Spectral Neighbor Analysis Potential for W-ZrC that can now be used to study these materials. In order to construct a potential suitable for large-scale atomistic simulations at fusion reactor temperatures, it is necessary to train on ab initio data generated for a diverse set of structures, chemical environments, and temperatures. Further accuracy and stability tests of the potential were achieved using objective functions for both material properties and high temperature stability. Validation of lattice parameters, surface energies, bulk moduli, and thermal expansion is confirmed on the optimized potential. Tensile tests of W/ZrC bicrystals show that although the W(110)-ZrC(111) C-terminated bicrystal has the highest ultimate tensile strength (UTS) at room temperature, observed strength decreases with increasing temperature. At 2500 K, the terminating C layer diffuses into the W, resulting in a weaker W-Zr interface. Meanwhile, the W(110)-ZrC(111) Zr-terminated bicrystal has the highest UTS at 2500 K.

Notes: This potential was optimized for bulk and surface structures for both W and ZrC. No optimization was performed on pure C structures, and no physical performance should be expected for pure C simulations. We expect the potential to perform well in the temperature range of 300 - 2500K. Primary optimization was performed on bulk modulus, (100) and (110) surface energies, thermal expansion, and several stability checks detailed in the publication.

LAMMPS pair_style hybrid/overlay zbl snap (2023--Sikorski-E-L--W-Zr-C--LAMMPS--ipr1)
See Computed Properties
Notes: These files were provided by Ember Sikorski on 28 March 2023. This potential can be used by adding "include in.pot_snapWZrC" to a LAMMPS input script. The zbl parameters in in.pot_snapWZrC must be included to achieve the accuracy and performance described in the publication.
File(s):
Date Created: October 5, 2010 | Last updated: April 12, 2023