Density functional theory study of hydrophobic zeolites for the removal of triclosan from aqueous solution

The chlorinated biphenyl ether triclosan (TCS), used as a disinfectant in health care settings and in various personal care products, is an emerging organic contaminant of significant concern. Adsorption-based methods have been proposed as one potential pathway for the removal of TCS from wastewaters. Hydrophobic high-silica zeolites could constitute suitable adsorbent materials for such applications. In order to gauge the impact of pore size, topology, and framework composition, the adsorption of TCS in six different all-silica zeolites (AFI, BEA, CFI, FAU, IFR, MOR frameworks) and two highly siliceous protonated zeolites (H-FAU, H-MOR) was investigated using dispersion-corrected density functional theory (DFT). While pore size was found to affect the interaction strength, the rather flexible TCS molecule can adjust to different pore shapes, resulting in very similar adsorption energies for most all-silica zeolites. Although the interaction with TCS is enhanced in protonated zeolites, the affinity towards water increases even more. In DFT-based molecular dynamics simulations of TCS and water co-adsorption, H 2 O molecules quickly replace TCS in the vicinity of the framework protons, deprotonating the framework and forming positively charged clusters. In addition to delivering atomic-level insights into TCS adsorption, the calculations indicate that a fine-tuning of pore size with a concurrent maximization of hydrophobicity should constitute a promising strategy to develop optimized zeolite adsorbents for TCS removal.


Introduction
8][9] Since potential risks to human health, including reproductive problems and enhanced risk of asthma and allergies, have been identified, the use of TCS in consumer products has been restricted in several countries. 7,9Given the persistence of TCS, which is stable to hydrolysis at near-neutral pH values.1][12][13] Incomplete removal of TCS by wastewater treatment facilities results in a significant discharge into the environment: In a global survey by the German Environment Agency (Umweltbundesamt), triclosan was listed among the 20 PPCP contaminants that were most frequently detected in surface water, groundwater, or drinking water. 14In the European Union, it has recently been proposed to add triclosan to a list of priority water pollutants requiring more stringent control measures (EU legislative procedure 2022/0344/COD).
A number of possible effects contribute to the environmental hazard potential of TCS.These include acute toxicity to aquatic organisms (however, at concentrations that are typically not reached in real-world scenarios 15 ), chronic toxicity, and endocrine-disrupting effects. 5reover, the antimicrobial TCS may contribute to selection processes that are responsible for the development of antibiotic-resistant bacteria. 5,8,16Due to the frequent presence of relatively low levels of antibiotics and disinfectants in wastewaters and wastewater treatment plants (WTTPs), these have been identified as "hot spots" for the development of microbial resistance. 17Hence, the frequent presence of TCS in such environments should be a cause for concern, not only due to its potential environmental impact, but also with regard to possible, more indirect consequences for human health.Non-negligible amounts of TCS can also be found in soils amended with biosolids from wastewater treatment plants.Here, TCS and its degradation product methyltriclosan, formed during the wastewater treatment process, may accumulate in terrestrial organisms like earthworms. 18Other relevant degradation products include highly toxic chlorinated dioxins and chlorophenols, which are produced through photolytic degradation of TCS. 5,19arting with the work of McAvoy et al., 10 several authors have measured TCS concentrations in WWTP influents and effluents, and the results have been compiled in review articles 3,12 and databases. 14,20While concentrations in effluents remain usually below 1 µg L -1 , a few studies reported values above 10 µg L -1 . 11The reported removal efficiencies of conventional WWTPs also vary widely, ranging from removal rates below 50% to essentially complete removal. 12,13,21veral advanced treatment methods that are under consideration for PPCP removal in general have also been studied for TCS removal in particular, including adsorption-based methods, advanced oxidation processes, and biodegradation.The advantages and drawbacks of individual methods have been discussed in recent reviews. 22,23e to the lipophilicity of TCS (log ( ) of 4.76 24 ), hydrophobic adsorbents appear as the most promising materials for the adsorption-based removal of TCS from aqueous solution.
Activated carbons 25,26 and other carbon-based adsorbents like biochars 27 have been shown to exhibit high TCS removal efficiencies.Although hydrophobic high-silica zeolites are more expensive to produce than these materials, they possess some advantages, such as high thermal and chemical stability, allowing for regeneration under fairly harsh conditions, 28 reduced co-adsorption of natural organic matter, 29 and high affinities towards species that fit well into the ordered pore system of the zeolite. 30ssner et al. studied the removal of a "cocktail" of 25 organic contaminants, including TCS with an initial concentration of about 0.6 µg L -1 , from spiked lake water by two high-silica zeolites. 25Whereas high-silica mordenite (MOR framework type 31 ) showed essentially complete removal of TCS (along with various other species), FAU-type zeolite Y removed only about 45%.Since the FAU-type sample should be more hydrophobic due to its higher Si/Al ratio (Si/Al ≈ 400 as compared to Si/Al ≈ 110 for MOR), this difference cannot be attributed to the zeolite composition, but is more likely related to differences in pore size and topology: Both zeolites have pore apertures outlined by 12-membered rings (12MRs) of tetrahedrally coordinated atoms, but the MOR framework possesses 1D channels with a diameter of about 6.5 Å, whereas the pore system of FAU consists of large supercages (diameter about 11 Å) connected by 12MR windows.More recently, Jiang et al. measured TCS adsorption isotherms (aqueous solution, concentrations ranging from ~200 to ~6,000 µg L -1 ) using zeolites with four different framework types: 32 In addition to FAU-and MOR-type zeolites, zeolite beta (BEA framework, 12MR pore openings) and ZSM-5 (MFI framework, 10MR pore openings) were considered, and samples with different Si/Al ratios were compared for FAU, BEA, and MFI.
Only negligible uptake was observed for MFI-type samples, indicating that the sizeable TCS molecules cannot diffuse through the 10MR windows.Among the other three zeolites, the highest TCS uptake of 378 mg g -1 was observed for the most Si-rich FAU sample (Si/Al ≈ 400), and both maximum uptake and affinity towards TCS, determined via a Langmuir-Freundlich fit, decreased with decreasing Si/Al ratio.While still being appreciable, uptakes and affinities of BEA-type zeolites (Si/Al ratios from ~80 to ~300) and a MOR-type sample (Si/Al ≈ 110) were lower.In addition to these investigations of synthetic zeolites, TCS adsorption experiments were also performed on natural zeolites (clinoptilolites) modified with organic surfactants, which could be cheaper to produce. 33,34It is, however, worth noting that maximal TCS loadings and affinities towards TCS reported in those works were considerably lower than for the bestperforming high-silica zeolites studied by Jiang et al. 32 Computational modelling has become a very widely used tool in zeolite science. 35Atomistic simulations at different levels of theory have been exploited for various purposes, e.g., to develop a deeper understanding of experimental observations, or to predict properties that were so far not characterized experimentally.With regard to PPCP removal, force field (FF) simulations were employed to compute the adsorption energies of 21 contaminants, including TCS, in all-silica MOR and FAU. 36Despite the simplistic description of the interatomic interactions, a good correlation between FF adsorption energies and experimental removal efficiencies reported by Rossner et al. 25 was observed.This indicates that such simulations could be employed to predict zeolite adsorbents having a high affinity towards a given PPCP prior to an experimental characterization.Subsequently, the adsorption energies were recomputed using periodic electronic structure calculations in the framework of dispersioncorrected density functional theory (DFT) with the purpose of identifying suitable DFT approaches for the study of functional organic molecules in zeolites. 37Several earlier DFT studies employed cluster models cut out from the zeolite structure in order to reduce the computational expense, with examples including investigations of sulfonamide antibiotics in all-silica FAU 38,39 and of salicylic acid, carbamazepine, and ciprofloxacin in cation-exchanged FAU. 40With ongoing hardware and software developments, periodic calculations are now routinely feasible for unit cells of typical zeolites.A recent DFT investigation studied the adsorption of carbamazepine (CBZ) in eleven all-silica zeolites. 41In addition to the optimization of CBZ@Zeo adsorption complexes, DFT-based ab initio molecular dynamics (AIMD) simulations were carried out for selected zeolites to analyze the influence of temperature.
Moreover, the role of guest-guest interactions was evaluated through an investigation of adsorbed CBZ dimers.
The present work addresses the adsorption of TCS in zeolites, employing, by and large, a similar approach as this previous study.As for CBZ, the investigation is restricted to zeolites having 12MR or 14MR pore openings, as experimental findings indicate that zeolites with smaller pores do not adsorb appreciable amounts of TCS. 32A somewhat smaller number of purely siliceous zeolite is considered here, including six, rather than eleven, frameworks.Apart from studying a different guest molecule, the present study goes beyond the previous work in the following respects: First, preliminary FF-based simulations are used to corroborate that TCS can diffuse through the pores of these zeolites.Second, the focus is no longer exclusively on all-silica zeolites, but models of highly siliceous, protonated zeolites are also considered for FAU and MOR topologies.Third, the co-adsorption of water, neglected in the previous study, is investigated for selected systems.In addition to providing atomic-level insights into TCS adsorption, the calculations allow predictions how pore size, pore shape, and framework composition affect the performance of hydrophobic zeolites as adsorbents for TCS removal.

