Further transformation of the primary ozonation products of tramadol- and venlafaxine N-oxide: Mechanistic and structural considerations.

Ozonation has been used to effectively remove micropollutants from the secondary effluent in several wastewater treatment plants. It is known that ozonation transforms tertiary amine compounds into their respective N-oxides, however in an earlier study a mass balance could not be closed at elevated ozone concentrations, leading to the assumption that more ozonation products are possible. This study was conducted to elucidate which (hitherto unknown) ozonation products can be formed from venlafaxine and tramadol when ozonating wastewater. Ozonation experiments were performed with tramadol and venlafaxine N-oxide in two different set-ups. Both tramadol- and venlafaxine N-oxide degraded during ozonation in pure (deionized) water in both semi-continuous and batch mode ozonation set-ups. 13 and 17 new transformation products were detected from tramadol- and venlafaxine N-oxide respectively, using high resolution mass spectrometry with ESI(+) ionization. Empirical chemical formulas were proposed based on the determination of the exact masses and interpretation of the product ion spectra. These transformation products result from the addition of one to three oxygen atoms and removal of C, -CH2, C2H2, C3H6, etc., from the parent molecule, respectively. Quenching experiments suggested that most of the transformation products originated from the direct reaction with ozone (eight for tramadol N-oxide and ten for venlafaxine N-oxide), whereas fewer products originated from the reaction with OH radicals (three for tramadol N-oxide and three for venlafaxine N-oxide). Reaction mechanisms and chemical structures of products are proposed, based on the available active sites and past literature on ozone reaction mechanisms. The experimental results are compared to theory and literature on ozone reactive sites and ozone reaction mechanisms. All in all this shows that there can be multiple ozonation products, and ozonation pathways can be complex, even if initially only one ozonation product is formed.

Toxic effects of substituted p-benzoquinones and hydroquinones in in vitro bioassays are altered by reactions with the cell assay medium.

Substituted para-benzoquinones and hydroquinones are ubiquitous transformation products that arise during oxidative water treatment of phenolic precursors, for example through ozonation or chlorination. The benzoquinone structural motive is associated with mutagenicity and carcinogenicity, and also with induction of the oxidative stress response through the Nrf2 pathway. For either endpoint, toxicological data for differently substituted compounds are scarce. In this study, oxidative stress response, as indicated by the AREc32 in vitro bioassay, was induced by differently substituted para-benzoquinones, but also by the corresponding hydroquinones. Bioassays that indicate defense against genotoxicity (p53RE-bla) and DNA repair activity (UmuC) were not activated by these compounds. Stability tests conducted under incubation conditions, but in the absence of cell lines, showed that tested para-benzoquinones reacted rapidly with constituents of the incubation medium. Compounds were abated already in phosphate buffer, but even faster in biological media, with reactions attributed to amino- and thiol-groups of peptides, proteins, and free amino acids. The products of these reactions were often the corresponding substituted hydroquinones. Conversely, differently substituted hydroquinones were quantitatively oxidized to p-benzoquinones over the course of the incubation. The observed induction of the oxidative stress response was attributed to hydroquinones that are presumably oxidized to benzoquinones inside the cells. Despite the instability of the tested compounds in the incubation medium, the AREc32 in vitro bioassay could be used as an unspecific sum parameter to detect para-benzoquinones and hydroquinones in oxidatively treated waters.

Micropollutant Oxidation Studied by Quantum Chemical Computations: Methodology and Applications to Thermodynamics, Kinetics, and Reaction Mechanisms.

