Halogen Bonding: A Halogen-Centered Noncovalent Interaction Yet to Be Understood

Hunting potential SIRT2 inhibitors towards cancer treatment: Drug repurposing, molecular dynamics, and binding energy computations

The histone deacetylase group of proteins, which includes the sirtuins, contributes to several cellular processes and is intimately involved in cancer development. Sirtuins type two (SIRT2) is a constituent of the human sirtuin family, which modulates a range of biological functions and is recognized as a potential biomarker for numerous cancers. The impact of SIRT2 knockout on tumorigenesis is debated and varies with the type of tumor; nonetheless, pharmacological inhibition of SIRT2 results exclusively in growth inhibition of diverse cancer cell lines. As a result, SIRT2 regulation is thought to be a viable protein for treating cancer. Herein, the DrugBank database, containing >14,000 drug molecules, was repurposed to find potential anticancer medications that have the capacity to inhibit the SIRT2 protein utilizing in-silico techniques. In light of the experimental findings, the capability of AutoDock Vina1.1.2 software to anticipate the docking scores and poses of the SIRT2 inhibitors was assessed. SirReal2, a potential SIRT2 inhibitor, was the controller for this study. Notably, drugs with docking scores less than SiReal2 were chosen and introduced to molecular dynamics (MD) simulations, accompanied by binding affinities estimations utilizing the MM-GBSA approach. Interestingly, MM-GBSA calculations demonstrated that five drugs, namely DB11526, DB11977, DB15133, DB04739, and DB04632, revealed potential affinities as SIRT2 inhibitors exhibiting ΔGbinding less than -50.0 kcal/mol. The post-MD analyses were inspected for DB11526, DB11977, DB15133, DB04739, and DB04632, indicating excellent steadiness of these drugs bound to SIRT2 protein throughout the 200 ns MD. The ADMET features were also examined and were acceptable. These findings suggested that more attention should be paid to DB11526, DB11977, DB15133, DB04739, and DB04632 as SIRT2 inhibitors utilizing in-vitro/in-vivo assays to treat cancer disease.

Control of Pt···Pt interactions in Pt(II) complex crystals as luminescence sensors for the detection of external stimuli: recent achievements and perspectives

This review summarizes recent advances in the control of the luminescent chromic behavior of Pt(II) complex crystals, focusing on the design strategies of their molecular and packing structures. Square-planar Pt(II) complexes exhibit unique chromic luminescence owing to flexible Pt···Pt interactions, enabling their use as highly sensitive optical sensors to visualize external stimuli and environmental changes. However, the energy of the Pt···Pt interaction is relatively weak, making it difficult to precisely design and control the chromic luminescence based on the Pt···Pt interactions. This review discusses the challenging control of Pt···Pt interactions based on molecular assembly and presents design strategies using hydrogen bonding, halogen interactions, π-π stacking, and ion-pairing and counterion effects, with specific examples for each approach. Future research will focus on expanding analytical approaches, integrating these materials into sensing devices, and exploring synergies between chromic luminescence and other physical properties such as magnetism and conductivity. These advancements are expected to lead to innovative sensing technologies and significant breakthroughs in analytical chemistry and materials science.

Different Stacking Types in a Single Hybrid Cocrystal System: π···π- and π-Hole-Based Organic-Inorganic Planar Assemblies

The planar bis-chelated complex [Pd(N∩O)2] (1; N∩O = 4-MeC5H3NNC(O)NMe2) exhibits two distinct stacking modes with electron-deficient aromatics: π···π stacking with hexafluorobenzene (C6F6) versus charge-transfer π-hole interactions with 1,2,4,5-tetracyanobenzene (TCB). Cocrystallization of the complex with C6F6 or TCB yields cocrystals 1·3(C6F6) and 1·2TCB, respectively, which display different colors and stacking patterns despite similar structural motifs. Comprehensive analysis using X-ray diffraction, combined with quantum theory of atoms-in-molecules (QTAIM), an independent gradient model based on Hirshfeld partition (IGMH), extended transition state natural orbital for chemical valence theory with charge displacement function (ETS-NOCV/CDF), many-body interaction analysis, and symmetry-adapted perturbation theory (SAPT), reveals fundamentally different interaction mechanisms. In 1·3(C6F6), the stacking is primarily governed by intermolecular polarization without significant charge transfer, with dispersion forces contributing approximately 70% of the attractive energy. In contrast, 1·2TCB exhibits pronounced charge transfer (35 me) and significant inductive components alongside dispersion forces, characteristic of π-hole interactions. This dichotomy in stacking behavior provides new insights into the nature of organic-inorganic planar assemblies and demonstrates that seemingly similar structural patterns can arise from distinctly different combinations of noncovalent forces, which is essential for rational crystal engineering of hybrid materials.

The Influence of the Solvation on the Bonding of Molecular Complexes of Diatomic Halogens With Nitrogen-Containing Donors and Their Stability With Respect to the Heterolytic Halogen-Halogen Bond Splitting

In the framework of SMD approach a systematic computational study of structural, electronic and thermodynamic properties of molecular complexes of Cl2, ICl and I2 with series of N-containing Lewis bases in solvents of different polarity was carried out. Results indicate that molecular complexes of Cl2 with strong and medium-strong LB undergo spontaneous ionization in the acetonitrile solution. The increase of the solvent polarity can change the nature of interaction in X'XLB systems from molecular X'X ← LB donor-acceptor complexes to 3-center 4-electron bound X'→X+ ← LB in solvents of medium polarity and to the contact ion pairs X'→[XLB]+ in polar solvents. Thus, the controlled generation of cationic [LB∙X]+ species is possible by varying the nature of LB, varying the nature of the solvent, and varying the nature of the halogen X. Molecular Cl2 has the greatest tendency to form ionic species in polar solvents. Spontaneous ionization of molecular nσ complexes of chlorine with strong LB in medium-polar solvents (starting from OEt2, ε = 4.24) should not be neglected and single point solvation energy computations on gas phase optimized geometries are not reliable for such systems.

Facile heterolytic bond splitting of molecular chlorine upon reactions with Lewis bases: Comparison with ICl and I2

Formation of molecular complexes and subsequent heterolytic halogen-halogen bond splitting upon reactions of molecular Cl2 with nitrogen-containing Lewis bases (LB) are computationally studied at M06-2X/def2-TZVPD and for selected compounds at CCSD(T)/aug-cc-pvtz//CCSD/aug-cc-pvtz levels of theory. Obtained results are compared with data for ICl and I2 molecules. Reaction pathways indicate, that in case of Cl2∙LB complexes the activation energies for the heterolytic Cl-Cl bond splitting are lower than the activation energies of the homolytic splitting of Cl2 molecule into chlorine radicals. The heterolytic halogen splitting of molecular complexes of X2∙Py with formation of [XPy2]+…X 3 - contact ion pairs in the gas phase is slightly endothermic in case of Cl2 and I2, but slightly exothermic in the case of ICl. Formation of {[ClPy2]+…Cl 3 - }2 dimers makes the overall process exothermic. Taking into account that polar solvents favor ionic species, generation of donor-stabilized Cl+ in the presence of the Lewis bases is expected to be favorable. Thus, in polar solvents the oxidation pathway via donor-stabilized Cl+ species is viable alternative to the homolytic Cl-Cl bond breaking.

