NMR Investigations of Noncovalent Carbon Tetrel Bonds. Computational Assessment and Initial Experimental Observation

C∙∙∙O and Si∙∙∙O Tetrel Bonds: Substituent Effects and Transfer of the SiF3 Group

The tetrel bond (TB) between 1,2-benzisothiazol-3-one-2-TF3-1,1-dioxide (T = C, Si) and the O atom of pyridine-1-oxide (PO) and its derivatives (PO-X, X = H, NO2, CN, F, CH3, OH, OCH3, NH2, and Li) is examined by quantum chemical means. The Si∙∙∙O TB is quite strong, with interaction energies approaching a maximum of nearly 70 kcal/mol, while the C∙∙∙O TB is an order of magnitude weaker, with interaction energies between 2.0 and 2.6 kcal/mol. An electron-withdrawing substituent on the Lewis base weakens this TB, while an electron-donating group has the opposite effect. The SiF3 group transfers roughly halfway between the N of the acid and the O of the base without the aid of cooperative effects from a third entity.

Intermolecular interactions between the heavy-atom analogues of acetylene T2H2 (T = Si, Ge, Sn, Pb) and HCN

METHODS

The intermolecular interactions between the heavy-atom analogues of acetylene T₂H₂ (T = Si, Ge, Sn, Pb) and HCN have been investigated by theoretical calculations at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVDZ level.

RESULTS

The global energy minimum of T₂H₂ is the butterfly structure A, and another energy minimum is the planar structure B. Both structures A and B exhibit the dual behavior when binding with HCN. The various hydrogen bond (HB), dihydrogen bond (DB) and tetrel bond (TB) complexes can be found according to the MEP maps of T₂H₂. One TB and three HB complexes formed between structure A and HCN can be located for Si₂H₂ and Ge₂H₂. One TB, two HB and one DB complexes formed between structure A and HCN can be located for Sn₂H₂ and Pb₂H₂. Four TB and one HB complexes formed between structure B and HCN can be located for all the T₂H₂. The geometries and binding strengths of the complexes are compared and analyzed.

CONCLUSIONS

The interactions in these complexes are generally weak, and the interaction energies of these complexes range from -0.53 to -8.23 kcal/mol. The interaction energies of the TB complexes are larger than those of the corresponding HB and DB complexes for structure A···HCN systems. The relative binding strength of the four TB complexes exhibits different order for different structure B···HCN systems, which is consistent with the MEP maps of the isolated monomers.

The importance of tetrel bonding interactions with carbon in two arrestive iso-structural Cd(ii)-Salen coordination complexes: a comprehensive DFT overview in crystal engineering

In this article, we describe the serendipitous synthesis of two remarkable iso-structural Cd(ii)-Salen complexes [L₂Cd₄(OAc)₂(NCS)₂] in the presence of H₂L and NaSCN {where L = L¹ (N,N'-bis(3-methoxysalicylidene)-1,2-diaminopropane) and L = L² (N,N'-bis(3-methoxysalicylidene)-ethylenediamine) in 1 and 2, respectively}. The complexes were characterized by using elemental analysis, SEM-EDX, PXRD, spectroscopy, and X-ray crystallography. The X-ray crystal structure revealed that both complexes crystallize in the orthorhombic space group Pbcn, with unit cell parameters: a = 20.758(6), b = 11.022(3), c = 21.396(6) Å, V = 4895(2) ų, and Z = 4. The inner N₂O₂ and outer O₄ compartments are essentially occupied by two different Cd(ii) metal ions resulting from the de-protonated form of the ligand (L^2-) with the Cd(1) metal ions adopting a capped octahedral geometry. At the same time, Cd(2) assumes a distorted trigonal prismatic geometry. The solid-state crystal structure involves various non-covalent supramolecular interactions delineated by Hirshfeld Surface and 2D fingerprint plot analysis. Noteworthily, interesting S⋯H, O⋯H, and N⋯H contacts were observed, which have identical percentages in both complexes. The sparse tetrel bonding interactions in the complex, involving the CH₃ group, were evaluated in a new dimension of DFT. We observed this privileged bonding landscape that leads to the formation of self-assembled dimers in the crystal complexes. DFT-based MEP, RDG surface, NBO, and QTAIM/NCI plot investigation quantified such unique tetrel bonding interactions.

