Correlating the physical-chemical properties of phthalate esters using the 'three solubility' approach

A quantitative structure-property relationship (QSPR) method for the correlation of physical-chemical properties and partition coefficients, namely the 'three solubility' approach, is described and applied to a group of 22 phthalate esters. The solubilities or 'apparent-solubilities' of these substances in the liquid state are compiled and correlated against Le Bas molar volume in the three primary media of air, water and octanol. From these solubilities the air-water (K(AW)), octanol-water (K(OW)) and octanol-air (K(OA)) partition coefficients are deduced. Estimated solubilities in water and octanol-water partition coefficients are shown to compare favourably with more recent accurate measurements. A set of selected values is presented, with error limits, which is recommended for use in modelling and assessment studies. Some environmental implications are discussed of the large range in property values for this series.     

Estimating K(oc) for persistent organic pollutants: limitations of correlations with K(ow)

The n-octanol/water partition coefficient (K(ow)) is commonly used to predict the soil or aquatic particle water partition coefficient normalized to organic carbon (K(oc)). Many correlations are available covering several chemical classes and ranges of hydrophobicity. This work indicates the K(ow) may not be a strong predictor for persistent organic pollutants (POPs) which are defined here as chemicals with logK(ow) > 5.0. In addition, the correlation developed in this work for POPs will still result in a predicted value which is of by a factor of 15. Accordingly, care must be taken when applying K(oc) estimations using K(ow) for POPs until more suitable correlations are developed.     

Selective removal of organochlorine pesticides (OCPs) from aqueous solution by triolein-embedded composite adsorbent

A novel composite adsorbent (CA-T) was used for the selective removal of organochlorine pesticides (OCPs) from aqueous solution. The adsorbent was composed of the supporting activated carbon and the surrounding triolein-embedded cellulose acetate membrane. Scanning electron microscopy (SEM), N2 adsorption isotherms and fluorescence methods were used to characterize the physicochemical properties of CA-T. Triolein was perfectly embedded in the cellulose acetate membrane and deposited on the surface of activated carbon. The adsorbent was stable in water and no triolein leakage was detected during the test periods. Some organochlorine pesticides (OCPs), such as dieldrin, endrin, aldrin, and heptachlor epoxide, were used as model contaminants and removed by CA-T in laboratory batch experiments. The adsorption isotherm followed the Freundlich equation and the kinetic data fitted well to the pseudo-second-order reaction model. Results also indicated that CA-T appeared to be a promising adsorbent with good selectivity and satisfactory removal rate for lipophilic OCPs from aqueous solutions when present in trace amounts. The adsorption rate and removal efficiency for lipophilic OCPs were positively related to their octanol-water partition coefficients (log K(ow)). Lower residual concentrations of OCPs were achieved when compared to granular activated carbon (GAC).     

Equilibrium partitioning of PCB congeners in the water column: Field measurements from the Hudson River

The environmental behavior of hydrophobic organic compounds in water is driven by partitioning between dissolved and sorbed phases. Partitioning behavior of a compound is often based on empirical relationships to other properties of the chemical, such as water solubility and octanol-water partition coefficients, but actual partitioning in the environment may differ significantly from such predictions. We conducted intensive studies of the distribution of PCBs in the fresh water portion of the Hudson River, using sensitive capillary-column gas chromatography methods to calibrate and resolve quantitations for 90 PCB congeners in 48 samples at 10 locations. A linear equilibrium model of PCB congener partitioning, when corrected for temperature and suspended-matter organic carbon content, provides a good representation of phase distribution. When particulate-phase concentrations are predicted from dissolved concentrations with a two-phase model, the predictions are unbiased for the majority of samples and the average percent difference between observed and predicted particulate concentrations is ±43 percent. Estimated in situ partition coefficients show systematic differences from partition coefficients predicted from octanol-water partitioning. Partitioning to colloids appears to be a significant component of total concentration for mono- and dichlorobiphenyls, but not for more highly chlorinated congeners.The colloidal fraction may still cause significant overestimation of the bioavailable fraction for more hydrophobic congeners when a two-component model is used.     

