Assessment of DFT Functionals for QTAIM Topological Analysis of Halogen Bonds with Benzene

Osme Bond: Geometric and Energetic Features in the Adducts between OsO4 and Lewis Bases

Adducts between OsO4 and Lewis bases exert a role in important oxidation processes such as epoxidation and dihydroxylation. It has been shown that the attractive interaction driving the formation of these adducts is a σ-hole bond involving the metal as the electrophilic species; the term Osme Bond (OmB) was proposed for designating it. Here some new adducts between OsO4 and various bases have been characterized through single crystal x-ray diffraction (XRD) and computational studies (density functional theory, DFT), confirming the existence of a robust correlation between σ-hole interaction energy and deformation of the tetrahedral geometry of OsO4. Also, some adducts formed by RuO4 with nucleophiles were investigated computationally.

Substituent-Guided Cluster Nuclearity for Tetranuclear Iron(III) Compounds with Flat {Fe4(μ3-O)2} Butterfly Core

The tetranuclear iron(III) compounds [Fe4(μ3-O)2(μ-LZ)4] (1–3) were obtained by reaction of FeCl3 with the shortened salen-type N2O2 tetradentate Schiff bases N,N’-bis(salicylidene)-o-Z-phenylmethanediamine H2LZ (Z = NO2, Cl and OMe, respectively), where the one-carbon bridge between the two iminic nitrogen donor atoms guide preferentially to the formation of oligonuclear species, and the ortho position of the substituent Z on the central phenyl ring selectively drives towards Fe4 bis-oxido clusters. All compounds show a flat almost-symmetric butterfly-like conformation of the {Fe4(μ3-O)2} core, surrounded by the four Schiff base ligands, as depicted by both the X-ray molecular structures of 1 and 2 and the optimized geometries of all derivatives as obtained by UM06/6-311G(d) DFT calculations. The strength of the antiferromagnetic exchange coupling constants between the iron(III) ions varies among the three derivatives, despite their magnetic cores remain structurally almost unvaried, as well as the coordination of the metal ions, with a distorted octahedral environment for the two-body iron ions, Feb, and a pentacoordination with trigonal bipyramidal geometry for the two-wing iron ions, Few. The different magnetic behavior within the series of examined compounds may be ascribed to the influence of the electronic features of Z on the electron density distribution (EDD) of the central {Fe4(μ3-O)2} core, substantiated by a Quantum Theory of Atoms In Molecules (QTAIM) topological analysis of the EDD, as obtained by UM06 calculations 1–3.

Principles Guiding the Square Bonding Motif Containing a Pair of Chalcogen Bonds between Chalcogenadiazoles

The bonding motif adopted by a dimer of chalcogenadiazole molecules is characterized by a pair of equivalent Ch···N chalcogen bonds. Quantum calculations show that the interaction energy is substantial, varying between 4 kcal/mol for Ch = S and 17 kcal/mol for Te. The interaction is cooperative in that the total bond strength is greater than either chalcogen bond individually. Neither the addition of a phenyl ring nor the addition of a pair of cyano substituents to the diazole ring has much influence on this binding. Removal of one N from the diazole weakens the binding, and addition of two nitrogens has little effect. The largest perturbation arises with three N atoms in each ring, for which the binding energy increases by some 25%. The ring size plays a minor role in most cases, although a near doubling of bond strength occurs if there are two N atoms present on a four-membered ring.

Influence of Substituents in the Benzene Ring on the Halogen Bond of Iodobenzene with Ammonia

The effects on the CI··N halogen bond between iodobenzene and NH3 of 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 NO2 were each placed ortho, meta, and para to the I. The depth of the σ-hole on I is deepened as R became more electron-withdrawing which is reflected in a strengthening of the halogen bond, which varied between 3.3 and 5.5 kcal/mol. 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.

Search for an exothermic halogen bond between anions

A halogen-bonded complex containing a pair of anions can be made more stable than the isolated anions if the Lewis acid is a long carbon chain, fully substituted by CN groups, with an I atom on one end and a COO− group on the other, with Cl− as base.

