[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Poncirus trifoliata Aqueous Extract Protects Cardiomyocytes against Doxorubicin-Induced Toxicity through Upregulation of NAD(P)H Dehydrogenase Quinone Acceptor Oxidoreductase 1
Next Article in Special Issue
Morphology Regulation of Zeolite MWW via Classical/Nonclassical Crystallization Pathways
Previous Article in Journal
Effect of Drying Methods on Chemical and Sensory Properties of Cannabis sativa Leaves
Previous Article in Special Issue
Efficient Selective Capture of Carbon Dioxide from Nitrogen and Methane Using a Metal-Organic Framework-Based Nanotrap
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 24
10.3390/molecules28248087
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulating Crystal Structure, Acidity, Proton Distribution, and IR Spectra of Acid Zeolite HSAPO-34: A High Accuracy Study

by
Xiaofang Chen
* and
Tie Yu
Institute of Molecular Sciences and Engineering, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, China
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(24), 8087; https://doi.org/10.3390/molecules28248087
Submission received: 6 November 2023 / Revised: 4 December 2023 / Accepted: 5 December 2023 / Published: 14 December 2023
(This article belongs to the Special Issue Zeolites and Porous Materials: Synthesis, Properties and Applications)
Browse Figures
Graphical abstract
"> Figure 1
<p>Single silicon-doped alumino-phosphates (HSAPO-34): (<b>a</b>) the single-unit cell (1UC) with rhombohedral representation, (<b>b</b>) the expanded lattice structure. Note that the proton is omitted for simplicity.</p> "> Figure 2
<p>The total energy depends on the lattice volume. The proton binding sites of HSAPO-34 are (<b>a</b>) O(1), (<b>b</b>) O(2), (<b>c</b>) O(3), and (<b>d</b>) O(4).</p> "> Figure 3
<p>Schematic diagram for four possible Brönsted acid sites of HSAPO-34. The proton binding sites are (<b>a</b>) O1, (<b>b</b>) O2, (<b>c</b>) O3, and (<b>d</b>) O4. Noting that Al in cycn, Si in blue, P in purple, O in red, and the number in figure is the numbering of O atom.</p> "> Figure 4
<p>The proton distribution of HSAPO-34 depends on the temperature.</p> "> Figure 5
<p>Bader charges on AlO<sub>3</sub>-O(H)-SiO<sub>3</sub> fragment from four possible Brönsted acid sites of HSAPO-34. The proton binding sites are (<b>a</b>) O(1), (<b>b</b>) O(2), (<b>c</b>) O(3), and (<b>d</b>) O(4). Noting that Al in cyan, Si in blue, and O in red.</p> "> Figure 6
<p>Calculated ensemble-averaged DPE (&lt;DPE&gt;) of HSAPO-34 depends on the temperature.</p> ">
Versions Notes

Abstract

:
It is a challenge to characterize the acid properties of microporous materials in either experiments or theory. This study presents the crystal structure, acid site, acid strength, proton siting, and IR spectra of HSAPO-34 from the SCAN + rVV10 method. The results indicate: the crystal structures of various acid sites of HSAPO-34 deviate from the space group of R 3 ¯ ; the acid strength inferred from the DPE value likely decreases with the proton binding sites at O(2), O(4), O(1),and O(3), contrary to the stability order in view of the internal energy; the calculated ensemble-averaged DPE is about 1525 kJ/mol at 673.15 K; and the proton siting and the proton distribution are distinctly influenced by the temperature: at low temperatures, the proton is predominantly located at O(3), while it prefers O(2) at high temperatures, and the proton at O(4) assumedly has the least distribution at 273.15–773.15 K. In line with the neutron diffraction experiment, a correction factor of 0.979 is needed to correct for the calculated hydroxyl stretching vibration (ν(O-H)) of HSAPO-34. It seems that the SCAN meta-GGA method, compensating for some drawbacks of the GGA method, could provide satisfying results regarding the acid properties of HSAPO-34.

Graphical Abstract">

Graphical Abstract

1. Introduction

Acidic silicon-doped silicoaluminophosphates (abbreviated as SAPOs) are developed from the famous microporous material silicoaluminophosphates (AlPOs), produced through the SII and SIII mechanisms, where some P atoms are substituted by the association of Si and H atoms, and P and Al atoms are substituted by two Si atoms, respectively [1,2,3,4,5]. Due to the excellent performance of the catalysis and the adsorption, SAPOs have widespread applications, especially in the field of methanol conversion into olefins (MTO), CO2 capture and conversion, NH3-SCR, and the DeNOx reaction [5,6,7,8,9,10,11,12,13,14,15,16].
The performance of zeolite is generally determined by the acid site, the acid amount, the acid strength, and accessibility to the acid site [17,18,19,20,21]. Most SAPOs are structurally constructed by isomorphous Si4+ displacing P5+, and the formed negative charge is effectively compensated for by the proton bound to the framework oxygen or metal cation [4,22] in the channels/cages. It was found that the P sites in SAPOs were substituted exclusively by Si sites [3]. The acid zeolite HSAPO-34 has the acidic center of SiO(H)Al. The compensatory proton is usually regarded as the intrinsic origin of the Brönsted acid of HSAPO-34. The incorporated Si atom directly changes the electrostatic potential of the adjacent O atoms, and then makes the O atom unequivalent. Thus, the acid properties of HSAPO-34 are largely determined by the nature of the incorporated Si species [23]. It seems feasible to identify the Brönsted acid site by examining the chemical environment of the oxygen atom chemically bound to the incorporated Si atom. The acidic single silicon-substituted SAPO-34 zeolite (HSAPO-34) has four unequivalent O atoms (i.e., O(1)-O(4)) chemically bound to the incorporated Si atom. It is expected that HSAPO-34 will have four possible Brönsted acid sites. The acid properties of HSAPO-34 are also influenced by the concentration of Si islands and the proton distribution [23,24].
Some acid properties of HSAPO-34 were studied by experimental technologies, such as neutron diffraction [2,3,25], probe molecule (for instance, NH3, CO, and C2H4) adsorb–desorption [4,11,19,26], IR [26,27] and NMR measurements [23,28], and transmission electron microscopes (TEMs) [11]. An isomorph of HSAPO-34, a chabazite specimen, was found to have the trigonal unit cell with the space group of R 3 ¯ m (or nearly so) and alattice parameter of a = 9.40 Å and 94.30° [25]. The infraredmass spectrometry/temperature-programmed desorption (IRMS-TPD) of ammonia indicated that HSAPO-34 had a slightly weaker acid strength than chabazite, and the primary Brönsted OH was on the edge of the Si domain (island) [4]. Among a series of hierarchical SAPO-34 samples, the nanostructure of the SAPO-34 sample with a lower acid concentration and weaker acid strength (i.e., a sample with a lower intensity in the NH3-TPD profiles) exhibits the appropriate catalysis performance for MTO reactions [11]. Two distinct peaks were identified as the signal of two bridging hydroxyl Brönsted sites of the crystalline HSAPO-34 catalyst [3,23], and the proton at O(4) likely had a greater acidity than that at O(2) [3]. It was assumed that the protons preferred to reside on the framework O(2) and O(4) atoms in the dehydrated [3] SAPO-34 sample, while the O(1) atom was coordinated by a hydroxonium ion (H3O+), and yet, the hydrogen bond formed between O(2) and a water molecule in the hydrated sample [2]. Meanwhile, three distinct acid sites of HSAPO-34 were further evidenced by FTIR spectra at 3631, 3617, and 3600 cm−1 [26]. The Si mole fraction and the particle size almost had no influence on the acid strength of HSAPO-34 [29], and the metal content had a strong influence on the concentration of acid sites [30].
The acid zeolite HSAPO-34 sample is usually synthesized under hydrothermal conditions. The high synthesis reactivity makes it difficult to monitor the crystal nucleus formation [31,32,33]. The proton exchanges also proceed fast between different adsorption sites [28]. This leads to the HSAPO-34 sample in the experiment being a mixture of various Brönsted acid sites. HSAPO-34 belongs to the small-pore materials, usually with the bore diameter of about 3–12 Å [5]. Some acid sites on the zeolitic surface are easily accessible, while the others, especially in the narrow channel or the inner cavity, are difficult to reach. The fast proton exchange, as well as the inaccessibility, also increases the challenge of characterizing the precise acid properties of SAPOs in the experiment.
The acid zeolite HSAPO-34 has also been investigated by diverse theories. The plane wave periodic gradient-corrected density functional theory (PDFT) is the favorite of researchers to simulate the structure, property, and reactivity of the bulk [34]. The general gradient approximation (GGA) PW91 functional was once used to calculate the exchange correlation energy; O(1) and O(3) were suggested to be two favored proton binding sites of HSAPO-34 [35,36]. In comparison with neutron diffraction results [2,3,28], it has been concluded that the GGA PW91 functional was capable of predicting the lattice parameter and the ion position at a certain level of precision, yet it failed to fully reproduce the acid properties of HSAPO-34 [36]. The divergences between the experiment and the theory are probably caused by the theoretical model’s simplicity, the experimental complexity, and the low energy difference between various Brönsted acid sites of HSAPO-34.
Besides the gradient-corrected correlation PW91 functional [37], PBE is another popular GGA functional [38]. As a simplification of PW91, PBE is capable of yielding numerical results almost identical to PW91 in most usual systems. The nonequivalence of PBE and PW91 has also been observed in the presence of surface effects; and PBE is generally superior to PW91 in treating these weak interactions [39]. During a plane wave theory simulation, the hierarchy of DFT approximations is adopted to resolve the Kohn–Sham equation through the displacement of the genuine potential by the relatively flexible pseudo potential [40,41,42,43,44,45]. The generalized gradient approximation is the second rung of Jacob’s ladder [46], which treats the density at a given point changeably [34]. It is desirable to use the suitable pseudo potential to maximally reproduce the crystal structure and the energy of the real bulk [47]. To compensate for the lack of kinetic energy, the meta generalized gradient approximation (meta-GGA) appears to be another wave plane approximation higher than GGA [48,49,50,51]. Recently, it was reported that the strongly constrained and appropriately normed (SCAN) approximation had the ability to evaluate solid stability and reactivity, although it demanded the long computation time and expensive calculation resources [48,49,50,51,52].
In this study, we attempt to use a higher-level approximation method (i.e., SCAN + rVV10) to simulate the fine crystal structure of all possible Brönsted acid sites of HSAPO-34 and analyze the acidity strength and the proton siting. The preferred proton binding site will be further analyzed in terms of the quantitative proton distribution of HSAPO-34. This study will shed light on the difficulties of accurately simulating the acid properties of porous materials.

2. Results and Discussion

2.1. Unit Cell

AlPO4-34 is structurally analogous to an isomorph, a chabazite. Both have CHA topology and yet have different element compositions. Following the SMIII mechanism, the all-silica chabazite is converted to AlPO4-34 by the alternative displacement of Si by Al and P atoms. The space group reduces from R 3 ¯ m in the all-silica CHA framework to R 3 ¯ in the AlPO4-34 framework. It was found that the silicon atom in SAPOs was substituted exclusively by P sites and the formed negative charge was compensated for by the proton, leading to the acid zeolite HSAPO-34 [3]. Chabazite has the primary building units (PBU) of [Si-O-Al] and [Si-O-Si] linkages, while SAPO-34 has [Si-O-Al], [Si-O-Si], and [Al-O-P] linkages. The rhombohedral representation of HSAPO-34 in a single-unit cell (1UC) has the chemical formula of HAl6P5Si1O24 (Figure 1). Some atoms in HSAPO-34 are structurally equivalent. The tetrahedral units of AlO4 and PO4 are coupled in HSAPO-34 crystalline. The Al-O-P bonds are formed and the [O-Al-O-P]n chemical network chains appear in the crystalline of HSAPO-34. The acid zeolite HSAPO-34 has a narrow eight-member ring (8MR) (3.8 Å × 3.8 Å) perpendicular to the crystalline ac or bc plane, and at welve-member ring (12MR) along the side face of the eight-member ring [53].
Apart from the newly added H atom, HSAPO-34 has a crystal structure similar to the AlPO4-34 framework. There are four unequivalent O atoms around the incorporated Si atom, denoted as O(x) (x = 1–4) (Figure 1). Any of the oxygen atoms which are chemically bound to the zeolitic substituted tetrahedral (T) sites could serve as the attached site of the Brönsted acid proton [54]. Each unequivalent O atom is shared by three windows: O(1), two 4MRs and one 8MR; O(2), one 8MR, one 6MR, and one 4MR; O(3), two 4MRs and one 6MR; O(4), two 8MRs and one 4MR. The O atom numbering in this study is the same as that mentioned by Smith et al. [3] and Jeanvoine et al. [36], which differs from the database of zeolite structures in the internal zeolite associate (IZA) [53] (Table 1). The four unequivalent O atoms are around the incorporated Si atom, resulting in four possible Brönsted acid sites of HSAPO-34, denoted as HSAPO-34-O(x), x = 1–4.

2.2. Calculated Crystal Structure

Generally, the Birch–Murnaghan equation [55] describes the total energy (E(V)) of the bulk dependence of the volume (V), the equilibrium volume (V0), the equilibrium electronic energy (E0), the bulk modulus (B0), and its derivative with respect to pressure ( B 0 ′ ) . On the basis of the E-V data, the total energy and the crystal volume in the equilibrium state could be deduced by fitting the E-V curve with the third order Birch–Murnaghan equation (Equation (1)).
E V = E 0 + 9 V 0 B 0 16 V 0 V 2 3 − 1 3 B 0 ′ + V 0 V 2 3 − 1 2 6 − 4 V 0 V 2 3
Scanning of the E-V curve is performed to search for the crystal structure in the equilibrium state. Figure 2 illustrates the total energy dependence of the lattice volume for four possible Brönsted acid sites of HSAPO-34. The constant volume selected is from 750 Å3 to 900 Å3 with an interval of 10 Å3. Thus, each E-V curve scan needs 16 structures at a constant volume, which are used to partially optimize. The partial optimization means that the crystal structure is optimized by relaxing the cell shape and the ion position but fixing the cell volume. The atomic force and the stress tensor are calculated. As shown in Figure 2, each curve corresponds to one Brönsted acid site of HSAPO-34, and all exhibit the parabola shape. The total energy of each curve at 750 Å3 is always lower than that at 900 Å3.
The fitted parameters of E0, V0, B0, and B 0 ′ for each Brönsted acid site of HSAPO-34 are provided in Table 2. When the proton is chemically bound to the O(1) or O(3) atom of the acidic zeolite HSAPO-34, the equilibrium volume (V0) of the Brönsted acid site is equal to 821.82 Å3 or 820.13 Å3. Previous neutron diffraction experiments detected that the lattice volume of HSAPO-34 was 822.24 Å3 for the dehydrated sample [3] and 820.20 Å3 for the hydrated sample [2]. In comparison, two Brönsted acid sites, HSAPO-34-O(1) and HSAPO-34-O(3), have lattice volumes that are the most close to the experimental value. When the proton rests on the O(2) or O(4) atom, the corresponding Brönsted acid site has a negative lattice volume deviation from either the dehydrated HSAPO-34 (−11.26 and −29.06 Å3) [3] or hydrated sample (−9.22 and −27.02 Å3) [2]. It is worthwhile noting that the lattice volume detected by the neutron diffraction experiment is generally related to the Brönsted acid site mixture of the acid zeolite HSAPO-34 in the equilibrium state.
To obtain the more precise crystal structure of HSAPO-34, full optimization is carried out by the period DFT, implemented with the strongly constrained and appropriately normed (SCAN) method. The atomic force and the stress tensor are calculated by relaxing the volume and the ion. Table 3 lists the calculated lattice parameters for four Brönsted acid sites of the acidic zeolite HSAPO-34. The dehydrated HSAPO-34 [3] is selected as the reference. For the detailed atomic coordinates of the calculated crystal structure of HSAPO-34, please refer to Supplementary Materials.
As shown in Table 3, three crystal edges (i.e., a, b and c) from each Brönsted acid site of HSAPO-34 are not completely equal to each other. Taking HSAPO-34-O(1) as an example, the crystal edge of a is calculated to be 9.40 Å, equal to the deduced experimental value of dehydrated HSAPO-34 [3], while the crystal edge of b increases by 0.09 Å and the crystal edge of c decreases by 0.14 Å. The crystal angle of α of HSAPO-34–O(1) is the calculated value the most close to the deduced experimental value. Among four Brönsted acid sites, HSAPO-34-O(1) has the least absolute values of both mean absolute deviation (MAD) and root mean square deviation (RMSD), followed by HSAPO-34-O(2). In total, the absolute value of MAD-abc is within 0.08 Å and the RMSD-abc ranges from 0.02 to 0.14. This implies that each Brönsted acid site of HSAPO-34 has already deviated from the space group of R 3 ¯ .
As for the lattice volume, one of the Brönsted acid sites with the proton at O(2) is the most close to the dehydrated HSAPO-34 sample reported by the neutron diffraction experiment [3], followed by the Brönsted acid site with the proton at O(1). Relative to the dehydrated HSAPO-34 sample, the absolute values of the lattice volume deviations are within 2.89 Å3 when the proton is located at O(1) or O(2), while it rises to 27.13 Å3 or 14.43 Å3 when the proton is chemically bound to O(4) or O(3).

2.3. Thermal Stability, Acidity Strength, and Vibrational Frequency of Brönsted Hydroxyl (ν(O-H))

2.3.1. Thermal Stability

The calculations indicated that four possible acid sites of HSAPO-34 had different internal energy (E), as shown in Table 4. Among them, the acid sites of HSAPO-34 with the proton at O(3) are at the bottom. When the proton switches to be chemically bound to the framework O(1), O(4), or O(2) atom, the corresponding Brönsted acid site’s internal energy increases by 1.025, 1.590, or 1.935 kcal/mol, respectively. Table 4 shows that the proton at O(2) has a higher internal energy than the others. This is probably because the proton at O(2) is lying at an intersection between two 8MR windows. The proton at O(3) has the lowest internal energy and can be regarded as the most stable acid site considering internal energy. The proton at O(3) is almost located in the plane of 6MR; as calculated, it is about 2.703, 2.838, 3.568, 2.572, and 2.864 Å away from the other five O atoms of 6MR. The hydrogen bond likely forms between the proton H atom and the other five O atoms of 6MR, respectively. This probably accounts for the proton at O(3) being the most stable acid site considering internal energy.
Generally, the lower the internal energy, the higher the stability is. The relatively stability of the Brönsted acid site of HSAPO-34 might increase in order of the proton binding sites of O(2), O(4), O(1), and O(3) when the temperature influence is not considered. Previous neutron diffraction experiments discovered that protons preferred to reside on the framework O(2) or O(4) atom in the dehydrated SAPO-34 sample [3], while the O(1) atom was coordinated by a hydroxonium ion (H3O+), and yet the hydrogen bond formed between O(2) and a water molecule in the hydrated [2] sample. Based on those neutron diffraction experiments, Jeanvoine et al. concluded that the thermal stability of HSAPO-34 increased in order of the proton binding sites of O(2), O(4), and O(1) [2,3]. In view of the relative internal energy, the stability order of various acid sites of HSAPO-34, predicted by the SCAN + rVV10 method, is almost fully confirmed by the experiment [2,3,28].

2.3.2. Vibrational Frequency of Brönsted Hydroxyl (ν(O-H))

The stretching frequency of Brönsted hydroxyl (ν(O-H)) of the acid zeolite HSAPO-34 has received wide attention [3,26,27,35,36,56,57,58]. Previously, Zubkov et al. [27] and Martins et al. [26] experimentally detected that the peaks at 3625/3631 and 3600 cm−1 were assigned to the fundamental stretching vibrations of two types of bridged hydroxyl groups at O(1) and O(2) in crystalline SAPO-34, respectively. Martins et al. [26] also detected a third distinct peak at 3617 cm−1. Yet, Smith et al. [3,56] used neutron diffraction technology to suggest that the ν(O-H) was 3625, 3601, and 3630 cm−1 for HSAPO-34 with protons at O(1), O(2), and O(4), respectively. On the basis of the fully optimized geometries in Section 2.2, phonon calculations of four acid sites of HSAPO-34 were carried out by using the SCAN + rVV10 method, providing some information on the vibrational frequencies. The ν(O-H) of the acid zeolite HSAPO-34 is calculated to be 3706, 3684, 3625, and 3700 cm−1 when the proton resides on the framework O(1), O(2), O(3), and O(4) atoms, respectively (Table 4). Our calculated ν(O-H) is very close to the experimental values [3,26,56], with an increment of 25–75 cm−1. The ν(O-H) of chabazite, calculated by the PDFT method with the GGA approximation of Perdew and Wang [59], was 3875, 3920, 3820, and 3815 cm−1 for the protons at the O(1), O(2), O(3), and O(4) atoms, respectively [58], which are about 100 cm−1 larger than our calculation results from HSAPO-34. It is common for almost all GGA functional to overestimate stretching vibrations [60]. In line with the neutron experiment [3,56], an average correction factor of about 0.979 might be introduced to our calculated vibrational frequencies of HSAPO-34 from the SCAN + rVV10 method.

2.3.3. Proton Siting and Proton Distribution

Proton siting. Proton siting and proton distribution play crucial roles in the hydrogen bond adsorption and the catalysis reaction [61]. Proton siting would be recognized by examining the calculated crystal structure of the Brönsted acid site [62]. As shown in Table 1 and Figure 3, HSAPO-34 has four types of protons: (1) the proton located at the framework O(1)atom (i.e., H(O1)) is almost at 8MR window, pointing to the positive direction of the crystal axis of b, and the H-O(1) bond is nearly parallel to the crystal plane of [010]; (2) the proton bound to O(2) (i.e., H(O2)) lies at an intersection between two narrow 8MR windows along the [010] and [100] planes, and H(O2) points toward the cavity of the 12-member ring; (3) H(O3) lies at an intersection between the 6MR and 4MR windows, and the H(O3) bond is almost parallel with the [001] plane, and (4) H(O4) is almost at the 8MR windows, and the H(O4) bond is almost parallel with the [100] plane.
Proton distribution. The experiments demonstrate that the proton siting is dependent on the synthesis conditions and is not random [63], but the distribution follows the law of thermodynamics under a given condition [29,64]. Therefore, the possible distribution of four possible Brönsted acid sites of HSAPO-34 could be estimated through Maxwell–Boltzmann statistics [65].
< N i >   = N · g i · e − G i ( T ) k b · T ∑ 1 i g i · e − G i ( T ) k b · T
G T = H − T · S
S = S e + S t + S r + S v
S = ∂ F ∂ T V , N
F = − k b · T · l n q N
Equation (2) shows the distribution of Maxwell–Boltzmann statistics. The natural abundance (<Ni>) of the ith configuration is related to the total number of particles (N), the free energy of the ith configuration (Gi), the temperature (T in Kelvin), and the Planck constant (kb). The Gibbs free energy (G(0)) at 0 K is given directly from the VASP calculation results, which are numerically equal to the enthalpy at 0 K. According to statistical mechanics, the Gibbs free energy (G(T)) at a given temperature is associated with the enthalpy (H), the temperature (T), and the entropy (S) (Equation (3)). Herein, the entropy (S) consists of the electric (Se), the transitional (St), the vibrational (Sv), and the rotational (Sr) contribution (Equation (4)). The entropy stems from the derivative of the Helmholtz free energy (F) divided by the temperature at the constant volume and the constant number of particles (Equation (5)). And the Helmholtz free energy (F) is calculated from the partition function (q) (Equation (6)).
The selected temperature ranges from 273.15 K to 773.15 K with a temperature interval of 50.00 K. The acidic single-silicon substituted SAPO-34 zeolite, HSAPO-34, has four possible Brönsted acid sites. The four calculated Brönsted acid sites are with the space group of P1, and the degeneracy gi in Equation (2) is cancellable to 1. When the numerical values of N, Gi, gi, kb, and T are introduced to the Maxwell–Planck statistics equation (Equation (2)), it initiates the distribution of the Brönsted acid site of zeolite HSAPO-34 and further leads to proton distribution.
Figure 4 illustrates the relationship between the proton distribution of HSAPO-34 and the temperature, as calculated by the periodic density functional theory with the SCAN + rVV10 approach. As shown, the proton distribution at either O(1) or O(4) varies weakly with the temperature change, while the temperature has a strong influence on the proton distribution at both O(2) and O(3). The HSAPO-34-O(3) distribution gradually declines with an increasing temperature, contrary to the other three Brönsted acid sites of HSAPO-34.
Predominate distribution. At 273.15–485.50 K, the predominant proton is assumedly located at the framework O(3) atom of HSAPO-34, and its contribution is about 67.85–35.28%. The proton at O(3) is almost in the 8MR window along the [100] plane. Searching for the largest proton distribution, the O(2) site competes with the O(3) site heavily. Above about 485.50 K, the proton at O(2) exceeds that at O(3) in abundance. At 485.50–773.15 K, the O(2) site of HSAPO-34 becomes the most dominant proton binding site, contributing to 35.28–49.36% of the proton. As shown in Table 1, the proton located at the framework O(2) atom mainly lies in an intersection of two 8MR windows. The proton siting with the dominant abundance at high temperatures probably contributes the experimental phenomenon that most catalysis reactions mainly take place in the specific 8MR window of HSAPO-34.
Secondary distribution. Figure 4 also shows that the O(1) site of HSAPO-34 competes with two other proton binding sites of O(2) and O(3) heavily. Below 348.46 K, the proton at O(1) is less than that at O(3) in abundance, and then the O(1) site becomes the secondary dominant Brönsted acid site. This is in line with the previous theoretical suggestion that the O(1) and O(3) sites are the two preferred binding sites of the proton [35]. When the temperature ranges from 348.46 K to 706.48 K, the proton at O(1) assumedly accounts for the third largest contribution, occupying 21.67–23.33%. Among four possible Brönsted acid sites, HSAPO-34 with a proton at O(4) assumedly has the least proton distribution, which is below 6.1% at 273.15–773.15K.
In sum, various Brönsted acid sites of HSAPO-34 are calculated to be dependent on the temperature in abundance, resulting in proton siting and its distribution change with the temperature.

2.3.4. Acid Strength

The charge distribution on the acid zeolites is of considerable interest. Figure 5 shows the calculated Bader charges [66] on the AlO3-O(H)-SiO3 fragment from four possible Brönsted acid sites of HSAPO-34. The Bader charge analysis is based on the full optimized geometry. All oxygen atoms around the incorporated Si atom appear to have a negative charge (about −1.580 e~−1.602 e), which is slightly lower than those (−1.530 e~−1.566 e) bound to the Al atom. The Bader charge on either the Si or Al atom is less than its respective valence electron; the former ranges from 3.148 e to 3.159 e, while the latter ranges from 2.478 e to 2.480 e. The proton charges from the four Brönsted acid sites are slightly different, calculated to be 0.677 e, 0.655 e, 0.659 e, and 0.660 e when the proton is at O(1), O(2), O(3), and O(4), respectively (Table 4 and Figure 5). In view of proton charge, the acid strength would gradually reduce in order of the proton binding site of O(1), O(4), O(3), and O(2).
DPE = E(H+) + E(Z−) − E(HZ)
The deprotonated energy (DPE) is capable of reflecting the acid strength of the solid acid catalyst. The deprotonated internal energy is defined as the internal energy difference between the deprotonated (Z−) plus the bare proton (H+) and the protonated (HZ) zeolites (Equation (7)). It is worth noting that the deprotonated zeolite has a negative charge. For the deprotonated HSAPO-34, its initial structure used for the E-V scan comes from any possible acid sites of HSAPO-34 by eliminating the Brönsted acid proton, and its precise structure would be obtained through a four-step procedure similar to that in Section 2.2. In comparison with the deprotonated HSAPO-34, the Si-O(H) bond length in each acid site of HSAPO-34 is elongated by about 0.134–0.138 Å in the presence of the proton. According to Equation (7), the calculated DPE value is 1526, 1522, 1530, and 1524 kJ/mol when HSAPO-34 is with the proton at O(1), O(2), O(3), and O(4), respectively (Table 4). In general, the less DPE, the stronger the acid strength is. The acid strength assumedly reduces in order of the proton binding site of O(2), O(4), O(1), and O(3). The DPE differences between one other are within 8 kJ/mol, and then those four acid sites have no vast difference in acid strength. Our calculated DPE value for each Brönsted acid site of HSAPO-34 is 20–29 kJ/mol lower than that for the zeolite chabazite previously reported by Jone and Iglesia [67]. It is worth noting that the DPE values are also dependent on the Si island [4].
< D P E > = E Z − + E H + − < E Z H >
Recently, Jone and Iglesia [67] put forward the ensemble-averaged DPE in view of the statistics (Equation (8)). Herein, <DPE> is the ensemble-averaged DPE value in the thermal equilibration state, < E Z H > is the ensemble-averaged energy of four possible Bronsted acid sites at each atom location of the acid source. It could be estimated by the statistical average over the internal energy and the free energy of each Brönsted acid site (Equation (9)). The item of free energy is related to the probability of the proton.
< E Z H >   = ∑ i = 1 4 E Z H , i · e − G Z H , i ( T ) k b · T ∑ i = 1 4 e − G Z H , i ( T ) k b · T
Figure 6 illustrates the calculated ensemble-averaged DPE (<DPE>) of HSAPO-34 dependent on the temperature. As we predicted, the <DPE> is about 1525 kJ/mol at a suitable methanol-to-olefins catalysis temperature of 673.15 K [6]. The <DPE> value gradually reduces as the temperature increases. The change magnitude of <DPE> is within 3.49 kJ/mol when the temperature ranges from 273.15 K to 773.15 K.

3. Materials and Methods

All calculations are performed by the Vienna Ab initio Simulation Package (VASP) [68,69,70,71] using the plane-wave periodic gradient-corrected density functional theory (DFT) methods. The strongly constrained and appropriately normed (SCAN) meta-generalized-gradient approximation (meta-GGA) [48,49,50,51] is used with PAW [72] Perdew–Burke–Ernzerhof (PBE) [38] pseudopotentials (version 52) and the rVV10 nonlocal correlation functional of Peng et al. [51]. This theoretical calculation approach is also called the SCAN + rVV10 method [51], implemented in VASP packages.
The initial structure of HSAPO-34 is constructed through the isomorphous substitution of silicon atom in the all-silica CHA structural topology by following the SMIII mechanism [73]. The CHA structure is obtained from the Database of Zeolite Structures in Internal Zeolite Association (IZA) [53]. The acidic single-silicon substituted SAPO-34 zeolites (HSAPO-34) with rhombohedral representation are selected as the calculation model. The model zeolite HSAPO-34 is the single-unit cell containing six Al atoms, five P atoms, one Si atom, twenty-four O atoms, and one H atom.
The lattice parameter is calculated from the crystal structure of zeolite in the equilibrium state [55]. Prior to the full optimization, scanning of the E-V curve is carried out by partially optimizing a series of crystal structures for each Brönsted acid site at the constant volume. The volume is selected from 750 Å3 to 900 Å3 with an interval of 10 Å3. Fitting the E-V curve with the Birch–Murnaghan equation would produce a rough crystal structure, which is used as the starting crystal structure for full relaxation of the cell shape, the lattice volume, and the ion position. The full optimized structure is further verified as the desirable structure of zeolite HSAPO-34 by phonon analysis. The energy cutoff is set to 570 eV for the E-V scan and 900 eV for the full optimization. The convergence threshold selected is 1 × 10−2 eV/Å for the atomic force and 1 × 10−7 eV for the energy. The K-point meshes are selected as 2 × 2 × 2 to sample the Brillouin zone for HSAPO-34. The density functional perturbation theory (DFPT) was used to calculate the second derivatives (Hessian matrix and phonon frequencies) without the symmetry restriction. Three small imaginary frequencies (below 6i cm−1) of each acid site of HSAPO-34 are present yet not easily removed, which would have an influence on the relative distribution of the proton. An analysis of each imaginary vibration mode indicates that the associated atoms oscillate almost along the same direction. To some extent, the associated nuclear motions are regarded as the translations/rotations of the species. The entropy contribution from 12 cm−1 harmonic vibration approximates is almost equal to that from a pseudotranslation/rotational degree of freedom [74,75,76]. In this study, the entropy contribution from those small imaginary frequencies was calculated after they were replaced by 12 cm−1 to reduce as much inaccuracy derived from the imaginary frequencies as possible.
The electronic charges on the atoms of each acid site of HSAPO-34 were analyzed by the Bader charge method developed by the Henkelman group [66] on the basis of the full optimized geometry of HSAPO-34.

4. Conclusions

Through a four-step procedure, the fine crystal structure of HSAPO-34 has been accurately simulated by the periodic density functional theory with the SCAN + rVV10 method. Four unequivalent framework O atoms around the incorporated Si atom of HSAPO-34 would produce four possible Brönsted acid sites, deviating from the perfect space group of R 3 ¯ . Under this high-level approximation (i.e., SCAN-meta-GGA), the calculated structural properties on the Brönsted acid sites of HSAPO-34 match the experiment well [2,3,28], including the lattice parameter, the internal ion position, the acid strength, the proton siting, and the relative thermal stability. It seems that the meta-GGA method could compensate for the some drawbacks of the GGA method [36] and then give rise to the satisfying calculation results of the acid zeolite HSAPO-34. The statistics analysis further indicates that both the proton distribution and the ensemble-averaged DPE of HSAPO-34 are distinctly altered by the temperature. It is worth noting that the numerical data mentioned in this study may fluctuate more or less with the selected method, especially on the natural abundance and the temperature threshold.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28248087/s1.

Author Contributions

X.C.: Conceptualization, methodology, software, validation, formal analysis, investigation, resources, data curation, writing—original draft preparation, writing—review and editing, visualization, supervision, project administration, funding acquisition. T.Y.: the IR spectra. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Basic Research Program (No. 2019GN018) of Shandong University.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Acknowledgments

We express great thanks to the support of Shandong University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wilson, S.T.; Lok, B.M.; Messina, C.A.; Cannan, T.R.; Flanigen, E.M. Alumniophosphate Molecular Sieves: A New Class of Microporous Crystalline Inorganic Solids. J. Am. Soc. Chem. 1982, 104, 1146–1147. [Google Scholar] [CrossRef]
  2. Smith, L.; Cheetham, A.K.; Morris, R.E.; Marchese, L.; Thomas, J.M.; Wright, P.A.; Chen, J. On the Nature of Water Bound to a Solide Acid Catalyst. Science 1996, 271, 799–802. [Google Scholar] [CrossRef]
  3. Smith, L.; Cheetham, A.K.; Marchese, L.; Thomas, J.M.; Wright, P.A.; Chen, J.; Gianotti, E. A Quantitative Description of the Active Sites in the Dehydrated Acid Catalyst HSAPO-34 for the Conversion of Methanol to Olefins. Catal. Lett. 1996, 41, 13–16. [Google Scholar] [CrossRef]
  4. Suzuki, K.; Nishio, T.; Katada, N.; Sastre, G.; Niwa, M. Ammonia IRMS-TPD measurements on Brønsted acidity of proton-formed SAPO-34. Phys. Chem. Chem. Phys. 2011, 13, 3311–3318. [Google Scholar] [CrossRef]
  5. Potter, M.E. Down the Microporous Rabbit Hole of Silicoaluminophosphates: Recent Developments on Synthesis, Characterization, and Catalytic Applications. ACS Catal. 2020, 10, 9758–9789. [Google Scholar] [CrossRef]
  6. Tian, P.; Wei, Y.; Ye, M.; Liu, Z. Methanol to Olefins (MTO): From Fundamentals to Commercialization. ACS Catal. 2015, 5, 1922–1938. [Google Scholar] [CrossRef]
  7. Rostami, R.B.; Lemraski, A.S.; Ghavipour, M.; Behbahani, R.M.; Shahraki, B.H.; Hamule, T. Kinetic modelling of methanol conversion to light olefins process over silicoaluminophosphate (SAPO-34) catalyst. Chem. Eng. Res. Des. 2016, 106, 347–355. [Google Scholar] [CrossRef]
  8. Dai, W.; Cao, G.; Yang, L.; Wu, G.; Dyballa, M.; Hunger, M.; Guan, N.; Li, L. Insights into the catalytic cycle and activity of methanol-to-olefin conversion over low-silica AlPO-34 zeolites with controllable Brønsted acid density. Catal. Sci. Technol. 2017, 7, 607–618. [Google Scholar] [CrossRef]
  9. Li, Z.; Wang, J.; Qu, Y.; Liu, H.; Tang, C.; Miao, S.; Feng, Z.; An, H.; Li, C. Highly Selective Conversion of Carbon Dioxide to Lower Olefins. ACS Catal. 2017, 7, 8544–8548. [Google Scholar] [CrossRef]
  10. Liu, Y.; Deng, D.; Bao, X. Catalysis for Selected C1 Chemistry. Chem 2020, 6, 2497–2514. [Google Scholar] [CrossRef]
  11. Han, J.; Yang, G.; Ding, H.; Chen, X. Revealing inherent factors of SAPO-34 zeolites etching towards the fabrication of hierarchical structure. Microporous Mesoporous Mater. 2021, 319, 111067. [Google Scholar] [CrossRef]
  12. Vallace, A.; Kester, G.; Casteel, W.; Lau, G.; Whitley, R.; Coe, C. A Study of Structural Defects in X- and Y-Type Zeolites and Their Effect on Their Transformation to Aluminum-Rich Chabazite. J. Phys. Chem. C 2021, 125, 12848–12856. [Google Scholar] [CrossRef]
  13. Usman, M.; Ghanem, A.S.; Shah, S.N.A.; Garba, M.D.; Khan, M.Y.; Khan, S.; Humayun, M.; Khan, A.L. A Review on SAPO-34 Zeolite Materials for CO2 Capture and Conversion. Chem. Rec. 2022, 22, e202200039. [Google Scholar] [CrossRef] [PubMed]
  14. Kheimi, M.; Salamah, S.K. Simulation of temperature swing adsorption process to purify hydrogen for fuel cell uses by SAPO34 as adsorbent. Chemosphere 2023, 338, 139454. [Google Scholar] [CrossRef] [PubMed]
  15. Yang, G.; Ran, J.; Du, X.; Wang, X.; Ran, Z.; Chen, Y.; Zhang, L.; Crittenden, J. Understanding the nature of NH3-coordinated active sites and the complete reaction schemes for NH3-SCR using Cu-SAPO-34 catalysts. Phys. Chem. Chem. Phys. 2021, 23, 4700. [Google Scholar] [CrossRef] [PubMed]
  16. Uzunova, E.L. Theoretical study of nitrogen dioxide and nitric oxide co-adsorption and DeNOx reaction on Cu-SAPO−34 and Cu-SSZ-13 in presence of Brønsted acid sites. Mol. Catal. 2018, 447, 47–55. [Google Scholar] [CrossRef]
  17. Wilson, S.; Barger, P. The characteristics of SAPO-34 which influence the conversion of methanol to light olefins. Microporous Mesoporous Mater. 1999, 29, 117–126. [Google Scholar] [CrossRef]
  18. Aguayo, A.T.; Gayubo, A.G.; Vivanco, R.; Olazar, M.; Bilbao, J. Role of acidity and microporous structure in alternative catalysts for the transformation of methanol into olefins. Appl. Catal. A Gen. 2005, 283, 197–207. [Google Scholar] [CrossRef]
  19. Wang, L.; Li, W.; Schmieg, S.J.; Weng, D. Role of Brønsted acidity in NH3 selective catalytic reduction reaction on Cu/SAPO-34 catalysts. J. Catal. 2015, 324, 98–106. [Google Scholar] [CrossRef]
  20. Li, X.; Shen, W.; Zheng, A. The influence of acid strength and pore size effect on propene elimination reaction over zeolites: A theoretical study. Microporous Mesoporous Mater. 2019, 278, 121–129. [Google Scholar] [CrossRef]
  21. Sierraalta, A.; Añez, R.; Coll, D.S.; Alejos, P. Conversion of methanol to dimethyl ether over silicoaluminophosphates: Isolated acid sites and the influence of silicon islands. A DFT-ONIOM study. Microporous Mesoporous Mater. 2020, 292, 109732. [Google Scholar] [CrossRef]
  22. Flanigen, E.M.; Lok, B.M.; Patton, R.L.; Wilson, S.T. Aluminophosphate molecular sieves and the periodic table. Pure Appl. Chem. 1986, 58, 1351–1358. [Google Scholar] [CrossRef]
  23. Shen, W.; Li, X.; Wei, Y.; Tian, P.; Deng, F.; Han, X.; Bao, X. A study of the acidity of SAPO-34 by solid-state NMR spectroscopy. Microporous Mesoporous Mater. 2012, 158, 19–25. [Google Scholar] [CrossRef]
  24. Sastre, G.; Lewis, D.W.; Catlow, C.R.A. Modeling of Silicon Substitution in SAPO-5 and SAPO-34 Molecular Sieves. J. Phys. Chem. B 1997, 101, 5249–5262. [Google Scholar] [CrossRef]
  25. Dent, L.S.; Smith, J.V. Crystal Structure of Chabazite, A Molecular Sieve. Nature 1958, 181, 1794–1796. [Google Scholar] [CrossRef]
  26. Martins, G.A.V.; Berlier, G.; Coluccia, S.; Pastore, H.O.; Superti, G.B.; Gatti, G.; Marchese, L. Revisiting the Nature of the Acidity in Chabazite-Related Silicoaluminophosphates: Combined FTIR and 29Si MAS NMR Study. J. Phys. Chem. C 2007, 111, 330–339. [Google Scholar] [CrossRef]
  27. ZubkoY, S.A.; Kustov, L.M.; Kazansky, V.B.; Girnus, I.; Fricke, R. Investigation of Hydroxyl Groups in Crystalline Silicoaluminophosphate SAPO-34 by Diffuse Reflectance Infrared Spectroscopy. J. Chem. Soc. Faraday Trans. 1991, 87, 897–900. [Google Scholar] [CrossRef]
  28. Pfeifer, H.; Freude, D.; Hunger, M. Nuclear magnetic resonance studies on the acidity of zeolites and related catalysts. Zeolites 1985, 5, 274–286. [Google Scholar] [CrossRef]
  29. Iwase, Y.; Motokura, K.; Koyama, T.-R.; Miyaji, A.; Baba, T. Influence of Si Distribution in Framework of SAPO-34 and its Particle Size on Propylene Selectivity and Production Rate for Conversion of Ethylene to Propylene. Phys. Chem. Chem. Phys. 2009, 11, 9268–9277. [Google Scholar] [CrossRef]
  30. Sena, F.C.; Souza, B.F.D.; Almeida, N.C.D.; Cardoso, J.S.; Fernandes, L.D. Influence of framework composition over SAPO-34 and MeAPSO-34 acidity. Appl. Catal. A Gen. 2011, 406, 59–62. [Google Scholar] [CrossRef]
  31. Zhang, X.; Yang, M.; Wang, L.; Han, J.; Lou, C.; Xu, S.; Zhang, Y.; Wu, R.A.; Tian, P.; Liu, Z. Recognizing the Minimum Structural Units Driving the Crystallization of SAPO-34 in a Top-Down Process. Chem. Eur. J. 2023, 29, e202203886. [Google Scholar] [CrossRef]
  32. Liu, Z.; Hua, Y.; Wang, J.; Dong, X.; Tian, Q.; Han, Y. Recent Progress in the Direct Synthesis of Hierarchical Zeolites: Synthetic Strategies and Characterization Methods. Mater. Chem. Front. 2017, 1, 2195–2212. [Google Scholar] [CrossRef]
  33. Izadbakhsh, A.; Farhadi, F.; Khorasheh, F.; Sahebdelfar, S.; Asadi, M.; Feng, Y.Z. Effect of SAPO-34′s composition on its physico-chemical properties and deactivation in MTO process. Appl. Catal. A Gen. 2009, 364, 48–56. [Google Scholar] [CrossRef]
  34. Skylaris, C.-K. A Benchmark for Materials Simulation-Material Properties Can Now Be Predicted Reliably from First-Principles Calculations. Science 2016, 351, 1394–1395. [Google Scholar] [CrossRef] [PubMed]
  35. Gale, J.D.; Shah, R.; Payne, M.C.; Stich, I.; Terakura, K. Methanol in microporous materials from first principles. Catal. Today 1999, 50, 525–532. [Google Scholar] [CrossRef]
  36. Jeanvoine, Y.; Ángyán, J.G.; Kresse, G.; Hafner, J. Brønsted Acid Sites in HSAPO-34 and Chabazite: An Ab Initio Structural Study. J. Phys. Chem. B 1998, 102, 5573–5580. [Google Scholar] [CrossRef]
  37. Perdew, J.P.; Wang, Y. Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B 1992, 45, 13244–13249. [Google Scholar] [CrossRef]
  38. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  39. Mattsson, A.E.; Armiento, R.; Schultz, P.A.; Mattsson, T.R. Nonequivalence of the generalized gradient approximations PBE and PW91. Phys. Rev. B 2006, 73, 195123. [Google Scholar] [CrossRef]
  40. Kohn, W. Nobel Lecture: Electronic Structure of Matter-Wave Functions and Density Functionals. Rev. Mod. Phys. 1999, 71, 1253–1266. [Google Scholar] [CrossRef]
  41. Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  42. Hellmann, H. A New Approximation Method in the Problem of Many Electrons. J. Chem. Phys. 1935, 3, 61. [Google Scholar] [CrossRef]
  43. Kohn, W.; Sham, L.J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133–A1138. [Google Scholar] [CrossRef]
  44. Lin-Chung, P.J. Limitation of the Pseudopotential Method. Phys. Rev. B 1973, 8, 4043–4045. [Google Scholar] [CrossRef]
  45. Ziesche, P.; Kurth, S.; Perdew, J.P. Density Functionals from LDA to GGA. Comput. Mater. Sci. 1998, 11, 122–127. [Google Scholar] [CrossRef]
  46. Perdew, J.P.; Schmidt, K. Jacob’s Ladder of Density Functional Approximations for the Exchange-Correlation Energy. AIP Conf. Proc. 2001, 577, 1–20. [Google Scholar]
  47. Mackrodt, W.C. A Note on An Aspect of Pseudopotential Theory. Theor. Chim. Acta 1973, 30, 119–126. [Google Scholar] [CrossRef]
  48. Sun, J.; Ruzsinszky, A.; Perdew, J.P. Strongly Constrained and Appropriately Normed Semilocal Density Functional. Phys. Rev. Lett. 2015, 115, 036402. [Google Scholar] [CrossRef]
  49. Sun, J.; Remsing, R.C.; Zhang, Y.; Sun, Z.; Ruzsinszky, A.; Peng, H.; Yang, Z.; Paul, A.; Waghmare, U.; Wu, X.; et al. Accurate first-principles structures and energies of diversely bonded systems from an efficient density functional. Nat. Chem. 2016, 8, 831–836. [Google Scholar] [CrossRef]
  50. Zhang, Y.; Kitchaev, D.A.; Yang, J.; Chen, T.; Dacek, S.T.; Samiento-Pérez, R.A.; Marques, M.A.L.; Peng, H.; Ceder, G.; Perdew, J.P.; et al. Efficient First-Principles Prediction of Solid Stability: Towards Chemical Accuracy. NPJ Comput. Mater. 2018, 4, 9. [Google Scholar] [CrossRef]
  51. Peng, H.; Yang, Z.-H.; Perdew, J.P.; Sun, J. Versatile van der Waals Density Functional Based on a Meta-Generalized Gradient Approximation. Phys. Rev. X 2016, 6, 041005. [Google Scholar] [CrossRef]
  52. Chen, X.-F. Periodic Density Functional Theory (PDFT) Simulating Crystal Structures with Microporous CHA Framework: An Accuracy and Efficiency Study. Inorganics 2023, 11, 215. [Google Scholar] [CrossRef]
  53. Baerlocher, C.; McCusker, L.B. Database of Zeolite Structures. Available online: http://www.iza-structure.org/databases (accessed on 5 November 2023).
  54. Ángyán, J.G.; Parsons, D.; Jeanvoine, Y. Ab initio simulations of zeolite reactivity. In Theoretical Aspects of Heterogeneous Catalysis; Nascimento, M.A.C., Ed.; Springer: Dordrecht, The Netherlands, 2001; pp. 77–108. Available online: https://link.springer.com/chapter/10.1007/0-306-47667-3_4 (accessed on 5 November 2023).
  55. Birch, F. Finite Elastic Strain of Cubic Crystals. Phys. Rev. 1947, 71, 809–824. [Google Scholar] [CrossRef]
  56. Smith, L.J.; Davidson, A.; Cheetham, A.K. Aneutron Diffraction and Infrared Spectroscopy Study of the Acid Form of the Aluminosilicate Zeolite, Chabazite (H-SSZ-13). Catal. Lett. 1997, 49, 143–146. [Google Scholar] [CrossRef]
  57. Marchese, L.; Chen, J.; Wright, P.A.; Thomas, J.M.J. Formation of hydronium at the Broensted site in SAPO-34 catalysts. J. Phys. Chem. 1993, 97, 8109–8112. [Google Scholar] [CrossRef]
  58. Shah, R.; Gale, J.D.; Payne, M.C. Methanol Adsorption in Zeolites: A First-Principles Study. J. Phys. Chem. 1996, 100, 11688–11697. [Google Scholar] [CrossRef]
  59. Perdew, J.P.; Chevary, J.A.; Vosko, S.H.; Jackson, K.A.; Pederson, M.R.; Singh, D.J.; Fiolhais, C. Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation. Phys. Rev. B 1992, 46, 6671–6687. [Google Scholar] [CrossRef] [PubMed]
  60. Farnworth, K.J.; O’Malley, P.J. A density functional study of acidic hydroxyl groups in zeolites and their interaction with carbon monoxide. J. Phys. Chem. 1996, 100, 1814–1819. [Google Scholar] [CrossRef]
  61. Kilburn, L.; DeLuca, M.; Hoffman, A.J.; Patel, S.; Hibbitts, D. Comparing alkene-mediated and formaldehyde-mediated diene formation routes in methanol-to-olefins catalysis in MFI and CHA. J. Catal. 2021, 400, 124–139. [Google Scholar] [CrossRef]
  62. Feng, P.; Chen, X.-F.; Li, X.-J.; Zhao, D.; Xie, S.-J.; Xu, L.-Y.; He, G.-Z. The distribution analysis on the proton siting and the acid strength of the zeolite ferrierite: A computational study. Microporous Mesoporous Mater. 2017, 239, 354–362. [Google Scholar] [CrossRef]
  63. Simperler, A.; Bell, R.G.; Foster, M.D.; Gray, A.E.; Lewis, D.W.; Anderson, M.W. Probing the acid strength of Brønsted acidic zeolites with acetonitrile: An atomistic and quantum chemical study. J. Phys. Chem. B 2004, 108, 7152–7161. [Google Scholar] [CrossRef]
  64. Bordiga, S.; Palomino, G.T.; Pazè, C.; Zecchina, A. Vibrational spectroscopy of H2, N2, CO and NO adsorbed on H, Li, Na, K-exchanged ferrierite. Microporous Mesoporous Mater. 2000, 34, 67–80. [Google Scholar] [CrossRef]
  65. Maxwell, J.C. Illustrations of the Dynamical Theory of Gases. Available online: https://www.worldscientific.com/doi/epdf/10.1142/9781848161337_0011 (accessed on 5 November 2023).
  66. Henkelman, G.; Arnaldson, A.; Jopsson, H. A fast and robust algorithm for Bader decomposition of charge density. Comput. Mater. Sci. 2006, 31, 354–360. [Google Scholar] [CrossRef]
  67. Jones, A.J.; Iglesia, E. The Strength of Brønsted Acid Sites in Microporous Aluminosilicates. ACS Catal. 2015, 5, 5741–5755. [Google Scholar] [CrossRef]
  68. Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 48, 13115–13118. [Google Scholar] [CrossRef]
  69. Kresse, G.; Hafner, J. Ab Initio Molecular-Dynamics Simulation of the Liquid-Metal Amorphous-Semiconductor Transition in Germanium. Phys. Rev. B Condens. Matter. Mater. Phys. 1994, 49, 14251–14269. [Google Scholar] [CrossRef] [PubMed]
  70. Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mat. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  71. Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-energy Calculations Using a Plane—Wave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef]
  72. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef]
  73. Flanigen, E.M. SMIII mechanism. Elsevier Sci. B 1986, 85, 653–685. [Google Scholar]
  74. Ghorbanpour, A.; Rimer, J.D.; Grabow, L.C. Computational Assessment of the Dominant Factors Governing the Mechanism of Methanol Dehydration over H-ZSM-5 with Heterogeneous Aluminum Distribution. ACS Catal. 2016, 6, 2287–2298. [Google Scholar] [CrossRef]
  75. Brogaard, R.Y.; Weckhuysen, B.M.; Nørskov, J.K. Guest–host interactions of arenes in H-ZSM-5 and their impact on methanol-to-hydrocarbons deactivation processes. J. Catal. 2013, 300, 235–241. [Google Scholar] [CrossRef]
  76. Brogaard, R.Y.; Henry, R.; Schuurman, Y.; Medford, A.J.; Moses, P.G.; Beato, P.; Svelle, S.; Nørskov, J.K.; Olsbye, U. Methanol-to-hydrocarbons conversion: The alkene methylation pathway. J. Catal. 2014, 314, 159–169. [Google Scholar] [CrossRef]