Physical origin for the nonlinear sorption of very hydrophobic organic chemicals in a membrane-like polymer film

Bioconcentration factor (BCF) is often assumed to be linearly associated with the octanol-water partition coefficient K(ow) for hydrophobic organic chemicals (HOCs). However, a large amount of data has suggested that the correlation between the logBCF and logK(ow) is curvilinear for HOCs. Similar curvilinear relationship has also been noticed for sorption of HOCs into poly(dimethyl)siloxane (PDMS), a polymer with cross-linked interior structures. So far no satisfactory explanation has been given to account for the deviation. In this study, we acquired additional experimental data to show that the curvilinear relationship between the log-based PDMS-coated fiber-water partition coefficient (logK(f)) and logK(ow) for polychlorinated biphenyls (PCBs) was indeed a reflection of the sorption process occurring in PDMS film other than experimental defects. The physical origin of the nonlinearity was pinpointed based on the theory of phase partitioning for HOCs. The linear relationship is observed if the solute molecule is considerably smaller than the size of a monomer unit of PDMS in that the Gibbs free energy required for cavity formation in PDMS is comparable to that in octanol. Higher free energy of cavity formation is needed to create sufficient free volume if the PCB molecular size is comparable to or larger than the monomer unit of PDMS. On the other hand, the free energy of cavity formation in octanol remains almost constant when this occurs, resulting in the observed curvilinear relationship. The proposed model adequately explains the observed data, as well as sheds lights into the physical origin of the steric interactions of large molecular size solute with the PDMS polymer network.     

Evaluating coal tar-water partitioning coefficient estimation methods and solute-solvent molecular interactions in tar phase

Equilibrium partitioning coefficients between an industrial coal tar sample and water (K(CT/w)) were determined for 41 polar and nonpolar solutes in batch systems. Together with literature values, 69 K(CT/w) data were analyzed using the following model approaches: Raoult's law, the single parameter linear free energy relationship (SPLFER) with octanol-water partitioning coefficients (K(ow)), the linear solvation energy relationships (LSERs), SPARC and COSMOtherm. Estimations by Raoult's law and the SPLFER agreed well with the experimental log K(CT/w) values for the investigated coal tar, with root mean square errors (RMSE) of 0.31 and 0.33, respectively. LSER resulted in as good estimations (RMSE=0.29) as the previous two. The LSER analysis revealed significant hydrogen (H)-bond acceptor properties of the studied coal tar phase. Using naphthalene as a surrogate solvent for the coal tar phase, SPARC and COSMOtherm provided fairly good predictions (RMSE of 0.63 and 0.65, respectively) of log K(CT/w), without any additional empirical parameter. Further calculations using SPARC and COSMOtherm for partitioning between water and other tar-components (e.g., benzofuran, phenol and quinoline) suggested that minor components in coal tar do not significantly influence K(CT/w) of nonpolar solutes, and that Raoult's law and the SPLFER thus may be generally applied to these types of solutes, e.g., polycyclic aromatic hydrocarbons and alkylbenzenes, regardless of coal tar compositions. In contrast, partitioning of H-bonding solutes (e.g., phenols) can significantly vary depending on the amount of polar tar-components such as N-heterocyclic aromatic compounds. Therefore, the presented successful applications of Raoult's law and SPLFER to the studied coal tar could be a special case, and these simple approaches may not provide reasonable estimations for partitioning of H-bonding solutes from compositionally different coal tars.     

Correlation relationships between physico-chemical properties and gas chromatographic retention index of polychlorinated-dibenzofurans

