# Twigs-derived activated carbons via H3PO4/ZnCl2 composite activation for methylene blue and congo red dyes removal

Aug 20, 2020

### Characteristics of activated carbon

#### Proximate analysis

Figure 1a shows the thermal degradation profile of twigs (TRM). There is a decrease in weight from 100 to 92% at temperature 105 °C due to the release of moisture and dehydration of water. The existence of thermally stable constituents can create a little weight loss in material at 300 °C7. The heat resistance of twigs material is due to semi-crystalline arrays of hemicellulose, which forms strong connected chains. Also, the polysaccharide chains in lignin can create a complex structure of phenolic and aromatic groups in the cell wall structure. These combinations render a durable and heat-resistant material8. At 375 °C, the material undergone major decomposition, wherein the weight drastically decreased from 90 to 20% due to the liberation of volatiles from semi-cellulose and lignin. Semi-cellulose is relatively easy to decompose because of linear polymer structure and short side chains9. The gradual weight loss at a temperature range of 400–900 °C can be attributed to the decomposition of heavy organic materials. The carbon-rich matrix is intensely gasified into CO and CO2, leaving behind the ash content of 10%. From the TGA analysis, the fixed carbon and volatile content were recorded as 10% and 72%, respectively. Also, a suitable temperature range for activation is around 400–600 °C.

Figure 1b shows the thermal degradation profiles of impregnated samples (imTRMs). Sample imTRM2 shows a weight loss of 10% at 100 °C due to the elimination of moisture content. The estimated value of moisture content is slightly higher than that of TRM because of the presence of ZnCl2 that is hygroscopic, i.e., able to adsorb more moisture from the surrounding. The maximum weight loss of imTRM2 at 200 °C indicates the decomposition of hemicellulose and cellulose. The degradation at 500 °C could be referred to the devolatilization of lignin and residual char. A rapid degradation of imTRM2 produces a high residue of 23%, which consists of fixed carbon and zinc element. This is because, ZnCl2 could inhibit the formation of tar, hence preserving the yield10. The evaporation of ZnCl2 could occur at a temperature range of 400–600 °C11. Aromatic condensation may also take place, resulting in the evolution of gaseous products due to the presence of hydro-aromatic structure in the raw material. Sample imTRM3 shows a weight loss nearly 5% at 95 °C, which is slightly lower than that of TRM due to H3PO4 activation. A weight loss of 5–20% between 150 and 250 °C is due to the elimination of volatile matters. The gradual decomposition of light organic constituents indicates that H3PO4 restricts the volatilization of carbon, leading to a high fixed carbon of 26%12. A 40% of ash content is likely from the phosphorus residues of activating agent. The moisture content of imTRM4 is quite similar with that of TRM. The combination of chemical reaction of ZnCl2 and H3PO4 in imTRM4 produces zinc phosphate cement, Zn3(PO4)2 as shown in Eq. (3). At 100–230 °C, some hydrated zinc phosphate salts (Zn3(PO4)2·4H2O) such as hopeite formed after the evolution of cement during the water evaporation13.

$$3{mathrm{ZnCl}}_{2} (mathrm{aq}) + 2{mathrm{H}}_{3}{mathrm{PO}}_{4} (mathrm{aq})rightarrow {mathrm{Zn}}_{3}{left({mathrm{PO}}_{4}right)}_{2} (mathrm{s}) + 6mathrm{HCl }(mathrm{aq})$$

(3)

Figure 1c shows the TGA profiles of activated carbons. The thermal degradation trends of TAC2 and TAC3 are almost similar. The moisture content due to the elimination of water for TAC2 and TAC3 are 23% and 28%, respectively. TAC3 shows excellent thermal stability that can withstand up to 600 °C, as compared to TAC2 which is only up to 480 °C. H3PO4 in TAC3 promotes the reducing of weight loss of carbon due to the formation of phosphate bonds which cross-linked with each other to form the polyphosphate bridges14. Consequently, it promotes the development of porous structure at wider activation temperature. The ash content of TAC2 and TAC3 are comparable at 38% and 35%, respectively. The TAC2, TAC3 and TAC4 gave the volatile matter of 4, 2 and 6%, respectively, reflecting the high stability of activated carbons at low temperature. The formation of hopeite in TAC4 reduces the melting point so that the mass of ZnCl2 remains in char. Thus, the ash content of TAC4 is slightly higher than that of TAC2 and TAC3, which is 76%. The values of volatile matter, fixed carbon, moisture content and ash content from the TGA data are given in Table S2.

