Bifurcated chalcogen bonds: A theoretical study on the structure, strength and bonding properties

Competition Between the Two σ-Holes in the Formation of a Chalcogen Bond

A chalcogen atom Y contains two separate σ-holes when in a R₁ YR₂ molecular bonding pattern. Quantum chemical calculations consider competition between these two σ-holes to engage in a chalcogen bond (ChB) with a NH₃ base. R groups considered include F, Br, I, and tert-butyl (tBu). Also examined is the situation where the Y lies within a chalcogenazole ring, where its neighbors are C and N. Both electron-withdrawing substituents R₁ and R₂ act cooperatively to deepen the two σ-holes, but the deeper of the two holes consistently lies opposite to the more electron-withdrawing group, and is also favored to form a stronger ChB. The formation of two simultaneous ChBs in a triad requires the Y atom to act as double electron acceptor, and so anti-cooperativity weakens each bond relative to the simple dyad. This effect is such that some of the shallower σ-holes are unable to form a ChB at all when a base occupies the other site.

Probing the Directionality of S···O/N Chalcogen Bond and Its Interplay with Weak C-H···O/N/S Hydrogen Bond Using Molecular Electrostatic Potential

The directionality of the chalcogen bond (Ch-bond) formed by S and its interplay with other weak interactions have important chemical and biological implications. Here, dimers made of CH₃-S-X and O/N containing nucleophiles are studied and found to be stabilized by coexisting S···O/N and C-H···O/N interactions. Based on experimentally accessible electron density and molecular electrostatic potentials (MESPs), we showed that reciprocity between S···O/N and C-H···O/N interactions in the stability of cumulative molecular interaction (ΔE) was dependent on the strength of the σ-hole on S (V_s,max). Direct correlation between ΔE of dimers with V_s,max of S supports the electrostatic nature of the Ch-bond. Such interplay of the Ch-bond is necessary for its directionality in complex nucleophiles (carbonyl groups) with multiple electron-rich centers, which is explained using MESP. A correlation between the MESP minima in the π-region and the strength of the S-π interaction explains the directional selectivity of the Ch-bond.

Supramolecular Chalcogen-Bonded Semiconducting Nanoribbons at Work in Lighting Devices

This work describes the design and synthesis of a π-conjugated telluro[3,2-β][1]-tellurophene-based synthon that, embodying pyridyl and haloaryl chalcogen-bonding acceptors, self-assembles into nanoribbons through chalcogen bonds. The ribbons π-stack in a multi-layered architecture both in single crystals and thin films. Theoretical studies of the electronic states of chalcogen-bonded material showed the presence of a local charge density between Te and N atoms. OTFT-based charge transport measurements showed hole-transport properties for this material. Its integration as a p-type semiconductor in multi-layered CuI -based light-emitting electrochemical cells (LECs) led to a 10-fold increase in stability (38 h vs. 3 h) compared to single-layered devices. Finally, using the reference tellurotellurophene congener bearing a C-H group instead of the pyridyl N atom, a herringbone solid-state assembly is formed without charge transport features, resulting in LECs with poor stabilities (<1 h).

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.

Participation of S and Se in hydrogen and chalcogen bonds

The heavier chalcogen atoms S, Se, and Te can each participate in a range of different noncovalent interactions. They can serve as both proton donor and acceptor in H-bonds. Each atom can also act as electron acceptor in a chalcogen bond.

Chalcogen Bonding "2S-2N Squares" versus Competing Interactions: Exploring the Recognition Properties of Sulfur

Chalcogen bonding (CB) is the focus of increased attention for its applications in medicinal chemistry, materials science, and crystal engineering. However, the origin of sulfur's recognition properties remains controversial, and experimental evidence for supporting theories is still emerging. Here, a comprehensive evaluation of sulfur CB interactions is presented by investigating 2,1,3-benzothiadiazole X-ray crystallographic structures gathered from the Cambridge Structure Database (CSD), Protein Data Bank (PDB), and own laboratory findings. Through the systematic analysis of substituent effects on a subset library of over thirty benzothiadiazole derivatives, the competing interactions have been categorized into four main classes, namely 2S-2N CB square, halogen bonding (XB), S⋅⋅⋅S, and hydrogen-bonding (HB). A geometric model is employed to characterize the 2S-2N CB square motifs and discuss the role of electrostatic, dipole, and orbital contributions toward the interaction.

Thienoisoindigo-Based Semiconductor Nanowires Assembled with 2-Bromobenzaldehyde via Both Halogen and Chalcogen Bonding

