Energy profiles of seven halogen-bonded complexes were analysed with the topological energy partitioning called Interacting Quantum Atoms (IQA) at MP4(SDQ)/6–31 + G(2d,2p) level of theory. Explicit interatomic electron correlation energies are included in the analysis. Four complexes combine X₂ (X = Cl or F) with HCN or NH₃, while the remaining three combine ClF with HCN, NH₃ or N₂. Each complex was systematically deformed by translating the constituent molecules along its central axis linking X and N, and reoptimising its remaining geometry. The Relative Energy Gradient (REG) method (Theor. Chem. Acc. 2017, 136, 86) then computes which IQA energies most correlate with the total energy during the process of complex formation and further compression beyond the respective equilibrium geometries. It turns out that the covalent energy (i.e., exchange) of the halogen bond, X…N, itself drives the complex formation. When the complexes are compressed from their equilibrium to shorter X…N distance then the intra-atomic energy of N is in charge. When the REG analysis is restricted to electron correlation then the interatomic correlation energy between X and N again drives the complex formation, and the complex compression is best described by the destabilisation of the through-space correlation energy between N and the “outer” halogen. Full article

In an effort to explain the experimentally observed variation of the photocatalytic activity of ${}^t$ Bu ${}_3$ P, ${}^n$ Bu ${}_3$ P and (MeO) ${}_3$ P in the blue-light regime [Helmecke et al., Org. Lett. 21 (2019) 7823], we have explored the absorption characteristics of several phosphine– and phosphite–IC ${}_4$ F ${}_9$ adducts by means of relativistic density functional theory and multireference configuration interaction methods. Based on the results of these computational and complementary experimental studies, we offer an explanation for the broad tailing of the absorption of ${}^t$ Bu ${}_3$ P-IC ${}_4$ F ${}_9$ and (MeO) ${}_3$ P-IC ${}_4$ F ${}_9$ into the visible-light region. Larger coordinate displacements of the ground and excited singlet potential energy wells in ${}^n$ Bu ${}_3$ P-IC ${}_4$ F ${}_9$ , in particular with regard to the P–I–C bending angle, reduce the Franck–Condon factors and thus the absorption probability compared to ${}^t$ Bu ${}_3$ P-IC ${}_4$ F ${}_9$ . Spectroscopic and computational evaluation of conformationally flexible and locked phosphites suggests that the reactivity of (MeO) ${}_3$ P may be the result of oxygen lone-pair participation and concomitant broadening of absorption. The proposed mechanism for the phosphine-catalyzed homolytic C–I cleavage of perfluorobutane iodide involves S₁ ← S₀ absorption of the adduct followed by intersystem crossing to the photochemically active T ${}_1$ state. Full article

Periodic local vibrational modes were calculated with the rev-vdW-DF2 density functional to quantify the intrinsic strength of the X-I⋯OA-type halogen bonding (X = I or Cl; OA: carbonyl, ether and N-oxide groups) in 32 model systems originating from 20 molecular crystals. We found that the halogen bonding between the donor dihalogen X-I and the wide collection of acceptor molecules OA features considerable variations of the local stretching force constants (0.1–0.8 mdyn/Å) for I⋯O halogen bonds, demonstrating its powerful tunability in bond strength. Strong correlations between bond length and local stretching force constant were observed in crystals for both the donor X-I bonds and I⋯O halogen bonds, extending for the first time the generalized Badger’s rule to crystals. It is demonstrated that the halogen atom X controlling the electrostatic attraction between the $\sigma$ -hole on atom I and the acceptor atom O dominates the intrinsic strength of I⋯O halogen bonds. Different oxygen-containing acceptor molecules OA and even subtle changes induced by substituents can tweak the $n\to {\sigma }^{\ast }$ (X-I) charge transfer character, which is the second important factor determining the I⋯O bond strength. In addition, the presence of the second halogen bond with atom X of the donor X-I bond in crystals can substantially weaken the target I⋯O halogen bond. In summary, this study performing the in situ measurement of halogen bonding strength in crystalline structures demonstrates the vast potential of the periodic local vibrational mode theory for characterizing and understanding non-covalent interactions in materials. Full article

