research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoSTRUCTURAL SCIENCE
CRYSTAL ENGINEERING
MATERIALS
ISSN: 2052-5206

Incommensurate structure of AlPO4-5 and its stability

crossmark logo

aGeochemical Research Center, Graduate School of Science, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo, Tokyo, 113-0033, Japan, bResearch Institute for Chemical Process Technology, National Institute of Advanced Industrial and Science and Technology, 4-2-1 Nigatake, Miyagino-ku, Sendai, Miyagi, 983-8551, Japan, and cResearch Institute for Chemical Process Technology, National Institute of Advanced Industrial and Science and Technology, 1-1-1 Higashi, Tsukuba, Ibaraki, 305-8565, Japan
*Correspondence e-mail: [email protected]

Edited by M. Dusek, Academy of Sciences of the Czech Republic, Czechia (Received 18 September 2025; accepted 28 April 2026; online 26 May 2026)

The incommensurate structure of AlPO4-5 and its stability have been investigated in detail. Under ambient conditions, single-crystal X-ray diffraction measurements reveal weak satellite reflections, which become significantly enhanced upon exposure to vacuum. Temperature-dependent X-ray diffraction and Raman scattering measurements of the O–H stretching region show that the appearance of satellite reflections is due to the desorption of water molecules under vacuum, and the satellite reflections disappear above 335 (5) K, whereas their intensity increases again upon cooling to room temperature. To elucidate the origin of the modulation, single-crystal X-ray diffraction experiments were performed on a dehydrated sample. Diffraction data including satellite reflections up to second order were successfully collected, with a modulation vector of q = 0.36710 (5)c*. Analysis of systematic absences up to second-order satellites indicate that the (3+1)D superspace group is P6(00γ)h. On the basis of these data, the first-order positional modulation amplitudes were refined. The results indicate that the modulation originates from correlated rotational displacements of AlO4 and PO4 tetrahedra. This behaviour suggests that phonon-like motions, as proposed in previous studies, persist over a long modulation period and give rise to the observed satellite reflections. These findings highlight a fundamental connection between the intrinsically spatial description of atomic positional modulations in the incommensurate structure and the underlying temporal dynamics.

1. Introduction

Aluminophosphate molecular sieves have been widely studied since their discovery in the 1980s (Wilson et al., 1983View full citation) owing to their unique adsorption properties and catalytic applications by the substitution of Al and P atoms with the heteroatoms of Si, Co, and so on (Kang et al., 2008View full citation; Li et al., 2025View full citation) . Numerous studies have focused on the physical and chemical properties of one of the aluminophosphate family, AlPO4-5, such as water adsorption (Li et al., 2023View full citation), membrane filtration (Khamis et al., 2023View full citation), and laser activation (Zhu et al., 2024View full citation). However, despite decades of study, its detailed crystal structure remains incompletely understood.

Earlier single-crystal X-ray diffraction studies of as-synthesized AlPO4-5 with organic structural directing agents (OSDAs) (Bennett et al., 1983View full citation, Qiu et al., 1989View full citation), or time-of-flight powder neutron diffraction of the calcined material (Richardson et al., 1987View full citation) have revealed its fundamental framework structure within the AFI-type zeolite family with space group P6cc. The structure is characterized by channels formed by 12-membered rings of AlO4 and PO4 tetrahedra stacked along the c axis, forming a one-dimensional channel of approximately 7.3 Å in diameter (Fig. 1[link]). Additionally, owing to the charge difference between Al and P, as well as the overall symmetry, the structure exhibits polarity. This polarity has been experimentally confirmed through pyroelectricity measurements (Klap et al., 1998View full citation).

[Figure 1]
Figure 1
(a) Polyhedral representation of the crystal structure of AlPO4-5. Green and purple tetrahedra depict AlO4 and PO4 tetrahedra, respectively. Oxygen atoms are not shown. (b) Ball-and-stick representation of the 12-membered ring in the crystal structure of AlPO4-5.

A subsequent powder X-ray diffraction study including fluoride anions and tropine as an OSDA showed a reversible phase transition from hexagonal to orthorhombic space group Ccc2 (Ohnishi et al., 1993View full citation), and further powder neutron and X-ray diffraction studies for the calcined sample also showed Pcc2 (Mora et al., 1996View full citation, Ikeda et al., 1999View full citation), a subgroup of Ccc2, as the most plausible space group. Subsequent computational studies of AlPO4-5 without OSDAs provided important suggestions that the symmetry reduction to P6 from P6cc resolved the unnatural linear Al–O–P angle in the P6cc structure (Henson et al., 1996View full citation), and the structure with the space group Pcc2 is very similar to the P6 structure (Ruiz-Salvador et al., 1996View full citation). Klap and co-authors conducted detailed structural analyses using single-crystal X-ray diffraction and proposed that the crystal structure of AlPO4-5 consists of three types of microdomains, each exhibiting P6 symmetry, and the averaged structure has the space group P6cc (Klap et al., 1998View full citation).

The crystal structure of AlPO4-5 has also attracted attention in terms of the structure and behaviour of guest water molecules. For example, Floquet et al. proposed a model in which water molecules inside the 12-membered rings form a double helix (Floquet et al., 2004View full citation), but this model has not been experimentally verified and may not be entirely plausible, considering the channel size. Computational studies have suggested that the guest structure varies depending on the number of water molecules present (Demontis et al., 2012View full citation). Additionally, extensive research has been conducted on the relationship between the polarity of guest molecules and that of the host framework (Klap et al., 1998View full citation), as well as the interaction between water dynamics and the channels (Gotoh et al., 2000View full citation).

