Free energy profiles of Al3+ and La3+ cation distribution in silica and soda silicate glasses

The factors that control the distribution of Al3+ and La3+ cations in silica and soda silicate glasses is examined by using molecular dynamics (MD) simulations. In particular, the response of the glass network to the presence of metal oxide is probed using a liquid state theory that treats the glass network as a solvent and the metal cation as a solute. MD simulations are used to obtain the mean solvent–solute and solute–solute force. The trajectory used to determine the free energy is analyzed to determine the stable configurations of a cation pair. Details of determining the potential of mean force for an Al cation pair in silica and silicate glass is presented. A comparison of these results with those previously calculated for a La cation pair in the same glass systems is given. The results reveal that there are differences in how the network accommodates the two different size cations. It is found that for the potential used here, based on the Vessal potential, the network wraps itself around the larger La cation forming a solvent shell, whereas, the smaller Al cation is incorporated into the network backbone. In silica and soda silicate glasses, La ion pairs cluster to form La–O–La linkages. In contrast, the glasses favor a separated state of the Al ion pair.