A computational study of the iodide discrimination by different neutral and cationic iodotriazole halogen bonding hosts was carried out by means of Density Functional Theory. The importance of the size of the scaffold was highlighted and its impact observed in the binding energies and intermolecular X⋯I distances. Larger scaffolds were found to reduce the electronic repulsion and increase the overlap between the halide electron lone pair and the corresponding I-C antibonding orbital, increasing the halogen bonding interactions. Additionally, the planarity plays an important role within the interaction, and can be tuned using hydroxyl to perform intramolecular hydrogen bonds (IMHB) between the scaffold and the halogen atoms. Structures with IMHB exhibit stronger halogen bond interactions, as evidenced by the shorter intramolecular distances, larger electron density values at the bond critical point and more negative binding energies. Full article

In the present work, a number of R–X⋯NH₃ (X = Cl, Br, and I) halogen bonded systems were theoretical studied by means of DFT calculations performed at the ωB97XD/6-31+G(d,p) level of theory in order to get insights on the effect of the electron-donating or electron-withdrawing character of the different R substituent groups (R = halogen, methyl, partially fluorinated methyl, perfluoro-methyl, ethyl, vinyl, and acetyl) on the stability of the halogen bond. The results indicate that the relative stability of the halogen bond follows the Cl < Br < I trend considering the same R substituent whereas the more electron-withdrawing character of the R substituent the more stable the halogen bond. Refinement of the latter results, performed at the MP2/6-31+G(d,p) level showed that the DFT and the MP2 binding energies correlate remarkably well, suggesting that the Grimme’s type dispersion-corrected functional produces reasonable structural and energetic features of halogen bond systems. DFT results were also observed to agree with more refined calculations performed at the CCSD(T) level. In a further stage, a more thorough analysis of the R–Br⋯NH₃ complexes was performed by means of a novel electron localization/delocalization tool, defined in terms of an Information Theory, IT, based quantity obtained from the conditional pair density. For the latter, our in-house developed C++/CUDA program, called KLD (acronym of Kullback–Leibler divergence), was employed. KLD results mapped onto the one-electron density plotted at a 0.04 a.u. isovalue, showed that (i) as expected, the localized electron depletion of the Br sigma-hole is largely affected by the electron-withdrawing character of the R substituent group and (ii) the R–X bond is significantly polarized due to the presence of the NH₃ molecule in the complexes. The afore-mentioned constitutes a clear indication of the dominant character of electrostatics on the stabilization of halogen bonds in agreement with a number of studies reported in the main literature. Finally, the cooperative effects on the [Br—CN]_n system (n = 1–8) was evaluated at the MP2/6-31+G(d,p) level, where it was observed that an increase of about ~14.2% on the complex stability is obtained when going from n = 2 to n = 8. The latter results were corroborated by the analysis of the changes on the Fermi-hole localization pattern on the halogen bond zones, which suggests an also important contribution of the electron correlation in the stabilization of these systems. Full article

Relativistic effects are found to be important for the estimation of NMR parameters in halogen-bonded complexes, mainly when they involve the heavier elements, iodine and astatine. A detailed study of 60 binary complexes formed between dihalogen molecules (XY with X, Y = F, Cl, Br, I and At) and four Lewis bases (NH₃, H₂O, PH₃ and SH₂) was carried out at the MP2/aug-cc-pVTZ/aug-cc-pVTZ-PP computational level to show the extent of these effects. The NMR parameters (shielding and nuclear quadrupolar coupling constants) were computed using the relativistic Hamiltonian ZORA and compared to the values obtained with a non-relativistic Hamiltonian. The results show a mixture of the importance of the relativistic corrections as both the size of the halogen atom and the proximity of this atom to the basic site of the Lewis base increase. Full article

