Abstract

Perchlorate (ClO4-) in drinking water is one of the most serious issues relating to drinking water contamination. This work reports both the experimental and computational investigation of wastewater perchlorate adsorption on surfactant-modified activated carbon cloth (ACC) in order to design an efficient water treatment system. Three different types of cationic surfactants including tetramethylammonium (TMA), benzyltrimethylammonium (BTMA), and dodecyltrimethylammonium (DDTMA) were employed in this study, and the ClO4- adsorption capacity of various adsorbents was determined in batch experiments for 24 hours. Among all, DDTMA-ACC exhibits the highest ClO4- adsorption capacity (4.5 times that of pristine ACC), the highest monolayer coverage capacity () value (0.59 mmol/g), an equilibrium adsorption constant () value (35.90 L/mmol), and the lowest value (-26.0 KJ/mol), revealing that ClO4- has a stronger affinity for ACC functionalized by DDTMA than ACC functionalized by other surfactants. Molecular dynamics (MD) simulations were also conducted to illustrate the interactions among ClO4-, surfactants, and a model graphite surface underwater environment. The relative number density profiles and cumulative number studies demonstrate that TMA/BTMA tends to form single-layered packing near the graphite surface, while DDTMA prefers to pack in a multilayered fashion. The dense packing behavior of DDTMA on the graphite surface can enhance the ClO4- adsorption, showing good agreement with experimental results.

1. Introduction

Due to population growth and pollution, the global demand for clean and drinkable water has steadily increased. Globally, water contamination by anionic pollutants (NO3-, NO2-, CN-, S2-, and Cl-) [1, 2] and heavy metals (Ni, Zn, Pb, Cr, Cu, Hg, As, and Cd) [3] is increasing, and it is essential to investigate the appropriate strategies for improving water quality and reducing water scarcity. The presence of perchlorate (ClO4-) in drinking water is one of the most critical issues about the contamination of source water [4]. Ingestion of perchlorate-containing water or food products would limit thyroidal iodine uptake because the size and shape of perchlorate are comparable to those of iodide (I-) [5]. The low level of iodine intake thus decreases the production of thyroid hormone and causes a detrimental impact on the metabolism and growth of newborns, hypothyroidism, and neurological malfunction associated with other thyroid disorders [6, 7].

Physicochemical methods like adsorption, membrane filtration, and ion exchange have been extensively investigated to remove perchlorate from contaminated water [4]. Among different strategies, adsorption is the most widely used technology due to its combined advantages of low cost, high treatment efficiency, and easy operation. Additionally, activated carbon (AC) is widely used as an adsorbent due to its high surface area and tunable surface morphology. Therefore, significant efforts have been devoted to the surface engineering of AC in order to improve the adsorption capacity [8].

For example, Chen et al. [9] functionalized a bituminous granular AC with ammonia gas and obtained 4 times higher ClO4- adsorption compared to a virgin counterpart. They found that the increase in ClO4- adsorption capacity was closely correlated to the positive surface charge, which was enhanced by increasing nitrogen content during the ammonia functionalization. The study by Cui and Atkinson [10] fabricated the AC from glycerol and acid-mediated dehydration, polymerization, and carbonization. Sulfuric acid and phosphoric acid were applied during the AC preparation. Their results demonstrated that sulfuric acid-tailored AC exhibited an extremely high surface area of 2470 m2 g-1. However, the highest adsorption capacity was obtained in phosphoric acid-tailored AC, which was attributed to the higher amount of oxygen surface functional groups and the formation of surface complexation. Mahmudov et al. utilized quaternary ammonium salts as cationic surfactants to functionalize the AC and reported that its capability for ClO4- adsorption was significantly improved [11, 12]. The removal of ClO4- was found to be directly related to the surfactant chain length, which induced the formation of more positive surface charges as evidenced by zeta potential measurements.

