The triel bond: a potential force for tuning anion–π interactions

Can H2 be Superacidic? A Computational Study of Triel-Bonded Brønsted Acids

The abundance of XIII group element compounds in science and industry together with their electron-deficient character gives rise to their influence on properties of the systems they interact with. This paper is an attempt to assess the strength, nature, and effect of formation of a triel bond on acidity. A wide set of Brønsted acids among others comprising hydrocarbons, halogen hydrides, and amines bonded with B, Al, and Ga trifluorides forming HX/TF3 was selected for the research. Various computational approaches (e.g., MP2, GFN2-xTB, SAPT2 + 3(CCD)δMP2, quantum theory of atoms in molecules analysis, and density overlap regions indicator) are used to describe the triel-bonded systems. Among other things, it was found that the electrostatics may not be the dominant contribution to the triel binding in some cases. Additionally, it was established that even weak Brønsted acids such as C2H2 or H2 may be superacidic if bonded to a Lewis acid (TF3) that is strong enough. The calculations indicate a significant covalent character of some of the studied HX/TF3 triel-bonded systems. Moreover, the effect of solvation of HX with TF3 as well as that of the reverse process on the acidity of the resulting system is thoroughly described.

Search for Osme Bonds with π Systems as Electron Donors

The Osme bond is defined as pairing a Group 8 metal atom as an electron acceptor in a noncovalent interaction with a nucleophile. DFT calculations with the ωB97XD functional consider MO4 (M = Ru, Os) as the Lewis acid, paired with a series of π electron donors C2H2, C2H4, C6H6, C4H5N, C4H4O, and C4H4S. The calculations establish interaction energies in the range between 9.5 and 26.4 kJ/mol. Os engages in stronger interactions than does Ru, and those involving more extensive π-systems within the aromatic rings form stronger bonds than do the smaller ethylene and acetylene. Extensive analysis questions the existence of a true Osme bond, as the bonding chiefly involves interactions with the three O atoms of MO4 that lie closest to the π-system, via π(C-C)→σ*(M-O) transfers. These interactions are supplemented by back donation from M-O bonds to the π*(CC) antibonding orbitals of the π-systems. Dispersion makes a large contribution to these interactions, higher than electrostatics and much greater than induction.

Ga···C Triel Bonds-Why They Are Not Strong Enough to Change Trigonal Configuration into Tetrahedral One: DFT Calculations on Dimers That Occur in Crystal Structures

Structures characterized by the trigonal coordination of the gallium center that interacts with electron rich carbon sites are described. These interactions may be classified as Ga···C triel bonds. Their properties are analyzed in this study since these interactions may be important in numerous chemical processes including catalytical activities; additionally, geometrical parameters of corresponding species are described. The Ga···C triel bonds discussed here, categorized also as the π-hole bonds, do not change the trigonal configuration of the gallium center into the tetrahedral one despite total interactions in dimers being strong; however, the main contribution to the stabilization of corresponding structures comes from the electrostatic forces. The systems analyzed theoretically here come from crystal structures since the Cambridge Structural Database, CSD, search was performed to find structures where the gallium center linked to CC bonds of Lewis base units occurs. The majority structures found in CSD are characterized by parallel, stacking-like arrangements of species containing the Ga-centers. The theoretical results show that interactions within dimers are not classified as the three-centers links as in a case of typical hydrogen bonds and numerous other interactions. The total interactions in dimers analyzed here consist of several local intermolecular atom-atom interactions; these are mainly the Ga···C links. The DFT results are supported in this study by calculations with the use of the quantum theory of atoms in molecules, QTAIM, the natural bond orbital, NBO, and the energy decomposition analysis, EDA, approaches.

Static and Dynamical Quantum Studies of CX3-AlX2 and CSiX3-BX2 (X = F, Cl, Br) Complexes with Hydrocyanic Acid: Unusual Behavior of Strong π-Hole at Triel Center

The set of TX₃-TrX₂ (T = C, Si, Ge; Tr = B, Al, Ga; X = F, Cl, Br) molecules offers a rather unique opportunity to study both σ-hole and π-hole dimerization on the tetrel and triel ends, respectively. According to the molecular electrostatic potential (MEP) distribution, the π-hole extrema (acidic sites) were more intense than their σ-hole counterparts. The molecules owning the most (CX₃-AlX₂) and least (SiX₃-BX₂) intense π-holes were chosen to evaluate their capacities to attract one and two HCN molecules (Lewis bases). We discovered that the energetic characteristics of π-hole dimers severely conflict with the monomers MEP pattern since the weakest π-hole monomer forms a dimer characterized by interaction energy compared to those created by the monomers with noticeably greater power in the π-hole region. This outcome is due to the deformation of the weakest π-hole donor. Furthermore, the MEP analysis for monomers in the geometry of respective dimers revealed a "residual π-hole" site that was able to drive second ligand attachment, giving rise to the two "unusual trimers" examined further by the NCI and QTAIM analyses. Apart from them, the π-hole/π-hole and σ-hole/π-hole trimers have also been obtained throughout this study and described using energetic and geometric parameters. The SAPT approach revealed details of the bonding in one of the "unusual trimers". Finally, Born-Oppenheimer Molecular Dynamics (BOMD) simulations were carried out to investigate the time evolution of the interatomic distances of the studied complexes as well as their stability.

