Charge-assisted triel bonding interactions in solid state chemistry: A combined computational and crystallographic study

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.

Osme Bond: Geometric and Energetic Features in the Adducts between OsO4 and Lewis Bases

Adducts between OsO4 and Lewis bases exert a role in important oxidation processes such as epoxidation and dihydroxylation. It has been shown that the attractive interaction driving the formation of these adducts is a σ-hole bond involving the metal as the electrophilic species; the term Osme Bond (OmB) was proposed for designating it. Here some new adducts between OsO4 and various bases have been characterized through single crystal x-ray diffraction (XRD) and computational studies (density functional theory, DFT), confirming the existence of a robust correlation between σ-hole interaction energy and deformation of the tetrahedral geometry of OsO4. Also, some adducts formed by RuO4 with nucleophiles were investigated computationally.

Resonance-assisted intramolecular triel bonds

The possibility that the intramolecular Tr⋯S triel bond is strengthened by resonance is examined by quantum chemical calculations within the planar five-membered ring of TrH₂-CRCR-CRS (Tr = Al, Ga, In; R = NO₂, CH₃). This internal bond is found to be rather short (2.4-2.7 Å) with a large bond energy between 12 and 21 kcal mol^-1. The pattern of bond length alternation and atomic charges within the ring is consistent with resonance involving the conjugated double bonds. This resonance enhances the triel bond strength by some 25%. The electron-withdrawing NO₂ group weakens the bond, but it is strengthened by the electron-donating CH₃ substituent. NICS analysis suggests the presence of a certain degree of aromaticity within the ring.

Principles Guiding the Square Bonding Motif Containing a Pair of Chalcogen Bonds between Chalcogenadiazoles

The bonding motif adopted by a dimer of chalcogenadiazole molecules is characterized by a pair of equivalent Ch···N chalcogen bonds. Quantum calculations show that the interaction energy is substantial, varying between 4 kcal/mol for Ch = S and 17 kcal/mol for Te. The interaction is cooperative in that the total bond strength is greater than either chalcogen bond individually. Neither the addition of a phenyl ring nor the addition of a pair of cyano substituents to the diazole ring has much influence on this binding. Removal of one N from the diazole weakens the binding, and addition of two nitrogens has little effect. The largest perturbation arises with three N atoms in each ring, for which the binding energy increases by some 25%. The ring size plays a minor role in most cases, although a near doubling of bond strength occurs if there are two N atoms present on a four-membered ring.

Search for an exothermic halogen bond between anions

A halogen-bonded complex containing a pair of anions can be made more stable than the isolated anions if the Lewis acid is a long carbon chain, fully substituted by CN groups, with an I atom on one end and a COO− group on the other, with Cl− as base.

On the Ability of Nitrogen to Serve as an Electron Acceptor in a Pnicogen Bond

Whereas pnicogen atoms like P and As have been shown repeatedly to act as electron acceptors in pnicogen bonds, the same is not true of the more electronegative first-row N atom. Quantum calculations assess whether N can serve in this capacity in such bonds and under what conditions. There is a positive π-hole belt that surrounds the central N atom in the linear arrangement of NNNF, NNN-CN, and NNO, which can engage a NH₃ base to form a pnicogen bond with binding energy between 3 and 5 kcal/mol. Within the context of a planar arrangement, the π-hole above the N in NO₂OF, N(CN)₃, and CF₃NO₂ is also capable of forming a pnicogen bond, the strongest of which amounts to 11 kcal/mol with NMe₃ as base. In their pyramidal geometry, NF₃ and N(NO₂)₃ engage with a base through the σ-hole on the central N, with variable binding energies between 2 and 9 kcal/mol. AIM and NBO provide somewhat different interpretations of the secondary interactions that occur in some of these complexes.

Maximal occupation by bases of π-hole bands surrounding linear molecules

Linear molecules such as CO₂ contain a positive π-hole ring that surrounds C on the molecule's equator. Quantum calculations examine the question as to how many bases can simultaneously bind to this ring. Linear molecules examined are TO₂ , where T = C, Si, Ge, Sn; bases are NCH and NH₃ . CO₂ engages in the weakest of the tetrel bonds, and can bind up to three NCH and two NH₃ . Unlike σ-hole tetrel bonds, Si forms the strongest tetrel bonds, with interaction energies as high as 43 kcal/mol with NH₃ . But like GeO₂ , SiO₂ can sustain only two bases in its equatorial ring. The π-hole ring of SnO₂ can engage in up to four tetrel bonds with either NCH or NH₃ , even though these bonds are weaker than those with GeO₂ or SiO₂ . As all of these complexes cast TO₂ in the role of multiple electron acceptor, the resulting negative cooperativity makes each successive bond weaker than its predecessor as bases are added, as well as reducing the magnitude of the central molecule's π-hole.

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.