On the Versatility of the sp-, sp 2-, and sp 3-Hybridized Chalcogen-Bearing Molecules To Engage in Type I Chalcogen···Chalcogen Interactions: A Quantum Mechanical Investigation of Like···Like and Unlike Complexes

The predilection of sp-, sp 2-, and sp 3-hybridized chalcogen-bearing molecules to engage in type I chalcogen···chalcogen interactions was comparatively unveiled in like···like/unlike CY···YC, YCY···YCY, and F2Y···YF2 (where Y = O, S, and Se) complexes, respectively. Upon the optimized monomers, a potential energy surface (PES) scan was conducted to pinpoint the most favorable complexes. The energetic findings unveiled the ability of the investigated systems to engage in the interactions under study with binding energy values ranging from -0.36 to -2.33 kcal/mol. Notably, binding energies were disclosed to align in the posterior sequence; sp 2- (i.e., YCY···YCY) > sp- (i.e., CY···YC) > sp 3- (i.e., F2Y···YF2) hybridized complexes, except the like···like oxygen-bearing complexes. Instead, the highest negative binding energy values were detected for the OCO···OCO followed by those of the F2O···OF2 and CO···OC complexes. Furthermore, the like···like selenium-bearing complexes demonstrated the most considerable binding energies compared to the other investigated complexes. Remarkably, the quantum theory of atoms in molecules and noncovalent interaction index analyses revealed the highly directional and closed-shell nature of the investigated chalcogen···chalcogen interactions. Symmetry adapted-perturbation theory findings outlined the dispersion forces as the commanding forces for all the studied complexes. These observations will provide convincing justifications for the nature of chalcogens within type I chalcogen···chalcogen interactions, leading to increased progress in various domains regarding drug design and materials science.

Cooperativity and halonium transfer in the ternary NCI···CH3I···-CN halogen-bonded complex: An ab initio gas phase study

CONTEXT

The strength and nature of the two halogen bonds in the NCI···CH3I···-CN halogen-bonded ternary complex are studied in the gas phase via ab initio calculations. Different indicators of halogen bond strength were employed to examine the interactions including geometries, complexation energies, Natural Bond Order (NBO) Wiberg bond indices, and Atoms in Molecules (AIM)-based charge density topological properties. The results show that the halogen bond is strong and partly covalent in nature when CH3I donates the halogen bond, but weak and noncovalent in nature when CH3I accepts the halogen bond. Significant halogen bond cooperativity emerges in the ternary complex relative to the corresponding heterodimer complexes, NCI···CH3I and CH3I···-CN, respectively. For example, the CCSD(T) complexation energy of the ternary complex (-18.27 kcal/mol) is about twice the sum of the complexation energies of the component dimers (-9.54 kcal/mol). The halonium transfer reaction that converts the ternary complex into an equivalent one was also investigated. The electronic barrier for the halonium transfer was calculated to be 6.70 kcal/mol at the CCSD(T) level. Although the MP2 level underestimates and the MP3 overestimates the barrier, their calculated MP2.5 average barrier (6.44 kcal/mol) is close to that of the more robust CCSD(T) level. Insights on the halonium ion transfer reaction was obtained by examining the reaction energy and force profiles along the intrinsic reaction coordinate, IRC. The corresponding evolution of other properties such as bond lengths, Wiberg bond indices, and Mulliken charges provides specific insight on the extent of structural rearrangements and electronic redistribution throughout the entire IRC space.

METHODS

The MP2 method was used for geometry optimizations. Energy calculations were performed using the CCSD(T) method. The aug-cc-pVTZ basis set was employed for all atoms other than iodine for which the aug-cc-pVTZ-PP basis set was used instead.

Recent advances in catalyst design, performance, and challenges of metal-heteroatom-co-doped biochar as peroxymonosulfate activator for environmental remediation

The escalation of global water pollution due to emerging pollutants has gained significant attention. To address this issue, catalytic peroxymonosulfate (PMS) activation technology has emerged as a promising treatment approach for effectively decontaminating a wide range of pollutants. Recently, modified biochar has become an increasingly attractive as PMS activator. Metal-heteroatom-co-doped biochar (MH-BC) has emerged as a promising catalyst that can provide enhanced performance over heteroatom-doped and metal-doped biochar due to the synergism between metal and heteroatom in promoting PMS activation. Therefore, this review aims to discuss the fabrication pathways (i.e., internal vs external doping and pre-vs post-modification) and key parameters (i.e., source of precursors, synthesis methods, and synthesis conditions) affecting the performance of MH-BC as PMS activator. Subsequently, an overview of all the possible PMS activation pathways by MH-BC is provided. Subsequently, Also, the detection, identification, and quantification of several reactive species (such as, •OH, SO4•-, O2•-, 1O2, and high valent oxo species) generated in the catalytic PMS system by MH-BC are also evaluated. Lastly, the underlying challenges associated with poor stability, the lack of understanding regarding the interaction between metal and heteroatom during PMS activation and quantification of radicals in multi-ROS system are also deliberated.

Halogen Bond via an Electrophilic π-Hole on Halogen in Molecules: Does It Exist?

This study reveals a new non-covalent interaction called a π-hole halogen bond, which is directional and potentially non-linear compared to its sister analog (σ-hole halogen bond). A π-hole is shown here to be observed on the surface of halogen in halogenated molecules, which can be tempered to display the aptness to form a π-hole halogen bond with a series of electron density-rich sites (Lewis bases) hosted individually by 32 other partner molecules. The [MP2/aug-cc-pVTZ] level characteristics of the π-hole halogen bonds in 33 binary complexes obtained from the charge density approaches (quantum theory of intramolecular atoms, molecular electrostatic surface potential, independent gradient model (IGM-δginter)), intermolecular geometries and energies, and second-order hyperconjugative charge transfer analyses are discussed, which are similar to other non-covalent interactions. That a π-hole can be observed on halogen in halogenated molecules is substantiated by experimentally reported crystals documented in the Cambridge Crystal Structure Database. The importance of the π-hole halogen bond in the design and growth of chemical systems in synthetic chemistry, crystallography, and crystal engineering is yet to be fully explicated.

2D Hybrid Perovskites Employing an Organic Cation Paired with a Neutral Molecule

Two-dimensional (2D) hybrid perovskites harness the chemical and structural versatility of organic compounds. Here, we explore 2D perovskites that incorporate both a first organic component, a primary ammonium cation, and a second neutral organic module. Through the experimental examination of 42 organic pairs with a range of functional groups and organic backbones, we identify five crystallization scenarios that occur upon mixing. Only one leads to the cointercalation of the organic modules with distinct and extended interlayer spacing, which is observed with the aid of X-ray diffraction (XRD) pattern analysis combined with cross-sectional transmission electron microscopy (TEM) and elemental analysis. We present a picture in which complementary pairs, capable of forming intermolecular bonds, cocrystallize with multiple structural arrangements. These arrangements are a function of the ratio of organic content, annealing temperature, and substrate surface characteristics. We highlight how noncovalent bonds, particularly hydrogen and halogen bonding, enable the influence over the organic sublattice in hybrid halide perovskites.

Hole interactions of aerogen oxides with Lewis bases: an insight into σ-hole and lone-pair-hole interactions

