Reactive simulations of the activation barrier to dissolution of amorphous silica in water

Michael Kagan, Glenn K. Lockwood, Steve Garofalini

Research output: Contribution to journalArticle

9 Citations (Scopus)

Abstract

Molecular dynamics simulations employing reactive potentials were used to determine the activation barriers to the dissolution of the amorphous SiO 2 surface in the presence of a 2 nm overlayer of water. The potential of mean force calculations of the reactions of water molecules with 15 different starting Q4 sites (Qi is the Si site with i bridging oxygen neighbors) to eventually form the dissolved Q0 site were used to obtain the barriers. Activation barriers for each step in the dissolution process, from the Q4 to Q3 to Q2 to Q1 to Q0 were obtained. Relaxation runs between each reaction step enabled redistribution of the water above the surface in response to the new Qi site configuration. The rate-limiting step observed in the simulations was in both the Q32 reaction (a Q3 site changing to a Q2 site) and the Q21 reaction, each with an average barrier of ∼14.1 kcal mol-1. However, the barrier for the overall reaction from the Q4 site to a Q0 site, averaged over the maximum barrier for each of the 15 samples, was 15.1 kcal mol-1. This result is within the lower end of the experimental data, which varies from 14-24 kcal mol-1, while ab initio calculations using small cluster models obtain values that vary from 18-39 kcal mol-1. Constraints between the oxygen bridges from the Si site and the connecting silica structure, the presence of pre-reaction strained siloxane bonds, and the location of the reacting Si site within slight concave surface contours all affected the overall activation barriers.

Original languageEnglish
Pages (from-to)9294-9301
Number of pages8
JournalPhysical Chemistry Chemical Physics
Volume16
Issue number20
DOIs
Publication statusPublished - May 28 2014

Fingerprint

Silicon Dioxide
dissolving
Dissolution
Chemical activation
activation
silicon dioxide
Water
water
Oxygen
Siloxanes
simulation
Molecular dynamics
Molecules
Computer simulation
siloxanes
oxygen
molecular dynamics
configurations
molecules

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Physics and Astronomy(all)

Cite this

Reactive simulations of the activation barrier to dissolution of amorphous silica in water. / Kagan, Michael; Lockwood, Glenn K.; Garofalini, Steve.

In: Physical Chemistry Chemical Physics, Vol. 16, No. 20, 28.05.2014, p. 9294-9301.

Research output: Contribution to journalArticle

@article{a391a5151c7542a79df20acb77e4d2ee,
title = "Reactive simulations of the activation barrier to dissolution of amorphous silica in water",
abstract = "Molecular dynamics simulations employing reactive potentials were used to determine the activation barriers to the dissolution of the amorphous SiO 2 surface in the presence of a 2 nm overlayer of water. The potential of mean force calculations of the reactions of water molecules with 15 different starting Q4 sites (Qi is the Si site with i bridging oxygen neighbors) to eventually form the dissolved Q0 site were used to obtain the barriers. Activation barriers for each step in the dissolution process, from the Q4 to Q3 to Q2 to Q1 to Q0 were obtained. Relaxation runs between each reaction step enabled redistribution of the water above the surface in response to the new Qi site configuration. The rate-limiting step observed in the simulations was in both the Q32 reaction (a Q3 site changing to a Q2 site) and the Q21 reaction, each with an average barrier of ∼14.1 kcal mol-1. However, the barrier for the overall reaction from the Q4 site to a Q0 site, averaged over the maximum barrier for each of the 15 samples, was 15.1 kcal mol-1. This result is within the lower end of the experimental data, which varies from 14-24 kcal mol-1, while ab initio calculations using small cluster models obtain values that vary from 18-39 kcal mol-1. Constraints between the oxygen bridges from the Si site and the connecting silica structure, the presence of pre-reaction strained siloxane bonds, and the location of the reacting Si site within slight concave surface contours all affected the overall activation barriers.",
author = "Michael Kagan and Lockwood, {Glenn K.} and Steve Garofalini",
year = "2014",
month = "5",
day = "28",
doi = "10.1039/c4cp00030g",
language = "English",
volume = "16",
pages = "9294--9301",
journal = "Physical Chemistry Chemical Physics",
issn = "1463-9076",
publisher = "Royal Society of Chemistry",
number = "20",

}

TY - JOUR

T1 - Reactive simulations of the activation barrier to dissolution of amorphous silica in water

AU - Kagan, Michael

AU - Lockwood, Glenn K.

AU - Garofalini, Steve

PY - 2014/5/28

Y1 - 2014/5/28

N2 - Molecular dynamics simulations employing reactive potentials were used to determine the activation barriers to the dissolution of the amorphous SiO 2 surface in the presence of a 2 nm overlayer of water. The potential of mean force calculations of the reactions of water molecules with 15 different starting Q4 sites (Qi is the Si site with i bridging oxygen neighbors) to eventually form the dissolved Q0 site were used to obtain the barriers. Activation barriers for each step in the dissolution process, from the Q4 to Q3 to Q2 to Q1 to Q0 were obtained. Relaxation runs between each reaction step enabled redistribution of the water above the surface in response to the new Qi site configuration. The rate-limiting step observed in the simulations was in both the Q32 reaction (a Q3 site changing to a Q2 site) and the Q21 reaction, each with an average barrier of ∼14.1 kcal mol-1. However, the barrier for the overall reaction from the Q4 site to a Q0 site, averaged over the maximum barrier for each of the 15 samples, was 15.1 kcal mol-1. This result is within the lower end of the experimental data, which varies from 14-24 kcal mol-1, while ab initio calculations using small cluster models obtain values that vary from 18-39 kcal mol-1. Constraints between the oxygen bridges from the Si site and the connecting silica structure, the presence of pre-reaction strained siloxane bonds, and the location of the reacting Si site within slight concave surface contours all affected the overall activation barriers.

AB - Molecular dynamics simulations employing reactive potentials were used to determine the activation barriers to the dissolution of the amorphous SiO 2 surface in the presence of a 2 nm overlayer of water. The potential of mean force calculations of the reactions of water molecules with 15 different starting Q4 sites (Qi is the Si site with i bridging oxygen neighbors) to eventually form the dissolved Q0 site were used to obtain the barriers. Activation barriers for each step in the dissolution process, from the Q4 to Q3 to Q2 to Q1 to Q0 were obtained. Relaxation runs between each reaction step enabled redistribution of the water above the surface in response to the new Qi site configuration. The rate-limiting step observed in the simulations was in both the Q32 reaction (a Q3 site changing to a Q2 site) and the Q21 reaction, each with an average barrier of ∼14.1 kcal mol-1. However, the barrier for the overall reaction from the Q4 site to a Q0 site, averaged over the maximum barrier for each of the 15 samples, was 15.1 kcal mol-1. This result is within the lower end of the experimental data, which varies from 14-24 kcal mol-1, while ab initio calculations using small cluster models obtain values that vary from 18-39 kcal mol-1. Constraints between the oxygen bridges from the Si site and the connecting silica structure, the presence of pre-reaction strained siloxane bonds, and the location of the reacting Si site within slight concave surface contours all affected the overall activation barriers.

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

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

U2 - 10.1039/c4cp00030g

DO - 10.1039/c4cp00030g

M3 - Article

VL - 16

SP - 9294

EP - 9301

JO - Physical Chemistry Chemical Physics

JF - Physical Chemistry Chemical Physics

SN - 1463-9076

IS - 20

ER -