The yield of TAC1 (char) is 21.3%, which is smaller than other TAC samples because there is no chemical impregnation involved, leading to high carbon and volatiles released. The yields of TAC2 and TAC3 are 41.3% and 50%, respectively. Zinc chloride acts as a dehydrating agent to enhance the removal of volatiles such as H2 and O2 from the carbon matrix. It exhibits a high boiling point of 732 °C, and highly soluble in water which accelerates the removal of water and volatile molecules. A high amount of ZnCl2 can accelerate the evaporation of volatiles through the breakdown of aromatic and aliphatic chains. Nevertheless, ZnCl2 can preserve the carbon content via the formation of long aromatic chain structure composed primarily of graphitic structure of carbon atoms, resulting in a high yield of TAC2 as compared to that of TAC1. It comprises of networks of pores with channels constructed within a rigid skeleton of disordered layers of carbon atoms, connected by chemical bonds and stacked unevenly. ZnCl2 also acts in pore drilling that enters the carbon matrix to form large pore volume and high surface area of activated carbon, while at the same time releasing some volatiles10.

The high percentage yield of TAC3 compared to TAC2 is due to the residual phosphorus element present after the activation process because. The phosphate residue from H3PO4 activation is water insoluble. The chemical reaction between H3PO4 and lignocellulosic materials like hemicellulose and lignin during activation creates weak hydrolyze glycosidic linkages, thus accelerating the decomposition process. However, cellulose is not easily broken down because it acts as acid hydrolysis resistant material, leading to the high yield of TAC315.

#### Elemental composition

Elemental composition of a material was determined through dynamic flash combustion technique. Table 2 shows the elemental composition of raw material and activated carbons. The presence of oxygen content in TACs can be attributed to the oxidation process which occurred during the activation and after the materials are exposed at room temperature. According to16, oxygen molecules are localized at the carbon surface and form electrostatic interactions with the positively charged functional groups. A smaller fraction of oxygen in TAC2 (35%) and TAC3 (18.8%) is due to rapid removal of volatiles during the activation at high temperature, thus leaving behind reactive oxygen element.

The elements of carbon, oxygen, and hydrogen make-up the building matrix of TRM which contains lignin, cellulose, and hemicellulose. Some of these constituents are easily decomposed at high temperature. As a result, the carbon content in char TAC1 decreased to 32.9%. Through suitable activation strategies, TRM could be a potential precursor for activated carbon production. Table S3 summarizes the comparison of elemental composition of TRM with other lignocellulosic raw materials of activated carbon in literature.

Figure 2 shows the EDX mappings and spectra of activated carbons. The distribution of carbon is dominant, with scattered elements of zinc, phosphorus and sulphur. Table 3 shows the elemental composition on the external surface of activated carbons. The surface of TAC2 consists of 81.1% carbon, 4% zinc and 0.47% sulphur. The presence of zinc and sulphur could be originated from ZnCl2 used in activation and sulphur content in organic precursor, respectively. The oxidation process also prompts the formation of zinc oxide, ZnO deposited on TAC2. The formation of ZnO is based on the reaction between ZnCl2 and oxygen released from the organic material as shown in Eq. (4). The ZnO produced acts as oxygen-soluble material, but insoluble in water.

$${mathrm{ZnCl}}_{2} (mathrm{s}) + 0.5{mathrm{O}}_{2} (mathrm{g}) rightleftharpoons mathrm{ZnO }(mathrm{s}) + {mathrm{Cl}}_{2} (mathrm{g})$$

(4)

TAC3 displays a higher carbon content of 95.3%. A high carbon composition is imperative in the formation of smooth surface and well-developed porous structure of activated carbon. The presence of phosphorus is attributed to the use of H3PO4 in activation17. The sodium (0.14%) residue could be resulted from the impurity during the washing step. TAC5 exhibits a smaller carbon content than TAC2 and TAC3. The formation of zinc phosphate increases the fractions of zinc and phosphorus elements, thereby decreasing the carbon content.

