# Adsorption of toluene and water over cationic-exchanged Y zeolites: A DFT exploration

Etienne P Hessou, Lucie A Bédé, Hicham Jabraoui, Abderrahmane Semmeq,

Michael Badawi, Valentin Valtchev

► To cite this version:

Etienne P Hessou, Lucie A Bédé, Hicham Jabraoui, Abderrahmane Semmeq, Michael Badawi, et al. Adsorption of toluene and water over cationic-exchanged Y zeolites: A DFT exploration. Molecules, 2021, 26 (18), pp.5486. 10.3390/molecules26185486. cea-03343141

HAL Id: cea-03343141

https://cea.hal.science/cea-03343141

Submitted on 14 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Adsorption of Toluene and Water over Cationic-Exchanged Y Zeolites: A DFT Exploration

Etienne P. Hessou <fnref n="1" />,* and Lucie A. Bédé <fnref n="2" />, Hicham Jabraoui <fnref n="3" />, Abderrahmane Semmeq <fnref n="1" />, Michael Badawi <fnref n="1" />
and Valentin Valtchev <fnref n="4" />

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Hessou, E.P.; Bédé, L.A.; Jabraoui, H.; Semmeq, A.; Badawi, M.; Valtchev, V. Adsorption of Toluene and Water over Cationic-Exchanged Y Zeolites: A DFT Exploration.
*Molecules* **2021**, *26*, 5486. https://doi.org/10.3390/molecules26185486

Academic Editors: Hristiyan A. Aleksandrov, Margarita D. Popova and T. Jean Daou

Received: 21 July 2021
Accepted: 2 September 2021
Published: 9 September 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** In this study, density functional theory (DFT) calculations have been performed to investigate the adsorption mechanisms of toluene and water onto various cationic forms of Y zeolite (LiY, NaY, KY, CsY, CuY and AgY). Our computational investigation revealed that toluene is mainly adsorbed via π-interactions on alkalis exchanged Y zeolites, where the adsorbed toluene moiety interacts with a single cation for all cases with the exception of CsY, where two cations can simultaneously contribute to the adsorption of the toluene, hence leading to the highest interaction observed among the series. Furthermore, we find that the interaction energies of toluene increase while moving down in the alkaline series where interaction energies are 87.8, 105.5, 97.8, and 114.4 kJ/mol for LiY, NaY, KY and CsY, respectively. For zeolites based on transition metals (CuY and AgY), our calculations reveal a different adsorption mode where only one cation interacts with toluene through two carbon atoms of the aromatic ring with interaction energies of 147.0 and 131.5 kJ/mol for CuY and AgY, respectively. More importantly, we show that water presents no inhibitory effect on the adsorption of toluene, where interaction energies of this latter were 10 kJ/mol (LiY) to 47 kJ/mol (CsY) higher than those of water. Our results point out that LiY would be less efficient for the toluene/water separation while CuY, AgY and CsY would be the ideal candidates for this application.

**Keywords:** adsorption; toluene; DFT; zeolite Y

## 1. Introduction

Volatile organic compounds (VOCs) are major air pollutants that are harmful to human health and to the environment [1–4]. These pollutants affect air quality and have an impact on global warming [5–8]. Among them, toluene is one of the worst pollutants classified by the Environmental Protection Agency of the United States due to its toxicity and significant emissions [9]. Toluene is used in many industrial products such as paints, thinners, glues, adhesives, cleaning products, and so forth [10]. Exposure or inhalation of toluene vapors can cause health issues such as neurotoxicity, cardiopulmonary problems, atrophy, encephalopathy, and eye irritation [11–18]. The severity of the symptoms is often correlated with the concentration of toluene vapor. It has been shown that exposure to high concentrations of toluene (>45,000 mg/m<fnref n="3" />) for about 12 h can be fatal [19,20]. Besides its toxic activity, toluene is also considered a biomarker [21–23]; abnormally high concentrations of toluene can be found in the breath of patients with lung cancer. Therefore, toluene sensors can be very helpful for the early-stage diagnosis of lung cancer.

Given the impact of toluene on the environment and on human health, several experimental and theoretical studies have focused on reducing the level of this pollutant using

<footnote n="1">¹ Laboratoire de Physique et Chimie Théoriques, Faculté des Sciences et Technologies, CNRS, Université de Lorraine, Boulevard des Aiguillettes, 54500 Vandoeuvre-lès-Nancy, France; abderrahmane.semmeq@univ-lorraine.fr (A.S.); michael.badawi@univ-lorraine.fr (M.B.)</footnote>

<footnote n="2">² Laboratoire de Constitution et Réaction de la Matière, Université Felix Houphouët-Boigny, 22 BP 582 Abidjan 22, Côte d'Ivoire; lucie.bede@univ-fhb.edu.ci</footnote>

<footnote n="3">³ Université Paris-Saclay, CEA, CNRS, NIMBE, 91191 Gif-sur-Yvette, France; hicham.jabraoui@cea.fr</footnote>

<footnote n="4">⁴ Laboratoire Catalyse et Spectrochimie, Normandie Université, ENSICAEN, CNRS, 6 Boulevard Maréchal Juin, 14050 Caen, France; valentin.valtchev@ensicaen.fr</footnote>

<footnote>* Correspondence: etienne.hessou@univ-lorraine.fr; Tel.: +33-372749867</footnote>
various adsorption methods [24–30]. One of the most promising techniques for remov-
ing VOCs from indoor or outdoor environments is to trap them in inorganic adsorbents
such as zeolites [26,31–34]. These materials have been the subject of many applications
in adsorption and catalysis processes due to their unique properties, such as the ability
to exchange cations, a high specific surface area and a relatively inexpensive production
process [35]. In addition, they have good performance in terms of VOC adsorption capac-
ity, good thermal stability, and good selectivity. These various properties make zeolites
extremely popular materials for the adsorption of VOCs [31,36]. Among the 200 types of
existing zeolites [36,37], mordenite (MOR), ZSM-5 (MFI) and faujasite (FAU) are the most
used zeolites in the treatment of VOCs [26,33,34]. In addition, they also find applications
in biomedical fields. Toluene sensors for lung cancer detection require a concentration
step prior to detection; adsorbing the molecule in highly porous material, such as zeolites,
greatly improves these sensors [21].

