Computational study of uniaxial deformations in silica aerogel using a coarse-grained model

Ferreiro-Rangel, Carlos A. and Gelb, Lev D. (2015) Computational study of uniaxial deformations in silica aerogel using a coarse-grained model. Journal of Physical Chemistry B, 119 (27). pp. 8640-8650. ISSN 1520-6106 (https://doi.org/10.1021/jp512998w)

[thumbnail of Ferreiro-Rangel-Gelb-JOPCB-2015-Computational-study-of-uniaxial-deformations-in-silica-aerogel]
Preview
Text. Filename: Ferreiro_Rangel_Gelb_JOPCB_2015_Computational_study_of_uniaxial_deformations_in_silica_aerogel.pdf
Accepted Author Manuscript

Download (1MB)| Preview

Abstract

Simulations of a flexible coarse-grained model are used to study silica aerogels. This model, introduced in a previous study (J. Phys. Chem. C 2007, 111, 15792), consists of spherical particles which interact through weak nonbonded forces and strong interparticle bonds that may form and break during the simulations. Small-deformation simulations are used to determine the elastic moduli of a wide range of material models, and large-deformation simulations are used to probe structural evolution and plastic deformation. Uniaxial deformation at constant transverse pressure is simulated using two methods: a hybrid Monte Carlo approach combining molecular dynamics for the motion of individual particles and stochastic moves for transverse stress equilibration, and isothermal molecular dynamics simulations at fixed Poisson ratio. Reasonable agreement on elastic moduli is obtained except at very low densities. The model aerogels exhibit Poisson ratios between 0.17 and 0.24, with higher-density gels clustered around 0.20, and Young's moduli that vary with aerogel density according to a power-law dependence with an exponent near 3.0. These results are in agreement with reported experimental values. The models are shown to satisfy the expected homogeneous isotropic linear-elastic relationship between bulk and Young's moduli at higher densities, but there are systematic deviations at the lowest densities. Simulations of large compressive and tensile strains indicate that these materials display a ductile-to-brittle transition as the density is increased, and that the tensile strength varies with density according to a power law, with an exponent in reasonable agreement with experiment. Auxetic behavior is observed at large tensile strains in some models. Finally, at maximum tensile stress very few broken bonds are found in the materials, in accord with the theory that only a small fraction of the material structure is actually load-bearing.