The abatement of organic micropollutants during oxidation processes has become an emerging issue for various urban water systems such as drinking water, wastewater, and water reuse. Reaction kinetics and mechanisms play an important role in terms of efficiency of these processes and the formation of transformation products, which are controlled by functional groups in the micropollutants and the applied oxidants. So far, the kinetic and mechanistic information on the underlying reactions was obtained by experimental studies; additionally, predictive quantitative structure-activity relationships (QSARs) were applied to determine reaction kinetics for the oxidation of emerging compounds. Since this experimental approach is very laborious and there are tens of thousands potential contaminants, alternative strategies need to be developed to predict the fate of micropollutants during oxidative water treatment. Due to significant developments in quantum chemical (QC) computations in recent years and increased computational capacity, QC-based methods have become an alternative or a supplement to the current experimental approach. This Account provides a critical assessment of the current state-of-the-art of QC-based methods for the assessment of oxidation of micropollutants. Starting from a given input structure, QC computations need to locate energetic minima on the potential energy surface (PES). Then, useful thermodynamic and kinetic information can be estimated by different approaches: Experimentally determined reaction mechanisms can be validated by identification of transition structures on the PES, which can be obtained for addition reactions, heavy atom transfer (Cl+, Br+, O·) and H atom transfer (simultaneous proton and electron transfer) reactions. However, transition structures in the PES cannot be obtained for e--transfer reactions. Second-order rate constants k for the reactions of micropollutants with chemical oxidants can be obtained by ab initio calculations or by QSARs with various QC descriptors. It has been demonstrated that second-order rate constants from ab initio calculations are within factors 3-750 of the measured values, whereas QSAR-based methods can achieve factors 2-4 compared to the experimental data. The orbital eigenvalue of the highest occupied molecular orbital ( EHOMO) is the most commonly used descriptor for QSAR-based computations of k-values. In combination with results from experimental studies, QC computations can also be applied to investigate reaction mechanisms for verification/understanding of oxidative mechanisms, calculation of branching ratios or regioselectivity, evaluation of the experimental product distribution and assessment of substitution effects. Furthermore, other important physical-chemical constants such as unknown equilibria for species, which are not measurable due to low concentrations, or p Ka values of reactive transient species can be estimated. With further development of QC-based methods, it will become possible to implement kinetic and mechanistic information from such computations in in silico models to predict oxidative transformation of micropollutants. Such predictions can then be complemented by tailored experimental studies to confirm/falsify the computations.

Gas-Phase Ozone Reactions with a Structurally Diverse Set of Molecules: Barrier Heights and Reaction Energies Evaluated by Coupled Cluster and Density Functional Theory Calculations.

Reactions with ozone transform organic and inorganic molecules in water treatment systems as well as in atmospheric chemistry, either in the aqueous phase, at gas/particle interfaces, or in the gas phase. Computed thermokinetic data can be used to estimate the reactivities of molecules toward ozone in cases where no experimental data are available. Although the gas-phase reactivity of olefins with ozone has been characterized extensively in the literature, this is not the case for the richer chemistry of ozone with polar molecules, which occurs in the aqueous phase or in microhydrated environments. Here, we selected a number of model reactions with small molecules (ethene, ethyne, hydrogen cyanide, hydrogen chloride, ammonia, bromide, and trimethylamine) to study the accuracy of different quantum chemical methods for describing the reactivities of these molecules with ozone. We calculated benchmark electronic energies of gas-phase reactions of these systems with single-reference coupled cluster (CC) theory. These benchmark results for the binding energy in the van der Waals complex, the energy of the transition structure, and the reaction energy were estimated to be accurate within 1-2 kcal mol-1. Singlet oxygen (1O2) is a common product of ozone reactions. Coupled cluster calculations with up to perturbative quadruples (CCSDT( Q)) were needed to obtain reaction energies accurate within 1 kcal mol-1 when this species was involved. In (micro)hydrated environments or at interfaces, coupled cluster methods are prohibitively expensive in most cases. We tested the suitability of some contemporary density functional theory (DFT) methods to reproduce the benchmark electronic energy differences. Range-separated functionals were found to be promising candidates to estimate forward barrier heights, with LC-ωPBE rivaling the accuracy of CCSD( T). For energies of reaction, however, DFT methods exhibited large systematic errors, depending on their fraction of orbital exchange. This was found to worsen when 1O2 is a product, and no safe recommendation can be given for DFT reaction energies in such cases.

Ozonation of Para-Substituted Phenolic Compounds Yields p-Benzoquinones, Other Cyclic α,ß-Unsaturated Ketones, and Substituted Catechols.

