Impact of high temperatures on aluminoceladonite studied by Mössbauer, Raman, X-ray diffraction and X-ray photoelectron spectroscopy

Mössbauer, Raman, X-ray diffraction and X-ray photoelectron spectroscopies were used to examine the effects of temperature on the structure of two aluminoceladonite samples. The process of oxidation of Fe2+ to Fe3+ ions started at about 350 °C for the sample richer in Al and at 300 °C for the sample somewhat lower Al-content. Mössbauer results show that this process may be associated with dehydroxylation or even initiate it. The first stage of dehydroxylation takes place at a temperature > 350 °C when the adjacent OH groups are replaced with a single residual oxygen atom. Up to ~500 °C, Fe ions do not migrate from cis-octahedra to trans-octahedra sites, but the coordination number of polyhedra changes from six to five. This temperature can be treated as the second stage of dehydroxylation. The temperature dependence on the integral intensity ratio between bands centered at ~590 and 705 cm−1 (I590/I705) clearly reflects the temperature at which six-coordinated polyhedra are transformed into five-coordinated polyhedra. X-ray photoelectron spectra obtained in the region of the Si2p, Al2p, Fe2p, K2p and O1s core levels, highlighted a route to identify the position of Si, Al, K and Fe cations in a structure of layered silicates with temperature. All the measurements show that the sample with a higher aluminum content and a lower iron content in octahedral sites starts to undergo a structural reorganization at a relatively higher temperature than the less aluminum-rich sample does. This suggests that iron may perform an important role in the initiation of the dehydroxylation of aluminoceladonites.


