Formulating Reduced Density Gradient Approaches for Noncovalent Interactions.

This work elucidates several forms of reduced electron density gradient (RDG) to describe noncovalent interactions (NCIs). By interpreting the RDG as a local moment function, we systematically leveraged Weizacker's and Fermi's local moments. This resulted in high-fidelity RDG representations consistent with the NCI analysis. In addition, the RDG version derived from the Lagrangian kinetic energy density is conveniently normalized. These results suggest the nonexistence of a particular RDG formulation when performing NCI analysis. Thus, an in-depth examination of the theoretical foundations connecting the RDG function with the nature of noncovalent interactions is necessary.

Revisiting the bonding evolution theory: a fresh perspective on the ammonia pyramidal inversion and bond dissociations in ethane and borazane.

This work offers a comprehensive and fresh perspective on the bonding evolution theory (BET) framework, originally proposed by Silvi and collaborators [X. Krokidis, S. Noury and B. Silvi, Characterization of elementary chemical processes by catastrophe theory, J. Phys. Chem. A, 1997, 101, 7277-7282]. By underscoring Thom's foundational work, we identify the parametric function characterizing bonding events along a reaction pathway through a three-step sequence to establish such association rigorously, namely: (a) computing the determinant of the Hessian matrix at all potentially degenerate critical points, (b) computing the relative distance between these points, and (c) assigning the unfolding based on these computations and considering the maximum number of critical points for each unfolding. In-depth examination of the ammonia inversion and the dissociation of ethane and ammonia borane molecules yields a striking discovery: no elliptic umbilic flag is detected along the reactive coordinate for any of the systems, contradicting previous reports. Our findings indicate that the core mechanisms of these chemical reactions can be understood using only two folds, the simplest polynomial of Thom's theory, leading to considerable simplification. In contrast to previous reports, no signatures of the elliptic umbilic unfolding were detected in any of the systems examined. This finding dramatically simplifies the topological rationalization of electron rearrangements within the BET framework, opening new approaches for investigating complex reactions.

A simple topology-based model for predicting the activation barriers of reactive processes at 0 K.

This work reveals an underlying correlation between the topology and energetic features of matter configurations/rearrangements by exploiting two topological concepts, namely, structural stability and persistency, leading thus to a model capable of predicting activation energies at 0 K. This finding provides some answers to the difficulties of applying Thom's functions for extracting energetic information of rate processes, which has been a limitation for exact, biological, and technological sciences. A linear relationship between the experimental barriers of 17 chemical reactions and both concepts was found by studying these systems' topography along the intrinsic reaction coordinate. Such a procedure led to the model , which accurately predicts the activation energy in reacting systems involving organic and organometallic compounds under different conditions, e.g., the gas-phase, solvent media, and temperature. This function was further recalibrated to enhance its predicting capabilities, generating the equation for this procedure, characterized by a squared Pearson correlation coefficient (r2 = 0.9774) 1.1 times higher. Surprisingly, no improvement was observed.

Elucidating the N-N and C-N Bond-breaking Mechanism in the Photoinduced Formation of Nitrile Imine.

In this study, we revealed the significance of chemical bonding for the photochemically induced mechanism of 2-phenyl tetrazole derivatives generating nitrile imines. The correlated electron localization function shows that the formation of imine nitrile involves two key bond events: (i) the heterolytic C-N breakage taking place in the T1 state and (ii) the homolytic N-N rupture occurring in the T2 excited state. In particular, a cation-radical specie results from the C-N cleavage, whereas the N-N rupture creates a biradical resonant form of imine nitrile. Additionally, we noticed that the substantial pair delocalization of the C-C-N bonded structure could play a significant role in the conversion of the biradical imine nitrile into both the propargylic and allenic forms via the T1 →S0 deactivation.

On the Notation of Catastrophes in the Framework of Bonding Evolution Theory: Case of Normal and Inverse Electron Demand Diels-Alder Reactions.

This paper generalizes very recent and unexpected findings [J. Phys. Chem. A, 2021, 125, 5152-5165] regarding the known "direct- and inverse-electron demand" Diels-Alder mechanisms. Application of bonding evolution theory indicates that the key electron rearrangement associated with significant chemical events (e. g., the breaking/forming processes of bonds) can be characterized via the simplest fold polynomial. For the CC bond formation, neither substituent position nor the type of electronic demand induces a measurable cusp-type signature. As opposed to the case of [4+2] cycloaddition between 1,3-butadiene and ethylene, where the two new CC single bonds occur beyond the transition state (TS) in the activated cases, the first CC bond occurs in the domain of structural stability featuring the TS, whereas the second one remains located in the deactivation path connecting the TS with the cycloadduct.

