A reactive interatomic potential based on an analytic bond-order scheme is developed for the ternary system W-C-H. The model combines Brenner's hydrocarbon potential with parameter sets for W-W, W-C and W-H interactions and is adjusted to materials properties of reference structures with different local atomic coordinations including tungsten carbide, W-H molecules as well as H dissolved in bulk W. The potential has been tested in various scenarios, like surface, defect, and melting properties, none of which were considered in the fitting. The intended area of application is simulations of hydrogen and hydrocarbon interactions with tungsten, that have a crucial role in fusion reactor plasma-walls. Furthermore, this study shows that the angular dependent bond-order scheme can be extended to second-nearest neighbor interactions, which are relevant in body-centered cubic metals. Moreover, it provides a possibly general route for modeling metal carbides.

The experimentally observed large difference in the depths of hydrogen and helium clusters formed in tungsten still lacks a fundamental explanation. Using density functional theory calculations, molecular dynamics simulations, and kinetic Monte Carlo calculations, we show that the fundamental mechanism behind the different clustering depths is significantly different behaviors of interstitial H and He atoms in W: H–H states are unstable for small interatomic distances whereas He–He states are strongly bound.

Stainless steels found in real-world applications usually have some C content in the base Fe-Cr alloy, resulting in hard and dislocation-pinning carbides-Fe3C (cementite) and Cr23C6-being present in the finished steel product. The higher complexity of the steel microstructure has implications, for example, for the elastic properties and the evolution of defects such as Frenkel pairs and dislocations. This makes it necessary to re-evaluate the effects of basic radiation phenomena and not simply to rely on results obtained from purely metallic Fe-Cr alloys. In this report, an analytical interatomic potential parameterization in the Abell-Brenner-Tersoff form for the entire Fe-Cr-C system is presented to enable such calculations. The potential reproduces, for example, the lattice parameter(s), formation energies and elastic properties of the principal Fe and Cr carbides (Fe3C, Fe5C2, Fe7C3, Cr3C2, Cr7C3, Cr23C6), the Fe-Cr mixing energy curve, formation energies of simple C point defects in Fe and Cr, and the martensite lattice anisotropy, with fair to excellent agreement with empirical results. Tests of the predictive power of the potential show, for example, that Fe-Cr nanowires and bulk samples become elastically stiffer with increasing Cr and C concentrations. High-concentration nanowires also fracture at shorter relative elongations than wires made of pure Fe. Also, tests with Fe3C inclusions show that these act as obstacles for edge dislocations moving through otherwise pure Fe.

No abstract

Analytical bond-order potentials for beryllium, beryllium carbide and beryllium hydride are presented. The reactive nature of the formalism makes the potentials suitable for simulations of non-equilibrium processes such as plasma-wall interactions in fusion reactors. The Be and Be-C potentials were fitted to ab initio calculations as well as to experimental data of several different atomic configurations and Be-H molecule and defect data were used in determining the Be-H parameter set. Among other tests, sputtering, melting and quenching simulations were performed in order to check the transferability of the potentials. The antifluorite Be(2)C structure is well described by the Be-C potential and the hydrocarbon interactions are modelled by the established Brenner potentials.

No abstract

An analytical bond-order interatomic potential has been developed for the iron-carbon system for use in molecular-dynamics and Monte Carlo simulations. The potential has been successfully fitted to cementite and Hägg carbide, which are most important crystalline polytypes among the many known metastable iron carbide phases. Predicted properties of other carbides and the simplest point defects are in good to reasonable agreement with available data from experiments and density-functional theory calculations. The potential correctly describes melting and recrystallization of cementite, making it useful for simulation of steels. We show that they correctly describe the metastability of cementite and can be used to model carbide growth and dissolution.

No abstract