We fabricated nanowires of a conjugated oligomer and applied them to organic field-effect transistors (OFETs). The supramolecular assemblies of a thienoisoindigo-based small molecular organic semiconductor (TIIG-Bz) were prepared by co-precipitation with 2-bromobenzaldehyde (2-BBA) via a combination of halogen bonding (XB) between the bromide in 2-BBA and electron-donor groups in TIIG-Bz, and chalcogen bonding (CB) between the aldehyde in 2-BBA and sulfur in TIIG-Bz. It was found that 2-BBA could be incorporated into the conjugated planes of TIIG-Bz via XB and CB pairs, thereby increasing the π − π stacking area between the conjugated planes. As a result, the driving force for one-dimensional growth of the supramolecular assemblies via π − π stacking was significantly enhanced. TIIG-Bz/2-BBA nanowires were used to fabricate OFETs, showing significantly enhanced charge transfer mobility compared to OFETs based on pure TIIG-Bz thin films and nanowires, which demonstrates the benefit of nanowire fabrication using 2-BBA.

Chalcogen bonding in synthesis, catalysis and design of materials

Chalcogen bonding is a type of noncovalent interaction in which a covalently bonded chalcogen atom (O, S, Se or Te) acts as an electrophilic species towards a nucleophilic (negative) region(s) in another or in the same molecule. In general, this interaction is strengthened by the presence of an electron-withdrawing group on the electron-acceptor chalcogen atom and upon moving down in the periodic table of elements, from O to Te. Following a short discussion of the phenomenon of chalcogen bonding, this Perspective presents some demonstrative experimental observations in which this bonding is crucial for synthetic transformations, crystal engineering, catalysis and design of materials as synthons/tectons.

The effect of anions on noncovalent interactions in model clusters of chalcogen-containing (CH3)2X (X = O, S, Se) molecules

A computational study of F⁻⋯(CH₃)₂O⋯CH₃F with F⁻ bound to the protons of the two methyl groups, found significant enhancement of the O⋯C interaction relative to the neutral (CH₃)₂O⋯CH₃F dyad.

The Origin of Chalcogen-Bonding Interactions

Favorable molecular interactions between group 16 elements have been implicated in catalysis, biological processes, and materials and medicinal chemistry. Such interactions have since become known as chalcogen bonds by analogy to hydrogen and halogen bonds. Although the prevalence and applications of chalcogen-bonding interactions continues to develop, debate still surrounds the energetic significance and physicochemical origins of this class of σ-hole interaction. Here, synthetic molecular balances were used to perform a quantitative experimental investigation of chalcogen-bonding interactions. Over 160 experimental conformational free energies were measured in 13 different solvents to examine the energetics of O···S, O···Se, S···S, O···HC, and S···HC contacts and the associated substituent and solvent effects. The strongest chalcogen-bonding interactions were found to be at least as strong as conventional H-bonds, but unlike H-bonds, surprisingly independent of the solvent. The independence of the conformational free energies on solvent polarity, polarizability, and H-bonding characteristics showed that electrostatic, solvophobic, and van der Waals dispersion forces did not account for the observed experimental trends. Instead, a quantitative relationship between the experimental conformational free energies and computed molecular orbital energies was consistent with the chalcogen-bonding interactions being dominated by n → σ* orbital delocalization between a lone pair (n) of a (thio)amide donor and the antibonding σ* orbital of an acceptor thiophene or selenophene. Interestingly, stabilization was manifested through the same acceptor molecular orbital irrespective of whether a direct chalcogen···chalcogen or chalcogen···H-C contact was made. Our results underline the importance of often-overlooked orbital delocalization effects in conformational control and molecular recognition phenomena.

The Many Facets of Chalcogen Bonding: Described by Vibrational Spectroscopy

A diverse set of 100 chalcogen-bonded complexes comprising neutral, cationic, anionic, divalent, and double bonded chalcogens has been investigated using ωB97X-D/aug-cc-pVTZ to determine geometries, binding energies, electron and energy density distributions, difference density distributions, vibrational frequencies, local stretching force constants, and associated bond strength orders. The accuracy of ωB97X-D was accessed by CCSD(T)/aug-cc-pVTZ calculations of a subset of 12 complexes and by the CCSD(T)/aug-cc-pVTZ //ωB97X-D binding energies of 95 complexes. Most of the weak chalcogen bonds can be rationalized on the basis of electrostatic contributions, but as the bond becomes stronger, covalent contributions can assume a primary role in the strength and geometry of the complexes. Covalency in chalcogen bonds involves the charge transfer from a lone pair orbital of a Lewis base into the σ* orbital of a divalent chalcogen or a π* orbital of a double bonded chalcogen. We describe for the first time a symmetric chalcogen-bonded homodimer stabilized by a charge transfer from a lone pair orbital into a π* orbital. New polymeric materials based on chalcogen bonds should take advantage of the extra stabilization granted by multiple chalcogen bonds, as is shown for 1,2,5-telluradiazole dimers.

Cooperativity between the hydrogen bonding and σ-hole interaction in linear NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 complexes (X = Cl, Br; Z = Ar, Kr): a comparative study