Structure of guest molecules
The molecular structure of TCS was taken from a previous X-ray diffraction study. 42In the crystal structure, the TCS hydroxyl groups simultaneously act as hydrogen bond donors and acceptors, resulting in one-dimensional chains of hydrogen-bonded TCS molecules.Moreover, the phenyl rings of neighboring molecules interact through π-π stacking interactions.Since it can be expected that the equilibrium molecular structure will be considerably different in the absence of intermolecular interactions, a conformational screening was carried out, using the Materials Studio (MS) Conformer module and the pcff force field 43 (see Electronic Supplementary Information [ESI], Table S1 for pcff parameters).Local energy minima were re-optimized using DFT calculations with the CP2K code (see subsection on DFT calculations for further details).In the lowest-energy conformer identified on the basis of these calculations, the hydroxyl group forms an intramolecular hydrogen bond to the central oxygen atom Oc (Figure 1, (H ••• O ) = 2.181 Å).A vibrational calculation confirmed the absence of imaginary frequencies.The DFT-calculated electrostatic potential (bottom panel of Figure 1) shows the comparatively strong positive polarization of the Hh proton of the hydroxyl group, which exceeds that of other positively polarized areas of the molecule (aromatic protons).This TCS conformer was used as reference for the calculation of adsorption energies and as starting point for the generation of TCS@Zeo adsorption complexes.
The structure of the H2O molecule was first optimized using pcff and then re-optimized with DFT, using the same settings as for TCS.As discussed in more detail in the Results section, there is little point in comparing adsorption energies obtained for single water molecules to those of one TCS molecule per cell due to the vast difference in molecular size.For the case of protonated zeolites, experimental results and prior DFT studies indicate that an adsorption of clusters of several H2O molecules, which are large enough to deprotonate the framework, is energetically favored over the adsorption of single H2O molecules at the framework protons. 44,45Therefore, calculations aimed at a comparison of TCS and water adsorption considered clusters of eight H2O molecules.To study energetic trends among (H2O)N clusters of different size, additional calculations were done for models with N = 1, 4, and 12. Structure figures were prepared using VESTA. 46

