Some Unmodified Household Adsorbents for the Adsorption of Benzalkonium Chloride—A Kinetic and Thermodynamic Case Study for Commercially Available Paper

The extensive use of biocide surfactant benzalkonium chloride (BAC) during the SARS-CoV-2 pandemic has led to the buildup of this hazardous chemical in waste, surface and groundwater. e The study aims to elucidate whether various low-cost household materials are suitable, in their unmodified and untreated form, to effectively adsorb BAC from its aqueous solutions.. Additionally, if a proper adsorbent is identified, a description of the kinetics and thermodynamics of the process is also targeted. From among the five tested materials, a commercially available white household paper towel was chosen to best satisfy the criteria of low price, large availability, and standardization degree, as well as high adsorption capacity within a fairly short time window needed until equilibrium. Batch experiments were carried out with a commercial mixture of BAC-12 and BAC-14 within a temperature range of 18-45 degrees Celsius, and a 25-100 mg/g BAC/adsorbent initial mass ratio range, respectively. The overall process follows a pseudo-second-order kinetic law, with an apparent activation energy of 73.35 KJ/mole. Both the Langmuir and the Redlich–Peterson isotherms describe the equilibrium data at 298 K well, with a Gibbs free energy of −20.64 KJ/mole. These findings are in agreement with previous reports and suggest a hybrid but chemisorption-dominated process.


