Does Chlorine in CH3Cl Behave as a Genuine Halogen Bond Donor?

A combined crystallographic and theoretical investigation of noncovalent interactions in 1,3,4-oxadiazole-2-thione-N-Mannich derivatives: in vitro bioactivity and molecular docking

Two 1,3,4-oxadiazole-2-thione-N-Mannich derivatives, specifically 5-(4-chlorophenyl)-3-[(2-trifluoromethylphenylamino)methyl]-1,3,4-oxadiazole-2(3H)-thione (1) and 5-(4-chlorophenyl)-3-[(2,5-difluorophenylamino)methyl]-1,3,4-oxadiazole-2(3H)-thione (2), were synthesized and then characterized by elemental analysis and NMR (1H and 13C) spectroscopy and the single crystal X-ray diffraction method. The formed weak intermolecular interactions in the solid-state structures of these derivatives were thoroughly investigated utilizing a variety of theoretical tools such as Hirshfeld surface analysis and quantum theory of atoms in molecules (QTAIM). Furthermore, the CLP-PIXEL and density functional theory calculations were used to study the energetics of molecular dimers. Numerous weak intermolecular interactions such as C-H⋯S/Cl/F/π interactions, a directional C-Cl⋯Cl halogen bond, π-stacking, type C-F⋯F-C contact and a short F⋯O interaction, help to stabilize the crystal structure of 1. Crystal structure 2 also stabilizes with several weak intermolecular contacts, including N-H⋯S, C-H⋯N//Cl/F interactions, a highly directional C1-Cl1⋯C(π) halogen bond and C(π)⋯C(π) interaction. In vitro antimicrobial potency of compounds 1 and 2 was assessed against various Gram-positive and Gram-negative bacterial strains and the pathogenic yeast-like Candida albicans. Both compounds showed marked activity against all tested Gram-positive bacteria and weak activity against Escherichia coli and lacked inhibitory activity against Pseudomonas aeruginosa. In addition, compounds 1 and 2 displayed good in vitro anti-proliferative activity against hepatocellular carcinoma (HepG-2) and mammary gland breast cancer (MCF-7) cancer cell lines. Molecular docking studies revealed the binding modes of title compounds at the active sites of prospective therapeutic targets.

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.

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.

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.

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.

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.

Combined crystallographic and computational investigation of the solvent disorder present in a new tipiracil hydrochloride methanol solvate–hydrate

We report a computational and crystallographic analysis of tipiracil hydrochloride (methanol solvate–hydrate). It shows a complex network of H-bonds with static discrete disorder in two out of the three symmetry-independent solvent molecules.

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

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

Anatomy of π-hole bonds: Linear systems

The list of σ-hole bonds is long and growing, encompassing both H-bonds and its closely related halogen, chalcogen, etc., sisters. These bonds rely on the asymmetric distribution of electron density, whose depletion along the extension of a covalent bond leaves a positive region of electrostatic potential from which these bonds derive their name. However, the density distributions of other molecules contain analogous positive regions that lie out of the molecular plane known as π-holes, which are likewise capable of engaging in noncovalent bonds. Quantum calculations are applied to study such π-hole bonds that involve linear molecules, whose positive region is a circular belt surrounding the molecule, rather than the more restricted area of a σ-hole. These bonds are examined in terms of their most fundamental elements arising from the spatial dispositions of their relevant molecular orbitals and the π-holes in both the total electron density and the electrostatic potential to which they lead. Systems examined comprise tetrel, chalcogen, aerogen, and triel bonds, as well as those involving group II elements, with atoms drawn from various rows of the Periodic Table. The π-hole bonds established by linear molecules tend to be weaker than those of comparable planar systems.

Yet another perspective on hole interactions

Hole interactions are known by different names depending on the key atom of the bond (halogen bond, chalcogen bond, hydrogen bond, etc.), and the geometry of the interaction (σ if in line, π if perpendicular to the Lewis acid plane). However, its origin starts with the creation of a Lewis acid by an underlying covalent bond, which forms an electrostatic depletion and a virtual antibonding orbital, which can create non-covalent interactions with Lewis bases. In this (maybe subjective) perspective, we will claim that hole interactions must be defined via the molecular orbital origin of the molecule. Under this premise we can better explore the richness of such bonding patterns. For that, we will study old, recent and new systems, trying to pinpoint some misinterpretations that are often associated with them. We will use as exemplars the triel bonds, a couple of metal complexes, a discussion on convergent σ-holes, and many cases of anti-electrostatic hole interactions.

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

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

σ- and π-Hole Interactions

Supramolecular chemistry is a very active research field that was initiated in the last century [...]

Local Vibrational Mode Analysis of π–Hole Interactions between Aryl Donors and Small Molecule Acceptors

11 aryl–lone pair and three aryl–anion π –hole interactions are investigated, along with the argon–benzene dimer and water dimer as reference compounds, utilizing the local vibrational mode theory, originally introduced by Konkoli and Cremer, to quantify the strength of the π –hole interaction in terms of a new local vibrational mode stretching force constant between the two engaged monomers, which can be conveniently used to compare different π –hole systems. Several factors have emerged which influence strength of the π –hole interactions, including aryl substituent effects, the chemical nature of atoms composing the aryl rings/ π –hole acceptors, and secondary bonding interactions between donors/acceptors. Substituent effects indirectly affect the π –hole interaction strength, where electronegative aryl-substituents moderately increase π –hole interaction strength. N-aryl members significantly increase π –hole interaction strength, and anion acceptors bind more strongly with the π –hole compared to charge neutral acceptors (lone–pair donors). Secondary bonding interactions between the acceptor and the atoms in the aryl ring can increase π –hole interaction strength, while hydrogen bonding between the π –hole acceptor/donor can significantly increase or decrease strength of the π –hole interaction depending on the directionality of hydrogen bond donation. Work is in progress expanding this research on aryl π –hole interactions to a large number of systems, including halides, CO, and OCH3− as acceptors, in order to derive a general design protocol for new members of this interesting class of compounds.