Hydride-Triel Bonds

Can H2 be Superacidic? A Computational Study of Triel-Bonded Brønsted Acids

The abundance of XIII group element compounds in science and industry together with their electron-deficient character gives rise to their influence on properties of the systems they interact with. This paper is an attempt to assess the strength, nature, and effect of formation of a triel bond on acidity. A wide set of Brønsted acids among others comprising hydrocarbons, halogen hydrides, and amines bonded with B, Al, and Ga trifluorides forming HX/TF3 was selected for the research. Various computational approaches (e.g., MP2, GFN2-xTB, SAPT2 + 3(CCD)δMP2, quantum theory of atoms in molecules analysis, and density overlap regions indicator) are used to describe the triel-bonded systems. Among other things, it was found that the electrostatics may not be the dominant contribution to the triel binding in some cases. Additionally, it was established that even weak Brønsted acids such as C2H2 or H2 may be superacidic if bonded to a Lewis acid (TF3) that is strong enough. The calculations indicate a significant covalent character of some of the studied HX/TF3 triel-bonded systems. Moreover, the effect of solvation of HX with TF3 as well as that of the reverse process on the acidity of the resulting system is thoroughly described.

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.

Static and Dynamical Quantum Studies of CX3-AlX2 and CSiX3-BX2 (X = F, Cl, Br) Complexes with Hydrocyanic Acid: Unusual Behavior of Strong π-Hole at Triel Center

The set of TX₃-TrX₂ (T = C, Si, Ge; Tr = B, Al, Ga; X = F, Cl, Br) molecules offers a rather unique opportunity to study both σ-hole and π-hole dimerization on the tetrel and triel ends, respectively. According to the molecular electrostatic potential (MEP) distribution, the π-hole extrema (acidic sites) were more intense than their σ-hole counterparts. The molecules owning the most (CX₃-AlX₂) and least (SiX₃-BX₂) intense π-holes were chosen to evaluate their capacities to attract one and two HCN molecules (Lewis bases). We discovered that the energetic characteristics of π-hole dimers severely conflict with the monomers MEP pattern since the weakest π-hole monomer forms a dimer characterized by interaction energy compared to those created by the monomers with noticeably greater power in the π-hole region. This outcome is due to the deformation of the weakest π-hole donor. Furthermore, the MEP analysis for monomers in the geometry of respective dimers revealed a "residual π-hole" site that was able to drive second ligand attachment, giving rise to the two "unusual trimers" examined further by the NCI and QTAIM analyses. Apart from them, the π-hole/π-hole and σ-hole/π-hole trimers have also been obtained throughout this study and described using energetic and geometric parameters. The SAPT approach revealed details of the bonding in one of the "unusual trimers". Finally, Born-Oppenheimer Molecular Dynamics (BOMD) simulations were carried out to investigate the time evolution of the interatomic distances of the studied complexes as well as their stability.

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.

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’.

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.

Exploring Intra- and Intermolecular Interactions in Selected N-Oxides-The Role of Hydrogen Bonds

Intra- and intermolecular interactions have been explored in selected N-oxide derivatives: 2-(N,N-dimethylamino-N-oxymethyl)-4,6-dimethylphenyl (1) and 5,5'-dibromo-3-diethylaminomethyl-2,2'-biphenol N-oxide (2). Both compounds possess intramolecular hydrogen bonding, which is classified as moderate in 1 and strong in 2, and resonance-assisted in both cases. Density Functional Theory (DFT) in its classical formulation as well as Time-Dependent extension (TD-DFT) were employed to study proton transfer phenomena. The simulations were performed in the gas phase and with implicit and explicit solvation models. The obtained structures of the studied N-oxides were compared with experimental data available. The proton reaction path was investigated using scan with an optimization method, and water molecule reorientation in the monohydrate of 1 was found upon the proton scan progress. It was found that spontaneous proton transfer phenomenon cannot occur in the electronic ground state of the compound 1. An opposite situation was noticed for the compound 2. The changes of nucleophilicity and electrophilicity upon the bridged proton migration were analyzed on the basis of Fukui functions in the case of 1. The interaction energy decomposition of dimers and microsolvation models was investigated using Symmetry-Adapted Perturbation Theory (SAPT). The simulations were performed in both phases to introduce polar environment influence on the interaction energies. The SAPT study showed rather minor role of induction in the formation of homodimers. However, it is worth noticing that the same induction term is responsible for the preference of water molecules' interaction with N-oxide hydrogen bond acceptor atoms in the microsolvation study. The Natural Bond Orbital (NBO) analysis was performed for the complexes with water to investigate the charge flow upon the polar environment introduction. Finally, the TD-DFT was applied for isolated molecules as well as for microsolvation models showing that the presence of solvent affects excited states, especially when the N-oxide acceptor atom is microsolvated.

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 Å).

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.

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.

