Aqueous phase hydration and hydrate acidity of perfluoroalkyl and n:2 fluorotelomer aldehydes.

The SPARC software program and comparative density functional theory (DFT) calculations were used to investigate the aqueous phase hydration equilibrium constants (Khyd) of perfluoroalkyl aldehydes (PFAlds) and n:2 fluorotelomer aldehydes (FTAlds). Both classes are degradation products of known industrial compounds and environmental contaminants such as fluorotelomer alcohols, iodides, acrylates, phosphate esters, and other derivatives, as well as hydrofluorocarbons and hydrochlorofluorocarbons. Prior studies have generally failed to consider the hydration, and subsequent potential hydrate acidity, of these compounds, resulting in incomplete and erroneous predictions as to their environmental behavior. In the current work, DFT calculations suggest that all PFAlds will be dominantly present as the hydrated form in aqueous solution. Both SPARC and DFT calculations suggest that FTAlds will not likely be substantially hydrated in aquatic systems or in vivo. PFAld hydrates are expected to have pKa values in the range of phenols (ca. 9 to 10), whereas n:2 FTAld hydrates are expected to have pKa values ca. 2 to 3 units higher (ca. 12 to 13). In order to avoid spurious modeling predictions and a fundamental misunderstanding of their fate, the molecular and/or dissociated hydrate forms of PFAlds and FTAlds need to be explicitly considered in environmental, toxicological, and waste treatment investigations. The results of the current study will facilitate a more complete examination of the environmental fate of PFAlds and FTAlds.

Performance of the major semiempirical, ab initio, and DFT methods for isomerization enthalpies of linear to branched heptanes.

The gas phase standard state (298.15 K, 1 atm) isomerization enthalpy (Δ(isom)H°(g)) prediction performance of the major semiempirical, ab initio, and density functional levels of theory for environmentally relevant transformations was investigated using the linear to branched heptanes as a representative case study. The M062X density functional, MP2 (and higher) levels of Moller-Plesset perturbation theory, and the CBS and Gaussian-n composite methods are well suited for investigating the thermodynamic properties of environmentally interesting isomerizations, although the M062X functional may not be appropriate for all systems. Where large molecular systems prohibit the use of higher levels of theory, the PM6 and PDDG semiempirical methods may offer an appropriate computational cost-accuracy compromise.

Air-water partition coefficients for a suite of polycyclic aromatic and other C10 through C20 unsaturated hydrocarbons.

The air-water partition coefficients (Kaw) for 86 large polycyclic aromatic hydrocarbons and their unsaturated relatives were estimated using high-level G4(MP2) gas and aqueous phase calculations with the SMD, IEFPCM-UFF, and CPCM solvation models. An extensive method validation effort was undertaken which involved confirming that, via comparisons to experimental enthalpies of formation, gas-phase energies at the G4(MP2) level for the compounds of interest were at or near thermochemical accuracy. Investigations of the three solvation models using a range of neutral and ionic compounds suggested that while no clear preferential solvation model could be chosen in advance for accurate Kaw estimates of the target compounds, the employment of increasingly higher levels of theory would result in lower Kaw errors. Subsequent calculations on the polycyclic aromatic and unsaturated hydrocarbons at the G4(MP2) level revealed excellent agreement for the IEFPCM-UFF and CPCM models against limited available experimental data. The IEFPCM-UFF-G4(MP2) and CPCM-G4(MP2) solvation energy calculation approaches are anticipated to give Kaw estimates within typical experimental ranges, each having general Kaw errors of less than 0.5 log10 units. When applied to other large organic compounds, the method should allow development of a broad and reliable Kaw database for multimedia environmental modeling efforts on various contaminants.

Use of the SPARC software program to calculate hydrolysis rate constants for the polymeric brominated flame retardants BC-58 and FR-1025.

The SPARC software program was used to estimate the acid-catalyzed, neutral, and base-catalyzed hydrolysis rate constants for the polymeric brominated flame retardants BC-58 and FR-1025. Relatively rapid hydrolysis of BC-58, producing 2,4,6-tribromophenol-and ultimately tetrabromobisphenol A-as the hydrolytically stable end products from all potential hydrolysis reactions, is expected in both environmental and biological systems with starting material hydrolytic half-lives (t(1/2,hydr)) ranging from less than 1 h in marine systems, several hours in cellular environments, and up to several weeks in slightly acid fresh waters. Hydrolysis of FR-1025 to give 2,3,4,5,6-pentabromobenzyl alcohol is expected to be slower (t(1/2,hydr) less than 0.5 years in marine systems up to several years in fresh waters) than BC-58, but is also expected to occur at rates that will contribute significantly to environmental and in vivo loadings of this compound.

Estimated pKa values for the environmentally relevant C1 through C8 perfluorinated sulfonic acid isomers.

In order to estimate isomer-specific acidity constants (pKa) for the perfluorinated sulfonic acid (PFSA) environmental contaminants, the parameterization method 6 (PM6) pKa prediction method was extensively validated against a wide range of carbon oxyacids and related sulfonic/sulfinic acids. Excellent pKa prediction performance was observed for the carbon oxyacids using the PM6 method, but this approach was found to have a severe positive bias for sulfonic/sulfinic acids. To overcome this obstacle, a correlation was developed between non-adjusted PM6 pKa values and the corresponding experimentally obtained/estimated acidity constants for a range of representative alkyl, aryl and halogen-substituted sulfonic acids. Application of this correction to the PM6 values allows for extension of this computational method to a new acid functional group. When used to estimate isomer-specific pKa values for the C1 through C8 PFSAs, the modified PM6 approach suggests an adjusted pKa range from -5.3 to -9.0, indicating that all members of this class of well-known environmental contaminants will be effectively completely dissociated in aquatic systems.

Thermodynamic properties of chloramine formation and related reactions during water treatment: a G4MP2, G4, and W1BD theoretical study.

A high-level gas and aqueous phase theoretical thermodynamic study was conducted on the primary and related chemical reactions which occur during chloramination for water treatment using the G4MP2, G4, and W1BD composite methods with the SMD, PCM, and CPCM solvation models. The standard state (298.15 K, 1 atm or 1M) formation of mono-, di-, and tri-chloramines from their precursors via hypochlorous acid chlorination is substantially exothermic and exergonic in both the gas and aqueous phases. The excellent agreement between experimental and theoretical values for a range of structural and thermodynamic calculations on a suite of calibration compounds suggests that the G4MP2, G4, and W1BD calculations meet or exceed criteria for thermochemical accuracy. The temperature influence on the thermodynamics of chloramine formation is projected to be negligible regardless of phase between 0 and 100°C. Additional thermodynamic calculations were undertaken on associated chloramination reactions involving the disproportionation of monochloramine, the decomposition of di- and tri-chloramine, and the reactions of trichloramine with ammonia and dichloramine. The results from these investigations not only provide a better understanding of the reaction thermodynamics, they also allow for a more rigorous interpretation of proposed chloramination mechanisms.

Prediction of the air-water partition coefficient for perfluoro-2-methyl-3-pentanone using high-level Gaussian-4 composite theoretical methods.

The air-water partition coefficient (Kaw) of perfluoro-2-methyl-3-pentanone (PFMP) was estimated using the G4MP2/G4 levels of theory and the SMD solvation model. A suite of 31 fluorinated compounds was employed to calibrate the theoretical method. Excellent agreement between experimental and directly calculated Kaw values was obtained for the calibration compounds. The PCM solvation model was found to yield unsatisfactory Kaw estimates for fluorinated compounds at both levels of theory. The HENRYWIN Kaw estimation program also exhibited poor Kaw prediction performance on the training set. Based on the resulting regression equation for the calibration compounds, the G4MP2-SMD method constrained the estimated Kaw of PFMP to the range 5-8 × 10(-6) M atm(-1). The magnitude of this Kaw range indicates almost all PFMP released into the atmosphere or near the land-atmosphere interface will reside in the gas phase, with only minor quantities dissolved in the aqueous phase as the parent compound and/or its hydrate/hydrate conjugate base. Following discharge into aqueous systems not at equilibrium with the atmosphere, significant quantities of PFMP will be present as the dissolved parent compound and/or its hydrate/hydrate conjugate base.

pK(a) values of the monohydroxylated polychlorinated biphenyls (OH-PCBs), polybrominated biphenyls (OH-PBBs), polychlorinated diphenyl ethers (OH-PCDEs), and polybrominated diphenyl ethers (OH-PBDEs).

