Halogen bonding and beyond: factors influencing the nature of CN–R and SiN–R complexes with F–Cl and Cl2

Quasi-atomic orbital analysis of halogen bonding interactions

A quasi-atomic orbital analysis of the halogen bonded NH3⋯XF complexes (X = F, Cl, Br, and I) is performed to gain insight into the electronic properties associated with these σ-hole interactions. It is shown that significant sharing of electrons between the nitrogen lone pair of the ammonia molecule and the XF molecule occurs, resulting in a weakening of the X-F bond. In addition, the N-X bond shows increasing covalent character as the size of the halogen atom X increases. While the Mulliken outer complex NH3⋯XF appears to be overall the main species, the strength of the covalent interaction of the N-X bond becomes increasingly similar to that of the N-X bond in the [NH3X]+ cation as the size of X increases.

Ability of Peripheral H Bonds to Strengthen a Halogen Bond

Quantum calculations study the manner in which the involvement of a halogen atom as a proton acceptor in one or more H bonds (HBs) affects the strength of the halogen bond (XB) it can form with a nucleophile aligned with the X σ-hole. A variety of Lewis acids wherein X = F, Cl, Br, and I are attached to a tetrel atom C or Ge engaged in a XB with nucleophile NH₃. One, two, and three HF molecules were positioned perpendicular to the XB axis so that they could form a HB to the X atom. Each such HB strengthened the XB by an increment of 1 kcal/mol or more that does not attenuate as each new HB is added, potentially increasing the interaction energy manyfold. Additionally, the presence of one or more HBs facilitates the formation of a XB by molecules which are reluctant to engage in such a bond in the absence of these auxiliary interactions. Even the F atom, which avoids such a XB, can be coaxed to participate in a XB of moderate strength by one or more of these external HBs.

Interaction and Polarization Energy Relationships in σ-Hole and π-Hole Bonding

We demonstrate that a wide range of σ- and π-hole interaction energies can be related to (a) the electrostatic potentials and electric fields of the σ- and π-hole molecules at the approximate positions of the negative sites and (b) the electrostatic potentials and polarizabilities of the latter. This is consistent with the Coulombic nature of these interactions, which should be understood to include both electrostatics and polarization. The energies associated with polarization were estimated and were shown to overall be greater for the stronger interactions; no new factors need be introduced to account for these. All of the interactions can be treated in the same manner.

Extended Assemblies of Ru(bpy)(CO)2X2 (X = Cl, Br, I) Molecules Linked by 1,4-Diiodotetrafluoro-Benzene (DITFB) Halogen Bond Donors

The ruthenium carbonyl compounds, Ru(bpy)(CO)2X2 (X = Cl, Br or I) act as neutral halogen bond (XB) acceptors when co-crystallized with 1,4-diiodotetrafluoro-benzene (DITFB). The halogen bonding strength of the Ru-X⋅⋅⋅I halogen bonds follow the nucleophilic character of the halido ligand. The strongest halogen bond occurs between the chlorido ligand and the iodide atoms of the DITFB. All three halogen bonded complexes form polymeric assemblies in the solid state. In Ru(bpy)(CO)2Cl2⋅DITFB (1) and in Ru(bpy)(CO)2Br2⋅DITFB (2) both halido ligands are halogen bonded to only one DITFB donor. In Ru(bpy)(CO)2I2⋅DITFB (3) only one of the halido ligands is involved in halogen bonding acting as ditopic center for two DITFB donors. The polymeric structures of 1 and 2 are isomorphic wave-like single chain systems, while the iodine complexes form pairs of linear chains attached together with weak F⋅⋅⋅O≡C interactions between the closest neighbors. The stronger polarization of the iodide ligand compared to the Cl or Br ligands favors nearly linear C-I⋅⋅⋅I angles between the XB donor and the metal complex supporting the linear arrangement of the halogen bonded chain.

An Overview of Strengths and Directionalities of Noncovalent Interactions: σ-Holes and π-Holes