Although surface modification on AC presents an effective approach to promoting its ClO4- adsorption capacity, most studies in the literature that discussed the interactions between substrate and surfactant as well as surfactant and ClO4- have been based on experimental observations such as kinetics and isotherm experiments that are inadequate for providing a fundamental understanding. To further support and guide experimental work, it is of interest to employ molecular simulations to examine the interaction among AC substrate, surfactants, and ClO4-. One of the most powerful techniques for modeling interaction behaviors between chemicals and substrates is molecular dynamics (MD) simulation [1318]. The MD simulation method is able to analyze the movement of molecules in a complex system at the atomic level by numerically solving Newton’s equation of motion [19]. Previous efforts have demonstrated that MD simulation can illustrate the adsorption state of nonionic and cationic surfactants on the low-rank coal surface [15, 20]. MD simulations can be used efficiently for the identification of adsorption sites, measurement of binding energies, elucidation of adsorption data, and optimization of adsorption processes [21]. In a recent study by Ouachtak et al. [22], MD simulation was used to investigate the adsorption behavior of methylene blue (MB) and crystal violet (CV) on green adsorbent montmorillonite @ activated carbon (Mt@AC). The RDF analysis and interaction energy calculations reveal that CV has a stronger affinity for the AC sheet and the montmorillonite framework than MB. Moreover, simulation results were found to be in good agreement with experimental adsorption isotherm measurements. MS simulations were used to evaluate the separation performance of porous C2N and g-C3N4 nanosheets for the separation of NO3- and NO2- pollutants from aqueous solutions. It was reported that the hydrophilic nature of these membranes (as confirmed by RDF between O2, H2, and N2 molecules) resulted in a threefold increase in self-diffusivity [23].

Attempts are being made to discover low-cost materials that can be used as green adsorbents from environmentally friendly resources without compromising efficiency. Low-cost green adsorbents can be made from a wide range of materials, such as clay, biomass waste, sawdust, and date palm charcoal [24]. Recently, orange peel waste has been utilized as a green biomass adsorbent for textile wastewater treatment and the production of biogas and biochar following regeneration and recycling. This adsorption system achieved a chemical oxygen demand (COD) of %, total dissolved solids (TDS) of %, and turbidity and color removal efficiencies of % and %, respectively [25]. Using NaOH chemically modified palm fiber (NaOH-CMDPF) waste biomass, Kumar et al. [26] achieved a high phenol removal efficiency from aquatic environments. It was reported that chemical modification with NaOH enhances phenol removal by 86% compared to raw date palm fiber (RDPF) with a removal efficiency of 81%. In another study, mesoporous-activated carbon (CSAC) prepared and activated by chitosan flakes (CS) and NaOH, respectively, was used for the adsorption of methylene blue (MB). A CSAC maximum adsorption capacity of 143.53 mg/L at 50°C was obtained by fitting to the Langmuir model [27]. A low-cost mesoporous zeolite-activated carbon (Z–AC) composite was used by Khanday et al. [28] as an effective adsorbent for dye removal. A maximum adsorption capacity of 285.71 mg/g of methylene blue was obtained at 50°C. Mesoporous activated carbon prepared from rattan (Lacosperma secundiflorum) hydrochar followed by NaOH activation and karanj (Pongamia pinnata) fruit hulls was used for adsorption of methylene blue (MB) with MB removal efficiency of 96% [29] and maximum adsorption capacities of 239.4 mg/gMB with hydrochar and karanj fruit hulls, respectively [30]. Khanday et al. [3133] utilize zeolite-hydroxyapatite oil palm ash (Z-HAP-AA), phosphoric acid-activated chitin (chitin-AC), and activated oil palm ash zeolite/chitosan composite (Z-AC/C) cross-linked beads for the adsorption of the antibiotic tetracycline, antibiotic cephalexin, and methylene blue (MB)/acid blue 29 (AB29), respectively. A maximum monolayer adsorption capacity of 186.09 mg/gMB, 245.19 mg/gMB, and 199.20 mg/gMB was obtained at 50°C for Z-HAP-AA, chitin-AC, and Z-AC/C, respectively. Marrakchi et al. [34] used cross-linked chitosan/sepiolite composite (CS50SP50) for methylene blue and reactive orange 16 (RO 16) adsorption. At an optimum adsorbent dosage of 0.2 g/100 mL, over a pH range of 3-11 and an initial dye concentration of 25-400 mg/l, a maximum monolayer adsorption capacity of 40.986 mg/g and 190.965 mg/g was calculated for CS50SP50 and RO 16, respectively. A pseudo-second-order kinetic was reported to have better adsorption than a pseudo-first-order kinetic.

