No abstract

We present a simple and computationally efficient classical atomistic model of silica in which the silicon and oxygen are simulated as hard spheres with four and two association sites, respectively. We have performed isobaric-isothermal Monte Carlo simulations to study the mechanical and phase behavior of this model. We have investigated solid phase structures of the model corresponding to quartz, cristobalite, and coesite, as well as some zeolite structures. For the model these phases are mechanically stable and highly incompressible. Ratios of zero-pressure bulk moduli and thermal expansion coefficients for alpha quartz, alpha cristobalite, and coesite are in quite good agreement with experimental values. The pressure-temperature phase diagram was constructed and shows three solid phases corresponding to cristobalite, quartz, and coesite, as well as a fluid or glass phase, behavior qualitatively similar to that seen for silica experimentally.

We have applied our previously reported model of silica based on low coordination and strong association ͓J. Chem. Phys. 121, 8415 ͑2004͔͒, to the calculation of phase stability of zeolite frameworks SOD, LTA, MFI, and FAU as silica polymorphs. We applied the method of Frenkel and Ladd for calculating free energies of these solids. Our model predicts that the MFI framework structure has a regime of thermodynamic stability at low pressures and above ϳ1400 K, relative to dense phases such as quartz. In contrast, our calculations predict that the less dense frameworks SOD, LTA, and FAU exhibit no regime of thermodynamic stability. We have also used our model to investigate whether templating extends the MFI regime of thermodynamic stability to lower temperatures, by considering templates with hard-sphere repulsions and mean-field attractions to silica. Within the assumptions of our model, we find that quartz remains the thermodynamically stable polymorph at zeolite synthesis temperatures ͑ϳ400 K͒ unless unphysically large template-silica attractions are assumed. These predictions suggest that some zeolites such as MFI may have regimes of thermodynamic stability even without template stabilization.

No abstract

No abstract