Quantum mechanics, through the Hellmann–Feynman theorem and the Schrödinger equation, show that noncovalent interactions are classically Coulombic in nature, which includes polarization as well as electrostatics. In the great majority of these interactions, the positive electrostatic potentials result from regions of low electronic density. These regions are of two types, designated as σ-holes and π-holes. They differ in directionality; in general, σ-holes are along the extensions of covalent bonds to atoms (or occasionally between such extensions), while π-holes are perpendicular to planar portions of molecules. The magnitudes and locations of the most positive electrostatic potentials associated with σ-holes and π-holes are often approximate guides to the strengths and directions of interactions with negative sites but should be used cautiously for this purpose since polarization is not being taken into account. Since these maximum positive potentials may not be in the immediate proximities of atoms, interatomic close contacts are not always reliable indicators of noncovalent interactions. This is demonstrated for some heterocyclic rings and cyclic polyketones. We briefly mention some problems associated with using Periodic Table Groups to label interactions resulting from σ-holes and π-holes; for example, the labels do not distinguish between these two possibilities with differing directionalities.

The thermal expansion properties of halogen bond containing 1,4 dioxane halogen complexes

Strong halogen bonds formed between 1,4 dioxane and dihalogens lead to minimum expansion in the direction of these bonds.

Influence of monomer deformation on the competition between two types of σ-holes in tetrel bonds

Competition between two competing sites on a tetrel atom is explained by balance between structural deformation and σ-hole intensity.

A perspective on quantum mechanics and chemical concepts in describing noncovalent interactions

Since quantum mechanical calculations do not typically lend themselves to chemical interpretation, analyses of bonding interactions depend largely upon models (the octet rule, resonance theory, charge transfer, etc.). This sometimes leads to a blurring of the distinction between mathematical modelling and physical reality. The issue of polarization vs. charge transfer is an example; energy decomposition analysis is another. The Hellmann-Feynman theorem at least partially bridges the gap between quantum mechanics and conceptual chemistry. It proceeds rigorously from the Schrödinger equation to demonstrating that the forces exerted upon the nuclei in molecules, complexes, etc., are entirely classically coulombic attractions with the electrons and repulsions with the other nuclei. In this paper, we discuss these issues in the context of noncovalent interactions. These can be fully explained in coulombic terms, electrostatics and polarization (which include electronic correlation and dispersion).

Computational study of the NO, SO2, and NH3 adsorptions on fragments of 3N-graphene and Al/3N graphene

The adsorption properties of common gas molecules (NO, NH₃, and SO₂) on the surface of 3N-graphene and Al/3N graphene fragments are investigated using density functional theory. The adsorption energies have been calculated for the most stable configurations of the molecules on the surface of 3N-graphene and Al/3N graphene fragments. The adsorption energies of Al/3N graphene-gas systems are -220.5 kJ mol^-1 for Al/3NG-NO, -111.9 kJ mol^-1 for Al/3NG-NH₃, and -347.7 kJ mol^-1 for Al/3NG-SO₂, respectively. Compared with the 3N-graphene fragment, the Al/3N graphene fragment has significant adsorption energy. Furthermore, the molecular orbital, density of states, and electron densities distribution were used to explore the interaction between these molecules and the surface. We found that orbital hybridization exists between these molecules and the Al/3N graphene surface, which indicates that doping Al significantly increases the interaction between the gas molecules and Al/3N graphene. In addition, compared with Li, Al can more powerfully enhance adsorption of the 3N-graphene fragment. The results indicate that Al/3N graphene can be viewed as a new nanomaterial adsorbent for NO, NH₃, and SO₂.

What kind of neutral halogen bonds can be modulated by solvent effects?

Halogen bonds with a large portion of polarization can be modulated by solvent effects.

Theoretical study of the CO, NO, and N2 adsorptions on Li-decorated graphene and boron-doped graphene

The adsorption properties of common gas molecules (CO, NO, and N2) on the surface of Li-decorated pristine graphene and Li-decorated boron doped graphene are investigated using density functional theory. The adsorption energy, charge transfer, and density of states of gas molecules on three surfaces have been calculated and discussed, respectively. The results show that Li-decorated pristine graphene has strong interaction with CO and N2. Compared with Li-decorated pristine graphene, Li-decorated boron doped graphene exhibit a comparable adsorption ability of CO and N2. Moreover, Li-decorated boron doped graphene have a more significant adsorption energy to NO than that of Li-decorated pristine graphene because of the chemical interaction of the NO gas molecule. The strong interaction between the NO molecule and substrate (Li-decorated boron doped graphene) induces dramatic changes to the electrical conductivity of Li-decorated boron doped graphene. The results indicate that Li-decorated boron doped graphene would be an excellent candidate for sensing NO gas.

Halogen Bonding Involving CO and CS with Carbon as the Electron Donor

