Competition between chalcogen bond and halogen bond interactions in YOX<Subscript>4</Subscript>:NH<Subscript>3</Subscript> (Y = S,Se; X = F,Cl, Br) complexes: An ab initio investigation

Using ab initio calculations, the geometries, interaction energies and bonding properties of chalcogen bond and halogen bond interactions between YOX₄ (Y = S, Se; X = F, Cl, Br) and NH₃ molecules are studied. These binary complexes are formed through the interaction of a positive electrostatic potential region (σ-hole) on the YOX₄ with the negative region in the NH₃. The ab initio calculations are carried out at the MP2/aug-cc-pVTZ level, through analysis of molecular electrostatic potentials, quantum theory of atoms in molecules and natural bond orbital methods. Our results indicate that even though the chalcogen and halogen bonds are mainly dominated by electrostatic effects, but the polarization and dispersion effects also make important contributions to the total interaction energy of these complexes. The examination of interaction energies suggests that the chalcogen bond is always favored over the halogen bond for all of the binary YOX₄:NH₃ complexes.     

The Prominent Enhancing Effect of the Cation–π Interaction on the Halogen–Hydride Halogen Bond in M1⋅⋅⋅C6H5X⋅⋅⋅HM2

We designed M¹???C₆H₅X???HM² (M¹=Li⁺, Na⁺; X=Cl, Br; M²=Li, Na, BeH, MgH) complexes to enhance halogen–hydride halogen bonding with a cation–π interaction. The interaction strength has been estimated mainly in terms of the binding distance and the interaction energy. The results show that halogen–hydride halogen bonding is strengthened greatly by a cation–π interaction. The interaction energy in the triads is two to six times as much as that in the dyads. The largest interaction energy is ?8.31 kcal mol^?1 for the halogen bond in the Li⁺???C₆H₅Br???HNa complex. The nature of the cation, the halogen donor, and the metal hydride influence the nature of the halogen bond. The enhancement effect of Li⁺ on the halogen bond is larger than that of Na⁺. The halogen bond in the Cl donor has a greater enhancement than that in the Br one. The metal hydride imposes its effect in the order HBeH<HMgH<HNa<HLi for the Cl complex and HBeH<HMgH<HLi<HNa for the Br complex. The large cooperative energy indicates that there is a strong interplay between the halogen–hydride halogen bonding and the cation–π interaction. Natural bond orbital and energy decomposition analyses indicate that the electrostatic interaction plays a dominate role in enhancing halogen bonding by a cation–π interaction.     

Theoretical study on intermolecular interactions in BrF/HnX adducts

Equilibrium geometries, interaction energies, and charge transfer for the intermolecular interactions between BrF and H_nX (HF, H₂O, and NH₃) were studied at the MP2/6-311++G(3d,3p) level. The halogen-bonded geometry and hydrogen-bonded geometry are observed in these interactions. The calculated interaction energies show that the halogen-bonded structures are more stable than the corresponding hydrogen-bonded structures. To study the nature of the intermolecular interactions, symmetry-adapted perturbation theory (SAPT) calculations were carried out and the results indicate that the halogen bonding interactions are dominantly inductive energy in nature, while electrostatic energy governs the hydrogen bonding interactions.     

Why are the Elemental Nonmetals (F2, Cl2, Br2, I2, S8, P4) of so Many Hues or of Any Hues and Where is the Chromophore? Insight into Intera-X–X Bonds

A unique approach is used to relate the HOMO-LUMO energy difference to the difference between the ionization potential (IP) and electron affinity (EA) to assist in deducing not only the colors, but also chromophores in elemental nonmetals. Our analysis focuses on compounds with lone pair electrons and σ electrons, namely X₂ (X = F, Cl, Br, I), S₈ and P₄. For the dihalogens, the [IP – EA] energies are found to be: F₂ (12.58 eV), Cl₂ (8.98 eV), Br₂ (7.90 eV), I₂ (6.78 eV). We suggest that the interahalogen X–X bond itself is the chromophore for these dihalogens, in which the light absorbed by the F₂, Cl₂, Br₂, I₂ leads to longer wavelengths in the visible by a π → σ* transition. Trace impurities are a likely case of cyclic S₈ which contains amounts of selenium leading to a yellow color, where the [IP – EA] energy of S₈ is found to be 7.02 eV. Elemental P₄ with an [IP – EA] energy of 9.09 eV contains a tetrahedral and σ aromatic structure. In future work, refinement of the analysis will be required for compounds with π electrons and σ electrons, such as polycyclic aromatic hydrocarbons (PAHs).     

