Bond paths between distant atoms do not necessarily indicate dominant interactions

σ-Hole Site-Based Interactions within Hypervalent Pnicogen, Halogen, and Aerogen-Bearing Molecules with Lewis Bases: A Comparative Study

σ-Hole site-based interactions in the trigonal bipyramidal geometrical structure of hypervalent pnicogen, halogen, and aerogen-bearing molecules with pyridine and NCH Lewis bases (LBs) were comparatively examined. In this respect, the ZF5···, XF3O2···, and AeF2O3···LB complexes (where Z = As, Sb; X = Br, I; Ae = Kr, Xe; and LB = pyridine and NCH) were investigated. The electrostatic potential (EP) analysis affirmations outlined the occurrence of σ-holes on the systems under consideration with disparate magnitudes that increased according to the following order: AeF2O3 < XF3O2 < ZF5. In line with EP outcomes, the proficiency of σ-hole site-based interactions increased as the atomic size of the central atom increased with a higher favorability for the pyridine-based complexes over NCH-based ones. The interaction energy showed the most favorable negative values of -35.97, -44.53, and -56.06 kcal/mol for the XeF2O3···, IF3O2···, and SbF5···pyridine complexes, respectively. The preferentiality pattern of the studied interactions could be explained as a consequence of (i) the dramatic rearrangement of ZF5 molecules from the trigonal bipyramid geometry to the square pyramidal one, (ii) the significant and tiny deformation energy in the case of the interaction of XF3O2 molecules with pyridine and NCH, respectively, and (iii) the absence of geometrical deformation within the AeF2O3···pyridine and ···NCH complexes other than the XeF2O3···pyridine one. Quantum theory of atoms in molecules and noncovalent interaction index findings reveal the partially covalent nature of most of the investigated interactions. Symmetry-adapted perturbation theory affirmations declared that the electrostatic component was the driving force beyond the occurrence of the considered interactions. The obtained findings will help in improving our understanding of the effect of geometrical deformation on intermolecular interactions.

Double Centrosymmetric Si···π Tetrel Bonds as New Synthons─Evidence from Crystal Structures and DFT Calculations

The crystal structure of bis((μ2-ethynylsilyloxo)-dichloro-aluminum), BEDCA, and a few related structures are characterized by the occurrence of tetrel bonds that link molecules. Particularly, centosymmetric dimers in such structures occur that are connected by two equivalent Si···π tetrel bonds. The dimer of BEDCA and dimers of other model species that similarly are linked by two equivalent Si···π tetrel bonds are analyzed theoretically. Some of the complexes calculated here are also characterized by the occurrence of triel bonds. Thus, ωB97XD/aug-cc-pVTZ calculations are performed and these DFT results are further supported by calculations with the use of other theoretical approaches: the quantum theory of atoms in molecules, QTAIM; the natural bond orbital, NBO; the energy decomposition analysis, EDA; and the noncovalent interactions method, NCI. The results show that the tetrel bonds analyzed here are rather weak, and they are not detected by the QTAIM approach; however, they are detected by other approaches, like NBO, for example. On the other hand, the triel bonds that occur in a few complexes discussed here are very strong and possess characteristics of covalent bonds.

Interplay of Hydrogen, Pnicogen, and Chalcogen Bonding in X(H2O)n=1-5 (X = NO, NO+, and NO-) Complexes: Energetics Insights via a Molecular Tailoring Approach

Nitric oxide (NO) and its redox congeners (NO+ and NO-), designated as X, play vital roles in various atmospheric and biological events. Understanding the interaction between X and water is inevitable to explain the different reactions that occur during these events. The present study is a unified attempt to explore the noncovalent interactions in microhydrated networks of X using the MP2/aug-cc-pVTZ//MP2/6-311++G(d,p) level of theory. The interactions between X and water have been probed by the molecular electrostatic potential (MESP) by exploiting the features of the most positive (Vmax) and most negative potential (Vmin) sites. The individual energy and cooperativity contributions of various types of noncovalent interactions present in X(H2O)n=1-5 complexes are estimated with the help of a molecular tailoring-based approach (MTA-based). The MTA-based analysis reveals that among various possible interactions in NO(H2O)n complexes, the water···water hydrogen bonds (HBs) are the strongest. Neutral NO can form hydrogen and pnicogen bonds (PBs) with water depending on the orientation; however, such HBs and PBs are the weakest. On the other hand, in the NO+(H2O)n complexes, the NO+···water interactions that occur through PBs are the strongest; the next one is the chalcogen bonding (CB), and the water···water HBs are the weakest. In the case of the NO-(H2O)n complexes, the HB interactions via both N and O atoms of NO- and water molecules are the strongest ones. The strength of water···water HB interactions is also seen to increase with the increase in the number of water molecules in NO-(H2O)n. The present study exemplifies the applicability of MTA-based calculations for quantifying various types of individual noncovalent interactions and their interplay in microhydrated networks of NO and its related ions.