Phenolic moieties are common functional groups in organic micropollutants and in dissolved organic matter, and are exposed to ozone during drinking water and wastewater ozonation. Although unsubstituted phenol is known to yield potentially genotoxic p-benzoquinone during ozonation, little is known about the effects of substitution of the phenol ring on transformation product formation. With batch experiments employing differing ozone/target compound ratios, it is shown that para-substituted phenols ( p-alkyl, p-halo, p-cyano, p-methoxy, p-formyl, p-carboxy) yield p-benzoquinones, p-substituted catechols, and 4-hydroxy-4-alkyl-cyclohexadien-1-ones as common ozonation products. Only in a few cases did para-substitution prevent the formation of these potentially harmful products. Quantum chemical calculations showed that different reaction mechanisms lead to p-benzoquinone, and that cyclohexadienone can be expected to form if no such pathway is possible. These products can thus be expected from most phenolic moieties. Kinetic considerations showed that substitution of the phenolic ring results in rather small changes of the apparent second order rate constants for phenol-ozone reactions at pH 7. Thus, in mixtures, most phenolic structures can be expected to react with ozone. However, redox cross-reactions between different transformation products, as well as hydrolysis, can be expected to further alter product distributions under realistic treatment scenarios.

Explicit solvent simulations of the aqueous oxidation potential and reorganization energy for neutral molecules: gas phase, linear solvent response, and non-linear response contributions.

First principles simulations were used to predict aqueous one-electron oxidation potentials (Eox) and associated half-cell reorganization energies (λaq) for aniline, phenol, methoxybenzene, imidazole, and dimethylsulfide. We employed quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations of the oxidized and reduced species in an explicit aqueous solvent, followed by EOM-IP-CCSD computations with effective fragment potentials for diabatic energy gaps of solvated clusters, and finally thermodynamic integration of the non-linear solvent response contribution using classical MD. A priori predicted Eox and λaq values exhibit mean absolute errors of 0.17 V and 0.06 eV, respectively, compared to experiment. We also disaggregate Eox into several well-defined free energy properties, including the gas phase adiabatic free energy of ionization (7.73 to 8.82 eV), the solvent-induced shift in the free energy of ionization due to linear solvent response (-2.01 to -2.73 eV), and the contribution from non-linear solvent response (-0.07 to -0.14 eV). The linear solvent response component is further apportioned into contributions from the solvent-induced shift in vertical ionization energy of the reduced species (ΔVIEaq) and the solvent-induced shift in negative vertical electron affinity of the ionized species (ΔNVEAaq). The simulated ΔVIEaq and ΔNVEAaq are found to contribute the principal sources of uncertainty in computational estimates of Eox and λaq. Trends in the magnitudes of disaggregated solvation properties are found to correlate with trends in structural and electronic features of the solute. Finally, conflicting approaches for evaluating the aqueous reorganization energy are contrasted and discussed, and concluding recommendations are given.

Exploring the aqueous vertical ionization of organic molecules by molecular simulation and liquid microjet photoelectron spectroscopy.

To study the influence of aqueous solvent on the electronic energy levels of dissolved organic molecules, we conducted liquid microjet photoelectron spectroscopy (PES) measurements of the aqueous vertical ionization energies (VIEaq) of aniline (7.49 eV), veratrole alcohol (7.68 eV), and imidazole (8.51 eV). We also reanalyzed previously reported experimental PES data for phenol, phenolate, thymidine, and protonated imidazolium cation. We then simulated PE spectra by means of QM/MM molecular dynamics and EOM-IP-CCSD calculations with effective fragment potentials, used to describe the aqueous vertical ionization energies for six molecules, including aniline, phenol, veratrole alcohol, imidazole, methoxybenzene, and dimethylsulfide. Experimental and computational data enable us to decompose the VIEaq into elementary processes. For neutral compounds, the shift in VIE upon solvation, ΔVIEaq, was found to range from ≈-0.5 to -0.91 eV. The ΔVIEaq was further explained in terms of the influence of deforming the gas phase solute into its solution phase conformation, the influence of solute hydrogen-bond donor and acceptor interactions with proximate solvent molecules, and the polarization of about 3000 outerlying solvent molecules. Among the neutral compounds, variability in ΔVIEaq appeared largely controlled by differences in solute-solvent hydrogen-bonding interactions. Detailed computational analysis of the flexible molecule veratrole alcohol reveals that the VIE is strongly dependent on molecular conformation in both gas and aqueous phases. Finally, aqueous reorganization energies of the oxidation half-cell ionization reaction were determined from experimental data or estimated from simulation for the six compounds aniline, phenol, phenolate, veratrole alcohol, dimethylsulfide, and methoxybenzene, revealing a surprising constancy of 2.06 to 2.35 eV.