Several intriguing recent findings come from Lotti et al. and Fischer et al. (2026View full citation). Lotti et al. (2016View full citation) observed satellite peaks at q ∼ 0.37c* in X-ray diffraction patterns using a synchrotron source, indicating the presence of an incommensurate modulated structure. A recent study by Fischer et al. (2026View full citation) has investigated the dynamic disorder in AlPO4-5 using single-crystal diffuse scattering in combination with ab initio molecular dynamics simulations. Their work revealed the presence of correlated atomic displacements within the framework and provided important insights into the origin of the observed diffuse scattering. In particular, the study suggested that the structural fluctuations can be interpreted in terms of collective motions resembling phonon-like behaviour. Despite these advances, the detailed structural analysis of the modulation in AlPO4-5 has remained unexplored. Furthermore, the stability of the modulated structure with respect to temperature, pressure, and guest water molecules has not been systematically investigated.

In this study, we utilized relatively large single crystals [0.1 × 0.07 × 0.07 (mm)] to perform X-ray diffraction analysis, allowing us to investigate the previously uncharacterized modulated structure of AlPO4-5. Furthermore, Raman spectroscopy under high-temperature and vacuum conditions was employed to explore the relationship between the modulated structure and guest water molecules, providing new insights into the structural dynamics of this material.

2. Experimental

2.1. Sample preparation

To synthesize single crystals of AlPO4-5, Alumisol-10 (10 wt%, Kawaken Fine Chemicals), orthophospho­ric acid (85 wt% Sigma-Aldrich), and tri­ethyl­amine [(C2H5)3N, +99%, Tokyo Chemical Industry] were used as Al and P sources, and an OSDA, respectively. Deionized water was used to dilute the ingredients. The synthetic solution was obtained by mixing these ingredients in the ratios of 1.0 Al2O3: 1.2 P2O5: 4.0 (C2H5)3N: 550 H2O. We referred to previous articles for the actual procedure to mix the ingredients (Demuth et al., 1995View full citation, Kodaira et al., 1999View full citation). The initial solution (pH = 7.15) was adjusted to pH = 2.80 by dropwise addition of concentrated sulfuric acid (98−99%, Fujifilm Wako). This solution was sealed in a Teflon-lined autoclave and was placed in a preheated oven at 463 K for 68 h. After the hydro­thermal reaction, AlPO4-5 crystals were collected, washed with deionized water, filtered, and dried at 333 K. Prior to the following measurements, the crystals were calcined at 1073 K for 30 h under a dried air flow of 0.5 L min−1 to remove the OSDA.

2.2. X-ray diffraction and Raman spectroscopy at high T

To ensure that both measurements of X-ray diffraction and Raman spectra were conducted under identical conditions on the same sample, an ad hoc sample chamber originally designed for low-temperature high-pressure measurements was used for ambient-pressure experiments (Fig. 2[link]). The single-crystal sample was placed on one of the diamonds in a diamond anvil cell (DAC) made of CuBe alloy [Fig. 2[link](b)], and the DAC was mounted onto a copper frame. A liquid nitro­gen cooling line and a ceramic heater were embedded in the copper frame, and the temperature was controlled using an external temperature controller (CHINO, KP-1000), while monitoring the temperature using a platinum resistance thermo­meter embedded within the frame. The flow rate of exhausted nitro­gen gas, evaporated from the liquid nitro­gen passing through the copper frame, was controlled using a mass flow controller. The entire copper frame was housed inside an aluminium sample chamber, allowing measurements to be performed under atmospheric pressure or vacuum conditions. This chamber was equipped with two windows, one made of a polyester film (Mylar) for optical observation and Raman spectra measurements, and the other was a polyimide film (Kapton) for X-ray diffraction.

[Figure 2]
Figure 2
(a) A photograph of instruments used in this study. The sample chamber is set on the goniometer of X-ray diffractometer with a microscope and a Raman probe. (b) Schematic drawings for the sample chamber and the enlarged side view around the sample.

The sample temperature was monitored using a K-type thermocouple fixed near the diamond on which the sample was placed. Within the experimental temperature range (298–340 K), the temperature difference between the thermocouple and the platinum resistance thermometer in the copper frame reached a steady state within approximately 10 min after changing the set temperature. In the steady state, the temperature difference was within 0.5 K, and given the excellent thermal conductivity of the diamond, the sample temperature was estimated to be accurately monitored within approximately ±1 K.

For X-ray diffraction measurements, a rotating anode X-ray source (Mo Kα, Rigaku, MicroMax-007) was used. The incident X-ray beam was focused onto the sample area using a confocal X-ray mirror (Rigaku, VariMax-Mo) and further collimated with a single-pinhole collimator with a hole radius of ϕ = 0.3 mm. Diffraction images were collected using an imaging plate detector (Rigaku, R-AXIS IV++). The sample chamber was mounted onto the goniometer of the X-ray diffractometer, and oscillation photographs were obtained by exposing the sample for 10 min while oscillating it within a φ range of −14° to −8° to observe the satellite peaks.

