Analysis of the simulation data was performed using the RevelsMD , site-analysis [119], ase [120], pymatgen , numpy , and scipy codes. The time-average fluorine density (Fig. ) was calculated using a linear combination of a conventional histogram with a triangular kernel and a force-extrapolated analogue, as described in Refs. .