Sulfate Kinetics and Adsorption Studies on a Zeolite/Polyammonium Cation Composite for Environmental Remediation

: Sulﬁde mineral mining produces highly sulfate-contaminated wastewater which needs to be treated before disposal. A composite material was made from natural zeolite (NZ) and Superﬂoc ® SC-581, a polyammonium cationic polymer. The resulting modiﬁed zeolite (MZ) demonstrated improved capacity for sulfate abatement from wastewater compared to NZ. Above pH 4.0, MZ retained positive surface charge while NZ remained negative. The effect of the ionic strength on the adsorption process was evaluated. Sulfate adsorption capacity was assessed and revealed MZ to be superior to NZ in all cases. Adsorption kinetics reached equilibrium after 10–12 h, with MZ adsorption being twice that of NZ; data ﬁtted a pseudo-second order kinetic model. Adsorption isotherms reﬂected the high capacity of MZ for sulfate adsorption with maximum of 3.1 mg g − 1 , while NZ only achieved 1.5 mg g − 1 . The process corresponds to heterogeneous partially reversible adsorption of ionic species over the solid adsorbent. Langmuir–Freundlich parameters revealed that adsorption over MZ corresponds to an interaction eight times stronger than that on NZ. The sulfate adsorption pattern changes with ionic strength. Taken together, the composite formed between natural zeolite and polyammonium represents an adsorbent that maintains the adsorption capacity of zeolite and proves suitable for anionic species removal. Further prospect considers the testing of the composite with other anionic pollutants (arsenate, phosphate, perchlorate, etc.)