Raman spectroscopic measurements were also conducted on the same sample inside this chamber used for X-ray diffraction measurements. A 50 mW DPSS laser (Airix) emitting green laser light at a wavelength of 532 nm was introduced into a confocal optical system placed beside the X-ray diffractometer via an optical fibre. The backscattered light from the sample passed through a narrow-band notch filter within the optical system and was then transmitted through an optical fibre to a spectrometer (Bruker Optics 250, grating of 1200 mm−1), where the spectra were recorded using a high-sensitivity CCD (Andor, iDus). Raman spectra were measured in the OH stretching region at 298 K under atmospheric pressure and at 300 K, 310 K, and 320 K under vacuum (≃ 10−1 Pa).

2.3. X-ray diffraction for detailed structure analysis at room temperature

The single-crystal sample was placed in a drying oven at 373 K for 1 h in order to completely dehydrate the sample. To prevent moisture absorption, the sample was immediately coated with an adhesive upon removal from the oven. After mounting the sample on a glass fibre, it was placed on the goniometer of an X-ray diffractometer (Rigaku, Synergy Custom). The X-ray source used was a rotating anode (Mo Kα, Rigaku, MicroMax-007), and the beam was focused onto the sample area using a confocal X-ray mirror (Rigaku, VariMax-MoHF) and collimated with a 0.3 mm collimator. Diffraction images were collected using a hybrid photon-counting X-ray detector (Rigaku, HyPix-6000HE). The obtained diffraction images were processed using CrysAlisPro (v. 1.171.42.70a; Rigaku Oxford Diffraction). The 3D averaged structure analysis based on the main reflections was performed with SHELXL (Sheldrick, 2015View full citation) implemented in Olex2-1.5 (Dolomanov et al., 2009View full citation), and the electron density distribution in the unit cell of AlPO4-5 was estimated using the maximum entropy method (MEM) implemented in the Dysnomia program (Momma et al., 2013View full citation). The crystal structures with the electron density were visualized using VESTA3 (Momma & Izumi, 2011View full citation). The analysis of the modulated structure, including satellite reflections, was performed using JANA2020 (Petříček et al., 2023View full citation). The initial model for the modulated structure was derived using the charge flipping algorithm implemented in SUPERFLIP (Palatinus & Chapuis, 2007View full citation).

3. Results and discussion

3.1. Stability of the incommensurate structure

During the initial screening to verify the crystallinity of the sample, we confirmed the presence of satellite reflections at approximately q ∼ 0.37c*, as reported by Lotti et al. (2016View full citation), even at ambient conditions. However, the intensity of these satellite reflections varied significantly between samples, with some exhibiting almost no observable satellites. Notably, samples that showed no satellite reflections had been exposed to air for extended periods since they were taken out of the sample bottle. This observation led us to hypothesize that moisture absorption from the atmosphere could influence the intensity of the satellite reflections.

To investigate this, we placed a single-crystal sample with weak satellite reflections under vacuum. As a result, faint satellite reflections became observable after evacuation [Fig. 3[link](a)]. Raman spectroscopy further revealed that placing the sample under vacuum alone caused a significant decrease in the intensity of peaks associated with OH stretching vibrations [Fig. 3[link](b)]. Moreover, as the temperature was gradually increased in 10 K increments, the OH stretching peak became nearly undetectable at 310 K, coinciding with an increase in the intensity of the satellite reflections [Fig. 3[link](c)]. With further heating, the intensity of the satellite reflections peaked at 320 K and then began to decrease, becoming almost undetectable at 340 K. Upon cooling, the satellite reflections reappeared at 330 K, and their intensity increased progressively as the temperature decreased. During the heating and cooling process, the positions of the satellite peaks did not change significantly.

[Figure 3]
Figure 3
(a) X-ray oscillation photographs of AlPO4-5 during a heating–cooling cycle. The upper-left panel shows the diffraction pattern at 298 K under atmospheric pressure together with its magnified view. The subsequent enlarged images to the right (upper row) were collected during heating under vacuum, while those in the lower row correspond to cooling under the same vacuum conditions. The lower-left panel shows the diffraction pattern after the heating–cooling cycle at 300 K. The powder rings observed in the vacuum measurements originate from scattering by the Kapton window and are absent in the ambient-pressure image. (b) Raman spectra measured at 298 K under atmospheric pressure and at 298 K, 310 K, and 320 K under vacuum. (c) Intensities of a satellite peak, 0211 reflection, recorded during the heating and cooling process.

These results indicate that the satellite reflections are most intense when the water molecules are fully evacuated from the structure. Additionally, we identified an incommensurate–disordered phase transition occurring at approximately 335 (5) K. This transition temperature is slightly lower than the transition at 380 K from hexagonal to orthorhombic symmetry observed in powder X-ray diffraction studies of AlPO4-5 containing OSDA molecules (Ohnishi et al., 1993View full citation). Furthermore, it closely matches the temperature at which the pyroelectric coefficient and spontaneous polarization begin to increase (Klap et al., 1998View full citation). They also reported that the temperature at which these properties change can vary depending on the presence and type of OSDA molecules. Although modulated structures have not been reported in AlPO4-5 containing OSDA molecules, the atomic displacements responsible for the long-period modulation, which will be discussed in the next section, strongly suggest a connection to the symmetry reduction observed in previous studies.

3.2. Analysis of 3D averaged structure

