skip to main content


Title: Characterizing the interplay of Pauli repulsion, electrostatics, dispersion and charge transfer in halogen bonding with energy decomposition analysis
The halogen bond is a class of non-covalent interaction that has attracted considerable attention recently. A widespread theory for describing them is the σ-hole concept, which predicts that the strength of the interaction is proportional to the size of the σ-hole, a region of positive electrostatic potential opposite a σ bond. Previous work shows that in the case of CX 3 I, with X equal to F, Cl, Br, and I, the σ-hole trend is exactly opposite to the trend in binding energy with common electron pair donors. Using energy decomposition analysis (EDA) applied to a potential energy scan as well as the recent adiabatic EDA technique, we show that the observed trend is a result of charge transfer. Therefore a picture of the halogen bond that excludes charge transfer cannot be complete, and permanent and induced electrostatics do not always provide the dominant stabilizing contributions to halogen bonds. Overall, three universally attractive factors, polarization, dispersion and charge transfer, together with permanent electrostatics, which is usually attractive, drive halogen bonding, against Pauli repulsion.  more » « less
Award ID(s):
1665315 1363342
NSF-PAR ID:
10066954
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
20
Issue:
2
ISSN:
1463-9076
Page Range / eLocation ID:
905 to 915
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The halogen bond formed by a series of Lewis acids TF 3 X (T = C, Si, Ge, Sn, Pb; X = Cl, Br, I) with NH 3 is studied by quantum chemical calculations. The interaction energy is closely mimicked by the depth of the σ-hole on the X atom as well as the full electrostatic energy. There is a first trend by which the hole is deepened if the T atom to which X is attached becomes more electron-withdrawing: C > Si > Ge > Sn > Pb. On the other hand, larger more polarizable T atoms are better able to transmit the electron-withdrawing power of the F substituents. The combination of these two opposing factors leaves PbF 3 X forming the strongest XBs, followed by CF 3 X, with SiF 3 X engaging in the weakest bonds. The charge transfer from the NH 3 lone pair into the σ*(TX) antibonding orbital tends to elongate the covalent TX bond, and this force is largest for the heavier X and T atoms. On the other hand, the contraction of this bond deepens the σ-hole at the X atom, which would enhance both the electrostatic component and the full interaction energy. This bond-shortening effect is greatest for the lighter X atoms. The combination of these two opposing forces leaves the T–X bond contracting for X = Cl and Br, but lengthening for I. 
    more » « less
  2. Abstract

    The replacement of a CH group of benzene by a triel (Tr) atom places a positive region of electrostatic potential near the Tr atom in the plane of the aromatic ring. This σ‐hole can interact with an X lone pair of XCCH (X=F, Cl, Br, and I) to form a triel bond (TrB). The interaction energy between C5H5Tr and FCCH lies in the range between 2.2 and 4.4 kcal/mol, in the order Tr=B+cation above the ring pulls density toward itself and thus magnifies the Tr σ‐hole. The TrB to the XCCH nucleophile is thereby magnified as is the strength of the TrB. This positive cooperativity is particularly large for Tr=B.

     
    more » « less
  3. Abstract

    The effects on the C−I⋅⋅N halogen bond between iodobenzene and NH3of placing various substituents on the phenyl ring are monitored by quantum calculations. Substituents R=N(CH3)2, NH2, CH3, OCH3, COCH3, Cl, F, COH, CN, and NO2were each placed ortho, meta, and para to the I. The depth of the σ‐hole on I is deepened as R becomes more electron‐withdrawing which is reflected in a strengthening of the halogen bond, which varied between 3.3 and 5.5 kcal mol−1. In most cases, the ortho placement yields the largest perturbation, followed by meta and then para, but this trend is not universal. Parallel to these substituent effects is a progressive lengthening of the covalent C−I bond. Formation of the halogen bond reduces the NMR chemical shielding of all three nuclei directly involved in the C−I⋅⋅N interaction. The deshielding of the electron donor N is most closely correlated with the strength of the bond, as is the coupling constant between I and N, so both have potential use as spectroscopic measures of halogen bond strength.

     
    more » « less
  4. null (Ed.)
    Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential. 
    more » « less
  5. The relationship between solid-state supramolecular interactions and crystal habits is highlighted based on experimental and computational analysis of the crystal structure of strong halogen-bonded (HaB) associations between iodine-containing dihalogens (ICl, IBr) with 1,4-diazabicyclo[2,2,2]octane (DABCO) as well as with substituted pyridines and phenazine. The pattern of the energy frameworks and the interplay of the attractive and repulsive interactions in the solid-state associations involving these HaB donors and acceptors directly correlated with their crystal habits. This correlation suggests that analysis of the energy framework serves as a useful tool (complementary to the earlier developed methods) to rationalize and predict the crystal habit. The X-ray structural analysis also revealed that the I⋯N distances in the complexes were in the 2.24–2.54 Å range, i.e. they were much closer to the I⋯N covalent bond length than to the van der Waals separation. The computational analysis of the nature of halogen bonding in these complexes showed delocalization of their molecular orbitals' between donor and acceptors resulting in a substantial charge transfer from the nucleophiles to dihalogens and elongation of the I⋯X bond. As a result, both I⋯N and I⋯X bonds in the strongest complexes ( e.g. , ICl with DABCO or 4-dimethylaminopyridine) are characterized by the comparable Mayer bonds orders of about 0.6, along with the electron and energy densities at their bond critical points of about 0.1 a.u. and −0.02 a.u., respectively. These data as well as the density overlap regions indicator (DORI) point to the covalency of the I⋯N bonding and suggest that the interaction within the IX complexes can be described as (unsymmetrical) hypervalent 3c/4e N⋯I⋯X bonding akin to that in trihalide or halonium ions. 
    more » « less