Dra. Durán-Guajardo, Rocío

In this research, we investigated the essential role of biogenic volatile organic compound emissions in regulating tropospheric ozone levels, atmospheric chemistry, and climate dynamics. We explored linalool ozonolysis and secondary organic aerosol formation mechanisms, providing key insights into atmospheric processes. Computational techniques, such as density functional theory calculations and molecular dynamics simulations, were employed for the analysis. Our study delves into the energetic and mechanistic aspects of the 1,3-dipolar cycloadditions involving linalool and its ozonolysis byproducts, known as Criegee intermediates. A total of 24 reactions were analyzed from the three possible Criegee intermediates formed, resulting from different reactant orientations and their endo/exo isomers. We found that only four of these reactions exhibit large rate constants that can compete with tropospheric reactions. This reactivity pattern was characterized by analyzing reactivity indices from conceptual density functional theory and determining that electron flux originates from linalool to the Criegee intermediates. Greater electrophilicity in the Criegee intermediates results in a lower reaction activation energy, confirmed by the global electrophilicity index. Furthermore, using the activation strain model and energy decomposition analysis, we found that differences in activation energies were primarily driven by nonorbital energy factors. Finally, molecular dynamics simulations showed that the final cycloaddition adducts of the most favorable 1,3-dipolar cycloaddition interact favorably with water molecules in an exergonic process, adsorbing up to 92% of the water molecules after 20 ns. Our findings provide insights that enhance our understanding of the interactions between natural emissions and atmospheric constituents.

Context: This study investigates the energetic and polarizability characteristics of three 1,3-dipolar cycloaddition reactions between diazene oxide and substituted ethylenes, focusing on the transition from synchronous to asynchronous mechanisms. Synchronicity analysis, using the reaction force constant, indicates that the bond evolution process becomes increasingly decoupled as the number of cyano groups increases. Polarizability analysis reveals that isotropic polarizability reaches its maximum near the transition state in all cases, while anisotropy of polarizability shifts from the transition state toward the product direction as asynchronicity increases. The larger the shift, the more asynchronous the mechanism, as refected by the weight of the transition region. A detailed examination of the parallel and perpendicular polarizability components to the newly formed sigma bonds shows that the evolution of the parallel component is closely aligned with the energetic changes along the reaction coordinate, particularly in the synchronous reaction. We have also identifed a relationship between the displacement in the maximum state of the parallel component from the transition state and the synchronicity of the mechanism. The larger the displacement, the more asynchronous the mechanism. These fndings suggest that asynchronous 1,3-dipolar cycloaddition mechanisms are characterized by a decoupling of isotropic and anisotropic polarizabilities and a shift in the maximum polarizability state of the parallel component toward the product direction. Methods: Density functional theory calculations were performed at the B3LYP/6–311+ +G(d,p)//B3LYP/6-31G(d,p) level of theory. The polarizability was calculated at each point of the reaction path, obtained using the intrinsic reaction coordinate method, as implemented in Gaussian 16.

Context: The negative of the Shannon entropy derivative is proposed to account for electron density contraction as the chemical bonds are breaking and forming during a chemical reaction. We called this property the electron density contraction index, EDC, which allows identifying stages in a reaction that are dominated by electron contraction or expansion. Four diferent reactions were analyzed to show how the EDC index changes along the reaction coordinate. The results indicate that the rate of change of Shannon entropy is directly related to the rate of change of the electron density at the bond critical points between all the atomic pairs in the molecular systems. It is expected that EDC will complement the detailed analysis of reaction mechanisms that can be performed with the theoretical tools available to date. Methods: Density functional theory calculations at the B3LYP/6-31G(d,p) level of theory were carried out using Gaussian 16 to analyze the reaction mechanisms of the four reactions studied. The reaction paths were obtained via the intrinsic reaction coordinate method, which served as the reaction coordinate to obtain the reaction force and the EDC profles in each case. Shannon entropy and electron density at the bond critical points were calculated using the Multiwfn 3.7 package.