Bader's Topological Bond Path Does Not Necessarily Indicate Stabilizing Interaction-Proof Studies Based on the Ng@[3n]cyclophane Endohedral Complexes

According to Bader's quantum theory of atoms in molecules (QTAIM), the simultaneous presence of a bond path and the corresponding bond critical point between any two atoms is both a necessary and sufficient condition for the atoms to be bonded to one another. In principle, this means that this pair of atoms should make a stabilizing contribution to the molecular system. However, the multitude of so-called counterintuitive bond paths strongly suggests that this statement is not necessarily true. Particularly 'troublesome' are endohedral complexes, in which encapsulation-enforced proximity between the trapped guest (e.g., an atom) and the host's cage system usually 'produces' many counterintuitive bond paths. In the author's opinion, the best evidence to demonstrate the repulsive nature of the intra-cage guest⋯host interaction is the use of some trapping systems containing small escape channels and then showing that the initially trapped entity spontaneously escapes outside the host's cage during geometry optimization of the initially built guest@host endohedral complex. For this purpose, a group of 24 Ng@[3n]cyclophane (3≤n≤6) endohedral complexes is used. As a result, arguments are presented showing that Bader's topological bond path does not necessarily indicate a stabilizing interaction.

The Ultrashort Spike-Ring Interaction in Substituted Iron Maiden Molecules

The in forms of molecular iron maidens are known for their unique ultrashort interaction between the apical hydrogen atom or its small substituent and the surface of the benzene ring. It is generally believed that this forced ultrashort X⋯π contact is associated with high steric hindrance, which is responsible for specific properties of iron maiden molecules. The main aim of this article is to investigate the influence of significant charge enrichment or depletion of the benzene ring on the characteristics of the ultrashort C-X⋯π contact in iron maiden molecules. For this purpose, three strongly electron-donating (-NH2) or strongly electron-withdrawing (-CN) groups were inserted into the benzene ring of in-[34,10][7]metacyclophane and its halogenated (X = F, Cl, Br) derivatives. It is shown that, despite such extremely electron-donating or electron-accepting properties, the considered iron maiden molecules surprisingly reveal quite high resistance to changes in electronic properties.

Dihydrogen Bonding-Seen through the Eyes of Vibrational Spectroscopy

In this work, we analyzed five groups of different dihydrogen bonding interactions and hydrogen clusters with an H3+ kernel utilizing the local vibrational mode theory, developed by our group, complemented with the Quantum Theory of Atoms-in-Molecules analysis to assess the strength and nature of the dihydrogen bonds in these systems. We could show that the intrinsic strength of the dihydrogen bonds investigated is primarily related to the protonic bond as opposed to the hydridic bond; thus, this should be the region of focus when designing dihydrogen bonded complexes with a particular strength. We could also show that the popular discussion of the blue/red shifts of dihydrogen bonding based on the normal mode frequencies is hampered from mode-mode coupling and that a blue/red shift discussion based on local mode frequencies is more meaningful. Based on the bond analysis of the H3+(H2)n systems, we conclude that the bond strength in these crystal-like structures makes them interesting for potential hydrogen storage applications.

Synthesis and Crystal Structure of Adamantylated 4,5,6,7-Tetrahalogeno-1H-benzimidazoles Novel Multi-Target Ligands (Potential CK2, M2 and SARS-CoV-2 Inhibitors); X-ray/DFT/QTAIM/Hirshfeld Surfaces/Molecular Docking Study

A series of new congeners, 1-[2-(1-adamantyl)ethyl]-1H-benzimidazole (AB) and 1-[2-(1-adamantyl)ethyl]-4,5,6,7-tetrahalogeno-1H-benzimidazole (Hal=Cl, Br, I; tClAB, tBrAB, tIAB), have been synthesized and studied. These novel multi-target ligands combine a benzimidazole ring known to show antitumor activity and an adamantyl moiety showing anti-influenza activity. Their crystal structures were determined by X-ray, while intermolecular interactions were studied using topological Bader's Quantum Theory of Atoms in Molecules, Hirshfeld Surfaces, CLP and PIXEL approaches. The newly synthesized compounds crystallize within two different space groups, P-1 (AB and tIAB) and P2₁/c (tClAB and tBrAB). A number of intramolecular hydrogen bonds, C-H⋯Hal (Hal=Cl, Br, I), were found in all halogen-containing congeners studied, but the intermolecular C-H⋯N hydrogen bond was detected only in AB and tIAB, while C-Hal⋯π only in tClAB and tBrAB. The interplay between C-H⋯N and C-H⋯Hal hydrogen bonds and a shift from the strong (C-H⋯Cl) to the very weak (C-H⋯I) attractive interactions upon Hal exchange, supplemented with Hal⋯Hal overlapping, determines the differences in the symmetry of crystalline packing and is crucial from the biological point of view. The hypothesis about the potential dual inhibitor role of the newly synthesized congeners was verified using molecular docking and the congeners were found to be pharmaceutically attractive as Human Casein Kinase 2, CK2, inhibitors, Membrane Matrix 2 Protein, M2, blockers and Severe Acute Respiratory Syndrome Coronavirus 2, SARS-CoV-2, inhibitors. The addition of adamantyl moiety seems to broaden and modify the therapeutic indices of the 4,5,6,7-tetrahalogeno-1H-benzimidazoles.

