Molecular Dynamics Simulation of the Structure and Hydroxylation of Silica Glass Surfaces

Abstract

The surface structure of silica glasses has been simulated using molecular dynamics. The surface hydroxyl concentration was estimated to be 4.5/nm 2 , based on surface defect statistics. Hydroxyl-silica potentials were developed and used to study the hydroxylation of silica surface. It is found that the energy of chemisorption of water declines in the sequence: three coordinated silicon (Si 3 ) and non-bridging oxygen (NBO) on separate sites, Si 3 and NBO on combined sites, two-and three-membered rings. Partial hydroxylation of the most reactive sites, which leads to an OH coverage of 2.5/nm 2 , was studied. Structural relaxation after hydroxylation was observed. J ournal

Silica Surface Features and Their Role in the Adsorption of Biomolecules: Computational Modeling and Experiments

“…In other words, the state of the fully hydroxylated silica surfaces, either crystalline or amorphous, the latter being the most common for real silica surfaces (see Section 2) has hardly been addressed with these force field (as exceptions see Ref. [462][463][464][465][466] ), as these studies were not concerned with the description of the silanol groups and their mutual H-bonding interactions. To specifically address the H-bonded features at the silica surface the FFSiOH force field has recently been proposed, demonstrating excellent performances for both the bulk and the surface of silica materials as regards structures, energetic and phonon spectra, including the OH stretching region.…”

“…The need for a much more extended set of parameters to deal with the interactions between the biomolecule and the surface has prevented this kind of studies from appearing in the literature. For the silica/water system (see Section 8 for a deeper discussion), Garofalini et al 464,465 Singer et al 462, 463 and Du et al 466 468 a reactive force field which was designed to model equilibrium structures as well as reaction intermediates and activated complexes in chemical reactions. The idea is to correlate bond lengths with bond orders and the latter with bond energies, following old insights from Linus Pauling's work.…”

UV-induced modification of fused silica: Insights from ReaxFF-based molecular dynamics simulations

“…two SiO 4 tetrahedrons share edges to form the smallest closed paths, which are apt to generate in silica surfaces through atom reorganization after cleavage. 22,23 Unlike those located in surfaces, wherein silicon remains regularly coordinated, a 2-membered ring in the irradiated bulk usually contains a Si 5 defect. This unique observation indicates that the formation of such 2-membered rings is caused by laser-induced Si 5 (not the native ones), which divided the network rings into smaller ones, resulting in small Si-O-Si bond angles and thus decreased Si-Si pair distances.…”

Temperature-dependent structural heterogeneity in calcium silicate liquids

“…This method has been successfully applied to a number of alkali silicate, rare-earth aluminate, and silicate simulations. 19,[22][23][24][25][26][27] The partial atomic charges and Buckingham potential parameters are shown in Table I.…”