Atomic diffusion in liquid gallium and gallium-nickel alloys probed by quasielastic neutron scattering and molecular dynamic simulations

Abstract

The atomic mobility in liquid pure gallium and a gallium-nickel alloy with 2 at% of nickel is studied experimentally by incoherent quasielastic neutron scattering. The integral diffusion coefficients for all-atom diffusion are derived from the experimental data at different temperatures. DFT-based ab-initio molecular dynamics (MD) is used to find numerically the diffusion coefficient of liquid gallium at different temperatures, and numerical theory results well agree with the experimental findings at temperatures below 500 K. Machine learning force fields derived from ab-initio molecular dynamics (AIMD) overestimate within a small 6% error the diffusion coefficient of pure gallium within the genuine AIMD. However, they better agree with experiment for pure gallium and enable the numerical finding of the diffusion coefficient of nickel in the considered melted alloy along with the diffusion coefficient of gallium and integral diffusion coefficient, that agrees with the corresponding experimental values within the error bars. The temperature dependence of the gallium diffusion coefficient follows the Arrhenius law experimentally for all studied temperatures and below 500 K also in the numerical simulations. However, can be well described alternatively by an Einstein-Stokes dependence with the metallic liquid viscosity following the Arrhenius law, especially for the MD simulation results at all studied temperatures. Moreover, a novel variant of the excess entropy scaling theory rationalized our findings for gallium diffusion. Obtained values of the Arrhenius activation energies are profoundly different in the competing theoretical descriptions, which is explained by different temperature-dependent prefactors in the corresponding theories. The diffusion coefficient of gallium is significantly reduced (at the same temperature) in a melted alloy with natural nickel, even at a tiny 2 at% concentration of nickel, as compared with its pure gallium value. This highly surprising behavior contradicts the existing excess entropy scaling theories and opens a venue for further research.