Introduction
The hydrous layer silicates, commonly known as clay minerals, are widely used in many industrial applications such as paper production, foundries and the chemical-and food industries (Hradila et al. 2003;Murad and Fabris 2010;Moretto et al. 2011). A 2:1 layer is the basic structural unit of various 2:1 phyllosilicates. It consists of an edge-sharing octahedral sheet sandwiched by two corner-sharing tetrahedral sheets. In general, the octahedral sheet of a 2:1 layer contains three symmetrically independent sites differing in the mutual arrangement of OH groups and oxygen atoms coordinating octahedral cations. The unit cell contains two symmetrically independent octahedral, i.e., one trans-ostahedra (M1) and one cis-octahedra (M2) differing in the positions of OH groups. The adjacent cis-octahedra share an edge formed by two OH groups whereas, in trans-octahedra, hydroxyls are located opposite to each other (Muller et al. 1999;Zviagina et al. 2007).
The thermal behavior and related process of dehydroxylation of dioctahedral silicates has been the subject of many investigations (Tsipursky et al. 1985;Muller et al. 2000aMuller et al. , 2000bMuller et al. , 2000c. Dehydroxylation slightly modifies the coordinates of atoms other than the residual oxygen which replaces Editorial handling: A. Beran two hydroxyl groups and is located in the plane of the octahedral cations midway between them. This replacement induces a change of the coordination of the octahedral sheet cations, becoming five instead of six. Also, it has been shown quantitatively that the dehydroxylation of Fe, Mg-rich dioctahedral micas, when heated, is accompanied by a migration of the octahedral sheet cations from cis to trans-sites (Drits et al. 1995;Muller et al. 1998Muller et al. , 2000cDainyak et al. 2006). Generally, on heating phyllosilicates, four different processes may occur separately or concurrently, namely, oxidation of cations, loss of water, migration of cations to different sites and disintegration of the structure. However, many questions still remain unanswered about the coordination of the octahedral cations after dehydroxylation, the details of the reaction mechanism and the fate of dehydroxylates formed in nature.
In this study, X-ray diffraction, X-ray photoelectron spectroscopy, Mössbauer and Raman analysis were used to examine the effects of heating on the structure of Febearing dioctahedral mica. The study tested conclusions about the structure of dehydroxylates and octahedralcation migration in phyllosilicates and attempts a global explanation of the role of the content of Fe in the octahedral sheet in the dehydroxylation process. Additionally, Xray photoelectron spectroscopy was used as a new tool to visualize cation arrangement in a structure of layer silicates at high temperatures.

Experimental details
Materials Several secondary minerals have been found in the altered volcanic rocks of intermediate to basaltic composition in Lubiechowa near Świerzawa (Sudety Mountains) from which two types of clay mineral were collected from amygdales, one emerald-green, the other of dark-green. Based on electronmicroprobe analyses, calculated formulas for these are K 0.629 Ca 0.041 (Mg 0.545 Fe 2+ 0.180 Mn 0.002 )(Fe 3+ 0.430 Al 0.950 )(S-i3.755Al 0.245 )O 10 (OH) 2 for emerald-green "celadonite" and K 0.766 Ca 0.031 (Mg 0.584 Fe 2+ 0.134 Mn 0.001 )(Fe 3+ 0.562 Al 0.655 )(S-i3.962Al 0.038 )O 10 (OH) 2 for dark-green "celadonite", where FeO and Fe 2 O 3 contents were determined by the use of the Fe 3+ /Fe 2+ values obtained from Mössbauer spectra. According to criteria for the classification of clay minerals (Rieder et al. 1998;Drits et al. 1997;Tischendorf et al. 2007;Zviagina et al. 2015), the samples cannot be considered celadonite, but aluminoceladonite. Here, they are named as aluminoceladonites AlCel1 and AlCel2, respectively. Several types of measurements were performed for the studied samples.

X-ray diffraction
In-situ, high-temperature X-ray diffraction (XRD) measurements were performed by CuK α radiation using a PANalytical, X'Pert Pro Multi-Purpose Diffractometer with a Philips PW3040/60 X-ray generator and an Anton Paar HTK-16 thermal chamber. A nickel filter was placed in the diffracted-beam path to reduce the K β radiation. XRD patterns were collected in the range 4 o 2θ to 90 o 2θ with a step size of 0.01 o 2θ and nominal time per step of 300 s, with the aid of a scanning X'Celerator detector (with an active length of 2.12 o 2θ).

Thermal analysis
The thermogravimetric analyses (TGA) and differential scanning calorimetry (DSC) were performed using the thermal analyzer NETZSCH Jupiter STA 449 F3. The curves were measured in the temperature range 25-800°C at a heating rate of 5°C/min. The samples (about 25-30 mg aliquots of powder) were measured in Al 2 O 3 crucibles in a protective argon atmosphere.

57
Fe Mössbauer spectroscopy 57 Fe Mössbauer transmission spectra were recorded at room temperature using a MS96 spectrometer and a linear arrangement of a 57 Co:Rh (25 mCi) source, a multichannel analyzer, an absorber and a detector. The spectrometer was calibrated at room temperature with a 30 μm thick α-Fe foil. The samples were pulverized in an agate mortar and heated for 8 h in air under static conditions in an electric furnace to a variety of temperatures in the range 100-600°C. After each heating treatment, the Mössbauer measurements were taken. Numerical analysis of the Mössbauer spectra was performed using the WMOSS program (Prisecaru 2009(Prisecaru -2016.

Raman spectroscopy
Raman spectra of the samples were obtained using a WITec confocal Raman microscope CRM alpha 300R equipped with a solid-state laser (λ = 532 nm) and a charge-coupled device (CCD) camera. The excitation laser radiation was coupled to the microscope through a single-mode optical fibre with a 50 μm diameter. An air Olympus MPLAN (100×/0.90NA) lens was used. The Raman spectra were initially laserexcited for a few hours to eliminate afterglow-type luminescence, and decrease the so-called Raman background effect. The spectrometer monochromator was calibrated using the emission lines of a Ne lamp. The temperature measurements were performed in-situ using a Linkam heating/cooling stage. The crystals of aluminosilicate were heated at a rate of 5°C/min with a temperature accuracy of 0.1°C in air. X-ray photoelectron spectroscopy XPS spectra were obtained using a PHI5700/660 Physical Electronics Photoelectron Spectrometer with monochromatic Al K α X-ray radiation (1486.6 eV). The energy of the electrons was measured with a hemispherical analyzer with resolution of about 0.3 eV. Photoelectron emission was measured from a surface area with a diameter of 800 μm and at a takeoff angle of 45°. Due to the occurrence of a charge effect, a neutralizer was used for nonconductive samples. To determine the binding energy (BE), the C1s component was set at 285 eV. The solid sample was mounted on a PTS CLEAV RES/C-K sample holder (Prevac Company), cleaved in UHV conditions in a vacuum of about 3 . 10 −9 Tr and transferred to the XPS chamber. The series of measurements in situ were performed at temperatures of 100, 300 and 550°C. The vacuum in the XPS chamber during the temperature measurements did not exceed 5 . 10 −8 Torr. The increase and stabilization of each temperature was executed over about 30 min using the Sample Heating Power Supply HEAT2-PS device and controlled by the use of EUROTHERM (Prevac Company). The HEAT2-PS worked in a resistive mode with a ramp of heating current of 1A/min.

X-ray diffraction
High-temperature XRD patterns (Fig. S1 in the Electronic Supplementary Material (ESM)) were collected every 50°C ranging from room temperature to 700°C and refined by the use of Rietveld method. Table 1 shows the unit-cell parameters of the investigated samples obtained for the temperature range < 700°C. The migration of octahedral cations from former cis-octahedra strongly decreases the repulsion between those cations bonded to residual oxygen anions and can be reflected by the temperature-dependence of the a and b parameter values (Table 1; Fig. 1). The decrease of a for AlCel1 is at 600°C. For AlCel2, no significant decrease of a can be observed in this temperature region. The minimum values of the b parameters are observed near 650°C and indicate that migration of octahedral cations takes place in these wide temperature ranges. Also, with the pure dehydroxylation increases the interlayer thickness, an increase which is reflected in the increase of the c parameter. This rapid increase of this parameter comes after heating at 550°C for AlCel1 and 500°C for AlCel2, which is also confirmed by the variation of the unitcell volume versus temperature (Fig. 1). The relative displacement of the adjacent layers along the a axis is characterized by dependence of the Ic . cosβ/aI value on temperature (Fig. 1). The deviation of the observed Ic . cosβ/aI from the ideal −0.333 value ) depends significantly on the difference in the sizes of the vacant and occupied octahedral site of the 2:1 layer. It is expected that increasing octahedral Al will be reflected in deviation decreases of the Ic . cosβ/aI values. The changes of this value from about 0.363 for the sample of AlCel1 heated at 550°C to 0.303 for that heated at 700°C, and from 0.367 for the sample of AlCel2 heated at 500°C to 0.313 for that heated at 700°C also indicate the transformation of cis into trans sites in these temperature ranges. More precisely, the coordination number of polyhedrons change from six into five. Detailed analysis of the unitcell parameters (Table 1; Fig. 1) reveal additional effects which take place near temperatures of 100°C and 250°C for AlCel2 and, to a much weaker degree, at 350°C for AlCel1. Changing values of the b and Ic . cosβ/aI parameters at 100°C can be a result of the loss of bound molecular waters in the structure of the AlCel2. The decrease in unit-cell volume and unit-cell parameters at 350°C for AlCel1, and at 250°C for AlCel2, can be a result of deprotonation via Fe oxidation, which reduces the octahedral thickness.

Thermal analysis
The dehydroxylation temperature can be determined from the DSC-TGA curves, as depicted in Fig. 2. These curves show some thermal movements that occur in the aluminoceladonite samples over a temperature range of~30°C to 800°C. TGA/ DSC experiments for the investigated aluminoceladonites show that a small loss of~0.84 wt% occurs in the temperature range 50-100°C for AlCel1 (Fig. 2a) and~0.65 wt% between 50 and 150°C for AlCel2 (Fig. 2a). The weight decrease, and the broad endothermic hump between the above temperature ranges, presumably result from the loss of loosely-bound molecular waters. For AlCel1, in the temperature ranges 100-550°C and 550-780°C, two areas of weight loss (3.15 wt% a n d 2 . 5 5 w t % ) a r e o b s e r v e d a n d i n d i c a t e t h e dehydroxylations. The endothermic events that occurred in the ranges 425-567°C and 567-675°C support this process. In the first region, the main sharp peak is seen at a temperature of 515°C, preceded by a small peak at 478°C. This region is deemed the first stage of dehydroxylation. It occurs in a wide temperature range, which is confirmed by a breadth of the endothermic hump which is about 142°C wide. In the second temperature range, another endothermic peak related to the second stage of dehydroxylation is visible at 625°C. For AlCel2, the dehydroxylations indicated by the mass change of <5.1% were detected from 200 to 660°C. On the DSC curve for AlCel2, endothermic peaks though visible are not as sharp as for AlCel1 and appear at different temperatures than for AlCel1. The first peak is visible at 408°C and a further two at 484°C and 510°C. For both samples, an endothermic peak was detected at 736°C. This is probably the temperature at w hich the decompositi on of the aluminoceladonites starts (Savko et al. 2015).

Fe Mössbauer spectroscopy
Mössbauer spectra of the initial samples were fitted as a superposition of six doublets, the hyperfine parameters (Table S1 in the ESM) of which indicate iron located only in the octahedral sheet. The initial sample of AlCel1 contains about 27% Fe 2+ ions. The Fe 3+ ions prefer the cis-M2 site where 70% of them are located and the rest are present in trans-M1. In AlCel2, Fe 2+ ions constitute 19% of total iron and are approximately equally divided between two octahedral sites with different local geometric environments around the Fe 2+ ions. The hyperfine parameters of the ferric doublets indicate that the majority of the Fe 3+ is located in the cis-M2 site in AlCel2. The remainder of these ions are in trans-M1 sites. It should be noted that, the ferrous doublets do not correspond to a specific crystallographic site (M1 or M2), but arise from local geometric and/or electronic distortion environments around the Fe atoms within the octahedral sheet of the micas (Rancourt 1994).
Selected Mössbauer spectra obtained after heating in air at definite temperatures are presented in Figs. 3 and 4. Figure 5 presents the changes in the hyperfine parameters for each component in the AlCel2 sample after heating at high temperatures. The fact that full line width at half maximum for all components connected with Fe 3+ was the same during the fitting procedure for all Mössbauer spectra allows tracking of the changes in hyperfine parameters each component in. The full line width of doublets connected with Fe 2+ changed, a change connected with visible changes in a shape of the Mössbauer spectrum and with processes taking place in the aluminoceladonite structure. Due to the dehydration stage, the escape of defects and changes that take place in neighboring octahedral sites, small changes in a shape of the Mössbauer Table 1 Temperature dependence of unit-cell parameters for the investigated aluminoseladonites calculated from the peak positions of the XRD patterns AlCel1 AlCel2  (8) 100.792(4) 5.2367(6) 9.0408 (7) 10.2187(9) 475.049(1) 100.918 (8)  450 5.2335(6) 9.0350(1) 10.1900(2) 473.329(4) 100.782 (7) 5.2375(9) 9.0342 (7) (7) 474.257 (7) 100.730(9) 5.2366(9) 9.0096 (2) (9) 5.2367 (7) 8.9922 (1) 10.4096(3) 483.256(5) 99.645 (4)  700 5.2215(2) 8.9859 (7) 10.3596 (4) 480.378 (5) 98.783(2) 5.2365(9) 8.9882 (8) 10.4021(1) 483.478(2) 99.068 (2) Fig . 1 The temperature-dependence of the a) b parameter, b) unit-cell volume, c) Ic . cosβ/aI values calculated from the XRD patterns for the investigated aluminoceladonites. The lines indicate the trends of the data and act as a guide to the eye spectra and hyperfine parameters of all components after heating at 100°C is seen. After heating at 200°C, an increase of abundance of Fe 2+ that is generated by the presence of structural water is visible. The oxidation of Fe 2+ to Fe 3+ started at~350°C for AlCel1 and at 300°C for AlCel2. This process took place over a wide temperature range and finished~150°C higher for each aluminoceladonite (Fig. 5). It may be assumed that at a temperature of~450°C for AlCel1 and 350°C for AlCel2, the first stage of dehydroxylation takes place. Up to a temperature of~550°C for AlCel1 and 500°C for AlCel2, the contribution of components connected with Fe 3+ ions in M1 sites increases rapidly (Fig. 5c). To be precise, the two new doublets represent Fe 3+ ions in M1 sites as a result of the decline of those connected with M2. The observed effects on the Mössbauer spectra at these temperatures can be treated as the second stage of dehydroxylation. It is also reflected by the changes in the contribution of Fe in M1 and M2 sites versus heating temperature ( Fig. 6a and

Raman spectroscopy
Raman spectra were collected every 50°C ranging from room temperature to 600°C. Post-processing analysis including baseline and luminescence correction as well as cosmic-ray removal was done using WitecProjectFour, and peak-fitting analysis using the GRAMS software package. Peak fitting was done by means of the Voigt function with a minimum number of components. The spectra of both investigated samples are presented in Fig. 6. In the hydroxyl-stretching range from 3400 to 3700 cm −1 , five bands at 3627, 3608, 3582, 3564 and 3541 cm −1 for AlCel1 are seen, and five at 3620, 3609, 3588, 3567 and 3542 cm −1 for AlCel2 ( Fig. 6; Table 2). The two bands with the highest Raman shifts at 3627 and 3620 cm −1 are interpreted as AlOHAl stretching vibrations.
In the light of previous analysis for clay minerals (Zviagina et al. 2004(Zviagina et al. , 2015(Zviagina et al. , 2020, other bands represent AlMgOH, 2MgOH, MgFe 3+ OH and/or AlFe 2+ OH with a possible contribution from AlFe 3+ OH, and 2Fe 3+ OH. The temperaturedependent experiment revealed only slight alteration in the band arrangement in the hydroxyl region (3500-3600 cm −1 ). The main changes were found to be connected to modification of the band intensity, whereas the band position linked to the vibration of hydroxyl groups at tetrahedral and octahedral units, especially at MgFe 2+ OH or AlFe 2+ OH, AlFe 3+ OH, 2Fe 3 + OH and MgFe 3 + OH, remained unchanged. Interestingly, temperatures >500°C (AlCel1) and > 400°C (AlCel2) practically deactivate vibration originating from the OH modes due to the first stage of dehydroxylation. This effect can also be explained as a result of molecular mobility which, in the case of aluminoceladonite, may lead to the loss of the hydrogen-bonding scheme. This observation is strictly correlated with the analysis of bands ascribed to the stretching vibration within (Al n Si 4-n O 10 ) 3− units, i.e., in the 1150-950 cm −1 region (Povarennykh 1978). Here, the intensity of the main bands centered at 1130 and 1083 cm −1 (AlCel1) and at 1144 and 1113 cm −1 (AlCel2) decrease up to 500°C and to 400°C, respectively. With both samples, those bands vanish completely above these temperatures. The temperature-dependent experiment also revealed a decline in hydrogen bonding strength that may be linked to the molecular re-arrangement within the tetrahedral sheets, including mutual rotation of silica tetrahedra. The reason for this structural reconfiguration is not completely clear, especially in the context of band analysis of this 1150-950 cm −1 region. Therefore, to shed light on this problem, the analyses of the bands located below 950 cm −1 were taken into consideration.  Here, huge differences were found in the position of bands and their intensities in analysis of the 750-500 cm −1 rangedifferences linked to the molecular vibration within tetrahedral and octahedral sheets (Fig. 6). For the record, the bands centred in between 710 and 690 cm −1 correspond to the overlapping signal of the stretching vibration of aluminum in tetrahedral and octahedral coordination as well as nonbridging and bridging oxygen in relation to the A g symmetry mode (McKeown et al. 1999). The appearance of the band between 690 and 600 cm −1 in turn may originate from the translational modes of i) nonbridging oxygen in the yz plane, ii) bridging oxygen in the yz plane as for the A g symmetry mode and iii) non-bridging oxygen in the xz plane as for the B g symmetry mode (McKeown et al.1999;Wang et al. 2015). Similarly, well-defined bands located at 600-500 cm −1 derived from the overlapping signal of ions occupying the tetrahedral as well as octahedral coordination. Tetrahedrally-coordinated ions are involved in the bending motion related to the xy plane (A g and B g symmetry mode), while octahedrally-coordinated ions bending motion in the translational modes in the xz plane (McKeown et al. 1999). Based on temperature-dependent data, bands located at 549 cm −1 (AlCel1) or 553 cm −1 (AlCel2) in the initial samples were shifted, respectively, towards lower frequencies and are found at 523 cm −1 and 544 cm −1 in the samples annealed at 400°C and 350°C, respectively. At the same time, an increase of full width at half maximum (FWHM) suggests strong molecular reorganization due to the heating. It is notable that those alterations correlate quite well with modifications within the hydroxyl group region and within tetrahedral sheets. This may suggest that the decrease in the hydrogen bonding scheme induced simultaneous rotation of tetrahedral-and octahedral units. Furthermore, the Fig. 6 Raman spectra of the initial sample of aluminoceladonites as well as collected at different temperatures of sample obtained for a) AlCel1 (left column) and for b) AlCel2 (right column); the zone of hydroxyl stretching are added as insert. Bands that are highlighted dark grey correspond to the overlapping signals of tetrahedral and octahedral units, and bands that are highlighted light grey represent vibrations of aluminum-iron octahedrons and tetrahedrons octahedral positions are not fully charge-balanced due to isomorphic substitutions occurring between Al 3+ , Fe 3+ , Mg 2+ and Fe 2+ that are atypical for celadonite micas. As a result, the temperature may additionally intensify variable structural effects by introducing structural point defects into the crystal structure, as is especially well seen in the case of AlCel1 where many new low intensity bands appear on the Raman spectrum (Fig. 6).

