Interatomic interactions for Mo and W based on the low-order moments of the density of states

On the binding of nanometric hydrogen-helium clusters in tungsten

In this work we developed an embedded atom method potential for large scale atomistic simulations in the ternary tungsten-hydrogen-helium (W-H-He) system, focusing on applications in the fusion research domain. Following available ab initio data, the potential reproduces key interactions between H, He and point defects in W and utilizes the most recent potential for matrix W. The potential is applied to assess the thermal stability of various H-He complexes of sizes too large for ab initio techniques. The results show that the dissociation of H-He clusters stabilized by vacancies will occur primarily by emission of hydrogen atoms and then by break-up of V-He complexes, indicating that H-He interaction does influence the release of hydrogen.

On the mobility of vacancy clusters in reduced activation steels: an atomistic study in the Fe-Cr-W model alloy

Reduced activation steels are considered as structural materials for future fusion reactors. Besides iron and the main alloying element chromium, these steels contain other minor alloying elements, typically tungsten, vanadium and tantalum. In this work we study the impact of chromium and tungsten, being major alloying elements of ferritic Fe-Cr-W-based steels, on the stability and mobility of vacancy defects, typically formed under irradiation in collision cascades. For this purpose, we perform ab initio calculations, develop a many-body interatomic potential (EAM formalism) for large-scale calculations, validate the potential and apply it using an atomistic kinetic Monte Carlo method to characterize the lifetime and diffusivity of vacancy clusters. To distinguish the role of Cr and W we perform atomistic kinetic Monte Carlo simulations in Fe-Cr, Fe-W and Fe-Cr-W alloys. Within the limitation of transferability of the potentials it is found that both Cr and W enhance the diffusivity of vacancy clusters, while only W strongly reduces their lifetime. The cluster lifetime reduction increases with W concentration and saturates at about 1-2 at.%. The obtained results imply that W acts as an efficient 'breaker' of small migrating vacancy clusters and therefore the short-term annealing process of cascade debris is modified by the presence of W, even in small concentrations.

Molybdenum at high pressure and temperature: melting from another solid phase

The Gibbs free energies of bcc and fcc Mo are calculated from first principles in the quasiharmonic approximation in the pressure range from 350 to 850 GPa at room temperatures up to 7500 K. It is found that Mo, stable in the bcc phase at low temperatures, has lower free energy in the fcc structure than in the bcc phase at elevated temperatures. Our density-functional-theory-based molecular dynamics simulations demonstrate that fcc melts at higher than bcc temperatures above 1.5 Mbar. Our calculated melting temperatures and bcc-fcc boundary are consistent with the Mo Hugoniot sound speed measurements. We find that melting occurs at temperatures significantly above the bcc-fcc boundary. This suggests an explanation of the recent diamond anvil cell experiments, which find a phase boundary in the vicinity of our extrapolated bcc-fcc boundary.

Angular-dependent matrix potentials for fast molecular-dynamics simulations of transition metals

The significance of the part played by the angular-dependent components of forces associated with d-d bonding between atoms in a transition metal has long remained a subject of debate. While almost all the large-scale molecular dynamics simulations of collision cascades and radiation damage in transition metals and alloys are currently performed using spherically symmetric many-body potentials, density functional calculations exhibit a highly anisotropic pattern of charge density deformation in and around the core of interstitial atom defects. This paper describes a fast second-order matrix recursion-based algorithm for including effects of angular anisotropy of d-d bonds in a large-scale molecular dynamics simulation.

High-pressure melting of molybdenum

The melting curve of the body-centered cubic (bcc) phase of Mo has been determined for a wide pressure range using both direct ab initio molecular dynamics simulations of melting as well as a phenomenological theory of melting. These two methods show very good agreement. The simulations are based on density functional theory within the generalized gradient approximation. Our calculated equation of state of bcc Mo is in excellent agreement with experimental data. However, our melting curve is substantially higher than the one determined in diamond anvil cell experiments up to a pressure of 100 GPa. An explanation is suggested for this discrepancy.