We revisited the crystal structure of AlPO4-5 using the observed main reflections, which represent the 3D averaged structure projected from the (3+1)D superstructure. We observed many strong hhl reflections (l = odd), which clearly violate the reflection conditions (hhl: l = even) of the previously proposed space group P6cc. We also observed weak but significant intensities for the 001 and 003 reflections (Fig. 4[link]), indicating that no systematic extinction conditions are present. This observation suggests that the true symmetry of AlPO4-5 is lower than that of the space group P6cc. Assuming that the true space group is a subgroup of P6cc and considering the absence of extinction conditions, P6 is a plausible candidate.

[Figure 4]
Figure 4
Pseudo-precession images of AlPO4-5 showing h0lm reflections. The right inset shows an enlarged reciprocal area around the 4000 reflection, which exhibits second-order satellites and reduced contrast was applied to clarify each satellite peak.

An initial crystal structure based on space group P6 was obtained by direct methods, and the resulting model is isotopological to the previously proposed structure with space group P6cc (Klap et al., 2000View full citation). Structure refinement based on this initial model, with anisotropic atomic displacement parameters (ADPs), resulted in extremely flattened, elongated, or even degenerate displacement ellipsoids for several O atoms. To account for this behaviour, we adopted a split-site model for these O atoms, assigning half or one-third occupancy depending on the number of split positions, with isotropic ADPs, as shown in Fig. 5[link]. The final R factors remained considerably larger (R1 = 0.0727) compared with the small Rint value of 0.0151 (see Table 1[link]).

Table 1
Experimental details

  3D Averaged structure (3+1)D Modulated structure
Crystal data    
Chemical formula AlO4P AlO4P
Mr 122 122
Crystal system, space group Hexagonal, P6 Hexagonal, P6(00γ)h
Temperature (K) 295 295
Modulation wavevectors q = 0.36710 (5)c*
a, c (Å) 13.76570 (17), 8.37040 (12) 13.7657 (17), 8.37040 (12)
V3) 1373.64 (3) 1373.64 (3)
Z 12 12
Radiation type Mo Kα Mo Kα
μ (mm−1) 0.67 0.67
Crystal size (mm) 0.1 × 0.07 × 0.07 0.1 × 0.07 × 0.07
 
Data collection
Diffractometer Rigaku, Synergy Custom, HyPix Rigaku, Synergy Custom, HyPix
Absorption correction Empirical Empirical
No. of measured, independent and observed [I > 2σ(I)] reflections 24151, 1617, 1614 103791, 8733, 8597
No. of parameters 111 180
Rint 0.015 0.057
(sin θ/λ)max−1) 0.624 0.625
h,k,l range −17 → 17, −17 → 17, −10 → 10 −17 → 17, −17 → 17, −10 → 10, −2 → 2
 
Refinement
Refinement method F2 F2
No. of reflections 1617 8733
No. of parameters 111 252
R[F2 > 2σ(F2)], wR(F2), S for overall reflections 0.073, 0.138, 1.15 0.082, 0.220, 7.54
R[F2 > 2σ(F2)], wR(F2), No. of main reflections 0.057, 0.184, 1733
R[F2 > 2σ(F2)], wR(F2), No. of first-order satellite reflections 0.066, 0.181, 3531
R[F2 > 2σ(F2)], wR(F2), No. of second-order satellite reflections 0.163, 0.312, 3469
Δρmax, Δρmin (e Å−3) 0.76, −0.79 0.66, −0.67
Symmetry operations: (1) x1, x2, x3, x4; (2) −x2, x1 − x2, x3, x4 + ⅓; (3) −x1 + x2, −x1, x3, x4 + Mathematical equation; (4) −x1, −x2, x3, x4 + ½; (5) x2, −x1 + x2, x3, x4 + Mathematical equation; (6) x1x2, x1, x3, x4 + Mathematical equation.
[Figure 5]
Figure 5
(a) 3D averaged and (b) (3+1)D crystal structure models with electron density distribution derived from the MEM for asymmetric unit of AlPO4-5 structure. Shaded surfaces depict averaged electron densities with isosurface level of 2.0 e Å−3 derived from the MEM using the main reflections. In (b), the (3+1)D structure models drawn while changing t (phase of modulation) in steps of 0.1 from 0 to 0.9 are superposed.

One possible reason for this discrepancy is that the complex distribution of atomic displacements in the (3+1)D superspace cannot be adequately modelled within split-site structural models. Supporting this interpretation, analysis using MEM, which allows greater flexibility in the electron density distribution than conventional structural models, yielded an excellent agreement factor of RF = 0.0243.

3.3. Analysis of incommensurate structure

The obtained single-crystal diffraction images for dried AlPO4-5 exhibit clear satellite peaks with the modulation vector of q = 0.36710 (5)c* up to second-order indices of m, which were first observed in this study (Fig. 4[link]). As shown in Fig. 4[link], the satellite reflections along the c* direction are absent, indicating that, for the 00lm reflections, those with m = 1 and 2 are systematically extinct. Based on space group P6, the possible (3+1)D superspace groups are therefore limited to two candidates: P6(00γ)h (No. 168.2; reflection condition: 00lm: m = 6n) and P6(00γ)t (No. 168.3; 00lm: m = 3n), according to International Tables for Crystallography, Vol. C, Table 9.8.3.5. Modulated structure models were then derived using the charge-flipping method based on these two superspace groups. As a result, a consistent modulated structure model was obtained only for P6(00γ)h. The modulation should be derived from the positional displacements of the host framework, not from occupancy modulation nor displacements of guest molecules, considering the results of the heating experiments under vacuum. Thus, the positional displacement parameters for the first-order harmonics were obtained using Superflip (Palatinus & Chapuis, 2007View full citation) and refined up to second-order harmonics using JANA2020 (Petříček et al., 2023View full citation).

