The thermodynamics of the Earth's interior is usually approached through the quasi-harmonic theory of solids, relying on the tacit assumption that explicit anharmonic contributions, which are difficult to calculate even for the most simple solids, are not important. The basic qualitative idea assumes that, as anharmonicity increases with temperature but decreases with pressure, in the Earth's interior the two effects are self-cancelling. A safer approach is to consider not the absolute but the homologous temperature, i.e. to take the melting temperature as a scaling factor. Since the geotherm lies rather close to melting, the relative importance of anharmonic terms in the thermodynamics of the Earth's interior is presumably high. At room pressure, the leading order anharmonic coefficients and their volume derivatives, as well as the quasi-harmonic terms from available thermodynamic data, can be estimated by fitting to experiment the theoretical expressions to the fourth order in the Hamiltonian of three functions of compressibility, heat capacity and thermal expansion, tailored to minimize the effect of experimental errors. The capability of this procedure to produce stable and accurate results has been previously tested through extensive applications to NaCl, and to several metals, chosen because of the large number of available independent data, which allowed the study in detail of the effect of external as well as internal errors. Comparatively accurate estimates were found to be possible. The procedure is here applied to forsterite, periclase, rutile and pyrope, allowing the estimation of the relative magnitude of the quasi-harmonic and leading order anharmonic terms. At temperatures of half the melting temperature the difference between the two is of the order of 20% on the Grüneisen parameter and just a few per cent on the heat capacity, which is nevertheless typically 10% from the classical 3 R limit.
Read full abstract