Introduction
Benzalkonium chloride (BAC) is a quaternary ammonium salt used mainly as a cationic surfactant with a broad anti-microbial spectrum.It may be observed from the structural formula illustrated in Scheme 1 that it contains a benzyl group, two methyl groups, and an alkyl group.The latter has a varying from 8 to 18, even-numbered carbon atom chain length.The commercially available formulation is usually a 50%wt aqueous BAC solution, containing a mixture of 12 and 14 carbon atom chains in different but usually equal ratios.
Toxics 2023, 11, x FOR PEER REVIEW 2 of 18 reached USD 715.064 million in 2020 and is expected to grow by 8.7% by 2027 [1].Although very effective for the purpose it has been designed for, BAC might accumulate in municipal or other wastewaters.Despite some efforts to inhibit its action or to biodegrade it [2,3], the usual corresponding wastewater treatment plants are not fully designed to deal with this species [4][5][6].Hence, it may reach various natural environments and build up in wastewater influents and/or effluents, surface and groundwater, respectively [6], and as such, raise ecological concerns [7][8][9][10][11][12].
Although there is no official threshold limit for the legally allowed benzalkonium chloride content of surface, waste, and/or groundwater, chemical companies mention 96 hour LC50 values for various fish in the range of ppm concentrations in the aquatic media [13,14]: 0.223 ÷ 0.46 mg/L for Lepomis macrochirus; 0.823 ÷ 1.61 mg/L for Oncorhynchus mykiss; 1.3 mg/L for Poecilia reticulate; and 2.4 mg/L for Oryzias latipes.Material safety sheets state that it is "very toxic to aquatic life" [14,15] and recommend that users avoid BAC discharge in lakes, streams, and sewers [15] and contact the environmental authorities in case of accidental leakages [15].Fuchsman et al. [7] express the opinion that even though BAC-12 and BAC-14 (BAC species in Scheme 1, with n = 12 and 14) are the most prevalent BAC compounds in wastewater influents, they generally "exhibit low bioaccumulation potential, tending to be retained on epithelial surfaces rather than crossing cell membranes".Hence, it is more likely for BAC-16 and BAC-18 (BAC species in Scheme 1, with n = 14 and 16) to be found in fish tissue [7].
Scheme 1.The structural formula of benzalkonium chloride (BAC), where n = 8 ÷ 18 stands for the even number of carbon atoms in the chain of the alkyl group.
It has been reported that BAC has a neurotoxic effect on zebrafish (Danio rerio).Sousa et al. [12] reported increasing overall swimming activity, thigmotaxis behavior, and erratic movements, whereas Sreevidya et al. [8] observed hatching delay or the inhibition of embryos, embryo enhanced mortality, and the prevalence of morphological malformations, Scheme 1.The structural formula of benzalkonium chloride (BAC), where n = 8 ÷ 18 stands for the even number of carbon atoms in the chain of the alkyl group.
One of the drivers of its increasing market importance since the recent SARS-CoV-2caused pandemic is the widespread use of BAC-containing hand sanitizers and surface disinfectants as alternatives to alcohol-based products.As a result, its market value reached USD 715.064 million in 2020 and is expected to grow by 8.7% by 2027 [1].Although very effective for the purpose it has been designed for, BAC might accumulate in municipal or other wastewaters.Despite some efforts to inhibit its action or to biodegrade it [2,3], the usual corresponding wastewater treatment plants are not fully designed to deal with this species [4][5][6].Hence, it may reach various natural environments and build up in wastewater influents and/or effluents, surface and groundwater, respectively [6], and as such, raise ecological concerns [7][8][9][10][11][12].
Although there is no official threshold limit for the legally allowed benzalkonium chloride content of surface, waste, and/or groundwater, chemical companies mention 96 h LC50 values for various fish in the range of ppm concentrations in the aquatic media [13,14]: 0.223 ÷ 0.46 mg/L for Lepomis macrochirus; 0.823 ÷ 1.61 mg/L for Oncorhynchus mykiss; 1.3 mg/L for Poecilia reticulate; and 2.4 mg/L for Oryzias latipes.Material safety sheets state that it is "very toxic to aquatic life" [14,15] and recommend that users avoid BAC discharge in lakes, streams, and sewers [15] and contact the environmental authorities in case of accidental leakages [15].Fuchsman et al. [7] express the opinion that even though BAC-12 and BAC-14 (BAC species in Scheme 1, with n = 12 and 14) are the most prevalent BAC compounds in wastewater influents, they generally "exhibit low bioaccumulation potential, tending to be retained on epithelial surfaces rather than crossing cell membranes".Hence, it is more likely for BAC-16 and BAC-18 (BAC species in Scheme 1, with n = 14 and 16) to be found in fish tissue [7].
It has been reported that BAC has a neurotoxic effect on zebrafish (Danio rerio).Sousa et al. [12] reported increasing overall swimming activity, thigmotaxis behavior, and erratic movements, whereas Sreevidya et al. [8] observed hatching delay or the inhibition of embryos, embryo enhanced mortality, and the prevalence of morphological malformations, with acute toxicity occurring at concentrations of hundreds of µg/L.Ikisa et al. [9] observed a direct positive relationship between O. niloticus juvenile mortality and BAC concentration in freshwater.BAC also induced relevant DNA damage, thus proving its genotoxicity, in both Daphnia magna and Ceriodaphnia dubia [10], with acute effects occurring at tens of µg/L in D. magna and hundreds of µg/L in C. dubia, respectively.On the other hand, comparable ranges of BAC concentrations of 15 to 200 µg/L did not significantly affect the physiological variables of Isochrysis galbana and Chaetoceros gracilis microalgae batch cultures, nor those of a coastal phytoplankton assemblage [11].
Tezel et al. [16] proved that at these low concentrations, counter-ions such as bromide, nitrate, and acetate did not affect the bioavailability, hence the toxicity of BAC, even if present in three-fold higher-order-of-magnitude concentrations.Yet, similar quantities of natural organic matter somewhat reduced the biocide's toxicity.Reported environmentally relevant surface and freshwater BAC concentrations may vary from units of µg/L [10] to hundreds of µg/L [11,12], with levels inducing adverse effects or acute toxicity, depending on the examined species, anywhere within this range.Therefore, environmental risk issues created by the wide and intensive use of benzalkonium chloride may not be ignored when striving for the preservation of a safe and healthy environment [6].
The removal of benzalkonium chloride from aqueous solutions by means of adsorption is also mentioned in the literature.One of the earliest studies reported that its adsorption on activated granular carbon [30] is controlled by intraparticle diffusion and characterized at equilibrium by a Freundlich isotherm.A decade ago, Zanini et al. [31] proposed the use of montmorillonite as an adsorbent and compared the adsorptions of BAC-12 and BAC-14, respectively.Both BAC variants exhibited two distinct stages in the isotherm, the one at lower pollutant concentrations corresponding almost exclusively to a cation exchange process, whereas the one at higher concentrations was considered to be mainly driven by lateral interactions between the surfactant BAC molecules.These studies were continued by adding triclosan, another biocide, to the pollutant's mixture, as well as by putting forward a Langmuir isotherm for the binary adsorption model [32].BAC's competitive adsorption with a herbicide on the same adsorbent was assessed by means of a theoretical cation exchange model [33].A laboratory-synthetized magnetite nanomaterial was also tried recently [34].Again, a Langmuir isotherm described very well the experimental data.
Other types of BAC adsorbents include polymers.Both an acrylic ester, and a styrenic polymer were used to remove it from aqueous streams in a single adsorbent column [35].Similar to previous findings [30], intraparticle mass transfer was considered to significantly affect the process rate.However, concerning the mechanism, the interaction between the adsorbent and adsorbate indicated slightly exothermic physical sorption [35].The adsorption of three BAC variants (BAC-12, BAC-14, and BAC-16) on polyethylene microplastics, as well as their combined toxic effects on Daphnia magna, were assessed by using batch experiments [36].Overall, pseudo-second-order kinetic laws were further proposed.
More detailed adsorption mechanisms involving various alkyl group chain length BACs were carried out just recently [37].The adsorption of BAC-12, BAC-14 and BAC-16, respectively, from aqueous media on powdered activated carbon showed that species with longer alkyl chains were adsorbed more effectively.Therefore, electrostatic interaction is more significant in BAC-12 adsorption, whereas for the other two species, van der Waals interaction plays a more important role.In all cases, pseudo-second-order kinetics and Langmuir isotherms proved to best fit the experimental data.
BAC adsorption was employed for other purposes: to probe the montmorillonite surface reactivity when encapsulated in a composite of alginate hydrogels [38], to decide whether zeolite-surfactant-drug composites are suitable drug deliverers when BAC is co-adsorbed with sulfamethoxazole on purified natural clinoptilolite [39], or to choose the best sterilizing grade filter membranes in the pharmaceutical industry, so that they retain the least BAC possible, and its concentration in the filtrate does not drop below its predetermined specification sterile limit [40,41].
It may be concluded that although there was some preoccupation for benzalkonium chloride removal from various water samples by means of adsorption, research efforts have mostly concentrated on finding the process conditions (temperature, pH, mixing rate, adsorbate/adsorbent mass ratio, ionic strength, process time, etc.) that ensure the most effective BAC elimination.Furthermore, various process mechanisms and isotherm models have been proposed, whereas detailed kinetic aspects have rarely been addressed, almost entirely since the SARS-CoV-2 pandemic.Except for the filtering membranes, which were acquired form commercial sources, all other adsorbents mentioned in the literature were subjected to some physical or chemical treatment prior to the adsorption experiments.
Therefore, this work aims to test the suitability of low-cost, unmodified adsorbents for BAC among various household items.The proposed material ought to be employed effectively in is available form, either from commercial sources or waste.Additionally, supplementary experiments are planned concerning the suitable candidate in order to conclude upon a more detailed kinetic model of the overall process, as well as to put forward an isotherm equation that best explains the adsorption mechanism under the employed conditions.