Photochemically Induced 1,3-Butadiene Ring-Closure from the Topological Analysis of the Electron Localization Function Viewpoint.

The electronic rearrangement featuring the photochemically-induced 1,3-cis-butadiene is discussed within a bonding evolution theory (BET) perspective based on the topological analysis of the electron localization function and Thom's catastrophe theory. The process involves the vertical singlet-singlet excitation S0 →S2 , and the subsequent deactivation implying the S2 /S1 and S1 /S0 conical intersection regions. BET results reveal that the new CC bond is finally formed on the S0 surface, as also recently found in the photochemical addition of two ethylenes [Phys. Chem. Chem. Phys. 23, 20598, (2021)].

Machine learning predictive classification models for the carcinogenic activity of activated metabolites derived from aromatic amines and nitroaromatics.

A 3D-QSAR study based on DFT descriptors and machine learning calculations is presented in this work. Our goal has been to build predictive models for classifying the carcinogenic activity of a set of aromatic amines (AA) and nitroaromatic (NA) compounds. As the main result, we stress that calculations must consider both the activated metabolites (derived from AA and NA species) and the water solvent to obtain reliable predictive classification models. We have obtained eight decision tree models that presented an accuracy of over 90% by using either Gázquez-Vela chemical potential (µ+) or the chemical hardness (η) of the activated metabolites in aqueous solvent.

On Electron Pair Rearrangements in Photochemical Reactions: 1,3-Cyclohexadiene Ring Opening.

1,3-Cyclohexadiene ring opening has been studied within the bonding evolution theory (BET) framework. We have focused on describing for the first time the electron pair rearrangements leading to the cis-1,3,5-hexatriene (HT) product from CHD. The nature of bonding in this process begins with the weakening of the double bonds in the Franck-Condon region. Along the 11B surface, the C-C sigma bond weakens. Meanwhile, its density redistributes toward the whole CHD ring, mainly over double bonds. Breaking of this bond occurs on the 21A surface due to the symmetrical splitting of pair density from this region. This density redistributes toward the reaction center once the pericyclic minimum is reached. The formation of the double bonds that characterize HT occurs gradually in the ground state. However, near the 21A/11A intersection, these bonds are partially established.

On the nature of bonding in the photochemical addition of two ethylenes: C-C bond formation in the excited state?

In this work, the 2s + 2s (face-to-face) prototypical example of a photochemical reaction has been re-examined to characterize the evolution of chemical bonding. The analysis of the electron localization function (as an indirect measure of the Pauli principle) along the minimum energy path provides strong evidence supporting that CC bond formation occurs not in the excited state but in the ground electronic state after crossing the rhombohedral S1/S0 conical intersection.

Unraveling the Bonding Nature Along the Photochemically Activated Paterno-Büchi Reaction Mechanism.

The photochemically activated Paterno-Büchi reaction mechanism following the singlet excited-state reaction path was analyzed based on a bonding evolution framework. The electronic rearrangements, which describe the mechanism of oxetane formation via carbon-oxygen attack (C-O), comprises of the electronic activation of formaldehyde and accumulation of pairing density on the O once the reaction system is approaching the conical intersection point. Our theoretical evidence based on the ELF topology shows that the C-O bond is formed in the ground-state surface (via C-O attack) returning from the S1 surface accompanied by 1,4-singlet diradical formation. Subsequently, the reaction center is fully activated near the transition state (TS), and the ring-closure (yielding oxetane) involves the C-C bond formation after the TS. For the carbon-carbon attack (C-C), both reactants (formaldehyde and ethylene) are activated, leading to C-C bond formation in the S1 excited state before reaching the conical intersection region. Finally, the C-O formation occurs in the ground-state surface, resulting from the pair density flowing primarily from the C to O atom.

A physiologic rise in cytoplasmic calcium ion signal increases pannexin1 channel activity via a C-terminus phosphorylation by CaMKII.