Halogen bonding: a theoretical study based on atomic multipoles derived from quantum theory of atoms in molecules

Atomic multipole moments derived from quantum theory of atoms in molecules are used to study halogen bonds in dihalogens (with general formula YX, in which X refers to the halogen directly interacted with the Lewis base) and some molecules containing C–X group. Multipole expansion is used to calculate the electrostatic potential in a vicinity of halogen atom (which is involved in halogen bonding) in terms of atomic monopole, dipole, and quadrupole moments. In all the cases, the zz component of atomic traceless quadrupole moments (where z axis taken along Y–X or C–X bonds) of the halogens plays a stabilizing role in halogen bond formation. The effects of atomic monopole and dipole moments on the formation of a halogen bond in YX molecules depend on Y and X atoms. In Br₂ and Cl₂, the monopole moment of halogens is zero and has no contribution in electrostatic potential and hence in halogen bonding, while in ClBr, FBr, and FCl it is positive and therefore stabilize the halogen bonds. On the other hand, the negative sign of dipole moment of X in all the YX molecules weakens the corresponding halogen bonds. In the C–X-containing molecules, monopole and dipole moments of X atom are negative and consequently destabilize the halogen bonds. So, in these molecules the quadrupole moment of X atom is the only electrostatic term which strengthens the halogen bonds. In addition, we found good linear correlations between halogen bonds strength and electrostatic potentials calculated from multipole expansion.     

The Role of Molecular Electrostatic Potentials in the Formation of a Halogen Bond in Furan⋅⋅⋅XY and Thiophene⋅⋅⋅XY Complexes

The halogen bonding of furan???XY and thiophene???XY (X=Cl, Br; Y=F, Cl, Br), involving σ‐ and π‐type interactions, was studied by using MP2 calculations and quantum theory of “atoms in molecules” (QTAIM) studies. The negative electrostatic potentials of furan and thiophene, as well as the most positive electrostatic potential (V_S,max) on the surface of the interacting X atom determined the geometries of the complexes. Linear relationships were found between interaction energy and V_S,max of the X atom, indicating that electrostatic interactions play an important role in these halogen‐bonding interactions. The halogen‐bonding interactions in furan???XY and thiophene???XY are weak, “closed‐shell” noncovalent interactions. The linear relationship of topological properties, energy properties, and the integration of interatomic surfaces versus V_S,max of atom X demonstrate the importance of the positive σ hole, as reflected by the computed V_S,max of atom X, in determining the topological properties of the halogen bonds.     

Theoretical Studies on the Halogen Bond between Propargyl Radical and Dihalogen Molecules

MP2/aug‐cc‐pVDZ calculations are carried out on the geometries, vibrational frequencies, interaction energies and topological properties for the π‐type halogen‐bonded complexes between propargyl radical and dihalogen molecules ClF, BrF and BrCl. There are two kinds of geometries: complex ( a ) involves the interaction between the X (X=Cl, Br) atom and the midpoint of C(1) –C(2) bond, complex ( b ) involves the interaction between the X atom and C(3) atom. The lengths of the halogen bond, the frequencies of the halogen bond, the elongation extent of the X–Y (XY=ClF, BrF, BrCl) bond, topological parameters at the BCPs of the halogen bond and X–Y bond are all well consistent with the interaction energies. The interaction of complex ( a ) is stronger than that of complex ( b ); the interaction of propargyl···BrF is stronger than that of propargyl···ClF and propargyl···BrCl. For the complexes ( a ) and ( b ), the charge transfer is observed from propargyl radical to XY, the atomic energy, the dipolar polarization, and the volume of the halogen atom X decrease upon complex formation.     

Cooperativity between the Halogen Bond and the Hydrogen Bond in H3N⋅⋅⋅XY⋅⋅⋅HF Complexes (X,Y=F,Cl, Br)