Although the apparent Laue class is 6/mmm, whereas space group P6 belongs to Laue class 6/m, the possibility of merohedral twinning should be considered. However, structure refinements taking into account the possible twin operations (−1, .m., .2.), as listed in Table 1.3.4.1 of International Tables for Crystallography, Vol. C (Prince, 2004View full citation), did not result in any significant improvement compared with refinements without considering twinning. Furthermore, the final R factor is sufficiently low (Table 1[link]), the anisotropic ADPs yielded physically reasonable ellipsoids, and the Al–O and P–O bond lengths fall within acceptable ranges without the application of any special restraints [Fig. 5[link](b)]. These observations suggest that the probability of crystal twinning is low. On the other hand, it should be noted that determination of the absolute structure is difficult in this study, because Mo Kα radiation was used, and the sample does not contain any heavy elements.

3.4. Origin of the modulation

The incommensurate structure observed in AlPO4-5 is characterized by rotations of the AlO4 and PO4 tetrahedra (Figs. 5[link] and 6[link]). As shown in Fig. 5[link], the positions of Al and P atoms exhibit minimal displacement with respect to changes in the modulation phase, t, whereas the O atoms either trace circular paths around their average positions or oscillate along linear trajectories. This distribution of atomic positions is in good agreement with the averaged electron density distribution derived from the MEM analysis of the main reflections [Fig. 5[link](b)].

[Figure 6]
Figure 6
(a) The ball-and-stick structure models of the 12-membered ring in the AlPO4-5 drawn while changing t (phase of modulation) in steps of 0.1 from 0 to 0.9. (b) Simply superposed image of the ten images in (a). (c) Superposed image of the ten images in (a), translated so that the basal O atoms of each AlO4 or PO4 tetrahedron overlap. The numbers 1–6 in (a) and (c) are assigned for convenience to the apical O atoms of the AlO4 tetrahedra. See also video 1 in supporting information.

Let us examine the structural changes with varying modulation phase, t. Focusing on the 12-membered ring, it is evident that the apical O atoms of the AlO4 tetrahedra—i.e. the O atoms in the Al–O–P bonds parallel to the c axis—rotate counterclockwise with increasing t [Fig. 6[link](a)]. Owing to the symmetry constraints of the superspace group, P6(00γ)h, the phases of rotation for apical O atoms located on opposite sides of the 12-membered ring are identical; that is, atoms 1–4, 2–5, and 3–6 in Fig. 6[link](a) rotate in phase, while adjacent apical O atoms exhibit phase shifts of 2π/3 with respect to each other. Such displacements are also suggested by rigid-unit-mode (RUM) calculations for SiO4 tetrahedra in SSZ-24, a silicate zeolite with an isotypic framework of AlPO4-5 (Liu et al., 2002View full citation). According to the RUM analysis, the B1 mode possesses a threefold rotational softening at k = (0, 0, 0.392). The softening of these phonon modes has been proposed to account for the satellite reflections observed at q = 0.37c* (Lotti et al., 2016View full citation) and q = 0.36710 (5)c* in this study.

When the structural models at different t values are simply superposed, the basal O atoms located in the ab plane also appear to move [Fig. 6[link](b)]. However, when the structures are superposed such that these basal O atoms coincide, they overlap almost perfectly [Fig. 6[link](c)]. This indicates that the basal O atoms within the 12-membered ring translate without changing their relative positions in the ab plane. This behaviour can be more intuitively appreciated in supplementary video 1.

On the other hand, when the modulated structure is viewed along a direction perpendicular to the c axis, the basal O atoms are seen to undergo undulating displacements along the direction parallel to the c axis as t varies, in concert with neighbouring O atoms (supporting video 2). During this modulation, the bond lengths and bond angles within the AlO4 and PO4 tetrahedra remain nearly unchanged—even though no geometrical restraints were imposed in the refinement—indicating that the tetrahedra behave essentially as rigid units. These characteristics of the structural variation with varying t apparently resemble the vibrational modes of phonons, as proposed for the origin of diffuse scattering in AlPO4-5 (Fischer et al., 2026View full citation).

It should be noted that the aforementioned rotation or vibration represents spatial phase shifts with respect to the modulation phase t and does not necessarily imply real dynamic temporal motion. However, in the case of AlPO4-5, dynamic motion may indeed play a role. The results from quasi-elastic neutron scattering and molecular dynamics simulations (Cortie et al., 2017View full citation) indicate that, at temperatures above at least 100 K, the O atoms bridging Al and P atoms exhibit motion along circular trajectories. In addition, diffuse scattering measurements at ambient conditions have been interpreted in terms of phonon-like correlated motions within the framework (Fischer et al., 2026View full citation), suggesting that such dynamics can give rise to long-range correlations that manifest as satellite reflections when stabilized as a long-period modulation. In the present study, high-temperature X-ray diffraction experiments reveal that the modulated structure disappears above 335 (5) K, whereas quasi-elastic neutron scattering (QENS) data indicate that O atoms are already mobile at 300 K. If this interpretation of the QENS data is correct, the tetrahedra must be dynamically vibrating while maintaining long-range phase correlations. The persistence of a long-period incommensurate structure in the presence of dynamically moving units, despite thermal vibrational perturbations, is an astonishing phenomenon. Furthermore, as shown by the appearance of satellite reflections upon dehydration in this study, it is reasonable to infer that such dynamic long-period structures are extremely sensitive to subtle changes in external conditions, such as the presence or absence of guest molecules. These findings highlight the intricate interplay between framework flexibility, modulation phenomena, and host–guest interactions in AlPO4-5.