Pannexin1 (Panx1) channels are ubiquitously expressed in vertebrate cells and are widely accepted as adenosine triphosphate (ATP)-releasing membrane channels. Activation of Panx1 has been associated with phosphorylation in a specific tyrosine residue or cleavage of its C-terminal domains. In the present work, we identified a residue (S394) as a putative phosphorylation site by Ca2+/calmodulin-dependent kinase II (CaMKII). In HeLa cells transfected with rat Panx1 (rPanx1), membrane stretch (MS)-induced activation-measured by changes in DAPI uptake rate-was drastically reduced by either knockdown of Piezo1 or pharmacological inhibition of calmodulin or CaMKII. By site-directed mutagenesis we generated rPanx1S394A-EGFP (enhanced green fluorescent protein), which lost its sensitivity to MS, and rPanx1S394D-EGFP, mimicking phosphorylation, which shows high DAPI uptake rate without MS stimulation or cleavage of the C terminus. Using whole-cell patch-clamp and outside-out excised patch configurations, we found that rPanx1-EGFP and rPanx1S394D-EGFP channels showed current at all voltages between ±100 mV, similar single channel currents with outward rectification, and unitary conductance (∼30 to 70 pS). However, using cell-attached configuration we found that rPanx1S394D-EGFP channels show increased spontaneous unitary events independent of MS stimulation. In silico studies revealed that phosphorylation of S394 caused conformational changes in the selectivity filter and increased the average volume of lateral tunnels, allowing ATP to be released via these conduits and DAPI uptake directly from the channel mouth to the cytoplasmic space. These results could explain one possible mechanism for activation of rPanx1 upon increase in cytoplasmic Ca2+ signal elicited by diverse physiological conditions in which the C-terminal domain is not cleaved.

Are There Only Fold Catastrophes in the Diels-Alder Reaction Between Ethylene and 1,3-Butadiene?

This work revisits the topological characterization of the Diels-Alder reaction between 1,3-butadiene and ethylene. In contrast to the currently accepted rationalization, we here provide strong evidence in support of a representation in terms of seven structural stability domains separated by a sequence of 10 elementary catastrophes, but all are only of the fold type (F and F†), that is, C4H6 + C2H4:1-7-[FF]F[F†F†][F†F†][FF]F†-0: C6H10. Such an unexpected finding provides fundamental new insights opening simplifying perspectives concerning the rationalization of the CC bond formation in pericyclic reactions in terms of the simplest Thom's elementary catastrophe, namely, the one-(state) variable, one-(control) parameter function.

Introduction to celebrating Latin American talent in chemistry.

In celebration of the excellence and breadth of Latin American research achievements across the chemical sciences, we are delighted to present an introduction to the themed collection, Celebrating Latin American talent in chemistry.

A Close Look to the Oxaphosphetane Formation along the Wittig Reaction: A [2+2] Cycloaddition?

The Wittig reaction between triphenylphosphine methylide and benzaldehyde has been studied both from conceptual and computational approaches. The supernucleophilic character of ylide accounts for the feasibility of the initial nucleophilic attack. The nature of bonding driving the formation of the first oxaphosphetane (OPA) intermediate in such a domino reaction is examined within a topological-based bonding evolution theory perspective. The sequence of the electronic flow associated to the changes in electron density supports a rationalization via two main electronic stages characterizing the single kinetic step: first, the C-C bond formation, which takes place via donation of electron density of the ylide carbon to the carbonyl carbon of benzaldehyde at a C-C distance of 2.02 Å, is formally associated to the transition state region; then, the P-O bond formation via the donation of electron density of the nonbonding region of the carbonyl oxygen to phosphorus at a P-O distance of 2.06 Å is located at the end of the reaction path. The detailed picture of bonding patterns suggests that the OPA formation in the Wittig mechanism can be better understood in terms of a two-stage one-step mechanism beyond molecular orbital considerations behind the traditionally accepted [2+2] cycloaddition proposal.

Are one-step aromatic nucleophilic substitutions of non-activated benzenes concerted processes?

Aromatic nucleophilic substitution (SNAr) reactions of non-electrophilically activated benzenes have been studied within the Molecular Electron Density Theory (MEDT) at the B3LYP/6-311+G(d) computational level. These reactions, taking place through a one-step mechanism, present a high activation Gibbs free energy, ΔG≠ = 31.0 kcal mol-1, which decreases to 22.1 kcal mol-1 in the intramolecular process. A topological analysis of the electron localisation function along the reaction paths permits establishing the non-concerted nature of these SNAr reactions. A series of unstable structures, with similar electronic structures to those of Meisenheimer intermediates, are characterised. The present MEDT study makes it possible to establish that even these one-step SNAr reactions involving only two single bonds are non-concerted processes.

On the electron flow sequence driving the hydrometallation of acetylene by lithium hydride.