σ-Hole and lone-pair (lp)-hole interactions of aerogen oxides with Lewis bases (LB) were comparatively inspected in terms of quantum mechanics calculations. The ZOn ⋯ LB complexes (where Z = Kr and Xe, n = 1, 2, 3 and 4, and LB = NH3 and NCH) showed favourable negative interaction energies. The complexation features were explained in light of σ-hole and lp-hole interactions within optimum distances lower than the sum of the respective van der Waals radii. The emerging findings outlined that σ-hole interaction energies generally enhanced according to the following order: KrO4 ⋯ < KrO⋯ < KrO3⋯ < KrO2⋯LB and XeO4⋯ < XeO⋯ < XeO2⋯ < XeO3⋯LB complexes with values ranging from -2.23 to -12.84 kcal mol-1. Lp-hole interactions with values up to -5.91 kcal mol-1 were shown. Symmetry-adapted perturbation theory findings revealed the significant contributions of electrostatic forces accounting for 50-65% of the total attractive forces within most of the ZOn⋯LB complexes. The obtained observations would be useful for the understanding of hole interactions, particularly for the aerogen oxides, with application in supramolecular chemistry and crystal engineering.

Proton tautomerism and stereoisomerism in 5-[(dimethylamino)methylidene]-4-[3/4-(trifluoromethylphenyl)amino]-1,3-thiazol-2(5H)-ones: synthesis, crystal structure and spectroscopic studies

5-[(Dimethylamino)methylidene]-4-{[3-(trifluoromethyl)phenyl]amino}-1,3-thiazol-2(5H)-one and the [4-(trifluoromethyl)phenyl]amino derivative, both C13H12F3N3OS, with the trifluoromethyl group substituted at the arene ring at the meta and para positions, were synthesized to study the structural changes associated with proton tautomerism of the amidine system. The studied compounds were found to be in the amine tautomeric form in both the solid and the liquid (dimethyl sulfoxide solutions) phase. In both isomers, the [(trifluoromethyl)phenyl]amino residue assumes a synperiplanar conformation with respect to the thiazolone system, while the 5-[(dimethylamino)methylidene] residue adopts the Z configuration. Density functional theory (DFT) calculations correctly predicted that the synperiplanar arrangement is favoured in both isomers. In the crystal, the whole independent molecule of the para compound is disordered over two alternative positions, with occupancy factors of 0.926 (3) and 0.074 (3).

Halogen Bonding Involving Isomeric Isocyanide/Nitrile Groups

2,3,5,6-Tetramethyl-1,4-diisocyanobenzene (1), 1,4-diisocyanobenzene (2), and 1,4-dicyanobenzene (3) were co-crystallized with 1,3,5-triiodotrifluorobenzene (1,3,5-FIB) to give three cocrystals, 1·1,3,5-FIB, 2·2(1,3,5-FIB), and 3·2(1,3,5-FIB), which were studied by X-ray diffraction. A common feature of the three structures is the presence of I···Cisocyanide or I···Nnitrile halogen bonds (HaBs), which occurs between an iodine σ-hole and the isocyanide C-(or the nitrile N-) atom. The diisocyanide and dinitrile cocrystals 2·2(1,3,5-FIB) and 3·2(1,3,5-FIB) are isostructural, thus providing a basis for accurate comparison of the two types of noncovalent linkages of C≡N/N≡C groups in the composition of structurally similar entities and in one crystal environment. The bonding situation was studied by a set of theoretical methods. Diisocyanides are more nucleophilic than the dinitrile and they exhibit stronger binding to 1,3,5-FIB. In all structures, the HaBs are mostly determined by the electrostatic interactions, but the dispersion and induction components also provide a noticeable contribution and make the HaBs attractive. Charge transfer has a small contribution (<5%) to the HaB and it is higher for the diisocyanide than for the dinitrile systems. At the same time, diisocyanide and dinitrile structures exhibit typical electron-donor and π-acceptor properties in relation to the HaB donor.

Catalytic activation via π-backbonding in halogen bonds

The role of halogen bonding (XB) in chemical catalysis has largely involved using XB donors as Lewis acid activators to modulate the reactivity of partner Lewis bases. We explore a more uncommon scenario, where a Lewis base modulates reactivity via a spectator halogen bond interaction. Our computational studies reveal that spectator halogen bonds may play an important role in modulating the rate of SN2 reactions. Most notably, π acceptors such as PF3 significantly decrease the barrier to substitution by decreasing electron density in the very electron rich transition state. Such π-backbonding represents an example of a heretofore unexplored situation in halogen bonding: the combination of both σ-donation and π-backdonation in this "non-covalent" interaction. The broader implications of this observation are discussed.

Homomeric chains of intermolecular bonds scaffold octahedral germanium perovskites

Perovskites with low ionic radii metal centres (for example, Ge perovskites) experience both geometrical constraints and a gain in electronic energy through distortion; for these reasons, synthetic attempts do not lead to octahedral [GeI6] perovskites, but rather, these crystallize into polar non-perovskite structures1-6. Here, inspired by the principles of supramolecular synthons7,8, we report the assembly of an organic scaffold within perovskite structures with the goal of influencing the geometric arrangement and electronic configuration of the crystal, resulting in the suppression of the lone pair expression of Ge and templating the symmetric octahedra. We find that, to produce extended homomeric non-covalent bonding, the organic motif needs to possess self-complementary properties implemented using distinct donor and acceptor sites. Compared with the non-perovskite structure, the resulting [GeI6]4- octahedra exhibit a direct bandgap with significant redshift (more than 0.5 eV, measured experimentally), 10 times lower octahedral distortion (inferred from measured single-crystal X-ray diffraction data) and 10 times higher electron and hole mobility (estimated by density functional theory). We show that the principle of this design is not limited to two-dimensional Ge perovskites; we implement it in the case of copper perovskite (also a low-radius metal centre), and we extend it to quasi-two-dimensional systems. We report photodiodes with Ge perovskites that outperform their non-octahedral and lead analogues. The construction of secondary sublattices that interlock with an inorganic framework within a crystal offers a new synthetic tool for templating hybrid lattices with controlled distortion and orbital arrangement, overcoming limitations in conventional perovskites.

Preferability of Molnupiravir, an Anti-COVID-19 Drug, toward Purine Nucleosides: A Quantum Mechanical Study

Structural aspects of molnupiravir complexed with the RNA of the SARS-CoV-2 virus have been recently resolved inside the RNA-dependent RNA polymerase (RdRp), demonstrating the interactions of molnupiravir with purine nucleosides. However, the preference of molnupiravir to interact with one purine nucleoside over another has not been clearly investigated. Herein, the complexation of molnupiravir in its active form with guanosine and adenosine was compared, using sundry density functional theory calculations. The plausible tautomeric structures of the molnupiravir drug in complex with guanosine/adenosine were minutely scrutinized. The relative energy findings outlined the favorability of amino-molnupiravir···keto-amino-guanosine and imino-molnupiravir···amino-adenosine optimized complexes. According to the interaction (Eint) and binding (Ebind) energy values, higher preferential base-pairing of molnupiravir with guanosine over the adenosine one was recognized with Eint/Ebind values of -31.16/-21.81 and -13.93/-12.83 kcal/mol, respectively. This could be interpreted by the presence of three and two hydrogen bonds within the former and latter complexes, respectively. Observable changes in the electronic properties and global indices of reactivity of the studied complexes also confirmed the preferential binding within the studied complexes. The findings from the quantum theory of atoms in molecules and the noncovalent interaction index also support the partially covalent nature of the investigated interactions. For both complexes, changes in thermodynamic parameters outlined the spontaneous, exothermic, and nonrandom states of the inspected interactions. Inspecting the solvent effect on the studied interactions outlined more observable amelioration within the water medium compared with the gas one. These results would be a durable ground for the forthcoming studies concerned with the interactions of the molnupiravir drug with purine nucleosides.

Quantifying the Impact of Halogenation on Intermolecular Interactions and Binding Modes of Aromatic Molecules

Halogenation of aromatic molecules is frequently used to modulate intermolecular interactions with ramifications for optoelectronic and mechanical properties. In this work, we accurately quantify and understand the nature of intermolecular interactions in perhalogenated benzene (PHB) clusters. Using benchmark binding energies from the fixed-node diffusion Monte Carlo (FN-DMC) method, we show that generalized Kohn-Sham semicanonical projected random phase approximation (GKS-spRPA) plus approximate exchange kernel (AKX) provides reliable interaction energies with mean absolute error (MAE) of 0.23 kcal/mol. Using the GKS-spRPA+AXK method, we quantify the interaction energies of several binding modes of PHB clusters ((C₆X₆)_n; X = F, Cl, Br, I; n = 2, 3). For a given binding mode, the interaction energies increase 3-4 times from X = F to X = I; the X-X binding modes have energies in the range of 2-4 kcal/mol, while the π-π binding mode has interaction energies in the range of 4-12 kcal/mol. SAPT-DFT-based energy decomposition analysis is then used to show that the equilibrium geometries are dictated primarily by the dispersion and exchange interactions. Finally, we test the accuracy of several dispersion-corrected density functional approximations and show that only the r2SCAN-D4 method has a low MAE and correct long-range behavior, which makes it suitable for large-scale simulations and for developing structure-function relationships of halogenated aromatic systems.

Methylammonium Tetrel Halide Perovskite Ion Pairs and Their Dimers: The Interplay between the Hydrogen-, Pnictogen- and Tetrel-Bonding Interactions

The structural stability of the extensively studied organic-inorganic hybrid methylammonium tetrel halide perovskite semiconductors, MATtX₃ (MA = CH₃NH₃⁺; Tt = Ge, Sn, Pb; X = Cl, Br, I), arises as a result of non-covalent interactions between an organic cation (CH₃NH₃⁺) and an inorganic anion (TtX₃^-). However, the basic understanding of the underlying chemical bonding interactions in these systems that link the ionic moieties together in complex configurations is still limited. In this study, ion pair models constituting the organic and inorganic ions were regarded as the repeating units of periodic crystal systems and density functional theory simulations were performed to elucidate the nature of the non-covalent interactions between them. It is demonstrated that not only the charge-assisted N-H···X and C-H···X hydrogen bonds but also the C-N···X pnictogen bonds interact to stabilize the ion pairs and to define their geometries in the gas phase. Similar interactions are also responsible for the formation of crystalline MATtX₃ in the low-temperature phase, some of which have been delineated in previous studies. In contrast, the Tt···X tetrel bonding interactions, which are hidden as coordinate bonds in the crystals, play a vital role in holding the inorganic anionic moieties (TtX₃^-) together. We have demonstrated that each Tt in each [CH₃NH₃⁺•TtX₃^-] ion pair has the capacity to donate three tetrel (σ-hole) bonds to the halides of three nearest neighbor TtX₃^- units, thus causing the emergence of an infinite array of 3D TtX₆^4- octahedra in the crystalline phase. The TtX₄^4- octahedra are corner-shared to form cage-like inorganic frameworks that host the organic cation, leading to the formation of functional tetrel halide perovskite materials that have outstanding optoelectronic properties in the solid state. We harnessed the results using the quantum theory of atoms in molecules, natural bond orbital, molecular electrostatic surface potential and independent gradient models to validate these conclusions.

Halogen bonds, chalcogen bonds, pnictogen bonds, tetrel bonds and other σ-hole interactions: a snapshot of current progress