Zeolite structures
The structures of the six all-silica zeolites (BEA, FAU, AFI, IFR, MOR, CFI frameworks, represented by the framework type codes [FTC] assigned in the structure database of the International Zeolite Associations IZA 31 ) are the same as in previous work (Table 1). 41Sources of the experimental structure data and the corresponding cell parameters are given in the ESI (Table S2).As described in more detail in previous work, short AIMD simulations were run for all zeolite structures (NVT ensemble, T = 298 K, timestep 1 fs, total duration 12.5 ps).The average structure from the production stage of the AIMD simulation (10 ps) was then subjected to a symmetry search and DFT-optimized.atom and one charge-balancing proton per unit cell/1×1×3 supercell, respectively.Due to the low amount of framework Al atoms, the same (experimental) cell parameters as for the allsilica systems were used. 49,53As there is only one type of T site in FAU, only one Al position was considered.Structure models with protons attached to each of the four surrounding oxygen atoms were constructed and DFT-optimized.Among the four models, the H@O2 case was found to be lowest in energy (Table S3, the labelling of T and O atoms in FAU and MOR follows the IZA database 31 ).However, the H@O2 proton points across a 6MR in an inaccessible area of the structure (the 6MR links a double six-membered ring and a sodalite cage, both of which are inaccessible for most guest molecules).In the H@O1 model, which is only 5 kJ mol -1 higher in energy (second lowest), the proton points into a 12MR, thus being accessible to guests accommodated in the supercages that are linked by these windows.This is also the most occupied proton site according to experimental investigations on protonated FAU samples with lower Si/Al ratios. 56,57The H-FAU(H@O1) model was hence used in all further calculations.In MOR, the increased complexity of the framework, with 4 non-equivalent T sites in the aristotype, is exacerbated by the symmetry reduction to space group 2 2 2 in the structure model obtained from the AIMD simulations. 41To limit the number of possible Al+H arrangements, only those cases where the proton points (approximately) into the 12MR channels were considered: Under this constraint, the framework proton can be attached to the O1, O4, O7, or O10 atoms, and the Al atom, which should be at a neighboring T site, can be located at the T1, T2, or T4 sites.Due to the aforementioned symmetry reduction, some of these sites are split into up to 4 non-equivalent sites, resulting in a total of 22 distinct arrangements.Models for these arrangements were constructed in a 1×1×3 supercell of the MOR unit cell and DFT-optimized.All models are listed in the ESI (Table S4 and Figure S2).
Among them, the model dubbed H-MOR(Al@T1_4,H@O4_4) was found to be lowest in energy and used in all following calculations on H-MOR.A previous DFT study of protonated MOR also delivered a model with Al@T1, H@O4 as lowest-energy case, in line with this result. 58As for the all-silica models, AIMD simulations (NVT ensemble, T = 298 K, timestep 0.5 fs, total duration 12.5 ps) were run for H-FAU and H-MOR.The average structures from the 10 ps production part of the trajectory were optimized again.These calculations converged to the same minima as the initial optimizations.

Force field simulations
In order to generate adsorption complexes of TCS, water ((H2O)N), or both TCS and water in different zeolites, a combination of Monte Carlo (MC) and molecular dynamics (MD) simulations was used, employing parameters from the pcff force field. 43,59The typical protocol consisted of fixed-loading MC simulations using the MS Sorption module, from which five lowenergy snapshots were extracted for each combination of guest and zeolite.To enhance the sampling, MD-based simulated annealing runs using the MS Forcite module were then carried out, starting from the five snapshots, with each run consisting of 20 annealing cycles.The configurations from the end of each annealing cycle were optimized.Out of the 100 configurations obtained via this procedure, five low-energy configurations were selected as starting points for the DFT optimizations.These configurations were labelled as "guest@Zeo, ConfigX", where "guest" represents TCS or H2O, "Zeo" corresponds to the zeolite framework type, and "X" is a running index from 1 to 5, assigned according to the pcff total energy.For co-adsorption of TCS and H2O, labels of the format "guest2(guest1@Zeo)" were assigned, where "guest1" is the species that was added first, and "guest2" the species that was added subsequently.Detailed settings of the force field simulations are provided in the ESI.
Additional MD simulations were carried out to assess whether TCS can diffuse through the pores of the six all-silica zeolites.These simulations, which started from low-energy TCS@Zeo configurations, also used pcff parameters.Unlike in the MC and simulated annealing simulations, the zeolite frameworks were treated as flexible.Self-diffusion coefficients  were determined from a least-square fitting to mean square displacements computed for the individual MD trajectories using the Einstein relation. 60

DFT calculations
DFT calculations and AIMD simulations used the electronic structure module Quickstep of the CP2K package, version 9.1, which uses the Gaussian and plane wave method. 61,62The rev-vdW-DF2 exchange-correlation functional was used in all calculations, 63 in line with earlier studies addressing the interaction of sizeable organics with zeolites. 37,41Core electrons were represented using Goedeker-Teter-Hutter (GTH) pseudopotentials devised by Krack. 64In the DFT structure optimizations, all atomic positions were relaxed using a BFGS algorithm, whereas cell parameters were held fixed at the experimental values compiled in Table S2.The optimizations employed Gaussian basis sets of double-zeta quality (DZVP-MOLOPT-SR 65 ) and a plane wave cutoff of 600 Ry.A maximal residual force of 5•10 -6 Ha bohr -1 and a maximal displacement of 2•10 -5 bohr between steps were used as convergence criteria.After successful optimizations, total energies were recomputed using larger triple-zeta basis sets (TZVP-MOLOPT) and a plane wave cutoff of 900 Ry.For a given configuration, the adsorption energy was then calculated as: The terms on the right-hand side represent the total energies of the adsorption complex, the guest-free zeolite framework, and the guest molecule in a cubic box with an edge length of 20 Å ( = 1 for TCS and  = 1, 4, 8, 12 for H2O; for the case of water, adsorption energies are either reported for the entire cluster or per H2O molecule, depending on context).Adsorption energies for systems containing co-adsorbed TCS and water were calculated accordingly, subtracting the total energies of both guests.
Single-point calculations using TZVP basis sets were also carried out for fragments of the optimized structures, permitting the calculation of the deformation energies: Here, the first term on the right-hand side corresponds to the total energy of the zeolite framework/TCS molecule extracted from the DFT-optimized adsorption complex, from which the total energy of the DFT-optimized system is subtracted.The deformation energy is always positive.For adsorption of a (H2O)8 cluster, an analogous calculation would correspond to: Here, the term  , contains both attractive H2O-H2O interactions and energy "penalties" arising from deformations of the adsorbed molecules.In the light of the ability of the H2O molecule to act as hydrogen bond donor and acceptor, it can be expected that the former contribution will typically outweigh the latter, and that the overall term will hence be negative.