On the nature of interactions of radicals with polar molecules.

Solvated radicals play an important role in many areas of chemistry, but to date, the nature of their interactions with polar solvent molecules lacks chemical interpretation. We present a computational quantum chemical analysis of the binding motives of binary complexes involving electron-poor and electron-rich radicals bound to water and hydrogen fluoride, considered here as model polar solvent molecules. By comparing the binding strengths of several open-shell and closed-shell complexes, in combination with natural localized molecular orbital analysis, we show that open-shell complexes can exhibit additional donor-acceptor interactions relative to analogous closed-shell systems. This may explain the unexpectedly large binding energies observed in some open-shell complexes. These exploratory results show that specific interactions in open-shell systems deserve more attention, and they imply that the quantum mechanical description of explicit solvent molecules needs to be considered carefully when designing simulation protocols for solvated radicals.

Aqueous oxidation of sulfonamide antibiotics: aromatic nucleophilic substitution of an aniline radical cation.

Sulfonamide antibiotics are an important class of organic micropollutants in the aquatic environment. For several, sulfur dioxide extrusion products have been previously reported upon photochemical or dark oxidation. Using quantum chemical modeling calculations and transient absorption spectroscopy, it is shown that single-electron oxidation from sulfadiazine produces the corresponding aniline radical cation. Density functional theory calculations indicate that this intermediate can exist in four protonation states. One species exhibits a low barrier for an intramolecular nucleophilic attack at the para position of the oxidized aniline ring, in which a pyrimidine nitrogen acts as a nucleophile. This attack can lead to a rearranged structure, which exhibits the same connectivity as the SO2 -extruded oxidation product that was previously observed in the aquatic environment and characterized by NMR spectroscopy. We report a detailed reaction mechanism for this intramolecular aromatic nucleophilic substitution, and we discuss the possibility of this reaction pathway for other sulfonamide drugs.

Reactivity of BrCl, Br2, BrOCl, Br2O, and HOBr toward dimethenamid in solutions of bromide + aqueous free chlorine.

HOBr, formed via oxidation of bromide by free available chlorine (FAC), is frequently assumed to be the sole species responsible for generating brominated disinfection byproducts (DBPs). Our studies reveal that BrCl, Br(2), BrOCl, and Br(2)O can also serve as brominating agents of the herbicide dimethenamid in solutions of bromide to which FAC was added. Conditions affecting bromine speciation (pH, total free bromine concentration ([HOBr](T)), [Cl(-)], and [FAC](o)) were systematically varied, and rates of dimethenamid bromination were measured. Reaction orders in [HOBr](T) ranged from 1.09 (±0.17) to 1.67 (±0.16), reaching a maximum near the pK(a) of HOBr. This complex dependence on [HOBr](T) implicates Br(2)O as an active brominating agent. That bromination rates increased with increasing [Cl(-)], [FAC](o) (at constant [HOBr](T)), and excess bromide (where [Br(-)](o)>[FAC](o)) implicate BrCl, BrOCl, and Br(2), respectively, as brominating agents. As equilibrium constants for the formation of Br(2)O and BrOCl (aq) have not been previously reported, we have calculated these values (and their gas-phase analogues) using benchmark-quality quantum chemical methods [CCSD(T) up to CCSDTQ calculations plus solvation effects]. The results allow us to compute bromine speciation and hence second-order rate constants. Intrinsic brominating reactivity increased in the order: HOBr ≪ Br(2)O < BrOCl ≈ Br(2) < BrCl. Our results indicate that species other than HOBr can influence bromination rates under conditions typical of drinking water and wastewater chlorination.