Description of non-covalent interactions in benzyl chalcocyanate crystals from smoothed Cromer-Mann electron density distribution functions

A well-known method to characterize non-covalent interactions consists in the topological analysis of electron density distribution (EDD) functions, complemented by the search for minima in the reduced density gradient (RDG) distributions. Here, we characterize intermolecular interactions occurring in crystals of benzyl chalcocyanate compounds through bond critical points (BCP) of the promolecular electron density (ED) built from the crystallographic Cromer-Mann parameters, at several smoothing levelst. The trajectories formed by thet-dependent BCP locations are interpreted in terms of the intermolecular interactions occurring within the crystal arrangements. Chalcogen…nitro BCPs are clearly present in the unsmoothed EDDs but are annihilated astincreases, while chalcogen…chalcogen BCPs appear and are among the only BCPs left at the highest smoothing level. The chalcogen bonds are differentiated from the other chalcogen interactions through the linear chalcogen…BCP…nitro geometry at low smoothing level and their more negative Laplacian values. The annihilation of CPs can be followed by the apparition of a RDG minimum, associated with a very weak interaction. Along the BCP trajectories, the Laplacian shows a progressive concentration of the ED in the intermolecular space within the crystals and adopts the most negative values at the shortest atom…atom separations. At the termination point of a BCP trajectory, the drastic increase of the ellipticity value illustrates the flattening of the EDD.

DFT–Assisted Structure Determination from Powder X-ray Diffraction Data of a New Zonisamide/ϵ-Caprolactam Cocrystal

The crystal structure of a new zonisamide cocrystal, an anticonvulsant drug used to treat the symptoms of epilepsy and Parkinson’s disease, with ϵ-caprolactam is reported herein. The structure has been solved by direct space methodologies from powder X-ray diffraction data. The refinement of the structure was conducted by the Rietveld method assisted by the dispersion-corrected density-functional theory (D-DFT) calculations and periodic boundary conditions. Further analysis of the structure reveals several H-bonded synthons and self–assembled dimers that have been further analyzed by DFT calculations and other computational tools such as molecular electrostatic potential (MEP) surfaces and the quantum theory of “atom-in-molecules” (QTAIM).

The physical nature of the ultrashort spike-ring interaction in iron maiden molecules

The so-called 'iron maiden' molecules belong to one of the most interesting subgroups of cyclophanes due to the presence of the ultrashort interaction between the CX apical bond and the benzene ring. This article presents an in-depth theoretical study of 16 'iron maiden' molecules, in which X = H, F, Cl or Br and the side chains are of various lengths and types: CSC, CSCC, CCC, and CCCC. It is shown that the H → F → Cl → Br substitution leads to a significant expansion of the 'iron maiden' molecule. Shorter chains lead to more pronounced effects, while insertion of sulfur atoms into the side chains lowers them. Structural changes are associated with an increase in energetic destabilization of X. Moreover, unlike for H, in the case of X = halogen, the out → in isomerization is energetically disadvantageous. The 'iron maiden' molecules are characterized by the presence of only three X⋯C_Ar bond paths. Particularly noteworthy are unusually large (even up to 32) values of the X⋯C_Ar bond ellipticity, which results from flat electron density distribution. The X⋯π interaction in each of the investigated 'iron maiden' molecule turned out to be multi-center, stabilizing and almost purely covalent in nature as indicated by the definitely dominant percentage (94.8%-101.6%) of the exchange-correlation energy. The spatial hindrance within the 'iron maiden' molecules appears to be not so much due to the X⋯π repulsion, but due to unfavorable steric interactions between X and the CC side bonds. It is also confirmed that some CH⋯HC interactions in aliphatic chains can be very weakly stabilizing.

Determining Repulsion in Cyclophane Cages

Superphane, i.e., [2.2.2.2.2.2](1,2,3,4,5,6)cyclophane, is a very convenient molecule in studying the nature of guest⋯host interactions in endohedral complexes. Nevertheless, the presence of as many as six ethylene bridges in the superphane molecule makes it practically impossible for the trapped entity to escape out of the superphane cage. Thus, in this article, I have implemented the idea of using the superphane derivatives with a reduced number of ethylene linkers, which leads to the [2n] cyclophanes where n<6. Seven such cyclophanes are then allowed to form endohedral complexes with noble gas (Ng) atoms (He, Ne, Ar, Kr). It is shown that in the vast majority of cases, the initially trapped Ng atom spontaneously escapes from the cyclophane cage, creating an exohedral complex. This is the best proof that the Ng⋯cyclophane interaction in endohedral complexes is indeed highly repulsive, i.e., destabilizing. Apart from the ‘sealed’ superphane molecule, endohedral complexes are only formed in the case of the smallest He atom. However, it has been shown that in these cases, the Ng⋯cyclophane interaction inside the cyclophane cage is nonbonding, i.e., repulsive. This highly energetically unfavorable effect causes the cyclophane molecule to ‘swell’.

