Understanding the coupling of non-metallic heteroatoms to CO2 from a Conceptual DFT perspective.

CONTEXT: A Conceptual DFT (CDFT) study has been carry out to analyse the coupling reactions of the simplest amine (CH3NH2), alcohol (CH3OH), and thiol (CH3SH) compounds with CO2 to form the corresponding adducts CH3NHCO2H, CH3OCO2H, and CH3SCO2H. The reaction mechanism takes place in a single step comprising two chemical events: nucleophilic attack of the non-metallic heteroatoms to CO2 followed by hydrogen atom transfer (HAT). According to our calculations, the participation of an additional nucleophilic molecule as HAT assistant entails important decreases in activation electronic energies. In such cases, the formation of a six-membered ring in the transition state (TS) reduces the angular stress with respect to the non-assisted paths, characterised by four-membered ring TSs. Through the analysis of the energy and reaction force profiles along the intrinsic reaction coordinate (IRC), the ratio of structural reorganisation and electronic rearrangement for both activation and relaxation energies has been computed. In addition, the analysis of the electronic chemical potential and reaction electronic flux profiles confirms that the highest electronic activity as well as their changes take place in the TS region. Finally, the distortion/interaction model using an energy decomposition scheme based on the electron density along the reaction coordinate has been carried out and the relative energy gradient (REG) method has been applied to identify the most important components associated to the barriers. METHODS: The theoretical calculation were performed with Gaussian-16 scientific program. The B3LYP-D3(BJ)/aug-cc-pVDZ level was used for optimization of the minima and TSs. IRC calculations has also been carried out connecting the TS with the associated minima. Conceptual-DFT (CDFT) calculations have been carried out with the Eyringpy program and in-house code. The distortion/interaction model along the reaction coordinate have used the decomposition scheme of Mandado et al. and the analysis of the importance of each components have been done with the relative energy gradient (REG) method.

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.

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.

(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.

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.

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).

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.

Two Shared Icosahedral Metallacarboranes through Iron: A Joint Experimental and Theoretical Refinement of Mössbauer Spectrum in [Fe(1,2-C2B9H11)2]Cs.

Mössbauer and X-ray photoelectron spectroscopies (XPS) are complemented with high-level quantum-chemical computations in the study of the geometric and electronic structure of the paramagnetic salt of the metallacarborane sandwich complex [Fe(1,2-C2B9H11)2]Cs = FeSanCs. Experimental 57Fe isomer shifts and quadrupole splitting parameters are compared with the theoretical prediction, with good agreement. The appearance of two sets of Cs(3d) doublets in the XPS spectrum, separated by 2 eV, indicates that Cs has two different chemical environments due to ease of the Cs(+) cation moving around the sandwich complex with low-energy barriers, as confirmed by quantum-chemical computations. Several minimum-energy geometries of the FeSanCs structure with the corresponding energies and Mössbauer parameters are discussed, in particular the atomic charges and spin population and the surroundings of the Fe atom in the complex. The Mössbauer spectra were taken at different temperatures showing the presence of a low-spin Fe atom with S = 1/2 and thus confirming a paramagnetic FeIII species.

Reactivity of a model of B3P3-doped nanographene with up to three CO2 molecules.

The reactivity of a B3P3-doped hexa-cata-hexabenzocoronene, as a model of nanographene (B3P3-NG), towards carbon dioxide was studied at the DFT M06-2X/6-311++G(3df,3pd)//M06-2X/6-31+G* level of theory. This compound can be classified as a poly-cyclic poly-Frustrated Lewis Pair (FLP) system, as it presents more than one Lewis Acid/Lewis Base pair on its surface, making the capture of several carbon dioxide molecules possible. Two scenarios were considered to fully characterize the capture of CO2 by this multi-FLP system: (i) fixation of three CO2 molecules sequentially one by one; and (ii) simultaneous contact of three CO2 molecules with the B3P3-NG surface. The resulting adducts were analyzed as function of activation barriers and the relative stability of the CO2 capture. A cooperativity effect due to the π-delocalization of the hexa-cata-hexabenzocoronene is observed. The fixation of a CO2 molecule modifies the electronic properties. It enhances the capture of additional CO2 molecules by changing the acidy and basicity of the rest of the boron and phosphorus atoms in the B3P3-NG system.

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.

Sequestration of Carbon Dioxide with Frustrated Lewis Pairs Based on N-Heterocycles with Silane/Germane Groups.

Frustrated Lewis pairs (FLPs) based on nitrogen heterocycles (pyridine, pyrazole, and imidazole) with a silane or germane group in the α-position of a nitrogen atom have been considered as potential molecules to sequestrate carbon dioxide. Three stationary points have been characterized in the reaction profile: a pre-reactive complex, an adduct minimum, and the transition state connecting them. The effect of external (solvent) or internal (hydroxyl group) electric fields in the reaction profile has been considered. In both cases, it is possible to improve the kinetics and thermodynamics of the complexation of CO2 by the FLP and favor the formation of adducts.