Michael
Fischer
ab
aFachgebiet Kristallographie, Fachbereich Geowissenschaften, Universität Bremen, Klagenfurter Straße, 28359 Bremen, Germany
bMAPEX Center for Materials and Processes, Universität Bremen, Germany. E-mail: michael.fischer@uni-bremen.de
First published on 17th May 2016
Porous aluminophosphates (AlPOs) and silicoaluminophosphates (SAPOs) with zeolite-like structures have received considerable attention as potential adsorbents for heat transformation applications using water adsorption/desorption cycles. Since a detailed experimental characterisation of the water adsorption properties has only been performed for some of these materials, such as AlPO-18 (AEI topology) and SAPO-34 (CHA topology), more systematic insights regarding the influence of the pore topology and (for SAPOs) the arrangement of the framework protons on the affinity towards water are lacking. To study the relationships between structure and properties in more detail, the interaction of water with six structurally different AlPOs (with AEI, AFX, CHA, ERI, GIS, RHO topologies) and their SAPO analogues was investigated using dispersion-corrected density-functional theory (DFT-D) calculations. Different possible locations of silicon atoms and charge-balancing protons were considered for the SAPO systems. The calculations for SAPOs at low water loadings (one H2O molecule per framework proton) revealed that the interaction energies exhibit a considerable variation, ranging from −75 to −100 kJ mol−1 (per water molecule). The differences in interaction energy were rationalised with the different structural environment of the framework protons at which the water molecules are adsorbed. At high water uptakes (near saturation), interaction energies in the range of −65 kJ mol−1 were obtained for all AlPOs, and there was no evidence for a marked influence of pore size and/or topology on the interaction strength. The interaction of water with SAPOs was found to be approximately 5 kJ mol−1 stronger than for AlPOs due to an increased contribution of electrostatic interactions. An analysis of the structural changes upon water adsorption revealed striking differences between the distinct topologies, with the materials with GIS and RHO topologies being distorted much more drastically than the systems based on double six-ring (d6r) units. Moreover, the direct coordination of water molecules to framework aluminium atoms occurs more frequently in these materials, an observation that points towards a reduced structural stability upon hydration.
AlPOs and SAPOs exhibit a significant structural variability: a review published in 2010 lists more than 25 AlPOs and 16 SAPOs with distinct topologies,7 and ongoing research efforts have led to the successful synthesis and structural characterisation of several new systems in the last couple of years.8–10 While some of the topologies of AlPOs/SAPOs are also found in zeolites, such as the chabazite (CHA) topology of AlPO-34/SAPO-34, others are unique to aluminophosphate-based systems, for example the AEI topology of AlPO-18/SAPO-18. In the absence of structural defects, AlPOs (composition AlPO4) correspond to a perfect alternating arrangement of Al and P atoms at the tetrahedrally coordinated sites (T sites). The tetrahedral framework of SAPOs has the general composition (SixAlyPz)O2, where x is typically in the range of 0.02 to 0.2. If only small amounts of silicon are incorporated in the structure, the Si atoms tend to replace phosphorus atoms at isolated T sites. In this case, y = 0.5 and z = 0.5 − x. The negative framework charge that arises from the replacement of P5+ by Si4+ is compensated by framework protons. At higher Si contents, silicon atoms may aggregate in larger assemblies (“silicon islands”), rather than forming isolated Si sites, and the occurrence of such heterogeneities is directly linked to the synthesis conditions.11
Sparked by the potential use of AlPOs and SAPOs in heat transformation applications, a number of researchers have investigated the adsorption of water in these systems. For example, water adsorption experiments using powder samples were performed for AlPO-5 (AFI topology),12,13 AlPO-17 (ERI),13 AlPO-18 (AEI),12–16 SAPO-34 (CHA),6,12–14,16 and a triclinically distorted CHA-type system termed AlPO-tric.16 These materials typically exhibit S-shaped water adsorption isotherms and isobars, a feature that is attractive because a large loading spread can be reached upon a moderate change in pressure and/or temperature (the loading spread corresponds to the difference in water uptake between adsorption and desorption conditions; a larger loading spread leads to a higher attainable energy density). The most promising systems, for example AlPO-18 and SAPO-34, exhibit a loading spread in the range of 300 g kg−1 when assuming H2O adsorption at 40 °C and desorption in the range of 100 °C, conditions that are relevant for the temporary storage of solar thermal energy.16 It is worth noting that, in cases where several systems were compared, fairly similar heats of adsorption were measured for different AlPOs and SAPOs, despite the topological (and, in some instances, compositional) differences among the systems studied. For example, Ristic et al. reported heats of adsorption in the range of 55 kJ mol−1 for AlPO-18, AlPO-tric, and SAPO-34.16 Typically, the heat of adsorption decreases considerably with increasing water loading in SAPO materials, since the first adsorbed water molecules can interact directly with the framework protons, whereas weaker interactions with the framework and other water molecules dominate at higher loadings.13 For AlPOs, where such strong interaction sites are absent, it can be expected that the heat of adsorption changes only moderately as a function of the water loading. In addition to the investigation of powder samples, composites of AlPOs or SAPOs and support materials, primarily aluminium sheets or foams,5,17–20 have been successfully prepared and characterised. The development of stable composites is pivotal for the actual application of these materials in adsorptive heat transformations. This is exemplified by a SAPO-34-coated heat exchanger, which was successfully integrated into a lab-scale adsorption chiller.20
Computational chemistry methods at different levels of theory have been employed to predict and understand the interaction of AlPOs/SAPOs with water at a microscopic level. Henninger and co-workers used grand-canonical Monte Carlo (GCMC) simulations with an empirical force field to predict water adsorption isobars of AlPO-18.15 Despite rather promising results, this method has not been widely employed for water in AlPOs or SAPOs, whereas there are numerous studies of water adsorption in aluminosilicate zeolites.21 A larger number of computational investigations relied on electronic structure methods, especially density-functional theory (DFT). Early periodic DFT studies addressed the interaction of water with the acid sites of SAPO-34.22,23 It was shown that the adsorption of a single water molecule does not lead to framework deprotonation, but that the interaction with several water molecules can induce the formation of H3O+(H2O)n clusters. A later study of the chabazite-type aluminosilicate SSZ-13 showed that the likelihood of framework deprotonation decreases with temperature.24 A number of DFT-based studies addressed the structure and, in some instances, dynamics of water in the pores of several aluminophosphates.25–28 These studies delivered insights into the interatomic interactions governing water adsorption: For example, the rapid filling of the unit cell of AlPO-34 could be explained with the formation of a collective network of hydrogen bonds,26 and the presence of octahedrally coordinated framework aluminium atoms in water-loaded AlPO-18 was found to be the energetically favoured scenario, in agreement with experimental observations.27 Furthermore, ab initio molecular dynamics (MD) calculations corroborated that the adsorbed molecules have a considerable freedom of motion at room temperature.25–27 Recent DFT-based studies of water in SAPO-34 have addressed various aspects: a series of investigations by Fjermestad et al. elucidated the role of water during framework desilication and silicon island formation.29–31 Van Speybroeck and co-workers used ab initio MD calculations to study the proton mobility and framework flexibility in water- and methanol-loaded SAPO-34, as well as the effect of water on methanol-to-olefin conversion reactions.32,33 Particular attention was paid to the changes in lattice parameters in guest-loaded structures, with the most pronounced change being a contraction of the structure along the c-axis upon water adsorption, in agreement with previous experimental observations.34 Finally, we studied the influence of local heterogeneities (silicon islands, SiAl domains) and defects on the interaction of SAPO-34 with water.35 While a significant effect was observed at low water loadings, the impact on the total interaction strength at high loadings was only modest.
The present work aims at a more systematic computational exploration of the influence of various structural factors (pore size, pore topology, environment of framework protons) on the interaction of AlPOs and SAPOs with water, covering both energetic and structural aspects. Dispersion-corrected density-functional theory (DFT-D) calculations are employed to study the adsorption in six structurally different AlPOs and their SAPO analogues. In addition to AlPO-34/SAPO-34,36,37 AlPO-17/SAPO-17,37–40 and AlPO-18/SAPO-18,38,41,42 pairs where either the AlPO or the SAPO system has already been proposed as adsorbent for heat transformation applications,5,6,12–16,20 three other systems are evaluated: AlPO-GIS/SAPO-43,40,43,44 AlPO-AFX/SAPO-56,45 and AlPO-RHO/SAPO-RHO.9,46,47 With the exception of the last two systems, which so far have only been reported as SAPO materials, not as pure aluminophosphates, all other materials have been successfully synthesised in both AlPO and SAPO form. In all six structures considered, the pores are connected through eight-ring windows. The different structure types were chosen in order to study the influence of pore size and pore geometry: while the main pore systems of AlPO-34, AlPO-17, and AlPO-AFX are formed by elongated cages, AlPO-GIS, AlPO-18, and AlPO-RHO contain more or less isometric pores of different size. The building unit constituting the six structures are visualised in Fig. 1, with one of the main cages being highlighted for each system.
The remainder of this article is organised as follows: after a description of the model systems and the computational methodology, the results of the DFT-D optimisations of guest-free AlPO and SAPO systems are presented. These are followed by a detailed discussion of the DFT-D results for water-containing systems: AlPOs typically exhibit a steep rise of the water adsorption isotherm at a certain relative pressure (corresponding to pore filling). Below this pressure, only very little water is adsorbed due to the absence of preferred adsorption sites. Therefore, only the interaction with large amounts of water (near saturation) is studied for these systems. For SAPOs, initial adsorption at low water pressures will occur at the framework protons, prior to a pore filling at higher relative pressures. Consequently, both the adsorption of small amounts of water (one H2O per framework proton) and large amounts (near saturation) are considered in the calculations for these materials. Besides reporting the DFT-D interaction energies (averaged over 5 snapshots for high water loadings), the structural changes upon water adsorption are also assessed.
In order to study the interaction with water, different strategies were pursued for AlPOs and SAPOs: while only high water loadings (near saturation) were considered for AlPOs, both low loadings (one H2O per framework proton) and near-saturation conditions were included for SAPOs. In the DFT-D calculations for water-loaded systems, all atomic coordinates and the lattice parameters were relaxed. The same constraints as for the guest-free models were applied for systems containing one water molecule per proton, whereas no constraints were invoked at high water loadings.
Starting structures for SAPOs with one molecule of H2O per proton were prepared by placing the guest molecule in the vicinity of the proton. To investigate AlPOs and SAPOs at high water loadings, preliminary grand-canonical Monte Carlo simulations were performed to determine the approximate saturation uptake of each system. These calculations used the Sorption module of the Accelrys “Materials Studio” suite, with each run comprising five million equilibration steps and five million production steps. The GCMC simulations employed parameters from the Consistent Valence Force Field (cvff).50 These parameters were also used in a GCMC study of water adsorption in AlPO-18 by Henninger and co-workers, where reasonable agreement with experimental adsorption isobars was observed.15 It is worth noting that several of the AlPOs/SAPOs contain smaller cavities bordered by six-rings in addition to the main pore system that is connected by eight-ring windows. As these cavities (d6r and can units shown in Fig. 1) are in principle large enough to accommodate a water molecule, it will depend on the ability of the molecules to pass through the six-ring windows whether the cavities contribute to the total water uptake. In a previous combined MD and experimental study of water in AlPO-34, a water molecule was observed at the centre of the d6r unit.26 On the other hand, these units seem to be inaccessible to water in AlPO-18.27 In any case, the largest part of the adsorbed water will be located in the larger pores, which is why the present study ignores the possibility of adsorption in these small cavities. These areas were blocked by non-interacting spheres in the preliminary GCMC simulations.
The input structures for the DFT-D calculations were generated from fixed-loading Monte Carlo simulations that also used the cvff parameters. The amount of water molecules used in these calculations was always somewhat lower than the saturation uptake predicted from the GCMC simulations (typically ∼15% below the saturation loading). For each system, the fixed-loading MC simulation consisted of at least one million equilibration steps and five million production steps. A total of 200 snapshots were stored throughout the production stage, and five snapshots were randomly selected from the whole production stage. These snapshots were then used as initial configurations in the DFT-D calculations.
One has to bear in mind that the present approach involves one key approximation: for the systems containing large amounts of adsorbed water, the actual interaction energy would correspond to a weighted average over a practically infinite number of possible arrangements of water molecules in the pores. In the approach used here, this is approximated by averaging the interaction energies calculated for a small number of optimised configurations. Admittedly, this is a rather simplistic procedure; however, the findings of our previous work,35 and the fact that the standard deviations arising from the averaging over five snapshots are fairly small (never exceeding ±2.0 kJ mol−1 and being smaller than ±1.0 kJ mol−1 in half of the cases, see below), indicate that one arrives at a reasonable estimate of the interaction energy, especially if the main aim is an identification of trends, rather than a highly accurate quantitative prediction. A more thorough sampling of configuration space (including temperature effects) could be obtained by running a series of ab initio MD calculations starting from different snapshots. However, in the view of the rather large estimated error bars of 10 kJ mol−1 reported in a recent ab initio MD study of water and methanol adsorbed in SAPO-34,33 it can be anticipated that a computational setup that reduces the error to an acceptable level would be prohibitively expensive, in particular since some of the systems studied have rather large unit cells.
The PBE-TS interaction energies reported in this study were calculated as follows:
Eint = EPBE-TS(adsorbent + nH2O) − EPBE-TS(adsorbent) − nEPBE-TS(H2O) |
The first term corresponds to the PBE-TS energy obtained for the adsorbent (AlPO or SAPO) with n water molecules adsorbed, the second term corresponds to the PBE-TS energy of the guest-free adsorbent, and the last term represents an isolated H2O molecule placed in a large box (edge length 20 Å). Throughout this work, all interaction energies are reported as negative values.
In order to directly compare the calculation results to experimental heats of water adsorption, it would be necessary to explicitly calculate vibrational contributions (zero-point vibrational energy ZPVE and temperature contributions). Since the water-containing AlPO and SAPO models studied in this work contain a large number of non-equivalent atoms in the unit cell, a full vibrational calculation would be computationally very expensive. However, we performed preliminary (PBE) calculations for SAPO-34 with one adsorbed water molecule, obtaining a correction term of ∼+7 kJ mol−1 for temperatures near room temperature. We observed previously that the PBE-TS interaction energies are typically 8 to 10 kJ mol−1 higher than experimental heats of adsorption (in absolute values – interaction energies are negative, whereas heats of adsorption are positive by definition).35 This magnitude is in good accordance with the PBE-based temperature correction. With few exceptions that are mentioned explicitly, only the uncorrected interaction energies are discussed throughout this work. As there is reason to assume that the contribution of vibrations will be similar for all AlPOs/SAPOs studied, the trends identified herein should not be affected by the neglect of vibrational effects.
Topology | Space group | a/Å | b/Å | c/Å | β/deg | |
---|---|---|---|---|---|---|
AlPO-34, DFT | CHA | R | 13.792 | 14.972 | ||
AlPO-34, Exp | 13.744 | 14.941 | ||||
AlPO-17, DFT | ERI | P63/m | 13.205 | 15.385 | ||
AlPO-17, Exp | 13.146 | 15.350 | ||||
AlPO-AFX, DFT | AFX | P1c | 13.774 | 20.048 | ||
AlPO-GIS, DFT | GIS | Fddd | 13.979 | 13.778 | 10.349 | |
AlPO-18, DFT | AEI | C2/c | 13.802 | 12.787 | 18.649 | 89.82 |
AlPO-18, Exp | 13.746 | 12.753 | 18.608 | 90.00 | ||
AlPO-RHO, DFT | RHO | I23 | 15.123 |
For those three systems for which experimental data for calcined, pure AlPO systems are available, the DFT-D optimised unit cell dimensions agree very well with the experimental cell parameters. There is a consistent tendency to overestimate the lattice parameters slightly, with relative deviations of not more than 0.5%.
After the DFT-D optimisations, geometric calculations using the PoreBlazer code were performed to determine the accessible volume fraction, the diameter of the largest sphere that can be included without overlapping with the framework (dLS), and the pore limiting diameter (dPLD).59 For the systems with elongated cages, the largest sphere diameter does not fully describe the pore dimensions, as its value depends only on the shortest axis of the cage. To account for this, the length of the longest axis of the cages of AlPO-34, AlPO-17, and AlPO-AFX was estimated from the distance between the centres of the six-rings at the top and the bottom of the cages, subtracting 2 Å to account for the van der Waals radii of the surrounding framework atoms. This quantity, termed dc because the direction of cage elongation is parallel to the crystallographic c-axis in these systems, is also given in Table 2.
V acc/cm3 g−1 | d LS/Å | d c/Å | d PLD/Å | |
---|---|---|---|---|
AlPO-34 | 0.376 | 6.93 | 9.8 | 3.57 |
AlPO-17 | 0.337 | 6.53 | 13.4 | 3.16 |
AlPO-AFX | 0.378 | 7.24 | 14.9 | 3.50 |
AlPO-GIS | 0.315 | 4.77 | — | 2.92 |
AlPO-18 | 0.378 | 7.24 | — | 3.62 |
AlPO-RHO | 0.412 | 10.32 | — | 3.76 |
The accessible volume fractions in the six structures fall in a reasonably narrow range, which leads us to expect that the maximal water uptake should also be similar. In the systems with elongated cages, the largest sphere diameter is fairly similar, which is not surprising given the structural similarities (Fig. 1). The degree of cage elongation in AlPO-34 (dc ≈ 1.4·dLS) is considerably smaller than in AlPO-17 and AlPO-AFX (dc ≈ 2.05·dLS). In the three systems with isometric cages, the largest sphere diameter increases in the order AlPO-GIS < AlPO-18 < AlPO-RHO. Regarding the pore limiting diameter, which describes the diameter of the largest sphere that can pass through the eight-ring windows, four of the systems have very similar values in the range of 3.5 to 3.6 Å. An inspection of the structures reveals that the eight-ring windows in these systems are almost circular, with typical oxygen–oxygen distances (measured across the cage) of at least 6.5 Å. The exceptions are AlPO-17 and AlPO-GIS: in these systems, the eight-ring windows exhibit an elliptical distortion, which is more pronounced in AlPO-GIS (shortest O–O distance across the cage of 5.72 Å) than in AlPO-17 (shortest O–O distance of 6.0 Å). Due to their small kinetic diameter (2.64 Å), water molecules will be able to diffuse through the windows of all these structures; however, it can be anticipated that an even more drastic elliptical distortion than in AlPO-GIS might impede the diffusion of water through eight-ring windows.
The SAPO models employed for each system are presented in detail in the ESI,† where the silicon and proton locations are also visualised. To distinguish the different models, labels of the form “SAPO-N_SiX_OY” are used, where X designates the location of the silicon atom, and Y labels the oxygen atom to which the framework proton is attached (OY is bonded to SiX). The label SiX is left out if there is only one Si site. Tables S2–S7 of the ESI† report the lattice parameters and relative energies for all models considered. It is apparent that there are many cases where several possible proton locations are very close in energy. To keep the total number of systems tractable, only the lowest-energy proton location for each silicon site was considered in the DFT-D calculations including adsorbed water molecules.
Table 3 gives an overview over the number of different models considered for each structure type, and reports the energetically most favourable scenarios. Due to a lack of available crystallographic data for calcined SAPOs, and the dependence of the DFT-optimised lattice parameters on the distribution of Si and H atoms in the framework, a detailed comparison to experimental data was not attempted for these systems. However, the good agreement for AlPOs found above, together with the previous observation that the PBE-TS optimised lattice parameters of different models of SAPO-34 agree well with experimental values,35 indicate that this functional should provide a realistic prediction of SAPO structures.
Space group | N(T) per u.c. | N(Si) per u.c. | x = N(Si)/N(T) | #(Si sites) | #(H sites) | Energetically preferred model(s) | |
---|---|---|---|---|---|---|---|
a Non-conventional setting of space group Cc. b Per primitive cell. | |||||||
SAPO-34 | P32 | 36 | 3 | 0.083 | 1 | 4 | SAPO-34_O1 |
SAPO-17 | P21 | 36 | 2 | 0.056 | 2 | 7 | SAPO-17_Si1_O3, SAPO-17_Si2_O5 |
SAPO-56 | Bn | 48 | 2 | 0.042 | 2 | 8 | SAPO-56_Si1_O5, SAPO-56_Si2_O4 |
SAPO-43 | P | 32 | 2 | 0.063 | 1 | 4 | SAPO-43_O12 |
SAPO-18 | Cn | 24b | 2b | 0.083 | 3 | 12 | SAPO-18_Si1_O12, SAPO-18_Si2_O11, SAPO-18_Si3_O31 |
SAPO-RHO | I2 | 24b | 2b | 0.083 | 1 | 4 | SAPO-RHO_O11 |
A clear trend regarding the energetically preferred proton positions can be identified in those structures in which d6r units are surrounded by four-rings and eight-rings (SAPO-34, SAPO-18, SAPO-56): in these systems, the proton is either bonded to an equatorial oxygen atom of the d6r unit, pointing into an eight-ring window, or it is located at the top of the d6r unit, pointing across a six-ring (ESI,† Fig. S1, S3 and S5). In many instances, the DFT-D energies of these two cases are so close together that a more or less statistical occupation of both positions appears likely for real systems. In SAPO-17, which also contains d6r units, both favoured proton positions are closely associated with six-rings: in SAPO-17_Si1, the proton points across one of the distorted six-ring windows of the can cage, forming a fairly short H⋯O contact of 2.2 Å across the window. If the silicon atom occupies position Si2, the most favourable proton location is located inside the isolated six-ring, which is practically planar (ESI,† Fig. S2). In SAPO-43 and SAPO-RHO, the energetically preferred proton locations are associated with eight-ring windows (ESI,† Fig. S4 and S6). For these systems, a comparison of the most favourable systems to other possible arrangements shows a tendency to maximise the distance between the framework protons. Scenarios in which two protons point into the same eight-ring window are particularly unfavourable, which seems plausible from simple electrostatic considerations.
N H2O,saturation (GCMC) | N H2O (DFT) | |
---|---|---|
AlPO-34/SAPO-34 | ∼45H2O per u.c. (∼370 g kg−1) | ∼40H2O per u.c. |
AlPO-17/SAPO-17 | ∼35H2O per u.c. (∼290 g kg−1) | ∼30H2O per u.c. |
AlPO-AFX/SAPO-56 | ∼58H2O per u.c. (∼360 g kg−1) | ∼50H2O per u.c. |
AlPO-GIS/SAPO-43 | ∼37H2O per u.c. (∼340 g kg−1) | ∼32H2O per u.c. |
AlPO-18/SAPO-18 | ∼29H2O per p.u.c. (∼360 g kg−1) | ∼25H2O per p.u.c. |
AlPO-RHO/SAPO-RHO | ∼30H2O per p.u.c. (∼370 g kg−1) | ∼25H2O per p.u.c. |
While the present study does not aim at an accurate quantitative prediction of the saturation uptake, it is nevertheless insightful to compare the results to available experimental data: Jänchen at co-workers reported a water uptake of 283 g kg−1 for AlPO-17 and 388 g kg−1 for AlPO-18 at room temperature and a relative pressure p/p0 of 0.3, in good agreement with the values predicted from the GCMC simulations.13 On the other hand, they obtained a lower uptake for SAPO-34 than the calculations (279 g kg−1). At a higher relative pressure (p/p0 = 0.7) and 40 °C, an uptake of ∼320 g kg−1 was reported in a subsequent study of SAPO-34, approaching the value predicted in the present work.16 Altogether, the reasonable correspondence with experimentally determined uptakes gives confidence that the GCMC simulations provide a sufficiently reliable estimate of the saturation uptake.
E int/kJ mol−1 | E disp/Eint | ΔV/% | E int/kJ mol−1 | E disp/Eint | ΔV/% | ||
---|---|---|---|---|---|---|---|
AlPO-34 | −65.3 ± 0.5 | 0.38 | −0.6 ± 0.1 | SAPO-34_O1 | −72.6 ± 0.8 | 0.34 | −1.2 ± 0.4 |
AlPO-17 | −64.0 ± 0.2 | 0.38 | −1.1 ± 0.2 | SAPO-17_Si1_O3 | −69.7 ± 0.8 | 0.34 | −0.8 ± 0.3 |
SAPO-17_Si2_O5 | −72.0 ± 1.2 | 0.33 | −1.2 ± 0.1 | ||||
AlPO-AFX | −66.0 ± 1.1 | 0.36 | −1.1 ± 0.2 | SAPO-56_Si1_O5 | −69.9 ± 0.6 | 0.34 | −1.2 ± 0.2 |
SAPO-56_Si2_O4 | −69.2 ± 0.7 | 0.35 | −1.3 ± 0.2 | ||||
AlPO-GIS | −65.2 ± 0.6 | 0.43 | 0.0 ± 0.4 | SAPO-43_O12 | −71.1 ± 1.6 | 0.39 | −0.3 ± 0.7 |
AlPO-18 | −66.0 ± 1.2 | 0.36 | −1.0 ± 0.3 | SAPO-18_Si1_O12 | −73.0 ± 0.8 | 0.33 | −1.5 ± 0.3 |
SAPO-18_Si2_O11 | −72.5 ± 2.0 | 0.33 | −1.5 ± 0.4 | ||||
SAPO-18_Si3_O31 | −73.3 ± 1.4 | 0.32 | −1.1 ± 0.4 | ||||
AlPO-RHO | −64.8 ± 1.3 | 0.34 | −2.6 ± 0.7 | SAPO-RHO_O11 | −70.7 ± 1.1 | 0.31 | −2.8 ± 0.4 |
An assessment of the changes in volume upon water adsorption (changes in individual lattice parameters are supplied in the ESI†), together with a closer inspection of the individual DFT-D optimised snapshots reveals a qualitatively different behaviour for all systems containing d6r units (AlPO-34, -17, -AFX, -18) on the one hand, and those systems that do not contain double six-rings (AlPO-GIS, -RHO) on the other hand. For the AlPOs with d6r units, a slight reduction in volume upon water adsorption is observed, which is typically in the range of −1%. In most cases, the changes in the lattice parameters differ considerably among different snapshots, meaning that no preferred directions of contraction/expansion can be identified. In AlPO-17, however, there is a relatively clear tendency to contract along the c-axis (average change −1.8%), which is partly compensated by an expansion along the a-axis (average change +0.4%). This deformation corresponds to a shortening of the long axis of the eri cage. In the four systems with d6r units, the direct coordination of water molecules to framework aluminium atoms, leading to five-coordinated Al (AlV), is not observed very frequently (the number of AlV atoms observed in each snapshot is supplied in the ESI†): while the individual number of AlV atoms varies between 0 and 3 among different snapshots, there is on average less than one case of AlV per unit cell. Similar observations were made in a previous study of AlPO-34 and SAPO-34, where models incorporating a somewhat smaller amount of water were studied (30 molecules per unit cell instead of 40).35
In AlPO-GIS, the volume remains practically constant. However, an evaluation of the changes in individual lattice parameters reveals a significant deformation of the unit cell, despite the negligible change in volume: an expansion by ∼1.8% along a and b is almost exactly compensated by a contraction of ∼−3.5% along the c-axis. As discussed above, the eight-ring windows in guest-free AlPO-GIS are elliptically distorted, with the elongated axis being approximately parallel to the c-axis. Upon water adsorption, this elliptical distortion is considerably reduced. The pronounced change in framework structure coincides with the frequent formation of five-coordinated Al atoms: between two and four of these atoms per unit cell are found in the different snapshots, and there is also one case of six-coordinated Al (Fig. 2).
For the case of AlPO-RHO, a pronounced volume reduction of ∼−2.6% is observed, with the volume changes of the snapshots varying between −1.9% and −3.6%. A closer inspection of the lattice parameters after transformation to the pseudo-cubic setting of the cell reveals that this contraction is not isotropic, but that the unit cell shrinks anisotropically along one or two directions (Table S13, ESI†). The contraction is associated with a distortion of the lta cages. Since a considerable degree of flexibility of the RHO structure has been observed and rationalised for the aluminosilicate form,60,61 it is not surprising that rather large structural changes upon water adsorption are predicted for the AlPO system. As in AlPO-GIS, the formation of five-coordinated framework aluminium atoms is observed quite frequently in AlPO-RHO, with between one and three occurrences (per primitive cell) in the different snapshots, and one occurrence of six-coordinated Al. In the majority of cases, the water molecules which coordinate to the framework Al atoms are located inside the double eight-ring units, as shown exemplarily in Fig. 2.
In those cases where the formation of five-/six-coordinated Al is observed, the Al–Owater distances range from 1.9 to 2.1 Å. The formation of the additional bond(s) causes an elongation of the intra-framework bonds around the aluminium atom to values between 1.75 and 1.9 Å (compared to ∼1.7 to 1.75 Å in guest-free AlPOs).
E int/kJ mol−1 | E disp/Eint | d(O–H)free/Å | d(O–H)ads/Å | d(H⋯OH2O)/Å | d(HH2O⋯O)/Å | H2O location | |
---|---|---|---|---|---|---|---|
SAPO-34_O1 | −93.4 | 0.25 | 0.974 | 1.034 | 1.531 | 2.21/2.71 | Inside 8-ring |
SAPO-17_Si1_O3 | −74.2 | 0.23 | 0.984 | 1.060 | 1.467 | 2.06/2.23 | Above 6-ring |
SAPO-17_Si2_O5 | −78.1 | 0.30 | 0.975 | 1.065 | 1.446 | 1.87/1.91 | Above 6-ring |
SAPO-56_Si1_O5 | −88.7 | 0.23 | 0.975 | 1.074 | 1.419 | 1.91/1.93 | Above 6-ring |
SAPO-56_Si2_O4 | −97.1 | 0.24 | 0.974 | 1.041 | 1.503 | 2.21/2.40 | Inside 8-ring |
SAPO-43_O12 | −101.5 | 0.24 | 0.974 | 1.060 | 1.450 | 1.88/2.23 | Inside 8-ring |
SAPO-18_Si1_O12 | −97.5 | 0.22 | 0.974 | 1.045 | 1.492 | 1.99/3.01 | Inside 8-ring |
SAPO-18_Si2_O11 | −98.5 | 0.23 | 0.974 | 1.043 | 1.499 | 2.11/2.76 | Inside 8-ring |
SAPO-18_Si3_O31 | −89.2 | 0.20 | 0.975 | 1.077 | 1.415 | 1.85/1.98 | Above 6-ring |
SAPO-RHO_O11 | −89.4 | 0.24 | 0.975 | 1.033 | 1.530 | 2.35/2.83 | Above 8-ring |
Very similar features are observed for all those systems where the framework proton is attached to an equatorial oxygen atom of a d6r unit, pointing into an eight-ring window (SAPO-34_O1, SAPO-56_Si2_O4, SAPO-18_Si1_O12, SAPO-18_Si2_O11): the interaction with the water molecule inside the eight-ring window is fairly strong, with Eint usually exceeding −95 kJ mol−1 (SAPO-34 is the only exception), while the length of the hydrogen bond d(H⋯OH2O) between the framework proton and the H2O oxygen atom is quite long, in the range of 1.5 Å. There are additional (“secondary”) contacts of the H2O hydrogen atoms to some oxygen atoms surrounding the eight-ring windows. Typically, one secondary contact is significantly shorter than the others (with d(HH2O⋯O) in the range of ∼2.0 to 2.2 Å), i.e. the molecule assumes an off-centre position inside the window to maximise the interaction with the framework. A representative example (SAPO-56_Si2_O4) is shown in Fig. 3 (top). In SAPO-43_O12, the water molecule is also located inside an eight-ring window. Here, the interaction is even stronger (−101 kJ mol−1) because the distortion of the eight-ring window permits two short secondary contacts (Fig. 3 middle). Of all systems where the water molecule is associated with an eight-ring window, the weakest interaction is observed in SAPO-RHO_O11 (−89 kJ mol−1). The decomposition of the total interaction energy shows that this is primarily due to a reduced non-dispersive contribution. A closer inspection of the local environment reveals that the water molecule is not located inside the eight-ring window, but above it (Fig. 3 bottom). This leads to fairly long secondary contacts, and thus weaker electrostatic interactions between the guest molecule and the framework.
In those cases where the proton is located at the top of a d6r unit (SAPO-56_Si1_O5, SAPO-18_Si3_O31), the water molecule assumes a tripod-like configuration above the d6r unit, with three relatively short hydrogen–oxygen contacts (Fig. 4, top illustrates this for SAPO-56_Si1_O5). Despite these short contacts, the total interaction is comparatively weak, in the range of −89 kJ mol−1, with both the non-dispersive term and the dispersion contribution being smaller than for systems where the water molecule is located inside an eight-ring window. In SAPO-17, the interaction with adsorbed water is weakest among all systems, with interaction energies in the range of −75 kJ mol−1. In both SAPO-17_Si1_O3 and SAPO-17_Si2_O5, the protons point across six-ring windows, with relatively short contacts to the oxygen atoms at the opposite side of the ring. Upon water adsorption, the proton is displaced considerably from its previous equilibrium position in both cases, and the weak overall interaction can be explained with the energetic “penalty” associated with this displacement (Fig. 4 middle and bottom).
In addition to the interatomic distances between the atoms of the water molecule and the closest framework atoms, the length of the intra-framework O–H bond d(O–H) before and after H2O adsorption is also given in Table 6. In the guest-free structures, the bond lengths d(O–H)free fall in a narrow range, the only exception being SAPO-17_Si1_O3, where the bond is somewhat elongated due to the interaction between the framework proton and the oxygen atom located at the opposite side of the distorted six-ring. Upon water adsorption, the O–H bonds are elongated to a different extent in different systems, and there is a dependence on the size of the ring with which the framework proton (and, correspondingly, the adsorbed molecule) is associated: in cases where water is adsorbed inside or above an eight-ring window, d(O–H)ads ranges between 1.03 Å and 1.06 Å, whereas it is increased to 1.06 to 1.08 Å in systems where the water adsorption site is associated with a six-ring. While there is no overall correlation between d(O–H)ads and the interaction energy, reasonable correlations are observed if systems where H2O is located inside eight-rings or above six-rings are treated separately (ESI,† Fig. S7). Within each group, an increased elongation corresponds to a stronger interaction, however, the elongation corresponding to a certain interaction energy is much larger if the water molecule is adsorbed at a proton that is associated with a six-ring. These observations provide a possible means for the spectroscopic distinction of different framework protons, as the frequency shift upon water adsorption should be more pronounced in systems where the proton is associated with a six-ring.
For completeness, we also summarise the changes in the unit cell dimensions upon adsorption of small amounts of water in SAPOs: the overall changes in volume are very modest, ranging between virtually no change and a slight volume contraction of up to −0.8%. The typical changes in individual lattice parameters are even smaller (Table S14, ESI†). For the case of SAPO-34, this observation is in accordance with the findings of the previous computational (ab initio MD) study of Van Speybroeck and co-workers, who observed only minimal changes in the lattice dimensions upon addition of small amounts of water.32
For those systems where different Si distributions were considered, particularly SAPO-56 and SAPO-18, the variations among the interaction energies obtained for different models are minimal, so it can be concluded that, in contrast to the pronounced differences at low water loadings, the actual location of silicon atoms and framework protons hardly affects the interaction with water at near-saturation conditions. The same observation was made in our earlier study of SAPO-34, where the interaction energies obtained for a model with eight-ring protons (SAPO-34_O1) was practically the same as for a model with six-ring protons (SAPO-34_O3).35 For this particular SAPO, the DFT-D interaction energy after applying an approximate correction for vibrational effects of +7 kJ mol−1 compares very favourably with experimentally measured heats of adsorption at high water loadings, which are in the range of 65 kJ mol−1.6
An analysis of the changes in volume (and in individual lattice parameters) upon water adsorption shows that the SAPOs behave very similarly to their AlPO counterparts, with the deformations being on average somewhat larger: SAPO-34, -56, and -18 exhibit a slight shrinkage upon water adsorption without preferential orientation, whereas a shortening of the longest axis of the eri cage is observed in SAPO-17. In SAPO-43, an expansion along a and b is offset by a contraction along c, and for the case of SAPO-RHO, a pronounced volume shrinkage in the range of ∼−3% is observed, which is often (though not always) fairly anisotropic.
An inspection of the individual snapshots of water-loaded SAPOs shows that the majority of the framework protons have left their initial position after the DFT-D optimisation, forming H3O+ ions or, less frequently, H3O+(H2O) clusters (Fig. 5 top). The charged species in the final configurations are often, but not always, located in the vicinity of the original position of the framework proton. It has been established in previous computational work that framework deprotonation does not occur when only one or a few water molecules are adsorbed per proton, while it becomes energetically favourable when larger water clusters (e.g. tetramers) interact with one proton.22–24 A recent ab initio MD study of water-loaded SAPO-34 emphasised the high mobility of the protons in this system at a temperature of 350 K.32 While the present work, which is based on static DFT optimisations, cannot quantify the effect of temperature on the mobility, the observation that framework deprotonation occurs in the large majority of cases is in accordance with these previous findings. The removal of framework protons from their initial positions is observed frequently in all SAPOs studied, regardless of the topology (or silicon content).
In the SAPOs, the formation of five-coordinated framework atoms occurs with a similar frequency as in their AlPO counterparts, with relatively few occurrences in SAPO-34, -17, -56, and -18, and a larger number of AlV atoms in SAPO-43 and SAPO-RHO. There are two examples where no AlV atoms are observed in any of the snapshots (SAPO-34_O1, SAPO-18_Si3_O31), however, this does not necessarily reflect a reduced tendency of these systems to form five-coordinated Al due to the limited number of configurations sampled (in this context, it is worth noting that AlV atoms were found in several snapshots of water-loaded SAPO-34_O1 in our previous study).35 As discussed previously, ab initio MD simulations would be very helpful to study the formation of five-coordinated Al in a more quantitative fashion, especially with regard to the influence of temperature. Another interesting observation from the present investigation is the preferential occurrence of five-coordinated aluminium in the direct vicinity of the silicon atoms, which is most pronounced in SAPO-43_O12 (Fig. 5 bottom).
(1) Due to the absence of strongly preferred adsorption sites in AlPOs, the interaction with small amounts of water was investigated only for SAPOs. The interaction is markedly stronger when the framework proton points into an eight-ring than in situations where the proton is located above a six-ring. If it was possible to fully control the proton locations in a real material, the interaction energy could be “tuned” in a certain range, with eight-ring protons leading to higher energy densities, and six-ring protons allowing for lower desorption temperatures. However, since the different proton sites are often extremely close in energy (e.g. in SAPO-18), such a control is unlikely to be possible in real materials. Moreover, the considerable mobility of the protons would lead to a rearrangement during water adsorption/desorption cycles. Nevertheless, further insights could be expected from a comparative study of zeotypes that contain eight-rings, but no six-rings on the one hand (e.g. SAPO-43), and systems with six-rings and larger rings, but no eight-rings, on the other hand (e.g. SAPO-5 [AFI topology, twelve-rings + six-rings] or SAPO-11 [AEL topology, ten-rings + six-rings]).40 It would be interesting to see whether the conclusions of the present study also hold for a wider range of structure types.
(2) At high water loadings, there is no evidence for a significant effect of the topology on the material's affinity towards water. When accounting for the contribution of thermal motion in an approximate fashion, heats of adsorption in the range of 55 to 60 kJ mol−1 can be expected for AlPOs, whereas somewhat higher values between 60 and 65 kJ mol−1 will be attained in SAPOs due to stronger electrostatic interactions. In SAPOs, the interaction with water can be enhanced by increasing the silicon content, which is directly correlated with the amount of framework protons as long as no silicon islands are formed. While no effect of the pore size on the interaction strength is observed, it has to be noted that a smaller pore diameter will lead to a pore filling at lower relative pressures. Therefore, the pore size may still determine the suitability of a certain adsorbent for given operating conditions.
(3) Unlike the adsorption energetics, the structural response of the adsorbents to water shows a considerable variation among the different topologies: while the systems that contain d6r units undergo only a moderate contraction, a much stronger deformation of the framework is observed in the zeotypes with GIS and RHO topologies (regardless of the composition). Moreover, the coordination of water molecules to framework aluminium atoms is observed more frequently in these materials. In this context, it is worth noting that a correlation between the strong binding of water to framework Al atoms and irreversible loss of crystallinity upon hydration of SAPOs has been found in a comparative MAS-NMR study of SAPO-34 (stable) and SAPO-37 (FAU topology, decomposes upon hydration at room temperature).62 Due to the severe structural changes observed for water-loaded GIS- and RHO-type AlPOs and SAPOs, it can be hypothesised that irreversible transformations might occur during water adsorption/desorption cycles. Such changes are likely to have a detrimental effect on the long-term performance of an adsorbent, rendering it unsuitable for heat transformation applications.
Altogether, the DFT-D computations permit us to conclude that the adsorption energetics in AlPOs and SAPOs are hardly influenced by the topology. There is thus little scope to “tune” the interaction strength – which is directly related to the amount of thermal energy that can be stored – by a judicious choice of a certain framework type. More promising strategies that can be envisaged are the variation of the silicon content in SAPOs and the incorporation of other elements in the framework (metal-containing aluminophosphates [MeAPOs]).40 On the other hand, the topology determines the structural response of the material to water, which is why rather large variations in long-term stability have to be expected for AlPOs and SAPOs having different framework types. In this context, the present results, and possible future studies using ab initio MD simulations, may provide useful information that will facilitate the identification of suitably stable adsorbents.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c6cp02289h |
This journal is © the Owner Societies 2016 |