Endo- and exohedral complexes of superphane with cations

Quite recently it has been shown in a previous study that superphane, that is, [2.2.2.2.2.2](1,2,3,4,5,6)cyclophane, is a very convenient molecule in the study of endohedral complexes and especially in the study of the influence of the caged entity (i.e., guest) on the structure of the host molecule. This advantage results from the presence of two parallel benzene rings joined together by six quite flexible ethylene bridges (spacers). This article examines the energetic and structural properties of endo- and exohedral complexes of superphane with the following cations: H⁺ , Li⁺ , Na⁺ , K⁺ , Be²⁺ , Mg²⁺ , Ca²⁺ , B³⁺ , Al³⁺ , Ga³⁺ . The stability of endohedral complexes has been shown to be strongly dependent on the charge and radius of the caged cation. The inclusion of the cation inside the superphane molecule causes its 'swelling', which is manifested by an increase in the distance between the benzene rings and elongations of the ring and spacer C-C bonds. In the case of exohedral complexes, three forms are investigated: with the cation above the benzene ring, with the cation interacting with the superphane window in the equatorial position, and with the cation interacting with the center of the C-C spacer bond. The first of these forms has been shown to be preferred. The cation⋯acceptor distance depends on the cation radius. Among the cations investigated, H⁺ and Be²⁺ are particularly reactive and predisposed to induce significant structural changes in the superphane molecule, forming C-H bond or C-Be-C bridges, respectively.

Does the Presence of a Bond Path Really Mean Interatomic Stabilization? The Case of the Ng@Superphane (Ng = He, Ne, Ar, and Kr) Endohedral Complexes

Using a fairly structurally flexible and, therefore, very suitable for this type of research, superphane molecule, we demonstrate that the inclusion of a noble gas atom (Ng = He, Ne, Ar, and Kr) inside it and, thus, the formation of the Ng@superphane endohedral complex, leads to its ‘swelling’. Positive values of both the binding and strain energies prove that encapsulation and in turn ‘swelling’ of the superphane molecule is energetically unfavorable and that the Ng⋯C interactions in the interior of the cage are destabilizing, i.e., repulsive. Additionally, negative Mayer Bond Orders indicate the antibonding nature of Ng⋯C contacts. This result in combination with the observed Ng⋯C bond paths shows that the presence of a bond path in the molecular graph does not necessarily prove interatomic stabilization. It is shown that the obtained conclusions do not depend on the computational methodology, i.e., the method and the basis set used. However, on the contrary, the number of bond paths may depend on the methodology. This is yet another disadvantageous finding that does not favor the treatment of bond paths on molecular graphs as indicators of chemical bonds. The Kr@superphane endohedral complex features one of the longest C–C bonds ever reported (1.753 Å).

Triel bonds within anionanion complexes

The ability of two anions to interact with one another is tested in the context of pairs of TrX₄^- homodimers, where Tr represents any of the triel atoms B, Al, Ga, In, or Tl, and X refers to a halogen substituent F, Cl, or Br. None of these pairs engage in a stable complex in the gas phase, but the situation reverses in water where the two monomers are held together by Tr⋯X triel bonds, complemented by stabilizing interactions between X atoms. Some of these bonds are quite strong, notably those involving TrF₄^-, with interaction energies surpassing 30 kcal mol^-1. Others are very much weaker, with scarcely exothermic binding energies. The highly repulsive electrostatic interactions are counteracted by large polarization energies.

Theoretical Study of N-Heterocyclic-Carbene-ZnX2 (X = H, Me, Et) Complexes

This article discusses the properties of as many as 30 carbene–ZnX2 (X = H, Me, Et) complexes featuring a zinc bond C⋯Zn. The group of carbenes is represented by imidazol-2-ylidene and its nine derivatives (labeled as IR), in which both hydrogen atoms of N-H bonds have been substituted by R groups with various spatial hindrances, from the smallest Me, iPr, tBu through Ph, Tol, and Xyl to the bulkiest Mes, Dipp, and Ad. The main goal is to study the relationship between type and size of R and X and both the strength of C⋯Zn and the torsional angle of the ZnX2 plane with respect to the plane of the imidazol-2-ylidene ring. Despite the considerable diversity of R and X, the range of dC⋯Zn is quite narrow: 2.12–2.20 Å. On the contrary, D0 is characterized by a fairly wide range of 18.5–27.4 kcal/mol. For the smallest carbenes, the ZnX2 molecule is either in the plane of the carbene or is only slightly twisted with respect to it. The twist angle becomes larger and more varied with the bulkier R. However, the value of this angle is not easy to predict because it results not only from the presence of steric effects but also from the possible presence of various interatomic interactions, such as dihydrogen bonds, tetrel bonds, agostic bonds, and hydrogen bonds. It has been shown that at least some of these interactions may have a non-negligible influence on the structure of the IR–ZnX2 complex. This fact should be taken into account in addition to the commonly discussed R⋯X steric repulsion.