Materials and Methods
A commercially available 50 wt% aqueous solution of benzalkonium chloride, containing a mixture of BAC-12 and BAC-14 in an unspecified ratio, was provided by SC Grande Gloria SRL, Galat , i, Romania and used as the source of a 5 g/L BAC stock solution.This was employed both in calibration (absorbance vs. BAC concentration) as well as in the adsorption experiments.All were carried out with deionized and double distilled water.
The following adsorbents were chosen from easily accessible sources: Expertto white household paper towel (SC Grande Gloria SRL, Galat , i, Romania), "Anna" Soft and Delicate light blue viscose household cloth (Stella Pack S.A., Lubartów, Poland), white tea filter paper (Demmer Gmbh, Wien, Austria), white cotton string (Românofir S.A., Tălmaciu, Romania), and a mixture of oak and beech sawdust from a private mountain farm in Braşov County, Romania, respectively.The adsorbent's relative humidity was determined via gravimetric measurements (Kern ABT220-5DNM analytical scale) by drying samples of 20 g (Memmert UN30 drying chamber) at 105 • C for 24 h and was assessed as a percentage of the mass before drying.The bulk density of dry sawdust was determined via the triplicate weighting of 200 mL solid.All adsorbents were employed in experiments in the form provided by their respective commercial sources, except the sawdust, which was dried first.No other physical or chemical treatment prior to the BAC adsorption experiments was employed.
Benzalkonium chloride adsorption experiments were carried out in batch arrangements by using 150 mL covered glass recipients in order to limit water evaporation during the course of the experiments.Each water-BAC-adsorbent mixture was prepared in triplicate in a 50 mL aqueous solution of 500 mg/L BAC.The start of the process was operationally defined by the moment at which a weighted amount of solid adsorbent was placed into the BAC solution.The variation of the initial adsorbate/adsorbent mass ratio in five steps, from 25 to 100 mg/g, was achieved by varying the mass of adsorbent placed into the recipients.The temperature was kept constant by pre-warming all of the employed materials and keeping the glasses in a thermostated drying chamber.
No stirring was employed during the experiments other than the rapid mixing of the sample when placing the adsorbent into it.The authors considered that such an approach might resemble the slow passing of wastewater through a filter.
At predetermined time intervals, varying from 30 min to 1-2 days, 3 mL aliquots of the liquid bulk solution were extracted in order to record their absorbance at 262 nm (Jasco V-650 UV-VIS spectrophotometer) by using quartz cuvettes of 1 cm optical path length.When 3 consecutive values did not vary for more than 2%, adsorption equilibrium was considered to be achieved.These aliquots were re-added to their corresponding mixtures in order to keep the liquid bulk's volume constant during the experiment.
The remnant BAC concentration was calculated from the absorbance vs. BAC concentration Lambert-Beer calibration line.This was obtained by using 10 triplicate aqueous biocide solutions in the concentration range of 100 to 2500 mg/L BAC.
Necessary calculus has been carried out by using Microsoft Office Excel and its Data Analysis toolpak, as well as the MatLab Software (Version 2021 A).

