Molecular dynamics simulation of secondary sorption behavior of montmorillonite modified by single chain quaternary ammonium cations

Organoclays synthesized from single chain quaternary ammonium cations (QAC) ((CH(3))(3)NR(+)) exhibit different mechanisms for the sorption of nonpolar organic compounds as the length of the carbon chain is increased. The interaction between a nonpolar sorbate and an organoclay intercalated with small QACs has been demonstrated to be surface adsorption, while partitioning is the dominant mechanism in clays intercalated with long chain surfactants. This study presents the results of a molecular dynamics (MD) simulation performed to examine the sorption mechanisms of benzene in the interlayer of three organoclays with chain lengths ranging from 1 to 16 carbons: tetramethylammonium (TMA) clay; decyltrimethylammonium (DTMA) clay; and hexadecyltrimethylammonium (HDTMA) clay. The basis of the overall simulation was a combined force field of ClayFF and CVFF. In the simulations, organic cations were intercalated and benzene molecules were introduced to the interlayer, followed by whole system NPT and NVT time integration. Trajectories of all the species were recorded after the system reached equilibrium and subsequently analyzed. Simulation results confirmed that the arrangement of the surfactants controlled the sorption mechanism of organoclays. Benzene molecules were observed to interact directly with the clay surface in the presence of TMA cations, but tended to interact with the aliphatic chain of the HDTMA cation in the interlayer. The simulation provided insight into the nature of the adsorption/partitioning mechanisms in organoclays, and explained experimental observations of decreased versus increased uptake capacities as a function of increasing total organic carbon (TOC) for TMA clay and HDTMA clay, respectively. The transition of sorption mechanisms was also quantified with simulation of DTMA clay, with a chain length between that of TMA and HDTMA. Furthermore, this study suggested that at the molecular level, the controlling factor for the ultimate sorption capacity is available surface sites in the case of TMA clay, and density of aliphatic chains within the interlayer space for HDTMA clay.     

Prediction of the temperature dependency of Henry's law constant from chemical structure

A new model to estimate the temperature dependency of Henry's law constant in water for organic compounds from the two-dimensional structure is presented. Air/water partition enthalpies of 456 chemicals were fitted to 46 substructural parameters with a squared correlation coefficient r2 of 0.81 and a standard error of 7.1 kJ/mol. The compound set covers various organic compound classes with the atom types C, H, N, O, F, Cl, Br, I, and S. Application of the model together with experimental data for 25 degrees C to a set of 462 compounds with 2119 experimental Henry's law constants at temperatures below 20 degrees C yields a predictive squared correlation coefficient q2 of 0.99 and a standard error of 0.21 logarithmic units. The prediction capability is further evaluated using cross validation and permutation.     

Enhanced sorption of polycyclic aromatic hydrocarbons to tetra-alkyl ammonium modified smectites via cation-pi interactions

The objective of this study was to characterize molecular sorptive interactions of polycyclic aromatic hydrocarbons (PAHs) by organoclays modified with quaternary ammonium cations. Three PAHs, naphthalene (NAPH), phenanthrene (PHEN), and pyrene (PYR), and three chlorobenzenes, 1,2-dichlorobenzene (DCB), 1,2,4,5-tetrachlorobenzene (TeCB), and pentachlorobenzene (PtCB), were sorbed from aqueous solution to reference montmorillonite clays (SWy-2) exchanged respectively with tetramethyl ammonium (TMA), tetraethyl ammonium (TEA), tetra-n-butyl ammonium (TBA), and hexadecyltrimethyl ammonium (HDTMA) cations. Solute hydrophobicities are compared between PAHs and chlorobenzenes using the solute n-octanol-water partition coefficient, n-hexadecane-water partition coefficient, and polyethylene-water distribution coefficient. The PAHs show several- to more than 10-fold greater sorption than the chlorobenzenes having close hydrophobicities but fewer delocalized pi electrons (NAPH/DCB, PHEN/TeCB, and PYR/ PtCB) by TEA-, TBA-, and HDTMA-clays. Furthermore, the PAHs show greater trends of solubility enhancement than the compared chlorobenzenes by TMA, TEA, and TBA in aqueous solution. The enhanced sorption and aqueous solubility of PAHs are best described by cation-pi interactions between ammonium cations and PAHs relative to chlorobenzenes that are incapable of such interactions. Cation-pi complexation between PAHs and tetra-alkyl ammonium cations in chloroform was verified by ring-current-induced upfield chemical shifts of the alkyl groups of cations in the 1H NMR spectrum.     