The sequence of the electronic flow driving the hydrometallation of acetylene by lithium hydride (and that of the opposite ß-hydride elimination reaction from the alkenyl metal intermediate), was examined within the perspective provided by the bonding evolution theory (BET). The analysis was based on the application of catastrophe theory to the changes of the electron localization function topology along the intrinsic reaction coordinate. The description of the electronic processes occurring on the process was represented in terms of topological structural stability domains (SSDs) and the associated elementary bifurcation catastrophes. Within such a framework of representation, the "evolution" of the system through the different SSDs reveals the key chemical events driving the transformation, including the large polarization effect as a consequence of Pauli repulsion between ions of the positive cationic metal on the hydride domain, the activation of the CC triple bond to attack the cationic center, and the agostic stabilizing interactions involving the hardest cationic metal, followed by the attack of the hydride center. These results contribute to emphasizing the intrinsic value and usefulness of using topological-based approaches and associated tools to increase our knowledge and understanding of the subtleties underlying the electronic flow as nuclei evolve along the reaction coordinate, providing detailed and complementary insights in comparison to other interpretative tool such those based on orbital-based representations, concerning the intimate nature of the electronic rearrangement of key mechanistic processes in chemistry. Graphical abstract The sequence of the electron flow (indicated by letters a and b) along the intrinsic reaction path for the hydrometallation of acetylene by lithium hydride to yield ethenyl lithium via a four-membered transition structure (TS), as determined within the bonding evolution theory to provide the key chemical events driven the changes in the key bonding patterns. Blue arrow Main event on the side HC ≡ CH + LiH → TS, red arrow the TS → HLiC=CH2 pathway, green arrows relative motion of nuclei along the imaginary frequency at the position of TS on the intrinsic reaction coordinate.

Reaction Electronic Flux Perspective on the Mechanism of the Zimmerman Di-π-methane Rearrangement.

The reaction electronic flux (REF) offers a powerful tool in the analysis of reaction mechanisms. Noteworthy, the relationship between aromaticity and REF can eventually reveal subtle electronic events associated with reactivity in aromatic systems. In this work, this relationship was studied for the triplet Zimmerman di-π-methane rearrangement. The aromaticity loss and gain taking place during the reaction is well acquainted by the REF, thus shedding light on the electronic nature of reactions involving dibenzobarrelenes.

Understanding the Highly Varying pKa of Arylamines. A Perspective from the Average Local Ionization Condensed-to-Atom Framework.

The highly varying experimental pKa values for 36 arylamines spanning 7 orders of magnitude is carefully examined. Within this framework, a valence condensed-to-atom model for the average ionization energy is introduced and tested. The theoretical approach is connected to orbital Fukui functions directly mapped into semilocal or regional site-specific responses. It is revealed that the average local ionization energies associated with the amino nitrogen atom is linearly correlated to the basicity of the substituted arylamines, properly reproducing the experimental ordering of basicity. The condensed-to-atom descriptor exhibits a high predictive power, providing a new direct reactivity evaluation of significant value.

Dual Emission of a Novel (P,N) Re(I) Complex: A Computational and Experimental Study on [P,N-{(C6H5)2(C5H4N)P}Re(CO)3Br].

The spectroscopic, electrochemical, and photophysical properties of the new complex [P,N-{(C6H5)2(C5H4N)P}Re(CO)3Br] are reported. The UV-vis spectrum in dichloromethane shows an absorption maximum centered at 315 nm and a shoulder at 350 nm. These absorption bands have been characterized to have MLCT character. Excitation at both wavelengths (maximum and shoulder) leads to an emission band centered at 550 nm. Cyclic voltammetry experiments show two ill-defined irreversible oxidation waves around +1.50 and 1.80 V that are assigned to Re(I)/Re(II) and Re(II)/Re(III) couples whereas an irreversible reduction signal centered at -1.80 V is likewise assigned to a ligand reduction process. These results support the proposal of the MLCT nature of the states implied by the emission of the complex. The luminescent decay fits to a biexponential function, where the lifetimes and emission quantum yields are dependent on the solvent polarity. DFT calculations suggest that dπ → π*pyridine and dπ → π*phenyl excited states may account for the existence of two decay lifetimes.

A computational and conceptual DFT study on the mechanism of hydrogen activation by novel frustrated Lewis pairs.

A computational and conceptual density functional theory (DFT) study on the mechanism of molecular hydrogen activation by a set of three frustrated Lewis pairs (FLPs) was performed at the ωB97X-D/6-311G(d,p) level of theory. A reduced model and other two prototypes derived from experimental data, based on the donor nitrogen and acceptor boron atoms, were used. Analysis based on the energy results, geometries and the global electron density transfer at the TSs made it possible to obtain some interesting conclusions: (i) despite the well-known very low reactivity of molecular hydrogen, the catalytic effectiveness of the three FLPs produces reactions with almost unappreciable activation energies; (ii) the reactions, being exothermic, follow a one-step mechanism via polarised TSs; (iii) there are neither substituent effects on the kinetics nor on the thermodynamics of these reactions; (iv) the activation of molecular hydrogen seems to be attained when the N-B distance in the FLP derivatives is around 2.74 Å; and (v) the proposed FLP model is consistent with the behaviour of the experimental prototypes. Finally, the ability of the three FLPs as efficient catalysts was evaluated studying the hydrogenation of acetylene to yield ethylene.