Three types of noncovalent interactions studied between pyrazine and XF

Three types noncovalent interactions (type I, II and III) between pyrazine (C₄H₄N₂) and XF (X = F, Cl, Br, and I) have been discovered at the MP2/aug-cc-pVTZ level. TypeI is σ-hole interaction between the positive site on the halogen X of XF and the negative site on one of the pyrazine nitrogens. Type II is counterintuitive σ-hole interaction driven by polarization between the positive site on the halogen X of XF and a portion of the pyrazine ring. Type III is an interaction between the lateral regions of the halogen X of XF and the position of the pyrazine ring. Through comparing the calculated interaction energy, we can know that the type II and type III interactions are weaker than the corresponding type I interactions, and type III interactions are weaker than the corresponding type II interactions in C₄H₄N₂-XF complexes. SAPT analysis shows that the electrostatic energy are the major source of the attraction for the type I (σ-hole) interactions while the type III interactions are mainly dispersion energy. For the type II (counterintuitive σ-hole) interactions in C₄H₄N₂-XF (X = F and Cl) complexes, electrostatic energy are the major source of the attraction, while in C₄H₄N₂-XF (X = Br and I) complexes, the electrostatic term, induction and dispersion play equally important role in the total attractive interaction. NBO analysis, AIM theory, and conceptual DFT are also being utilized.

On the Ability of Nitrogen to Serve as an Electron Acceptor in a Pnicogen Bond

Whereas pnicogen atoms like P and As have been shown repeatedly to act as electron acceptors in pnicogen bonds, the same is not true of the more electronegative first-row N atom. Quantum calculations assess whether N can serve in this capacity in such bonds and under what conditions. There is a positive π-hole belt that surrounds the central N atom in the linear arrangement of NNNF, NNN-CN, and NNO, which can engage a NH₃ base to form a pnicogen bond with binding energy between 3 and 5 kcal/mol. Within the context of a planar arrangement, the π-hole above the N in NO₂OF, N(CN)₃, and CF₃NO₂ is also capable of forming a pnicogen bond, the strongest of which amounts to 11 kcal/mol with NMe₃ as base. In their pyramidal geometry, NF₃ and N(NO₂)₃ engage with a base through the σ-hole on the central N, with variable binding energies between 2 and 9 kcal/mol. AIM and NBO provide somewhat different interpretations of the secondary interactions that occur in some of these complexes.

Maximal occupation by bases of π-hole bands surrounding linear molecules

Linear molecules such as CO₂ contain a positive π-hole ring that surrounds C on the molecule's equator. Quantum calculations examine the question as to how many bases can simultaneously bind to this ring. Linear molecules examined are TO₂ , where T = C, Si, Ge, Sn; bases are NCH and NH₃ . CO₂ engages in the weakest of the tetrel bonds, and can bind up to three NCH and two NH₃ . Unlike σ-hole tetrel bonds, Si forms the strongest tetrel bonds, with interaction energies as high as 43 kcal/mol with NH₃ . But like GeO₂ , SiO₂ can sustain only two bases in its equatorial ring. The π-hole ring of SnO₂ can engage in up to four tetrel bonds with either NCH or NH₃ , even though these bonds are weaker than those with GeO₂ or SiO₂ . As all of these complexes cast TO₂ in the role of multiple electron acceptor, the resulting negative cooperativity makes each successive bond weaker than its predecessor as bases are added, as well as reducing the magnitude of the central molecule's π-hole.

Exploring Orthogonality between Halogen and Hydrogen Bonding Involving Benzene

The concept of orthogonality between halogen and hydrogen bonding, brought out by Ho and coworkers some years ago, has become a widely accepted idea within the chemists' community. While the original work was based on a common carbonyl oxygen as acceptor for both interactions, we explore here, by means of M06-2X, M11, ωB97X, and ωB97XD/aug-cc-PVTZ DFT calculations, the interdependence of halogen and hydrogen bonding with a shared π-electron system of benzene. The donor groups (specifically NCBr and H2O) were placed on either or the same side of the ring, according to a double T-shaped or a perpendicular geometry, respectively. The results demonstrate that the two interactions with benzene are not strictly independent on each other, therefore outlining that the orthogonality between halogen and hydrogen bonding, intended as energetical independence between the two interactions, should be carefully evaluated according to the specific acceptor group.