Quantum chemical calculations are performed to investigate the cooperativity of hydrogen bonding with halogen or aerogen bonding interactions in linear NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 clusters, where X = Cl, Br and Z = Ar, Kr. To understand the cooperativity mechanism in these systems, the corresponding binary NCX···NCH and O3Z···NCH complexes are also considered. The binding distances, interaction energies, and bonding properties of the NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 clusters are analyzed in detail. It is found that the cooperative effects in the hydrogen bonding tend to strengthen X···N and Z···N interactions. For both NCX···(NCH)n and O3Z···(NCH)n clusters, a small bond shrinkage is observed from n = 4 to n = 5, which suggests that the cooperativity effects are almost saturated in the larger clusters (n > 5). As the size of the X or Z atom is increased, the magnitude of the cooperative energy in these systems is also increased, which is mainly ascribed to changes in electrostatic potentials and orbital interactions. Our results indicate that the cooperative effects lead to a substantial change in the 14N nuclear quadrupole coupling constants of the NCH molecule.

Modulating of the pnicogen-bonding by a H⋯π interaction: An ab initio study

An ab initio study of the cooperativity in XH₂P⋯NCH⋯Z and XH₂P⋯CNH⋯Z complexes (X=F, Cl, Br, CN, NC; Z=C₂H₂,C₆H₆) connected by pnicogen-bonding and H⋯π interactions is carried out by means of MP2 computational method. A detailed analysis of the structures, interaction energies and bonding properties is performed on these systems. For each set of the complexes considered, a favorable cooperativity is observed, especially in X=F and CN complexes. However, for a given X or Z, the amount of cooperativity effects in XH₂P⋯CNH⋯Z complexes are more important than XH₂P⋯NCH⋯Z counterparts. Besides, the influence of a H⋯π interaction on a P⋯N (C) bond is more pronounced than that of a P⋯N (C) bond on a H⋯π bond. The quantum theory of atoms in molecules shows that ternary complexes have increased electron densities at their bond critical points relative to the corresponding binary systems. The results also indicate that the strength of the P⋯N(C) and H⋯π interactions increases in the presence of the solvent.

Difference in Energy between Two Distinct Types of Chalcogen Bonds Drives Regioisomerization of Binuclear (Diaminocarbene)PdII Complexes

The reaction of cis-[PdCl₂(CNXyl)₂] (Xyl = 2,6-Me₂C₆H₃) with various 1,3-thiazol- and 1,3,4-thiadiazol-2-amines in chloroform gives a mixture of two regioisomeric binuclear diaminocarbene complexes. For 1,3-thiazol-2-amines the isomeric ratio depends on the reaction conditions and kinetically (KRs) or thermodynamically (TRs) controlled regioisomers were obtained at room temperature and on heating, respectively. In CHCl₃ solutions, the isomers are subject to reversible isomerization accompanied by the cleavage of Pd-N and C-N bonds in the carbene fragment XylNCN(R)Xyl. Results of DFT calculations followed by the topological analysis of the electron density distribution within the formalism of Bader's theory (AIM method) reveal that in CHCl₃ solution the relative stability of the regioisomers (ΔG_exp = 1.2 kcal/mol; ΔG_calcd = 3.2 kcal/mol) is determined by the energy difference between two types of the intramolecular chalcogen bonds, viz. S···Cl in KRs (2.8-3.0 kcal/mol) and S···N in TRs (4.6-5.3 kcal/mol). In the case of the 1,3,4-thiadiazol-2-amines, the regioisomers are formed in approximately equal amounts and, accordingly, the energy difference between these species is only 0.1 kcal/mol in terms of ΔG_exp (ΔG_calcd = 2.1 kcal/mol). The regioisomers were characterized by elemental analyses (C, H, N), HRESI⁺-MS and FTIR, 1D (¹H, ¹³C{¹H}) and 2D (¹H,¹H-COSY, ¹H,¹H-NOESY, ¹H,¹³C-HSQC, ¹H,¹³C-HMBC) NMR spectroscopies, and structures of six complexes (three KRs and three TRs) were elucidated by single-crystal X-ray diffraction.

The mutual influence of Y···N and H···H interactions in XHY···NCH···HM complexes (X = F, Cl, Br; Y = S, Se; M = Li, Na, BeH, MgH): tuning of the chalcogen bond by dihydrogen bond interaction

The equilibrium structures, interaction energies, and bonding properties of ternary XHY···NCH···HM complexes are studied by ab initio calculations, where X = F, Cl, Br, Y = S, Se, and M = Li, Na, BeH, MgH. The ab initio calculations are carried out at the MP2/aug-cc-pVTZ level. The results indicate that all optimized Y···N and H···H binding distances in the ternary complexes are smaller than the corresponding values in the binary systems. The calculated cooperative energies (Ecoop) are between −0.20 kcal/mol in BrHS···NCH···HBeH and −3.29 kcal/mol in FHSe···NCH···HNa. For a given Y and M, the estimated Ecoop values increase as X = F > Cl > Br. In addition, the selenium-bonded complexes exibit larger Ecoop values than those of the sulfur-bonded counterparts. The cooperativity between Y···N and H···H interactions is further analyzed by quantum theory of atoms in molecules and natural bond orbital methods. Cooperative effects make an increase in the J(Y–N) and J(H–H) spin–spin coupling constants of the ternary complexes with respect to the binary systems.