Dra. Durán-Guajardo, Rocío

We present a clear relationship between hydride affinity and the condensed-to-boron electrophilicity index for 27 intramolecular aminoborane frustrated Lewis pairs (FLPs), employing Density Functional Theory calculations. First, the Gibbs free energy change's dependence of H2 splitting on the hydride affinity of the FLPs is explored, identifying strong linear correlations for aromatic and aliphatic nitrogen substituents. Additionally, we found that the condensed-to-boron electrophilicity index's inverse and logarithm show strong linear relationships with hydride affinity across all FLP systems, generating a direct connection between a purely theoretical reactivity index with an experimental property like the hydride affinity. Using these models, we applied our findings to an experimentally tested FLP system, predicting hydride affinity values that yield ΔG (ΔH) for H2 splitting of −0.2 (−8.2) and 0.9 (−9.4) kcal mol−1 with the inverse and logarithmic relationships, respectively, which are close to the DFT calculated values: 1.3 (−6.0) kcal mol−1 and closely matching the experimental value of −1.7 ± 4.0 (−8 ± 1) kcal mol−1. These results confirm the predictive power of our models, which can be used to guide the rational design of new aminoborane FLPs to be employed in reversible H2 activation and metal-free hydrogenation of unsaturated compounds.

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.

The dimerization of intramolecular aminoborane and aminoalane frustrated Lewis pairs was investigated using density functional theory. We systematically varied the substituents to gradually increase their bulkiness, including H, CH3, t-Bu, Ph, and Mes groups. Starting from the most stable conformer of the monomers, a frustrated Lewis pair or classic Lewis adduct, we studied the dimerization process for all systems, revealing significant variations in the Gibbs free energy. Dimerization was favored in four aminoboranes and six aminoalanes, depending on the specific combinations of substituents. Applying an energy decomposition analysis, we found that the preparation energy of the monomers and the non-orbital interactions between them are the primary contributors to the observed energetic differences, showing a clear linear relationship. Additionally, we analyzed the electronic effects by increasing the acidity of the Lewis acid, observing a shift toward endergonic and exergonic directions in aminoboranes and aminoalanes, respectively. This shift was attributed to the stabilization of a classic Lewis adduct. This study underscores three crucial factors influencing dimer formation: (i) substituent size, (ii) stabilization of the classic Lewis adduct conformation, and (iii) covalent radii of the Lewis centers. Understanding these factors is essential for designing FLPs and preventing unwanted dimerization that could affect their catalytic performance in H2 activation processes.

In this article, we present a comprehensive computational investigation into the reaction mechanism of N-arylation of substituted aryl halides through Ullmann-type coupling reactions. Our computational findings, obtained through DFT ωB97X-D/6-311G(d,p) and ωB97X-D/LanL2DZ calculations, reveal a direct relation between the previously reported experimental reaction yields and the activation energy of haloarene activation, which constitutes the rate-limiting step in the overall coupling process. A detailed analysis of the reaction mechanism employing the Activation Strain Model indicates that the strain in the substituted iodoanilines is the primary contributor to the energy barrier, representing an average of 80% of the total strain energy. Additional analysis based on conceptual Density Functional Theory (DFT) suggests that the nucleophilicity of the nitrogen in the lactam is directly linked to the activation energies. These results provide valuable insights into the factors influencing energetic barriers and, consequently, reaction yields. These insights enable the rational modification of reactants to optimize the N-arylation process.

The mechanistic paradigm in which the Schmittel cyclization transitions from one-step to stepwise has been investigated through the stabilization of a full hidden intermediate in the framework of the Diabatic Model of Intermediate Stabilization. Hidden intermediate activation was studied in silico employing quasi-classical trajectories and the Electron Localization Function. The stabilization of hidden intermediates achieved by substituting enyne–allenes with cyano and nitro groups generates the appearance of a partially hidden and an explicit intermediate, leading to one-step asynchronous biradical and stepwise biradical/zwitterionic mechanisms, respectively. The mechanistic feature associated with the activation level of the hidden intermediate arises from the Thornton effect and non-RRKM dynamics, where in the case of the CN-substituted system, despite having a single transition state, 54% of the effective trajectories remain in the intermediate zone after 540 fs, indicating that a mixture of mechanisms is observed.