Choice of Adsorbent
The five different adsorbent candidates were compared by means of the following criteria in order to choose the most suitable one for further studies: 1.
Low price and easy availability (commercial sources); 2.
Good standardization (already ensured by the producer); 3.
An approximatelly zero value for the 262 nm absorbance of the adsorbent-water mixture, so that it does not interfere with the values provided by BAC; 4.
A high equilibrium adsorption yield coupled with a low total S-L contact time required until equilibrium.
The first three criteria are summarized in Table 1.Prices were estimated by using the Romanian local currency conversion rate into EUR as of 30 October 2023.For the outer surface, only the apparent surface of both sides of a cloth or a sheet of paper is taken into account.It is not to be mistaken for the total surface available for adsorption, which also includes the inner surface of the pores.As such, the values in Table 1 were assessed via calculation from the geometrical size and weight measurements of one cloth, one sheet of household paper, and one tea filter, respectively.The parameter has no practical value except for being a helpful tool in assessing the standardization quality ensured by the producer.The same is true for the relative humidity.The listed absorbance at 262 nm represents the average value registered for five distinct mixtures within the range of 2.5-25 g/L solid adsorbent in water.Figure 1 sheds light on the choice of the use of a 262 nm wavelength to monitor the absorbance of the liquid phase in order to assess its remaining BAC content.Curve (a) corresponds to the UV spectrum of a 500 mg/L aqueous BAC solution.It has three peaks, at 268, 262 and 257 nm, respectively, with the one at 262 nm being the highest.Therefore, it was chosen for the purpose of quantifying BAC concentrations from the absorbance measurements as long as the adsorbent does not add to the A value.This is exactly the case with the white household paper towel.Its UV spectrum is illustrated in curve (b), and shows the negligible value of A = 0.0136 at 262 nm.Curve (c) illustrates the spectrum of the aqueous bulk phase of an initially 500 mg/L BAC solution containing 10 g/L paper at adsorption equilibrium (after 48 h of total S-L contact time).It may be observed that (i) the shape of the spectrum is similar to that of curve (a), demonstrating that the adsorbent does not affect A values, and (ii) the peak at 262 nm can be correlated with the remaining BAC content in the liquid phase.As a result, the difference of absorbance A at 262 nm between curves (a) and (c) is directly correlated with the adsorbed BAC.Methods relying on spectrophotometric determination of BAC use the same wavelength of 262 nm in UV [37] or a close one of 215 nm in the far-UV [41].
Toxics 2023, 11, 950 6 of 17 at adsorption equilibrium (after 48 h of total S-L contact time).It may be observed that (i) the shape of the spectrum is similar to that of curve (a), demonstrating that the adsorbent does not affect A values, and (ii) the peak at 262 nm can be correlated with the remaining BAC content in the liquid phase.As a result, the difference of absorbance A at 262 nm between curves (a) and (c) is directly correlated with the adsorbed BAC.Methods relying on spectrophotometric determination of BAC use the same wavelength of 262 nm in UV [37] or a close one of 215 nm in the far-UV [41].The data in Table 1 indicate sawdust as the first choice in terms of price.However, it presents the disadvantage of needing standardization as well as having a significant absorbance signal at 262 nm.This is due to some extraction-dissolution processes of organic species from sawdust [42,43].Hence, it might interfere with the determination of BAC concentrations via spectrophotometric measurements and could potentially interfere with the BAC adsorption process.The next best option is the household paper towel.It has a reasonable price, the producer ensures a good standardization quality, and no other species interferes with BAC concentration determinations at 262 nm or with its adsorption.
Criterion 4 in the list above was assessed by using the equilibrium 262 nm absorbance values of the aqueous bulk solutions, Aeq, of the water-BAC-adsorbent mixtures, as well as the corresponding A0 values of BAC aqueous solutions of the same concentration as the tests' initial ones.In this study, 10 g/L of each adsorbent was employed for C 0 BAC = 1000 mg/L initial BAC concentrations at 18.0 ± 0.1 °C.The adsorption equilibrium yield ηeq (%) was calculated by means of Equation (1).Equilibrium was considered to be achieved when three consecutive absorbance measurements (at different total S-L contact times) did not vary for more than 2%.Figure 2a,b illustrate the results registered during the assessment of criterion 4 for the choice of the most suitable BAC adsorbent.
The calculus for ηeq was carried out with absorbance values instead of BAC concentrations because not all adsorbents have negligible A262 values in aqueous media, see sawdust, for example.The first four adsorbents in Table 1 do not add to the overall absorbance value at 262 nm, and the A values in Equation ( 1) are only due to the presence of BAC.Hence, the corresponding ηeq values of these adsorbents in Figure 2a  The data in Table 1 indicate sawdust as the first choice in terms of price.However, it presents the disadvantage of needing standardization as well as having a significant absorbance signal at 262 nm.This is due to some extraction-dissolution processes of organic species from sawdust [42,43].Hence, it might interfere with the determination of BAC concentrations via spectrophotometric measurements and could potentially interfere with the BAC adsorption process.The next best option is the household paper towel.It has a reasonable price, the producer ensures a good standardization quality, and no other species interferes with BAC concentration determinations at 262 nm or with its adsorption.
Criterion 4 in the list above was assessed by using the equilibrium 262 nm absorbance values of the aqueous bulk solutions, A eq , of the water-BAC-adsorbent mixtures, as well as the corresponding A 0 values of BAC aqueous solutions of the same concentration as the tests' initial ones.In this study, 10 g/L of each adsorbent was employed for C 0 BAC = 1000 mg/L initial BAC concentrations at 18.0 ± 0.1 • C. The adsorption equilibrium yield η eq (%) was calculated by means of Equation (1).Equilibrium was considered to be achieved when three consecutive absorbance measurements (at different total S-L contact times) did not vary for more than 2%.Figure 2a,b illustrate the results registered during the assessment of criterion 4 for the choice of the most suitable BAC adsorbent.
Toxics 2023, 11, x FOR PEER REVIEW 7 of 18 On the other hand, since the sawdust itself contributes to any A262 value registered for the bulk liquid phase, the Aeq value of the sawdust-BAC mixtures contains contributions from both species.As a result, the numerator (A0 − Aeq) in Equation ( 1) is not quantitatively linked with the remaining liquid phase BAC concentration at equilibrium.Hence, the value of ηeq for sawdust in Figure 2a ought to be regarded with reservations.According to the data presented in Figure 2, the most suitable adsorbent that satisfies criterion 4 above is tea filter paper.In 4 h, it reaches an approximately 80% adsorption yield, and in 24, almost 100%.The next suitable choices would be cotton string, household paper towel and household viscose cloth, in this order.Yet, the cheapest among these is The calculus for η eq was carried out with absorbance values instead of BAC concentrations because not all adsorbents have negligible A 262 values in aqueous media, see sawdust, for example.The first four adsorbents in Table 1 do not add to the overall absorbance value at 262 nm, and the A values in Equation ( 1) are only due to the presence of BAC.Hence, the corresponding η eq values of these adsorbents in Figure 2a are pertinent.
On the other hand, since the sawdust itself contributes to any A 262 value registered for the bulk liquid phase, the A eq value of the sawdust-BAC mixtures contains contributions from both species.As a result, the numerator (A 0 − A eq ) in Equation ( 1) is not quantitatively linked with the remaining liquid phase BAC concentration at equilibrium.Hence, the value of η eq for sawdust in Figure 2a ought to be regarded with reservations.
According to the data presented in Figure 2, the most suitable adsorbent that satisfies criterion 4 above is tea filter paper.In 4 h, it reaches an approximately 80% adsorption yield, and in 24, almost 100%.The next suitable choices would be cotton string, household paper towel and household viscose cloth, in this order.Yet, the cheapest among these is the towel.According to the data in Table 1, it also would be a good choice.
Hence, the white household paper towel provided by Expertto (SC Grande Gloria SRL, Galat , i, Romania) was chosen for further BAC batch adsorption experiments for the purpose of establishing the kinetic and thermodynamic aspects of this process.