Importance of adsorption (hole-filling) mechanism for hydrophobic organic contaminants on an aquifer kerogen isolate

Sorption and desorption behaviors of four hydrophobic organic compounds (HOCs) were investigated for an isolated kerogen material from Borden aquifer material with total organic carbon of 0.021%. The solubility-normalized modified Freundlich equation and the combined linear and Polanyi-Dubinin (PD) equation can quite well describe the sorption or desorption isotherms. The partition component is estimated and compared using desorption data, dual-mode modeling, and the reported partition coefficients. The result suggests that the dual-mode modeling and the combined linear and PD modeling may overestimate the partitioning component. The partition component is not so important as assumed before in sorption of HOCs for the studied sorbent. As the fitted PD equation has an exponent parameter b' approaching 1, it is equivalent to the modified Freundlich equation. The small molecules 1,2-dichlorobenzene (DCB) and naphthalene (Naph) have higher adsorption volumes. The lower adsorption volumes for 1,3,5-trichlorobenzene (TCB) and phenanthrene (Phen) suggest that accessibility to the holes of kerogen by large HOC molecules is reduced. The desorption hysteresis is approximately constant for DCB when the relative aqueous concentration ranges from 0.0007 to 0.6, but for Phen is only obvious at higher relative aqueous concentrations. The varied sorption and desorption behaviors for DCB and Phen are satisfactorily explained by an adsorption/ hole filling mechanism and entrapment of some adsorbates in the kerogen matrix and by possible pore deformation mechanism at high concentrations.     

Detailed sorption isotherms of polar and apolar compounds in a high-organic soil

Sorption isotherms of 13 apolar liquids and solids and polar solids-six in unprecedented detail-are used to evaluate a polymer-based model for natural organic matter. While all isotherms are nonlinear, the "running" Freundlich exponent n varies markedly with concentration. The isotherms show linear-scale inflection consistent with the presence of flexible (deformable) porosity as predicted by the glassy polymer-based Extended Dual-Mode Model (EDMM). The EDMM assumes dissolution and hole-filling domains in the organic solid, with provision for sorbate-caused plasticization of the solid and "melting" of the holes. Features of the EDMM are illustrated for chlorinated benzenes in poly(vinyl chloride). The solutes fall into categories of "hard" (aliphatics and 2,4-dichlorophenol) and "soft" (chlorinated benzenes, 2-chloronitrobenzene) according to their ability to plasticize organic matter. Comparison of domain coefficients at infinite dilution reveals that organic solutes have a modestly greater affinity for holes than dissolution sites (by 0.1-0.6 log unit), as expected by the polymer model. Sorption of CHCl3 shows time-dependent hysteresis diminished at high concentrations by the plasticizing effect. Sorption of CHCl3 also shows a type of hysteresis for glassy solids known as the "conditioning effect" in which high loading of sorbate increases hole population upon its removal and thus leads to enhanced uptake and nonlinearity when a second sorption is performed. A Polanyi-based, fixed-pore filling model applied to the adsorption component of the isotherms gave widely variant volumetric pore capacity, contrary to its own stipulations, and could not explain the hysteresis.     

Organic sorbate-organoclay interactions in aqueous and hydrophobic environments: sorbate-water competition

Sorption of nitrobenzene, phenol, and m-nitrophenol from water and n-hexadecane was measured on Na-montmorillonite and organoclays in which 41 and 90% of the exchange capacity of the Na-clay was occupied by hexadecyltrimethylammonium. The strength of sorbate-sorbent interactions in n-hexadecane for all three sorbents was in the following order: nitrobenzene < phenol < m-nitrophenol. The magnitude of the distribution coefficients suggests that the contribution to solute uptake of partitioning between n-hexadecane and the organic pseudophase of the dried organoclays is minor, whereas the major contribution is from adsorptive sorbate-sorbent interactions. Sorption isotherms obtained in different solvents were compared using a sorbate activity scale. In the organoclays, the stronger the tendency of a sorbate to interact with sorption sites, the less pronounced is the reduction in the activity-based sorption due to competition with water. The order of this reduction for the different sorbates is nitrobenzene > phenol > m-nitrophenol. The weakening of sorbate-sorbent interactions resulting from water-sorbate competition might be mitigated by interaction between the organic sorbate and sorbed water molecules. Since the more strongly interacting organic compounds are less susceptible to suppression of sorption in the presence of water, hydrating organoclays may result in an increased differentiation between "weakly" and "strongly" interacting ("nonpolar" and "polar") compounds in the organoclay phase.     