MP2/aug'-cc-pVTZ calculations have been carried out to investigate the halogen-bonded complexes formed when CO and CS act as electron-pair donors through C to ClF, ClNC, ClCl, ClOH, ClCN, ClCCH, and ClNH₂. CO forms only complexes stabilized by traditional halogen bonds, and all ClY molecules form traditional halogen-bonded complexes with SC, except ClF which forms only an ion-pair complex. Ion-pair complexes are also found on the SC:ClNC and SC:ClCl surfaces. SC:ClY complexes stabilized by traditional halogen bonds have greater binding energies than the corresponding OC:ClY complexes. The largest binding energies are found for the ion-pair SC-Cl⁺:-Y complexes. The transition structures which connect the complex and the ion pair on SC:ClNC and SC:ClCl potential surfaces provide the barriers for inter-converting these structures. Charge-transfer from the lone pair on C to the σ-hole on Cl is the primary charge-transfer interaction stabilizing OC:ClY and SC:ClY complexes with traditional halogen bonds. A secondary charge-transfer occurs from the lone pairs on Cl to the in-plane and out-of-plane π antibonding orbitals of ClY. This secondary interaction assumes increased importance in the SC:ClNH₂ complex, and is a factor leading to its unusual structure. C-O and C-S stretching frequencies and 13C chemical shieldings increase upon complex formation with ClY molecules. These two spectroscopic properties clearly differentiate between SC:ClY complexes and SC-Cl⁺:-Y ion pairs. Spin-spin coupling constants 1xJ(C-Cl) for OC:ClY complexes increase with decreasing distance. As a function of the C-Cl distance, 1xJ(C-Cl) and ¹J(C-Cl) provide a fingerprint of the evolution of the halogen bond from a traditional halogen bond in the complexes, to a chlorine-shared halogen bond in the transition structures, to a covalent bond in the ion pairs.

Comparison of Various Means of Evaluating Molecular Electrostatic Potentials for Noncovalent Interactions

The various heterodimers formed by a series of Lewis acids with NH₃ as Lewis base are identified. Lewis acids include those that can form chalcogen (HSF and HSBr), pnicogen (H₂ PF and H₂ PBr), and tetrel (H₃ SiF and H₃ SiBr) bonds, as well as H-bonds and halogen bonds. The molecular electrostatic potential (MEP) of each Lewis acid is considered in a number of ways. Pictorial versions show broad regions of positive and negative MEP, on surfaces that vary with respect to either the value of the chosen isopotential, or their distance from the nuclei. Specific points are identified where the MEP reaches a maximum on a particular isodensity surface (V_s,max ). The locations and values of V_s,max were evaluated on different isodensity surfaces, and compared to the stabilities of the various equilibrium geometries. As the chosen isodensity is decreased, and the MEP maxima drift away from the molecule, some points maintain their angular positions with respect to the molecule, whereas others undergo a reorientation. The lowering isodensity also causes some of the maxima to disappear. In general, there is a fairly good correlation between the energetic ordering of the equilibrium structures and the values of V_s,max . A number of possible Lewis acid sites on the heteroaromatic imidazole ring were also considered and presents some cautions about application of V_s,max as the principal criterion for predicting equilibrium geometries. © 2017 Wiley Periodicals, Inc.

σ-holes and π-holes: Similarities and differences

σ-Holes and π-holes are regions of molecules with electronic densities lower than their surroundings. There are often positive electrostatic potentials associated with them. Through these potentials, the molecule can interact attractively with negative sites, such as lone pairs, π electrons, and anions. Such noncovalent interactions, "σ-hole bonding" and "π-hole bonding," are increasingly recognized as being important in a number of different areas. In this article, we discuss and compare the natures and characteristics of σ-holes and π-holes, and factors that influence the strengths and locations of the resulting electrostatic potentials. © 2017 Wiley Periodicals, Inc.

Theoretical insights into the π-hole interactions in the complexes containing triphosphorus hydride (P3H3) and its derivatives

The π-hole of triphosphorus hydride (P₃H₃) and its derivatives Z₃X₃ (Z = P, As; X = H, F, Cl, Br) was discovered and analyzed. MP2/aug-cc-pVDZ calculations were performed on the π-hole interactions in the HCN...Z₃X₃ complexes and the mutual influence between π-hole interactions and the hydrogen bond in the HCN...HCN...Z₃X₃ and HCN...Z₃X₃...HCN complexes studied. The π-hole interaction belongs to the typical closed-shell noncovalent interaction. The linear relationship was found between the most positive electrostatic potential of the π-hole (V_S,max) and the interaction energy. Moreover, the V_S,max of the π-hole was also found to be linearly correlated to the electrostatic energy term, indicating the important contribution of the electrostatic energy term to the π-hole interaction. There is positive cooperativity between the π-hole interaction and the hydrogen bond in the termolecular complexes. The π-hole interaction has a greater influence on the hydrogen bond than vice versa. The mutual enhancing effect between the π-hole interaction and the hydrogen bond in the HCN...HCN...Z₃X₃ complexes is greater than that in the HCN...Z₃X₃...HCN complexes.