Introduction
Mining activities worldwide represent an important economic activity, with reported annual incomes near USD 600 billion during 2017 from the top 40 mining companies [1]. Mining has inherent environmental consequences, as intense mineral processing produces large volumes of industrial waste [2,3]. Mineral extraction from deposits, metal purification with leaching, and production of smelter wastes result in wastewater of varied composition depending on the source and local geological composition. Several reports indicate that typical wastewater from Chilean mines [4] contains high levels of sulfate (ca. 14,400 ppm). The release of sulfate-polluted wastewaters into natural water bodies presents significant environmental consequences [5] such as sulfide liberation through microorganism anaerobic respiration from sulfates [6] and eutrophication [7]. Therefore, treatment of wastewater generated from mining activities is critical for sustainable development [8][9][10].
Sulfate abatement from wastewater remains a technological challenge. Commonly employed methods include chemical precipitation [11], biological treatment [12] and adsorption techniques [13,14]. Chemical precipitation, biological or adsorption treatments are associated with high installation and operational costs that may contain harmful substances with greater health and environmental risks than sulfate [15]. A suitable sulfate remediation method would ideally involve naturally abundant, low-cost materials able to reduce water sulfate, but with environmentally friendly products.
Notably, Chile has enormous deposits of zeolite, a renowned adsorbent material used in a range of applications from adsorption to catalysis [16]. Natural zeolite is an aluminosilicate mineral that typically presents permanent negative surface charge due to isomorphic substitution of its Si 4+ ions with Al 3+ [17]. Due to its high cation exchange ability and molecular sieve properties, natural zeolite (NZ) has been widely used as a metal adsorbent. [18,19]. The adsorption characteristic of zeolites is influenced by their high specific surface area, with active sites that have the capability to adsorb a large variety of compounds. Some authors have applied zeolite to remove ammonium from wastewater, recovering nitrogen nutrients and preventing eutrophication of natural water bodies [20][21][22]. These properties can be modulated by chemical treatments to promote anion adsorption efficiency of natural clays [23][24][25][26][27]. For instance, cationic surfactants [13,[28][29][30] polycationic chitosan [31], and cationic polyacrylamide [32] can be used for this purpose. Due to typical cationic surfactant behavior, these compounds are adsorbed over zeolite through the positively charged head, leaving the non-polar tails in solution. Then, a second layer of cationic surfactant deposits through their lipophilic tails, leaving the positively charged head in contact with the bulk solution and free to adsorb anions, as explained by Kamble et al. [33].
In this work, we describe for the first time the preparation of a composite made from Chilean natural zeolite samples and Superfloc ® SC-581, a polyamine polycation (or cationic polymer) of commercial origin. Superfloc ® SC-581 is a co-polymer between neutral and cationic polyacrylamide commonly used as an organic flocculant for wastewater treatment in mining activities [34,35]. The modified zeolite (MZ) composite is assayed as adsorbent for sulfates in a synthetic sulfate-polluted water solution. Surface charge analysis and sulfate adsorption experiments were performed to optimize the modification of NZ and assess the sulfate affinity of the resulting composite. Since mining wastewater contains a high ionic concentration, the ionic strength effect [36,37] on adsorption is also studied. Finally, adsorption kinetics and isotherms were performed and modeled to evaluate the adsorption capabilities of the composite material.

Materials
Natural zeolite (NZ) samples were obtained from a source located in the centralsouth zone of Chile (geographical coordinates 36 • 16'S 71 • 40'W). Before adsorption studies, natural zeolite was characterized to check its suitability as adsorbent material. Properties of natural zeolite (NZ) are presented in Table 1. Several characteristics of NZ reflect its suitability to environmental remediation applications. First, the pH value of NZ is close to neutral (6.6), suggesting that the application of this mineral in remediation technologies would not result in major pH shifts to water bodies. Moreover, the Si/Al ratio as determined though X-Ray diffractometry (XRD) analysis corresponds to a medium level (5.5) according to the Van Bekkum et al. [38] proposed scale, indicating considerable isomorphic substitution of Si 4+ atoms by Al 3+ to produce a mineral with a permanent negative surface charge. According to the isoelectric point (IEP) value, this negative surface charge is developed from a relatively low pH (3.1). This is important for synthesizing the composite: SC-581 is a permanent positive charge water-stable polymer (or polycation); thus, formation of the composite between the negatively charged NZ and SC-581 is feasible simply through electrostatic attraction. XRD analysis reveals that NZ corresponds to a rehydrated mordenite-type zeolite that possesses exchangeable calcium in its structure as represented by its general formula, Ca 3.4 Al 7.4 Si 40.6 O 96 (H 2 O) 31 [39]. Finally, approximately 75% of NZ is 20 µm or less in size. This property is promising, since particles in this size range present the highest adsorption capacities due to their elevated contact surface [39][40][41]. Collectively, these characteristics make NZ a suitable starting material for further development of composites to be used for adsorption of anionic species such as sulfates. Then, modified zeolite was prepared by combination with SC-581.

Methods
The pH was measured in the supernatant of the solid: double distilled water in a mass ratio of 1.0:2.5 after 2 h of continuous orbital stirring and 3 h of rest.
The total specific surface area (TSSA) was determined gravimetrically, applying the ethylene glycol monoethyl ether (EGME) technique, proposed by Carter et al. [43] to determine the TSSA of the soils.
External surface area was determined using a 0.3 g sample. The sample was firstly degassed and adsorption-desorption isotherms of N 2 at −195 • C (77 K) were acquired, changing the relative pressure (P/Po) of the gas and recording the volume adsorbed on the solid's surface. The specific surface area was calculated from the amount of N 2 adsorbed employing the Brunauer-Emmett-Teller or BET equation [44].
The silicon/aluminum (Si/Al) ratio and the chemical formula of NZ was determined by X-ray diffraction experiments run on a Siemens D5000 X-Ray diffractometer (Germany) with Bragg-Brentano geometry and an X-ray tube with copper anode (wavelength 1.54 Å). Data analysis was carried out using the Direct Plus EVA 15 (2009) database PDF-2 (Powder Diffraction File, 2003). The particle size distribution was measured for each adsorbent through methodology based on Stoke's law relating sedimentation speed with particle diameter [45]. Exchangeable cations were determined for NZ samples. After hominization with KCl 0.05 mol L −1 , samples were centrifuged and the quantity of Ca and Mg at the supernatant were determined through Inductively Coupled Plasma Atomic Emission Spectroscopy analysis.
Isoelectric point (IEP) values were determined by micro-electrophoresis directly on NZ and MZ particles with reported particle size distribution. Electrophoretic mobilities were measured with a zeta meter ZM-77 apparatus (USA). Dilute dispersions of adsorbents (0.05 g L −1 ) were prepared in 10 −3 mol L −1 of KCl and pH was adjusted with 10 −2 mol L −1 of HCl or NaOH. The mobilities were averaged and the zeta potential was calculated using the Helmholtz-Smoluchowski equation [46]. Experimental points obtained showed a 10% coefficient of variation in zeta potential.

Modification Process
Zeolite was firstly washed by adding 1 g in 50 mL centrifuge tubes with 25 mL of double distilled water. The mix was orbitally agitated at 300 rpm for 30 min and centrifuged at 6000 rpm for 5 min and the supernatant was eliminated. This washing process was repeted three times. The washed zeolite was then homoionized. Homoionization of the zeolite was identical to the wash steps but used 0.1 mol L −1 of KCl solution. After the third elimination of supernatan, the homoionized zeolite was freeze-dried at low presure.
Modified zeolite was prepared by placing 1 g of homoionized NZ in 50 mL centrifuge tubes with 25 mL of SC-581 at a specific concentration (0.1 to 12.0 mg L −1 ) followed by orbital agitation for 3 h. The mixture was then separated through centrifugation (6000 rpm per 10 min). The supernatant was discarded, and the solids were stored for further analyses, including measurement of pH, surface area, particle size distribution, zeta-potential, and isoelectric point.

Ionic Strength Effect
Adsorption isotherms of sulfate onto KCl-homoionized zeolite were obtained through batch experiments at different ionic strengths using 0.001, 0.050 and 0.100 mol L −1 of KCl. The concentration of sulfate at equilibrium ranged from 100 to 500 mg L −1 . Sulfate in solution at equilibrium was determined through turbidimetry.

Adsorption Kinetics
In 50 mL polypropylene centrifuge tubes, 20 mL of 1000 ppm SO 4 2− solution in 0.05 mol L −1 of KCl was mixed with a 5% w/v solid ratio of adsorbent (NZ or MZ). The tubes were stirred (300 rpm) in a reciprocating shaker using time intervals of 1, 3, 6, 9, 12, 16, 20 and 24 h. At the end of each time interval, the tubes were centrifuged at 6000 rpm for 10 min. The supernatant was filtered and used for the quantification of sulfate ions using the turbidimetric method. The experimental results were used to calculate the amount of adsorbed sulfate as a function of time, C t (mg g −1 ), determined by mass balance, according to the equation: where C S0 and C St are the sulfate initial and final concentrations (mg L −1 ), respectively, in the supernatant; V (L) is the total volume of solution in the tubes; and m is the mass (g) of adsorbent employed. Adsorption kinetic data were fitted employing pseudo-first and -second order models. The pseudo-first order kinetic model assumes that a solute ion is adsorbed on a surface site of the adsorbent [47,48]. The integrated form of the pseudo-first order kinetic rate equation is: where C t is the adsorbed quantity at any time (t), C m is the maximum experimental adsorbed amount (obtained from the kinetic curve), k 1 is a combination of the adsorption and the desorption rate constants, and C m-cal corresponds to the maximum adsorbed amount as determined by the application of the model. C m-cal and C m will agree if the model succeeds in giving a perfect description of the experimental data. The pseudo-second order kinetic model postulates that solute adsorption occurs on two available surface sites of the adsorbent. The integrated form of the pseudo-second order kinetic rate equation is [47][48][49]: The observed rate constant of the pseudo-second order equation (k 2 ) is a complex function of the initial solute concentration. This linear plot may be described as a chemisorption process.

Adsorption Isotherms
Experimental points were obtained in a similar manner to those for adsorption kinetics but with the sulfate concentration ranging from to 0 to 700 mg L −1 of SO 4 2− in 0.05 mol L −1 of KCl. Quantification and mass balance were also made the same way. Adsorption isotherms data were fitted using Freundlich [50] Langmuir [47], and the Langmuir-Freundlich [51,52] adsorption models. The Freundlich model does not have a mechanistic interpretation; it only represents an empirical approach to predict the species distribution between a solid/solution phase, and described by: (4) where K F and n fads are adjustable parameters with 0 < 1/n fads < 1 (usually), C s is the concentration of adsorbed species (mg g −1 ), and C Se is the concentration in the supernatant solution (mg L −1 ). K F corresponds to the empirical Freundlich adsorption coefficient and n fads is a linearity factor.
The Langmuir model assumes that adsorption of the solute on a solid takes place through formation of a monolayer. It considers ideal solutions where adsorption is independent of the coverage of the surface; the sites on which the adsorption takes place are uniformly distributed on the adsorbent surface; energy due to the interaction is equal at all the interaction sites; and a fixed number of exchange sites exists. The Langmuir model is applicable to L-type isotherms, characterized by a high affinity between adsorbate and adsorbent, commonly implying a chemical adsorption with the formation of inner sphere complexes [53]. The Langmuir model is described by: where C s and C Se have the same meaning as in the Freundlich model, C m-cal is the maximum adsorption capacity (mg g −1 ) calculated by the model and K L is the empirical affinity Langmuir coefficient (L mg −1 ). The combined Langmuir-Freundlich model (Sips model) [54] is expressed by: where K d is the apparent equilibrium affinity constant for heterogeneous solids, which includes k 2 /k 1 contributions; k 1 is the adsorption rate constant; k 2 is the dissociation rate constant; and "n" is the Langmuir-Freundlich heterogeneity parameter. This "n" value can also be employed as an empirical coefficient, representing the type and extent of cooperativity present in the binding interaction [52].

Turbidimetric Analysis of Sulfate
1 mL of sulfate-containing samples was placed in an Eppendorf tube with 0.05 mL of conditioning solution (prepared by mixing 30 mL of concentrated HCl, 300 mL of double distilled water, 100 mL of ethanol 96%, 75 g of NaCl, and 50 mL of glycerol). Then 0.01 g of BaCl 2 crystals were added and the mix was stirred continuously for 1 min. Immediately after that, sample was transferred into a cuvette and absorbance was measured at 500 nm of wavelength [55] using a Thermo Electron Spectronic Helios Alpha Beta UV-Visible double beam spectrophotometer (USA). In parallel, a calibration curve (sulfate concentration 0-40 mg mL −1 ) was prepared following same procedure as described for samples. Sulfate concentration was determined using the calibration curve equation [55].

Preparation of Modified Zeolite
Modified zeolite (MZ) synthesis was optimized through measurement of zeta potential across increasing quantities of SC-581. This optimization is shown in Figure 1. pH values for each experimental point were on average 6.6 ± 0.1.

Preparation of Modified Zeolite
Modified zeolite (MZ) synthesis was optimized through measurement of zeta potential across increasing quantities of SC-581. This optimization is shown in Figure 1. pH values for each experimental point were on average 6.6 ± 0.1. Optimization of MZ synthesis revealed that low amounts of polycation effectively change the surface charge of NZ. As shown in Figure 1, 25 mL of a 0.6 mg L −1 of SC-581 effectively neutralizes the negative charge of 1 g of NZ. As zeta potential reached a plateau of 45.2 mV after the addition of 4.8 mg L −1 of the polycation, this concentration was employed for preparation of the composite between NZ and SC-581 to produce MZ. IEP determination for NZ and MZ is presented in Figure 2. Optimization of MZ synthesis revealed that low amounts of polycation effectively change the surface charge of NZ. As shown in Figure 1, 25 mL of a 0.6 mg L −1 of SC-581 effectively neutralizes the negative charge of 1 g of NZ. As zeta potential reached a plateau of 45.2 mV after the addition of 4.8 mg L −1 of the polycation, this concentration was employed for preparation of the composite between NZ and SC-581 to produce MZ. IEP determination for NZ and MZ is presented in Figure 2.

Preparation of Modified Zeolite
Modified zeolite (MZ) synthesis was optimized through measurement of zeta potential across increasing quantities of SC-581. This optimization is shown in Figure 1. pH values for each experimental point were on average 6.6 ± 0.1. Optimization of MZ synthesis revealed that low amounts of polycation effectively change the surface charge of NZ. As shown in Figure 1, 25 mL of a 0.6 mg L −1 of SC-581 effectively neutralizes the negative charge of 1 g of NZ. As zeta potential reached a plateau of 45.2 mV after the addition of 4.8 mg L −1 of the polycation, this concentration was employed for preparation of the composite between NZ and SC-581 to produce MZ. IEP determination for NZ and MZ is presented in Figure 2. At first glance, these adsorbents exhibit marked variation in the curve of the slope. NZ develops negative surface charge at pH values greater than 3.1 (NZ isoelectric point), while MZ displays positive surface charge from pH values greater than its IEP (4.2). Interestingly, the curve for IEP determination of MZ always maintains a positive slope. When determining IEP through electrophoretic mobility studies, pH elevation typically potentiates negative charge over adsorbents. However, if multivalent cationic species (divalent: M 2+ ; polyvalent: M n+ , with n > 1) are abundant, the formation of surface complexes may result in the development of positive charge with increasing pH. Such behavior has been observed before for different systems. Escudey and Gil-Llambias [56] reported positive slope curves when determining the IEP of MoO 3 /γ-Al 2 O 3 catalyst if divalent cations (Ca 2+ and Ba 2+ ) were present in the background electrolyte. More recently, Arancibia-Miranda et al. [57] studied the effect of cations on the IEP of imogolite samples, observing a positive slope for zeta potential vs. pH curves in the presence of divalent cations (Ca 2+ and Mg 2+ ). In the interaction of SC-581 with NZ, the polymer acts analogously to multivalent cations. SC-581 consists of a co-polymer between neutral and cationic polyacrylamide. Smets and Hesbain [58] studied the effect of pH on the structure and stability of the polyacrylamide polymer, proving that its structure was extremely sensitive to low pH conditions. This work demonstrated that the polymer breaks down into its monomeric components at pH values less than 5.0, with lower pH values accelerating this process. Moreover, the product sheet for Superfloc C-581 [34,35] states the adequate working pH window to be 4.0-7.0. Indeed, at pH values below 4.0, SC-581 polycation structures break down into their components (2-chloromethyl oxirane, ethane-1,2-diamine and N-methylmethanamine) [59] which may still form multivalent oligomeric species. These polymer fragments may stablish surface complexes over the NZ surface, elevating surface charge to positive values with increasing pH. Therefore, sulfate adsorption experiments should be carried out at pH values ≥ 4.5 to ascertain that zeolite remains effectively modified by the polycation. In this pH range, NZ and MZ adsorbents should present negative and positive surface charge, respectively ( Figure 2). Regarding external surface area (BET) and pore size for the composite, it can be presumed that was maintained nearly the same as those for the natural zeolite since Dickson et al. (2020) [30] reported that zeolite containing similar Si/Al ratios did not modify these properties greatly after being modified with cationic surfactant hexadecyltrimethylammonium bromide. Also, if the hydrated radii of sulfate are reported to be 3.79 Å [60], then the interaction with the adsorbents may take place in the external surface or inside the adsorbent pores; sulfate would probably lose its solvated water molecules before. Finally, regarding the interaction between sulfate and the composite, according to zeta potential vs. pH diagrams at pH values of 4.2 or higher, the positive charge of the composite ensures the electrostatic interaction between sulfate and the solid, explaining in part the interaction.

Adsorption
The effect of ionic strength on sulfate adsorption is presented in Figure 3.  Fitted curves for ionic strength of 0.001 and 0.050 mol L −1 of KCl are quite similar. Both curves fit the Langmuir adsorption model well with a maximum of adsorption (C m-cal ) of 1.421 and 1.439 mg g −1 , respectively. Interestingly, increasing the KCl concentration to 0.100 mol L −1 changed the adsorption isotherm curve pattern from L-type to S-type. The L-type isotherm curve describes a high-affinity adsorption between sulfate and NZ, indicating chemisorption, while the S-type isotherm describes an interaction between sulfate species over the NZ and/or the adsorption of sulfate on NZ through solution ligands [47]. The difference in the adsorption curve pattern indicates that the sulfate adsorption process is modified as ionic strength increases, thus, maintaining low concentrations of KCl would be preferential to avoid shifting the adsorption mechanism. This could be explained by the formation of outer sphere complexes, in which electrostatic forces are responsible for adsorption. Therefore, these complexes may be easily affected by variations in ionic strength [61,62]. Consequently, it is possible that at KCl concentrations higher than 0.050 mol L −1 , chloride ions begin to compete with sulfate for adsorption sites and consequently reduce sulfate adsorption.

Adsorption Kinetics
Adsorption kinetics experiments determine both adsorption rates, the time needed to reach equilibrium, the type of adsorption taking place as well as the route and mechanism of the reaction [63]. Results from adsorption kinetics experiments for sulfate on NZ and MZ are in Figure 4.  Adsorbents became saturated after 10-12 h of equilibration time (Figure 4). Analyzing the plateaus, maximum adsorption capacity for kinetic studies on MZ is almost 2.2 times higher than on NZ. This preliminary result indicates that SC-581 effectively modified the surface of NZ to enhance sulfate adsorption capacity. Fitting parameters for pseudo-first-and pseudo-second order kinetic models are presented in Table 2. For both adsorbents, the pseudo-second order model better fits the kinetic data according to the correlation coefficient (R 2 ). Other authors have made the same observation for sulfate adsorption on modified aluminosilicate-based adsorbents with increased positive surface charge, resulting in the enhancement of their anion adsorption capacities. Chen and Liu [64] determined pseudo-second order kinetics when studying sulfate adsorption over natural zeolite modified with cationic surfactant cetyltrimethylammonium bromide (CTAB), reaching saturation after only 2 h. Moreover, Matusik [65] determined pseudo-second order kinetics for the adsorption of sulfate,  Adsorbents became saturated after 10-12 h of equilibration time (Figure 4). Analyzing the plateaus, maximum adsorption capacity for kinetic studies on MZ is almost 2.2 times higher than on NZ. This preliminary result indicates that SC-581 effectively modified the surface of NZ to enhance sulfate adsorption capacity. Fitting parameters for pseudo-firstand pseudo-second order kinetic models are presented in Table 2. For both adsorbents, the pseudo-second order model better fits the kinetic data according to the correlation coefficient (R 2 ). Other authors have made the same observation for sulfate adsorption on modified aluminosilicate-based adsorbents with increased positive surface charge, resulting in the enhancement of their anion adsorption capacities. Chen and Liu [64] determined pseudo-second order kinetics when studying sulfate adsorption over natural zeolite modified with cationic surfactant cetyltrimethylammonium bromide (CTAB), reaching saturation after only 2 h. Moreover, Matusik [65] determined pseudo-second order kinetics for the adsorption of sulfate, nitrate, orthophosphate, and arsenate onto N-methyl-diethanolamine-grafted kaolinite samples, achieving equilibrium time swiftly after 20 min. In the present work, adjustment parameters reveal that MZ adsorbs sulfate from solution 3.7 times faster (k 2 ) and displays 1.7 times greater adsorption capacity (C m ) compared to NZ. Therefore, SC-581 modification of NZ effectively increases adsorption rate of sulfate. To further understand adsorption dynamics, we produced adsorption isotherms by exposing the adsorbents to different concentrations of sulfate until equilibrium is reached (after 12 h). The adsorption isotherms are in Figure 5. As Figure 5 illustrates, the amount of adsorbed sulfate, Cs, increases rapidly when the concentration of sulfate at equilibrium, CSe, increases. As experiments were carried out at pH 4.8 ± 0.1, these data indicate that the composite structure of MZ was not threatened after the adsorption process. The isotherm shape of NZ resembles L-type isotherms, which represent a relatively strong adsorption interaction [53]. On the other hand, the MZ shape is more similar to H-type isotherms, implying a stronger interaction between analyte and adsorbents [53]. Experimental sulfate maximum adsorption capacity on NZ and MZ over is about 1.5 and 3.1 mg g −1 , respectively. Therefore, as expected from the IEP and adsorption kinetic data, sulfate adsorption capacities increased after modification with SC-581. Comparable results have been attained by other authors, including 2.5 mg g −1 of sulfate maximum adsorption on hexadecyltrimethylammonium-modified zeolite [66]. Moreover, maximum sulfate adsorption of 5.02 mg g −1 was reported for methyldiethanolamine-grafted kaolinite samples by Matusik [65]. Other studies have achieved higher maximum adsorption capacities, but with different experimental conditions. For instance, Oliveira and Rubio [67] modified zeolite samples (average diameter 25.4 μm) to produce a flocculated barium-zeolite composite, reporting sulfate maximum adsorption as high as 54 mg g −1 . Despite achieving this impressive level of sulfate adsorption, this approach produces highly toxic barium residues. On the other hand, Chen and Liu [64] made use of zeolite (sieved with 75 μm standard sieves) modified with CTAB, reaching 38 mg g −1 of sulfate maximum adsorption, but with isotherm construction carried out at 40 °C. Therefore, it is critical to consider experimental details when comparing results across different studies. To obtain further insight into these systems, Freundlich (F), Langmuir (L) and Langmuir-Freundlich (LF) models were applied to determine which As Figure 5 illustrates, the amount of adsorbed sulfate, C s , increases rapidly when the concentration of sulfate at equilibrium, C Se , increases. As experiments were carried out at pH 4.8 ± 0.1, these data indicate that the composite structure of MZ was not threatened after the adsorption process. The isotherm shape of NZ resembles L-type isotherms, which represent a relatively strong adsorption interaction [53]. On the other hand, the MZ shape is more similar to H-type isotherms, implying a stronger interaction between analyte and adsorbents [53]. Experimental sulfate maximum adsorption capacity on NZ and MZ over is about 1.5 and 3.1 mg g −1 , respectively. Therefore, as expected from the IEP and adsorption kinetic data, sulfate adsorption capacities increased after modification with SC-581. Comparable results have been attained by other authors, including 2.5 mg g −1 of sulfate maximum adsorption on hexadecyltrimethylammonium-modified zeolite [66]. Moreover, maximum sulfate adsorption of 5.02 mg g −1 was reported for methyl-diethanolamine-grafted kaolinite samples by Matusik [65]. Other studies have achieved higher maximum adsorption capacities, but with different experimental conditions. For instance, Oliveira and Rubio [67] modified zeolite samples (average diameter 25.4 µm) to produce a flocculated bariumzeolite composite, reporting sulfate maximum adsorption as high as 54 mg g −1 . Despite achieving this impressive level of sulfate adsorption, this approach produces highly toxic barium residues. On the other hand, Chen and Liu [64] made use of zeolite (sieved with 75 µm standard sieves) modified with CTAB, reaching 38 mg g −1 of sulfate maximum adsorption, but with isotherm construction carried out at 40 • C. Therefore, it is critical to consider experimental details when comparing results across different studies. To obtain further insight into these systems, Freundlich (F), Langmuir (L) and Langmuir-Freundlich (LF) models were applied to determine which model better represented the adsorption mechanism involved ( Figure 5). Fitting parameters are presented in Table 3. The behavior of the isotherms fit all three adsorption models, as evidenced by the marginal differences between correlation coefficients. This indicates that the relevant adsorption processes do not belong to a single type of process (pure physical or chemical adsorption, for instance). When comparing the R 2 values, the Langmuir-Freundlich model (LF) best fit the adsorption isotherm data. However, maximum adsorption capacity calculated through Langmuir model of 3.08 mg g −1 is closer to the experimental value of 3.1 mg g −1 , and the model tendency at higher C Se values seems more realistic. Therefore, there are good arguments for selecting both models for the description. In this sense, most sulfate adsorption isotherms adjustment reports found in the literature using similar adsorbents agree that the Langmuir model better represented the experimental data [65][66][67] compared to the Freundlich model [64]. Authors do not necessarily attempt to use the Langmuir-Freundlich model for fitting their adsorption data, in spite of its applicability to heterogeneous immobilization of ions on solid adsorbents, which involves both Langmuir and Freundlich processes [68].
Analyzing LF parameters, maximum adsorption capacity (C m-cal ) on MZ is two times higher than for NZ. As K d corresponds to the model dissociation constant, its value is a measure if the instability of the interaction between analyte and adsorbent; thus, adsorption interaction over MZ is about eight times stronger than on NZ. Finally, the LF coefficient number, n, represents the cooperativity extent of the process. Since n values are 0 < n < 1 on both adsorbents, this points to negative cooperativity, meaning that both zeolite-based adsorbents have low affinity for sulfate immobilization and the interaction sites act independently from each other [51,52]. The fact that the three models adjust well to the experimental data means that the adsorption process is a combination of chemical and physical processes. Taking into account that other studies [65] made on surfactant modified zeolites described a reversible adsorption process, it is reasonable to infer that sulfate adsorption on MZ is a partially reversible process. A partially reversible adsorption process is advantageous from the point of view that the composite will maintain the adsorbed sulfate retained while the sulfate concentration and pH values are constant. However, if these conditions change, the adsorbent will liberate sulfate, recovering the adsorbing surface activity for further remediation cycles.
Thus, as explained after the adsorption kinetic results, SO 4 2− adsorption capacities on the adsorbents in the present study seem to strongly relate to the surface charge differences observed at the working pH. MZ corresponds to a zeolite-based adsorbent of sulfate, acting faster and displaying superior adsorption capacities. It should be noted that its adsorption efficiency is lower than other reported materials for sulfate removal [69,70]. However, many of these materials are synthetic (mainly polymers) and not friendly to the environment, so in this sense the product developed in this work has a composition equivalent to soils and sediments with a marginal amount of polycation, reaching a maximum of 0.1 mg per gram of zeolite. MZ represents a low-cost, environmentally friendly, scalable material capable of adsorbing (remediating) sulfate from aqueous solutions with potential applications in the field of copper mining wastewater treatment.

Conclusions
The proposed composite between natural zeolite (NZ) samples and a commercialorigin polycationic flocculant (SC-581) resulted in a modified zeolite (MZ) surface with enhanced capacities for sulfate adsorption. Combination of 1 g of zeolite with 0.015 mg of of SC-581 in 25 mL is sufficient to achieve high charge reversion. Evaluation of the ionic strength effect on the adsorption process revealed that KCl concentration should not surpass 0.050 mol L −1 to avoid changes in the adsorption process. MZ exhibited double the sulfate adsorption capacities of NZ during adsorption kinetics and isotherms construction. Isotherm data were modeled with three different adsorption models achieving good fitting parameters in all cases. The best correlation is achieved by the Langmuir-Freundlich model, but also seems that the maximum adsorption capacity and the tendency is better described by the Langmuir model. LF parameters corroborated the enhanced adsorption capacities of MZ compared to NZ, with two times higher maximum adsorption capacities (C m-cal ) and an eight times lower dissociation constant (K d ). The differences observed between the adsorption capacity of MZ and NZ can be attributed to distinct surface charge development at working pH (around 5.0). Most of the adsorption takes place by an electrostatic interaction, defining a partially reversible sulfate adsorbent suitable for recycling. Collectively, these data support the further development of MZ for sulfate remediation technologies in wastewater. The methods are simple, low-cost and present low health and environmental risk. Future prospect for this technology is to optimize the adsorption of further anionic pollutants (phosphate, nitrate, arsenate, perchlorate, among others) considering different zeolite-type and polycation combinations.