Supporting information


Computing details top

(I) top
Crystal data top
AlO4PF(000) = 720
Mr = 122Dx = 1.769 Mg m3
Hexagonal, P6(00γ)hMo Kα radiation, λ = 0.71073 Å
q = 0.367100c*Cell parameters from 17280 reflections
a = 13.76570 (17) Åθ = 2.9–26.3°
c = 8.37040 (12) ŵ = 0.67 mm1
V = 1373.64 (3) Å3T = 295 K
Z = 120.10 × 0.08 × 0.07 mm
† Symmetry operations: (1) x1, x2, x3, x4; (2) −x2, x1x2, x3, x4+1/3; (3) −x1+x2, −x1, x3, x4+2/3; (4) −x1, −x2, x3, x4+1/2; (5) x2, −x1+x2, x3, x4+5/6; (6) x1x2, x1, x3, x4+1/6.

Data collection top
ROD, SynergyCustom system, HyPix
diffractometer
8733 independent reflections
Radiation source: Rotating-anode X-ray tube8597 reflections with I > 3σ(I)
Mirror monochromatorRint = 0.057
Detector resolution: 10.0000 pixels mm-1θmax = 26.4°, θmin = 1.7°
ω scansh = 1717
Absorption correction: multi-scan
CrysAlisPro 1.171.44.103a (Rigaku Oxford Diffraction, 2025) Empirical absorption correction using spherical harmonics, implemented in SCALE3 ABSPACK scaling algorithm.
k = 1717
Tmin = 0.81, Tmax = 1l = 1010
103791 measured reflections
Refinement top
Refinement on F22 constraints
R[F2 > 2σ(F2)] = 0.082Weighting scheme based on measured s.u.'s w = 1/[σ2(Fo2) + (0.02P)2]
where P = (Fo2 + 2Fc2)/3
wR(F2) = 0.221(Δ/σ)max = 0.020
S = 7.49Δρmax = 0.66 e Å3
8733 reflectionsΔρmin = 0.67 e Å3
252 parametersAbsolute structure: 0 of Friedel pairs used in the refinement
0 restraints
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
P10.45425 (7)0.33235 (7)0.92319 (11)0.0221 (4)
P20.33328 (7)0.45551 (7)0.42294 (11)0.0217 (4)
Al10.33220 (7)0.45424 (7)0.79830 (11)0.0123 (4)
Al20.45537 (7)0.33311 (7)0.29851 (11)0.0119 (4)
O10.42184 (18)0.2104 (2)0.8603 (4)0.0353 (11)
O70.2102 (2)0.42125 (17)0.3610 (4)0.0336 (11)
O40.3630 (2)0.3636 (2)0.8602 (4)0.0375 (11)
O60.4230 (3)0.5794 (3)0.3619 (4)0.0390 (12)
O30.5779 (3)0.4238 (3)0.8613 (4)0.0405 (12)
O80.3688 (2)0.3682 (2)0.3618 (4)0.0397 (12)
O50.3333 (2)0.4553 (2)0.6121 (4)0.0471 (14)
O20.4551 (2)0.3335 (2)1.1129 (4)0.0477 (14)
Atomic displacement parameters (Å2) top
U11U22U33U12U13U23
P10.0222 (5)0.0184 (4)0.0268 (6)0.0109 (4)0.0019 (3)0.0015 (3)
P20.0194 (5)0.0225 (5)0.0265 (6)0.0130 (4)0.0010 (3)0.0011 (3)
Al10.0100 (4)0.0141 (4)0.0136 (5)0.0067 (4)0.0021 (3)0.0021 (3)
Al20.0139 (5)0.0116 (4)0.0131 (5)0.0087 (4)0.0007 (3)0.0008 (3)
O10.0296 (13)0.0223 (13)0.0556 (15)0.0143 (12)0.0002 (13)0.0017 (11)
O70.0241 (13)0.0259 (12)0.0505 (15)0.0122 (12)0.0013 (12)0.0001 (13)
O40.0362 (15)0.0363 (15)0.0463 (13)0.0229 (12)0.0022 (13)0.0001 (13)
O60.0387 (16)0.0364 (16)0.0485 (14)0.0238 (12)0.0061 (12)0.0047 (12)
O30.0362 (15)0.0377 (16)0.0517 (15)0.0216 (12)0.0063 (12)0.0068 (12)
O80.0355 (15)0.0351 (15)0.0519 (15)0.0202 (12)0.0028 (13)0.0055 (13)
O50.0500 (18)0.0547 (18)0.0385 (17)0.0275 (15)0.0041 (13)0.0012 (13)
O20.0511 (18)0.0500 (18)0.0417 (18)0.0250 (15)0.0008 (14)0.0040 (13)
Geometric parameters (Å, º) top
AverageMinimumMaximum
P1—O11.615 (9)1.602 (9)1.635 (9)
P1—O41.622 (10)1.615 (10)1.632 (10)
P1—O31.628 (7)1.621 (8)1.633 (8)
P1—O21.622 (10)1.602 (10)1.640 (10)
P2—O71.610 (8)1.596 (9)1.622 (9)
P2—O61.632 (7)1.622 (8)1.643 (8)
P2—O81.615 (11)1.597 (11)1.632 (11)
P2—O51.621 (9)1.605 (10)1.640 (10)
Al1—O1i1.598 (11)1.584 (12)1.613 (12)
Al1—O41.610 (10)1.600 (10)1.623 (10)
Al1—O3ii1.613 (8)1.602 (8)1.621 (8)
Al1—O51.593 (9)1.580 (10)1.608 (10)
Al2—O7iii1.600 (11)1.590 (11)1.610 (11)
Al2—O6ii1.613 (7)1.603 (8)1.620 (8)
Al2—O81.609 (11)1.599 (11)1.632 (11)
Al2—O2iv1.591 (9)1.565 (10)1.607 (10)
O1—P1—O4109.2 (4)107.8 (4)111.2 (4)
O1—P1—O3111.0 (5)101.5 (4)122.1 (5)
O1—P1—O2110.5 (5)107.2 (6)112.0 (5)
O4—P1—O3108.9 (4)104.0 (4)112.6 (5)
O4—P1—O2108.6 (6)105.1 (5)111.7 (6)
O3—P1—O2108.7 (4)106.0 (5)110.7 (4)
O7—P2—O6109.6 (4)108.0 (4)111.1 (5)
O7—P2—O8110.4 (4)104.4 (4)117.7 (4)
O7—P2—O5109.2 (5)107.1 (5)112.8 (5)
O6—P2—O8109.6 (4)96.7 (4)122.9 (5)
O6—P2—O5107.9 (4)104.8 (4)109.8 (4)
O8—P2—O5109.6 (6)107.3 (6)113.1 (5)
O1i—Al1—O4109.0 (4)107.2 (4)111.7 (4)
O1i—Al1—O3ii110.6 (5)98.8 (4)123.5 (5)
O1i—Al1—O5110.2 (5)106.0 (6)113.4 (5)
O4—Al1—O3ii108.7 (4)103.2 (4)112.9 (5)
O4—Al1—O5108.5 (6)105.6 (5)111.5 (6)
O3ii—Al1—O5109.3 (4)107.8 (4)110.8 (5)
O7iii—Al2—O6ii108.7 (4)104.9 (4)112.5 (5)
O7iii—Al2—O8110.4 (4)100.3 (4)122.3 (4)
O7iii—Al2—O2iv109.4 (5)107.7 (5)111.1 (6)
O6ii—Al2—O8108.3 (4)99.5 (4)115.8 (5)
O6ii—Al2—O2iv109.1 (4)106.4 (4)112.8 (5)
O8—Al2—O2iv110.1 (6)106.8 (7)113.5 (6)
Symmetry codes: (i) x1x2, x1, x3, x4+1/6; (ii) x1+1, x2+1, x3, x4+1/2; (iii) x2, x1+x2, x3, x4+5/6; (iv) x1, x2, x31, x4.
 