DFT-based AIMD simulations
Like the DFT calculations described in the preceding subsection, DFT-based AIMD simulations made use of the CP2K package, also using the rev-vdW-DF2 exchange-correlation functional, GTH pseudopotentials, and double-zeta basis sets.The AIMD simulations were carried out in the  ensemble for T = 298 K, using a Nosé-Hoover thermostat 66,67 (time constant 50 fs).
Trajectories were calculated for a total duration of 12.5 ps, using a time step of 0.5 fs (25,000 steps).Unless otherwise noted, the final 10 ps of each AIMD trajectory were taken as production part that was used in the analysis.The AIMD simulations usually started from DFToptimized structures (exception: simulations of TCS and H2O co-adsorption in protonated zeolites, discussed in more detail in the Results section).For each structure considered, three independent 12.5 ps simulations were run.After completing AIMD simulations for a given TCS@Zeo configuration, the zeolite framework alone, and TCS in a box, the internal energy of adsorption was calculated as: Here, the terms in brackets correspond to the averages over the total energies (sum of potential and kinetic energy) obtained from the production stages of the AIMD simulations.The index "3tr" highlights that these averages were computed over the three independent trajectories.To have a measure of the uncertainty in the internal energies of adsorption, ∆ , was recalculated using only subsets of two trajectories in the calculations of the averages in equation ( 5), considering all possible permutations.The resulting estimated boundary values of ∆ , are reported as ∆ and ∆ , respectively.For relevant combinations of atoms/group of atoms, radial distribution functions (RDFs) were calculated using the VMD code, version 1.9.4. 68With one exception, discussed in the Results section, the presented RDFs always correspond to averages over three trajectories.Furthermore, angles between the planes defined by the two TCS phenyl rings  were measured using the TRAVIS code, both for lowest-energy configurations and as averages over the AIMD trajectories. 69,70esults and Discussion

Preliminary investigation: Diffusion of TCS in all-silica zeolites
As a first step, FF-based MD simulations were used to evaluate whether TCS could diffuse through the pores of the six all-silica zeolites chosen for this study.For each zeolite, five trajectories were calculated, starting from the TCS@Zeo configurations that were also used as starting points for the DFT optimizations.,550/573•10 -8 cm 2 s -1 computed for benzene/para-xylene in FAU at 300 K 71,72 ).While diffusion is distinctly slower in IFR, it can be inferred that the TCS molecule is still able to move through the pores.No diffusion on the nanosecond timescale occurs for BEA, where TCS remains in the vicinity of its initial location during all 50 ns simulations.This finding disagrees with the experimental observation of a significant TCS uptake from aqueous solution in zeolite beta samples having different Si/Al ratios. 32Although the difference between the idealized all-silica BEA model used in the calculations and real zeolite beta samples, which exhibit stacking disorder of the tetragonal beta layers along the [001] direction, 48 could be responsible for this discrepancy, other explanations cannot be ruled out.Altogether, the findings for BEA indicate that the MD simulations make a rather conservative estimate of the diffusion behavior, giving confidence in the prediction that TCS can diffuse through the pore systems of the other zeolites.Adsorption energies of TCS@Zeo complexes

All-silica zeolites
For the six all-silica zeolites, Boltzmann-weighted adsorption energies ∆ , as well as adsorption and deformation energies of the lowest-energy TCS@Zeo configurations are compiled in Table 3 (results for all configurations are given in Table S5).With the exception of FAU, the adsorption energies of the other five all-silica zeolites fall in a relatively narrow range from −153 to −160 kJ mol -1 .While the weaker interaction in FAU can be explained with the much larger pore diameter, resulting in reduced dispersion interactions between TCS and the pore walls, there is no apparent correlation of the adsorption energy with pore diameter and/or framework density for the remaining systems.The deformation energies are always small, with the sum of the two ∆ terms amounting to less than 5% of the total adsorption energy in all zeolites except MOR (where it reaches 6.5%).Rather different observations were made in the earlier DFT study of CBZ adsorption, where the sum of the deformation energies exceeded 10% of the total adsorption energies for some zeolites, notably BEA and MOR. 41ereas the CBZ molecule is rather rigid, containing a tricyclic aromatic system, TCS can adjust much more easily to different pore environments because the two phenyl rings can rotate essentially independently.The flexibility manifests in the variation of the angle  , which amounts to 89.0° in the crystal structure of triclosan at 150 K. 42 This angle is reduced to 75.9° in the DFT-optimized structure of isolated TCS, and  values varying from 53.9 to 92.5° occur in the lowest-energy TCS@Zeo complexes.There is no apparent correlation with the magnitude of the ∆ , term (Table 3).
As discussed above, the DFT-optimized lowest-energy TCS conformer possesses an intramolecular hydrogen bond from the hydroxyl group to the central oxygen atom.This bond remains present in all of the lowest-energy TCS@Zeo complexes, and an evaluation of the hydrogen bond distances reveals a certain shortening of the distance (H As the adsorption energies are not correlated with the hydrogen bond distance, it can be concluded that hydrogen bonds make only a minuscule contribution to the total zeolite-TCS interaction.To illustrate typical adsorption complexes with and without hydrogen bonds, the lowest-energy TCS@AFI and TCS@IFR configurations are visualized in Figure 2.  S5).The adsorption energies are distinctly more negative than for the corresponding all-silica zeolites, with the differences amounting to −31 kJ mol -1 (corresponding to a change by −25%) for FAU and to −36 kJ mol -1 (−23%) for MOR.This increase in interaction strength can be straightforwardly attributed to the formation of a hydrogen bond from the framework proton to the TCS Oh atom, visualized in the lowestenergy configurations shown in Due to steric reasons, the hydroxyl group of TCS does not form any additional hydrogen bonds to framework oxygen atoms in the lowest-energy TCS@H-FAU/H-MOR complexes.

