Potential Application of Machine-Learning-Based Quantum Chemical Methods in Environmental Chemistry

It is an important topic in environmental sciences to understand the behavior and toxicology of chemical pollutants. Quantum chemical methodologies have served as useful tools for probing behavior and toxicology of chemical pollutants in recent decades. In recent years, machine learning (ML) techniques have brought revolutionary developments to the field of quantum chemistry, which may be beneficial for investigating environmental behavior and toxicology of chemical pollutants. However, the ML-based quantum chemical methods (ML-QCMs) have only scarcely been used in environmental chemical studies so far. To promote applications of the promising methods, this Perspective summarizes recent progress in the ML-QCMs and focuses on their potential applications in environmental chemical studies that could hardly be achieved by the conventional quantum chemical methods. Potential applications and challenges of the ML-QCMs in predicting degradation networks of chemical pollutants, searching global minima for atmospheric nanoclusters, discovering heterogeneous or photochemical transformation pathways of pollutants, as well as predicting environmentally relevant end points with wave functions as descriptors are introduced and discussed.

Atmospheric Chemistry of Allylic Radicals from Isoprene: A Successive Cyclization-Driven Autoxidation Mechanism

The atmospheric chemistry of isoprene has broad implications for regional air quality and the global climate. Allylic radicals, taking 13-17% yield in the isoprene oxidation by ^•Cl, can contribute as much as 3.6-4.9% to all possible formed intermediates in local regions at daytime. Considering the large quantity of isoprene emission, the chemistry of the allylic radicals is therefore highly desirable. Here, we investigated the atmospheric oxidation mechanism of the allylic radicals using quantum chemical calculations and kinetics modeling. The results indicate that the allylic radicals can barrierlessly combine with O₂ to form peroxy radicals (RO₂^•). Under ≤100 ppt NO and ≤50 ppt HO₂^• conditions, the formed RO₂^• mainly undergo two times "successive cyclization and O₂ addition" to finally form the product fragments 2-alkoxy-acetaldehyde (C₂H₃O₂^•) and 3-hydroperoxy-2-oxopropanal (C₃H₄O₄). The presented reaction illustrates a novel successive cyclization-driven autoxidation mechanism. The formed 3-hydroperoxy-2-oxopropanal product is a new isomer of the atmospheric C₃H₄O₄ family and a potential aqueous-phase secondary organic aerosol precursor. Under >100 ppt NO condition, NO can mediate the cyclization-driven autoxidation process to form C₅H₇NO₃, C₅H₇NO₇, and alkoxy radical-related products. The proposed novel autoxidation mechanism advances our current understanding of the atmospheric chemistry of both isoprene and RO₂^•.

Formation of Low-Volatile Products and Unexpected High Formaldehyde Yield from the Atmospheric Oxidation of Methylsiloxanes

With stricter regulation of atmospheric volatile organic compounds (VOCs) originating from fossil fuel-based vehicles and industries, the use of volatile chemical products (VCPs) and the transformation mechanism of VCPs have become increasingly important to quantify air quality. Volatile methylsiloxanes (VMS) are an important class of VCPs and high-production chemicals. Using quantum chemical calculations and kinetics modeling, we investigated the reaction mechanism of peroxy radicals of VMS, which are key intermediates in determining the atmospheric chemistry of VMS. L2-RSiCH₂O₂^• and D3-RSiCH₂O₂^• derived from hexamethyldisiloxane and hexamethylcyclotrisiloxane, respectively, were selected as representative model systems. The results indicated that L2-RSiCH₂O₂^• and D3-RSiCH₂O₂^• follow a novel Si-C-O rearrangement-driven autoxidation mechanism, leading to the formation of low volatile silanols and high yield of formaldehyde at low NO/HO₂^• conditions. At high NO/HO₂^• conditions, L2-RSiCH₂O₂^• and D3-RSiCH₂O₂^• react with NO/HO₂^• to form organic nitrate, hydroperoxide, and active alkoxy radicals. The alkoxy radicals further follow a Si-C-O rearrangement step to finally form formate esters. The novel Si-C-O rearrangement mechanism of both peroxy and alkoxy radicals are supported by available experimental studies on the oxidation of VMS. Notably, the high yield of formaldehyde is estimated to significantly contribute to formaldehyde pollution in the indoor environment, especially during indoor cleaning.

Mechanism and predictive model development of reaction rate constants for N-center radicals with O2

Atmospheric oxidation of NH_x-containing (x = 1, 2) compounds can produce N-center radicals, a precursor of toxic nitrosamines. The reaction rate constant (k_O2) with O₂ has been considered as an important parameter to determine the nitrosamines yield in the subsequent reactions of N-center radicals. However, available k_O2 values of N-center radicals are limited. Here, a three-step scheme including mechanistic analysis and kinetics calculation of the reactions of 28 various N-center radicals with O₂, and model development was taken to solve the k_O2 data shortage. Mainly employed tools include highly cost-expensive coupled-cluster theory (CCSD(T)), kinetic model and statistics. The results indicate that the direct H-abstraction pathway is the most favorable for the reactions of all considered N-center radicals with O₂. The specific molecular conformation and the C-H bond energy of the N-center radicals are two important factors to determine k_O2 values. Based on the mechanistic understanding of k_O2 values, a quantitative structure-activity relationship (QSAR) model of k_O2 values was developed. The model has satisfactory goodness-of-fit, robustness and predictive ability. The determined k_O2 values and the in silico methods provide a scientific base for assessing formation risk of toxic nitrosamines in the atmosphere.