Experimental Data and Their Processing Method
The BAC adsorption experiments from its synthetic aqueous solutions on solid household paper towels were carried out, as described in Section 2. The one variable at a time approach was used.Firstly, the temperature varied stepwise at constant initial BAC and adsorbent contents in the mixtures (constant adsorbate/adsorbent initial mass ratio).Then, the temperature was kept constant at 25.0 ± 0.1 • C, and the adsorbate/adsorbent initial mass ratio varied stepwise.These data ought to elucidate the process thermodynamics and serve for the calculus of the equilibrium constant at 25 • C, the standard temperature.
Since, in all cases, the absorbance values were recorded at predetermined values of time, more often after the initiation of the process and less often towards the equilibrium, all data serve for the elucidation of process kinetics, as well as for the calculus of overall rate coefficients.
A calibration line absorbance vs. C BAC was obtained at 262 nm in order to link the absorbance of the liquid phase bulk with the remaining, un-adsorbed BAC concentration.Its mathematical expression is provided in Equation ( 2), where C BAC is expressed in mg/L.
The R 2 adjusted = 0.9974 value demonstrates the good quality of the calibration.The slope permits the calculus of the molar absorption coefficient for known BAC variants.Because BAC solutions have been prepared from a commercial source with an unspecified BAC-12 and BAC-14 ratio, concentrations are expressed in "overall" mg/L BAC instead of mole/L.Therefore, the slope of Equation ( 2) is expressed in (mg/L) −1 •cm −1 .
For correct C BAC calculus, the A 262 values used in Equation ( 2) were corrected by subtracting the adsorbent's own absorbance (of an average value of 0.015 as listed in Table 1) from the one of the test's mixture.
Further, C BAC values served to calculate the adsorption capacity q (mg/g), according to Equation (3).C 0 BAC (mg/L) stands for the initial concentration of BAC, whereas C BAC (mg/L) stands for the one at a certain moment t after process initiation, respectively.C adsorbent (g/L) is the adsorbent content of the aqueous BAC-solid adsorbent mixture.Hence, q stands for the mass (mg) of BAC adsorbed by 1 g of adsorbent after a certain S-L contact time t.
Toxics 2023, 11, 950 8 of 17 The equilibrium value of q is represented by q eq and was determined experimentally in each case.Figure 3a presents a set of typical results of calculus with Equation (3) in the form of the adsorption capacity q vs. total S-L contact time profile at 25.0 ± 0.1 • C and 50 mg/g initial BAC/household paper towel mass ratio.The star symbols stand for all of the data derived from the triplicate experiments.Negligible differences among the triplicates may be observed (RSD < 5%).The data show that approximately 90% of the process occurs within the first 5 h, with dramatic q value increases.Yet, equilibrium is attained far more slowly, taking about three or more days.These findings are in agreement with the data reported in the literature [31,36,37].
Cadsorbent (g/L) is the adsorbent content of the aqueous BAC-solid adsorbent mixture.Hence, q stands for the mass (mg) of BAC adsorbed by 1 g of adsorbent after a certain S-L contact time t.
The equilibrium value of q is represented by qeq and was determined experimentally in each case.Figure 3a presents a set of typical results of calculus with Equation (3) in the form of the adsorption capacity q vs. total S-L contact time profile at 25.0 ± 0.1 °C and 50 mg/g initial BAC/household paper towel mass ratio.The star symbols stand for all of the data derived from the triplicate experiments.Negligible differences among the triplicates may be observed (RSD < 5%).The data show that approximately 90% of the process occurs within the first 5 h, with dramatic q value increases.Yet, equilibrium is attained far more slowly, taking about three or more days.These findings are in agreement with the data reported in the literature [31,36,37].The first line of experiments was carried out at 18, 25, 35 and 45 • C, respectively, by keeping the starting BAC/paper mass ratio at constant 50 mg/g.The second line of trials kept the temperature constant at 25 • C, and varied the initial adsorbate/adsorbent mass ratio by giving it values of 25, 37.5, 50, 75 and 100 mg/g, respectively.In all cases, similar curves to the one depicted in Figure 3a were obtained, with significantly shorter time frames necessary until equilibrium at higher temperatures.
Single-factor ANOVA statistical testing was employed in order to assess whether the q = f(time) profiles like the one in Figure 3a were statistically significantly affected by the variation of either the temperature or the initial adsorbate/adsorbent mass ratio.Calculus was performed by using the Data Analysis toolpak of Microsoft Office Excel.
Statistical testing based on analysis of variance (ANOVA) is a technique that allows the isolation and estimation of variations between various data sets.It further classifies the differences as statistically "significant" or "insignificant", depending on the value of the Fisher criterion (the F value).If its value exceeds the F significance threshold, then the differences are statistically significant [44] (pp.59-65).Otherwise, they are statistically insignificant.The threshold value depends on the number of degrees of freedom, which, in turn, is related to the number of experimental data subjected to testing.
For the effect of temperature, the value of the Fisher criterion was found to be equal to 17.62, higher than the threshold F significance = 2.68.Similar results were obtained for the effect of the initial BAC/paper mass ratio, with F = 1874.4being much higher than F significance = 2.41.The results show that both parameters affect statistically significantly the course and outcome of the adsorption process, yet with a far more dramatic influence by initial conditions than the temperature (at least within the range of employed process conditions).