AIMD simulations
AIMD simulations for all-silica FAU and MOR as well as H-FAU and H-MOR started from the lowest-energy TCS@Zeo configurations.Internal energies of adsorption were computed according to equation ( 5), using the results from three independent trajectories computed for the TCS@Zeo complex, the zeolite alone, and TCS in a box.The first column of Table 5 gives the internal energies of adsorption as well as the upper and lower boundary values obtained by using averages over two, rather than three, trajectories.On this basis, the uncertainties in the reported values can be estimated to be on the order of ±8 kJ mol -1 for the all-silica zeolites and ±13 kJ mol -1 for the protonated zeolites.As the AIMD simulations used DZVP basis sets, the ∆ , values are not directly comparable to the Boltzmann-weighted adsorption energies ∆ , computed with TZVP basis sets.To allow such a comparison, the ∆ , values were rescaled using scaling factor  , which was calculated as the ratio of the Boltzmann-weighted TZVP and DZVP adsorption energies, i.e.,  = ∆ , ()/ ∆ , ().The scaling factors fall in a fairly narrow range (Table 5), indicating that this is a suitable strategy to extrapolate the DZVP internal energies of adsorption to larger basis sets.Comparing the scaled internal energies of adsorption ∆ , , to the ∆ , values, it is clear that the AIMD-computed values are systematically less negative.The relative difference can be evaluated through the ratio of both quantities, given in the last column of Table 5.The largest deviation of about 20% occurs for all-silica FAU, followed by H-FAU (16%), whereas the differences for all-silica MOR (10%) and H-MOR (7%) are significantly smaller.As observed in the previous study of CBZ adsorption, where FAU, IFR, and MOR were compared, 41 the deviation between "dynamic" and "static" results is directly correlated with the guest molecule's freedom of motion in the zeolite pores: Much more pronounced oscillations are possible in the large cages of FAU than in the comparatively narrow channels of MOR.This results in a larger impact of thermal motion on the internal energy of adsorption in the former system.In the protonated zeolites, the formation of a relatively strong hydrogen bond further reduces the freedom of motion.
The AIMD trajectories were further analyzed to investigate the evolution of the hydrogen bonds during the simulation.Although this analysis remains limited to the picosecond timescale, it can permit some conclusions on the stability/lability of these bonds.For MOR, H-FAU, and H-MOR, average RDFs over the three trajectories were used as basis for the analysis, because differences among the individual trajectories were primarily of a statistical nature.For FAU, however, a distinct behavior was observed in one of the three trajectories (labelled as "Traj1"), which will be discussed in more depth at the end of this section.Therefore, the average RDFs for FAU were calculated only over the other two trajectories.Figure 4a shows the RDFs of the intramolecular hydrogen bond of TCS, including the cumulative RDFs that correspond to the number of interatomic contacts within a given distance range.Both in free TCS and in the adsorption complexes, the intramolecular hydrogen bond persists during the AIMD simulations, with 97.5 to 99% of the distances (H ••• O ) remaining below 2.5 Å.While the maxima in the RDFs of free TCS and of the TCS@FAU and TCS@MOR complexes are very close to the DFT-optimized bond distances tabulated in Table 2, those computed for TCS in protonated zeolites are shifted towards somewhat higher distances.S11), do not exhibit any systematic trends.