Sorption of Ethyl Butyrate and Octanal Constituents of Orange Essence by Polymeric Adsorbents

Sorption characteristics, rates and capacities of polymeric adsorbents (XAD-4, XAD-7, XAD-16, Duolite ES-865, Duolite S-761, Porapak-Q and XUS-43436) for ethyl butyrate and octanal were determined using model solutions. Sorption of these two principal components of aqueous orange essence, was evaluated utilizing the Freundlich isotherm model. Breakthrough curves were determined using XAD-16 as adsorbent. The capacity of XAD-16 for ethyl butyrate in the column system was 426±212 mg/g. There was no column breakthrough for octanal even after 130 bed volumes. Most (91.4%) of the adsorbed ethyl butyrate was eluted from XAD-16 resin by 95% ethanol elution. Recovery of adsorbed octanal from the resin was 66.4%.     

Sorption of ionic surfactants to estuarine sediment and their influence on the sequestration of phenanthrene

The sorption of an anionic surfactant (sodium dodecyl sulfate; SDS) and a cationic surfactant (hexadecyl trimethylammonium bromide; HDTMA) to estuarine sediment has been studied in river water and seawater. Sorption isotherms for SDS were essentially linear in both waters, suggesting a nonspecific, hydrophobic interaction between the SDS tail and particle surface. Sorption of HDTMA was considerably greater, more nonlinear, and more sensitive to water composition. These observations were attributed to a combination of both electrostatic and hydrophobic interactions between the surfactant and particle surface, the formation of admicelles, and salinity-induced structural alteration of the hydrophobic tail of the HDTMA molecule. Presence of SDS caused a reduction in the sorption of phenanthrene to estuarine sediment because of the competitive effects of the surfactant tail for hydrophobic sorption sites on the particle surface. Conversely, the presence of HDTMA caused significant enhancement in phenanthrene sequestration because of head-on sorption of surfactant molecules and a resulting, more hydrophobic particle surface. The most persistent feature of our results was an inverse dependence of unit sorption on particle concentration, and an empirical algorithm defining the effect was used to calculate the sediment-water fractionation of realistic concentrations of reactants in the estuarine water column. The results of these calculations, and the more general findings of this study, significantly improve our understanding of both the transport and fate of ionic surfactants in the estuarine environment, and the effects that these surfactants have on the partitioning of hydrophobic organic micropollutants.     

Interaction of the macrolide antimicrobial clarithromycin with dissolved humic acid

