Introduction

Extreme temperature fluctuations due to daily illumination conditions affect the remote sensing response of minerals that build up Solar-System solid surfaces. Mercury, for instance, is a slow-rotation planet very close the Sun that undergoes to surface temperature ranges of hundreds of degrees1. Similarly, several minor rocky bodies like Near Earth asteroids experience large temperature variations during their inbound and outbound journeys around the Sun because of the great eccentricity of their orbits2. Such variations, which strongly affect the crystal structure of minerals, make the interpretation of remote sensing thermal infrared (TIR) spectra much more challenging. It is known, indeed, that in the middle portion of the electromagnetic spectrum (i.e. 100–5000 cm−1) absorption features arise by cation–anion and lattice vibrations of the crystal structure3. The way in which any crystal structure responds to temperature variations is by expanding or contracting4, therefore modifying bond lengths and energies absorbed during atom vibrations, and affecting shape and position of TIR absorption features.

The influence that the thermal expansion exerts on TIR spectra of minerals has been recently explored in support of the future employ of the BepiColombo Mercury Radiometer and Thermal Infrared Spectrometer (MERTIS)5, which was launched on October 20, 2018. Extended laboratory studies have been conducted on silicates6,7,8,9, which are the main constituents of Solar System rocky surfaces, in order to unveil their spectral behavior as a function of the daily variation of Mercury’s temperatures. In particular, it has been shown7 that, increasing the sample temperature from 298 to 773 K, the Si–O stretching vibration absorption bands of a Mg1.84Fe0.16SiO4 olivine reach markedly shifts toward lower wavenumbers, up to 8 cm−1. Concerning olivine, previous works investigated the effect of temperature on the cell size10,11 and on IR spectra12,13,14,15 as separate consequences. The authors7 investigated the unit-cell volume changes of the same Mg-rich olivine in situ, from 298 to 773 K, to correlate the olivine thermal expansion with the TIR spectra collected within the same range of temperatures. The authors confirmed that the unit-cell expansion couples with the spectrum shift. As this process is completely reversible (i.e. the volume expands upon heating and reduces back to the starting value after cooling) also the spectrum will show no signs of the heating process that has occurred. On the other hand, the same shift can be obtained without temperature variations but with a change of the Mg/Fe ratio16. Olivine, indeed, is an ideal solid solution between forsterite (Fo, Mg2SiO4), and fayalite (Fa, Fe2SiO4), and thus the unit-cell volume along such binary join changes linearly17 substituting Mg with Fe2+. As reported by these authors, the unit-cell volume increases linearly with increasing the Fe content through the following equation:

$${\text{V }}\left( {{\AA}^{{\text{3}}} } \right){\text{ }} = {\text{ 3}}0{\text{8}}.{\text{56 }}{-}{\text{ }}0.{\text{18}}0{\text{1 }} \times {\text{ }}\% {\text{Fo}}$$
(1)

Based on Eq. (1), we could observe an increase in the unit-cell volume because of the increasing Fe-content and/or because of an increase in temperature. An increase in temperature by about 420 K causes the same increase in volume caused by an enrichment in iron from Fo92 to Fo71 (i.e. from Mg1.84Fe0.16SiO4 to Mg1.42Fe0.68SiO4), without any temperature variations7. Therefore, spectral variations due to thermal expansion provide a temporary mismatch between remote sensing data and those of traditional room-temperature databases, affecting the correct interpretation of the surface composition as much as irreversible transformations (e.g., spectral flattening due to shock impacts) affect possible structural inferences (e.g., particle size). Knowing the effect of thermal expansion on TIR spectra becomes crucial to understand what happens to the unit-cell volume of olivine in the complete low–high temperature range.

In this work, we investigate the thermal expansion of the Mg1.84Fe0.16SiO4 olivine sample7 between 298 and 20 K by means of precise and accurate X-ray powder diffraction technique. These data, combined with those provided by7, cover the complete range of temperatures attainable in the inner Solar System and suggest which would be the best ranges to acquire remote sensing thermal infrared spectra from bodies characterized by extremely variable surface temperatures.

In addition, this study aims to find the best temperature range in which the thermal infrared spectra of olivine should be collected by spectrometers, whether they are remote or lander instruments, on a cold silicate surface of the Solar System in order to obtain a reliable mineralogical composition. Indeed, several cold surfaces of minor bodies, as well as of Mars and the Giant Planets’ moons, retain or could potentially retain silicates.

