Applications of Mössbauer Spectroscopy in Meteoritical and Planetary Science, Part II: Differentiated Meteorites, Moon, and Mars

: Mössbauer (nuclear γ -resonance) spectroscopy is a powerful technique which is actively used in various ﬁelds from physics and chemistry to biology and medicine. Rudolf L. Mössbauer, who observed nuclear γ -resonance and published his results in 1958, got a Nobel Prize in physics in 1961 for this discovery. 57 Fe is the most widely used nucleus in Mössbauer spectroscopy. Therefore, a large variety of compounds containing iron can be studied by Mössbauer spectroscopy. It is well known that planetary matter contains various iron-bearing phases and minerals. Therefore, the extraterrestrial material from different meteorites, asteroids, and planets can be studied using 57 Fe Mössbauer spectroscopy as an additional powerful technique. Two parts of this review consider the results of more than 50 years of experience of Mössbauer spectroscopy applied for the studies of various meteorites, soils, and rocks from the Moon and a recent investigation of the Martian surface using two rovers equipped with miniaturized Mössbauer spectrometers. Part I considered the results of Mössbauer spectroscopy of undifferentiated meteorites. Part II discusses the results of Mössbauer spectroscopy of differentiated meteorites formed in asteroids and protoplanets due to matter differentiation, as well as Lunar and Martian matter. phase, kamacite ( δ = − 0.002 mm/s, H eff = 332.2 kOe, A = 35%), ( ii ) γ -Fe(Ni) phase, taenite ( δ = − 0.04 mm/s, H eff = 306 kOe, A = 3%), ( iii ) troilite ( δ = 0.76 mm/s, H eff = 302 kOe, A = 2%), and ( iv ) cohenite (Fe, Ni, Co) 3 C ( δ = 0.11 mm/s, H eff = 199 kOe, A = 2%). Three quadrupole doublets were associated with ( i ) olivine ( δ = 1.151 mm/s, ∆ E Q = 2.983 mm/s, A = 32%), ( ii ) pyroxene, which was assigned to unresolved contributions from orthopyroxene and pigeonite ( δ = 1.147 mm/s, ∆ E Q = 2.139 mm/s, A = 24%), and ( iii ) an unknown Fe 3+ compound or pyrite ( δ = 0.52 mm/s, ∆ E Q = 0.67 mm/s, A = 2%). consider the meaning of a small peak at around +3 mm/s. Similar results were obtained for the Piplia Kalan spectra measured for two different lithologies A and B in [19]. The Mössbauer parameters for both lithologies were in fact the same and are as follows for the lithology A: δ = 1.13 mm/s, ∆ E Q = 2.40 mm/s, A = 22% (M1 sites), δ = 1.08 mm/s, ∆ E Q = 1.94 mm/s, A = 73% (M2 sites) and δ = 0.26 mm/s, ∆ E Q = 0.93 mm/s, A = 5% (ferric compound). the ratios of Fe 2+ occupancies among the and sites in silicate crystals, as well as the of K D and T eq , were using XRD and Mössbauer data. These ratios were the following: A M1 /A M2 ≈ 0.20 and X FeM1 /X FeM2 ≈ 0.18 for orthopyroxene, A M1 /A M2 ≈ 0.66 and X FeM1 /X FeM2 ≈ 0.66 for Ca-poor clinopyroxene, and A M1 /A M2 ≈ 2.40 and X FeM1 /X FeM2 ≈ 1.84 for Ca-rich clinopyroxene. These ratios agreed for the two independent techniques. The Fe 2+ and Mg 2+ cation equilibrium distributions among the M1 and M2 sites in orthopyroxene in Sariçiçek howardite related to its thermal history were K D = 0.12 and T eq = 886 K (XRD) and K D = 0.12 and T eq = 878 K (Mössbauer spectroscopy). The value of Fs was 28% (from MBD). The K D and T eq values also agreed for the two independent techniques. sextets one paramagnetic singlet. The 57 Fe parameters for sextets (singlet). sextets assigned with in the local iron microenvironment.


Introduction
Part I of this review [1] briefly presented the main features and applications of 57 Fe Mössbauer spectroscopy, together with the various studies on undifferentiated meteorites. The main parameters which can be obtained from 57 Fe Mössbauer spectroscopy are the following: (i) isomer shift δ, which is determined by the electron density on the 57 Fe nucleus and, therefore, is related to the iron valence/spin state; (ii) quadrupole splitting ∆E Q (quadrupole shift ε for magnetically split spectra, 2ε = ∆E Q ), which is related to the electric field gradient on the 57 Fe nucleus and reflects any tiny variations in the 57 Fe local microenvironment (both δ and ∆E Q determine the iron electron structure); (iii) magnetic hyperfine field H eff on the 57 Fe nucleus; (iv) the line width at a half maximum Γ; (v) the relative spectrum/component area A, which is proportional to the product of the Mössbauer effect probability (f -factor) and the number of 57 Fe nuclei in the corresponding compound.
Modern classification divides meteorites into differentiated and undifferentiated [2]. This section (Part II) presents the investigations conducted on differentiated meteorites, which consist of two main types: the primitive achondrites, which are supposed to be a result of low-degree partial melting of chondritic material, and achondrites, resulting from a high degree of chondrite melting in asteroids and protoplanets. The primitive achondrite

Primitive Achondrites
The first study of ureilites, ultramafic achondrites, the majority of which consist of olivine and uninverted pigeonite, i.e., a clinopyroxene subgroup with 5-25% of Ca 2+ fraction (see MBD), by Mössbauer spectroscopy was done in [6]. The authors studied two falls: Novo Urei (this meteorite gave the name "ureilites" for the group of primitive achondrites) and Haverö, and seven finds in Antarctica: ALHA84136 (not mentioned in MBD), ALHA81101, Pecora Escarpment (PCA) 82506, Elephant Moraine (EET) 87517, EET 87511, Meteorite Hills (MET) A78008, and ALHA77257. The room temperature Mössbauer spectra of Novo Urei, Haverö, ALHA84136, and ALHA81101 meteorites are shown in Figure 1. All Mössbauer spectra consist of one or two magnetic sextets assigned to Fe-Ni-Co alloys and Fe 3+ compounds (hematite or maghemite and goethite, ferrihydrite, akaganéite, or nanophase goethite as the authors supposed), two quadrupole doublets related to olivine and pyroxene, and one quadrupole doublet associated with the paramagnetic Fe 3+ compound. The δ and ∆E Q values of olivine are in the ranges~1.14-1.21 mm/s and~2.89-2.95 mm/s, respectively, while those for pyroxene (pigeonite) are in the ranges~1.12-1.18 mm/s and~1.88-2.35 mm/s, respectively. The δ and ∆E Q values for the paramagnetic ferric compounds are in the ranges~0.36-0.43 mm/s and~0.59-0.71 mm/s, respectively. The relative fractions of Fe were found (i) in metallic alloy in the range~0-47%, (ii) in olivine in the range~31-90%, (iii) in pyroxene in the range~10-69%, and (iv) in ferric compounds in the range~3-69%.

Achondrites (Stony Meteorites Except HED and Martian Meteorites)
The Zakłodzie meteorite, an ungrouped enstatite-rich achondrite, was studied by Mössbauer spectroscopy in [8]. The room temperature Mössbauer spectrum of this meteorite is shown in Figure 3. The authors used three fits with different components and compared their results with those obtained for the Abee EH4 enstatite chondrite in [9] (see Figure 8 in Part I [1]). The second and the third fits revealed the same components 1-6, while components 7 and 8 were different (the second fit is shown in Figure 3). The first six components, which are common for the two fits, were assigned by the authors to (1)
Mössbauer spectroscopy was applied to the study of D'Orbigny angrite, a relat rare type of basaltic achondrite, in [11,12]. The authors extracted and analyzed firs bulk and glass material [11] and then the druse clinopyroxene [12]. Glasses may be fo as a result of impact processes in the parent body or other factors. A comparison o room temperature Mössbauer spectra of the bulk and glass materials from D'Orbign grite is shown in Figure 4. Indicated components are the result of the authors' second fit (see [8] and the text above). Adapted from [8].
However, the results of these fits lead to many questions for the following reasons: the δ value for the γ-Fe(Ni) phase (2) is unusually large; the parameters of two components associated with pyrrhotite are very different and also differ from Mössbauer parameters for pyrrhotite obtained in [10]; plagioclase (6) is a mineral that does not contain iron and cannot give the Mössbauer spectrum, while the 57 Fe hyperfine parameters for component 6 are close to olivine; if singlets 7 and 8 in the third fit are considered as a quadrupole doublet, they can correspond to the paramagnetic ferric compound. Therefore, the Mössbauer parameters obtained in [8] should be carefully verified.
Mössbauer spectroscopy was applied to the study of D'Orbigny angrite, a relatively rare type of basaltic achondrite, in [11,12]. The authors extracted and analyzed first the bulk and glass material [11] and then the druse clinopyroxene [12]. Glasses may be formed as a result of impact processes in the parent body or other factors. A comparison of the room temperature Mössbauer spectra of the bulk and glass materials from D'Orbigny angrite is shown in Figure 4.
The Mössbauer spectrum of the bulk material from D'Orbigny angrite was fitted by the authors of [11] using four quadrupole doublets, two outer doublets of which were related to olivine while the two inner doublets were assigned to pyroxene. These four quadrupole doublets were associated with the M1 and M2 sites in both olivine and pyroxene. In contrast, the spectrum of glass material was decomposed with three quadrupole doublets only. Two doublets with the larger ∆E Q values were assigned to the M1 and M2 sites in pyroxene. The 57 Fe hyperfine parameters for these sites in pyroxene from the bulk and glass materials were close: (i) δ = 1.15 mm/s, ∆E Q = 2.56 mm/s (M1) and δ = 1.16 mm/s, ∆E Q = 2.09 mm/s (M2) for the bulk material and (ii) δ = 1.12 mm/s, ∆E Q = 2.70 mm/s, A = 17% (M1) and δ = 1.08 mm/s, ∆E Q = 2.07 mm/s, A = 40% (M2) for the glass material. The parameters for the third quadrupole doublet in the glass spectrum were determined as δ = 1.03 mm/s, ∆E Q = 1.53 mm/s, A = 43%. The latter component was not identified. However, there is a question concerning the fit of the glass spectrum because the third quadrupole doublet has a broad line width Γ = 0.74 mm/s which is twice larger than that for the M1 component (Γ = 0.36 mm/s).
The room temperature Mössbauer spectra of the druse clinopyroxene extracted from the D'Orbigny angrite represent an asymmetrical quadrupole doublet with a small admixture of the paramagnetic ferric component in both single crystal and powder forms [12]. The authors fitted their spectra using the quadrupole splitting distribution with a Voigt line shape instead of individual quadrupole doublets with a Lorentzian line shape. Therefore, the obtained Mössbauer parameters seem questionable and require new decomposition with individual quadrupole doublets which can be better related to the M1 and M2 sites in clinopyroxene.

