A molecular dynamics examination of the relationship between self-diffusion and viscosity in liquid metals

Non-Arrhenius behaviour of nickel self-diffusion in liquid Ni77Si23

Nickel self-diffusion was measured for a Ni₇₇Si₂₃alloy in the liquid state over a temperature range of about 400 K through quasielastic neutron scattering. At the two lowest temperature points the derived diffusion coefficients deviate from a high-temperature Arrhenius-type behaviour and indicate a change in dynamics above the liquidus temperature. A fit with a power-law temperature dependence as predicted by the mode coupling theory for the liquid to glass transition can describe the diffusion coefficients quite well over the whole measured temperature range. The obtained results agree with predictions from a classical molecular dynamics (MD)-simulation, which evidenced an increasing glass forming ability with increasing silicon content. A crossover to a super-Arrhenius behaviour was reported for metallic glass formers above the liquidus temperature and the here investigated NiSi alloy demonstrates the same signature.

Probing the self-diffusion process in Aluminium

An in-depth understanding of the diffusion process in liquid metals is a key to design and engineering new high-performance materials. In this study, using molecular dynamics simulations supplemented with the embedded atom potential, we investigate/compare the self-diffusion process in liquid Aluminium. To understand the self-diffusion process, we analyse the radial distribution functions, velocity distributions, mean square displacements, and self-diffusion coefficients at various temperatures well above the melting temperature of Aluminium in the temperature range of 1000 K to 1800 K. As a key result, in both the [Formula: see text] and [Formula: see text] phases, the self-diffusion coefficients show a non-linear variation with rise in temperatures in the range of 1000 K to 1200 K. From 1300 K to 1800 K, the self-diffusion coefficients increase more or less monotonically with rise in temperature. We found that a higher temperature in the range of 1300 K to 1800 K leads to a greater self-diffusion coefficient, suggesting the more violent movement of the atoms around their equilibrium positions. The results presented in this work can help to understand the differences in the self-diffusion process in the technologically relevant Al phases.

Revisiting the Stokes-Einstein relation for glass-forming melts

Molecular dynamics simulations of Ni36Zr64, Cu65Zr35 and Ni80Al20 were carried out over a broad range of temperature (900-3000 K) to investigate the Stokes-Einstein (SE) relation for glass-forming melts. Our results reproduce experimental structural and transport properties. Results show that the breakdown temperature of the SE relation (TSE) equals the dynamical crossover temperature (TA) and both are roughly twice the glass-transition temperature (Tg) for the three glass-forming melts (TSE = TA ≈ 2.0Tg). The product of the individual component self-diffusion coefficient and viscosity Dαη can be roughly regarded as a constant at the transition zone (a small temperature range around TSE) in which the temperature behaviors of self-diffusion coefficient and viscosity switch from high-temperature Arrhenius to a low-temperature VFT behavior. Below TSE, the decoupling of component diffusion coefficients was found. In particular, the decoupling of component diffusion coefficients can be ascribed to the decoupling of the partial pair structural correlation of components, which can be clearly reflected by the intersection of the high-temperature and low-temperature behaviors of the ratio between the partial pair correlation entropy of components (Sβ2/Sα2). Furthermore, the ratio between the partial pair correlation entropy of components may be used to predict the validity of the SE relation, in the absence of both transport coefficients and atomic coordinates.

Heterogeneous diffusion, viscosity, and the Stokes-Einstein relation in binary liquids

We investigate the origin of the breakdown of the Stokes-Einstein relation (SER) between diffusivity and viscosity in undercooled melts. A binary Lennard-Jones system, as a model for a metallic melt, is studied by molecular dynamics. A weak breakdown at high temperatures can be understood from the collectivization of motion, seen in the isotope effect. The strong breakdown at lower temperatures is connected to an increase in dynamic heterogeneity. On relevant time scales some particles diffuse much faster than the average or than predicted by the SER. The van Hove self-correlation function allows one to unambiguously identify slow particles. Their diffusivity is even less than predicted by the SER. The time span of these particles being slow particles, before their first conversion to be a fast one, is larger than the decay time of the stress correlation. The contribution of the slow particles to the viscosity rises rapidly upon cooling. Not only the diffusion but also the viscosity shows a dynamically heterogeneous scenario. We can define a "slow" viscosity. The SER is recovered as the relation between slow diffusivity and slow viscosity.