Process Kinetics
The q = f(time) profiles obtained via the above-described means have been tested against the pseudo-first, pseudo-second, and intra-particle diffusion models described by Yagub et al. [21], as well as against a more detailed kinetic model consisting of two parallel competing first-order individual stages [45].The only fit was the one describing a pseudosecond-order overall adsorption process.According to it, the value of q follows a hyperbolic profile until it reaches the equilibrium value q eq , see Figure 3a, with a pseudo-second-order kinetic rate coefficient k expressed in (g/mg•h), see Equation (4).
Its linear form (5) permits the calculus of both the rate coefficient k and the equilibrium adsorption capacity q eq , from the intercept and slope of the linear plot of the (t/q) values against t. Figure 3b presents this for the data shown in Figure 3a.Good linearization is observed, proving the suitability of the chosen kinetic model.
The plot in Figure 3b enabled the calculation of q eq = 23.53mg/g, which is practically equal to the experimental 23.56 mg/g value (see also Table 3).Similar matches were found for the other employed experimental conditions.This fact is a further argument in favor of the kinetic model described in Equation ( 4).Such findings are in agreement with previously reported ones [36,37].Equation ( 5) was used to determine adsorption rate coefficients for three more adsorbents of Table 1: household viscose cloth, tea filter paper, and cotton string, respectively, at 18.0 ± 0.1 • C and 50 mg/g adsorbate/adsorbent initial mass ratio.For these species, the calculus algorithm, described in Section 3.2, was also suitable because the adsorbent did not contribute to the overall absorbance value of the bulk liquid phase at 262 nm.The results are summarized in Figure 4 and are in good correlation with those presented in Figure 2b, supporting the decision favoring the use of this paper towel as a suitable BAC absorbent.
Table 2 contains the calculated pseudo-second-order rate coefficients k as a function of the employed temperature, T, as well as the statistical R 2 adjusted values of the linear regressions from which they were calculated, such as the one in Figure 3b.The last column proves the good quality of the (i) experimental data and (ii) linearizations, hence model fit.
species, the calculus algorithm, described in Section 3.2., was also suitable because the adsorbent did not contribute to the overall absorbance value of the bulk liquid phase at 262 nm.The results are summarized in Figure 4 and are in good correlation with those presented in Figure 2b, supporting the decision favoring the use of this paper towel as a suitable BAC absorbent.Table 2 contains the calculated pseudo-second-order rate coefficients k as a function of the employed temperature, T, as well as the statistical R 2 adjusted values of the linear regressions from which they were calculated, such as the one in Figure 3b.The last column proves the good quality of the (i) experimental data and (ii) linearizations, hence model fit. is a strong statistical tool when comparing models with different numbers of independent variables m and data points n.Therefore, it is a better indicator of model quality than the simple R 2 , because it shows the variability in the target variable due only to independent variables that actually affect it [44] (pp.92-108), [46,47].The relationship between these two statistical parameters is described in Equation ( 6).It may be observed that, for the same m = 1 (the case here), if n differs, so will R 2 , and hence R 2 adjusted.The value of the latter being close to 1 demonstrates that the data fit very well with the proposed equation, in this case, with the second-order kinetic law.
The rate coefficient values in Table 2 have been used to calculate the overall activation energy Ea and the pre-exponential factor Z by means of the linearized Arrhenius law.It fits fairly well (R 2 adjusted = 0.7855) with the four experimental ln(k) vs. 1/T values.Hence, Ea = 73.35± 21.26 KJ/mole and ln(Z) = 27.01 ± 8.43 were obtained, pointing to the chemisorption of BAC as the main adsorption mechanism.Previously, similar [31] and opposite (physical sorption) findings [35] have been previously reported, mainly  adjusted is a strong statistical tool when comparing models with different numbers of independent variables m and data points n.Therefore, it is a better indicator of model quality than the simple R 2 , because it shows the variability in the target variable due only to independent variables that actually affect it [44] (pp.92-108), [46,47].The relationship between these two statistical parameters is described in Equation ( 6).It may be observed that, for the same m = 1 (the case here), if n differs, so will R 2 , and hence R 2 adjusted .The value of the latter being close to 1 demonstrates that the data fit very well with the proposed equation, in this case, with the second-order kinetic law.
The rate coefficient values in Table 2 have been used to calculate the overall activation energy E a and the pre-exponential factor Z by means of the linearized Arrhenius law.It fits fairly well (R 2 adjusted = 0.7855) with the four experimental ln(k) vs. 1/T values.Hence, E a = 73.35± 21.26 KJ/mole and ln(Z) = 27.01 ± 8.43 were obtained, pointing to the chemisorption of BAC as the main adsorption mechanism.Previously, similar [31] and opposite (physical sorption) findings [35] have been previously reported, mainly depending on the nature of the adsorbent but also on the initial adsorbate/adsorbent mass ratio.
The experiments carried out at a constant 25 • C and variable initial BAC/paper mass ratio generated the rate coefficient values presented in Figure 5.It shows that the overall adsorption process is favored by both low and high ratios.At low ones, the adsorption might be favored by the large number of available empty sorption sites, thus occurring faster.On the other hand, at high mass ratios, the concentration gradient of BAC between the liquid bulk phase and the adsorbent's surface increases, favoring diffusion, hence raising the overall rate constant.concentration in the liquid bulk phase, Ceq (mg/L).Both experimental and calculated qeq values are provided.The latter are obtained from plots described through the use of Equation ( 5) and exemplified by Figure 3b.There is a very good match between the qeq value pairs.Hence, (i) the suitability of the second-order kinetic model is reinforced since the slope of its linearized form (5) serves for the determination of qeq calculated and (ii) the good quality of the experimental data is proven (see also the corresponding excellent R 2 adjusted values).34,37] and Freundlich [30] isotherms are mentioned in the literature, the fitting quality of eight different adsorption isotherms [47][48][49][50] over these data was evaluated in order to decide upon one that might elucidate the