A deeper computational mechanistic study of an environmentally friendly metal-free CO2 reduction process towards obtaining methanol is presented, employing a previously tested kind of intramolecular frustrated Lewis pair (2-[bis(R)boryl]-N,N-dimethylaniline) as the catalyst and H2 as the reducing agent. The Lewis acid strength of the electrophilic boron atom was adjusted to facilitate hydride release by changing the R group, using electron-donating groups (EDGs) based on methylated aryls (Mes and Mes′) and electron-withdrawing groups (EWGs) based on fluorinated alkyls (CF3 and PFtB) and aryls (FMes and C6F5), to analyse its effect on both the H2 splitting and CO2 hydrogenation processes. The acidity of boron was measured from the local electrophilicity index obtained using conceptual density functional theory, where an excellent correlation with the Gibbs free energy of the H2 splitting process was found (R2 = 0.95), indicating that the higher the acid power of boron is, the more exergonic the H2 activation process is. The reversibility of H2 activation directly impacts the CO2 and formic acid hydrogenations, where the less exergonic the H2 splitting process is, the lower the activation energies for these hydrogenation processes are. To obtain methanol at the end, methanediol dehydration forming formaldehyde is crucial, because methanediol has a high energetic barrier, hindering the catalytic cycle from being more efficient. Conversely, formaldehyde can be easily hydrogenated to methanol in the same way as CO2 and formic acid. Finally, the catalytic activity in each case was analysed in terms of the energetic span model, where the local electrophilicity index condensed on boron shows a good linear correlation with the logarithm of the relative turnover frequency (R2 = 0.91), indicating that this reactivity index can be employed to guide the design of optimal catalytic systems to increase its catalytic activity, opening new routes directing future experiments in the field.

A large set of intramolecular aminoborane-based FLPs was studied employing density functional theory in the H2 activation process to analyze how the acidity and basicity of boron and nitrogen atoms, respectively, affect the reversibility of the process. Three different linkers were employed, keeping the C–C nature in the connection between both Lewis centers: –CH2–CH2–, –CH[double bond, length as m-dash]CH–, and –C6H4–. The results show that significant differences in the Gibbs free energy of the process are found by considering all the combinations of substituents. Of the 75 systems studied, only 9 showed the ability to carry out the process reversibly (ΔGH2 in the range of −3.5 to 2.0 kcal mol−1), where combinations of alkyl/aryl or aryl/alkyl in boron/nitrogen generate systems capable of reaching reversibility. If the alkyl/alkyl or aryl/aryl combination is employed, highly exergonic (non-reversible H2 activation) and endergonic (unfeasible H2 activation) reactions are found, respectively. No appreciable differences in the linker were found, allowing us to continue the analysis with the most entropically favorable linker, the –C6H4– linker. From this, 25 different FLP systems of type 2-[bis(X)boryl]-(Y)aniline (X: H, CF3, C6F5, PFtB, FMes and Y: H, CH3, t-but, Ph, Mes) can be formed. By analyzing the electronic properties of each system, we have found that the condensed-to-boron electrophilicity index ωB+ is inversely related to the ΔGH2. Interestingly, two relationships were found; the first is for alkyl groups (Y: CH3 and t-but) and the second for aryl groups (Y: H, Ph, and Mes), which is intimately related to the proton affinity of each aniline. In addition, it is quite interesting when the frustration degree, given by B⋯N distance dB–N, is brought together with ωB+, since the Image ID:d3sc03992g-t1.gif quotient has unit energy/length corresponding to unit force; concomitantly, a measure of the FLP strength in H–H bond activation can be defined. With this finding, a rational design of this kind of FLP can be performed by analyzing the acidity of boron through condensed-to-boron electrophilicity and knowing the nature of the substituent of nitrogen according to whether the Y is alkyl or aryl, optimizing the H2 reversible activation in a rational way, which is crucial to improve the catalytic performance.

In this paper, we will study the reactivity along with substituent changes in the OH insertion reaction in copper carbenoids. To this end, we have used M06-2X functional with cc-pVDZ for light atoms and LanL2DZ for copper. We have studied the IRC insertion profiles and analysed reactivity indexes such as electrophilicity (ω) and pKa calculations. We have found that with R1 substitutions phenyl group, R2 substitutions amide group lower the reaction barrier considerably. Concerning the substrate reactions, the most favoured substituent is NO2 in para position.