Influence of Substituents in the Benzene Ring on the Halogen Bond of Iodobenzene with Ammonia

The effects on the CI··N halogen bond between iodobenzene and NH3 of placing various substituents on the phenyl ring are monitored by quantum calculations. Substituents R = N(CH3)2, NH2, CH3, OCH3, COCH3, Cl, F, COH, CN, and NO2 were each placed ortho, meta, and para to the I. The depth of the σ-hole on I is deepened as R became more electron-withdrawing which is reflected in a strengthening of the halogen bond, which varied between 3.3 and 5.5 kcal/mol. In most cases, the ortho placement yields the largest perturbation, followed by meta and then para, but this trend is not universal. Parallel to these substituent effects is a progressive lengthening of the covalent C-I bond. Formation of the halogen bond reduces the NMR chemical shielding of all three nuclei directly involved in the C-I··N interaction. The deshielding of the electron donor N is most closely correlated with the strength of the bond, as is the coupling constant between I and N, so both have potential use as spectroscopic measures of halogen bond strength.

S66x8 noncovalent interactions revisited: new benchmark and performance of composite localized coupled-cluster methods

The S66x8 noncovalent interactions benchmark has been re-evaluated at the “sterling silver” level. Against this, a selection of computationally more economical alternatives has been assayed, ranging from localized CC to double hybrids and SAPT(DFT).

Frustrated Lewis Pairs Based on Carbon⋅⋅⋅Carbon+ Tetrel Bonds: A DFT Study

The ability of Triangulenium (TA⁺ ) compounds to form Frustrated Lewis Pairs (FLPs) with N-HeteroCycle Carbenes (NHCs) is analysed in this manuscript at the PBE0-D3/def2-TZVP level of theory. We have used six TA⁺ -based moieties, three presenting similar bridging groups (O (trioxo), -CH₂ (triaryl) and -NH (triaza)) and another three mixing, O, -CH₂ and NH moieties. In addition, several aryl-substituted NHCs have been used as electron donor moieties to undergo carbon⋅⋅⋅carbon⁺ tetrel bonds with the TA⁺ derivatives. More precisely, -Me,-iPr, -tBu and -Ph groups were used. Finally, we have used Bader's quantum theory of "atoms in molecules" (QTAIM) and Natural Bonding Analysis (NBO) to characterize the carbon⋅⋅⋅carbon⁺ tetrel bonds described herein. We expect the results gathered herein will be useful for further exploitation of carbon⋅⋅⋅carbon⁺ bonds in the formation of FLPs as well as to expand the current knowledge of tetrel bonds to the fields of synthetic chemistry and catalysis.

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

The tetrel bond formed by HC≡CX, H₂C=CHX, and H₃CCH₂X (X=F, Cl, Br, I) as an electron donor and TH₃F (T=C, Si, Ge) was explored by ab initio calculations. The tetrel bond formed by H₃CCH₂X is the strongest, as high as -3.45 kcal/mol for the H₃CCH₂F···GeH₃F dimer, followed by H₂C=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 H₃CCH₂X 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 H₂C=CHX complexes.

Intermolecular interactions between the heavy alkenes H2Si = TH2 (T = C, Si, Ge, Sn, Pb) and acetylene

