Discovering trends in the Lewis acidity of beryllium and magnesium hydrides and fluorides with increasing clusters size.

We have reported in the last years the strong effect that Be- and Mg-containing Lewis acids have on the intrinsic properties of typical bases, which become acids upon complexation. In an effort to investigate these changes when the Be and Mg derivatives form clusters of increasing size, we have examined the behavior of the (MX2)n (M = Be, Mg; X = H, F; n = 1, 2, 3) clusters when they interact with ammonia, methanimine, hydrogen cyanide and pyridine, and with their corresponding deprotonated forms. The complexes obtained at the M06-2X/aug-cc-pVTZ level were analyzed using the MBIE energy decomposition formalism, in parallel with QTAIM, ELF, NCIPLOT and AdNDP analyses of their electron density. For n = 1 the interaction enthalpy for the different families of monomers, Be (Mg) hydrides and Be (Mg) fluorides, follows the same trend as the intrinsic basicity of the base that interacts with them. This interaction is greatly reinforced after the deprotonation of the base, resulting in a significant enhancement of the intrinsic acidity of the corresponding MX2-Base complex. For (MX2)2 clusters a further reinforcement of the interaction with the base is observed, this reinforcement being again larger for the deprotonated complexes. However, the concomitant increase of their intrinsic acidity is one order of magnitude larger for hydrides than for fluorides. Unexpectedly, the cyclic conformers (MX2)3, which are more unstable than the linear ones, become the global minima after association with the base and the same is true for the deprotonated complex. Accordingly, a further increase of the intrinsic acidity of the (MX2)3-Base complexes with respect to the (MX2)2-Base ones is observed. This effect is maximum for (MgF2)3 clusters, to the point that the (MgF2)3-Base complexes become more acidic than nitric acid, the extreme case being the cluster (MgF2)3-NCH, whose acidity is higher than that of perchloric acid.

Striking Borane Planarization in the Thermal Rearrangement (η5-C5H5)Fe(η3-B5H10)→(η5-C5H5)Fe(η5-B5H10).

In 1977 Weiss and Grimes, by means of mass spectrometry and 1H and 11B NMR spectroscopy, proposed two structures (I and II) for the ferraborane (η5-C5H5)Fe(B5H10), isoelectronic with ferrocene. In this work, by means of high-level quantum-chemical computations, we confirm the experimental structures of the two isomers with their corresponding energies, and assign the reported 1H and 11B NMR chemical shifts. A striking result from this study is the planarization (3D→2D) of the B5H10 - ligand - an unknown isolated anion, isoelectronic with aromatic cyclopentadienyl anion C5H5 - - when attached to the (η5-C5H5)Fe+ moiety, thus resulting in a more stable ferraborane isomer II.

An ab Initio Study of High-Energetic Trioxatriazinane Derivatives: Conformational Analysis, Nitrogen Inversion, Heats of Formation, Dissociation Reaction, and Dimerization.

High-energetic materials belong to two main classes: propellants and explosives. The still rather unexplored family of 1,3,5,2,4,6-trioxatriazinanes, N3O3R3, has a representative of each class. We have selected three compounds, R = H, R = CH3 and R = NO2, this last compound being known as TNTOTA, "trinitro-trioxa-triazinane". Of these compounds we have studied the conformational analysis, the nitrogen inversion, the heats of formation, and the dissociation reaction into the three monomers. In addition, the corresponding 1,3,2,4-dioxadiazetidines (N2O2R2) have also been studied.

Low valence triel (I) systems as hydrogen bond acceptors and their stability with respect to triel (III) compounds.

A theoretical study of the complexes formed by carbene like Al(I), Ga(I), In(I) and Tl(I) compounds with hydrogen bond donors (HBD), XH (HCCH, HSH ,HOH, HCN, HCl, HBr, HF, and HNC) have been carried out at MP2 computational level. The isolated triel(I) compounds show a negative region of the molecular electrostatic potential region associated with the triel atom suitable to interact with electron deficient groups. This region is associated to a lone pair based on the ELF analysis and to the location of the HOMO orbital. The complexes are similar to those found in nitrogen heterocyclic carbenes (NHC) with HBD. In addition, the oxidative addition reactions of those complexes to yield the corresponding valence III compounds have been characterized. The Al(III) compounds are much more stable than the corresponding Al(I) complexes. However, the stability of the triel(III) compounds decreases with the size of the triel atom and for the thallium derivatives, the Tl(I) complexes are more stable than the Tl(III) compounds in accordance with the number of the structures found in the CSD. The barrier of the TS connecting the triel(I) and triel(III) systems increases with the size of the triel atoms.