Methyl groups as widespread Lewis bases in noncovalent interactions

It is well known that, under certain conditions, C(sp³) atoms behave, via their σ-hole, as Lewis acids in tetrel bonding. Here, we show that methyl groups, when bound to atoms less electronegative than carbon, can counterintuitively participate in noncovalent interactions as electron density donors. Thousands of experimental structures are found in which methyl groups behave as Lewis bases to establish alkaline, alkaline earth, triel, tetrel, pnictogen, chalcogen and halogen bonds. Theoretical calculations confirm the high directionality and significant strength of the interactions that arise from a common pattern based on the electron density holes model. Moreover, despite the absence of lone pairs, methyl groups are able to transfer charge from σ bonding orbitals into empty orbitals of the electrophile to reinforce the attractive interaction.

Molecular Tailoring Approach for the Estimation of Intramolecular Hydrogen Bond Energy

Hydrogen bonds (HBs) play a crucial role in many physicochemical and biological processes. Theoretical methods can reliably estimate the intermolecular HB energies. However, the methods for the quantification of intramolecular HB (IHB) energy available in the literature are mostly empirical or indirect and limited only to evaluating the energy of a single HB. During the past decade, the authors have developed a direct procedure for the IHB energy estimation based on the molecular tailoring approach (MTA), a fragmentation method. This MTA-based method can yield a reliable estimate of individual IHB energy in a system containing multiple H-bonds. After explaining and illustrating the methodology of MTA, we present its use for the IHB energy estimation in molecules and clusters. We also discuss the use of this method by other researchers as a standard, state-of-the-art method for estimating IHB energy as well as those of other noncovalent interactions.

Focal Point Evaluation of Energies for Tautomers and Isomers for 3-Hydroxy-2-Butenamide: Evaluation of Competing Internal Hydrogen Bonds of Types -OH…O=, -OH…N, -NH…O=, and CH…X (X=O and N)

The title compound is a small molecule with many structural variations; it can illustrate a variety of internal hydrogen bonds, among other noncovalent interactions. Here we examine structures displaying hydrogen bonding between carbonyl oxygen and hydroxyl H; between carbonyl oxygen and amino H; hydroxyl H and amino N; hydroxyl O and amino H. We also consider H-bonding in its tautomer 2-oxopropanamide. By extrapolation algorithms applied to Hartree-Fock and correlation energies as estimated in HF, MP2, and CCSD calculations using the cc-pVNZ correlation-consistent basis sets (N = 2, 3, and 4) we obtain reliable relative energies of the isomeric forms. Assuming that such energy differences may be attributed to the presence of the various types of hydrogen bonding, we attempt to infer relative strengths of types of H-bonding. The Atoms in Molecules theory of Bader and the Local Vibrational Modes analysis of Cremer and Kraka are applied to this task. Hydrogen bonds are ranked by relative strength as measured by local stretching force constants, with the stronger =O…HO- > NH…O= > -OH…N well separated from a cluster > NH…O= ≈ >NH…OH ≈ CH…O= of comparable and intermediate strength. Weaker but still significant interactions are of type CH…N which is stronger than CH…OH.

Study of Beryllium, Magnesium, and Spodium Bonds to Carbenes and Carbodiphosphoranes

The aim of this article is to present results of theoretical study on the properties of C⋯M bonds, where C is either a carbene or carbodiphosphorane carbon atom and M is an acidic center of MX2 (M = Be, Mg, Zn). Due to the rarity of theoretical data regarding the C⋯Zn bond (i.e., the zinc bond), the main focus is placed on comparing the characteristics of this interaction with C⋯Be (beryllium bond) and C⋯Mg (magnesium bond). For this purpose, theoretical studies (ωB97X-D/6-311++G(2df,2p)) have been performed for a large group of dimers formed by MX2 (X = H, F, Cl, Br, Me) and either a carbene ((NH2)2C, imidazol-2-ylidene, imidazolidin-2-ylidene, tetrahydropyrymid-2-ylidene, cyclopropenylidene) or carbodiphosphorane ((PH3)2C, (NH3)2C) molecule. The investigated dimers are characterized by a very strong charge transfer effect from either the carbene or carbodiphosphorane molecule to the MX2 one. This may even be over six times as strong as in the water dimer. According to the QTAIM and NCI method, the zinc bond is not very different than the beryllium bond, with both featuring a significant covalent contribution. However, the zinc bond should be definitely stronger if delocalization index is considered.

Towards developing a criterion to characterize non-covalent bonds: a quantum mechanical study