The intermolecular interactions between the heavy alkenes H₂Si = TH₂ (T = C, Si, Ge, Sn, Pb) and C₂H₂ have been calculated at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVDZ level, and the nature of these complexes has been investigated by natural bond orbital. The four types (type-A, type-B, type-C and type-D) of complexes can be located for H₂Si = TH₂···C₂H₂ system. The complexes involving H₂Si = TH₂···C₂F₂ and H₂Si = TH₂···C₂(CN)₂ have also been examined to explore the substituent effects. Some complexes which are stable for H₂Si = TH₂···C₂H₂ system become unstable for H₂Si = TH₂···C₂F₂ or H₂Si = TH₂···C₂(CN)₂ system, while other complexes which are unstable for H₂Si = TH₂···C₂H₂ system become stable for H₂Si = TH₂···C₂F₂ or H₂Si = TH₂···C₂(CN)₂ system.

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.

Intermolecular Interactions Involving Heavy Alkenes H2Si=TH2 (T = C, Si, Ge, Sn, Pb) with H2O and HCl: Tetrel Bond and Hydrogen Bond

The intermolecular interactions between the heavy alkenes H₂Si=TH₂ (T = C, Si, Ge, Sn, Pb) and H₂O or HCl have been explored at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVDZ level. The various hydrogen bond (HB) and tetrel bond (TB) complexes can be located on the basis of molecular electrostatic potential maps of the isolated monomers. The competition between TB and HB interactions has been investigated through the relaxed potential energy surface scan. The results indicate that the HB complexes become more and more unstable relative to the TB complexes with the increase of the T atomic number, and cannot even retain as a minimum in some cases, for H₂Si=TH2···H₂O systems. In contrast, the HB complexes are generally more stable than TB complexes, and the TB complexes exhibit rather weak binding strength, for H₂Si=TH₂···HCl systems. The majority of the TB complexes formed between H₂Si=TH₂ and H₂O possesses very strong binding strength with covalent characteristics. The noncovalent TB complexes can be divided into two types on the basis of the orbital interactions: π-hole complexes, with binding angles ranging from 91 to 111°, and hybrid σ/π-hole complexes, with binding angles ranging from 130 to 165°. The interplay between different molecular interactions has been explored, and an interesting result is that the covalent TB interaction is significantly abated and becomes noncovalent because of the competitive effect.

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.

Theoretical and Spectroscopic Characterization of API-Related Azoles in Solution and in Solid State

Azoles are a family of five-membered azacyclic compounds with relevant biological and pharmacological activity. Different subclasses of azoles are defined depending on the atomic arrangement and the number of nitrogen atoms present in the ring: pyrazoles, indazoles, imidazoles, benzimidazoles, triazoles, benzotriazoles, tetrazoles and pentazoles. The complete characterization of their structure and the knowledge about their crystal packing and physical and chemical properties are of vital importance for the advancement in the design of new azole-containing drugs. In this review, we report the latest recent contributions to azole chemistry, in particular, those in which theoretical studies have been performed.

Influence of alkali substituents on the strength, properties, and nature of tetrel bond between TH3F and pyridine

Ab initio calculations have been performed for the complexes of TH₃F (T=C, Si, and Ge) with pyridine and its alkali derivatives to study the influence of an alkali substituent on the strength, properties, and nature of tetrel bond. The introduction of an alkali atom into the electron donor has a prominent enhancing effect on the strength of tetrel bond, which depends on the T atom as well as the alkali atom and its substitution position. The enhancing effect becomes larger in the C < Ge < Si, Li < Na < K, and para- < meta- < ortho- patterns. The interaction energy varies in a wide range from 2 to 40 kcal/mol. Both electrostatic and polarization including charge transfer are responsible for the enhancing effect of an alkali atom. The formation of a tetrel bond results in an elongation of F-T bond and a red shift of F-T stretch vibration, which is big enough to be detected with infrared spectroscopy. Electrostatic interaction is dominant in all complexes, while polarization is smaller or larger than dispersion in the complexes of CH₃F or TH₃F(T=Si and Ge).

A GIPAW versus GIAO-ZORA-SO study of 13C and 15N CPMAS NMR chemical shifts of aromatic and heterocyclic bromo derivatives