The present study is aimed at exploring the adsorption behavior of cationic surfactants with distinct moieties on the surface of AC substrate and further investigating the correlation between surface functionalization and adsorption of anionic pollutants, i.e., ClO4-, from water. Most significantly, a joint experiment and MD simulation effort was carried out to study the interaction mechanisms among three species, which are rarely discussed in the literature. Different quaternary ammonium salts including tetramethylammonium (TMA), benzyltrimethylammonium (BTMA), and dodecyltrimethylammonium (DDTMA) were chosen as cationic surfactants in this study (Figure 1). Alkyl chains are bonded to the N center in TMA, which has a symmetric structure. In BTMA, a CH2 group connects the benzyl ring to the ammonium quaternary group, while in DDTMA’s ammonium quaternary group, it is attached to a 12-alkyl chain [35]. The examination of these three cationic surfactants through our integrated simulation-experimental approach will reveal the impacts of long alkyl chains and the presence of benzyl rings on adsorption behaviors onto the ACC surface, which provide a transformative impact on material design and enable the development of systems that greatly surpass the performance of conventional adsorbents for perchlorate removal in wastewater.

2. Experimental

2.1. Materials

All cationic surfactants used in this study were procured from Alfa Aesar Company (United States) and utilized without further purification. The aqueous solution was prepared using deionized water with a resistivity greater than 18 ΩM cm and distilled water. The standard perchlorate (ClO4-) solution was acquired from SPEX CertiPrep (USA). Stock solutions of sodium perchlorate (>98%, Sigma-Aldrich, USA) and sodium hydroxide (>97%, Fisher Scientific, USA) were prepared by individually dissolving metal salts into deionized water. In ion chromatography analysis, a 50% sodium hydroxide solution from Fisher Scientific (USA) was used as the mobile phase solution.

2.2. Surfactant Adsorption

Surfactant stock solutions were prepared at predetermined concentrations of TMA, BTMA, and DDTMA in distilled deionized water. Batch adsorption experiments were conducted in a 50 mL amber glass bottle. To achieve a final ionic concentration of 10-2 M, 20 mL of deionized water, followed by 0.5 mL of a 1.0 M NaCl solution, was added to each bottle. Then, 0.1 g of ACC was added to each bottle and kept hydrating for one hour. To attain the desired concentrations, a predetermined quantity of surfactant stock solution was added, and pH was adjusted using 0.1 M, HCl, and NaOH solutions while maintaining the final sample volume at 50 mL. The adsorption reaction reached equilibrium within 8 h. The samples were then agitated on the reciprocal shaker for 10 h, and the residual surfactant content was determined in terms of total organic carbon content (TOC) by filtering 10 mL of solution aliquots with a Millex 5-m PVDF syringe-type filter (Millipore P/N SLSV025LS). The surfactants adsorbed on the ACC surface were calculated by dividing the difference between the final and initial concentrations by the ACC concentration.

2.3. Perchlorate Adsorption