Correlation relationships between physico-chemical properties including vapor pressures (P), water solubilities (S), Henry's law constants (H(c)), n-octanol-water partition coefficients (K(ow)), sediment-water partition coefficient (K(pw)) and biotic lipid-water partition coefficient (K(bw), bioconcentration factor) of polychlorinated-dibenzofurans (PCDFs) and their gas chromatographic retention indices (GC-RIs) were established. A model equation between GC-RIs (= RI) and these physico-chemical properties (K) of PCDFs was in a form of log K = aRI2 + bRI + c with correlation coefficients (R2) greater than 0.94, except H(c). These equations were derived from six experimental data (five experimental data for log K(bw)) in each physico-chemical properties of PCDFs reported previously. The values of log P, log S, log H(c), log K(ow), log K(pw) and log K(bw) of PCDFs predicted by these equations based on their GC-RIs in the present study derviated from those calculated by the solubility parameters for fate analysis method in a previous study by 0.49, 0.32, 0.11, 0.34, 0.14 and 0.22 log units, respectively.     

Acute toxicity of substituted phenols to Rana japonica tadpoles and mechanism-based quantitative structure-activity relationship (QSAR) study.

Acute 12 h and 24 h lethal toxicity (12 h-LC50 and 24 h-LC50) of 31 substituted phenols to Rana japonica tadpoles was determined. Results indicate that toxicity of phenols to tadpoles varied only slightly with length of exposure and the 12-h test could serve as surrogate of the 24-h test. A mechanism-based quantitative structure-activity relationship (QSAR) method was employed and 1-octanol/water partition coefficient (log K(ow))-dependent models were developed to study different modes of toxic action. Most phenols elicited their response via a polar narcotic mechanism and an excellent logK(ow)-dependent model was obtained. Soft electrophilicity and pro-electrophilicity were observed for some phenols and a good log K(ow)-dependent model was also achieved. Additionally, the significant dissociation of carboxyl on benzoic acid derivatives sharply reduced their toxicity. A statistically robust QSAR model was developed for all studied compounds with the combined application of log K(ow), energy of lowest unoccupied orbital (E(lumo)), heat of formation (HOF) and the first-order path molecular connectivity dices (1chi(p)).     

Physicochemical properties of selected polybrominated diphenyl ethers and extension of the UNIFAC model to brominated aromatic compounds

The aqueous solubilities (S(w)) at various temperatures from 283 K to 308 K and 1-octanol/water partition coefficients (K(ow)) for four polybrominated diphenyl ethers (PBDEs: 4,4'-dibromodiphenyl ether (BDE-15), 2,2',4,4'-tetrabromodiphenyl ether (BDE-47), 2,2',4,4',5-pentabromodiphenyl ether (BDE-99), and 2,2',4,4',5,5'-hexabromodiphenyl ether (BDE-153)) were measured by the generator column method. The S(w) and K(ow) data revealed the effect of bromine substitution and basic structure on S(w) and K(ow). To estimate the infinite dilution activity coefficients (gamma(i)(w,infinity)) of the PBDEs in water from the S(w) data, enthalpies of fusion and melting points for those compounds were measured with a differential scanning calorimeter. Henry's Law constants (H(w)) of the PBDEs were derived from the determined gamma(i)(w,infinity) and literature vapor pressure data. Some physicochemical characteristics of PBDEs were also suggested by comparing the present property data with that of polychlorinated dibenzo-p-dioxins, brominated phenols and brominated benzenes in past studies. Furthermore, in order to represent different phase equilibria including solubility and partition equilibrium for other brominated aromatic compounds using the UNIFAC model, a pair of UNIFAC group interaction parameters between the bromine and water group were determined from the S(w) and K(ow) data of PBDEs and brominated benzenes. The ability of the determined parameters to represent both properties of brominated aromatics was evaluated.     

Sorption behaviour of acetochlor,atrazine, carbendazim,diazinon, imidacloprid and isoproturon on Hungarian agricultural soil