Theoretical simulation of NMR parameters in compounds bearing heavy atoms generally requires the application of relativistic corrections. We report herein the theoretical characterization of ¹³C and ¹⁵N CPMAS NMR of known bromo-derivative crystals by using both the GIPAW and the combined GIAO-ZORA-SO approximation methods. Several statistical analyses were performed to compare both approaches, with non-relativistic GIPAW method being more useful to predict the ¹³C and ¹⁵N chemical shifts. The problem of applying GIPAW to crystal structures showing static or dynamic crystalline disorder of the special class resulting in half-protons will be discussed in detail.

Growth Pattern, Stability, and Properties of Complexes of C2H5OH and nCO2 (n = 1-5) Molecules: A Theoretical Study

This work is dedicated to theoretically investigate the formation process of C2H5OH···nCO2 (n = 1-5) complexes and to shed light on the nature of interactions formed under the variation of CO2 concentration. It is found that CO2 molecules tend to locate around the polarized -OH group to interact with the lone pairs of the O atom. The interaction of ethanol with three CO2 molecules (C2H5OH···3CO2) induces the most stable structure in the sequence considered. The atoms in molecules (AIM), NCIplot, and natural bond orbital (NBO) analyses point out that the Oethanol···CCO2 tetrel bond overcomes hydrogen, chalcogen, and CO2···CO2 tetrel-bonded interactions and mainly contributes to the strength of C2H5OH···nCO2 (n = 1-5) complexes. All intermolecular interactions in the examined complexes are weakly noncovalent, and their positive cooperativity is evaluated to be slightly weaker than that of CO2 pure systems. SAPT2+ and molecular electrostatic potential (MEP) calculations indicate that the electrostatic force is the main factor underlying the attractive interplay in the complexes of C2H5OH and CO2.

Double Chalcogen Bonds: Crystal Engineering Stratagems via Diffraction and Multinuclear Solid-State Magnetic Resonance Spectroscopy

Group 16 chalcogens potentially provide Lewis-acidic σ-holes, which are able to form attractive supramolecular interactions with electron rich partners through chalcogen bonds. Here, a multifaceted experimental and computational study of a large series of novel chalcogen-bonded cocrystals, prepared using the principles of crystal engineering, is presented. Single-crystal X-ray diffraction studies reveal that dicyanoselenadiazole and dicyanotelluradiazole derivatives work as promising supramolecular synthons with the ability to form double chalcogen bonds with a wide range of electron donors including halides and oxygen- and nitrogen-containing heterocycles. Extensive ⁷⁷ Se and ¹²⁵ Te solid-state nuclear magnetic resonance spectroscopic investigations of cocrystals establish correlations between the NMR parameters of selenium and tellurium and the local chalcogen bonding geometry. The relationships between the electronic environment of the chalcogen bond and the ⁷⁷ Se and ¹²⁵ Te chemical shift tensors were elucidated through a natural localized molecular orbital density functional theory analysis. This systematic study of chalcogen-bond-based crystal engineering lays the foundations for the preparation of the various multicomponent systems and establishes solid-state NMR protocols to detect these interactions in powdered materials.

Direct investigation of chalcogen bonds by multinuclear solid-state magnetic resonance and vibrational spectroscopy

We report a multifaceted experimental and computational study of three self-complementary chalcogen-bond donors as well as a series of seven chalcogen bonded cocrystals.

Conformational control of N-methyl-N,N′-diacylhydrazines by noncovalent carbon bonding in solution

Unusual N(amide)⋯C–X noncovalent carbon bonding interactions stabilize the trans–cis (t–c) amide bond rotamers of N-methyl-N,N′-diacylhydrazines over the expected trans–trans (t–t) rotamers in solution.

Binding motif of ebselen in solution: chalcogen and hydrogen bonds team up