The functionalized ACC was prepared by dispersing 0.03 g of ACC in 50 mL of a predetermined ClO4- concentration in perchlorate adsorption experiments. After 24 h of adsorption, the ClO4- concentration was recorded periodically using an ion chromatograph (Dionex 500) by withdrawing small aliquots (2 mL) of perchlorate at varied time intervals. The functionalized ACC was then rinsed three times with distilled deionized water to remove superficially adsorbed surfactants and dried in the oven at 38°C for 24 h. The as-prepared samples were cooled down and then placed in a desiccator before use. Adsorptions of perchlorate on functionalized ACC were carried out in batch experiments using the same methodology as the adsorption of cationic surfactants. All the experiments with the adsorption of fluoride by either pristine or functionalized ACC were performed in duplicate to assure reproducibility, and the adsorption values shown in this work are the average of two samples at different conditions.

2.4. Characterization

The Brunauer-Emmett-Teller (BET) model is used to determine the specific surface areas of all samples from surface measurements performed with the ASAP 2020 (Micromeritics), and the pore size distribution was calculated from the desorption branch using the Barrett-Joyner-Halenda (BJH) technique.

3. Computational Methods

3.1. Adsorption Energies between Single Surfactant Molecule and Graphite Model

The adsorption energies between surfactant molecules and the ACC substrate are calculated using semiempirical quantum mechanics (SQM). The molecular adsorption is modeled as a barrier-less exothermic reaction, and the adsorption energy is obtained from calculating the binding energy using the following equation: where is the minimal potential energy of the complex in vacuum mum, and and represent the energies of the substrate and the surfactant separated at an infinite distance.

The calculation of adsorption energy was based on the procedure reported previously [37]. First, the graphite model went through energy minimization in a vacuum with the -coordinate of all graphite atoms fixed at zero. Then, a single surfactant molecule was placed at the center of the graphite model. After optimizing the geometry of the surfactant with all graphite atoms fixed in space, was obtained. and were then determined as single-point energies by separating the surfactant and graphite models at an infinite distance. The adsorption energies at the SQM level were calculated using the GFN2-xTB method [38] and a squared graphite model including 504 carbon atoms with hydrogen atoms to saturate the valence of edge carbons, as shown in Figure 2(a).

3.2. Simulations of Surfactant-Perchlorate Mixtures on a Graphite Surface

Molecular dynamics (MD) simulations were performed to understand the mechanism of different surfactant-ClO4- adsorption behaviors at the atomic level. The system setup is shown in Figures 2(b) and 2(c). The squared graphite surface is comprised of 1056 neutral carbon atoms with the dimensions  Å. Initially, positively charged surfactant molecules with random orientation were added within 5 Å above the graphite surface. Then the water molecules and counter-ions were added to the system to make the concentrations of 0.275/0.55/0.825 M. The actual numbers of each species are listed in Table 1. Initial structures for simulations were prepared using Packmol [39]. A liquid-vapor interface is formed ~50 Å above the graphite surface depending on the surfactant and concentration of the solution. Periodic boundary conditions were applied, but the -dimension was extended to 100 Å (see Figure 2(c)); therefore, only one liquid-solid interface is formed in each system.

All simulations were conducted in the NVT ensemble using NAMD 2.14 [40] with the CHARMM 36 and CGenFF force field parameters [41, 42]. Carbon atoms in the rigid graphite surface are modeled as neutral Lennard–Jones (LJ) spheres with the atom type CG2R61 (aromatic carbon) and are fixed throughout the simulations. We use TIP3P [43] water with SHAKE [44] to enforce the rigid water geometry. The particle mesh Ewald (PME) [45] method with a grid of 1.0 Å resolution and the Langevin method with a damping coefficient of 5 ps-1 were used to treat long-range electrostatic interactions and to fix the temperature at 300 K. The Lennard–Jones interactions were switched off gradually over the range of 10 and 11 Å. The equations of motion were integrated with a time step of 2.0 fs.

For each surfactant-ClO4- system and each concentration, we have at least 5 independent simulation runs with different initial structures and different seed numbers. Each simulation run is 10 ns with snapshots saved every 10 ps for further analyses. The first 5 ns of each run are treated as equilibration. The total sampling time for each system was between 25 and 50 ns.