Dissection of the Origin of π-Holes and the Noncovalent Bonds in Which They Engage

Accompanying the rapidly growing list of σ-hole bonds has come the acknowledgment of parallel sorts of noncovalent bonds which owe their stability in large part to a deficiency of electron density in the area above the molecular plane, known as a π-hole. The origins of these π-holes are probed for a wide series of molecules, comprising halogen, chalcogen, pnicogen, tetrel, aerogen, and spodium bonds. Much like in the case of their σ-hole counterparts, formation of the internal covalent π-bond in the Lewis acid molecule pulls density toward the bond midpoint and away from its extremities. This depletion of density above the central atom is amplified by an electron-withdrawing substituent. At the same time, the amplitude of the π*-orbital is enhanced in the region of the density-depleted π-hole, facilitating a better overlap with the nucleophile's lone pair orbital and a stabilizing n → π* charge transfer. The presence of lone pairs on the central atom acts to attenuate the π-hole and shift its position somewhat, resulting in an overall weakening of the π-hole bond. There is a tendency for π-hole bonds to include a higher fraction of induction energy than σ-bonds with proportionately smaller electrostatic and dispersion components, but this distinction is less a product of the σ- or π-character and more a function of the overall bond strength.

Proximity Effects of Substituents on Halogen Bond Strength

It is well known that the presence of an electron-withdrawing substituent (EWS) placed near the halogen (X) atom on a Lewis acid molecule amplifies the ability of this unit to engage in a halogen bond with a base. Quantum calculations are applied to examine how quickly these effects fade as the EWS is moved further and further from the X atom. Conjugated alkene and alkyne chains of varying lengths with a terminal C-I first facilitate analysis as to how the number of these multiple bonds affects the strength of CI··N XB to NH₃. Then, electron-withdrawing F and C≡N substituents are placed on the opposite end of the chain, and their effects on the XB properties are monitored as a function of their distance from I. These same EWSs are added to the ortho, meta, and para positions of aromatic iodobenzene. It is found that the XB grows in strength as more triple bonds are added to the alkyne, but there is little change caused by elongating an alkene. The cyano group has a much stronger effect than does F. While F strengthens the XB, its effects are quickly attenuated as it is moved further from I. The consequences of C≡N substitution are stronger and extend over a longer distance. Placement of an EWS on the phenyl ring diminishes with distance: o > m > p, and the effects of disubstitution are nearly additive. These trends apply not only to energetics but also to geometries, properties of the wave function, σ-hole depth, and NMR shielding.

Relative Strengths of a Pnicogen and a Tetrel Bond and Their Mutual Effects upon One Another

The ability of the T and Z atoms of TR₃ZR₂ to engage in a noncovalent interaction with NH₃ is assessed by DFT calculations, where the T atom refers to C, Si, and Ge; Z = As, Sb, and P; and substituents R = H and F. In most instances, the tetrel bond (TB) is both stronger and shorter than the pnicogen bond (ZB). These two bond strengths can be equalized, or preference shifted to the ZB, if F substituents are placed on the Z and H on the T atoms. Employing C as the T atom results in a very weak TB, with the ZB clearly favored energetically. The simultaneous formation of both TB and ZB weakens both, particularly the latter, but both bonds survive intact. Geometric and spectroscopic perturbations of the subunits reflect the two types of noncovalent bonds.

Understanding noncovalent bonds and their controlling forces

The fundamental underpinnings of noncovalent bonds are presented, focusing on the σ-hole interactions that are closely related to the H-bond. Different means of assessing their strength and the factors that control it are discussed. The establishment of a noncovalent bond is monitored as the two subunits are brought together, allowing the electrostatic, charge redistribution, and other effects to slowly take hold. Methods are discussed that permit prediction as to which site an approaching nucleophile will be drawn, and the maximum number of bonds around a central atom in its normal or hypervalent states is assessed. The manner in which a pair of anions can be held together despite an overall Coulombic repulsion is explained. The possibility that first-row atoms can participate in such bonds is discussed, along with the introduction of a tetrel analog of the dihydrogen bond.

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.

F-Halogen Bond: Conditions for Its Existence