Adsorption of water and co-adsorption of TCS and water
The calculations presented up to this point made the simplifying assumption that no coadsorption of water occurs when adsorbing TCS from aqueous solution.While this simplification might appear reasonably justified for all-silica zeolites, which are known to be highly hydrophobic, 74,75 the introduction of framework Al atoms and charge-balancing protons will increase the hydrophilicity.The relative affinity towards water and organic molecules such as phenol has sometimes been evaluated by comparing interaction energies computed for adsorption complexes of individual molecules. 76However, when looking at molecules of vastly different size, adsorption energies obtained for single molecules will (in almost all cases) be much less negative for water, simply because the contribution of dispersion interactions is much larger for a bulky organic molecule like TCS.For this reason, it appears to be more plausible to compare TCS adsorption energies to those obtained for a cluster of water molecules that is similar in size to TCS.An estimation of the molecular volumes based on the crystal structures delivered molecular volumes of 32.5 Å 3 for H2O (hexagonal ice 77 ) and 302 Å 3 for TCS, 42 corresponding to a ratio of about 1:9.3.Moreover, experimental investigations on protonated zeolites have shown that the typical size of water clusters in the vicinity of the acid sites (framework protons) is in the range of eight H2O molecules (for MFI-type ZSM-5 45 ).For these reasons, an adsorption of eight water molecules per simulation cell was assumed in calculations aimed at a comparison of TCS and H2O adsorption energies and in studies of TCS+H2O co-adsorption.Additional results obtained for 1, 4, and 12 H2O molecules per cell were evaluated to investigate the impact of cluster size on the adsorption energy.
For the FAU-and MOR-type all-silica zeolites as well as H-FAU and H-MOR, five configurations containing different arrangements of eight H2O molecules per cell were DFToptimized, and the Boltzmann-weighted adsorption energies were computed in an analogous fashion as for TCS. Figure 6 shows the total adsorption energies and compares them to the TCS adsorption energies (full results for individual configurations are given in Table S6).Even for the two all-silica zeolites, the total adsorption energies for the cluster are much more negative than for TCS, amounting to −389 and −427 kJ mol -1 , respectively.However, a calculation of the guest-guest contribution shows that larger part of the adsorption energy stems from H2O-H2O interactions, which are on the order of −295 kJ mol -1 for both zeolites.
The remaining host-guest contribution is about 20 to 30 kJ mol -1 smaller than the TCS adsorption energies, indicating a preference for TCS over H2O.On a per molecule basis, the H2O adsorption energies correspond to −48.7/−53.3kJ mol -1 for all-silica FAU/MOR.As guestguest interactions are virtually identical, the more negative value for MOR can be attributed to stronger interactions with the framework in the narrower channels.It has to be noted that these values are not directly comparable to thermodynamically measurable quantities that represent the intermolecular interactions in liquid water, specifically the enthalpy of vaporization (44 kJ mol -1 at 298 K).For a semi-quantitative estimation, however, they can be compared to the intermolecular contribution to the total energy obtained in DFT optimizations of a box of water molecules using the same approach.Calculations for four different boxes containing 256 H2O molecules (DFT-optimized density ≈ 1.06 g cm -3 ) delivered energies on the order of −60 kJ mol -1 per H2O molecule.In other words, the DFT calculations correctly predict the interaction to be weaker in the all-silica zeolites than in liquid water, in line with the  This finding agrees with experimental observations 45 and previous DFT results. 44The framework deprotonation precludes a simple calculation of the contribution of H2O-H2O interactions as done for the all-silica zeolites.However, it is apparent from Figure 6 that the increase in interaction strength when moving from an all-silica zeolite to its protonated counterpart is much more pronounced for water than for TCS.The occurrence of framework deprotonation upon water adsorption, and its absence during TCS adsorption, also confirm the notion of a stronger interaction with water.
It is clear that calculations addressing the individual adsorption of either TCS molecules or H2O clusters can only deliver indirect insights into the relative affinity towards the two species.In particular, the magnitude of the adsorption energy of a (H2O)N cluster that competes with TCS will depend on the cluster size, limiting the quantitative interpretation.To obtain insights into the co-adsorption of TCS and water, the simultaneous interaction of both species with FAUand MOR-type zeolites was investigated.For all-silica FAU and MOR, an approach based on DFT optimizations was employed: Starting from DFT-optimized low-energy TCS@Zeo and 8H2O@Zeo configurations, the other species (eight H2O molecules/one TCS molecule) was added using MC simulations.After a simulated annealing of the MC snapshots, five configurations for each case were extracted and DFT-optimized. Figure 7 shows the Boltzmann-averaged adsorption energies for the 8H2O@(TCS@Zeo) and TCS@(8H2O@Zeo) complexes.First of all, it is worth noting that the adsorption energies are very close, regardless of whether H2O molecules were added to a TCS@Zeo complex or vice versa.Second, the total adsorption energies are more negative than the sum of the individual adsorption energies computed for TCS and water, indicating that H2O-TCS interactions during co-adsorption of both species are, overall, attractive.The difference is larger for FAU than for MOR, which can be understood when looking at exemplary low-energy configurations (Figure 8): In the spacious cages of FAU, the water molecules not only form hydrogen bonds to each other and to framework oxygen atoms, but also to the hydroxyl group of TCS, which simultaneously acts as hydrogen bond donor and acceptor.In MOR, the water cluster occupies a different channel section than the TCS molecule, limiting the possibility of guest-guest interactions.For the protonated zeolites, it is most interesting to elucidate what happens to a TCS molecule that is adsorbed at a framework proton in the presence of water at finite temperature.To investigate this, 8H2O@(TCS@Zeo) configurations were generated for H-FAU and H-MOR using MC simulations.Subsequently, three AIMD simulations starting from different MCgenerated configurations were run for each zeolite.In these cases, the first 2.5 ps of the trajectories were not discarded, but included in the analysis, since the most significant structural changes were typically found to occur during the very initial stages of the simulations.
An inspection of the temperature evolution of the individual simulations showed that no significant "overshoot" in temperature occurred during this initial phase.Although a qualitative interpretation could be based on the visualization of individual frames from the AIMD trajectories alone, a more comprehensive picture can be obtained by looking at the evolution of relevant interatomic distances over time.As analogous observations were made for all three trajectories for H-FAU and H-MOR, a representative example trajectory is discussed here for each zeolite, and additional results are shown in the ESI (Figures S7 to S12).
Since the behavior of H-MOR is somewhat simpler, the results for this zeolite are presented first.9d).This hydrogen bond is formed to a framework oxygen atom that does not neighbor the Al atom (Figure 9g).A somewhat different chain of events unfolds in the AIMD simulations for 8H2O@(TCS@H-FAU).Here, two concerted proton transfers occur during the very initial stages of the AIMD simulation (completed within 0.1 to 0.4 ps): The HZeo proton leaves its position at the framework oxygen atom and moves to the Oh atom of TCS, whereas the Hh proton, initially bonded to Oh, is transferred to a water molecule, forming a hydronium ion (Figure 10a  To compute the adsorption energy in a way that is comparable to the calculations for all-silica zeolites, described above, the following approach was used: For both H-FAU and H-MOR, the last 5 ps of each trajectory were analyzed to identify the frame having the lowest potential energy.This frame was then optimized with DZVP basis sets, and the total energy recomputed using TZVP basis sets.For 8H2O@(TCS@H-FAU), a Boltzmann-averaged adsorption energy of −613.5 kJ mol -1 was obtained.This value is 17.5 kJ mol -1 less negative than the sum of the adsorption energies computed for TCS and for eight H2O molecules, which amount to −153.5 and −478.5 kJ mol -1 , respectively.For H-MOR, on the other hand, the ∆ , value computed from three optimized AIMD frames of −686.8 kJ mol -1 is about 2.5 kJ mol -1 more negative than the sum of the individual values (TCS: −190.6 kJ mol -1 , (H2O)8: −493.7 kJ mol -1 ).
Given the limited sampling of only a few frames, these observations should not be overinterpreted in a quantitative way.Qualitatively, however, it seems reasonable to surmise that the contribution of attractive TCS-H2O interactions is offset by a reduced ability of TCS to interact with the framework proton in the presence of water as compared to the water-free system.

Conclusions
Calculations of the adsorption energies for TCS in all-silica zeolites delivered very similar adsorption energies in the range of −153 to −160 kJ mol -1 for five of the six frameworks, with a weaker interaction being observed only for FAU, which possesses large supercages.This is a notably different finding compared to the previous DFT study of CBZ adsorption, where adsorption energies obtained for these five zeolites varied between −122 kJ mol -1 for BEA and −172 kJ mol -1 for IFR. 41This qualitative difference can be attributed to the higher flexibility of TCS, which allows it to adjust in a way that the interaction with the pore walls in differently shaped pores is maximized.At first glance, it might appear that the low affinity of FAU disagrees with the findings of Jiang et al., who observed the highest affinity towards TCS in a FAU-type sample. 32However, it has to be taken into account that this sample also had the highest Si/Al ratio among the zeolites with 12MR pore systems, thus being most hydrophobic.
Moreover, all isotherm points reported in that study correspond to loadings of at least eight TCS molecules per unit cell (the maximal reported loading is close to two TCS in every supercage, i.e., 16 per FAU unit cell).At such high loadings, attractive guest-guest interactions will contribute to the overall affinity.Such interactions can be expected to be stronger in the supercages of FAU, where the adsorbed molecules have space to reorient in a way that maximizes guest-guest interactions, in comparison to frameworks with narrower pores or channels. 41For a more direct comparison between experiment and DFT, it would be required to measure TCS adsorption in the Henry regime.
The stronger interaction of TCS with protonated zeolites can be attributed to the formation of hydrogen bonds between the framework proton and the Oh atom of TCS, which, unlike the weak Hh•••OZeo hydrogen bonds found in all-silica zeolites, are too stable to be broken at room temperature.However, calculations including the co-adsorption of water showed that the stronger interaction would not necessarily translate into higher removal efficiencies: The HZeo•••Oh hydrogen bond is not strong enough to remain stable in the presence of a cluster of water molecules, and framework deprotonation with concurrent formation of a hydronium ion in the pores is energetically favored.Thus, it appears unlikely that a bonding of TCS to framework protons could be realized in real systems, where these more hydrophilic areas of the structure will be occupied by clusters of H2O molecules.In such a material, TCS could still be adsorbed in the more hydrophobic areas of the pore system, meaning that the presence of some framework Al atoms and charge-balancing protons will not necessarily have a completely detrimental effect on the removal efficiency.However, it can be inferred that the removal efficiency will deteriorate with decreasing Si/Al ratio, as indeed observed experimentally. 32tting together the findings outlined above, it appears that the hydrophobicity should be maximized in order to develop adsorbents with maximal TCS removal efficiency.Since smaller pores afford stronger dispersion interactions with the pore walls, the pores of an "ideal" adsorbent should, of course, be large enough to accommodate TCS and allow its diffusion, but not much larger.However, there appears to be little scope to optimize specific features of the pore shape.Although the present work looked exclusively at high-silica zeolites, it seems likely that these general criteria could also be applied to adsorbents from other classes.Furthermore, similar relationships should be valid for organic molecules of similar hydrophobicity and flexibility.In contrast, the shape and size of the pores have a much larger impact on the affinity towards more rigid molecules like CBZ.This difference highlights that guest molecule flexibility should be taken into account when searching for suitable adsorbents.

Conflicts of interest
There are no conflicts of interest to declare.

Acknowledgments
This research was funded by the Deutsche Forschungsgemeinschaft (German Research Foundation, DFG) through a Heisenberg fellowship (project no.455871835).The author gratefully acknowledges the computing time granted by the Resource Allocation Board and provided on the supercomputer Lise and Emmy at NHR@ZIB and NHR@Göttingen as part of the NHR infrastructure.The calculations for this research were conducted with computing resources under the project hbc00062.The author is grateful to Jakob Brauer for comments on different versions of the manuscript.

Figure 1 :
Figure 1: Top: DFT-optimized structure of triclosan (color code: red = O, grey = C, green = Cl, white = H).Labels of atoms that are of particular relevance in the analysis of hydrogen bonds are given, and the intramolecular hydrogen bond is shown as blue line.Bottom: DFTcomputed electrostatic potential plotted on a DFT electron density isosurface.The TCS molecule is rotated with respect to the top panel to better visualize key features of the ESP.Structure figures were prepared using VESTA.46

Figure 3 .
With a bond length of 1.49/1.63Å in H-FAU/H-MOR, these bonds are much shorter than the Hh•••OZeo hydrogen bonds that are formed in all-silica zeolites.As a result, the OZeo−HZeo bond is elongated significantly, explaining the increased ∆ , term.Moreover, the HZeo•••Oh hydrogen bond withdraws some electron density from the Oh atom, leading to a weakening of the Oh−Hh bond and, concurrently, a shortening of the intramolecular Hh•••Oc hydrogen bond.These changes in the molecular structure can serve to explain the relatively large TCS deformation energies.For TCS@H-MOR, the pronounced rotation of the rings ( = 37.4°) probably also contributes to ∆ , .

Figure 2 :
Figure 2: Lowest-energy TCS@AFI (top) and TCS@IFR (bottom) configurations.Atoms of the all-silica zeolite frameworks are shown in pale yellow (Si) and pale red (O).Panels on the left show the full channel in a projection along the channel axis, those on the right show only the vicinity of the adsorbed TCS molecule.

Figure 3 :
Figure 3: Lowest-energy configurations of TCS in protonated zeolites H-FAU and H-MOR.The framework Al atoms are shown as turquoise spheres.

Figure 5 :
Figure 5: Radial distribution functions (RDFs) of HZeo•••Oh hydrogen bonds in protonated zeolites.The top panel shows the cumulative RDFs.
experimentally observed hydrophobicity of the materials.The hydrophobic nature is confirmed when adsorption energies computed for clusters of different size (1, 4, 8, 12 H2O molecules) are compared (FigureS3): The interaction with a single H2O molecule is very weak, with adsorption energies of −25.1/−38.2kJ mol -1 for FAU/MOR.The computed adsorption energies per molecule become increasingly negative with increasing cluster size, with the most prominent change occurring when moving from one to four molecules per cell.Differences between (H2O)8 and (H2O)12 clusters are relatively modest, amounting to 3.5/1.5 kJ mol -1 per molecule.This indicates that the stability of adsorbed water clusters in this size range does not depend strongly on the exact number of molecules involved.

Figure 7 :
Figure 7: Boltzmann-averaged adsorption energies obtained for all-silica FAU and MOR with co-adsorbed TCS and water.Sums of the adsorption energies of individual guest molecules are shown for comparison (top columns).

Figure 8 :
Figure 8: Lowest-energy configurations of co-adsorbed TCS and water in all-silica zeolites FAU and MOR.Oxygen atoms of H2O molecules are shown in blue.
Figure 9a and 9b shows the time evolution of the OZeo−HZeo and HZeo•••Oh distances during the first 3 ps of the simulation.For a time of about 1.2 ps, the framework proton remains bonded to the framework, and hydrogen-bonded to the TCS molecule.This situation is depicted in Figure 9e, which shows a snapshot from the initial phase of the simulation.After this period, a breaking of the HZeo•••Oh hydrogen bond occurs, and a transient hydrogen bond to a water molecule is formed (Figure 9f).After about 1.6 ps, the covalent OZeo−HZeo bond is broken, with the distance (O − H ) quickly increasing to values above 2 Å.A visualization of a simulation frame obtained after this period (Figure 9g) shows that the proton is now bonded to the H2O molecule that participated in the transient hydrogen bond.It remains bonded to the same H2O oxygen atom for the remainder of the simulation, as is visible from the time evolution of the corresponding HZeo•••OH2O distance (while Figure 9c shows only the first 3 ps, the full 12.5 ps are visualized in Figure S10).After the framework deprotonation, further proton transfers between water molecules occur, as evidenced by the fact that the hydronium ion that is present after 3 ps is not the same molecule that deprotonated the framework (Figure 9g).The TCS molecule moves away from the framework Al atom after the HZeo•••Oh hydrogen bond is broken, now acting as a hydrogen bond donor through formation of a new Hh•••OZeo hydrogen bond (Figure

Figure 9 :
Figure 9: a) to d) Evolution of selected interatomic distances during the first 3 ps of one AIMD trajectory (Traj1) computed for 8H2O@(TCS@H-MOR), e) to g) environment of the framework Al atom at different stages of the simulation.
to 10d).Although both transfers happen almost at the same time, the time evolution of the distances indicates that the framework deprotonation precedes the TCS  H2O proton transfer.After moving towards the TCS molecule, the HZeo proton may form a hydrogen bond to the OZeo atom to which it was originally attached or to a water molecule in the vicinity.While Figure 10e shows a simulation frame at the very beginning of the simulation (after 0.05 ps), prior to deprotonation of the framework, Figure 10e (after 1.0 ps) and 10f (after 2.5 ps) exemplify the consecutive formation of HZeo•••OH2O and HZeo•••OZeo hydrogen bonds after completion the proton transfer.Similar alternations can be observed in the other two trajectories.Furthermore, the Hh proton, now part of a hydronium ion, forms a hydrogen bond to the Oh atom of TCS.

