The structure and dynamics of silica polymorphs using a two‐body effective potential model

Anomalous surface states modify the size-dependent mechanical properties and fracture of silica nanowires

Molecular dynamics simulations of amorphous silica nanowires under tension were analyzed for size and surface stress effects on mechanical properties and for structural modifications via bond angle distributions. Their fracture behavior was also investigated beyond the elastic limit. The Young's moduli of silica nanowires were predicted to be about 75-100 GPa, depending on the nanowire size. The ultimate strength was calculated to be ∼10 GPa, depending on the diameter, which is in excellent agreement with the experiments. The dependence of the Young's modulus on nanowire diameter is explained in terms of surface compressive stress effects. The fracture behavior of nanowires was also found to be influenced by surface compressive stresses. Bond angle distribution analysis of various nanowires reveals significant compressive surface states, as evidenced by the appearance of a secondary peak in the Si-O-Si bond angle distribution at ∼97°, which is absent in bulk silica. The strain rate was found to have a negligible effect on the Young's modulus of the silica nanowires, but it has a critical role in determining their fracture mode.

A molecular dynamics study of ambient and high pressure phases of silica: structure and enthalpy variation with molar volume

Extensive molecular dynamics studies of 13 different silica polymorphs are reported in the isothermal-isobaric ensemble with the Parrinello-Rahman variable shape simulation cell. The van Beest-Kramer-van Santen (BKS) potential is shown to predict lattice parameters for most phases within 2%-3% accuracy, as well as the relative stabilities of different polymorphs in agreement with experiment. Enthalpies of high-density polymorphs - CaCl2-type, α-PbO2-type, and pyrite-type - for which no experimental data are available as yet, are predicted here. Further, the calculated enthalpies exhibit two distinct regimes as a function of molar volume-for low and medium-density polymorphs, it is almost independent of volume, while for high-pressure phases a steep dependence is seen. A detailed analysis indicates that the increased short-range contributions to enthalpy in the high-density phases arise not only from an increased coordination number of silicon but also shorter Si-O bond lengths. Our results indicate that amorphous phases of silica exhibit better optimization of short-range interactions than crystalline phases at the same density while the magnitude of Coulombic contributions is lower in the amorphous phase.

Comparison of model potentials for molecular-dynamics simulations of silica

Structural, thermomechanical, and dynamic properties of pure silica SiO2 are calculated with three different model potentials, namely, the potential suggested by van Beest, Kramer, and van Santen (BKS) [Phys. Rev. Lett. 64, 1955 (1990)], the fluctuating-charge potential with a Morse stretch term for the short-range interactions proposed by Demiralp, Cagin, and Goddard (DCG)[Phys. Rev. Lett. 82, 1708 (1999)], and a polarizable force field proposed by Tangney and Scandolo (TS) [J. Chem. Phys. 117, 8898 (2002)]. The DCG potential had to be modified due to flaws in the original treatment. While BKS reproduces many thermomechanical properties of different polymorphs rather accurately, it also shows qualitatively wrong trends concerning the phononic density of states, an absence of the experimentally observed anomaly in the c/a ratio at the quartz alpha-beta transition, pathological instabilities in the beta-cristobalite phase, and a vastly overestimated transition pressure for the stishovite I --> II transition. These shortcomings are only partially remedied by the modified DCG potential but greatly improved by the TS potential. DCG and TS both reproduce a pressure-induced transition from alpha-quartz to quartz II, predicted theoretically based on the BKS potential.

Infrared and Raman spectra of silica polymorphs from anab initioparametrized polarizable force field

The general aim of this study is to test the reliability of polarizable model potentials for the prediction of vibrational (infrared and Raman) spectra in highly anharmonic systems such as high temperature crystalline phases. By using an ab initio parametrized interatomic potential for SiO2 and molecular dynamics simulations, we calculate the infrared and Raman spectra for quartz, cristobalite, and stishovite at various thermodynamic conditions. The model is found to perform very well in the prediction of infrared spectra. Raman peak positions are also reproduced very well by the model; however, Raman intensities calculated by explicitly taking the derivative of the polarizability with respect to the atomic displacements are found to be in poorer agreement than intensities calculated using a parametrized "bond polarizability" model. Calculated spectra for the high temperature beta phases, where the role of dynamical disorder and anharmonicities is predominant, are found to be in excellent agreement with experiments. For the octahedral phases, our simulations are able to reproduce changes in the Raman spectra across the rutile-to-CaCl2 transition around 50 GPa, including the observed phonon softening.