X-ray photoelectron spectroscopy
One of the most powerful direct methods for investigating electronic structure, chemical bonds, and the nearest surroundings of atoms in substances is X-ray photoelectron spectroscopy. As XPS signals arise from within about 1-10 nm of the outermost surface layers of investigated samples (Elmi et al. 2016), this method is generally treated as a surface sensitivity method. Discrepancies, e.g., apparent shifts of the photoelectron peaks, can be explained as due to selective attenuation of some signals (Bhattacharyya 1993). Shifts due to surface charging do not depend on the thickness of the sample, but on whether or not there is an adsorbed layer on the mica surface that reduces the charging and, hence, the shift is also reduced. Moreover, many solid systems can exhibit different compositions at their surfaces due to oxidation, layering, or other surface forces. However, the surface analysis of clay minerals is indeed representative of the overall bulk composition. This is to be expected since clay mineral layers should be of similar composition whether at the surface or in the center of an aggregate. Therefore, the interpretation of the XPS spectra of the investigated aluminoceladonite can provide information on the molecular environment, i.e., oxidation state, multiplet structure, chemical bonding, etc. with an error of <10%. This can enable visualization of cation arrangement in a layersilicate structure at high temperatures. Reference survey spectra obtained in a wide range of energies, and high-resolution  Ospitali (2008) and Zviagina (2020) RT = room temperature core lines of C1s, O1s, Si2p, Al2p, Fe2p and K2p, were recorded at room temperature. Then, the series of measurements in situ were performed at temperatures of 100, 300 and 550°C. Quantification of XPS spectra utilizing peak area and peak height sensitivity factor was possible thanks to the Multipak Physical Electronics application. The XPS core level spectra were fitted with Gaussian-Lorentzian functions. XPS spectra obtained in the region of the Si2p, Al2p, Fe2p, K2p and O1s core levels for the aluminoceladonite AlCel2 are shown on Fig. 8. The Si2p spectrum should show a doublet branching ratio of 2:1 for Si 2p 3/2 to Si 2p 1/2 with a spin-orbit splitting of about 1 eV difference. Values of Si2p binding energies show Si 4+ in tetrahedral coordination (Bhattacharyya 1993;Biino and Gröning 1998;Elmi et al. 2014Elmi et al. , 2016. These spectra in Fig. 8 present two Si2p (3/2, 1/2) doublets. The silicon binding energy shows that different populations of silicon correlate with substitutions in the tetrahedral sheet. A lower Si2p binding energy may result from a negative charge increase on the silicate framework as a result of Si 4+ replacement by Al 3+ . Our results indicate that the tetrahedral sheet contains only Si in tetrahedral coordination with some Al 3+ substitution. On the basis of the area ratio of Si2p components, the ratio of Si 4+ in tetrahedral coordination to Al 3+ can be calculated; for AlCel2 this ratio is 72/28. Because, the Al2p spectra have evident asymmetry in peak shape, they were decomposed into two Al2p doublets (Fig. 8). The obtained Al binding energy value is very close to both fourfold and sixfold coordination. However, an octahedral sheet must compensate for the charge of two relatively negative and adjacent tetrahedral sheets. The positive nature of the octahedral sheet coupled with the relatively negative silicate sheet results in more ionic Al-O bonds. The enhanced ionicity of the Al-O bond may be reflected in the relatively higher Al2p binding energy value. Also, the presence of OH and/or F enhances the positive character of the Al resulting in a moderate increase in the Al2p binding energy (Barr et al. (1997)). That is why the Al2p 3/2 line located at a binding energy of 74.3 eV relates to Al present in octahedral coordination, while the other line with a BE of 72.8 eV arises from Al in tetrahedral coordination. The ratio of Al ions in tetrahedral to those in octahedral coordination is 27/73 and is very close to those obtained for Si2p, supporting the correctness of our interpretation. Heating at 100°C causes the decrease of the intensity lines located at BE of 102.8 eV for the Si2p spectrum and 74.3 eV for the Al2p spectrum. The calculated ratio of Si 4+ in tetrahedral coordination to Al 3+ which substituted Si 4+ is 85/15 and the ratio of Al ions in tetrahedral coordination to these in octahedral coordination is 14/86. Heating at 300°C does not lead to significant changes in the Si2p and Al2p spectra but, after heating at a temperature of 550°C, these spectra become almost symmetrical. The ratios of the intensity lines mentioned above for Si are 91/9 at 300°C and 97/3 at 550°C. For the Al2p spectra, the obtained ratios amount to 12/88 at 300°C and 97/3 at 550°C. Potassium K2p has clearly spaced spin-orbit components (Fig. 8). Because, the K2p spectra have asymmetric peak shapes, they were decomposed into two K2p doublets. The K2p 3/2 binding energies at 294.4 eV and 292.8 eV are connected with the K + ions occupying interlayer sites between two 2:1 layers. The area ratio for the K2p doublets is 80/20 for AlCel2, which results from the ratio of K + and Ca 2+ , respectively. Obtained the calcium content is too high and it may suggest that not only Ca 2+ is substituting the K + in interlayer sites. This explanation is supported by K2p spectra obtained at 100°C and 300°C. As is seen on Fig. 8, the intensity line of K2p 3/2 with a binding energy of 294.4 eV decreases and the area ratio for the doublets of K2p is 88/12 at 100°C and 95/5 at 300°C. This is also supported by XPS spectra of O1s regions obtained at the studied temperatures (Fig. 8). The O1s core-level line of the initial sample of dark-green celadonite comprises three separate peaks located at binding energies of 530.2, 531.7 and 533 eV, which can be related to the existence of different forms of oxygen bondings. The most intense line at a binding energy of 531.7 eV is typical of lattice oxygen O 2− (Kloprogge and Wood, 2018) and the peak positioned at the highest binding energy of 533 eV corresponds to hydroxyl oxygen OH − (Schingaro et al. 2013). The lowest intensity line at 530.2 eV is due to the presence of O-Me bonds or H 2 O or the sum of both. The intensity of this line decreases after heating at 100°C, which confirms the presence of a small amount of structural water in interlayer sites in the initial sample. Heating at 100°C also causes an increase of the intensity line corresponding to hydroxyl oxygen at a binding energy of 533 eV; this can be an effect of reaction of the water molecule H 2 O with oxygen O 2− creating OH − and OH groups. At 300°C, the ratio of the peaks located at binding energies of 531.7 and 533 eV is almost the same but, at 550°C, the intensity of the hydroxyl oxygen line increases again.
The Fe2p spectra are characterised by doublet structures relating to the spin-orbit splitting between 2p 1/2 and 2p 3/2 states. In addition to the multiplet structures, shake-up satellite peaks arising from intrinsic energy losses when the photoelectron leaves the hosting atom, are present. Deconvolution of the Fe2p peak (Fig. 8) shows three main components. The first, a weaker component related to the Fe 2+ valence state, is a doublet with binding energies Fe2p 3/2 = 709.5 eV and Fe2p 1/2 = 723 eV. The second, the highest component, is characteristic of the Fe 3+ state observed as a multiplet with binding energies of the main component Fe2p 3/2 = 712 eV and Fe2p 1/2 = 725.5 eV. The third component at a binding energy = 714.4 eV is associated with the Fe 3+ state. The Fe 2+ /Fe 3+ ratio was determined to be 18/82. The area ratio of the peaks at binding energies = 712 eV and 714.2 eV is 62/20 and reflects the M2/M1 ratio of Fe 3+ ions in octahedral positions. Thẽ 2 eV distance between these components can reflect dehydroxylated sites. Small changes of the Fe2p spectra are observed at 300°C where the obtained Fe 2+ /Fe 3+ ratio is 22/ 78 and the ratio M2/M1 of Fe 3+ ions is 56/22, based on the area ratio of the Fe2p lines. The main changes of the Fe2p spectrum are observe at 550°C where the binding energies of the peaks and their intensities change. All of the doublets observed are connected with Fe 3+ ions. The two peaks located at 710.4 and 712.5 eV can be associated with Fe in M1 sites and that with a binding energy of 714.2 eV with a small contribution of iron in dimeric pyramid sites. The area ratio of peaks at bonding energies = 710.4 eV and 712.5 to that at 714.2 eV is 83/17, which is similar to the ratio 87/13 obtained from the Mössbauer studies at this temperature.

Discussion
TGA-DSC analysis reveal a small weight loss in the temperature range 50-150°C that can result of the loss of bound molecular waters in the structure of the investigated aluminoceladonites. The endothermic events, areas of weight loss and the temperature-dependence of the unit-cell parameters indicate two-stage process of the dehydroxylation. The first stage of dehydroxylation takes place at a temperature of about 450°C for AlCel1 and at a 350°C for AlCel2, when the adjacent OH groups are replaced with a single residual oxygen atom. In turn, temperature about 550°C for AlCel1 and 500°C for AlCel2 can be treated as the temperature where the second stage of dehydroxylation started. Finally, sample decomposition process starts at 736°C.
The crystal-chemical models obtained due to 57 Fe Mössbauer spectroscopy (Tsipursky and Drits (1984); Dainyak et al. 2004Dainyak et al. , 2006Dainyak et al. , 2009Dainyak et al. , 2013Zviagina et al. 2015) have shown that the values of the fitted quadrupole doublets Fe 3+ and Fe 2+ (visible on the Mössbauer spectra) may correspond to specific local cationic arrangements around these ions. Such cation distribution reconstructions can indicate random or ordered distribution of octahedral cations and can also display Fe clustering or domain structure. Due to above, the inner doublets with quadrupole splitting of 0.17 mm s −1 and 0.19 mm s −1 is assigned to Fe 3+ with the three nearest Fe 3+ ions. The quadrupole splitting of 0.49 mm s −1 and 0.54 mm s −1 reflects Fe 3+ in the arrangement of MgFe 2+ Fe 3+ and 2AlFe 3+ , respectively. The value of QS, which is 0.81 mm s −1 , is close to the theoretical value responsible for Fe 3+ in the 3Al arrangement. Though AlCel1 contains more Al ions than AlCel2, no doublets reflecting Fe 3+ ions with 3Al arrangement are observed on the Mössbauer spectra, suggesting that Al ions are steadily ordered in octahedral sites. The large values of quadrupole splitting, i.e., 1.01, 1.29 and 1.44 mm s −1 , associated with the trans-OH (M1) sites, are the consequence of Fe 3+ in arrangements AlMgFe 3+ , 2Fe 3+ Mg and 2AlMg, respectively. For AlCel1, a higher proportion of components related to the Al environment than in AlCel2 is observed due to the higher Al content. A similar interpretation of the quadrupole splitting can be made for the local cation arrangements around Fe 2+ ions. The smallest value of QS (1.66 mm s −1 ) is connected with the arrangement of 2MgFe 2+ ions and the highest value (2.77 mm s −1 ) with 2Fe 3+ Fe 2+ . Intermediate QS reflect an almost similar arrangement, i.e., 2.23 mm s −1 is the result of 2Fe 2+ Al and 2.66 mm s −1 is the result of AlMgFe 2+ local cation distributions. The above interpretation of Mössbauer spectra suggests that Fe 3+ ions are not randomly located in the structure of the studied samples but in some position within the preferred neighbourhood.
According to so-performed analysis, a determined temperatures of the dehydroxylation by analysis of Mössbauer spectra turned out to be in good agreement with the XRD and DSC-TGA results. However, Mössbauer spectroscopy results indicate additional factors responsible for initiation of the proc e s s o f d e h y d r o x y l a t i o n i n t h e i n v e s t i g a t e d aluminoceladonite. Namely, at temperature~350°C for AlCel1 and at 300°C for AlCel2 the oxidation of Fe 2+ to Fe 3 started. In the oxidation process, the charge balance can be provided by the loss of protons; the occurrence of proton jumps between neighboring OH − groups and O 2− ions may be postulated. A proton jump or proton exchange causes a breakdown of the ordering of the OH − groups and the reduction of bond strength between Fe and hydroxyl. In facilitating loss of one of these groups, oxidation of Fe 2+ may thus be associated with dehydroxylation or even initiate the process. The first stage of dehydroxylation takes place as the adjacent OH groups are replaced by a single residual oxygen atom. In this stage of dehydroxylation Fe ions do not migrate from M2 to M1 sites, but coordination numbers undergo change. This replacement of OH by O induces a change in the coordination of the octahedral sheet cations, becoming five instead of six. Additionally, above temperature 550°C for AlCel1 and 500°C for AlCel2, where second stage of dehydroxylation take place, on Mössbauer spectra appear doublet characteristic for Fe 3+ in tetrahedral coordination. This can be the effect of the remodelling of some M1 coordination polyhedrons. The changing of coordination number from five to four is associated with detachment of an OH group or of O with the result that a dimeric pyramid can be formed.
An interesting supply of the previously mentioned techniques provides temperature-dependent Raman and XPS data. First of all, Raman spectra reveal huge differences in the position of bands and their intensities in the 750-500 cm −1 range were found (Fig. 6). More precisely, a gradual decrease of the intensity of the band centered at 706 cm −1 (AlCel1) and 708 cm −1 (AlCel2) was correlated with an increase of the band at 567 cm −1 (AlCel1) and 575 cm −1 (AlCel2). Additionally, the peak-fitting analysis was used to estimate the integrated intensity of two bands at~705 cm −1 and at~590 cm −1 and calculation of the integrated intensity ratio I 590 /I 705 (Fig. 7c). It is worth noting that this approach was purely hypothetical and not known from the literature. An interesting difference emerged, namely, the sample with higher aluminum started to undergo a structural reorganization at a higher temperature (400°C) than the less aluminum-rich sample did (~300°C). It may be speculated that the difference may reflect some change in the structural stability of the celadonite micas. What is more, a plot of temperature-dependent integral intensity revealed a jump of the I 590 /I 705 values at temperatures of 500°C may indicate a phase transition of celadonite into ferri-phlogopite or ferri-biotite. However, the appearance of potentially new phases was countered by the XRD and Mössbauer results. The correlation of the Raman data with the Mössbauer spectroscopy outcomes does raise the possibility of analyzing the contribution of Fe in octahedral sites in relation to temperature. The Mössbauer spectroscopy data suggest that the I 590 /I 705 pattern may reflect the increase in occupancy of the M1 site with increasing temperature (Fig.  7). The observed maximum of the intensity ratio of I 590 /I 705 at 550°C for AlCel1 and 500°C for AlCel2 indicates that the polyhedra have been transformed from six-coordinated into five-coordinated polyhedra. Afterwards, a slow decrease in the curves with temperature ( Fig. 7c) reflects the migration of cations from cis-sites to trans-sites. Those data are in good agreement with results obtained from the XRD analysis and reflect the conclusions taken from the Mössbauer data.
Secondly, X-ray photoelectron spectroscopy were applied to look more precisely on the chemical environment and follow the impact of temperature on the modification of chemical states. According to XPS data, the two Si2p (3/2, 1/2) doublets indicate that the tetrahedral sheet contains only Si in tetrahedral coordination with some Al 3+ substitution. On the basis of the area ratio of Si2p components, the ratio of Si 4+ in tetrahedral coordination to Al 3+ is 72/28. On the other hand, the Al2p spectrum indicate that Al ions are present in octahedral and in tetrahedral coordination in structure of AlCel2. The ratio of Al ions in tetrahedral to those in octahedral coordination is 27/73 and is very close to those obtained for Si2p. Heating the sample at 100°C causes the change of the intensity lines for the Si2p and for the Al2p. This observation confirms the XRD results which indicated the loss of bound molecular waters in the structure of the AlCel2 close to this temperature. After heating at a temperature of 550°C these spectra become almost symmetrical. The ratio of the intensity lines mentioned above for Si is 97/3 and for Al is 97/3, at 550°C. These results suggest that at higher temperatures, the Al cations are mostly located in octahedral sites what confirm conclusions from XRD, Raman and Mössbauer analysis.
The K2p spectrum of initial sample AlCel2 were decomposed into two doublets. The origin of the extra K2p doublet is not well understood. On the basis of the chemical analysis, it might be concluded that it is the result of calcium substitution modifying the chemical environment in the interlayer configuration of the micas. The area ratio for the mentioned above components is 80/20, which results from the ratio of K + and Ca 2+ , respectively. However, the calcium content is too high compared to the results of the chemical analysis. It may suggest that not only Ca 2+ is substituting the K + in interlayer sites but that in these sites also, a small amount of structural water is present. This explanation is supported by K2p spectra obtained at 100°C and 300°C. The intensity of extra doublet decreases during the heating and the area ratio for the doublets of K2p is 88/12 at 100°C and 95/5 at 300°C. This is the result of dehydration of interlayer water from inner surfaces with the increase in temperature. This is also supported by XPS spectra of O1s regions obtained at the studied temperatures.
The Fe 2+ /Fe 3+ ratio, based on the area ratio of Fe2p, was determined to be 18/82, in very good agreement with the chemical analysis and with the Mössbauer results. Also, the area ratio of the peaks connected with Fe 3+ is 62/20 and reflects the M2/M1 ratio of Fe 3+ ions in octahedral positions. The main changes of the Fe2p spectrum are observe at 550°C where all of the doublets observed are connected with Fe 3+ Fig. 7 Plots of temperature dependence of average fraction of Fe species in M1 and M2 sites obtained from Mössbauer measurements for a) AlCel1 and b) AlCel2 and c) integral intensity ratio between bands centred at~590 and~705 cm −1 versus temperature obtained from Raman results. Note: values of fraction of Fe species in octahedral sites were obtained after cooling the sample from high to room temperature and the lines indicate only the trends of the data ions. The two peaks can be associated with Fe in M1 sites and one, with a small contribution, of iron in dimeric pyramid sites. Relationship of these Fe2p peaks with specific Fe sites in structure of aluminoceladonite or more generally in structure of layer silicates was showed first time. The area ratio two first peaks to the last one is 83/17, which is similar to the ratio 87/13 obtained from the Mössbauer studies at this temperature.

Conclusions
We show that all of the measured data show that the aluminoceladonite with a higher aluminum content started to structurally reorganize at a higher temperature range than the less aluminum-rich sample. The oxidation of Fe 2+ to Fe 3+ ions started at~350°C for the Al-richer sample and at 300°C for the other. This process may be associated with, or even initiate, dehydroxylation. Up to a temperature of~500°C, cation migration in the octahedral sheet is not seen. Based on Mössbauer results, visible changes in the coordination numbers of polyhedrons occurred; octahedrally-coordinated Fe 3+ occupying M2 sites becomes 5-coordinated. This temperature can be viewed as the second stage of dehydroxylation. Based on Raman data, the proposed temperature-dependent integral intensity ratio between bands centred at~590 and 750 cm −1 (I 590 / I 705 ) clearly indicates the temperature at which sixcoordinated polyhedra are transformed to five-coordinated. Above~500°C, cation migration from former cis-sites to former trans-sites take place. Changing of coordination number from five into four is associated with the detachment of OH or O and, as a result, dimetric pyramids can be formed. Electronic supplementary material The online version of this article (https://doi.org/10.1007/s00710-021-00753-z) contains supplementary material, which is available to authorized users.