Non-linear equation Parameters Freundlich ${q}_{e}={K}_{F}{C}_{e}^{1/n}$ KF (mg∙g−1)(L∙mg−1)1/n: Freundlich adsorption constant n (dimensionless): empirical parameter representing the energetic heterogeneity of surface qe (mg∙g−1): equilibrium adsorbed quantity Langmuir ${q}_{e}=\frac{{q}_{m}{K}_{L}{C}_{e}}{1+{K}_{L}{C}_{e}}$ qm (mg∙g−1): monolayer adsorption capacity Ce (mg∙L−1): equilibrium concentration KL (L∙mg−1): equilibium adsorption constant Dubinin-Radushkevich ${q}_{e}={q}_{s}\mathrm{exp}\left(-{k}_{ad}{\left[RT\mathrm{ln}\left(1+\frac{1}{{C}_{e}}\right)\right]}^{2}\right)$ qs: maximum adsorbed amount R (8.314 J∙mol−1K−1): gas constant T (K): temperature Pseudo-first-order ${q}_{t}={q}_{e}\left[1-\mathrm{exp}\left(-{k}_{1}t\right)\right]$ qt (mg∙g−1): adsorbed quantity at time t t (min): contact time k1 (min−1): pseudo-fist order rate constant Pseudo-second-order ${q}_{t}=\frac{{q}_{e}^{2}{k}_{2}t}{1+{q}_{e}{k}_{2}t}$ k2 (g.mg−1∙min−1): pseudo-fist order rate constant Intraparticular diffusion ${q}_{t}={k}_{ip}{t}^{2}+{C}_{i}$ Kip (mg∙g−1∙min1/2): intraparticle diffusion rate constant Ci: constant value depicting the boundary layer effect