RMC modelling

The model consists in 1950 atoms in a cubic box (30.12 Å length for corresponding to the number density Rho = 0.07135 as determined from the glass density). The initial positions were generated from a random filling of the box by the M atoms first, then the Na, the Pb and the F atoms successively were inserted. Positions at this filling stage were accepted if minimal predefined interatomic distances were respected. Then, during the RMC simulation runs, distances (Å) of closest approach were applied too : Na-Na = 3.15, Na-Pb = 3.25, Na-Fe = 3.05, Na-F = 2.05, Pb-Pb = 3.35, Pb-Fe = 3.25, Pb-F = 2.15, Fe-Fe = 3.05, Fe-F = 1.73 and F-F = 2.28. The MF6 coordination was constrained to occur with a maximum M-F distance equal to 2.15 Å. The neutron diffraction data of NaPbM2F9 (M = Fe, V) glasses and the X-ray data for M = Fe were simulated as F(Q) data, the total coherent scattering functions (TSF) :

F(Q) = [Icoh(Q)-<f2>] / <f>2,

where the <f2> and <f>2 are the usual mean diffusion factors, depending on Q (X-ray) or not (neutron). About 16 million accepted moves were completed before the sixfold constraint was fulfilled, the last of the 300 octahedra was particularly long to be obtained. The whole process took 720 hours CPU on a DEC-ALPHA AXP 4620.

The R factors calculated as 100*|Iobs-kIcalc|/|Iobs| (%) (according to the definition I(Q) = F(Q)+1, k being a scale factor) are 1.42, 1.29 and 1.72 % respectively for the Fe and V-based neutron data and the Fe-based X-ray data. The observed and calculated F(Q) curves are shown on the figure 1. A global view of the octahedra network is shown on the figure 2.

Fig. 1. Experimental (dots) and RMC simulated (solid line) TSF of NaPbM2F9; (a) neutron data for M = Fe, (b) neutron data for M = V, (c) X-ray data for M ) Fe.

Fig. 2. Projection of the RMC cubic box showing the 300 [MF6] polyhedra network.

Armel Le Bail - June 1997