skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Crack propagation in silica from reactive classical molecular dynamics simulations

ORCiD logo [1]; ORCiD logo [2]; ORCiD logo [1]
  1. Department of Geochemistry, Sandia National Laboratories, Albuquerque NM
  2. Department of Mechanics of Materials, Sandia National Laboratories, Livermore CA
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of the American Ceramic Society
Additional Journal Information:
Journal Volume: 101; Journal Issue: 4; Related Information: CHORUS Timestamp: 2018-02-02 03:01:06; Journal ID: ISSN 0002-7820
Country of Publication:
United States

Citation Formats

Rimsza, Jessica M., Jones, Reese E., and Criscenti, Louise J.. Crack propagation in silica from reactive classical molecular dynamics simulations. United States: N. p., 2017. Web. doi:10.1111/jace.15292.
Rimsza, Jessica M., Jones, Reese E., & Criscenti, Louise J.. Crack propagation in silica from reactive classical molecular dynamics simulations. United States. doi:10.1111/jace.15292.
Rimsza, Jessica M., Jones, Reese E., and Criscenti, Louise J.. 2017. "Crack propagation in silica from reactive classical molecular dynamics simulations". United States. doi:10.1111/jace.15292.
title = {Crack propagation in silica from reactive classical molecular dynamics simulations},
author = {Rimsza, Jessica M. and Jones, Reese E. and Criscenti, Louise J.},
abstractNote = {},
doi = {10.1111/jace.15292},
journal = {Journal of the American Ceramic Society},
number = 4,
volume = 101,
place = {United States},
year = 2017,
month =

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on November 8, 2018
Publisher's Accepted Manuscript

Save / Share:
  • We have used molecular dynamics simulations based on a three-body potential with charge transfer to generate nanoporous silica aerogels. Care was taken to reproduce the sol-gel condensation reaction that forms the gel backbone as realistically as possible and to thereby produce credible gel structures. The self-similarity of aerogel structures was investigated by evaluating their fractal dimension from geometric correlations. For comparison, we have also generated porous silica glasses by rupturing dense silica and computed their fractal dimension. The fractal dimension of the porous silica structures was found to be process dependent. Finally, we have determined that the effect of supercriticalmore » drying on the fractal nature of condensed silica gels is not appreciable.« less
  • Molecular dynamics simulations of heat transfer in gases are computationally expensive when the wall molecules are explicitly modeled. To save computational time, an implicit boundary function is often used. Steele's potential has been used in studies of fluid–solid interface for a long time. In this work, the conceptual idea of Steele's potential was extended in order to simulate water–silicon and water–silica interfaces. A new wall potential model is developed by using the electronegativity-equalization method (EEM), a ReaxFF empirical force field and a non-reactive molecular dynamics package PumMa. Contact angle simulations were performed in order to validate the wall potential model.more » Contact angle simulations with the resulting tabulated wall potentials gave a silicon–water contact angle of 129°, a quartz–water contact angle of 0°, and a cristobalite–water contact angle of 40°, which are in reasonable agreement with experimental values.« less
  • We have investigated the transition from the atomistic to the macroscopic impact mechanism by simulating large Argon cluster impacts on amorphous silica. The transition occurs at cluster sizes less than 50 000 atoms at hypervelocity regime (22 km/s). After that, the crater volume increases linearly with the cluster size opposite to the nonlinear scaling typical of small cluster impacts. The simulations demonstrate that the molecular dynamics method can be used to explore atomistic mechanisms that lead to damage formation in small particle impacts, for example, in impacts of micrometeorites on spacecraft.
  • Ab initio molecular dynamics (AIMD) simulations of the hydration shells surrounding the Zn2+ ion are reported for temperatures near 300oC. Simulations using a combined ab initio and classical molecular dynamics (AIMD/MM) approach are also carried out. Both simulations are done with 64 solvating water molecules (~15 ps). The hydration structure predicted from both simulations is found to agree very well with known results from X-ray data. The 1st hydration shell contains six water molecules in an octahedral structure with the hydrogen atoms oriented away from the Zn2+ ion. The six waters in the 1st shell are located at an averagemore » distance of 2.44Å. A 2nd hydration shell is observed at 4.59Å. Beyond these shells, the bonding pattern substantially returns to the tetrahedral structure of bulk water. No exchanges are seen between the 1st and 2nd hydrations shells, however many water transfers between the 2nd and outer hydrations shells are observed to occur on a picosecond (ps) time scale via dissociative and associative mechanisms. In general, it is found that the AIMD and AIMD/MM simulations give nearly identical results for structural parameters, EXAFS spectra, and exchange dynamics. These results suggest that AIMD/MM can be used to extend the particle scale and time scale of AIMD simulations of highly charged ions in solution.« less