A force field for the modelling of neodymium fluorobritholites, having the general composition Ca 1 0 - y Nd y (SiO 4 ) y (PO 4 ) 6 - y F 2 , is derived. The validity of the established analytical potentials is tested by calculating the crystallographic parameters of the compositions for 0=<y=<6 and comparing to the experimental values. The validity is further tested by calculating the variation of the phonon spectra when transitioning from the fluoroapatite Ca 1 0 (PO 4 ) 6 F 2 to the full silicate neodymium fluorobritholite Ca 4 Nd 6 (SiO 4 ) 6 F 2 . Trivalent plutonium ions are introduced into the mono-silicate neodymium fluorobritholite using transferable potentials in order to obtain the composition Ca 9 Pu 0 . 5 Nd 0 . 5 (SiO 4 )(PO 4 ) 5 F 2 , which corresponds to a Pu concentration of about 10 wt%. The thermodynamic properties of the resulting structure are calculated, and the threshold displacement energies of various sublattices are estimated by applying the sudden approximation (SA) method. The activation energies for Frenkel defect annealing are also calculated using the rational function optimisation (RFO) procedure for the saddle point research.