Directional Proton Transfer in the Reaction of the Simplest Criegee Intermediate with Water Involving the Formation of Transient H3O. J Phys Chem Lett 2021;12:3379-3386. [PMID: 33784110 DOI: 10.1021/acs.jpclett.1c00448]

Barrierless reactions of C2 Criegee intermediates with H2SO4 and their implication to oligomers and new particle formation

The formation of oligomeric hydrogen peroxide triggered by Criegee intermediate maybe contributes significantly to the formation and growth of secondary organic aerosol (SOA). However, to date, the reactivity of C2 Criegee intermediates (CH3CHOO) in areas contaminated with acidic gas remains poorly understood. Herein, high-level quantum chemical calculations and Born-Oppenheimer molecular dynamics (BOMD) simulations are used to explore the reaction of CH3CHOO and H2SO4 both in the gas phase and at the air-water interface. In the gas phase, the addition reaction of CH3CHOO with H2SO4 to generate CH3HC(OOH)OSO3H (HPES) is near-barrierless, regardless of the presence of water molecules. BOMD simulations show that the reaction at the air-water interface is even faster than that in the gas phase. Further calculations reveal that the HPES has a tendency to aggregate with sulfuric acids, ammonias, and water molecules to form stable clusters, meanwhile the oligomerization reaction of CH3CHOO with HPES in the gas phase is both thermochemically and kinetically favored. Also, it is noted that the interfacial HPES- ion can attract H2SO4, NH3, (COOH)2 and HNO3 for particle formation from the gas phase to the water surface. Thus, the results of this work not only elucidate the high atmospheric reactivity of C2 Criegee intermediates in polluted regions, but also deepen our understanding of the formation process of atmospheric SOA induced by Criegee intermediates.

Ab initio molecular dynamics of solvation effects and reactivity at the interface between water and ascorbic acid covered anatase TiO2 (101)

In this work, we present a detailed study of the interaction between ascorbic acid (L-asc) and anatase TiO2 (101) surface both in gas phase and in contact with water by using density functional theory and ab initio molecular dynamics simulations. In gas phase, L-asc strongly binds the TiO2 (101) surface as a dianion (L-asc2- ), adopting a bridging bidentate coordination mode (BB), with the two acid protons transferred to two surface 2-fold bridging oxygens (O2c). AIMD simulations show that the interaction between the organic ligand and the anatase surface is stable and comparable to the vacuum one despite the possible solvent effects and/or possible structural distortions of the ligand. In addition, during the AIMD simulations hydroxylation phenomena occur forming transient H3 O+ ions at the solid-liquid interface. For the first time, our results provide insight into the role of the ascorbic acid on the electronic properties of the TiO2 (101), the influence of the water environment on the ligand-surface interaction and the nature of the solid-liquid interface.

Machine-learning assisted scheduling optimization and its application in quantum chemical calculations

Easy and effective usage of computational resources is crucial for scientific calculations, both from the perspectives of timeliness and economic efficiency. This work proposes a bi-level optimization framework to optimize the computational sequences. Machine-learning (ML) assisted static load-balancing, and different dynamic load-balancing algorithms can be integrated. Consequently, the computational and scheduling engine of the ParaEngine is developed to invoke optimized quantum chemical (QC) calculations. Illustrated benchmark calculations include high-throughput drug suit, solvent model, P38 protein, and SARS-CoV-2 systems. The results show that the usage rate of given computational resources for high throughput and large-scale fragmentation QC calculations can primarily profit, and faster accomplishing computational tasks can be expected when employing high-performance computing (HPC) clusters.

Multiple evaluations of atmospheric behavior between Criegee intermediates and HCHO: Gas-phase and air-water interface reaction