The SPARC software program aqueous pK(a) prediction module was validated against corresponding experimental acidity constants for chlorinated and brominated phenols and the limited experimental aqueous pK(a) data sets for monohydroxylated polychlorinated biphenyls (OH-PCBs), polychlorinated diphenyl ethers (OH-PCDEs), and polybrominated diphenyl ethers (OH-PBDEs). pK(a) values were then estimated for all 837 monohydroxylated mono- through nona-halogenated congeners in each of the OH-PCB, OH-PCDE, and OH-PBDE classes, as well as for the monohydroxylated polybrominated biphenyls (OH-PBBs), giving a total of 3348 compounds. Large intrahomolog pK(a) variation by up to six units is expected within each contaminant class, with pK(a) values ranging from about 4 to 11 dependent on the degree and pattern of halogenation. Increasing halogenation generally decreased the average pK(a) within each homolog group. Significant intrahomolog differences in pK(a) values exist between OH-PCB, OH-PBB, OH-PCDE, and OH-PBDE congeners, including large acidity constant variation between isomers with equivalent halogenation patterns but varying location of the hydroxy moiety. Congener specific pH dependent investigations into the partitioning and degradation behaviors of these compounds are necessary, including greater consideration of analyte ionization effects during their extraction and analysis.

Dow and Kaw,eff vs. Kow and Kaw degrees: acid/base ionization effects on partitioning properties and screening commercial chemicals for long-range transport and bioaccumulation potential.

A set of 543 ionizable commercial organic compounds with various acid/base functionalities and experimental octanol-water partitioning coefficients (log Kow) were obtained from the Canadian Domestic Substances List. Corresponding pH-dependent octanol-water distribution coefficients (log Dow) and air-water partitioning coefficients (log Kaw,eff) were estimated using the SPARC software program, as were log Kow and log Kaw degrees values for the neutral forms of each chemical. Significant ionization dependent effects on chemical screening results at various pH values were obtained using established criteria for bioaccumulation potential (BAP) in aquatic organisms, terrestrial animals, and humans, as well as for atmospheric long range transport potential (LRTP). Future modelling efforts for environmental and toxicological screening of commercial chemicals should therefore explicitly include the influence of ionization for both weak and strong organic acids and bases on bioavailability and air-water mobility within the respective regulatory frameworks. Functional group specific sorption of both ionizable and neutral compounds to particulate and dissolved inorganic and organic matter will also affect chemical screening results for BAP and LRTP. More complex sorption related modelling in various types of representative aquatic systems also appears necessary to achieve reliable chemical screening results for commercial organic compounds.

Quantitative structure-activity relationship (QSAR) studies for predicting activation of the ryanodine receptor type 1 channel complex (RyR1) by polychlorinated biphenyl (PCB) congeners.

A quantitative structure-activity relationship (QSAR) was developed to predict the congener specific ryanodine receptor type RyR1 activity of all 209 polychlorinated biphenyl (PCB) congeners. A three-variable QSAR equation was obtained via stepwise forward linear regression on an unsupervised forward selection reduced data set from an initial database. Application of the QSAR towards predicting EC(2x) values for all 209 PCB congeners indicated good agreement in substitution pattern trends between the experimental and estimated data sets. The QSAR model predicts a less than two-fold increase in maximal potency among all congeners outside the experimental database, and it appears that no high-potency PCB congeners with EC(2x) values much less than 0.2 microM exist. Increasing RyR1-neuro toxicity equivalents with increasing homologue number and Aroclor chlorination likely reflect indirect molecular controls on toxicity, since congeners with multiple ortho substituents-the primary structural feature controlling a lack of coplanarity and resulting neurotoxicity-are more likely to be found in higher homologues.

Modeling the hydrolysis of perfluorinated compounds containing carboxylic and phosphoric acid ester functions and sulfonamide groups.