In particular, faujasite zeolites are very often used in the separation of aromatics and hydrocarbons in petroleum processing due to the size of their pores going up to 7.4 Å [38]. They have been synthesized and exchanged with many different extra framework cations such as Na⁺, K⁺, Rb⁺, Cs⁺, Ag⁺, Ca²⁺ and Ba²⁺ [39,40]. In addition, faujasites are often used as a molecular sieve, an application that exploits their adsorption properties, and in the selective catalytic reduction process [26,33,41]. Accordingly, for many years now, faujasites have been used in the refining industry as the catalytically active component in cracking catalysts [42–44]. Previous works on the adsorption of xylene on faujasite have shown that the selectivity towards different mesomers of this molecule can be linked to the nature of cations [45–47]. Deng et al. [48] performed experimental and *ab initio* adsorption studies to compare the toluene adsorption capacity on FAU, BEA, and MFI zeolites where they reported superior selectivity of FAU USY towards toluene compared to water and dichloromethane. While Serra et al. [49,50] investigated the adsorption of toluene on different zeolite samples with various Si/Al ratios that were modified with varying loadings of exchangeable cations (Na⁺, H⁺ or Cs⁺). They showed that the thermal stability of the adsorbed toluene was mostly related to the type of the exchanged cation, which controls the overall basicity, rather than the structure of the zeolite. It arises that the most important factor is indeed the chemical nature of the chosen cations.

Wang and co-authors [51] conducted an experimental study to investigate the capture of VOCs by MFI-type zeolites in humid conditions, and they reported the competitive adsorption of toluene and water vapor. Liu et al. [52] showed that the presence of water could substantially influence the adsorption of toluene, stressing the importance of selective adsorbents of VOCs. Indeed, ambient air is a carrier gas, where a typical relative humidity (hr) of 60% corresponds to a water content of approximately 12.3 gm⁻³ or 1.66 vol-% at room temperature (T = 23 °C) [53]. However, to date, studies on VOCs-water mixtures' adsorption remain scarce [48].

In this study, we investigate the adsorption mechanisms of toluene and water in some
monovalent cation exchanged Y zeolites; simultaneously, we systematically identify the
cations that give the best adsorption capacity for the target molecule (i.e., toluene). We
use state-of-the-art Density Functional Theory (DFT) calculations, including dispersion
corrections. This approach has been proven to provide reliable results to elucidate the
mechanisms of adsorption in zeolite materials at the atomistic level [54–57]. In fact, nowa-
days, the chemical accuracy of DFT simulations is high enough to predict the enthalpies
of the adsorption of various molecules in zeolites, which agrees remarkably with the data
obtained in calorimetric experiments [55–59]. Periodic DFT approaches have also been
used to study the capture of various pollutants such as carbon dioxide, methane [60], or
iodine species from nuclear streams [61–63].

The paper is organized as follows: we present and discuss the results obtained for the type and strength of interactions between toluene/water and FAU Y exchanged with Li⁺, Na⁺, K⁺, Cs⁺, Cu⁺ and Ag⁺ cations. Then, a section dedicated to the calculation procedure is presented. The main conclusions drawn from this work are given in the last section.
## 2. Results

Using the DFT methodology as detailed in Section 3, we investigated the adsorption of toluene and water on different cation-exchanged Y faujasites. The preferential adsorption modes, as well as the interaction energies, of both molecules (toluene and water) over LiY, NaY, KY, CsY, CuY and AgY are presented and discussed below. The CONTAR files of most stable configuration calculated for toluene and water on different cation-exchanged Y faujasites are available in supplementary materials.

### 2.1. Toluene Adsorption

As can be seen from Figure 1, where we presented the most stable configurations of adsorbed toluene, the flat adsorption is the most stable mode found for toluene molecules adsorbed onto alkaline exchanged zeolites. The preference of this layout is explained by the fact that it allows π-interactions, which are known to be the most favorable adsorption mode between aromatic rings and cations in zeolites [64–66]. In addition, with alkalis, the cation acts as Lewis acid and, hence, it is expected to polarize the electron cloud of the aromatic ring [27]. Toluene is then adsorbed onto LiY, NaY and KY over one cation of Li, Na and K, respectively.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** The most stable adsorption modes calculated for toluene over LiY, NaY, KY, CsY, CuY and AgY. For the sake of clarity, only the 6-MR involved in the adsorption process was displayed (One 6-MR of LiY, NaY, KY, CuY and AgY; two 6-MR for CsY).</figcaption></figure>

In the case of LiY, the carbon atoms of the aromatic ring are located in the range of 2.57–3.12 Å (in average 2.84 Å) with respect to Li⁺, while the aromatic center is about 2.48 Å of the Li⁺ (see Table 1). In the case of NaY, the calculated distance between carbon atoms of the aromatic ring of the toluene and the interacting Na⁺ cation are in the range of 2.82–2.95 Å (in average 2.87 Å), and the aromatic center is about 2.5 Å of the Na⁺ cation.
For KY, the distances (K to the carbon of the aromatic ring or K to the aromatic center) are respectively 3.10–3.20 Å and 2.82 Å. One can observe that the increase of these average distances is proportional to increasing ionic radii of Li⁺, Na⁺ and K⁺ and, as can be expected, this same observation naturally extends to Cs⁺. In fact, in CsY, the increasing of the ionic radius allows the toluene molecule to interact with two Cs cations simultaneously while it only interacts with a single cation in the case of LiY, NaY, and KY. Toluene establishes a π-interaction with both Cs cations, where the average distance between the first Cs cation and the carbon atoms of the aromatic ring is 3.88 Å. The second cesium, diametrically opposed to the first, interacts with toluene at a distance of 3.4 Å from the aromatic center of toluene and an average of 3.68 Å from the carbon atoms of the aromatic ring.

<table><caption>**Table 1.** Distance between the aromatic ring of toluene and sorbents upon adsorption. C corresponds to the aromatic center of toluene; M = Li, Na, K, Cs1, Cs2, Cu, and Ag.</caption><thead><tr><th></th><th>LiY</th><th>NaY</th><th>KY</th><th>CsY</th><th>CuY</th><th>AgY</th></tr></thead><tbody><tr><td>M-C1</td><td>3.125</td><td>2.949</td><td>3.192</td><td>3.706</td><td>3.998</td><td>2.859</td></tr><tr><td>M-C2</td><td>3.095</td><td>2.896</td><td>3.132</td><td>3.729</td><td>3.840</td><td>2.052</td></tr><tr><td>M-C3</td><td>2.837</td><td>2.841</td><td>3.101</td><td>3.709</td><td>3.738</td><td>2.064</td></tr><tr><td>M-C4</td><td>2.569</td><td>2.821</td><td>3.116</td><td>3.657</td><td>3.782</td><td>2.866</td></tr><tr><td>M-C5</td><td>2.566</td><td>2.834</td><td>3.156</td><td>3.620</td><td>3.918</td><td>3.474</td></tr><tr><td>M-C6</td><td>2.844</td><td>2.891</td><td>3.189</td><td>3.639</td><td>4.020</td><td>3.460</td></tr><tr><td>M-C</td><td>2.478</td><td>2.504</td><td>2.819</td><td>3.396</td><td>3.619</td><td>2.481</td></tr></tbody></table>