Adsorption Isotherms and Mechanism
The experiments carried out until equilibrium under the conditions described in Figure 5 led to the data presented in Table 3.It contains the equilibrium adsorption capacities q eq (mg/g), which are quantitatively linked to the BAC retained on the adsorbent's surface as a function of their corresponding equilibrium remnant concentration in the liquid bulk phase, C eq (mg/L).Both experimental and calculated q eq values are provided.The latter are obtained from plots described through the use of Equation ( 5) and exemplified by Figure 3b.There is a very good match between the q eq value pairs.Hence, (i) the suitability of the second-order kinetic model is reinforced since the slope of its linearizted form (5) serves for the determination of q eq calculated and (ii) the good quality of the experimental data is proven (see also the corresponding excellent R 2 adjusted values).
Table 3. Equilibrium data: BAC adsorption capacities q eq (mg/g) vs. their corresponding equilibrium remnant concentration in the liquid bulk phase, C eq (mg/L).Experimental and calculated data, at 25 ± 0.1 • C.
The Langmuir model describes the ideal situation for which all adsorption sites are homogeneous in structure and, therefore, energetically equivalent.The adsorption is assumed to take place in a monolayer, one site at a time, only through chemisorption, without any consequences on the neighboring positions of an occupied site.Its equation is presented in Table 4. Its linearized form in Equation ( 7) permits the easy calculus of both the maximum adsorption capacity of a single layer at a given temperature, q max (mg/g), and the adsorption equilibrium constant at that temperature, K L (L/mg).
C eq q eq = C eq q max + 1 K L q max (7) The data presented in Table 3 fitted well Equation ( 7), with R 2 adjsuted = 0.9562.Hence, at the standard temperature of 298 K, the following were obtained: q max = 35.09mg/g and K L = 11.76L/g, respectively.
However, despite being easy, the use of simple linear regressions to evaluate the fitting quality and to calculate the parameters of a model can generate significant errors.Moreover, it is limited to models that allow linearization.Therefore, the parameters of the eight isotherms were calculated by means of an optimization algorithm in MatLab Software.This works with the normalized mean square error (NMSE) as an objective function and strives to find its minimum.Its expression is provided in Equation (8).q exp(i) and q calc(i) stand for the i th experimental and calculated q values, respectively.The denominator represents the average value of q exp , and n is the number of experiments (here n = 5).The decision variables of the objective function are the sought-after model parameters.
After finding the parameters of each adsorption isotherm, three more error functions were calculated as simultaneous tools to be used when deciding upon the best model: (1) RMSD function (root mean squared deviations), see Equation ( 9), (2) HYBRID function (hybrid fractional error deviation), see Equation (10) in which p is the number of parameters of the model [42,50], and (3) R 2 adjusted , see Equation (6).
Table 4 lists the calculus results for the eight isotherms taken into consideration: the value of their parameters as well as of the three above-mentioned error functions.The footnotes indicate the units of the parameters for each isotherm [49,50].A good quality fit simultaneously satisfies the following requests: R 2 adjusted as close to 1 as possible, whereas RMSD and HYBRID ought to be as close to zero as possible.
The analysis of the values in Table 4 suggests that the models fitting the best data in Table 3 are the Redlich-Peterson (RP) isotherm (lowest RMSD, relatively low HYBRID, very good R 2 adjusted of 0.988), seconded by the Langmuir (L) isotherm (best R 2 adjusted , yet weakest HYBRID), respectively.A comparison between the experimental values of q eq and those calculated by these two models is presented in Figure 6.Both seem to do a good job (see also results of Table 4); however, the L-isotherm predicts values that are scattered among the experimental ones, whereas the RP-equation predicts lower q eq values, regardless of the C eq level.Hence, despite the better statistical error functions of the RP model, the Langmuir one still holds best.Similar conclusions have been previously reported [32,34,37].The RP isotherm combines the characteristics of both the Langmuir and Freundlich (F) ones [48,49].Unlike the L-model, it does not assume ideal monolayer adsorption.The equation in Table 4 shows that when the parameter g = 1, it will reduce to the Langmuir relationship.On the other hand, at high liquid-phase adsorbate concentrations or low g values, the model reduces to the F-isotherm.In this case, the value of g = 0.780 shows a stronger similarity with the L-model.This fact also supports the conclusions reached in The RP isotherm combines the characteristics of both the Langmuir and Freundlich (F) ones [48,49].Unlike the L-model, it does not assume ideal monolayer adsorption.The equation in Table 4 shows that when the parameter g = 1, it will reduce to the Langmuir relationship.On the other hand, at high liquid-phase adsorbate concentrations or low g values, the model reduces to the F-isotherm.In this case, the value of g = 0.780 shows a stronger similarity with the L-model.This fact also supports the conclusions reached in Figure 6.
Kim et al. [37] mention somewhat different sorption mechanisms for various BAC molecules: electrostatic interaction seems to be more significant in the case of BAC-12, whereas van der Waals bonding is more relevant for BAC molecules with longer alkyl chains.Since this study was carried out with a mixture of BAC-12 and BAC-14 in an unknown ratio, the equilibrium data describe a rather overall process, which seems to mainly follow chemisorption, with some physisorption (condensation) going on at high BAC concentrations: after completion of the monolayer, further BAC molecules can still attach themselves to the adsorbent's surface via Van der Waals bonding to the previous layers.
These conclusions are supported by the findings of Zanini et al. [31], who studied the synergistic effect of BAC-12 and BAC-14 when simultaneously adsorbed from mixtures.An isotherm with two well-differentiated steps was proposed, each corresponding to a different adsorption mechanism.These impregnate a letter S-shaped q eq vs. C eq plot, each plateau (leveling off) corresponding to a step: the first being explained by chemical and the second by physical processes.The experimental data in Figure 6 exhibit such a shape somewhat, with a leveling off in the 100-320 mg/L C eq domain.However, there are too few points to correctly apply a double L-model (a sum of two terms, like the Langmuir equation presented in Table 4) [31], with four parameters to be determined from only five q eq -C eq pairs.It is possible to estimate the Gibbs free energy of adsorption ∆G 0 by using the Van't Hoff relationship (11) and the value of the equilibrium constant; however, the latter ought to be dimensionless [51].
∆G 0 = −RTln K 0 e (11) Yet, the K L = 11.76L/g provided by the Langmuir model is not dimensionless.Equation ( 12) [51] describes a way of computing the dimensionless equilibrium constant K 0 e in such cases: The [adsorbate]/[adsorbate] 0 = 353.5 g/mole is the dimensionless standardized average molar mass of the employed BAC-12 and BAC-14 mixture (the ratio is unspecified by the provider; however, for the sake of this calculus, a 1:1 value was considered); oγ is the activity coefficient, which was here considered to be equal to 1.The obtained K 0 e value enabled the calculus of ∆G 0 with the help of (11) (R = 8.314 J/mole•K, T = 298 K).The resulting value of the ∆G 0 = −20.64KJ/mole indicates a spontaneous process at 298 K.The magnitude of this adsorption free energy is comparable to those reported for other cationic adsorbates [52], yet it is approximately five-fold higher than values reported for BAC physisorption on powdered activated carbon [37].Hence, again, a somewhat hybrid chemi-and physisorption mechanism may be assumed for the adsorption of BAC on the studied household paper.
This study proposes a low-cost, unmodified adsorbent for benzalkonium chloride from its aqueous solutions prepared by using a commercial mixture of BAC-12 and BAC-14.Both are readily accessible to the general public, both in terms of availability as well as price.Even so, the achieved adsorption capacities at 25 • C are comparable to the ones mentioned previously for another low-cost adsorbent: 11 ÷ 22 mg/g on polyethylene microplastics [40].On the other hand, up to one order of magnitude higher values of (100 to 500 mg/g) may be targeted when using activated carbon [30,37] or montmorillonite [31,33,36,38].

