Molecular dynamics simulation for ternary lithium-aluminium-phosphate glass

P2O5-Al2O3-Li2O vitreous system structure was simulated by molecular dynamic (MD) method. The knowledge of analytical formula of the interaction forces enables, for a limited number of ions, to solve the Newtonian equations, which describe the particles movement. The simulation gives atomic system configuration images, at several temperatures, starting at 3000K, down to 1000K, 700K and 300K temperatures. Information about glass structure was determined by calculation of radial distribution function. The maxima from the radial distribution function represent the most likely distances between two ionic species.