Eur. Phys. J. E (2013) 36: 60
https://doi.org/10.1140/epje/i2013-13060-9

Regular Article

Computer simulation of diffusion in silica liquid under temperature and pressure

Department of Computational Physics, Hanoi University of Science and Technology, 1 Dai Co Viet, Hanoi, Viet Nam

^* e-mail: pkhung@fpt.vn

Received: 2 October 2012
Accepted: 21 May 2013
Published online: 20 June 2013

Abstract

We have studied the diffusion mechanism in silica liquid following a new approach where the diffusion rate is estimated via the rate of SiO _x SiO _x±1 and the mean square displacement of Si particles per SiO _x SiO _x±1 . Molecular dynamics simulation has been conducted for a model consisting of 1998 particles over a wide range of temperatures (3000-4500K) and pressure (from 0 to 25.75GPa). Our results show that the rate of SiO _x SiO _x±1 increases either with increasing the temperature or pressure. Further, we find that SiO _x SiO _x±1 is heterogeneously distributed through the network structure of the liquid. In particular, it is concentrated on a small section of Si particles in a low-temperature regime and at ambient pressure. The spatial localisation of SiO _x SiO _x±1 originates from the fact that the stable unit in low- and high-pressure regime is SiO₄ and SiO₆ , respectively. The major change in the diffusion mechanism under pressure or temperature concerns the change in the distribution of SiO _x SiO _x±1 through the network structure. It is finally shown that the spatial localisation of SiO _x SiO _x±1 is responsible for the dynamics heterogeneity and the diffusion anomaly for silica liquid. This finding supports the concept that as the temperature approaches the glass transition point, SiO _x SiO _x±1 spatially localises such that the diffusivity drops and the dynamics are anomalously slow.