Figure 1 .
Figure 1.UV absorbance spectrum of the liquid phase for (a) a 500 mg/L aqueous BAC solution, (b) a water-paper mixture containing 10 g/L paper, and (c) a water-BAC-paper mixture containing initially 500 mg/L BAC and 10 g/L paper, at adsorption equilibrium and at 18.0 ± 0.1 °C, respectively.

Figure 1 .
Figure 1.UV absorbance spectrum of the liquid phase for (a) a 500 mg/L aqueous BAC solution, (b) a water-paper mixture containing 10 g/L paper, and (c) a water-BAC-paper mixture containing initially 500 mg/L BAC and 10 g/L paper, at adsorption equilibrium and at 18.0 ± 0.1 • C, respectively.

Figure 2 .
Figure 2. (a) Values of the equilibrium adsorption yield; (b) total S-L contact time until equilibrium.Data describe the 5 different employed adsorbents, with corresponding error bars for triplicate determinations, at 50 mg/g BAC/adsorbent initial mass ratio and 18.0 ± 0.1 °C, respectively.

Figure 2 .
Figure 2. (a) Values of the equilibrium adsorption yield; (b) total S-L contact time until equilibrium.Data describe the 5 different employed adsorbents, with corresponding error bars for triplicate determinations, at 50 mg/g BAC/adsorbent initial mass ratio and 18.0 ± 0.1 • C, respectively.

Figure 4 .
Figure 4. Pseudo-second-order rate coefficients k for various adsorbents at 18.0 ± 0.1 °C and 50 mg/g initial BAC/adsorbent mass ratio.

Figure 4 .
Figure 4. Pseudo-second-order rate coefficients k for various adsorbents at 18.0 ± 0.1 • C and 50 mg/g initial BAC/adsorbent mass ratio.

Figure 5 .
Figure 5.The dependence of overall pseudo-second-order rate coefficients k on the BAC/paper initial mass ratio, at 25.0 ± 0.1 °C.

Figure 5 .
Figure 5.The dependence of overall pseudo-second-order rate coefficients k on the BAC/paper initial mass ratio, at 25.0 ± 0.1 • C.

Table 2 .
Pseudo-second-order rate coefficients k as a function of temperature, at 50 mg/g initial BAC/household paper towel mass ratio.

Table 2 .
Pseudo-second-order rate coefficients k as a function of temperature, at 50 mg/g initial BAC/household paper towel mass ratio.

Table 4 .
Equations, parameters and error function values for the tested adsorption isotherm models.