Beryllium as a base: complexes of Be(CO)3 with HX (X = F, Cl, Br, CN, NC, CCH, OH).

Beryllium chemistry is typically governed by its electron deficient character, but in some compounds it can act as a base. In order to understand better the unusual basicity of Be, we have systematically explored the complexes of one such compound, Be(CO)3, towards several hydrogen bond donors HX (X = F, Cl, Br, CN, NC, CCH, OH). For all complexes we find three different minima, two hydrogen bonded minima (to the Be or O atoms), and one weak beryllium bonded minimum. Further characterization of the interactions using a topological analysis of the electron density and Symmetry Adapted Perturbation Theory (SAPT) provide insight into the nature of these interactions. Overall these results highlight the capability of certain beryllium compounds to act as either a weak Lewis acid or, unconventionally, a Lewis base whose basicity towards hydrogen bonding is comparable to that of π systems.

Towards 2D Borane Chemistry in Hexagonal Cyclic Compounds.

A comprehensive comparison between known benzene mono-substituted compounds R-Ph and the corresponding isoelectronic unknown R-cyclohexaborane(12) molecules is carried out from a geometric and electronic structure point of view, with R={H, BH2, CH3, NH2, OH, F ; AlH2, SiH3, PH2, SH, Cl ; NO2, OCH3}. We suggest new chemical names for the 2D borane compounds and analyze the geometric and electronic structure carbon vs. boron comparatives by means of HOMO-LUMO gaps, bonding schemes, electron density topological properties and predicted NMR chemical shifts. The predictions on the properties in planar hexagonal cyclic boranes may help in the design of synthesis procedures for these yet-unkown compounds.

(Pyridin-2-ylmethyl)triel Derivatives as Masked Frustrated Lewis Pairs. Interactions and CO2 -Sequestration.

The isolated (pyridin-2-ylmethyl)triel derivatives (triel=B, Al and Ga) show an intramolecular N⋅⋅⋅Tr triel bond as shown by compounds found in the Cambridge Structural Database and DFT calculations. The possibility to use them as masked frustrated Lewis pairs (mFLP) has been explored theoretically concerning their reaction with CO2 . The adduct formation proceeds in two steps. In the first one, the (pyridin-2-ylmethyl)triel derivatives break the intramolecular N⋅⋅⋅Tr bond assisted by CO2 and in the second step the adduct is formed with Tr-O and N-C covalent bonds. The corresponding energy minima and transition states (TS) of the reaction have been characterized and analyzed. The distortion/interaction model analysis of the stationary points indicates that the whole process can be divided in two parts: reorganization of the mFLP in the first steps of the reaction while the reaction with CO2 (associated to the distortion of this molecule) is more important in the formation of the final adduct. In all cases studied, the final products are more stable than the starting molecules that combine with reasonable TS energies indicating that these reactions can occur.

A multi-FLP approach for CO2 capture: investigating nitrogen, boron, phosphorus and aluminium doped nanographenes and the influence of a sodium cation.

The reactivity of B3N3-doped hexa-cata-hexabenzocoronene (B3N3-NG), Al3N3-NG, B3P3-NG and Al3P3-NG, models of doped nanographenes (NGs), towards carbon dioxide was studied with density functional theory (DFT) calculations at the M06-2X/6-311++G(3df,3pd)//M06-2X/6-31+G* level of theory. The NG systems exhibit a poly-cyclic poly-frustrated Lewis pair (FLP) nature, featuring multiple Lewis acid/Lewis base pairs on their surface enabling the capture of several CO2 molecules. The capture of CO2 by these systems was investigated within two scenarios: (A) sequential capture of up to three CO2 molecules and (B) capture of CO2 molecules in the presence of a sodium cation. The resulting adducts were analyzed in terms of the activation barriers and relative stabilities. The presence of aluminium atoms changes the asynchrony of the reaction favoring the aluminium-oxygen bond and influences the regioselectivity of the multi-capture. A cooperative effect is predicted due to π-electron delocalization, with the sodium cation stabilizing the stationary points and favoring the addition of CO2 to the NGs.