The calculation of the electron density profiles confirms the geometrical analysis; where the double contribution of the Cs⁺ cations is clearly shown through the electron density map (see Figure 2a). This observation is totally different from what we found for Li⁺, Na⁺ and K⁺ cations, where the computed electron density shows that there is only one cation that can contribute to the adsorption of toluene in FAU, due to a shorter radius of these cations.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** The difference in electron density (Δρ) induced by the adsorption of toluene in CsY (a), LiY (b), and AgY (c) faujasites, respectively. The blue (green) zones indicate density increase (decrease). Where cyan balls represent Cs, light green balls represent Li, grey balls represent Ag, blue balls represent Si, and red balls represent O.</figcaption></figure>

The interaction energies of toluene with LiY, NaY, KY and CsY are 87.8/105.5/97.8/ 114.4 kJ/mol, respectively (see Table 2). An increase of the interaction energy is observed as we go down in the alkaline series (Li to Cs) since the electronegativity decreases considerably from Na to Cs. Although one would expect that the interaction energy of toluene with CsY would be weaker than that with NaY, the possibility of bidentate interaction explains the fact that toluene is 10 kJ/mol more adsorbed on CsY than NaY. A previous study of benzene adsorption on faujasite Y exhibits the same energy trends [67].
<table><caption>**Table 2.** Calculated interaction energies (kJ/mol) of toluene and water onto LiY, NaY, KY, CsY, CuY and AgY.</caption><thead><tr><th rowspan="2"></th><th colspan="2">ΔE<sub>int</sub> (kJ/mol)</th></tr><tr><th>Toluene</th><th>Water</th></tr></thead><tbody><tr><td>LiY</td><td>87.8</td><td>74.3</td></tr><tr><td>NaY</td><td>105.5</td><td>70.2</td></tr><tr><td>KY</td><td>97.8</td><td>68.8</td></tr><tr><td>CsY</td><td>114.4</td><td>67.3</td></tr><tr><td>CuY</td><td>147.0</td><td>79.1</td></tr><tr><td>AgY</td><td>131.5</td><td>76.6</td></tr></tbody></table>

Furthermore, we investigate the adsorption of toluene with CuY and AgY where transition metallic cations are embedded. The distances between the cation and the aromatic ring vary from 2.05 to 3.47 Å for CuY and from 2.37 to 3.57 Å in the case of AgY. Unlike alkaline cations, the interaction of toluene is established rather by two of the carbon atoms of the aromatic ring. The two distances M–C (M = Cu and Ag) are 2.05/2.06 Å for CuY and 2.37/2.38 Å for AgY. The interaction energies of toluene are 147.0/131.5 kJ/mol, respectively, for CuY and AgY. These two systems give the highest interaction energy of toluene. Here, we can see that only transition metal cation that can contribute to the adsorption of toluene in FAU. This is supported by the electron density map in which the charge of the second metal was not affected by the adsorption of toluene, evidencing the absence of any interaction (see Figure 2c).

## 2.2. Water Adsorption

Adsorption of water has been investigated on LiY, NaY, KY, CsY, CuY, and AgY. The most stable adsorption modes are displayed in Figure 3. We can see that the preferential adsorption mode of water is through its oxygen atom for all cases. In the case of systems with alkali metals, LY, NaY, KY, and CsY, the M–O distance between the metal and oxygen atom of water is 1.95/2.31/2.75/3.09 Å, respectively. This distance increases in proportion to the increase in the ionic radius from Li to Cs. A hydrogen bond is established during water adsorption between one hydrogen atom of water and the nearest oxygen atom of the 6MR. The H-bond distances are 2.04/1.90/1.86/1.92 Å, respectively, for LiY, NaY, KY, and CsY.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** The most stable adsorption modes calculated for water over LiY, NaY, KY, CsY, CuY and AgY. For the sake of clarity, only the 6-MR involved in the adsorption process was displayed. Color codes: silicon in blue and oxygen in red.</figcaption></figure>

The interaction energy calculated for water is 74.3 kJ/mol in LiY, 70.2 kJ/mol in NaY, 68.8 kJ/mol in KY, and 67.3 kJ/mol in CsY (see Table 2). The highest interaction is obtained
with LiY, and the smallest with CsY. These results are in agreement with experimental ones, as Dzhigit et al. [68] reported that the adsorption enthalpy of water at 23 °C is about 67 kJ/mol in LiNaY (44Li, 6Na); 66 kJ/mol in NaY and about 65 kJ/mol in KNaY (55 K, 2Na). The interaction energy decreases slightly from Li to Cs in line with the decrease in electronegativity along the alkali metal series. This observation can be explained by the HSAB theory of Pearson [69–71] and was described in our previous work with a higher Si/Al ratio [72,73]. In addition, the experimental work of Jentys et al. [74] has shown that in alkali-metal-exchanged ZSM-5, the adsorption of water molecules at low coverage on alkali-metal cations are more favored for Li-, Na-, and K-ZSM5 than for Cs- and Rb-ZSM5.

In the case of CuY and AgY, the interactions of the metallic atom (M) and the oxygen atom in water are characterized with M–O distances of 1.91 Å and 2.26 Å, respectively. In comparison, one water H atom is involved in a hydrogen bond with an oxygen atom of the 6MR at 1.75 Å and 1.81 Å for CuY and AgY, respectively. The interaction energy trends show that both structures yield similar adsorption behavior where interaction energies are 79.1 kJ/mol and 76.6 kJ/mol for CuY and AgY, respectively. In previous work, we have found a similar trend with an Si/Al ratio of 47 [72], where water adsorption underlines close interaction energies with an only slight preference for CuY and AgY compared to alkali ones (LiY, NaY, KY and CsY). We point out that this contrasts with the results reported by Kumar et al. [67]. Their damped Car–Parrinello molecular dynamics results showed that the adsorption energies of water in CuY and AgY are about 65 and 70 kJ/mol, respectively. According to their results, AgY is expected to adsorb water more efficiently than CuY. However, these differences remain close to the expected chemical accuracy of the employed methods.