Ebselen, a compound active against SARS-CoV-2, forms a bifurcated supramolecular synthon thanks to chalcogen bond and hydrogen bond cooperation.

Theoretical Studies of IR and NMR Spectral Changes Induced by Sigma-Hole Hydrogen, Halogen, Chalcogen, Pnicogen, and Tetrel Bonds in a Model Protein Environment

Various types of σ-hole bond complexes were formed with FX, HFY, H₂FZ, and H₃FT (X = Cl, Br, I; Y = S, Se, Te; Z = P, As, Sb; T = Si, Ge, Sn) as Lewis acid. In order to examine their interactions with a protein, N-methylacetamide (NMA), a model of the peptide linkage was used as the base. These noncovalent bonds were compared by computational means with H-bonds formed by NMA with XH molecules (X = F, Cl, Br, I). In all cases, the A-F bond, which lies opposite the base and is responsible for the σ-hole on the A atom (A refers to the bridging atom), elongates and its stretching frequency undergoes a shift to the red with a band intensification, much as what occurs for the X-H bond in a H-bond (HB). Unlike the NMR shielding decrease seen in the bridging proton of a H-bond, the shielding of the bridging A atom is increased. The spectroscopic changes within NMA are similar for H-bonds and the other noncovalent bonds. The C=O bond of the amide is lengthened and its stretching frequency red-shifted and intensified. The amide II band shifts to higher frequency and undergoes a small band weakening. The NMR shielding of the O atom directly involved in the bond rises, whereas the C and N atoms both undergo a shielding decrease. The frequency shifts of the amide I and II bands of the base as well as the shielding changes of the three pertinent NMA atoms correlate well with the strength of the noncovalent bond.

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.

Nitropyridine-1-Oxides as Excellent π-Hole Donors: Interplay between σ-Hole (Halogen, Hydrogen, Triel, and Coordination Bonds) and π-Hole Interactions

In this manuscript, we use the primary source of geometrical information, i.e., Cambridge Structural Database (CSD), combined with density functional theory (DFT) calculations (PBE0-D3/def2-TZVP level of theory) to demonstrate the relevance of π-hole interactions in para-nitro substituted pyridine-1-oxides. More importantly, we show that the molecular electrostatic potential (MEP) value above and below the π–hole of the nitro group is largely influenced by the participation of the N-oxide group in several interactions like hydrogen-bonding (HB) halogen-bonding (XB), triel bonding (TrB), and finally, coordination-bonding (CB) (N⁺–O⁻ coordinated to a transition metal). The CSD search discloses that p-nitro-pyridine-1-oxide derivatives have a strong propensity to participate in π-hole interactions via the nitro group and, concurrently, N-oxide group participates in a series of interactions as electron donor. Remarkably, the DFT calculations show from strong to moderate cooperativity effects between π–hole and HB/XB/TrB/CB interactions (σ-bonding). The synergistic effects between π-hole and σ-hole bonding interactions are studied in terms of cooperativity energies, using MEP surface analysis and the Bader’s quantum theory of atoms in molecules (QTAIM).

Single-Crystal NMR Characterization of Halogen Bonds

Oxygen-17-enriched triphenylphosphine oxide and three of its halogen-bonded cocrystals featuring 1,4-diiodotetrafluorobenzene and 1,3,5-trifluoro-2,4,6-triiodobenzene as halogen bond donors have been characterized by ³¹P and ¹⁷O single-crystal NMR spectroscopy. Single-crystal NMR allows for the measurement of not only the magnitudes of various NMR interaction tensors, but also their orientations relative to the crystal lattice and therefore relative to the halogen bonds themselves. ³¹P chemical shift tensors, ¹⁷O chemical shift tensors, ¹⁷O quadrupolar coupling tensors, and ³¹P-¹⁷O indirect nuclear spin-spin (J) coupling tensors are reported here for P═O···I halogen bonds. The angular deviations in the directions of the pseudo-unique components of the ³¹P chemical shift tensors, the ¹⁷O chemical shift tensors, and the ¹⁷O quadrupolar coupling tensors from the direction of the oxygen-iodine halogen bond correlate with the deviations in linearity of the P═O···I halogen bond. There is also a clear decrease in anisotropy and an increase in asymmetry of the J(³¹P,¹⁷O) coupling tensors attributable to the formation of iodine-oxygen halogen bonds. The small but quantifiable changes in the tensors are consistent with the weak nature of these halogen bonds relative to the P═O motif. Overall, this work establishes single-crystal NMR as a novel probe of halogen bonds in solids. Analysis of the results has provided insights into the correlations between the magnitude and orientation of various NMR interaction tensors and the local geometry of the halogen bond. Gauge-including projector-augmented wave computations corroborate the experimental findings.

