Explaining the High Catalytic Activity in Bis(indenyl)methyl Zirconium Cation Using Combined EDA-NOCV/QTAIM Approach.

The main purpose of this study is to elucidate some discrepancies already observed in the catalytic activity values of some zirconocene methyl cations. The EDA-NOCV scheme was employed for a theoretical description of the interactions between an ethylene molecule and five catalysts of zirconocene methyl cation. The nature of the chemical interactions has been elucidated through the QTAIM topological analysis. The steric hindrance due to the ligands was evaluated qualitatively by means of an IRI-based analysis and quantitively through Fisher information. The findings prove that the indenyl ligand seems to favor the orbital interaction between the ethylene molecule and the metal centre of zirconocene methyl cation. Both electrostatic and orbital contributions play a crucial role in stabilising the studied complexes. Based on the NOCV deformation density contributions, the strongest orbital interaction is reached with the bis(indenyl)methyl zirconium cation, which is the only one exhibiting covalent interactions. Especially, the strong contribution of π-back donation (occurring from the occupied orbitals of the zirconium atom to the π* anti-bonding orbital of ethylene) may be a key to understand why this catalyst has a higher polymerisation yield than the other studied catalysts. This work suggests a perspective for predicting values of catalytic activity when theoretically designing novel catalysts of zirconocene type.

Analysis in silico of chemical reactivity employing the local hyper-softness in some classic aromatic compounds, boron aromatic clusters and all-metal aromatic clusters.

In the current work, the authors analyzed and compared the chemical behavior of some (anti)aromatic compounds. The species selected are benzene and cyclobutadiene as the aromatic and antiaromatic classical examples, respectively. Next, the anion Al 4 2 - , which is the first all-metal molecule catalogued as aromatic and its non-metallic isoelectronic analog, B 4 2 - were also analyzed. The antiaromatic clusters Al 4 4 - and B 4 4 - were studied in form of lithium salts. And the end, the non-planar B 12 boron cluster and its isoelectronic analogs ( B 11 Be - and B 11 C + ) were considered for being analyzed under the same criterium. The analysis was realized employing the following descriptors: molecular electrostatic potential and local hyper-softness to get insights concerning local reactivity when facing reagents leading to ionic or covalent interactions. The results show that all the molecules analyzed presented some specific variations in their respective local reactivity despite being labeled as aromatic compounds. This analysis provides a notion that the local reactivity is more based on the nature of the atoms and the molecular geometry than the aromatic character by itself.

Physicochemical and Theoretical Characterization of a New Small Non-Metal Schiff Base with a Differential Antimicrobial Effect against Gram-Positive Bacteria.

Searching for adequate and effective compounds displaying antimicrobial activities, especially against Gram-positive bacteria, is an important research area due to the high hospitalization and mortality rates of these bacterial infections in both the human and veterinary fields. In this work, we explored (E)-4-amino-3-((3,5-di-tert-butyl-2-hydroxybenzylidene)amino) benzoic acid (SB-1, harboring an intramolecular hydrogen bond) and (E)-2-((4-nitrobenzilidene)amino)aniline (SB-2), two Schiff bases derivatives. Results demonstrated that SB-1 showed an antibacterial activity determined by the minimal inhibitory concentration (MIC) against Staphylococcus aureus, Enterococcus faecalis, and Bacillus cereus (Gram-positive bacteria involved in human and animal diseases such as skin infections, pneumonia, diarrheal syndrome, and urinary tract infections, among others), which was similar to that shown by the classical antibiotic chloramphenicol. By contrast, this compound showed no effect against Gram-negative bacteria (Klebsiella pneumoniae, Escherichia coli, and Salmonella enterica). Furthermore, we provide a comprehensive physicochemical and theoretical characterization of SB-1 (as well as several analyses for SB-2), including elemental analysis, ESMS, 1H and 13C NMR (assigned by 1D and 2D techniques), DEPT, UV-Vis, FTIR, and cyclic voltammetry. We also performed a computational study through the DFT theory level, including geometry optimization, TD-DFT, NBO, and global and local reactivity analyses.