#### Morphology and textural properties

Figure 3 shows the morphology of activated carbons at different magnifications. Figure 3a,b shows the SEM images of TAC2. TAC2 exhibits a surface with large pore volume. The development of pores within the matrix of TAC2 is associated to the evaporation of volatiles in the presence of dehydrating agent, ZnCl2. Nevertheless, insoluble oxides can remain on TAC2 surface due to oxidation at activation temperature that is below the boiling point of ZnCl2 of 732 °C, thus blocking some pathways10.

The smooth surface with more structured and ordered pores arrangement could be observed from the morphology of TAC3 as shown in Fig. 3c,d. Although insoluble phosphorus deposits may present, TAC3 demonstrates the smooth development of pores compared to TAC2 and TAC5. This is due to the evaporation of phosphoric acid during carbonization process, leaving behind the space formerly occupied by the activating agent.

Figure 3e,f shows the SEM images of TAC5. Based on the topography view, the insoluble Zn3(PO4)2 deposits scatter and clog the pore openings. The high melting point of zinc phosphate (900 °C) as opposed to the activation temperature at 600 °C could be the reason for its distribution over the carbon surface This situation can hamper the diffusion of molecules in adsorption due to inaccessible pores.

The yield, pH, pHpzc and textural properties of activated carbons are summarized in Table 4. Figure S1 shows the N2 adsorption–desorption isotherm and pore size distribution of TACs. The decreasing order of BET surface area is TAC3 (1,445 m2/g) > TAC2 (907 m2/g) > TAC5 (498 m2/g) > TAC6 (394 m2/g) > TAC4 (344 m2/g) > TAC1 (190 m2/g). In the absence of activation, char (TAC1) displays a smaller surface area. A relatively smaller surface area of TAC4, TAC5 and TAC6 as opposed to TAC2 and TAC3 is due to the formation of ZnO, hopeite, and parahopeite minerals that partly clog the pores.

The Langmuir surface area and total pore volume signify the pore development of activated carbon. The total pore volume of TAC3 and TAC2 are 0.730 cm3/g and 0.450 cm3/g, respectively. A high total pore volume can be attributed to the formation of pores and pathways within the matrix of activated carbon. Furthermore, H3PO4 in TAC3 preparation has two important roles, (1) to increase the rate of thermal decomposition of raw material, and (2) to promote the formation of cross-linked bonding18. Table S4 shows the comparison between TAC3 and other activated carbons from the viewpoint of textural properties.

The adsorption behaviour of activated carbon depends not only on its inner surface area but also the pore size. TAC2, TAC3 and TAC4 are microporous with pore width of 1.98 nm, 2.0 nm 1.96 nm, respectively. Meanwhile, TAC1, TAC5 and TAC6 are mesoporous with pore width of 4.13 nm, 2.61 nm and 2.77 nm, respectively. It shows that H3PO4 instigates the development of micropores in TAC3 and TAC4 and mesopores in TAC5 and TAC6. The underdeveloped porous texture of TAC1 indicates that it is not suitable for adsorption because of the restricted access and limited active sites, even though it is highly mesoporous.

#### Surface chemistry and functional groups

The FTIR spectra of raw material and activated carbons are shown in Fig. 4. The peak at 3,850–3,700 cm−1 is due to physisorbed H2O molecules19. All samples display identical band of medium intensity around 3,500 cm−1 which is attributed to the vibration of hydroxyl groups (OH) such as alcohol and phenol due to hydrogen bonding from water absorption, and primary amines (NH) stretching20. The presence of hydroxyl groups is associated with the lignocellulosic structures in the parent material. The band of 2,370 cm−1 indicates the stretching of carboxylic groups (COOH). However, it disappeared in TAC2 spectrum because of strong dehydrating effect of ZnCl2 during activation.

