A Simple Model for Halogen Bond Interaction Energies

Solid-State 19F NMR Chemical Shift in Square-Planar Nickel-Fluoride Complexes Linked by Halogen Bonds

The halogen bond (XB) is a highly directional class of noncovalent interactions widely explored by experimental and computational studies. However, the NMR signature of the XB has attracted limited attention. The prediction and analysis of the solid-state NMR (SSNMR) chemical shift tensor provide useful strategies to better understand XB interactions. In this work, we employ a computational protocol for modeling and analyzing the ¹⁹F SSNMR chemical shifts previously measured in a family of square-planar trans Ni^II-L₂-iodoaryl-fluoride (L = PEt₃) complexes capable of forming self-complementary networks held by a NiF···I(C) halogen bond [Thangavadivale, V.; Chem. Sci. 2018, 9, 3767-3781]. To understand how the ¹⁹F NMR resonances of the nickel-bonded fluoride are affected by the XB, we investigate the origin of the shielding in trans-[NiF(2,3,5,6-C₆F₄I)(PEt₃)₂], trans-[NiF(2,3,4,5-C₆F₄I)(PEt₃)₂], and trans-[NiF(C₆F₅)(PEt₃)₂] in the solid state, where a XB is present in the two former systems but not in the last. We perform the ¹⁹F NMR chemical shift calculations both in periodic and molecular models. The results show that the crystal packing has little influence on the NMR signatures of the XB, and the NMR can be modeled successfully with a pair of molecules interacting via the XB. Thus, the observed difference in chemical shift between solid-state and solution NMR can be essentially attributed to the XB interaction. The very high shielding of the fluoride and its driving contributor, the most shielded component of the chemical shift tensor, are well reproduced at the 2c-ZORA level. Analysis of the factors controlling the shielding shows how the highest occupied Ni/F orbitals shield the fluoride in the directions perpendicular to the Ni-F bond and specifically perpendicular to the coordination plane. This shielding arises from the magnetic coupling of the Ni(3d)/F(2p lone pair) orbitals with the vacant σ_Ni-F^* orbital, thereby rationalizing the very highly upfield (shielded) resonance of the component (δ₃₃) along this direction. We show that these features are characteristic of square-planar nickel-fluoride complexes. The deshielding of the fluoride in the halogen-bonded systems is attributed to an increase in the energy gap between the occupied and vacant orbitals that are mostly responsible for the paramagnetic terms, notably along the most shielded direction.

Strength of the [Z-I···Hal]- and [Z-Hal···I]- Halogen Bonds: Electron Density Properties and Halogen Bond Length as Estimators of Interaction Energy

Bond energy is the main characteristic of chemical bonds in general and of non-covalent interactions in particular. Simple methods of express estimates of the interaction energy, Eint, using relationships between Eint and a property which is easily accessible from experiment is of great importance for the characterization of non-covalent interactions. In this work, practically important relationships between Eint and electron density, its Laplacian, curvature, potential, kinetic, and total energy densities at the bond critical point as well as bond length were derived for the structures of the [Z-I···Hal]- and [Z-Hal···I]- types bearing halogen bonds and involving iodine as interacting atom(s) (totally 412 structures). The mean absolute deviations for the correlations found were 2.06-4.76 kcal/mol.

The HOX⋯SO2 (X=F, Cl, Br, I) Binary Complexes: Implications for Atmospheric Chemistry

Sulfur dioxide and hypohalous acids (HOX, X=F, Cl, Br, I) are ubiquitous molecules in the atmosphere that are central to important processes like seasonal ozone depletion, acid rain, and cloud nucleation. We present the first theoretical examination of the HOX⋯SO₂ binary complexes and the associated trends due to halogen substitution. Reliable geometries were optimized at the CCSD(T)/aug-cc-pV(T+d)Z level of theory for HOF and HOCl complexes. The HOBr and HOI complexes were optimized at the CCSD(T)/aug-cc-pV(D+d)Z level of theory with the exception of the Br and I atoms which were modeled with an aug-cc-pwCVDZ-PP pseudopotential. 27 HOX⋯SO₂ complexes were characterized and the focal point method was employed to produce CCSDT(Q)/CBS interaction energies. Natural Bond Orbital analysis and Symmetry Adapted Perturbation Theory were used to classify the nature of each principle interaction. The interaction energies of all HOX⋯SO₂ complexes in this study ranged from 1.35 to 3.81 kcal mol^-1 . The single-interaction hydrogen bonded complexes spanned a range of 2.62 to 3.07 kcal mol^-1 , while the single-interaction halogen bonded complexes were far more sensitive to halogen substitution ranging from 1.35 to 3.06 kcal mol^-1 , indicating that the two types of interactions are extremely competitive for heavier halogens. Our results provide insight into the interactions between HOX and SO₂ which may guide further research of related systems.

Computational Investigations of Dispersion Interactions between Small Molecules and Graphene-like Flakes

We investigate dispersion interactions in a selection of atomic, molecular, and molecule-surface systems, comparing high-level correlated methods with empirically corrected density functional theory (DFT). We assess the efficacy of functionals commonly used for surface-based calculations, with and without the D3 correction of Grimme. We find that the inclusion of the correction is essential to get meaningful results, but there is otherwise little to distinguish between the functionals. We also present coupled-cluster quality interaction curves for H₂, NO₂, H₂O, and Ar interacting with large carbon flakes, acting as models for graphene surfaces, using novel absolutely localized molecular orbital based methods. These calculations demonstrate that the problems with empirically corrected DFT when investigating dispersion appear to compound as the system size increases, with important implications for future computational studies of molecule-surface interactions.

Relativistic Effects on NMR Parameters of Halogen-Bonded Complexes

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.

Relationships between Interaction Energy and Electron Density Properties for Homo Halogen Bonds of the [(A)nY-X···X-Z(B)m] Type (X = Cl, Br, I)

Relationships between interaction energy (Eint) 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 Eint is the kinetic energy density (Gb), which reasonably approximates the whole set of the structures as -Eint = 0.128Gb2 - 0.82Gb + 1.66 (R2 = 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 Eint for the individual series Cl···Cl, Br···Br, and I···I. A number of the Eint(property) correlations are recommended for the practical application in the express estimates of the strength of the homo-halogen bonds.