Achondrites (HED)
HED meteorites, whose parent body is considered the asteroid (4) Vesta [3], consist of the main iron bearing phases such as Fe-rich pyroxene in eucrites, orthopyroxene in diogenite, and pyroxene (mainly orthopyroxene) in howardites which are a polymict-breccia, i.e., a mixture of eucrites and diogenites (see MBD). The first Mössbauer studies of HED meteorites were performed in [13][14][15][16][17][18]. The Mössbauer spectra of ALHA77256 diogenite and Stannern eucrite measured in [13] were different. One quadrupole doublet corresponding to pyroxene was observed in the spectra of the former meteorite samples and assigned to the M2 sites occupied by Fe 2+ . Four quadrupole doublets were revealed in the spectrum of the latter sample. These quadrupole doublets were related to the M1 and M2 sites and to the M1 sites in orthopyroxenes (probably pigeonite), as well as to ilmenite. Similar results were obtained for the paramagnetic components observed in the Stannern eucrite Mössbauer spectrum in [14]; however, in addition, the authors revealed the magnetic sextet with parameters corresponding to troilite. The room temperature Mössbauer spectrum of Ibitira eucrite consisted of three quadrupole doublets that were assigned to the M1 and M2 sites in pigeonite and to ilmenite [15]. Furthermore, the authors estimated the ratio of Fe 2+ populations in the M1 and M2 sites in pigeonite from Ibitira as A M2 /A M1 = 1.92 (in the 55 K Ibitira spectrum, the authors of [15,17] revealed two quadrupole doublets only related to the M1 and M2 sites in pigeonite with A M2 /A M1 = 2.11). Similar results were obtained in [16] for the Mössbauer spectra of Juvinas, ALHA80102, and EET 83232 eucrites with two quadrupole doublets related to the M1 and M2 sites in pyroxene and one quadrupole doublet assigned to ilmenite, as well as an additional quadrupole doublet associated with ferric compound. The Mössbauer spectra of Tatahouine, ALHA77256, and EETA79002 diogenites measured in [16] consisted of one quadrupole doublet related to the M2 sites in pyroxene and an additional quadrupole doublet assigned to ferric compounds. However, the room temperature Mössbauer spectrum of Kapoeta howardite was fitted using two quadrupole doublets related to the M1 and M2 sites in pyroxene [18].
The Mössbauer spectra of Cachari and Piplia Kalan eucrites (monomict breccias) measured in [11,19] are shown in Figure 5. In the spectrum of Cachari, the authors of [11] revealed three quadrupole doublets related to the M1 and M2 sites in pyroxenes and paramagnetic ferric compound. The 57 Fe hyperfine parameters for pyroxene were δ = 1.16 mm/s, ∆E Q = 2.46 mm/s (M1 sites) and δ = 1.13 mm/s, ∆E Q = 2.00 mm/s (M2 sites). However, the authors did not consider the meaning of a small peak at around +3 mm/s. Similar results were obtained for the Piplia Kalan spectra measured for two different lithologies A and B in [19]. The Mössbauer parameters for both lithologies were in fact the same and are as follows for the lithology A: δ = 1.13 mm/s, ∆E Q = 2.40 mm/s, A = 22% (M1 sites), δ = 1.08 mm/s, ∆E Q = 1.94 mm/s, A = 73% (M2 sites) and δ = 0.26 mm/s, ∆E Q = 0.93 mm/s, A = 5% (ferric compound).
inerals 2021, 11, x FOR PEER REVIEW 6 of 50 Mössbauer spectrum of Ibitira eucrite consisted of three quadrupole doublets that were assigned to the M1 and M2 sites in pigeonite and to ilmenite [15]. Furthermore, the authors estimated the ratio of Fe 2+ populations in the M1 and M2 sites in pigeonite from Ibitira as A M2 /A M1 = 1.92 (in the 55 K Ibitira spectrum, the authors of [15,17] revealed two quadrupole doublets only related to the M1 and M2 sites in pigeonite with A M2 /A M1 = 2.11). Similar results were obtained in [16] for the Mössbauer spectra of Juvinas, ALHA80102, and EET 83232 eucrites with two quadrupole doublets related to the M1 and M2 sites in pyroxene and one quadrupole doublet assigned to ilmenite, as well as an additional quadrupole doublet associated with ferric compound. The Mössbauer spectra of Tatahouine ALHA77256, and EETA79002 diogenites measured in [16] consisted of one quadrupole doublet related to the M2 sites in pyroxene and an additional quadrupole doublet assigned to ferric compounds. However, the room temperature Mössbauer spectrum of Kapoeta howardite was fitted using two quadrupole doublets related to the M1 and M2 sites in pyroxene [18]. The Mössbauer spectra of Cachari and Piplia Kalan eucrites (monomict breccias) measured in [11,19] are shown in Figure 5. In the spectrum of Cachari, the authors of [11] revealed three quadrupole doublets related to the M1 and M2 sites in pyroxenes and paramagnetic ferric compound. The 57 Fe hyperfine parameters for pyroxene were δ = 1.16 mm/s, ΔEQ = 2.46 mm/s (M1 sites) and δ = 1.13 mm/s, ΔEQ = 2.00 mm/s (M2 sites). However the authors did not consider the meaning of a small peak at around +3 mm/s. Similar results were obtained for the Piplia Kalan spectra measured for two different lithologies A and B in [19]. The Mössbauer parameters for both lithologies were in fact the same and are as follows for the lithology A: δ = 1.13 mm/s, ΔEQ = 2.40 mm/s, A = 22% (M1 sites), δ = 1.08 mm/s, ΔEQ = 1.94 mm/s, A = 73% (M2 sites) and δ = 0.26 mm/s, ΔEQ = 0.93 mm/s, A = 5% (ferric compound).
(a) (b) Figure 5. Comparison of the room temperature Mössbauer spectra of Cachari (a) and Piplia Kalan, lithology A (b) monomict eucrites. Indicated spectral components are the results of the fits. Adapted from [11,19].
A comparison of five cumulative and basaltic eucrites and one diogenite was done by Mössbauer spectroscopy in [20]. Cumulate eucrites and diogenites seem to have formed at depth in the asteroid (4) Vesta and crystallized quite slowly, while basaltic eucrites apparently formed at or near Vesta's surface and cooled relatively fast. The authors of [20] decomposed their Mössbauer spectra with different numbers of components: (i) one quadrupole doublet for cumulate eucrites and diogenite and (ii) two quadrupole doublets for basaltic eucrites as shown in Figure 6. Comparison of the room temperature Mössbauer spectra of Cachari (a) and Piplia Kalan, lithology A (b) monomict eucrites. Indicated spectral components are the results of the fits. Adapted from [11,19].
A comparison of five cumulative and basaltic eucrites and one diogenite was done by Mössbauer spectroscopy in [20]. Cumulate eucrites and diogenites seem to have formed at depth in the asteroid (4) Vesta and crystallized quite slowly, while basaltic eucrites apparently formed at or near Vesta's surface and cooled relatively fast. The authors of [20] decomposed their Mössbauer spectra with different numbers of components: (i) one quadrupole doublet for cumulate eucrites and diogenite and (ii) two quadrupole doublets for basaltic eucrites as shown in Figure 6. The Mössbauer parameters for two quadrupole doublets in the spectra of basaltic eucrites were related to the M1 and M2 sites in pyroxene: δ = 1.15 mm/s, ΔEQ = 2.42 mm/s, A = 35% (M1 sites), δ = 1.10 mm/s, ΔEQ = 2.01 mm/s, A = 65% (M2 sites) for Sioux County and δ = 1.14 mm/s, ΔEQ = 2.49 mm/s, A = 27% (M1 sites), δ = 1.10 mm/s, ΔEQ = 2.01 mm/s, A = 73% (M2 sites) for Millbillillie. The 57 Fe hyperfine parameters for the M2 sites in pyroxene in cumulate eucrites and diogenite were almost the same (δ = 1.10-1.11 mm/s, ΔEQ = 2.09-2.12 mm/s). The authors of [19] focused on the asymmetry of peak intensities in the Mössbauer spectra of cumulate eucrites which were fitted using one quadrupole doublet, while similar asymmetry in the spectra of basaltic eucrites was fitted using a superposition of two quadrupole doublets. The authors suggested that this asymmetry may be related to the pyroxene crystal orientation in the former group of meteorites. However, this should be verified by new studies, because the Mössbauer spectrum of Ibbenbüren diogenite demonstrated symmetrical peaks in one quadrupole doublet.
Earlier, one quadrupole doublet was applied to fit the Mössbauer spectra of three diogenites: ALHA77256, Tatahouine, and Bilanga [21]. These spectra are shown in Figure  7. The 57 Fe hyperfine parameters assigned to the M2 sites in pyroxene were almost the same: δ = 1.10-1.11 mm/s, ΔEQ = 2.10-2.12 mm/s. The effect of high pressure on Lohawat howardite and Piplia Kalan eucrite was studied by Mössbauer spectroscopy in [22,23]. The authors observed a decrease in the δ and ΔEQ values for the M2 sites in pyroxene in both meteorites with increasing pressure as shown in Figure 8. These results may be useful for the study of meteorites and their parent body collisions. The Mössbauer parameters for two quadrupole doublets in the spectra of basaltic eucrites were related to the M1 and M2 sites in pyroxene: δ = 1.15 mm/s, ∆E Q = 2.42 mm/s, A = 35% (M1 sites), δ = 1.10 mm/s, ∆E Q = 2.01 mm/s, A = 65% (M2 sites) for Sioux County and δ = 1.14 mm/s, ∆E Q = 2.49 mm/s, A = 27% (M1 sites), δ = 1.10 mm/s, ∆E Q = 2.01 mm/s, A = 73% (M2 sites) for Millbillillie. The 57 Fe hyperfine parameters for the M2 sites in pyroxene in cumulate eucrites and diogenite were almost the same (δ = 1.10-1.11 mm/s, ∆E Q = 2.09-2.12 mm/s). The authors of [19] focused on the asymmetry of peak intensities in the Mössbauer spectra of cumulate eucrites which were fitted using one quadrupole doublet, while similar asymmetry in the spectra of basaltic eucrites was fitted using a superposition of two quadrupole doublets. The authors suggested that this asymmetry may be related to the pyroxene crystal orientation in the former group of meteorites. However, this should be verified by new studies, because the Mössbauer spectrum of Ibbenbüren diogenite demonstrated symmetrical peaks in one quadrupole doublet.
Earlier, one quadrupole doublet was applied to fit the Mössbauer spectra of three diogenites: ALHA77256, Tatahouine, and Bilanga [21]. These spectra are shown in Figure 7. The 57 Fe hyperfine parameters assigned to the M2 sites in pyroxene were almost the same: δ = 1.10-1.11 mm/s, ∆E Q = 2.10-2.12 mm/s. The Mössbauer parameters for two quadrupole doublets in the spectra of basaltic eucrites were related to the M1 and M2 sites in pyroxene: δ = 1.15 mm/s, ΔEQ = 2.42 mm/s, A = 35% (M1 sites), δ = 1.10 mm/s, ΔEQ = 2.01 mm/s, A = 65% (M2 sites) for Sioux County and δ = 1.14 mm/s, ΔEQ = 2.49 mm/s, A = 27% (M1 sites), δ = 1.10 mm/s, ΔEQ = 2.01 mm/s, A = 73% (M2 sites) for Millbillillie. The 57 Fe hyperfine parameters for the M2 sites in pyroxene in cumulate eucrites and diogenite were almost the same (δ = 1.10-1.11 mm/s, ΔEQ = 2.09-2.12 mm/s). The authors of [19] focused on the asymmetry of peak intensities in the Mössbauer spectra of cumulate eucrites which were fitted using one quadrupole doublet, while similar asymmetry in the spectra of basaltic eucrites was fitted using a superposition of two quadrupole doublets. The authors suggested that this asymmetry may be related to the pyroxene crystal orientation in the former group of meteorites. However, this should be verified by new studies, because the Mössbauer spectrum of Ibbenbüren diogenite demonstrated symmetrical peaks in one quadrupole doublet.
Earlier, one quadrupole doublet was applied to fit the Mössbauer spectra of three diogenites: ALHA77256, Tatahouine, and Bilanga [21]. These spectra are shown in Figure  7. The 57 Fe hyperfine parameters assigned to the M2 sites in pyroxene were almost the same: δ = 1.10-1.11 mm/s, ΔEQ = 2.10-2.12 mm/s. The effect of high pressure on Lohawat howardite and Piplia Kalan eucrite was studied by Mössbauer spectroscopy in [22,23]. The authors observed a decrease in the δ and ΔEQ values for the M2 sites in pyroxene in both meteorites with increasing pressure as shown in Figure 8. These results may be useful for the study of meteorites and their parent body collisions. The effect of high pressure on Lohawat howardite and Piplia Kalan eucrite was studied by Mössbauer spectroscopy in [22,23]. The authors observed a decrease in the δ and ∆E Q values for the M2 sites in pyroxene in both meteorites with increasing pressure as shown in Figure 8. These results may be useful for the study of meteorites and their parent body collisions. The Kapoeta howardite was studied by both XRD and Mössbauer spectroscopy in order to analyze its thermal history in [24]. The calculated values of the cation distribution coefficient KD and closure temperature TCl (see Part I, Section 6.7 [1]) for orthopyroxene in Kapoeta were different for XRD and Mössbauer data, e.g., TCl = 684 K and 630 K from XRD and TCl = 1017 K from Mössbauer data. The obtained values of A M1 and A M2 for the M1 and M2 sites in orthopyroxene were ~20% and ~80%. However, these values disagreed with data A M1 = ~41% and A M2 = ~59% known from a previous Mössbauer study of Kapoeta in [18], which may be not correct.
The Kapoeta howardite was studied by both XRD and Mössbauer spectroscopy in order to analyze its thermal history in [24]. The calculated values of the cation distribution coefficient K D and closure temperature T Cl (see Part I, Section 6.7 [1]) for orthopyroxene in Kapoeta were different for XRD and Mössbauer data, e.g., T Cl = 684 K and 630 K from XRD and T Cl = 1017 K from Mössbauer data. The obtained values of A M1 and A M2 for the M1 and M2 sites in orthopyroxene were~20% and~80%. However, these values disagreed with data A M1 =~41% and A M2 =~59% known from a previous Mössbauer study of Kapoeta in [18], which may be not correct.
Mössbauer spectroscopy with a high velocity resolution was used to study the bulk interior and the fusion crust of the Sariçiçek howardite in [25]. These Mössbauer spectra are shown in Figure 9. Two magnetic sextets, nine quadrupole doublets, and one paramagnetic singlet were revealed in the spectrum of the bulk interior. On the base of the 57 Fe hyperfine parameters, these spectral components were assigned to the following phases: (i) α-Fe(Ni, Co), δ = −0.002 mm/s, H eff = 333. 8  close to that obtained for Kapoeta in [24]. Using approaches considered in Part 1, Sections 6.6 and 6.7 [1], the ratios of Fe 2+ occupancies among the M1 and M2 sites in silicate crystals, as well as the values of K D and T eq , were calculated using both XRD and Mössbauer data. These ratios were the following: A M1 /A M2 ≈ 0.20 and X Fe M1 /X Fe M2 ≈ 0.18 for orthopyroxene, A M1 /A M2 ≈ 0.66 and X Fe M1 /X Fe M2 ≈ 0.66 for Ca-poor clinopyroxene, and A M1 /A M2 ≈ 2.40 and X Fe M1 /X Fe M2 ≈ 1.84 for Ca-rich clinopyroxene. These ratios agreed for the two independent techniques. The Fe 2+ and Mg 2+ cation equilibrium distributions among the M1 and M2 sites in orthopyroxene in Sariçiçek howardite related to its thermal history were K D = 0.12 and T eq = 886 K (XRD) and K D = 0.12 and T eq = 878 K (Mössbauer spectroscopy). The value of Fs was 28% (from MBD). The K D and T eq values also agreed for the two independent techniques. Sariçiçek howardite related to its thermal history were KD = 0.12 and Teq = 886 K (XRD) and KD = 0.12 and Teq = 878 K (Mössbauer spectros copy). The value of Fs was 28% (from MBD). The KD and Teq values also agreed for the two independent techniques.
(a) (b) Figure 9. Room temperature Mössbauer spectra of Sariçiçek howardite measured with a high-velocity resolution: the bulk interior (a) and the fusion crust (b). Indicated components are the results of the best fits. The differential spectra are shown on the bottom. Adapted from [25].
The first study of the fusion crust by Mössbauer spectroscopy was carried out for Stannern eucrite in [14]. The authors used a 1:1 mixture of the bulk interior and the fusion crust and observed in the Mössbauer spectrum two quadrupole doublets related to the M1 and M2 sites in pyroxene, two quadrupole doublets assigned to ilmenite and ferric compound, and two magnetic sextets associated with troilite and Fe-Ni-Co alloy. In the study of Sariçiçek fusion crust [25], the sample contained only the fusion crust material. I was a glass-like matter as demonstrated by XRD. The Mössbauer spectrum of Sariçiçek fusion crust (see Figure 9b) did not contain any magnetic sextets but consisted of nine quadrupole doublets revealed from the fit with known components of orthopyroxene and unknown ferrous and ferric compounds. Two quadrupole doublets with the largest quad rupole splitting were related to the M1 and M2 sites in orthopyroxene (δ = 1.088 mm/s ΔEQ = 2.699 mm/s, A = ~7.3% for the M1 sites; δ = 1.139 mm/s, ΔEQ = 2.099 mm/s, A = ~21.7% for the M2 sites). The ratio of Fe 2+ occupancies among the M1 and M2 sites in orthopyroxene was A M1 /A M2 ≈ 0.34, which is larger than that obtained for the bulk interior The values of KD and TFC calculated for Sariçiçek fusion crust were 0.23 and 1177 K. The latter indicates that the fast cooling and solidification of the molten fusion crust was started in the temperature range 1100-1200 K.

Achondrites (Lunar and Martian Meteorites)
Lunar and Martian stony achondrites are rocks which were ejected from the Moon The first study of the fusion crust by Mössbauer spectroscopy was carried out for Stannern eucrite in [14]. The authors used a 1:1 mixture of the bulk interior and the fusion crust and observed in the Mössbauer spectrum two quadrupole doublets related to the M1 and M2 sites in pyroxene, two quadrupole doublets assigned to ilmenite and ferric compound, and two magnetic sextets associated with troilite and Fe-Ni-Co alloy. In the study of Sariçiçek fusion crust [25], the sample contained only the fusion crust material. It was a glass-like matter as demonstrated by XRD. The Mössbauer spectrum of Sariçiçek fusion crust (see Figure 9b) did not contain any magnetic sextets but consisted of nine quadrupole doublets revealed from the fit with known components of orthopyroxene and unknown ferrous and ferric compounds. Two quadrupole doublets with the largest quadrupole splitting were related to the M1 and M2 sites in orthopyroxene (δ = 1.088 mm/s, ∆E Q = 2.699 mm/s, A =~7.3% for the M1 sites; δ = 1.139 mm/s, ∆E Q = 2.099 mm/s, A =~21.7% for the M2 sites). The ratio of Fe 2+ occupancies among the M1 and M2 sites in orthopyroxene was A M1 /A M2 ≈ 0.34, which is larger than that obtained for the bulk interior. The values of K D and T FC calculated for Sariçiçek fusion crust were 0.23 and 1177 K. The latter indicates that the fast cooling and solidification of the molten fusion crust was started in the temperature range 1100-1200 K.