Results

The good agreement between the low-T (present study) and high-T7 data allowed us to reliably combine the two datasets into a single one. The evolution of unit-cell volumes of our sample as a function of temperature from 298 to 20 K is reported in Table 1, and shown in Fig. 1 together with data measured from 298 to 773 K7. The unit-cell volumes decrease non linearly from 298 to 80 K, whereas from this temperature down to 20 K we do not observe any volume decrease higher than the experimental uncertainties. The volume variation between 298 and 80 K (a temperature variation equal to about 220 K) is only 1.13 Å3 (i.e. 0.39%). The high-temperature study7 reported between 298 and 498 K (a temperature variation equal to 200 K) an almost double volume variation. Those measurements show how different are the thermal expansion behaviors at high and low temperature conditions for olivine.

Table 1 Unit-cell parameters as a function of temperature.
Figure 1
figure 1

Evolution of the unit-cell volume of olivine as a function of temperature. The unit-cell volume of olivine with composition Mg1.84Fe0.16SiO4 is plotted versus temperature between 20 and 773 K. Triangles are from7 whereas circles correspond to the data collected in this study. The blue curve represents the fit of all data (see text for more details).

From this extended dataset, we can determine the elastic behavior of olivine in the entire range of temperature occurring on Mercury surface (i.e. from about 700 K to about 70 K). On the other hand, our dataset can be easily applied to any Solar System surface that shows extreme temperatures. However, processing such extended dataset requires a careful data analysis in order to choose the correct equation of state (hereafter EoS) to be used. Following recent reviews18, fitting of the extended dataset was performed using a Kroll-type EoS, available on EoSFit7 software19,20 to guarantee robustness of the fittings. The resulting volume thermal expansion coefficient (αv), room-T unit-cell volume (V0) and the Einstein temperature (θE) values are reported in Table 2 (linear thermal expansion coefficients are reported in Table 3). The advantage of this model is that the only additional parameter is a characteristic Einstein temperature that can be estimated from entropy measurements21 or obtained by fitting V-T data18. In our study, we used the value of K′ equal to 4.47, which is the best estimate obtained in the recent review of the elastic properties of olivine18. Similar fitting strategies were used in19,22,23. The θE value of 450 K reported in Table 2 seems to be consistent with the values of 497(7) K for Fo9118, and with the estimate of 531 K for Fo100 from the measured entropy21 and with the decrease in θE with increasing Fe content10. The fit is reported in Fig. 1 along the entire low–high temperature data set7. Our volume thermal expansion values are in perfect agreement with previously published data at room-T and high-temperature10 and with a recent review on mantle olivine24. The value of volume thermal expansion coefficient, αV,298 = 2.673 10−5 K−1 (see Table 2) for the olivine studied in this work, is similar to that of other planetary rock-forming minerals like plagioclases, which show values in αV,298 between 1.6 and 2.8 10−5 K−125, and pyroxenes for which most values of αV is between 2.3 and 2.9 10−5 K−18,26,22 (see Table 4). Such similarities allow us to extend the following discussion carried out for olivine to most of the silicates considered the main constituting minerals of the surface of most Solar System bodies.

Table 2 Thermoelastic parameters for olivine Mg1.84Fe0.16SiO4.
Table 3 Linear thermal expansion coefficients of olivine Mg1.84Fe0.16SiO4.
Table 4 Volume thermal expansion coefficients.

Discussion

Our results (data between 20 and 298 K) combined with those by7 (data between 300 and 773 K) in Fig. 1 show how the volume of olivine displays larger volume changes at high-temperature values and smaller volume variations at low-temperature regimes. Since the shifts of the position of the TIR absorption bands strongly depend on unit-cell volume changes, greater shifts will occur at higher temperatures as well as negligible shifts will occur at very low temperatures. Such observation is extremely important in terms of spectral interpretation: it is likely that TIR spectra acquired on low-temperature surfaces would not show any observable shift between 20 and 120 K (Fig. 1, Table 1).