Tetrel Interactions from an Interacting Quantum Atoms Perspective

Tetrel bonds, the purportedly non-covalent interaction between a molecule that contains an atom of group 14 and an anion or (more generally) an atom or molecule with lone electron pairs, are under intense scrutiny. In this work, we perform an interacting quantum atoms (IQA) analysis of several simple complexes formed between an electrophilic fragment (A) (CH₃F, CH₄, CO₂, CS₂, SiO₂, SiH₃F, SiH₄, GeH₃F, GeO₂, and GeH₄) and an electron-pair-rich system (B) (NCH, NCO-, OCN-, F-, Br-, CN-, CO, CS, Kr, NC-, NH₃, OC, OH₂, SH-, and N₃-) at the aug-cc-pvtz coupled cluster singles and doubles (CCSD) level of calculation. The binding energy ( E bind AB ) is separated into intrafragment and inter-fragment components, and the latter in turn split into classical and covalent contributions. It is shown that the three terms are important in determining E bind AB , with absolute values that increase in passing from electrophilic fragments containing C, Ge, and Si. The degree of covalency between A and B is measured through the real space bond order known as the delocalization index ( δ AB ). Finally, a good linear correlation is found between δ AB and E xc AB , the exchange correlation (xc) or covalent contribution to E bind AB .

121/123Sb Nuclear Quadrupole Resonance Spectroscopy: Characterization of Non-Covalent Pnictogen Bonds and NQR Crystallography

Pnictogen (or pnicogen) bonding is an attractive interaction between the electrophilic region of group 15 elements (N, P, As, Sb, Bi) and a nucleophile. This interaction for which unique applications in catalysis have recently been uncovered continues to gain popularity. Here, we investigate a series of pnictogen-bonded cocrystals based on SbF₃ and SbCl₃, prepared via mechanochemical ball milling, with ^121/123Sb ( I = ⁵/₂ and ⁷/₂, respectively) nuclear quadrupole resonance (NQR) spectroscopy. Observed NQR frequency shifts upon cocrystallization are on the order of 0.1 to 10 MHz and are clearly diagnostic of the formation of pnictogen bonds to antimony. Further evidence for pnictogen bonding is obtained by complementary ¹³C cross-polarization magic-angle spinning solid-state NMR experiments. DFT calculations of NMR parameters as well as natural localized molecular orbital analyses support the experimental findings and elucidate the electronic origins of the experimental NQR frequency shifts. This work provides insights into the changes in the antimony quadrupolar coupling constant upon pnictogen bonding: strikingly, the decreases noted here parallel those known for hydrogen bonds, but contrast with the increases reported for halogen bonds. The utility of the observed antimony nuclear quadrupolar coupling constants in constraining structural models of cocrystals for which diffraction-based structures are unavailable, i.e., a rudimentary implementation of NQR crystallography, is established. Overall, this work offers a new approach to understand emerging classes of electrophilic interactions and to contextualize them in the broader landscape of established chemical bonding paradigms.

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.

Unexpected chalcogen bonds in tetravalent sulfur compounds