Mechanism and kinetics of the oxidation of dimethyl carbonate by hydroxyl radical in the atmosphere

The mechanism and kinetics for the reaction of dimethyl carbonate (DMC) with OH radical have been studied by using quantum chemical methods. Four reaction pathways were identified for the initial reaction. In the first two pathways, hydrogen atom abstraction is taking place and alkyl radical intermediate is formed with the energy barrier of 6.4 and 7.9 kcal/mol. In the third pathway, OH addition reaction to the carbonyl carbon (C2) atom of DMC and intermediate, I2, is formed with an energy barrier of 11.9 kcal/mol. In the fourth pathway, along with CH₃O^●, methyl hydrogen carbonate is formed. For this C-O bond breaking and O-H addition reaction, the energy barrier is 27 kcal/mol. The calculated enthalpy and Gibbs energy values show that the studied initial reactions are exothermic and exoergic except the OH addition reaction. For the initial reactions, the rate constants were calculated by using canonical variational transition state theory (CVT) with small curvature tunneling (SCT) correction over the temperature range of 278-1200 K. At 298 K, the calculated rate coefficient for the in-plane and out-of-plane hydrogen atom abstraction reaction pathway is 2.30 × 10^-13 and 0.02 × 10^-13 cm³ molecule^-1 s^-1. Further, the reaction between alkyl radical intermediate formed from the first pathway and O₂ is studied. The reaction of alkyl peroxy radical intermediate with atmospheric oxidants, HO₂, NO, and NO₂ is also studied. It was found that the formic (methyl carbonic) anhydride is the end product formed from the atmospheric oxidation and secondary reactions of DMC.

Phototransformation of 2,3-Dibromopropyl-2,4,6-tribromophenyl ether (DPTE) in Natural Waters: Important Roles of Dissolved Organic Matter and Chloride Ion

Novel brominated flame retardants (NBFRs) have become ubiquitous emerging organic pollutants. However, little is known about their transformation in natural waters. In this study, aquatic photochemical behavior of a representative NBFR, 2,3-dibromopropyl-2,4,6-tribromophenyl ether (DPTE), was investigated by simulated sunlight irradiation experiment. Results show that DPTE can undergo direct photolysis (apparent quantum yield 0.008 ± 0.001) and hydroxyl radical (·OH) initiated oxidation (second order reaction rate constant 2.4 × 10⁹ M^-1·s^-1). Dissolved organic matter (DOM) promotes the photodegradation due to generation of excited triplet DOM and ·OH. Two chlorinated intermediates were identified in the photodegradation of DPTE in seawaters. Density functional theory calculation showed that ·Cl or ·Cl₂^- addition reactions on C-Br sites of the phenyl group and H-abstraction reactions from the propyl group are main reaction pathways of DPTE with the chlorine radicals. The ·Cl or ·Cl₂^- addition proceeds via a replacement mechanism to form chlorinated intermediates. Environmental half-lives of DPTE relevant with photodegradation are estimated to be 6.5-1153.9 days in waters of the Yellow River estuarine region. This study provides valuable insights into the phototransformation behavior of DPTE in natural waters, which is helpful for persistence assessment of the NBFRs.

Atmospheric Oxidation of Piperazine Initiated by ·Cl: Unexpected High Nitrosamine Yield

Chlorine radicals (·Cl) initiated amine oxidation plays an important role for the formation of carcinogenic nitrosamine in the atmosphere. Piperazine (PZ) is considered as a potential atmospheric pollutant since it is an alternative solvent to monoethanolamine (MEA), a benchmark solvent in a leading CO₂ capture technology. Here, we employed quantum chemical methods and kinetics modeling to investigate ·Cl-initiated atmospheric oxidation of PZ, particularly concerning the potential of PZ to form nitrosamine compared to MEA. Results showed that the ·Cl-initiated PZ reaction exclusively leads to N-center radicals (PZ-N) that mainly react with NO to produce nitrosamine in their further reaction with O₂/NO. Together with the PZ + ·OH reaction, the PZ-N yield from PZ oxidation is still lower than that of the corresponding MEA reactions. However, the nitrosamine yield of PZ is higher than the reported value for MEA when [NO] is <5 ppb, a concentration commonly encountered in a polluted urban atmosphere. The unexpected high nitrosamine yield from PZ compared to MEA results from a more favorable reaction of N-center radicals with NO compared to O₂. These findings show that the yield of N-center radicals cannot directly be used as a metric for the yield of the corresponding carcinogenic nitrosamine.