Figure 10 :
Figure 10: a) to d) Evolution of selected interatomic distances during the first 3 ps of one AIMD trajectory (Traj1) computed for 8H2O@(TCS@H-FAU), e) to g) environment of the framework Al atom at different stages of the simulation.Not all eight H2O molecules are visible in the portion of the cage shown in the figure.

Table 1 :
31olites considered in this study.Cell content refers to the supercell used in the calculations.correspondsto the framework density.The last column reports the diameters of the largest included sphere  and the largest diffusing sphere  , taken from the IZA database31.

Table 2
gives the average self-diffusion coefficients  , (averaged over five trajectories) as well as the minimal and maximal individual values of  .As it should be expected, the trend in the self-diffusion coefficients follows the evolution of diameters of the largest diffusing sphere.A fairly large scatter in individual values, which span almost one order of magnitude in the most extreme case, is evident.It is, however, not the purpose of these calculations to deliver accurate predictions of  , but rather to determine whether a diffusion of TCS can be expected at all.This is clearly the case for FAU, AFI, MOR, and CFI, with self-diffusion coefficients on the order of 20 to 2,000•10 -8 cm 2 s -1 .Values in this range are typical for sizeable organic molecules in zeolites with 12 MR windows (e.g., 9

Table 2 :
Self-diffusion coefficients averaged over five trajectories  , ••• O ) with respect to free TCS in all zeolites.This shortening is most pronounced in MOR and IFR, the two zeolites where the ∆ , term is largest.The elongation of the covalent Oh−Hh bond remains negligible.Additional hydrogen bonds to framework oxygen atoms can be found in all zeolites except IFR, with distances (H ••• O ) varying from 2.11 Å in AFI to 2.48 Å in MOR.

Table 3 :
Results for TCS in all-silica zeolites.Top: Boltzmann-averaged adsorption energies ∆ ,, adsorption energies and deformation energies of lowest-energy configurations.Bottom: Selected interatomic distances and angles of lowest-energy configurations.

Table 4
compiles the results for protonated zeolites H-FAU and H-MOR in a largely analogous fashion to Table3(full results are given in Table

Table 4 :
Results for TCS in protonated zeolites.

Table 5 :
Results of AIMD simulations: Unscaled and scaled internal energies of adsorption, including boundary values obtained by considering subsets of trajectories.Boltzmannweighted adsorption energies ∆ , are given for comparison.