The density polarization reveals directions of electron displacements due to the substituent effect: Analysis performed on a metal-organic Mo-Oxo catalyst.

Some Mo-oxo complexes bearing pyridine rings have the capability for dihydrogen production from water. However, energy barrier and overall energy vary depending on the effect exerted by several substituent groups located at different positions around one or more pyridine rings which are ligands of these compounds. Based on the Karunadasa and coworkers investigation where the para-position was experimentally tested in compounds derivatised from the 2,6-bis[1,1-bis(2-pyridil)ethyl]-pyridine oxo-molybdenum complex synthesized (Karunadasa et al., Nature, 2010, 464, 1329), we tested the combined effect of electron-withdrawing and electron-donating groups simulated as perturbations represented by point-charges. Then, we used the density polarization concept, δρ(r), a local reactivity descriptor corresponding to the partially integrated linear response function, χ(r, r') (a non-local reactivity descriptor), which is able to reveal different displacements of π-electrons on molecular structures. We perturbed the para-positions in the pentadentate ligand 2,6-bis[1,1-bis(2-pyridil)ethyl]-pyridine in the Mo-based complex by means of point-charges. They were located in three different configurations of the organic ligand (trans, geminal, and cis) which could help to explain energy barriers and overall energy of reactions catalyzed by this type of Mo-complexes. Our results indicate that the trans configuration of point-charges induces the most amount of fraction of electron shifted on the complex. A Mo-based complex bearing the same trans configuration for electron-withdrawing and electron-donating substituent groups (cyano and amino, respectively), leads to a kinetically more favorable H2 release than the cis or geminal configuration of the substituent groups aforementioned.

1,3-Dipolar Cycloadditions by a Unified Perspective Based on Conceptual and Thermodynamics Models of Chemical Reactivity.

The main aim in the present report is to gain a deeper understanding of typical 1,3-dipolar cycloadditions by means of three chemical reactivity models in a unified perspective: conceptual density functional theory, distortion/interaction, and reaction force analysis. The focus is to explore the information provided by each reactivity model and how they complement or reinforce each other. Our results showed that the Bell-Evans-Polanyi (BEP) relationship is fulfilled, which is consistent with the Hammond-Leffler postulate. The electronic chemical potential based analysis classifies the reactions as HOMO-, HOMO/LUMO-, and LUMO-controlled reactions as the activation energy increases. It seems likely that HOMO-controlled reaction shifts into LUMO-controlled one as the transition state (TS) position does from early into late. Therefore, the transition from HOMO- (and early TS) into LUMO-controlled (and late TS) is paid by shifting the overall energy change into an endothermic direction, thus supporting the fulfillment of the BEP principle. While thermodynamic models unveil that the distortion or structural rearrangements mainly drive the activation barriers rather than interaction or electronic rearrangements in accord with the distortion/interaction and reaction force analysis, respectively. It is also found that both models are consistent when energy associated with structural and electronic reordering from reaction force analysis is respectively confronted with destabilizing (distortion and Pauli repulsion) and stabilizing (electrostatic and orbital interactions) contributions from the distortion/interaction model, which, on the other hand, increases as low activation barrier and high exothermicity are converted into the high barrier and low exothermicity along with the BEP relation. Finally, the reaction force constant reveals that all 1,3-dipolar cycloaddition reactions proceed by a synchronous single-step mechanism, unveiling that the degree of synchronicity is quite the same in all reactions, confirming the statement that BEP is fulfilled for similar reactions proceeding by a quite alike degree of synchronicity.

A statistical thermodynamics view of electron density polarisation: application to chemical selectivity.

A fundamental link between conceptual density functional theory and statistical thermodynamics is herein drawn, showing that intermolecular electrostatic interactions can be understood in terms of effective work and heat exchange. From a more detailed analysis of the heat exchange in a perturbation theory framework, an associated entropy can be subsequently derived, which appears to be a suitable descriptor for the local polarisability of the electron density. A general rule of thumb is evidenced: the more the perturbation can be spread, both through space and among the excited states, the larger the heat exchange and entropy.

Scrutinizing the substituent effect on Mo-based electrocatalysts for molecular hydrogen release through axial-equatorial decomposition: a DFT study.

Based on the experimental precedent discovered by Kuranadasa and coworkers [H. I. Karunadasa et al., Nature, 2010, 464, 1329] to produce dihydrogen from water electrocatalyzed by 2,6-bis[1,1-bis(2-pyridyl)ethyl]-pyridine oxo-molybdenum complexes, we performed an extensive analysis to study the substituent group effect of derivatised compounds coming from the before mentioned Mo-based metal-organic cations in terms of two kinds of substitutions: axial and equatorial at the para-position of pyridine rings; several conceptual tools were used to back up our conclusions. We found that each type of substituent group (electron-withdrawing and electron-donating ones) exerts an independent influence on energetic parameters (energy barrier and overall energy). This opens the chance to search for a synergistic effect by combining these opposite behaviours of these substituents located in the equatorial and axial para-positions of pyridine rings to computationally modulate the aforementioned energetic parameters. This procedure will make easier the proposal of new catalysts to favour either kinetically or thermodynamically or in both ways the production of dihydrogen from water. Additionally, we encompassed a key point: the number of solvent molecules, so that including their presence in further investigations is mandatory.

Effect of the exchange-correlation functional on the synchronicity/nonsynchronicity in bond formation in Diels-Alder reactions: a reaction force constant analysis.

In this paper, we assess the performance of 24 density functional theory (DFT) based methods classified into 5 categories (GGA, MGGA, HGGA, HMGGA and DHGGA) in predicting reaction energetics, transition state geometries, and the degree of synchronicity/nonsynchronicity in the formation of two new C-C single-bonds in three Diels-Alder reactions between symmetrically and unsymmetrically substituted cyanoethylenes and cyclopentadiene, which gradually proceed from fully synchronous to highly asynchronous concerted mechanisms. This important concept in reaction mechanisms is revealed by the fine structure of the reaction force constant κ(ξ) along the transition region. Some wave function theory (WFT) based methods are also assessed against the CCSD(T) and CCSD benchmarks for the energy and geometry, respectively. The results and the statistical analysis of the errors confirm the robustness of SCS-MP2 (a WFT-based method) as one of the most reliable computational approaches. Regarding DFT-based methods, hybrid exchange-correlation functionals combined with medium-range electron correlation effects or long-range corrected exchange appear as the best performing methods, highlighting both M11 and M06-2X, since a certain percentage of exact Hartree-Fock exchange could counterbalance the delocalization errors that affect pure functionals. Thus, they reliably describe energetics, geometries and the degree of synchronicity in the formation of new C-C single bonds in Diels-Alder reactions. Noticeably, moderate performance for double hybrid functionals and poor performance for the most popular B3LYP method were found as well.

A generalized operational formula based on total electronic densities to obtain 3D pictures of the dual descriptor to reveal nucleophilic and electrophilic sites accurately on closed-shell molecules.

By means of the conceptual density functional theory, the so-called dual descriptor (DD) has been adapted to be used in any closed-shell molecule that presents degeneracy in its frontier molecular orbitals. The latter is of paramount importance because a correct description of local reactivity will allow to predict the most favorable sites on a molecule to undergo nucleophilic or electrophilic attacks; on the contrary, an incomplete description of local reactivity might have serio us consequences, particularly for those experimental chemists that have the need of getting an insight about reactivity of chemical reagents before using them in synthesis to obtain a new compound. In the present work, the old approach based only on electronic densities of frontier molecular orbitals is replaced by the most accurate procedure that implies the use of total electronic densities thus keeping consistency with the essential principle of the DFT in which the electronic density is the fundamental variable and not the molecular orbitals. As a result of the present work, the DD will be able to properly describe local reactivities only in terms of total electronic densities. To test the proposed operational formula, 12 very common molecules were selected as the original definition of the DD was not able to describe their local reactivities properly. The ethylene molecule was additionally used to test the capability of the proposed operational formula to reveal a correct local reactivity even in absence of degeneracy in frontier molecular orbitals. © 2016 Wiley Periodicals, Inc.

The effect of {O,N}=X⋯M={Ti,Zr,Hf} interactions on the sensitivity of CNO2 trigger bonds in FOX-7: Approach based on the QTAIM/EDA-NOCV analysis.

The local chemical reactivity of FOX-7 (1,1-diamino-2,2-nitroethylene, also known as DADNE from DiAminoDiNitroEthylene) was elucidated through a quantitative study of the electrostatic potential on the molecular surface, topological analysis based on Bader's theory, and the EDA-NOCV method. Unlike (O2N)2CC(NH2)H2N⋯Cp2MCH3+ complexes, which exhibit both σ-donor and π-acceptor features, the situation is different concerning the (H2N)2CC(NO2)(O)NO⋯Cp2MCH3+ complexes, where both charge transfers correspond to the σ-donation. The two charge transfers reinforce each other, resulting in increased stability for (H2N)2CC(NO2)(O)NO⋯Cp2MCH3+. This seems to strengthen the (H2N)2CC(NO2)(O)NO⋯M={Ti,Zr,Hf} bond, which may explain the high stability of (H2N)2CC(NO2)(O)NO⋯Cp2MCH3+ compared to (O2N)2CC(NH2)-H2N⋯Cp2MCH3+. Results from topological analysis revealed that the decreased sensitivity to decomposition of CNO2 bonds depends on the chemical nature of the interacting metal, and the best achievements are obtained for the Hf-based complex. Our results demonstrate that the interaction of M={Ti,Zr,Hf} with CNO2 is more favourable than that with CNH2, this specific action on the trigger bond may support the use of Metallocene Methyl Cations (MMC) as possible neutralisers.

Aromaticity of six-membered nitro energetic compounds through molecular electrostatic potential, magnetic, electronic delocalization and reactivity-based indices.

The electron density depletion associated with π-hole at the ring center typical of energetic compounds was clearly revealed by the molecular electrostatic potential (ESP). In addition, the spatial arrangement of NO2 groups appears to affect the ESP value in the ring center, and therefore sensitivity to detonation. Indeed, for monocyclic nitrobenzene compounds with the same number of NO2 groups, the ESP value in the ring center decreases as the NO2 groups are more closely spaced. As expected, the central rings become less aromatic as NO2 groups are added. The MCI, PDI, PLR, NICSzz(1), FLU indices are all strongly correlated with the ESP values observed in the ring center of the set of nitrobenzenes. Aromaticity indices based on electron delocalization criteria appear to be very sensitive to small variations in aromaticity. Among magnetic-based indices, only NICSzz(1) is capable to predict small changes in aromaticity. The PLR index derived from conceptual DFT is quite relevant for predicting small variations in aromaticity. According to our results, the most suitable aromaticity index is not based on a single criterion, and that selecting it is more subtle. Therefore, it is important to combine information from several criteria to obtain a more complete description of the aromaticity of the studied compounds.

Theoretical investigation of the effect of O⋯M={Ti,Zr,Hf} interactions on the sensitivity of energetic N-nitro compounds.

This paper outlines the role of intermolecular interactions involving group 4 transition metals in stabilising the N-NO2 trigger bonds. Minimising sensitivity is the foremost priority in designing energetic compounds. A quantitative analysis with Molecular Electrostatic Potential (MEP) evidenced anomalies arising from the marked depletion of negative charge distribution of RDX and HMX. The Energy Decomposition Analysis with Natural Orbitals for Chemical Valence (EDA-NOCV) results reveal that the electrostatic and orbital contributions are the dominant factors driving the assembly of the M={Ti,Zr,Hf}-based complexes. Sensitivity of the N-NO2 trigger bonds is investigated by using the Quantum Theory of Atoms in Molecules (QTAIM). The QTAIM topological analysis showed that the O⋯M={Ti,Zr,Hf} interaction strengthens these trigger bonds, revealing an increased stability to decomposition. This effect is more marked in the Hf- and Zr-based complexes. Finally, the results based on Interaction Indicator Region (IRI) are fully consistent with those generated from QTAIM analysis.

The Dual Descriptor Reveals the Janus-Faced Behaviour of Diiodine.

The Janus-faced ligand behavior of diiodine (I2) was evidenced after applying the dual descriptor (DD or second-order Fukui function), thus providing additional support to the work performed by Rogachev and Hoffmann in 2013. Along with its capacity to reveal sites susceptible to undergo attacks simultaneously of nucleophilic and electrophilic types, another advantage of DD lies in being an orbital-free descriptor. That means it is based only upon total electron densities when written in its most accurate operational formula. This quality is not exclusive of DD because when Fukui functions are written in terms of electron densities instead of densities of frontier molecular orbitals, they become orbital-free descriptors too. Furthermore, the present work is an application of the generalized operational formula of the dual descriptor published in 2016 that takes into account any possible degeneracy in frontier molecular orbitals. As a proof about capabilities of DD, the possible sites for a favorable interaction between I2 with two organometallic compounds [Rh2(O2CCF3)4] and [(C8H11N2)Pt (CH3)] were correctly revealed by overlapping the biggest lobe for receiving nucleophilic attacks of one molecule with the biggest lobe for receiving electrophilic attacks of the other molecule, so allowing to predict the same coordination modes as experimentally known: linear "end-on" for the [(C8H11N2)Pt (CH3)]…I2, and bent "end-on" for the [Rh2(O2CCF3)4]…I2 interactions.

A combined QTAIM/IRI topological analysis of the effect of axial/equatorial positions of NH2 and CN substituents in the [(PY5Me2)MoO]+ complex.

By means of the Interaction Region Indicator (IRI) and Quantum Theory of Atoms in Molecules (QTAIM), the influence exerted by NH2 (amino) and CN (cyano) as electron donor and electron acceptor substituent groups, respectively, located at para-positions of axial and equatorial pyridine rings of derivatized complexes coming from the [(PY5Me2)MoO]+ complex during the hydrogen molecular release in the gas phase was analyzed. In any case, a H-H covalent bond is forming at the transition state, with a strengthening of the electron density of 5.5% when the substituent group involved is NH2 at the para-position of the axial pyridine ring. However, there was no difference between NH2 and CN when these substituent groups are located at the para-positions of the equatorial pyridine rings. The topological properties of electron densities from the QTAIM are not perturbed by the electron donor and electron acceptor nature of the substituents, even when these substituent groups are located at the axial or equatorial pyridine rings of the Mo-based complex.

KID Procedure Applied on the [(PY5Me2)MoO]+ Complex.

The KID (Koopmans in DFT) protocol usually applies in organic molecules of the closed-shell type. We used the KID procedure on an open-shell Mo-based system for the first time to choose the most suitable density functional to compute global and local reactivity descriptors obtained from the conceptual density-functional theory (DFT). From a set of 18 density functionals, spread from the second until the fourth rung of Jacob's ladder: BLYP, BP86, B97-D, MN12-L, MN15-L, M06-L, M11-L, CAM-B3LYP, PBE0, B3LYP, N12-SX, M06-2X, MN15, MN12-SX, ωB97X-D, M11, LC-ωHPBE, and APFD, we concluded that CAM-B3LYP provides the best outcome, and in the second place, M06-2X. Because the vertical first ionization potential and vertical first electron affinity in the ground state (gs) are defined as follows I = E gs(N - 1) - E gs(N) and A = E gs(N) - E gs(N + 1), where E gs(N - 1), E gs(N), and E gs(N + 1) correspond to energies of the system bearing N, N + 1, and N - 1 electrons, along with Koopmans' theorem (KT) given by I ≈ -εHOMO (εHOMO, highest occupied molecular orbital energy) and A ≈ -εLUMO (εLUMO, lowest unoccupied molecular orbital energy), the deviation from the KT was performed by the use of the index, such that J I = |E gs(N - 1) - E gs(N) + εHOMO| and J A = |E gs(N) - E gs(N + 1) + εLUMO|, which are absolute deviations from the perspective of I and A, respectively. Furthermore, the εSOMO (SOMO: singly-occupied molecular orbital energy) leads us to another index given by |ΔSL| = |εSOMO - εLUMO|. Therefore, J HL and |ΔSL| are indexes defined to evaluate the quality of the KT when employed within the context of quantum chemical calculations based on DFT and not the Hartree-Fock theory. We propose the index that could be more suitable to choose the most proper density functional because the J HL and |ΔSL| are independent indexes.

Assessment of ten density functionals through the use of local hyper-softness to get insights about the catalytic activity : Iron-based organometallic compounds for ethylene polymerization as testing molecules.

Ten functionals were used to assess their capability to compute a local reactivity descriptor coming from the Conceptual Density Functional Theory on a group of iron-based organometallic compounds that have been synthesized by Zohuri, G.H. et al. in 2010; these compounds bear the following substituent groups: H-, O2N- and CH3O- at the para position of the pyridine ring and their catalytic activities were experimentally measured by these authors. The present work involved a theoretical analysis applied on the aforementioned iron-based compounds thus leading to suggest a new 2,6-bis(imino)pyridine catalyst based on iron(II) bearing a fluorine atom whose possible catalytic activity is suggested to be near the catalytic activity of the complex bearing a hydrogen atom as a substituent group by means of the so called local hyper-softness (LHS) thus opening a chance to estimate a possible value of catalytic activity for a new catalyst that has not been synthesized yet without simulating the entire process of ethylene polymerization. Since Conceptual DFT is not a predictive theory, but rather interpretative, an analysis of the used reactivity descriptor and its dependence upon the level of theory was carried in the present work, thus revealing that care should be taken when DFT calculations are used for these purposes.

A 3D visualization of the substituent effect : A brief analysis of two components of the operational formula of dual descriptor for open-shell systems.

Six organometallic compounds coming from a basic Mo-based complex were analyzed from the perspective of the dual descriptor in order to detect subtle influences that a substituent group could exert on the reactive core at a long range. Since the aforementioned complexes are open-shell systems, the used operational formula for the dual descriptor is that one defined for those aforementioned systems, which was then compared with spin density. In addition, dual descriptor was decomposed into two terms, each of which was also applied on every molecular system. The obtained results indicated that components of dual descriptor could become more useful than the operational formula of dual descriptor because differences exerted by the substituents at the para position were better detected by components of dual descriptor rather than the dual descriptor by itself.

Solvent effect on the degree of (a)synchronicity in polar Diels-Alder reactions from the perspective of the reaction force constant analysis.

In this work, we computationally evaluated the influence of six different molecular solvents, described as a polarizable continuum model at the M06-2X/6-31+G(d,p) level, on the activation barrier/reaction rate, overall energy change, TS geometry, and degree of (a)synchronicity of two concerted Diels-Alder cycloadditions of acrolein (R1) and its complex with Lewis acid acrolein···BH3 (R2) to cyclopentadiene. In gas-phase, we found that both exothermicity and activation barrier are only reduced by about 2.0 kcal mol-1, and the asynchronicity character of the mechanism is accentuated when BH3 is included. An increment in the solvent's polarity lowers the activation energy of R1 by 1.3 kcal mol-1, while for R2 the reaction rate is enhanced by more than 2000 times at room temperature (i.e., the activation energy decreases by 4.5 kcal mol-1) if the highest polar media is employed. Therefore, a synergistic effect is achieved when both external agents, i.e., Lewis acid catalyst and polar solvent, are included together. This effect was ascribed to the ability of the solvent to favor the encounter between cyclopentadiene and acrolein···BH3. This was validated by the asymmetry of the TS which becomes highly pronounced when either both or just BH3 is considered or the solvent's polarity is increased. Finally, the reaction force constant κ(ξ) reveals that an increment in the solvent's polarity is able to turn a moderate asynchronous mechanism of the formation of the new C-C σ-bonds into a highly asynchronous one. Graphical abstract A synergistic effect is achieved when both external agents, i.e., Lewis acid catalyst and polar solvent, are included together: lowered energy barriers and increased asynchronicities.