The association of the cationic macrolide antimicrobial clarithromycin (CLA) with dissolved Elliot soil humic acid (ESHA) was investigated as a function of solution chemistry. CLA-ESHA association was strongly pH-dependent, reaching a maximum near pH 6.5, and was modeled successfully using FITEQL as a 1:1 complexation reaction between CLA+ and discreet deprotonated acidic functional groups with pKa values of 4 and 6. Approximate order of magnitude increases in ionic strength produced approximately 10-fold decreases in CLA+ -ESHA association. Coefficients for CLA+ -ESHA association were significantly smaller in the presence of K+ vs Na+. Sorption data were well-fit by the Freundlich model; the Freundlich exponent was <1, suggesting CLA+ interacted with sites having a range of binding energies. Sorption appeared largely reversible; little sorption-desorption hysteresis was observed. The affinities of erythromycin and CLA+ for ESHA association sites were indistinguishable, suggesting interaction with specific sorption sites. Comparison of experimentally determined CLA-ESHA association coefficients with those predicted from single-parameter linear free energy relationships based on log Kow suggested limited contribution of hydrophobic interactions to CLA-ESHA association at environmentally relevant pH values. CLA-ESHA association constants were similar in magnitude (10(3.9)-10(4.6) to those of many nonpolar organic contaminants, and macroscopic binding data were consistent with cation exchange dominating CLA+ -ESHA association.     

Sorption behavior of nonylphenol in terrestrial soils

Nonylphenol (NP) as an intermediate from anaerobic degradation of widely used nonionic surfactants occurs widespread in the environment. Partition behavior of this toxic and endocrine-disrupting chemical between soil and water was not examined until yet. The objective of this investigation was to quantify sorption and desorption behavior of 4-nonyl[14C]phenol in a set of 51 soils using the batch equilibrium approach. Kinetic studies indicated apparent equilibrium within 20 h. Sorption was influenced by sorbate structure as could be shown with branched 4-nonyl[14C]phenol and the linear 4-n-NP, respectively. Linear 4-n-NP behaves differently from the branched isomers of 4-NP. Sorption of 4-nonyl[14C]phenol tested with five different initial concentrations resulted in linearly fitted isotherms that provided calculation of sorption partition coefficients (KP). Desorption partition coefficients (KP-des) revealed hysteresis independent of soil properties but decreasing with decreasing initial NP concentrations. KP values were correlated with organic carbon content of the soils yielding a log KOC of 3.97.     

Sorption of very hydrophobic organic compounds (VHOCs) on dissolved humic organic matter (DOM). 2. Measurement of sorption and application of a Flory-Huggins concept to interpret the data

Sorption phenomena of very hydrophobic compounds (VHOCs, log K(OW) > 5) on dissolved humic organic matter (DOM) are overwhelmingly based on partitioning processes. In this respect, DOM is very similar to "rubbery" soil/sediment OM. To exclude system adsorption effects, the DOM sorption coefficients (K(DOM)) of VHOCs were determined using a dynamic approach based on the VHOCs' aqueous solubility enhancement in the presence of DOM. Partition coefficients are strongly correlated to the analytes' Kow across the alkane, PAH, and PCB groups under study. These three "families" are regarded to be good models of hydrophobic partitioning. On the basis of a uniform one-parameter concept characterizing sorption on amorphous polymers, Hildebrand solubility parameters of amorphous polymeric sorbents, including DOM, and of sorbates can be calculated on the basis of partition coefficients. Likewise, partition coefficients can be estimated using Hildebrand solubility parameters. Literature-based partition coefficients on DOM fit very well in this universal one-parameter concept. On using our own sorption data of PAHs, PCBs, and alkanes on DOM, an almost identical solubility parameter for the DOM polymer under study is obtained. The concept is also very useful in understanding both waterborne and airborne bioconcentration processes, which are considered to be partitioning phenomena.     

Influence of rhamnolipids and triton X-100 on the desorption of pesticides from soils

A rhamnolipid biosurfactant mixture produced by P. aeruginosa UG2 and the surfactant Triton X-100 were tested for their effectiveness of enhancing the desorption of trifluralin, atrazine, and coumaphos from soils. Sorption of both surfactants by the soils was significant and adequately described by the Langmuir-type isotherm. Values of maximum sorption capacity (Qmax) and Langmuir constant (Klang) did not correlate with the amount of soil organic matter. Our results indicate that clay surfaces play an important role in the sorption of surfactants. When surfactant dosages were high enough to reach soil saturation and maintain an aqueous micellar phase, pesticide desorption was only enhanced. At dosages below soil saturation, surfactants sorbed onto soil, increasing its hydrophobicity and enhancing the sorption of the pesticides by a factor of 2. Similar values of water-soil partition coefficients (Ksol*) for aged and fresh added pesticides to soils indicate that the aging process used did not significantly after the capability of either surfactant to desorb the pesticides. A model able to estimate equilibrium distributions of organic compounds in soil-aqueous-micellar systems was tested against experimental results. The determined organic carbon partition coefficients, Koc values, indicate that, on a carbon normalized basis, sorbed Rh-mix is a much better sorbent of pesticides than TX-100 or soil organic matter. These results have significant implications on determining the effectiveness of surfactants to aid soil remediation technologies.     

How do the partitioning properties of polyhalogenated POPs change when chlorine is replaced with bromine?

Information about the mobility and the partitioning properties of brominated persistent organic pollutants, the environmental levels of which are sometimes higher than those of the chlorinated analogues, is limited. We estimated n-octanol/ water (log K(OW)), n-octanol/air (log K(OA)), and air/water (log K(AW)) partition coefficients for 1436 chloro- and bromo-analogues of dibenzo-p-dioxins, dibenzofurans, biphenyls, naphthalenes, diphenyl ethers, and benzenes by employing quantitative structure-property relationship (QSPR) techniques. The searches for similar partitioning patterns were performed by means of two-dimensional cluster analysis. Five classes of compounds were identified. Each of the class is characterized by similar partition coefficients and, in consequence, similar environmental properties. Finally the data was fitted into a simple multimedia model involving the partitioning map. In addition, we found thatthe changes in the partition coefficients upon the replacement of chlorine with bromine were constant: 0.11, 0.31, and -0.21 per bromine atom for log K(OW), log K(OA), and log K(AW), respectively. On the basis of this observation, a method for rapid estimation of changes in the partition coefficient upon chlorine-bromine substitution was proposed.     

Sorption of a diverse set of organic vapors to urban aerosols

Sorption to urban aerosols is a key process in determining the transport and fate of organic pollutants in the atmosphere. The sorption properties of two urban aerosol samples have been determined using aerosol/air partition coefficients measured for a large set of diverse organic vapors. The dominant sorption process could be identified for both samples with two complementary methods: (a) by applying poly-parameter linear free energy relationships (LFERs) to the data sets, and (b) by evaluating the specific surface area, the elemental carbon (EC) content, and the organic matter (OM) content of the aerosols in combination with various sorbent-air partition coefficients from the literature. This revealed that sorption to the two urban aerosols was dominated by absorption into OM and that the diverse data set could be evaluated with an absorption model. The data further revealed that neither EC nor OM was fully available for sorption. The latter leads to the hypothesis that aerosol OM in urban samples has characteristics comparable to those of glassy polymers.     

Black carbon-inclusive modeling approaches for estimating the aquatic fate of dibenzo-p-dioxins and dibenzofurans

A novel black carbon (BC) inclusive modeling tool is applied to estimate the distribution and long-term fate of dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) in the Norwegian Grenland Fjords. Three versions of the model were developed in which sediment-water partitioning was described using (i) an amorphous organic carbon (AOC) partitioning sorption model without BC sorption, (ii) a combined AOC and BC sorption model based on the Freundlich isotherm, and (iii) a combined BC-AOC model based on the Langmuir isotherm. The predictive ability of the three different models was evaluated for 17 PCDD/Fs by comparison of model predictions with observed organic carbon normalized sediment-water partition coefficients (K(TOC)) and with measured concentrations. All three versions of the model were able to predict concentrations that were in reasonable agreement with measured particulate concentrations (i.e., within a factor of 4 of median values). Estimated particulate concentrations were less sensitive to the model choice because the majority of the mass of these hydrophobic chemicals is associated with particulates regardless. However, for estimation of K(TOC) or dissolved water concentrations, both versions of the combined AOC and BC sorption models provided greatly improved estimates compared to the AOC-only model.     

