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.

Photosensitization mechanisms at the air-water interface of aqueous aerosols

Photosensitization reactions are believed to provide a key contribution to the overall oxidation chemistry of the Earth's atmosphere. Generally, these processes take place on the surface of aqueous aerosols, where organic surfactants accumulate and react, either directly or indirectly, with the activated photosensitizer. However, the mechanisms involved in these important interfacial phenomena are still poorly known. This work sheds light on the reaction mechanisms of the photosensitizer imidazole-2-carboxaldehyde through ab initio (QM/MM) molecular dynamics simulations and high-level ab initio calculations. The nature of the lowest excited states of the system (singlets and triplets) is described in detail for the first time in the gas phase, in bulk water, and at the air–water interface, and possible intersystem crossing mechanisms leading to the reactive triplet state are analyzed. Moreover, the reactive triplet state is shown to be unstable at the air–water surface in a pure water aerosol. The combination of this finding with the results obtained for simple surfactant-photosensitizer models, together with experimental data from the literature, suggests that photosensitization reactions assisted by imidazole-2-carboxaldehyde at the surface of aqueous droplets can only occur in the presence of surfactant species, such as fatty acids, that stabilize the photoactivated triplet at the interface. These findings should help the interpretation of field measurements and the design of new laboratory experiments to better understand atmospheric photosensitization processes.

First-principles molecular dynamics simulations of imidazole-2-carboxaldehyde at the air–water interface highlight the role of surfactants in stabilising the reactive triplet state involved in photosensitisation reactions in aqueous aerosols.

Multiscale Modeling Strategy of 2D Covalent Organic Frameworks Confined at an Air-Water Interface

Two-dimensional covalent organic frameworks (2D COFs) have attracted attention as versatile active materials in many applications. Recent advances have demonstrated the synthesis of monolayer 2D COF via an air-water interface. However, the interfacial 2D polymerization mechanism has been elusive. In this work, we have used a multiscale modeling strategy to study dimethylmethylene-bridged triphenylamine building blocks confined at the air-water interface to form a 2D COF via Schiff-base reaction. A synergy between the computational investigations and experiments allowed the synthesis of a 2D-COF with one of the linkers considered. Our simulations complement the experimental characterization and show the preference of the building blocks to be at the interface with a favorable orientation for the polymerization. The air-water interface is shown to be a key factor to stabilize a flat conformation when a dimer molecule is considered. The structural and electronic properties of the monolayer COFs based on the two monomers are calculated and show a semiconducting nature with direct bandgaps. Our strategy provides a first step toward the in silico polymerization of 2D COFs at air-water interfaces capturing the initial steps of the synthesis up to the prediction of electronic properties of the 2D material.

Multiconfigurational Quantum Chemistry Determinations of Absorption Cross Sections (σ) in the Gas Phase and Molar Extinction Coefficients (ε) in Aqueous Solution and Air-Water Interface

Theoretical determinations of absorption cross sections (σ) in the gas phase and molar extinction coefficients (ε) in condensed phases (water solution, interfaces or surfaces, protein or nucleic acids embeddings, etc.) are of interest when rates of photochemical processes, J = ∫ ϕ(λ) σ(λ) I(λ) dλ, are needed, where ϕ(λ) and I(λ) are the quantum yield of the process and the irradiance of the light source, respectively, as functions of the wavelength λ. Efficient computational strategies based on single-reference quantum-chemistry methods have been developed enabling determinations of line shapes or, in some cases, achieving rovibrational resolution. Developments are however lacking for strongly correlated problems, with many excited states, high-order excitations, and/or near degeneracies between states of the same and different spin multiplicities. In this work, we define and compare the performance of distinct computational strategies using multiconfigurational quantum chemistry, nuclear sampling of the chromophore (by means of molecular dynamics, ab initio molecular dynamics, or Wigner sampling), and conformational and statistical sampling of the environment (by means of molecular dynamics). A new mathematical approach revisiting previous absolute orientation algorithms is also developed to improve alignments of geometries. These approaches are benchmarked through the nπ* band of acrolein not only in the gas phase and water solution but also in a gas-phase/water interface, a common situation for instance in atmospheric chemistry. Subsequently, the best strategy is used to compute the absorption band for the adduct formed upon addition of an OH radical to the C6 position of uracil and compared with the available experimental data. Overall, quantum Wigner sampling of the chromophore with molecular dynamics sampling of the environment with CASPT2 electronic-structure determinations arise as a powerful methodology to predict meaningful σ(λ) and ε(λ) band line shapes with accurate absolute intensities.