4. Results and Discussion

4.1. Surfactant Adsorption on Activated Carbon Cloth

The adsorption isotherm of different surfactants on the activated carbon cloth (ACC) is provided in Figure 3. It shows that the adsorption capacity is significantly influenced by the chemical moieties of surfactants. All the isotherms in Figure 3 exhibit a sharp increase in adsorption followed by a single plateau that can be fitted to the Langmuir model [46] using the following equation: where is the amount of adsorbed surfactants (mmol/g), represents the monolayer coverage capacity (mmol/g), is the equilibrium constant (L/mmol), and is the equilibrium concentration in solution (mM). Table 2 enlists the fitted parameters of the Langmuir equation. The values of and are in the order of DDTMA > BTMA > TMA. Results show that extended alkyl chains and aromatic rings facilitate the adsorption of surfactants on the ACC, particularly for DDTMA. The higher suggests that the 12 alkyl chains in DDTMA can increase the affinity of surfactant towards the ACC surface, which is potentially contributed by the strong hydrophobicity of DDTMA [11, 47]. Moreover, a high value of is also obtained in DDTMA (Table 2). This result indicates that the DDTMA surfactants are adsorbed more strongly on the surface of ACC compared to BTMA and TMA, which supports the findings that hydrophobic forces are more dominant for surfactant adsorption.

Additionally, the adsorption-free energy, , was calculated based on using Eq. (3) [48, 49]. where and are the temperature and ideal gas constant, respectively, and the corresponding values of adsorption-free energy are reported in Table 2. Results demonstrated that all three surfactants provided a negative value, indicating a spontaneous adsorption process onto the ACC surface under experimental conditions. Moreover, DDTMA exhibited the most negative , revealing a stronger bonding tendency.

4.2. Perchlorate Adsorption on Functionalized Activated Carbon Cloth

In this work, the surfactant loading on functionalized ACCs for ClO4- adsorption is fixed at ca. 0.4 mmol/g. This specific mass loading is employed to prevent the influence of micelle aggregation during the surfactant loading process. The critical micelle concentration (CMC) for TMA [50], BTMA [51], and DTMA [52] is 5.4-5.7 mM, 110 mM, and 15 mM, respectively. It was reported that when the concentration of the surfactant is lower than its CMC, the surfactant molecules adsorbing onto the ACC surface exist primarily as monomers [53]. Therefore, the surfactant loading of 0.4 mmol/g on ACC presents an ideal platform to study the interactions among ACC, ClO4-, and various chemical moieties of cationic surfactants.

The ClO4- adsorption isotherms of pristine and functionalized ACCs are shown in Figure 4, and the corresponding fitted parameters for the Langmuir model are reported in Table 3. All three modified ACCs exhibit an increase in adsorption capacities compared to the pristine ACC. The maximum adsorption capacity of ClO4- is obtained in DDTMA-ACC, which is about 4.5 times higher than the performance of pristine ACC. Additionally, DDTMA-ACC presents the largest equilibrium adsorption constant and the most negative , indicating a stronger affinity for ClO4- towards ACC functionalized by DDTMA than that by other surfactants.

To explore the influence of surface modification on the performance of ClO4- adsorption, the surface properties of pristine and functionalized ACCs were examined. The BET-specific surface area is detailed in Table 3. The functionalization process is found to slightly decrease the surface area of ACC. The nitrogen adsorption–desorption isotherms of all ACCs are demonstrated in Figure S1. The obtained hysteresis loop exhibits a type-I pattern, suggesting the ACCs are mainly microporous solids [54]. The pore volume of samples decreased from 0.052 cm3/g for pristine ACC to 0.04 cm3/g for BTMA-ACC and DDTMA-ACC. The reduction of pore volume along with decreased surface area suggests that the surfactants were successfully adsorbed on the ACC surface but in different packing behaviors, thereby influencing their pore structure.