Acknowledgements

We thank Vaclav Petricek for his kind advice on the usage of JANA2020. Crystal structure models in Figs. 1[link] and 5[link] were drawn using the VESTA program (Momma & Izumi, 2011).

Funding information

The following funding is acknowledged: Japan Society for the Promotion of Science (grant No. 25K01724 to Kazuki Komatsu; grant No. 21K18154 to Kazuki Komatsu; grant No. 20310066 to Tetsuya Kodaira; grant No. 23350107 to Tetsuya Kodaira; grant No. 17K19000 to Tetsuya Kodaira; grant No. 25K01582 to Takuji Ikeda; grant No. 25K01579 to Takuji Ikeda).

References

Return to citationBennett, J. M., Cohen, J. P., Flanigen, E. M., Pluth, J. J. & Smith, J. V. (1983). Intrazeolite Chemistry, Vol. 218, edited by G. D. Stucky and F. G. Dwyer, pp. 109–118. American Chemical Society.  Google Scholar
Return to citationCortie, D. L., McBride, B. R., Narayanan, N., de Souza, N. R., Avdeev, M., Mole, R. A., McIntyre, G. J., Kearley, G. J., Withers, R., Yu, D. H. & Liu, Y. (2017). J. Phys. Chem. C 121, 18762–18770.  CrossRef CAS Google Scholar
Return to citationDemontis, P., Gulín-González, J. & Suffritti, G. B. (2012). J. Phys. Chem. C 116, 11100–11109.  CrossRef CAS Google Scholar
Return to citationDemuth, D., Stucky, G. D., Unger, K. K. & Schüth, F. (1995). Microporous Mater. 3, 473–487.  CrossRef CAS Google Scholar
Return to citationDolomanov, O. V., Bourhis, L. J., Gildea, R. J., Howard, J. A. K. & Puschmann, H. (2009). J. Appl. Cryst. 42, 339–341.  Web of Science CrossRef CAS IUCr Journals Google Scholar
Return to citationFischer, M., Chaney, D. A., Uwumwonse, O. A., Neumann, E., Spieß, I., Wiggers, C. & Schmidt, E. M. (2026). Microporous Mesoporous Mater. 401, 113912.  CrossRef Google Scholar
Return to citationFloquet, N., Coulomb, J. P., Dufau, N. & Andre, G. (2004). J. Phys. Chem. B 108, 13107–13115.  CrossRef CAS Google Scholar
Return to citationGotoh, K., Ishimaru, S. & Ikeda, R. (2000). Phys. Chem. Chem. Phys. 2, 1865–1869.  CrossRef CAS Google Scholar
Return to citationHenson, N. J., Cheetham, A. K. & Gale, J. D. (1996). Chem. Mater. 8, 664–670.  CrossRef CAS Google Scholar
Return to citationIkeda, T., Miyazawa, K., Izumi, F., Huang, Q. & Santoro, A. (1999). J. Phys. Chem. Solids 60, 1531–1535.  CrossRef CAS Google Scholar
Return to citationKang, L., Zhang, T., Liu, Z. & Han, K. (2008). J. Phys. Chem. C 112, 5526–5532.   CrossRef CAS Google Scholar
Return to citationKhamis, F., Kumar, M., Koutsoukos, P. G., Banat, F., Arafat, H. A. & Karanikolos, G. N. (2023). J. Water Process. Eng. 54, 104038.  CrossRef Google Scholar
Return to citationKlap, G. J., van Klooster, S. M., Wübbenhorst, M., Jansen, J. C., van Bekkum, H. & van Turnhout, J. (1998). J. Phys. Chem. B 102, 9518–9524.  CrossRef CAS Google Scholar
Return to citationKlap, G. J., van Koningsveld, H., Graafsma, H. & Schreurs, A. M. M. (2000). Microporous Mesoporous Mater. 38, 403–412.  Web of Science CrossRef CAS Google Scholar
Return to citationKodaira, T., Miyazawa, K., Ikeda, T. & Kiyozumi, Y. (1999). Microporous Mesoporous Mater. 29, 329–337.  CrossRef CAS Google Scholar
Return to citationLi, X., Li, X., Liu, X., Song, X., Wang, B., Xu, L. & Tao, S. (2025). Microporous Mesoporous Mater. 382, 113396.   Google Scholar
Return to citationLi, X., Wang, J., Chen, S., Lv, H., Wang, Q., Dong, E. & Yang, X. (2023). J. Mol. Struct. 1292, 136245.  CrossRef Google Scholar
Return to citationLiu, Z., Fujita, N., Terasaki, O., Ohsuna, T., Hiraga, K., Camblor, M. A., Díaz-Cabañas, M.-J. & Cheetham, A. K. (2002). Chem. Eur. J. 8, 4549–4556.  CrossRef PubMed CAS Google Scholar
Return to citationLotti, P., Gatta, G. D., Comboni, D., Merlini, M., Pastero, L. & Hanfland, M. (2016). Microporous Mesoporous Mater. 228, 158–167.  CrossRef CAS Google Scholar
Return to citationMomma, K., Ikeda, T., Belik, A. A. & Izumi, F. (2013). Powder Diffr. 28, 184–193.  Web of Science CrossRef CAS Google Scholar
Return to citationMomma, K. & Izumi, F. (2011). J. Appl. Cryst. 44, 1272–1276.  Web of Science CrossRef CAS IUCr Journals Google Scholar
Return to citationMora, A. J., Fitch, A. N., Cole, M., Goyal, R., Jones, R. H., Jobic, H. & Carr, S. W. (1996). J. Mater. Chem. 6, 1831–1835.  CrossRef CAS Google Scholar
Return to citationOhnishi, N., Qiu, S., Terasaki, O., Kajitani, T. & Hiraga, K. (1993). Microporous Mater. 2, 73–74.  CrossRef CAS Google Scholar
Return to citationPalatinus, L. & Chapuis, G. (2007). J. Appl. Cryst. 40, 786–790.  Web of Science CrossRef CAS IUCr Journals Google Scholar
Return to citationPetříček, V., Palatinus, L., Plášil, J. & Dušek, M. (2023). Z. Kristallogr. Cryst. Mater. 238, 271–282.  Google Scholar
Return to citationPrince, E. (2004). International Tables for Crystallography, Vol. C, Mathematical, physical and chemical tables. Springer.  Google Scholar
Return to citationQiu, S., Pang, W., Kessler, H. & Guth, J.-L. (1989). Zeolites 9, 440–444.  CrossRef CAS Google Scholar
Return to citationRichardson, J. W. Jnr, Pluth, J. J. & Smith, J. V. (1987). Acta Cryst. C43, 1469–1472.  CrossRef CAS IUCr Journals Google Scholar
Return to citationRuiz-Salvador, A. R., Sastre, G., Lewis, D. W. & Catlow, C. R. A. (1996). J. Mater. Chem. 6, 1837–1842.  CAS Google Scholar
Return to citationSheldrick, G. M. (2015). Acta Cryst. C71, 3–8.  Web of Science CrossRef IUCr Journals Google Scholar
Return to citationWilson, S. T., Lok, B. M., Messina, C. A., Cannan, T. R. & Flanigen, E. M. (1983). Intrazeolite Chemistry, Vol. 218, edited by G. D. Stucky and F. G. Dwyer, pp. 79–106. American Chemical Society.  Google Scholar
Return to citationZhu, J., Chen, X., Sang, X. & Yang, G. (2024). Appl. Catal. Gen. 669, 119486.  CrossRef Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoSTRUCTURAL SCIENCE
CRYSTAL ENGINEERING
MATERIALS
ISSN: 2052-5206
Follow Acta Cryst. B
Sign up for e-alerts
Follow Acta Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds