Assembly of Effective Halide Receptors from Components. Comparing Hydrogen, Halogen, and Tetrel Bonds

Cooperativity between Intermolecular Hydrogen and Carbon Bonds in ZY···CH3CN/CH3NC···HX Trimers (ZY = H2O, H2S, HF, HCl, HBr, NH3, and H2CO; HX = HF, HCl, and HBr)

Hydrogen-bonding and carbon-bonding interactions are widespread in nature. We studied the cooperativity between these interactions in 42 trimeric complexes ZY···CH3CN/CH3NC···HX, where ZY molecules are H2O, H2S, HF, HCl, HBr, NH3, and H2CO, and HX molecules are HF, HCl, and HBr. Acetonitrile (CH3CN) and isoacetonitrile (CH3NC) act as hydrogen bond acceptors as well as carbon bond donors in these trimers. Various theoretical methods, such as electronic structure calculations, quantum theory of atoms in molecule (QTAIM), natural bond orbital (NBO), and reduced density gradient analysis, are employed to study these trimers, and the results are compared with the corresponding ZY···CH3CN/CH3NC and CH3CN/CH3NC···HX dimers. Electronic structure calculations are performed at the second-order Mo̷ller-Plesset perturbation theory using the 6-311++G(2d,2p) basis set. We show that both the interactions act synergistically in these trimers leading to an increase in their bond strength as compared to the strength in the individual dimers. The cooperative energies for these trimers are in the range of 0.69 to 3.22 kJ/mol. It is seen that the carbon bonds benefit more from the cooperativity than the hydrogen bonds. The trends of cooperativity and correlations of interaction energies and cooperative energies with relevant QTAIM and NBO parameters are reported.

Competition between Binding to Various Sites of Substituted Imidazoliums

The imidazolium cation has a number of different sites that can interact with a nucleophile. Adding a halogen atom (X) or a chalcogen (YH) group introduces the possibility of an NX···nuc halogen or NY···nuc chalcogen bond, which competes against the various H-bonds (NH and CH donors) as well as the lone pair···π interaction wherein the nucleophile lies above the plane of the cation. Substituted imidazoliums are paired with the NH₃ base, and the various different complexes are evaluated by density functional theory (DFT) calculations. The strength of XB and YB increases quickly along with the size and polarizability of the X/Y atom, and this sort of bond is the strongest for the heavier Br, I, Se, and Te atoms, followed by the NH···N H-bond, but this order reverses for Cl and S. The various CH···N H-bonds are comparable to one another and to the lone pair···π bond, all with interaction energies of 10-13 kcal/mol, values which show very little dependence upon the substituent placed on the imidazolium.

Selectivity Rule of Cryptands for Anions: Molecular Rigidity and Bonding Site

Cryptands utilize inside CH or NH groups as hydrogen bond (H-bond) donors to capture anions such as halides. In this work, the nature and selectivity of confined hydrogen bonds inside cryptands were computationally analyzed with the energy decomposition scheme based on the block-localized wavefunction method (BLW-ED), aiming at an elucidation of governing factors in the binding between cryptands and anions. It was revealed that the intrinsic strengths of inward hydrogen bonds are dominated by the electrostatic attraction, while the anion preferences (selectivity) of inner CH and NH hydrogen bonds are governed by the Pauli exchange repulsion and electrostatic interaction, respectively. Typical conformers of cages are classified into two groups, including the C_3(h) -symmetrical conformers, in which all halide anions are located near the centroids of cages, and the "semi-open" conformers, which exhibit shifted bonding sites for different halide anions. Accordingly, the difference in governing factors of selectivity is attributed to either the rigidity of cages or the binding site of anions for these two groups. In details, the C₃ conformers of NH cryptands can be enlarged more remarkably than the C_3(h) -symmetrical conformers of CH cryptands as the size of anion (ionic radius) increases, resulting in the relaxation of the Pauli repulsion and a dramatic reduction in electrostatic attraction, which eventually rules the selectivity of NH cryptands for halide anions. By contrary, the CH cryptands are more rigid and cannot effectively reduce the Pauli repulsion, which subsequently governs the anion preference. Unlike C₃ conformers whose rigidity determines the selectivity, semi-open conformers exhibit different binding sites for different anions. From F^- to I^- , the bonding site shifts toward the outside end of the pocket inside the semi-open NH cryptand, leading to the significant reduction of the electrostatic interaction that dominates the anion preference. Differently, binding sites are much less affected by the size of anion inside the semi-open CH cryptand, in which the Pauli exchange repulsion remains the key factor for the selectivity of inner hydrogen bonds.

Promotion of TH3 (T = Si and Ge) group transfer within a tetrel bond by a cation-π interaction

The possibility of the transfer of the TH₃ group across a tetrel bond is considered by ab initio calculations. The TB is constructed by pairing PhTH₃ (Ph = phenyl; T = Si and Ge) with bases NH₃, NHCH₂, and the C₃N₂H₄ carbene. The TH₃ moves toward the base but only by a small amount in these dimers. However, when a Be²⁺ or Mg²⁺ dication is placed above the phenyl ring, the tetrel bond strength is greatly magnified reaching up to nearly 100 kcal mol^-1. This dication also induces a much higher degree of transfer which can be best categorized as half-transfer for the two N-bases and a near complete transfer for the carbene.

Enhancement of the Tetrel Bond by the Effects of Substituents, Cooperativity, and Electric Field: Transition from Noncovalent to Covalent Bond

The T⋅⋅⋅N tetrel bond (TB) formed between TX₃ OH (T=C, Si, Ge; X=H, F) and the Lewis base N≡CM (M=H, Li, Na) is studied by ab initio calculations at the MP2/aug-cc-pVTZ level. Complexes involving TH₃ OH contain a conventional TB with interaction energy less than 10 kcal/mol. This bond is substantially strengthened, approaching 35 kcal/mol and covalent character, when fluorosubstituted TF₃ OH is combined with NCLi or NCNa. Along with this enhanced binding comes a near equalization of the TB T⋅⋅⋅N and the internal T-O bond lengths, and the associated structure acquires a trigonal bipyramidal shape, despite a high internal deformation energy. This structural transformation becomes more complete, and the TB is further strengthened upon adding an electron acceptor BeCl₂ to the Lewis acid and a base to the NCM unit. This same TB strengthening can be accomplished also by imposition of an external electric field.

Comparison for Electron Donor Capability of Carbon-Bound Halogens in Tetrel Bonds

The tetrel bond formed by HC≡CX, H2C=CHX, and H3CCH2X (X=F, Cl, Br, I) as an electron donor and TH3F (T=C, Si, Ge) was explored by ab initio calculations. The tetrel bond formed by H3CCH2X is the strongest, as high as -3.45 kcal/mol for the H3CCH2F···GeH3F dimer, followed by H2C=CHX, and the weakest bond is from HC≡CX, where the tetrel bond can be as small as -0.8 kcal/mol. The strength of the tetrel bond increases in the order of C < Si < Ge. For the H3CCH2X and HC≡CX complexes, the tetrel bond strength shows a similar increasing tendency with the decrease of the electronegativity of the halogen atom. Electrostatic interaction plays the largest role in the stronger tetrel bonds, while dispersion interaction makes an important contribution to the H2C=CHX complexes.

Evaluation of Electron Density Shifts in Noncovalent Interactions

In the present paper, we report the quantitative evaluation of the electron density shift (EDS) maps within different complexes. Values associated with the total EDS maps exhibited good correlation with different quantities such as interaction energies, Eint, intermolecular distances, bond critical points, and LMOEDA energy decomposition terms. Besides, EDS maps at different cutoffs were also evaluated and related with the interaction energies values. Finally, EDS maps and their corresponding values are found to correlate with Eint within systems with cooperative effects. To our knowledge, this is the first time that the EDS has been quanitatively evaluated.

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 bond between tetrel π-hole and hydride

The π-hole above the plane of the X2T'Y molecule (T' = Si, Ge, Sn; X = F, Cl, H; Y = O, S) was allowed to interact with the TH hydride of TH(CH3)3 (T = Si, Ge, Sn). The resulting THT' tetrel bond is quite strong, with interaction energies exceeding 30 kcal mol-1. F2T'O engages in the strongest such bonds, as compared to F2T'S, Cl2T'O, or Cl2T'S. The bond weakens as T' grows larger as in Si > Ge > Sn, despite the opposite trend in the depth of the π-hole. The reverse pattern of stronger tetrel bond with larger T is observed for the Lewis base TH(CH3)3, even though the minimum in the electrostatic potential around the H is nearly independent of T. The THT' arrangement is nonlinear which can be understood on the basis of the positions of the extrema in the molecular electrostatic potentials of the monomers. The tetrel bond is weakened when H2O forms an OT' tetrel bond with the second π-hole of F2T'O, and strengthened if H2O participates in an OHO H-bond.

Comparison of Bifurcated Halogen with Hydrogen Bonds

Bifurcated halogen bonds are constructed with FBr and FI as Lewis acids, paired with NH₃ and NCH bases. The first type considered places two bases together with a single acid, while the reverse case of two acids sharing a single base constitutes the second type. These bifurcated systems are compared with the analogous H-bonds wherein FH serves as the acid. In most cases, a bifurcated system is energetically inferior to a single linear bond. There is a larger energetic cost to forcing the single σ-hole of an acid to interact with a pair of bases, than the other way around where two acids engage with the lone pair of a single base. In comparison to FBr and FI, the H-bonding FH acid is better able to participate in a bifurcated sharing with two bases. This behavior is traced to the properties of the monomers, in particular the specific shape of the molecular electrostatic potential, the anisotropy of the orbitals of the acid and base that interact directly with one another, and the angular extent of the total electron density of the two molecules.

Origins and properties of the tetrel bond

The tetrel bond (TB) recruits an element drawn from the C, Si, Ge, Sn, Pb family as electron acceptor in an interaction with a partner Lewis base. The underlying principles that explain this attractive interaction are described in terms of occupied and vacant orbitals, total electron density, and electrostatic potential. These principles facilitate a delineation of the factors that feed into a strong TB. The geometric deformation that occurs within the tetrel-bearing Lewis acid monomer is a particularly important issue, with both primary and secondary effects. As a first-row atom of low polarizability, C is a reluctant participant in TBs, but its preponderance in organic and biochemistry make it extremely important that its potential in this regard be thoroughly understood. The IR and NMR manifestations of tetrel bonding are explored as spectroscopy offers a bridge to experimental examination of this phenomenon. In addition to the most common σ-hole type TBs, discussion is provided of π-hole interactions which are a result of a common alternate covalent bonding pattern of tetrel atoms.

Reliable Comparison of Pnicogen, Chalcogen, and Halogen Bonds in Complexes of 6-OXF2-Fulvene (X = As, Sb, Se, Te, Be, I) With Three Electron Donors

The pnicogen, chalcogen, and halogen bonds between 6-OXF₂-fulvene (X = As, Sb, Se, Te, Br, and I) and three nitrogen-containing bases (FCN, HCN, and NH₃) are compared. For each nitrogen base, the halogen bond is strongest, followed by the pnicogen bond, and the chalcogen bond is weakest. For each type of bond, the binding increases in the FCN < HCN < NH₃ pattern. Both FCN and HCN engage in a bond with comparable strengths and the interaction energies of most bonds are < -6 kcal/mol. However, the strongest base NH₃ forms a much more stable complex, particularly for the halogen bond with the interaction energy going up to -18 kcal/mol. For the same type of interaction, its strength increases as the mass of the central X atom increases. These bonds are different in strength, but all of them are dominated by the electrostatic interaction, with the polarization contribution important for the stronger interaction. The presence of these bonds changes the geometries of 6-OXF₂-fulvene, particularly for the halogen bond formed by NH₃, where the F-X-F arrangement is almost vertical to the fulvene ring.

Interplay between σ Holes, Anion···H-C, and Cation-π Interactions in Dibromo[2,2]paracyclophane Complexes

Theoretical calculations were performed to investigate the interplay between σ-hole, anion-HC and cation-π interactions in the complexes of dibromo[2,2]paracyclophane (DBr[2,2]PCP) with alkali (Li⁺, Na⁺, K⁺), alkaline earth metal cations (Be²⁺, Mg²⁺, and Ca²⁺), and halogen anions (F^-, Cl^-, and Br^-) using the wave function (MP2) and density functional theory (M06-2X and B3LYP) methods with the 6-311++G(d,p) basis set. The study reveals that DBr[2,2]PCP behaves as amphoteric molecule with a predominance of basic character. It prefers to interact with hard cations and hard anions such as Be²⁺ and F^- through cation-π and anion···HC interactions, respectively. Substitution of Br by F and Cl atoms in DBr[2,2]PCP decreases slightly the interaction energies of DX[2,2]PCP-halogen complexes (X = F, Cl, and Br) by 2.0 and 0.3 kcal/mol (M06-2X), respectively. The anion-HC interactions in DBr[2,2]PCP complexes are ∼10 kcal/mol stronger (B3LYP; ∼15 kcal/mol at M06-2X and 7 kcal/mol at MP2) than the σ-hole interactions.

How Many Pnicogen Bonds can be Formed to a Central Atom Simultaneously?

A central ZF₃ molecule (Z = P, As, Sb, Bi) is allowed to interact with a number of nucleophiles exemplified by NCH, NH₃, and NC^– anion. The Z···N pnicogen bond (ZB) to a single base grows stronger for heavier Z atom: P < A < Sb < Bi and follows the NCH < NH₃ < NC^– order for the three bases. The maximum number of ZBs depends on both the nature of the base and pnicogen atom. PF₃ and AsF₃ can pnicogen bond with only a single CN^–; SbF₃ and BiF₃ can interact with two anions but only weakly. The weak NCH nucleophile can engage in a maximum of two ZBs, while three ZBs occur for NH₃. The latter NH₃ maximum can be extended to four ZBs but only for BiF₃. The fourth ZB is somewhat longer and weaker than the others, and the entire (H₃N)₄···BiF₃ complex relies partially on secondary interactions for its stability.

Anion Recognition by Neutral and Cationic Iodotriazole Halogen Bonding Scaffolds

A computational study of the iodide discrimination by different neutral and cationic iodotriazole halogen bonding hosts was carried out by means of Density Functional Theory. The importance of the size of the scaffold was highlighted and its impact observed in the binding energies and intermolecular X···I distances. Larger scaffolds were found to reduce the electronic repulsion and increase the overlap between the halide electron lone pair and the corresponding I-C antibonding orbital, increasing the halogen bonding interactions. Additionally, the planarity plays an important role within the interaction, and can be tuned using hydroxyl to perform intramolecular hydrogen bonds (IMHB) between the scaffold and the halogen atoms. Structures with IMHB exhibit stronger halogen bond interactions, as evidenced by the shorter intramolecular distances, larger electron density values at the bond critical point and more negative binding energies.

Definition of the chalcogen bond (IUPAC Recommendations 2019)

This recommendation proposes a definition for the term “chalcogen bond”; it is recommended the term is used to designate the specific subset of inter- and intramolecular interactions formed by chalcogen atoms wherein the Group 16 element is the electrophilic site.

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.

Dual Geometry Schemes in Tetrel Bonds: Complexes between TF₄ (T = Si, Ge, Sn) and Pyridine Derivatives

When an N-base approaches the tetrel atom of TF₄ (T = Si, Ge, Sn) the latter molecule deforms from a tetrahedral structure in the monomer to a trigonal bipyramid. The base can situate itself at either an axial or equatorial position, leading to two different equilibrium geometries. The interaction energies are considerably larger for the equatorial structures, up around 50 kcal/mol, which also have a shorter R(T··N) separation. On the other hand, the energy needed to deform the tetrahedral monomer into the equatorial structure is much higher than the equivalent deformation energy in the axial dimer. When these two opposite trends are combined, it is the axial geometry which is somewhat more stable than the equatorial, yielding binding energies in the 8–34 kcal/mol range. There is a clear trend of increasing interaction energy as the tetrel atom grows larger: Si < Ge < Sn, a pattern which is accentuated for the binding energies.

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.

Influence of monomer deformation on the competition between two types of σ-holes in tetrel bonds

Competition between two competing sites on a tetrel atom is explained by balance between structural deformation and σ-hole intensity.

Calculation of VS,max and Its Use as a Descriptor for the Theoretical Calculation of pKa Values for Carboxylic Acids

The theoretical calculation of pKa values for Brønsted acids is a challenging task that involves sophisticated and time-consuming methods. Therefore, heuristic approaches are efficient and appealing methodologies to approximate these values. Herein, we used the maximum surface electrostatic potential (V_S,max) on the acidic hydrogen atoms of carboxylic acids to describe the H-bond interaction with water (the same descriptor that is used to characterize σ-bonded complexes) and correlate the results with experimental pKa values to obtain a predictive model for other carboxylic acids. We benchmarked six different methods, all including an implicit solvation model (water): Five density functionals and the Møller–Plesset second order perturbation theory in combination with six different basis sets for a total of thirty-six levels of theory. The ωB97X-D/cc-pVDZ level of theory stood out as the best one for consistently reproducing the reported pKa values, with a predictive power of 98% correlation in a test set of ten other carboxylic acids.

Tetrel Bond between 6-OTX₃-Fulvene and NH₃: Substituents and Aromaticity

Carbon bonding is a weak interaction, particularly when a neutral molecule acts as an electron donor. Thus, there is an interesting question of how to enhance carbon bonding. In this paper, we found that the –OCH₃ group at the exocyclic carbon of fulvene can form a moderate carbon bond with NH₃ with an interaction energy of about −10 kJ/mol. The –OSiH₃ group engages in a stronger tetrel bond than does the –OGeH₃ group, while a reverse result is found for both –OSiF₃ and –OGeF₃ groups. The abnormal order in the former is mainly due to the stronger orbital interaction in the –OSiH₃ complex, which has a larger deformation energy. The cyano groups adjoined to the fulvene ring not only cause a change in the interaction type, from vdW interactions in the unsubstituted system of –OCF₃ to carbon bonding, but also greatly strengthen tetrel bonding. The formation of tetrel bonding has an enhancing effect on the aromaticity of the fulvene ring.

Crystallographic and Computational Characterization of Methyl Tetrel Bonding in S-Adenosylmethionine-Dependent Methyltransferases

Tetrel bonds represent a category of non-bonding interaction wherein an electronegative atom donates a lone pair of electrons into the sigma antibonding orbital of an atom in the carbon group of the periodic table. Prior computational studies have implicated tetrel bonding in the stabilization of a preliminary state that precedes the transition state in S_N2 reactions, including methyl transfer. Notably, the angles between the tetrel bond donor and acceptor atoms coincide with the prerequisite geometry for the S_N2 reaction. Prompted by these findings, we surveyed crystal structures of methyltransferases in the Protein Data Bank and discovered multiple instances of carbon tetrel bonding between the methyl group of the substrate S-adenosylmethionine (AdoMet) and electronegative atoms of small molecule inhibitors, ions, and solvent molecules. The majority of these interactions involve oxygen atoms as the Lewis base, with the exception of one structure in which a chlorine atom of an inhibitor functions as the electron donor. Quantum mechanical analyses of a representative subset of the methyltransferase structures from the survey revealed that the calculated interaction energies and spectral properties are consistent with the values for bona fide carbon tetrel bonds. The discovery of methyl tetrel bonding offers new insights into the mechanism underlying the S_N2 reaction catalyzed by AdoMet-dependent methyltransferases. These findings highlight the potential of exploiting these interactions in developing new methyltransferase inhibitors.

Strong Tetrel Bonds: Theoretical Aspects and Experimental Evidence

In recent years, noncovalent interactions involving group-14 elements of the periodic table acting as a Lewis acid center (or tetrel-bonding interactions) have attracted considerable attention due to their potential applications in supramolecular chemistry, material science and so on. The aim of the present study is to characterize the geometry, strength and bonding properties of strong tetrel-bond interactions in some charge-assisted tetrel-bonded complexes. Ab initio calculations are performed, and the results are supported by the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) approaches. The interaction energies of the anionic tetrel-bonded complexes formed between XF₃M molecule (X=F, CN; M=Si, Ge and Sn) and A- anions (A-=F-, Cl-, Br-, CN-, NC- and N₃-) vary between -16.35 and -96.30 kcal/mol. The M atom in these complexes is generally characterized by pentavalency, i.e., is hypervalent. Moreover, the QTAIM analysis confirms that the anionic tetrel-bonding interaction in these systems could be classified as a strong interaction with some covalent character. On the other hand, it is found that the tetrel-bond interactions in cationic tetrel-bonded [p-NH₃(C₆H₄)MH₃]⁺···Z and [p-NH₃(C₆F₄)MH₃]⁺···Z complexes (M=Si, Ge, Sn and Z=NH₃, NH₂CH₃, NH₂OH and NH₂NH₂) are characterized by a strong orbital interaction between the filled lone-pair orbital of the Lewis base and empty BD*M-C orbital of the Lewis base. The substitution of the F atoms in the benzene ring provides a strong orbital interaction, and hence improved tetrel-bond interaction. For all charge-assisted tetrel-bonded complexes, it is seen that the formation of tetrel-bond interaction is accompanied bysignificant electron density redistribution over the interacting subunits. Finally, we provide some experimental evidence for the existence of such charge-assisted tetrel-bond interactions in crystalline phase.

Comparative Strengths of Tetrel, Pnicogen, Chalcogen, and Halogen Bonds and Contributing Factors

Ab initio calculations are employed to assess the relative strengths of various noncovalent bonds. Tetrel, pnicogen, chalcogen, and halogen atoms are represented by third-row atoms Ge, As, Se, and Br, respectively. Each atom was placed in a series of molecular bonding situations, beginning with all H atoms, then progressing to methyl substitutions, and F substituents placed in various locations around the central atom. Each Lewis acid was allowed to engage in a complex with NH₃ as a common nucleophile, and the strength and other aspects of the dimer were assessed. In the context of fully hydrogenated acids, the strengths of the various bonds varied in the pattern of chalcogen > halogen > pnicogen ≈ tetrel. Methyl substitution weakened all bonds, but not in a uniform manner, resulting in a greatly weakened halogen bond. Fluorosubstitution strengthened the interactions, increasing its effect as the number of F atoms rises. The effect was strongest when the F atom lay directly opposite the base, resulting in a halogen > chalcogen > pnicogen > tetrel order of bond strength. Replacing third-row atoms by their second-row counterparts weakened the bonds, but not uniformly. Tetrel bonds were weakest for the fully hydrogenated acids and surpassed pnicogen bonds when F had been added to the acid.

Tetrel Bonding as a Vehicle for Strong and Selective Anion Binding

Tetrel atoms T (T = Si, Ge, Sn, and Pb) can engage in very strong noncovalent interactions with nucleophiles, which are commonly referred to as tetrel bonds. The ability of such bonds to bind various anions is assessed with a goal of designing an optimal receptor. The Sn atom seems to form the strongest bonds within the tetrel family. It is most effective in the context of a -SnF₃ group and a further enhancement is observed when a positive charge is placed on the receptor. Connection of the -SnF₃ group to either an imidazolium or triazolium provides a strong halide receptor, which can be improved if its point of attachment is changed from the C to an N atom of either ring. Aromaticity of the ring offers no advantage nor is a cyclic system superior to a simple alkyl amine of any chain length. Placing a pair of -SnF₃ groups on a single molecule to form a bipodal dicationic receptor with two tetrel bonds enhances the binding, but falls short of a simple doubling. These two tetrel groups can be placed on opposite ends of an alkyl diamine chain of any length although SnF₃⁺NH₂(CH₂)nNH₂SnF₃⁺ with n between 2 and 4 seems to offer the strongest halide binding. Of the various anions tested, OH− binds most strongly: OH− > F− > Cl− > Br− > I−. The binding energy of the larger NO₃− and HCO₃− anions is more dependent upon the charge of the receptor. This pattern translates into very strong selectivity of binding one anion over another. The tetrel-bonding receptors bind far more strongly to each anion than an equivalent number of K⁺ counterions, which leads to equilibrium ratios in favor of the former of many orders of magnitude.

Implications of monomer deformation for tetrel and pnicogen bonds

Monomer rearrangement raises the interaction energy by up to 20 kcal mol⁻¹and intensifies its σ-hole by a factor of 1.5–2.9.

Halogen Bonds Formed between Substituted Imidazoliums and N Bases of Varying N-Hybridization

Heterodimers are constructed containing imidazolium and its halogen-substituted derivatives as Lewis acid. N in its sp³, sp² and sp hybridizations is taken as the electron-donating base. The halogen bond is strengthened in the Cl < Br < I order, with the H-bond generally similar in magnitude to the Br-bond. Methyl substitution on the N electron donor enhances the binding energy. Very little perturbation arises if the imidazolium is attached to a phenyl ring. The energetics are not sensitive to the hybridization of the N atom. More regular patterns appear in the individual phenomena. Charge transfer diminishes uniformly on going from amine to imine to nitrile, a pattern that is echoed by the elongation of the C-Z (Z=H, Cl, Br, I) bond in the Lewis acid. These trends are also evident in the Atoms in Molecules topography of the electron density. Molecular electrostatic potentials are not entirely consistent with energetics. Although I of the Lewis acid engages in a stronger bond than does H, it is the potential of the latter which is much more positive. The minimum on the potential of the base is most negative for the nitrile even though acetonitrile does not form the strongest bonds. Placing the systems in dichloromethane solvent reduces the binding energies but leaves intact most of the trends observed in vacuo; the same can be said of ∆G in solution.