The characterization results clearly demonstrate that the functionalization of ACCs with various surfactants changes their surface properties including specific surface area, pore volume, and pore diameter. However, these features, such as reduced surface area, are difficult to explain the promoted adsorption capability of modified ACC and to reflect the impact of different surfactants on the adsorption performance [55]. To better understand the adsorption mechanisms of ClO4- in these ACC materials, MD simulations are employed to investigate the interactions among ClO4-, ACC, and various surfactants used in this work.

4.3. Adsorption Behavior of Surfactants Using In Silico Methods
4.3.1. Adsorption Energy between a Single Surfactant and the Graphite Surface

Initially, the adsorption energies among surfactants, perchlorate, and a model graphite surface were analyzed, and the results for TMA, BTMA, and DDTMA are shown in Table 4. The adsorption energies computed using DFT, force field, and GFN2-xTB techniques for other molecules in the supplementary information (Table S1) are benchmarks. For short-chain alkanes (methane and ethane), all three methods are within 1 kcal/mol differences from the experimental values. However, for aromatic compounds (benzene, naphthalene, and o-dichlorobenzene), differences of all methods are still within 2 kcal/mol, with the force field method showing a larger deviation from the experimental values. With regard to the gas/solvent molecules (methanol, ammonia, and water), all methods significantly underestimate the adsorption energies. The possible reasons have been discussed elsewhere [32] and are out of the scope of this manuscript. Based on the values reported for these molecules, we conclude that the performance of the GFN2-xTB SQM method has a good balance between efficiency and accuracy for adsorption energy calculations.

The calculated adsorption energies for TMA, BTMA, and DDTMA qualitatively reproduce the trend of the experimentally fitted in Table 2. For BTMA, the optimized structure shows an orientational preference for the aromatic ring to be nearly parallel to the graphite surface with a ~90 kcal/mol difference from the dihedral angle scanning of the CH2-aromatic ring bond (shown in Figure S2). As for DDTMA, there are more conformations than TMA and BTMA due to its longer chain length. The conformation with the lowest energy is extended with a distance of 15.3 Å between the nitrogen atom and terminal CH3 (as shown in Figure 1(c)). We use the vector from nitrogen to the terminal CH3 () to define the angle with the graphite surface and examine the interaction energies at different angles, as shown in Figure S3. The interaction energy shows a minimum when DDTMA is parallel to the graphite surface.

In practice, we rarely have only one surfactant molecule interacting with the graphite surface as in adsorption energy calculations. While multiple surfactant molecules have interactions with each other, they also interact with the perchlorate, which may lead to different adsorption behaviors.

4.3.2. Number Density Profile of Perchlorate and Surfactant Ions with the Graphite Surface

We use atomistic molecular dynamics simulations to mimic multibody interactions and analyze the adsorption behaviors of different surfactant-perchlorate systems. The number density profiles of cations/surfactants and perchlorate along the -dimension relative to the graphite surface are shown in Figure 5. The number density is defined as where and represent the length of - and -dimensions in the simulations. represents the unit length in -dimension of each volume slab, with denoting the number of species within this slab. In each panel, the total area under each curve/species equals 1.

We first focus on the number density profiles of Na+/surfactants as . We use the positions of the nitrogen atoms in the surfactant molecules to represent the corresponding positions in these calculations since the nitrogen atoms carry most of the positive charges. Na+ shows no preference near the graphite surface in all concentrations (Figures 5(a)–5(c)). TMA shows two peaks at and 7.5 Å (Figures 5(d)–5(f)), with the peaks both becoming more evident as the concentration increases. BTMA shows a strong first peak at  Å, with a weak second peak only shown in the highest concentration (Figures 5(g)–5(i)). In the case of DDTMA, there is only one strong peak at 0.275 M (Figure 5(j)); as the concentration elevates, the second peak becomes more evident with an extremely long tail (Figures 5(k) and 5(l)).