Anatomy of π-hole bonds: Linear systems

The list of σ-hole bonds is long and growing, encompassing both H-bonds and its closely related halogen, chalcogen, etc., sisters. These bonds rely on the asymmetric distribution of electron density, whose depletion along the extension of a covalent bond leaves a positive region of electrostatic potential from which these bonds derive their name. However, the density distributions of other molecules contain analogous positive regions that lie out of the molecular plane known as π-holes, which are likewise capable of engaging in noncovalent bonds. Quantum calculations are applied to study such π-hole bonds that involve linear molecules, whose positive region is a circular belt surrounding the molecule, rather than the more restricted area of a σ-hole. These bonds are examined in terms of their most fundamental elements arising from the spatial dispositions of their relevant molecular orbitals and the π-holes in both the total electron density and the electrostatic potential to which they lead. Systems examined comprise tetrel, chalcogen, aerogen, and triel bonds, as well as those involving group II elements, with atoms drawn from various rows of the Periodic Table. The π-hole bonds established by linear molecules tend to be weaker than those of comparable planar systems.

Dissection of the Origin of π-Holes and the Noncovalent Bonds in Which They Engage

Accompanying the rapidly growing list of σ-hole bonds has come the acknowledgment of parallel sorts of noncovalent bonds which owe their stability in large part to a deficiency of electron density in the area above the molecular plane, known as a π-hole. The origins of these π-holes are probed for a wide series of molecules, comprising halogen, chalcogen, pnicogen, tetrel, aerogen, and spodium bonds. Much like in the case of their σ-hole counterparts, formation of the internal covalent π-bond in the Lewis acid molecule pulls density toward the bond midpoint and away from its extremities. This depletion of density above the central atom is amplified by an electron-withdrawing substituent. At the same time, the amplitude of the π*-orbital is enhanced in the region of the density-depleted π-hole, facilitating a better overlap with the nucleophile's lone pair orbital and a stabilizing n → π* charge transfer. The presence of lone pairs on the central atom acts to attenuate the π-hole and shift its position somewhat, resulting in an overall weakening of the π-hole bond. There is a tendency for π-hole bonds to include a higher fraction of induction energy than σ-bonds with proportionately smaller electrostatic and dispersion components, but this distinction is less a product of the σ- or π-character and more a function of the overall bond strength.

Proximity Effects of Substituents on Halogen Bond Strength

It is well known that the presence of an electron-withdrawing substituent (EWS) placed near the halogen (X) atom on a Lewis acid molecule amplifies the ability of this unit to engage in a halogen bond with a base. Quantum calculations are applied to examine how quickly these effects fade as the EWS is moved further and further from the X atom. Conjugated alkene and alkyne chains of varying lengths with a terminal C-I first facilitate analysis as to how the number of these multiple bonds affects the strength of CI··N XB to NH₃. Then, electron-withdrawing F and C≡N substituents are placed on the opposite end of the chain, and their effects on the XB properties are monitored as a function of their distance from I. These same EWSs are added to the ortho, meta, and para positions of aromatic iodobenzene. It is found that the XB grows in strength as more triple bonds are added to the alkyne, but there is little change caused by elongating an alkene. The cyano group has a much stronger effect than does F. While F strengthens the XB, its effects are quickly attenuated as it is moved further from I. The consequences of C≡N substitution are stronger and extend over a longer distance. Placement of an EWS on the phenyl ring diminishes with distance: o > m > p, and the effects of disubstitution are nearly additive. These trends apply not only to energetics but also to geometries, properties of the wave function, σ-hole depth, and NMR shielding.

Relative Strengths of a Pnicogen and a Tetrel Bond and Their Mutual Effects upon One Another