Ab initio calculations are used to provide information on H₃N???XY???HF triads (X, Y=F, Cl, Br) each having a halogen bond and a hydrogen bond. The investigated triads include H₃N???Br₂‐HF, H₃N???Cl₂???HF, H₃N???BrCI???HF, H₃N???BrF???HF, and H₃N???ClF???HF. To understand the properties of the systems better, the corresponding dyads are also investigated. Molecular geometries, binding energies, and infrared spectra of monomers, dyads, and triads are studied at the MP2 level of theory with the 6‐311++G(d,p) basis set. Because the primary aim of this study is to examine cooperative effects, particular attention is given to parameters such as cooperative energies, many‐body interaction energies, and cooperativity factors. The cooperative energy ranges from ?1.45 to ?4.64 kcal mol^?1, the three‐body interaction energy from ?2.17 to ?6.71 kcal mol^?1, and the cooperativity factor from 1.27 to 4.35. These results indicate significant cooperativity between the halogen and hydrogen bonds in these complexes. This cooperativity is much greater than that between hydrogen bonds. The effect of a halogen bond on a hydrogen bond is more pronounced than that of a hydrogen bond on a halogen bond.     

H2O, H2S与双卤分子间相互作用的电子密度拓扑研究

运用量子化学微扰理论MP2和密度泛函B3LYP方法, 采用6-311＋＋G(d,p)基组, 对H₂O, H₂S与双卤分子XY (XY＝F₂, Cl₂, Br₂, ClF, BrF, BrCl)形成的卤键复合物进行构型全优化, 并计算得到了这些体系的分子间相互作用能. 利用电子密度拓扑分析方法对卤键复合物的拓扑性质进行了分析研究, 探讨了该类分子间卤键的作用本质. 结果表明, 形成卤键后, 作为电子受体的双卤分子X—Y键长增长, 振动频率减小. 复合物体系中的卤键介于共价键与离子键之间, 偏于静电作用成分为主.     

Source of Cooperativity in Halogen‐Bonded Haloamine Tetramers

Inspired by the isostructural motif in α‐bromoacetophenone oxime crystals, we investigated halogen–halogen bonding in haloamine quartets. Our Kohn–Sham molecular orbital and energy decomposition analysis reveal a synergy that can be traced to a charge‐transfer interaction in the halogen‐bonded tetramers. The halogen lone‐pair orbital on one monomer donates electrons into the unoccupied σ*_N?X orbital on the perpendicular N?X bond of the neighboring monomer. This interaction has local σ symmetry. Interestingly, we discovered a second, somewhat weaker donor–acceptor interaction of local π symmetry, which partially counteracts the aforementioned regular σ‐symmetric halogen‐bonding orbital interaction. The halogen–halogen interaction in haloamines is the first known example of a halogen bond in which back donation takes place. We also find that this cooperativity in halogen bonds results from the reduction of the donor–acceptor orbital‐energy gap that occurs every time a monomer is added to the aggregate.     

Halogen Bonding or Hydrogen Bonding between 2,2,6,6-Tetramethyl-piperidine-noxyl Radical and Trihalomethanes CHX3 (X=Cl,Br, I)

The halogen and hydrogen bonding complexes between 2,2,6,6-tetramethylpiperidine-noxyl and trihalomethanes (CHX₃, X=Cl, Br, I) are simulated by computational quantum chem-istry. The molecular electrostatic potentials, geometrical parameters and interaction energy of halogen and hydrogen bonding complexes combined with natural bond orbital analysis are obtained. The results indicate that both halogen and hydrogen bonding interactions obey the order Cl相似文献   

Substituent effects on noncovalent halogen/π interactions: Theoretical study

Noncovalent halogen/π interactions of FCl with substituted benzenes have been investigated using ab initio calculations. It was shown that the predicted maximum interaction energy gap between the substituted and unsubstituted systems amounts to 1.14 kcal/mol, and therefore substituents on benzene have a pronounced effect on the strength of halogen/π interactions. While the presence of electron‐donating groups (NH₂, CH₃, and OH) on benzene enhances the interaction energy appreciably, an opposite effect is observed for electron‐accepting groups (NO₂, CN, Br, Cl, and F). The large gain of the attraction by electron correlation illustrates that the stabilities of the systems considered arise primarily from the dispersion interaction. Beside the dispersion interaction, the charge‐transfer interaction also plays an important role in halogen/π interactions, as a charge density analysis suggested. To provide more insight into the nature of halogen/π interactions, topological analysis of the electron density distribution and properties of bond critical points were determined in terms of the atoms in molecules (AIM) theory. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