Given the high abundance of water in the atmosphere, the reaction of Criegee intermediates (CIs) with (H₂O)₂ is considered to be the predominant removal pathway for CIs. However, recent experimental findings reported that the reactions of CIs with organic acids and carbonyls are faster than expected. At the same time, the interface behavior between CIs and carbonyls has not been reported so far. Here, the gas-phase and air-water interface behavior between Criegee intermediates and HCHO were explored by adopting high-level quantum chemical calculations and Born-Oppenheimer molecular dynamics (BOMD) simulations. Quantum chemical calculations evidence that the gas-phase reactions of CIs + HCHO are submerged energy or low energy barriers processes. The rate ratios speculate that the HCHO could be not only a significant tropospheric scavenger of CIs, but also an inhibitor in the oxidizing ability of CIs on SO_x in dry and highly polluted areas with abundant HCHO concentration. The reactions of CH₂OO with HCHO at the droplet's surface follow a loop structure mechanism to produce i) SOZ (), ii) BHMP (HOCH₂OOCH₂OH), and iii) HMHP (HOCH₂OOH). Considering the harsh reaction conditions between CIs and HCHO at the interface (i.e., the two molecules must be sufficiently close to each other), the hydration of CIs is still their main atmospheric loss pathway. These results could help us get a better interpretation of the underlying CIs-aldehydes chemical processes in the global polluted urban atmospheres.

Water Charge Transfer Accelerates Criegee Intermediate Reaction with H2O- Radical Anion at the Aqueous Interface

Criegee intermediates (CIs) are important carbonyl oxides that may react with atmospheric trace chemicals and impact the global climate. The CI reaction with water has been widely studied and is a main channel for trapping CIs in the troposphere. Previous experimental and computational reports have largely focused on reaction kinetic processes in various CI-water reactions. The molecular-level origin of CI's interfacial reactivity at the water microdroplet surface (e.g., as found in aerosols and clouds) is unclear. In this study, by employing the quantum mechanical/molecular mechanical (QM/MM) Born-Oppenheimer molecular dynamics with the local second-order Møller-Plesset perturbation theory, our computational results reveal a substantial water charge transfer up to ∼20% per water, which creates the surface H₂O⁺/H₂O^- radical pairs to enhance the CH₂OO and anti-CH₃CHOO reactivity with water: the resulting strong CI-H₂O^- electrostatic attraction at the microdroplet surface facilitates the nucleophilic attack to the CI carbonyl by water, which may counteract the apolar hindrance of the substituent to accelerate the CI-water reaction. Our statistical analysis of the molecular dynamics trajectories further resolves a relatively long-lived bound CI(H₂O^-) intermediate state at the air/water interface, which has not been observed in gaseous CI reactions. This work provides insights into what may alter the oxidizing power of the troposphere by the next larger CIs than simple CH₂OO and implicates a new perspective on the role of interfacial water charge transfer in accelerating molecular reactions at aqueous interfaces.

Theoretical Study on the Gas Phase and Gas-Liquid Interface Reaction Mechanism of Criegee Intermediates with Glycolic Acid Sulfate

Criegee intermediates (CIs) are important zwitterionic oxidants in the atmosphere, which affect the budget of OH radicals, amines, alcohols, organic/inorganic acids, etc. In this study, quantum chemical calculation and Born-Oppenheimer molecular dynamic (BOMD) simulation were performed to show the reaction mechanisms of C2 CIs with glycolic acid sulfate (GAS) at the gas-phase and gas-liquid interface, respectively. The results indicate that CIs can react with COOH and OSO₃H groups of GAS and generate hydroperoxide products. Intramolecular proton transfer reactions occurred in the simulations. Moreover, GAS acts as a proton donor and participates in the hydration of CIs, during which the intramolecular proton transfer also occurs. As GAS widely exists in atmospheric particulate matter, the reaction with GAS is one of the sink pathways of CIs in areas polluted by particulate matter.

Organic Peroxides in Aerosol: Key Reactive Intermediates for Multiphase Processes in the Atmosphere