Effect of fluorotelomer alcohol chain length on aqueous solubility and sorption by soils

Fluorotelomer alcohols (FTOHs) are a group of polyfluorinated alkyl chemicals that have been widely studied as precursorsto perfluorocarboxylates such as perfluorooctanoic acid and for which knowledge on their fate in soils is sparse. The solubility and sorption by soil of the homologous 4:2 to 10:2 FTOHs were measured in water or cosolvent/ water solutions. For the smaller 4:2 and 6:2 FTOHs, solubility and sorption could be measured adequately in aqueous systems although transformation was apparent even in gamma-irradiated and autoclaved systems. Sorption coefficients estimated by measuring both sorbed and solution-phase concentrations were not significantly affected by the biotransformation process. The use of cosolvents was employed for probing the behavior of the longer-chain FTOHs with limited aqueous solubility. A single log-linear correlation between aqueous solubility and modified McGowan molar volumes resulted for the n-alkanols and FTOHs. Soil organic carbon (OC) consistently appeared to be the key soil property influencing sorption of the FTOHs while the perfluorocarbon chain length was the dominant structural feature influencing solubility and sorption. Each CF2 moiety decreased the aqueous solubility by -0.78 log units (compared to 0.60 log units for each CH2 addition in hydrogenated primary alcohols), and increased OC-normalized sorption coefficients (Koc) by -0.87 log units. Good log-log linear correlations between Koc and both octanol-water partition coefficients and solubility were observed for the FTOHs.     