The ability of the T and Z atoms of TR₃ZR₂ to engage in a noncovalent interaction with NH₃ is assessed by DFT calculations, where the T atom refers to C, Si, and Ge; Z = As, Sb, and P; and substituents R = H and F. In most instances, the tetrel bond (TB) is both stronger and shorter than the pnicogen bond (ZB). These two bond strengths can be equalized, or preference shifted to the ZB, if F substituents are placed on the Z and H on the T atoms. Employing C as the T atom results in a very weak TB, with the ZB clearly favored energetically. The simultaneous formation of both TB and ZB weakens both, particularly the latter, but both bonds survive intact. Geometric and spectroscopic perturbations of the subunits reflect the two types of noncovalent bonds.

Halogen Complexes of Anionic N-Heterocyclic Carbenes

The lithium complexes [(WCA-NHC)Li(toluene)] of anionic N-heterocyclic carbenes with a weakly coordinating anionic borate moiety (WCA-NHC) reacted with iodine, bromine, or CCl4 to afford the zwitterionic 2-halogenoimidazolium borates (WCA-NHC)X (X=I, Br, Cl; WCA=B(C6 F5 )3 , B{3,5-C6 H3 (CF3 )2 }3 ; NHC=IDipp=1,3-bis(2,6-diisopropylphenyl)imidazolin-2-ylidene, or NHC=IMes=1,3-bis(2,4,6-trimethylphenyl)imidazolin-2-ylidene). The iodine derivative (WCA-IDipp)I (WCA=B(C6 F5 )3 ) formed several complexes of the type (WCA-IDipp)I⋅L (L=C6 H5 Cl, C6 H5 Me, CH3 CN, THF, ONMe3 ), revealing its ability to act as an efficient halogen bond donor, which was also exploited for the preparation of hypervalent bis(carbene)iodine(I) complexes of the type [(WCA-IDipp)I(NHC)] and [PPh4 ][(WCA-IDipp)I(WCA-NHC)] (NHC=IDipp, IMes). The corresponding bromine complex [PPh4 ][(WCA-IDipp)2 Br] was isolated as a rare example of a hypervalent (10-Br-2) system. DFT calculations reveal that London dispersion contributes significantly to the stability of the bis(carbene)halogen(I) complexes, and the bonding was further analyzed by quantum theory of atoms in molecules (QTAIM) analysis.

Polycentric binding in complexes of trimethylamine-N-oxide with dihalogens

Dihalogens readily interact with trimethylamine-N-oxide under ambient conditions. Accordingly, herein, stable 1 : 1 adducts were obtained in the case of iodine chloride and iodine bromide. The crystal and molecular structure of the trimethylamine-N-oxide-iodine chloride adduct was solved. Furthermore, the geometry and electronic structure of the trimethylamine-N-oxide-dihalogen complexes were studied computationally. Only molecular ensembles were found in the global minimum for the 1 : 1 stoichiometry. The O⋯X-Y halogen bond is the main factor for the thermodynamic stability of these complexes. Arguments for electrostatic interactions as the driving force for this noncovalent interaction were discussed. Also, the equilibrium structures are additionally stabilised by weak C-H⋯X hydrogen bonds. Consequently, formally monodentate ligands are bound in a polycentric manner.

Comparison of Bifurcated Halogen with Hydrogen Bonds

Bifurcated halogen bonds are constructed with FBr and FI as Lewis acids, paired with NH₃ and NCH bases. The first type considered places two bases together with a single acid, while the reverse case of two acids sharing a single base constitutes the second type. These bifurcated systems are compared with the analogous H-bonds wherein FH serves as the acid. In most cases, a bifurcated system is energetically inferior to a single linear bond. There is a larger energetic cost to forcing the single σ-hole of an acid to interact with a pair of bases, than the other way around where two acids engage with the lone pair of a single base. In comparison to FBr and FI, the H-bonding FH acid is better able to participate in a bifurcated sharing with two bases. This behavior is traced to the properties of the monomers, in particular the specific shape of the molecular electrostatic potential, the anisotropy of the orbitals of the acid and base that interact directly with one another, and the angular extent of the total electron density of the two molecules.

Ag(i) and Cu(i) cyclic-triimidazole coordination polymers: revealing different deactivation channels for multiple room temperature phosphorescences