Structural, vibrational and electronic properties of some tetrel-bonded complexes of the fluorinated methanes methyl fluoride, difluoromethane and fluoroform: an ab initio study

A search has been conducted, by means of ab initio molecular orbital theory, for potential tetrel-bonded complexes formed between the fluorinated methanes methyl fluoride, difluoromethane and fluoroform, and the related hydrides ammonia, water, hydrogen fluoride, phosphine, hydrogen sulphide and hydrogen chloride. Eleven such complexes have been identified, six containing CH₃F and five CH₂F₂. The complexes are typically less strongly bound than their hydrogen-bonded counterparts, and the interaction energies vary in a consistent way with the periodic trend of the electron donors. The intermolecular separations and changes of the relevant intramolecular bond lengths, the wavenumber shifts of the critical vibrational modes and the extents of charge transfer correlate, by and large, with the strengths of interaction.

Benzene, an Unexpected Binding Unit in Anion–π Recognition: The Critical Role of CH/π Interactions

We report high-level ab initio calculations (CCSD(T)(full)/CBS//SCS-RI-MP2(full)/aug-cc-pwCVTZ) that demonstrate the importance of cooperativity effects when Anion–π and CH/π interactions are simultaneously established with benzene as the π-system. In fact, most of the complexes exhibit high cooperativity energies that range from 17% to 25.3% of the total interaction energy, which is indicative of the strong influence of the CH/π on the Anion–π interaction and vice versa. Moreover, the symmetry-adapted perturbation theory (SAPT) partition scheme was used to study the different energy contributions to the interaction energies and to investigate the physical nature of the interplay between both interactions. Furthermore, the Atoms in Molecules (AIM) theory and the Non-Covalent Interaction (NCI) approach were used to analyze the two interactions further. Finally, a few examples from the Protein Data Bank (PDB) are shown. All results stress that the concurrent formation of both interactions may play an important role in biological systems due to the ubiquity of CH bonds, phenyl rings, and anions in biomolecules.

Triel bonds within anionanion complexes

The ability of two anions to interact with one another is tested in the context of pairs of TrX₄^- homodimers, where Tr represents any of the triel atoms B, Al, Ga, In, or Tl, and X refers to a halogen substituent F, Cl, or Br. None of these pairs engage in a stable complex in the gas phase, but the situation reverses in water where the two monomers are held together by Tr⋯X triel bonds, complemented by stabilizing interactions between X atoms. Some of these bonds are quite strong, notably those involving TrF₄^-, with interaction energies surpassing 30 kcal mol^-1. Others are very much weaker, with scarcely exothermic binding energies. The highly repulsive electrostatic interactions are counteracted by large polarization energies.

Anatomy of π-hole bonds: Linear systems

The list of σ-hole bonds is long and growing, encompassing both H-bonds and its closely related halogen, chalcogen, etc., sisters. These bonds rely on the asymmetric distribution of electron density, whose depletion along the extension of a covalent bond leaves a positive region of electrostatic potential from which these bonds derive their name. However, the density distributions of other molecules contain analogous positive regions that lie out of the molecular plane known as π-holes, which are likewise capable of engaging in noncovalent bonds. Quantum calculations are applied to study such π-hole bonds that involve linear molecules, whose positive region is a circular belt surrounding the molecule, rather than the more restricted area of a σ-hole. These bonds are examined in terms of their most fundamental elements arising from the spatial dispositions of their relevant molecular orbitals and the π-holes in both the total electron density and the electrostatic potential to which they lead. Systems examined comprise tetrel, chalcogen, aerogen, and triel bonds, as well as those involving group II elements, with atoms drawn from various rows of the Periodic Table. The π-hole bonds established by linear molecules tend to be weaker than those of comparable planar systems.

Yet another perspective on hole interactions

Hole interactions are known by different names depending on the key atom of the bond (halogen bond, chalcogen bond, hydrogen bond, etc.), and the geometry of the interaction (σ if in line, π if perpendicular to the Lewis acid plane). However, its origin starts with the creation of a Lewis acid by an underlying covalent bond, which forms an electrostatic depletion and a virtual antibonding orbital, which can create non-covalent interactions with Lewis bases. In this (maybe subjective) perspective, we will claim that hole interactions must be defined via the molecular orbital origin of the molecule. Under this premise we can better explore the richness of such bonding patterns. For that, we will study old, recent and new systems, trying to pinpoint some misinterpretations that are often associated with them. We will use as exemplars the triel bonds, a couple of metal complexes, a discussion on convergent σ-holes, and many cases of anti-electrostatic hole interactions.

Competition between Inter and Intramolecular Tetrel Bonds: Theoretical Studies Complemented by CSD Survey

Crystal structures document the ability of a TF₃ group (T=Si, Ge, Sn, Pb) situated on a naphthalene system to engage in an intramolecular tetrel bond (TB) with an amino group on the adjoining ring. Ab initio calculations evaluate the strength of this bond and evaluate whether it can influence the ability of the T atom to engage in a second, intermolecular TB with another nucleophile. A very strong CN^- anionic base can approach the T either along the extension of a T-C or T-F bond and form a strong TB with an interaction energy approaching 100 kcal/mol, although this bond is weakened a bit by the presence of the internal T⋅⋅⋅N bond. The much less potent NCH base engages in a correspondingly longer and weaker TB, less than 10 kcal/mol. Such an intermolecular TB is weakened by the presence of the internal TB, to the point that it only occurs for the two heavier tetrel atoms Sn and Pb.

Noncovalent Bonds through Sigma and Pi-Hole Located on the Same Molecule. Guiding Principles and Comparisons

Over the last years, scientific interest in noncovalent interactions based on the presence of electron-depleted regions called σ-holes or π-holes has markedly accelerated. Their high directionality and strength, comparable to hydrogen bonds, has been documented in many fields of modern chemistry. The current review gathers and digests recent results concerning these bonds, with a focus on those systems where both σ and π-holes are present on the same molecule. The underlying principles guiding the bonding in both sorts of interactions are discussed, and the trends that emerge from recent work offer a guide as to how one might design systems that allow multiple noncovalent bonds to occur simultaneously, or that prefer one bond type over another.

Anion-anion and anion-neutral triel bonds

The ability of a TrCl4- anion (Tr = Al, Ga, In, Tl) to engage in a triel bond with both a neutral NH3 and CN- anion is assessed by ab initio quantum calculations in both the gas phase and in aqueous medium. Despite the absence of a positive σ or π-hole on the Lewis acid, strong triel bonds can be formed with either base. The complexation involves an internal restructuring of the tetrahedral TrCl4- monomer into a trigonal bipyramid shape, where the base can occupy either an axial or equatorial position. Although this rearrangement requires a substantial investment of energy, it aids the complexation by imparting a much more positive MEP to the site that is to be occupied by the base. Complexation with the neutral base is exothermic in the gas phase and even more so in water where interaction energies can exceed 30 kcal mol-1. Despite the long-range coulombic repulsion between any pair of anions, CN- can also engage in a strong triel bond with TrCl4-. In the gas phase, complexation is endothermic, but dissociation of the metastable dimer is obstructed by an energy barrier. The situation is entirely different in solution, with large negative interaction energies of as much as -50 kcal mol-1. The complexation remains an exothermic process even after the large monomer deformation energy is factored in.

In the search for ditriel B⋯Al non-covalent bonding

The ditriel B⋯Al interaction has been characterised using SAPT, AIM and ELF.

Crystal engineering strategies towards halogen-bonded metal–organic multi-component solids: salts, cocrystals and salt cocrystals

This highlight presents an overview of the current advances in the preparation of halogen bonded metal–organic multi-component solids, including salts and cocrystals comprising neutral and ionic constituents.

Dominance of unique Pπ phosphorus bonding with π donors: evidence using matrix isolation infrared spectroscopy and computational methodology

Albeit the first account of hypervalentπ interactions has been reported with halogenπ interactions, the feasibility of their extension to other hypervalent atoms as possible Lewis acids is still open. In this work, the role of phosphorus as an acceptor from the π electron cloud (Pπ pnicogen or phosphorus bonding) in PCl3-C2H2 and PCl3-C2H4 heterodimers is explored, by combining matrix isolation infrared spectroscopy with ab initio and DFT computational methodologies. The respective potential energy surfaces of the PCl3-C2H2 and PCl3-C2H4 heterodimers reveal unique minima stabilized by a concert of reasonably strong to weak interactions, of which Pπ phosphorus bonding was energetically dominant. Heterodimers, trimers and tetramers bound primarily by this unique phosphorus bond were generated at low temperatures. The dominance of phosphorus bonding in the PCl3-C2H2 and PCl3-C2H4 heterodimers over other interactions (such as Hπ, HCl, HP, Clπ and lone pair-π interactions) was confirmed and substantiated using extended quantum theory of atoms in molecules, natural bond orbital, electrostatic potential mapping and energy decomposition analyses. The following inferences in correlation with results from non-covalent-interaction analysis offer a complete understanding of the nature of the Pπ phosphorus bonding interactions. The significance of electrostatic forces kinetically favoring the formation of phosphorus bonded heterodimers, in addition to thermodynamic stabilization, is demonstrated experimentally.

Competition between Intra and Intermolecular Triel Bonds. Complexes between Naphthalene Derivatives and Neutral or Anionic Lewis Bases

: A TrF2 group (Tr = B, Al, Ga, In, Tl) is placed on one of the α positions of naphthalene, and its ability to engage in a triel bond (TrB) with a weak (NCH) and strong (NC-) nucleophile is assessed by ab initio calculations. As a competitor, an NH2 group is placed on the neighboring Cα, from which point it forms an intramolecular TrB with the TrF2 group. The latter internal TrB reduces the intensity of the π-hole on the Tr atom, decreasing its ability to engage in a second external TrB. The intermolecular TrB is weakened by a factor of about two for the smaller Tr atoms but is less severe for the larger Tl. The external TrB can be quite strong nonetheless; it varies from a minimum of 8 kcal/mol for the weak NCH base, up to as much as 70 kcal/mol for CN-. Likewise, the appearance of an external TrB to a strong base like CN- lessens the ability of the Tr to engage in an internal TrB, to the point where such an intramolecular TrB becomes questionable.

α-methylation and α-fluorination electronic effects on the regioselectivity of carbonyl groups of uracil by H and triel bonds in the interaction of U, T and 5FU with HCl and TrH3 (Tr = B, Al)

Quantum chemical calculations at the ωB97XD/6-311++G(d,p) level of theory have been executed to investigate the effect of substituents via hydrogen-bonded and triel-bonded complexes between uracil (U), thymine (T) and 5-fluorouracil (5FU) with HCl for the former complexes, and with BH₃ and AlH₃ for the latter complexes. These calculations are supported by single-point energy calculations at MP2/6-311++G(d,p) and CCSD/6-31 + G(d,p) levels of theory, Natural Bond Orbital (NBO) and Molecular Electrostatic Potentials (MEPs) analyses, and global/local reactivity descriptors. The results reveal that triel-bonded complexes are strongly bounded than hydrogen-bonded ones, and Al-containing dimers stronger than B-containing ones. In addition, as the central triel atom grows in size, B-containing dimers (B-O triel bond) are accompanied by weak B-H⋯O unconventional H-bonds. According to local reactivity descriptors, the B-O triel bond is hard-hard interaction that indicates that the association is primarily charge controlled, while the Al-O triel bond is soft-soft interaction that is primarily orbital controlled. In both Hydrogen as well as triel-bonded complexes, the α-methylation slightly overestimates the binding strength of U, while the α-fluorination exerts the opposite role by underestimating the binding strength of U. In overall, the effect of substituents on the bond strength and thus on the regioselectivity is very small, suggesting a competition between the two carbonyl groups in terms of structures and binding energies.

Cooperative effects between π-hole triel and π-hole chalcogen bonds

MP2/aug-cc-pVTZ calculations have been performed on π-hole triel- and chalcogen-bonded complexes involving a heteroaromatic compound. These complexes are very stable with large interaction energy up to −47 kcal mol⁻¹. The sp²-hybridized nitrogen atom engages in a stronger π-hole bond than the sp-hybridized species although the former has smaller negative electrostatic potential. The sp²-hybridized oxygen atom in 1,4-benzoquinone is a weaker electron donor in the π-hole bond than the sp²-hybridized nitrogen atom. The π-hole triel bond is stronger than the π-hole chalcogen bond. A clear structural deformation is found for the triel or chalcogen donor molecule in these π-hole-bonded complexes. The triel bond exhibits partially covalent interaction, whereas the chalcogen bond exhibits covalent interaction in the SO₃ complexes of pyrazine and pyridine derivatives with a sp²-hybridized nitrogen atom. Intermolecular charge transfer (>0.2e) occurs to a considerable extent in these complexes. In ternary complexes involving an aromatic compound, wherein a triel bond and a chalcogen bond coexist, both the interactions are weakened or strengthened when the central aromatic molecule acts as a double Lewis base or plays a dual role of both a base and an acid. Both electrostatic and charge transfer effects have important contributions toward changes in the strength of both interactions.

MP2/aug-cc-pVTZ calculations have been performed on π-hole triel- and chalcogen-bonded complexes involving a heteroaromatic compound. Both interactions exhibit cooperative/diminutive effect, depending on the role of the central heteroaromatic compound.