
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39237487
51788
10.1038/s41467-024-51788-7
Article
Quantitative localisation of titanium in the framework of titanium silicalite-1 using anomalous X-ray powder diffraction
http://orcid.org/0000-0003-3185-3535
Rzepka Przemyslaw 123
http://orcid.org/0000-0003-0521-3702
Signorile Matteo 4
http://orcid.org/0000-0003-4174-8958
Huthwelker Thomas 5
http://orcid.org/0000-0003-0499-4885
Checchia Stefano 6
http://orcid.org/0000-0001-9980-5551
Rosso Francesca 4
http://orcid.org/0000-0003-2371-4156
Bordiga Silvia silvia.bordiga@unito.it

4
http://orcid.org/0000-0002-4166-2284
van Bokhoven Jeroen A. jeroen.vanbokhoven@chem.ethz.ch

23
1 grid.425073.7 0000 0004 0633 9822 J. Heyrovsky Institute of Physical Chemistry Dolejškova 2155/3, 182 23 Prague 8, Czech Republic
2 https://ror.org/05a28rw58 grid.5801.c 0000 0001 2156 2780 Institute for Chemical and Bioengineering, ETH Zurich, 8093 Zurich, Switzerland
3 https://ror.org/03eh3y714 grid.5991.4 0000 0001 1090 7501 Paul Scherrer Institute, Center for Energy and Environmental Sciences, PSI, 5232 Villigen, Switzerland
4 grid.7605.4 0000 0001 2336 6580 Department of Chemistry, NIS and INSTM Reference Centre, Università di Torino, Via G. Quarello 15, I-10135 and Via P. Giuria 7, I-10125 Torino, Italy
5 https://ror.org/04xrb7c55 Swiss Light Source, PSI, 5232 Villigen, Switzerland
6 grid.5398.7 0000 0004 0641 6373 Beamline ID15A, European Synchrotron Radiation Facility 71 Avenue des Martyrs, 38000 Grenoble, France
5 9 2024
5 9 2024
2024
15 775719 5 2023
14 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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/.
One of the biggest obstacles to developing better zeolite-based catalysts is the lack of methods for quantitatively locating light heteroatoms on the T-sites in zeolites. Titanium silicalite-1 (TS-1) is a Ti-bearing zeolite-type catalyst commonly used in partial oxidation reactions with H2O2, such as aromatic hydroxylation and olefin epoxidation. The reaction mechanism is controlled by the configuration of titanium sites replacing silicon in the zeolite framework, but these sites remain unknown, hindering a fundamental understanding of the reaction. This study quantitatively determines heteroatoms within the zeolite-type framework using anomalous X-ray powder diffraction (AXRD) and the changes in the titanium X-ray scattering factor near the Ti K-edge (4.96 keV). Two TS-1 samples, each with approximately 2 Ti atoms per unit cell, were examined. Half of the titanium atoms are primarily split between sites T3 and T9, with the remainder dispersed among various T-sites within both MFI-type frameworks. One structure showed significant non-framework titanium in the micropores of a more distorted lattice. In both samples, isolated titanium atoms were more prevalent than dinuclear species, which could only potentially arise at site T9, but with a significant energy penalty and were not detected.

Here authors use anomalous X-ray powder diffraction to locate titanium atoms in two TS-1 catalysts. Titanium appears as isolated species at T3 and T9 sites within the frameworks of both catalysts and as non-framework species in one structure.

Subject terms

Heterogeneous catalysis
Porous materials
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Titanium silicalite-1 (TS-1) is a catalytically active zeolite-type material1, widely used in several partial oxidation reactions2–4. Its unique abilities to epoxidize olefins with hydrogen peroxide (H2O2) and to release only water as a by-product are of clear interest from an industrial perspective, in particular concerning the epoxidation of propylene via the so-called HPPO (hydrogen peroxide to propylene oxide) process5. Owing to the outstanding properties in catalytically active processes, TS-1 is also one of the best-characterized materials. However, the mechanism that governs those processes still remains unknown mostly because of the high dilution of titanium in the TS-1 siliceous framework (ca. 3 wt% TiO2), hence hampering a precise assessment of the exact position of the metal within the possible crystallographic sites by conventional structural and spectroscopic techniques.

TS-1 structure belongs to MFI-type topology6 and, like any zeolitic material, consists of tetrahedrally coordinated silicon atoms interconnected by covalently bonded oxygen bridges arranged into a porous construction. The complicated orthorhombic structure (Pnma space group) is constructed by two intersecting perpendicular systems of 10-membered ring channels (straight and sinusoidal) and 12 different crystallographic positions of silicon (T-sites). The outstanding catalytic performance of TS-1 in partial oxidation reactions presumably arises from the active sites that emerged by the isomorphous substitution of silicon for titanium. The minor concentration of tetrahedral titanium (up to ca. 3 wt% TiO2) corresponds to around up to 2.5 Ti atoms per unit cell that can replace only a fraction of 8 symmetry-equivalent silicon positions at each of 12 T-sites (96 per unit cell in total). The fractional occupancy of T-sites and the intrinsic complexity of the TS-1 structure blurs the image of titanium distribution, although the material has been studied for years and by many different methods (NMR7–9, other spectroscopic methods10–14, and theory15–19). It remains unknown whether the active species responsible for the formation of Ti-peroxo-species are comprised of framework titanium monomers, di-titanium species, or higher nuclearity clusters and to what extent they reside within the crystalline framework4,7. Since the geometry of the sites is defined by the vicinity of the Ti positions in the framework, pinpointing titanium to specific T-sites is essential to understand their character. Despite the long-standing discussion, diffraction studies using synchrotron X-ray20,21, and neutron radiation22,23, also exploiting Ti isotope exchanged24, have not managed to unambiguously elucidate the siting of titanium on the T-sites, partially due to its low concentration in the framework. Hence, to date, no definitive picture of the nature of the active site(s) in TS-1 has been established.

The primary objective of this study was the unequivocal identification of titanium siting within the TS-1 structure and determining whether titanium emerges as isolated species or dinuclear sites. Though X-ray atomic scattering factors f=f0+f1+i⋅f2 are primarily functions reflecting the number of electrons in elements (f0) that should enable the discrimination of Ti (Z = 22) from Si (Z = 14) by using the conventional diffraction methods, the low concentration of Ti requires a technique that is more element sensitive, such as anomalous X-ray diffraction (AXRD) at the Ti K-edge. AXRD is a method that exploits the resonant scattering of atoms near their absorption edges25,26. Many, if not most, elements have absorption edges within energy ranges suitable for anomalous diffraction experiments, e.g., Br (13.5 keV), Rb (15.2 keV), Se (12.66 keV), Zn (9.66 keV), and Fe (7.11 keV)27,28. Recently, aluminum was successfully pinpointed to the crystallographic T-sites in zeolite ferrierite29, and cations on the extra-framework positions have been located30–32. Multiwavelength anomalous diffraction is broadly used to facilitate the structure determination of biological macromolecules33.