The photophysics of isostructural Ag(i) and Cu(i) 1D and 3D coordination polymers based on cyclic triimidazole reveals excitation wavelength-dependent multiple emissions with different radiative paths according to the coordinated metal.

Halogen Bond of Halonium Ions: Benchmarking DFT Methods for the Description of NMR Chemical Shifts

Because of their anisotropic electron distribution and electron deficiency, halonium ions are unusually strong halogen-bond donors that form strong and directional three-center, four-electron halogen bonds. These halogen bonds have received considerable attention owing to their applicability in supramolecular and synthetic chemistry and have been intensely studied using spectroscopic and crystallographic techniques over the past decade. Their computational treatment faces different challenges to those of conventional weak and neutral halogen bonds. Literature studies have used a variety of wave functions and DFT functionals for prediction of their geometries and NMR chemical shifts, however, without any systematic evaluation of the accuracy of these methods being available. In order to provide guidance for future studies, we present the assessment of the accuracy of 12 common DFT functionals along with the Hartree-Fock (HF) and the second-order Møller-Plesset perturbation theory (MP2) methods, selected from an initial set of 36 prescreened functionals, for the prediction of 1H, 13C, and 15N NMR chemical shifts of [N-X-N]+ halogen-bond complexes, where X = F, Cl, Br, and I. Using a benchmark set of 14 complexes, providing 170 high-quality experimental chemical shifts, we show that the choice of the DFT functional is more important than that of the basis set. The M06 functional in combination with the aug-cc-pVTZ basis set is demonstrated to provide the overall most accurate NMR chemical shifts, whereas LC-ωPBE, ωB97X-D, LC-TPSS, CAM-B3LYP, and B3LYP to show acceptable performance. Our results are expected to provide a guideline to facilitate future developments and applications of the [N-X-N]+ halogen bond.

Coordination of a Central Atom by Multiple Intramolecular Pnicogen Bonds

A central pnicogen Z atom (Z = Sb, As) is covalently attached to the O atom of three -O(CH₂)_nX chains where X represents either an aldehyde or amine group. The chain can fold around so that the basic X group can engage in a noncovalent pnicogen bond with the central Z. The formation of up to three pnicogen bonds is energetically favored. The amine appears to engage in stronger pnicogen bonds than does the aldehyde, and bonds to Sb are favored over As, but there is little dependence on the length of the chain. The formation of each successive pnicogen bond reduces the magnitude of the σ-holes surrounding the Z atom, which tends to weaken the attraction for the basic end of the chain.

Halogen bonding (HaB) in E–I⋯X–M systems: influence of the halogen donor on the HaB nature

A study on how the covalent character of the halogen bond can be modulated by using different halogen donors.

The ditetrel bond: noncovalent bond between neutral tetrel atoms

The ability of a tetrel atom to serve in the capacity of electron donor in a σ-hole noncovalent bond is tested by quantum calculations.

The mechanism of monochloramine disproportionation under acidic conditions

Monochloramine is a widely employed agent in water treatment technologies. However, its utilization has some drawbacks like the transformation of the active species into the undesired dichloramine. Although it is more pronounced in acidic solutions, the features of this reaction have still remained largely unexplored in the pH < 4 region. In this study the decomposition of monochloramine is examined under such conditions by using kinetic and computational methods. Fast kinetics measurements have convincingly showed that the disproportion into dicloramine is relatively fast and can be studied without any interference from side reactions. By varying the pH, the deprotonation constant of monochloramine has been determined by UV spectroscopy (K_a = 0.023 ± 0.005 M for I = 1.0 M NaClO₄, and T = 25.0 °C). Dichloramine formation via monochloramine disproportion was found to follow second-order kinetics. The computations have provided the reaction mechanism and its free energy profile in accord with the proposed kinetic model. This involves the reaction between the protonated and unprotonated forms of monochloramine, with a rate constant k = 335.3 ± 11.8 M^-1 s^-1, corresponding to an activation free energy barrier of 14.1 kcal mol^-1. The simulations predicted a barrier of 14.9 kcal mol^-1 and revealed a key short-lived chlorine-bridged intermediate which yields dichloroamine and ammonium ion through a deprotonation-coupled chlorine shift.