Combined CSD analysis and theoretical calculations show the importance of the polarizability in chalcogen bonding interactions. We provide evidence that the Lewis base has a preference in some cases for the σ-hole that is opposite to the more polarizable group instead of the more electron withdrawing one.

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.

S⋅⋅⋅Sn Tetrel Bonds in the Activation of Peroxisome Proliferator-Activated Receptors (PPARs) by Organotin Molecules

In this study, a PDB (Protein Data Bank) analysis and theoretical calculations (PBE0-D3/def2-TZVP level of theory) were combined to analyze the impact of S⋅⋅⋅Sn tetrel-bonding interactions in the activation mechanism of peroxisome proliferator-activated receptors (PPARs) by two organotin derivatives, triphenyltin (TPT) and tributyltin (TBT). The presence of a covalently bonded CYS285 to the organotin molecule was found to be key to enhance the σ-hole-donor ability of the tin atom, thus strengthening the tetrel-bonding interaction with a sulfur atom belonging to a vicinal methionine residue (MET364).

Quantitative Assessment of Tetrel Bonding Utilizing Vibrational Spectroscopy

A set of 35 representative neutral and charged tetrel complexes was investigated with the objective of finding the factors that influence the strength of tetrel bonding involving single bonded C, Si, and Ge donors and double bonded C or Si donors. For the first time, we introduced an intrinsic bond strength measure for tetrel bonding, derived from calculated vibrational spectroscopy data obtained at the CCSD(T)/aug-cc-pVTZ level of theory and used this measure to rationalize and order the tetrel bonds. Our study revealed that the strength of tetrel bonds is affected by several factors, such as the magnitude of the σ-hole in the tetrel atom, the negative electrostatic potential at the lone pair of the tetrel-acceptor, the positive charge at the peripheral hydrogen of the tetrel-donor, the exchange-repulsion between the lone pair orbitals of the peripheral atoms of the tetrel-donor and the heteroatom of the tetrel-acceptor, and the stabilization brought about by electron delocalization. Thus, focusing on just one or two of these factors, in particular, the σ-hole description can only lead to an incomplete picture. Tetrel bonding covers a range of -1.4 to -26 kcal/mol, which can be strengthened by substituting the peripheral ligands with electron-withdrawing substituents and by positively charged tetrel-donors or negatively charged tetrel-acceptors.

Tetrel Bonding Interactions in Perchlorinated Cyclopenta- and Cyclohexatetrelanes: A Combined DFT and CSD Study

In this manuscript, we combined DFT calculations (PBE0-D3/def2-TZVP level of theory) and a Cambridge Structural Database (CSD) survey to evaluate the ability of perchlorinated cyclopenta- and cyclohexatetrelanes in establishing tetrel bonding interactions. For this purpose, we used Tr₅Cl₁₀ and Tr₆Cl₁₂ (Tr = Si and Ge) and HCN, HF, OH⁻ and Cl⁻ as electron donor entities. Furthermore, we performed an Atoms in Molecules (AIM) analysis to further describe and characterize the interactions studied herein. A survey of crystal structures in the CSD reveals that close contacts between Si and lone-pair-possessing atoms are quite common and oriented along the extension of the covalent bond formed by the silicon with the halogen atom.

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.

A theoretical NMR study of selected benzazoles: Comparison of GIPAW and GIAO-PCM (DMSO) calculations

This paper compares the absolute shieldings obtained by gauge-including-projected-augmented-wave (GIPAW) to those obtained by gauge-invariant atomic orbital/Becke, 3-parameter, Lee-Yang-Parr (GIAO/B3LYP)/6-311++G(d,p)-polarizable continuum model (PCM, dimethyl sulfoxide) for nine benzazoles (benzimidazoles, indazoles, and benzotriazoles) recorded in the solid-state. Three nuclei were explored, ¹³ C, ¹⁵ N, and ¹⁹ F, and the gauge-including-projected-augmented-wave approach only proved better for ¹⁵ N MAS NMR.