From 352 to 773 K7, the main absorption band of olivine occurring at about 940 cm−1 undergoes a total shift toward higher wavenumbers of 8 cm−1. For such temperature range, the unit-cell volume of the same sample undergoes an increase from 292.71 to 296.58 Å3, i.e. a positive increase of 3.87 Å3, which corresponds to an average increase of 0.0091 Å3/K if taken as linear. Instead, our new results on olivine (see Table 1) measured between 80 and 298 K show that the average increase is 0.0052 Å3/K (linear calculation) and becomes an even smaller volume variation between 160 and 80 K being 0.0030 Å3/K.

How these results could affect the TIR spectra of olivine at very low temperature values remain still unknown. However, assuming a linear behavior of the band shifts8,16, we could predict a proportional shift of 4.1 cm−1 between 298 and 80 K. We know that between 80 and 298 K the volume increase is almost halved with respect to the 352 to 773 K range, and thus it would be reasonable to expect a much smaller band shift. Based on our results, we could predict that such shift would be very small for temperature ranges between 160 and 298 K, and almost negligible between 80 and 160 K. If such prediction will be confirmed by the low-temperature TIR spectra dataset, we will be able to comprehend the thermal behavior of olivine and of any other silicates on Mercury and other Solar System bodies in order to retrieve their correct chemical composition bypassing any temperature influence.

A number of Solar System solid surfaces, indeed, have relatively low surface temperatures (e.g., Fig. 2) and, as well as the Mercury’s one, some of them are going to be explored by international missions. For instance, the Near Earth asteroid 1999 RQ36/101955 Bennu were recently globally scanned by the OSIRIS-REx Thermal Emission Spectrometer (OTES)23. Such orbit acquisitions reasonably depend on the signal-to-noise ratio, which partly depends on the body temperature. Acquisitions, indeed, were scheduled on the sunlit portion of the body in a well-calibrated range of temperatures, which are higher than those suggested in this work for negligible band position shifts (i.e., 80–160 K). The surface of Bennu can undergo midday temperatures between 200 and 420 K depending on the solar distance, and Planck functions having temperatures between 150 and 380 K were used to obtain OTES emissivity spectra. The average TIR spectrum of Bennu describes a radiative surface volumetrically dominated by silicates24, but no conclusive matches have been found yet among available aqueously-altered meteorites. Nevertheless, considering the best resolution in cm−1 achievable by hyperspectral instruments, we can be confident that any possible chemical composition derivable between 80 and 300 K will be not affected by thermal expansion effects.

Figure 2
figure 2

Surface temperature ranges of several Solar System bodies. The light cyan window encloses the values of volume expansion discussed in this work; orange bars show surface temperature range of several Solar System bodies that match the range discussed in this work; black circles indicate the evolution of volume as a function of temperature for olivine Mg1.84Fe0.16SiO4 here investigated.

Our analysis suggests that high-resolution TIR spectra of Main Belt asteroids, giant planet’s moons and Mars—of which surfaces do not exceed 300 K—provide reliable information about not only the silicate mineral identification but also about their chemical composition. Such information becomes critical in determining the geological environment of body evolution. On the other hand, Solar System hot surfaces like those of Mercury and most of the Near Earth asteroids would benefit of sunset TIR acquisitions as well as suitable high temperature spectra databases.

Materials and methods

In this study, we have investigated a natural olivine7 with composition Mg1.84Fe0.16SiO4. The olivine sample is from the high-pressure ultramafic nodules from the Newer Volcanics of Mt. Leura, Victoria, Australia25,26. The sample has been fully characterized by single-crystal X-ray diffraction (SCXRD) and X-ray powder diffraction (XRPD)7,25,26 and by Electron Microprobe Analysis (EMPA)29.

Precise measurements of unit-cell lattice parameters and volumes were performed using step-scan powder X-ray diffraction measurements in the temperature range 20 K—300 K. Each measurement was performed in a 2θ range between 10 and 100° (continuous scan) on a Philips X'Pert diffractometer (CuKα1,2 radiation, fixed 0.125° divergence slit, primary and secondary 0.04 rad Soller slits, X'Celerator detector) using a Janis CCS-250 cryostat at the Institute of Crystallography, RWTH Aachen, Germany. Temperature stability was achieved to better than 1 K. Lattice parameters were obtained from whole diffraction pattern refinement using the FULLPROF-suite31 and are reported at each temperature in Table 1.