Halogen Bonding Interactions: Revised Benchmarks and a New Assessment of Exchange vs Dispersion

A Simple Model for Halogen Bond Interaction Energies

Halogen bonds are prevalent in many areas of chemistry, physics, and biology. We present a statistical model for the interaction energies of halogen-bonded systems at equilibrium based on high-accuracy ab initio benchmark calculations for a range of complexes. Remarkably, the resulting model requires only two fitted parameters, X and B—one for each molecule—and optionally the equilibrium separation, R e , between them, taking the simple form E = X B / R e n . For n = 4 , it gives negligible root-mean-squared deviations of 0.14 and 0.28 kcal mol - 1 over separate fitting and validation data sets of 60 and 74 systems, respectively. The simple model is shown to outperform some of the best density functionals for non-covalent interactions, once parameters are available, at essentially zero computational cost. Additionally, we demonstrate how it can be transferred to completely new, much larger complexes and still achieve accuracy within 0.5 kcal mol - 1 . Using a principal component analysis and symmetry-adapted perturbation theory, we further show how the model can be used to predict the physical nature of a halogen bond, providing an efficient way to gain insight into the behavior of halogen-bonded systems. This means that the model can be used to highlight cases where induction or dispersion significantly affect the underlying nature of the interaction.

Is the Fluorine in Molecules Dispersive? Is Molecular Electrostatic Potential a Valid Property to Explore Fluorine-Centered Non-Covalent Interactions?

Can two sites of positive electrostatic potential localized on the outer surfaces of two halogen atoms (and especially fluorine) in different molecular domains attract each other to form a non-covalent engagement? The answer, perhaps counterintuitive, is yes as shown here using the electronic structures and binding energies of the interactions for a series of 22 binary complexes formed between identical or different atomic domains in similar or related halogen-substituted molecules containing fluorine. These were obtained using various computational approaches, including density functional and ab initio first-principles theories with M06-2X, RHF, MP2 and CCSD(T). The physical chemistry of non-covalent bonding interactions in these complexes was explored using both Quantum Theory of Atoms in Molecules and Symmetry Adapted Perturbation Theories. The surface reactivity of the 17 monomers was examined using the Molecular Electrostatic Surface Potential approach. We have demonstrated inter alia that the dispersion term, the significance of which is not always appreciated, which emerges either from an energy decomposition analysis, or from a correlated calculation, plays a structure-determining role, although other contributions arising from electrostatic, exchange-repulsion and polarization effects are also important. The 0.0010 a.u. isodensity envelope, often used for mapping the electrostatic potential is found to provide incorrect information about the complete nature of the surface reactive sites on some of the isolated monomers, and can lead to a misinterpretation of the results obtained.

Halogen-bonded cocrystallization with phosphorus, arsenic and antimony acceptors

The formation of non-covalent directional interactions, such as hydrogen or halogen bonds, is a central concept of materials design, which hinges on using small compact atoms of the 2nd period, notably nitrogen and oxygen, as acceptors. Heavier atoms are much less prominent in that context, and mostly limited to sulfur. Here, we report the experimental observation and theoretical study of halogen bonds to phosphorus, arsenic and antimony in the solid state. Combining 1,3,5-trifluoro-2,4,6-triiodobenzene with triphenylphosphine, -arsine, and -stibine provides cocrystals based on I···P, I···As and I···Sb halogen bonds. The demonstration that increasingly metallic pnictogens form halogen bonds sufficiently strong to enable cocrystal formation is an advance in supramolecular chemistry which opens up opportunities in materials science, as shown by colossal thermal expansion of the cocrystal involving I···Sb halogen bonds.

Halogen bonding can be exploited for the design of functional supramolecular materials, but heavier elements that are known to accept a halogen bond remain limited. Here, the authors demonstrate the formation of two-component cocrystals based on halogen bonds with phosphorus, arsenic and antimony.

A comparison between hydrogen and halogen bonding: the hypohalous acid–water dimers, HOX⋯H2O (X = F, Cl, Br)

Hypohalous acids (HOX) are a class of molecules that play a key role in the atmospheric seasonal depletion of ozone and have the ability to form both hydrogen and halogen bonds.

A DFT assessment of some physical properties of iodine-centered halogen bonding and other non-covalent interactions in some experimentally reported crystal geometries

A set of six binary complexes that feature iodine-centered halogen bonding, extracted from structures deposited in the Cambridge Structure Database, has been examined computationally using density functional theory calculations with the M06-2X global hybrid, and dispersion corrected B3LYP-D3 and B97-D3, to determine their equilibrium geometries, binding energies and electronic properties. The results show that gas phase calculations are very informative in evaluating what occurs in the solid state, even though these calculations ignore the importance of lattice packing and counter ion effects. The calculated binding energies for the non-covalent interactions responsible for these complexes lie between -4.15 and -7.48 kcal mol-1 (M06-2X), which enables us to characterize them as weak-to-moderate in strength. The basis set superposition error energies are calculated to vary between 0.60 and 2.42 kcal mol-1 for all the complexes examined, even though an all-electron QZP basis set used in the analysis was of quadrupole-ζ (plus polarization) quality. Dispersion is found to have a profound effect on the binding energy of some of these complexes, and was estimated to be as large as 5.0 kcal mol-1. For one complex, the crystal geometry could not be precisely reproduced using a gas phase calculation. While both halogen- and hydrogen-bonding interactions were found competitive, they cooperate with each other to determine the stable configuration of the binary complex. The molecular electrostatic surface potential, quantum theory of atoms in molecules, and reduced density gradient non-covalent Interaction models were utilized to arrive at a fundamental understanding of the various inter- and intra-molecular molecular interactions involved, as well as some other previously-overlooked non-covalent interactions that emerge in the modelling.

Revealing Factors Influencing the Fluorine-Centered Non-Covalent Interactions in Some Fluorine-Substituted Molecular Complexes: Insights from First-Principles Studies

We examine the equilibrium structure and properties of six fully or partially fluorinated hydrocarbons and several of their binary complexes using computational methods. In the monomers, the electrostatic surface of the fluorine is predicted to be either entirely negative or weakly positive. However, its lateral sites are always negative. This enables the fluorine to display an anisotropic distribution of charge density on its electrostatic surface. While this is the electrostatic surface scenario of the fluorine atom, its negative sites in some of these monomers are shown to have the potential to engage in attractive engagements with the negative site(s) on the same atom in another molecule of the same type, or a molecule of a different type, to form bimolecular complexes. This is revealed by analyzing the results of current state-of-the-art computational approaches such as DFT, together with those obtained from the quantum theory of atoms in molecules, molecular electrostatic surface potential and symmetry adapted perturbation theories. We demonstrate that the intermolecular interaction energy arising in part from the universal London dispersion, which has been underappreciated for decades, is an essential factor in explaining the attraction between the negative sites, although energy arising from polarization strengthens the extent of the intermolecular interactions in these complexes.