The question as to whether the F atom can engage in a halogen bond (XB) remains unsettled. This issue is addressed via density functional theory calculations which pair a wide range of organic and inorganic F-acids with various sorts of Lewis bases. From an energetic perspective, perfluorinated hydrocarbons with sp, sp², or sp³ C-hybridization are unable to form an XB with an N-base, but a very weak bond can be formed if electron-withdrawing C≡N substituents are added to the acid. There is little improvement for inorganic acids O₂NF, FOF, ClF, BrF, SiF₄, or GeF₄, but F₂ is capable of a stronger XB of up to 5 kcal/mol. These results are consistent with a geometric criterion, which compares the intermolecular equilibrium distance with the sum of atomic van der Waals radii. The intensity of the σ-hole on the F atom has predictive value in that a V_s,max of at least 10-15 kcal/mol is required for XB formation. Adding a positive charge to the Lewis acid enhances the strength of any XB and even more so if the base is anionic. The acid-base interaction induces a contraction of the r(AF) covalent bond in the acid in most cases and a deshielding of the NMR signal of the F nucleus.

Coordination of a Central Atom by Multiple Intramolecular Pnicogen Bonds

A central pnicogen Z atom (Z = Sb, As) is covalently attached to the O atom of three -O(CH₂)_nX chains where X represents either an aldehyde or amine group. The chain can fold around so that the basic X group can engage in a noncovalent pnicogen bond with the central Z. The formation of up to three pnicogen bonds is energetically favored. The amine appears to engage in stronger pnicogen bonds than does the aldehyde, and bonds to Sb are favored over As, but there is little dependence on the length of the chain. The formation of each successive pnicogen bond reduces the magnitude of the σ-holes surrounding the Z atom, which tends to weaken the attraction for the basic end of the chain.

Not Only Hydrogen Bonds: Other Noncovalent Interactions

In this review, we provide a consistent description of noncovalent interactions, covering most groups of the Periodic Table. Different types of bonds are discussed using their trivial names. Moreover, the new name “Spodium bonds” is proposed for group 12 since noncovalent interactions involving this group of elements as electron acceptors have not yet been named. Excluding hydrogen bonds, the following noncovalent interactions will be discussed: alkali, alkaline earth, regium, spodium, triel, tetrel, pnictogen, chalcogen, halogen, and aerogen, which almost covers the Periodic Table entirely. Other interactions, such as orthogonal interactions and π-π stacking, will also be considered. Research and applications of σ-hole and π-hole interactions involving the p-block element is growing exponentially. The important applications include supramolecular chemistry, crystal engineering, catalysis, enzymatic chemistry molecular machines, membrane ion transport, etc. Despite the fact that this review is not intended to be comprehensive, a number of representative works for each type of interaction is provided. The possibility of modeling the dissociation energies of the complexes using different models (HSAB, ECW, Alkorta-Legon) was analyzed. Finally, the extension of Cahn-Ingold-Prelog priority rules to noncovalent is proposed.

The ditetrel bond: noncovalent bond between neutral tetrel atoms

The ability of a tetrel atom to serve in the capacity of electron donor in a σ-hole noncovalent bond is tested by quantum calculations.

Effects of Halogen, Chalcogen, Pnicogen, and Tetrel Bonds on IR and NMR Spectra

Complexes were formed pairing FX, FHY, FH2Z, and FH3T (X = Cl, Br, I; Y = S, Se, Te; Z = P, As, Sb; T = Si, Ge, Sn) with NH3 in order to form an A⋯N noncovalent bond, where A refers to the central atom. Geometries, energetics, atomic charges, and spectroscopic characteristics of these complexes were evaluated via DFT calculations. In all cases, the A-F bond, which is located opposite the base and is responsible for the σ-hole on the A atom, elongates and its stretching frequency undergoes a shift to the red. This shift varies from 42 to 175 cm-1 and is largest for the halogen bonds, followed by chalcogen, tetrel, and then pnicogen. The shift also decreases as the central A atom is enlarged. The NMR chemical shielding of the A atom is increased while that of the F and electron donor N atom are lowered. Unlike the IR frequency shifts, it is the third-row A atoms that undergo the largest change in NMR shielding. The change in shielding of A is highly variable, ranging from negligible for FSnH3 all the way up to 1675 ppm for FBr, while those of the F atom lie in the 55-422 ppm range. Although smaller in magnitude, the changes in the N shielding are still easily detectable, between 7 and 27 ppm.

A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts

The founding in 1965 of what is now called the Cambridge Structural Database (CSD) has reaped dividends in numerous and diverse areas of chemical research. Each of the million or so crystal structures in the database was solved for its own particular reason, but collected together, the structures can be reused to address a multitude of new problems. In this Review, which is focused mainly on the last 10 years, we chronicle the contribution of the CSD to research into molecular geometries, molecular interactions, and molecular assemblies and demonstrate its value in the design of biologically active molecules and the solid forms in which they are delivered. Its potential in other commercially relevant areas is described, including gas storage and delivery, thin films, and (opto)electronics. The CSD also aids the solution of new crystal structures. Because no scientific instrument is without shortcomings, the limitations of CSD research are assessed. We emphasize the importance of maintaining database quality: notwithstanding the arrival of big data and machine learning, it remains perilous to ignore the principle of garbage in, garbage out. Finally, we explain why the CSD must evolve with the world around it to ensure it remains fit for purpose in the years ahead.

Differential Binding of Tetrel-Bonding Bipodal Receptors to Monatomic and Polyatomic Anions

Previous work has demonstrated that a bidentate receptor containing a pair of Sn atoms can engage in very strong interactions with halide ions via tetrel bonds. The question that is addressed here concerns the possibility that a receptor of this type might be designed that would preferentially bind a polyatomic over a monatomic anion since the former might better span the distance between the two Sn atoms. The binding of Cl⁻ was thus compared to that of HCOO⁻, HSO₄⁻, and H₂PO₄⁻ with a wide variety of bidentate receptors. A pair of SnFH₂ groups, as strong tetrel-binding agents, were first added to a phenyl ring in ortho, meta, and para arrangements. These same groups were also added in 1,3 and 1,4 positions of an aliphatic cyclohexyl ring. The tetrel-bonding groups were placed at the termini of (-C≡C-)_n (n = 1,2) extending arms so as to further separate the two Sn atoms. Finally, the Sn atoms were incorporated directly into an eight-membered ring, rather than as appendages. The ordering of the binding energetics follows the HCO₂⁻ > Cl⁻ > H₂PO₄⁻ > HSO₄⁻ general pattern, with some variations in selected systems. The tetrel bonding is strong enough that in most cases, it engenders internal deformations within the receptors that allow them to engage in bidentate bonding, even for the monatomic chloride, which mutes any effects of a long Sn···Sn distance within the receptor.

Dual function of the boron center of BH(CO)2/BH(N2)2 in halogen- and triel-bonded complexes with hypervalent halogens

The complexes between BH(CO)₂/BH(N₂)₂ and XF₃/XF₅ are stabilized by a halogen bond and a triel bond. The MEP analyses of BH(CO)₂/BH(N₂)₂ indicate that there are both a region with negative MEPs on the B atom in the vertical direction of the molecular plane and a σ-hole at the B-H bond end. Therefore, the boron atom in BH(CO)₂/BH(N₂)₂ plays a dual role of a Lewis base and an acid in the halogen bond and triel bond, respectively. The halogen and triel bonds are stronger in order of IF₃< BrF₃< ClF₃, IF₅< BrF₅< ClF₅, and BH(CO)₂< BH(N₂)₂ in most complexes. These complexes have large stability since the interaction energy varies from -5 to -115 kcal/mol. The halogen bond belongs to a covalent interaction or a partially covalent interaction in most complexes. The subsystems in these complexes have prominent deformation, accompanied with big charge transfer and large polarization.

Concurrent aerogen bonding and lone pair/anion-π interactions in the stability of organoxenon derivatives: a combined CSD and ab initio study

In this manuscript the ability of organoxenon fluorides to establish concurrent π-hole aerogen bonding and lone pair/anion-π interactions has been studied at the RI-MP2/def2-TZVP level of theory. The presence of both an aromatic system (benzene, trifluorobenzene and pentafluorobenzene) and a xenon atom makes these molecules suitable for simultaneously establishing both interactions. In this regard, we have used CH₃CN, NH₃, O(CH₃)₂, Cl^-, CN^- and BF₄^- as neutral and charged electron donors, respectively. Moreover, the NBO analysis showed that orbital effects contribute to the global stabilization of the complexes studied. Furthermore, we have used Bader's theory of "atoms in molecules" to analyse and characterize the noncovalent interactions described herein from a charge-density perspective. Finally, several examples retrieved from the CSD are also included, highlighting the impact of these interactions in the solid state chemistry of Xe.

Alkyl groups as electron density donors in π-hole bonding

A combined structural and computational analysis has demonstrated that alkyl groups can act as Lewis bases in π-hole bonding.