Sign Crossover in All Maxwell-Stefan Diffusivities for Molten Salt LiF-BeF2: A Molecular Dynamics Study

Applying Green-Kubo formalism and equilibrium molecular dynamics (MD) simulations, we have studied for the first time the dynamic correlation, Onsager coefficients, and Maxwell-Stefan (MS) diffusivities of molten salt LiF-BeF2, which is a potential candidate for a coolant in a high temperature reactor. We observe an unusual composition dependence and strikingly a crossover in sign for all the MS diffusivities at a composition of around 7% of LiF where the MS diffusivity between cation-anion pair (Đ(BeF) and Đ(LiF)) jumps from positive to negative value while the MS diffusivity between cation-cation pair (Đ(LiBe)) becomes positive from a negative value. Even though the negative MS diffusivities have been observed for electrolyte solutions between cation-cation pair, here we report negative MS diffusivity between cation-anion pair where Đ(BeF) shows a sharp rise around 66% of BeF2, reaches maximum value at 70% of BeF2, and then decreases almost exponentially with a sign change for BeF2 around 93%. For low mole fraction of LiF, Đ(BeF) follows the Debye-Huckel theory and rises with the square root of LiF mole fraction similar to the MS diffusivity between cation-anion pair in aqueous solution of electrolyte salt. Negative MS diffusivities while unusual are, however, shown to satisfy the non-negative entropy constraints at all thermodynamic states as required by the second law of thermodynamics. We have established a strong correlation between the structure and dynamics and predict that the formation of flouride polyanion network between Be and F ions and coulomb interaction is responsible for sharp variation of the MS diffusivities which controls the multicomponent diffusion phenomenon in LiF-BeF2 which has a strong impact on the performance of the reactor.