Trapping CO2 by Adduct Formation with Nitrogen Heterocyclic Carbenes (NHCs): A Theoretical Study

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.

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.

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.

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

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

The reactivity of a B₃P₃-doped hexa-cata-hexabenzocoronene, as a model of nanographene (B₃P₃-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 CO₂ by this multi-FLP system: (i) fixation of three CO₂ molecules sequentially one by one; and (ii) simultaneous contact of three CO₂ molecules with the B₃P₃-NG surface. The resulting adducts were analyzed as function of activation barriers and the relative stability of the CO₂ capture. A cooperativity effect due to the π-delocalization of the hexa-cata-hexabenzocoronene is observed. The fixation of a CO₂ molecule modifies the electronic properties. It enhances the capture of additional CO₂ molecules by changing the acidy and basicity of the rest of the boron and phosphorus atoms in the B₃P₃-NG system.

An in silico study of the selective adsorption and separation of CO2 from a flue gas mixture (CH4, CO2, N2) by ZnLi5+ clusters

Due to the increasing concentration of CO₂ in the atmosphere and its negative effect on the environment, selective adsorption of CO₂ from flue gas has become significantly important. In this study, we have considered a Zn-doped lithium cluster, ZnLi₅⁺ cluster, featuring a planar pentacoordinate Zn centre, as a potential candidate for selective CO₂ capture and separation from a flue gas mixture (CH₄, CO₂, N₂). The binding energy calculation and non-covalent interaction study showed that CO₂ molecules bind relatively strongly as compared to N₂ and CH₄ molecules. The metal cluster can bind five CO₂, five CH₄, and four N₂ molecules with average binding energies of -9.2, -4.4, and -6.1 kcal mol^-1, respectively. Decomposition of the binding energy through symmetry-adapted perturbation theory analysis reveals that the electrostatic component plays a major role. The cationic cluster may be a promising candidate for selective CO₂ capture and can be used as a pollution-controlling agent. The calculated adsorption energy of H₂S is quite closer to that of CO₂, suggesting competitive adsorption between CO₂ and H₂S. The adsorption energies of H₂O and NH₃ are higher compared to CO₂, indicating that these gases may be a potential threat to CO₂ capture.

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.

Interacting Quantum Atoms Analysis of the Reaction Force: A Tool to Analyze Driving and Retarding Forces in Chemical Reactions

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.

pKa Scale for Cyclopropenium Ions with Applications in CO2 Capture

Molecular acid-base properties are core to understanding chemical systems and the prediction of reactivity. This axiom holds for cyclopropenium ions in terms of their broad use as (organo)catalysts, ligands, redox-flow batteries, and applications in materials sciences. In view of this significant status, and with it, the critical importance of acidity, we disclose in this report the first comprehensive computational study of the pK_a values of cyclopropenium ions employing a subset of 70 structurally diverse cyclopropenium derivatives, density functional computations, and Hammett linear free-energy relationships. Capitalizing upon these computed findings, and with an eye toward greenhouse gas trapping, we further document the timely use of a cyclopropenium-cyclopropenylidene coupled platform for CO₂ capture and light-triggered release.

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.

Theoretical study on the noncovalent interactions involving triplet diphenylcarbene

The properties of some types of noncovalent interactions formed by triplet diphenylcarbene (DPC³) have been investigated by means of density functional theory (DFT) calculations and quantum theory of atoms in molecule (QTAIM) studies. The DPC³···LA (LA = AlF₃, SiF₄, PF₅, SF₂, ClF) complexes have been analyzed from their equilibrium geometries, binding energies, and properties of electron density. The triel bond in the DPC³···AlF₃ complex exhibits a partially covalent nature, with the binding energy - 65.7 kJ/mol. The tetrel bond, pnicogen bond, chalcogen bond, and halogen bond in the DPC³···LA (LA = SiF₄, PF₅, SF₂, ClF) complexes show the character of a weak closed-shell noncovalent interaction. Polarization plays an important role in the formation of the studied complexes. The strength of intermolecular interaction decreases in the order LA = AlF₃ > ClF > SF₂ > SiF₄ > PF₅. The electron spin density transfers from the radical DPC³ to ClF and SF₂ in the formation of halogen bond and chalcogen bond, but for the DPC³···AlF₃/SiF₄/PF₅ complexes, the transfer of electron spin density is minimal.

Competitive Reactivity of SO2 and NO2 with N-Heterocyclic Carbene: A Mechanistic Study

Recent DFT based molecular engineering to obtain stable oxathiirane S-oxide derivatives evokes the recommencement of the use of carbenes for the sequestering of SO₂, which has been kept separate so far. Carbene is one of the key chemicals for the sequestering of various premier greenhouse gases like CO₂, CO, N₂O, etc. In this respect, a comparative study of the reactivity of carbenes with variant greenhouse gases is highly demanding. The present investigation is engrossed in the comparative reactivity of SO₂ and NO₂ with carbenes. All three selected carbenes are highly susceptible to SO₂ and NO₂. Through an immaculate mechanistic study, we are able to corroborate that the end product of the carbene-SO₂ reaction is an adduct which has a preferable structure having a six-membered ring with hydrogen bonding instead of ketone and SO with higher thermodynamic stability than the corresponding oxathiirane S-oxide derivative. Carbene reacts with NO₂ to form a stable carbene N, N-dioxide derivative which forms vibrationally excited oxaziridine N-oxide which rapidly dissociates to form a ketone derivative. The formation of carbene S, S-dioxide and carbene N, N-dioxide is a barrierless process. The dissociation of oxaziridene N-oxide is also a barrierless process.

Interacting Quantum Atoms-A Review

The aim of this review is threefold. On the one hand, we intend it to serve as a gentle introduction to the Interacting Quantum Atoms (IQA) methodology for those unfamiliar with it. Second, we expect it to act as an up-to-date reference of recent developments related to IQA. Finally, we want it to highlight a non-exhaustive, yet representative set of showcase examples about how to use IQA to shed light in different chemical problems. To accomplish this, we start by providing a brief context to justify the development of IQA as a real space alternative to other existent energy partition schemes of the non-relativistic energy of molecules. We then introduce a self-contained algebraic derivation of the methodological IQA ecosystem as well as an overview of how these formulations vary with the level of theory employed to obtain the molecular wavefunction upon which the IQA procedure relies. Finally, we review the several applications of IQA as examined by different research groups worldwide to investigate a wide variety of chemical problems.

Hydrogen bonds and halogen bonds in complexes of carbones L→C←L as electron donors to HF and ClF, for L = CO, N2, HNC, PH3, and SH2

MP2 and EOM-CCSD calculations have been carried out to determine the structures, binding energies, and spin-spin coupling constants of carbone complexes L→C←L with the carbone the electron donor to HF or ClF, for L = CO, N₂, HNC, PH₃, and SH₂.

Formation of Carbon Dioxide Attached Fragment Ions in the Fragmentation of Deprotonated Tolfenpyrad and Tebufenpyrad

The in-source collision-induced dissociation (CID) and MS/MS mass spectra of deprotonated tolfenpyrad and tebufenpyrad both showed an unusual fragment ion at m/z 187, but its fragmentation pattern and structure could not be explained by logical neutral losses. Accurate mass measurement indicated that the mass difference between this fragment ion and the dominant fragment ion at m/z 143 equaled to a carbon dioxide (CO₂) molecule. The isolation of the fragment ion m/z 143 in the mass analyzer could spontaneously give rise to the ion m/z 187. The Gibbs free energy of carbon dioxide addition to deprotonated pyrazole ion was significantly negative from the computational results. According to these results, we derived a proposal for the formation and structure of the ion m/z 187, which was an attachment of molecular carbon dioxide to the fragment ion m/z 143 to produce a carboxylate anion. The trace carbon dioxide was speculated to be derived from the residual atmosphere or collision gas in the instrument. This study is valuable for the qualitative and quantitative mass spectrometry analysis of pesticides containing the pyrazole functional group.

Lennard-Jones Intermolecular Potentials for the Description of 6-Membered Aromatic Heterocycles Interacting with the Isoelectronic CO2 and CS2

We have generated Lennard-Jones potentials for the interaction between CX₂ (X = O, S) and 11 nitrogen-doped benzene derivatives in different orientations at the M06-2X/def2-tzvpp level as tools to parametrize accurate force fields and to better understand the interaction of these greenhouse gases with heterocyclic building blocks used in the design of capture and detection systems. We find that the most favorable interactions are found between the carbon in CO₂ and the main heterocycle in the ring in a parallel orientation, whereas the preferred interaction mode of CS₂ is established between sulfur and the π density of the aromatic ring. The fact that the preferences for interaction sites and orientations of CO₂ and CS₂ are most of the times opposite helps in terms of ensuring the selectivity of these systems in front of these two isoelectronic compounds. The existence of very good linear correlations ( R² values very close to one) between the number of nitrogen atoms in the heterocyclic ring and the depth of the interaction potential wells opens the door to the use of these results in generating coarse-grained potentials or models with predictive power for use in the design of larger systems.

Comparison of pnicogen and tetrel bonds in complexes containing CX2 carbenes (X = F, Cl, Br, OH, OMe, NH2, and NMe2)

The similarities and differences of pnicogen and tetrel bonds formed by carbenes CX₂ with H₃AsO and H₃SiCN were investigated by carrying out ab initio calculations in association with topological analysis of electron density.

A relative energy gradient (REG) study of the planar and perpendicular torsional energy barriers in biphenyl

Biphenyl is a prototype molecule, the study of which is important for a proper understanding of stereo-electronic effects. In the gas phase it has an equilibrium central torsion angle of ~ 45° and shows both a planar (0°) and a perpendicular (90°) torsional energy barrier. The latter is analysed for the first time. We use the newly proposed REG method, which is an exhaustive procedure that automatically ranks atomic energy contributions according to their importance in explaining the energy profile of a total system. Here, the REG method operates on energy contributions computed by the interacting quantum atoms method. This method is minimal in architecture and provides a crisp picture of well-defined and well-separated electrostatic, steric and exchange (covalent) energies at atomistic level. It is shown that the bond critical point occurring between the ortho-hydrogens in the planar geometry has been wrongly interpreted as a sign of repulsive interaction. A convenient metaphor of analysing football matches is introduced to clarify the role of a REG analysis.

Solvent and Substituent Effects on the Phosphine + CO2 Reaction

A theoretical study of the substituent and solvent effects on the reaction of phosphines with CO2 has been carried out by means of Møller-Plesset (MP2) computational level calculations and continuum polarizable method (PCM) solvent models. Three stationary points along the reaction coordinate have been characterized, a pre-transition state (TS) assembly in which a pnicogen bond or tetrel bond is established between the phosphine and the CO2 molecule, followed by a transition state, and leading finally to the adduct in which the P–C bond has been formed. The solvent effects on the stability and geometry of the stationary points are different. Thus, the pnicogen bonded complexes are destabilized as the dielectric constant of the solvent increases while the opposite happens within the adducts with the P–C bond and the TSs trend. A combination of the substituents and solvents can be used to control the most stable minimum.

Using the Relative Energy Gradient Method with Interacting Quantum Atoms to Determine the Reaction Mechanism and Catalytic Effects in the Peptide Hydrolysis in HIV-1 Protease

The reaction mechanism in an active site is of the utmost importance when trying to understand the role that an enzyme plays in biological processes. In a recently published paper [Theor. Chem. Acc. 2017, 136, 86], we formalised the Relative Energy Gradient (REG) method for automating an Interacting Quantum Atoms (IQA) analysis. Here, the REG method is utilised to determine the mechanism of peptide hydrolysis in the aspartic active site of the enzyme HIV-1 Protease. Using the REG method along with the IQA approach we determine the mechanism of peptide hydrolysis without employing any arbitrary parameters and with remarkable ease (albeit at large computational cost: the system contains 133 atoms, which means that there are 17 689 individual IQA terms to be calculated). When REG and IQA work together it is possible to determine a reaction mechanism at atomistic resolution from data directly derived from quantum calculations, without arbitrary parameters. Moreover, the mechanism determined by this novel method gives concrete insight into how the active site residues catalyse peptide hydrolysis.

Complexes of O=C=S with Nitrogen Bases: Chalcogen Bonds, Tetrel Bonds, and Other Secondary Interactions

Ab initio MP2/aug'-cc-pVTZ calculations have been carried out to investigate chalcogen-bond formation through the σ-hole at S and tetrel-bond formation through the π-hole at C in complexes of OCS with a series of nitrogen bases. The binding energies of chalcogen- and tetrel-bonded complexes with the sp-hybridized bases correlate exponentially with the N-S and N-C distances, respectively. The presence of secondary interactions between an N-H or C-H group of an sp2 -hybridized base and OCS in chalcogen-bonded complexes decreases the correlation between binding energies and the N-S distance. These secondary interactions are stronger in the tetrel-bonded complexes with the sp2 bases, particularly in the isomers of OCS:imidazole and OCS : N2 H2 , where they may be described as distorted N-H⋅⋅⋅O or N-H⋅⋅⋅S hydrogen bonds. Charge-transfer interactions are consistent with the nature of the primary and secondary interactions in these complexes. The in-plane OCS bending frequencies are blue-shift in the chalcogen-bonded complexes, and red-shifted in the tetrel-bonded complexes. EOM-CCSD spin-spin coupling constants 1c J(N4-S) across chalcogen bonds have absolute values less than 9.0 Hz, while the two-bond coupling constants 2c J(N4-C) do not exceed 4.0 Hz. These are greater in absolute value that the one-bond coupling constants 1t J(N4-C) across tetrel bonds that are less than 0.5 Hz at much shorter N-C distances.

Fluorine Gauche Effect Explained by Electrostatic Polarization Instead of Hyperconjugation: An Interacting Quantum Atoms (IQA) and Relative Energy Gradient (REG) Study

We present an interacting quantum atoms (IQA) study of the gauche effect by comparing 1,2-difluoroethane, 1,2-dichloroethane, and three conformers of 1,2,3,4,5,6-hexafluorocyclohexane. In the 1,2-difluoroethane, the gauche effect is observed in that the gauche conformation is more stable than the anti, whereas in 1,2-dichloroethane the opposite is true. The analysis performed here is exhaustive and unbiased thanks to using the recently introduced relative energy gradient (REG) method [ Thacker , J. C. R. ; Popelier , P. L. A. Theor. Chem. Acc . 2017 , 136 , 86 ], as implemented in the in-house program ANANKE. We challenge the common explanation that hyperconjugation is responsible for the gauche stability in 1,2-difluoroethane and instead present electrostatics as the cause of gauche stability. Our explanation of the gauche effect is also is seen in other molecules displaying local gauche conformations, such as the recently synthesized "all-cis" hexafluorocyclohexane and its conformers where all the fluorine atoms are in the equatorial positions. Using our extension of the traditional IQA methodology that allows for the partitioning of electrostatic terms into polarization and charge transfer, we propose that the cause of gauche stability is 1,3 C···F electrostatic polarization interactions. In other words, if a number of fluorine atoms are aligned, then the stability due to polarization of nearby carbon atoms is increased.

Binding indirect greenhouse gases OCS and CS2by nitrogen heterocyclic carbenes (NHCs)

Carbon disulfide (CS₂) and carbonyl sulfide (OCS) are indirect greenhouse gases that can be effectively trapped by classical, abnormal and remote nitrogen heterocyclic carbenes (NHCs), according to high levelab initiocalculations.

The ANANKE relative energy gradient (REG) method to automate IQA analysis over configurational change

Much chemical insight ultimately comes down to finding out which fragment of a total system behaves like the total system, in terms of an energy profile. A simple example is that of the water dimer, where this system is regarded as held together by a hydrogen bond. The hydrogen bond consists of two atoms (H···O), which energetically behave similarly to the total system (H₂O)₂. However, from a quantum mechanical point of view, each atom in the total system interacts with any other atom. Thus, the view that the hydrogen bond by itself governs the energetic stability of the water dimer needs rigorous justification. In this work, we propose a method that provides such a justification, in general, but only illustrated on the water dimer here. This method is based on the topological energy partitioning method called interacting quantum atoms (IQA). The method is implemented in the program ANANKE, which calculates correlations between the energy profile of the total system and those of subsystems (or fragments). ANANKE acts on the IQA energy contributions obtained for a sequence of full-system geometries controlled by a coordinate of interest (e.g. the O···H distance in the water dimer). Although applied only for the water dimer in this work, the method is general and able to explain the gauche effect, the torsional barrier in biphenyl, the arrow-pushing scheme of an enzymatic reaction (peptide hydrolysis in the HIV-1 Protease active site), and halogen-alkane nucleophilic substitution (S_N2) reactions. Those applications will appear elsewhere as separate and elaborated case studies; here we focus on the details of the ANANKE method and its justification, using the water dimer as a concrete case.