3.1. Material Characterization
The morphological characteristics of examined samples were analyzed by scanning electron microscopy (
Figure 1). The main morphological characteristic of cotton fibers within the structure of Cott and Cott/PES (
Figure 1a and 1d, respectively) samples is the appearance of a spirally twisted ribbon. Along the rough surface of cotton fibers, the longitudinal cracks are visible [
21,
22]. The straight and smooth fiber visible in the structure of the Cott/PES sample (
Figure 1d) represents a PES component. After carbonization cotton and cotton/polyester yarns retain a fibrous structure which is visible in SEM photographs (
Figure 1b,e). Carbonized cotton fibers retain a spirally twisted structure with more visible microfibrils that contribute to an increase in surface roughness. The surface of carbonized cotton fibers within the structure of the Cott/PES
c sample (
Figure 1e) is additionally decorated by evenly distributed particles. During the pyrolysis, depolymerization of the PES component occurs, leading to the tearing of the PES filament and the creation of shorter linear segments of fibers with a significantly lower degree of polymerization. Also, the pyrolysis products obtained as a consequence of PES decomposition are evenly deposited on the surface of the cotton component remaining as condensed particles after carbonization.
Unlike the carbonized samples that retain the fibrous form of the precursor, the samples obtained after activation become breakable and lose their fibrous form, acquiring the form of needle powder (which is visible from the optical photographs in
Figure 1c,f). However, the SEM photographs (
Figure 1c,f) show that, at the micro level, these samples retain the precursor structure. Additionally, on the SEM photograph of sample Cott/PES
ac (
Figure 1f) only the presence of the activated cotton component is noticeable, while the presence of the activated PES component is not visible. It may be assumed that two subsequent pyrolysis processes were too aggressive for the PES component, leading to its decomposition and disappearance.
Textural characteristics of carbonized and activated samples are given in
Table 1. Although non-carbonized samples contain rough surfaces with longitudinal cracks, their specific surface area is not developed and cannot be measured. The carbonization of cotton and cotton/PES yarn does not bring a considerable increase in specific surface area and porosity. The Cott
c sample showed very low values of specific surface area and pore volumes, while the specific surface area and porosity of sample Cott/PES
c were below the measurement limit. This immeasurable surface area may be the consequence of the presence of condensed material, visible in
Figure 1e, which blocks the pores on the surface of the carbonized cotton component. As can be seen from
Table 1, high values of specific surface area and pore volume are obtained for activated samples. It is evident that starting from cotton and cotton/polyester yarn as carbon precursors, microporous materials with a highly developed specific surface area can be obtained by carbonization and activation in the presence of KOH. This double pyrolysis procedure of cotton/polyester yarn leads to some higher S
BET and S
micro values in the case of Cott/PES
ac, due to the decomposition of the synthetic component (the disappearance of the PES component is confirmed by SEM photographs in
Figure 1f) with the opening of the pores upon activation in the presence of KOH.
The surface charge of non-carbonized samples is determined through the values of the isoelectric point: for Cott pHIEP = 2.25 and Cott/PES pHIEP = 1.34. The isoelectric point indicates that the surface of examined material is positively charged below, and negatively charged above the pHIEP value. For the carbonized samples point of zero charge is in the basic pH region indicating that the surface of carbonized samples is negatively charged in the surrounding solution having a pH value above 8.71 for Cottc and 8.19 for Cott/PESc. On the other hand, activated samples have a more acidic surface with the pHPZC 5.11 for Cottac and 6.43 for Cott/PESac.
FTIR analysis (
Figure 2) was used to examine the type of surface functional groups of unmodified and thermally modified yarn samples. FTIR spectra of untreated samples (Cott and Cott/PES) show a broad band between 3350 and 3250 cm
−1 that originates from cellulose in the structure of cotton component: peaks at 3340 and 3279 cm
−1, corresponds to the intramolecular O(3)H…O(5) and intermolecular O(6)H…O(3) hydrogen bonds, respectively [
23], while shoulder at 3410 cm
−1 can be assigned to the O(2)H…O(6) intramolecular hydrogen bonds [
23]. Symmetrical and asymmetrical vibrations of the C-H bond in methyl and methylene groups of cellulose give two peaks at wavenumbers of 2850 cm
−1 and 2920 cm
−1 [
24]. The broad peak around 1630 cm
−1, can be assigned to the OH bending of adsorbed water [
23], aromatic skeletal vibration, or C=O stretching vibrations in carbonyl groups of hemicelluloses [
24].
The band at 1365 cm
−1, observed in
Figure 2a, may be attributed to C-H bending vibrations in cellulose and hemicellulose, while the bands centered at 1313–1317 cm
−1 and 1336 cm
−1 are assigned to CH
2 wagging, characteristic of crystalline cellulose, and the C–O–H in-plane bending from amorphous cellulose, respectively [
25]. The C-O stretching around 1162 cm
−1 represents asymmetric bridge stretching of C-O-C groups in cellulose, while the band at 1110 cm
−1 corresponds to asymmetric glucose ring stretching. The low-intensity peak at 894 cm
−1 indicates the presence of
β-glycosidic linkages between monosaccharides [
26], and the peak at 668 cm
−1 can be assigned to C–OH out-of-plane bending [
23]. Compared to cotton, the FTIR spectrum of cotton/polyester yarns showed some additional peaks originating from polyester. The presence of ester groups in a polyester component is confirmed by the intense peaks at 1710 cm
−1 and 1240 cm
−1. Also, the peak at 1505 cm
−1 is assigned to the skeletal vibrations of the aromatic systems in polyester chains [
27], while out-of-plane bending vibrations of the benzene ring in the polyester appear at 870 cm
−1 (C-C) and 720 cm
−1 (C-H and C=O) [
27]. The FTIR spectra of the carbonized and activated samples also show a broad band in the region of 3450-3650 cm
−1, which originates from stretching vibrations of the O-H bond in carboxyl or hydroxyl groups. Peaks at 2855 cm
−1 and 2925 cm
−1 originate from the symmetrical and asymmetrical vibrations of the C-H bond in methyl and methylene groups [
28], while the doublet present at 2360 and 2340 cm
−1 originate from the carbon dioxide. Also, a band appearing at 1640 cm
−1 corresponds to the bending vibrations of the O-H bond or the stretching vibrations of the aliphatic C=C bond, while the peak at 1384 cm
−1 belongs to the deformational vibration of the C-O bond in the carboxyl group [
29].
3.2. Adsorption Experiments
The initial pH value of the adsorbate solution influences the adsorption efficiency by dictating the form of dye molecule in the solution, as well as, the charge of the adsorbent surface depending on pH
PZC. In the adsorbate solution with a pH below 4, the molecular form of MB and the zwitterion form of MO are predominating. On the other hand, with a pH increase above 4, the cationic form of MB and anionic form of MO become dominant. The influence of the initial pH value of the MB and MO solutions on the adsorption efficiency of examined samples is shown in
Figure 3.
The adsorption efficiency of non-carbonized samples increases with pH since the surfaces of these samples (pHPZC around 2) are negatively charged in the examined pH range and attract MB cations. On the other hand, these samples repel the anionic MO dye showing low adsorption efficiency. Both carbonized samples show the lowest adsorption efficiency for MB adsorption, which is slightly increased with the initial pH value. According to the pHPZC values, the surface of carbonized samples is positively charged in the examined pH range, repelling the MB cations from the surface. As the pH value of the solution approaches the pHPZC value, the adsorption efficiency slightly increases. Carbonized samples show similar behavior for MO adsorption as non-carbonized samples, having the highest adsorption efficiency in pH range 4-6. Activated samples showed the highest adsorption efficiencies for all examined pH values of solution for both adsorption of MO and MB, which is the consequence of developed specific surface area. The high adsorption efficiency of activated samples does not depend on the initial pH, and the sample Cottac completely removes both dyes from the water solution. Based on the obtained results, all following adsorption experiments were performed with the initial pH value adjusted to 8 for MB, and for MO solution without adjustment (pH 5.5).
The dynamics of MB and MO adsorption onto the surface of examined samples is given in
Figure 4, as the influence of contact time on adsorption capacities. The pseudo-first and pseudo-second order kinetic, Elovich, and intraparticle diffusion models were selected to test the adsorption dynamics and determine the kinetic parameters (
Tables S1 and S2).
For most examined samples, adsorption of both dyes can be described by pseudo-second order kinetic model, while experimental data obtained for MB adsorption onto the activated cotton sample, fit better with pseudo-first order kinetic (
Figure 4 and
Table S1). The pseudo-second and pseudo-first order models are the most commonly used models, which can fit kinetic data originating from systems limited by the surface reaction and by diffusion, thus not associated with just one adsorption mechanism [
30]. However, according to Mita et al. 2017 [
31] a good correlation of experimental data with the pseudo-second order model suggests that chemisorption is most likely the rate-limiting step in the process of adsorption. On the other hand, a good correlation with the pseudo-first order model assumes that the rate of occupation of sorption sites is proportional to the number of unoccupied sites [
32].
High values of correlation coefficients obtained for the Elovich model (
Table S2) indicate that this model can describe dye adsorption on all examined samples, except on sample Cott
ac. Applied thermal treatment of cotton and cotton/polyester yarns led to differences in surface coverage (
β), and an increase in the initial rate of adsorption (
α). Adsorption of MB and MO is a relatively fast process, reaching the equilibrium after 60 minutes for all samples, except for sample Cott
ac, which shows the fastest adsorption, with the equilibrium at 10 minutes for MB and 30 minutes for MO (
Figure 4). These findings are in agreement with the value of Elovich constant
α, which shows that the initial rate of adsorption is the highest for activated samples, especially for Cott
ac. The values of R
E approaching the equilibrium parameter based on the Elovich equation [
33], given in
Table S2, also show that dependence qt-t is mild rising (0.3 > R
E > 0.1) for almost all samples, and rapid rising (0.1 > R
E > 0.02) for adsorption of MB on carbonized and activated samples. For MB adsorption, adsorption capacities decrease in the following order: activated, unmodified yarn, and carbonized samples, while for MO adsorption unmodified and carbonized samples show similar, extremely low adsorption capacities, and the highest adsorption capacity (up to 100 mg/g) was obtained for Cott
ac.
The intraparticle diffusion model was used to evaluate the influence of diffusion on the adsorption process (
Figure 5). By applying Equation (8) on experimental data, the values of intraparticle diffusion rate constant,
kid, and constant
C were obtained and given in
Table S2.
Intraparticle diffusion plots of unmodified and activated cotton, and cotton/polyester samples are generally given as multi-linear plots of
qt vs.
t1/2, consisting of three or two (
Figure 5) consecutive steps of the process of dye adsorption. For all carbonized samples and MB adsorption on the Cott/PES
ac sample, the intraparticle diffusion plot was fitted with one straight line that does not pass through the origin. Therefore, MB and MO adsorption onto these samples is not controlled by the intraparticle diffusion step of the adsorption process. Carbonized samples do not have developed specific surface area and porosity, and the overall adsorption process, which occurs on the external surface of the material, is not controlled by intraparticle diffusion. Additionally, experimental data obtained for Cott/PES
ac (sample with a developed specific surface area and microporosity) fitted by intraparticle diffusion model gave one straight line, confirming that the process of dye adsorption occurs only on the external surface without diffusion of dye molecules into micropores of material. Similarly, the adsorption of MB and MO onto other activated samples occurs through two steps: fast adsorption onto the external surface with the highest
kid value, and equilibrium adsorption with the lowest values of
kid (
Table S2). For MB adsorption onto Cott and MO adsorption onto Cott and Cott/PES samples between the initial external adsorption and final equilibrium process, the intraparticle diffusion step is observed as a second step in the overall process of adsorption. This step is moderately fast and controlled by intraparticle diffusion of dye molecules through the macropores, cracks, and cavities of the cotton component in the yarn structure, while adsorption onto the smooth surface of the polyester component is less influenced by the intraparticle diffusion.
According to the values obtained for interparticle diffusion rate constants, kid, the initial steps of adsorption are the fastest, especially for activated samples, which is in agreement with the values of Elovich constant α. As the adsorption process proceeds, kid values decrease while values of constant C increase, and the highest are at equilibrium.
To analyze the adsorption process at the equilibrium, as well as the influence of initial dye concentration on adsorption, experimental data were fitted with Langmuir and Freundlich isotherm models (
Figure 6). Although adsorption capacity increases with the initial concentration, there is no characteristic plot on the
qe-
ce dependence for adsorption on activated samples, which indicates the lack of surface saturation in an examined concentration range, especially in the case of Cott
ac.
Isotherm parameters obtained by fitting adsorption equilibrium data with isotherm models are given in
Table 2. Correlation coefficient values (R
2) indicate that equilibrium data obtained for MB and MO adsorption fit better with the Langmuir isotherm model.
The highest values of
Q0 and
Kf (
Table 2) follow the experimentally obtained highest capacities of activated samples, which are the consequence of the developed specific surface area. However, the heterogeneity factor (1/
n) values were less than a unit, implying that the surfaces of examined materials are relatively homogeneous and that the adsorption of selected dyes is a chemical process occurring on the surface functional groups as active sites [
34].
To examine the influence of surrounding temperature on the adsorption capacities, adsorption experiments were performed at 20, 35, and 45
oC (
Figure 7).
It can be observed that adsorption capacities decrease with the temperature increase, except for MO adsorption on samples Cott/PES
c and Cott/PES
ac. This decrease in adsorption capacities as temperature increases may be the consequence of the increase in dye solubility and weakening of the physical bonds between the adsorbate and adsorbent followed by the partial removal of dye molecules from the adsorbent surface [
35,
36]. Thermodynamic parameters obtained from experimental results are shown in
Table 3. Obtained negative values for Δ
H0 and Δ
S0 for all samples, except for adsorption of MO on Cott/PES
c and Cott/PES
ac, indicate that adsorption of selected dyes is an exothermic process with decreased randomness at the solid/solute adsorption systems. According to the values of Δ
G0, adsorption of MB and MO is a feasible and spontaneous process only on samples Cott, Cott
ac, and Cott/PES
ac.
Optimal parameters of adsorption obtained in previous experiments were applied for the adsorption of methylene blue from real wastewater samples using Cott
ac and Cott/PES
ac as the most efficient adsorbents (
Figure 8). Along with the adsorption efficiency obtained for wastewater samples,
Figure 8 shows the adsorption efficiency obtained for spiked distilled water to demonstrate the effect of the matrix on adsorption. The wastewater matrix affects the adsorption of MB onto Cott/PES
ac and decreases adsorption efficiency by approximately 20 %. On the other hand, the adsorption efficiency of Cott
ac is not influenced by wastewater matrix, since this adsorbent completely removes MB from wastewater.