Results

Ti anomalous X-ray scattering factors

To pinpoint titanium to T-sites in the TS-1 structure, the X-ray energy was set close to and far away from the Ti K-edge, near- and off-resonance. In the former case, the complex contributions of the scattering factor f=f0+f1+if2 become significant and the scattering power of the element changes dramatically while those from the other elements, silicon and oxygen, remain unchanged. Any changes to the pattern, therefore, relate to the location of titanium in the crystalline structure. Since the anomalous scattering on a coordinated titanium can differ from the one of the isolated atom as given in the LBL database34, f1 and f2 values were derived from Ti K-edge XANES spectrum of TS-1 structure by mean of Kramers–Kronig relation (Fig. 1A)35,36.Fig. 1 Structure factors for TS-1 structure model.

A Top: Ti K-edge XANES recorded from TS-1A structure13. Bottom: The derived f1 and f2 contributions to titanium X-ray scattering factor near the Ti K-edge (4.96 keV). The diffraction data were collected at the near- and on-resonance energies represented by vertical lines. B |ΔFhkl/Fhkl | = |(Fhkl,4.600keV - Fhkl,4.960keV)/Fhkl,4.600keV|: Relative differences between the selected structure factors (ΔFhkl > 1), Fhkl,E, calculated at off-resonance (E = 4.600 keV) and near-resonance (E = 4.960 keV) energies from a TS-1 structure model holding 0.5 Ti atoms per unit cell concentrated at T3 site.

The different intensities of the Bragg reflections measured at off- and near-resonance energies vary and enable the most unequivocal assignment of the titanium to the T-atom in the framework structure. In contrast to short-range spectroscopic techniques, such as X-ray absorption spectroscopy, which probe only the local environments of Ti atoms, AXRD searches through long-range order and can quantify titanium at different T-sites and pinpoint their positions unambiguously in the crystal lattice.

Figure 1B depicts the variations in structure factors (Fhkl) for the TS-1 model, holding half of the Ti atom in the unit cell concentrated at the T3 site when energy is switched from off- (4.600 keV) to near-resonance (4.960 keV). The structural information on titanium distribution related to the specific lattice planes (hkl) can be quantified by the intensity change of the corresponding reflection (Ihkl) when the energy approaches the Ti K-edge. Fhkl values are altered for the Ti-bearing model, whereas those for Ti-free S-1 structure remain unchanged when changing the photon energy. ΔFhkl values vary for different reflections and enable the titanium distribution to be determined.

Location of titanium in TS-1 structures

Two industrial TS-1 samples (provided by Evonik Industries AG), originating from different synthetic routes and marked as TS-1A and TS-1B, were experimentally examined (basic characterization data available in Supplementary Fig. 1 and discussed therein). Both structures have similar chemical composition Si96−xTixO192, where x ~ 2 (ca. 2.5 wt% TiO2), and demonstrated similar catalytic activity in propylene epoxidation14. TS-1B contains some fraction of segregated TiO2 (anatase polymorph) and hexacoordinated (i.e., distorted octahedral) Ti embedded in the MFI framework at Si-defect positions13, that seem to slightly affect selectivity in HPPO reaction by increasing the concentration of reaction by-products (alcoholysis and hydrolysis products) at long reaction times14.

