The structural, thermochemical, and thermophysical properties of the AnCl₄, and NaCl-AnCl₄ (An = Th, U) melts were investigated using molecular dynamics (MD) simulations based on the Polarisable Ion Model (PIM). New force fields were proposed and used to compute key properties including density, thermal expansion, enthalpy of mixing, heat capacity, as well as the local structure and chemical speciation in the molten (Na, An)Cl_x (An = Th, U) salts. Thermodynamic models were then developed based on the CALPHAD method, using both PIM-MD and experimental data as input. Employing the modified quasichemical formalism in the quadruplet approximation for the liquid solution, the models account for the chemical speciation in the melt as calculated by MD simulations, and reproduce well phase equilibria in those systems. In particular, the models included monomeric and dimeric species to represent the physical nature of the ionic melt, which shows progressive oligomerisation with increasing AnCl₄ fraction. Our studies confirm that the melt becomes highly volatile at high AnCl₄ fractions, which is discussed in light of the results obtained herein.

The short-range structures of LiF-ThF4, NaF-AnF4, KF-AnF4, and Cs-AnF4 (An = Th, U), were probed using in situ high temperature Extended X-ray Absorption Fine Structure (EXAFS) spectroscopy. Signally, the EXAFS spectra of pure molten ThF4 and UF4 were measured for the first time. The data were interpreted with the aid of Molecular Dynamics (MD) and standard fitting of the EXAFS equation. As in related studies, a speciation distribution dominated by [AnFx]4-x (x = 7, 8, 9) coordination complexes was observed. The average coordination number was found to decrease with the increasing size of the alkali cation, and increase with AnF4 content. An average coordination number close to 6, which had not been detected before in melts of alkali actinide fluorides, was seen when CsF was used as solvent. This journal is

The structural, thermochemical, and thermophysical properties of the NaF-ThF4 fuel system were studied with experimental methods and molecular dynamics (MD) simulations. Equilibrium MD (EMD) simulations using the polarizable ion model were performed to calculate the density, molar volume, thermal expansion, mixing enthalpy, heat capacity, and distribution of [ThFn]m- complexes in the (Na,Th)Fx melt over the full concentration range at various temperatures. The phase equilibria in the 10-50 mol % ThF4 and 85-95 mol % ThF4 regions of the NaF-ThF4 phase diagram were measured using differential scanning calorimetry, as were the mixing enthalpies at 1266 K of (NaF/ThF4) = (0.8:0.2), (0.7:0.3) mixtures. Furthermore, the β-Na2ThF6 and NaTh2F9 compounds were synthesized and subsequently analyzed with the use of X-ray diffraction. The heat capacities of both compounds were measured in the temperature ranges (2-271 K) and (2-294 K), respectively, by thermal relaxation calorimetry. Finally, a CALPHAD model coupling the structural and thermodynamic data was developed using both EMD and experimental data as input and a quasichemical formalism in the quadruplet approximation. Here, 7- and 8-coordinated Th4+ cations were introduced on the cationic sublattice alongside a 13-coordinated dimeric species to reproduce the chemical speciation, as calculated by EMD simulations and to provide a physical description of the melt.

The (Li,Be)F_x fluoride salt is an ionic liquid with complex non-ideal thermodynamic behaviour due to the formation of short-range order. In this work, we explore the relationship between local structure, thermo-physical and thermodynamic properties in this system using a multidisciplinary approach that couples molecular dynamics simulations using the Polarizable Ion Model (PIM) and thermodynamic modelling assessment using the CALPHAD method. The density, thermal expansion, viscosity, thermal conductivity, molar and mixing enthalpies and heat capacity of the (Li,Be)F_x melt are extracted from the polarizable ionic interaction potentials and investigated across a wide range of compositions and temperatures. The agreement with the available experimental data is generally very good. The local structure is also examined in detail, in particular the transition between a molecular liquid with Li⁺, BeF₄ ²⁻ and F⁻ predominant species at low BeF₂ content, and a polymeric liquid at high BeF₂ content, with the formation of polymers (Be₂F₇ ³⁻, Be₃F₁₀ ⁴⁻, Be₄F₁₃ ⁵⁻, etc.), and finally of a three-dimensional network of corner-sharing tetrahedrally coordinated Be²⁺ cations for pure BeF₂. Based on the available experimental information and the output of the MD simulations, we moreover develop for the first time a coupled structural-thermodynamic model for the LiF-BeF₂ system based on the quasi-chemical formalism in the quadruplet approximation, that provides a physical description of the melt and reproduces (in addition to the thermodynamic data) the chemical speciation of beryllium polymeric species predicted from the simulations.

The optimized Gibbs energies for the second-nearest neighbours exchange reactions of the liquid solution were not reported correctly in Eqs. (17)–(19) in our recent work [1]. The corresponding equations should read: [Formula presented] [Formula presented] [Formula presented]