**3. Discussion and Conclusions**

In the previous sections, we have presented individual details on the adsorption of toluene (Section 2.1) and water (Section 2.2) in some cationic-exchanged forms of zeolite Y (LiY, NaY, KY, CsY, CuY and AgY). At the base of these results, it is now necessary to be interested in the probable effect of water on the adsorption of the target molecule (toluene).

As shown in Table 2 and Figure 4, toluene is more adsorbed than water on all the systems studied. However, the energy gap between the adsorption of toluene and that of water is a crucial factor that may well inform the choice of the systems.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Variation of the interaction energy (kJ/mol) of toluene and water depending on the adsorbent LiY, NaY, KY, CsY, CuY, and AgY.</figcaption></figure>
Therefore, by comparing the interaction energies calculated for toluene and water, the following observations can be made:

- LiY adsorbs more toluene than water with a difference of more than 13 kJ/mol;

- NaY adsorbs more toluene than water with a difference of more than 35 kJ/mol;

- KY adsorbs toluene a little less than NaY, but the toluene/water difference is 29kJ/mol;

- CsY adsorbs toluene more than water with a difference of 47 kJ/mol, thus putting CsY at the top of the list of Y zeolites exchanged for alkali cations;

- CuY and AgY adsorb toluene more than water with 68 and 58 kJ/mol, respectively.

One can therefore classify the investigated systems as follows, CuY > AgY > CsY > NaY > KY > LiY, and conclude that CuY, AgY, CsY and KY are, in decreasing order, the best options for effective selectivity in the adsorption of toluene in the presence of water vapor.

In order to provide a fundamental understanding governing the interaction between toluene/water and several formulations of FAU, we limit ourselves only to the adsorption of individual molecules without considering the impact of one molecule on another. This opens new perspectives for in-depth studies on the adsorption of these two molecules, including experimental investigations or co-adsorption molecular dynamics simulations.

# 4. Computational Methods

## 4.1. Computation Details

We investigated the toluene and water adsorption onto Y zeolites using the Vienna Ab-initio Simulation Package (VASP) [75]. DFT calculations employed the PBE functional and projector augmented plane wave (PAW) method [76,77] with plane-wave cutoff energy of 450 eV. Kohn-Sham equations were solved iteratively until the energy difference between cycles became lower than $10^{-6}$ eV, while Gaussian smearing was fixed to 0.1 eV. The Brillouin zone was sampled at the Γ-point only. FAU structures were fully optimized until all the forces fell below 0.01 eV/Å per atom. In order to account for the van der Waals (vdW) interactions, the TS/HI method [78] was used. The TS/HI scheme is an improved version of the TS (Tkatchenko-Scheffler) [79] dispersion correction; it is based on an Iterative Hirshfeld partitioning leading to an accurate description of ionic solids such as cationic zeolites [78–80].

To describe the adsorption phenomena, the interaction energy between the molecule (Toluene/Water) and the Y zeolite formulations were evaluated at 0 K using the following equation [67,81]:

$$ \Delta E_{\text{int}} = E_{\text{zeolite}} + E_{\text{molecule}} - E_{\text{zeolite-molecule}} \quad (1) $$

where $E_{\text{zeolite-molecule}}$ is the energy of the Y zeolite with adsorbed molecules, $E_{\text{zeolite}}$ is the energy of the empty zeolite, and $E_{\text{molecule}}$ is the energy of the isolated molecule in the gaseous phase. According to this equation, a positive value of $\Delta E_{\text{int}}$ corresponds to an exothermic process.

In addition to adsorption investigations in terms of the total energy interaction, we are using the charge density difference ($\Delta\rho$) to further enhance our understanding of the adsorption of molecules into the considered faujasite formulations. To visualize the charge density difference ($\Delta\rho$) introduced by the adsorption of toluene, we combine the three charge densities: (1) the density of the complex ($\rho_{\text{Fau-molecule}}$); (2) the density of the isolated molecule ($\rho_{\text{molecule}}$); and (3) the density of the clean faujasite $\rho_{\text{FAU}}$, which is formulated as follows:

$$ \Delta\rho = \rho_{\text{FAU-molecule}} - \rho_{\text{FAU}} - \rho_{\text{molecule}}. \quad (2) $$

## 4.2. Structural Model

The siliceous structure of faujasite crystallizes within the Fd3m symmetry space group [82]. The FAU-type framework consists of cuboctahedral sodalite units connected by hexagonal prisms (D6R), forming large empty cavities called supercages. These supercages are interconnected by their 12-membered ring windows, also called hexagonal windows (12MR), thereby forming the porous network. The standard cell has a cubic structure with
the lattice parameters *a* = *b* = *c* = 25.028 Å (576 atoms, Si<sub>192</sub>O<sub>384</sub>) and α = β = γ = 90° [83,84]. However, to reduce the computational effort, we have used a primitive rhombohedral cell (two supercages and eight hexagonal windows connecting the sodalite cages to the supercages) containing 144 atoms. In this structure, 14 atoms of silicon Si(+IV) are replaced by 14 atoms of aluminum Al(+III) to thus obtain a Y faujasite with an Si/Al ratio equal to 2.43. In order to once again reach the neutrality of the system lost due to the introduction of Al(+III), 14 monovalent cations of the same atom are added to obtain LiY, NaY, KY, CsY, CuY and AgY, respectively.

**Supplementary Materials:** The CONTCAR files are available online.

**Author Contributions:** Conceptualization, E.P.H., L.A.B. and M.B.; methodology, E.P.H., H.J. and A.S.; software, E.P.H. and L.A.B.; validation, E.P.H., M.B. and V.V.; formal analysis, E.P.H., L.A.B., H.J. and A.S.; investigation, E.P.H., L.A.B. and H.J.; resources, M.B. and V.V.; data curation, E.P.H., L.A.B. and H.J.; writing—original draft preparation, E.P.H., A.S. and H.J.; writing—review and editing, E.P.H., L.A.B., H.J., A.S., M.B. and V.V.; visualization, E.P.H.; supervision, E.P.H., M.B. and V.V.; project administration, E.P.H., M.B. and V.V.; funding acquisition, M.B. and V.V. All authors have read and agreed to the published version of the manuscript.

**Funding:** H.J. and M.B. acknowledge financial support through the COMETE project (COception in silico de Matériaux pour l'Environnement et l'Energie) co-funded by the European Union under the program "FEDER-FSE Lorraine et Massif des Vosges 2014-2020".

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available in supplementary material.