Capture of CO2 by Melamine Derivatives: A DFT Study Combining the Relative Energy Gradient Method with an Interaction Energy Partitioning Scheme.

A theoretical study of the interaction between melamine and CO2 was carried out using density functional theory (DFT) with the B3LYP-D3(BJ)/aug-cc-pVTZ level of theory. The presence of anions interacting with melamine transforms the weakly bonded tetrel complexes into adducts. Thus, melamine acts as an FLP (frustrated Lewis pair) with acid groups (NHs as hydrogen bond donors) and a base group (N of the triazine ring). The application of the relative energy gradient formalism (REG) along the reaction coordinate has demonstrated that the ability of the melamine-anion systems to capture CO2 is linked to its capacity to polarize the CO2 molecule. These results have been confirmed by placing the melamine:CO2 complex in a uniform electric field with different strengths.

Dismantlement of ammonia upon interaction with Ben (n ≤ 10) clusters.

The interaction of ammonia with Ben (n < 1-10) clusters has been investigated by density functional theory and ab initio calculations. The main conclusion is that, regardless of the size of the Be cluster, neither the structure of ammonia nor that of the Be clusters are preserved due to a systematic dissociation of its NH bonds and a spontaneous H-shift toward the available Be atoms. This H migration not only leads to rather stable BeH bonds, but dramatically enhances the strength of the BeN bonds as well. Accordingly, the maximum stability is found for the interaction with the beryllium trimer, leading to a complex with three NBe and three BeH bonds. Another maximum in stability, although lower than that reached for n = 3, is found for the Be heptamer, since from n = 6, a new NBe bond is formed, so that complexes from n = 6 to n = 10 are characterized by the formation of a NBe4 moiety, whose stability reaches a maximum at n = 7. The bonding characteristics of the different species formed are analyzed by means of AIM, NBO, ELF and AdNDP approaches.

Production of Dihydrogen Using Ammonia Borane as Reagent and Pyrazole as Catalyst.