Achondrites (Lunar and Martian Meteorites)
Lunar and Martian stony achondrites are rocks which were ejected from the Moon and Mars by impacts and later fell to the Earth as meteorites (see MBD). Martian meteorites are divided into three main groups: shergottites (named after the Shergotty meteorite), nakhlites (named after the Nakhla meteorite), and chassignites (named after the Chassigny meteorite). These groups were named as SNC meteorites as an abbreviation using the first letters of the above-mentioned meteorite names. Until now, no studies on Lunar meteorites using Mössbauer spectroscopy have been carried out. On the contrary, Mössbauer investigations of Martian meteorites and their extracted iron-bearing phases have been carried out since mid-1980s (see, e.g., [6,16,26]). In the first study [26] of three SNC achondrites (Zagami and EETA79001 shergottites and Nakhla) the authors identified different iron-bearing phases in these meteorites. The spectrum of Zagami consisted of two quadrupole doublets related to the M1 and M2 sites in pyroxene (δ = 1.16 mm/s, ∆E Q = 2.47 mm/s and δ = 1.14 mm/s, ∆E Q = 2.01 mm/s, respectively). The Mössbauer spectrum of EETA79001 was decomposed into three quadrupole doublets, two of which were assigned to the M1 and M2 sites in pyroxene (δ = 1.16 mm/s, ∆E Q = 2.40 mm/s and δ = 1.14 mm/s, ∆E Q = 2.27 mm/s, respectively) while one was assigned to the M1 site in olivine (δ = 1.15 mm/s, ∆E Q = 2.88 mm/s). The spectrum of Nakhla was much more complex and consisted of magnetic and paramagnetic components. Therefore, the authors also measured the paramagnetic and magnetic separation. The Mössbauer spectrum of the paramagnetic separate from Nakhla was decomposed into four quadrupole doublets: two different M1 sites and M2 sites in pyroxene (δ = 1.17 mm/s, ∆E Q = 2.60 mm/s; δ = 1.16 mm/s, ∆E Q = 2.03 mm/s and δ = 1.16 mm/s, ∆E Q = 1.83 mm/s, respectively) and M1 sites in olivine (δ = 1.17 mm/s, ∆E Q = 2.90 mm/s). The spectrum of magnetic separate from Nakhla was fitted using two magnetic sextets for the (A) and [B] sites in magnetite (δ = 0.27 mm/s, H eff = 475 kOe and δ = 0.68 mm/s, H eff = 463 kOe, respectively) and four quadrupole doublets assigned to ilmenite (δ = 1.06 mm/s, ∆E Q = 0.68 mm/s), ulvöspinel Fe 2 TiO 4 (δ = 1.16 mm/s, ∆E Q = 1.52 mm/s), and residual pyroxene. Investigations of [6,16] were dedicated to the SNC meteorite weathering and detection of ferric compounds in order to consider the possible preterrestial oxidation, i.e., oxidation on Mars.
Ten various SNC meteorites with different lithologies and extracted olivine and pyroxene were studied by Mössbauer spectroscopy in [28] in order to analyze ferric compounds. The Mössbauer spectra of the bulk material from selected SNC meteorites are shown in Figure 11. Ten various SNC meteorites with different lithologies and extracted olivine and pyroxene were studied by Mössbauer spectroscopy in [28] in order to analyze ferric compounds. The Mössbauer spectra of the bulk material from selected SNC meteorites are shown in Figure 11. These spectra show different fractions of olivine and pyroxene in these meteorites and different spectral components in the spectra of the bulk materials and separated olivine and pyroxene. The presence of ferric compounds was observed in the Mössbauer spectra of some SNC meteorites. These spectra were fitted using the quadrupole splitting distribution that lead to some questions because, in the case of the thin absorber with well- Ten various SNC meteorites with different lithologies and extracted olivine and p roxene were studied by Mössbauer spectroscopy in [28] in order to analyze ferric co pounds. The Mössbauer spectra of the bulk material from selected SNC meteorites a shown in Figure 11. These spectra show different fractions of olivine and pyroxene in these meteori and different spectral components in the spectra of the bulk materials and separated o vine and pyroxene. The presence of ferric compounds was observed in the Mössbau spectra of some SNC meteorites. These spectra were fitted using the quadrupole splitti distribution that lead to some questions because, in the case of the thin absorber with we These spectra show different fractions of olivine and pyroxene in these meteorites and different spectral components in the spectra of the bulk materials and separated olivine and pyroxene. The presence of ferric compounds was observed in the Mössbauer spectra of some SNC meteorites. These spectra were fitted using the quadrupole splitting distribution that lead to some questions because, in the case of the thin absorber with welldefined crystal sites for the 57 Fe nuclei, the fit with individual quadrupole doublets directly related to these sites seems much more reasonable than the fit by the model-independent distribution function. It is noteworthy that the Mössbauer spectra of Zagami shergottite measured in [27,28] were different.
Investigations of MIL 03346 nakhlite with dominant clinopyroxene content and extracted clinopyroxene, as well as NWA 2737 chassignite, which is rich with olivine and extracted olivine, were studied by Mössbauer spectroscopy in [29][30][31], respectively. The Mössbauer spectra of the bulk materials and extracted silicates from the MIL 03346 and NWA 2737 meteorites are shown in Figure 12.
Minerals 2021, 11, x FOR PEER REVIEW 12 of 50 defined crystal sites for the 57 Fe nuclei, the fit with individual quadrupole doublets directly related to these sites seems much more reasonable than the fit by the model-independent distribution function. It is noteworthy that the Mössbauer spectra of Zagami shergottite measured in [27] and [28] were different. Investigations of MIL 03346 nakhlite with dominant clinopyroxene content and extracted clinopyroxene, as well as NWA 2737 chassignite, which is rich with olivine and extracted olivine, were studied by Mössbauer spectroscopy in [29] and [30,31], respectively. The Mössbauer spectra of the bulk materials and extracted silicates from the MIL 03346 and NWA 2737 meteorites are shown in Figure 12. The spectra of the bulk materials were similar to those of the extracted silicates, confirming the dominant content of clinopyroxene in MIL 03346 and olivine in NWA 2737. Measuring the low temperature Mössbauer spectra of NWA 2737 samples did not reveal any magnetically split components. The room temperature Mössbauer parameters obtained for MIL 03346 samples were as follows: (i) Fe 2+ (δ = 1.14 mm/s, ΔEQ = 2.60 mm/s, A = 23.1% for the bulk material and δ = 1.15 mm/s, ΔEQ = 2.68 mm/s, A = 19.8% for the ex-  The spectra of the bulk materials were similar to those of the extracted silicates, confirming the dominant content of clinopyroxene in MIL 03346 and olivine in NWA 2737. Measuring the low temperature Mössbauer spectra of NWA 2737 samples did not reveal any magnetically split components. The room temperature Mössbauer parameters obtained for MIL 03346 samples were as follows: (i) Fe 2+ (δ = 1.14 mm/s, ∆E Q = 2.60 mm/s, A = 23.1% for the bulk material and δ = 1.15 mm/s, ∆E Q = 2.68 mm/s, A = 19.8% for the extracted clinopyroxene), (ii) Fe 2+ (δ = 1.13 mm/s, ∆E Q = 1.95 mm/s, A = 54.2% for the bulk material and δ = 1.15 mm/s, ∆E Q = 1.97 mm/s, A = 56.1% for the extracted clinopyroxene), and (iii) Fe 3+ compounds (δ = 0.42 mm/s, ∆E Q = 0.66 mm/s, A = 22.8% for the bulk material and δ = 0.27 mm/s, ∆E Q = 0.61 mm/s, A = 24.2% for the extracted clinopyroxene). The authors of [29] revealed a large content of Fe 3+ compounds and assigned them to augite. The room temperature Mössbauer parameters which were obtained for NWA 2737 samples by the authors of [30] were the following: (i) Fe 2+ in olivine (δ = 1.18 mm/s, ∆E Q = 3.07 mm/s, A = 28% for the bulk material and δ = 1.08 mm/s, ∆E Q = 2.91 mm/s, A = 23% for extracted olivine), (ii) Fe 2+ in olivine (δ = 1.15 mm/s, ∆E Q = 2.84 mm/s, A = 38% for the bulk material and δ = 1.16 mm/s, ∆E Q = 2.98 mm/s, A = 55% for the extracted olivine), (iii) Fe 2+ in pyroxene (δ = 1.16 mm/s, ∆E Q = 2.64 mm/s, A = 21% for the bulk material only), (iv) Fe 2+ in chromite (δ = 0.99 mm/s, ∆E Q = 1.72 mm/s, A = 3% for the bulk material and δ = 0.99 mm/s, ∆E Q = 1.97 mm/s, A = 6% for the extracted olivine), and (v) Fe 3+ in olivine, chromite (δ = 0.48 mm/s, ∆E Q = 0.82 mm/s, A = 10% for the bulk material and δ = 0.48 mm/s, ∆E Q = 0.84 mm/s, A = 15% for the extracted olivine). It should be noted that these results for chromite are quite different from those for the chromite component observed in the room temperature Mössbauer spectra of some ordinary chondrites (see Figures 29 and 34 in Part I [1]) and Sariçiçek howardite (see Figure 9a above), which is a paramagnetic singlet (see, e.g., [32] and references therein and [33]). However, the presence of Fe 2+ and Fe 3+ in the tetrahedral or octahedral sites, respectively, is possible in the synthetic spinels Fe 1+x Cr 2−x O 4 if x = 0 or in the natural Cr-bearing spinels, such as magnesiochromite containing Fe and Al, by the observation of several quadrupole doublets in the Mössbauer spectra (see [33,34]). In general, the term "chromite" is used for Cr-bearing spinels [34]. Therefore, this general term mixes with the term meaning the stoichiometric or close to stoichiometric chromite with small Cr substitutions by Al, Mg, Ti, and Mn atoms. The Mössbauer spectra of the latter chromites demonstrate a single peak with some variations in δ and Γ values. The authors of [29]  These results for Cr-bearing spinels can be compared with those for chromite separate from the ALH 84001 Martian meteorite which does not belong to the SNC group, and it is considered as an orthopyroxene-rich Martian meteorite (OPX). The Mössbauer spectrum of chromite extracted from ALH 84001 OPX and measured in [35] is shown in Figure 13. This spectrum shows two quadrupole doublets corresponding to the ferrous and ferric compounds. It is possible that the Cr-bearing spinel extracted from ALH 84001 OPX also has a complex formula containing Fe 2+ in tetrahedral and Fe 3+ in octahedral sites and some substitutions of Cr by Fe, Al, Ti, and other metals presented in this spinel. if x ≠ 0 or in the natural Cr-bearing spinels, such as magnesi Al, by the observation of several quadrupole doublets in [33,34]). In general, the term "chromite" is used for Cr-bearin general term mixes with the term meaning the stoichiomet chromite with small Cr substitutions by Al, Mg, Ti, and Mn a of the latter chromites demonstrate a single peak with some The authors of [29] determined the NWA 2737 chromite fo (Fe 2+ 0.05Ti0.05Fe 3+ 0.10Al0.34Cr1.46)O4. These results for Cr-bearing s those for chromite separate from the ALH 84001 Martian met to the SNC group, and it is considered as an orthopyroxene-r The Mössbauer spectrum of chromite extracted from ALH [35] is shown in Figure 13. This spectrum shows two quadru to the ferrous and ferric compounds. It is possible that the from ALH 84001 OPX also has a complex formula containin in octahedral sites and some substitutions of Cr by Fe, Al, T in this spinel. The authors used two fragments of Y-984028 from the bul including the fusion crust. Moreover, they additionally ext from both chips. The room temperature Mössbauer spectra as olivine and pyroxene extracts from the interior and ext shown in Figure 14. These spectra were measured in 2048 therefore, can be considered as measured with a high veloci tion of the velocity reference signal on 2 10 steps. A detailed study of Yamato 984028 (Y-984028), a lherzolitic shergottite from Antarctica, was carried out with various techniques including Mössbauer spectroscopy in [36]. The authors used two fragments of Y-984028 from the bulk interior and exterior chips including the fusion crust. Moreover, they additionally extracted olivine and pyroxene from both chips. The room temperature Mössbauer spectra of the bulk material, as well as olivine and pyroxene extracts from the interior and exterior chips of Y-984028, are shown in Figure 14. These spectra were measured in 2048 channels before folding and, therefore, can be considered as measured with a high velocity resolution with discretization of the velocity reference signal on 2 10 steps. The 57 Fe hyperfine parameters for the same spectral components in the studied samples were very similar while the relative areas (the relative fractions) of these components in some cases were different. Mössbauer parameters for the bulk samples were as follows: (i) Fe 2+ in pyroxene (δ = 1.12 mm/s, ΔEQ = 2.04 mm/s, A = 22% for the exterior material and δ = 1.10 mm/s, ΔEQ = 2.25 mm/s, A = 31% for the interior material), (ii) Fe 2+ in olivine (δ = 1.16 mm/s, ΔEQ = 2.99 mm/s, A = 19% for the exterior material and δ = 1.17 mm/s, ΔEQ = 2.99 mm/s, A = 22% for the interior material), (iii) Fe 2+ in olivine ± pyroxene (δ = 1.11 mm/s, ΔEQ = 2.82 mm/s, A = 48% for the exterior material and δ = 1.11 mm/s, ΔEQ = 2.78 mm/s, A = 39% for the interior material), and (iv) Fe 3+ (δ = 0.17 mm/s, ΔEQ = 0.55 mm/s, A = 11% for the exterior material and δ = 0.51 mm/s, ΔEQ = 0.73 mm/s, A = 9% for the interior material). Mössbauer parameters for the olivine separates were as follows: (i) Fe 2+ in olivine (δ = 1.18 mm/s, ΔEQ = 2.98 mm/s, A = 24% for the exterior material and δ = 1.16 mm/s, ΔEQ = 3.01 mm/s, A = 39% for the interior material), (ii) Fe 2+ in olivine ± pyroxene (δ = 1.12 mm/s, ΔEQ = 2.74 mm/s, A = 34% for the exterior material and δ = 1.12 mm/s, ΔEQ = 2.82 mm/s, A = 48% for the interior material), and (iii) Fe 3+ (δ = 0.29 mm/s, ΔEQ = 0.71 mm/s, A = 42% for the exterior material and δ = 0.36 mm/s, ΔEQ = 0.61 mm/s, A = 13% for the interior material). Mössbauer parameters for the pyroxene separates were as follows: (i) Fe 2+ in pyroxene (δ = 1.13 mm/s, ΔEQ = 2.17 mm/s, A = 25% for the exterior material and δ = 1.13 mm/s, ΔEQ = 2.17 mm/s, A = 56% for the interior material), (ii) Fe 2+ in olivine ± pyroxene (δ = 1.14 mm/s, ΔEQ = 2.83 mm/s, A = 49% for the exterior material and δ = 1.14 mm/s, ΔEQ = 2.83 mm/s, A = 32% for the interior material), (iii) Fe 3+ (δ = 0.15 mm/s, ΔEQ = 0.55 mm/s, A = 14% for the exterior material and δ = 0.11 mm/s, ΔEQ = 0.42 mm/s, A = 12% for the interior material), and (iv) Fe 3+ (δ = 0.28 mm/s, ΔEQ = 0.96 mm/s, A = 12% for the exterior material). The spectral component marked by the authors as olivine ± pyroxene indicated a quadrupole doublet with probably unresolved contributions from both olivine and pyroxene. The exterior olivine separate contained more Fe 3+ fraction than the interior olivine separate. This component was related to phyllosilicate.
Mössbauer spectroscopy was used for the determination of Fe 3+ in the augite (clinopyroxene) core-crystal from MIL 03346 nakhlite in [37]. The room temperature Mössbauer spectra of the MIL 03346 augite single crystal with two different orientations are shown in Figure 15. The authors determined that the content of ferric compound was around 7%, which is significantly smaller than the ferric compound fraction obtained in [29] for clinopyroxene separate from MIL 03346.  Figure 15. The authors determined that the content of ferric compound was around 7%, which is significantly smaller than the ferric compound fraction obtained in [29] for clinopyroxene separate from MIL 03346.
(a) (b) Figure 15. Room temperature Mössbauer spectra of MIL 03346 augite (clinopyroxene) single crystal at two different orientations (a,b). Indicated components are the results of the fits. The ferric component is shaded in red. The differential spectra are shown on the top. Adapted from [37].

