MODELING OF THE THERMAL-DEPENDENCE OF STRUCTURAL AND ELASTIC PROPERTIES OF CALCITE, CACO3

A computational method, based on the quasiharmonic approximation, has been computer-coded to calculate the temperature dependence of elastic constants and structural features of crystals. The model is applied t o calcite, CaCO3; an interatomic potential based on a C-O Morse functi on and Ca-O and O-O Born-type interactions, including a shell model fo r O, has been used. Equilibrations in the range 300-800 K reproduce th e experimental unit-cell edges and bond lengths within 1%. The simulat ed thermal expansion coefficients are 22.3 (parallel to c) and 2.6 (pe rpendicular to c), against 25.5 and -3.7 x 10(-6) K-1 experimental val ues, respectively. The thermal coefficients of elastic constants tend to be underestimated; for the bulk modulus, -2.3 against -3.7 x 10(-4) K-1 is obtained.