Sorption Mechanisms and Enhancement of Selected or
Sorption Mechanisms and Enhancement of Selected or
Journal of Chemistry
Volume 2024, Article ID 4008315, 19 pages
https://doi.org/10.1155/2024/4008315
Research Article
Sorption Mechanisms and Enhancement of Selected
Organochlorine Pollutants in Water on Zeolites
Received 14 April 2023; Revised 29 November 2023; Accepted 23 December 2023; Published 8 January 2024
Copyright © 2024 Timothy Manda et al. Tis is an open access article distributed under the Creative Commons Attribution
License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly
cited.
Tis study investigates the adsorption capacities of selected organochlorines on zeolites, focusing on hexachlorobenzene
(HCB), hexachlorotetradecane (HCTD), hexachlorodecane (HCD), hexachlorocyclohexane (HCH), heptachlorodecane
(HPCD), octachlorodecane (OCD), dichlorodiphenyltrichloroethane (DDT), and octachlorotetradecane (OCTD). Te
structures of the organochlorines were optimized and their Frontier molecular orbitals were calculated. Te analysis of
HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) energies provided insights
into the molecules’ electron-donating and -accepting capabilities. Te present research identifed the universal force feld as
suitable for the investigation and used it to evaluate the adsorption capacities of the pollutants on various zeolites. It was
found that CLO (a cubic microporous gallophosphate) demonstrated the highest adsorption capacity for HCB among 245
zeolites, with a loading capacity of 65.84 wt%. In terms of molecules adsorbed per cell, CLO remained the highest with 120
molecules per cell for HCB, 113 molecules per cell for HCH, 43 molecules per cell for DDT, 21 molecules per cell for HCTD,
19 molecules per cell for OCTD, 47 molecules per cell for HCD, 30 molecules per cell for HPCD, and 22 molecules per cell for
OCD. Te analysis revealed correlations between the structural parameters of zeolites (mass, density, HVF, APV, VSA, GSA,
DPS, and Di) and their adsorption capacities. Te investigation delved into cluster models to understand the interaction of
organochlorines with the zeolite framework. Te study explored the impact of doping CLO zeolite with diferent atoms (Al,
Si, and Na) on adsorption capacity. Te results showed that doping with aluminum improved both loading capacity and
adsorption energy and dissociate the chlorinated compounds during adsorption. Quantum chemical calculations show that
hydrogen-based bonding of the organochlorides on the CLO is thermodynamically favorable compared to dissociative
adsorption. In addition, oxygen atoms in the zeolites provide active adsorption sites. In the present work, laboratory
adsorption experiments were performed, treating zeolites with heat at 400°C. Surprisingly, untreated zeolites outperformed
treated ones, adsorbing up to 91% of HCB, while treated zeolites reached saturation after the third run. Te study attributed
the better performance of untreated zeolites to the presence of interstitial water and hydrogen atoms, which are critical for
electrostatic interactions with organic compounds. In general, this research provides a comprehensive analysis of the
adsorption capacities of organochlorines on zeolites, combining computational simulations and laboratory experiments.
2 Journal of Chemistry
Tis work’s distinctive quality is its methodology that combines molecular simulations, experimental verifcation, doping, and
interstitial water efects. Te fndings emphasize the importance of zeolite (a high-porosity nanostructured material) structure,
composition, and treatment methods in determining their efectiveness as adsorbents for environmental pollutants.
processes play crucial roles in environmental remediation, with Table 1: List of pollutants with their calculated zeolite-guest in-
adsorption providing immediate pollutant removal and deg- teraction energies for zeolite MOR and FAU obtained from
radation ofering a long-term solution by breaking down literature.
pollutants into less toxic forms. Te choice of materials, their MOR FAU
properties, and the specifc pollutants targeted are critical factors Pollutant
Ezg (kJ mol−1) Ezg (kJ mol−1)
infuencing the success of these processes. Acetaminophen −6 −12
Te strength and stability of the adsorption are de- Atrazine 43 2
termined by the zeolite-guest interaction energy, which is Cafeine 12 5
a critical factor in the adsorption process (Table 1). A Carbamazepine 40 11
stronger afnity between the zeolite and the organic pol- DEET 97 6
lutant is indicated by a higher interaction energy, which Diazepam 17 5
leads to a more efcient adsorption process. Te zeolite’s Diclofenac −15 −2
chemical makeup, the framework’s structure, the size and Dilantin 14 1
shape of the organic pollutant molecule, and the type of Estrone 100 35
chemical bonds that connect the zeolite and the pollutant Fluoxetine 100 98
Gemfbrozil 98 6
can all have an impact on the interaction energy. Designing
Hydrocodone 23 26
efective adsorbents for the removal of organic contaminants Ibuprofen 98 6
and optimizing adsorption procedures depend on an un- Meprobamate 97 7
derstanding of the zeolite-guest interaction energy [21]. Naproxen 82 2
Researchers can increase the adsorption capacity and se- Oxybenzone 99 47
lectivity of zeolites by adjusting the zeolite characteristics Pentoxifylline 21 3
and managing the interaction energy, leading to more ef- Sulfamethoxazole 13 0
fective methods for water purifcation and environmental TCEP 21 7
remediation. Triclosan 99 45
Te sorption and desorption processes are pivotal as- Trimethoprim 46 5
pects of the studied zeolite-based organochlorine adsorption
system. Te investigation into sorption involves un- implementation of zeolite-based remediation strategies.
derstanding how zeolites attract and retain organochlorine Table 1 shows a summary of various organic pollutants with
compounds, shedding light on the mechanisms governing their associated interaction energy as reported in the liter-
the initial uptake. Te zeolite’s microporous structure and ature for two zeolite types, MOR and FAU.
high surface area play a crucial role, providing ample sites for Molecular mechanical techniques enable atomic-level
the adsorption of organochlorines. Te simulations and analysis of materials and determination of bulk properties
experimental data collectively suggest that the adsorption using force felds [22]. Te interactions in a molecular system
process is infuenced by factors such as the type of zeolite, are determined by potential functions. A specifc instance of an
structural modifcations, and the specifc characteristics of adsorption process is represented by each model statement in
the organochlorine compounds. On the other hand, the statistical physics modeling. Tese models are developed
exploration of desorption is equally signifcant, as it ad- progressively, starting with the simplest Langmuir model (one
dresses the practical applicability and sustainability of site energy, one site occupation, and one layer) and progressing
zeolite-based remediation. Te ability to desorb and release to more complicated adsorption systems such as the Hill model
organochlorines from zeolites is vital for the regeneration (monolayer, monoenergy, and multioccupancy), double Hill
and reusability of the adsorbent material. Te study’s em- (monolayer, diferent site energies, and mono-occupancy), or
phasis on investigating desorption characteristics provides BET (multilayer, multisite energies, mono-, or multi-
valuable insights into the reversibility of the adsorption occupancy) [23, 24]. In order to choose the optimal material for
process. Factors afecting desorption efciency, such as adsorption of organochlorines (Figure 1) from water, this
temperature and pore structure, are crucial considerations research ofers a simulation-based method to forecast the
for optimizing zeolite regeneration. Te comprehensive loading ability of several zeolites. In order to determine the
nature of the study, incorporating both simulation and signifcance of the many descriptors that have an impact on the
experimental results, allows for a nuanced understanding of elimination of these organic pollutants from water, a data-
sorption and desorption dynamics. Te structural modif- based approach was used. Insights into the sorption mecha-
cations, including silicon, aluminum, or sodium doping, nisms were provided by quantum chemical simulations. Te
contribute to tailoring zeolite properties for enhanced impact of doping the zeolites with particular elements on the
sorption and desorption capabilities. Tis strategic approach loading capacity was then thoroughly examined. Laboratory
aligns with the goal of developing efcient and sustainable experiments were performed using adsorption columns, and
environmental remediation technologies. In essence, the the role of interstitial water was established.
discussion on sorption and desorption in this study high- Tis comprehensive study on the adsorption capacities
lights the multifaceted nature of zeolite-organochlorine of various zeolites for selected organochlorines, Figure 1,
interactions. It underscores the importance of not only ef- holds signifcant importance in the feld of environmental
fective adsorption but also the ability to release contami- science and material engineering. Te research not only
nants, providing a holistic perspective on the practical provides detailed insights into the structural characteristics
4 Journal of Chemistry
Figure 1: Examples of chlorinated organic pollutants that have been reported in surface and underground waters.
and performance of diferent zeolites but also establishes trials, operations of translation, rotation, regeneration,
a correlation between their properties and adsorption ca- CFCMC interchange of “fractional” particles, and identities
pacities. Te identifcation of the universal force feld as the adjustments. Whole adsorption isotherms were calculated
most suitable for predicting adsorption behaviors enhances for each chlorinated compound. Material Studio software
the understanding of molecular interactions. Moreover, the was used for all simulations [32]. Te IRASPA software [33]
investigation into the impact of zeolite doping on adsorption was used to quantify the characteristic features of zeolites,
capacity, supported by thorough analysis and charge dis- including mass (M), density (D), helium void fraction
tribution studies, adds a valuable dimension to the study. (HVF), accessible pore volume (APV), gravimetric surface
Te laboratory experiments, particularly the comparison area (GSA), and dimensionality of the pore system (Di). To
between treated and untreated zeolites, reveal practical investigate the efect of silicon, aluminum, or sodium doping
implications for zeolite application in pollutant removal. Te in the CLO zeolite, fve phosphorus atoms were replaced
fndings emphasize the critical role of oxygen atoms in with silicon, aluminum, or sodium. In each case, an extra
creating active adsorption sites, shedding light on a key framework was created in order to balance the charges. After
structural feature infuencing adsorption efciency. Gen- this, new charges were assigned on all the atoms both in the
erally, this research contributes to the optimization of zeolite zeolite structure and the extra framework created.
selection and treatment methods for environmental re- Clinoptilolite and faujasite zeolites were obtained
mediation, ofering valuable insights for designing efective commercially for experimental use. Before sorption exper-
adsorbents in pollution control applications. iments, the zeolites were activated by heating at high
temperatures of 400°C under vacuum pressure for 24 hrs.
2. Materials and Methods Te experimental setup was meticulously designed,
employing columns as the primary vessels for adsorption.
Te International Zeolites Association (IZA) database Tese columns were flled with various adsorbents, which
provided all of the zeolites used in this simulation work [25]. were compacted using distilled water [34] to ensure con-
Te structures for the organochlorines were obtained from sistency in the packing density. Two diferent types of zeolite
PubChem [26] and optimized using the Dmol3 module adsorbent materials were used, faujasite zeolite (zeolite x)
program [27] in material studio software [28] using the and clinoptilolite, for comparison purposes. To enhance the
GGA/PBE functional [29] and DND basis set [30]. Grand- hydrophobic nature and Si/Al ratio of zeolite x, a mixture
canonical Monte Carlo (GCMC) simulations were carried with silica gel in a specifc ratio (5 : 3) was prepared. Tis
out utilizing the continuous fractional component approach careful manipulation aimed to optimize the adsorption ef-
(CFCMC) to explore adsorption [31]. Mostly in the grand- fciency of zeolite x when utilized in the columns. Tis was
canonical ensemble, the chemical potential, volume, and done for both treated and untreated zeolite. Simultaneously,
temperature remain constant. Te frugality coefcients were clinoptilolite, known for its Si/Al ratio ranging between 4.0
taken to be 1. Every location on the adsorption isotherms and 5.3 [35], was investigated independently. Silica gel,
was created using production cycles 5 × 105. Each Monte a material with potential adsorption capabilities for organic
Carlo cycle was performed using the random Monte Carlo pollutants, was also included in the experimental design for
Journal of Chemistry 5
comparison purposes. Te experimental apparatus com- organometallic systems and other systems for which other
prised fve columns: one contains silica gel, another with force felds do not have parameters [40]. As a result, we
treated zeolite x, a third with untreated zeolite x, the fourth concentrated on the universal force feld in the primary text.
with treated clinoptilolite, and the ffth with untreated cli- Table S1 (supporting information) displays the structural
noptilolite zeolite. A meticulous calibration process was characteristics of 245 diferent types of zeolites along with
implemented using various standard solutions derived from the loading of HCB, including mass (M), density (D), helium
a 30 ppm stock solution of the hexachlorobenzene pollutant. void fraction (HVF), accessible pore volume (APV),
Each column was flled with a predetermined amount of gravimetric surface area (GSA), and dimensionality of the
water to establish the initial pollutant concentration. Te pore system (Di). Out of the 245 zeolites that were studied,
pollutant solution, with a known concentration of 15 ppm, 71 zeolites demonstrated the ability to adsorb hexa-
was systematically passed through each column. Filtrate chlorobenzene with CLO, a cubic microporous gallophos-
concentrations were confrmed through multiple iterations, phate with three-dimensional channels surrounded by 20-
providing a comprehensive dataset for determining equi- membered rings in the shape of a four-leaf clover, having the
librium concentrations and evaluating the reusability of each highest capacity of 65.84 wt%. CLO showed a stepwise
zeolite under diferent conditions. multilayer adsorption isotherm that only occurs when the
sample surface has various adsorption site types with en-
3. Results and Discussion ergetically diverse properties (Figure 4).
Conversely, the loading of all the other zeolites was
3.1. Simulations of Organochlorine Adsorption Mechanisms to relatively low with nearly a horizontal plateau isotherm. Tis
Zeolites. Te structures of selected organochlorines (Fig- is attributed to the limited active sites for HCB adsorption.
ure 1) were optimized and their Frontier molecular orbitals CLO has a mass of 17728 g/mol, density of 1129 kg/m3,
were calculated (Figure 2). Te global hardness, calculated helium void fraction of 0.55, a dimensionality of the pore
using the HOMO energy of the orbital (EHOMO ), is fre- system of 3, a largest overall cavity diameter of 15.33, and
quently connected with a molecule’s capability of donating a larger volumetric surface area of 1070 m2/cm3. On the
electrons, while the LUMO energy of the orbital (ELUMO ), other hand, RWY has a lower mass and density (2884 g/mol
depends on the molecule’s capability to accept electrons. Te and 1129 kg/m3, respectively) and a higher helium void
diference in energy between HOMO and LUMO also fraction, accessible pore volume, and gravimetric surface
contributes to the explanation of the electrical and chemical area (0.66 and 0.77 cm3/g, respectively) compared to CLO. It
properties of the particles. Hexachlorodecane possesses the has a comparable volumetric surface area with respect to
smallest HOMO energy (−0.2558 eV), but dichlor- CLO. It has 54.23 wt% of hexachlorobenzene and is ranked
odiphenyltrichloroethane has a considerably greater HOMO the second after CLO in terms of wt%. ITV is ranked the
energy (−0.1832 eV), as determined by analyses of the third after RWY. Despite ITV having a larger mass, density
electron orbitals energy levels. Te HOMO-LUMO energy and VSA than RWY, it possesses a lower HVF, APV, GSA,
gaps are 0.1456, 0.1952, 0.1841, 0.1777, 0.1912, 0.1761, 0.1406, and Di compared to RWY. Tese diferences may be used to
0.1696, and 0.1977 eV for hexachlorobenzene, hexa- explain why it has 33.65 wt% compared to 54.23 wt% of
chlorotetradecane, hexachlorodecane hexachloro RWY. ITV is followed closely by IFU which has adsorption
cyclohexane, heptachlorodecane, octachlorodecane, capacity of 32.41% wt. IFU has a small mass, HVF, APV,
dichlorodiphenyltrichloroethane, heptachlorotetradecane, VSA, and GSA compared to ITV. Tis may be the reason
and octachlorotetradecane, respectively. Terefore, octa- why it has a lower capacity than ITV. However, IFU has
chlorotetradecane has a higher stability. It has relatively a larger Di than ITV and this explains why it ranked after
a higher chemical hardness of 0.1977 eV, while dichlor- ITV. LTA has almost similar loading capacity as IFU
odiphenyltrichloroethane is less stable and has a relatively (32.31 wt% vs. 32.41 wt%, respectively). Tis closeness may
low chemical harness of 0.1406 eV. be attributed the closeness in VSA and Di between the two
Considering the structure of HCB (a phenyl ring with six zeolites. However, LTA has a lower mass than HVF, APV,
chlorine atoms attached on each of the carbon atoms), we and GSA and this explains why it has a lower capacity for
frst evaluated its adsorption capacities on the zeolites. HCB than IFU. On the other hand, zeolite LTA has a higher
Materials Studio software has a variety of inbuilt force felds density than IFU and this may be an attribute as to why it
which are designed for the adsorption of diferent categories ranked ffth after IFU. Zeolite IRR has a higher mass than
of compounds. Tus, we frst compared the outcomes of HVF, APV, GSA, and Di. It is expected that it has a higher
several force felds (Figure 3), in order to choose a force feld capacity than LTA. However, this is not the case; instead, it
that is appropriate for our systems. has a lower wt% capacity than LTA. It has 31.31 wt%
In contrast to compass [36], valency-consistent force compared to 32.31 wt% of LTA. Tis scenario might have
feld, CVFF [37], and polymer-consistent force feld, PCFF occurred due to the fact that zeolite IRR has a lower density
[38], we discovered that the universal force feld [39] was and VSA compared to LTA. Zeolites FAU, SBT, and VFI are
best for organic compounds used in this investigation. Since among those with similar HCB adsorption capacity, with
the geometries and conformational energy diferences of 28.32 wt%, and it may be attributed to their negligible dif-
organic molecules, main-group inorganics, and metal ferences in HVF, APV, and Di. In addition, they have almost
complexes have been found to be predicted with moderate equal densities and VSA. On the other hand, zeolite TSC has
accuracy by the universal force feld, it is suggested for the largest mass, density, and Di and yet it has a lower wt%
6 Journal of Chemistry
-0.2341 -0.2558
HOMO -0.2294 ΔE = 0.1952 ev
ΔE = 0.1841 ev
ΔE = 0.1456 ev -0.0389 -0.0717
-0.0838
LUMO
HOMO -0.2495
-0.2553 -0.2427 ΔE = 0.1761 ev
ΔE = 0.1777 ev
ΔE = 0.1912 ev
-0.0776 -0.0734
LUMO -0.0515
(d) (e) (f )
HOMO
-0. 1832
-0.2388 -0.2395
ΔE = 0.1406 ev
ΔE = 0.1696 ev ΔE = 0.1977 ev
-0.0777
LUMO -0.0692 -0.0418
Figure 2: Optimized structures of hexachlorobenzene (a), hexachlorotetradecane (b), hexachlorodecane (c) hexachlorocyclohexane
(d), heptachlorodecane (e), octachlorodecane (f ), dichlorodiphenyltrichloroethane (g), heptachlorotetradecane (h), and octa-
chlorotetradecane (i), respectively. Te colors gray, white, and green, respectively, represent the atoms of carbon, hydrogen, and chlorine.
Positive and negative orbital values are represented by the yellow and blue sections, respectively.
HCB capacity. In comparison with CLO, TSC would be adsorption capacity less than 1 wt%. All the other zeolites
expected to have a higher percentage capacity considering among the 71 had capacities above 5 wt%. A correlation
mass, density, and Di. However, the opposite is seen for HCB analysis was done in order to get more understanding of how
loading. Tis may be attributed to low HVF, APV GSA, and various structural features afect the wt% (Figure 5). Te
VSA in TSC zeolite. Among the 71 zeolites that showed the results showed that the mass of the zeolites and density were
possibility of adsorbing HCB compounds, eight of them had negatively correlated to adsorption capacity. Tis means that
an adsorption capacity less than 5 wt% with two having an a heavier material would have less uptake capacity compared
Journal of Chemistry 7
17
120
16
100
15
Loading of HCB on FAU
13 60
12
40
11
20
10
9 0
200 400 600 800 1000 200 400 600 800 1000
Fugacity Fugacity
Figure 3: Loading of hexachlorobenzene (HCB) molecules on FAU and CLO zeolites using various force felds.
60 60
40 40
30 30
20 20
10 10
0 0
20000 40000 60000 80000 800 1000 1200 1400 1600 1800
Mass (g/mol ) Density (Kg/m3)
(a) (b)
60 60
50 50
40 40
30 30
20 20
10 10
0
0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Helium void fraction (HVF) Accessible pore volume (APV)
(c) (d)
60 60
Adsorption capacity (wt%)
Loading capacity (wt%)
50 50
40 40
30 30
20 20
10 10
800 1000 1200 1400 1600 1800 2000 2200 500 1000 1500 2000 2500
Volumetric surface area (VSA) Gravimetric surface area (GSA)
(e) (f )
60 60
Adsorption capacity (wt%)
50 50
40 40
30 30
20
20
10
10
0
0 2 0 2 4 6 8 10 12 14 16
Dimensionality of pore system (DPS) Overall largest cavity diameter (Di)
(g) (h)
Figure 5: Correlation between zeolites structural parameters; M, D, HVF, APV, VSA, GSA DPS, Di, and adsorption capacity (wt%).
Journal of Chemistry 9
O +
P
Ga
ΔE = 0.02 Ha
0.00
ΔE = -0.04 Ha
Cl
Figure 6: Structures of optimized CLO cluster and adsorbed HCB molecules on the CLO cluster.
doping CLO zeolite provides improved results in terms of (Table S2). Te best zeolite for adsorption of hexa-
both loading and wt%, whereby aluminum provided the best chlorodecane is CLO which has 47 molecules per cell an
results (Table 3). Te energy diference between the reactants equivalent of 46 wt%. Tis is followed by TSC, FAU, ITV,
adsorbed states upon doping was also investigated (Fig- IFT, RWY, IRY, SBT, and SBS, respectively. Finally, CLO is
ure 7). Tis was done by substituting phosphorus atom with still the best adsorbent for heptachlorodecane adsorption
a single aluminum atom in the CLO cluster. Contrary to the with 30 molecules per cell. Te adsorption isotherms rep-
undoped case, it is thermodynamically favorable for HCB to resenting this information are shown (Figure 8). Te shapes
adsorb through dissociative adsorption as opposed to of these fgures can be used to inform us that the zeolites
through hydrogen bonding. It is also worth noting that studied have diferent distribution of active sites for ad-
doping increases the system energy, i.e., 0.02 verses −12.25 sorption of the selected organic compounds. However, it is
and −0.04 verses −12.63 Ha for the undoped and doped evident that CLO zeolite has more active sites as seen from
zeolites, respectively. its adsorption isotherm. It is important to note, from these
Te efect of Al doping on Muliken charge distribution isotherms, that there are diferent sites on each zeolite with
was investigated (Table 4). Te charge on each atom before diferent energy requirement leading us to stepwise iso-
and after adsorption was calculated as 0.060 on each carbon therm. Tis applies to all zeolites studied (See Figure 8). Tis
atom and −0.060 on each chlorine atom before adsorption, fnding led us into investigation to fnd out about these active
while −0.001, −0.001, 0.000, −0.009, −0.036, and −0.002 for sites. Adsorption locator module was thus used to identify
carbon atoms and 0.072, 0.031, 0.033, 0.028, 0.024, and 0.028 the adsorption sites as seen in Figure 9. It is noted that the
for chlorine atoms after adsorption. While HCB has a zero organic compounds are adsorbed through electrostatic in-
charge before adsorption, after adsorption, HCB molecule teraction within the zeolite structure. However, in this case,
inside the zeolite becomes positively charged with a net oxygen atoms provide active interaction sites. It is indicated
charge of 0.167. On the other hand, a doped CLO cluster is that the pores surrounded by more oxygen atoms are the
positively charged with a charge of 0.001 before adsorption. most preferred. Tus, afnity for organic pollutants can be
After adsorption, it acquires a −0.096 net charge, making it increased by increasing the number of oxygen atoms.
negatively charged. According to these Mulliken charge Te performance of CLO zeolite over other zeolite
values, charge transfer from HCB to the CLO zeolite takes materials enabled an investigation to fnd out possible
place, similar to the undoped scenario, and thus it is noted features responsible for its best performance. Previously, it
that doping does not alter the charge transfer pattern. has been reported that various general structural features
For comparison, we also looked into the adsorption of have an impact on sorption ability of zeolites [43-44].
hexachlorocyclohexane (Table S2), dichlorodiphenyltri- However, few reports have focused on how sorption per-
chloroethane (Table S3), hexachlorotetradecane (Table S4), formance is infuenced by the structural composition of
heptachlorodecane (Table S5), octachlorodecane (Table S6), zeolites [44]. We therefore investigated the distribution of
hexachlorodecane (Table S7), heptachlorodecane (Table S8), active adsorption sites on the surfaces of the zeolites which
and octachlorobenzene (Table S9). It is observed that for facilitated the afnity for organic compounds. Since most of
hexachlorotetradecane, CLO zeolite was also the best with the zeolites have oxygen atoms in their framework, it was
a loading capacity of 21 molecules per cell. Tis is equivalent very necessary to fnd out the infuence of these atoms to
to 31 wt%. Tis is followed by ITV, TSC, FAU, IFT, IRY, adsorption capacity. Interestingly, it was discovered that
RWY, SBT, and SYT, respectively. For hepta- oxygen atoms played a key role in active sites provision. Tis
chlorotetradecane, CLO is still the best adsorbent with 22 implies that the more the oxygen atoms, the more the active
molecules adsorbed per cell, which corresponds to 36 wt%. sites on the zeolite surface, and hence more adsorption
Tis is followed by ITV, FAU, IRY, SBT, IFT, RWY, SBS, and capacity. Table 5 shows a summary of the adsorption ca-
SYT. Octachlorotetradecane also exhibits CLO as the ad- pacity versus the number of oxygen atoms.
sorbent with 17 molecules per cell. Tis is equivalent to 38 wt
%. It is followed by TSC, ITV, FAU, RWY, IRY, IFT, SBT,
DFO, and EMT, respectively. Dichlorodiphenyltrichloro- 3.2. Laboratory Adsorption Experiments of Organochlorines
ethane still exhibits CLO as the best candidate for its ad- on Zeolites. Zeolites treatment is one of the key practices
sorption with an uptake capacity of 43 molecules per cell. performed in order to activate them prior to any application
After CLO, TSC, ITV, FAU, RWY, IRY, IFT, SBT, DFO, and [54]. One of the reasons why zeolites are treated is to remove
EMT show good adsorption capacities. On the other hand, interstitial water from their framework. Tere are several
hexachlorocyclohexane adsorbs on CLO with 113 molecules zeolites treatment methods that have been employed in
per cell which is equivalent to 54 wt%, followed by TSC, ITV, previous studies [54-56]. Tese include heat treatment at
IFT, FAU, IRY, RWY, DFO, SYT, and LTA, respectively high temperatures and acid treatment methods among
Journal of Chemistry 11
O
Al
Ga P
0.00
E = -12.2513 Ha
Cl E = -12.6285 Ha
Figure 7: Optimized structures of Al-doped CLO cluster and HCB adsorbed on the doped cluster.
20 20
Loading per cell
0 0
200 400 600 800 1000 200 400 600 800 1000
Fugacity (KPa) Fugacity (KPa)
20
Loading (per cell) 40
Loading per cell
10
20
0
0
200 400 600 800 1000 200 400 600 800 1000
Fugacity (KPa) Fugacity (KPa)
100 40
Loading (per cell)
Loading (per cell)
80
60
20
40
20
0
0
200 400 600 800 1000
200 400 600 800 1000
Fugacity (KPa)
Fugacity (KPa)
CLO FAU DFO CLO ITV IRY
TSC IRY SYT TSC IFT SBT
ITV RWY LTA FAU RWY SBS
IFT
(e) (f )
Figure 8: Continued.
Journal of Chemistry 13
30
20
25
Loading (per cell)
15 10
10
5
0
0
200 400 600 800 1000 200 400 600 800 1000
Fugacity (KPa) Fugacity (KPa)
(a) (b)
Figure 9: (a) CLO zeolite framework. (b) Adsorption sites for hexachlorobenzene within the CLO zeolite.
others. In the current study, the efect of heat treatment at obtained from the column for analysis. On the other hand,
400°C to adsorption of chlorinated compounds was the sorption results for silica gel were not pronounced. Tis
investigated. therefore showed that silica gel could not adsorb organic
Since columns were used for adsorption experiments, pollutants but increased hydrophilic nature of zeolites
packing was done using distilled water by frst preparing (Table 7). Tis result closely agreed with the present liter-
slurry to facilitate compact packing. Te results showed that ature on the use of silica gel in adsorption [57].
the selected zeolites could not make slurry since they were Basing on the results obtained, it was discovered that the
highly hydrophobic. Tis hydrophobicity is attributed to untreated zeolite performed better than the treated one.
a low Si/Al ratio. Interestingly, it was discovered that the Si/ From Table 7 above, the treated zeolite adsorbed 3.2 ppm,
Al ratio could be improved by adding silica gel powder to the 0.3 ppm, 0.1 ppm, 0 ppm, and 0.2 ppm for the fve time runs,
zeolite in the ratio 3 : 5 (silica: zeolites). Tis was discovered respectively. Tis trend indicates that the treated zeolite
after several trials as shown in Table 6. could reach saturation point only after the third run, after
Te sorption results showed that clinoptilolite zeolite which no more pollutant could be adsorbed. Tis means that
could not be used in columns for adsorption since it caused if one were to re-use the adsorbent, then it could be used
logging in the columns. As such no solution could be three times before recovery is done. On the other hand,
14
Table 6: Determination of appropriate zeolite silica mixing ratio in Table 8: Te remaining hexachlorobenzene compound concen-
order to improve Si/Al ratio in the zeolite. tration after adsorption.
Amount Concentration
Amount of silica (g) Comment Material
of zeolite (g) 1 2 3 4 5
0.0 Hydrophobic Treated zeolite x 2.7149 2.3756 0.2262 0.2262 0.1131
0.1 Hydrophobic Untreated zeolite x 1.4706 0.9615 0.2262 0 ND
0.2 Hydrophobic Silica gel 4.4118 4.4118 4.2986 4.4118 4.2986
0.5
0.3 Hydrophilic
0.4 Hydrophilic
0.5 Hydrophilic zeolites for organic pollutant adsorption heavily depended
on the number of oxygen atoms present in the structure
(Table 5). Integrating these three simulation fndings with
Table 7: Hexachlorobenzene removal by faujasite zeolite and silica the experimental results, we found out that, heating the
gel for fve runs into columns. zeolite, apart from removing inorganic impurities, also
Amount removed (ppm) removes interstitial water from the zeolite framework. Tis
Material interstitial water within the zeolite structure contains the
1 2 3 4 5
hydroxyl radicals which are reactive with organic molecules.
Treated faujasite zeolite 3.2851 0.3393 0.149 ND 0.2149 Tey contain a negative charge which reacts with the positive
Untreated faujasite
5.3294 0.5091 0.7353 0.2262 ND charges in the compound (ion exchange). Tus, heating the
zeolite
Silica gel 0.5882 ND 0.1132 ND ND zeolite material reduces the value of negative hydroxyl ions
and therefore the reaction between the zeolite and the
ND � not detected.
positive organic compound reduces. Tis causes poor per-
formance of the heated zeolite and the better performance of
untreated zeolite adsorbed 5.3 ppm, 0.5 ppm, 0.7 ppm, and the untreated zeolite material. Secondly, since the in-
0.2 ppm for the fve run times, respectively. Tis trend is teraction mechanism involved between the zeolite and the
quite interesting in the sense that the material could not hexachlorobenzene is electrostatic through hydrogen
reach saturation after the fve re-uses. Tis is due to the fact bonding, heating reduces the number of hydrogen atoms in
that the remaining concentrations after the third run could the zeolite framework which could facilitate the interaction.
not be detected, implying that most the compounds had Hence, the untreated zeolite with more hydrogen atoms
been adsorbed and the material could still adsorb more. performs better. Finally, the simulation results indicated that
Table 8 shows the remaining concentrations after every run the active adsorption sites for organic pollutants under study
through the columns. In addition, the amount of compound were provided by the oxygen atoms. Since water molecules
adsorbed after every run through the columns was much compose of oxygen atoms, heating would mean removing
higher than that of treated zeolite. Previous researchers have more oxygen atoms from the framework. Terefore, there
explained the efect of various treatment methods to ad- would be less oxygen atoms which could bond with the
sorption of methylene blue. Hor and co-authors [58] re- chlorine atoms in the organochloride compounds (Figures 6
ported that both the physico-chemical methods of acid and 7). Tus, the treated zeolite adsorbs less than the un-
treatment and base treatment were able to boost the per- treated one. It is worth noting that the present fndings, both
formance of natural zeolite from 41% to 98.8%, and 52.2%, experimental and computational, have strongly agreed as
respectively. However, their report showed that combining discussed.
reduced performance would be observed when any of the Comparing our study with the previous study by Kul-
methods were combined. Tis shows that diferent treatment deyev et al. [59] on the adsorption of heavy metals on ze-
methods favor diferent adsorbates. Our present work olites, it noted that our fnding on the zeolite treatment
therefore reveals that thermal treatment as a method of methods contradicts what was suggested in their study. Te
zeolite treatment lowers the zeolite performance, as com- study focused on enhancing the adsorption capabilities of
pared to untreated zeolites. natural zeolites from the Shankanay deposit in Kazakhstan
Te percentage removal calculation reveals that un- for heavy metal ions in water. Te researchers employed
treated zeolite removes up to 91% of hexachlorobenzene. thermal methods to modify the zeolites and increase their
Comparing this experimental work with the simulated re- specifc surface area and porosity. Tree stages of activation
sults, we were able to fnd out that for organic pollutants, were conducted, involving cleaning the zeolites with water
untreated zeolites performs better than treated zeolite. Tis vapor, demineralization through acid treatment, and high
is based on the fact that, one, our computations of the thermal activation in an inert atmosphere. Te results
Mulliken charge analysis performed (Table 2) revealed that showed that thermal activation signifcantly improved the
there is charge transfer between the zeolite and organic adsorption capacity of the zeolites for various heavy metals.
pollutants during adsorption. Secondly, the interaction Te optimal condition was found to be a thermal treatment
mechanisms between the zeolite and hexachlorobenzene at 550°C for approximately 2 hours. Te heat-treated zeolites
involve electrostatic interaction through hydrogen bonding demonstrated high adsorption capacities, with percentages
(Figure 6), and thirdly, the number of active sites in the ranging from 83% to 100% for copper, cadmium, lead, and
16 Journal of Chemistry
nickel. Te adsorption process was found to be rapid, with helps us create more efcient methods for removing organ-
intense adsorption occurring within the frst 10 minutes. Te ochlorine pollutants from contaminated sites. Despite its
residual concentration of heavy metals in the fltered water value, the present study has a number of shortcomings that
remained within the permissible limits for drinking pur- need to be acknowledged. First of, the study was carried out
poses. Te study concluded that natural zeolites from the in a controlled laboratory setting, which might not accurately
Shankanay deposit have high potential as cost-efective refect the complexity of real-world settings. In addition, the
sorbents for the treatment of heavy-metal-polluted water, investigation of the kinetics of the process, which is essential
particularly in wastewater treatment and water softening for developing efective remediation solutions, was not done
applications in the metallurgical industry. Te fndings also as much in the study as it was largely on equilibrium ad-
suggested that thermal treatment may be more efcient than sorption. Te analysis of individual organochlorine com-
chemical treatment in enhancing the adsorption capacity of pounds also fails to take into account the occurrence of
zeolites, but further research is needed to confrm this. combinations in real-world situations, where many pollutants
However, this contradiction may be as a result of diferent may interact and compete for adsorption sites. Te re-
interaction mechanisms between the zeolite and the strictions on mass transfer, reactor design, and the behavior of
material. zeolite materials at greater scales may all provide difculties
for scaling up the fndings to broader uses. In addition, not
much study has been done on the long-term stability and
3.3. Factors Tat Afect Adsorption of Organochlorine Pol-
durability of zeolite materials, which is crucial for assessing the
lutants on Zeolites. Te outcomes of the current study may
viability and sustainability of zeolite-based remediation so-
be infuenced by a number of factors. Te zeolite’s own
lutions. By addressing these restrictions, future research will
characteristics, such as its surface area, pore size, surface
improve our knowledge and speed the creation of more ef-
charge, and composition, are important. Te adsorption
fcient and useful environmental cleanup technologies.
behavior can be afected by changes in the solution’s pH and
temperature, as well as the presence of other chemicals or
pollutants. Te molecular size and hydrophobicity of the 4. Conclusions
organochlorine compounds, as well as other physico- In this work, the investigation into organochlorine ad-
chemical characteristics, afect their adsorption afnity. Te sorption using zeolites, combining simulation and experi-
adsorption kinetics and achievement of equilibrium are mental approaches, has delivered impactful fndings.
infuenced by the contact time and agitation level. Important Simulation work utilizing the IZA database and PubChem,
factors to take into account also include the presence of rival along with GCMC simulations, deepened our understanding
chemicals and the capacity to regenerate and re-use zeolite of adsorption mechanisms. Material Studio and IRASPA
materials. For precise and dependable fndings in research software facilitated quantitative zeolite analysis. Experi-
involving the adsorption of organochlorine chemicals uti- mental validation with clinoptilolite and faujasite zeolites,
lizing zeolite materials, it is important to comprehend and subjected to strategic treatments, aligned with simulation
regulate these aspects. results, emphasizing the study’s robustness. Structural
modifcations, including silicon, aluminum, or sodium
3.4. Signifcance of Tis Study. Tis study has signifcant doping, and variations in the Si/Al ratio revealed their
consequences for environmental cleanup initiatives. It clarifes impact on adsorption efciency. Te inclusion of diferent
the adsorption mechanisms at work and shows the critical role zeolite types and silica gels for comparison enhanced
played by electrostatic interactions, particularly hydrogen practical relevance. Rigorous column adsorption tests and
bonding. Researchers can create more efective adsorbents for calibration curves provided quantitative insights. Tis study
removing organochlorines from contaminated settings based efectively bridges simulation and experimentation, high-
on this understanding. In addition, in order to help with the lighting zeolites’ adaptability in pollutant remediation. Te
selection of zeolites for particular remediation projects, the results afrm zeolites’ potential as versatile and sustainable
study compares the adsorption capacities of various zeolite materials for environmental applications, advancing the
materials and identifes certain zeolites with high adsorption understanding of their efcacy in water purifcation and
capacities for particular pollutants. Understanding how zeolite environmental protection.
structural characteristics afect adsorption performance might
help in the design and optimization of zeolite materials with Data Availability
improved qualities. Te study examines the efects of doping
zeolites with various elements and fnds that doping, especially Te research data used to support the fndings of this study
with aluminum doping, can increase loading capacity. Lastly, are included within the article and supporting information
the experiments show how crucial it is to take into account thereof attached.
how zeolite materials are treated, with untreated zeolite
exhibiting greater adsorption capacity and reusability. Col- Conflicts of Interest
lectively, these efects help us better understand how organ-
ochlorine adsorption operates, direct the design and selection Te authors declare that they have no conficts of interest
of zeolites, and ofer helpful tips for enhancing adsorption that could have appeared to infuence the work reported in
procedures in environmental remediation. Ultimately, this this paper.
Journal of Chemistry 17
simulations,” in Forcefelds for Atomistic-Scale Simulations: zeolites,” Journal of Molecular Structure, vol. 957, no. 1–3,
Materials and Applications, pp. 299–313, Springer, Berlin, pp. 47–54, 2010.
Germany, 2022. [43] N. M. Mahmoodi and M. H. Safar-Dastgerdi, “Zeolite
[29] L. Song, W. Liu, F. Xin, and Y. Li, “Materials studio simulation nanoparticle as a superior adsorbent with high capacity:
study of the adsorption and polymerization mechanism of synthesis, surface modifcation and pollutant adsorption
sodium silicate on active silica surface at diferent tempera- ability from wastewater,” Microchemical Journal, vol. 145,
tures,” International Journal of Metalcasting, vol. 15, no. 3, pp. 74–83, 2019.
pp. 1091–1098, 2021. [44] F. S. Wanyonyi, T. T. Fidelis, G. K. Mutua, F. Orata, and
[30] S. F. Zaman and K. J. Smith, “A study of synthesis gas A. M. Pembere, “Role of pore chemistry and topology in the
conversion to methane and methanol over a Mo6p3 cluster heavy metal sorption by zeolites: from molecular simulation
using density functional theory,” Molecular Simulation, to machine learning,” Computational Materials Science,
vol. 34, no. 10-15, pp. 1073–1084, 2008. vol. 195, Article ID 110519, 2021.
[31] W. Shi and E. J. Maginn, “Continuous fractional component [45] M. Estermann, L. McCusker, C. Baerlocher, A. Merrouche,
Monte Carlo: an adaptive biasing method for open system and H. Kessler, “A synthetic gallophosphate molecular sieve
atomistic simulations,” Journal of Chemical Teory and with a 20-tetrahedral-atom pore opening,” Nature, vol. 352,
Computation, vol. 3, no. 4, pp. 1451–1463, 2007. no. 6333, pp. 320–323, 1991.
[32] D. S. Biovia, Biovia Pipeline Pilot, Dassault Systèmes, San [46] N. Zheng, X. Bu, B. Wang, and P. Feng, “Microporous and
Diego, BW, USA, 2017. photoluminescent chalcogenide zeolite analogs,” Science,
[33] D. Dubbeldam, S. Calero, and T. J. Vlugt, “Iraspa: gpu- vol. 298, no. 5602, pp. 2366–2369, 2002.
accelerated visualization software for materials scientists,” [47] J. Sun, C. Bonneau, Á. Cantı́n et al., “Te itq-37 mesoporous
Molecular Simulation, vol. 44, no. 8, pp. 653–676, 2018. chiral zeolite,” Nature, vol. 458, no. 7242, pp. 1154–1157, 2009.
[34] F. Gritti and M. F. Wahab, “Understanding the science behind [48] J. Jiang, Y. Yun, X. Zou, J. L. Jorda, and A. Corma, “ITQ-54:
packing high-efciency columns and capillaries: facts, fun- a multi-dimensional extra-large pore zeolite with 20 × 14 × 12-
damentals, challenges, and future directions,” LCGC North ring channels,” Chemical Science, vol. 6, no. 1, pp. 480–485,
America, vol. 36, pp. 82–98, 2018. 2015.
[35] C. Wang, S. Leng, H. Guo et al., “Quantitative arrangement of [49] J. Antúnez-Garcı́a, D. Galván, V. Petranovskii et al., “Te
Si/Al ratio of natural zeolite using acid treatment,” Applied
efect of chemical composition on the properties of LTA
Surface Science, vol. 498, Article ID 143874, 2019.
zeolite: a theoretical study,” Computational Materials Science,
[36] H. Sun, “Compass: an ab initio force-feld optimized for
vol. 196, Article ID 110557, 2021.
condensed-phase applications overview with details on alkane
[50] F. Sifuna Wanyonyi, T. T. Fidelis, H. Louis et al., “Simulation
and benzene compounds,” Te Journal of Physical Chemistry
guided prediction of zeolites for the sorption of selected
B, vol. 102, no. 38, pp. 7338–7364, 1998.
anions from water: machine learning predictors for enhanced
[37] I. Yungerman, I. Starodumov, A. Fulati, K. Uto, M. Ebara, and
loading,” Journal of Molecular Liquids, vol. 355, Article ID
Y. J. T. J. O. P. C. B. Moskovitz, “Full-atomistic optimized
potentials for liquid simulations and polymer consistent force 118913, 2022.
feld models for biocompatible shape-memory poly (Ε-Cap- [51] D. T. Bregante, J. Z. Tan, A. Sutrisno, and
rolactone),” Te Journal of Physical Chemistry B, vol. 126, D. W. J. C. S. Flaherty, “Heteroatom substituted zeolite fau
no. 21, pp. 3961–3972, 2022. with ultralow Al contents for liquid-phase oxidation cataly-
[38] J. Blomqvist, L. Ahjopalo, B. Mannfors, and L.-O. Pietilä, sis,” Catalysis Science and Technology, vol. 10, no. 3,
“Studies on aliphatic polyesters I: ab initio, density functional pp. 635–647, 2020.
and force feld studies of esters with one carboxyl group,” [52] Y. Yun, M. Hernández, W. Wan et al., “Te frst zeolite with
Journal of Molecular Structure, vol. 488, no. 1-3, pp. 247–262, a tri-directional extra-large 14-ring pore system derived using
1999. a phosphonium-based organic molecule,” Chemical Com-
[39] S. L. Mayo, B. D. Olafson, and W. A. Goddard, “Dreiding: munications, vol. 51, no. 36, pp. 7602–7605, 2015.
a generic force feld for molecular simulations,” Journal of [53] Y. Fu, Y. Liu, Z. Li et al., “Insights into adsorption separation
Physical Chemistry, vol. 94, no. 26, pp. 8897–8909, 1990. of N2/O2 mixture on FAU zeolites under plateau special
[40] A. Golchoobi and H. Pahlavanzadeh, “Extra-framework conditions: a molecular simulation study,” Separation and
charge and impurities efect, grand canonical Monte Carlo Purifcation Technology, vol. 251, Article ID 117405, 2020.
and volumetric measurements of Co2/ch4/N2 uptake on nax [54] H. Maghsoudi, “Defects of zeolite membranes: character-
molecular sieve,” Separation Science and Technology, vol. 52, ization, modifcation and post-treatment techniques,” Sepa-
no. 16, pp. 2499–2512, 2017. ration and Purifcation Reviews, vol. 45, no. 3, pp. 169–192,
[41] P. Sherwood, A. H. De Vries, S. J. Collins et al., “Computer 2016.
simulation of zeolite structure and reactivity using embedded [55] R. Mahdavi Far, B. Van der Bruggen, A. Verliefde, and
cluster methods,” Faraday Discussions, vol. 106, pp. 79–92, E. J. J. O. C. T. Cornelissen, “A review of zeolite materials used
1997. in membranes for water purifcation: history, applications,
[42] D. K. Papayannis and A. M. Kosmas, “A computational study challenges and future trends,” Journal of Chemical Technology
of ethyl chloride conversion reactions catalyzed by acidic and Biotechnology, vol. 97, no. 3, pp. 575–596, 2022.
Journal of Chemistry 19