Reactivity of Undissociated Molecular Nitric Acid at the Air-Water Interface

Recent experiments and theoretical calculations have shown that HNO₃ may exist in molecular form in aqueous environments, where in principle one would expect this strong acid to be completely dissociated. Much effort has been devoted to understanding this fact, which has huge environmental relevance since nitric acid is a component of acid rain and also contributes to renoxification processes in the atmosphere. Although the importance of heterogeneous processes such as oxidation and photolysis have been evidenced by experiments, most theoretical studies on hydrated molecular HNO₃ have focused on the acid dissociation mechanism. In the present work, we carry out calculations at various levels of theory to obtain insight into the properties of molecular nitric acid at the surface of liquid water (the air-water interface). Through multi-nanosecond combined quantum-classical molecular dynamics simulations, we analyze the interface affinity of nitric acid and provide an order of magnitude for its lifetime with regard to acid dissociation, which is close to the value deduced using thermodynamic data in the literature (∼0.3 ns). Moreover, we study the electronic absorption spectrum and calculate the rate constant for the photolytic process HNO₃ + hν → NO₂ + OH, leading to 2 × 10^-6 s^-1, about twice the value in the gas phase. Finally, we describe the reaction HNO₃ + OH → NO₃ + H₂O using a cluster model containing 21 water molecules with the help of high-level ab initio calculations. A large number of reaction paths are explored, and our study leads to the conclusion that the most favorable mechanism involves the formation of a pre-reactive complex (HNO₃)(OH) from which product are obtained through a coupled proton-electron transfer mechanism that has a free-energy barrier of 6.65 kcal·mol^-1. Kinetic calculations predict a rate constant increase by ∼4 orders of magnitude relative to the gas phase, and we conclude that at the air-water interface, a lower limit for the rate constant is k = 1.2 × 10^-9 cm³·molecule^-1·s^-1. The atmospheric significance of all these results is discussed.

Molecular Mechanisms and Atmospheric Implications of Criegee Intermediate-Alcohol Chemistry in the Gas Phase and Aqueous Surface Environments

Criegee intermediates and alcohols are important species in the atmosphere. In this study, we use quantum chemistry and Born-Oppenheimer molecular dynamics (BOMD) simulations to investigate the reaction between methanol/ethanol and Criegee intermediates (anti- or syn-CH₃CHOO) in the gas phase and at the air-water interface. Reactions at the interface are found to be much faster than those in the gas phase. When water molecules are available, loop structures can be formed to facilitate the reaction. In addition, nonloop reaction pathways characterized by the formation of hydrated protons, although with a low possibility, are also identified at the air-water interface. Implications of our results on the fate of Criegee intermediates in the atmosphere are discussed, which deepen our understanding of Criegee intermediate-alcohol chemistry in humid environments.

Reaching multi-nanosecond timescales in combined QM/MM molecular dynamics simulations through parallel horsetail sampling

We report an enhanced sampling technique that allows to reach the multi-nanosecond timescale in quantum mechanics/molecular mechanics molecular dynamics simulations. The proposed technique, called horsetail sampling, is a specific type of multiple molecular dynamics approach exhibiting high parallel efficiency. It couples a main simulation with a large number of shorter trajectories launched on independent processors at periodic time intervals. The technique is applied to study hydrogen peroxide at the water liquid-vapor interface, a system of considerable atmospheric relevance. A total simulation time of a little more than 6 ns has been attained for a total CPU time of 5.1 years representing only about 20 days of wall-clock time. The discussion of the results highlights the strong influence of the solvation effects at the interface on the structure and the electronic properties of the solute. © 2017 Wiley Periodicals, Inc.

A polarizable continuum model for molecules at spherical diffuse interfaces

We present an extension of the Polarizable Continuum Model (PCM) to simulate solvent effects at diffuse interfaces with spherical symmetry, such as nanodroplets and micelles. We derive the form of the Green's function for a spatially varying dielectric permittivity with spherical symmetry and exploit the integral equation formalism of the PCM for general dielectric environments to recast the solvation problem into a continuum solvation framework. This allows the investigation of the solvation of ions and molecules in nonuniform dielectric environments, such as liquid droplets, micelles or membranes, while maintaining the computationally appealing characteristics of continuum solvation models. We describe in detail our implementation, both for the calculation of the Green's function and for its subsequent use in the PCM electrostatic problem. The model is then applied on a few test systems, mainly to analyze the effect of interface curvature on solvation energetics.