Noncovalent Bonds through Sigma and Pi-Hole Located on the Same Molecule. Guiding Principles and Comparisons

Over the last years, scientific interest in noncovalent interactions based on the presence of electron-depleted regions called σ-holes or π-holes has markedly accelerated. Their high directionality and strength, comparable to hydrogen bonds, has been documented in many fields of modern chemistry. The current review gathers and digests recent results concerning these bonds, with a focus on those systems where both σ and π-holes are present on the same molecule. The underlying principles guiding the bonding in both sorts of interactions are discussed, and the trends that emerge from recent work offer a guide as to how one might design systems that allow multiple noncovalent bonds to occur simultaneously, or that prefer one bond type over another.

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.

On differences in substituent effects in substituted ethene and acetylene derivatives and their boranyl analogs

This article is the first attempt to present different influence of substituent effects on double and triple bonds and, conversely, to present the impact of these bonds on the electronic structure of substituents. For this purpose, quantum-mechanical calculations were made for X-substituted derivatives of ethene and acetylene with 27 diverse substituents representing a wide spectrum of electronic properties, from strongly electron-accepting to strongly electron-donating ones. In addition to these systems, their boranyl derivatives are also investigated. It turns out that the Hammett substituent constants do not correctly describe changes in the CC bond length in any of the considered family of systems. However, the relationships with the CB bond length are significantly better. It is shown that the triple bond in acetylene derivatives is much more resistant to external perturbations than the double bond in the analogs containing an ethene unit. As a consequence, in acetylene derivatives, the substituent effects on CC bond length are about half of the substituent effects in ethene derivatives. We suggest that the observed lack of a clear linear correlation between the length of the CC triple bond in acetylene derivatives and the value of electron density on this bond is due to the disturbing additional interaction between the structure of the X substituent in the xy plane and the π bond being in the same plane in the acetylene unit—on the contrary, this interaction is not possible in ethene analogs.

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.

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.

Influence of N-Base and O-Base Hybridization on Triel Bonds

The complexes of TrR₃ (Tr = B and Al; R = H, F, Cl, and Br) with three N-bases (NH₃, CH₂NH, and HCN) and three O-bases (CH₃OH, H₂CO, and CO) are utilized to explore the hybridization effect of N and O atoms on the strength, properties, and nature of the triel bond. The sp-hybridized O and N atoms form the weakest triel bond, followed by the sp²-hybridized O atom or the sp³-hybridized N atom, and the sp³-hybridized O atom or the sp²-hybridized N atom engages in the strongest triel bond. The hybridization effect is also related to the substituent of TrR₃. Most complexes are dominated by electrostatic, with increasing polarization contribution from sp to sp² to sp³. Although the CO oxygen engages in a weaker triel bond, its carbon atom is a better electron donor and the interaction energy even amounts to -37 kcal/mol in the BH₃ complex.

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.

Competition between Intra and Intermolecular Triel Bonds. Complexes between Naphthalene Derivatives and Neutral or Anionic Lewis Bases

: A TrF2 group (Tr = B, Al, Ga, In, Tl) is placed on one of the α positions of naphthalene, and its ability to engage in a triel bond (TrB) with a weak (NCH) and strong (NC-) nucleophile is assessed by ab initio calculations. As a competitor, an NH2 group is placed on the neighboring Cα, from which point it forms an intramolecular TrB with the TrF2 group. The latter internal TrB reduces the intensity of the π-hole on the Tr atom, decreasing its ability to engage in a second external TrB. The intermolecular TrB is weakened by a factor of about two for the smaller Tr atoms but is less severe for the larger Tl. The external TrB can be quite strong nonetheless; it varies from a minimum of 8 kcal/mol for the weak NCH base, up to as much as 70 kcal/mol for CN-. Likewise, the appearance of an external TrB to a strong base like CN- lessens the ability of the Tr to engage in an internal TrB, to the point where such an intramolecular TrB becomes questionable.

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.