Ion-Pair Halogen Bonds in 2-Halo-Functionalized Imidazolium Chloride Receptors: Substituent and Solvent Effects

The interaction of 2-halo-functionalized imidazolium derivatives (n-X⁺ ; X=Cl, Br, I) with a chloride anion through ion-pair halogen bonds (n-X⋅Cl) was studied by means of DFT and ab initio calculations. A method benchmark was performed on 2-bromo-1H-imidazol-3-ium in association with chloride (1-Br⋅Cl); MP2 yielded the best results when compared with CCSD(T) calculations. The interaction energies (ΔE) in the gas phase are high and, although the electrostatic interaction is strong owing to the ion-pair nature of the system, large X⋅⋅⋅Cl^- Wiberg bond orders and contributions from charge transfer (nCl- →σ*C-X) are obtained. These values drop considerably in chloroform and water; this shows that solvent plays a role in modulating the interaction and that gas-phase calculations are particularly unrealistic for experimental applications. The introduction of electron-withdrawing groups in the 4,5-positions of the imidazolium (e.g., -NO₂ , -F) increases the halogen-bond strength in both the gas phase and solvent, including water. The effect of the substituents on the 1,3-positions (N-H groups) also depends on the solvent. The variation of ΔE can be predicted through a two-parameter linear regression that optimizes the weights of charge-transfer and electrostatic interactions, which are different in vacuum and in solvent (chloroform and water). These results could be used in the rational design of efficient chloride receptors based on halogen bonds that work in solution, in particular, in an aqueous environment.

Structural organization in the trimethylamine iodine monochloride complex

The combination of a strong N⋯I–Cl halogen bond and a weak C–H⋯Cl hydrogen bond lead to the formation of dimeric species in the solid.

Quantitative Assessment of Halogen Bonding Utilizing Vibrational Spectroscopy

A total of 202 halogen-bonded complexes have been studied using a dual-level approach: ωB97XD/aug-cc-pVTZ was used to determine geometries, natural bond order charges, charge transfer, dipole moments, electron and energy density distributions, vibrational frequencies, local stretching force constants, and relative bond strength orders n. The accuracy of these calculations was checked for a subset of complexes at the CCSD(T)/aug-cc-pVTZ level of theory. Apart from this, all binding energies were verified at the CCSD(T) level. A total of 10 different electronic effects have been identified that contribute to halogen bonding and explain the variation in its intrinsic strength. Strong halogen bonds are found for systems with three-center-four-electron (3c-4e) bonding such as chlorine donors in interaction with substituted phosphines. If halogen bonding is supported by hydrogen bonding, genuine 3c-4e bonding can be realized. Perfluorinated diiodobenzenes form relatively strong halogen bonds with alkylamines as they gain stability due to increased electrostatic interactions.

The protonated 2-halogenated imidazolium cation as the noncovalent interaction donor: the σ-hole and π-hole interactions

The σ-hole and π-hole of the protonated 2-halogenated imidazolium cation (XC₃H₄N₂⁺; X = F, Cl, Br, I) were investigated and analyzed. The monomers of (CH₃)₃SiY(Y=F, Cl, Br, I), considered as the Lewis base, were combined with the σ-hole and π-hole of XC₃H₄N₂⁺ to form the σ-hole and π-hole interactions in the bimolecular complexes (CH₃)₃SiY · · · XC₃H₄N₂⁺ and (CH₃)₃SiY · · · C₃(X)H₄N₂⁺(X/Y=F, Cl, Br, I), respectively. For both the σ-hole and π-hole interactions, the equilibrium geometries of complexes show regular changes according to the sequence of heavy sequence of the noncovalent interaction acceptors and donors. The electrostatic energy is the main contribution in the formation of both kinds of interactions, it has linear relations with the V _S,max values of σ-hole and the V' _S,max values of π-hole. Both the σ-hole and π-hole interactions belong to the closed-shell and noncovalent interactions. The π-hole interactions are stronger than the σ-hole interactions. For the π-hole interactions, the contribution percents of the dispersion energies are somewhat greater than those of the σ-hole interactions, while it is contrary for the polarization energy. Graphical Abstract The protonated 2-halogenated imidazolium cation as the noncovalent interaction donor: the σ-hole and π-hole interactionsᅟ.