Binding in Radical-Solvent Binary Complexes: Benchmark Energies and Performance of Approximate Methods.

In many situations, weak interactions between radicals and their environment potentially influence their properties and reactivity. We computed benchmark binding energies of 12 binary complexes involving radicals, using basis set extrapolated coupled cluster theory with up to CCSDT(Q) excitations plus corrections for core correlation and relativistic effects. The set was comprised of both electron-rich and electron-poor small radicals which were either neutral or positively charged. The radicals were complexed with the closed-shell polar (model) solvent molecules H2O and HF. On the basis of these accurate ab initio binding energies, we assess the performance of many modern DFT functionals for these radical-solvent molecule interactions. Radical hydrogen bonded complexes are well-described by most DFT methods, but two-center-three-electron interactions are at least slightly overbound by most functionals evaluated here, including range-separated functionals. No such systematic error was found for electron-rich metal-water complexes. None of the functionals tested yield chemical accuracy for all types of complexes.

Geometries and Vibrational Frequencies of Small Radicals: Performance of Coupled Cluster and More Approximate Methods.

We generated a new set of reference geometries of small radicals using experimental equilibrium structures, as well as a benchmark-quality coupled cluster additivity scheme including up to quadruples excitations (CCSDTQ). Using these geometries and a set of experimental vibrational frequencies of open shell diatomics, we evaluated the performance of various coupled cluster methods based mainly on unrestricted references, using Dunning basis sets both with and without core correlation. Contrary to previous results, we found that UCCSD(T) and ROCCSD(T) perform equally well for geometries, better than CCSD, and close to their performances for closed shell systems. No improvement over CCSD(T) was achieved by using a Brueckner reference (BD(T)) or full triples (CCSDT). For frequencies, ROCCSD(T), BD(T), and CCSDT improve upon UCCSD(T), especially for the troublesome NO and CO(+) cases. EOMIP-CCSD yields geometries and harmonic frequencies similar to CCSD, and qualitatively correct anharmonic (VPT2) contributions in all cases, like the RO-CC methods. The double hybrid DFT functional B2PLYP-D yields geometries and frequencies of similar quality to that of CCSD but at a much reduced cost. The meta hybrid functionals M06-2X, M06-HF, and BMK perform worse than CCSD, and worse than B3LYP, on average.

Carbon, hydrogen, and nitrogen isotope fractionation during light-induced transformations of atrazine.

The 13C, 2H, and 15N fractionation associated with light-induced transformations of N-containing pesticides in surface waters was investigated using atrazine as a model compound. In laboratory model systems, bulk isotope enrichment factors epsilonC, epsilonH, and epsilonN were determined during the photooxidation of atrazine by excited triplet states of 4-carboxybenzophenone ((3)4-CBBP*), by OH radicals, and during direct photolysis at 254 nm. Moderately large 2H fractionations, quantified by EH values of -51.2 +/- 2.5% per hundred and -25.3 +/- 1.7% per hundred, were found for the transformation of atrazine by (3)4-CBBP* and OH radicals, respectively. 13C and 15N enrichment factors were rather small (-0.3% per hundred > epsilon(C, N) > -1.7% per hundred). The combined delta(13)C, delta(2)H, and delta(15)N analysis suggests that isotope effects are most likely due to H abstraction at the N-H and C-H bonds of the N-alkyl side chains. Direct photolysis of atrazine yielding hydroxyatrazine as main product was characterized by inverse 13C and 15N fractionation (epsilonC = 4.6 +/- 0.3% per hundred, epsilonN = 4.9 +/- 0.2% per hundred) and no detectable 2H fractionation. We hypothesize that isotope effects from photophysical processes involving the excited states of atrazine as well as magnetic isotope effect originating from the magnetic interactions of spin-carrying C and N nuclei have contributed to the observed inverse fractionation. Our study illustrates how compound-specific isotope analysis can be used to differentiate between important direct and indirect phototransformation pathways of agrochemicals in the environment.