The AXRD measurements were performed at the PHOENIX beamline at the Swiss Light Source synchrotron facility (SLS) operating at the tender energy ranges (https://www.psi.ch/en/sls/phoenix) using the previously developed set-up29. Since the 4.960 keV X-ray beam is totally absorbed after penetrating a few tens of micrometers of a zeolite sample, the experiment was carried out in reflection geometry inside the vacuum chamber. The preferred orientation of the zeolite crystallites in the vertical direction was prevented by the special sample preparation presented in the “Methods” section and depicted by SEM images (Supplementary Fig. 7). The silicon drift diode detector37 enabled energy-specific detection of single photons. The discrimination of inelastic scattering suppressed the background elevated by the fluorescence emission released after an edge and enhanced the signal-to-noise ratio. Ti K-edge XANES scans were measured for each sample to determine the positions of the absorption edge. Further technical details regarding the experiments are summarized in the “Methods” section.

The analysis followed the protocol that was previously established29. The formulated data analysis strategy involved the single structure and profile parameters model refined with the Topas software38 concurrently to fit all different energy diffraction datasets. The assignment of titanium to T-sites was preceded by the careful refinement of the orthorhombic TS-1 structure (Pnma space group) parameters against conventional diffraction data collected in transmission mode (0.5 mm capillary) at wavelength λ = 0.70860 Å (E = 17.5 keV) on the Materials Science beamline at the SLS39. All diffraction data observed and those calculated from the refined models are shown in Supplementary Figs. 3, 4. There was no evidence of monoclinic distortions, as could be expected at low (<1 Ti per unit cell) loading40. A model that converged to the conventional diffraction data was simultaneously refined against the datasets collected across the Ti K-edge (Fig. 1A), at 4.600 keV (off-resonance), 4.960 keV, and 4.972 keV (both near-resonance) (Supplementary Figs. 3, 4). The dataset collected at 4.968 keV was excluded from the further evaluation since the Kramers–Kronig transformation of XANES data fails to decouple f1 and f2 values due to strong changes in the absorption coefficient for the on-resonance condition (Fig. 1A). At first, the profile parameters for AXRD data (scale factor, peak profile, specimen displacement, and axial model) were fitted. The initial scale factors were estimated using a previously established approach29, which involved considering normalized μE—attenuation coefficients for the sample composition at the specified energies and I0—the incoming photon flux (calculated from the energy-dependent diode response function). Then the scales were allowed to be refined individually.

Once the profile functions had converged, the structure parameters were refined. The limit for observed diffraction at 4.960 keV corresponded to a d-spacing of 1.27 Å and enabled refining the atomic positions in the structure models against data across the absorption edge. Finally, the content of titanium, initially set as 2.5 Ti atoms equally divided between 12 T-sites, was refined while other parameters remained fixed at the values obtained before. The uncertainties of Ti occupancies on individual T-sites were determined by Bootstrap methods41,42. The details on the refinements (Supplementary Tables 1 and 2) and corresponding cif files are provided in the Supporting Information.

The refinement of TS-1A revealed 0.6(1) and 0.5(1) titanium per unit cell that occupy sites T3 and T9, respectively (Fig. 2A). This finding corroborates the theoretical calculations on this system that suggested those sites as the most stable titanium siting15. The tetrahedral titanium signal was also highlighted at other T-atoms of orthorhombic MFI framework6, but the individual contributions at each of the sites remained insignificant. The refinement of total titanium content converged to 2.2 Ti atoms per unit cell in TS-1A. The tetrahedral Ti loading in the structure was reported to be 1.814. Thus, all titanium is identified as a framework species.Fig. 2 Ti distribution across T-sites of TS-1 framework structure.

A Refined TS-1A framework structure model. T3 and T9 sites are represented as magenta and green spheres, remaining Si and O atoms as gray sticks. Within a unit cell, 0.6 and 0.5 Ti atoms are distributed across T3 and T9 sites (both navy blue). B Computed (PBE/DZVP) relative energy ΔE of all possible unique Ti pair models (including SDA) featuring the T3 and/or T9 site(s) as a function of the minimal Ti-Ti distance. The red circle points to the only model featuring the Ti-O-Ti moiety. The yellow box highlights the optimal range of Ti-Ti distances (11.8 ± 1.2 Å), giving the most populated structures, identified as the average of the Ti-Ti distances typical of the six most stable (ΔE < 15 kJ/mol) structures.

Whereas the TS-1B structure displays less distinct Bragg reflections than TS-1A (Supplementary Figs. 3, 4), it also shows a titanium population of 1.8 primarily hosted by sites T3 (0.5(1) Ti atom) and T9 (0.4(2) Ti atom), with the remaining titanium signal spread over various T-sites. Since the reported total Ti occurrence in the sample (2.89 wt% TiO2)14 would correspond to 2.1 Ti atoms per unit cell, the differential 0.3 titanium is non-framework. This stands in agreement with the previous findings on TS-1B that suggested over 10% contributions from higher coordination and non-framework titanium in the sample based on combined ICP, vibrational spectroscopy, and XAS studies13,14.

The robustness of the refinement is further corroborated by the relations between the lattice parameters and Ti loadings in the samples as is shown in Fig. 3. The refined lattice parameters of orthorhombic TS-1B structure a = 20.127(5), b = 19.947(5), and c = 13.417(3) stay distinctly contracted in comparison to those revealed by TS-1A, respectively a = 20.140(2) Å, b = 19.966(2) Å, and c = 13.428(1) Å (Fig. 3). One can see the linear correlation between Ti/(Ti+Si) ratios in TS-1 and the lattice parameters when compared to those from Ti-free silicalite-1 (a = 20.089(3), b = 19.895(3), and c = 13.384(2)). Linear regression coefficients of unit cell expansion (Supplementary Table 3) corroborate with those reported by ref. 21 and ref. 40.Fig. 3 TS-1 lattice expansion as a function of the Ti content.

Linear fits of A lattice parameters and B volume of the unit cell vs. the titanium content in TS-1 structures (Pnma space group). Ti-free S-1 was refined in monoclinic P21/n11 space group with α = 90.055(3)°.

0.6(1) Ti atoms occupancy at a T3 site in the unit cell of TS-1A structure occurs as an isolated species that is equivalent to the T15 site in the lower monoclinic symmetry possible for the MFI framework, found as a probably occupied one15. On the other hand, the two nearest symmetry-equivalent positions ascribed to the T9 site, possibly occupied by titanium are located adjacent to each other resulting in four possible Ti pairs per unit cell (Fig. 2A). The oxygen bridge interconnecting those two titanium atoms points towards the intersection between straight and sinusoidal 10 ring channels and enables them to act as potential dinuclear titanium site, whose existence and reactivity in TS-1 is debated7,13. Differential Pair Distribution Function (dPDF) comparing the reference pattern of the defect-free Ti-free silicalite-1 (S-1) with G(r) of TS-1 structures (see SI for further details on sample preparation) provided no evidence supporting such preferential pairing (Supplementary Fig. 5), however, such measurement does not rule out that titanium dimers on T9 site emerge during synthesis and remain undetected due to extremely low occurrence of such ordering, additionally suppressed by the associated presence of silicon defect sites. Another possible structure with a quite short Ti-Ti distance may occur with a simultaneous occupation of a T3 and a T9 sites on the same five-membered ring and separated by a single siloxane bridge (i.e., generating a Ti-O-Si-O-Ti moiety). Though not properly a dimeric Ti-O-Ti moiety, due to the metal centers' proximity, it may potentially act similarly under reaction conditions.

Computational models of titanium pairs

To pinpoint the exact reciprocal position of T3 and/or T9 sites(s) within the unit cell, a set of periodic DFT calculations (PBE/DZVP, see Methods section for full computational details) were performed on the 22 possible unique pair models generated from the two sites by assuming a loading of two Ti atoms per unit cell (main computational results in Supplementary Table 4). Since the Ti distribution in the framework takes place during the hydrothermal crystallization of the zeolite framework, the structure directing agent (SDA), i.e., tetrapropylammonium hydroxide (TPA+OH−), has been explicitly included in the models. Previous attempts to explicitly consider the role of SDA in the simulation are rare and did not provide relevant support to the interpretation of experimental data23. The probability of having a given structure was evaluated on the basis of the relative energy of the different models. This is shown as a function of the minimal distance between the two Ti atoms in Fig. 2B. Known the uncertainty of DFT methods, also related to the selected computational setup15,43, that can be in the order of 10–15 kJ/mol for the formation energy of siliceous zeolites44, it clearly emerges that few structural models account for the larger fraction of sites (six models falls below a threshold of 15 kJ/mol). Furthermore, these structures are characterized by a very specific Ti-Ti distance, on average being 11.8 ± 1.2 Å. This site distribution significantly differs from that computed in the absence of the SDA, which shows much smaller differences in relative energies (16 models below the 15 kJ/mol threshold, see Supplementary Table 5), preventing a clear discrimination of siting.

Discussion

A dimeric titanium model, capable of generating a bridging μ-η2-peroxo species upon interaction with H2O2, was proposed on the basis of analogy of its NMR signatures with those of reference metalorganic compounds, and modeled by DFT considering a Ti-O-Ti moiety across the T7 and T11 positions (both located on the walls of the straight channel of MFI along the b-axis)7. The Ti-Ti distance for such a structure, evaluated by reoptimizing the model with the method adopted in this work (without including SDA), is 3.357 Å, comparable to the Ti-Ti distance for a Ti-O-Ti moiety where both Ti atoms sit on adjacent T9 sites (3.259 Å, in the absence of SDA as well). On the other hand, assuming the statistical distribution of 0.5(1) Ti atom across equivalent positions of the T9 site in the unit cell of TS-1 structure, there is only P = 1/28 probability to find two Ti atoms in the same unit cell (P = 1/4) on two neighboring positions (P = 1/7). The probability further lowers considering the poor relative stability for this structure (56 kJ/mol over the most stable structure modeled by DFT), making the formation of dimeric Ti pairs at T9 sites unlikely on a thermodynamic basis. Thus, while dinuclear sites may theoretically emerge from bare geometrical considerations, their formation is accompanied by a severe energy penalty. The absence of distinctive features for titanium dimers, further corroborated by AXRD and dPDF and in combination with isolated siting at T3, suggests that monomeric titanium prevails in our materials.

In the mixed T3–T9 model, the Ti pair displays a Ti-Ti distance of 4.492 Å in absence of SDA, which appears too long to enable the formation of the μ2-η2:η2-peroxo species on geometrical basis (see also Supplementary Fig. 8). It could be possible that these titanium sites can still act together as a dimeric site, but hypothetically assuming a μ2- η1:η1-peroxo bridging coordination obtained through a different reaction mechanism. Nonetheless, such a Ti pair is also improbable on a thermodynamic basis due to its poor stability (36 kJ/mol) compared to the most probable models. Instead, structures with an “optimal” Ti-Ti distance (11.8 ± 1.2 Å) are much more stable. One can clearly see from Fig. 2B, as a general trend, that the relative stability of the considered models increases as Ti-Ti spacing increases. Within the six most stable structures, titanium will react with H2O2 as isolated Ti sites generating monomeric η2-peroxo species, as broadly accepted in the literature4. The graphical representation for all commented models is provided in Supplementary Fig. 6.

In conclusion, TS-1A exclusively holds tetrahedral titanium, with more than half of it preferentially distributed across sites T3 and T9. The site T3 can host only isolated titanium species. Although dinuclear titanium sites could potentially emerge at adjacent T9 sites, such a model exhibits very low stability compared to more separated Ti sites and could only occur in the case of highly preferred localization experienced during sample synthesis. The framework titanium distribution in TS-1B resembles one in TS-1A, however ca. 15% titanium presence remains in non-framework locations. These two TS-1 catalysts of similar chemical composition Si96−xTixO192 (x ~ 2) exhibit significant differences in the distribution of tetra- and octahedral titanium, despite both structures displaying elevated titanium concentration on similar T-sites within the MFI-type framework. The quantitative determination of heteroatoms in the TS-1 catalysts shows that titanium ordering in their structures follows an individual fashion of specific samples controlled by the catalysts’ origins and can differ due to various synthesis conditions. Nonetheless, provided that Ti preferentially sits at namely isolated tetrahedral sites, TS-1 can be consistently described as a heterogeneous single-site catalyst, as broadly accepted in the literature43. Anomalous X-ray diffraction (AXRD) at the Ti K-edge constitutes the tool that allows to locate titanium of concentration even below 2.5 wt% TiO2 in the TS-1 catalysts’ structures of different origin.

Methods

The AXRD experiments at PHOENIX beamline

The experiments were performed on the PHOENIX (PHotons for the Exploration of Nature by Imaging and XAFS) undulator beamline at the Swiss Light Source (SLS) in Villigen, Switzerland that uses an elliptical undulator as the source. An in-vacuum diffractometer with Bragg-Brentano geometry was installed at the PHOENIX I branch that covers the energy range from 0.8–8 keV using a double crystal monochromator. A Si 111 crystal generated monochromatic light. High harmonics was rejected by reflections on three mirrors with a reflection angle of 0.4 degree. The diffractometer chamber, separated by a 500 nm thick silicon nitride window from the beamline, provided a residual pressure lower than 10−4 mbar by using a turbomolecular pump. The slit, made of simple razorblades, placed at a fixed distance, shaped the incoming beam into a horizontal line of 2 mm × 200 μm. Two motorized goniometers rotated the sample holder and the detector to collect the diffracted photons in Θ-2Θ geometry. The sample holder was continuously rotating to compensate for the possible horizontal preferred orientation of the crystallites in the sample. The powder samples were dispersed in isopropanol and allowed to be set on a flat zero-diffraction Si wafer allowing all orientations of crystallites within a few micrometers of the sample layer. The specimen layer was also ensured to be settled flat to avoid micro-absorption effects on a rough surface. An energy-dispersive silicon drift diode (SDD, Viamp diode 50 mm2, with MOXYEC AP3 window, manufacturer Ketek, Germany) set ∼120 mm from the sample was used as the detector. Further details on the detector and the set-up exploited for AXRD experiments can be found in ref. 29.

To ensure the high statistical quality of all registered Bragg reflections, the weaker reflections were integrated for longer times to collect more photons compared to more intense reflections. All data points shown in the diffractograms are in arbitrary units proportional to the photon flux and normalized to the dwell time. The data are also corrected for the detector dead time.

Computational details

DFT calculations were performed with the CP2K package (2022.1 version)45,46 by employing the PBE GGA functional47, including Grimme’s D3 dispersion correction. Atomic cores were represented using Goedecker-Teter-Hutter (GTH) pseudopotentials48. The GPW method was adopted49, being the Kohn-Sham orbitals expanded in the MOLOPT-DZVP basis set50, while an auxiliary plane-wave basis (cutoff of 550 Ry) was used to represent the density. The energy in the SCF cycles was converged to 10−7 hartree, using the orbital transformation method51. Geometry optimization was conducted on atomic coordinates and cell parameters simultaneously, with the default convergence criteria (please reference to the CP2K online manual for further details: https://manual.cp2k.org/). Initially, a fully siliceous Silicalite-1 model including the TPA+OH- SDA in the channels intersection (derived from the crystallographic structure by ref. 52) was optimized, then two Si occupying T3 and/or T9 sites were substituted by Ti. Only the 22 unique Ti pair models generated by substitution were considered, as obtained from the CONFCNT routine53 implemented in the CRYSTAL17 code54. As Ti was inserted, an OH− anion from the nearest SDA molecule was manually coordinated to each Ti site, which is in line with the Parker and Millini hypothesis8. These structures were further optimized, and their relative energies ΔE, with respect to the most stable optimized model, were computed as from Eq. (1):1 ΔE=Ei−Emin

where Ei is the energy of the ith pair model, Emin is the lowest energy among those of all possible Ti pair models. The main geometrical parameters, the relative energies, and the derived Boltzmann probabilities for the 22 pairs featuring two Ti at T3/T9 sites are given in Supplementary Table 4.

A set of models where the SDA was removed has been considered too. As the geometrical features of these structures (Supplementary Table 5) are obviously comparable with refinement data presented in this work, this second dataset was exploited for the sake of the computational method validation. The adopted method and models reproduced quite well the structural parameters experimentally found for TS-1. As an indication, the average RMSD for the calculated cell parameters with respect to the refined ones for TS-1B (having the closest Ti loading) is 0.13 Å. Also, the average Ti-O bond length is in line with previous literature (see ref. 15 and references herein). RESP charges for Ti atoms were computed for all models and their average (+1.88) was adopted to rescale the value of f1 obtained via Kramers–Kronig transformation of XANES data.

Pair distribution function

Total scattering data for X-ray pair distribution function (PDF) analysis was collected at ESRF (European Synchrotron Radiation Facility, Grenoble, France), on the ID15A beamline55. An X-ray wavelength of 0.2480 Å (50 keV) ensured contrast between Ti and Si (~5 barns/atom) while enabling scattering measurement up to a maximum scattering vector Qmax = 22 Å−1. Powders were contained in 0.7 mm diameter polyimide tubes and spun during measurement to ensure powder averaging. The scattering of air and empty polyimide tube in the same conditions was collected for subsequent background subtraction. The scattered radiation was collected by a Dectris Pilatus 2 M CdTe detector sitting at 206 mm from the sample. Raw images were azimuthally integrated, accounting for detector flat-field, spatial distortion, and beam polarization using pyFAI56. The resulting total scattering patterns were background-subtracted and converted to G(r) using PDFgetx357. Fits to the experimental G(r) were done using Topas Academic v7 using the refined AXRD model as input structures.

Supplementary information

Supplementary Information for 'Quantitative locating titanium in the framework of titanium silicalite-1 by exploiting anomalous X-ray powder diffraction'

Peer Review File

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-51788-7.

Acknowledgements

P.R. thanks Dr. Ola Gjonnes Grendal from ID22 at ESRF for help with the Kramers–Kronig transformation of XANES data and Dr. Michal Andrzejewski from PSI for support with the ex situ capillary measurements. M.S., F.R., and S.B. acknowledge Evonik Resource Efficiency for providing samples investigated in this paper. M.S., F.R., and S.B. acknowledge Prof. Francesca Bonino and Prof. Valentina Crocellà for the fruitful discussion. M.S., F.R., and S.B. acknowledge support from the Project CH4.0 under the MUR program “Dipartimenti di Eccellenza 2023-2027” (CUP: D13C22003520001). M.S. thanks the C3S consortium for granting computational resources on the OCCAM cluster, funded by the Compagnia di San Paolo. P.R. and J.A.v.B. acknowledge the Paul Scherrer Institute, Villigen, Switzerland, for the provision of synchrotron radiation beamtime on the PHOENIX and MS beamlines. S.C. acknowledges the ESRF for the provision of in-house research beamtime under proposal number IH-HC-4004 (10.15151/ESRF-ES-1555168862). All authors thank the reviewers, Dr. Roberto Millini and Prof. Ettore Fois, for their inputs, which were instrumental in enhancing the robustness and credibility of this publication.

Author contributions

J.A.v.B. conceived the project, which was further developed by P.R., M.S., and S.B. P.R., M.S., and T.H. developed the methodology. P.R., M.S., T.H., S.C., and F.R. conducted the experiments. P.R., M.S., and S.C. performed the data analysis. All authors contributed to the data interpretation and writing the manuscript.

Peer review

Peer review information

Nature Communications thanks Ettore Fois, Roberto Millini and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Funding

Open access funding provided by Swiss Federal Institute of Technology Zurich.

Data availability

Data supporting the findings of this study are available within the article and the Supplementary Information file. Source data are provided as a Source Data file. All raw data generated during the study can be obtained from P.R. and M.S. Source data are provided with this paper.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Przemyslaw Rzepka, Matteo Signorile.
==== Refs
References

1. Taramasso, M., Perego, G. & Notari, B. Preparation of porous crystalline synthetic material comprised of silicon and titanium oxides. US patent 4,410,501 (1983).
2. Roffia, P., Padovan, M. & De Alberti, G. Catalytic process for preparing cyclohexanone-oxime. US patent 4,745,221 (1988).
3. Neri, C., Anfossi, B., Esposito, A. & Buonomo, F. Process for the epoxidation of olefinic compounds. US patent 4,833,260 21 (1984).
4. Bordiga S The structure of the peroxo species in the TS-1 catalyst as investigated by resonant Raman spectroscopy Angew. Chem. Int. Ed. 2002 41 4734 4737 10.1002/anie.200290032
Bordiga, S. et al. The structure of the peroxo species in the TS-1 catalyst as investigated by resonant Raman spectroscopy. Angew. Chem. Int. Ed. 41, 4734–4737 (2002).10.1002/anie.200290032
5. Cavani F Liquid phase oxidation via heterogeneous catalysis. Organic synthesis and industrial applications. Edited by Mario G. Clerici and Oxana A. Kholdeeva. Angew. Chem. Int. Ed. 2014 53 7707 7707 10.1002/anie.201404813
Cavani, F. Liquid phase oxidation via heterogeneous catalysis. Organic synthesis and industrial applications. Edited by Mario G. Clerici and Oxana A. Kholdeeva. Angew. Chem. Int. Ed. 53, 7707–7707 (2014).10.1002/anie.201404813
6. Baerlocher, C. & McCusker, L. B. Database of zeolite structures. http://www.iza-structure.org/databases/ (2008).
7. Gordon CP Efficient epoxidation over dinuclear sites in titanium silicalite-1 Nature 2020 586 708 713 10.1038/s41586-020-2826-3 33116285
Gordon, C. P. et al. Efficient epoxidation over dinuclear sites in titanium silicalite-1. Nature 586, 708–713 (2020).33116285 10.1038/s41586-020-2826-3
8. Parker WO Millini R Ti coordination in titanium silicalite-1 J. Am. Chem. Soc. 2006 128 1450 1451 10.1021/ja0576785 16448106
Parker, W. O. & Millini, R. Ti coordination in titanium silicalite-1. J. Am. Chem. Soc. 128, 1450–1451 (2006).16448106 10.1021/ja0576785
9. Lätsch L NMR signatures and electronic structure of Ti sites in titanosilicalite-1 from solid-state 47/49Ti NMR spectroscopy J. Am. Chem. Soc. 2023 145 15018 15023 10.1021/jacs.2c09867 37418311
Lätsch, L. et al. NMR signatures and electronic structure of Ti sites in titanosilicalite-1 from solid-state 47/49Ti NMR spectroscopy. J. Am. Chem. Soc. 145, 15018–15023 (2023).37418311 10.1021/jacs.2c09867
10. Ricchiardi G Vibrational structure of titanium silicate catalysts. A spectroscopic and theoretical study J. Am. Chem. Soc. 2001 123 11409 11419 10.1021/ja010607v 11707118
Ricchiardi, G. et al. Vibrational structure of titanium silicate catalysts. A spectroscopic and theoretical study. J. Am. Chem. Soc. 123, 11409–11419 (2001).11707118 10.1021/ja010607v
11. Su J Amorphous Ti species in titanium silicalite-1: Structural features, chemical properties, and inactivation with sulfosalt J. Catal. 2012 288 1 7 10.1016/j.jcat.2011.12.006
Su, J. et al. Amorphous Ti species in titanium silicalite-1: Structural features, chemical properties, and inactivation with sulfosalt. J. Catal. 288, 1–7 (2012).10.1016/j.jcat.2011.12.006
12. Guo Q Feng Z Li G Fan F Li C Finding the ‘missing components’ during the synthesis of TS-1 zeolite by UV resonance raman spectroscopy J. Phys. Chem. C. 2013 117 2844 2848 10.1021/jp310900a
Guo, Q., Feng, Z., Li, G., Fan, F. & Li, C. Finding the ‘missing components’ during the synthesis of TS-1 zeolite by UV resonance raman spectroscopy. J. Phys. Chem. C. 117, 2844–2848 (2013).10.1021/jp310900a
13. Signorile M Titanium defective sites in TS-1: structural insights by combining spectroscopy and simulation Angew. Chem. Int. Ed. 2020 59 18145 18150 10.1002/anie.202005841
Signorile, M. et al. Titanium defective sites in TS-1: structural insights by combining spectroscopy and simulation. Angew. Chem. Int. Ed. 59, 18145–18150 (2020).10.1002/anie.202005841
14. Signorile M Effect of Ti speciation on catalytic performance of TS-1 in the hydrogen peroxide to propylene oxide reaction J. Phys. Chem. C. 2018 122 9021 9034 10.1021/acs.jpcc.8b01401
Signorile, M. et al. Effect of Ti speciation on catalytic performance of TS-1 in the hydrogen peroxide to propylene oxide reaction. J. Phys. Chem. C. 122, 9021–9034 (2018).10.1021/acs.jpcc.8b01401
15. Signorile M Computational assessment of relative sites stabilities and site-specific adsorptive properties of titanium silicalite-1 J. Phys. Chem. C. 2018 122 1612 1621 10.1021/acs.jpcc.7b10104
Signorile, M. et al. Computational assessment of relative sites stabilities and site-specific adsorptive properties of titanium silicalite-1. J. Phys. Chem. C. 122, 1612–1621 (2018).10.1021/acs.jpcc.7b10104
16. Dong J Toward a unified identification of Ti location in the MFI framework of high-Ti-loaded TS-1: combined EXAFS, XANES, and DFT study J. Phys. Chem. C. 2016 120 20114 20124 10.1021/acs.jpcc.6b05087
Dong, J. et al. Toward a unified identification of Ti location in the MFI framework of high-Ti-loaded TS-1: combined EXAFS, XANES, and DFT study. J. Phys. Chem. C. 120, 20114–20124 (2016).10.1021/acs.jpcc.6b05087
17. Ricchiardi G De Man A Sauer J The effect of hydration on structure and location of Ti-sites in Ti- silicalite catalysts. A computational study Phys. Chem. Chem. Phys. 2000 2 2195 2204 10.1039/a909617e
Ricchiardi, G., De Man, A. & Sauer, J. The effect of hydration on structure and location of Ti-sites in Ti- silicalite catalysts. A computational study. Phys. Chem. Chem. Phys. 2, 2195–2204 (2000).10.1039/a909617e
18. Deka RC Comparison of all sites for Ti substitution in zeolite TS-1 by an accurate embedded-cluster method J. Phys. Chem. B 2005 109 24304 24310 10.1021/jp050056l 16375428
Deka, R. C. et al. Comparison of all sites for Ti substitution in zeolite TS-1 by an accurate embedded-cluster method. J. Phys. Chem. B 109, 24304–24310 (2005).16375428 10.1021/jp050056l
19. Gale JD A periodic density functional study of the location of titanium within TS-1 Solid State Sci. 2006 8 234 240 10.1016/j.solidstatesciences.2006.02.011
Gale, J. D. A periodic density functional study of the location of titanium within TS-1. Solid State Sci. 8, 234–240 (2006).10.1016/j.solidstatesciences.2006.02.011
20. Marra GL Artioli G Fitch AN Milanesio M Lamberti C Orthorhombic to monoclinic phase transition in high-Ti-loaded TS-1: an attempt to locate Ti in the MFI framework by low temperature XRD Microporous Mesoporous Mater. 2000 40 85 94 10.1016/S1387-1811(00)00244-4
Marra, G. L., Artioli, G., Fitch, A. N., Milanesio, M. & Lamberti, C. Orthorhombic to monoclinic phase transition in high-Ti-loaded TS-1: an attempt to locate Ti in the MFI framework by low temperature XRD. Microporous Mesoporous Mater. 40, 85–94 (2000).10.1016/S1387-1811(00)00244-4
21. Lamberti C Structural characterization of Ti-Silicalite-1: a synchrotron radiation X-ray powder diffraction study J. Catal. 1999 183 222 231 10.1006/jcat.1999.2403
Lamberti, C. et al. Structural characterization of Ti-Silicalite-1: a synchrotron radiation X-ray powder diffraction study. J. Catal. 183, 222–231 (1999).10.1006/jcat.1999.2403
22. Lamberti C Ti location in the MFI framework of Ti-silicalite-1: a neutron powder diffraction study J. Am. Chem. Soc. 2001 123 2204 2212 10.1021/ja003657t 11456866
Lamberti, C. et al. Ti location in the MFI framework of Ti-silicalite-1: a neutron powder diffraction study. J. Am. Chem. Soc. 123, 2204–2212 (2001).11456866 10.1021/ja003657t
23. Hijar CA The siting of Ti in TS-1 is non-random. Powder neutron diffraction studies and theoretical calculations of TS-1 and FeS-1 J. Phys. Chem. B 2000 104 12157 12164 10.1021/jp002167k
Hijar, C. A. et al. The siting of Ti in TS-1 is non-random. Powder neutron diffraction studies and theoretical calculations of TS-1 and FeS-1. J. Phys. Chem. B 104, 12157–12164 (2000).10.1021/jp002167k
24. Henry PF Weller MT Wilson CC Structural investigation of TS-1: determination of the true nonrandom titanium framework substitution and silicon vacancy distribution from powder neutron diffraction studies using isotopes J. Phys. Chem. B 2001 105 7452 7458 10.1021/jp0107715
Henry, P. F., Weller, M. T. & Wilson, C. C. Structural investigation of TS-1: determination of the true nonrandom titanium framework substitution and silicon vacancy distribution from powder neutron diffraction studies using isotopes. J. Phys. Chem. B 105, 7452–7458 (2001).10.1021/jp0107715
25. Cromer DT Liberman D Relativistic calculation of anomalous scattering factors for X rays J. Chem. Phys. 1970 53 1891 1898 10.1063/1.1674266
Cromer, D. T. & Liberman, D. Relativistic calculation of anomalous scattering factors for X rays. J. Chem. Phys. 53, 1891–1898 (1970).10.1063/1.1674266
26. Grenier S Joly Y Basics of resonant elastic X-ray scattering theory J. Phys. Conf. Ser. 2014 519 12001 10.1088/1742-6596/519/1/012001
Grenier, S. & Joly, Y. Basics of resonant elastic X-ray scattering theory. J. Phys. Conf. Ser. 519, 12001 (2014).10.1088/1742-6596/519/1/012001
27. Hendrickson WA Anomalous diffraction in crystallographic phase evaluation Q. Rev. Biophys. 2014 47 49 93 10.1017/S0033583514000018 24726017
Hendrickson, W. A. Anomalous diffraction in crystallographic phase evaluation. Q. Rev. Biophys. 47, 49–93 (2014).24726017 10.1017/S0033583514000018
28. Finkelstein KD Shen Q Shastri S Resonant X-ray diffraction near the iron K edge in hematite (α-Fe2O3) Phys. Rev. Lett. 1992 69 1612 1615 10.1103/PhysRevLett.69.1612 10046266
Finkelstein, K. D., Shen, Q. & Shastri, S. Resonant X-ray diffraction near the iron K edge in hematite (α-Fe2O3). Phys. Rev. Lett. 69, 1612–1615 (1992).10046266 10.1103/PhysRevLett.69.1612
29. Pinar AB Pinpointing and quantifying the aluminum distribution in zeolite catalysts using anomalous scattering at the Al absorption edge J. Am. Chem. Soc. 2021 143 17926 17930 10.1021/jacs.1c06925 34695360
Pinar, A. B. et al. Pinpointing and quantifying the aluminum distribution in zeolite catalysts using anomalous scattering at the Al absorption edge. J. Am. Chem. Soc. 143, 17926–17930 (2021).34695360 10.1021/jacs.1c06925
30. Goldbach, A., Saboungi, M. L., Iton, L. & Price, D. L. Stabilization of selenium in zeolites: an anomalous X-ray scattering study. Chem. Commun. 10.1039/a900253g(1999).
31. Pichon C In situ characterisation by anomalous X-ray diffraction of the cationic distribution of dehydrated SrRbX zeolite Stud. Surf. Sci. Catal. 2004 154 1641 1648 10.1016/S0167-2991(04)80689-X
Pichon, C. et al. In situ characterisation by anomalous X-ray diffraction of the cationic distribution of dehydrated SrRbX zeolite. Stud. Surf. Sci. Catal. 154, 1641–1648 (2004).10.1016/S0167-2991(04)80689-X
32. Jones RH Lightfoot P Ormerod RM The location of sorbed Kr in silicalite using resonant X-ray diffraction J. Phys. Chem. Solids 1995 56 1377 1381 10.1016/0022-3697(95)00073-9
Jones, R. H., Lightfoot, P. & Ormerod, R. M. The location of sorbed Kr in silicalite using resonant X-ray diffraction. J. Phys. Chem. Solids 56, 1377–1381 (1995).10.1016/0022-3697(95)00073-9
33. Hendrickson WA Determination of macromolecular structures from anomalous diffraction of synchrotron radiation Science 1991 254 51 58 10.1126/science.1925561 1925561
Hendrickson, W. A. Determination of macromolecular structures from anomalous diffraction of synchrotron radiation. Science 254, 51–58 (1991).1925561 10.1126/science.1925561
34. Henke BL Gullikson EM Davis JC X-ray interactions: photoabsorption, scattering, transmission, and reflection at E = 50-30, 000 eV, Z = 1-92. At. Data Nucl. Data Tables 1993 54 181 342 10.1006/adnd.1993.1013
Henke, B. L., Gullikson, E. M. & Davis, J. C. X-ray interactions: photoabsorption, scattering, transmission, and reflection at E = 50-30, 000 eV, Z = 1-92.At. Data Nucl. Data Tables 54, 181–342 (1993).10.1006/adnd.1993.1013
35. Evans G Pettifer RF CHOOCH: a program for deriving anomalous-scattering factors from X-ray fluorescence spectra J. Appl. Crystallogr. 2001 34 82 86 10.1107/S0021889800014655
Evans, G. & Pettifer, R. F. CHOOCH: a program for deriving anomalous-scattering factors from X-ray fluorescence spectra. J. Appl. Crystallogr. 34, 82–86 (2001).10.1107/S0021889800014655
36. Beek Untangling diffraction intensity: modulation enhanced diffraction on ZrO2 powder J. Appl. Crystallogr. 2012 45 738 747 10.1107/S0021889812018109
Beek et al. Untangling diffraction intensity: modulation enhanced diffraction on ZrO2 powder. J. Appl. Crystallogr. 45, 738–747 (2012).10.1107/S0021889812018109
37. Owen RL Holton JM Schulze-Briese C Garman EF Determination of X-ray flux using silicon pin diodes J. Synchrotron Radiat. 2009 16 143 151 10.1107/S0909049508040429 19240326
Owen, R. L., Holton, J. M., Schulze-Briese, C. & Garman, E. F. Determination of X-ray flux using silicon pin diodes. J. Synchrotron Radiat. 16, 143–151 (2009).19240326 10.1107/S0909049508040429
38. Coelho AA TOPAS and TOPAS-academic: an optimization program integrating computer algebra and crystallographic objects written in C++ J. Appl. Crystallogr. 2018 51 210 218 10.1107/S1600576718000183
Coelho, A. A. TOPAS and TOPAS-academic: an optimization program integrating computer algebra and crystallographic objects written in C++. J. Appl. Crystallogr. 51, 210–218 (2018).10.1107/S1600576718000183
39. Willmott PR The Materials Science beamline upgrade at the Swiss light source J. Synchrotron Radiat. 2013 20 667 682 10.1107/S0909049513018475 23955029
Willmott, P. R. et al. The Materials Science beamline upgrade at the Swiss light source. J. Synchrotron Radiat. 20, 667–682 (2013).23955029 10.1107/S0909049513018475
40. Millini R Previde Massara E Perego G Bellussi G Framework composition of titanium silicalite-1 J. Catal. 1992 137 497 503 10.1016/0021-9517(92)90176-I
Millini, R., Previde Massara, E., Perego, G. & Bellussi, G. Framework composition of titanium silicalite-1. J. Catal. 137, 497–503 (1992).10.1016/0021-9517(92)90176-I
41. Efron B Tibshirani R Bootstrap methods for standard errors, confidence intervals, and other measures of statistical accuracy Stat. Sci. 1986 1 54 75
Efron, B. & Tibshirani, R. Bootstrap methods for standard errors, confidence intervals, and other measures of statistical accuracy. Stat. Sci. 1, 54–75 (1986).
42. DiCiccio TJ Efron B Bootstrap confidence intervals Stat. Sci. 1996 11 189 212 10.1214/ss/1032280214
DiCiccio, T. J. & Efron, B. Bootstrap confidence intervals. Stat. Sci. 11, 189–212 (1996).10.1214/ss/1032280214
43. Gamba A Tabacchi G Fois E TS-1 from first principles J. Phys. Chem. A 2009 113 15006 15015 10.1021/jp905110s 19785451
Gamba, A., Tabacchi, G. & Fois, E. TS-1 from first principles. J. Phys. Chem. A 113, 15006–15015 (2009).19785451 10.1021/jp905110s
44. Román-Román EI Zicovich-Wilson CM The role of long-range van der Waals forces in the relative stability of SiO2-zeolites Chem. Phys. Lett. 2015 619 109 114 10.1016/j.cplett.2014.11.044
Román-Román, E. I. & Zicovich-Wilson, C. M. The role of long-range van der Waals forces in the relative stability of SiO2-zeolites. Chem. Phys. Lett. 619, 109–114 (2015).10.1016/j.cplett.2014.11.044
45. Zhang KYJ Main P Histogram matching as a new density modification technique for phase refinement and extension of protein molecules Acta Crystallogr. Sect. A 1990 46 41 46 10.1107/S0108767389009311
Zhang, K. Y. J. & Main, P. Histogram matching as a new density modification technique for phase refinement and extension of protein molecules. Acta Crystallogr. Sect. A 46, 41–46 (1990).10.1107/S0108767389009311
46. Hutter J Iannuzzi M Schiffmann F Vandevondele J Cp2k: Atomistic simulations of condensed matter systems Wiley Interdiscip. Rev. Comput. Mol. Sci. 2014 4 15 25 10.1002/wcms.1159
Hutter, J., Iannuzzi, M., Schiffmann, F. & Vandevondele, J. Cp2k: Atomistic simulations of condensed matter systems. Wiley Interdiscip. Rev. Comput. Mol. Sci. 4, 15–25 (2014).10.1002/wcms.1159
47. Perdew JP Burke K Ernzerhof M Generalized gradient approximation made simple Phys. Rev. Lett. 1996 77 3865 3868 10.1103/PhysRevLett.77.3865 10062328
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).10062328 10.1103/PhysRevLett.77.3865
48. Goedecker S Teter M Separable dual-space Gaussian pseudopotentials Phys. Rev. B - Condens. Matter Mater. Phys. 1996 54 1703 1710 10.1103/PhysRevB.54.1703
Goedecker, S. & Teter, M. Separable dual-space Gaussian pseudopotentials. Phys. Rev. B - Condens. Matter Mater. Phys. 54, 1703–1710 (1996).10.1103/PhysRevB.54.1703
49. Lippert G Hutter J Parrinello M A hybrid Gaussian and plane wave density functional scheme Mol. Phys. 1997 92 477 488 10.1080/00268979709482119
Lippert, G., Hutter, J. & Parrinello, M. A hybrid Gaussian and plane wave density functional scheme. Mol. Phys. 92, 477–488 (1997).10.1080/00268979709482119
50. VandeVondele J Hutter J Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases J. Chem. Phys. 2007 127 114105 10.1063/1.2770708 17887826
VandeVondele, J. & Hutter, J. Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. J. Chem. Phys. 127, 114105 (2007).17887826 10.1063/1.2770708
51. VandeVondele J Hutter J An efficient orbital transformation method for electronic structure calculations J. Chem. Phys. 2003 118 4365 4369 10.1063/1.1543154
VandeVondele, J. & Hutter, J. An efficient orbital transformation method for electronic structure calculations. J. Chem. Phys. 118, 4365–4369 (2003).10.1063/1.1543154
52. Van Koningsveld H Van Bekkum H Jansen JC On the location and disorder of the tetrapropylammonium (TPA) ion in zeolite ZSM‐5 with improved framework accuracy Acta Crystallogr. Sect. B 1987 43 127 132 10.1107/S0108768187098173
Van Koningsveld, H., Van Bekkum, H. & Jansen, J. C. On the location and disorder of the tetrapropylammonium (TPA) ion in zeolite ZSM‐5 with improved framework accuracy. Acta Crystallogr. Sect. B 43, 127–132 (1987).10.1107/S0108768187098173
53. Mustapha S On the use of symmetry in configurational analysis for the simulation of disordered solids J. Phys. Condens. Matter 2013 25 105401 105417 10.1088/0953-8984/25/10/105401 23388579
Mustapha, S. et al. On the use of symmetry in configurational analysis for the simulation of disordered solids. J. Phys. Condens. Matter 25, 105401–105417 (2013).23388579 10.1088/0953-8984/25/10/105401
54. Dovesi R Quantum-mechanical condensed matter simulations with CRYSTAL Wiley Interdiscip. Rev. Comput. Mol. Sci. 2018 8 e1360 10.1002/wcms.1360
Dovesi, R. et al. Quantum-mechanical condensed matter simulations with CRYSTAL. Wiley Interdiscip. Rev. Comput. Mol. Sci. 8, e1360 (2018).10.1002/wcms.1360
55. Vaughan GBM ID15A at the ESRF-a beamline for high speed operando X-ray diffraction, diffraction tomography and total scattering J. Synchrotron Rad. 2020 27 515 528 10.1107/S1600577519016813
Vaughan, G. B. M. et al. ID15A at the ESRF-a beamline for high speed operando X-ray diffraction, diffraction tomography and total scattering. J. Synchrotron Rad. 27, 515–528 (2020).10.1107/S1600577519016813
56. Kieffer J Valls V Blanc N Hennig C New tools for calibrating diffraction setups J. Synchrotron Rad. 2020 27 558 566 10.1107/S1600577520000776
Kieffer, J., Valls, V., Blanc, N. & Hennig, C. New tools for calibrating diffraction setups. J. Synchrotron Rad. 27, 558–566 (2020).10.1107/S1600577520000776
57. Juhás P Davis T Farrow CL Billinge SJL PDFgetX3: a rapid and highly automatable program for processing powder diffraction data into total scattering pair distribution functions J. Appl. Crystallogr. 2013 46 560 566 10.1107/S0021889813005190
Juhás, P., Davis, T., Farrow, C. L. & Billinge, S. J. L. PDFgetX3: a rapid and highly automatable program for processing powder diffraction data into total scattering pair distribution functions. J. Appl. Crystallogr. 46, 560–566 (2013).10.1107/S0021889813005190