We report here on the status of research on halogen bonds and other σ-hole interactions involving p-block elements in Lewis acidic roles, such as chalcogen bonds, pnictogen bonds and tetrel bonds. A brief overview of the available literature in this area is provided via a survey of the many review articles that address this field. Our focus has been to collect together most review articles published since 2013 to provide an easy entry into the extensive literature in this area. A snapshot of current research in the area is provided by an introduction to the virtual special issue compiled in this journal, comprising 11 articles and entitled `Halogen, chalcogen, pnictogen and tetrel bonds: structural chemistry and beyond.'

Halogen Bond to Experimentally Significant N-Heterocyclic Carbenes (I, IMe2, IiPr2, ItBu2, IPh2, IMes2, IDipp2, IAd2; I = Imidazol-2-ylidene)

The subjects of the article are halogen bonds between either XCN or XCCH (X = Cl, Br, I) and the carbene carbon atom in imidazol-2-ylidene (I) or its derivatives (IR2) with experimentally significant and systematically increased R substituents at both nitrogen atoms: methyl = Me, iso-propyl = iPr, tert-butyl = tBu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad. It is shown that the halogen bond strength increases in the order Cl < Br < I and the XCN molecule forms stronger complexes than XCCH. Of all the carbenes considered, IMes2 forms the strongest and also the shortest halogen bonds with an apogee for complex IMes2⋯ICN for which D0 = 18.71 kcal/mol and dC⋯I = 2.541 Å. In many cases, IDipp2 forms as strong halogen bonds as IMes2. Quite the opposite, although characterized by the greatest nucleophilicity, ItBu2 forms the weakest complexes (and the longest halogen bonds) if X ≠ Cl. While this finding can easily be attributed to the steric hindrance exerted by the highly branched tert-butyl groups, it appears that the presence of the four C-H⋯X hydrogen bonds may also be of importance here. Similar situation occurs in the case of complexes with IAd2.

The Tetrel Bond and Tetrel Halide Perovskite Semiconductors

The ion pairs [Cs+•TtX3-] (Tt = Pb, Sn, Ge; X = I, Br, Cl) are the building blocks of all-inorganic cesium tetrel halide perovskites in 3D, CsTtX3, that are widely regarded as blockbuster materials for optoelectronic applications such as in solar cells. The 3D structures consist of an anionic inorganic tetrel halide framework stabilized by the cesium cations (Cs+). We use computational methods to show that the geometrical connectivity between the inorganic monoanions, [TtX3-]∞, that leads to the formation of the TtX64- octahedra and the 3D inorganic perovskite architecture is the result of the joint effect of polarization and coulombic forces driven by alkali and tetrel bonds. Depending on the nature and temperature phase of these perovskite systems, the Tt···X tetrel bonds are either indistinguishable or somehow distinguishable from Tt-X coordinate bonds. The calculation of the potential on the electrostatic surface of the Tt atom in molecular [Cs+•TtX3-] provides physical insight into why the negative anions [TtX3-] attract each other when in close proximity, leading to the formation of the CsTtX3 tetrel halide perovskites in the solid state. The inter-molecular (and inter-ionic) geometries, binding energies, and charge density-based topological properties of sixteen [Cs+•TtX3-] ion pairs, as well as some selected oligomers [Cs+•PbI3-]n (n = 2, 3, 4), are discussed.

A New Motif in Halogen Bonding: Cooperative Intermolecular S-Br⋅⋅⋅O, O⋅⋅⋅F, and F⋅⋅⋅F Associations in the Crystal Packing of α,ω-Di(sulfonyl bromide) Perfluoroalkanes

We recently reported the first examples of S-Cl⋅⋅⋅O halogen bonding complemented by short F⋅⋅⋅F contacts between neighboring chains that resulted in stabilized crystals of ClSO₂ (CF₂ )₄ SO₂ Cl and ClSO₂ (CF₂ )₆ SO₂ Cl. More recently, other researchers studied our crystallographic data further using an Independent Gradient Model (IGM), and they suggested if one goes beyond IUPAC's proposed 'less than the sum of the van der Waals radii' criterion that even more noncovalent interactions between fluorine atoms on neighboring chains as well as Cl⋅⋅⋅Cl, Cl⋅⋅⋅S, O⋅⋅⋅F, and O⋅⋅⋅S attractive interactions can be found. With that said, we have prepared samples of the related BrSO₂ (CF₂ )_n SO₂ Br derivatives (where n=4, 6, 8, and others), which give rise to even stronger S-Br⋅⋅⋅O halogen bonding interactions complemented minimally by O⋅⋅⋅F and F⋅⋅⋅F intermolecular interactions as shown by X-ray crystallography and computational chemistry using IGM isosurface plots. Additional spectroscopic characterization (multinuclear NMR, FT-IR, and MS) of the disulfonyl bromide derivatives BrSO₂ (CF₂ )₄ SO₂ Br, BrSO₂ (CF₂ )₆ SO₂ Br, and BrSO₂ (CF₂ )₈ SO₂ Br has also been obtained as well as some preliminary spectroscopic evidence for BrSO₂ (CF₂ )₂ SO₂ Br and BrSO₂ CF₂ O(CF₂ )₂ OCF₂ SO₂ Br. The implication of these results toward the preparation of the corresponding disulfonyl iodides is discussed.

Boat and Chair Shaped Hexahalogen Synthons

Non-covalent halogen bonding interactions are quintessential in crystal engineering for the construction of distinctive supramolecular synthons. Here, we report the first crystalline evidences of unique boat and chair shaped cyclic hexahalogen synthons in the crystal structures of α,α,α',α',4-pentabromo-o-xylene (PBX) and α,α,α',α',4,5-hexabromo-o-xylene (HBX) respectively. Nature and stability of constituent interactions in the supramolecular synthons are scrutinized with the help of quantum-chemical calculations. Pendás' interacting quantum atoms approach confirmed the stability of Br⋅⋅⋅Br interactions leading to boat and chair shaped synthons with major contribution from exchange-correlation. Although both the molecules are achiral in nature, the packing forces guide PBX to crystallize in the chiral space group P2₁ with a helix-like orientation while HBX packs in a centrosymmetric P2₁ /n space group. The extended furcations in the pentabromo derivative construct a molecular framework consisting of macrocycles realized through halogen bonding.

Understanding the impact of halogen functional group (Br, Cl, F, OH) in amprenavir ligand of the HIV protease

We focused our attention towards the most dreadful disease that threatens the mankind of 20th century - Acquired immunodeficiency syndrome (AIDS), caused through the human immunodeficiency virus (HIV) and a sexually transmitted infection (STI). In this study, our foremost interest was to identify the potency and stability of HIV ligand- Amprenavir (APV) and its modelled functional group (Br, Cl, F, CF3, CH3, NH2) ligands through halogen and hydrogen bond contact, which will have a clear portrait on the structure activity of protein ligand interactions. This will assist chemist in synthesizing novel APV ligands, which are expected to inhibit the activity of HIV-1 protease enzyme. The binding strength of Amprenavir ligand with interacting hinge region amino acid side chains: Isoleucine (ILE 147, 150, 184), Valine (VAL 82), Alanine (ALA 28), Aspartic acid (25, 30, 125, 130) and Glycine (GLY 127, 149) were understood through interaction energy calculations at HF, B3LYP, M052X, MP2 level of theories for different basis set (6-311 G**, LANL2DZ). The present work will reveal an understandable picture about the halogen and hydrogen bond interaction that grip the contact of ligand and amino acids in the hinge region. Overall the Halogen atom (Br, Cl, F) functional groups improved the binding strength of APV in HIV protease; which provide a new novel path for the functional group preference on the ligand that enclose perfectly with the amino acid in the hinge region.Communicated by Ramaswamy H. Sarma.

Copper(II) complexes with 3,5-dihalogeno-salicylaldehydes: Synthesis, structure and interaction with DNA and albumins

Eight copper(II) complexes of 3,5-dichloro-salicyladehyde or 3,5-dibromo-salicyladehyde (3,5-diX-saloH, X = Br or Cl) were synthesized in the absence or presence of a N,N'-donor co-ligand such as 2,2'-bipyridylamine, 1,10-phenanthroline, or 2,2'-bipyridine. The resultant compounds were formulated as [Cu(3,5-diX-salo)₂(MeOH)₂] (1-2) and [Cu(3,5-diX-salo)(N,N'-donor)Cl] (3-8) and were characterized by diverse techniques. The crystal structures of three complexes were determined by single-crystal X-ray crystallography. Diverse techniques were employed in order to investigate the interaction of the complexes with calf-thymus DNA which showed intercalation as the most possible mode of their interaction. The affinity of the complexes for bovine serum albumin and human serum albumin was evaluated by fluorescence emission spectroscopy in order to calculate the binding constants which suggested a tight and reversible binding. SYNOPSIS: A series of copper(II) complexes with 3,5-dihalogen-substituted salicylaldehydes as ligands were isolated and characterized. In vitro biological studies showed the intercalation of the compounds with calf-thymus DNA and their tight and reversible binding with serum albumins.

Tetrel Bonding in Anion Recognition: A First Principles Investigation

Twenty-five molecule-anion complex systems [I₄Tt···X^-] (Tt = C, Si, Ge, Sn and Pb; X = F, Cl, Br, I and At) were examined using density functional theory (ωB97X-D) and ab initio (MP2 and CCSD) methods to demonstrate the ability of the tetrel atoms in molecular entities, I₄Tt, to recognize the halide anions when in close proximity. The tetrel bond strength for the [I₄C···X^-] series and [I₄Tt···X^-] (Tt = Si, Sn; X = I, At), was weak-to-moderate, whereas that in the remaining 16 complexes was dative tetrel bond type with very large interaction energies and short Tt···X close contact distances. The basis set superposition error corrected interaction energies calculated with the highest-level theory applied, [CCSD(T)/def2-TZVPPD], ranged from -3.0 to -112.2 kcal mol^-1. The significant variation in interaction energies was realized as a result of different levels of tetrel bonding environment between the interacting partners at the equilibrium geometries of the complex systems. Although the ωB97X-D computed intermolecular geometries and interaction energies of a majority of the [I₄Tt···X^-] complexes were close to those predicted by the highest level of theory, the MP2 results were shown to be misleading for some of these systems. To provide insight into the nature of the intermolecular chemical bonding environment in the 25 molecule-anion complexes investigated, we discussed the charge-density-based topological and isosurface features that emanated from the application of the quantum theory of atoms in molecules and independent gradient model approaches, respectively.

Enhancing Effects of the Cyano Group on the C-X∙∙∙N Hydrogen or Halogen Bond in Complexes of X-Cyanomethanes with Trimethyl Amine: CH3−n(CN)nX∙∙∙NMe3, (n = 0–3; X = H, Cl, Br, I)

In this paper, density functional theory and wave function theory calculations are carried out to investigate the strength and nature of the intermolecular C-X∙∙∙N bond interaction as a function of the number of cyano groups, CN, in the X-bond donor while maintaining the X-bond acceptor as fixed. Specifically, complexes of X-cyanomethanes with trimethyl amine CH_3−n(CN)_nX∙∙∙NMe₃ (n = 0–3; X = H, Cl, Br, I) are used as model systems. Geometrical parameters and vibrational C-X-stretching frequencies as well as interaction energies are used as relevant indicators to gauge hydrogen or halogen bond strength in the complexes. Additional characteristics of interactions that link these complexes, i.e., hydrogen or halogen bonds, are calculated with the use of the following theoretical tools: the atoms in molecules (AIM) approach, the natural bond orbital (NBO) method, and energy decomposition analysis (EDA). The results show that, for the specified X-center, the strength of C-X∙∙∙N interaction increases significantly and in a non-additive fashion with the number of CN groups. Moreover, the nature (noncovalent or partly covalent) of the interactions is revealed via the AIM approach.

Definition of the Pnictogen Bond: A Perspective

This article proposes a definition for the term “pnictogen bond” and lists its donors, acceptors, and characteristic features. These may be invoked to identify this specific subset of the inter- and intramolecular interactions formed by elements of Group 15 which possess an electrophilic site in a molecular entity.

The Pnictogen Bond, Together with Other Non-Covalent Interactions, in the Rational Design of One-, Two- and Three-Dimensional Organic-Inorganic Hybrid Metal Halide Perovskite Semiconducting Materials, and Beyond

The pnictogen bond, a somewhat overlooked supramolecular chemical synthon known since the middle of the last century, is one of the promising types of non-covalent interactions yet to be fully understood by recognizing and exploiting its properties for the rational design of novel functional materials. Its bonding modes, energy profiles, vibrational structures and charge density topologies, among others, have yet to be comprehensively delineated, both theoretically and experimentally. In this overview, attention is largely centered on the nature of nitrogen-centered pnictogen bonds found in organic-inorganic hybrid metal halide perovskites and closely related structures deposited in the Cambridge Structural Database (CSD) and the Inorganic Chemistry Structural Database (ICSD). Focusing on well-characterized structures, it is shown that it is not merely charge-assisted hydrogen bonds that stabilize the inorganic frameworks, as widely assumed and well-documented, but simultaneously nitrogen-centered pnictogen bonding, and, depending on the atomic constituents of the organic cation, other non-covalent interactions such as halogen bonding and/or tetrel bonding, are also contributors to the stabilizing of a variety of materials in the solid state. We have shown that competition between pnictogen bonding and other interactions plays an important role in determining the tilting of the MX6 (X = a halogen) octahedra of metal halide perovskites in one, two and three-dimensions. The pnictogen interactions are identified to be directional even in zero-dimensional crystals, a structural feature in many engineered ordered materials; hence an interplay between them and other non-covalent interactions drives the structure and the functional properties of perovskite materials and enabling their application in, for example, photovoltaics and optoelectronics. We have demonstrated that nitrogen in ammonium and its derivatives in many chemical systems acts as a pnictogen bond donor and contributes to conferring stability, and hence functionality, to crystalline perovskite systems. The significance of these non-covalent interactions should not be overlooked, especially when the focus is centered on the rationale design and discovery of such highly-valued materials.

The Pnictogen Bond: The Covalently Bound Arsenic Atom in Molecular Entities in Crystals as a Pnictogen Bond Donor

In chemical systems, the arsenic-centered pnictogen bond, or simply the arsenic bond, occurs when there is evidence of a net attractive interaction between the electrophilic region associated with a covalently or coordinately bound arsenic atom in a molecular entity and a nucleophile in another or the same molecular entity. It is the third member of the family of pnictogen bonds formed by the third atom of the pnictogen family, Group 15 of the periodic table, and is an inter- or intramolecular noncovalent interaction. In this overview, we present several illustrative crystal structures deposited into the Cambridge Structure Database (CSD) and the Inorganic Chemistry Structural Database (ICSD) during the last and current centuries to demonstrate that the arsenic atom in molecular entities has a significant ability to act as an electrophilic agent to make an attractive engagement with nucleophiles when in close vicinity, thereby forming σ-hole or π-hole interactions, and hence driving (in part, at least) the overall stability of the system’s crystalline phase. This overview does not include results from theoretical simulations reported by others as none of them address the signatory details of As-centered pnictogen bonds. Rather, we aimed at highlighting the interaction modes of arsenic-centered σ- and π-holes in the rationale design of crystal lattices to demonstrate that such interactions are abundant in crystalline materials, but care has to be taken to identify them as is usually done with the much more widely known noncovalent interactions in chemical systems, halogen bonding and hydrogen bonding. We also demonstrate that As-centered pnictogen bonds are usually accompanied by other primary and secondary interactions, which reinforce their occurrence and strength in most of the crystal structures illustrated. A statistical analysis of structures deposited into the CSD was performed for each interaction type As···D (D = N, O, S, Se, Te, F, Cl, Br, I, arene’s π system), thus providing insight into the typical nature of As···D interaction distances and ∠R–As···D bond angles of these interactions in crystals, where R is the remainder of the molecular entity.

Surface and Void Space Analysis of the Crystal Structures of Two Lithium Bis(pentafluoroethanesulfonyl)imide Salts

Analysis of two crystal structures of lithium bis(pentafluoroethanesulfonyl)imide is presented. Two orientations of the anion, that is a cis and trans orientation, are observed. Both structures exhibit unique interactions leading to the formation of discrete fluorous domains in the solid-state. A notable difference in the F···F interactions is seen when contrasting the two orientations wherein the trans geometry has a higher percentage of fluorine interactions than the cis orientation. The inclusion of water molecules in one of the structures also leads to the formation of a polar domain formed through a series of cyclical hydrogen bonding rings. The two structures allow for a detailed examination of the bond distances and angles involved in the formation of the two structures. Analysis of the void space in the two structures leads to the observation that the trans conformation exhibits notably higher void space as compared with the cis orientation. Hirshfeld surface analysis is used to help rationalize the interactions leading to unique changes in geometries and structure.

External Electric Field Effect on the Strength of σ-Hole Interactions: A Theoretical Perspective in Like⋯Like Carbon-Containing Complexes

For the first time, σ-hole interactions within like⋯like carbon-containing complexes were investigated, in both the absence and presence of the external electric field (EEF). The effects of the directionality and strength of the utilized EEF were thoroughly unveiled in the (F-C-F₃)₂, (F-C-H₃)₂, and (H-C-F₃)₂ complexes. In the absence of the EEF, favorable interaction energies, with negative values, are denoted for the (F-C-F₃)₂ and (H-C-F₃)₂ complexes, whereas the (F-C-H₃)₂ complex exhibits unfavorable interactions. Remarkably, the strength of the applied EEF exhibits a prominent role in turning the repulsive forces within the latter complex into attractive ones. The symmetrical nature of the considered like⋯like carbon-containing complexes eradicated the effect of directionality of the EEF. The quantum theory of atoms in molecules (QTAIM), and the noncovalent interaction (NCI) index, ensured the occurrence of the attractive forces, and also outlined the substantial contributions of the three coplanar atoms to the total strength of the studied complexes. Symmetry-adapted perturbation theory (SAPT) results show the dispersion-driven nature of the interactions.

The Stibium Bond or the Antimony-Centered Pnictogen Bond: The Covalently Bound Antimony Atom in Molecular Entities in Crystal Lattices as a Pnictogen Bond Donor

A stibium bond, i.e., a non-covalent interaction formed by covalently or coordinately bound antimony, occurs in chemical systems when there is evidence of a net attractive interaction between the electrophilic region associated with an antimony atom and a nucleophile in another, or the same molecular entity. This is a pnictogen bond and are likely formed by the elements of the pnictogen family, Group 15, of the periodic table, and is an inter- or intra-molecular non-covalent interaction. This overview describes a set of illustrative crystal systems that were stabilized (at least partially) by means of stibium bonds, together with other non-covalent interactions (such as hydrogen bonds and halogen bonds), retrieved from either the Cambridge Structure Database (CSD) or the Inorganic Crystal Structure Database (ICSD). We demonstrate that these databases contain hundreds of crystal structures of various dimensions in which covalently or coordinately bound antimony atoms in molecular entities feature positive sites that productively interact with various Lewis bases containing O, N, F, Cl, Br, and I atoms in the same or different molecular entities, leading to the formation of stibium bonds, and hence, being partially responsible for the stability of the crystals. The geometric features, pro-molecular charge density isosurface topologies, and extrema of the molecular electrostatic potential model were collectively examined in some instances to illustrate the presence of Sb-centered pnictogen bonding in the representative crystal systems considered.

σ-Hole and LP-Hole Interactions of Pnicogen···Pnicogen Homodimers under the External Electric Field Effect: A Quantum Mechanical Study

σ-Hole and lone-pair (lp)-hole interactions within σ-hole···σ-hole, σ-hole···lp-hole, and lp-hole···lp-hole configurations were comparatively investigated on the pnicogen···pnicogen homodimers (PCl3)2, for the first time, under field-free conditions and the influence of the external electric field (EEF). The electrostatic potential calculations emphasized the impressive versatility of the examined PCl3 monomers to participate in σ-hole and lp-hole pnicogen interactions. Crucially, the sizes of σ-hole and lp-hole were enlarged under the influence of the positively directed EEF and decreased in the case of reverse direction. Interestingly, the energetic quantities unveiled more favorability of the σ-hole···lp-hole configuration of the pnicogen···pnicogen homodimers, with significant negative interaction energies, than σ-hole···σ-hole and lp-hole···lp-hole configurations. Quantum theory of atoms in molecules and noncovalent interaction index analyses were adopted to elucidate the nature and origin of the considered interactions, ensuring their closed shell nature and the occurrence of attractive forces within the studied homodimers. Symmetry-adapted perturbation theory-based energy decomposition analysis alluded to the dispersion force as the main physical component beyond the occurrence of the examined interactions. The obtained findings would be considered as a fundamental underpinning for forthcoming studies pertinent to chemistry, materials science, and crystal engineering.

Type I-IV Halogen⋯Halogen Interactions: A Comparative Theoretical Study in Halobenzene⋯Halobenzene Homodimers

In the current study, unexplored type IV halogen⋯halogen interaction was thoroughly elucidated, for the first time, and compared to the well-established types I−III interactions by means of the second-order Møller−Plesset (MP2) method. For this aim, the halobenzene⋯halobenzene homodimers (where halogen = Cl, Br, and I) were designed into four different types, parodying the considered interactions. From the energetic perspective, the preference of scouted homodimers was ascribed to type II interactions (i.e., highest binding energy), whereas the lowest binding energies were discerned in type III interactions. Generally, binding energies of the studied interactions were observed to decline with the decrease in the σ-hole size in the order, C6H5I⋯IC6H5 > C6H5Br⋯BrC6H5 > C6H5Cl⋯ClC6H5 homodimers and the reverse was noticed in the case of type IV interactions. Such peculiar observations were relevant to the ample contributions of negative-belt⋯negative-belt interactions within the C6H5Cl⋯ClC6H5 homodimer. Further, type IV torsional trans → cis interconversion of C6H5X⋯XC6H5 homodimers was investigated to quantify the π⋯π contributions into the total binding energies. Evidently, the energetic features illustrated the amelioration of the considered homodimers (i.e., more negative binding energy) along the prolonged scope of torsional trans → cis interconversion. In turn, these findings outlined the efficiency of the cis configuration over the trans analog. Generally, symmetry-adapted perturbation theory-based energy decomposition analysis (SAPT-EDA) demonstrated the predominance of all the scouted homodimers by the dispersion forces. The obtained results would be beneficial for the omnipresent studies relevant to the applications of halogen bonds in the fields of materials science and crystal engineering.

The Phosphorus Bond, or the Phosphorus-Centered Pnictogen Bond: The Covalently Bound Phosphorus Atom in Molecular Entities and Crystals as a Pnictogen Bond Donor

The phosphorus bond in chemical systems, which is an inter- or intramolecular noncovalent interaction, occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a covalently or coordinately bonded phosphorus atom in a molecular entity and a nucleophile in another, or the same, molecular entity. It is the second member of the family of pnictogen bonds, formed by the second member of the pnictogen family of the periodic table. In this overview, we provide the reader with a snapshot of the nature, and possible occurrences, of phosphorus-centered pnictogen bonding in illustrative chemical crystal systems drawn from the ICSD (Inorganic Crystal Structure Database) and CSD (Cambridge Structural Database) databases, some of which date back to the latter part of the last century. The illustrative systems discussed are expected to assist as a guide to researchers in rationalizing phosphorus-centered pnictogen bonding in the rational design of molecular complexes, crystals, and materials and their subsequent characterization.

Chalcogen Bonding in the Molecular Dimers of WCh2 (Ch = S, Se, Te): On the Basic Understanding of the Local Interfacial and Interlayer Bonding Environment in 2D Layered Tungsten Dichalcogenides

Layered two-dimensional transition metal dichalcogenides and their heterostructures are of current interest, owing to the diversity of their applications in many areas of materials nanoscience and technologies. With this in mind, we have examined the three molecular dimers of the tungsten dichalcogenide series, (WCh₂)₂ (Ch = S, Se, Te), using density functional theory to provide insight into which interactions, and their specific characteristics, are responsible for the interfacial/interlayer region in the room temperature 2H phase of WCh₂ crystals. Our calculations at various levels of theory suggested that the Te···Te chalcogen bonding in (WTe₂)₂ is weak, whereas the Se···Se and S···S bonding interactions in (WSe₂)₂ and (WS₂)₂, respectively, are of the van der Waals type. The presence and character of Ch···Ch chalcogen bonding interactions in the dimers of (WCh₂)₂ are examined with a number of theoretical approaches and discussed, including charge-density-based approaches, such as the quantum theory of atoms in molecules, interaction region indicator, independent gradient model, and reduced density gradient non-covalent index approaches. The charge-density-based topological features are shown to be concordant with the results that originate from the extrema of potential on the electrostatic surfaces of WCh₂ monomers. A natural bond orbital analysis has enabled us to suggest a number of weak hyperconjugative charge transfer interactions between the interacting monomers that are responsible for the geometry of the (WCh₂)₂ dimers at equilibrium. In addition to other features, we demonstrate that there is no so-called van der Waals gap between the monolayers in two-dimensional layered transition metal tungsten dichalcogenides, which are gapless, and that the (WCh₂)₂ dimers may be prototypes for a basic understanding of the physical chemistry of the chemical bonding environments associated with the local interfacial/interlayer regions in layered 2H-WCh₂ nanoscale systems.

Chalcogen···Chalcogen Bonding in Molybdenum Disulfide, Molybdenum Diselenide and Molybdenum Ditelluride Dimers as Prototypes for a Basic Understanding of the Local Interfacial Chemical Bonding Environment in 2D Layered Transition Metal Dichalcogenides

An attempt was made, using computational methods, to understand whether the intermolecular interactions in the dimers of molybdenum dichalcogenides MoCh2 (Ch = chalcogen, element of group 16, especially S, Se and Te) and similar mixed-chalcogenide derivatives resemble the room temperature experimentally observed interactions in the interfacial regions of molybdenites and their other mixed-chalcogen derivatives. To this end, MP2(Full)/def2-TVZPPD level electronic structure calculations on nine dimer systems, including (MoCh2)2 and (MoChCh′2)2 (Ch, Ch′ = S, Se and Te), were carried out not only to demonstrate the energetic stability of these systems in the gas phase, but also to reproduce the intermolecular geometrical properties that resemble the interfacial geometries of 2D layered MoCh2 systems reported in the crystalline phase. Among the six DFT functionals (single and double hybrids) benchmarked against MP2(full), it was found that the double hybrid functional B2PLYPD3 has some ability to reproduce the intermolecular geometries and binding energies. The intermolecular geometries and binding energies of all nine dimers are discussed, together with the charge density topological aspects of the chemical bonding interactions that emerge from the application of the quantum theory of atoms in molecules (QTAIM), the isosurface topology of the reduced density gradient noncovalent index, interaction region indicator and independent gradient model (IGM) approaches. While the electrostatic surface potential model fails to explain the origin of the S···S interaction in the (MoS2)2 dimer, we show that the intermolecular bonding interactions in all nine dimers examined are a result of hyperconjugative charge transfer delocalizations between the lone-pair on (Ch/Ch′) and/or the π-orbitals of a Mo–Ch/Ch′ bond of one monomer and the dπ* anti-bonding orbitals of the same Mo–Ch/Ch′ bond in the second monomer during dimer formation, and vice versa. The HOMO–LUMO gaps calculated with the MN12-L functional were 0.9, 1.0, and 1.1 eV for MoTe2, MoSe2 and MoS2, respectively, which match very well with the solid-state theoretical (SCAN-rVV10)/experimental band gaps of 0.75/0.88, 0.90/1.09 and 0.93/1.23 eV of the corresponding systems, respectively. We observed that the gas phase dimers examined are perhaps prototypical for a basic understanding of the interfacial/inter-layer interactions in molybdenum-based dichalcogenides and their derivatives.

Crystal structure prediction of N-halide phthalimide compounds: Halogen bond synthon as a touchstone

We address crystal structure prediction problem by combining evolutionary algorithm USPEX (used to predict sets of low-energy crystal structures) and synthon approach (extracting preferable supramolecular synthons from Cambridge Structural Database,...

Exploring halogen⋯halogen interactions in supramolecular self-assemblies of BODIPY networks

In this study, the efficiency of halogen⋯halogen interactions to control supramolecular assemblies of boron dipyrromethene (BODIPY) (B1–B5) derivatives was explored.

Influence of Fluorine Substitution on Nonbonding Interactions in Selected Para-Halogeno Anilines

A series of 4-halogeno aniline derivatives was studied employing combined theoretical and experimental methods (i. e. crystal structure analysis and vibrational spectroscopies). This simplified model system was selected to shed light on the impact of fluorine substitution on the formation of noncovalent interactions such as halogen bonds (XBs) and hydrogen bonds (HBs), which are key interactions in fluorinated/halogenated drug-protein complex formation. Comparative analysis of three previously reported and five newly determined crystal structures indicated that, in most cases, 2-fluoro and 2,6-difluoro substitution of 4-X anilines increases the ability of adjacent amine to form strong N-H⋅⋅⋅N HBs. Additionally, fluorine substituents in the difluorinated derivatives are competitive and attractive HB and XB acceptors and increase the probability of halogen-halogen contacts. A peculiar observation was made for 4-iodoaniline and 2,6-difluoro-4-iodoaniline, which form distinct interaction patterns compared to the corresponding 4-Cl and 4-Br analogs. The observed intramolecular N-H⋅⋅⋅F interactions lead to additional NH bands in the FT-IR spectra.

Metal-Involving Chalcogen Bond. The Case of Platinum(II) Interaction with Se/Te-Based σ-Hole Donors

Platinum(II) complexes exhibiting an expressed d_z²-nucleophilicity of the positively charged metal centers, namely, [Pt(ppy)(acac)] (1; acacH is acetylacetone; ppyH is 2-Ph-pyridine) and [Pt(ppy)(tmhd)] (2; tmhdH is 2,2,6,6-tetramethylheptanedione-3,5), were cocrystallized with the chalcogen bond donors (4-NC₅F₄)₂Ch (Ch = Se, Te) to form two isostructural cocrystals 1·¹/₂(4-NC₅F₄)₂Ch, and 2·²/₃(4-NC₅F₄)₂Se and 2·(4-NC₅F₄)₂Te. The X-ray data for these cocrystals and appropriate theoretical DFT calculations (PBE0-D3BJ) allowed the recognition of the metal-involving chalcogen bond, namely, Ch···d_z²-Pt^II (its energy spans from -7 to -12 kcal/mol). In 1·¹/₂(4-NC₅F₄)₂Ch, Ch···d_z²-Pt^II bonding is accompanied by the C···d_z²-Pt^II interaction, representing a three-center bifurcate, whereas in 2·(4-NC₅F₄)₂Te the chalcogen bond Te···d_z²-Pt^II is purely two-centered and is stronger than that in 1·¹/₂(4-NC₅F₄)₂Ch because of more efficient orbital overlap. The association of 2 with (4-NC₅F₄)₂Te and the structure of the formed adduct in CDCl₃ solutions was studied by using ¹H, ¹³C, ¹⁹F, ¹⁹⁵Pt, ¹²⁵Te NMR, ¹⁹F-¹H HOESY, and diffusion NMR methods. The ¹⁹⁵Pt and ¹²⁵Te NMR titration and the isothermal titration calorimetry results revealed a 1:1 association of 2 with (4-NC₅F₄)₂Te.

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.

Halogen Interactions in Halogenated Oxindoles: Crystallographic and Computational Investigations of Intermolecular Interactions

X-ray structural determinations and computational studies were used to investigate halogen interactions in two halogenated oxindoles. Comparative analyses of the interaction energy and the interaction properties were carried out for Br···Br, C-H···Br, C-H···O and N-H···O interactions. Employing Møller-Plesset second-order perturbation theory (MP2) and density functional theory (DFT), the basis set superposition error (BSSE) corrected interaction energy (E_int(BSSE)) was determined using a supramolecular approach. The E_int(BSSE) results were compared with interaction energies obtained by Quantum Theory of Atoms in Molecules (QTAIM)-based methods. Reduced Density Gradient (RDG), QTAIM and Natural bond orbital (NBO) calculations provided insight into possible pathways for the intermolecular interactions examined. Comparative analysis employing the electron density at the bond critical points (BCP) and molecular electrostatic potential (MEP) showed that the interaction energies and the relative orientations of the monomers in the dimers may in part be understood in light of charge redistribution in these two compounds.