Organic peroxides (POs) are organic molecules with one or more peroxide (-O-O-) functional groups. POs are commonly regarded as chemically labile termination products from gas-phase radical chemistry and therefore serve as temporary reservoirs for oxidative radicals (HOx and ROx) in the atmosphere. Owing to their ubiquity, active gas-particle partitioning behavior, and reactivity, POs are key reactive intermediates in atmospheric multiphase processes determining the life cycle (formation, growth, and aging), climate, and health impacts of aerosol. However, there remain substantial gaps in the origin, molecular diversity, and fate of POs due to their complex nature and dynamic behavior. Here, we summarize the current understanding on atmospheric POs, with a focus on their identification and quantification, state-of-the-art analytical developments, molecular-level formation mechanisms, multiphase chemical transformation pathways, as well as environmental and health impacts. We find that interactions with SO2 and transition metal ions are generally the fast PO transformation pathways in atmospheric liquid water, with lifetimes estimated to be minutes to hours, while hydrolysis is particularly important for α-substituted hydroperoxides. Meanwhile, photolysis and thermolysis are likely minor sinks for POs. These multiphase PO transformation pathways are distinctly different from their gas-phase fates, such as photolysis and reaction with OH radicals, which highlights the need to understand the multiphase partitioning of POs. By summarizing the current advances and remaining challenges for the investigation of POs, we propose future research priorities regarding their origin, fate, and impacts in the atmosphere.

Toward High-level Machine Learning Potential for Water Based on Quantum Fragmentation and Neural Networks

Accurate and efficient simulation of liquids, such as water and salt solutions, using high-level wave function theories is still a formidable task for computational chemists owing to the high computational costs. In this study, we develop a deep machine learning potential based on fragment-based second-order Møller-Plesset perturbation theory (DP-MP2) for water through neural networks. We show that the DP-MP2 potential predicts the structural, dynamical, and thermodynamic properties of liquid water in better agreement with the experimental data than previous studies based on density functional theory (DFT). The nuclear quantum effects (NQEs) on the properties of liquid water are also examined, which are noticeable in affecting the structural and dynamical properties of liquid water under ambient conditions. This work provides a general framework for quantitative predictions of the properties of condensed-phase systems with the accuracy of high-level wave function theory while achieving significant computational savings compared to ab initio simulations.

Ab initio multi-level layered elongation method and its application to local interaction analysis between DNA bulge and ligand molecules

A multi-level layered elongation method was developed for efficiently analyzing the electronic states of local structures in large bio/nano-systems at the full ab initio level of theory. The original elongation method developed during the last three decades in our group has focused on the system in one direction from one terminal to the other terminal to sequentially construct the electronic states of a polymer, called a theoretical synthesis of polymers. In this study, an important region termed the central (C) part is targeted in a large polymer and the remainder are terminal (T) parts. The electronic structures along with polymer elongation are calculated repeatedly from both end T parts to the C central part at the same time. The important C part is treated with large basis sets (high level) and the other regions are treated with small basis sets (low level) in the ab initio theoretical framework. The electronic structures besides the C part can be reused for other systems with different structures at the C part, which renders the method computationally efficient. This multi-level layered elongation method was applied to the investigation on DNA single bulge recognition of small molecules (ligands). The reliability and validity of our approach were examined in comparison with the results obtained by direct calculations using a conventional quantum chemical method for the entire system. Furthermore, stabilization energies by the formation of the complex of bulge DNA and a ligand were estimated with basis set superposition error corrections incorporated into the elongation method.

Hydrogen Bonds and H3O+ Formation at the Water Interface with Formic Acid Covered Anatase TiO2

Carboxylic acid-modified TiO₂ surfaces in aqueous environment are of widespread interest, yet atomic-scale understanding of their structure is limited. We here investigate formic acid (FA) on anatase TiO₂ (101) (A-101) in contact with water using density functional theory (DFT) and ab initio molecular dynamics (AIMD). Isolated FA molecules adsorbed in a deprotonated bridging bidentate (BD) form on A-101 are found to remain stable at the interface with water, with the acid proton transferred to a surface oxygen to form a surface bridging hydroxyl (O_brH). With increasing FA coverage, adsorbed monolayers of only BD and successively of alternating monodentate (MD) and BD species give rise to a higher concentration of surface O_brH's. Simulations of these adsorbed monolayers in water environment show that some protons are released from the surface O_brH's to water resulting in a negatively charged surface with nearby solvated H₃O⁺ ions. These results provide insight into the complex acid-base equilibrium between an oxide surface, adsorbates and water and can also help obtain a better understanding of the wetting properties of chemically modified TiO₂ surfaces.