The alkene (C=C) stretching vibration is the characteristic of peak at 1,640 cm−1 in TRM spectrum. The vibration probably occurred in olefinic and acyclic compounds, and conjugates. The profile is closely related to lignocellulosic-rich mustard card as reported in21. Yet, the peak is obviously diminished in all heated samples due to the conversion of aliphatic long chain to graphene (aromatic) sheets. The multiple bands between 1,500 and 1,400 cm−1 could be assigned to C=C stretching in aromatic groups or benzene rings and the bending vibration of methylene (alkane) groups.

The peaks at wavenumber below 650–450 cm−1 are generally attributed to C-X bonds where X could be C, O and N atoms for volatile matters such as sulphur, silicone, and halogenated compounds. The weak peaks around 950–800 cm−1 signify the CH vibration of alkenes and compounds derived from benzene that is associated with the out-of-plane deformation from the bending of CH and hydroxyl groups. The strong intensity and broad peak around 1,100–1,010 cm−1 in all spectra indicates strong absorption of CO groups of primary alcohols (R-OH), phenol groups and asymmetric vibration of Si–O–Si. Specifically, the peaks at 950–800 cm−1 are related to asymmetric vibration of Si–OH and symmetric stretching vibration of Si–O–Si.

Table 5 shows the concentrations of acidic and basic groups in activated carbons by Boehm titration. According to22, the basic nature is mainly due to oxygen-containing groups like carbonyl, ether, diketone and chromene, delocalized of π-electrons of carbon basal planes, and inorganic minerals. While, the acidic nature is due to the presence of carboxylic, lactonic and phenolic groups.

TAC3 shows a total acidic concentration of 1.15 mmol/g. TAC4 and TAC5 exhibit the absence of basic sites, while TAC1 endows a 0.761 mmol/g of basic groups. The presence of acidic and basic groups in activated carbons are in accordance with the peaks in FTIR spectra. Basically, the density of acidic and basic sites is due to the nature of Pterocarpus indicus twigs and chemical activation strategies. The presence of surface functional groups is also the contributing factor in adsorption performance. The physicochemical properties of raw material are transformed due to the use activating agents. The increasing pattern of total acidic sites (TAC3 > TAC5 > TAC4 > TAC6) is reflected from the increase of H3PO4 concentration during impregnation. Besides, the concentration of carboxylic groups in TAC3 is the highest among the samples at 0.70 mmol/g. The absence of lactonic groups in TAC2, TAC3, TAC5 and TAC6 indicates that the entire of lactonic groups in twigs have been converted into carboxylic and phenolic groups during chemical activation19. TAC1 shows higher basic sites than the other TACs because of the inherent nature of the plant material. Also, there is no chemical agent involved to promote CO2 gasification during twigs carbonization at 600 °C. The presence of basic sites in TAC3 and TAC6 is presumably governed in part by anionic phosphate minerals23.

Figure S2 shows the pH drift method to determine pHpzc of activated carbons. The interception of coloured lines with the 45° diagonal line represent the pHpzc of activated carbons. The respective values are summarized in Table 4. The decreasing pattern of pHpzc is often associated with the increasing concentration of acidic groups on the adsorbent surface. TAC3 shows a pHpzc of 3.6, while TAC1 exhibits a basic-rich surface with pHpzc of 7.9. The pHpzc of TAC2 and TAC6 are slightly acidic, and close to neutral (pH 7.0), which correlates well with the Boehm titration data. The high value of pHpzc indicates that the surface is negatively charged and may favour the removal of positively charged molecules19.

Figure 5a illustrates the equilibrium adsorption (Qe) of methylene blue onto the activated carbons. For all activated carbons, the equilibrium adsorption increases with increasing the initial concentration (Ce) until the saturation point is reached. At equilibrium, the adsorption rate is equal to desorption rate. TAC2 and TAC3 show excellent removal performance at 401 mg/g and 413 mg/g, respectively. The two TACs possess high surface area to encourage interaction probabilities with methylene blue molecules. A slightly bigger methylene blue capacity by TAC3 is due to the high degree of carboxylic groups and total acidic sites on the carbon surface.