Pallasites
Pallasites are stony-iron meteorites that formed as a result of the solidification of the mixture of molten Fe-Ni-Co alloy and olivine-rich rock fragments. This process may take place in differentiated asteroids and protoplanets, which contain the molten metallic core and olivine-rich mantle or a molten metallic layer between the solid metallic core and the olivine-rich mantel. After destructive impact with another large body in space, this mixture can be formed and slowly cooled. Therefore, pallasites consist of Fe-Ni-Co alloy matrix with stony fragments. There are several groups of pallasites, the largest of which is named the main group pallasites (PMG).
The first study of olivine extracted from Omolon PMG and Seymchan PMG using Mössbauer spectroscopy with a high velocity resolution was performed in [38][39][40]. The Mössbauer spectra of olivine separates from both pallasites measured at 295 and 90 K are shown in Figure 16. These spectra demonstrate well-known peak asymmetry, which becomes inverse at low temperature. These spectra were decomposed with two main quadrupole doublets marked 1 and 2 in Figure 16 and related to the M1 and M2 sites in olivine. All spectra also contain the minor quadrupole doublet 3 assigned to unknown ferrous compound X, while the 295 K spectrum of Omolon only shows a very small quadrupole doublet 4 associated with ferric compound. The room temperature Mössbauer parameters for both olivines were the following: (1) δ = 1.156 mm/s, ΔEQ = 3.020 mm/s, A = 51% for Omolon and δ = 1.158 mm/s, ΔEQ = 3.058 mm/s, A = 50% for Seymchan; (2) δ = 1.117 mm/s, ΔEQ = 2.849 mm/s, A = 43% for Omolon and δ = 1.137 mm/s, ΔEQ = 2.876 mm/s, A = 45% for Seymchan; (3) δ = 1.007 mm/s, ΔEQ = 3.108 mm/s, A = 5% for Omolon and δ = 0.952 mm/s, ΔEQ = 3.032 mm/s, A = 5% for Seymchan; (4) δ = 0.389 mm/s, ΔEQ = 0.570 mm/s, A = 1% for Omolon. Mössbauer parameters at 90 K were also similar for both olivines. It was observed that the δ values for the 57 Fe in the M1 and M2 sites in both meteorites showed different tendencies of decrease with the increase in temperature. The decrease in the isomer shift was related to the second-order Doppler shift (see details in, e.g., [41]), while the differences in the observed tendencies indicated the different Mössbauer temperatures for the 57 Fe in the M1 and M2 sites in olivines from Omolon and Seymchan, resulting from the different 57 Fe mean square velocities in these sites. Using the approach described in Part I, Section 6.7, for the ordinary chondrites [1], the values of KD and Teq were estimated