Temperature-dependent rate constants were estimated for the acid- and base-catalyzed and neutral hydrolysis reactions of perfluorinated telomer acrylates (FTAcrs) and phosphate esters (FTPEs), and the S(N)1 and S(N)2 hydrolysis reactions of fluorotelomer iodides (FTIs). Under some environmental conditions, hydrolysis of monomeric FTAcrs could be rapid (half-lives of several years in marine systems and as low as several days in some landfills) and represent a dominant portion of their overall degradation. Abiotic hydrolysis of monomeric FTAcrs may be a significant contributor to current environmental loadings of fluorotelomer alcohols (FTOHs) and perfluoroalkyl carboxylic acids (PFCAs). Polymeric FTAcrs are expected to be hydrolyzed more slowly, with estimated half-lives in soil and natural waters ranging between several centuries to several millenia absent additional surface area limitations on reactivity. Poor agreement was found between the limited experimental data on FTPE hydrolysis and computational estimates, requiring more detailed experimental data before any further modeling can occur on these compounds or their perfluoroalkyl sulfonamidoethanol phosphate ester (PFSamPE) analogs. FTIs are expected to have hydrolytic half-lives of about 130 days in most natural waters, suggesting they may be contributing to substantial FTOH and PFCA inputs in aquatic systems. Perfluoroalkyl sulfonamides (PFSams) appear unlikely to undergo abiotic hydrolysis at the S-N, C-S, or N-C linkages under environmentally relevant conditions, although potentially facile S-N hydrolysis via intramolecular catalysis by ethanol and acetic acid amide substituents warrants further investigation.

Congener-specific organic carbon-normalized soil and sediment-water partitioning coefficients for the C1 through C8 perfluoroalkyl carboxylic and sulfonic acids.

Organic carbon-normalized soil and sediment-water partitioning coefficients (K(oc)) were estimated for all C(1) through C(8) perfluoroalkyl carboxylic (PFCA) and sulfonic (PFSA) acid congeners. The limited experimental K(oc) data set for the straight chain C(7) through C(10) PFCAs and C(8) and C(10) PFSAs was correlated to SPARC and ALOGPS computationally estimated octanol-water partitioning/distribution constants and used to predict K(oc) values for both branched and linear C(1) through C(8) isomers. Branched and linear congeners in this homologue range are generally expected to have K(oc) values > 1, leading to their accumulation in organic matter on sediments and soils, retardation during ground and pore water flow, and the preferential association with dissolved organic matter in aquatic systems. Both increasing perfluoroalkyl chain length and linearity increase K(oc) values with substantial intra- and inter-homologue variation and interhomologue mixing. Variability in K(oc) values among the PFCA and PFSA congeners will likely lead to an enrichment of more linear and longer-chain isomers in organic matter fractions, resulting in aqueous phases fractionated towards shorter-chain branched congeners. The expected magnitude of fractionation will require inclusion in source apportionment models and risk assessments. A comparison of representative established quantitative structure property relationships for estimating K(oc) values from octanol-water partitioning constants suggests that these equilibrium partitioning frameworks may be applicable towards modeling PFCA and PFSA environmental fate processes.

A new class of perfluorinated acid contaminants: primary and secondary substituted perfluoroalkyl sulfonamides are acidic at environmentally and toxicologically relevant pH values.

The SPARC software program was validated for nitrogen-hydrogen acidity constant estimation of primary and secondary sulfonamides against a broad suite of substituted derivatives with experimental datasets in water and dimethylsulfoxide solvent systems and across a wide pK(a) range. Following validation, amidic proton pK(a) values were estimated for all C(1) through C(8) congeners of five major perfluoroalkyl sulfonamide classes: unsubstituted sulfonamides, N-methyl and N-ethyl sulfonamides, sulfonamidoethanols, and sulfonamidoacetates. Branching of the perfluoroalkyl chain is expected to have substantial impacts on amide moiety acidity in these contaminant groups, with intrahomologue variability of up to four pK(a) units and increasing pK(a) values with both increasing chain branching and greater proximity of the chain branching to the sulfonamide head group. Perfluoroalkyl chain length is not predicted to have a substantial influence on sulfonamide acidity. The predicted pK(a) values and variability are anticipated to have substantial impacts on the environmental partitioning and degradation of these compounds, as well as the modes and magnitudes of toxicological effects. Substantial pH dependent isomeric fractionation of perfluoroalkyl sulfonamides is expected both in situ and in vivo, necessitating the incorporation of amide group acidities in multimedia environmental models and pharmacokinetic studies.