Laboratory studies were conducted to determine the sorption behaviour of six commonly used pesticides (acetochlor, atrazine, carbendazim, diazinon, imidacloprid and isoproturon) on Hungarian brown forest soil with clay alluviation (Luvisol) using the batch equilibrium technique. The sorption isotherms could be described by the Freundlich equation in non-linear form (n < 1) for all compounds, however in case of diazinon using the extended Freundlich equation proved to be a better approach. The adsorption constant related soil organic carbon content (Koc) calculated from Freundlich equation were 314 for acetochlor, 133 for atrazine, 2805 for carbendazim, 1589 for diazinon, 210 for imidacloprid and 174 for isoproturon. The octanol-water partition coefficients (Pow), which can be a useful parameter to predict of adsorption behaviour of a chemical on soil, and dissociation coefficients of these pesticides were calculated based on the chemical structure of them using a computerized expert system. The octanol-water partition coefficients were determined experimentally from high performance liquid chromatographic parameters as well. Good agreement was observed between experimental and the computer expert system estimated data. Computer estimated log Pow values ranged 0.5 and 3.86 for the examined pesticides, with imidacloprid and diazinon being the least and most hydrophobic respectively. Experimentally determined logPow ranged between 0.92 and 3.81 with the same tendency. It can be concluded that the Freundlich adsorption constants (Kf) are slightly related to the octanol-water partition coefficients of investigated chemicals, nevertheless no close correlation could be established because of the influence of further characteristics of solutes and soil.     

Modeling transport and deposition of contaminants to ecosystems of concern: a case study for the Laurentian Great Lakes

Transfer efficiency (TE) is introduced as a model output that can be used to characterize the relative ability of chemicals to be transported in the environment and deposited to specific target ecosystems. We illustrate this concept by applying the Berkeley-Trent North American contaminant fate model (BETR North America) to identify organic chemicals with properties that result in efficient atmospheric transport and deposition to the Laurentian Great Lakes. By systematically applying the model to hypothetical organic chemicals that span a wide range of environmental partitioning properties, we identify combinations of properties that favor efficient transport and deposition to the Lakes. Five classes of chemicals are identified based on dominant transport and deposition pathways, and specific examples of chemicals in each class are identified and discussed. The role of vegetation in scavenging chemicals from the atmosphere is assessed, and found to have a negligible influence on transfer efficiency to the Great Lakes. Results indicate chemicals with octanol-water (K(ow)) and air-water (K(aw)) partition coefficients in the range of 10(5)-10(7) and 10(-4)-10(-1) combine efficient transport and deposition to the Great Lakes with potential for biaccumulation in the aquatic food web once they are deposited. A method of estimating the time scale for atmospheric transport and deposition process is suggested, and the effects of degrading reactions in the atmosphere and meteorological conditions on transport efficiency of different classes of chemicals are discussed. In total, this approach provides a method of identifying chemicals that are subject to long-range transport and deposition to specific target ecosystems as a result of their partitioning and persistence characteristics. Supported by an appropriate contaminant fate model, the approach can be applied to any target ecosystem of concern.     

Comparison of nonideal sorption formulations in modeling the transport of phthalate esters through packed soil columns

Sorption of dimethyl phthalate (DMP), diethyl phthalate (DEP) and dipropyl phthalate (DPP) to two soil materials that vary in organic matter content was investigated using miscible displacement experiments under saturated flow conditions. Generated breakthrough curves (BTCs) were inversely simulated using linear, equilibrium sorption (LE), nonlinear, equilibrium sorption (NL), linear, first-order nonequilibrium sorption (LFO), linear, radial diffusion (LRD), and nonlinear, first-order nonequilibrium sorption (NFO) models. The Akaike information criterion was utilized to determine the preferred model. The LE model could not adequately describe phthalate ester (PE) BTCs in higher organic matter soil or for more hydrophobic PEs. The LFO and LRD models adequately described the BTCs but a slight improvement in curve-fitting was gained in some cases when the NFO model was used. However, none of the models could properly describe the desorptive tail of DPP for the high organic matter soil. Transport of DPP through this soil was adequately predicted when degradation or sorption hysteresis was considered. Using the optimized parameter values along with values reported by others it was shown that the organic carbon distribution coefficient (K(oc)) of PEs correlates well with the octanol/water partition coefficient (K(ow)). Also, a strong relationship was found between the first-order sorption rate coefficient normalized to injection pulse size and compound residence time. A similar trend of timescale dependence was found for the rate parameter in the radial diffusion model. Results also revealed that the fraction of instantaneous sorption sites is dependent on K(ow) and appears to decrease with the increase in the sorption rate parameter.     