The influence of the rate of bond stretching on the differential cross sections for ion-pair formation in K + XY collisions

Experimental differential cross sections are presented for ion-pair formation in collisions between K and I₂, IBr, Br₂, ICl, Cl₂. In this sequence the vibrational frequency of the halogen molecule is increased gradually from one system to the next. Therefore this is a suitable set of cros sections for analysing the influence of the rate at which the halogen bond stretches after transition to the ionic state. The experimental data are compared with trajectory calculations on diabatic and on adiabatic potential energy surfaces. The results of the two types of calculations differ significantly; the adiabatic model gives the better agreement with experiment. The essential difference between the two models is that in the diabatic model stretching of the halogen bond starts only when the crossing seam R_c between the ionic and covalent configurations is reached, while in adiabatic scattering, pre-stretching of the halogen bond occurs as the alkali atom approaches the crossing seam. As a result of pre-stretching the electron affinity is increased, which in turn affects the endoergicity and the total and differential cross sections for ion-pair formation. At energie below the threshold for ion-pair formation, the cross section for reactive scattering is also affected.     

Computer‐Aided Design of Ionic Liquids as CO2 Absorbents

Ionic liquids (ILs), vary strongly in their interaction with CO₂. We suggest simple theoretical approach to predict the CO₂ absorption behavior of ILs. Strong interaction of the CO₂ with the IL anions corresponds to chemical absorption whereas weak interaction indicates physical absorption. A predictive estimate with a clear distinction between physical and chemical absorption can be simply obtained according to geometries optimized in the presence of a solvation model instead of optimizing it only in gas phase as has been done to date. The resulting Gibbs free energies compare very well with experimental values and the energies were correlated with experimental capacities. Promising anions, for ionic liquids with reversible CO₂ absorption properties can be defined by a reaction Gibbs free energy of absorption in the range of ?30 to 16 kJ mol^?1.     

Theoretical study of the interaction of fluorinated dimethyl ethers and the ClF and HF molecules. Comparison between halogen and hydrogen bonds

A theoretical study of the halogen‐bonded complexes formed between fluorinated dimethyl ethers (n_F = 0–4) and ClF is carried out using the wB97XD method combined with the 6‐311++G(d,p) basis set. The properties of the complexes are compared with the corresponding properties of the hydrogen‐bonded complexes formed between the same electron donors and HF. The optimized geometries, the interaction energies, relevant natural bonding orbital characteristics along with some vibrational data are calculated. The analyzed properties also include the symmetry adapted perturbation theory decomposition of the energies along with the atoms‐in molecule analysis. For both the halogen and hydrogen bonds, the interaction energies are ruled by the intermolecular hyperconjugation energies. In contrast, the correlations between the binding energies and the basic properties of the ethers or the charge transfer are different for the halogen and hydrogen bonds. The applicability of the Bent's rule to these systems is discussed. © 2016 Wiley Periodicals, Inc.     

How Dihalogens Catalyze Michael Addition Reactions

We have quantum chemically analyzed the catalytic effect of dihalogen molecules (X₂=F₂, Cl₂, Br₂, and I₂) on the aza‐Michael addition of pyrrolidine and methyl acrylate using relativistic density functional theory and coupled‐cluster theory. Our state‐of‐the‐art computations reveal that activation barriers systematically decrease as one goes to heavier dihalogens, from 9.4 kcal mol^?1 for F₂ to 5.7 kcal mol^?1 for I₂. Activation strain and bonding analyses identify an unexpected physical factor that controls the computed reactivity trends, namely, Pauli repulsion between the nucleophile and Michael acceptor. Thus, dihalogens do not accelerate Michael additions by the commonly accepted mechanism of an enhanced donor–acceptor [HOMO(nucleophile)–LUMO(Michael acceptor)] interaction, but instead through a diminished Pauli repulsion between the lone‐pair of the nucleophile and the Michael acceptor's π‐electron system.     

Molecular Dynamics Simulation Study of ClF in Water: Halogen Bonding Interaction in Liquid