**Acknowledgments:** This work was granted access to the HPC resources of TGCC under the allocations 2020-A0080810169 and 2021-A010081433 by GENCI-EDARI project. E.P.H. and M.B. acknowledge financial support through the COMETE project (COception in silico de Matériaux pour l'Environnement et l'Energie) co-funded by the European Union under the program "FEDER-FSE Lorraine et Massif des Vosges 2014-2020".

**Conflicts of Interest:** The authors declare no conflict of interest.

**Sample Availability:** Zeolite models are available in Supplementary Materials.

## References

1. Chen, X.; Qian, W.; Kong, L.; Xiong, Y.; Tian, S. Performance of a suspended biofilter as a new bioreactor for removal of toluene. *Biochem. Eng. J.* **2015**, *98*, 56–62. [CrossRef]

2. Liang, W.-J.; Ma, L.; Liu, H.; Li, J. Toluene degradation by non-thermal plasma combined with a ferroelectric catalyst. *Chemosphere* **2013**, *92*, 1390–1395. [CrossRef]

3. Yeom, S.H.; Daugulis, A.J.; Nielsen, D.R. A strategic approach for the design and operation of two-phase partitioning bioscrubbers for the treatment of volatile organic compounds. *Biotechnol. Prog.* **2010**, *26*, 1777–1786. [CrossRef]

4. Banerjee, T.; Srivastava, R.K. Evaluation of environmental impacts of Integrated Industrial Estate—Pantnagar through application of air and water quality indices. *Environ. Monit. Assess.* **2010**, *172*, 547–560. [CrossRef] [PubMed]

5. Borbon, A.; Gilman, J.B.; Kuster, W.; Grand, N.; Chevaillier, S.; Colomb, A.; Dolgorouky, C.; Gros, V.; Lopez, M.; Sarda-Esteve, R.; et al. Emission ratios of anthropogenic volatile organic compounds in northern mid-latitude megacities: Observations versus emission inventories in Los Angeles and Paris. *J. Geophys. Res. Atmos.* **2013**, *118*, 2041–2057. [CrossRef]

6. Abtahi, M.; Fakhri, Y.; Conti, G.O.; Ferrante, M.; Taghavi, M.; Tavakoli, J.; Heshmati, A.; Keramati, H.; Moradi, B.; Amanidaz, N.; et al. The Concentration of BTEX in the Air of Tehran: A Systematic Review-Meta Analysis and Risk Assessment. *Int. J. Environ. Res. Public Health* **2018**, *15*, 1837. [CrossRef] [PubMed]

7. Guerra, L.R.; Souza, A.; Côrtes, J.A.; Lione, V.D.O.F.; Castro, H.C.; Alves, G.G. Assessment of predictivity of volatile organic compounds carcinogenicity and mutagenicity by freeware in silico models. *Regul. Toxicol. Pharmacol.* **2017**, *91*, 1–8. [CrossRef]

8. Rene, E.R.; Kar, S.; Krishnan, J.; Pakshirajan, K.; López, M.E.; Murthy, D.; Swaminathan, T. Start-up, performance and optimization of a compost biofilter treating gas-phase mixture of benzene and toluene. *Bioresour. Technol.* **2015**, *190*, 529–535. [CrossRef]

9. Yeom, S.-H.; Daugulis, A.J. Development of a novel bioreactor system for treatment of gaseous benzene. *Biotechnol. Bioeng.* **2001**, *72*, 156–165. [CrossRef]

10. Cámara-Lemarroy, C.R.; González-Moreno, E.I.; Rodriguez-Gutierrez, R.; González-González, J.G. Clinical presentation and management in acute toluene intoxication: A case series. *Inhal. Toxicol.* **2012**, *24*, 434–438. [CrossRef]
11. Knox, J.W.; Nelson, J.R. Permanent Encephalopathy from Toluene Inhalation. *N. Engl. J. Med.* **1966**, *275*, 1494–1496. [CrossRef] [PubMed]

12. Streicher, H.Z.; Gabow, P.A.; Moss, A.H.; Kono, D.; Kaehny, W.D. Syndromes of Toluene Sniffing in Adults. *Ann. Intern. Med.* **1981**, *94*, 758–762. [CrossRef] [PubMed]

13. Carlisle, E.J.; Donnelly, S.M.; Vasuvattakul, S.; Kamel, K.S.; Tobe, S.; Halperin, M.L. Glue-sniffing and distal renal tubular acidosis: Sticking to the facts. *J. Am. Soc. Nephrol.* **1991**, *1*, 1019–1027. [CrossRef]

14. Ran, J.; Qiu, H.; Sun, S.; Tian, L. Short-term effects of ambient benzene and TEX (toluene, ethylbenzene, and xylene combined) on cardiorespiratory mortality in Hong Kong. *Environ. Int.* **2018**, *117*, 91–98. [CrossRef] [PubMed]

15. Benignus, V.A.; Boyes, W.K.; Kenyon, E.; Bushnell, P.J. Quantitative Comparisons of the Acute Neurotoxicity of Toluene in Rats and Humans. *Toxicol. Sci.* **2007**, *100*, 146–155. [CrossRef] [PubMed]

16. Ran, J.; Sun, S.; Yang, A.; Yang, L.; Han, L.; Mason, T.G.; Chan, K.-P.; Li, J.; Tian, L. Effects of ambient benzene and toluene on emergency COPD hospitalizations: A time series study in Hong Kong. *Sci. Total Environ.* **2019**, *657*, 28–35. [CrossRef] [PubMed]

17. Lin, C.-M.; Liu, C.-K. Reversible Cerebral Periventricular White Matter Changes with Corpus Callosum Involvement in Acute Toluene-Poisoning. *J. Neuroimaging* **2014**, *25*, 497–500. [CrossRef]

18. Kao, H.-W.; Pare, L.; Kim, R.; Hasso, A.N. Toxic leukoencephalopathy with atypical MRI features following a lacquer thinner fire. *J. Clin. Neurosci.* **2014**, *21*, 878–880. [CrossRef]

19. Tsatsakis, A.; Dolapsakis, G.; Troulakis, G.; Christodoulou, P.; Relakis, K.; Trikilis, N.; Michalodimitrakis, M. Fatal and non-fatal outcome by accidental intoxication with paint thinner. *J. Clin. Forensic Med.* **1997**, *4*, 133–137. [CrossRef]