Estimation of selected physicochemical properties for methylated naphthalene compounds

Liquid aqueous solubility (S(w,L)), octanol/water partition coefficients (K(ow)), liquid vapor pressure (P(v,L)), and Henry's law constants (H(c)) were estimated for 20 methylated naphthalenes ranging from monomethyl to tetramethylnaphthalenes. Chromatographic methods were used for the estimation. Chromatographic retention measurements were conducted for 11 reference compounds and regressions were fit between the retention indices and the physicochemical properties. HPLC octadecylsilyl column with acetonitrile/water eluent was used for the estimation of S(w,L) and K(ow). Two GC columns, HP5-MS and a more hydrophobic HP-1, were tested for the estimation of P(v,L). Measured retention indices for the methylated naphthalenes were entered to the regression equations to calculate the physicochemical properties for these compounds. Literature values, where available, were used to validate the calculated values. The method accurately estimated the physicochemical properties. Estimated S(w,L) and P(v,L) decreased with the number of methyl groups. K(ow) increased with the number of methyl groups. There was no obvious relation between H(c) and the number of methyl groups. Log S(w,L) ranged from 0.885 for 1,2,5,6-tetramethylnaphthalene to 2.269 for 1-methylnaphthalene (mmol/m(3)). Log K(ow) varied from 3.89 for 1-methylnaphthalene to 4.95 for 1,2,5,6-tetramethylnaphthalene. Log P(v,L) ranged from -0.983 for 1,2,5,6-tetramethylnaphthalene to 0.789 for 2-methylnaphthalene (Pa). Log H(c) varied from 1.03 for 1,4,5-trimethylnaphthalene to 1.73 for 2,6-dimethylnaphthalene (Pa m(3)/mol). There were no apparent effects of GC column hydrophobicity on the accuracy of the results. Estimation of S(w,L) and K(ow) based on GC retention indices was not as accurate as with HPLC. Comparison of the estimated values with values predicted by EPIWIN indicated that EPIWIN is useful in giving order-of-magnitude prediction of physicochemical properties.     

The relationship between elimination rates and partition coefficients

Rate constants for uptake and elimination of chemicals in organisms are often related to partition coefficients (typically the octanol-water partition coefficient). We show that the well-mixed one-compartment model for toxico-kinetics implies that the elimination rate is inversely proportional to the square root of the partition coefficient. When chemical exchange is limited by diffusion in the boundary layers adjacent to the interface, two-film models are appropriate, which have more complex implications for the relationships between the exchange rates and the partition coefficient. We also show that the popular steady-flux approximation of the two-film model is not a conceptual generalization of the one-compartment model, although it shares the first-order kinetics. We compare the kinetics of a series of models with an increasing number of well-mixed compartments for exchange, such that the two-film model results for an infinite number of compartments. The latter model formulation in terms of partial differential equations, and more in particular its boundary condition at the interface of the two media, is believed to be new. In the steady-flux approximation and in the model with single well-mixed boundary layers and low diffusivities, the elimination rate depends hyperbolically on the partition coefficient. The available data for abiotic systems (SPME fibers) supports a hyperbolic relationship, whereas the data for aquatic biota are less discriminating between a hyperbolic or a square root relationship with the partition coefficient. The daphnia data showed less scatter than the fish data, possibly due to the small variance in body sizes, since elimination rates are inversely proportional to body length. The square root relationship fitted these data best.     

生物滴滤塔净化苯乙烯废气的实验研究