Chemical bonds are central to chemistry, biology, and allied fields, but still, the criterion to characterize an interaction as a non-covalent bond has not been studied rigorously.

A Critical Overview of Current Theoretical Methods of Estimating the Energy of Intramolecular Interactions

This article is probably the first such comprehensive review of theoretical methods for estimating the energy of intramolecular hydrogen bonds or other interactions that are frequently the subject of scientific research. Rather than on a plethora of numerical data, the main focus is on discussing the theoretical rationale of each method. Additionally, attention is paid to the fact that it is very often possible to use several variants of a particular method. Both of the methods themselves and their variants often give wide ranges of the obtained estimates. Attention is drawn to the fact that the applicability of a particular method may be significantly limited by various factors that disturb the reliability of the estimation, such as considerable structural changes or new important interactions in the reference system.

Understanding noncovalent bonds and their controlling forces

The fundamental underpinnings of noncovalent bonds are presented, focusing on the σ-hole interactions that are closely related to the H-bond. Different means of assessing their strength and the factors that control it are discussed. The establishment of a noncovalent bond is monitored as the two subunits are brought together, allowing the electrostatic, charge redistribution, and other effects to slowly take hold. Methods are discussed that permit prediction as to which site an approaching nucleophile will be drawn, and the maximum number of bonds around a central atom in its normal or hypervalent states is assessed. The manner in which a pair of anions can be held together despite an overall Coulombic repulsion is explained. The possibility that first-row atoms can participate in such bonds is discussed, along with the introduction of a tetrel analog of the dihydrogen bond.

Cooperativity in a cycloalkane-1,2/1,3-polyol corona: Topological hydrogen bonding in 1,2-diol motifs

A corona, consisting of 18 carbon atoms bearing 12 hydroxy groups in a continuous hydrogen-bonded chain, is built up by alternating degenerate conformations of alternating alkane-1,2-diol and 1,3-diol motifs. Geometries, proton nuclear magnetic resonance shifts and interaction energies for the dodecahydroxycyclo-octadecane and selected fragments are determined by density functional calculations at the B3LYP/6-311+G(d,p) level. Cooperative effects of O-H⋯O-H bonding are evident from the simple juxtaposition of these two motifs with a common OH group in butane-1,2,4-triol conformers. Bracketing a 1,2-diol motif with two 1,3-diol motifs in hexane-1,3,4,6-tetrol leads to a structure in which the 1,2-diol motif displays a bond critical point for hydrogen bonding. This is associated with enhancement of the shift of the hydrogen-bonded OH proton and of the corresponding H⋯O interaction energy. The full corona has a complete outer ring of O-H⋯O-H bond paths, and an inner ring of bond paths, due to C-H⋯H-C hydrogen-hydrogen bonding, which result in a central ring critical point. The topological O-H⋯O-H hydrogen bond, never seen in simple alkane-1,2-diols, is associated with cooperative enhancement of the H⋯O interaction energy, but this is not a necessary condition for a bond path: values for topological C-H⋯H-C hydrogen-hydrogen bonds can be as low as -0.4 kcal mol^-1 .

Versatility of the Cyano Group in Intermolecular Interactions

Several cyano groups are added to an alkane, alkene, and alkyne group so as to construct a Lewis acid molecule with a positive region of electrostatic potential in the area adjoining these substituents. Although each individual cyano group produces only a weak π-hole, when two or more such groups are properly situated, they can pool their π-holes into one much more intense positive region that is located midway between them. A NH3 base is attracted to this site, where it forms a strong noncovalent bond to the Lewis acid, amounting to as much as 13.6 kcal/mol. The precise nature of the bonding varies a bit from one complex to the next but typically contains a tetrel bond to the C atoms of the cyano groups or the C atoms of the linkage connecting the C≡N substituents. The placement of the cyano groups on a cyclic system like cyclopropane or cyclobutane has a mild weakening effect upon the binding. Although F is comparable to C≡N in terms of electron-withdrawing power, the replacement of cyano by F substituents substantially weakens the binding with NH3.

F-Halogen Bond: Conditions for Its Existence

The question as to whether the F atom can engage in a halogen bond (XB) remains unsettled. This issue is addressed via density functional theory calculations which pair a wide range of organic and inorganic F-acids with various sorts of Lewis bases. From an energetic perspective, perfluorinated hydrocarbons with sp, sp², or sp³ C-hybridization are unable to form an XB with an N-base, but a very weak bond can be formed if electron-withdrawing C≡N substituents are added to the acid. There is little improvement for inorganic acids O₂NF, FOF, ClF, BrF, SiF₄, or GeF₄, but F₂ is capable of a stronger XB of up to 5 kcal/mol. These results are consistent with a geometric criterion, which compares the intermolecular equilibrium distance with the sum of atomic van der Waals radii. The intensity of the σ-hole on the F atom has predictive value in that a V_s,max of at least 10-15 kcal/mol is required for XB formation. Adding a positive charge to the Lewis acid enhances the strength of any XB and even more so if the base is anionic. The acid-base interaction induces a contraction of the r(AF) covalent bond in the acid in most cases and a deshielding of the NMR signal of the F nucleus.