20. Park, J.T.; Jung, H.H.; Choi, K.H. Cerebellar infarction following accidental inhalation of toluene-mixed paint. *Am. J. Emerg. Med.* **2018**, *36*, 1924.e5–1924.e7. [CrossRef]

21. Gregis, G.; Schaefer, S.; Sanchez, J.-B.; Fierro, V.; Berger, F.; Bezverkhyy, I.; Weber, G.; Bellat, J.-P.; Celzard, A. Characterization of materials toward toluene traces detection for air quality monitoring and lung cancer diagnosis. *Mater. Chem. Phys.* **2017**, *192*, 374–382. [CrossRef]

22. Kim, N.-H.; Choi, S.-J.; Yang, D.-J.; Bae, J.; Park, J.; Kim, I.-D. Highly sensitive and selective hydrogen sulfide and toluene sensors using Pd functionalized WO3 nanofibers for potential diagnosis of halitosis and lung cancer. *Sens. Actuators B Chem.* **2014**, *193*, 574–581. [CrossRef]

23. Gregis, G.; Sanchez, J.-B.; Bezverkhyy, I.; Guy, W.; Berger, F.; Fierro, V.; Bellat, J.-P.; Celzard, A. Detection and quantification of lung cancer biomarkers by a micro-analytical device using a single metal oxide-based gas sensor. *Sens. Actuators B Chem.* **2018**, *255*, 391–400. [CrossRef]

24. Bal'Zhinimaev, B.S.; Paukshtis, E.A.; Toktarev, A.V.; Kovalyov, E.V.; Yaranova, M.A.; Smirnov, A.E.; Stompel, S. Effect of water on toluene adsorption over high silica zeolites. *Microporous Mesoporous Mater.* **2019**, *277*, 70–77. [CrossRef]

25. Zaitan, H.; Mohamed, E.F.; Valdés, H.; Nawdali, M.; Rafqah, S.; Manero, M.H. Toluene, Methanol and Benzaldehyde Removal from Gas Streams by Adsorption onto Natural Clay and Faujasite-Y type Zeolite. *Acta Chim. Slov.* **2016**, *63*, 798–808. [CrossRef]

26. Romero, D.; Chlala, D.; Labaki, M.; Royer, S.; Bellat, J.-P.; Bezverkhyy, I.; Giraudon, J.-M.; Lamonier, J.-F.; Chlala, D. Removal of Toluene over NaX Zeolite Exchanged with Cu²⁺. *Catalysts* **2015**, *5*, 1479–1497. [CrossRef]

27. Simperler, A.; Bell, R.G.; Philippou, A.; Anderson, M. Theoretical Study of Toluene Adsorbed on Zeolites X and Y: Calculation of 13C NMR Parameters. *J. Phys. Chem. B* **2002**, *106*, 10944–10954. [CrossRef]

28. Lee, H.C.; Kim, K.M.; Choi, S.I.; Kim, Y.H.; Woo, H.C.; Won, Y.S. A Computational Study on the Adsorption Characteristics of Hydrocarbons (Propylene, n-Butane and Toluene) by Uing Cation-Exchanged ZSM-5 Zeo-lites. *Korean Chem. Eng. Res.* **2018**, *56*, 909–913. [CrossRef]

29. Sacchetto, V.; Bisio, C.; Olivera, D.F.O.; Paul, G.; Gatti, G.; Braschi, I.; Berlier, G.; Cossi, M.; Marchese, L. Interactions of Toluene and n-Hexane on High Silica Zeolites: An Experimental and Computational Model Study. *J. Phys. Chem. C* **2015**, *119*, 24875–24886. [CrossRef]

30. Oh, G.-Y.; Ju, Y.-W.; Kim, M.-Y.; Jung, H.-R.; Kim, H.J.; Lee, W.-J. Adsorption of toluene on carbon nanofibers prepared by electrospinning. *Sci. Total Environ.* **2008**, *393*, 341–347. [CrossRef]

31. Guisnet, M.; Gilson, J.-P. (Eds.) *Zeolites for Cleaner Technologies, Catalytic Science Series*; Imperial College Press: London, UK, 2002; ISBN 978-1-86094-329-4.

32. Dyer, A. Bish D. L. & Ming D. W. (Editors) Natural Zeolites: Occurrences, Properties, Applications.: Reviews in Mineralogy and Geochemistry, 45. 2001, 662 pp. Price US$32. ISBN 1529-6466. *Clay Minerals* **2002**, *37*, 733.

33. de Fonseca, B.; Rossignol, J.; Bezverkhyy, I.; Bellat, J.; Stuerga, D.; Pribetich, P. Detection of VOCs by microwave transduction using dealuminated faujasite DAY zeolites as gas sensitive materials. *Sens. Actuators B Chem.* **2015**, *213*, 558–565. [CrossRef]

34. Brunchi, C.C.; Sanchez, J.M.C.; Stankiewicz, A.I.; Kramer, H.J.M.; Vlugt, T.J.H. Adsorption of Volatile Organic Compounds. Experimental and Theoretical Study. *Ind. Eng. Chem. Res.* **2012**, *51*, 16697–16708. [CrossRef]

35. Góra-Marek, K.; Palomares, A.; Glanowska, A.; Sadowska, K.; Datka, J. Copper Sites in Zeolites—Quantitative IR studies. *Microporous Mesoporous Mater.* **2012**, *162*, 175–180. [CrossRef]

36. Jacobs, P.A.; Flanigen, E.M.; Jansen, J.C.; van Bekkum, H. *Introduction to Zeolite Science and Practice*; Elsevier: Amsterdam, The Netherlands, 2001; ISBN 0-444-82421-9.

37. Auerbach, S.M.; Carrado, K.A.; Dutta, P.K. *Handbook of Zeolite Science and Technology*; CRC Press: Boca Raton, FL, USA, 2003.
38. Khalil, I.; Thomas, K.; Jabraoui, H.; Bazin, P.; Maugé, F. Selective elimination of phenol from hydrocarbons by zeolites and silica-based adsorbents—Impact of the textural and acidic properties. *J. Hazard. Mater.* **2020**, *384*, 121397. [CrossRef] [PubMed]

39. Kulprathipanja, S. Process for Adsorptive Separation of Metaxylene from Xylene Mixtures. U.S. Patent No. 5,900,523, 4 May 1999.

40. Roeseler, C.M.; Kulprathipanja, S.; Rekoske, J.E. Adsorptive Separation Process for Recovery of Para-Xylene. U.S. Patent No 6,706,938, 16 March 2004.

