In this paper, we present a new finding, the basis electronic activity (BEA) of molecular systems; it corresponds to the significant, although nonreactive, vibrationally induced electronic activity that takes place in any molecular system. Although the molecule's BEA is composed of an equal number of local contributions as the vibrational degrees of freedom, our results indicate that only stretching modes contribute to it. To account for this electronic activity, a new descriptor, the bond electronic flux (BEF), is introduced. The BEF combined with the force constant of the potential well hosting the electronic activity gives rise to the effective bond reactivity index (EBR), which turns out to be the first density functional theory-based descriptor that simultaneously accounts for structural and electronic effects. Besides quantifying the bond reactivity, EBR provides a basis to compare the reactivities of bonds inserted in different chemical environments and paves the way for the exertion of selective control to enhance or inhibit their reactivities. The new concepts formulated in this paper and the associated computational tools are illustrated with characterization of the BEA of a set of representative molecules. In all cases, the BEFs follow the same linear pattern, whose slopes indicate the intensity of the electronic activity and quantify the reactivity of chemical bonds.
The chiral-induced spin selectivity effect (CISS) is a breakthrough phenomenon that has revolutionized the field of electrocatalysis. We report the first study on the electron spin-dependent electrocatalysis for the oxygen reduction reaction, ORR, using iron phthalocyanine, FePc, a well-known molecular catalyst for this reaction. The FePc complex belongs to the non-precious catalysts group, whose active site, FeN4, emulates catalytic centers of biocatalysts such as Cytochrome c. This study presents an experimental platform involving FePc self-assembled to a gold electrode surface using chiral peptides (L and D enantiomers), i.e., chiro-self-assembled FePc systems (CSAFePc). The chiral peptides behave as spin filters axial ligands of the FePc. One of the main findings is that the peptides' handedness and length in CSAFePc can optimize the kinetics and thermodynamic factors governing ORR. Moreover, the D-enantiomer promotes the highest electrocatalytic activity of FePc for ORR, shifting the onset potential up to 1.01â V vs. RHE in an alkaline medium, a potential close to the reversible potential of the O2 /H2 O couple. Therefore, this work has exciting implications for developing highly efficient and bioinspired catalysts, considering that, in biological organisms, biocatalysts that promote O2 reduction to water comprise L-enantiomers.
The critical and vanishing points of the reaction force F(ξ) = -dV(ξ)/dξ yield five important coordinates (ξR, ξR* , ξTS, ξP* , ξP) along the intrinsic reaction coordinate (IRC) for a given concerted reaction or reaction step. These points partition the IRC into three well-defined regions, reactants (ξRâξR* ), transition state (ξR* âξP* ), and products (ξP* âξP), with traditional roles of mostly structural changes associated with the reactants and products regions and mostly electronic activity associated with the transition state (TS) region. Following the evolution of chemical bonding along the IRC using formal descriptors of synchronicity, reaction electron flux, Wiberg bond orders, and their derivatives (or, more precisely, the intensity of the electron activity) unambiguously indicates that for nonsynchronous reactions, electron activity transcends the TS region and takes place well into the reactants and products regions. Under these circumstances, an extension of the TS region toward the reactants and products regions may occur.
The synthesis of novel aluminyl anion complexes has been well exploited in recent years. Moreover, the elucidation of the structure and reactivity of these complexes opens the path toward a new understanding of low-valent aluminum complexes and their chemistry. This work computationally treats the substituent effect on aluminyl anions to discover suitable alternatives for H2 activation at a high level of theory utilizing coupled-cluster techniques extrapolated to the complete basis set. The results reveal that the simplest AlH2- system is the most reactive toward the activation of H2, but due to the low steric demand, severe difficulty in the stabilization of this system makes its use nonviable. However, the results indicate that, in principle, aluminyl systems with -C, -CN, -NC, and -N chelating centers would be the best choices of ligand toward the activation of molecular hydrogen by taking care of suitable steric demand to prevent dimerization of the catalysts. Furthermore, computations show that monosubstitution (besides -H) in aluminyl anions is preferred over disubstitution. So our predictions show that bidentate ligands may yield less reactive aluminyl anions to activate H2 than monodentate ones.
Hydrogen atom transfer (HAT) is a crucial step in the physiological conversion of dopamine into norepinephrine catalyzed by dopamine ß-monooxygenase. The way the reaction takes place is unclear, and a rational explanation on how the electronic activity drives the HAT seems to be necessary. In this work, we answer this question using the reaction electronic flux (REF), a DFT-based descriptor of electronic activity. Two reaction mechanisms will be analyzed using the REF's decomposition in polarization and electron transfer effects. Results show that both mechanisms proceed as follows: (1) polarization effects initiate the reactions producing structural distortions; (2) electron transfer processes take over near the transition states, triggering specific chemical events such as bond forming and breaking which are responsible to push the reactions toward the products; (3) after passing the transition state, polarization shows up again and drives the relaxation process toward the product. Similar polarization effects were observed in both reactions, but they present an opposite behavior of the electronic transfer flux disclosing the fact that electron transfer phenomena govern the reaction mechanisms.
Dopamine , Norepinephrine , Electron Transport , Electronics , Hydrogen/chemistry
Aluminyl anions are low-valent aluminum species bearing a lone pair of electrons and a negative charge. These systems have drawn recent synthetic interest for their nucleophilic nature, allowing for the activation of σ-bonds, and have been proposed as a pathway to hydrogen energy storage. In this research, we provide high-level ab initio geometries and energies for both the simplest aluminyl anion (AlH2-) and several substituted derivatives. Geometries are reported using the gold-standard CCSD(T)/aug-cc-pV(T+d)Z level of theory. Energies were extrapolated to the complete basis set limit through the focal point approach, utilizing coupled-cluster methods through perturbative quadruples and basis sets up to five-ζ quality. Geometries were rationalized using electrostatic, steric, and orbital donation effects. The donation from substituents to Al is accompanied by back-donation effects, a property traditionally thought of in transition-metal systems. Stereoelectronic effects through the secondary orbital interaction play a fundamental role in stabilizing these low-valent aluminum compounds and would likely also affect the feasibility of their use within several industrial applications. The energetic analysis of the formation of each substituted anion is rationalized as the result of three energetic schemes. The effectiveness of these schemes for determining the relative formation energies is discussed.
Aluminyl anions are low-valent, anionic, and carbenoid aluminum species commonly found stabilized with potassium cations from the reaction of Al-halogen precursors and alkali compounds. These systems are very reactive toward the activation of σ-bonds and in reactions with electrophiles. Various research groups have detected that the potassium atoms play a stabilization role via electrostatic and cation ⯠π interactions with nearby (aromatic)-carbocyclic rings from both the ligand and from the reaction with unsaturated substrates. Since stabilizing Kâ¯H bonds are witnessed in the activation of this class of molecules, we aim to unveil the role of these metals in the activation of the smaller and less polarizable H2 molecule, together with a comprehensive characterization of the reaction mechanism. In this work, the activation of H2 utilizing a NON-xanthene-Al dimer, [K{Al(NON)}]2 (D) and monomeric, [Al(NON)]- (M) complexes are studied using density functional theory and high-level coupled-cluster theory to reveal the potential role of K+ atoms during the activation of this gas. Furthermore, we aim to reveal whether D is more reactive than M (or vice versa), or if complicity between the two monomer units exits within the D complex toward the activation of H2 . The results suggest that activation energies using the dimeric and monomeric complexes were found to be very close (around 33â kcal mol-1 ). However, a partition of activation energies unveiled that the nature of the energy barriers for the monomeric and dimeric complexes are inherently different. The former is dominated by a more substantial distortion of the reactants (and increased interaction energies between them). Interestingly, during the oxidative addition, the distortion of the Al complex is minimal, while H2 distorts the most, usually over 0.77 Δ E d i s t ≠ . Overall, it is found here that electrostatic and induction energies between the complexes and H2 are the main stabilizing components up to the respective transition states. The results suggest that the K+ atoms act as stabilizers of the dimeric structure, and their cooperative role on the reaction mechanism may be negligible, acting as mere spectators in the activation of H2 . Cooperation between the two monomers in D is lacking, and therefore the subsequent activation of H2 is wholly disengaged.
The analysis of the reaction force and its topology has provided a wide range of fruitful concepts in the theory of chemical reactivity over the years, allowing to identify chemically relevant regions along a reaction profile. The reaction force (RF), a projection of the Hellmann-Feynman forces acting on the nuclei of a molecular system onto a suitable reaction coordinate, is partitioned using the interacting quantum atoms approach (IQA). The exact IQA molecular energy decomposition is now shown to open a unique window to identify and quantify the chemical entities that drive or retard a chemical reaction. The RF/IQA coupling offers an extraordinarily detailed view of the type and number of elementary processes that take reactants into products, as tested on two sets of simple reactions.
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.
Quantum chemical calculations were used to study the mechanism of Diels-Alder reactions involving chiral anthracenes as dienes and a series of dienophiles. The reaction force analysis was employed to obtain a detailed scrutiny of the reaction mechanisms, it has been found that thermodynamics and kinetics of the reactions are quite consistent: the lower the activation energy, the lower the reaction energy, thus following the Bell-Evans-Polanyi principle. It has been found that activation energies are mostly due to structural rearrangements that in most cases represented more than 70% of the activation energy. Electronic activity mostly due to changes in σ and π bonding were revealed by the reaction electronic flux (REF), this property helps identify whether changes on σ or π bonding drive the reaction. Additionally, new global indexes describing the behavior of the electronic activity were introduced and then used to classify the reactions in terms of the spontaneity of their electronic activity. Local natural bond order electronic population analysis was used to check consistency with global REF through the characterization of specific changes in the electronic density that might be responsible for the activity already detected by the REF. Results show that reactions involving acetoxy lactones are driven by spontaneous electronic activity coming from bond forming/strengthening processes; in the case of maleic anhydrides and maleimides it appears that both spontaneous and non-spontaneous electronic activity are quite active in driving the reactions.
Qualitative relationships between the reaction force constant κ(ξ), the second derivative of the potential energy V(ξ), and the reactive vibrational mode that drives the reaction in the transition state region have been used in the past to measure the synchronicity of key chemical events that lead a chemical reaction. In this work, we provide a formal demonstration that κ(ξ) can be expressed in terms of the frequencies of normal modes at each point of the reaction path. This produce a decomposition of κ(ξ) that is used to analyze few representatives chemical reactions, a series of intramolecular proton transfer on formic, thioformic and dithioformic acids, and an intermolecular double proton transfer in the HNS2:H2O complex. It has been found that this partitioning allows to identify unambigously the reactive mode that drives the reaction at each point along the reaction coordinate thus giving relevant and detailed information about the mechanism of the chemical reactions under study.
The non-dissociative and dissociative adsorptions of zwitterionic Gly on graphene oxide (GO) was studied in the framework of DFT using a cluster model approach. In this work, the interaction with an epoxy group of GO basal plane was mainly considered. As a comparison, the non-dissociative and dissociative adsorptions of neutral Gly were also taken into account. The non-dissociative adsorption modes for zwitterionic and neutral Gly conformers show binding energies of 12.2 and 14.4 kcal mol-1, respectively. These molecules are thought to remain over the GO surface due to attractive noncovalent interactions. Two dissociative adsorption modes, for Z-Gly and N-Gly, show smaller binding energies of 7.2 and 8.4 kcal mol-1, where the deprotonated species links strongly through a C-O or C-N covalent bond to the GO surface. The results obtained in the present theoretical approach to the glycine/graphene oxide system support the fact that glycine can be attached to epoxy groups of graphene oxide basal planes in addition to the anchoring on edge oxidation groups. In summary, we conclude that glycine can be used as a reducing agent as well as a functionalizer of GO sheets.
Glycine/chemistry , Graphite/chemistry , Models, Chemical , Adsorption
The chemistry of low-valent main group elements has attracted much attention in the past decade. These species are relevant because they have been able to mimic transition metal behavior in catalytic applications, with decreased material costs and diminished toxicity. In this contribution, we study the L1EH catalysts (E = Si(II), Ge(II), Sn(II), and Pb(II); L1 = [ArNC(Me)CHC(Me)NAr] with Ar = 2,6-iPr2C6H3) for the formation of formic acid derivatives through hydroboration of CO2. Detailed characterization of relevant structures on the potential energy surface enabled us to rationalize different paths for the hydroboration of CO2. Interestingly, it was found that according to the activation energies for the whole catalytic cycle, the process of transformation of CO2 becomes more favored going down group 14. However, an effective energetic decrease for the process (taking as the reference the uncatalyzed reaction between CO2 and HBpin) is evidenced just from the germanium analogue. The trend in reactivity found in the present study is a direct consequence of the change in the central main group element, enabling enhanced polar character of the E-H (L1EH in the CO2 activation step) and E-O (metal formates in the hydroboration step) bonds as the atomic radius increases. The transient stabilization of reaction intermediates found in the hydroboration step was rationalized through the non-covalent interaction index (NCI) and symmetry-adapted perturbation theory (SAPT). This computational study highlights the reactivity trends in group-14-based hydride catalysts in hydrometalation and posterior hydroboration to form formic acid intermediates. We hope that this study will motivate further experimental work in low-valent lead chemistry.
The insertion of CO into hydrogenated boron-amidine 1 yielding five-membered diazaborolone (5DAB) 3 has been studied within the molecular electron density theory (MEDT) at the DFT ωB97X-D/6-311G(d,p) level. This is a domino process comprised of two consecutive reactions: (i) the dehydrogenation of 1 yielding the frustrated Lewis pair (FLP) boron-amidine 4, which quickly equilibrates with four-membered diazaborolone (4DAB) 2; and (ii) the addition of CO into FLP 4, yielding the final 5DAB 3. Analysis of the Gibbs free energies indicates that the extrusion of H2 demands a high ΔG≠ of 28.6 kcal·mol-1, being endergonic by 6.7 kcal·mol-1. The subsequent addition of CO into FLP 4 presents a low ΔG≠ of 15.0 kcal·mol-1; formation of 5DAB 3 being exergonic by -5.7 kcal·mol-1 from hydrogenated boron-amidine 1. An analysis of the bonding changes along the insertion of CO in a smaller FLP model indicates that this reaction can be considered a [4 + 1] cycloaddition reaction taking place via a five-membered pseudocyclic transition state associated with a two-stage one-step mechanism. Analysis of the conceptual DFT reactivity indices suggests that the initial attack of CO on FLP 4 is an acid/base process in which the carbenoid carbonyl character allows CO to participate as a Lewis base, rather than a nucleophilic/electrophilic interaction. The results arising from the analysis of the Parr functions, however, coincide with this behaviour.
Fullerenes, in particular C60, are important molecular entities in many areas, ranging from material science to medicinal chemistry. However, chemical transformations have to be done in order to transform C60 in added-value compounds with increased applicability. The most common procedure corresponds to the classical Diels-Alder cycloaddition reaction. In this research, a comprehensive study of the electronic activity that takes place in the cycloaddition between C60 and cyclopentadiene toward the [5,6] and [6,6] reaction pathways is presented. These are competitive reaction mechanisms dominated by σ and π fluctuating activity. To better understand the electronic activity at each stage of the mechanism, the reaction force (RF) and the symmetry-adapted reaction electronic flux (SA-REF, JΓi(ξ)) have been used to elucidate whether π or σ bonding changes drive the reaction. Since the studied cycloaddition reaction proceeds through a Cs symmetry reaction path, two SA-REF emerge: JA'(ξ) and JA''(ξ). In particular, JA'(ξ) mainly accounts for bond transformations associated with π bonds, while JA''(ξ) is sensitive toward σ bonding changes. It was found that the [6,6] path is highly favored over the [5,6] with respect to activation energies. This difference is primarily due to the less intensive electronic reordering of the σ electrons in the [6,6] path, as a result of the pyramidalization of carbon atoms in C60 (sp2 â sp3 transition). Interestingly, no substantial differences in the π electronic activity from the reactant complex to the transition state structure were found when comparing the [5,6] and [6,6] paths. Partition of the kinetic energy into its symmetry contributions indicates that when a bond is being weakened/broken (formed/strengthened) non-spontaneous (spontaneous) changes in the electronic activity occur, thus prompting an increase (decrease) of the kinetic energy. Therefore, contraction (expansion) of the electronic density in the vicinity of the bonding change is expected to take place.
A computational DFT study of the reaction mechanism of hydrogenation and hydration of carbon dioxide is presented. It has been found that hydrogenation and hydration are endoenergetic reactions that are carried out in two steps, passing by a stable intermediate that is surrounded by energy barriers of 70 kcal/mol and 10 kcal/mol for hydrogenation and 50 kcal/mol and 10 kcal/mol for hydration. Using the reaction force analysis, we were able to characterize the physical nature of the activation barriers and found that activation energies are mostly due to structural rearrangements. An interesting difference in the reaction mechanisms disclosed by the reaction force and electronic flux analyses is that while in the hydrogenation reaction the mechanisms is conditioned by the H2 cleavage with a high energy barrier, in the hydration reaction the formation of a transient four member ring structure favoring an attractive local hydrogen bond interaction pushes the reaction toward the product with a considerably lower energy barrier.
The mechanism of ethylene polymerization by means of neutral methallyl-Ni(ii) complexes has been studied by quantum chemical calculations. Two isomer complexes having different ligand functionalization at the ortho or para position, and co-activated with trispentafluorephenylborane [B(C6F5)3], were studied according to the Cossee-Arlman's mechanism. Comparison of the reaction mechanism of both isomers shows that energy barriers strongly depend on ligand-functionalization and are mostly due to structural rearrangements. In addition, it was found that para-functionalization can be distinguished by favorable σ-donation whereas ortho-functionalization is more prone to a π back-donation process. Our results concerning the polymerization process for ortho and para isomers not only provide a theoretical perspective of available experimental data, but also explain the experimentally observed higher molecular weight of the methallyl-Ni(ii) ortho isomer co-catalyzed by B(C6F5)3, revealing the role of ligand-functionalization in polyethylene production.
An analytic potential function consistent with the Marcus equation for activation energy is formulated and used to reveal new insights into the activation process in chemical reactions. As for the Marcus equation, the new potential function depends only on two parameters, the reaction energy and the activation energy (or the so-called Marcus intrinsic activation energy). Combination of the Marcus potential with the reaction force analysis provides two-parameter analytic expressions for the reaction force, reaction force constant, and reaction works. Moreover, since the parameters necessary to define the Marcus potential energy function can be obtained experimentally, the present model may produce experimental analytic potentials allowing for new and interesting applications, thus emerging as a powerful tool to characterize activation processes in chemical reactions.
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.
The mechanism of ground and excited state double proton transfer reaction in formic acid dimer has been analyzed with the help of reaction force and the reaction electronic flux. The separation of reaction electronic flux in terms of electronic activity and reactivity, NBO, and dual descriptor lends additional support for the mechanism. Interestingly we found that the ground state double proton transfer mechanism is concerted synchronic, whereas the excited state double proton transfer is concerted asynchronic in nature.
