Abstract. -A fitting scheme is proposed to obtain effective potentials from Car-Parrinello molecular dynamics (CPMD) simulations. It is used to parameterize a new pair potential for silica. MD simulations with this new potential are done to determine structural and dynamic properties and to compare these properties to those obtained from CPMD and a MD simulation using the so-called BKS potential. The new potential reproduces accurately the liquid structure generated by the CPMD trajectories, the experimental activation energies for the self-diffusion constants and the experimental density of amorphous silica. Also lattice parameters and elastic constants of α−quartz are well-reproduced, showing the transferability of the new potential.

We show that finite-range alternatives to the standard long-range BKS pair potential for silica might be used in molecular dynamics simulations. We study two such models that can be efficiently simulated since no Ewald summation is required. We first consider the Wolf method, where the Coulomb interactions are truncated at a cutoff distance rc such that the requirement of charge neutrality holds. Various static and dynamic quantities are computed and compared to results from simulations using Ewald summations. We find very good agreement for rc ≈ 10Å. For lower values of rc, the long-range structure is affected which is accompanied by a slight acceleration of dynamic properties. In a second approach, the Coulomb interaction is replaced by an effective Yukawa interaction with two new parameters determined by a force fitting procedure. The same trend as for the Wolf method is seen. However, slightly larger cutoffs have to be used in order to obtain the same accuracy with respect to static and dynamic quantities as for the Wolf method.