Theoretical study of the structures and first hyperpolarizabilities of C60Cl n and Li@C60Cl n (n = 4, 6, 8, 10)

We recently reported (Song Y-D et al., 2016, J Mol Model 22:50) that doping with Li greatly affects the static first hyperpolarizability of C60Cl2. In this work, with a view to creating nonlinear optical materials with high thermodynamic stability and wide transparent regions, a series of Li@C60Cl n (n = 4, 6, 8, 10) were designed. The structures, electrostatic potentials, electronic structures, absorption spectra, and linear and nonlinear optical properties of C60Cl n and Li@C60Cl n were systematically investigated using density functional theory (DFT) methods. The results of our calculations indicated that the stability of these molecules decreases in the order Li@C60Cl10 > Li@C60Cl8 > Li@C60Cl6 > Li@C60Cl4. It is clear that the number of Cl atoms greatly influences the stability of Li@C60Cl n . Li@C60Cl n showed greater thermodynamic stability than Li@C60Cl2. We also investigated the first hyperpolarizabilities of Li@C60Cl n and found them to be 2973, 3640, 4307, and 2627 au for n = 4, 6, 8, and 10, respectively-higher than that of Li@C60Cl2. Finally, we noted that the transparent region could be modulated by increasing the number of Cl atoms: Li@C60Cl n possess wider transparent regions than that of Li@C60Cl2. We therefore believe that this study has highlighted an effective method for designing nonlinear optical materials with high thermodynamic stability and wide transparent regions.

Comparative studies on group III σ-hole and π-hole interactions

The σ-hole of M2 H6 (M = Al, Ga, In) and π-hole of MH3 (M = Al, Ga, In) were discovered and analyzed, the bimolecular complexes M2 H6 ···NH3 and MH3 ···N2 P2 F4 (M = Al, Ga, In) were constructed to carry out comparative studies on the group III σ-hole interactions and π-hole interactions. The two types of interactions are all partial-covalent interactions; the π-hole interactions are stronger than σ-hole interactions. The electrostatic energy is the largest contribution for forming the σ-hole and π-hole interaction, the polarization energy is also an important factor to form the M···N interaction. The electrostatic energy contributions to the interaction energy of the σ-hole interactions are somewhat greater than those of the π-hole interactions. However, the polarization contributions for the π-hole interactions are somewhat greater than those for the σ-hole interactions.

The encapsulated lithium effect on the first hyperpolarizability of C60Cl2 and C60F2

In this paper, we report a study on the structure and first hyperpolarizability of C60Cl2 and C60F2. The calculation results show that the first hyperpolarizabilities of C60Cl2 and C60F2 were 172 au and 249 au, respectively. Compared with the fullerenes, the first hyperpolarizability of C60Cl2 increased from 0 au to 172 au, while the first hyperpolarizability of C60F2 increased from 0 au to 249 au. In order to further increase the first hyperpolarizability of C60Cl2 and C60F2, Li@C60Cl2 and Li@C60F2 were obtained by introducing a lithium atom to C60Cl2 and C60F2. The first hyperpolarizabilities of Li@C60Cl2 and Li@C60F2 were 2589 au and 985 au, representing a 15-fold and 3.9-fold increase, respectively, over those of C60Cl2 and C60F2. The transition energies of four molecules (C60Cl2, Li@C60Cl2, C60F2, Li@C60F2) were calculated, and were found to be 0.17866 au, 0.05229 au, 0.18385 au, and 0.05212 au, respectively. A two-level model explains why the first hyperpolarizability increases for Li@C60Cl2 and Li@C60F2.

Influence of the protonation of pyridine nitrogen on pnicogen bonding: competition and cooperativity

Ab initio MP2/aug-cc-pVTZ calculations were performed to investigate the pnicogen-bonded complexes of PyZX₂ (Py = pyridine, Z = P and As, X = H and F) and their protonated analogues.

Computational evidence that hyperconjugative orbital interactions are responsible for the stability of intramolecular Te⋯O/Te⋯S non-covalent interactions and comparable to hydrogen bonds in quasi-cyclic systems

The intramolecular secondary bonding interactions involving quasi-cyclic tellurium are comparable to H-bond strength and partially governed by orbital interactions.