Molecular Modelling of the Deformation Mechanisms Acting in Auxetic Silica

Molecular mechanics simulations have been performed to undertake a systematic investigation into the structure and mechanical properties of -cristobalite undergoing uniaxial loading along each of the 3 mutually orthogonal principal directions and also hydrostatic pressure loading. Simulations were performed using both the BKS and Burchart force-fields. The simulations indicate that pressure loading and uniaxial loading along the x 3 direction leads to uniform variation of the four independent Si-O-Si intertetrahedral angles, indicative of cooperative tetrahedral rotation about tetrahedral axes which transform the -cristobalite structure into the 'idealised' -cristobalite structure. Uniaxial loading along either of the transverse directions (x 1 and x 2 ) leads to a divergence of the intertetrahedral angles, consistent with tetrahedral rotation about the tetrahedral axes which transform the idealised -cristobalite structure into the 'ordered' -cristobalite structure. The data also indicate that a phase transition to one of the proposed phases may be induced by a negative hydrostatic pressure or tensile stress along x 3 . The phase transition is accompanied by a change in sign of some of the Poisson's ratios (i.e. from positive to negative). A negative hydrostatic pressure is also predicted to lead to conversion of initially positive to negative Poisson's ratio values (within the same phase).