Effects of Halogen, Chalcogen, Pnicogen, and Tetrel Bonds on IR and NMR Spectra

Complexes were formed pairing FX, FHY, FH2Z, and FH3T (X = Cl, Br, I; Y = S, Se, Te; Z = P, As, Sb; T = Si, Ge, Sn) with NH3 in order to form an A⋯N noncovalent bond, where A refers to the central atom. Geometries, energetics, atomic charges, and spectroscopic characteristics of these complexes were evaluated via DFT calculations. In all cases, the A-F bond, which is located opposite the base and is responsible for the σ-hole on the A atom, elongates and its stretching frequency undergoes a shift to the red. This shift varies from 42 to 175 cm-1 and is largest for the halogen bonds, followed by chalcogen, tetrel, and then pnicogen. The shift also decreases as the central A atom is enlarged. The NMR chemical shielding of the A atom is increased while that of the F and electron donor N atom are lowered. Unlike the IR frequency shifts, it is the third-row A atoms that undergo the largest change in NMR shielding. The change in shielding of A is highly variable, ranging from negligible for FSnH3 all the way up to 1675 ppm for FBr, while those of the F atom lie in the 55-422 ppm range. Although smaller in magnitude, the changes in the N shielding are still easily detectable, between 7 and 27 ppm.

Differential Binding of Tetrel-Bonding Bipodal Receptors to Monatomic and Polyatomic Anions

Previous work has demonstrated that a bidentate receptor containing a pair of Sn atoms can engage in very strong interactions with halide ions via tetrel bonds. The question that is addressed here concerns the possibility that a receptor of this type might be designed that would preferentially bind a polyatomic over a monatomic anion since the former might better span the distance between the two Sn atoms. The binding of Cl⁻ was thus compared to that of HCOO⁻, HSO₄⁻, and H₂PO₄⁻ with a wide variety of bidentate receptors. A pair of SnFH₂ groups, as strong tetrel-binding agents, were first added to a phenyl ring in ortho, meta, and para arrangements. These same groups were also added in 1,3 and 1,4 positions of an aliphatic cyclohexyl ring. The tetrel-bonding groups were placed at the termini of (-C≡C-)_n (n = 1,2) extending arms so as to further separate the two Sn atoms. Finally, the Sn atoms were incorporated directly into an eight-membered ring, rather than as appendages. The ordering of the binding energetics follows the HCO₂⁻ > Cl⁻ > H₂PO₄⁻ > HSO₄⁻ general pattern, with some variations in selected systems. The tetrel bonding is strong enough that in most cases, it engenders internal deformations within the receptors that allow them to engage in bidentate bonding, even for the monatomic chloride, which mutes any effects of a long Sn···Sn distance within the receptor.

Halogen Bonds Formed between Substituted Imidazoliums and N Bases of Varying N-Hybridization

Heterodimers are constructed containing imidazolium and its halogen-substituted derivatives as Lewis acid. N in its sp³, sp² and sp hybridizations is taken as the electron-donating base. The halogen bond is strengthened in the Cl < Br < I order, with the H-bond generally similar in magnitude to the Br-bond. Methyl substitution on the N electron donor enhances the binding energy. Very little perturbation arises if the imidazolium is attached to a phenyl ring. The energetics are not sensitive to the hybridization of the N atom. More regular patterns appear in the individual phenomena. Charge transfer diminishes uniformly on going from amine to imine to nitrile, a pattern that is echoed by the elongation of the C-Z (Z=H, Cl, Br, I) bond in the Lewis acid. These trends are also evident in the Atoms in Molecules topography of the electron density. Molecular electrostatic potentials are not entirely consistent with energetics. Although I of the Lewis acid engages in a stronger bond than does H, it is the potential of the latter which is much more positive. The minimum on the potential of the base is most negative for the nitrile even though acetonitrile does not form the strongest bonds. Placing the systems in dichloromethane solvent reduces the binding energies but leaves intact most of the trends observed in vacuo; the same can be said of ∆G in solution.