The high adsorption capacity signifies the smooth diffusion of methylene blue molecules from bulk solution to the surface of activated carbon. A steep and high adsorption capacity at low concentration indicates a stronger carbon affinity towards dye molecules. Char (TAC1) shows the smallest methylene blue capacity of 16.9 mg/g, which agrees with its inferior surface area as compared with other TACs. The methylene blue capacities of TAC4, TAC5 and TAC6 are 130, 264 and 221 mg/g, respectively. TAC5 exhibits a slightly greater capacity than TAC6 because of higher surface area and smaller pHpzc. This signifies the primary role of surface area (pore volume) and pHpzc in methylene blue adsorption by TACs. Generally, the positively charged dye shows a stronger electrostatic attraction with the negatively charged carbon surface, which correlates well with the pHpzc and Boehm titration data.

Figure 5b shows the equilibrium adsorption of congo red onto activated carbons. A small adsorption indicates that the adsorption is concentration-dependent due to limited number of dye molecules in solution at low concentration. The highest congo red capacity is 214 mg/g (TAC5), which is about half that of methylene blue (in mass basis). TAC6 endows a 190 mg/g capacity of congo red. It implies the role of zinc phosphate complex in the performance of TAC5 and TAC6. The formation of complex can be represented as ligand that forms strong covalent bonds with anionic molecules24. A better performance of TAC5 than that of TAC6 is due to high H3PO4 concentration used in activation that results in a higher surface area of 498 m2/g.

The negatively charged surface of TAC3 creates electrostatic repulsion that weakens the dispersion forces, hence decreasing the adsorption capacity25. Nevertheless, the removal performance of TAC2 (161.1 mg/g) and TAC3 (176.2 mg/g) can still be considered high when compared to that of TAC1 and TAC4. This can be attributed to the high surface area and pore volume of the former. Furthermore, the high adsorption capacity of TAC3 is due to the presence of H3PO4 that increases the chemisorbed oxygen on the carbon surface26 . TAC1 shows a 90.3 mg/g capacity of congo red. TAC4 also possesses a low adsorption capacity (119.1 mg/g) than other treated activated carbons due to poor pore development.

Table 6 show the isotherm constants for methylene blue and congo red adsorption by activated carbons. The isotherm parameters were solved by non-linear regression to provide a mathematically rigorous method from the original equation. The adsorption capacity at molecular level was evaluated in molar basis. TAC3 exhibits a methylene blue capacity of 1.29 mmol/g. From the perspective of molar basis, methylene blue adsorption is greater than congo red adsorption. This could be governed by the smaller molecular size of the former (0.95 nm)27 than the latter (2.30 nm)28. Thus, methylene blue molecules could easily enter and diffuse through the pore channels of activated carbon. The large molecular size of congo red affect the performance of TACs, whereby the dye molecules can only interact with the mesoporous surface of activated carbon.

From Table 6, the equilibrium data fitted well with Langmuir model with R2 in the range of 0.9031 to 0.9978, and 0.9960 to 0.9996 for methylene blue and congo red adsorption, respectively. It can be concluded that the adsorption forms monolayer of molecules onto homogenous carbon surface19. It is supported by the smaller sum-of-squared error (SSE) values. The predicted capacities (Qm) of 438 mg/g (TAC3) and 217 mg/g (TAC5) for methylene blue and congo red, respectively are closely tallied with the experimental values.

The Langmuir constant, b indicates the affinity of carbon active sites towards dye molecules. From Table 6, TAC1 displays a high b value (0.521 L/mg) suggesting a high methylene blue adsorption at low concentration probably due to basic surface that instigates electrostatic interaction. The low b but high Qe reflects the other contributing factors in adsorption. On the contrary, TAC1 shows a small b value for congo red removal. The high constant b for congo red adsorption is credited to the chemical modification which enhances the textural characteristics of activated carbon. However, the predominance of carboxylic and phenolic groups on carbon surface can create a high degree of H+ dissociation towards the repulsion of anionic molecules.

The favourability of dye adsorption was predicted by separation factor, RL a feature of Langmuir isotherm as shown in Fig. S3. All adsorption data obeyed 0 < RL < 1, implying that the process is favourable. While, RL > 1, RL = 1, and RL = 0 are attributed to unfavourable, linear, and irreversible adsorption, respectively29.

