Tysonite structure fluorides doped with divalent cations, represented by Ce 0.95 Ca 0.05 F 2.95 , are a class of good F − ion conductors together with fluoritestructured compounds. Computational understanding of the F − conduction process is difficult because of the complicated interactions between three symmetrically distinct F sites and the experimentally observed change in the F diffusion mechanism slightly above room temperature, effectively making first-principles molecular dynamics (FP-MD) simulations, which are often conducted well above the transition temperature, useless when analyzing behavior below the transition point. Neural network potential (NNP) MD simulations showed that the F diffusion coefficient is higher when the divalent dopant cation size is similar to the trivalent cation size. The diffusion behavior of F in different sites changes at roughly 500 K in Ce 0.95 Ca 0.05 F 2.95 because only the F1 site sublattice contributes to F diffusion below this temperature, but the remaining F2 and F3 sublattices become gradually active above this temperature. The paradox of higher diffusion coefficients in CeF 3 -based compounds than similar LaF 3 -based compounds even though the lattice parameters are larger in the latter may be caused by a shallower potential of Ce and F in CeF 3 compared to the LaF 3 counterparts.