Interatomic exchange-correlation interaction energy from a measure of quantum theory of atoms in molecules topological bonding: A diatomic case

The potential relations between the measure of topological interatomic bonding-integrals of electron density with respect to internuclear axis over the corresponding quantum theory of atoms in molecules (QTAIM)-defined interatomic surface (IAS)-and interatomic exchange-correlation contributions from the interacting quantum atoms approach are discussed. The quantum chemical computations of 38 equilibrium diatomic systems at different levels of theory (HF, MP2, MP4SDQ, and CCSD) are invoked to support abstract considerations. Parameters of excellent correlations between IAS integrals and interatomic exchange-correlation energy are found by the optimization. The performance of these trends depends on the accuracy of the electronic correlation treatment. The resulting trends are a unique feature of equilibrium states, whereas more complicated dependencies are explored for several systems at non-equilibrium conditions. The relations of established trends with other IAS-based estimations of strength of bonding interactions between topological atoms and issues explored for multiatomic systems are briefly discussed.

Naphthyl-Fused Phosphepines: Luminescent Contorted Polycyclic P-Heterocycles

This article presents the synthesis of a new family of naphthyl-fused phosphepines through Ni-mediated C-C coupling. Interestingly, the chlorophosphine oxide intermediate shows strong resistance toward oxidation/hydrolysis owing to a combination of steric hindrance and pnictogen interactions. However, it can undergo substitution reactions under specific conditions. The optical/redox properties and the electronic structure of these new π-systems were studied experimentally (UV/Vis absorption, emission, cyclic voltammetry) and computationally (TD-DFT calculations, NICS investigation). Taking advantage of the luminescence of these derivatives, a blue-emitting OLED has been prepared, highlighting that these novel π-conjugated P-heterocycles appear to be promising building blocks for solid-state lighting applications.

Study of the influence of intermolecular interaction on classical and reverse substituent effects in para-substituted phenylboranes

The substituent effect and the reverse substituent effect in para-substituted phenylboranes and the influence of the intermolecular interaction of H⋯B type with either silane or methylsilane on the latter of these effects are extensively studied.

The ditetrel bond: noncovalent bond between neutral tetrel atoms

The ability of a tetrel atom to serve in the capacity of electron donor in a σ-hole noncovalent bond is tested by quantum calculations.

Metal-Halogen Bonding Seen through the Eyes of Vibrational Spectroscopy

Incorporation of a metal center into halogen-bonded materials can efficiently fine-tune the strength of the halogen bonds and introduce new electronic functionalities. The metal atom can adopt two possible roles: serving as halogen acceptor or polarizing the halogen donor and acceptor groups. We investigated both scenarios for 23 metal–halogen dimers trans-M(Y₂)(NC₅H₄X-3)₂ with M = Pd(II), Pt(II); Y = F, Cl, Br; X = Cl, Br, I; and NC₅H₄X-3 = 3-halopyridine. As a new tool for the quantitative assessment of metal–halogen bonding, we introduced our local vibrational mode analysis, complemented by energy and electron density analyses and electrostatic potential studies at the density functional theory (DFT) and coupled-cluster single, double, and perturbative triple excitations (CCSD(T)) levels of theory. We could for the first time quantify the various attractive contacts and their contribution to the dimer stability and clarify the special role of halogen bonding in these systems. The largest contribution to the stability of the dimers is either due to halogen bonding or nonspecific interactions. Hydrogen bonding plays only a secondary role. The metal can only act as halogen acceptor when the monomer adopts a (quasi-)planar geometry. The best strategy to accomplish this is to substitute the halo-pyridine ring with a halo-diazole ring, which considerably strengthens halogen bonding. Our findings based on the local mode analysis provide a solid platform for fine-tuning of existing and for design of new metal–halogen-bonded materials.

Bifurcated Triel Bonds—Hydrides and Halides of 1,2-Bis(Dichloroboryl)Benzene and 1,8-Bis(Dichloroboryl)Naphthalene

MP2/aug-cc-pVTZ calculations were performed on hydrides, fluorides, and chlorides of 1,8-bis(dichloroboryl)naphthalene and 1,2-bis(dichloroboryl)benzene. The theoretical analysis of BHB−, BFB−, and BClB− arrangements occurring in these complexes and classified as bifurcated triel bonds was partly based on decomposition of the energy of interaction. The latter was carried out for structures optimized using the DFT method. The complexes analyzed were characterized by a partly covalent character of the links to the hydride and halide anions; these anions strongly influenced the geometry of the complexes. The boron centers’ links for the neutral 1,8-bis(dichloroboryl)naphthalene and 1,2-bis(dichloroboryl)benzene molecules were characterized by approximately trigonal and planar configurations, while for anionic complexes, tetrahedral configurations were observed. The crystal structures of compounds related to species calculated here were found in the Cambridge Structural Database (CSD).