Sorption and displacement of pyrene in soils and sediments

Sorption isotherms of pyrene on soils and sediments were examined to understand its sorption behavior. All systems examined exhibited nonlinear sorption. Sorption nonlinearity was found to be a function of the polarity index of soil/sediment organic matter (SOM), suggesting that the degree of condensation of SOM, characterized by its polarity index, was correlated with the sorption behavior of pyrene. The polarity index of SOM could be a new factor for explaining the sorption nonlinearity. The sorption affinity of two soils and two sediments for pyrene increased with decreasing SOM polarity. A higher sorption affinity in the two soils was associated with a higher degree of condensation of SOM compared to that of the two sediments. A displacement test was performed after pyrene sorption using phenanthrene as a displacer. Pyrene was displaced in all systems examined, and nonlinearity became less pronounced after displacement. Such an increase in isotherm linearity implied that sorption site energies became more homogeneous after displacement. Furthermore, the site energy distribution IE*) derived from the Freundlich model parameters showed that energy reduction of high-energy sites was more significant than that of low-energy sites after displacement. In addition, a decrease in sorption capacity after displacement could be ascribed to the partial depletion of sorption sites by the displacer. The displacement data indicated that the cocontaminant can have potential effects on the fate and bioavailability of anthropogenic organic pollutants sorbed in soils and sediments, thus affecting their exposure risks.     

Solubility-normalized combined adsorption-partitioning sorption isotherms for organic pollutants

Equilibrium sorption isotherms were measured for five different low-polarity organic compounds (benzene, trichloroethene, 1,2- and 1,4-dichlorobenzene, and phenanthrene) over a wide concentration range. The investigated sorbents can be grouped into the following three classes: (1) humic soil organic matter, which shows linear sorption isotherms (solely partitioning, as observed in the peat sample); (2) carbon materials, which were thermally altered (due to their natural history or industrial production) and thus contain a high specific surface area and exhibit nonlinear isotherms, and (3) pure engineered microporous materials (e.g., zeolites and activated carbon), where adsorption is solely due to a pore-filling process. Sorption of all compounds was fitted very well by the Polanyi-Dubinin-Manes (PDM) model, which for sorbents containing humic organic matter (e.g., peat) was combined with linear partitioning. Both the partitioning and the Polanyi-Dubinin-Manes model predict unique sorption isotherms of similar compounds if the solubility-normalized aqueous concentration is used. In addition, an inverse linear relationship between the distribution coefficient (Kd) and water solubility, which was very well confirmed by the data, is obtained. This also leads to unit-equivalent Freundlich sorption isotherms and explains the often observed apparent correlation between sorption capacity at a given concentration (e.g., Freundlich coefficient) and sorption nonlinearity (Freundlich exponent).     

山楂汁中有机酸和总黄酮吸附动力学研究

为了探讨山楂汁中有机酸和总黄酮的同时吸附的动力学规律,实验设计含有不同浓度有机酸和总黄酮的山楂果汁,与一定量的树脂(D941树脂或NKA树脂)混合,搅拌吸附,用Langmuir等温模型和Freundlich等温模型分析吸附平衡数据。结果表明,D941树脂和NKA树脂对山楂汁中的有机酸和总黄酮的吸附平衡数据符合Langmuir等温模型。根据Langmuir等温模型,求出这两种树脂分别对有机酸和总黄酮的最大吸附,与大孔吸附树脂NKA相比,发现大孔弱碱性阴离子交换树脂D941吸附较多的有机酸和较少的总黄酮,可以作为山楂果汁降低有机酸而保留总黄酮研究中较为理想的树脂。D941树脂和NKA树脂对山楂汁中的有机酸的吸附能力低于对纯柠檬酸溶液的吸附能力,证明总黄酮的吸附会阻碍有机酸的吸附。