41. Hessou, E.; Jabraoui, H.; Chebbi, M.; Chibani, S.; Cantrel, L.; Badawi, M. Evaluation of the Inhibiting Effect of Organic Compounds on the Adsorption of Iodine Compounds in Cation-Exchanged Zeolites: A DFT Study. In Proceedings of the Recent Advances in Environmental Science from the Euro-Mediterranean and Surrounding Regions, Sousse, Tunisia, 22–25 November 2017; Kallel, A., Ksibi, M., Ben Dhia, H., Khélifi, N., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 107–109.

42. Pine, L. Prediction of cracking catalyst behavior by a zeolite unit cell size model. *J. Catal.* **1984**, *85*, 466–476. [CrossRef]

43. Verboekend, D.; Nuttens, N.; Locus, R.; Van Aelst, J.; Verolme, P.; Groen, J.C.; Pérez-Ramírez, J.; Sels, B.F. Synthesis, characterization, and catalytic evaluation of hierarchical faujasite zeolites: Milestones, challenges, and future directions. *Chem. Soc. Rev.* **2015**, *45*, 3331–3352. [CrossRef]

44. Kulprathipanja, S. *Zeolites in Industrial Separation and Catalysis*; John Wiley & Sons: Hoboken, NJ, USA, 2010; ISBN 978-3-527-62957-2.

45. Bellat, J.-P.; Simonot-Grange, M.-H. Adsorption of gaseous p-xylene and m-xylene on NaY, KY, and BaY zeolites. Part 2: Modeling. Enthalpies and entropies of adsorption. *Zeolites* **1995**, *15*, 219–227. [CrossRef]

46. Cottier, V.; Bellat, J.-P.; Simonot-Grange, M.-H.; Méthivier, A. Adsorption of p-Xylene/m-Xylene Gas Mixtures on BaY and NaY Zeolites. Coadsorption Equilibria and Selectivities. *J. Phys. Chem. B* **1997**, *101*, 4798–4802. [CrossRef]

47. Hulme, R.; Rosensweig, R.E.; Ruthven, D.M. Binary and ternary equilibria for C8 aromatics on KY faujasite. *Ind. Eng. Chem. Res.* **1991**, *30*, 752–760. [CrossRef]

48. Deng, H.; Pan, T.; Zhang, Y.; Wang, L.; Wu, Q.; Ma, J.; Shan, W.; He, H. Adsorptive removal of toluene and dichloromethane from humid exhaust on MFI, BEA and FAU zeolites: An experimental and theoretical study. *Chem. Eng. J.* **2020**, *394*, 124986. [CrossRef]

49. Serra, R.M.; Miró, E.E.; Boix, A.V. FTIR study of toluene adsorption on Cs-exchanged mordenites. *Microporous Mesoporous Mater.* **2010**, *127*, 182–189. [CrossRef]

50. Serra, R.M.; Miró, E.E.; Sapag, M.K.; Boix, A.V. Adsorption and diffusion of toluene on Na and Cs mordenites for hydrocarbon traps. *Microporous Mesoporous Mater.* **2011**, *138*, 102–109. [CrossRef]

51. Wang, S.; Bai, P.; Wei, Y.; Liu, W.; Ren, X.; Bai, J.; Lu, Z.; Yan, W.; Yu, J. Three-Dimensional-Printed Core–Shell Structured MFI-Type Zeolite Monoliths for Volatile Organic Compound Capture under Humid Conditions. *ACS Appl. Mater. Interfaces* **2019**, *11*, 38955–38963. [CrossRef]

52. Liu, S.; Peng, Y.; Chen, J.; Shi, W.; Yan, T.; Li, B.; Zhang, Y.; Li, J. Engineering surface functional groups on mesoporous silica: Towards a humidity-resistant hydrophobic adsorbent. *J. Mater. Chem. A* **2018**, *6*, 13769–13777. [CrossRef]

53. Kraus, M.; Trommler, U.; Holzer, F.; Kopinke, F.-D.; Roland, U. Competing adsorption of toluene and water on various zeolites. *Chem. Eng. J.* **2018**, *351*, 356–363. [CrossRef]

54. Van Speybroeck, V.; Hemelsoet, K.; Joos, L.; Waroquier, M.; Bell, R.G.; Catlow, C.R.A. Advances in theory and their application within the field of zeolite chemistry. *Chem. Soc. Rev.* **2015**, *44*, 7044–7111. [CrossRef] [PubMed]

55. De Moor, B.A.; Reyniers, M.-F.; Gobin, O.C.; Lercher, J.A.; Marin, G.B. Adsorption of C2–C8 n-Alkanes in Zeolites. *J. Phys. Chem. C* **2011**, *115*, 1204–1219. [CrossRef]

56. Göltl, F.; Grüneis, A.; Bučko, T.; Hafner, J. Van der Waals interactions between hydrocarbon molecules and zeolites: Periodic calculations at different levels of theory, from density functional theory to the random phase approximation and Møller-Plesset perturbation theory. *J. Chem. Phys.* **2012**, *137*, 114111. [CrossRef] [PubMed]

57. Piccini, G.; Alessio, M.; Sauer, J.; Zhi, Y.; Liu, Y.; Kolvenbach, R.; Jentys, A.; Lercher, J.A. Accurate Adsorption Thermodynamics of Small Alkanes in Zeolites. Ab initio Theory and Experiment for H-Chabazite. *J. Phys. Chem. C* **2015**, *119*, 6128–6137. [CrossRef]

58. Khalil, I.; Jabraoui, H.; Maurin, G.; Lebègue, S.; Badawi, M.; Thomas, K.; Maugé, F. Selective Capture of Phenol from Biofuel Using Protonated Faujasite Zeolites with Different Si/Al Ratios. *J. Phys. Chem. C* **2018**, *122*, 26419–26429. [CrossRef]

59. Khalil, I.; Jabraoui, H.; Lebègue, S.; Kim, W.J.; Aguilera, L.-J.; Thomas, K.; Maugé, F.; Badawi, M. Biofuel purification: Coupling experimental and theoretical investigations for efficient separation of phenol from aromatics by zeolites. *Chem. Eng. J.* **2020**, *402*, 126264. [CrossRef]

60. Fischer, M.; Bell, R.G. A dispersion-corrected density-functional theory study of small molecules adsorbed in alkali-exchanged chabazites. *Z. Krist. Cryst. Mater.* **2013**, *228*, 124–133. [CrossRef]