Physical nature of silane⋯carbene dimers revealed by state-of-the-art ab initio calculations

Using the SAPT2 + 3(CCD)δMP2 method in complete basis set (CBS) limit, it is shown that the interactions in the recently studied silane⋯carbene dimers are mainly dispersive in nature. Consequently, slow convergence of dispersion energy also forces slow convergence of the interaction energy. Therefore, obtaining very accurate values requires extrapolation of the correlation part to the CBS limit. The most accurate values obtained at the CCSD(T)/CBS level of theory show that the studied silane⋯carbene dimers are rather weakly bound, with interaction energies ranging from about -1.9 to -1.3 kcal/mol. Comparing to CCSD(T)/CBS, it will be shown that SCS-MP2 and MP2C methods clearly underestimate and methods based on SAPT2+ and having some third-order corrections, as well as the MP2 method, overestimate values of interaction energies. Popular SAPT(DFT) method performs better than SCS-MP2 and MP2C; nevertheless, underestimation is still considerable. The underestimation is slightly quenched if third-order dispersion energy and its exchange counterpart is added to the SAPT(DFT). The closest value of CCSD(T)/CBS has been given by the SAPT2 + (3)(CCD)δMP2 method in quadruple-ζ basis set. © 2019 Wiley Periodicals, Inc.

Triel bonds in RZH2···NH3: hybridization, solvation, and substitution

The influence of hybridization, substitution, and solvation on the triel bond has been investigated in the complexes of RZH₂···NH₃ (Z = B and Al). The magnitude of the π-hole on the triel atom is related to the nature of the Z atom and the hybridization of R. CH₃BH₂ has the largest π-hole among RBH₂, while for RAlH₂ the largest π-hole is found in CH≡CAlH₂. The interaction energy is partly inconsistent with the magnitude of the π-hole on the triel atom and the orbital interaction from the N lone pair of NH₃ into the empty p orbital of the triel atom. The strongest B···N triel bond is found in CH≡CBH₂···NH₃, while the weakest Al···N triel bond is in CH₃AlH₂···NH₃. The strength of the triel bond is increased in solvents, and its enhancement is prominent with the increase of solvent polarity. Solvents also change the nature of the Al···N triel bond from an electrostatic interaction to a partially covalent one. The F substituent in the triel donor strengthens the triel bond, depending on the substitution position and number. Graphical Abstract The π-hole triel bonded complexes between RZH2 (Z =B and Al) and NH3 have been investigated. We focused on the effects of hybridization, solvent, and substitution on the strength and nature of π-hole triel bond.

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.

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.

Bond paths between distant atoms do not necessarily indicate dominant interactions

The goal of the article is to revive discussion on the interpretation of bond paths linking distant atoms, particularly tracing weak interactions in dimers. According to the Pendás' concept of privileged exchange channel, a bond path is formed between this pair of competing atoms, which is associated with larger value of the exchange energy. We point out that, due to the short-range nature of the exchange energy, bond paths linking distant atoms clearly become doubtful indicators of dominant intermolecular interactions, particularly if some other characteristics (geometric, spectroscopic, based on electrostatic parameters, etc.) indicate other intermolecular interactions as dominant. Several such cases are thoroughly investigated. We show that electrostatic parameters are much more reliable indicators of dominant intermolecular interactions than bond paths. Then, we pay attention that the presence of ("unexpected", i.e., not necessarily indicating dominant intermolecular interactions) bond paths between pairs of atoms featuring highly expanded charge distributions can be easily explained by visual exploration of isodensity contour plots. As always pointing in the direction of the steepest increase, the gradient vector of the electron density favors areas of its high values gaining higher exchange energy, yet being blind to highly electron deficient areas nearby, which, however, can quite often be involved in dominant intermolecular interactions as strongly suggested by many other bonding analysis. We also suggest that an interatomic component of Hellmann-Feynman force would most likely be the most reliable indicator of attractive or repulsive character of individual interatomic interaction. © 2018 Wiley Periodicals, Inc.