Trapping of Li(+) Ions by [ThFn](4-n) Clusters Leading to Oscillating Maxwell-Stefan Diffusivity in the Molten Salt LiF-ThF4

J Phys Chem B. 2016 Aug 18;120(32):7926-36. doi: 10.1021/acs.jpcb.6b05173. Epub 2016 Aug 5.

Abstract

A molten salt mixture of lithium fluoride and thorium fluoride (LiF-ThF4) serves as a fuel as well as a coolant in the most sophisticated molten salt reactor (MSR). Here, we report for the first time dynamic correlations, Onsager coefficients, Maxwell-Stefan (MS) diffusivities, and the concentration dependence of density and enthalpy of the molten salt mixture LiF-ThF4 at 1200 K in the composition range of 2-45% ThF4 and also at eutectic composition in the temperature range of 1123-1600 K using Green-Kubo formalism and equilibrium molecular dynamics simulations. We have observed an interesting oscillating pattern for the MS diffusivity for the cation-cation pair, in which ĐLi-Th oscillates between positive and negative values with the amplitude of the oscillation reducing as the system becomes rich in ThF4. Through the velocity autocorrelation function, vibrational density of states, radial distribution function analysis, and structural snapshots, we establish an interplay between the local structure and multicomponent dynamics and predict that formation of negatively charged [ThFn](4-n) clusters at a higher ThF4 mole % makes positively charged Li(+) ions oscillate between different clusters, with their range of motion reducing as the number of [ThFn](4-n) clusters increases, and finally Li(+) ions almost get trapped at a higher ThF4% when the electrostatic force on Li(+) exerted by various surrounding clusters gets balanced. Although reports on variations of density and enthalpy with temperature exist in the literature, for the first time we report variations of the density and enthalpy of LiF-ThF4 with the concentration of ThF4 (mole %) and fit them with the square root function of ThF4 concentration, which will be very useful for experimentalists to obtain data over a range of concentrations from fitting the formula for design purposes. The formation of [ThFn](4-n) clusters and the reduction in the diffusivity of the ions at a higher ThF4% may limit the percentage of ThF4 that can be used in the MSR to optimize the neutron economy.