On the ability of pnicogen atoms to engage in both σ and π-hole complexes. Heterodimers of ZF2C6H5 (Z = P, As, Sb, Bi) and NH3

When bound to a pair of F atoms and a phenyl ring, a pyramidal pnicogen (Z) atom can form a pnicogen bond wherein an NH₃ base lies opposite one F atom. In addition to this σ-hole complex, the ZF₂C₆H₅ molecule can distort in such a way that the NH₃ approaches on the opposite side to the lone pair on Z, where there is a so-called π-hole. The interaction energies of these π-hole dimers are roughly 30 kcal/mol, much larger than the equivalent quantities for the σ-hole complexes, which are only 4-13 kcal/mol. On the other hand, this large interaction energy is countered by the considerable deformation energy required for the Lewis acid to adopt the geometry necessary to form the π-hole complex. The overall energetics of the complexation reaction are thus more exothermic for the σ-hole dimers than for the π-hole dimers.

Nonbonding interaction analyses on PVDF/[BMIM][BF4] complex system in gas and solution phase

The present study provides a detailed quantum chemical description of the physicochemical interactions between poly-vinylidene fluoride (PVDF) and 1-butyl-3-methyl-imidazolium tetrafluoro borate ([BMIM][BF₄]) ionic liquid (IL). Geometry optimization and frequency calculations are carried out for four monomer units of α- and β-PVDF, [BMIM][BF₄], and PVDF/[BMIM][BF₄] using dispersion corrected density functional theory. The effects of solvation on the systems under study are demonstrated for three polar aprotic solvents, namely tetra-hydrofuran (THF), acetone, and n,n-dimethyl formamide (DMF) using the integral equation formalism polarizable continuum model (IEFPCM). Calculated negative solvation free energy values suggest solution phase stability of the systems under study. Binding and interaction energies for β-PVDF/IL are found higher in magnitude than those for α-PVDF/IL. The nonbonding interaction phenomenon of β-PVDF/[BMIM][BF₄] is elucidated on the basis of natural bond orbital (NBO), Bader's quantum theory of atoms in molecules (QTAIM), delocalization indices, Hirshfeld surface, and reduced density gradient (RDG) analyses. Both anions and cations of ionic liquids are found to show weak van der Waals interaction with PVDF molecule but the anion ([BF₄]^-)/PVDF interaction is found to be stronger than cation ([BMIM]⁺)/PVDF interaction. Inter-unit C-H⋯F type hydrogen bonds are found to show improper (causing blue shifts in vibrational frequencies) nature. Frontier molecular orbital analysis is carried out, and different chemical parameters like electronegativity, chemical potential, chemical hardness and softness, and electrophilicity index are calculated using Koopmans' theorem. Thermochemical calculations are also performed, and the variation in different standard thermodynamic parameters with temperature is formulated. Graphical abstract (a) Hirshfeld surface mapped onto electron density and (b) NCI isosurfaces showing inter-unit interactions of β-PVDF/[BMIM][BF₄].

On the Uselessness of Bond Paths Linking Distant Atoms and on the Violation of the Concept of Privileged Exchange Channels

We refer to frequently used determinants suggesting dominant interactions between distant atoms in various dimers. First of all, we show, against the still-prevailling opinion, that, in general, bond paths have nothing in common with dominant intermolecular interactions and therefore they are useless in such cases. Quite the contrary, reliable information about dominant intermolecular interactions can be obtained by means of electrostatic potential maps, which very convincingly explain mutual orientation of molecules in a dimer. For the first time, numerous examples of interactions that violate both the concept of privileged exchange channels proposed by Pendás and his collaborators as well as inequalities obtained by Tognetti and Joubert for the β parameter related to secondary interactions are presented. The possible cause of this violation is suggested. We also show that the so-called counterintuitive bond paths result from quite natural behavior of the electron density gradient vector, i. e. searching for those areas of space that are characterized by large values of electron density or the most expanded its distributions.

Cooperative Effects in Weak Interactions: Enhancement of Tetrel Bonds by Intramolecular Hydrogen Bonds

A series of silyl and germanium complexes containing halogen atoms (fluorine and chlorine atoms) and exhibiting tetrel bonds with Lewis bases were analyzed by means of Møller-Plesset computational theory. Binding energies of germanium derivatives were more negative than silicon ones. Amongst the different Lewis bases utilized, ammonia produced the strongest tetrel bonded complexes in both Ge and Si cases, and substitution of the F atom by Cl led to stronger complexes with an ethylene backbone. However, with phenyl backbones, the fluorosilyl complexes were shown to be less stable than the chlorosilyl ones, but the opposite occurred for halogermanium complexes. In all the cases studied, the presence of a hydroxyl group enhanced the tetrel bond. That effect becomes more remarkable when an intramolecular hydrogen bond between the halogen and the hydrogen atom of the hydroxyl group takes places.