Molecules 28 08087 g001
Figure 1. Single silicon-doped alumino-phosphates (HSAPO-34): (a) the single-unit cell (1UC) with rhombohedral representation, (b) the expanded lattice structure. Note that the proton is omitted for simplicity.
Figure 1. Single silicon-doped alumino-phosphates (HSAPO-34): (a) the single-unit cell (1UC) with rhombohedral representation, (b) the expanded lattice structure. Note that the proton is omitted for simplicity.
Molecules 28 08087 g001
Molecules 28 08087 g002
Figure 2. The total energy depends on the lattice volume. The proton binding sites of HSAPO-34 are (a) O(1), (b) O(2), (c) O(3), and (d) O(4).
Figure 2. The total energy depends on the lattice volume. The proton binding sites of HSAPO-34 are (a) O(1), (b) O(2), (c) O(3), and (d) O(4).
Molecules 28 08087 g002
Molecules 28 08087 g003
Figure 3. Schematic diagram for four possible Brönsted acid sites of HSAPO-34. The proton binding sites are (a) O1, (b) O2, (c) O3, and (d) O4. Noting that Al in cycn, Si in blue, P in purple, O in red, and the number in figure is the numbering of O atom.
Figure 3. Schematic diagram for four possible Brönsted acid sites of HSAPO-34. The proton binding sites are (a) O1, (b) O2, (c) O3, and (d) O4. Noting that Al in cycn, Si in blue, P in purple, O in red, and the number in figure is the numbering of O atom.
Molecules 28 08087 g003
Molecules 28 08087 g004
Figure 4. The proton distribution of HSAPO-34 depends on the temperature.
Figure 4. The proton distribution of HSAPO-34 depends on the temperature.
Molecules 28 08087 g004
Molecules 28 08087 g005
Figure 5. Bader charges on AlO3-O(H)-SiO3 fragment from four possible Brönsted acid sites of HSAPO-34. The proton binding sites are (a) O(1), (b) O(2), (c) O(3), and (d) O(4). Noting that Al in cyan, Si in blue, and O in red.
Figure 5. Bader charges on AlO3-O(H)-SiO3 fragment from four possible Brönsted acid sites of HSAPO-34. The proton binding sites are (a) O(1), (b) O(2), (c) O(3), and (d) O(4). Noting that Al in cyan, Si in blue, and O in red.
Molecules 28 08087 g005
Molecules 28 08087 g006
Figure 6. Calculated ensemble-averaged DPE (<DPE>) of HSAPO-34 depends on the temperature.
Figure 6. Calculated ensemble-averaged DPE (<DPE>) of HSAPO-34 depends on the temperature.
Molecules 28 08087 g006
Table 1. The atomic numbering relationship of AlPO4-34 framework between this work (same as refs. [3,36]) and the database of zeolite structures in the internal zeolite associates (IZA), as well as the location of O atom.
Table 1. The atomic numbering relationship of AlPO4-34 framework between this work (same as refs. [3,36]) and the database of zeolite structures in the internal zeolite associates (IZA), as well as the location of O atom.
This workO(1)O(2)O(3)O(4)
IZAO(4)O(3)O(1)O(2)
Location4MR(2) α
8MR(1,+) α,β		4MR(1) α
6MR(1) α,β
8MR(1,+)		4MR(2) α
6MR(1,+) α,β		4MR(1) α
8MR(2,+) α,β
α The number in parentheses is the number of window ring; β The symbol “+” in parentheses is the rough location that the Brönsted acid proton points toward.
Table 2. The fitted parametersof the Birch–Murnaghan equation.
Table 2. The fitted parametersof the Birch–Murnaghan equation.
ParametersBrönsted Acid Sites of HSAPO-34
O(1)O(2)O(3)O(4)
V0/Å3821.82810.98820.13806.62
E0/eV−229.33−229.33−229.34−229.31
B00.130.080.120.01
B 0 ′ −27.37−43.85−30.35−400.24
Table 3. The calculated lattice parameters of HSAPO-34 with the rhombohedral presentation.
Table 3. The calculated lattice parameters of HSAPO-34 with the rhombohedral presentation.
ParametersBrönsted Acid Sites of HSAPO-34
O(1)O(2)O(3)O(4)Exp [3]
a (Å)9.409.429.409.199.40
b (Å)9.499.399.529.389.40
c (Å)9.269.429.329.399.40
α (°)94.1294.7993.5895.7794.27
β (°)94.1394.3694.1696.1294.27
γ (°)94.0394.6795.5896.0494.27
V0 (Å3)819.50824.46807.96795.26822.39
MAD-abc−0.020.01−0.05−0.08-
RMSD-abc0.100.020.140.12-
MAD-αβγ−0.180.34−0.701.71-
RMSD-αβγ0.180.380.851.71-
Table 4. The relative internal energy (ΔE), the deprotonated internal energy (DPE), and the stretching frequency of hydroxyl (ν(O-H)) for the acid zeolite HSAPO-34.
Table 4. The relative internal energy (ΔE), the deprotonated internal energy (DPE), and the stretching frequency of hydroxyl (ν(O-H)) for the acid zeolite HSAPO-34.
ParametersBrönsted Acid Sites of HSAPO-34
O(1)O(2)O(3)O(4)
ΔE/(kcal/mol)1.0251.93501.59
DPE/(kJ/mol)1526152215301524
proton charge (e)0.6770.6550.6690.66
R(O-H)(Å)0.9770.9790.9810.978
R(Si-O(H))1.7661.7541.7741.751
ν(O-H)/cm−13706368436253700
ν(O-H)/cm−1 [3,56]36253601-3630
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, X.; Yu, T. Simulating Crystal Structure, Acidity, Proton Distribution, and IR Spectra of Acid Zeolite HSAPO-34: A High Accuracy Study. Molecules 2023, 28, 8087. https://doi.org/10.3390/molecules28248087

AMA Style

Chen X, Yu T. Simulating Crystal Structure, Acidity, Proton Distribution, and IR Spectra of Acid Zeolite HSAPO-34: A High Accuracy Study. Molecules. 2023; 28(24):8087. https://doi.org/10.3390/molecules28248087

Chicago/Turabian Style

Chen, Xiaofang, and Tie Yu. 2023. "Simulating Crystal Structure, Acidity, Proton Distribution, and IR Spectra of Acid Zeolite HSAPO-34: A High Accuracy Study" Molecules 28, no. 24: 8087. https://doi.org/10.3390/molecules28248087

APA Style

Chen, X., & Yu, T. (2023). Simulating Crystal Structure, Acidity, Proton Distribution, and IR Spectra of Acid Zeolite HSAPO-34: A High Accuracy Study. Molecules, 28(24), 8087. https://doi.org/10.3390/molecules28248087

Article Metrics

Back to TopTop