Freundlich model showed the least fitting model for equilibrium adsorption data of methylene blue and congo red by TACs. However, the 1/n value lesser than unity agreed well with Langmuir model, while 1/n > 1 reflects the adsorption activity between readily adsorbed molecules and free molecules in bulk solution. The Redlich–Peterson model is a hybrid of Langmuir and Freundlich models30. At g < 1, the large molecules normally restrict the adsorption due to pore clogging, while g = 1 simplifies the model into Langmuir equation. A linear Henry’s Law is when g = 0. From Table 6, the K values are in the range of 15.5–66.1 (methylene blue) and 2.31–11 (congo red). The Redlich–Peterson model supports Freundlich equation when 1/K → 0. Conversely, 1/K ≠ 0 signifies the dominance of monolayer coverage of dye molecules in adsorption. Table S5 summarizes the comparison of methylene blue and congo red adsorption by wood-based activated carbons. TACs demonstrate comparable and substantial performance, although the surface area (pore volume) may not the sole contributing factor in dyes adsorption.

Figure 6 illustrates the kinetics of methylene blue and congo red adsorption by TACs. Generally, all curves exhibit a similar convex upward shape of Qt against t but at varying magnitudes at equilibrium. The capacity rapidly increases at the beginning due to plentiful vacant sites on the carbon surface that renders concentration gradient as driving force for adsorption. Thereafter, the adsorption rate subsides because the sites are gradually occupied when approaching equilibrium. Desorption may also take place between the already adsorbed molecules and the free-moving molecules in bulk at prolonged contact time.

From Fig. 6a, TAC2 and TAC3 demonstrate a rapid equilibrium of 3 h, followed by TAC5 and TAC6 (4 h), and TAC4 (5 h) for methylene blue concentration of 20 mg/L. The time taken to attain equilibrium when in contact with 100 mg/L dye solution increases to about 6 h for TAC2 and TAC3, as displayed in Fig. 6b. In Fig. 6c, TAC5 and TAC6 require 5 h to achieve equilibrium for adsorption in 20 mg/L congo red, while a longer contact time was recorded for TAC, TAC3 and TAC4. In Fig. 6d, a 9 h contact time was needed for TAC5 and TAC6 to achieve equilibrium in 100 mg/L congo red. A steep adsorption rate of methylene blue by TAC2 and TAC3 at short contact time is probably due to instantaneous electrostatic interaction of dye molecules with surface acidic functional groups. This effect is more prevalent because methylene blue molecules can easily diffuse through the pores within the carbon matrix.

Table 7 summarizes the kinetics constants of methylene blue and congo red adsorption by TACs. The models demonstrate the rate-limiting step in adsorption. From Table 7, the pseudo-second-order model shows a good correlation with the kinetics data. The close agreement between the calculated Qe and experimental value shows the applicability of the model to explain the mechanisms of dyes adsorption by TACs. Dyes removal can be described as chemisorption through the exchange of electrons or electrons sharing between the dye molecules and activated carbon31. Additionally, the rate-limiting step for dye adsorption is surface adsorption.

The pseudo-first-order equation did not obey the kinetics data for the entire range of contact time but may only applicable for the initial adsorption stage31. A high dye concentration induces stronger molecular repulsion and competition for active sites. Consequently, the rate constant, k2 and adsorption rate decreased, resulting in a longer contact time to reach equilibrium. A high k2 value is also associated with the developed porous structure of activated carbon.

Table 8 summarizes the kinetics constants for intraparticle diffusion (Webber-Morris) and Boyd’s models to elucidate the adsorption diffusion mechanism. The adsorption is often controlled by a sequence of four stages: (1) Bulk diffusion—molecules diffuse onto the activated carbon surface from the bulk solution; (2) Film diffusion—molecules diffuse through the boundary layer; (3) Intraparticle diffusion—molecules diffuse from the film surface into the pores of activated carbon; (4) Adsorption—active sites of activated carbon is occupied by molecules32. The intercept of intraparticle diffusion plot, C signifies that the pore or intraparticle diffusion is not the only rate-limiting step in adsorption. It indicates the boundary layer (film) thickness. A greater C implies a thicker film, or film diffusion as dominant rate-limiting step in adsorption19.

