Empirical three-body potential for vitreous silica

B. P. Feuston, Steve Garofalini

Research output: Contribution to journalArticle

358 Citations (Scopus)

Abstract

A three-body potential suitable for molecular dynamics (MD) simulations has been developed for vitreous silica by adding three-body interactions to the Born-Mayer-Huggins (BMH) pair potential. Previous MD simulations with the BMH potential have formed glassy SiO2 through the melt-quench method with some success. Though bond lengths were found to be in fair agreement with experiment, the distribution of tetrahedral angles was too broad and the model glass contained 6%-8% bond defects. This is indicative of a lack of the local order that is present in the laboratory glass. The nature of the short range order is expected to play an important role when investigating defect formation, surface reconstruction, or surface reactivities. An attempt has been made to increase the local order in the simulated glass by including a directional dependent term in the effective potential to model the partial covalency of the Si-O bond. The vitreous state obtained through MD simulation with this modified BMH potential shows an increase in the short range order with a narrow O-Si-O angle distribution peaked about the tetrahedral angle and a low concentration of bond defects, typically ∼1%-2%. The static structure factor S(q) is calculated and found to be in good agreement with neutron scattering results. Intermediate range order is also discussed in reference to the distribution of ring sizes.

Original languageEnglish
Pages (from-to)5818-5824
Number of pages7
JournalJournal of Chemical Physics
Volume89
Issue number9
Publication statusPublished - 1988

Fingerprint

Fused silica
Molecular dynamics
silicon dioxide
Glass
Defects
Computer simulation
molecular dynamics
Surface reconstruction
Bond length
glass
Neutron scattering
defects
simulation
low concentrations
neutron scattering
reactivity
Experiments
rings
interactions

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Empirical three-body potential for vitreous silica. / Feuston, B. P.; Garofalini, Steve.

In: Journal of Chemical Physics, Vol. 89, No. 9, 1988, p. 5818-5824.

Research output: Contribution to journalArticle

@article{f9d12d659cb440648276f67af5c6832b,
title = "Empirical three-body potential for vitreous silica",
abstract = "A three-body potential suitable for molecular dynamics (MD) simulations has been developed for vitreous silica by adding three-body interactions to the Born-Mayer-Huggins (BMH) pair potential. Previous MD simulations with the BMH potential have formed glassy SiO2 through the melt-quench method with some success. Though bond lengths were found to be in fair agreement with experiment, the distribution of tetrahedral angles was too broad and the model glass contained 6{\%}-8{\%} bond defects. This is indicative of a lack of the local order that is present in the laboratory glass. The nature of the short range order is expected to play an important role when investigating defect formation, surface reconstruction, or surface reactivities. An attempt has been made to increase the local order in the simulated glass by including a directional dependent term in the effective potential to model the partial covalency of the Si-O bond. The vitreous state obtained through MD simulation with this modified BMH potential shows an increase in the short range order with a narrow O-Si-O angle distribution peaked about the tetrahedral angle and a low concentration of bond defects, typically ∼1{\%}-2{\%}. The static structure factor S(q) is calculated and found to be in good agreement with neutron scattering results. Intermediate range order is also discussed in reference to the distribution of ring sizes.",
author = "Feuston, {B. P.} and Steve Garofalini",
year = "1988",
language = "English",
volume = "89",
pages = "5818--5824",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "9",

}

TY - JOUR

T1 - Empirical three-body potential for vitreous silica

AU - Feuston, B. P.

AU - Garofalini, Steve

PY - 1988

Y1 - 1988

N2 - A three-body potential suitable for molecular dynamics (MD) simulations has been developed for vitreous silica by adding three-body interactions to the Born-Mayer-Huggins (BMH) pair potential. Previous MD simulations with the BMH potential have formed glassy SiO2 through the melt-quench method with some success. Though bond lengths were found to be in fair agreement with experiment, the distribution of tetrahedral angles was too broad and the model glass contained 6%-8% bond defects. This is indicative of a lack of the local order that is present in the laboratory glass. The nature of the short range order is expected to play an important role when investigating defect formation, surface reconstruction, or surface reactivities. An attempt has been made to increase the local order in the simulated glass by including a directional dependent term in the effective potential to model the partial covalency of the Si-O bond. The vitreous state obtained through MD simulation with this modified BMH potential shows an increase in the short range order with a narrow O-Si-O angle distribution peaked about the tetrahedral angle and a low concentration of bond defects, typically ∼1%-2%. The static structure factor S(q) is calculated and found to be in good agreement with neutron scattering results. Intermediate range order is also discussed in reference to the distribution of ring sizes.

AB - A three-body potential suitable for molecular dynamics (MD) simulations has been developed for vitreous silica by adding three-body interactions to the Born-Mayer-Huggins (BMH) pair potential. Previous MD simulations with the BMH potential have formed glassy SiO2 through the melt-quench method with some success. Though bond lengths were found to be in fair agreement with experiment, the distribution of tetrahedral angles was too broad and the model glass contained 6%-8% bond defects. This is indicative of a lack of the local order that is present in the laboratory glass. The nature of the short range order is expected to play an important role when investigating defect formation, surface reconstruction, or surface reactivities. An attempt has been made to increase the local order in the simulated glass by including a directional dependent term in the effective potential to model the partial covalency of the Si-O bond. The vitreous state obtained through MD simulation with this modified BMH potential shows an increase in the short range order with a narrow O-Si-O angle distribution peaked about the tetrahedral angle and a low concentration of bond defects, typically ∼1%-2%. The static structure factor S(q) is calculated and found to be in good agreement with neutron scattering results. Intermediate range order is also discussed in reference to the distribution of ring sizes.

UR - http://www.scopus.com/inward/record.url?scp=33646268111&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=33646268111&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:33646268111

VL - 89

SP - 5818

EP - 5824

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 9

ER -