Theoretical chemistry (DLPNO-CCSD(T)/def2-TZVP//M06-2x/aug-cc-pVDZ) was used to design a system based on ammonia boranes catalyzed by pyrazoles with the aim of producing dihydrogen, nowadays of high interest as clean fuel. The reactivity of ammonia borane and cyclotriborazane were investigated, including catalytic activation through 1H-pyrazole, 4-methoxy-1H-pyrazole, and 4-nitro-1H-pyrazole. The results point toward a catalytic cycle by which, at the same time, ammonia borane can initially store and then, through catalysis, produce dihydrogen and amino borane. Subsequently, amino borane can trimerize to form cyclotriborazane that, in presence of the same catalyst, can also produce dihydrogen. This study proposes therefore a consistent progress in using environmentally sustainable (metal free) catalysts to efficiently extract dihydrogen from small B-N bonded molecules.

Metastable Charged Dimers in Organometallic Species: A Look into Hydrogen Bonding between Metallocene Derivatives.

Multiply charged complexes bound by noncovalent interactions have been previously described in the literature, although they were mostly focused on organic and main group inorganic systems. In this work, we show that similar complexes can also be found for organometallic systems containing transition metals and deepen in the reasons behind the existence of these species. We have studied the structures, binding energies, and dissociation profiles in the gas phase of a series of charged hydrogen-bonded dimers of metallocene (Ru, Co, Rh, and Mn) derivatives isoelectronic with the ferrocene dimer. Our results indicate that the carboxylic acid-containing dimers are more strongly bonded and present larger barriers to dissociation than the amide ones and that the cationic complexes tend to be more stable than the anionic ones. Additionally, we describe for the first time the symmetric proton transfer that can occur while in the metastable phase. Finally, we use a density-based energy decomposition analysis to shine light on the nature of the interaction between the dimers.

A theoretical study of the reaction of borata derivatives of benzene, anthracene and pentacene with CO2.

A theoretical study of the reaction between several borataacenes (1-methylboratabenzene, 9-methyl-9-borataanthracene and cis and trans diboratapentacene) and CO2 has been carried out at the M06-2X computational level. The influence of a counterion (potassium cation), the cation complexation by 18-crown-6-ether and solvent effects have been explored. The computational results predict anti/syn selectivity as found experimentally in the cis- and trans-diboratapentacene reaction with CO2 (Baker et al., J. Am. Chem. Soc., 2023, 145, 2028).

Reply to the 'Comment on "Theoretical study of the NO3 radical reaction with CH2ClBr, CH2ICl, CH2BrI, CHCl2Br, and CHClBr2"' by C. J. Nielsen and Y. Tang, Phys. Chem. Chem. Phys., 2022, 24, DOI: 10.1039/D2CP03013F.

In this Reply, we answer the main argument raised in the Comment about the energy of the NO3 radical and its influence in the reaction profiles of the reaction of the NO3 radical with CH2ClBr, CH2ICl, CH2BrI, CHCl2Br, and CHClBr2 by C. J. Nielsen and Y. Tang. The optimized geometry of the NO3 radical has been obtained using 49 DFT functionals: 26 functionals predict a minimum with D3h symmetry and 23 with C2v symmetry. The former functionals have been used to calculate the thermodynamic values of three reactions (X + HNO3 → XH + NO3, X= OH, CH3 and CCl3) and compared with experimental data. Those functionals with smaller errors have been used to recalculate the barriers of the reaction of NO3 with CH2ClBr, CH2ICl, CH2BrI, CHCl2Br, and CHClBr2. The results show differences of 10.5 kJ mol-1 when compared to those obtained with the M08HX functional.

Hydrogen-Bond Dissociation Energies from the Properties of Isolated Monomers.

The strength of binding, as measured by the equilibrium dissociation energy De of an isolated hydrogen-bonded complex B···HX, where B is a simple Lewis base and X = F, Cl, Br, I, CN, CCH, or CP, can be determined from the properties of the infinitely separated components B and HX. The properties in question are the maximum and minimum values σmax(HX) and σmin(B) of the molecular electrostatic surface potentials on the 0.001 e/bohr3 iso-surfaces of HX and B, respectively, and two recently defined quantities: the reduced electrophilicity ΞHX of HX and the reduced nucleophilicity ИB of B. It is shown that De is given by the expression De = {σmax(HX)σmin(B)} ИB ΞHX. This is tested by comparing De calculated ab initio at the CCSD(T)(F12c)/cc-pVDZ-F12 level of theory with that obtained from the equation. A large number of complexes (203) falling into four categories involving different types of hydrogen-bonded complex B···HX are investigated: those in which the hydrogen-bond acceptor atom of B is either oxygen or nitrogen, or carbon or boron. The comparison reveals that the proposed equation leads to De values in good agreement in general with those calculated ab initio.

Linking the Interatomic Exchange-Correlation Energy to Experimental J-Coupling Constants.

The main aim of the current work is to find an experimental connection to the interatomic exchange-correlation energy as defined by the energy decomposition method Interacting Quantum Atoms (IQA). A suitable candidate as (essentially) experimental quantity is the nuclear magnetic resonance (NMR) J-coupling constant denoted 3J(H,H'), which a number of previous studies showed to correlate well with QTAIM's delocalization index (DI), which is essentially a bond order. Inspired by Karplus equations, here, we investigate correlations between 3J(H,H') and a relevant dihedral angle in six simple initial compounds of the shape H3C-YHn (Y = C, N, O, Si, P, and S), N-methylacetamide (as prototype of the peptide bond), and five peptide-capped amino acids (Gly, Ala, Val, Ile, and Leu) because of the protein direction of the force field FFLUX. In conclusion, except for methanol, the inter-hydrogen exchange-correlation energy Vxc(H,H') makes the best contact with experiment, through 3J(H,H'), when multiplied with the internuclear distance RHH'.

Dispersion, Rehybridization, and Pentacoordination: Keys to Understand Clustering of Boron and Aluminum Hydrides and Halides.

The structure, stability, and bonding characteristics of dimers and trimers involving BX3 and AlX3 (X = H, F, Cl) in the gas phase, many of them explored for the first time, were investigated using different DFT (B3LYP, B3LYP/D3BJ, and M06-2X) and ab initio (MP2 and G4) methods together with different energy decomposition formalisms, namely, many-body interaction-energy and localized molecular orbital energy decomposition analysis. The electron density of the clusters investigated was analyzed with QTAIM, electron localization function, NCIPLOT, and adaptive natural density partitioning approaches. Our results for triel hydride dimers and Al2X6 (X = F, Cl) clusters are in good agreement with previous studies in the literature, but in contrast with the general accepted idea that B2F6 and B2Cl6 do not exist, we have found that they are predicted to be weakly bound systems if dispersion interactions are conveniently accounted for in the theoretical schemes used. Dispersion interactions are also dominant in both homo- and heterotrimers involving boron halide monomers. Surprisingly, B3F9 and B3Cl9 C3v cyclic trimers, in spite of exhibiting rather strong B-X (X = F, Cl) interactions, were found to be unstable with respect to the isolated monomers due to the high energetic cost of the rehybridization of the B atom, which is larger than the two- and three-body stabilization contributions when the cyclic is formed. Another important feature is the enhanced stability of both homo- and heterotrimers in which Al is the central atom because Al is systematically pentacoordinated, whereas this is not the case when the central atom is B, which is only tri- or tetra-coordinated.

A Holistic View of the Interactions between Electron-Deficient Systems: Clustering of Beryllium and Magnesium Hydrides and Halides.

In the search for common bonding patterns in pure and mixed clusters of beryllium and magnesium derivatives, the most stable dimers and trimers involving BeX2 and MgX2 (X = H, F, Cl) have been studied in the gas phase using B3LYP and M06-2X DFT methods and the G4 ab initio composite procedure. To obtain some insight into their structure, stability, and bonding characteristics, we have used two different energy decomposition formalisms, namely MBIE and LMO-EDA, in parallel with the analysis of the electron density with the help of QTAIM, ELF, NCIPLOT, and AdNDP approaches. Some interesting differences are already observed in the dimers, where the stability sequence observed for the hydrides differs entirely from that of the fluorides and chlorides. Trimers also show some peculiarities associated with the presence of compact trigonal cyclic structures that compete in stability with the more conventional hexagonal and linear forms. As observed for dimers, the stability of the trimers changes significantly from hydrides to fluorides or chlorides. Although some of these clusters were previously explored in the literature, the novelty of this work is to provide a holistic approach to the entire series of compounds by using chemical bonding tools, allowing us to understand the stability trends in detail and providing insights for a significant number of new, unexplored structures.

Use of 5,10-Disubstituted Dibenzoazaborines and Dibenzophosphaborines as Cyclic Supports of Frustrated Lewis Pairs for the Capture of CO2.

The reactivity of 5,10-disubstituted dibenzoazaborines and dibenzophosphaborines towards carbon dioxide was studied at the DFT, M06-2X/def2-TZVP, computational level. The profile of this reaction comprises of three stationary points: the pre-reactive complex and adduct minima and the transition state(TS) linking both minima. Initial results show that dibenzoazaborines derivatives are less suitable to form adducts with CO2 than dibenzophosphaborine systems. The influence of the basicity on the P atom and the acidity on the B center of the dibenzophosphaborine in the reaction with CO2 was also explored. Thus, an equation was developed relating the properties (acidity, basicity and boron hybridization) of the isolated dibenzophosphaborine derivatives with the adduct energy. We found that modulation of the boron acidity allows to obtain more stable adducts than the pre-reactive complexes and isolated monomers.

Stand up for Electrostatics: The Disiloxane Case.

The basicity of the simplest silicone, disiloxane (H3 Si-O-SiH3 ), is strongly affected by the Si-O-Si angle (α). We use high-level ab initio MP2/aug'-cc-pVTZ calculations and the molecular electrostatic potential (MEP) to analyze the relationship between the increase in basicity and the reduction of α. Our results clearly point out that this increase can be explained through the MEP, as the interactions between oxygen from disiloxane and the acceptors are mostly electrostatic. Furthermore, the effect of α on the tetrel bond between disiloxane and several Lewis bases can again be rationalized using the MEP. Finally, we explore the cooperativity throughout α for ternary complexes where disiloxane simultaneously interacts with a Lewis acid and a Lewis base. Both non-covalent interactions remain cooperative for all α values, although the largest cooperativity effects are not always those maximizing the binding energy in the binary complexes. Overall, the MEP remains a powerful predictor for noncovalent interactions.