1. Introduction
In recent times, the world has witnessed a rapid increase in the global population, industrialization, agricultural activities, and unplanned urbanization, along with the excessive use of chemicals. Unfortunately, these factors have collectively contributed to a concerning issue of environmental pollution. One significant consequence is the release of colored wastewaters containing dyes from various industries such as textiles, paper, leather, printing, plastic, and food processing [
1].
The presence of dyes in textile wastewater poses a serious environmental challenge due to their inherent characteristics like resistance, high visibility, non-biodegradability, and considerable toxicity. Moreover, there is a potential risk of these dyes transforming into carcinogenic, teratogenic, and mutagenic substances, amplifying the concern even further [
2]. For centuries, dyes have found extensive use in applications such as paint and textiles due to their ease of application, vibrant colors, and water-fastness [
3]. However, their usage comes with detrimental consequences, including adverse health effects on humans, leading to issues like cancer, jaundice, tumors, heart defects, skin irritation, allergies, and mutations.
Furthermore, the discharge of these dyes into environmental water bodies gives rise to grave problems for aquatic and terrestrial ecosystems. The impact includes disrupting photosynthetic processes in aquatic plants, reducing oxygen levels in water, and, in severe cases, causing the suffocation of aquatic fauna and flora [
3]. Addressing these issues has become a critical priority for safeguarding the environment and human health.
Commercial dyes are typically classified according to various factors, such as their color, chemical structure, chemical nature, and application method. Additionally, based on the charge they carry upon dissolution in an aqueous medium, dyes are categorized as cationic, anionic or nonionic. Azo dyes, in particular, pose a significant threat due to the presence of amine groups in the effluent. These dyes are known to be toxic and highly hazardous to both environmental life and human health. Consequently, dye industries have faced mounting pressure to effectively eliminate or minimize the discharge of dyes into water streams before wastewater is released [
4]. Such measures are essential to safeguard the environment and protect human well-being from the potential harmful effects of these dyes.
As a result, it becomes imperative to carefully choose an appropriate treatment method to effectively remove dyes from wastewater and enhance the quality of treated water released into the environment. Consequently, numerous techniques have been explored and applied to tackle the challenge of dye wastewater treatment [
5].
Among the various methods employed, some common approaches include adsorption on activated carbons, coagulation and flocculation, chemical oxidation, reverse osmosis, bacterial action, activated sludge, ozonation, membrane filtration, ion exchange, and electrochemical techniques, all aimed at the removal of dyes from wastewater [
6]. However, these methods have shown certain limitations, such as high operational costs, limited effectiveness, generation of excess sludge, strict environmental requirements, or the possibility of producing more toxic byproducts [
6].
Of all these techniques, adsorption stands out as a highly effective approach for dye removal due to its straightforward design, high efficiency, affordability, production of non-toxic byproducts, rapid adsorption rate, and versatile applicability. Recently, there has been a growing focus on materials based on natural polymers, which offer the advantage of removing pollutants from contaminated water at a reduced cost [
7]. By exploring and utilizing these natural polymer-based materials, progress can be made in addressing the challenges associated with dye wastewater treatment in a more sustainable and efficient manner.
Chitosan has emerged as a promising alternative adsorbent for conventional wastewater treatment processes, garnering significant attention due to its unique properties. Notably, it exhibits a high adsorption capacity, macromolecular structure, abundance, and cost-effectiveness [
8]. As part of the new wave of adsorption techniques, biopolymers like chitin and chitosan have been utilized to effectively remove heavy metal ions and dyes from wastewater [
9].
Chitosan, also known as poly-(1-4)-2-amino-2-deoxy-β-D-glucose, ranks as the second most abundant polysaccharide globally. It is synthesized through the deacetylation of chitin, a substance abundantly found in the shells of shellfish (crabs, lobsters, prawns, crayfish), fungi, insects, and other crustaceans. This feature makes chitosan a sustainable and eco-friendly choice [
8,
10].
The biopolymer, chitosan, possesses several noteworthy characteristics: it is heterogeneous, linear, cationic, and boasts a high molecular weight. Its hydrophobic nature contributes to its exceptional adsorption properties. Additionally, chitosan is non-toxic, biodegradable, hydrophilic, and biocompatible, making it even more desirable for various applications, including wastewater treatment [
8]. These properties collectively position chitosan as a versatile and effective material for addressing pollution challenges while being environmentally friendly.
However, the application of adsorbent materials in wastewater treatment is hindered by the challenge of separating them from the aqueous solution. To address this limitation, magnetic nanoparticles (Fe
3O
4) have been widely adopted in environmental applications due to their unique magnetic properties, enabling a quick and effective solid-liquid separation technique when an external magnetic field is applied during wastewater treatment [
5,
11]. Consequently, chitosan magnetic beads were developed, leveraging the advantages of magnetic nanoparticles in wastewater treatment.
The primary objective of this study was to analyze the adsorption behavior of the Direct Red 83:1 dye on the magnetic adsorbent. To achieve this, experimental data were fitted to various kinetic and isotherm models, providing insights into the adsorption characteristics of magnetic polymers. Additionally, a comprehensive characterization of the adsorbent material was carried out. This research aims to contribute to the understanding and improvement of wastewater treatment processes using magnetic adsorbents.
2. Results
2.1. Effect of Contact Time
The effect of contact time on the uptake of Direct Red 83:1 was investigated across various concentrations of the dye, ranging from 50 to 300 mg/L. The adsorption capacity data at different contact times are presented in
Figure 1. All experiments were conducted under a constant pH of 7.0, utilizing 1g of the polymer, and employing continuous stirring at 500 rpm.
In
Figure 1, it can be observed that the adsorption capacity increased steadily until the point where the adsorption of dye on the polymer reached equilibrium at each concentration. At equilibrium, there exists a dynamic balance between the adsorbed and desorbed dye inside and outside the polymer. The time required to achieve this equilibrium is known as the equilibrium time, and the amount of dye removed by the polymer at this point represents the maximum adsorption capacity [
12].
Within the concentration range analyzed (50–300 mg/L), various adsorption phases can be distinguished. Initially, the adsorption process exhibited a rapid increase during the early stages of contact between Direct Red 83:1 and the adsorbent. For the chitosan magnetic polymer adsorbing Direct Red 83:1 (
Figure 1), equilibrium was rapidly attained at a concentration of 50 mg/L, with only 20 minutes of contact time. At a concentration of 100 mg/L, the adsorption was still fast, although slower than at lower concentrations, reaching equilibrium after 40 minutes of contact time.
As the concentration of the dye increased to 150 mg/L, the adsorption process became asymptotic after approximately 60 minutes of contact time. However, for higher concentrations of the dye (200 and 300 mg/L), the adsorption curves did not display the typical asymptotic form. Instead, the equilibrium time extended with increasing concentrations of Direct Red 83:1. These findings reveal important insights into the dynamics of the adsorption process and the varying behavior of the adsorbent at different dye concentrations.
2.2. Adsorption Kinetics
To determine the equilibration time for achieving maximum adsorption and to understand the rate-determining step of the adsorption process, both time effect and kinetics were investigated for the adsorption of the Direct Red dye on the prepared adsorbents. To accomplish this, three kinetic models were employed in this study: the pseudo first order, pseudo second order, and the intraparticle diffusion model.
The pseudo first order rate of Lagergren [
13] is given by the Equation (1):
This model assumes that the rate of adsorption is directly proportional to the difference between the equilibrium adsorption capacity (qe) and the adsorption capacity at any given time (qt). By fitting the experimental data to this model, the rate constant (k1) and other kinetic parameters can be determined, providing insights into the adsorption mechanism and the rate-limiting step of the process.
The adsorption kinetics can be effectively described by a pseudo-second-order model, as proposed by Ho [
14]. This model is represented by Equation (2):
The pseudo-second-order model assumes that the adsorption rate is influenced by the interaction between the adsorbate (dye) and the adsorbent (the prepared adsorbents). Unlike the pseudo-first-order model, the pseudo-second-order model considers that the rate of adsorption is directly proportional to the square of the equilibrium adsorption capacity (qe), and inversely proportional to the adsorption capacity at any given time (qt).
Fitting the experimental data to this model allows us to determine the rate constant (k2) and further understand the adsorption mechanism, providing valuable information about the kinetics of the adsorption process.
According to the intraparticle diffusion model proposed by Weber and Morris [
15], the root time dependence can be expressed by the following Equation (3):
The intraparticle diffusion model is used to investigate the intraparticle diffusion mechanism during the adsorption process. It suggests that the adsorption occurs in multiple steps, including external surface adsorption, intraparticle diffusion, and equilibrium. The root time dependence in Equation (3) indicates that the rate of intraparticle diffusion is proportional to the square root of the contact time (t0.5).
Fitting the experimental data to this model allows us to determine the intraparticle diffusion rate constant (ki) and evaluate the contribution of intraparticle diffusion to the overall adsorption process. Additionally, the intercept (C) provides insights into the boundary layer effect, which represents the contribution of external surface adsorption to the overall adsorption mechanism.
2.2.1. Pseudo First Order Model
The adsorption kinetic data of Direct Red 83:1 was analyzed using three different kinetics models: the pseudo-first-order model, the pseudo-second-order model, and the intra-particle diffusion model. The aim was to investigate the adsorption kinetics and mechanisms involved in the adsorption process of the dye onto chitosan magnetic polymers. To determine the most appropriate model, the linear determination coefficient (R2) was used as a criterion for selection.
The plots for the pseudo-first-order, pseudo-second-order, and intra-particle diffusion models for the adsorption of Direct Red 83:1 onto the chitosan magnetic adsorbent are presented in Figures 2, 3 and 4, respectively. These plots showcase the fitting of the experimental data to the respective kinetic models, enabling a visual assessment of their agreement with the observed adsorption behavior.
The results and corresponding kinetic parameters, as well as the correlation coefficients (R
2), are reported in
Table 1. By comparing the R
2 values for each model, the best-fitting model can be identified, providing valuable insights into the dominant adsorption mechanisms and the overall kinetics of the adsorption process for Direct Red 83:1 onto the chitosan magnetic adsorbent.
2.2.2. Pseudo Second Order Model
As shown in
Figure 2, the plot of t/qt versus t resulted in straight lines throughout the entire range of measurement. However, despite this linearity, the R
2 values obtained from fitting the experimental data to the pseudo-first-order model varied, and the calculated q
e (equilibrium adsorption capacity) values were notably different from the experimental q
e values. As a consequence, the pseudo-second-order kinetic model emerged as the most suitable model to describe the adsorption kinetics of the dye onto the adsorbent surface.
The application of the pseudo-second-order model indicates that the adsorption process is primarily governed by chemisorption or chemical adsorption, and it is the rate-limiting step that controls the overall adsorption process. This finding aligns with the observation of similar kinetics in various other studies involving different adsorbent systems, such as the adsorption of Malachite Green using magnetic-cyclodextrin-graphene oxide nanocomposites [
16], the adsorption of methyl orange and Pb (II) on β-CD and polyethyleneimine bi-functionalized magnetic nano-adsorbent [
17], the removal of Eu (III) onto Fe
3O
4 cyclodextrin magnetic composite [
18], and the removal of organic pollutants on magnetic β-cyclodextrin porous polymer nano-spheres [
19].
These collective findings provide strong support for the pseudo-second-order kinetic model as a reliable and consistent approach to understand and characterize the adsorption kinetics of different adsorbent systems for various pollutants.
The adsorption process is a multi-step phenomenon involving several stages. It begins with the transport of solute molecules (dye) from the aqueous solution to the surface of the adsorbent (polymer). Subsequently, the dye molecules diffuse onto the adsorbent, leading to their adsorption. To gain a comprehensive understanding of the adsorption of Direct Red 83:1 onto magnetic polymers, the kinetics of the adsorption process were analyzed using the intraparticle diffusion model. This analysis aims to determine whether intraparticle diffusion plays a significant role and if it is the rate-limiting step in the overall adsorption process.
The intraparticle diffusion model provides insights into the role of intraparticle diffusion in governing the adsorption process. By plotting the amount of Direct Red 83:1 dye absorbed against the square root of time (
Figure 3), the influence of intraparticle diffusion on the adsorption kinetics can be assessed. This approach follows the work of Weber in 1963 [
15] and Crini in 2007 [
20], who employed this model to study similar adsorption phenomena.
Analyzing the adsorption process using the intraparticle diffusion model allows researchers to ascertain the contribution of intraparticle diffusion and its significance in determining the rate of the adsorption process. Understanding this aspect is vital in comprehending the overall kinetics and mechanisms of adsorption of Direct Red 83:1 onto the magnetic polymers.
2.2.3. Intraparticle Diffusion Model
Figure 4 illustrates the plot of qt (amount of Direct Red 83:1 dye adsorbed) versus t
0.5 (square root of time) for the intraparticle diffusion of Direct Red 83:1 onto chitosan at different concentrations of the dye. For the chitosan magnetic polymer, two distinct regions were observed at high concentrations (200 and 300 mg/L), while a single linear portion was seen at low concentrations (50, 100, and 150 mg/L).
At low concentrations, the single linear portion indicates that the adsorption process inside the polymer is primarily controlled by the diffusion of dye molecules to the polymer surface. This suggests that the intraparticle diffusion process is the rate-limiting step for adsorption at these lower concentrations.
However, at high concentrations, the graph displays two zones. The initial curved portion is attributed to a rapid diffusion adsorption stage, which is influenced by the boundary layer effect. This effect arises due to the high concentration of dye in the bulk solution, leading to a faster initial adsorption rate. The subsequent linear stage reflects the gradual adsorption due to intraparticle diffusion, indicating the further penetration of dye molecules into the porous structure of the polymer.
The summary of kinetic parameters obtained from these intraparticle diffusion plots is presented in
Table 1, providing valuable insights into the different stages and mechanisms involved in the adsorption process of Direct Red 83:1 onto the chitosan magnetic polymer at various dye concentrations.
The similar results reported for the adsorption of Direct Blue 78 onto β-CDs-EPI and chitosan-NaOH, as well as for the adsorption of Direct Red dye using CDs-EPI polymers, indicate that the adsorption behaviour and kinetics of these dye-polymer systems share common characteristics. These studies have observed comparable adsorption patterns, suggesting that the adsorption mechanisms and rate-limiting steps are consistent across different adsorbent materials and dye molecules.
The findings from these studies provide valuable insights into the effectiveness of CDs-EPI and chitosan-NaOH polymers as potential adsorbents for the removal of Direct Blue 78 and Direct Red dye from aqueous solutions [
5,
12,
21]. The consistency of the results across multiple studies strengthens the understanding of the adsorption processes and highlights the potential applicability of these adsorbents in water treatment and wastewater remediation [
22].
2.3. Adsorption Equilibrium
The adsorption isotherm plays a crucial role in describing how the dye interacts with the adsorbent and provides an understanding of the adsorption capacity of the adsorbent [
23]. The Freundlich isotherm model is commonly used for adsorption on heterogeneous surfaces and in cases where multilayer sorption occurs. The application of this model also suggests that the energy decreases as the sorption centers of the adsorbent become saturated [
24].
The linearized form of the Freundlich isotherm model is expressed by the following Equation (4):
where
qe is the equilibrium dye concentration on absorbent (mg/g),
Ce the equilibrium dye concentration in solution (mg/L),
KF the Freundlich constant (L/g) and
1/nF is the heterogeneity factor. The plot of ln
qe versus ln
Ce was used to determine the intercept value of
KF and the slope of
1/nF.
By analyzing the experimental data and fitting them to the Freundlich isotherm model, valuable information about the adsorption capacity, intensity, and favorability of the adsorption process can be obtained, providing insights into the interaction between the dye and the adsorbent.
The linearized form of the Langmuir model is given by the following Equation (5):
where
Ce (mg/L) and
qe (mg/g) are the liquid phase concentration and solid phase concentration of adsorbate at equilibrium respectively.
KL (L/g) and
aL (L/mg) are the Langmuir isotherm constants. Plotting
Ce/
qe versus
Ce is possible to know the value of
KL from the intercept (1/
KL) and the value of
aL from the slope (
aL/
KL);
qmax is the maximum adsorption capacity of the polymer and is defined by
KL/aL.
The most important characteristic of this isotherm can be described using a dimensionless constant (
RL) which is called separation factor and is described by the following Equation (6):
The value of the separation factor (RL) provides insight into the nature of adsorption as follows:
RL >1: Indicates unfavorable adsorption, meaning the adsorption capacity decreases with increasing equilibrium concentration. This suggests that the adsorption process is less efficient at higher dye concentrations.
RL =1: Indicates linear adsorption, implying a constant adsorption capacity regardless of the equilibrium concentration.
0<RL<1: Indicates favorable adsorption, where the adsorption capacity increases with increasing equilibrium concentration. This suggests that the adsorption process is more efficient at higher dye concentrations.
The separation factor (RL) is a valuable tool for assessing the feasibility and efficiency of the adsorption process and helps to understand the type of adsorption occurring on the adsorbent surface.
Temkin isotherm is an alternative model that assumes a linear decrease in the heat of sorption as the adsorption process progresses, in contrast to the logarithmic relationship in the Freundlich equation. The Temkin isotherm takes into account the interactions between the adsorbate and the adsorbent, considering the adsorbent’s finite surface area and the heterogeneity of the surface [
25].
The Temkin isotherm equation suggests a linear decrease in sorption energy as the adsorption sites on the adsorbent become saturated. This implies that the heat of sorption gradually decreases as the adsorption progresses.
The linearized form of the Temkin isotherm model is expressed by the following Equation (7):
where
bT is the Temkin constant, related to the heat of adsorption (kJ/mol),
aT is the constant of Temkin isotherm (L/g),
R the universal gas constant (8.314 J/mol K) and
T is the absolute temperature in Kelvin.
2.3.1. Freundlich Isotherm
The Freundlich isotherm is applicable to both monolayer (chemisorption) and multilayer adsorption (physisorption) processes and is based on the assumption that the adsorbate adsorbs onto the heterogeneous surface of the adsorbent. The Freundlich isotherm [
25] provides two important constants, K
F and n
F, which reflect the adsorption capacity and intensity, respectively. For a favorable adsorption, the Freundlich constant n
F typically falls within the range of 1 to 10. A higher value of n
F (equivalently, a smaller value of 1/n
F) indicates a stronger and more effective interaction between the adsorbent and the adsorbate. When 1/n
F is less than 1, it corresponds to a normal L-type isotherm, while 1/n
F greater than 1 indicates cooperative sorption.
According to the results presented in
Table 2 and
Figure 5, the K
F value for the magnetic chitosan polymer was found to be 1.28. The most crucial parameter related to this isotherm is n
F, which represents the heterogeneity factor. In the case of the chitosan magnetic polymer, the value of n
F was determined to be 1.68, falling within the favorable adsorption range of 1 to 10. This suggests that the adsorption process is favored and that there is an effective interaction between the adsorbent and the adsorbate. The experimental data were fitted to the Freundlich model, resulting in a determination coefficient of 0.844 for the chitosan polymer. This coefficient represents the goodness of fit of the experimental data to the Freundlich isotherm model for the adsorption process.
2.3.2. Langmuir Isotherm
The Langmuir adsorption isotherm is based on the assumption of homogeneous adsorption, where adsorption occurs in a monolayer coverage on the adsorbent surface. This isotherm assumes that all adsorption sites are equivalent, leading to uniform surface coverage. According to the Langmuir model, the ability of a molecule to adsorb at a specific site is independent of the occupation of neighboring sites.
For the Langmuir isotherm, it is essential to analyze the q
max value, which represents the maximum adsorption capacity of the adsorbent under specific experimental conditions. In the case of the magnetic chitosan polymer, the q
max value was found to be 17.46 (
Table 2,
Figure 6A). This value provides valuable information about the adsorption capacity of the adsorbent and the maximum amount of dye that can be adsorbed per unit weight of the adsorbent.
2.3.3. Temkin Isotherm
The experimental data were also analyzed using the Temkin isotherm. The best fit for the magnetic chitosan polymer was achieved using this model with a determination coefficient of 0.946, indicating that the adsorption occurred on heterogeneous surfaces [
24]. The Temkin isotherm model assumes that heat of adsorption (function of temperature) of all molecules in the layer would decrease linearly rather than logarithmic with coverage. The Temkin isotherm model mainly describes the chemical adsorption process as electrostatic interaction. In the present study, the Temkin isotherm fitted the observed results well, hence, electrostatic interaction is an important mechanism affecting the interaction be-tween the chitosan adsorbent and the pollutant [
25].
Figure 7.
Temkin isotherm plot for chitosan magnetic polymer.
Figure 7.
Temkin isotherm plot for chitosan magnetic polymer.
2.3. Polymer Characterization
2.3.1. Thermogravimetric Analysis
Different pattern of decomposition was obtained for the original chitosan and the iron oxide coated polymer (
Figure 8).
Together with a first mass loss, probably due to water hydration (below 150° C) a single decomposition path was observed for chitosan. Common elimination of volatile compounds at similar temperatures was also recorded for the iron polymer, followed by another four not completely defined steps leading to a total mass loss of 36.05%. No defined major compounds or stable intermediates were identified within the MS fragments, with multiple masses between 30 and 73 amu, probably corresponding to the decomposition series of oxygen and nitrogen single bonded compounds Cn H2n+1 O and Cn H2n+2 N.
The major weight loss in chitosan above 300 degrees, where most of the organic matter is lost, was not observed in the iron containing polymer. However, the second and third losses, occurring between 150 and 400 degrees are very similar to those suffered by the iron particles although to a higher extent (8.42 and 4.09% in the polymer towards 1.56 and 0.85% in the iron particles).
A last significant loss of 7.65% by weight above 700 degrees was registered leaving a total mass loss almost halved (from 67.91% to 36.05%) than original chitosan with about 64% of matter remaining after the constant heating up to 800 °C (
Figure 8B).
From thermogravimetry, the weight percentage of chitosan in the Chitosan-Fe polymer can be calculated as follows [
26]:
where the Ws was the weight loss (%) from 120 to 800 °C for the chitosan-Fe samples, and Wn was the weight loss (%) from 120 to 800 °C for the pure Fe
3O
4 nanoparticles. After calculation, the weight percentage of chitosan in chitosan-Fe reached 31.5 %. In the case of Jaafari’s chitosan-Fe nanoparticles [
24] they obtain 16.3%. The double content in organic matter obtained for the in situ generation of iron particles could suggest a higher efficiency in the incorporation into the structure of the polymer.
2.3.2. DRX Powder X-ray Diffraction
A degree of crystallinity between 68-70% was recorded in the material for the independently prepared iron particles (PI) and for the chitosan polymer formed in the presence of iron particles (CHI-Fe). The characteristic face-centred cubic lattice Fd-3m of magnetite (ICDD database, standard sample PDF 04-006-6551) observed for independently PI particles diffraction pattern is majorly found in the final chitosan iron polymer (
Figure 9).
However, novel diffractions were also observed together with those of magnetite which could be attributed to a distortion of the crystalline structure or to a substitution of magnetite by haematite in some lattices. This last process has been previously described by a partial oxidative substitution of Fe
2+ or through a redox-independent dissolution-reprecipitation reaction in the acidic media in which they are prepared [
27].
2.3.3. IRTF
Similar profiles were observed in the spectrum of commercial chitosan and iron containing chitosan polymer (
Figure 10).
The absorption bands at 2916 and 2868 cm−1 attributed to C-H symmetric and asymmetric stretching, respectively, were also recorded for the polymer with a lower definition. The stretching vibration of the double bound C=O (1654 cm-1) with minor intensity than the N–H (3355.42 cm-1), typical for chitosan and related polysaccharides, was also found for both compounds although at slightly higher frequencies in the case of the polymer.
This shift of the bands to lower frequencies (C=O bands 1654.16 and 1678.75 and NH at 1589 to 1560 cm-1 can be attributed to the establishment of new hydrogen bridge interactions that weaken the bonds vibrating in the structure. On the other hand, both the C-N stretching vibration bands, originally located at 1419.85 cm-1, 1378.63 and 1319.67 cm-1 and the C-O bands at 1062.12, 1022.52 and 984.01 cm-1 are shifted to higher frequencies 1459.14, 1407, and 1057.84, 1018.19, which is an evidence of an increasement in the strain of the structure when forming the coating for the iron particles.
Finally, a significant band corresponding to the vibration of the iron-oxygen bonds at 534 cm-1 was also identified in the spectrum.
2.3.4. SEM and EDX
Several representative polymeric macrostructures were analyzed by field emission scanning electron microscopy to characterize the morphology of iron coated chitosan polymer. Highly regular spherical structures were observed for the polymer (
Figure 11A). However, some small dark areas dotted the structure suggesting the presence of zones where the coating of iron nanoparticles seems to be not completed (
Figure 11B).
The EDX analysis revealed a huge difference in iron content for both structures (
Figure 12), being virtually absent in clear areas and the majority in dark ones, and thus supporting the previous assumption.
2.3.5. Molecular Models
The interaction of chitosan with iron was eventually assessed by using density functional theory (DFT). Simulations were performed with divalent cations and glucosamine monomers, a model system by Gomes and co-workers that correctly mimics the stability of chitosan complexes with metallic ions [
28]. As discussed by these authors, metal centres can be coordinate to chitosan by hydroxyl and/or amino groups at five available sites (labelled as C1 – C5).
Figure 13 illustrates the optimized complexes with metals at these positions. Values in parenthesis stand for the computed relative energies in kcal/mol, where the global minimum is taken as reference. As noted, the C3 coordination leads to the most stable complex, while all other coordination patterns are associated to less favourable (more positive) relative energies, which lie in the range of 13.14 – 29.21 kcal/mol. Indeed, C3 outperforms both the anchoring to the amino groups of the chitosan entity (C4 and C5) and the alternative C1 and C2 sites, even if the former established three interactions (two hydroxyl contacts + the heterocyclic oxygen). A close inspection of
Figure 13 reveals that the large stability at C3 arises from a concomitant proton transfer equilibrium from the hydroxyl to the amino group. Consequently, theory predicts that iron coordination induces a protonation transfer that leads to a deprotonated hydroxyl, which in turn enhances the coordination. This computational outcome agrees the recorded IR spectrum, where Fe – O bands are observed while Fe – N peaks were not present.
3. Materials and Methods
3.1. Chemicals and Reagents
Chitosan, iron (III) chloride hexahydrate (FeCl3·6H2O), ammonium hydroxide, and ethanol were procured from Sigma-Aldrich (Spain). Meanwhile, iron (II) chloride tetrahydrate (FeCl2·4H2O) was obtained from Fluka (Spain). The Direct Red 83:1 dye used in the experiments was supplied by AITEX (Asociación de Investigación de la Industria Textil, Alcoy, Spain).
3.2. Iron Nanoparticles Preparation
In order to compare possible differences in the structure and efficacy of the incorporation into the polymer between iron oxide nanoparticles performed by classical co-precipitation under alkaline conditions and in situ formation in acidic media during polymerization of chitosan, magnetic nanoparticles were also prepared following a standard procedure as follows. 3.49 g of FeCl2·4H2O and 9.5 g of FeCl3·6H2O were sequentially dissolved in 100 mL of water at room temperature. After complete solution, 30 mL of ammonium hydroxide was added drop by drop and the mixture was subjected to magnetic stirring for 30 minutes at 80 °C to promote proper reaction and particle formation. The mixture was then centrifuged at 4000 rpm for 10 minutes to separate the solid precipitate from the liquid solution. The resulting precipitate was then washed twice with a mixture of ethanol and water in a 1:1 ratio to remove basic medium and unreacted reagents. Finally, the washed precipitate was left to dry overnight, resulting in the formation of magnetic nanoparticles. These magnetic nanoparticles are now ready for use in further experiments or applications.
3.3. Chitosan-Fe Polymer Preparation
To prepare the chitosan magnetic beads, the following procedure was followed. 5 g of chitosan powder was dissolved in 250 mL of a 5% acetic acid aqueous solution. The mixture was stirred for 2 hours at 50 °C to obtain a homogeneous gel. After achieving a uniform gel, 16 g of FeCl2·4H2O and 32 g of FeCl3·6H2O were added to the gel. The mixture was further stirred for 1 hour at 50 °C to facilitate the incorporation of iron ions into the gel. To form the polymer beads, the homogeneous gel containing chitosan and iron ions was carefully dropped into a bath containing a 30% aqueous solution of NaOH. Upon contact with the NaOH solution, the gel rapidly transformed into spherical particles, thus forming the chitosan magnetic beads. The newly formed beads were then washed with water to remove any residual chemicals or impurities. Finally, the chitosan magnetic beads were left to dry at room temperature, resulting in the final product, ready for use in the experiments.
3.4. Dye solution Preparation
To conduct the adsorption experiments, various concentrations of Direct Red 83:1 were prepared. Direct Red 83:1 is an azo dye with the chemical formula C33H20N6Na4O17S4 and a molecular weight of 992.77 g/mol.
The concentrations of Direct Red 83:1 used in the adsorption tests for both adsorbents were as follows: 50 mg/L, 100 mg/L, 150 mg/L, 200 mg/L, and 300 mg/L. These different concentrations were employed to investigate and evaluate the adsorption properties of the respective adsorbents under varying dye concentrations. By analyzing the adsorption behavior at different dye concentrations, a comprehensive understanding of the adsorption process and efficiency can be obtained.
3.5. Analyses and Data Evaluation
In the adsorption experiments, the concentration of the dye in the supernatant was measured using a spectrophotometer, specifically the Shimadzu UV-1700 model (Kyoto, Japan). The measurements were taken after reaching the equilibrium state between the dye and the adsorbents. The absorbance of the dye solution was measured both before and after the treatment with the polymers at the wavelength of maximum absorbance for the dye, which is λmax = 526 nm. The molar absorptivity (extinction coefficient) of the dye at this wavelength is ε526 = 1065 M-1cm-1.
By comparing the absorbance values before and after the treatment, the extent of dye adsorption onto the polymers can be determined. This analysis helps in evaluating the effectiveness of the adsorbents in removing the dye from the solution at different concentrations.
3.6. Adsorption Experiments
The adsorption tests were conducted at a constant temperature of 25 °C, utilizing solutions containing various concentrations of the dye ranging from 50 to 300 mg/L. For each experiment, 1g of the polymer was mixed with 50 mL of the dye solution. The mixture was continuously stirred at a fixed speed of 500 rpm throughout the experiment to ensure uniform contact between the dye and the adsorbent.
To determine the amount of dye that was not retained by the polymer at specific time intervals (every 10 minutes), external magnets were employed. The mixture was subjected to the influence of these magnets for 5 minutes, facilitating the separation of the magnetic polymer beads along with the adsorbed dye from the solution.
Once the separation was accomplished, the concentration of the remaining dye in the solution was measured using the spectrophotometer, which allows for precise and accurate quantification of the dye’s concentration. This procedure was repeated for each concentration of dye tested, and the data obtained from the spectrophotometer measurements were used to assess the adsorption efficiency of the polymer under varying dye concentrations and contact times.
The amount of dye adsorbed on the polymers (q
t) in mg/g was determined using a mass balance approach, as described by Renard et al. in 1997 [
29]. The equation used for this determination is given by Equation (9):
where:
qt is the amount of dye adsorbed on the polymer at time t (mg/g).
Co is the initial concentration of dye in the solution (mg/L).
Ct is the concentration of dye in the solution at time t (mg/L).
V is the volume of the dye solution used (L).
m is the mass of the polymer used (g)
3.7. Adsorption Kinetics
The study of kinetics holds significant importance in adsorption research as it allows us to predict the rate at which pollutants are eliminated from aqueous solutions and provides valuable insights into the mechanism of sorption reactions. To gain a deeper understanding of the mechanism behind dye adsorption onto the newly developed adsorbents, several kinetic models, including the pseudo first order, pseudo second order, and intra-particle diffusion models, were carefully examined using the experimental data. These models were utilized to analyze and characterize the adsorption process, enabling a comprehensive investigation of the interactions between the adsorbents and the dye molecules.
3.8. Isotherm Analysis
The adsorption isotherm plays a critical role in understanding how adsorbing molecules are distributed between the liquid and solid phases once the adsorption process reaches an equilibrium state. To optimize the design of an adsorption system for efficiently removing Direct Red 83:1 from solutions, it becomes crucial to identify the most suitable correlation for the equilibrium studies. In this study, three different isotherms, namely Freundlich, Langmuir, and Temkin isotherms, were employed and fitted to the experimental data.
By utilizing these isotherms, the adsorption behavior of Direct Red 83:1 onto the adsorbents can be thoroughly investigated, and the appropriate model can be selected based on its best fit to the experimental data. This choice will aid in obtaining valuable information about the adsorption capacity, surface characteristics, and interactions between the adsorbent and the dye molecules, ultimately contributing to the optimal design and performance of the adsorption system.
3.9. Polymer Characterization
Thermogravimetric analysis was performed on a thermobalance model TGA/DSC 1 HT (Mettler-Toledo GmbH, Schwerzenbach, Switzerland) coupled to a Balzers Thermostar mass spectrometer (Pfeiffer Vacuum, Asslar, Germany) for gas analysis. Scan bargraph cycles were performed in the range 15–94 m/z in the quadrupole mass spectrometer (QMS 200 M3), with a dwell time of 2 s for every ion and cathode voltage in the ion source of 65 V.
About 10 mg of sample were used to study the thermal decomposition from 30 to 800 °C at a heating rate of 10 °C/min under N2 atmosphere with a flow rate of 50 mL/min. During the experiment, the sample was placed in an uncovered sapphire crucible of 70 μL capacity. The result was corrected with a blank curve. For each sample, the m/z-ratio profiles in the range 15-94 amu, the experimental mass loss (TG curve) and its first derivative (DTG curve) were recorded.
The crystalline phases in pre-formed iron nanoparticles and within the polymer were also analysed. An x-ray powder diffractometer in mode θ-θ (Bruker D8 Advance, Bruker Corporation, Billerica, MA, USA), using Cu Kα radiation at 40 kV and 30 mA and one dimensional detector with a 2° window. The primary optics consisted of a 2° Soller slit, a 1 mm incidence slit and an air anti-scattering screen. The secondary optics included an 8 mm anti-scattering slit, a Ni filter and a 2.5° Soller slit.
Approximately 1 g of powder reduced iron coated macrostructures was placed on a back-loaded sample holder and analysed at steps of 0.05º (dwell time of 1 s/step) over an angular range of 10-70° and 30 rpm of angular velocity. The powder diffraction file was evaluated with the linked software (DIFFRAC.EVA 5.2, Bruker AXS, 2020, MA, USA) and a crystalline powder database (PDF-4+ 2021, ICDD).
The infrared study was performed on a Nicolet 5700 spectrometer (Nicolet, Madison, WI, USA) equipped with a Ge/KBr beam splitter, a DTGS-KBr detector and a ceramic infrared source. The sampling module was the Smart Orbit diamond attenuated total reflectance attenuator accessory. The samples were analyzed finely divided and pre-dried at 60 °C. Spectral analysis was performed in the range 4000 to 400 cm-1, with a resolution of 4 cm-1.
FESEM uncovered samples mounted with the aid of conductive carbon cement on aluminum stubs. FESEM-EDX analyses were performed with a FE-SEM (ApreoS Lovac IML Thermofisher, MA, USA) and coupled with an Octane Plus EDAX microanalyzer (AMETEK, PA, USA). Sample analysis were performed between 1 - 2 kV, and 20 kV for EDX analysis.
3.10. Computational Calculations
All DFT calculations were performed with the Gaussian 16 suite of programs [
30]. The molecular structures were fully optimized with the dispersion-corrected version of the B3LYP functional (B3LYP-D3) [
31,
32]. This functional has been successfully used for simulating chitosan-metal complexes [
33]. The Def2SVP basis set was used for all atoms except Fe, which was treated with the relative effective core potentials (Def2-ECP) [
34,
35]. Frequencies were also computed to confirm that the localized structures correctly correspond to real minima in the potential energy surface.
Figure 1.
Influence of contact time on chitosan adsorption capacity at various dye concentrations: 50 mg/L (•), 100 mg/L (■), 150 mg/L (♦), 200 mg/L (▲), and 300 mg/L (○).
Figure 1.
Influence of contact time on chitosan adsorption capacity at various dye concentrations: 50 mg/L (•), 100 mg/L (■), 150 mg/L (♦), 200 mg/L (▲), and 300 mg/L (○).
Figure 2.
Pseudo first order model plots for chitosan at different dye concentrations 50 mg/L (•), 100 mg/L (■), 150 mg/L (♦), 200 mg/L (▲) and 300 mg/L (○).
Figure 2.
Pseudo first order model plots for chitosan at different dye concentrations 50 mg/L (•), 100 mg/L (■), 150 mg/L (♦), 200 mg/L (▲) and 300 mg/L (○).
Figure 3.
Pseudo second order model plots for chitosan at different dye concentrations 50 mg/L (•), 100 mg/L (■), 150 mg/L (♦), 200 mg/L (▲) and 300 mg/L (○).
Figure 3.
Pseudo second order model plots for chitosan at different dye concentrations 50 mg/L (•), 100 mg/L (■), 150 mg/L (♦), 200 mg/L (▲) and 300 mg/L (○).
Figure 4.
Intraparticle diffusion model plots for chitosan at different dye concentrations 50 mg/L (•), 100 mg/L (■), 150 mg/L (♦), 200 mg/L (▲) and 300 mg/L (○).
Figure 4.
Intraparticle diffusion model plots for chitosan at different dye concentrations 50 mg/L (•), 100 mg/L (■), 150 mg/L (♦), 200 mg/L (▲) and 300 mg/L (○).
Figure 5.
Freundlich isotherm plot for chitosan magnetic adsorbent.
Figure 5.
Freundlich isotherm plot for chitosan magnetic adsorbent.
Figure 6.
(A). Langmuir isotherm plot. (B) Separation factor.
Figure 6.
(A). Langmuir isotherm plot. (B) Separation factor.
Figure 8.
A TG of chitosan-Fe (Blue line), chitosan (Grey Line) and Fe (Orange line). B. DTG of chitosan-Fe (Blue line), chitosan (Grey Line) and Fe (Orange line).
Figure 8.
A TG of chitosan-Fe (Blue line), chitosan (Grey Line) and Fe (Orange line). B. DTG of chitosan-Fe (Blue line), chitosan (Grey Line) and Fe (Orange line).
Figure 9.
DRX powder X-Ray Diffraction of Chitosan-Fe (Blue line) and Fe (Orange line).
Figure 9.
DRX powder X-Ray Diffraction of Chitosan-Fe (Blue line) and Fe (Orange line).
Figure 10.
IRTF of Chitosan (Grey line), Chitosan-Fe (Blue line) and Fe (Orange Line).
Figure 10.
IRTF of Chitosan (Grey line), Chitosan-Fe (Blue line) and Fe (Orange Line).
Figure 11.
A: Chitosan Polymer. B: Detail of small dark areas of the chitosan polymer.
Figure 11.
A: Chitosan Polymer. B: Detail of small dark areas of the chitosan polymer.
Figure 12.
EDX of Chitosan (A) and Chitosan-Fe (B).
Figure 12.
EDX of Chitosan (A) and Chitosan-Fe (B).
Figure 13.
Optimized Chitosan-Fe(II) model systems. Colour atom code: grey: carbon, red: oxygen, blue: nitrogen, cyan: iron.
Figure 13.
Optimized Chitosan-Fe(II) model systems. Colour atom code: grey: carbon, red: oxygen, blue: nitrogen, cyan: iron.
Table 1.
Kinetic Parameters for Adsorption of Direct Red 83:1 onto Chitosan Magnetic Polymer.
Table 1.
Kinetic Parameters for Adsorption of Direct Red 83:1 onto Chitosan Magnetic Polymer.
PFOM 1
|
Chitosan-Fe |
C0 (mg/L) |
qeexp
|
qecal
|
K1(min-1) |
R2
|
50 |
2.27 |
1.05 |
0.033 |
0.987 |
100 |
4.69 |
3.11 |
0.040 |
0.881 |
150 |
6.97 |
5.31 |
0.031 |
0.960 |
200 |
8.60 |
9.36 |
0.027 |
0.950 |
300 |
12.47 |
18.29 |
0.029 |
0.772 |
PSOM 2 |
Chitosan-Fe |
C0 (mg/L) |
qeexp
|
qecal
|
K2(min-1) |
R2
|
50 |
2.27 |
2.27 |
- |
0.999 |
100 |
4.69 |
4.69 |
- |
0.999 |
150 |
6.97 |
6.97 |
- |
0.999 |
200 |
8.60 |
8.60 |
- |
0.999 |
300 |
12.47 |
12.47 |
- |
0.999 |
IDM 3 |
Chitosan-Fe |
C0 (mg/L) |
qeexp
|
qecal
|
Ki (mg/g min1/2) |
R2
|
50 |
2.27 |
1.06 |
0.11 |
0.930 |
100 |
4.69 |
1.50 |
0.32 |
0.662 |
150 |
6.97 |
2.08 |
0.46 |
0.834 |
200 |
8.60 |
0.63 |
0.71 |
0.978 |
300 |
12.47 |
-1.37 |
1.20 |
0.977 |
Table 2.
Freundlich, Langmuir and Temkin parameters.
Table 2.
Freundlich, Langmuir and Temkin parameters.
Isotherm |
Parameter |
Chitosan-Fe |
Freundlich |
KF
|
1.28 |
nF
|
1.68 |
R2
|
0.844 |
Langmuir |
qmax
|
17.46 |
KL
|
0.76 |
aL
|
0.043 |
R2
|
0.884 |
RL
|
0.13-0.024 |
Temkin |
aT
|
0.48 |
bT
|
0.67 |
R2
|
0.946 |