Pallasites
Pallasites are stony-iron meteorites that formed as a result of the solidification of the mixture of molten Fe-Ni-Co alloy and olivine-rich rock fragments. This process may take place in differentiated asteroids and protoplanets, which contain the molten metallic core and olivine-rich mantle or a molten metallic layer between the solid metallic core and the olivine-rich mantel. After destructive impact with another large body in space, this mixture can be formed and slowly cooled. Therefore, pallasites consist of Fe-Ni-Co alloy matrix with stony fragments. There are several groups of pallasites, the largest of which is named the main group pallasites (PMG).
The first study of olivine extracted from Omolon PMG and Seymchan PMG using Mössbauer spectroscopy with a high velocity resolution was performed in [38][39][40]. The Mössbauer spectra of olivine separates from both pallasites measured at 295 and 90 K are shown in Figure 16. These spectra demonstrate well-known peak asymmetry, which becomes inverse at low temperature. These spectra were decomposed with two main quadrupole doublets marked 1 and 2 in Figure 16 and related to the M1 and M2 sites in olivine. All spectra also contain the minor quadrupole doublet 3 assigned to unknown ferrous compound X, while the 295 K spectrum of Omolon only shows a very small quadrupole doublet 4 associated with ferric compound. The room temperature Mössbauer parameters for both olivines were the following: (1) δ = 1.156 mm/s, ∆E Q = 3.020 mm/s, A = 51% for Omolon and δ = 1.158 mm/s, ∆E Q = 3.058 mm/s, A = 50% for Seymchan; (2) δ = 1.117 mm/s, ∆E Q = 2.849 mm/s, A = 43% for Omolon and δ = 1.137 mm/s, ∆E Q = 2.876 mm/s, A = 45% for Seymchan; (3) δ = 1.007 mm/s, ∆E Q = 3.108 mm/s, A = 5% for Omolon and δ = 0.952 mm/s, ∆E Q = 3.032 mm/s, A = 5% for Seymchan; (4) δ = 0.389 mm/s, ∆E Q = 0.570 mm/s, A = 1% for Omolon. Mössbauer parameters at 90 K were also similar for both olivines. It was observed that the δ values for the 57 Fe in the M1 and M2 sites in both meteorites showed different tendencies of decrease with the increase in temperature. The decrease in the isomer shift was related to the second-order Doppler shift (see details in, e.g., [41]), while the differences in the observed tendencies indicated the different Mössbauer temperatures for the 57 Fe in the M1 and M2 sites in olivines from Omolon and Seymchan, resulting from the different 57 Fe mean square velocities in these sites. Using the approach described in Part I, Section 6.7, for the ordinary chondrites [1], the values of K D and T eq were estimated for olivine extracted from Omolon PMG on the basis of XRD and Mössbauer data. For the room temperature measurements, the values of K D and T eq were 1.39 and 795 K, respectively (XRD). In the case of the room temperature Mössbauer data, these values were unrealistic until the authors considered the value of A M1 + A X instead of A M1 , supposing the possibility of the presence of some distorted M1 sites associated with the X spectral component. In the latter case, the values of K D and T eq were 1.40 and 743 K, respectively, which agree with the values calculated from XRD data. For the room temperature Mössbauer data for olivine from Seymchan PMG, the values of K D and T eq were 1.25 and 1118 K, respectively. The same calculations for the 90 K olivine Mössbauer data demonstrated K D = 1.29, T eq = 997 K for Omolon and K D = 1.31, T eq = 938 K for Seymchan. The results calculated for the 90 K Mössbauer spectra were considered the most reliable [39]. It was noteworthy that the ∆E Q values for the 57 Fe in both M1 and M2 sites were slightly different for olivines extracted from Omolon and Seymchan PMG meteorites at 295 K and at 90 K. This result indicates that the 57 Fe local microenvironments in the M1 sites and in the M2 sites in the two olivines were slightly different. Moreover, a comparison of these parameters for olivine in the studied pallasites and earlier studied ordinary chondrites Farmington L5 and Tsarev L5 showed some small differences, indicating possible variations in the 57 Fe local microenvironments in the M1 and M2 sites in olivine from different meteorites [40].  57 Fe local microenvironments in the M1 sites and in the M2 sites in the two olivines were slightly different. Moreover, a comparison of these parameters for olivine in the studied pallasites and earlier studied ordinary chondrites Farmington L5 and Tsarev L5 showed some small differences, indicating possible variations in the 57 Fe local microenvironments in the M1 and M2 sites in olivine from different meteorites [40]. Investigations of the Fe-Ni-Co alloy in the slices with a thickness of ~100 μm prepared from Krasnojarsk PMG-an ("an" means with anomalous properties), Imilac PMG, and Brenham PMG-an were done by Mössbauer spectroscopy in [42]. The room temperature Mössbauer spectrum of Imilac PMG is shown in Figure 17. This spectrum was fitted with Investigations of the Fe-Ni-Co alloy in the slices with a thickness of~100 µm prepared from Krasnojarsk PMG-an ("an" means with anomalous properties), Imilac PMG, and Brenham PMG-an were done by Mössbauer spectroscopy in [42]. The room temperature Mössbauer spectrum of Imilac PMG is shown in Figure 17. This spectrum was fitted with two magnetic sextets and one paramagnetic singlet, which were assigned to the following metal phases: kamacite, the α-Fe(Ni, Co) phase (δ = 0.02 mm/s, H eff = 337 kOe, A = 95%), tetrataenite, the γ-FeNi phase (δ = 0.04 mm/s, H eff = 286 kOe, A = 4%), and "antitaenite" (the incorrect term used by the authors of [42]; see Part I, Figure 15c,d, and further explanation in the text [1]), which is in fact the paramagnetic γ-Fe(Ni, Co) phase (δ = 0.14 mm/s, A = 1%). The same spectral components were revealed for the Mössbauer spectra of Krasnojarsk PMG-an and Brenham PMG-an: the α-Fe(Ni, Co) phase (δ = 0.02 mm/s, H eff = 337 kOe, A = 94% for Krasnojarsk and δ = 0.02 mm/s, H eff = 337 kOe, A = 97% for Brenham), the γ-FeNi phase (δ = 0.04 mm/s, H eff = 286 kOe, A = 4% for Krasnojarsk and δ = 0.02 mm/s, H eff = 280 kOe, A = 1% for Brenham), and the paramagnetic γ-Fe(Ni, Co) phase (δ = −0.05 mm/s, A = 2% for Krasnojarsk and δ = −0.07 mm/s, A = 2% for Brenham). rals 2021, 11, x FOR PEER REVIEW two magnetic sextets and one paramagnetic singlet, which were a metal phases: kamacite, the α-Fe(Ni, Co) phase (δ = 0.02 mm/s, tetrataenite, the γ-FeNi phase (δ = 0.04 mm/s, Heff = 286 kOe, A (the incorrect term used by the authors of [42]; see Part I, Figure  nation in the text [1]), which is in fact the paramagnetic γ-Fe(Ni, C A = 1%). The same spectral components were revealed for the Mö nojarsk PMG-an and Brenham PMG-an: the α-Fe(Ni, Co) phase kOe, A = 94% for Krasnojarsk and δ = 0.02 mm/s, Heff = 337 kOe the γ-FeNi phase (δ = 0.04 mm/s, Heff = 286 kOe, A = 4% for Krasn Heff = 280 kOe, A = 1% for Brenham), and the paramagnetic γ-Fe mm/s, A = 2% for Krasnojarsk and δ = −0.07 mm/s, A = 2% for Br It should be noted that the authors of [42] used a Mössbau sinusoidal velocity reference signal that increases the nonlinear The sample thickness was substantially higher than the limit of th for the Lorentzian line shape. The values of Γ for sextets were in th for kamacite and 0.30-0.34 mm/s for tetrataenite, indicating that th The unique studies of Esquel PMG were carried out using bauer spectroscopy with a spatial resolution of 10-20 μm [43,44 meteorite was polished to obtain a thin section with a thickness o ted the authors to measure the Mössbauer spectra in the local phases, their mixtures (plessite structure α-Fe(Ni, Co)/α2-Fe(Ni, It should be noted that the authors of [42] used a Mössbauer spectrometer with the sinusoidal velocity reference signal that increases the nonlinearity of the velocity scale. The sample thickness was substantially higher than the limit of thin absorber and the limit for the Lorentzian line shape. The values of Γ for sextets were in the ranges 0.33-0.37 mm/s for kamacite and 0.30-0.34 mm/s for tetrataenite, indicating that the lines were not narrow. The unique studies of Esquel PMG were carried out using 57 Fe synchrotron Mössbauer spectroscopy with a spatial resolution of 10-20 µm [43,44]. The sample of Esquel meteorite was polished to obtain a thin section with a thickness of 20-40 µm. This permitted the authors to measure the Mössbauer spectra in the local areas of different metal phases, their mixtures (plessite structure α-Fe(Ni, Co)/α 2 -Fe(Ni, Co) + γ-Fe(Ni, Co), and the cloudy zone as a transition region between kamacite and plessite), and inclusions such as schreibersite. The authors chose two line-profile transects, (1) the kamacite-cloudy zone-plessite transition (300 µm) and (2) the schreibersite grain close to kamacite/taenite interface (140 µm), for measuring 18 and 15 spectra, respectively. Some of these Mössbauer spectra are shown in Figure 18. These spectra demonstrate the presence of two magnetic sextets related to kamacite and tetrataenite, a paramagnetic singlet (the authors assigned this singlet to questionable "antitaenite", while this component should be considered as the paramagnetic γ-Fe(Ni, Co) phase), and a quadrupole doublet for schreibersite. The values of H eff for all sextets assigned to the α-Fe(Ni, Co) phase were the same within the error (~334-338 kOe), while those for all sextets related to the γ-FeNi phase were in the range~281-288 kOe. separately and investigated these separates. The room temperature Mössbauer spectra of Seymchan metallic part in comparison with the reference bcc α-Fe foil with a thickness of 7 μm are shown in Figure 19. The α-Fe foil spectrum is a symmetrical sextet with respect to the sextet center (the slightly higher intensity of the second and the fifth absorption lines than the ratio 3:2:1 for intensities of the first (sixth) to the second (fifth) and to the third (fourth) lines is the result of the texture effect after 7 μm foil preparation). In contrast, the spectrum of metallic part extracted from Seymchan PMG demonstrated an asymmetrical six-line pattern, which was more pronounced than that in the spectra of the Imilac metallic part (see Figure 17) and Esquel plessite (see Figure 18c). This six-line pattern was fitted using eight magnetic sextets and one quadrupole doublet.  A new fragment of Seymchan PMG was studied by Mössbauer spectroscopy with a high velocity resolution in [45]. The authors extracted the metallic and the stony parts separately and investigated these separates. The room temperature Mössbauer spectra of Seymchan metallic part in comparison with the reference bcc α-Fe foil with a thickness of 7 µm are shown in Figure 19. The α-Fe foil spectrum is a symmetrical sextet with respect to the sextet center (the slightly higher intensity of the second and the fifth absorption lines than the ratio 3:2:1 for intensities of the first (sixth) to the second (fifth) and to the third (fourth) lines is the result of the texture effect after 7 µm foil preparation). In contrast, the spectrum of metallic part extracted from Seymchan PMG demonstrated an asymmetrical six-line pattern, which was more pronounced than that in the spectra of the Imilac metallic part (see Figure 17) and Esquel plessite (see Figure 18c). This six-line pattern was fitted using eight magnetic sextets and one quadrupole doublet. Minerals 2021, 11, x FOR PEER REVIEW 19 of 50 (a) (b) Figure 19. Room temperature Mössbauer spectra of Seymchan PMG metallic part measured in 4096 channels and presented in 1024 channels (a) and reference α-Fe foil measured in 4096 channels (b). Indicated components are the results of the best fits. The differential spectra are shown on the bottom. Adapted from [45].
The Seymchan PMG matrix alloy consisted of α-Fe(Ni, Co) phase, γ-Fe(Ni, Co) phase, and plessite structure as shown by metallography. Therefore, the following phases were assigned to revealed magnetic sextets:  [45], was different than that in [43]. It is possible that the increase in the number of spectral points in the complex Mössbauer spectra permits decomposing these spectra with a larger number of spectral components, which can be assigned to slightly different local 57 Fe microenvironments.
The Mössbauer spectra of the stony part of Seymchan PMG were measured in the large and small velocity ranges in order to verify the presence of the magnetically split components first and then to obtain the spectrum with a better resolution in the small velocity range (see Figure 20).  The Seymchan PMG matrix alloy consisted of α-Fe(Ni, Co) phase, γ-Fe(Ni, Co) phase, and plessite structure as shown by metallography. Therefore, the following phases were assigned to revealed magnetic sextets:  [45], was different than that in [43]. It is possible that the increase in the number of spectral points in the complex Mössbauer spectra permits decomposing these spectra with a larger number of spectral components, which can be assigned to slightly different local 57 Fe microenvironments.
The Mössbauer spectra of the stony part of Seymchan PMG were measured in the large and small velocity ranges in order to verify the presence of the magnetically split components first and then to obtain the spectrum with a better resolution in the small velocity range (see Figure 20).
The stony part of Seymchan PMG did not contain Fe-Ni-Co alloy and orthopyroxene; the absence of these phases was confirmed by XRD and scanning electron microscopy with energy dispersive spectroscopy. These techniques and Mössbauer spectroscopy showed the presence of olivine as the expected main mineral (δ = 1.188 mm/s, ∆E Q = 3.024 mm/s, A =~51.5% for the M1 sites and δ = 1. 106  was different than that in [43]. It is possible that the increase in the number of spectral points in the complex Mössbauer spectra permits decomposing these spectra with a larger number of spectral components, which can be assigned to slightly different local 57 Fe microenvironments. The Mössbauer spectra of the stony part of Seymchan PMG were measured in the large and small velocity ranges in order to verify the presence of the magnetically split components first and then to obtain the spectrum with a better resolution in the small velocity range (see Figure 20).

Mesosiderites
Mesosiderites are a rare type of stony-iron meteorite, which appear as brecciated rocks and contain subequal metal (Fe-Ni-Co alloy) and silicate components (olivine, pyroxenes); the silicates are dominantly igneous rock fragments. Mesosiderites are divided into different groups according to the petrologic class (A, B, C) and metamorphic grade (1, 2, 3, etc.) (see MBD). Only two articles have been published on investigations of mesosiderites using Mössbauer spectroscopy. The authors of [42]  The stony part of Seymchan PMG did not contain Fe-Ni-Co alloy and orthopyroxene; the absence of these phases was confirmed by XRD and scanning electron microscopy with energy dispersive spectroscopy. These techniques and Mössbauer spectroscopy showed the presence of olivine as the expected main mineral (δ = 1.188 mm/s, ΔEQ = 3.024 mm/s, A = ~51.5% for the M1 sites and δ = 1.

Mesosiderites
Mesosiderites are a rare type of stony-iron meteorite, which appear as brecciated rocks and contain subequal metal (Fe-Ni-Co alloy) and silicate components (olivine, pyroxenes); the silicates are dominantly igneous rock fragments. Mesosiderites are divided into different groups according to the petrologic class (A, B, C) and metamorphic grade   [42] did not reveal components which could be assigned to the ferromagnetic γ-Fe(Ni, Co) phase (taenite).
The authors of [44] applied synchrotron Mössbauer spectroscopy to study individual metal grains in Estherville Mes-A3/4. They measured the room temperature spectra of two Fe-Ni-Co alloy grains with α-phase and with a mixture of α-phase and γ-phase (see Figure 22). The authors observed one magnetic sextet corresponding to kamacite (Figure 22a), two magnetic sextets assigned to kamacite and tetrataenite, and one paramagnetic singlet considered as "antitaenite" which is, in fact, the paramagnetic taenite ( Figure 22b). The authors of [44] applied synchrotron Mössbauer spectroscopy to study individual metal grains in Estherville Mes-A3/4. They measured the room temperature spectra of two Fe-Ni-Co alloy grains with α-phase and with a mixture of α-phase and γ-phase (see Figure  22). The authors observed one magnetic sextet corresponding to kamacite (Figure 22a), two magnetic sextets assigned to kamacite and tetrataenite, and one paramagnetic singlet considered as "antitaenite" which is, in fact, the paramagnetic taenite (Figure 22b).

Iron Meteorites
Iron meteorites consist of Fe-Ni-Co alloy with small inclusions of minerals such as troilite, schreibersite and rhabdites, daubréelite, and silicate phases. The main feature of iron meteorites is the Widmanstätten structures that can be seen on the polished slice. Iron meteorites are structurally divided in octahedrites O (fine Of, finest Off, medium Om, coarse Og, coarsest Ogg, and plessitic Opl), hexahedrites H, and ataxites D (see, e.g., [2]). Furthermore, chemical classification is used for iron meteorites, and they can be divided into 13 groups denoted using Roman numbers I-IV and one or two Latin letters A, B, C, D, E, F, and G: IAB complex iron, IC, IIAB, IIC, IID, IIE, IIF, IIG, IIIAB, IIIE, IIIF, IVA, and IVB  (see MBD). Additionally, in the cases of anomalous or ungrouped iron meteorites, in the group name, the letters "-an" or "-ung", as well as "MG" for the main group, are added.
Further study of Fe-Ni alloys in the Antarctic Ni-rich ataxite Y-791694 IAB meteorite in comparison with non-Antarctic Ni-rich ataxites Santa Catharina IAB-ung, Twin City IABung, and some other meteorites was done in [67]. The room temperature Mössbauer spectra of thin slices of Y-791694 and Twin City are shown in Figure 23. The spectrum of Y-791694 demonstrated two magnetic sextets indicating the presence of two Fe-Ni phases with the main component associated by the authors with fully disordered Fe-Ni alloy in spite of the Ni content being~40 at.%. This may be a result of the fast cooling of Y-791694. In contrast, the spectra of Santa Catharina and Twin City (see Figure 23b) demonstrated the presence of a magnetic sextet and paramagnetic singlet. The magnetic sextet corresponded to the ordered γ-FeNi phases or tetrataenite (50 at.% Fe-50 at.% Ni) with L1 0 superstructure. This indicates that the decomposition of the phases in these meteorites was started at low temperature.

Iron Meteorites
Iron meteorites consist of Fe-Ni-Co alloy with small inclusions of minerals such as troilite, schreibersite and rhabdites, daubréelite, and silicate phases. The main feature of iron meteorites is the Widmanstätten structures that can be seen on the polished slice. Iron meteorites are structurally divided in octahedrites O (fine Of, finest Off, medium Om, coarse Og, coarsest Ogg, and plessitic Opl), hexahedrites H, and ataxites D (see, e.g., [2]). Furthermore, chemical classification is used for iron meteorites, and they can be divided into 13 groups denoted using Roman numbers I-IV and one or two Latin letters A, B, C, D, E, F, and G: IAB complex iron, IC, IIAB, IIC, IID, IIE, IIF, IIG, IIIAB, IIIE, IIIF, IVA, and  IVB (see MBD). Additionally, in the cases of anomalous or ungrouped iron meteorites, in the group name, the letters "-an" or "-ung", as well as "MG" for the main group, are added.
Further study of Fe-Ni alloys in the Antarctic Ni-rich ataxite Y-791694 IAB meteorite in comparison with non-Antarctic Ni-rich ataxites Santa Catharina IAB-ung, Twin City IAB-ung, and some other meteorites was done in [67]. The room temperature Mössbauer spectra of thin slices of Y-791694 and Twin City are shown in Figure 23. The spectrum of Y-791694 demonstrated two magnetic sextets indicating the presence of two Fe-Ni phases with the main component associated by the authors with fully disordered Fe-Ni alloy in spite of the Ni content being ~40 at.%. This may be a result of the fast cooling of Y-791694. In contrast, the spectra of Santa Catharina and Twin City (see Figure 23b) demonstrated the presence of a magnetic sextet and paramagnetic singlet. The magnetic sextet corresponded to the ordered γ-FeNi phases or tetrataenite (50 at.% Fe-50 at.% Ni) with L10 superstructure. This indicates that the decomposition of the phases in these meteorites was started at low temperature.  A quite different Mössbauer spectrum was measured for Henbury IIIAB iron meteorite in [68] (see Figure 24a). This meteorite consisted of~91.9 at.% Fe and~7.5 at.% Ni. This alloy was a homogeneous bcc α-phase. Different fits of the measured spectrum showed the same 57 Fe hyperfine parameters: δ =~0.02 mm/s, H eff = 337.7 kOe. These values corresponded to the α-Fe(Ni, Co) phase (kamacite). In contrast, the Mössbauer spectrum of Santa Catharina IAB-ung measured in [69] was similar to the spectrum of Twin City (compare Figures 23b  and 24b). Santa Catharina ataxite contained the highest known Ni content (35 wt.% Ni). The Mössbauer spectrum was fitted using two sextets: one symmetrical sextet with larger H eff and broad lines, which was assigned to disordered tetrataenite (in fact this is taenite because tetrataenite is always ordered), and one asymmetrical sextet with smaller H eff and narrow lines, which was associated with ordered tetrataenite, as well as one paramagnetic singlet, which the authors of [69] denoted as "antitaenite" similar to [70] and wrote about the experimental proof of their observation. A discussion about the questionable term "antitaenite" was already reported in Part I, Section 6.1 [1]. Moreover, the authors of [64], who also observed the paramagnetic singlet in the Mössbauer spectrum of Santa Catharina, determined the Ni content of 28 at.% for this component that corresponded to the wellknown paramagnetic γ-Fe(Ni) phase (paramagnetic taenite). [70] and wrote about the experimental proof of their observation. A discussion ab questionable term "antitaenite" was already reported in Part I, Section 6.1 [1]. Mo the authors of [64], who also observed the paramagnetic singlet in the Mössbaue trum of Santa Catharina, determined the Ni content of 28 at.% for this compone corresponded to the well-known paramagnetic γ-Fe(Ni) phase (paramagnetic taen A comparison of the different areas of the coarse octahedrite Soledade IAB-Mössbauer spectroscopy was carried out in [71]. The Fe-Ni alloy contained about 6 Ni, indicating the α-Fe(Ni, Co) phase. The Mössbauer spectra of samples obtaine the interior and oxidized surface of Soledade measured at room temperature are in Figure 25. These spectra demonstrate asymmetrical six-line patterns with an add quadrupole doublet for the oxidized sample.
A similar study was done for the interior and surface of hexahedrite from the del Cielo meteorite shower in [72] (the name of this meteorite is unknown, whil were many registered iron meteorites found from this shower). This meteorite con 94.6 wt.% Fe and 5.4 wt.% Ni. The room temperature Mössbauer spectra of the tw ples are shown in Figure 26. The interior spectrum demonstrated an asymmetric pattern similar to the other meteoritic Fe-Ni-Co alloys with low Ni content. In co the surface spectrum mainly showed the presence of ferric compounds in the m and paramagnetic states.  A comparison of the different areas of the coarse octahedrite Soledade IAB-MG by Mössbauer spectroscopy was carried out in [71]. The Fe-Ni alloy contained about 6.7 at.% Ni, indicating the α-Fe(Ni, Co) phase. The Mössbauer spectra of samples obtained from the interior and oxidized surface of Soledade measured at room temperature are shown in Figure 25. These spectra demonstrate asymmetrical six-line patterns with an additional quadrupole doublet for the oxidized sample.
A similar study was done for the interior and surface of hexahedrite from the Campo del Cielo meteorite shower in [72] (the name of this meteorite is unknown, while there were many registered iron meteorites found from this shower). This meteorite contained 94.6 wt.% Fe and 5.4 wt.% Ni. The room temperature Mössbauer spectra of the two samples are shown in Figure 26. The interior spectrum demonstrated an asymmetric six-line pattern similar to the other meteoritic Fe-Ni-Co alloys with low Ni content. In contrast, the surface spectrum mainly showed the presence of ferric compounds in the magnetic and paramagnetic states. The study of Sikhote-Alin IIAB, Bilibino IIAB, Chinga and Dronino (both are u grouped iron meteorites), and products of their weathering was carried out by Mössbau spectroscopy in [73][74][75]. The measured Mössbauer spectra demonstrated asymmetri six-line patterns which were fitted with different numbers of magnetic sextets. The spec of the corrosion (weathered) products showed the presence of various ferric oxides a oxyhydroxides depending on the place of meteorite weathering. These Mössbauer spec were measured using an SM-2201 spectrometer with a high velocity resolution but fit in the presentation in 512-channel spectra only. Furthermore, these meteorites were measured to obtain the high velocity resolution Mössbauer spectra. A comparison of Mössbauer spectra of Chinga iron-ung foil presented in 512 and 4096 channels in [76 shown in Figure 27. The first spectrum was decomposed using three magnetic sextets fit an asymmetrical six-line pattern, while the second spectrum was fitted with six m netic sextets. The values of Heff for these six sextets were in the range ~354-305 kOe. T Mössbauer spectra of Chinga metal fine shaving and metal powder were also fitted us six magnetic sextets with Heff in the ranges ~343-296 kOe and ~354-295 kOe, respective These differences can be attributed to the effect of sample preparation on the phase co position. The study of Sikhote-Alin IIAB, Bilibino IIAB, Chinga and Dronino (both are u grouped iron meteorites), and products of their weathering was carried out by Mössbau spectroscopy in [73][74][75]. The measured Mössbauer spectra demonstrated asymmetri six-line patterns which were fitted with different numbers of magnetic sextets. The spec of the corrosion (weathered) products showed the presence of various ferric oxides a oxyhydroxides depending on the place of meteorite weathering. These Mössbauer spec were measured using an SM-2201 spectrometer with a high velocity resolution but fitt in the presentation in 512-channel spectra only. Furthermore, these meteorites were measured to obtain the high velocity resolution Mössbauer spectra. A comparison of t Mössbauer spectra of Chinga iron-ung foil presented in 512 and 4096 channels in [76] shown in Figure 27. The first spectrum was decomposed using three magnetic sextets fit an asymmetrical six-line pattern, while the second spectrum was fitted with six ma netic sextets. The values of Heff for these six sextets were in the range ~354-305 kOe. T Mössbauer spectra of Chinga metal fine shaving and metal powder were also fitted usi six magnetic sextets with Heff in the ranges ~343-296 kOe and ~354-295 kOe, respective These differences can be attributed to the effect of sample preparation on the phase co position. The study of Sikhote-Alin IIAB, Bilibino IIAB, Chinga and Dronino (both are ungrouped iron meteorites), and products of their weathering was carried out by Mössbauer spectroscopy in [73][74][75]. The measured Mössbauer spectra demonstrated asymmetrical six-line patterns which were fitted with different numbers of magnetic sextets. The spectra of the corrosion (weathered) products showed the presence of various ferric oxides and oxyhydroxides depending on the place of meteorite weathering. These Mössbauer spectra were measured using an SM-2201 spectrometer with a high velocity resolution but fitted in the presentation in 512-channel spectra only. Furthermore, these meteorites were re-measured to obtain the high velocity resolution Mössbauer spectra. A comparison of the Mössbauer spectra of Chinga iron-ung foil presented in 512 and 4096 channels in [76] is shown in Figure 27. The first spectrum was decomposed using three magnetic sextets to fit an asymmetrical six-line pattern, while the second spectrum was fitted with six magnetic sextets. The values of H eff for these six sextets were in the range~354-305 kOe. The Mössbauer spectra of Chinga metal fine shaving and metal powder were also fitted using six magnetic sextets with H eff in the ranges~343-296 kOe and~354-295 kOe, respectively. These differences can be attributed to the effect of sample preparation on the phase composition. An investigation of an iron meteorite found in Colombia named "Gaspar" was carried out by Mössbauer spectroscopy in [77]. The name "Gaspar" was not registered, but this fragment was considered as a part of the Santa Rosa IC meteorite. The authors measured the Mössbauer spectra of two samples from different sides of "Gaspar" (fragments A and B), which are shown in Figure 28. Although these spectra were asymmetrical sixline patterns, they were fitted using one magnetic sextet only with the same value of Heff = 337 Oe. However, these fits cannot be considered as final and adequate because the enveloped lines could not fit the experimental points well, while the authors did not use differential spectra for quality control. Samples of Sikhote-Alin IIAB were studied by Mössbauer spectroscopy in [78] and [79]. The measured room temperature Mössbauer spectra in both works are shown in Figure 29 for comparison. Both spectra showed similar asymmetric six-line patterns which were fitted in different ways. The authors of [78] used a superposition of multiple magnetic sextets with different Heff and hyperfine field distribution as a result of different local 57 Fe microenvironments due to a distribution of Ni and Co atoms in the first coordination sphere. The average value of Heff was 335.7 kOe. On the other hand, the authors of [79] fitted their spectrum using three magnetic sextets and one paramagnetic singlet. The 57 Fe hyperfine parameters for these sextets were δ = 0.03 mm/s, Heff = 346.5 kOe, A = 25% (sextet 1), δ = 0.02 mm/s, Heff = 336.5 kOe, A = 50% (sextet 2), δ = 0.03 mm/s, Heff = 326.7 kOe, A = Figure 27. Room temperature Mössbauer spectra of Chinga iron-ung meteorite: the spectrum presented in 512 channels (a) and the spectrum presented in 4096 channels (b). Indicated components are the results of the best fits. The differential spectrum is shown on the bottom in (b). Adapted from [76].
An investigation of an iron meteorite found in Colombia named "Gaspar" was carried out by Mössbauer spectroscopy in [77]. The name "Gaspar" was not registered, but this fragment was considered as a part of the Santa Rosa IC meteorite. The authors measured the Mössbauer spectra of two samples from different sides of "Gaspar" (fragments A and B), which are shown in Figure 28. Although these spectra were asymmetrical six-line patterns, they were fitted using one magnetic sextet only with the same value of H eff = 337 Oe. However, these fits cannot be considered as final and adequate because the enveloped lines could not fit the experimental points well, while the authors did not use differential spectra for quality control. An investigation of an iron meteorite found in Colombia named "Gaspar" was ca ried out by Mössbauer spectroscopy in [77]. The name "Gaspar" was not registered, b this fragment was considered as a part of the Santa Rosa IC meteorite. The authors mea ured the Mössbauer spectra of two samples from different sides of "Gaspar" (fragmen A and B), which are shown in Figure 28. Although these spectra were asymmetrical si line patterns, they were fitted using one magnetic sextet only with the same value of H = 337 Oe. However, these fits cannot be considered as final and adequate because the e veloped lines could not fit the experimental points well, while the authors did not u differential spectra for quality control. Samples of Sikhote-Alin IIAB were studied by Mössbauer spectroscopy in [78] an [79]. The measured room temperature Mössbauer spectra in both works are shown in Fi ure 29 for comparison. Both spectra showed similar asymmetric six-line patterns whic were fitted in different ways. The authors of [78] used a superposition of multiple ma netic sextets with different Heff and hyperfine field distribution as a result of different loc 57 Fe microenvironments due to a distribution of Ni and Co atoms in the first coordinatio sphere. The average value of Heff was 335.7 kOe. On the other hand, the authors of [7 fitted their spectrum using three magnetic sextets and one paramagnetic singlet. The 57 F hyperfine parameters for these sextets were δ = 0.03 mm/s, Heff = 346.5 kOe, A = 25% (sext 1), δ = 0.02 mm/s, Heff = 336.5 kOe, A = 50% (sextet 2), δ = 0.03 mm/s, Heff = 326.7 kOe, A Samples of Sikhote-Alin IIAB were studied by Mössbauer spectroscopy in [78,79]. The measured room temperature Mössbauer spectra in both works are shown in Figure 29 for comparison. Both spectra showed similar asymmetric six-line patterns which were fitted in different ways. The authors of [78] used a superposition of multiple magnetic sextets with different H eff and hyperfine field distribution as a result of different local 57 Fe microenvironments due to a distribution of Ni and Co atoms in the first coordination sphere. The average value of H eff was 335.7 kOe. On the other hand, the authors of [79] fitted their spectrum using three magnetic sextets and one paramagnetic singlet. The 57 Fe hyperfine parameters for these sextets were δ = 0.03 mm/s, H eff = 346.5 kOe, A = 25% (sextet 1), δ = 0.02 mm/s, H eff = 336.5 kOe, A = 50% (sextet 2), δ = 0.03 mm/s, H eff = 326.7 kOe, A = 23% (sextet 3), and δ = 0 mm/s, A = 2% (singlet). These sextets were assigned to α-Fe(Ni, Co) with different Ni concentrations in the local iron microenvironment. (a) (b) Figure 29. Comparison of the room temperature Mössbauer spectra of Sikhote-Alin IIAB measured in [78] (a) and in [79] (b). Indicated spectral components are the results of the fits. The differential spectra are shown on the bottom. Adapted from [78,79].
The Ni contents in Gebel Kamil and NWA 6259 were 19.8 wt.% and 42.2 wt.%, respectively (MBD). Therefore, the authors of [42] revealed 85% of α-Fe(Ni, Co) phase, 12% of γ-Fe(Ni, Co) phase, and 3% of paramagnetic γ-Fe(Ni, Co) phase (which the authors called "antitaenite") in the Mössbauer spectrum of Gebel Kamil. In contrast, these authors obtained 12% of γ-Fe(Ni, Co) phase, 74% of γ-FeNi phase (tetrataenite), and 14% of paramagnetic γ-Fe(Ni, Co) phase for the Mössbauer spectrum of NWA 6259. Toluca and São João Nepomuceno meteorites contained~8 wt.% Ni. Therefore, their Mössbauer spectra were fitted using one magnetic sextet assigned to the α-Fe(Ni, Co) phase. However, the Mössbauer spectra of Ni-enriched metal from both meteorites demonstrated the presence of 72% of γ-FeNi phase and 28% of paramagnetic γ-Fe(Ni, Co) phase for Toluca and 16% of α-Fe(Ni, Co) phase, 60% of γ-Fe(Ni, Co) phase, and 24% of paramagnetic γ-Fe(Ni, Co) phase for São João Nepomuceno. For the studied iron meteorites, except San Cristobal, the authors of [42] obtained the following ranges of H eff : 331-339 kOe for α-Fe(Ni, Co) phase, 297-312 kOe for γ-Fe(Ni, Co) phase, and 284-295 kOe for γ-FeNi phase. As for San Cristobal, this meteorite contained 25.2 wt.% Ni. Therefore, the value of H eff = 349 kOe may be assigned to the α 2 -Fe(Ni, Co) phase which contained Ni in the range of 8-25 at.%, while H eff = 333 kOe should be related to the α-Fe(Ni, Co) phase (the authors of [42] assigned these values of H eff to kamacite and taenite, respectively, which seems to be incorrect). The δ values for the paramagnetic γ-Fe(Ni, Co) were in the range −0.08 to +0.02 mm/s. The comparison of powdered samples of one fragment of Chinga iron-ung, in which visually different areas (1, 2, 3, 4, and 5) of the saw-cut surface were observed, was done using Mössbauer spectroscopy with a high velocity resolution (4096 channels) in [80]. The Mössbauer spectra of these samples obtained from four areas are shown in Figure 31. These spectra are six-line patterns with asymmetry, unlike that observed for Chinga foil (Figure 27b) in [76]. The new spectra were measured with larger signal-to-noise ratios and decomposed using nine magnetic sextets. The obtained distributions of the magnetic hyperfine fields shown in the histograms in Figure 32 demonstrated some differences that could be related to structural variations in different areas of the studied fragment of Chinga. Visual differences in the absorption line shape of the most negative peaks of sextets are also shown in Figure 31. These differences can be seen due to a large number of spectral points in the 4096-channel Mössbauer spectra.
A comparative study of iron meteorites Anyujskij IIAB, Sikhote-Alin IIAB, Aliskerovo IIIE-an, and Sterlitamak IIIAB was carried out using Mössbauer spectroscopy with a high velocity resolution in [81,82]. The 1024-channel Mössbauer spectra are shown in Figure 33. These spectra demonstrate similar asymmetric six-line patterns. However, these spectra were decomposed with different numbers of components; three magnetic sextets were revealed in the spectra of IIAB iron meteorites, while five magnetic sextets and one paramagnetic singlet were revealed in the spectra of IIIE-an and IIIAB iron meteorites. This may be related to different Ni contents and phase compositions in these meteorites: Anyujskij IIAB and Sikhote-Alin IIAB consisted of an α-Fe(Ni, Co) phase with ~5 at.% Ni, while Aliskerovo IIIE-an contained an α-Fe(Ni, Co) phase with ~4-7 at.% Ni, as well as probably an α2-Fe(Ni, Co) phase with ~7-9 at.% Ni, γ-Fe(Ni, Co) phase with ~32-35 at.% Ni, and γ-FeNi phase with 50 at.% Ni, and Sterlitamak IIIAB consisted of an α-Fe(Ni, Co) The comparison of powdered samples of one fragment of Chinga iron-ung, in which visually different areas (1, 2, 3, 4, and 5) of the saw-cut surface were observed, was done using Mössbauer spectroscopy with a high velocity resolution (4096 channels) in [80]. The Mössbauer spectra of these samples obtained from four areas are shown in Figure 31. These spectra are six-line patterns with asymmetry, unlike that observed for Chinga foil (Figure 27b) in [76]. The new spectra were measured with larger signal-to-noise ratios and decomposed using nine magnetic sextets. The obtained distributions of the magnetic hyperfine fields shown in the histograms in Figure 32 demonstrated some differences that could be related to structural variations in different areas of the studied fragment of Chinga. Visual differences in the absorption line shape of the most negative peaks of sextets are also shown in Figure 31. These differences can be seen due to a large number of spectral points in the 4096-channel Mössbauer spectra.
A comparative study of iron meteorites Anyujskij IIAB, Sikhote-Alin IIAB, Aliskerovo IIIE-an, and Sterlitamak IIIAB was carried out using Mössbauer spectroscopy with a high velocity resolution in [81,82]. The 1024-channel Mössbauer spectra are shown in Figure 33. These spectra demonstrate similar asymmetric six-line patterns. However, these spectra were decomposed with different numbers of components; three magnetic sextets were revealed in the spectra of IIAB iron meteorites, while five magnetic sextets and one paramagnetic singlet were revealed in the spectra of IIIE-an and IIIAB iron meteorites. This may be related to different Ni contents and phase compositions in these meteorites: Anyujskij IIAB and Sikhote-Alin IIAB consisted of an α-Fe(Ni, Co) phase with~5 at.% Ni, while Aliskerovo IIIE-an contained an α-Fe(Ni, Co) phase with~4-7 at.% Ni, as well as probably an α 2 -Fe(Ni, Co) phase with~7-9 at.% Ni, γ-Fe(Ni, Co) phase with~32-35 at.% Ni, and γ-FeNi phase with 50 at.% Ni, and Sterlitamak IIIAB consisted of an α-Fe(Ni, Co) phase with~6.5-7 at.% Ni, as well as probably an α 2 -Fe(Ni, Co) phase with~7-8 at.% Ni and γ-Fe(Ni, Co) phase with~30-34 at.% Ni.
inerals 2021, 11, x FOR PEER REVIEW 28 of 50 Figure 31. Comparison of the room temperature Mössbauer spectra of different powdered samples of the Chinga ironung meteorite obtained from areas 1, 2, 3, and 5; the most negative peaks of these spectra are enlarged. The differential spectra are shown on the bottom. Adapted from [80].  Figure 31. Adapted from Ref. [80].
Therefore, the 57 Fe hyperfine parameters for revealed magnetic sextets were assigned in [81,82] 1, 2, 3, and 5; the most negative peaks of these spectra are enlarged. The differential spectra are shown on the bottom. Adapted from [80].
Minerals 2021, 11, x FOR PEER REVIEW 28 of Figure 31. Comparison of the room temperature Mössbauer spectra of different powdered samples of the Chinga ironung meteorite obtained from areas 1, 2, 3, and 5; the most negative peaks of these spectra are enlarged. The differential spectra are shown on the bottom. Adapted from [80].  Figure 31. Adapted from Ref. [80].
Synchrotron Mössbauer spectroscopy was applied to study the Tazewell IAB-sLH iron meteorite with a spatial resolution of~10-20 µm in [43,44]. Therefore, the authors were able to measure the spectra of selected local areas with different phase composition. Examples of the room temperature Mössbauer spectra of the cloudy zone and plessite area of Tazewell are shown in Figure 35. The spectrum of the cloudy zone demonstrated the presence of two magnetic sextets related to kamacite (with a larger relative area) and tetrataenite (with a smaller relative area) and a paramagnetic singlet associated by the authors with "antitaenite" (in fact, the paramagnetic γ-Fe(Ni, Co) phase). On the other hand, the spectrum of plessite area surprisingly consisted of only two components: the main magnetic sextet assigned to kamacite and a small paramagnetic singlet which should be associated with the paramagnetic taenite (because the plessite structure is a mixture of α-Fe(Ni, Co)/α 2 -Fe(Ni, Co) + γ-Fe(Ni, Co) phases). The values of H eff for kamacite were in the range 335.4-339.9 kOe, while those for tetrataenite varied between 284 and 289 kOe. The peak positions δ of the paramagnetic taenite were in the range between −0.004 and −0.08 mm/s. three magnetic sextets and one paramagnetic singlet. The Ni concentration in Muonionalusta IVA iron meteorite was determined to be about 7.9%. Therefore, magnetic sextets with parameters (i) δ = 0.04 mm/s, Heff = 346 kOe, A = 34%, (ii) δ = 0.01 mm/s, Heff = 337 kOe, A = 40%, and (iii) δ = 0.03 mm/s, Heff = 328 kOe, A = 24% were related to the α-Fe(Ni, Co) phase with variations of Ni in the iron local microenvironment. A small singlet with parameters δ = 0.02 mm/s and A = 2% was assigned to the paramagnetic γ-Fe(Ni, Co) phase.
(a) (b) Figure 35. Room temperature synchrotron Mössbauer spectra of the Tazewell IAB-sLH iron meteorite selected areas: cloudy zone (a) and plessite area (b). Indicated spectral components are the results of the fits. The differential spectra are shown on the top. Adapted from [43]. It is well known that iron meteorites contain various iron-bearing inclusions such as troilite, schreibersite, and some other compounds, in addition to the α-and γ-phases. Therefore, several papers were dedicated to the investigation of the extracted or separated compounds observed in iron meteorites. The ferromagnetic phases in graphite nodule from Canyon Diablo IAB-MG were studied by Mössbauer spectroscopy in [86]. The authors measured two Mössbauer spectra of the extracted graphite nodule material at room temperature and at 15 K (see Figure 37). The first spectrum showed the presence of three magnetic sextets and two quadrupole doublets. The magnetic sextets marked MA and MB with Heff values of 490 kOe and 458 kOe, respectively, were related to the (A) and [B] sites in magnetite while the sextet K with δ = 0.03 mm/s and Heff = 335 kOe was assigned to Figure 35. Room temperature synchrotron Mössbauer spectra of the Tazewell IAB-sLH iron meteorite selected areas: cloudy zone (a) and plessite area (b). Indicated spectral components are the results of the fits. The differential spectra are shown on the top. Adapted from [43].
A sample of Muonionalusta IVA iron meteorite was studied by conversion electron Mössbauer spectroscopy in [85]. The room temperature CEMS spectrum is shown in Figure 36. The authors of [85] decomposed their CEMS spectrum using four components: three magnetic sextets and one paramagnetic singlet. The Ni concentration in Muonionalusta IVA iron meteorite was determined to be about 7.9%. Therefore, magnetic sextets with parameters (i) δ = 0.04 mm/s, H eff = 346 kOe, A = 34%, (ii) δ = 0.01 mm/s, H eff = 337 kOe, A = 40%, and (iii) δ = 0.03 mm/s, H eff = 328 kOe, A = 24% were related to the α-Fe(Ni, Co) phase with variations of Ni in the iron local microenvironment. A small singlet with parameters δ = 0.02 mm/s and A = 2% was assigned to the paramagnetic γ-Fe(Ni, Co) phase.
It is well known that iron meteorites contain various iron-bearing inclusions such as troilite, schreibersite, and some other compounds, in addition to the αand γ-phases. Therefore, several papers were dedicated to the investigation of the extracted or separated compounds observed in iron meteorites. The ferromagnetic phases in graphite nodule from Canyon Diablo IAB-MG were studied by Mössbauer spectroscopy in [86]. The authors measured two Mössbauer spectra of the extracted graphite nodule material at room temperature and at 15 K (see Figure 37). The first spectrum showed the presence of three magnetic sextets and two quadrupole doublets. The magnetic sextets marked M A and M B with H eff values of 490 kOe and 458 kOe, respectively, were related to the (A) and [B] sites in magnetite while the sextet K with δ = 0.03 mm/s and H eff = 335 kOe was assigned to kamacite. One quadrupole doublet marked A with δ = 0.36 mm/s and ∆E Q = 0.69 mm/s was associated with akaganéite, while the other quadrupole doublet E with δ = 1.17 mm/s and ∆E Q = 2.00 mm/s was related to iron in the M2 site of enstatite. The 15 K Mössbauer spectrum consisted of magnetically split components with H eff > 300 kOe.
(a) (b) Figure 35. Room temperature synchrotron Mössbauer spectra of the Tazewell IAB-sLH iron meteorite selected areas: cloudy zone (a) and plessite area (b). Indicated spectral components are the results of the fits. The differential spectra are shown on the top. Adapted from [43]. It is well known that iron meteorites contain various iron-bearing inclusions such as troilite, schreibersite, and some other compounds, in addition to the α-and γ-phases. Therefore, several papers were dedicated to the investigation of the extracted or separated compounds observed in iron meteorites. The ferromagnetic phases in graphite nodule from Canyon Diablo IAB-MG were studied by Mössbauer spectroscopy in [86]. The authors measured two Mössbauer spectra of the extracted graphite nodule material at room temperature and at 15 K (see Figure 37). The first spectrum showed the presence of three magnetic sextets and two quadrupole doublets. The magnetic sextets marked MA and MB with Heff values of 490 kOe and 458 kOe, respectively, were related to the (A) and [B] sites in magnetite while the sextet K with δ = 0.03 mm/s and Heff = 335 kOe was assigned to kamacite. One quadrupole doublet marked A with δ = 0.36 mm/s and ΔEQ = 0.69 mm/s was associated with akaganéite, while the other quadrupole doublet E with δ = 1.17 mm/s and ΔEQ = 2.00 mm/s was related to iron in the M2 site of enstatite. The 15 K Mössbauer spectrum consisted of magnetically split components with Heff > 300 kOe. Troilite extracted from the iron meteorites Cape York IIIAB, Nantan IAB-MG, and Sikhote-Alin IIAB was studied by Mössbauer spectroscopy in [87][88][89][90][91]. The Mössbauer spectra of troilite extracted from Cape York IIIAB (called as Agpalilik) were measured in [87] at various temperatures from 77 K up to 645 K (in the temperature range 400-500 K troilite has two structural transitions). The spectrum at 645 K represented one quadrupole doublet, while the spectra at 499 K and 409 K were a six-line pattern and a superposition of three magnetic sextets, respectively. These magnetic sextets were fitted using the full static Hamiltonian. Other studies of troilite extracted from Cape York IIIAB and from Nantan IAB-MG by Mössbauer spectroscopy were carried out at very low temperatures with external magnetic field application [88][89][90]. These spectra were fitted without the full static Hamiltonian. In contrast, the Mössbauer spectra of troilite extracted from Sikhote-Alin IIAB measured with a high velocity resolution at 295 K and 90 K were fitted using the full static Hamiltonian [91]. These authors also revealed the spectral components assigned to nonstoichiometric troilite Fe1−xS (which were fitted using the hyperfine field distribution) and a paramagnetic singlet with δ = ~0.65 mm/s at 295 K, which was related to daubréelite. The authors of [92] studied samples of troilite, kamacite, and taenite extracted from the Morasko IAB-MG iron meteorite by Mössbauer spectroscopy. They measured the room temperature spectra of these samples; however, the results obtained were partly questionable. For example, in the spectrum of extracted troilite, the authors revealed a magnetic sextet with parameters δ = 0.00 mm/s and Heff = 210 kOe, which was assigned to daubréelite. However, the room temperature Mössbauer spectrum of daubréelite is a paramagnetic singlet (see, e.g., [93], in which the FeCr2S4 Mössbauer spectrum has a singlet shape even at 170 K, and [91]). Therefore, this sextet could not be related to daubréelite, while its value of Heff was the same as that obtained by the authors for cohenite in the sample of extracted taenite (e.g., Heff = 199 kOe was assigned to cohenite in [7]). Moreover, the authors of [92] assigned the magnetic sextet with Heff = 340 kOe to disordered taenite in the spectrum of extracted taenite, which was also incorrect, because the values of Heff larger than ~329 kOe should be related to α-Fe(Ni, Co) and α2-Fe(Ni, Co) phases, as shown above. Thus, the authors of [92] did not correctly identify the mineralogical phases.
The separates of (Fe, Ni)3P inclusions from Sikhote-Alin IIAB in the forms of massive schreibersite and rhabdite microcrystals were studied by Mössbauer spectroscopy in [94][95][96][97]. These iron nickel phosphides had a crystal structure with three non-equivalent sites M1, M2, and M3 for the metal atoms, which presented different probabilities of occupation by Fe and Ni, as well as differences in these probabilities between schreibersite and Troilite extracted from the iron meteorites Cape York IIIAB, Nantan IAB-MG, and Sikhote-Alin IIAB was studied by Mössbauer spectroscopy in [87][88][89][90][91]. The Mössbauer spectra of troilite extracted from Cape York IIIAB (called as Agpalilik) were measured in [87] at various temperatures from 77 K up to 645 K (in the temperature range 400-500 K troilite has two structural transitions). The spectrum at 645 K represented one quadrupole doublet, while the spectra at 499 K and 409 K were a six-line pattern and a superposition of three magnetic sextets, respectively. These magnetic sextets were fitted using the full static Hamiltonian. Other studies of troilite extracted from Cape York IIIAB and from Nantan IAB-MG by Mössbauer spectroscopy were carried out at very low temperatures with external magnetic field application [88][89][90]. These spectra were fitted without the full static Hamiltonian. In contrast, the Mössbauer spectra of troilite extracted from Sikhote-Alin IIAB measured with a high velocity resolution at 295 K and 90 K were fitted using the full static Hamiltonian [91]. These authors also revealed the spectral components assigned to nonstoichiometric troilite Fe 1−x S (which were fitted using the hyperfine field distribution) and a paramagnetic singlet with δ =~0.65 mm/s at 295 K, which was related to daubréelite.
The authors of [92] studied samples of troilite, kamacite, and taenite extracted from the Morasko IAB-MG iron meteorite by Mössbauer spectroscopy. They measured the room temperature spectra of these samples; however, the results obtained were partly questionable. For example, in the spectrum of extracted troilite, the authors revealed a magnetic sextet with parameters δ = 0.00 mm/s and H eff = 210 kOe, which was assigned to daubréelite. However, the room temperature Mössbauer spectrum of daubréelite is a paramagnetic singlet (see, e.g., [93], in which the FeCr 2 S 4 Mössbauer spectrum has a singlet shape even at 170 K, and [91]). Therefore, this sextet could not be related to daubréelite, while its value of H eff was the same as that obtained by the authors for cohenite in the sample of extracted taenite (e.g., H eff = 199 kOe was assigned to cohenite in [7]). Moreover, the authors of [92] assigned the magnetic sextet with H eff = 340 kOe to disordered taenite in the spectrum of extracted taenite, which was also incorrect, because the values of H eff larger than~329 kOe should be related to α-Fe(Ni, Co) and α 2 -Fe(Ni, Co) phases, as shown above. Thus, the authors of [92] did not correctly identify the mineralogical phases.
The separates of (Fe, Ni) 3 P inclusions from Sikhote-Alin IIAB in the forms of massive schreibersite and rhabdite microcrystals were studied by Mössbauer spectroscopy in [94][95][96][97]. These iron nickel phosphides had a crystal structure with three non-equivalent sites M1, M2, and M3 for the metal atoms, which presented different probabilities of occupation by Fe and Ni, as well as differences in these probabilities between schreibersite and rhabdite. The Mössbauer spectra of these iron nickel phosphides extracted from Sikhote-Alin IIAB iron meteorite and measured with a high velocity resolution at 295 K and at 90 K are shown in Figure 38. rhabdite. The Mössbauer spectra of these iron nickel phosphides extracted from Sikhote-Alin IIAB iron meteorite and measured with a high velocity resolution at 295 K and at 90 K are shown in Figure 38.
The room temperature spectra of schreibersite and rhabdite samples demonstrated a clear difference due to the superparamagnetic state of rhabdite microcrystals, as the spectra of rhabdite microcrystals indicated magnetically split components at lower temperatures. The magnetically split Mössbauer spectra of schreibersite and rhabdite microcrystals were better fitted using six magnetic sextets and one paramagnetic doublet. These six sextets were related to the 57 Fe nuclei in the M1, M2, and M3 sites in iron nickel phosphides, with two sextets per site. Using this fit, the authors of [94][95][96][97] estimated the Fe and Ni occupancies of the M1, M2, and M3 sites in schreibersite and rhabdite from Sikhote-Alin iron meteorites, which appeared to be different for the two phosphides. São João Nepomuceno IVA-an is an anomalous iron meteorite containing silicate inclusions. The authors of [98,99] studied the extracted orthopyroxene using Mössbauer spectroscopy in order to analyze the thermal history of meteorite by determining the Fe 2+ populations in the M1 and M2 sites in orthopyroxene. The room temperature spectrum of extracted orthopyroxene is shown in Figure 39. This spectrum demonstrated a superposition of two quadrupole doublets related to the 57  The differential spectra are shown on the bottom. Adapted from [96].
The room temperature spectra of schreibersite and rhabdite samples demonstrated a clear difference due to the superparamagnetic state of rhabdite microcrystals, as the spectra of rhabdite microcrystals indicated magnetically split components at lower temperatures. The magnetically split Mössbauer spectra of schreibersite and rhabdite microcrystals were better fitted using six magnetic sextets and one paramagnetic doublet. These six sextets were related to the 57 Fe nuclei in the M1, M2, and M3 sites in iron nickel phosphides, with two sextets per site. Using this fit, the authors of [94][95][96][97] estimated the Fe and Ni occupancies of the M1, M2, and M3 sites in schreibersite and rhabdite from Sikhote-Alin iron meteorites, which appeared to be different for the two phosphides.
São João Nepomuceno IVA-an is an anomalous iron meteorite containing silicate inclusions. The authors of [98,99] studied the extracted orthopyroxene using Mössbauer spectroscopy in order to analyze the thermal history of meteorite by determining the Fe 2+ populations in the M1 and M2 sites in orthopyroxene. The room temperature spectrum of extracted orthopyroxene is shown in Figure 39. This spectrum demonstrated a superposition of two quadrupole doublets related to the 57

Lunar Soils and Rocks
Investigation of Lunar material became possible due to successful return missions to the Moon in 1969-1976 by two space programs: Apollo (USA) and Luna (the former USSR). Apollo missions were crewed, and the first two men appeared on the Moon on July 21, 1969 (Apollo 11). Astronauts collected Lunar soils and rocks from different parts of the Moon near the places of six Apollo mission landings, with a total weight of 380 kg, and delivered Lunar matter to the Earth. Luna space missions to the Moon were uncrewed and equipped with return modules which delivered Lunar material to the Earth with a total weight of about 320 g. Therefore, intensive studies of delivered Lunar material by various techniques including Mössbauer spectroscopy were started in 1969 by various research groups (some of these studies were reviewed in [100]).
The backscattered Mössbauer spectra of Lunar gabbroic rock and vesicular basalt measured at room temperature are shown in Figure 41. These spectra also demonstrated a large contribution of ilmenite (two peaks in the middle of the spectrum with a small quadrupole doublet), in addition to silicates and probably glasses [102]. However, the relative part of ilmenite in gabbroic rock was smaller than that in vesicular basalt.
The transmission Mössbauer spectra of two rock samples delivered by the Apollo 12

Lunar Soils and Rocks
Investigation of Lunar material became possible due to successful return missions to the Moon in 1969-1976 by two space programs: Apollo (USA) and Luna (the former USSR). Apollo missions were crewed, and the first two men appeared on the Moon on July 21, 1969 (Apollo 11). Astronauts collected Lunar soils and rocks from different parts of the Moon near the places of six Apollo mission landings, with a total weight of 380 kg, and delivered Lunar matter to the Earth. Luna space missions to the Moon were uncrewed and equipped with return modules which delivered Lunar material to the Earth with a total weight of about 320 g. Therefore, intensive studies of delivered Lunar material by various techniques including Mössbauer spectroscopy were started in 1969 by various research groups (some of these studies were reviewed in [100]).
The backscattered Mössbauer spectra of Lunar gabbroic rock and vesicular basalt measured at room temperature are shown in Figure 41. These spectra also demonstrated a large contribution of ilmenite (two peaks in the middle of the spectrum with a small quadrupole doublet), in addition to silicates and probably glasses [102]. However, the relative part of ilmenite in gabbroic rock was smaller than that in vesicular basalt.  The transmission Mössbauer spectra of two rock samples delivered by the Apollo 12 mission are shown in Figure 42. These spectra showed that the components in both rocks were similar: silicate phases and ilmenite, metallic iron alloy, and troilite. However, the relative contents of these phases were different for specimens 12002 and 12063 [105]. For example, ilmenite contents were~7% (12002) and~14% (12063), while olivine contents were~25% (12002) and~14% (12063). On the other hand, the contents of metallic iron alloy and of troilite are very close, respectively, in specimens 12002 and 12063, while pyroxene content differed by about 5% only. The authors compared the Mössbauer spectra of several rocks from the Apollo 12 mission and observed some similarities and differences. Moreover, these spectra were compared with those of rocks from the Apollo 11 mission, and the authors found some differences.
Some other investigations of Lunar soil and rocks delivered by the Apollo 11 mission using Mössbauer spectroscopy were done in [106,107]. The authors of [106] studied Lunar fines (specimen 10084,89) and light and heavy dust separated from the fines. The room temperature Mössbauer spectra of these samples are shown in Figure 43. All these spectra demonstrated the presence of large contents of pyroxene (components related to the M1 and M2 sites were identified) and ilmenite (components marked 1, 2, and 3 in Figure 43). A component corresponding to a metallic iron alloy was found in fines and light dust, while troilite was found in the spectrum of light dust only. The authors explained the absence of metallic iron alloy in the heavy dust spectrum as being due to removal during sample preparation. , and heavy dust from fines (c). Fe is metallic iron alloy, FeS is troilite, 1 represents the M1 sites in pyroxene, 2 represents the M2 sites in pyroxene, and 3 is ilmenite. Adapted from [106].
Furthermore, the large number of Lunar samples delivered by Apollo missions 11,12,14,15,16, and 17 were studied by Mössbauer spectroscopy in order to determine the iron-bearing phases and ratios of their contents in [107]. The room temperature Mössbauer spectra of basalt, specimen 70017, and poikilitic rock, specimen 60315, are shown in Figure  44. These spectra demonstrated that the basalt sample contained a small amount of troilite and large amounts of pyroxene (components related to the M1 and M2 sites) and ilmenite, while poikilitic rock additionally contained metallic iron alloy and olivine. However, the content of ilmenite was substantially smaller in the latter sample.
The results of various other studies of Lunar soils, rocks, fines, glasses, etc. from , and heavy dust from fines (c). Fe is metallic iron alloy, FeS is troilite, 1 represents the M1 sites in pyroxene, 2 represents the M2 sites in pyroxene, and 3 is ilmenite. Adapted from [106].
The authors of [106] measured Mössbauer spectra of their various Lunar samples at low temperatures down to liquid helium temperature. They deduced similar 57 Fe hyperfine parameters for pyroxene and for ilmenite in the studied sample at both 295 K and 77 K. The room temperature parameters were as follows: δ = 1.14 mm/s, ∆E Q = 2.64 mm/s (the M1 sites in pyroxene), δ = 1.10 mm/s, ∆E Q = 1.96 mm/s (the M2 sites in pyroxene), and δ = 1.07 mm/s, ∆E Q = 0.74 mm/s (ilmenite) for Lunar fines and δ = 1.11 mm/s, ∆E Q = 2.58 mm/s (the M1 sites in pyroxene), δ = 1.08 mm/s, ∆E Q = 1.98 mm/s (the M2 sites in pyroxene), and δ = 1.02 mm/s, ∆E Q = 0.72 mm/s (ilmenite) for Lunar rocks.
Furthermore, the large number of Lunar samples delivered by Apollo missions 11,12,14,15,16, and 17 were studied by Mössbauer spectroscopy in order to determine the iron-bearing phases and ratios of their contents in [107]. The room temperature Mössbauer spectra of basalt, specimen 70017, and poikilitic rock, specimen 60315, are shown in Figure 44. These spectra demonstrated that the basalt sample contained a small amount of troilite and large amounts of pyroxene (components related to the M1 and M2 sites) and ilmenite, while poikilitic rock additionally contained metallic iron alloy and olivine. However, the content of ilmenite was substantially smaller in the latter sample.
Furthermore, the large number of Lunar samples delivered by Apollo missions 11,12,14,15,16, and 17 were studied by Mössbauer spectroscopy in order to determine the iron-bearing phases and ratios of their contents in [107]. The room temperature Mössbauer spectra of basalt, specimen 70017, and poikilitic rock, specimen 60315, are shown in Figure  44. These spectra demonstrated that the basalt sample contained a small amount of troilite and large amounts of pyroxene (components related to the M1 and M2 sites) and ilmenite, while poikilitic rock additionally contained metallic iron alloy and olivine. However, the content of ilmenite was substantially smaller in the latter sample.
(a) (b) Figure 44. Mössbauer spectra of Lunar samples measured at room temperature: basalt, specimen 70017 (a) and poikilitic rock, specimen 60315 (b). F is metallic iron alloy, T is troilite, O is olivine, M1 represents the M1 sites in pyroxene, M2 represents the M2 sites in pyroxene, and I is ilmenite. Adapted from [107].
A comparison of the Mössbauer spectra of Lunar fines delivered by Luna-16 and Luna-20 return missions was carried out in [121]. The room temperature spectra of Lunar Figure 44. Mössbauer spectra of Lunar samples measured at room temperature: basalt, specimen 70017 (a) and poikilitic rock, specimen 60315 (b). F is metallic iron alloy, T is troilite, O is olivine, M1 represents the M1 sites in pyroxene, M2 represents the M2 sites in pyroxene, and I is ilmenite. Adapted from [107].
A comparison of the Mössbauer spectra of Lunar fines delivered by Luna-16 and Luna-20 return missions was carried out in [121]. The room temperature spectra of Lunar fines specimens L1627,2, L1627,9, and L2015,2 measured in the large velocity range are shown in Figure 45a. These spectra demonstrated the presence of magnetic sextet with H eff ≈ 330 kOe in specimens L1627,2 and L1627,9, which is characteristic for ferromagnetic Fe-Ni alloy. However, this magnetic sextet was not revealed in the spectrum of specimen L2015,2. The central (paramagnetic) parts of these spectra were slightly different for the studied fines. Therefore, the authors measured the Mössbauer spectra of specimens L1627,2 and L2015,2 in comparison with the Lunar soil, specimen 15601,104, from the Apollo 15 mission, in the small velocity range (Figure 45b). The central part of these spectra indicated the presence of at least a superposition of three quadrupole doublets with different relative areas (in fact, with different iron content in these phases) but with the same isomer shift and quadrupole splitting for the studied samples. The authors of [121] believed that these quadrupole doublets could be assigned to orthopyroxenes, clinopyroxenes, olivine, ilmenite, glassy phases, and superparamagnetic iron. Parameters for ilmenite were δ = 0.70 mm/s and ∆E Q = 1.06 mm/s. Different contents of ilmenite in specimens L1627,2, L2015,2, and 15601,104 were visible in the spectra (shown in Figure 45b) due to different contributions of the quadrupole doublet related to ilmenite. On the other hand, the authors did not reveal any spectral component which could be related to ferric compound. roxenes, olivine, ilmenite, glassy phases, and superparamagnetic iron. Parameters for ilmenite were δ = 0.70 mm/s and ΔEQ = 1.06 mm/s. Different contents of ilmenite in specimens L1627,2, L2015,2, and 15601,104 were visible in the spectra (shown in Figure 45b) due to different contributions of the quadrupole doublet related to ilmenite. On the other hand, the authors did not reveal any spectral component which could be related to ferric compound.

Martian Surface
The study of the surface of Mars by means of Mössbauer spectroscopy became possible with development of the miniaturized spectrometer MIMOS II (see [133][134][135]). Two robotic rovers of NASA's Mars Exploration Rover Missions named Spirit and Opportunity were equipped with MIMOS II spectrometers that recorded the backscattered

Martian Surface
The study of the surface of Mars by means of Mössbauer spectroscopy became possible with development of the miniaturized spectrometer MIMOS II (see [133][134][135]). Two robotic rovers of NASA's Mars Exploration Rover Missions named Spirit and Opportunity were equipped with MIMOS II spectrometers that recorded the backscattered Mössbauer spectra of the surface targets. These rovers landed on Mars at two different locations: on January4, 2004 at the Gusev Crater (Spirit) and on January 25, 2004 at the Miridiani Planum (Opportunity). The measured Mössbauer spectra transferred to the Earth were fitted and demonstrated the presence of spectral components that could be related to olivine and possibly pyroxene, the paramagnetic ferric compounds, and magnetically split ferric compounds [136][137][138].
The first Mössbauer spectrum measured on the Mars surface (Gusev Crater) by Spirit rover is shown in Figure 46a. This spectrum demonstrated the presence of three quadrupole doublets that were assigned to olivine, probably pyroxene and ferric compound. Another spectrum of Adirondack rock (also Gusev Crater) demonstrated the same three quadrupole doublets and an additional two magnetic sextets that were assigned to the tetrahedral (A) and octahedral [B] sites in magnetite (Figure 46b). These Mössbauer spectra indicated that the Mars surface and rocks contained olivine and pyroxene, as well as some ferric compounds. measured the Mössbauer spectrum of outcropped target "El Capitan", which was different from the spectra measured at the Gusev Crater (see Figure 47). This spectrum demonstrated the presence of basaltic components (olivine and pyroxene), as well as jarosite (K, Na)Fe3(SO4)2(OH)6 and paramagnetic and magnetic ferric compounds (the latter was associated with hematite). Observation of jarosite was the most important because this mineral contains hydroxyl groups, which was a direct confirmation of aqueous alteration processes and the presence of water at this site of Mars in the past.  On the other side of Mars, at the Miridiani Planum, the second rover Opportunity measured the Mössbauer spectrum of outcropped target "El Capitan", which was different from the spectra measured at the Gusev Crater (see Figure 47). This spectrum demonstrated the presence of basaltic components (olivine and pyroxene), as well as jarosite (K, Na)Fe 3 (SO 4 ) 2 (OH) 6 and paramagnetic and magnetic ferric compounds (the latter was associated with hematite). Observation of jarosite was the most important because this mineral contains hydroxyl groups, which was a direct confirmation of aqueous alteration processes and the presence of water at this site of Mars in the past.
The first Mössbauer spectrum measured on the Mars surface (Gusev Crater) by Spirit rover is shown in Figure 46a. This spectrum demonstrated the presence of three quadrupole doublets that were assigned to olivine, probably pyroxene and ferric compound. Another spectrum of Adirondack rock (also Gusev Crater) demonstrated the same three quadrupole doublets and an additional two magnetic sextets that were assigned to the tetrahedral (A) and octahedral [B] sites in magnetite (Figure 46b). These Mössbauer spectra indicated that the Mars surface and rocks contained olivine and pyroxene, as well as some ferric compounds.
On the other side of Mars, at the Miridiani Planum, the second rover Opportunity measured the Mössbauer spectrum of outcropped target "El Capitan", which was different from the spectra measured at the Gusev Crater (see Figure 47). This spectrum demonstrated the presence of basaltic components (olivine and pyroxene), as well as jarosite (K, Na)Fe3(SO4)2(OH)6 and paramagnetic and magnetic ferric compounds (the latter was associated with hematite). Observation of jarosite was the most important because this mineral contains hydroxyl groups, which was a direct confirmation of aqueous alteration processes and the presence of water at this site of Mars in the past. (a) (b) Figure 47. Mössbauer spectra of outcropped target "El Capitan" from Miridiani Planum, Mars measured in the large velocity range (a) and in the small velocity range (b). Adapted from [138]. Figure 47. Mössbauer spectra of outcropped target "El Capitan" from Miridiani Planum, Mars measured in the large velocity range (a) and in the small velocity range (b). Adapted from [138].
In addition to jarosite, goethite α-FeOOH was found in the Mössbauer spectra of rock "Clovis" in the Columbia Hills at the Gusev Crater [139]. Two Mössbauer spectra measured at slightly different temperature ranges are shown in Figure 48. These spectra demonstrated the presence of goethite in the magnetic and paramagnetic states, thus testifying to the presence of water on Mars in the past.
Various Mössbauer measurements of soils and rocks by two rovers at the Gusev Crater and Miridiani Planum were considered, e.g., in [140][141][142][143][144]. The authors determined the iron-bearing minerals in these soils and rocks and obtained the modal phase analyses of the studied materials. The iron-bearing phases detected at the Gusev Crater by Spirit rover Mössbauer spectrometer were olivine, pyroxene, ilmenite, magnetite, chromite, nanophase ferric oxide, Fe 3+ -sulfate, hematite, goethite, and pyrite/marcasite. The ironbearing phases identified by Mössbauer spectroscopy with the Opportunity rover at the Meridiani Planum were jarosite, hematite, olivine, pyroxene, magnetite, nanophase ferric oxides, an unassigned ferric phase, and metallic iron alloy (kamacite). Selected Mössbauer spectra of various soils and rocks from the Gusev Crater are shown in Figure 49.
In addition to jarosite, goethite α-FeOOH was found in the Mössbauer spectra of rock "Clovis" in the Columbia Hills at the Gusev Crater [139]. Two Mössbauer spectra measured at slightly different temperature ranges are shown in Figure 48. These spectra demonstrated the presence of goethite in the magnetic and paramagnetic states, thus testifying to the presence of water on Mars in the past.
(a) (b) Figure 48. Mössbauer spectra of rock "Clovis" in the Columbia Hills at the Gusev Crater, Mars, measured at two temperature ranges: 210-220 K (a) and 230-240 K (b). Indicated spectral components are the results of the fits. Adapted from [139].
Various Mössbauer measurements of soils and rocks by two rovers at the Gusev Crater and Miridiani Planum were considered, e.g., in [140][141][142][143][144]. The authors determined the iron-bearing minerals in these soils and rocks and obtained the modal phase analyses of the studied materials. The iron-bearing phases detected at the Gusev Crater by Spirit rover Mössbauer spectrometer were olivine, pyroxene, ilmenite, magnetite, chromite, nanophase ferric oxide, Fe 3+ -sulfate, hematite, goethite, and pyrite/marcasite. The ironbearing phases identified by Mössbauer spectroscopy with the Opportunity rover at the Meridiani Planum were jarosite, hematite, olivine, pyroxene, magnetite, nanophase ferric oxides, an unassigned ferric phase, and metallic iron alloy (kamacite). Selected Mössbauer spectra of various soils and rocks from the Gusev Crater are shown in Figure 49.
Some particular studies by Mössbauer spectroscopy were carried out for investigations of Martian hematite [145,146], mineralogy of several Martian cobbles [147,148], shergottite-like basalt on Mars [149], and Martian olivine and pyroxene [150]. It is interesting to observe that the chemical and mineralogical data for Martian shergottite-like basalt were similar to those for Martian meteorites shergottites found on Earth, especially for Elephant Moraine (EETA) 79001 (see Mössbauer data in Section 2.4 and [26]) and Queen Alexandra Range 94201.
Moreover, Opportunity rover discovered the first iron meteorite at the Meridiani Planum, Mars which was named Maridiani Planum. The Mössbauer spectrum of the brushed part of Meridiani Planum is shown in Figure 52 [151] and demonstrates the presence of a magnetic sextet and paramagnetic doublet with the following parameters (for the temperature range 220-230 K): (i) δ = 0.01 mm/s, Heff = 343 kOe, A = 91% and (ii) δ = 0.34 mm/s, ΔEQ = 0.75 mm/s, A = 9%. The magnetic sextet was assigned to kamacite (the Ni content was ~7 wt.%). The paramagnetic doublet was associated with paramagnetic ferric compound. Moreover, Opportunity rover discovered the first iron meteorite at the Meridiani Planum, Mars which was named Maridiani Planum. The Mössbauer spectrum of the brushed part of Meridiani Planum is shown in Figure 52 [151] and demonstrates the presence of a magnetic sextet and paramagnetic doublet with the following parameters (for the temperature range 220-230 K): (i) δ = 0.01 mm/s, H eff = 343 kOe, A = 91% and (ii) δ = 0.34 mm/s, ∆E Q = 0.75 mm/s, A = 9%. The magnetic sextet was assigned to kamacite (the Ni content was~7 wt.%). The paramagnetic doublet was associated with paramagnetic ferric compound.

Conclusions
About 60 years of successful applications of 57 Fe Mössbauer spectroscopy in the investigation of extraterrestrial materials have demonstrated the great utility of this technique in the study of various iron-bearing minerals and phases. Almost all groups/clans of meteorites have been investigated by Mössbauer spectroscopy. Lunar soils and rocks were also studied by Mössbauer spectroscopy after delivery to the Earth by returned Lunar space missions. In the case of Mars, two rovers Spirit and Opportunity equipped with miniaturized Mössbauer spectrometers measured the Mössbauer spectra in situ of Martian soils and rocks, as well as meteorites that fell on Mars. The application of synchrotron Mössbauer spectroscopy demonstrated the unique possibility of measuring the selective spectra of thin slices of meteorite with a spatial resolution ~10-20 μm. Mössbauer spectroscopy with a high velocity resolution demonstrates the possibility to excavate spectral components including minor ones which cannot be extracted correctly from the spectra

Conclusions
About 60 years of successful applications of 57 Fe Mössbauer spectroscopy in the investigation of extraterrestrial materials have demonstrated the great utility of this technique in the study of various iron-bearing minerals and phases. Almost all groups/clans of meteorites have been investigated by Mössbauer spectroscopy. Lunar soils and rocks were also studied by Mössbauer spectroscopy after delivery to the Earth by returned Lunar space missions. In the case of Mars, two rovers Spirit and Opportunity equipped with miniaturized Mössbauer spectrometers measured the Mössbauer spectra in situ of Martian soils and rocks, as well as meteorites that fell on Mars. The application of synchrotron Mössbauer spectroscopy demonstrated the unique possibility of measuring the selective spectra of thin slices of meteorite with a spatial resolution~10-20 µm. Mössbauer spectroscopy with a high velocity resolution demonstrates the possibility to excavate spectral components including minor ones which cannot be extracted correctly from the spectra measured using conventional spectrometers.
The main results which can be obtained using Mössbauer spectroscopy of extraterrestrial materials include the following: (i) decomposition of the Mössbauer spectra of various materials with determination of the 57 Fe hyperfine and other Mössbauer parameters; (ii) determination of the iron-bearing phases/minerals using the 57 Fe hyperfine parameters; (iii) comparison of the 57 Fe hyperfine parameters for the same phases/minerals in different meteorites/planetary materials in order to analyze the microstructural variations in the iron local microenvironments and evaluate similarities and differences in their formation and evolution; (iv) modal phase analysis of the iron-bearing phases/minerals and its comparison for different samples; (v) determination of the Fe 2+ partitioning among the M1 and M2 sites in silicate crystals (olivine, orthopyroxene and clinopyroxene), as well as iron occupancies of the M1, M2, and M3 sites in iron nickel phosphides (and in any other extraterrestrial crystals consisting of non-equivalent sites for iron); (vi) estimation of the distribution coefficient and the temperature of equilibrium cation distribution for olivines and orthopyroxenes in order to analyze their thermal history (with a thermal history of extraterrestrial materials containing these silicates); (vii) investigation of the processes of terrestrial and extraterrestrial weathering of meteorites and planetary materials, as well as the fusion crust formation; (viii) estimation of the terrestrial age of meteorites; (ix) classification of meteorites using Mössbauer parameters. Therefore, Mössbauer spectroscopy will be further used for the study of meteorites and planetary materials. However, to obtain reliable results, it is necessary to use the most appropriate Mossbauer spectrometers and carry out measurements with high quality and precision that can produce the high-quality spectra. Moreover, the most careful studies should be done using other powerful complementary techniques such as X-ray diffraction, scanning electron microscopy with energy/wavelength dispersive spectroscopy, electron microprobe analysis, magnetization measurements, and optical microscopy (metallography). A multidisciplinary approach to study extraterrestrial materials that also includes Mössbauer spectroscopy is most successful because there is the possibility of verifying the quality of the data by comparing the analyses performed with different techniques.