采用生物滴滤(BTF)系统对含苯乙烯的有机废气进行了生物净化实验并研究该系统VOCs生物降解性能。实验表明,苯乙烯进气浓度低于20 mg/m3时BTF去除效率可达92%以上,出口苯乙烯浓度低于1.6 mg/m3,达到GB14554-1993中规定的排放标准;该BTF装置对苯乙烯的去除负荷在2.0 g/(m3.h)左右;系统稳定运行时循环液COD、浊度和pH等都保持稳定,无脱落生物膜积累现象;生物滴滤塔系统适宜的气液比为300;系统总压降约100 Pa,鲍尔环填料和聚氨酯发泡填料混合装填方式可以降低系统压降并有利于微生物挂膜。     

Accumulation of organochlorine pesticides by semipermeable membrane devices using composite complex

Semipermeable membrane devices (SPMDs) were developed for passive in situ monitoring of organochlorine pesticides (OCPs) in aqueous solution in both laboratory and field (Pearl River Delta, China) studies. The device consisted of a thin film of neutral lipid triolein, enclosed in thin-walled tubing made of composite cellulose acetate membrane (CA) supported by linear low density polyethylene (LLDPE) (CAPE). Results from the laboratory and field application indicated that triolein-CAPE (TCAPE) could quickly and efficiently accumulate hydrophobic OCPs in water and uptake equilibrium could reached within 20h in the laboratory. Some mathematical relationships of TCAPE-water partition coefficient (logK(sw)), triolein-water partition coefficient (logK(tw)) and octanol-water partition coefficient (logK(ow)) were developed under the laboratory conditions. A good correlation of accumulation in TCAPE with r(2) values ranging from 0.55 to 0.86 for individual OCPs (n=8) and an excellent correlation of logK(sw) and logK(ow) was also obtained under the field conditions. The average OCPs concentration in the surface water could be estimated by measuring OCPs concentration in the device under the field conditions.     

A screening method for preliminary assessment of risk to groundwater from land-applied chemicals

A simple mathematical model for initial screening is presented that can aid in evaluating the relative risk to groundwater from applying nonpolar synthetic organic chemicals to soil. The basic premise is that the magnitude of the quotient of the chemical concentration of the water entering the aquifer and the maximum allowable concentration (as established by EPA or Health Departments) represents the health risk of a chemical. The chemical concentration of the soil water is estimated based on conservative, simplifying assumptions and requires only readily available data such as: basic soil properties (organic matter and saturated hydraulic conductivity), organic chemical properties (octanol-water partition coefficient and degradation rate) and environmental factors (recharge rate and depth to groundwater).The methodology was applied to assess the relative risk of organic chemicals in municipal sewage sludge and pesticides applied to agricultural land. The results are realistic.     

Structure-activity relationships and response-surface analysis of nitroaromatics toxicity to the yeast (Saccharomyces cerevisiae)

Inhibition of growth of the yeast Saccharomyces cerevisiae (Cmiz, the minimum concentration that produced a clear inhibition zone within 12 h) for 24 nitroaromatic compounds was investigated and a quantitative structure-activity relationship (QSAR) developed based on hydrophobicity expressed as the l-octanol/water partition coefficient in logarithm form, log K(ow), electrophilicity based on the energy of the lowest unoccupied orbital (E(lumo)). All nitrobenzene derivatives exhibited enhanced reactive toxicity than baseline. The toxicities of mono-nitrobenzenes and di-nitrobenzenes were elicited by different mechanisms of toxic action. For mono-nitro-derivatives, both significant log K(ow) based and strong E(lumo)-dependent relationships were observed indicating that their toxicities were affected both by the penetration process and the interaction with target sites of interaction. The toxicities of di-nitrobenzenes were greater than mono-nitrobenzenes and no log K(ow)-dependent but highly significant E(lumo)-based relationship was obtained. This suggests that toxicity of di-nitrobenzenes was highly electrophilic and involved mainly their in vivo electrophilic interaction with biomacromolecules. In an effort to model the elevated toxicity of all nitrobenzenes, a response-surface analysis was performed and this resulted in a highly predictive two-variable QSAR without reference to their exact mechanisms (Cmiz = 0.41 log K(ow) - 0.89 E(lumo) - 0.46, r2 = 0.87, Q2 = 0.86, n = 24).