Does the halogen bonding interaction co-exist in liquid when it competes with the hydrogen bonding interaction? The classical molecular dynamics simulations for the solvation properties of ClF molecule in water are performed with the Lennard-Jones plus Coulomb electrostatic potential parameters that are optimized with ab initio interaction energy calculations for the pre-reactive H₂O…ClF complex. We find that the halogen bonding interactions occur between O and Cl atoms and have the comparable strength and population with respect to the hydrogen bonding interactions of Cl…H.     

Comparative Study of Halogen‐ and Hydrogen‐Bond Interactions between Benzene Derivatives and Dimethyl Sulfoxide

The halogen bond, similar to the hydrogen bond, is an important noncovalent interaction and plays important roles in diverse chemistry‐related fields. Herein, bromine‐ and iodine‐based halogen‐bonding interactions between two benzene derivatives (C₆F₅Br and C₆F₅I) and dimethyl sulfoxide (DMSO) are investigated by using IR and NMR spectroscopy and ab initio calculations. The results are compared with those of interactions between C₆F₅Cl/C₆F₅H and DMSO. First, the interaction energy of the hydrogen bond is stronger than those of bromine‐ and chlorine‐based halogen bonds, but weaker than iodine‐based halogen bond. Second, attractive energies depend on 1/rⁿ, in which n is between three and four for both hydrogen and halogen bonds, whereas all repulsive energies are found to depend on 1/r^8.5. Third, the directionality of halogen bonds is greater than that of the hydrogen bond. The bromine‐ and iodine‐based halogen bonds are strict in this regard and the chlorine‐based halogen bond only slightly deviates from 180°. The directional order is iodine‐based halogen bond>bromine‐based halogen bond>chlorine‐based halogen bond>hydrogen bond. Fourth, upon the formation of hydrogen and halogen bonds, charge transfers from DMSO to the hydrogen‐ and halogen‐bond donors. The CH₃ group contributes positively to stabilization of the complexes.     

Halogen Bonding Propensity in Solution: Direct Observation and Computational Prediction

Halogen-bonded complexes are often designed by consideration of electrostatic potential (ESP) predictions. ESP predictions do not capture the myriad variables associated with halogen bond (XB) donors and acceptors; thus, binding interaction cannot be quantitatively predicted. Here, a discrepancy between predictions based on ESP energy difference (ΔV_s) and computed gas phase binding energy (ΔE_bind) motivated the experimental determination of the relative strength of halogen bonding interactions in solution by Raman spectroscopic observation of complexes formed from interacting five iodobenzene-derived XB donors and four pyridine XB acceptors. Evaluation of ΔE_bind coupled with absolutely-localized molecular orbital energy decomposition analysis (ALMO-EDA) deconvolutes halogen bonding energy contributions and reveals a prominent role for charge transfer (CT) interactions. Raman spectra reveal ΔE_bind accurately predicts stronger interactions within iodopentafluorobenzene (IPFB) complexes than with 1-iodo-3,5-dinitrobenzene (IDNB) complexes even though IPFB has similar electrostatics to IDNB and contains a smaller σ-hole.     

Regular/abnormal variation in the strength and nature of the halogen bond between H2Te and the dihalogens: Prominent effect of methyl substituents

The halogen-bonded complexes between H₂Te/Me₂Te and the dihalogen molecules XY (XY = F₂, Cl₂, Br₂, I₂, ClF, ClBr, BrF, BrCl, BrI, IF, ICl, IBr) have been studied to investigate the dependence of its strength and nature on the halogen donor X and its adjoining atom Y, as well as the methyl groups in the electron donor. The interaction energy varies between −1.7 and − 43.5 kcal/mol, indicating that the Te atom in H₂Te/Me₂Te has a strong affinity for the dihalogen molecules. For the H₂Te-XY complex, the halogen bond is stronger for the heavier halogen donor X atom and the strong electron-withdrawing group Y. However, for Me₂Te-XY, the halogen bond is stronger for the lighter halogen donor X atom. The H₂Te/Me₂Te-F₂ complex has the largest interaction energy, although the σ-hole on F₂ is the smallest in magnitude. In most of the complexes, the electrostatic and polarization contributions to the binding strength are similar in magnitude. However, for H₂Te/Me₂Te-F₂, the polarization contribution is much larger than the electrostatic contribution, with a significant contribution from charge transfer.