3.1. Properties of different sized SARC
The SEM images of S-SARC, M-SARC, and L-SARC are shown in Fig. 1 (c), (d), and (e), which exhibit different surface morphology after NaOH modification. L-SARC samples have uneven surfaces with unique equiaxed pores in cross directions, thus exponentially amplifying their adsorption capacity. Their internal pore structure is clearer and orderly, and the surface pores are more open with thinner walls. These thin walls may increase the surface area to some degree, however, making them vulnerable to collapse. As shown in Fig. 1 (e), a certain number of destructed lamellar pores are visible upon the corrosion of NaOH. In some areas, adjacent macropores are merging into super large pores. One of the possible reactions may be attributed to NaOH corrosion was proposed by Tryk et al. 7:
\(\text{C}+6{\text{O}\text{H}}^{-}= {\text{C}\text{O}}_{3}^{2-}+3{\text{H}}_{2}\text{O}+4{\text{e}}^{-} ; {{\Delta }\text{G}}^{^\circ }=-70.72 \text{k}\text{c}\text{a}\text{l}/\text{m}\text{o}\text{l}\text{e}\) Eq. (3)
It is also known that:
\(4{\text{H}}^{+}+{\text{O}}_{2}+4{\text{e}}^{-}=2{\text{H}}_{2}\text{O}; {{\Delta }\text{G}}^{^\circ }=-113.38 \text{k}\text{c}\text{a}\text{l}/\text{m}\text{o}\text{l}\text{e}\) Eq. (4)
Combining both equations:
\(\text{C}+2{\text{O}\text{H}}^{-}+{\text{O}}_{2}= {\text{C}\text{O}}_{3}^{2-}+{\text{H}}_{2}\text{O} ; {{\Delta }\text{G}}^{^\circ }=-184.1 \text{k}\text{c}\text{a}\text{l}/\text{m}\text{o}\text{l}\text{e}\) Eq. (5)
The OH− ions from NaOH are consumed during alkali modification, and react with carbon with the aid of oxygen. This process generates ionized carbon (CO32−), resulting in the collapse of carbon walls meanwhile generating more internal pores for L-SARC. On the other hand, the medium and small sized SARC particles (M-SARC and S-SARC) have relatively thicker pore walls, which can resist the excessive corrosion collapse caused by NaOH solution. As shown in Fig. 1 (d), M-SARC particles possess well defined micropattern in their porous channels. The honeycomb structure is visible among the larger M-SARC particles. Most M-SARC particles are crosslinked with each other, possibly due to the strong interactions among their remaining acidic groups. These groups are mostly oxygen-containing groups, e.g. carboxyl, lactone and hydroxyl groups, which can generate hydrogen bonding or electrostatic interactions between each other. While for S-SARC samples (Fig. 1 (c), their sizes are much smaller than those of the others, indicating that S-SARC consist of more micropores. These pores are composed of crystalline structure, usually featuring internal cavities of variable shapes and diameters which are able to adsorb VOCs molecules via capillary and partitioning mechanisms 8. The individual particles are evenly dispersed inside of the domain and have very minimum interference between particles. As shown in the size distribution histogram (Fig. 1 (f)), L-SARC particles are relatively uniform in their size with a narrow distribution, and dav is about 175.1 µm. In contrast, M-SARC has a heterogenous size distribution with an average diameter of 80.5 µm. The majority of S-SARC particles have fairly homogeneous texture and size distribution. The average diameter of S-SARC is 28.6 µm, which is the smallest among the three SARC samples.
The micro- and meso-pore size distributions of three SARC samples were determined with BJH model. As shown in Fig. 1 (a), the pore size distribution was similar for all samples and relatively intensive. The average pore diameter increased as the size of SARC increased, and the peak values were 1.71 nm, 1.93 nm and 2.45 nm for S-SARC, M-SARC and L-SARC, respectively. According to the IUPAC (International Union of Pure and Applied Chemistry) classification, both S-SARC and M-SARC were microporous (< 2 nm), while L-SARC were mesoporous (2–50 nm)9.
With high resolution energy-dispersive X-ray spectroscopy (EDS), the elemental compositions of different sized SARC were determined and the analysis results are shown in Fig. 1 (b). Carbon (C) and oxygen (O) are the main elements for all three samples with total concentration of more than 99% (99.17% for S-SARC; 99.21 for M-SARC; 99.22 for L-SARC). As the size of SARC increases, the carbon content decreases, while the oxygen content increases. The majority of oxygen elements are sourced from acidic functional groups, and ratio of O/C can be used to characterize the degree of polarity and hydrophobicity. The high O/C ratio may result in a weakened affinity towards organic VOCs, thus negatively affecting the adsorption capacity. Alkali modification is able to reduce the oxygen-containing groups. As shown in Figure. 1(b), the smaller the diameter of SARC, the less the number of oxygen element and the ratio of O/C. This is because the smaller sized SARC particles are more accessible by NaOH solution due to their larger specific area, resulting in more acidic groups being modified.
3.2. Molecular dynamic simulation
Particle size, as an important factor, can significantly influence the migration of SARC particles inside of LCCs, thus affecting their adsorption towards VOCs molecules. Figure 2 (a-c) shows the molecular models of LCCs matrix consisting of different sized SARC particles in PP domain. The figures (a), (b), and (c) correspond to S-SARC, M-SARC, and L-SARC, respectively. Due to the arbitrary scales of MD simulation, magnitudes of these three SARC samples were relative. It can be observed that as the size decreases, the distribution of SARC becomes more homogeneous, suggesting better contact between SARC particles and polymer chains.
The migration behavior of SARC particles inside the domain was simulated by LAMMPS-VMD software. In the user-defined coordinate system, the elapsed time (fs) and the x- and y-coordinates of the individual particles were recorded for the calculation of particles’ migration velocity:
\(v=\frac{\sqrt{{\left({x}_{s}-{x}_{i}\right)}^{2}+{\left({y}_{s}-{y}_{i}\right)}^{2}}}{{t}_{s}-{t}_{i}}\) Eq. (6)
where v is the particles’ velocity, xs is the particle’s x-coordinate in the stopping position, xi is the particle’s x-coordinate in the initial position, ys is the particle’s y-coordinate in the in the stopping position, yi is the particle’s y-coordinate in the initial position, ts is the time at the stopping position, and ti is the time at the initial position.
Figure 2 (d), (e), and (f) are the migration profiles of S-SARC, M-SARC, and L-SARC, respectively, in 3D view (Ortho, Front, Left and Perspective). The flow of these particles may be driven by shear stress, Coulomb friction as well as inter-particle electrostatic interactions. The relaxed confining stress present at the surface of particles would allow them to readily arrange during compounding. The results were plotted in terms of particles’ velocity and elapse time. There was a significant decrease in particles’ velocity as their size increased, and the M-SARC and L-SARC particles migrated with velocities of only 10–40% of those of S-SARC particles. This may be because S-SARC meandered more easily through the fortuitous channels inside of the domain, thus having higher velocities. Moreover, the large fluctuation of velocity trajectory indicates S-SARC particles’ travel may be impeded by frictional interactions between particles, thus migrating in multiple modes such as shear flow, rolling and sliding with different speed. In contrast, a relative narrow range of velocities was observed for both M-SARC and L-SARC. This may be attributed to their larger size and the viscous PP matrix, which would have fewer particle-particle interactions to affect their rearrangement. The higher migration velocity of particles allows them to contact the domain more frequently and distribute homogeneously, thus being more effective in capturing VOCs.
3.3. Calculation of adsorption energy
Shown in Fig. 3 (a), (b), and (c) are the adsorption isotherm of N2 on the three clusters of SARC (insert), which can be categorized as Type I (Langmuir isotherm), where nitrogen molecules are adsorbed physically by micropore filling, according to IUPAC classification (IUPAC, 1985). The adsorption energy is further calculated based on the isotherm profiles. As is known, adsorption is driven by the adsorption force field on the surface of SARC. This unique force field generates an adsorption potential that is determined by the adsorbate, the adsorbent, and the distance between the adsorbate molecules and adsorbent surface. The Dubinin-Radushkevich (D-R) equation is able to adequately characterize these parameters, and the adsorption energy of SARC can be calculated by D-R equation via linearly fitting the adsorption isotherm. The D-R equation can be expressed as follows 10,11:
\(W={W}_{0}exp{\left[\frac{-A}{{\beta E}_{0}}\right]}^{2}={W}_{0}exp{\left[\frac{-A}{E}\right]}^{2}\) Eq. (7)
\(lnW=ln{W}_{0}-2.303{\left(\frac{RT}{E}\right)}^{2}\) Eq. (8)
where W is the equilibrium adsorption amount under P/P0, mL/g; W0 is the micropore adsorbed volume, mL/g; A is the adsorption potential, which equals to RTln (P0/P); β is an affinity coefficient; E0 and E are the normal N2 adsorption energy and characteristic adsorption energy, respectively; P0 and P are the saturation vapor pressure and equilibrium pressure, respectively; R is the conventional gas constant; and T is the absolute temperature.
The linear correlation of the N2 adsorption isotherm for the three clusters of SARC was determined according to the D-R equation, and the relationship between ln2(P0/P) and lnW was obtained. The slope of the fitting curve in Fig. 3 (a-c) is 2.303(RT/βE0)2, and the affinity coefficient of nitrogen is 0.33. Thus, the adsorption energies of S-SARC, M-SARC, and L-SARC were calculated as 25.29 kJ/mol, 24.80 kJ/mol, and 23.98 kJ/mol, respectively.
Acidity and basicity of three SARC samples are presented in Fig. 3 (d). Boehm titration confirmed the lowering of total acidic groups as the size of SARC decreased, with values of 0.178 mmol/g, 0.165mmol/g, and 0.125 mmol/g for L-SARC, M-SARC, and S-SARC, respectively. This result is in good agreement with the trend of oxygen element and the ratio of O/C in the previous section of 3.1. Properties of different sized SARC. While for total alkaline groups, an opposite trend is presented, with S-SARC having the highest alkaline groups on its surface, followed by M-SARC and L-SARC.
3.4. Identification and quantification of VOCs by GC/MS
GC/MS spectroscopy was applied for measuring the odorous VOCs. Both lignocellulosic fibers and PP can emit VOCs due to the thermal degradation of lignocellulosic fibers, and the reaction of PP with oxygen present in the extruder barrel. However, at the temperature range for LCCs compounding (180–220°C), lignocellulosic fibers are the main producers for VOCs. Figure 4 presents the various emission levels of VOCs from LCCs embedded with S-SARC, M-SARC, and L-SARC, respectively. A complicated composition of VOCs can be observed, which includes more than 120 compounds. Most of these compounds possess unpleasant odors, and some of them are carcinogens. Therefore, it is necessary to adsorb the volatile compounds and prevent them emitting into the environment. Worth mentioning, these volatiles can also be useful for other purposes, e.g., pesticide or pharmacology, upon proper collection and treatment. Moreover, as the size of SARC decreases, both the numbers and intensity counts of the GC peaks - which stand for the compounds escaped from the adsorption, are reduced, indicating that the smaller sized SARC particles are more efficient in the capture of VOCs.
MS analysis was applied to identify the GC peaks, and the detected compounds were categorized in Table 2. The concentration of compounds was calculated based on the integrated chromatographic area of GC peaks. The main types of VOCs were phenols, aldehydes, ketones, acid derivatives, furfural and furan derivatives, hydrocarbons, and others. Among which, the most abundant compounds were acid derivatives. These acids were mainly from fiber-pyrolysates, and their high content indicated that a large amount of cellulose and hemicellulose were degraded during the LCCs compounding. Similarly, the furfural and furan derivatives were produced from lignocellulosic fibers which were rich in pentosan, and phenolic volatiles were generated from the pyrogenation of lignin. Aldehyde, ketones, and other molecules from wood-resin, were among the least emitted volatiles, however, usually had the strongest odor. For these identified as hydrocarbons, e.g., Toluene and Benzene, the minority of them were from PP, while the majority ones were generated by lignocellulosic fibers. This result is in good agreement with previous studies which concluded that wood sawdust was the main pollution source because the unpleasant odor was mostly generated from sawdust rather than from PP during LCCs manufacturing 12.
Table 2
GC/MS analysis of VOCs and their chemical nature.
Compounds | Contents ((µg/L) | Acidity | Hydrophobicity | Polarity |
S-SARC | M-MARC | L-SARC |
Phenols | 21.09 | 28.57 | 34.44 | Weak acid | Mostly hydrophilic | Polar |
Aldehydes | 12.15 | 11.41 | 12.55 | Neutral | Hydrophilic | Polar |
Ketones | 13.37 | 11.83 | 15.31 | Neutral | Hydro-neutral | Polar |
Acid derivatives | 36.25 | 41.65 | 42.99 | Acid | Hydrophilic | Polar |
Furfural and furan derivatives | 19.70 | 22.08 | 22.25 | Weak base | Hydrophobic | Polar |
Hydrocarbons | 16.98 | 18.01 | 17.40 | Neutral | Hydrophobic | Nonpolar |
Others* | 15.47 | 8.37 | 2.15 | Mixed | Mixed | Mixed |
*Others include all detected minor compounds with or without odor.
The detected compounds were further classified for their acidity, hydrophobicity and polarity. The emitted VOCs of S-SARC embedded LCCs had the lowest amount of acidic, hydrophilic, and polar compounds (about 42%, 51%, and 61%), followed by M-SARC (about 49%, 58%, and 66%), and L-SARC (about 53%, 61%, and 72%). This indicated S-SARC had a good adsorption capacity towards acidic, hydrophilic and polar VOCs. The reason is, as discussed in section 3.1. Properties of different sized SARC, S-SARC were more easily to be infiltrated by alkali (NaOH) and had more alkaline groups being introduced, thus exhibiting higher affinity towards VOCs with opposite nature. Similar results were reported in the previous studies upon investigating the influence of chemical properties of VOCs on their adsorption onto carbons with different surface functionalities 13,14. The modification which weakened the surface acidity of carbons was able to enhance their hydrophobicity and affinity with VOCs of various polarity 15.
The Langmuir isotherm is valid for the adsorption with uniform energies on a surface containing a finite number of identical sites, and only one molecule can be adsorbed at each site. The adsorption fitted Langmuir isotherm is considered as monolayer adsorption with no transmigration of adsorbates in the plane of the adsorbent surface. Figure 5 (a) shows the fitting curve between Ce/qe versus Ce, with a slope of 1/Qm and an intercept of 1/KLQm. As the particle size of the samples decreased, the maximum coverage capacity Qm increased. Among the three samples, S-SARC had the best fitting of R2 = 0.8847, and the highest VOCs adsorption with Qm=0.058 mg/g. Similar influences of adsorbents’ particle size on their adsorption capacity have been observed previously, e.g., Alkherraz et al. studied the biosorption behavior of activated carbon towards metal ions including Pb(II), Zn(II), Cu(II), and Cd(II). Their activated carbon was prepared from olive branches with different particle sizes of 125, 500 and 800 µm, and the results showed an increase in the uptake of metal ions as the size of activated carbon decreased 16.
To compare, the adsorption data was further analyzed by Freundlich isotherm, which is applicable for multilayer adsorption on heterogeneous sites. A linear Freundlich isotherm model was applied by plotting LogQe against LogCe, due to the fact that adsorptions isotherm at low adsorbate concentration, e.g., VOCs emitted from LCCs, are often linear. In comparison with Langmuir model, regression analysis of the data from all three samples fitted Freundlich adsorption isotherm better (Fig. 5 (b)), with R2 of 0.9919, 0.9694, and 0.9253 for S-SARC, M-SARC and L-SARC, respectively. The intercept of the fitting curve is log(-Kf) and the slope is 1/n. The constant Kf is an approximate indicator of adsorption affinity for VOCs. The larger the Kf, the higher the maximum adsorption capacity. As for n, its value is a function of adsorption strength, and a favorable adsorption tends to have a n between 1 and 10. Of which, n = 1 implies a linear adsorption leading to identical adsorption energies for all sites, while the larger value of n (smaller value of 1/n) corresponds to a more heterogeneous adsorption, which usually results in an increased non-linearity. As shown in Fig. 5 (b), all SARC samples have n > 1; meanwhile S-SARC has the minimum n = 1.491 and highest Kf=0.106 mg/g, implies its strong linear Freundlich adsorption towards VOCs.
3.6. Influence of physical and chemical properties on VOCs removal capacity
The alkali (NaOH) modification provided renewable carbon with high porous structure as well as rich surface functionalities including both acidic and alkaline groups, e.g., oxygen-containing, nitrogen-containing, and aromatic groups. These unique physical and chemical properties allowed the prepared SARC particles adsorbing various VOCs as well as thermally degraded gases during LCCs manufacturing.
Figure 6 (a-b) presents the physical properties of SARC samples and their influence on VOCs removal capacity. The NaOH solution erodes the carbon surface, and the reaction between carbon and NaOH may form more pores from micro-size to macro-size (as shown in Eq. (3), (4) and (5)). These newly generated pores can increase the BET surface area (SBET) and migration channels for VOCs. Specifically, S-SARC had the highest SBET after NaOH modification, followed by M-SARC and L-SARC. The high SBET was proposed to have stronger adsorption towards adsorbates due to higher surface energy in previous work 17,18, which is in good agreement with the results in this study. Similarly, the total pore volume and micropore volume have the following order: S-SARC>M-SARC>L-SARC. Moreover, S-SARC has the maximum ratio of micropore/total pore volume (78.41%). There is a general agreement that micropores are favorable in the adsorption of small guest molecules, e.g., VOCs, while meso- and macropores provide fast and unobstructed channels for transporting the VOCs to the active sites in micropores 19. As the diameter of SARC increases, the ratio of micropore/pore volume decreases, which negatively affects the adsorption of VOCs as shown in Figure 6 (b).
Linear regressions were also performed on the chemical properties versus removal capacity. As shown in Fig. 6 (d), the effect of total alkaline groups on the removal capacity demonstrated a similar positive trend with that of physical properties (SBET and pore/micropore volume), while the total acidic groups presented an opposite trend (Fig. 6 (c)). Compared with physical properties (R2 = 0.9746 for SBET, 0.9796 for total pore, and 0.9921 for micropore), the adsorption capacity of SARC was more sensitive to the changes in its surface functionalities (R2 = 0.9845 for total acidic, and 0.9977 for total alkaline), indicating its superior chemical adsorption (Freundlich Isotherm). Similar results were report with other carbon-based materials, e.g. Lin et al. found that activated-carbon fiber (ACF) with less surface oxygen groups (< 900 µmol/g) tended to adsorb more nonpolar VOCs (benzene, toluene) rather than polar VOCs (acetaldehyde, acetone) 20. Yu et al. compared the adsorption behavior of GO and rGO towards benzene and toluene at room temperature. Their results showed rGO exhibited higher adsorption capacities of benzene and toluene (276.4 and 304.4 mg/g) than that of GO (216.2 and 240.6 mg/g), possibly due to the more hydrophobic nature of rGO after the plentiful oxygen groups were removed, which favored the adsorption of nonpolar or weak polar VOCs 21. Li et al. compared the performance of various acids and bases for their surface reduction modification of coconut shell-derived granular activated carbon (AC), in order to improve the adsorption of hydrophobic carbon-based gases. Their results showed that alkali modification increased the specific surface areas and pore volumes of AC, and decreased its total content of oxygen functional groups. The modified AC exhibited a good adsorption towards o-xylene; in contrast, the acid-modified AC exhibited a weaker adsorption as its total surface oxygen-containing groups increased 22.
Also noticed quite a few literatures reported an opposite result that the acid-modified carbon had higher VOCs adsorption, which was reasonable considering the polarity, hydrophilicity and the limited types of compounds in their investigated system 23. Moreover, there existed competitive adsorptions, in which weakly affinitive VOCs would be displaced by the strongly affinitive ones, though they had similar volatility and chemical nature 24. In the present study, the VOCs from LCCs were a wide variety of organic mixtures (as shown in section 3.4. Identification and quantification of VOCs by GC/MS), and alkali treatment was demonstrated as a more effective approach in modifying SARC to enhance its removal capacity of VOCs.
The synergistic effect of physical characteristics and surface chemistry on the adsorption dynamic was illustrated in Fig. 7. During LCCs compounding, lignocellulosic fibers were thermally degraded into various volatiles, meanwhile generating certain amount of residual ash particles. These ash particles could adhere to the external surface of SARC, and some particles could even block the entrance of macropores once their size matched, thus hindering the diffusion of VOCs from surface into the internal pores. Therefore, the micropores connected with them could not adsorb VOCs and decrease the adsorption rate. As shown in the SEM images, macropore were more easily being blocked, while meso- and micropores were majorly structured with open pores. S-SARC contained a lower proportion of macropores, thus was less vulnerable to ash blockage. More open pores, higher micropore volume as well as larger surface area determined the superior performance of S-SARC towards VOCs adsorption.
As for the chemical properties, it was reported that the surface functionalities of adsorbent can make certain contributions to total surface area, thus enhancing adsorption as well 25. On top of that, the interactions between functional groups and VOCs molecules can generate chemical bonds or electrostatic attractions, allowing the volatile molecules to be efficiently captured. The oxygen-containing groups are the among the most dominant species on the surface of pre-modified carbon, which can be divided into acidic, neutral, and alkaline groups. While for nitrogen-containing groups, the extra p-electrons of pyrollic and quaternary nitrogen at high energy state would promote the oxidation reaction by forming hydrophilic superoxide ions, thus generating abundant active adsorption sites 26. The majority of these two groups prefer to adsorb the polar VOCs. On the other hand, aromatic and other minor groups are mostly neutral or alkaline, and have more affinity towards nonpolar VOCs. After NaOH modification, more alkaline groups on SARC were introduced, thus allowing it to adsob more VOCs from LCCs, which were majorly composed of nonpolar organic compounds.