Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 12, Article number: 15934 (2022 ) Cite this article
For the first time, the structural investigation of a Pb-exchanged zeolite (Pb13.4(OH)10Al17.4Si54.6O144 ∙38H2O) with STI framework type, revealed a highly unusual and intriguing sudden volume increase under continuous heating. Understanding the fundamental mechanisms leading to such an unusual behaviour is essential for technological applications and interpretation of chemical bonding in zeolites. The dehydration was tracked in situ from 25 to 450 °C by single crystal X-ray diffraction, infrared and X-ray absorption spectroscopy. Further interpretation of the experimental observations was supported by ab initio molecular dynamics simulations. Initially, Pb-STI unit-cell volume contracts (ΔV = − 3.5%) from 25 to 100 °C. This agrees with the trend observed in STI zeolites. Surprisingly, at 125 °C, the framework expanded (ΔV = + 2%), adopting a configuration, which resembles that of the room temperature structure. Upon heating, the structure loses H2O but no de-hydroxylation occurred. The key mechanism leading to the sudden volume increase was found to be the formation of Pbx(OH)y clusters, which prevent the shrinking of the channels, rupture of the tetrahedral bonds and occlusion of the pores. This zeolite has therefore an increased thermal stability with respect to other STI metal-exchanged zeolites, with important consequences on its applications.
Zeolites are an important class of crystalline porous materials, which occur in nature and are extensively synthesized in laboratory. Thanks to their specific properties, they find applications in several fields like petro-chemical industry, pollution control, environmental remediation, energy storage, etc.1,2.
The practical use of zeolites depends on numerous factors, which are mainly related to their porous structure. The latter consists of a three-dimensional interconnected framework of tetrahedra (TO4) (T mainly Si,Al), which form cavities and channels where a large variety of extraframework (EF) species and H2O are host3,4.
The porous structure is very responsive to changes in temperature and water vapour pressure and can, therefore, easily release or absorb H2O. Of particular interest are the structural modifications occurring as a function of increasing temperature since many zeolites find use only after thermal treatment5,6. The extent of such modifications can vary from minor to severe. For instance, some zeolites experience upon dehydration the rupture of the tetrahedral bonds T-O-T of the framework7 with subsequent collapse of the structure. These changes influence the microporous properties and therefore the performance of the material. Thus, their understanding and ultimately the prediction of the thermal stability is of paramount importance.
In contrast to dense solid phases, most zeolites undergo a negative thermal expansion (NTE) and typically their volume contracts with increasing temperature6. This process is caused by the release of the structural H2O and the subsequent migration of the EF cations toward the wall of the framework in order to reach a more balanced coordination environment. In contrast, positive thermal expansion has been rarely observed8. It must be noticed that a slight expansion of the unit-cell volume at the early stage of heating is commonly due to weakening of hydrogen bonds, but not to H2O loss.
In general, the extent of the transformations and therefore the modification of the structural topologies upon heating, depend on the framework type and the kind of EF species that occupy the pores6,7.
In previous studies, we reported the effect of heavy-metals incorporation on the thermal stability of a microporous zeolite with STI framework type9,10,11. The latter has symmetry Fmmm that can be lowered depending on the type and distribution of the extraframework occupants. The topology at room temperature (phase A) consists of two systems of interconnected channels: The first is composed by 8-membered rings of tetrahedra, running parallel to [001], and the second is formed by ten-membered rings and runs parallel to [100]12. The interconnection of tetrahedra originates the t-sti-1* cavity13, where the EF cations and H2O are hosted. Zeolites with STI framework type experience upon heating severe framework modifications, which lead to: (i) the rupture of specific T-O-T connections of the framework; (ii) the formation of new structural topologies (namely phase B14,15,16,17, B’11, D18, D’10) and, (iii) the occlusion of the ten-membered ring channels with subsequent loss of porosity9,10,18. The aforementioned structural changes lead to significant contraction of the pores accompanied by decrease of the unit-cell volume.
The extraframework cations in STI-type zeolites (Fig. 1) influence: (1) the structural modifications (i.e. phase B, B’, D, D’) that form upon heating, (2) the unit-cell volume trend, and (3) the upper temperature limit of the structural stability.
Typical temperature-dependent behavior and structural transformations observed in zeolites with STI framework type, containing different EF cations. (a) The graph reports the unit-cell volume as a function of increasing temperature for Ag10-, Ca9-, Na9-, and Cd11-stellerite obtained by single crystal X-ray diffraction. The capital letters (A, B, B’, D’, D) refer to the name of the corresponding phases observed as a function of increasing temperature for each metal-exchanged form. The contraction of the unit-cell volume (ΔV) is reported as a percentage with respect to the one measured at room temperature. (b) Representative structures for the A, B, and D-phase topology. The aluminosilicate framework is shown as blue tetrahedra. Cyan, yellow and red tetrahedra represent those tetrahedra involved in the T-O-T breaking process for each phase.
In this study, we report on the behaviour of Pb-exchanged stellerite (Pb-STI) upon heating. We monitor in situ the structural changes with increasing temperature and we demonstrate that the dehydration behaviour as well as the associate structural transformations are anomalous compared to, not only those of STI framework type zeolites, but also to the commonly reported unit-cell volume decrease observed in zeolitic materials. To the best of our knowledge, it is the first time that a sudden volume increase at elevated temperature, is reported for a zeolite structure. This increase is accompanied by a reversal of the structural transformations, which finally result in a high-temperature phase with the same framework configuration as the RT structure. We will show that this phenomenon can be explained by the formation of Pbx(OH)y clusters during the dehydration process. Most important, it leads to increased thermal stability and has therefore important consequences for the applications of these materials.
The unit-cell volume of Pb-STI, with corresponding framework modifications, was determined by in situ Single Crystal X-ray Diffraction (SC-XRD) from room temperature (RT) to 450 °C (Fig. 2). Three different transformation stages can be distinguished:
From RT to 100 °C, the unit-cell volume gradually decreases by 3.5% with respect to that measured at RT19. This change agrees with the common trend observed for STI framework type zeolites (Fig. 1) and can be associated with an initial dehydration process9,10,11. The structural transformation in this temperature range (Fig. 2a–c) is characterized by the shrinking and deformation of the ten-membered ring channels parallel to [100] that become more elliptical. This structural modification is at 50 °C accompanied, by the lowering of the crystal symmetry from Fmmm to A2/m (Table 1).
From 100 to 150 °C, in contrast to the common behaviours of STI type zeolites, the unit-cell volume increases by approximately 2%. The space group at 125 °C turns into orthorhombic (space group Amma) and the channels expand, adopting the original roundish shape, characteristic to the RT structure (Table 1, Fig. 2d).
From 175 to 450 °C, the unit-cell volume does not vary significantly. The structure gradually transforms from Amma to Fmmm space group (Table 1, Fig. 2e).
(a) Evolution of the unit-cell volume of Pb-STI as a function of temperature obtained by SC-XRD. (b–e) Corresponding structural changes in STI framework of Pb-stellerite observed in situ at different temperatures.
Thus, the thermal behaviour of Pb-STI is characterized by an initial volume contraction and a subsequent expansion, which leads to a structural topology (phase A) equal to that observed at RT. This mechanism is anomalous not only for STI framework-type zeolites, but in general for zeotype materials, which are normally characterized by a negative thermal expansion. Most impressive is the sudden volume increase at a given temperature accompanied by the expansion of the contracted framework, which is usually observed only after rehydration. Moreover, the thermogravimetric analysis19 indicated a gradual and continuous release of (most probably) H2O upon heating that should be associated with a channel contraction and corresponding decrease of the unit-cell volume.
Because of the strong disorder of the EF occupants (Pb2+, H2O and OH−), their positions as well as the eventual diffusion and release of the EF species could not be unequivocally monitored by SC-XRD. From the structural refinements, it was possible, however, to follow the variation of the electron-density distribution associated with the EF content. The latter was refined as several low occupied sites, close to each other. At RT, the EF species were mainly concentrated at the centre of the ten-membered ring channels and in the middle of the eight-membered ring window between two t-sti-1* cavities19. With the increase of temperature and subsequent shrinking of these channels, Pb and eventually OH groups came closer to the wall of the framework. When the structure expanded and the channels adopted again a roundish configuration, the EF species remained close to the framework wall (Fig. 3).
Polyhedral representation of the Pb-STI structure at RT, 50 °C and 400 °C, refined from SC-XRD data. Partially colored dark-grey spheres represent partially occupied Pb sites. Red arrows indicate the direction of EF cations movement upon heating. For sake of clarity, only the Pb sites are shown.
Spectroscopic analyses were performed to obtain more information about the dehydration process. A distinction between H2O and OH− by infrared (IR) analysis was complicated by the fact that these molecules vibrate at similar frequencies. Vibrational Density of States (VDOS) curves were calculated for hydrogen atoms (H) based on the DFT-equilibrated Pb-stellerite structure at room temperature from our previous study19. Different bands in the IR spectra were assigned to specific molecular groups (Supplementary material S1): vibrations of H2O are found as a sharp peak between 1500 and 1800 cm−1 and a broader peak between 3000 and 3800 cm−1, while only one rather sharp vibrational peak of OH− is found between 3500 and 3800 cm−1 (Fig. 4a,b).
(a) VDOS curves calculated for hydrogen belonging to OH− bonded to 1 Pb (orange), hydrogen belonging to OH− bonded to 2 Pb (red), hydrogen belonging to “free” H2O molecules (light blue) and hydrogen belonging to H2O bonded to 1 Pb (dark blue), based on the theoretical structural model of Pb-STI at RT19. (b) VDOS curves calculated for hydrogen belonging to H2O (blue) and hydrogen belonging to OH− groups (red). (c) Selected IR absorption spectra of Pb-STI at RT (blue), after ex situ thermal treatment (red) and after rehydration (purple).
Experimentally measured IR absorption spectra of Pb-stellerite are reported in Fig. 4c and Fig. S1. The region of the H2O bending vibrations (1500–1800 cm−1) is strongly influenced by the background measurements due to thickness variations in the quartz glass capillary (quartz glass has an IR absorption peak at 1635 cm−1) and was therefore not included in the interpretation. Instead, another absorption peak from 5100 to 5300 cm−1 was detected in the spectrum that results from the bend and stretch combination mode of the H2O molecule20. This peak is present in the spectra obtained at RT (Fig. 4c, Fig. S1a), indicating that the structure is hydrated.
IR spectra collected after the ex situ thermal treatment up to 430 °C are displayed in Fig. 4c, Fig. S1b. The peak between 5100 and 5300 cm−1, as well as the left shoulder of the broad, flat absorption band between ca. 2700 and 3700 observed at RT, disappeared. The only remaining absorption band in the spectra after the thermal treatment is sharp and at ca. 3500 cm−1, similar to the calculated VDOS curve for H belonging to OH− groups (Fig. 4b). Thus, it can be concluded that Pb-STI, at high temperatures, contains OH− groups but is H2O free.
After the ex situ thermal treatment, Pb-STI crystals were exposed to high humidity conditions (relative humidity > 90%) in order to allow reabsorption of H2O. IR absorption curves of the crystals obtained after 15 days equilibration time are shown in Fig. 4c, Fig. S1c. The absorption on the left shoulder of the broad peak between 2700 and 3700 cm−1 increased and so did the absorption of the H2O bend and stretch combination peak, illustrated by purple-coloured arrows in Fig. 4c. This indicates that the structure of Pb-STI is able to reabsorb H2O and that the dehydration process is reversible.
The reversibility of the dehydration process was also investigated by SC-XRD. After an equilibration time of 15 days under high humidity conditions (relative humidity > 90%), new data were collected at RT on a previously dehydrated crystal. The rehydrated structure has Fmmm symmetry and lattice parameters a = 13.6119(7), b = 18.1834(9), c = 17.8460 (10) Å and V = 4417.1(4) Å3. The latter are very similar to those of the RT structure. The deviation of the unit cell volume with respect to the RT Pb-STI is smaller than 0.34%. The positions of the EF sites do not differ from those in RT Pb-STI19, despite the disorder (Fig. S2, Table S1).
The pronounced similarity between RT Pb-STI and rehydrated Pb-STI in terms of both structural and spectroscopy analyses consistently indicates that Pb-STI de- and rehydration as well as its structural modifications upon heating are reversible processes.
The experimental results showed that the fully exchanged Pb-STI has a higher thermal stability with respect to the pristine material (Ca-STI9) and the other metal-exchanged forms (Na-STI9, Ag-STI10, Cd-STI11). In particular, no formation of the structural modifications B, B’, C and D, D’ occurs upon heating and the framework does not experience the breaking of the tetrahedral bonds T-O-T. Consequently, it does not lose its microporous structure at high temperatures. The most surprising observation is the expansion of the zeolitic framework after the initial contraction, which is induced by dehydration. According to TGA19 and IR analysis, H2O is lost as a function of increasing temperature whereas, as demonstrated by IR spectra collected at higher temperatures (430 °C), the OH− groups are retained in the dehydrated structure.
The interpretation of this unique behaviour was complicated by the disorder of the EF species. Thus, different hypotheses were tested by means of MD simulations. A detailed description of the set up and the length of trajectories for each structural model are given in the Methods session.
According to SC-XRD data, two main structural configurations form upon heating: (i) a contracted STI topology, with space group A2/m characterized by elliptical channels, and (ii) an “expanded” structure, compared to that observed at 75 °C, with Fmmm space group and roundish channels, which resemble those of the RT phase.
The first structural configuration was well reproduced by the theoretical model containing 12 Pb p.f.u. and 0.8 H2O per Pb, simulated at 75 °C (Pb12-75-10 W) (Fig. S3, Table S2). Thus, we can assume that at 75 °C, Pb-STI structure is partially dehydrated.
To reveal the mechanism, which brings to the formation of the second structural topology, and to understand the evolution of the EF species upon heating, three potential scenarios (M1, M2, M3) were considered:
M1: Dehydration of the structure without any further reaction involved. In this case, the reference model consisted of a fully dehydrated Pb-STI (Pb12OH8Si56Al16O144).
M2: Partial H2O dissociation at elevated temperatures. Dissociation of H2O molecules based on the hydrolysis reaction H2O ⇌ OH− + H+ has been reported in literature for zeolite L21.
M3: Oxidation of some Pb2+ ions to Pb4+. The reaction Pb2+ + 2 H2O ⇌ Pb4+ + 2OH− + H2(g) was hypothesized for Pb-exchanged zeolite A treated at high temperatures22. The produced OH− groups coordinate Pb and the H2, a gas, is assumed to leave the crystal structure.
The analysis of the tested models showed that the dehydrated Pb-STI (M1) has an extremely contracted unit- cell volume with the resulting structural configuration similar to that of the other dehydrated forms of STI zeolites9,10,11 (Fig. S4a, Table S3). Therefore, this model was not further considered. The “hydrolysis” hypothesis (M2) was ruled out based on the obtained equilibrated structures (Fig. S4b, Table S3), the configuration of which has also contracted channels, not in agreement with the experimental structure measured at 400 °C. Moreover, 60% of OH− groups re-protonated, indicating that the deprotonated form is less stable than the H2O molecules. The theoretical model that better matched the experimental data, in terms of unit-cell volume and framework configuration, was the M3 scenario, i.e. a Pb-STI with 17% oxidized Pb4+ (Fig. S4c, Table S3). However, XANES spectra collected in situ as a function of temperature, did not show any evidence of a change of the oxidation state of Pb2+ (Fig. S4d), not supporting this hypothesis.
Finally, a series of structural models (M4) containing different PbxOy and Pbx(OH)y clusters (1 < x,y < 4) was tested (Table 2). Among the structural configurations, the one that better reproduced the experimental data was a dehydrated model with a composition Pb12OH8Si56Al16O144 where 33% and 16% of Pb formed Pb4OH4 and Pb2OH2, respectively (reference model name Pb12-400-C4). The framework of the equilibrated structure well reproduced that of the refined Pb-stellerite at 400 °C, characterized by the roundish channels (Fig. 5a) parallel to [100]. The size of the horizontal and vertical maximum length was 7.25 and 8.68 Å, respectively, in good agreement with the experimental values (7.01 and 8.99 Å).
(a) Comparison between the experimental structure measured at 400 °C by SC-XRD and the simulated one corresponding to the model Pb12-400-C4. The framework is displayed by polyhedral representation. (b–e) Fragment of the simulated structure (average atomic coordinates) showing the different arrangement of PbOH species in the t-sti-1* cavities of Pb-STI. Dark grey spheres represent Pb atoms. Oxygen and hydrogen are depicted as red and white spheres, respectively. The framework is shown as blue sticks.
The simulation suggests the existence of different [Pbx(OH)z]y+ complexes in the zeolitic cavities: i) [Pbx(OH)x]y+ (x = 2, 4) clusters in the middle of the cavities; ii) as single cation Pb2+ bonded to the oxygen-atoms of the framework, and iii) as [Pb2(OH)]3+ species bonded to the oxygen atoms of the framework. The Pb atoms in the Pb4(OH)4 clusters adopt an almost tetrahedral shape (Fig. 5b) with average Pb-OH bond distances equal to 2.35 Å. Pb atoms are on average at 3.82 Å far from each other, with the shortest Pb-Pb contact found at 3.64 Å. Similarly, the Pb2(OH)2 clusters are found close to the wall of the framework, where they bond with the tetrahedral oxygen atoms (Fig. 5c). The corresponding Pb-OH (average = 2.27 Å) and Pb–O distances (average = 2.63 Å) are shorter than those in Pb4(OH)4.
The single Pb2+ cations are mainly found at the wall of the six-membered ring apertures of the t-sti-1* cavity, with average Pb–O distances = 2.5 Å (Fig. 5d). Finally, the [Pb2(OH)]+3 species consist of two Pb atoms, which share an OH group (with average Pb-OH distances = 2.26 Å) and bond to further oxygen atoms of the framework at longer distances (Pb–O = 2.64 Å) (Fig. 5e).
Our structural characterisation of Pb-STI reveals for the first time a highly unusual thermal behaviour in a microporous zeolite. The mechanism is unique in two ways:
The sudden volume increase at ca. 100 °C is accompanied by the reversal of the structural modification. The structure transforms back to its original framework configuration. Commonly, a slight expansion of the unit-cell volume is observed in some zeolites, but only after complete dehydration, before the structural collapse23,24. Such minor volume changes do not lead however to the transformation back to the RT structure, unless rehydration occurs.
The expansion apparently prevents an immediate structural collapse commonly observed in zeolites upon dehydration. The thermal stability of Pb-STI increased significantly with respect the other STI zeolites. The diffraction data collected at high temperature indicate good crystallinity of the sample up to 450 °C.
This unprecedented behaviour is of high relevance because Pb-STI is the only exchanged form maintaining the RT framework configuration (and therefore its microporous properties) up to 450 °C. To the best of our knowledge, this is the first time that such an anomalous thermal behaviour is reported in zeolites.
The sudden increase of the unit-cell volume is associated with the formation of Pb(OH) clusters, which prevent the further shrinking of the channels. The MD simulations demonstrated that a dehydrated Pb-STI without clustering of Pb2+ and OH−, results in a contracted topology. Therefore, the formation of Pb4(OH)4 clusters must be the key mechanism, which leads to the channel expansion between 125 and 150 °C.
The unique dehydration behaviour of Pb-STI can be explained according to the following process: at first, the water loss is accompanied by the unit-cell volume contraction, in line with the general negative thermal expansion observed in zeolitic materials. At 75 °C, the structure is partially hydrated and between 100 and 125 °C, in correspondence of the increase of the unit-cell volume, the cluster formation starts. Residual water might be retained up to 250–270 °C, and in a subsequent step completely released, as indicated by the slight unit-cell volume decrease between 250 and 450 °C and by TGA curve19 (Fig. S5). The different thermal behavior of STI zeolites exchanged with different EF cations (Ca, Na, Ag, Cd) was related not only to the size and charge of the counter cation9,10,11, but also to the electronic structure and strength of chemical bonding of the EF cations with H2O and framework oxygen atoms11. In light of the results obtained on Pb-STI, we can conclude that the nature of the cationic species (cation, hydroxyl, etc.) which are up taken during the exchange process as well those forming upon heating, have an impact on the dehydration path. In Pb-STI, the clustering of Pb(OH)+ species during the dehydration prevents the collapse of the structure because no enough space is left in the channels to accommodate the “flipping” of the tetrahedra, which would be involved in the breaking of T–O–T bonds, typical of zeolites with STI framework type.
The formation of metal clusters in zeolites under thermal treatment was at first hypothesized by Uytterhoeven et al.25. The authors stated that, “H2O are dissociated by the strong electrostatic fields of the Me2+ ions” leading to the formation of metal clusters of the type Me+–O–Me+ or Me+–OH. Baekelant et al.26 used the same hypothesis to explain the behaviour of Pb-LTA zeolite upon heating. They reported the formation of Pb4(OH2)x (x = 2–4) and subsequently, of Pb4Oz (z = 2–4) clusters upon increasing temperature, suggesting that the OH− groups dehydroxylate and transform to extraframework oxygen. Pb4(O,OH)4 clusters were also found in the mineral maricopaite, the first natural zeolite with Pb as a dominant extraframework cation27. This natural zeolite contains, at room temperature, H2O and OH− groups. Pb atoms form Pb4 tetrahedra capped by OH or H2O. At date, no investigation of its dehydration behaviour has been conducted.
In our study, the process leading to the clusters formation is different from that described by Uytterhoeven et al.25 and Baekelant et al.26 for Pb-LTA zeolite. The PbxOHx clusters do not form because of water dissociation, since the OH species are already present in the hydrated structure at RT, i.e. the STI zeolite takes up Pb in form of Pb2+ and Pb(OH)+ species19 during the exchange procedure. In addition, no dehydroxylation takes place up to 450 °C, as demonstrated by the IR spectra of the thermally treated sample. The latter process would indeed imply the presence of protons attached to framework oxygen atoms. Last hypothesis was tested by the MD simulations (Pb12-400-C2, Pb12-400-C4O) and could be rejected on the basis of the obtained results.
These results suggest that the incorporation of Pb in form of (PbOH)+ species into zeolites can increase the thermal stability and widen the range of temperatures for which the structural changes are reversible. Further studies must be conducted to ascertain if this unusual trend can be generalized to all Pb-zeolites containing Pb(OH)+ and H2O at RT, and if and how this mechanism might be tuned for specific purposes.
The sample under investigation was a Pb-stellerite (Pb-STI) produced and characterized at room temperature in our previous study19. This Pb-STI was obtained by cation exchange of a natural stellerite28, previously exchanged with Na (Na-STI), in 0.5 M Pb-acetate solution for 4 weeks at 95 °C. The electron microprobe analyses coupled with X-ray Absorption Spectroscopy19 indicated the presence of [Pb(OH)]n+ species within the channels and resulted in the chemical composition: Pb13.4(OH)10Al17.4Si54.6O144 ∙38H2O.
A single crystal with dimensions 0.070 × 0.100 × 0.200 mm was glued on the tip of a glass fibre and mounted on a goniometer head. This was the same crystal used for the RT data collection in our previous study19. The diffraction data were collected on a Bruker APEX II diffractometer equipped with a CCD detector, MoKα fine focus tube and an in house N2 blower. The high temperature (HT) data were obtained by heating the crystal stepwise from 50 to 450 °C (maximum allowed temperature of our equipment), in steps of 25 °C. Before each data collection, the crystal was equilibrated for at least 30 min. The same experimental set up was used for the investigation of the thermal stability of the other cationic forms of stellerite (i.e. Ca-, Na-, Ag-, and Cd-STI)9,10,11 and is regarded as dry experimental conditions.
The data were integrated and corrected for absorption using the software package Apex3.v.2018.7-229. Structure solution (ShelxS30) for the data set at 50, 75 and 100 °C pointed to the monoclinic space group C2/m, transformed to the non-standard setting A2/m for better comparison with literature data. From 125 to 150 °C, the structures were solved and refined (ShelxL31) in the Amma space group, and from 200 to 450 °C in Fmmm.
The reversibility of the dehydration process (rehydration) was also investigated by SC-XRD. Another Pb-STI crystal was selected, because the one used during the first set of high-temperature measurements was lost after the experiment. After a stepwise thermal treatment up to 200 °C for approximately 6 h, the crystal was exposed to high humidity conditions (relative humidity > 90%) for 15 days. Subsequently, diffraction data were collected at RT with data collection strategy similar to that used for the first crystal.
Crystal data and refined parameters of the Pb-STI measured at 50 (PbSTI 50), 125 (PbSTI 125), and 400 °C (PbSTI 400) are reported in Table 1.
Structural data of Pb-STI obtained at different temperatures have been deposited in form of Crystallographic Information File (Cif). All structural drawings were produced with the software VESTA32.
Ab initio simulations were performed with the CP2K package33. The interatomic forces were obtained based on the Density Functional Theory (DFT)34. The electron exchange and correlation energies were approximated with the Perdew-Burke-Ernzerhof (PBE) functional35. Dispersion correction was applied using the semi-empirical pair potential type DFT-D336. The Gaussian and plane waves basis sets were used to represent electron wave function and the electron density, respectively37. Valence electron wave functions were described by the Gaussian short range Double-Zeta Valence Polarized basis set (DZVP-MOLOPT)38. Core electrons were represented by dual space norm conserving pseudopotentials39. MD simulations were performed in an isothermal-isobaric (NPT) ensemble with a flexible cell. The temperature was controlled by the Canonical Sampling Velocity Rescaling (CSVR) thermostat40. The integration time step was set to 0.5 fs. Initial structures for MD simulations were obtained by geometry optimisation of the hypothetical structural models.
Different structural models of Pb-stellerite in terms of H2O-, Pb-content and temperature were tested. Here we report those that are discussed in the Results session.
All structures consisted of a fully flexible 2 × 1 × 1 Pb-STI supercell. The starting coordinates of the framework were taken from the refined structures obtained from SC-XRD data. The following configurations were considered:
Partially dehydrated Pb-stellerite at 75 °C (Pb12-75-10 W) (Table S2); The initial coordinates were obtained starting from the equilibrated structure of RT Pb-stellerite, with chemical composition Pb12(OH)8Si56Al16O144 ∙34H2O19 by removing 14 H2O.
Dehydrated Pb-stellerite at 400 °C (M1: Pb12-400-0W) (Table S3); The dehydrated structure, with chemical composition Pb12(OH)8Si56Al16O144 was simulated by using two different approaches. At first, the water was removed stepwise by increasing the temperature and the MD trajectories were run at 400 °C. A second approach consisted in using the atomic coordinates of the framework of the refined structure at 400 °C, and distributing the Pb and OH- groups as indicated by the electron density distribution of the XRD data, i.e. close to the wall of the framework. No significant differences were noticed between the two models; therefore, we reported here the results from the second approach.
Pb-stellerite with hydrolysed water 400 °C (M2: Pb12-400-hy) (Table S3); This model had chemical composition Pb12(OH)28H20Si56Al16O144, and corresponds to a structure where 20 H2O (estimated as the number of H2O in the partially hydrated structure) underwent the hydrolysis reaction: H2O ⇌ OH− + H+. The structure contains additional 20 OH− groups p.f.u., bonded to Pb, and 20 H+ p.f.u., attached to framework oxygen atoms.
Pb-stellerite with oxidized Pb4+ at 400 °C (M3: Pb12-400-10Pb4 and Pb12-400-4Pb4) (Table S3). In this case, the oxidation of Pb2+ to Pb4+ was hypothesized. Two MD simulations were run at 400 °C, one assuming 5 out of 12 Pb p.f.u. (42%) and the other one 2 out of 12 Pb p.f.u. (17%) in the oxidized form. Additional 10 OH− groups p.f.u. were added to the structural model hypothesizing 42% of all Pb to be oxidized whereas 4 additional OH− groups were added assuming an oxidation of 17% of all Pb ions.
Dehydrated Pb-stellerite containing Pbx(OH)y and PbxOy clusters, with initial values x = y = 2, 4 (M4: Pb12-400-Cx) (Table 2). The clusters were distributed in the t-sti-1* cavities and the remaining Pb2+ and OH− species were located, in agreement with the electron density distribution revealed by XRD, close to the wall of the framework.
Vibrational frequencies of hydrogen (H) atom in different structural environment were derived from MD simulations. The analysis of the 15 ps long equilibrated trajectory of the RT Pb-STI19 was performed with the program TRAVIS41. The total and partial VDOS was calculated for: (1) all H atoms, (2) H atoms belonging to H2O (H(water)) and H atoms belonging to OH− groups (H(hydroxyl)) and (3) H atoms belonging to H2O molecules that are not bonded to Pb (distance Pb–O > 3 Å) (H(water)0), H atoms belonging to H2O molecules that are bonded to Pb (distance Pb–O < 3 Å) (H(water)1), H atoms belonging to OH− groups that are bonded to 1 Pb (H(hydroxyl)1) and H atoms belonging to OH− groups that are bonded to 2 Pb (H(hydroxyl)2). The spectrum was calculated up to wavenumber 5000 cm−1 and the correlation depth of the ACF (autocorrelation function) was set to 256 for the calculation of all curves.
Transmission IR spectroscopy analyses were carried out on a Bruker Tensor II spectrometer equipped with a Hyperion 3000 FT-IR microscope and an MCT detector which is cooled with liquid nitrogen. The side length of the aperture square was approximately 40 µm. Each sample and background measurement consisted of 64 scans. The resolution was set to 4 cm−1.
Spectra were collected for several Pb-STI crystals, placed inside a Perspex chamber. Moreover, Pb-STI crystals, with dimensions ranging from ca. 20–60 µm, were placed inside a quartz glass capillary with 100 µm diameter and a wall thickness of 10–20 µm. Spectra of several crystals, inserted in the capillary, were measured at RT. The background spectrum was measured on the empty glass capillary. Afterwards, the crystals were ex situ thermally treated: the temperature was increased from RT to 430 °C with a heating rate of 40 °C/h. Then, it was kept constant for 8 h at 430 °C, before cooling down slowly. As soon as a temperature of 250 °C was reached, the capillary with the crystals was removed from the oven and the capillary was closed immediately with wax to prevent rehydration of the crystals. After cooling down to RT, the IR spectra of several crystals were measured again. The capillary was then opened and exposed to high humidity conditions (relative humidity higher than 90%). After 15 days, IR absorption spectra of several rehydrated crystals were measured.
Pb-STI crystals were gently grounded and the obtained powder was investigated by Pb-L3 X-ray Absorption Near-Edge Structure (XANES) spectroscopy at the Rossendorf Beamline (BM2042) at the European Synchrotron Radiation Facility (Grenoble, France). A double crystal monochromator mounted with a Si (111) crystal was used. For each XANES measurement, the spectra of a Pb reference foil were systematically collected for post-treatment of the energy calibration. The sample was heated using a hot gas generator from Cyberstar controlled by a Eurotherm with a precision of ± 1 °C. Data were collected in situ from 25 to 400 °C.
The datasets generated during and/or analysed during the current study are available in the Crystallography Open Database (COD) repository (entries number: 3000391, 3000392, 300093, 300094) and/or included in the article as supplementary information files.
Millini, R., & Bellussi, G. Zeolite science and perspectives in Zeolites in Catalysis: Properties and Applications (eds. Čejka, J., Morris, R. E., Nachtigall, P.) 1–36 (The royal Society of Chemistry 2017).
Bish, D. L., & Ming, D. W. Natural zeolites: Occurrence, properties and applications. Rev. Mineral. Geochem. 45, 654.
Smeets, S., & Zou, X. Zeolite structures in Zeolites in Catalysis: Properties and Applications (eds. Čejka, J., Morris, R. E., Nachtigall, P.) 37–72 (The royal Society of Chemistry 2017).
Armbruster, T., Gunter, M. Crystal structures of natural zeolites in Natural zeolites: occurrence, properties and applications, Reviews in Mineralogy and Geochemistry, Vol. 45, (eds. Bish, D. L., Ming, D. W.) 1–67 (Mineralogical Society of America, 2001).
Bish, D. L., Carey, J. W. Thermal properties of natural zeolites in Natural zeolites: occurrence, properties and applications, Reviews in Mineralogy and Geochemistry, Vol. 45, (eds. Bish, D. L., Ming, D. W.) 403–452 (Mineralogical Society of America, 2001).
Cruciani, G. Zeolites upon heating: Factors governing their thermal stability and structural changes. J. Phys. Chem. Solids 67, 1973–1994 (2006).
Alberti, A. & Martucci, A. Reconstructive phase transitions in microporous materials: Rules and factors affecting them. Microporous and Mesoporous Mater. 141, 192–198 (2011).
Carey, T. et al. The effect of extra framework species on the intrinsic negative thermal expansion property of zeolites with the LTA topology. Chem. Commun. 48, 5829–5831 (2012).
Cametti, G., Fisch, M. & Armbruster, T. Thermal behavior of stilbite and stellerite revisited and dehydration of their Na-exchanged forms: Considerations on the memory effect of the STI framework type. Micropor. Mesopor. Mat. 253, 239–250 (2017).
Cametti, G., Scheinost, A. C., Giordani, M. & Churakov, S. V. Framework modifications and dehydration path of a Ag+-modified zeolite with STI framework type. J. Phys. Chem. C 123, 13651–13663 (2019).
Cametti, G., Scheinost, A. C. & Churakov, S. V. Structural modifications and thermal stability of Cd2+-exchanged stellerite, a zeolite with STI framework type. J. Phys. Chem. C 123, 25236–25245 (2019).
Baerlocher, C., Meier, W. M., Olson, D. H. Atlas of zeolite framework types, structure commission of the IZA, fifth ed., Elsevier (2001).
Smith, V. J. Topochemistry of zeolites and related materials. 1. Topology and geometry. Chem. Rev. 88, 149–182 (1988).
Alberti, A., Rinaldi, R. & Vezzalini, G. Dynamics of dehydration in stilbite-type structures Stellerite phase B. Phys. Chem. Miner. 2, 365–375 (1978).
Cruciani, G., Artioli, G., Gualtieri, A., Ståhl, K. & Hanson, J. C. Dehydration dynamics of stilbite using synchrotron X-ray powder diffraction. Am. Mineral. 82, 729–739 (1997).
Ori, S., Mazzucato, E. & Vezzalini, G. Dehydration dynamics of barreirte: An in situ synchrotron XRPD study. Am. Mineral. 94, 64–73 (2009).
Arletti, R., Mazzucato, E. & Vezzalini, G. Influence of dehydration kinetics on T-O-T bridge breaking in zeolites with framework type STI: The case of stellerite. Am. Mineral. 91, 628–634 (2006).
Sacerdoti, M. The crystal structure of zeolite barrerite dehydrated in air at 400–450 °C. Micorpor. Mesopor. Mat. 102, 299–303 (2007).
Roos, D. P., Scheinost, A. C., Churakov, S. V., Nagashima, M. & Cametti, G. On the nature of Pb species in Pb-exchanged zeolite stellerite (STI): A combined experimental and theoretical study. Micropor. Mesopor. Mater. 327, 111444 (2021).
Aines, R., D., Rossman, G. R. Water in minerals? A peak in the infrared. J. Geophys. Res. 89, 4059–4071 (1984).
Albuquerque, R. Q. & Calzaferri, G. Proton activity inside the channels of zeolite L. Chem. A Eur. J. 13, 8939–8952 (2007).
Ronay, C., & Seff, K. Crystal structure of Pb6-A and Pb9(OH)8(H2O)3-A. Zeolite A ion exchanged with Pb2+ at pH 4.3 and 6.0 and evacuated. J. Phys. Chem. 89, 1965–1970 (1985).
Wadoski-Romeijn, E. & Armbruster, T. Topotactic transformation and dehydration of the zeolite gismondine to a novel Ca feldspar structure. Am. Mineral. 98, 1988–1997 (2013).
Cametti, G., Danisi, R. M., Armbruster, T. & Nagashima, M. De- and re-hydration of scolecite revisited: An in situ single-crystal X-ray study under low and high humidity conditions. Micropor. Mesopor. Mater. 208, 171–180 (2015).
Uytterhoeven, J. B., & Schoonheydt, R. Studies of the hydrogen held by solids. XVI Infrared spectroscopy of X- and Y-type zeolites containing univalent and divalent cations. J. Catal. 13, 425–434 (1969).
Baekelant, W. et al. Confinement of highly luminescent lead clusters in zeolite A. J. Phys. Chem. C 122, 13953–13961 (2018).
Rouse, R. C. & Peacor, D. R. Maricopaite, an unusual lead calcium zeolite with an interrupted mordenite-like framework and intrachannel Pb4 tetrahedral clusters. Am. Mineral. 79, 175–184 (1994).
Armbruster, T. et al. The zeolite, fluorite, quartz assemblage of the fissures at Gibelsbach, Fiesch (Valais, Switzerland): Crystal chemistry, REE patterns and genetic speculations. Mineral. Petrogr. Mitt. 76, 131–146 (1996).
Sheldrick, G. M. A short history of SHELX. Acta Crystallogr. A 64, 112–122 (2008).
Sheldrick, G. M. Crystal structure refinement with SHELXL. Acta Crystallogr. C 71, 3–8 (2015).
Momma, K. & Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Crystallogr. 44, 1272–1276 (2011).
Kühne, T. et al. CP2K: An electronic structure and molecular dynamics software package—Quickstep: Efficient and accurate electronic structure calculations. J. Chem. Phys. 152, 194103 (2020).
Kohn, W. & Sham, L. J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 140, 1133–1138 (1965).
ADS MathSciNet Article Google Scholar
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).
Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 132, 1–19 (2010).
VandeVondele, J. et al. Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach. Comput. Phys. Comm. 167, 103–128 (2005).
VandeVondele, J. & Hutter, J. Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. J. Chem. Phys. 127, 114105 (2007).
Hartwigsen, C., Goedecker, S. & Hutter, J. Relativistic separable dual-space Gaussian pseudopotentials from H to Rn. Phys. Rev. 58, 3641–3662 (1998).
Bussi, G., Donadio, D. & Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 126, 014101 (2007).
Brehm, M. & Kirchner, B. TRAVIS - A free analyzer and visualizer for monte carlo and molecular dynamics trajectories. J. Chem. Inf. Model. 51, 2007–2023 (2011).
Scheinost, A. C. et al. ROBL-II at ESRF: A synchrotron toolbox for actinide research. J. Synchrotron Radiat. 28, 333–349 (2021).
Jorg Hermann is acknowledged for his assistance during the IR measurements. We are thankful to Thomas Armbruster for fruitful comments and discussion. We acknowledge access to the Swiss National Supercomputing Center (CSCS) and UBELIX HPC cluster at the University of Bern. GC acknowledges the Swiss National Science Foundation (SNSF) for the grant nr. PZ00P2 173997.
Institute of Geological Sciences, Bern University, Baltzerstrasse 1+3, 3012, Bern, Switzerland
Georgia Cametti, Diana P. Roos & Sergey V. Churakov
The Rossendorf Beamline at the European Synchrotron Radiation Facility (ESRF), Avenue des Martyrs 71, 38043, Grenoble, France
Damien Prieur & Andreas C. Scheinost
Institute of Resource Ecology, Helmholtz Zentrum Dresden Rossendorf, Bautzner Landstrasse 400, 01328, Dresden, Germany
Damien Prieur & Andreas C. Scheinost
Paul Scherrer Institut, Forschungstrasse 111, 5232, Villingen PSI, Switzerland
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
G.C. conceived and supervised the study, performed the XRD experiments and MD simulations, analysed the data and wrote the manuscript. D.R. performed the XRD and IR experiments, set up the MD simulations, wrote the manuscript. S.C. supervised and validated the MD simulations, reviewed and edited the manuscript. D.P. and A.S. performed the XANES measurements, analysed the data and reviewed the manuscript.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Cametti, G., Roos, D.P., Prieur, D. et al. Pbx(OH)y cluster formation and anomalous thermal behaviour in STI framework-type zeolites. Sci Rep 12, 15934 (2022). https://doi.org/10.1038/s41598-022-20317-1
DOI: https://doi.org/10.1038/s41598-022-20317-1
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.