With regard to the number density profiles of perchlorate , when the counter-ion is Na+, ClO4- still shows some surface activity with a peak at  Å. When the surfactant is TMA, the ClO4- adsorption peaks are enhanced relative to the cases of Na-ClO4 at all concentrations. This adsorption peak enhancement is more evident in the BTMA-ClO4 systems. In the case of DDTMA as the surfactant, the ClO4- adsorption peak is slightly enhanced compared to the Na-ClO4 system at low concentration; as the concentration increases, the ClO4- adsorption peak becomes even weaker than the Na-ClO4 system. While other surfactants such as TMA/BTMA tend to show strong peaks/single-layered packing behavior near the graphite surface, DDTMA shows relatively flattened (or multilayered, contrary to the single-layered) adsorption behaviors from the graphite surface. We further investigate this effect by plotting cumulative numbers of ClO4- and Na+/surfactants as functions of relative -positions in Figure 6. The cumulative number is defined as follows: where is the total number of anion-cation pairs at each concentration; therefore, when equals infinity. For all species in our systems, reaches at  Å. We use  Å as the distance cutoff to inspect each , as listed in Table 5. When  M, and in both BTMA-ClO4 and DDTMA-ClO4 systems are similar (~15 and ~20). At  M, we start to see the amount of DDTMA () to be significantly larger than BTMA (), while both are ~32. Both and in DDTMA-ClO4 are larger than and in BTMA-ClO4 at the highest concentration. Despite being the largest surfactant in this study, DDTMA molecules show the stronger capability to pack themselves within identical volumes ( Å in our case). For surfactants like TMA/BTMA, they can only form a single layer near the graphite surface, while the rest of the surfactants are further away from the surface, which consequently reduces the adsorption of perchlorate anions. On the other hand, most DDTMA molecules reside near the graphite surface, thus inducing more perchlorate anion adsorption.

5. Conclusions

Activated carbon cloth is successfully functionalized using surfactant adsorption for perchlorate removal. A constant surfactant loading of 0.4 mmol/g is used to prevent micelle aggregation and to analyze the interaction among ACC, ClO4-, and various cationic surfactants. The adsorption of the surfactants on the surface of the ACC is confirmed by a decrease in pore volume and an increase in surface area. Regardless of the nature of the surfactants, high functionality and perchlorate removal efficiency compared with pristine ACC are obtained. Of the most significant, DDTMA-ACC demonstrates about 350% improvement in the ClO4- adsorption compared to pristine ACC. DDTMA-ACC also outperforms BTMA-ACC and TMA-ACC with a 30% and 120% increase in , respectively. Further, MD simulation studies reveal that the enhanced ClO4- adsorption capability of DDTMA-ACC is attributable to the closely packed or multilayered packing features of DDTMA on the ACC surface. It is found that minimal interaction energy occurs when DDTMA is parallel with the model graphite surface. Relative number density profiles and cumulative number studies further demonstrate that the majority of DDTMA molecules reside near the graphite surface, whereas TMA/BTMA can only form a single layer, which deteriorates the adsorption performance. Our joint experiment-simulation study offers a new route to elucidate the underlying mechanism of surfactant-based surface modification and enable a transformative impact on material design to surpass the performance of conventional perchlorate adsorbents for wastewater treatment.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

The authors acknowledge the financial support from the National Science and Technology Council of Taiwan with grant numbers 109-2222-E-131-002-MY3, 110-2222-E-182-002-MY3, 112-2628-E-182-001-MY3, and 112-2221-E-131-007. W.-F. K. and C.-L. C. also express gratitude for the financial support provided by Chang Gung University (URRPD2N0041).

Supplementary Materials

Figure S1: nitrogen adsorption–desorption isotherms of various adsorbents. Figure S2: relative energy profile for different orientations between BTMA and graphite surface. Figure S3: relative energy profile for different orientations of DDTMA. Table S1: calculated adsorption energies to graphite using CGenFF force field, DFT, and semiempirical GFN2-xTB methods. (Supplementary Materials)