Figure S4 shows the intraparticle diffusion curves for methylene blue and congo red removal by TACs at different concentrations. The first region indicates a strong electrostatic interaction between dye molecules and carbon surface; the second region is where the adsorption takes place via molecular diffusion in pores; and the last region signifies the diffusion rate decelerates to a point of equilibrium. In Boyd’s model, the rate-limiting step is attributed to pore diffusion. Effective diffusion coefficient is expressed as, B = πDi/r2, where, B is a linear line gradient from Bt versus t (s), Di (cm2/s) is the effective diffusion coefficient, and r (cm) is the activated carbon radius. According toss33, the intraparticle diffusion is the rate-limiting step for Di in the range of 10–6–10–12 cm2/s. Hence, all TACs studied possess the same rate-limiting step in the adsorption, which is intraparticle diffusion.

Figure S5 shows the equilibrium adsorption for (a) methylene blue and (b) congo red onto activated carbons at 30–60 °C. Dye molecules easily diffused from bulk solution into the internal pores at high temperature due to thin film thickness (intraparticle transport pore diffusion)34. The adsorption capacity of dyes studied increased with increasing temperature, implying an endothermic nature of adsorption. At high temperature, the solution becomes less viscous and dye molecules gained more kinetics energy to vibrate and surpass the mass transfer resistance, thus accelerating the migration of dye molecules towards the internal pores of carbon matrix35. For high dye concentration, the increase in solution temperature aids in weakening the retarding forces against molecular migration. The mobility of dye molecules stimulates more interaction probabilities with carbon sites.

Figure S6 shows the van’t Hoff plots for (a) methylene blue (b) congo red. Table 9 summarizes the thermodynamics parameters of methylene blue and congo red adsorption by activated carbons. The Gibb’s free energy, ΔG° becomes more negative at high temperature, suggesting a feasible and spontaneous adsorption process19. The positive entropy, ΔS° describes the orientation disorder and randomness between dye molecules and carbon surface36. From Table 9, the ΔS° values are in the range of 112–333 J/mol K (methylene blue) and 66.3–470 J/mol K (congo red). At higher concentration, the dye solution becomes more dispersed, rendering the molecules more freedom of arrangement, hence increasing the entropy.

The possible mechanisms of methylene blue adsorption by TAC3 consist of ion-exchange, hydrogen bonding, π-π conjugation, electrostatic attraction, and physical adsorption. The π–π conjugation is a lateral attraction between aromatic ring of dye molecule and graphitic (aromatic) layer of carbon37, while physisorption is generally tied-up with van der Waals forces because of the well-developed porous structure of activated carbon. From the equilibrium and kinetics data, TAC3 exhibits chemisorption by electron withdrawn containing π-electron acceptors such as ketone, aldehyde, and benzene groups. In addition, N–H group of methylene blue can establish hydrogen bonding with carboxylic (-COOH) and hydroxyl groups (-OH) in TAC3 surface. Also, the carboxylic groups dissociate the protons and become negatively charged to offer electrostatic attraction with methylene blue molecules. Accordingly, the solution pH slightly decreased due to the liberation of H+ ions38. Electrostatic interaction also exists between acidic carbon sites and basic methylene blue dye when the solution is higher than pHpzc, at which the surface is surrounded by excess OH. Besides, the Si–OH group present in TAC3 possibly prompts the n-π interaction with methylene blue. The presence of Si–OH group in TAC3 (Fig. 4) enhances the removal of methylene blue in basic environment39.

Complex interaction of Zn3(PO4)2·4H2O with congo red is the main mechanism for TAC5. The hydrated salt that holds a positively charged Zn2+ allows electrostatic interaction with anionic dye. According to40, ion–dipole interaction may take place from the interaction of Zn2+ with polar molecules (partial negatively charged oxygen in congo red). Additionally, amine group in congo red can also form hydrogen bonding with phosphate ion (PO43−). These two interactions are attributed to hydrophilic adsorption. Similar mechanism has been proposed for congo red removal by calcium hydroxyapatite41.