We show that several structural and dynamic properties of silica, like the IR spectrum, distribution functions, mechanical properties, and negative thermal expansion computed using the MZHB potential agree well with experimental data. 
THERM
A potential with low point charges for pure siliceous zeolites.A modified CHARMM force-field (ZHB potential) with low point charges for silica was previously proposed by Zimmerman et al. (J. Chem. Theory Comput. 2011, 7, 1695). 
CHARMM