61. Chebbi, M.; Chibani, S.; Paul, J.-F.; Cantrel, L.; Badawi, M. Evaluation of volatile iodine trapping in presence of contaminants: A periodic DFT study on cation exchanged-faujasite. *Microporous Mesoporous Mater.* **2017**, *239*, 111–122. [CrossRef]

62. Chibani, S.; Chebbi, M.; Lebègue, S.; Bučko, T.; Badawi, M. A DFT investigation of the adsorption of iodine compounds and water in H-, Na-, Ag-, and Cu- mordenite. *J. Chem. Phys.* **2016**, *144*, 244705. [CrossRef] [PubMed]

63. Chibani, S.; Medlej, I.; Lebègue, S.; Ángyán, J.G.; Cantrel, L.; Badawi, M. Performance of CuII-, PbII-, and HgII-Exchanged Mordenite in the Adsorption of I2, ICH3, H2O, CO, ClCH3, and Cl2: A Density Functional Study. *ChemPhysChem* **2017**, *18*, 1642–1652. [CrossRef]
64. Jabraoui, H.; Hessou, E.; Chibani, S.; Cantrel, L.; Lebègue, S.; Badawi, M. Adsorption of volatile organic and iodine compounds over silver-exchanged mordenites: A comparative periodic DFT study for several silver loadings. *Appl. Surf. Sci.* **2019**, *485*, 56–63.
[CrossRef]

65. Hessou, E.P.; Ponce-Vargas, M.; Mensah, J.-B.; Tielens, F.; Santos, J.C.; Badawi, M. Dibenzyl Disulfide Adsorption on Cationic Exchanged Faujasites: A DFT Study. *Nanomaterials* **2019**, *9*, 715. [CrossRef]

66. Hessou, E.P.; Jabraoui, H.; Hounguè, M.T.A.K.; Mensah, J.-B.; Pastore, M.; Badawi, M. A first principle evaluation of the adsorption mechanism and stability of volatile organic compounds into NaY zeolite. *Z. Krist. Cryst. Mater.* **2019**, *234*, 469–482. [CrossRef]

67. Hessou, E.P.; Jabraoui, H.; Khalil, I.; Dziurla, M.-A.; Badawi, M. Ab initio screening of zeolite Y formulations for efficient adsorption of thiophene in presence of benzene. *Appl. Surf. Sci.* **2021**, *541*, 148515. [CrossRef]

68. Dzhigit, O.M.; Kiselev, A.V.; Rachmanova, T.A.; Zhdanov, S.P. Influence of Li⁺, Na⁺ and K⁺ cation concentrations in X and Y zeolites on isotherms and heats of adsorption of propane and water. *J. Chem. Soc. Faraday Trans. 1 Phys. Chem. Condens. Phases* **1979**, *75*, 2662–2677. [CrossRef]

69. Parr, R.G.; Pearson, R.G. Absolute hardness: Companion parameter to absolute electronegativity. *J. Am. Chem. Soc.* **1983**, *105*, 7512–7516. [CrossRef]

70. Pearson, R.G. Hard and soft acids and bases, HSAB, part 1: Fundamental principles. *J. Chem. Educ.* **1968**, *45*, 581–587. [CrossRef]

71. Pearson, R.G. Absolute electronegativity and hardness: Application to inorganic chemistry. *Inorg. Chem.* **1988**, *27*, 734–740. [CrossRef]

72. Hessou, E.P.; Kanhounnon, W.G.; Rocca, D.; Monnier, H.; Vallières, C.; Lebègue, S.; Badawi, M. Adsorption of NO, NO₂, CO, H₂O and CO₂ over isolated monovalent cations in faujasite zeolite: A periodic DFT investigation. *Theor. Chem. Acc.* **2018**, *137*, 161. [CrossRef]

73. Jabraoui, H.; Khalil, I.; Lebègue, S.; Badawi, M. Ab initio screening of cation-exchanged zeolites for biofuel purification. *Mol. Syst. Des. Eng.* **2019**, *4*, 882–892. [CrossRef]

74. Jentys, A.; Warecka, G.; Derewinski, M.; Lercher, J.A. Adsorption of water on ZSM 5 zeolites. *J. Phys. Chem.* **1989**, *93*, 4837–4843. [CrossRef]

75. Sun, G.; Kürti, J.; Rajczy, P.; Kertesz, M.; Hafner, J.; Kresse, G. Performance of the Vienna ab initio simulation package (VASP) in chemical applications. *J. Mol. Struct. THEOCHEM* **2003**, *624*, 37–45. [CrossRef]

76. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. *Phys. Rev. B* **1999**, *59*, 1758. [CrossRef]

77. Blöchl, P.E. Projector augmented-wave method. *Phys. Rev. B* **1994**, *50*, 17953–17979. [CrossRef]

78. Bucko, T.; Lebègue, S.; Ángyán, J.G.; Hafner, J. Extending the applicability of the Tkatchenko-Scheffler dispersion correction via iterative Hirshfeld partitioning. *J. Chem. Phys.* **2014**, *141*, 034114. [CrossRef]

79. Tkatchenko, A.; Scheffler, M. Accurate Molecular Van Der Waals Interactions from Ground-State Electron Density and Free-Atom Reference Data. *Phys. Rev. Lett.* **2009**, *102*, 073005. [CrossRef] [PubMed]

80. Bučko, T.; Lebègue, S.; Hafner, J.; Ángyán, J.G. Improved Density Dependent Correction for the Description of London Dispersion Forces. *J. Chem. Theory Comput.* **2013**, *9*, 4293–4299. [CrossRef] [PubMed]

81. Jabraoui, H.; Charpentier, T.; Gin, S.; Delaye, J.-M.; Pollet, R. Atomic Insights into the Events Governing the Borosilicate Glass–Water Interface. *J. Phys. Chem. C* **2021**, *125*, 7919–7931. [CrossRef]

82. Baur, W.H. On the Cation and Water Positions in Faujasite. *Am. Mineral.* *J. Earth Planet. Mater.* **1964**, *49*, 697–704.

83. Dempsey, E.; Kuehl, G.H.; Olson, D.H. Variation of the lattice parameter with aluminum content in synthetic sodium faujasites. Evidence for ordering of the framework ions. *J. Phys. Chem.* **1969**, *73*, 387–390. [CrossRef]

84. Porcher, F.; Souhassou, M.; Dusausoy, Y.; Lecomte, C. The crystal structure of a low-silica dehydrated NaX zeolite. *Eur. J. Miner.* **1999**, *11*, 333–344. [CrossRef]