Next-generation quantum theory of atoms in molecules (QTAIM) was used to investigate the competition between hydrogen bonding and halogen bonding for the recently proposed (Y = Br, I, At)/halogenabenzene/NH₃ complex. Differences between using the SR-ZORA Hamiltonian and effective core potentials (ECPs) to account for relativistic effects with increased atomic mass demonstrated that next-generation QTAIM is a much more responsive tool than conventional QTAIM. Subtle details of the competition between halogen bonding and hydrogen bonding were observed, indicating a mixed chemical character shown in the 3-D paths constructed from the bond-path framework set B. In addition, the use of SR-ZORA reduced or entirely removed spurious features of B on the site of the halogen atoms. Full article

Relationships between interaction energy (E_int) and electron density properties at the X···X bond critical point or the d(X···X) distance were established for the large set of structures [(A)_nY–X···X–Z(B)_m] bearing the halogen bonds Cl···Cl, Br···Br, and I···I (640 structures in total). The best estimator of E_int is the kinetic energy density (G_b), which reasonably approximates the whole set of the structures as −E_int = 0.128G_b² − 0.82G_b + 1.66 (R² = 0.91, mean absolute deviation 0.39 kcal/mol) and demonstrates low dispersion. The potential and kinetic energy densities, electron density, and the d(X···X) distance behave similarly as estimators of E_int for the individual series Cl···Cl, Br···Br, and I···I. A number of the E_int(property) correlations are recommended for the practical application in the express estimates of the strength of the homo-halogen bonds. Full article

Iodothyronine deiodinases (Dios) are involved in the regioselective removal of iodine from thyroid hormones (THs). Deiodination is essential to maintain TH homeostasis, and disruption can have detrimental effects. Halogen bonding (XB) to the selenium of the selenocysteine (Sec) residue in the Dio active site has been proposed to contribute to the mechanism for iodine removal. Polybrominated diphenyl ethers (PBDEs) and polychlorinated biphenyls (PCBs) are known disruptors of various pathways of the endocrine system. Experimental evidence shows PBDEs and their hydroxylated metabolites (OH-BDEs) can inhibit Dio, while data regarding PCB inhibition are limited. These xenobiotics could inhibit Dio activity by competitively binding to the active site Sec through XB to prevent deiodination. XB interactions calculated using density functional theory (DFT) of THs, PBDEs, and PCBs to a methyl selenolate (MeSe⁻) arrange XB strengths in the order THs > PBDEs > PCBs in agreement with known XB trends. THs have the lowest energy C–X*-type unoccupied orbitals and overlap with the Se lp donor leads to high donor-acceptor energies and the greatest activation of the C–X bond. The higher energy C–Br* and C–Cl* orbitals similarly result in weaker donor-acceptor complexes and less activation of the C–X bond. Comparison of the I···Se interactions for the TH group suggest that a threshold XB strength may be required for dehalogenation. Only highly brominated PBDEs have binding energies in the same range as THs, suggesting that these compounds may inhibit Dio and undergo debromination. While these small models provide insight on the I···Se XB interaction itself, interactions with other active site residues are governed by regioselective preferences observed in Dios. Full article

The strong, specific, and directional halogen bond (XB) is an ideal supramolecular synthon in crystal engineering, as well as rational catalyst and drug design. These attributes attracted strong growing interest in halogen bonding in the past decade and led to a wide range of applications in materials, biological, and catalysis applications. Recently, various research groups exploited the XB mode of activation in designing halogen-based Lewis acids in effecting organic transformation, and there is continual growth in this promising area. In addition to the rapid advancements in methodology development, computational investigations are well suited for mechanistic understanding, rational XB catalyst design, and the study of intermediates that are unstable when observed experimentally. In this review, we highlight recent computational studies of XB organocatalytic reactions, which provide valuable insights into the XB mode of activation, competing reaction pathways, effects of solvent and counterions, and design of novel XB catalysts. Full article